MeteoInfoLab脚本示例:中尺度气旋散点图

全球长时间序列中尺度气旋数据(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])

原文地址:https://www.cnblogs.com/yaqiang/p/4625452.html