netCDF格式是很灵活的,程序自动读出数据的空间维需要数据遵循相应的约定(convensions),MeteoInfo桌面软件目前支持CF convensions,不过主要针对格点数据。这个是站点数据,有空了我看看这个站点数据是不是也遵循了CF convensions,如果是的话可以考虑在MeteoInfo桌面软件中支持。
目前可以用MeteoInfoLab脚本来处理,比如读取站点经纬度和能见度数组并绘图:
 - fn = 'D:/Temp/nc/Surface_Synoptic_20170221_0000.nc'
- f = addfile(fn)
- lon = f['Lon'][:]
- lat = f['Lat'][:]
- stname = f['stnName'][:,:]
- vis = f['VIS'][:]
- vis[vis<0] = nan
- #Plot
- axesm()
- lworld = shaperead('D:/Temp/map/country1.shp')
- geoshow(lworld, edgecolor='k')
- levs = arange(2, 50, 2)
- layer = scatterm(lon, lat, vis, levs, edge=False, size=3)
- colorbar(layer)
- title('Synoptic stations visibility observation')
|