“闪电数据的统计处理”程序的使用一例

chengxf

网友@vanni010提问如何使用我曾经发表的闪电统计程序，现根据他提供的数据实际操作如下：
先用记事本打开数据文件看看，原来这些数据是国家气象局ICS 07.060行业标准，根据粗略估计，总共有83567行记录。下图：

再来看看成都市的底图，确定一下它的范围，就是最大最小的经纬度范围值是多少。下图：

运行闪电统计程序，按图中的1、2、3步说明，设置好数据类型和范围参数。下图：

然后调入原始文件：

可以一次调入多个文件，只要按照Windows的多选操作办法，选中多个文件即可。这个例子中只有一个文件。
由于原始数据中包含了闪电强度数据，因此将有提示窗口如下：

调入数据以后会发现，实际的记录数与原来我们估计的记录数不同，少了一些，下面再次打开原始数据，可以发现序号不是连续的，其中有跳越的号码。如下图：

点击导入数据按钮后，程序排除了相当一部分范围之外的数据记录。这是因为我们事先要求排除范围以外的记录。这样做的好处是可以减少要处理的数据量，避免内存溢出故障。当然，你如果知道这些数据都在范围之内，就没有必要打勾，经纬度范围值也没有必要填写。

现在调整经度方向的格点数，程序将自动计算纬度方向的格点数，以及格点之间的间距，单位有度和米两种。调整到你希望的数值，例如1公里见方的格子。按照下图的步骤做：

现在保存一些文件，它们都是今后有用的，参看下图，选择一些你需要的文件来保存：

现在开始作图，按照下图的步骤来做，没什么困难的，但是，做出来的图并不理想：

还是回到 Surfer 中，首先去除网格中的负值(有相关帖子可以参考：Surfer作图如何消除等值线中的负值)：

只是做了一点点工作，图形还是不很满意。余下的工作，希望@vanni010自己去完善。

补充一张，真正的散点图。如果认真数一数，应该有13150个点，这是在矩形框范围之内的。真正落在成都市的闪电还要少，可惜这个程序还没有办法将成都边界图作为范围，排除掉多边形以外的散点。我正在想如何去掉多边形外的各个点，像白化等值线图那样。2014.11.13补充：这个功能已经做到，请看这里： 散点图：剔除多边形外的数据点


闪电统计程序也作了小小的修改，这是最新的(没有病毒)：
闪电统计程序(20140928).rar (1.67 MB, 下载次数: 235)

2014-10-1 07:50 上传

点击文件名下载附件
阅读权限: 10

为适应安徽省DL-II闪电定位仪的不同数据，又修改了程序，新的程序如下：
LightningStatistics20160516.rar (1.05 MB, 下载次数: 68)

2016-5-19 09:47 上传

点击文件名下载附件
阅读权限: 10

kongfeng0824

谢谢程老师的力作

kongfeng0824

抓紧研读一下

mofangbao

程老师辛苦了！
我们不得已删除了一些用户上传的观测数据，希望程老师能理解！

林夕

特别喜欢程老师每一个教程，每一个细节都很清楚，详细的解释说明对于我们这些初学者帮助特别大，谢谢程老师

aa22556688

程老师好厉害，非常感谢能够这样帮助新手。

vanni010

谢谢程老师{:eb303:}{:eb303:}{:eb303:}{:eb303:}

ltf991

学习了，真是好人呀

vanni010

程老师您好，如果依靠闪电统计程序加surfer做一“24h降水量分布等值线图”叠加“当天正负地闪分布的散点图”，怎样来画呢？我设想第一步先依靠降水资料绘制等值线图，第二步利用闪电统计程序处理过经纬度的原始资料的经纬度绘制散点图，这里根据强度区分正负地闪后，正负地闪怎样表现在surfer作图中呢？

chengxf

vanni010 发表于 2014-9-26 09:41
程老师您好，如果依靠闪电统计程序加surfer做一“24h降水量分布等值线图”叠加“当天正负地闪分布的散点图 ...

要区分正负地闪，可以采用分类散点图。除了经纬度数据外，还有一个闪电强度数据，它有正负值。你将第三列强度这一列作为分类的依据，在分类规则的对话框中，简单地分成>0 和 <=0 两类，你就可以对这两类进行属性设置，包括它的符号样式、大小、颜色等等。没错，你可以将它和雨量等值线图叠加。
现在的问题是如何将成都区域外的散点都去掉，只留多边形区域之内的散点。这要有一段程序来判断这些点是否落在多边形之内。