求助帮忙编程！急用！不甚感激！

dl凋零

求热心人帮忙，我需要计算地闪密度。把一片区域分割成经纬度为0.1*0.1的网格，然后判断每次闪电是否发生在这个网格里，如果是，就把发生在这个网格里的闪电次数累加起来，再除以面积，就得到每一个网格的地闪密度。我需要用程序实现它，求大神们帮帮我，不甚感激啊！毕业论文用的。感觉我好像没讲清楚，我的算法和下面我从一篇文献里的算法基本一样，如下：
  格子的中心位置为评估点所在经纬度。每个格子边长为３ｋｍ，每边共有９ 个格子。算法具体步骤如下：
评估点附近网格分布图：网格是9*9的网格
（１）在程序中声明４个二维数组：①ＡｒｒａｙＮｇ，②ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ，③ＡｒｒａｙＢｏｘｓＬｎｇ，④Ａｒｒａｙ－ＢｏｘｓＬａｔ。分别用来存放图１中每个格子的①Ｎｇ值，②雷电流平均强度，③格子中心点经度，④格子中心点纬度。
（２）已知用户输入的评估点经纬度，而且根据地球上任意两点的经纬度又可以得到两点间距离。计算公式如下：
ｄ（ｘ１，ｙ１，ｘ２，ｙ２）＝ｒａｒｃｃｏｓ［ｓｉｎ（ｘ１）ｓｉｎ（ｘ２）＋ｃｏｓ（ｘ１）ｃｏｓ（ｘ２）ｃｏｓ（ｙ１－ｙ２）］
其中：ｘ１，ｙ１是纬度、经度的弧度单位，ｒ 为地球半径。
据此，可以经过计算得出每个格子的经纬度范围，ＡｒｒａｙＢｏｘｓＬｎｇ，ＡｒｒａｙＢｏｘｓＬａｔ能够通过循环成功赋值。
（３）逐行读取经过处理的闪电定位系统（ＬＬＳ）数据，取闪电的经纬度存入局部变量ｌｎｇＴｍｐ、ｌａｔＴｍｐ。
使用ｉｆ语句判断该行闪电定位数据是否属于图１中所示９×９＝８１个格子。如果属于，则ＡｒｒａｙＮｇ数组中第ｉ行ｊ列格子数据执行自加运算（ＡｒｒａｙＮｇ［ｉ，ｊ］＋＋；）；同时，地闪强度也累加（ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ［ｉ，ｊ］＋＝ ＬｉｇｈｔｎｉｎｇＣｕｒｒｅｎｔ；）。
（４）对数据源循环一遍后，使用ｆｏｒ循环语句重新对８１个格子进行循环赋值。使用公式为：平均电流强度＝ 总电流强度／闪击次数
                     Ｎｇ ＝ 总地闪次数／资料年数／格子面积
经过以上４步，即可得出每个格子的Ｎｇ值和雷

1912347158

同求。很想知道。

coco4046

来学习1一下！

chengxf

coco4046 发表于 2014-9-30 10:46
来学习1一下！

楼主发了这个帖子后应该有一年时间没来论坛了。其实他的问题我比他早几天发布的“闪电统计程序”就可以解决，而且非常直观。有现成的程序，又是免费的，为何不用？