全球长时间序列中尺度气旋数据(http://cioss.coas.oregonstate.edu/eddies/)有netCDF格式,散点数据类型,只有一个很大的维Nobs = 2590938。尝试读取了部分经纬度和track数据,并绘制散度图。
脚本程序:
fn = 'D:/Temp/nc/tracks.20130125.nc' f = addfile(fn) n = 1000 lon = f['lon'][0:n] lat = f['lat'][0:n] track = f['track'][0:n] axesm() mlayer = shaperead('D:/Temp/map/country1.shp') geoshow(mlayer, edgecolor=(0,0,255)) layer = scatterm(lon, lat, track) title('Eddies') colorbar(layer, orientation='horizontal') axism([-180, 180, -90, 90])