DEM的建立及其在林業(yè)上的應(yīng)用
2017-09-06 06:19吳顯橋
現(xiàn)代農(nóng)業(yè)科技 2017年14期
吳顯橋
摘要 利用高程矢量數(shù)據(jù)，在Arcgis地理信息系統(tǒng)中建立不規(guī)則三角網(wǎng)TIN及數(shù)字高程模型DEM，運(yùn)用DEM進(jìn)行三維地理信息表達(dá)，制作不同類型的專題圖，獲取相關(guān)的信息和數(shù)據(jù)，滿足林業(yè)調(diào)查規(guī)劃設(shè)計、生態(tài)旅游景點開發(fā)等生產(chǎn)、經(jīng)營活動中反映地形地貌以及地面物體的形狀和逼真性的需要。文中較為完整地闡述了數(shù)字高程模型DEM從數(shù)據(jù)源準(zhǔn)備到建立再到應(yīng)用的全過程，并著重展示了在林業(yè)方面的應(yīng)用。
關(guān)鍵詞 DEM；矢量化；數(shù)字高程模型；不規(guī)則三角網(wǎng)；三維地形可視化
中圖分類號 S718.5 文獻(xiàn)標(biāo)識碼 A 文章編號 1007-5739（2017）14-0149-04
Abstract Using elevation vector data，TIN triangulated irregular network and digital elevation model DEM were established in Arcgis geographic information system.Using DEM 3D geographic information expression，thematic map of different types of access were made to relevant information and data，met the forestry investigation and planning design，ecological tourist attractions development production，reflected the needs of topography and the shapes of ground objects and realistic operating activities.This paper is a complete exposition of the digital elevation model DEM from the data source to the establishment to the whole process of application，and mainly presents the application in forestry.
Key words DEM；vectorization；DTM；TIN；3D terrain visualization
三維地理信息的產(chǎn)生和發(fā)展，是GIS技術(shù)及其應(yīng)用發(fā)展到一定水平的必然要求，它不僅局限于利用計算機(jī)技術(shù)手段對地理信息進(jìn)行可視化表達(dá)及其空間查詢，而且滿足了自動、全面顯示地理信息的要求[1]。在三維地理信息的表達(dá)中，DEM是一種獲取相對容易、具有較高精度的有效方式。DEM（digital elevation model）是用一組有序數(shù)值陣列形式表示地面高程的一種實體地面模型，形成柵格結(jié)構(gòu)的數(shù)據(jù)集，是地形曲面的數(shù)字化表達(dá)[2]，是數(shù)字地面模型DTM（digital terrain model）的一種。DEM包含規(guī)則格網(wǎng)Grid和不規(guī)則三角網(wǎng)TIN等表面模型，TIN指用一系列無重疊的三角形近似模擬陸地表面，從而構(gòu)成不規(guī)則的三角網(wǎng)[3]，使用TIN能更好地體現(xiàn)地形起伏變化特征。DEM直觀反映了現(xiàn)實的地貌情況，是現(xiàn)實世界山川、河流、地面起伏在計算機(jī)中的數(shù)字化表達(dá)。由于DEM描述的是地面高程信息，它在測繪、水文、氣象、地質(zhì)、地貌、土壤、工程建設(shè)、通訊、軍事等國民經(jīng)濟(jì)和國防建設(shè)以及人文和自然科學(xué)領(lǐng)域有著廣泛的應(yīng)用[4]。
本文通過對掃描地形圖等高線矢量化，進(jìn)而生成不規(guī)則三角網(wǎng)TIN及數(shù)字高程模型DEM，用實例說明了利用DEM制作三維景觀圖、電子林相沙盤和提取調(diào)查小班平均海拔、坡度、坡向值，旨在探討、優(yōu)化DEM的獲取途徑和創(chuàng)建方法，并對DEM在林業(yè)上的應(yīng)用進(jìn)行探索與展望。
1 獲取高程矢量數(shù)據(jù)
1.1 DEM獲取的途徑
目前，數(shù)字高程模型DEM的獲取主要有以下3種方式[5]：一是通過航空攝影測量獲取地面影像，并借助于干涉雷達(dá)和激光掃描儀等新型傳感器快速獲取高精度、高分辨率的DEM數(shù)據(jù)源；二是利用掃描地形圖獲取DEM數(shù)據(jù)源；三是通過全站儀、GPS等儀器配合計算機(jī)在野外進(jìn)行觀測獲取地面點數(shù)據(jù)，然后進(jìn)行轉(zhuǎn)換處理獲得DEM數(shù)據(jù)源。第1種方式需要擁有先進(jìn)的設(shè)備，一般的單位和個人不具備這個條件；第3種方式工作量大，費(fèi)用高，不適合用于大規(guī)模的數(shù)據(jù)采集任務(wù)；第2種方式獲取比較容易，又能滿足精度要求，適合用于一般的單位和個人。本文DEM數(shù)據(jù)的獲取采用第2種方式。
1.2 DEM數(shù)據(jù)源的提取
DEM數(shù)據(jù)源就是高程矢量數(shù)據(jù)，對于如何從一幅掃描地形圖中獲取等高線矢量數(shù)據(jù)的方法較多，有基于AutoC-AD、MAPGis、CASSCAn、Geoway、R2V、vpstudio、Arcgis（Arcs-can模塊）等軟件的矢量化方法，這些矢量化軟件各有優(yōu)缺點。R2V（Raster2Vector）是Windows環(huán)境下一款高級光柵圖矢量化軟件系統(tǒng)，具有快速矢量化、高程批量賦值和坐標(biāo)校正功能，以及良好的適應(yīng)性與高精確度，非常適合于GIS、地形圖、CAD及科學(xué)計算等應(yīng)用[6]。
1.3 等高線矢量化
1.3.1 裁剪柵格地形圖。柵格地形圖的裁剪需要使用Arcgis地理信息系統(tǒng)，ArcGis 是美國環(huán)境系統(tǒng)研究所（ESRI）研制的一套桌面地理信息系統(tǒng)軟件，ArcGis Desktop主要包含Arcmap、ArcCatalog、ArcToolbox三大桌面應(yīng)用，還有許多功能擴(kuò)展模塊。本文Arcgis的操作均以10.2中文版為例。
打開Arcmap10.2，在ArcToolbox中點擊數(shù)據(jù)管理工具→柵格→柵格處理→裁剪，用預(yù)定的一定范圍的面狀矢量數(shù)據(jù)裁剪出需要的柵格地形圖（一般采用5 m或10 m等高距的地形圖效果較好），并導(dǎo)出為tif格式圖像備用（本文以湖南省鳳凰縣南華山國家森林公園為例）。endprint
1.3.2 柵格地形圖的處理。將導(dǎo)出的tif格式地形圖用photoshop軟件進(jìn)行處理，去掉周邊的黑色斑塊，同時可以將等高線以外不必要的點狀、面狀、線狀地物以及文字、數(shù)字等標(biāo)注清除，以提高在R2V中等高線矢量化的速度和質(zhì)量。1.3.3 R2V矢量化。啟動R2V，加載經(jīng)過處理的柵格地形圖，用光標(biāo)拉動編輯窗口邊框和滾動條可移動到圖像的不同位置，按住鍵盤數(shù)字區(qū)的“+”、“-”鍵可放大、縮小圖像。如果掃描圖像質(zhì)量好，且內(nèi)容單一（以線條為主），可以選擇 “自動矢量化”；如果圖像比較復(fù)雜，有各種圖像元素混和在一起，只能使用R2V的交互跟蹤功能（即系統(tǒng)自動跟蹤配合手工處理）。若有多個編輯圖層，則需在“編輯”菜單中先定義編輯圖層。線條開始跟蹤時，點擊菜單“編輯”→“編輯線段”→“編輯線開/關(guān)”，激活線編輯器，用鼠標(biāo)在需要跟蹤矢量化的等高線上的任意位置點擊，即可進(jìn)入編輯跟蹤狀態(tài)。當(dāng)一條線跟蹤矢量化完成后，按
等高線賦值完成后，為了將生成的矢量數(shù)據(jù)轉(zhuǎn)換到特定的投影坐標(biāo)系統(tǒng)中（與裁剪柵格地形圖的坐標(biāo)系保持一致），點擊菜單“文件”→“載入控制點/TFW”，選擇裁剪柵格地形圖中后綴為tfw的文件作為控制點文件。最后將矢量數(shù)據(jù)導(dǎo)出為shp格式文件備用[6]。如果有高程點需要矢量化，可以建立點圖層進(jìn)行編輯，操作步驟與線圖層編輯大致相同。
2 建立不規(guī)則三角網(wǎng)TIN及數(shù)字高程模型DEM
2.1 創(chuàng)建不規(guī)則三角網(wǎng)TIN
啟動Arcmap，加載從R2V中導(dǎo)出的高程矢量數(shù)據(jù)并進(jìn)行拓?fù)錂z查，主要是消除數(shù)據(jù)中存在的懸掛點等問題[8]。拓?fù)錂z查處理后，重新加載高程矢量數(shù)據(jù)以及裁剪柵格地形圖，在統(tǒng)一的投影坐標(biāo)系中兩者疊加后位置基本對應(yīng)、吻合，能夠滿足精度要求。如果發(fā)現(xiàn)矢量化線條與柵格地形圖等高線的偏移過大不符合精度要求，可在Arcmap中用空間校正功能選擇地形圖與矢量圖同一位置的點作為控制點進(jìn)行校正。
在ArcToolbox中點擊3DAnalyst工具→數(shù)據(jù)管理→Tin→創(chuàng)建Tin（如3DAnalyst未能正常啟動，可先點擊Arc-map菜單“自定義”→“擴(kuò)展模塊”→勾選3DAnalyst前面的方框），在彈出的對話框中“輸入要素類”選擇等高線矢量數(shù)據(jù)，“高度字段”“Tag Field”均選擇高程字段，“SF Type”選擇“Hard_Line”（選擇“硬”隔斷線能更好地模擬地表形態(tài)），選擇原用投影坐標(biāo)系，確定Tin保存位置，點擊“確定”按鈕即可生成不規(guī)則三角網(wǎng)TIN[9]，如圖2所示。在“編輯Tin”中還可根據(jù)需要用面狀矢量數(shù)據(jù)（shp格式）對生成的Tin進(jìn)行裁剪。
2.2 創(chuàng)建數(shù)字高程模型DEM
在ArcToolbox中點擊3DAnalyst工具→轉(zhuǎn)換→由Tin轉(zhuǎn)出→Tin轉(zhuǎn)柵格，在彈出的對話框中“輸入Tin”選擇已生成的Tin，“輸出數(shù)據(jù)類型”選擇“FLOAT”，“方法”“采樣距離”“Z因子”可以選擇默認(rèn)設(shè)置（如果想進(jìn)一步提高DEM的精度，可在“方法”選項中選擇“NATURAL_NEIGHBOBS”，“采樣距離”選擇“CELSIZE…”，像元值大小與裁剪柵格地形圖的像元值大小保持一致，“Z因子”用默認(rèn)設(shè)置），確定DEM保存位置，點擊“確定”按鈕即可生成數(shù)字高程模型DEM[9]，如圖3所示。在ArcToolbox中點擊spatial Analyst工具→表面分析→山體陰影，在彈出的對話框中“輸入柵格”選擇已生成的DEM，在“輸出柵格”中選擇山體投影后DEM的保存位置，“方位角”“高度角”“Z因子”選擇默認(rèn)設(shè)置，點擊“確定”按鈕還可以生成山體投影DEM[7]。
對生成DEM精度的檢驗，在缺乏有效理論支持的情況下，采用回放DEM等高線與原柵格地形圖等高線疊置法是目前國際上較為廣泛采用的一種評價方法[10]。
3 數(shù)字高程模型DEM的應(yīng)用
以數(shù)字高程模型DEM為基礎(chǔ)，可以制作不同類型的專題圖。用DEM制作的三維景觀圖可以很好地展示地形地貌以及標(biāo)志性景點的位置、形狀等，制作的電子林相沙盤具有直觀、逼真、通視的效果，用DEM提取調(diào)查小班立地因子的數(shù)值快捷、高效、準(zhǔn)確率高，下面用實例進(jìn)行說明。
3.1 三維景觀圖的制作
啟動Arcscene10.2，加載已生成的南華山DEM以及道路、河流、景點、TIN裁剪范圍等在Arcmap中編輯好的矢量數(shù)據(jù)（其中裁剪范圍數(shù)據(jù)必須是線狀數(shù)據(jù)），右鍵南華山DEM→屬性，在彈出的對話框中“基本高度”選項“從表面獲取的高程”欄下點選“在自定義表面上浮動”并選擇南華山TIN，在“從要素獲取的高程”欄下選擇“沒有基于要素的高度”，“用于將圖層高程值轉(zhuǎn)換為場景單位的系數(shù)”后選擇“自定義”，變換系數(shù)值可根據(jù)需要的效果輸入具體數(shù)值，然后點擊“確定”按鈕。按照以上操作步驟，對道路、河流、景點、裁剪范圍線等矢量數(shù)據(jù)進(jìn)行同樣的設(shè)置，尤其是變換系數(shù)值必須保持一致。同時，對裁剪范圍的線狀矢量數(shù)據(jù)進(jìn)行如下設(shè)置：右鍵線狀矢量數(shù)據(jù)→屬性，在彈出的對話框中選擇“拉伸”選項，勾選“拉伸圖層中的要素”，“拉伸值或表達(dá)式”可根據(jù)需求效果輸入具體數(shù)值或選擇線圖層中的字段值，“拉伸方式”選擇“將其用作要素的拉伸數(shù)值”，這樣可形成剖斷面，同時對景點進(jìn)行符號編輯啟用3D標(biāo)記符號并選擇下載好的3D模型，最后生成的三維景觀圖[11]如圖4所示。如果疊加裁剪的相應(yīng)范圍衛(wèi)星影像圖還可以生成更接近現(xiàn)實表面的三維景觀圖。通過Arcscene的場景導(dǎo)航和飛行功能還可以獲取不同視角和高度的景觀效果。
3.2 電子林相沙盤的制作
在Arcmap中對林相圖的點、線、面矢量數(shù)據(jù)按彩色林相圖的要求進(jìn)行著色和標(biāo)注，然后導(dǎo)出為jpg或tif圖像（范圍與生成的Tin保持一致）。對導(dǎo)出的圖像進(jìn)行坐標(biāo)配準(zhǔn)，使之與Tin坐標(biāo)系保持一致。啟動Arcscene，加載Tin和配準(zhǔn)好的林相圖，右鍵林相圖→屬性，在彈出的對話框中“基本高度”選項“從表面獲取的高程”欄下點選“在自定義表面上浮動”并選擇“Tin”，在“從要素獲取的高程”欄下選擇“沒有基于要素的高度”，“用于將圖層高程值轉(zhuǎn)換為場景單位的系數(shù)”后選擇“自定義”，變換系數(shù)值可根據(jù)需要的效果輸入具體數(shù)值。在“渲染”選項中“效果”欄下勾選“相對于場景的光照位置為面要素創(chuàng)建陰影”和“如果可能使用平滑陰影”，同時將林相圖的優(yōu)先顯示級別設(shè)置為1，點擊“確定”按鈕。對Tin進(jìn)行同樣的設(shè)置，優(yōu)先顯示級別設(shè)置為2，點擊“確定”按鈕即可制作出具有林相分布的地理模型的三維圖像，稱之為電子林相沙盤[12]，如圖5所示。
3.3 提取調(diào)查小班平均海拔、坡度、坡向值
數(shù)字高程模型DEM 可以自動提取地形因子，包括高程、高程差、平均高程、等高線、坡度、坡向、坡度變化率、地面形態(tài)、地形剖面、地表粗糙度、脊谷線、溝谷密度、溝谷深度、太陽輻射強(qiáng)度及可視區(qū)等信息[13]。在森林資源調(diào)查中可用DEM來提取調(diào)查小班的平均海拔、坡度、坡向值。提取小班平均海拔的步驟：①打開ArcMap，加載DEM 數(shù)據(jù)和小班面shp數(shù)據(jù)；②點擊ArcToolbox→Spatial Analyst工具→區(qū)域分析→以表格顯示分區(qū)統(tǒng)計，在對話框中“輸入柵格數(shù)據(jù)或要素區(qū)域數(shù)據(jù)”欄選擇小班面shp 數(shù)據(jù)，在“區(qū)域字段”欄中選擇小班屬性表中的小班關(guān)鍵字（如小班號），在“輸入賦值柵格”欄中選擇相應(yīng)的DEM數(shù)據(jù)，在“輸出表”欄指定統(tǒng)計表輸出名稱及路徑，在“統(tǒng)計類型”欄中選擇MEAN（平均值）。設(shè)置完成后，點擊“確定”按鈕，生成小班平均海拔表；③分區(qū)統(tǒng)計后生成的小班平均海拔表是一個獨(dú)立的表，還需要建立基于公共字段的屬性連接，將小班平均海拔賦值給小班面的屬性表。操作步驟：點擊 ArcToolbox→數(shù)據(jù)管理工具→連接→連接字段，在“輸入表”欄中選擇小班面，在“輸入連接字段”欄中選擇小班面中確定的公共字段，在“連接表”欄中選擇小班平均海拔表，在“輸出連接字段”欄中選擇小班平均海拔表中的公共字段，在“連接字段”欄中選擇“平均海拔”字段，點擊“確定”按鈕，即可將小班海拔表中的平均海拔自動寫入小班面屬性表中。
提取平均坡度、坡向，先要通過ArcToolbox→Spatial Analyst工具→表面分析→坡度（“輸出測量單位”選擇DEGREE）和ArcToolbox→Spatial Analyst工具→表面分析→坡向，分別生成坡度、坡向柵格圖，在“區(qū)域分析”中按平均值生成小班平均坡度表和坡向表，然后再通過掛接和賦值將提取的數(shù)值直接導(dǎo)入小班面的屬性表中，方法與平均海拔的提取大致相同。最后在小班面的屬性表中還需要將坡向值按坡向柵格圖數(shù)值范圍對應(yīng)的方向替換成文字。
通過數(shù)字高程模型DEM提取調(diào)查小班的平均海拔、坡度、坡向值，準(zhǔn)確率達(dá)到85%以上，符合森林資源二類調(diào)查精度要求，可以適當(dāng)減少內(nèi)業(yè)工作量，提高工作效率[2]。用上述方法對南華山國家森林公園4個小班的平均海拔、坡度、坡向值提取結(jié)果如圖6、表1所示。
4 結(jié)語
數(shù)字高程模型DEM數(shù)據(jù)的獲取有多種途徑，柵格地形圖等高線矢量化方法也多種多樣，在實際應(yīng)用中應(yīng)根據(jù)本單位和個人的現(xiàn)有條件及經(jīng)濟(jì)狀況來決定采用哪種方式方法。數(shù)字高程模型DEM的應(yīng)用主要表現(xiàn)在4個方面[10]：地形分析方面；體積、面積計算方面；選址與規(guī)劃方面；管理應(yīng)用方面。該模型DEM在林業(yè)上具有廣泛的用途，除上面提到的應(yīng)用實例外，在林區(qū)建設(shè)填挖土石方計算以及水土保持、水源涵養(yǎng)、更新造林規(guī)劃設(shè)計等方面也有較高的應(yīng)用價值，通過可視域分析建立森林防火可視域效果圖[14]，能夠為森林防火撲救和指揮、決策提供科學(xué)依據(jù)。隨著虛擬現(xiàn)實技術(shù)VR的不斷發(fā)展，應(yīng)用三維可視化技術(shù)實現(xiàn)動態(tài)仿真變成可能[15]，將DEM描述的場景坐標(biāo)系轉(zhuǎn)換成為視點坐標(biāo)系，建立立體地形圖，結(jié)合數(shù)字林相圖的著色和配準(zhǔn)以及樹種處理，還可以合成虛擬林相圖，相當(dāng)逼真地表達(dá)出地形地貌和森林、活立木的景觀特點[16]。數(shù)字高程模型DEM對真實場景的表現(xiàn)畢竟有其局限性，通過疊加數(shù)字正射影像DOM（經(jīng)過處理糾正的航片或衛(wèi)星影像圖）、矢量圖層和導(dǎo)入Arcgis軟件自帶或SketchUp、3dsmax等軟件創(chuàng)建的三維模型，可進(jìn)一步增強(qiáng)地形可視化的空間立體感[17]。基于DEM的三維空間分析，由于受到分析窗口“近視”的制約，無法完整展示特定比例下的地形全景，在一定程度上限制了其應(yīng)用的廣度與深度，亟待從技術(shù)上解決、突破[10]。
5 參考文獻(xiàn)
[1] 康紅霞.基于Arcgis的三維景觀建模建設(shè)研究[D].西安：西安科技大學(xué)，2006.
[2] 李希飛.基于 DEM 的森林資源二類調(diào)查中小班立地因子提取方法的探討[J].內(nèi)蒙古林業(yè)調(diào)查設(shè)計，2016，39（2）：13-15.
[3] KANG-TSUNG C.地理信息系統(tǒng)導(dǎo)論[M].3版.陳健飛，譯.北京：清華大學(xué)出版社，2009.
[4] 苗險峰.淺談DEM的制作及應(yīng)用[J].中國地名，2013（3）：73-74.
[5] 劉志平.基于Arcgis的三維地貌影像圖生成方法與應(yīng)用研究：以大冶市為例[D].武漢：湖北大學(xué)，2010.
[6] 百度文庫.12 圖形矢量化及R2V使用方法簡介[EB/OL].[2017-03-20].http：//www.docin.com/p-1582677049.html.
[7] 許沖，莫艷合，劉燦.使用R2V、AutoCAD與ArcGIS建立DEM的方法研究與實踐[J].工程地質(zhì)計算機(jī)應(yīng)用，2008（2）：1-8.
[8] 郝紅科，田江博，陳永貴.用地形圖生成DEM的新方法研究[J].陜西林業(yè)科技，2005（1）：35-37.
[9] 劉志平，張素華，杜啟勝，等.基于ArcGIS的DEM生成方法及應(yīng)用[J].地理空間信息，2009，7（5）：69-71.
[10] 高英杰，潘建君，劉世峰.1：10000大比例尺DEM的制作及其地學(xué)應(yīng)用研究[J].安徽農(nóng)業(yè)科學(xué)，2009，37（2）：894-898.
[11] 湯國安，楊昕.ARCGIS地理信息系統(tǒng)空間分析實驗教程[M].北京：科學(xué)出版社，2006.
[12] 唐允森，王萍，李玉堂.利用 ArcScene 制作電子林相沙盤的初步實踐[C]//首屆吉林省林業(yè)學(xué)術(shù)大會論文集.長春：首屆吉林省林業(yè)學(xué)術(shù)大會，2006：90-92.
[13] 付曉，周維.數(shù)字高程模型及其在森林資源調(diào)查中的應(yīng)用[J].廣西林業(yè)科學(xué)，2002，31（1）：12-15.
[14] 王佳璆，張貴，肖化順.ArcGIS可視域分析在瞭望臺管理中的應(yīng)用[J].湖南林業(yè)科技，2005，32（2）：24-26.
[15] 楊光華，陳濤.基于ERDAS的三維地形可視化及其應(yīng)用[J].西華師范大學(xué)學(xué)報（自然科學(xué)版），2008，29（3）：307-312.
[16] 羅傳文.帽兒山虛擬林相的構(gòu)建[J].東北林業(yè)大學(xué)學(xué)報，2002，30（5）：40-43.
[17] 劉陽，張培松，韋仕川，等.基于ArcGIS三維地形可視化及其應(yīng)用研究：以陽江農(nóng)場為例[J].廣西農(nóng)業(yè)科學(xué)，2009，40（6）：772-776.endprint
