MeteoInfoLab脚本示例：闪电位置图

MeteoInfo

wuwei2163 发表于 2016-10-19 10:53
嗯，是的，王老师。我上面的脚本就是先对河南区域内的闪电判断，然后再统计变化。但运行程序后不出图或者 ...

参考下面的脚本：

1. #Read lighting data
   
2. fn = 'D:/Temp/ascii/lighting/2009_06_06.txt'
   
3. tf = open(fn)
   
4. m_henan = shaperead('D:/Temp/Map/henan.shp')
   
5. times = []
   
6. nums = []
   
7. for h in range(0, 24):
   
8.     times.append(datetime.datetime(2009,6,6,h))
   
9.     nums.append(0)
   
10. for aline in tf:   
    
11.     datalist = aline.split()
    
12.     lat = float(datalist[3].split('=')[1])
    
13.     lon = float(datalist[4].split('=')[1])
    
14.     isin = inpolygon(lon, lat, m_henan.shapes()[0])
    
15.     if isin:
    
16.         h = int(datalist[2].split(':')[0])
    
17.         nums[h] = nums[h] + 1
    
18. 
    
19. #Plot
    
20. plot(times, nums, '-ro')
    
21. xaxis(axistype='time')
    
22. ylabel('Lighting number')
    
23. xlim(times[0], times[-1])
    
24. title('Lighting number hourly variation in ' + times[0].strftime('%Y/%m/%d'))

#Read lighting data<br /> fn = 'D:/Temp/ascii/lighting/2009_06_06.txt'<br /> tf = open(fn)<br /> m_henan = shaperead('D:/Temp/Map/henan.shp')<br /> times = []<br /> nums = []<br /> for h in range(0, 24):<br />     times.append(datetime.datetime(2009,6,6,h))<br />     nums.append(0)<br /> for aline in tf:    <br />     datalist = aline.split()<br />     lat = float(datalist[3].split('=')[1])<br />     lon = float(datalist[4].split('=')[1])<br />     isin = inpolygon(lon, lat, m_henan.shapes()[0])<br />     if isin:<br />         h = int(datalist[2].split(':')[0])<br />         nums[h] = nums[h] + 1<br /> <br /> #Plot<br /> plot(times, nums, '-ro')<br /> xaxis(axistype='time')<br /> ylabel('Lighting number')<br /> xlim(times[0], times[-1])<br /> title('Lighting number hourly variation in ' + times[0].strftime('%Y/%m/%d'))

wuwei2163

MeteoInfo 发表于 2016-10-19 11:54
参考下面的脚本：

非常感谢王老师。。。根据您的参考脚本，对照着一点一点调试自己的脚本程序，找到自己的错误原因了，原来是把时间的循环放到后面，导致这个for循环是个空的。同时通过这个，自己确实收获不少，准备开始学习教程了~

qdsonglinyy

王老师好厉害 谢谢

lu19320

王老师，你好！由于编程基础比较弱，你编写的脚本我只能勉强看懂，自己动手写脚本还比较难。现在部分省市已经开始上三维闪电定位仪了，闪电的类型多了IC云闪（CG为云地闪），三维闪电数据格式（csv格式）跟二维的有所不同，如附件所示，如果要做三维的闪电统计，数据该如何读取呢，烦请指点。

MeteoInfo

lu19320 发表于 2017-4-7 16:21
王老师，你好！由于编程基础比较弱，你编写的脚本我只能勉强看懂，自己动手写脚本还比较难。现在部分省市已 ...

1. fn = 'D:/Temp/ascii/lighting/1.csv'
   
2. table = readtable(fn, delimiter=',', format='%i%{yyyy/M/d H:mm}D%s%4f%3s%4f', \
   
3.     encoding='gb2312')
   
4. lon = table[u'经度']
   
5. lat = table[u'纬度']
   
6. v = table[u'强度(千安)']
   
7. ltype = table[u'类型']
   
8. #Plot
   
9. axesm()
   
10. mlayer = shaperead('D:/Temp/map/bou2_4p.shp')
    
11. geoshow(mlayer)
    
12. ss = makesymbolspec('point', {'value':(-10000,0), 'color':'b', 'marker':'m', 'size':6, 'caption':'Negative'}, \
    
13.     {'value':(0,10000), 'color':'r', 'marker':'+', 'size':6, 'caption':'Positive'})
    
14. layer = scatterm(lon, lat, v, symbolspec=ss)
    
15. legend(legend=layer.legend(), loc='lower left')
    
16. xlim(100, 120)
    
17. ylim(20, 30)
    
18. title('Lighting locations')

fn = 'D:/Temp/ascii/lighting/1.csv'<br /> table = readtable(fn, delimiter=',', format='%i%{yyyy/M/d H:mm}D%s%4f%3s%4f', \<br />     encoding='gb2312')<br /> lon = table[u'经度']<br /> lat = table[u'纬度']<br /> v = table[u'强度(千安)']<br /> ltype = table[u'类型']<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(mlayer)<br /> ss = makesymbolspec('point', {'value':(-10000,0), 'color':'b', 'marker':'m', 'size':6, 'caption':'Negative'}, \<br />     {'value':(0,10000), 'color':'r', 'marker':'+', 'size':6, 'caption':'Positive'})<br /> layer = scatterm(lon, lat, v, symbolspec=ss)<br /> legend(legend=layer.legend(), loc='lower left')<br /> xlim(100, 120)<br /> ylim(20, 30) <br /> title('Lighting locations')

lu19320

MeteoInfo 发表于 2017-4-7 22:34

非常感谢王老师!不过我把老师脚本修改路径运行后，出现了如下错误：
Traceback (most recent call last):
  File "<iostream>", line 4, in <module>
  File "F:\MeteoInfo_Java_1.4.1_Files\MeteoInfo\pylib\mipylib\numeric\mitable.py", line 29, in __getitem__
    coldata = self.data.getColumnData(key)
        at org.meteoinfo.table.ColumnData.<init>(ColumnData.java:33)
        at org.meteoinfo.table.DataTable.getColumnData(DataTable.java:550)
        at org.meteoinfo.table.DataTable.getColumnData(DataTable.java:529)
        at org.meteoinfo.data.TableData.getColumnData(TableData.java:385)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)

java.lang.NullPointerException: java.lang.NullPointerException

MeteoInfo

lu19320 发表于 2017-4-8 14:25
非常感谢王老师!不过我把老师脚本修改路径运行后，出现了如下错误：
Traceback (most recent call last) ...

确保你的文件路径是对的

lu19320

MeteoInfo 发表于 2017-4-8 15:48
确保你的文件路径是对的

找到原因了，数据里有两列的顺序反了。还有一个问题，三维闪电类型有云闪（只有正闪）和云地闪（正地闪和负地闪），怎样通过控制闪电类型、高度值、强度值，把云闪、正地闪、和负地闪三种类型一起画出来，正地闪里又不包括云闪。

MeteoInfo

lu19320 发表于 2017-4-8 16:28
找到原因了，数据里有两列的顺序反了。还有一个问题，三维闪电类型有云闪（只有正闪）和云地闪（正地闪和 ...

我这里没有问题，要不你下载最新版本试试。

xuzheyong

赞赞赞{:eb502:}