| 查看: 819 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助编程序，算雷击大地密度，用网格法。
|
||
|
求热心人帮忙，我需要计算地闪密度。把一片区域分割成经纬度为0.1*0.1的网格，然后判断每次闪电是否发生在这个网格里，如果是，就把发生在这个网格里的闪电次数累加起来，再除以面积，就得到每一个网格的地闪密度。我需要用程序实现它，求大神们帮帮我，不甚感激啊！毕业论文用的。 感觉我好像没讲清楚，我的算法和下面我从一篇文献里的算法基本一样，如下： 格子的中心位置为评估点所在经纬度。每个格子边长为３ｋｍ，每边共有９ 个格子。算法具体步骤如下： 评估点附近网格分布图：网格是9*9的网格 （１）在程序中声明４个二维数组：①ＡｒｒａｙＮｇ，②ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ，③ＡｒｒａｙＢｏｘｓＬｎｇ，④Ａｒｒａｙ－ＢｏｘｓＬａｔ。分别用来存放图１中每个格子的①Ｎｇ值，②雷电流平均强度，③格子中心点经度，④格子中心点纬度。 （２）已知用户输入的评估点经纬度，而且根据地球上任意两点的经纬度又可以得到两点间距离。计算公式如下： ｄ（ｘ１，ｙ１，ｘ２，ｙ２）＝ｒａｒｃｃｏｓ［ｓｉｎ（ｘ１）ｓｉｎ（ｘ２）＋ｃｏｓ（ｘ１）ｃｏｓ（ｘ２）ｃｏｓ（ｙ１－ｙ２）］ 其中：ｘ１，ｙ１是纬度、经度的弧度单位，ｒ 为地球半径。 据此，可以经过计算得出每个格子的经纬度范围，ＡｒｒａｙＢｏｘｓＬｎｇ，ＡｒｒａｙＢｏｘｓＬａｔ能够通过循环成功赋值。 （３）逐行读取经过处理的闪电定位系统（ＬＬＳ）数据，取闪电的经纬度存入局部变量ｌｎｇＴｍｐ、ｌａｔＴｍｐ。 使用ｉｆ语句判断该行闪电定位数据是否属于图１中所示９×９＝８１个格子。如果属于，则ＡｒｒａｙＮｇ数组中第ｉ行ｊ列格子数据执行自加运算（ＡｒｒａｙＮｇ［ｉ，ｊ］＋＋；）；同时，地闪强度也累加（ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ［ｉ，ｊ］＋＝ ＬｉｇｈｔｎｉｎｇＣｕｒｒｅｎｔ；）。 （４）对数据源循环一遍后，使用ｆｏｒ循环语句重新对８１个格子进行循环赋值。使用公式为：平均电流强度＝ 总电流强度／闪击次数 Ｎｇ ＝ 总地闪次数／资料年数／格子面积 经过以上４步，即可得出每个格子的Ｎｇ值和雷 |
回复此楼» 猜你喜欢
今天基金会出结果吗？20260819
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有118人回复
-
今天放榜吗？
已经有15人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
filecode，4个jtjc了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
3楼2013-05-16 13:29:43
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xzhdty: 金币+1, 感谢参与谢谢骠骑将军 2013-05-15 22:52:50
ben_ladeng: 专家考核 2013-05-16 13:22:05
nirvana_fwl: 金币+6, ★有帮助 2013-05-21 13:24:44
感谢参与，应助指数 +1
xzhdty: 金币+1, 感谢参与谢谢骠骑将军 2013-05-15 22:52:50
ben_ladeng: 专家考核 2013-05-16 13:22:05
nirvana_fwl: 金币+6, ★有帮助 2013-05-21 13:24:44
|
你看的论文挺早的吧,看起来像纯c程序,这个设计很不OOP, c++的做法是弄一个类叫Point,保存经纬度信息,2个属性,加一个方法计算距离 再弄一个类叫Grid,保存Ng,雷电强度,与一个Point变量(中心经纬度),加一个方法地闪加1 不过看你的描述,感觉格子光保存中心点经纬度不够方便,因为有宽度,虽然每次可以计算,但是由于格子一旦初始化完成就固定了,最好还是把每个格子的经纬度范围保存起来作为Grid的属性,用[x,y,w,h]这种顶点宽高法也行,用[x1,y1,x2,y2]这种两点法也行 这样程序就只需要用一个Gfid的二维数组来保存,9*9 |
2楼2013-05-15 21:59:54
