 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 671 | 回復: 3 | |||
nirvana_fwl新蟲 (初入文壇)
|
[求助]
求助編程序，算雷擊大地密度，用網(wǎng)格法。
|
|
求熱心人幫忙，我需要計算地閃密度。把一片區(qū)域分割成經(jīng)緯度為0.1*0.1的網(wǎng)格，然后判斷每次閃電是否發(fā)生在這個網(wǎng)格里，如果是，就把發(fā)生在這個網(wǎng)格里的閃電次數(shù)累加起來，再除以面積，就得到每一個網(wǎng)格的地閃密度。我需要用程序實現(xiàn)它，求大神們幫幫我，不甚感激��！畢業(yè)論文用的。 感覺我好像沒講清楚，我的算法和下面我從一篇文獻里的算法基本一樣，如下： 格子的中心位置為評估點所在經(jīng)緯度。每個格子邊長為３ｋｍ，每邊共有９ 個格子。算法具體步驟如下： 評估點附近網(wǎng)格分布圖：網(wǎng)格是9*9的網(wǎng)格 （１）在程序中聲明４個二維數(shù)組：①ＡｒｒａｙＮｇ，②ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ，③ＡｒｒａｙＢｏｘｓＬｎｇ，④Ａｒｒａｙ－ＢｏｘｓＬａｔ。分別用來存放圖１中每個格子的①Ｎｇ值，②雷電流平均強度，③格子中心點經(jīng)度，④格子中心點緯度。 （２）已知用戶輸入的評估點經(jīng)緯度，而且根據(jù)地球上任意兩點的經(jīng)緯度又可以得到兩點間距離。計算公式如下： ｄ（ｘ１，ｙ１，ｘ２，ｙ２）＝ｒａｒｃｃｏｓ［ｓｉｎ（ｘ１）ｓｉｎ（ｘ２）＋ｃｏｓ（ｘ１）ｃｏｓ（ｘ２）ｃｏｓ（ｙ１－ｙ２）］ 其中：ｘ１，ｙ１是緯度、經(jīng)度的弧度單位，ｒ 為地球半徑。 據(jù)此，可以經(jīng)過計算得出每個格子的經(jīng)緯度范圍，ＡｒｒａｙＢｏｘｓＬｎｇ，ＡｒｒａｙＢｏｘｓＬａｔ能夠通過循環(huán)成功賦值。 （３）逐行讀取經(jīng)過處理的閃電定位系統(tǒng)（ＬＬＳ）數(shù)據(jù)，取閃電的經(jīng)緯度存入局部變量ｌｎｇＴｍｐ、ｌａｔＴｍｐ。 使用ｉｆ語句判斷該行閃電定位數(shù)據(jù)是否屬于圖１中所示９×９＝８１個格子。如果屬于，則ＡｒｒａｙＮｇ數(shù)組中第ｉ行ｊ列格子數(shù)據(jù)執(zhí)行自加運算（ＡｒｒａｙＮｇ［ｉ，ｊ］＋＋；）；同時，地閃強度也累加（ＡｒｒａｙＡｖｇＬｉｇＣｕｒｔ［ｉ，ｊ］＋＝ ＬｉｇｈｔｎｉｎｇＣｕｒｒｅｎｔ；）。 （４）對數(shù)據(jù)源循環(huán)一遍后，使用ｆｏｒ循環(huán)語句重新對８１個格子進行循環(huán)賦值。使用公式為：平均電流強度＝ 總電流強度／閃擊次數(shù) Ｎｇ ＝ 總地閃次數(shù)／資料年數(shù)／格子面積 經(jīng)過以上４步，即可得出每個格子的Ｎｇ值和雷 |
回復此樓» 猜你喜歡
數(shù)一英一271專碩（085401）求調(diào)劑，可跨
已經(jīng)有8人回復
-
295求調(diào)劑
已經(jīng)有4人回復
-
0856求調(diào)劑
已經(jīng)有5人回復
-
材料學碩333求調(diào)劑
已經(jīng)有10人回復
-
299求調(diào)劑
已經(jīng)有6人回復
-
一志愿鄭州大學，080500學碩，總分317分求調(diào)劑
已經(jīng)有6人回復
-
289求調(diào)劑
已經(jīng)有12人回復
-
070305高分子化學與物理 304分求調(diào)劑
已經(jīng)有6人回復
-
0703化學調(diào)劑，求導師收
已經(jīng)有10人回復
-
各位老師好，我的一志愿為北京科技大學085601材料專碩
已經(jīng)有5人回復
libralibra
至尊木蟲 (著名寫手)
驃騎將軍
- 程序強帖: 40
- 應助: 817 (博后)
- 金幣: 12914.1
- 紅花: 64
- 帖子: 2238
- 在線: 287.3小時
- 蟲號: 696514
- 注冊: 2009-02-05
- 專業(yè): 計算機軟件
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +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
感謝參與，應助指數(shù) +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,保存經(jīng)緯度信息,2個屬性,加一個方法計算距離 再弄一個類叫Grid,保存Ng,雷電強度,與一個Point變量(中心經(jīng)緯度),加一個方法地閃加1 不過看你的描述,感覺格子光保存中心點經(jīng)緯度不夠方便,因為有寬度,雖然每次可以計算,但是由于格子一旦初始化完成就固定了,最好還是把每個格子的經(jīng)緯度范圍保存起來作為Grid的屬性,用[x,y,w,h]這種頂點寬高法也行,用[x1,y1,x2,y2]這種兩點法也行 這樣程序就只需要用一個Gfid的二維數(shù)組來保存,9*9 |
2樓2013-05-15 21:59:54
nirvana_fwl
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 10.3
- 帖子: 20
- 在線: 3.4小時
- 蟲號: 2163332
- 注冊: 2012-12-03
- 專業(yè): 大氣遙感和大氣探測
3樓2013-05-16 13:29:43
deephill
鐵桿木蟲 (職業(yè)作家)
- 應助: 59 (初中生)
- 金幣: 5262.3
- 散金: 100
- 紅花: 16
- 帖子: 4425
- 在線: 425小時
- 蟲號: 305318
- 注冊: 2006-12-09
- 性別: GG
- 專業(yè): 微量元素地球化學
4樓2013-05-16 17:52:00
相關版塊跳轉
我要訂閱樓主 nirvana_fwl 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+6
|
nnnnnnn5 2026-03-25 | 6/300 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學 070300 學碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 生物技術與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|