Meteoinfolab脚本经验分享：读取轨迹数据（聚类图层）叠加在地形图层

︶ㄣ安定■丶

Meteoinfolab脚本经验分享：读取轨迹数据（聚类图层）叠加在地形图层+添加聚类图层legend（1）根据Meteoinfo软件算出郑州（随便选的地方）2013年12月的轨迹，通过软件聚类生成cluster图层
（2）根据Meteoinfolab脚本编辑器读取轨迹图层+聚类图层
（3）叠加地形图层（地形数据来自于: 之前分享的中国地形数据:http://bbs.06climate.com/forum.php?mod=viewthread&tid=49302)
具体脚本如下：
[

1. ]#add topographic data
   
2. fn = 'H:/Meteoinfo/dixing/dixing-china.nc'
   
3. f1 = addfile(fn)
   
4. tdata1 = f1['z'][:,:]
   
5. 
   
6. #plot
   
7. axesm(tickfontsize=18)
   
8. lchina = shaperead('D:/Temp/map/bou2_4p.shp')
   
9. cluster = shaperead(r'H:\test\zhengzhou-test\2013-12\tgs\cluster\cluster.shp')
   
10. tra = shaperead(r'H:\test\zhengzhou-test\2013-12\tgs\2013-12-zhengzhou.shp')
    
11. geoshow(lchina,edgecolor='k',size=0.1)
    
12. ss = makesymbolspec('line', {'value':'1', 'color':'r','size':3,'caption':'30.65%'  }, \
    
13.     {'value':'2', 'color':'g', 'size':3,'caption':'33.87%'},\
    
14.     {'value':'3', 'color':'y', 'size':3,'caption':'19.35%'},\
    
15.     {'value':'4', 'color':'b', 'size':3,'caption':'16.13%'},field='Cluster'
    
16.     )
    
17. geoshow(tra,color='k')
    
18. geoshow(cluster,fontsize=17,bold='True',yoffset=13, avoidcoll=False,symbolspec=ss)
    
19. legend(legend=cluster.legend(),loc='upper right',fontsize=18,markerscale=2)
    
20. text(90.5, 48, 'Cluster (Zhengzhou)', fontname='Arial', fontsize=18, bold=True)
    
21. 
    
22. #add dixing
    
23. levels = arange(0, 6000, 100)
    
24. layer = imshowm(tdata1,levels,cmap='MPL_terrain')
    
25. xlim(90,120)
    
26. ylim(32,50)
    
27. xlabel('Longitude',fontsize=18)
    
28. ylabel('Latitude',fontsize=18)
    
29. #add a new colorbar
    
30. axes([0.16,0.12,0.7,0.8], newaxes=True, frameon=False, axis=False)
    
31. colorbar(layer,shrink=1,label='Height(m)',fontsize=18)

]#add topographic data<br /> fn = 'H:/Meteoinfo/dixing/dixing-china.nc'<br /> f1 = addfile(fn)<br /> tdata1 = f1['z'][:,:]<br /> <br /> #plot<br /> axesm(tickfontsize=18)<br /> lchina = shaperead('D:/Temp/map/bou2_4p.shp')<br /> cluster = shaperead(r'H:\test\zhengzhou-test\2013-12\tgs\cluster\cluster.shp')<br /> tra = shaperead(r'H:\test\zhengzhou-test\2013-12\tgs\2013-12-zhengzhou.shp')<br /> geoshow(lchina,edgecolor='k',size=0.1)<br /> ss = makesymbolspec('line', {'value':'1', 'color':'r','size':3,'caption':'30.65%'  }, \<br />     {'value':'2', 'color':'g', 'size':3,'caption':'33.87%'},\ <br />     {'value':'3', 'color':'y', 'size':3,'caption':'19.35%'},\<br />     {'value':'4', 'color':'b', 'size':3,'caption':'16.13%'},field='Cluster'<br />     )<br /> geoshow(tra,color='k')<br /> geoshow(cluster,fontsize=17,bold='True',yoffset=13, avoidcoll=False,symbolspec=ss)<br /> legend(legend=cluster.legend(),loc='upper right',fontsize=18,markerscale=2)<br /> text(90.5, 48, 'Cluster (Zhengzhou)', fontname='Arial', fontsize=18, bold=True)<br /> <br /> #add dixing<br /> levels = arange(0, 6000, 100)<br /> layer = imshowm(tdata1,levels,cmap='MPL_terrain')<br /> xlim(90,120)<br /> ylim(32,50)<br /> xlabel('Longitude',fontsize=18)<br /> ylabel('Latitude',fontsize=18)<br /> #add a new colorbar<br /> axes([0.16,0.12,0.7,0.8], newaxes=True, frameon=False, axis=False)<br /> colorbar(layer,shrink=1,label='Height(m)',fontsize=18)

出图效果如下：


轨迹加聚类文件： tgs.rar (126.32 KB, 下载次数: 26)

2017-3-24 10:29 上传

点击文件名下载附件
PS：Meteoinfolab越来越强大了。。强烈推荐大家使用

goddess豆豆

优秀

lgz

优秀

黑色青春

桂博

梁YX

{:eb502:}

liuxiunuli

虽然没看懂,点赞

静言_GRMC

解决了问题

天边的彩虹

好厉害，我想问这个地形图的色标怎么修改？

︶ㄣ安定■丶

天边的彩虹 发表于 2019-1-21 15:19
好厉害，我想问这个地形图的色标怎么修改？

修改imshowm函数里面的cmap参数，把色标改为你需要的就行。

︶ㄣ安定■丶

天边的彩虹 发表于 2019-1-21 15:19
好厉害，我想问这个地形图的色标怎么修改？

修改imshowm函数里面的cmap参数，把色标改为你需要的就行