登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 MeteoInfo 于 2018-12-25 12:28 编辑
数据来源见此贴:http://bbs.06climate.com/forum.php?mod=viewthread&tid=38961,需要说明的此青藏高原的范围只包含中国国内的青藏高原部分。
脚本程序:
- f = addfile('D:/Temp/nc/elev.0.25-deg.nc')
- elev = f['data'][0,'25:42','72:105']
- axesm()
- ltibet = shaperead('D:/GeoData/DBATP/DBATP_polygon.shp')
- geoshow('country', edgecolor='k')
- geoshow('cn_province', edgecolor='b', size=2)
- geoshow(ltibet, edgecolor='r', size=2)
- #layer = imshowm(elev, 20)
- layer = contourfm(elev, 20, cmap='MPL_terrain')
- title('Tibet Plateau - China')
- colorbar(layer, orientation='horizontal')
|