| 4 | 1/1 | 返回列表 |
| 查看: 668 | 回復(fù): 3 | |||
nirvana_fwl新蟲 (初入文壇)
|
[求助]
求助編程序,算雷擊大地密度,用網(wǎng)格法。
|
|
求熱心人幫忙,我需要計(jì)算地閃密度。把一片區(qū)域分割成經(jīng)緯度為0.1*0.1的網(wǎng)格,然后判斷每次閃電是否發(fā)生在這個(gè)網(wǎng)格里,如果是,就把發(fā)生在這個(gè)網(wǎng)格里的閃電次數(shù)累加起來,再除以面積,就得到每一個(gè)網(wǎng)格的地閃密度。我需要用程序?qū)崿F(xiàn)它,求大神們幫幫我,不甚感激。‘厴I(yè)論文用的。 感覺我好像沒講清楚,我的算法和下面我從一篇文獻(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。分別用來存放圖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)過計(jì)算得出每個(gè)格子的經(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個(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)語句重新對(duì)81個(gè)格子進(jìn)行循環(huán)賦值。使用公式為:平均電流強(qiáng)度= 總電流強(qiáng)度/閃擊次數(shù) Ng = 總地閃次數(shù)/資料年數(shù)/格子面積 經(jīng)過以上4步,即可得出每個(gè)格子的Ng值和雷 |
至尊木蟲 (著名寫手)
驃騎將軍
|
你看的論文挺早的吧,看起來像純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 不過看你的描述,感覺格子光保存中心點(diǎn)經(jīng)緯度不夠方便,因?yàn)橛袑挾?雖然每次可以計(jì)算,但是由于格子一旦初始化完成就固定了,最好還是把每個(gè)格子的經(jīng)緯度范圍保存起來作為Grid的屬性,用[x,y,w,h]這種頂點(diǎn)寬高法也行,用[x1,y1,x2,y2]這種兩點(diǎn)法也行 這樣程序就只需要用一個(gè)Gfid的二維數(shù)組來保存,9*9 |

新蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +9 | 路癡小琪 2026-03-27 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +3 | 丟掉懶惰 2026-03-27 | 6/300 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +7 | 大饞小子 2026-03-27 | 8/400 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 求調(diào)劑 +5 | 蘆lty 2026-03-25 | 6/300 |
|
|
[考研] 286求調(diào)劑 +13 | Faune 2026-03-21 | 13/650 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +6 | 谷雨上岸 2026-03-23 | 7/350 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 打過很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|