| 2 | 1/1 | 返回列表 |
| 查看: 560 | 回復(fù): 1 | |||
askding金蟲(chóng) (小有名氣)
|
[求助]
環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)及作圖請(qǐng)教
|
|
想做一個(gè)環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)的圖,使用VBA編程看的不是很懂,有人能幫助下嗎? 我試著編了下,求高人幫我看下編的對(duì)嗎? 如下:在評(píng)價(jià)區(qū)域內(nèi),將區(qū)域污染源可能影響區(qū)域14km×14km的范圍網(wǎng)格化,取步長(zhǎng)為500 m,形成29行29列的網(wǎng)格。在每個(gè)網(wǎng)格對(duì)7個(gè)源的擴(kuò)散后的風(fēng)險(xiǎn)值進(jìn)行迭加得到總和r。然后在Excel中將網(wǎng)格的數(shù)據(jù)轉(zhuǎn)成X/Y/Z格式。具體編程如下:Sub 環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)() Dim k As Integer, j As Integer, i As Integer, x As Integer, y As Integer, m As Integer, n As Integer, d1 As Integer, d2 As Integer, d3 As Integer, d4 As Integer, d5 As Integer, d6 As Integer, d7 As Integer, r1 As Double, r2 As Double, r3 As Double, r4 As Double, r5 As Double, r6 As Double, r7 As Double , rr As Double k=0'定義X軸最小值 For j= 2 To 30 Sheet1.Cells(1,j)=k k=k+50 '間隔為500m(將實(shí)際距離換算成相對(duì)坐標(biāo)系下的距離長(zhǎng)度,實(shí)際為500m,坐標(biāo)下為50m) Next k=1400 'Y軸最大值 For j=2 To 30 Sheet1.Cells(j,1)=k k=k-50 Next x=0 For m=2 To 30 y=1400 For n=2 To 30 Sheet1.Cells(n,m)=Sqr((705-x)^2+(652-y)^2)'輸入源1坐標(biāo) Sheet1.Cells(n+30,m)=Sqr((690-x)^2+(523-y)^2)'輸入源2坐標(biāo) Sheet1.Cells(n+60,m)=Sqr((769-x)^2+(416-y)^2)'輸入源3坐標(biāo) Sheet1.Cells(n+90,m)=Sqr((552-x)^2+(682-y)^2)'輸入源4坐標(biāo) Sheet1.Cells(n+120,m)=Sqr((758-x)^2+(815-y)^2)'輸入源5坐標(biāo) Sheet1.Cells(n+150,m)=Sqr((894-x)^2+(983-y)^2)'輸入源6坐標(biāo) Sheet1.Cells(n+180,m)=Sqr((667-x)^2+(381-y)^2)'輸入源7坐標(biāo) y=y-50 Next x=x+50 Next For m=2 To 30 For n=2 To 30 d1=Sheet1.Cells(n,m)'表示距離 d2=Sheet1.Cells(n+30,m) d3=Sheet1.Cells(n+60,m) d4=Sheet1.Cells(n+90,m) d5=Sheet1.Cells(n+120,m) d6=Sheet1.Cells(n+150,m) d7=Sheet1.Cells(n+180,m) r1=(400-d1)/350*8.33*0.00001'風(fēng)險(xiǎn)值計(jì)算,需輸入影響半徑m r2=(400-d2)/350*8.33*0.00001 r3=(400-d3)/ 350*8.33*0.00001 r4=(417-d4)/367*8.08*0.0001 r5=(250-d5)/200*9.16*0.00001 r6=(333-d6)/283*5.7*0.00001 r7=(400-d7)/350*3.72*0.001 Sheet2.Cells(n,m)=r1 If r1<0 Then r1=0 End If Sheet2.Cells(n+30,m)=r2 If r2<0 Then r2=0 End If Sheet2.Cells(n+60,m)=r3 If r3<0 Then r3=0 End If Sheet2.Cells(n+90,m)=r4 If r4<0 Then r4=0 End If Sheet2.Cells(n+120,m)=r5 If r5<0 Then r5=0 End If Sheet2.Cells(n+150,m)=r6 If r6<0 Then r6=0 End If Sheet2.Cells(n+180,m)=r7 If r7<0 Then r7=0 End If Sheet3.Cells(n,m)=rl+r2+r3+r4+r5+r6+r7 Next Next For m=2 To 30 For n=2 To 30 rr=Sheet3.Cells(n,m) If rr>0 Then Sheet5.Cells(n,m)=Int(Log(rr)/Log(10)+8)'轉(zhuǎn)化為小數(shù)形式,并存入表單5 Else Sheet5.Cells(n,m)=1 End If Next Next For m=2 To 30 For n=2 To 30 Sheet4.Cells((m-2)*29+n,3)=Sheet5.Cells(n,m)'以YXZ形式存入表單4 Next Next For m=2 To 30 For n=2 To 30 Sheet4.Cells((m-2)*29+n,2)=1400-50*(n-2) Next Next For m=2 To 30 For n=2 To 30 Sheet4.Cells((m-2)*29+n,l)=0+(m-2)*50 Next Next End Sub |
鐜 |
金蟲(chóng) (小有名氣)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚(yú)567 2026-03-24 | 4/200 |
|
|
[考研] 359求調(diào)劑 +3 | 王了個(gè)楠 2026-03-25 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 +9 | A-哆啦Z夢(mèng) 2026-03-23 | 15/750 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | 林之夕 2026-03-24 | 5/250 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 315分,誠(chéng)求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|