MeteoInfoLab脚本示例：闪电位置图

任我行

你下载MeteoInfo最新版本试试

呼小喵喵

谢谢王老师

呼小喵喵

老师，我现在数据格式跟您用的格式是一致的，但不知道为什么读不了闪电数据，显示错误。

MeteoInfo

呼小喵喵 发表于 2017-5-15 17:38
老师，我现在数据格式跟您用的格式是一致的，但不知道为什么读不了闪电数据，显示错误。

是不是冒号弄成中文的了？

呼小喵喵

MeteoInfo 发表于 2017-5-15 22:49
是不是冒号弄成中文的了？

检查了是英文的，但运行就提示这个错误，不知道为什么。数据格式我也检查了一遍，还是不对。

MeteoInfo

呼小喵喵 发表于 2017-5-16 09:19
检查了是英文的，但运行就提示这个错误，不知道为什么。数据格式我也检查了一遍，还是不对。

你把你的脚本程序作为附件贴出来，看看问题在哪里

呼小喵喵

MeteoInfo 发表于 2017-5-16 09:28
你把你的脚本程序作为附件贴出来，看看问题在哪里

王老师，这是代码。

1. fn = 'F:/Meteinfo/20150514.txt'
   
2. tf = open(fn)
   
3. lats = []
   
4. lons = []
   
5. vs = []
   
6. for aline in tf:  
   
7.     datalist = aline.split()
   
8.     lat = float(datalist[3].split('=')[1])
   
9.     lon = float(datalist[4].split('=')[1])
   
10.     v = float(datalist[5].split('=')[1])
    
11.     lats.append(lat)
    
12.     lons.append(lon)
    
13.     vs.append(v)
    
14. lon = array(lons)
    
15. lat = array(lats)  
    
16. v = array(vs)
    
17. axesm()
    
18. mlayer = shaperead('F:/Meteinfo/MeteoInfo/CHN_adm/hubei-1.shp')
    
19. geoshow(mlayer)
    
20. ss = makesymbolspec('point', {'value':(-10000,0), 'color':'b', 'marker':'m', 'size':6}, \
    
21.     {'value':(0,10000), 'color':'r', 'marker':'+', 'size':6})
    
22. layer = scatterm(lon, lat, v, symbolspec=ss)
    
23. xlim(108, 117)
    
24. ylim(29, 34)
    
25. title(u'闪电位置图', fontname=u'黑体', fontsize=18)

复制代码

MeteoInfo

呼小喵喵 发表于 2017-5-16 09:47
王老师，这是代码。

代码应该没问题，你把20150514.txt文件附上，看看文件的内容该怎么读。

呼小喵喵

20150514.txt (1.76 MB, 下载次数: 4)

2017-5-16 11:27 上传

点击文件名下载附件
王老师，这是数据附件。

MeteoInfo

呼小喵喵 发表于 2017-5-16 11:27
王老师，这是数据附件。

数据中有空白行，加一条判断语句即可：
for aline in tf:   
    if aline.isspace():
        continue
    datalist = aline.split()
    lat = float(datalist[3].split('=')[1])
    lon = float(datalist[4].split('=')[1])
    v = float(datalist[5].split('=')[1])
    lats.append(lat)
    lons.append(lon)
    vs.append(v)