标签:
这里演示从micaps第一类数据(地面全要素观测)中读取一个变量(用DimDataFile类的stationdata方法),然后maskout掉中国区域之外的数据,利用scatterm函数绘制散点图。
脚本程序:
f = addfile_micaps(‘D:/Temp/micaps/10101414.000‘) pr = f.stationdata(‘Visibility‘) layer = shaperead(‘D:/Temp/map/china.shp‘) pr = pr.maskout(layer) axesm() mlayer = shaperead(‘D:/Temp/map/country1.shp‘) geoshow(mlayer, edgecolor=(0,0,255)) layer = scatterm(pr) title(‘Visibility‘) xlim(70, 140) ylim(15, 55) colorbar(layer, orientation=‘horizontal‘)
标签:
原文地址:http://www.cnblogs.com/yaqiang/p/4602131.html