| 4 | 1/1 | 返回列表 |
| 查看: 667 | 回復(fù): 3 | ||
nirvana_fwl新蟲 (初入文壇)
|
[求助]
求助編程序,算雷擊大地密度,用網(wǎng)格法。
|
|
求熱心人幫忙,我需要計算地閃密度。把一片區(qū)域分割成經(jīng)緯度為0.1*0.1的網(wǎng)格,然后判斷每次閃電是否發(fā)生在這個網(wǎng)格里,如果是,就把發(fā)生在這個網(wǎng)格里的閃電次數(shù)累加起來,再除以面積,就得到每一個網(wǎng)格的地閃密度。我需要用程序?qū)崿F(xiàn)它,求大神們幫幫我,不甚感激!畢業(yè)論文用的。 感覺我好像沒講清楚,我的算法和下面我從一篇文獻里的算法基本一樣,如下: 格子的中心位置為評估點所在經(jīng)緯度。每個格子邊長為3km,每邊共有9 個格子。算法具體步驟如下: 評估點附近網(wǎng)格分布圖:網(wǎng)格是9*9的網(wǎng)格 (1)在程序中聲明4個二維數(shù)組:①ArrayNg,②ArrayAvgLigCurt,③ArrayBoxsLng,④Array-BoxsLat。分別用來存放圖1中每個格子的①Ng值,②雷電流平均強度,③格子中心點經(jīng)度,④格子中心點緯度。 (2)已知用戶輸入的評估點經(jīng)緯度,而且根據(jù)地球上任意兩點的經(jīng)緯度又可以得到兩點間距離。計算公式如下: d(x1,y1,x2,y2)=rarccos[sin(x1)sin(x2)+cos(x1)cos(x2)cos(y1-y2)] 其中:x1,y1是緯度、經(jīng)度的弧度單位,r 為地球半徑。 據(jù)此,可以經(jīng)過計算得出每個格子的經(jīng)緯度范圍,ArrayBoxsLng,ArrayBoxsLat能夠通過循環(huán)成功賦值。 (3)逐行讀取經(jīng)過處理的閃電定位系統(tǒng)(LLS)數(shù)據(jù),取閃電的經(jīng)緯度存入局部變量lngTmp、latTmp。 使用if語句判斷該行閃電定位數(shù)據(jù)是否屬于圖1中所示9×9=81個格子。如果屬于,則ArrayNg數(shù)組中第i行j列格子數(shù)據(jù)執(zhí)行自加運算(ArrayNg[i,j]++;);同時,地閃強度也累加(ArrayAvgLigCurt[i,j]+= LightningCurrent;)。 (4)對數(shù)據(jù)源循環(huán)一遍后,使用for循環(huán)語句重新對81個格子進行循環(huán)賦值。使用公式為:平均電流強度= 總電流強度/閃擊次數(shù) Ng = 總地閃次數(shù)/資料年數(shù)/格子面積 經(jīng)過以上4步,即可得出每個格子的Ng值和雷 |
至尊木蟲 (著名寫手)
驃騎將軍
|
你看的論文挺早的吧,看起來像純c程序,這個設(shè)計很不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 |

新蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +4 | 腳滑的守法公民 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|
[考研] 329求調(diào)劑 +4 | 星野? 2026-03-26 | 4/200 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +4 | 舉個栗子oi 2026-03-24 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強芒? 2026-03-21 | 8/400 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|