登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
天气现象符号分布图实际就是散点图,可以用scatterm函数绘制,但之前需要创建天气符号图例,用weatherspec函数。如果只需要绘制某些天气现象(比如雾、霾),可以给出相应的天气符号序号列表(可以在网上查找)。
脚本程序:
- #Set data folders
- basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
- datadir = os.path.join(basedir, 'sample/MICAPS')
- mapdir = os.path.join(basedir, 'map')
- #Read shape files
- lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
- lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
- #Read station data
- f = addfile_micaps(os.path.join(datadir, '10101414.000'))
- data = f.stationdata('WeatherNow')
- #Plot
- axesm(bgcolor=(204,255,255))
- geoshow(lworld, edgecolor='k', facecolor=(255,251,195))
- geoshow(lchina, edgecolor='k')
- ls = weatherspec()
- layer = scatterm(data, symbolspec=ls)
- #colorbar(layer)
- xlim(72, 136)
- ylim(16, 55)
|