绘制闪电密度图

MeteoInfo · 发表于 2013-5-7 20:03:57

koikaze 发表于 2013-5-7 17:02

登录/注册后可看大图

了解，就是经纬度坐标数组。
还有一点疑问需要请教：
代码中：test.grd

test.grd是已经计算了闪电次数的格点数据。

把次数也读成一个变量，在给格点数据加闪电次数时用这个变量就可以了。

koikaze · 发表于 2013-5-8 09:16:21

MeteoInfo 发表于 2013-5-7 20:03

登录/注册后可看大图

test.grd是已经计算了闪电次数的格点数据。

把次数也读成一个变量，在给格点数据加闪电次数时用这个变 ...

多谢，数据统计还需要另外第二方软件fortran编写代码，感激不尽。

koikaze · 发表于 2013-5-8 15:16:04

本帖最后由 koikaze 于 2013-5-8 15:49 编辑

MeteoInfo 发表于 2013-5-7 20:03

登录/注册后可看大图

test.grd是已经计算了闪电次数的格点数据。

把次数也读成一个变量，在给格点数据加闪电次数时用这个变 ...

王老师，您好。
脚本中：
gData.Data[yIdx, xIdx] += 1
为何要把yIdx放在数组第一位，yIdx应该是纬度数据来的，是否应该为： gData.Data[xIdx, yIdx] += 1

脚本中：
aLS.ImportFromXMLFile(dataDir + "test.lgs")
test.lgs文件是如何生成的呢？

万望指教，谢谢。

MeteoInfo · 发表于 2013-5-8 15:50:34

koikaze 发表于 2013-5-8 15:16

登录/注册后可看大图

王老师，您好。
脚本中：
gData.Data[yIdx, xIdx] += 1

没错，数组元素的顺序是由西向东、由南向北。你仔细体会一下就明白了。

koikaze · 发表于 2013-5-8 16:04:49

MeteoInfo 发表于 2013-5-8 15:50

登录/注册后可看大图

没错，数组元素的顺序是由西向东、由南向北。你仔细体会一下就明白了。

刚才，调试脚本的时候发现如下提示：
Define grid data...
Read lighting data file...
Traceback (most recent call last):
File "<string>", line 107, in <module>
IOError: [Errno 2] 未能找到文件“D:\Temp\ascii\test.lgs”。
是否test.lgs出错了？

MeteoInfo · 发表于 2013-5-8 16:57:45

koikaze 发表于 2013-5-8 16:04

登录/注册后可看大图

刚才，调试脚本的时候发现如下提示：
Define grid data...
Read lighting data file...

脚本用的图例文件：

test.lgs (4.46 KB, 下载次数: 9)

chengxf · 发表于 2013-5-8 21:03:17

王老师，我照那些闪电资料用Surfer做出来的图显得很粗糙，和重庆的底图边缘吻合很差。因为是直接利用的格点的数据作为GRD文件(ASCII格式)，不是经过数据文件插值出来的网格文件，因此也就没有办法调整插值格点数。它总共只有50x41=2050个数据格点，粗糙也是难免的。有什么办法解决呢？

TempLig.jpg (63.08 KB, 下载次数: 2)

MeteoInfo · 发表于 2013-5-8 22:06:08

chengxf 发表于 2013-5-8 21:03

登录/注册后可看大图

王老师，我照那些闪电资料用Surfer做出来的图显得很粗糙，和重庆的底图边缘吻合很差。因为是直接利用的格点 ...

其实次数这种变量并不适合做等值线图，因为是和网格面积相关的，严格的讲应该计算出单位面积发生的闪电次数然后再来做等值线图。目前的做法网格的划分不同的话结果会有很大的区别，比如网格加密的话每个网格的闪电次数必然会减少。

我在此贴（求助！如何做珠江三角洲的底图？
http://bbs.06climate.com/forum.p ... 995&fromuid=106
）的35楼讲到了GrADS和MeteoInfo的maskout机制，Surfer我并不了解，不过从你给出的结果来看应该也是先maskout格点数据然后在绘制等值线，不过在边界追踪上并未做矩形处理。

chengxf · 发表于 2013-5-9 13:45:10

MeteoInfo 发表于 2013-5-8 22:06

登录/注册后可看大图

其实次数这种变量并不适合做等值线图，因为是和网格面积相关的，严格的讲应该计算出单位面积发生的闪电次 ...

谢谢你的回复。我也体会到“次数这种变量并不适合做等值线图”，用Surfer做出的“闪电密度分布”等值线图实在很勉强。倒是MeteoInfo做出来的分布图好得多。那边35楼的回帖也在第一时间拜读过，maskout的先后影响很大。如果先做等值线着色，然后再白化maskout，那样吻合度就会很好，如果是先白化以后再进行画等值线和着色，估计Surfer就是这样干的，因此边缘的吻合度就显得很粗糙。再次谢谢指教，对白化过程有了更深入的了解。

变换 · 发表于 2013-5-11 19:06:23

谢谢楼主！

		自动登录	找回密码
密码			立即注册