这里演示从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