MeteoInfoLab脚本示例：站点数据绘制等值线

MeteoInfo

wuwei2163 发表于 2017-1-9 11:40
王老师，您好：
      之前一直通过站点数据作图并mask本省数据，没出现问题；今天上午想mask自定义地图， ...

作为mask的shape文件必须是polygon的，这样才能有闭合区域。

wuwei2163

MeteoInfo 发表于 2017-1-9 12:24
作为mask的shape文件必须是polygon的，这样才能有闭合区域。

王老师，这样的文件可以通过MeteoInfo做出来吗？

MeteoInfo

wuwei2163 发表于 2017-1-9 12:28
王老师，这样的文件可以通过MeteoInfo做出来吗？

可以呀

wuwei2163

还是根据菜单中输入地图保存成shp文件做出来吗，王老师

wuwei2163

MeteoInfo 发表于 2017-1-9 12:29
可以呀

王老师，打扰您了，问题已经解决了，之前自定义地图是很早版本MeteoInfo的c#做出来，刚刚用新版作出的地图shp完全没问题，谢谢王老师~

crhj

辛苦了，学习中

呼小喵喵

王老师，我在画城市的累积6小时降水，想试试脚本的运行，多学习学习。运行后提示错误，是等值线颜色那提示错误，不太明白。

1. #Set data folders
   
2. basedir = 'F:/Meteinfo/MeteoInfo_Java_1.4.2_Files/MeteoInfo'
   
3. datadir = os.path.join(basedir, 'sample/MICAPS/surface20150515/surface/plot')
   
4. mapdir = os.path.join(basedir, 'map')
   
5. #Read shape files
   
6. bou2_layer = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
   
7. bou1_layer = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
   
8. hubei_layer = shaperead(os.path.join(mapdir, 'hubei-1.shp'))
   
9. #Read station data
   
10. f = addfile_micaps(os.path.join(datadir, '15051502.000'))
    
11. pr = f.stationdata('Precipitation6h')
    
12. #griddata function - interpolate
    
13. x = arange(108, 116, 0.5)
    
14. y = arange(29, 34, 0.5)
    
15. prg = pr.griddata((x, y), method='idw', radius=3)
    
16. #Plot
    
17. axesm()
    
18. geoshow(bou2_layer, edgecolor='lightgray')
    
19. geoshow(bou1_layer, facecolor=(0,0,255))
    
20. geoshow(hubei_layer, facecolor='r', size=4, labelfield='NAME2', fontname=u'Times New Roman', fontsize=16, yoffset=15)
    
21. geoshow(hubei_layer, visible=False)
    
22. levs = [0.1, 1, 2, 5, 10, 20, 25, 50, 100]
    
23. cols = [(255,255,255),(170,240,255),(120,230,240),(200,220,50),(240,220,20),(255,120,10),(255,90,10), \
    
24.     (240,40,0),(180,10,0),(120,10,0)]
    
25. layer = contourfm(prg, levs, colors=cols)
    
26. masklayer(hubei_layer, [layer])
    
27. colorbar(layer)
    
28. xlim(107, 117)
    
29. ylim(28, 35)
    
30. text(95, 52, u'Accumulated Precipitatio 6h', fontname=u'黑体', fontsize=16)
    
31. text(95, 50, u'(2015-05-14 09:00 至 2015-05-15 02:00)', fontname=u'黑体', fontsize=14)
    

复制代码

提示错误如下：
>>> run script...
  File "<iostream>", line 24
    &#194;&#160; &#194;&#160; (240,40,0),(180,10,0),(120,10,0)]
    ^
SyntaxError: no viable alternative at character '&#160; 15051502.000 (774.45 KB, 下载次数: 0)

2017-5-19 10:28 上传

点击文件名下载附件

MeteoInfo

呼小喵喵 发表于 2017-5-19 10:28
王老师，我在画城市的累积6小时降水，想试试脚本的运行，多学习学习。运行后提示错误，是等值线颜色那提示 ...

你把23，24行合并写在一行试试

MeteoInfo

气象人xyy 发表于 2017-11-17 13:56
请问示例的数据在哪里呢？想试着学习学习

软件Sample目录中有