| 5 | 1/1 | 返回列表 |
| 查看: 670 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
nirvana_fwl新蟲 (初入文壇)
|
[求助]
求助編程序,算雷擊大地密度,用網(wǎng)格法。
|
|
|
求熱心人幫忙,我需要計(jì)算地閃密度。把一片區(qū)域分割成經(jīng)緯度為0.1*0.1的網(wǎng)格,然后判斷每次閃電是否發(fā)生在這個(gè)網(wǎng)格里,如果是,就把發(fā)生在這個(gè)網(wǎng)格里的閃電次數(shù)累加起來(lái),再除以面積,就得到每一個(gè)網(wǎng)格的地閃密度。我需要用程序?qū)崿F(xiàn)它,求大神們幫幫我,不甚感激。‘厴I(yè)論文用的。 感覺(jué)我好像沒(méi)講清楚,我的算法和下面我從一篇文獻(xiàn)里的算法基本一樣,如下: 格子的中心位置為評(píng)估點(diǎn)所在經(jīng)緯度。每個(gè)格子邊長(zhǎng)為3km,每邊共有9 個(gè)格子。算法具體步驟如下: 評(píng)估點(diǎn)附近網(wǎng)格分布圖:網(wǎng)格是9*9的網(wǎng)格 (1)在程序中聲明4個(gè)二維數(shù)組:①ArrayNg,②ArrayAvgLigCurt,③ArrayBoxsLng,④Array-BoxsLat。分別用來(lái)存放圖1中每個(gè)格子的①Ng值,②雷電流平均強(qiáng)度,③格子中心點(diǎn)經(jīng)度,④格子中心點(diǎn)緯度。 (2)已知用戶輸入的評(píng)估點(diǎn)經(jīng)緯度,而且根據(jù)地球上任意兩點(diǎn)的經(jīng)緯度又可以得到兩點(diǎn)間距離。計(jì)算公式如下: 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)過(guò)計(jì)算得出每個(gè)格子的經(jīng)緯度范圍,ArrayBoxsLng,ArrayBoxsLat能夠通過(guò)循環(huán)成功賦值。 (3)逐行讀取經(jīng)過(guò)處理的閃電定位系統(tǒng)(LLS)數(shù)據(jù),取閃電的經(jīng)緯度存入局部變量lngTmp、latTmp。 使用if語(yǔ)句判斷該行閃電定位數(shù)據(jù)是否屬于圖1中所示9×9=81個(gè)格子。如果屬于,則ArrayNg數(shù)組中第i行j列格子數(shù)據(jù)執(zhí)行自加運(yùn)算(ArrayNg[i,j]++;);同時(shí),地閃強(qiáng)度也累加(ArrayAvgLigCurt[i,j]+= LightningCurrent;)。 (4)對(duì)數(shù)據(jù)源循環(huán)一遍后,使用for循環(huán)語(yǔ)句重新對(duì)81個(gè)格子進(jìn)行循環(huán)賦值。使用公式為:平均電流強(qiáng)度= 總電流強(qiáng)度/閃擊次數(shù) Ng = 總地閃次數(shù)/資料年數(shù)/格子面積 經(jīng)過(guò)以上4步,即可得出每個(gè)格子的Ng值和雷 |
新蟲 (初入文壇)
至尊木蟲 (著名寫手)
驃騎將軍
|
你看的論文挺早的吧,看起來(lái)像純c程序,這個(gè)設(shè)計(jì)很不OOP, c++的做法是弄一個(gè)類叫Point,保存經(jīng)緯度信息,2個(gè)屬性,加一個(gè)方法計(jì)算距離 再弄一個(gè)類叫Grid,保存Ng,雷電強(qiáng)度,與一個(gè)Point變量(中心經(jīng)緯度),加一個(gè)方法地閃加1 不過(guò)看你的描述,感覺(jué)格子光保存中心點(diǎn)經(jīng)緯度不夠方便,因?yàn)橛袑挾?雖然每次可以計(jì)算,但是由于格子一旦初始化完成就固定了,最好還是把每個(gè)格子的經(jīng)緯度范圍保存起來(lái)作為Grid的屬性,用[x,y,w,h]這種頂點(diǎn)寬高法也行,用[x1,y1,x2,y2]這種兩點(diǎn)法也行 這樣程序就只需要用一個(gè)Gfid的二維數(shù)組來(lái)保存,9*9 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +12 | minht 2026-03-26 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 085602 307分 求調(diào)劑 +6 | 不知道叫什么! 2026-03-26 | 6/300 |
|
|
[考研] 315分求調(diào)劑 +6 | 26考研上岸版26 2026-03-26 | 6/300 |
|
|
[考研] 286求調(diào)劑 +8 | PolarBear11 2026-03-26 | 8/400 |
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 085600材料與化工306 +10 | z1z2z3879 2026-03-21 | 11/550 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
[考研] 296求調(diào)劑 +4 | 汪!?! 2026-03-25 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開(kāi)開(kāi)心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|