| 查看: 677 | 回复: 3 | ||
[求助]
求助编程序，算雷击大地密度，用网格法。
|
|
求热心人帮忙，我需要计算地闪密度。把一片区域分割成经纬度为0.1*0.1的网格，然后判断每次闪电是否发生在这个网格里，如果是，就把发生在这个网格里的闪电次数累加起来，再除以面积，就得到每一个网格的地闪密度。我需要用程序实现它，求大神们帮帮我，不甚感激啊！毕业论文用的。 感觉我好像没讲清楚，我的算法和下面我从一篇文献里的算法基本一样，如下： 格子的中心位置为评估点所在经纬度。每个格子边长为３ｋｍ，每边共有９ 个格子。算法具体步骤如下： 评估点附近网格分布图：网格是9*9的网格 （１）在程序中声明４个二维数组：①ＡｒｒａｙＮｇ，②ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ，③ＡｒｒａｙＢｏｘｓＬｎｇ，④Ａｒｒａｙ－ＢｏｘｓＬａｔ。分别用来存放图１中每个格子的①Ｎｇ值，②雷电流平均强度，③格子中心点经度，④格子中心点纬度。 （２）已知用户输入的评估点经纬度，而且根据地球上任意两点的经纬度又可以得到两点间距离。计算公式如下： ｄ（ｘ１，ｙ１，ｘ２，ｙ２）＝ｒａｒｃｃｏｓ［ｓｉｎ（ｘ１）ｓｉｎ（ｘ２）＋ｃｏｓ（ｘ１）ｃｏｓ（ｘ２）ｃｏｓ（ｙ１－ｙ２）］ 其中：ｘ１，ｙ１是纬度、经度的弧度单位，ｒ 为地球半径。 据此，可以经过计算得出每个格子的经纬度范围，ＡｒｒａｙＢｏｘｓＬｎｇ，ＡｒｒａｙＢｏｘｓＬａｔ能够通过循环成功赋值。 （３）逐行读取经过处理的闪电定位系统（ＬＬＳ）数据，取闪电的经纬度存入局部变量ｌｎｇＴｍｐ、ｌａｔＴｍｐ。 使用ｉｆ语句判断该行闪电定位数据是否属于图１中所示９×９＝８１个格子。如果属于，则ＡｒｒａｙＮｇ数组中第ｉ行ｊ列格子数据执行自加运算（ＡｒｒａｙＮｇ［ｉ，ｊ］＋＋；）；同时，地闪强度也累加（ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ［ｉ，ｊ］＋＝ ＬｉｇｈｔｎｉｎｇＣｕｒｒｅｎｔ；）。 （４）对数据源循环一遍后，使用ｆｏｒ循环语句重新对８１个格子进行循环赋值。使用公式为：平均电流强度＝ 总电流强度／闪击次数 Ｎｇ ＝ 总地闪次数／资料年数／格子面积 经过以上４步，即可得出每个格子的Ｎｇ值和雷 |
回复此楼» 猜你喜欢
生物学308求调剂（一志愿华东师大）
已经有6人回复
-
调剂求收留
已经有7人回复
-
272分材料子求调剂
已经有36人回复
-
275求调剂
已经有8人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有20人回复
-
070300化学学硕311分求调剂
已经有19人回复
-
材料与化工调剂
已经有13人回复
-
材料与化工调剂
已经有33人回复
-
复试调剂
已经有7人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有17人回复
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
3楼2013-05-16 13:29:43
deephill
铁杆木虫 (职业作家)
- 应助: 59 (初中生)
- 金币: 5262.3
- 散金: 100
- 红花: 16
- 帖子: 4425
- 在线: 425小时
- 虫号: 305318
- 注册: 2006-12-09
- 性别: GG
- 专业: 微量元素地球化学
4楼2013-05-16 17:52:00
