MeteoInfoLab脚本示例：闪电位置图

呼小喵喵

谢谢王老师，威武~~{:eb502:}{:eb502:}{:eb502:}{:eb502:}

candylu

王老师您好，根据您提供的脚本画闪电的位置分布图，在14行出现了NameError: name 'array' is not defined的报错，请问怎么解决？谢谢老师

MeteoInfo

candylu 发表于 2017-6-26 10:03
王老师您好，根据您提供的脚本画闪电的位置分布图，在14行出现了NameError: name 'array' is not defined的 ...

要在MeteoInfoLab中运行此脚本，参考这里：http://bbs.06climate.com/forum.p ... &extra=page%3D1

candylu

MeteoInfo 发表于 2017-6-26 10:06
要在MeteoInfoLab中运行此脚本，参考这里：http://bbs.06climate.com/forum.php?mod=viewthread&tid=3615 ...

谢谢老师指导，问题已解决。

qdsonglinyy

王老师您好，我拷贝您的脚本运行我的数据，总是提示这样的错误，完全用您的数据尝试，也是一样？
fn = 'F:/grads/lightning6.17/1600test.txt'
tf = open(fn)
lats = []
lons = []
vs = []
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)
lon = array(lons)
lat = array(lats)  
v = array(vs)
axesm()
mlayer = shaperead('F:/grads/lightning6.17/test.shp')
geoshow(mlayer)
ss = makesymbolspec('point', {'value':(-10000,0), 'color':'b', 'marker':'m', 'size':6}, \
    {'value':(0,10000), 'color':'r', 'marker':'+', 'size':6})
layer = scatterm(lon, lat, v, symbolspec=ss)
xlim(119, 123)
ylim(35, 38)
title(u'闪电位置图', fontname=u'黑体', fontsize=18)

qdsonglinyy

王老师，另外想请教，想把雷达的组合反射率作为底图图层，在其上画闪电位置，有类似的脚本么？要怎么样设置？

qdsonglinyy

qdsonglinyy 发表于 2017-7-20 15:37
王老师您好，我拷贝您的脚本运行我的数据，总是提示这样的错误，完全用您的数据尝试，也是一样？
fn = 'F: ...

我的txt格式数据文件

MeteoInfo

qdsonglinyy 发表于 2017-7-20 15:42
我的txt格式数据文件

脚本里读文本文件用的是Python自有的函数，具体用法上网搜搜就知道了。

有事秦天

王老师您好，我看了上面的计算1°*1°的闪电次数的脚本。如果是要计算每平方千米上闪电次数，请问王老师脚本该怎么改?

MeteoInfo

有事秦天 发表于 2017-9-27 07:22
王老师您好，我看了上面的计算1°*1°的闪电次数的脚本。如果是要计算每平方千米上闪电次数，请问王老师脚 ...

可以，需要进行投影转换。