登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
MeteoInfoLab的georead函数提供了读取shape文件、image文件(JPG、PNG等,需要有相应的地理定位文件)文件生成图层的功能(事实上shaperead也是同样的功能,不过函数名太过具体)。对于Polygon图层可以设置其颜色透明度,可以用makecolors函数生成颜色列表,函数中设置参数alpha(255是不透明,0是完全透明)。在contourfm函数中设置edgecolor参数可以绘制颜色区域的边框。
脚本程序:
- f = addfile('D:/Temp/GrADS/model.ctl')
- ps = f['PS'][0,[10,60],[70,140]]
- lrelief = georead('D:/Temp/Map/GLOBALeb3colshade.jpg')
- world = georead('D:/Temp/Map/country1.shp')
- axesm()
- geoshow(lrelief)
- geoshow(world)
- cols = makecolors(20, alpha=128)
- layer = contourfm(ps, colors=cols, edgecolor='gray', order=1)
- title('Pressure')
- colorbar(layer, aspect=12)
|