經(jīng)過觀察可以發(fā)現(xiàn)某種珍貴植物在山區(qū)的某個區(qū)域位置生長狀況明顯要比其他區(qū)域好很多，通過研究了解到這種植物生長具有嚴格的生長條件。為了能更好地保護該種植物的生長環(huán)境，現(xiàn)在需要使用GIS空間分析方法，將適合植物生長的區(qū)域找出來。
需要的數(shù)據(jù)包括該山區(qū)的等高線數(shù)據(jù)和山區(qū)若干觀測點采集的年平均溫度和年總降水數(shù)據(jù)。植物生長條件主要包括4個：
(1)這種植物一般生長在河流兩側(cè)500m的區(qū)域內(nèi)；
(2)這種植物喜陰；
(3)這種植物適宜生長氣候環(huán)境為年平均溫度為9°—11°；
(4)這種植物適宜的年總降水量為580—710mm。
下面小編說一說該選址的整體工作思路哈，為了簡潔直觀明了，小編直接上技術(shù)流程圖：
一、提取陰坡數(shù)據(jù)
（1）提取坡向數(shù)據(jù)。選擇【Spatial Analyst工具】|【表面分析】|【坡向】工具，【輸入柵格】選擇“DEM”，【輸出柵格】設(shè)置為“aspect”，根據(jù)需要，設(shè)置其他參數(shù)，點擊【確定】，生成坡向數(shù)據(jù)“aspect”，如下圖所示：
（2）提取陰坡數(shù)據(jù)。選擇【Spatial Analyst 工具】|【重分類】|【重分類】工具,【輸入柵格】選擇“aspect”，【值字段】選擇“Value”，點擊【分類】，類別選擇3，分類方法選擇“手動”，在右邊的框輸入90，270。
（3）這時候由于我們要提取的是陰坡（坡向為270°—360°和0°—90°），所以小編把90°—270°陽坡的范圍新值設(shè)置為“NoData”，其他的新值都設(shè)置成1，如下圖所示：
（4）【輸出柵格】設(shè)置為“yinpo”，點擊確定，并且將最終輸出結(jié)果疊加到“DEM”數(shù)據(jù)上，如下圖所示：
二、提取河流緩沖區(qū)
（1）填洼。選擇【Spatial Analyst工具】|【水文分析】|【填洼】工具，【輸入表面柵格數(shù)據(jù)】選擇“DEM”，【輸出表面柵格】設(shè)置為“fill”，【Z限制】意思是要填充的匯與其傾瀉點之間的最大高程差。默認情況下將填充所有匯(不考慮深度)。
（2）提取水流流向。選擇【Spatial Analyst工具】|【水文分析】|【流向】工具，【輸入表面柵格數(shù)據(jù)】選擇“fill”，【輸出表面柵格】設(shè)置為“dir”，其他參數(shù)視具體情況來設(shè)置，本次小編都默認哦，結(jié)果圖層如下圖所示：
（3）提取水流流量。選擇【Spatial Analyst工具】|【水文分析】|【流量】工具，【輸入流向柵格數(shù)據(jù)】選擇“fill”，【輸出儲積柵格數(shù)據(jù)】設(shè)置為“sum”，其他參數(shù)默認即可。
（4）提取河流。這時候可以看見河流的模樣了，可是怎么提取出來呢？選擇【Spatial Analyst 工具】|【地圖代數(shù)】|【柵格計算器】，輸入一個判斷語句：“sum>300”，輸出文件名為“stream”，有的小伙伴就會問了，為什么是大于300？不是大于250或者大于1000，這個其實是看你研究情況來定的，各位小伙伴可以通過符號化，漸變分級把想要的河網(wǎng)給盡可能顯示出來，然后通過使用識別工具（藍色感嘆號）比較多個河流最上游的像元值大小，選擇最小的像元值作為判斷的閾值,即可把你想要的河流提取出來，結(jié)果如下圖所示（值為1的都是河流，背景值為0）：
（5）對河流做緩沖區(qū)分析。右鍵剛才提取的河流數(shù)據(jù)“stream”選擇屬性表，選擇“Value”為1的行，也就是選中了河流。這時候選擇【轉(zhuǎn)換工具】|【自柵格轉(zhuǎn)換】|【柵格轉(zhuǎn)線】，將剛才的數(shù)據(jù)轉(zhuǎn)換成矢量線，輸出結(jié)果數(shù)據(jù)為“Buffer”，接著使用菜單欄“地理處理”下面的緩沖區(qū)工具，相關(guān)系數(shù)設(shè)置如下圖所示：
三、提取合適的溫度和降水區(qū)域
（1）加載觀測點。由于觀測點為文本信息，首先將其轉(zhuǎn)換成要素，以便后續(xù)分析。在ArcMap主菜單上選擇【文件】【添加數(shù)據(jù)】【添加XY數(shù)據(jù)】，【X字段】選擇“X”，【Y字段】選擇“Y”，點擊確定完成導(dǎo)入后，內(nèi)容列表將新增一個圖層，此時該圖層為事件，還不是要素類，許多高級編輯處理功能不能使用，因此要轉(zhuǎn)換成要素類。右擊新圖層，選擇【數(shù)據(jù)】【導(dǎo)出數(shù)據(jù)】，將觀測點數(shù)據(jù)保存為“dot”。
（2）對觀測點進行插值分析。由于只有局部點數(shù)據(jù)，要生成表面數(shù)據(jù)需要使用插值方法內(nèi)插生成，本次選擇內(nèi)插方法“普通克里金”（記得使用插值方法前在環(huán)境里面設(shè)置處理范圍和柵格分析范圍與“DEM”數(shù)據(jù)相同）。選擇【3D Analyst 工具】|【柵格插值】|【克里金】工具，輸入點要素選擇“dot”，Z值字段分別選擇“溫度”和“降水”，輸出柵格分別設(shè)置為“wendu”和“jiangshui”，
其他參數(shù)設(shè)置默認即可。
溫度插值結(jié)果：
降水插值結(jié)果：
（3）根據(jù)植物生長條件提取合適的溫度和降水區(qū)域。打開【Spatial Analyst 工具】|【地圖代數(shù)】|【柵格計算器】，輸入一個判斷語句：("jiangyu">580)&("jiangyu"9)&("wendu"
上圖中黑色的是不符合條件的柵格區(qū)域，這時候我們希望只要藍色的，選擇【Spatial Analyst 工具】|【重分類】|【重分類】工具,【輸入柵格】選擇“aspect”，【值字段】選擇“Value”，將舊值為0的新值改成“NoData”，輸出文件名為“area”，點擊確定即可得到我們想要的藍色區(qū)域，如下圖所示：
（4）這時候還要疊加上500米的緩沖區(qū)。使用【Spatial Analyst Tools】|【提取分析】|【按掩膜提取】，第一欄輸入“area”，第二欄輸入數(shù)據(jù)“Buffer”，輸出文件名為“area1”，點擊確定得到結(jié)果圖層：
（5）最終我們?yōu)榱酥庇^地展示出該植物的生長區(qū)域，小編將其疊加到DEM數(shù)據(jù)上面，如下圖所示：
練習(xí)數(shù)據(jù)下載：植物生長區(qū)評估
本文鏈接：http://www.lajichun.cn/blog/128.html
本文標簽：DEM
