衛(wèi)星遙感在貴州喀斯特地域研究中的應(yīng)用
——以興義萬(wàn)峰林為例

楊 娟，宋善海，劉 蕓

(貴州省生態(tài)氣象和衛(wèi)星遙感中心，貴州 貴陽(yáng) 550002)

1 引言

萬(wàn)峰林位于貴州省興義市境內(nèi)，是中國(guó)西南三大喀斯特地貌之一[1]。景區(qū)長(zhǎng)200多km，寬30～50 km，從海拔2000多m的七捧高原邊緣延伸，與萬(wàn)峰湖北岸、黃泥河?xùn)|岸成扇形展開，連綿至安龍、貞豐等地[2]。萬(wàn)峰林景區(qū)內(nèi)峰巒疊嶂，峰林中錯(cuò)落分布著田園小河、布依小寨，人與自然和諧共融。萬(wàn)峰林的神奇與秀美，形成一道天下罕見的峰林畫廊，被不少專家和游人譽(yù)為“天下奇觀”[3]。

2 研究區(qū)域概況

萬(wàn)峰林夏無(wú)酷暑，冬無(wú)嚴(yán)寒，雨量充沛，日照長(zhǎng)，屬于低緯度高海拔地帶，具有亞熱帶季風(fēng)氣候特征。1月份平均氣溫4.5 ℃，7月份平均氣溫26.8 ℃，年平均氣溫14～19 ℃，降雨量1300～1600 mm，無(wú)霜期300 d左右。

圖1 萬(wàn)峰林風(fēng)貌

萬(wàn)峰林景區(qū)山峰密集奇特，氣勢(shì)宏大壯闊，整體造型完美，北部為峰林盆地，中南部為峰林洼地和峰叢山地[4]。從地質(zhì)學(xué)的角度看，萬(wàn)峰林峰林、峰叢大多呈錐形，部分為鐘狀、平頂狀和馬鞍狀，堪稱一座“中國(guó)錐狀喀斯特博物館”[1]。萬(wàn)峰林景區(qū)內(nèi)含有1個(gè)鄉(xiāng)、2個(gè)街道辦事處、5個(gè)鎮(zhèn)(表1)，研究區(qū)將景區(qū)范圍進(jìn)行適當(dāng)擴(kuò)展，納入研究區(qū)的面積為717.18 km2(圖2)。

圖2 萬(wàn)峰林景區(qū)及研究區(qū)范圍

表1 研究區(qū)域范圍及面積 km2

3 研究思路

3.1 數(shù)據(jù)來(lái)源

本研究所采用的數(shù)據(jù)為高分二號(hào)1 m空間分辨率的多光譜遙感影像數(shù)據(jù)和資源三號(hào)立體像對(duì)提取的柵格分辨率為10 m數(shù)字高程模型(圖3、表2)，利用ArcGIS軟件，采用空間數(shù)據(jù)分析、鄰域分析等分析方法，完成了基于規(guī)則格網(wǎng)DEM的山峰頂點(diǎn)的快速提取方法研究，并分析了各影像因素之間的相互關(guān)系，完成了研究區(qū)內(nèi)山峰個(gè)數(shù)的計(jì)算和分析評(píng)估。

表2 采用的遙感影像數(shù)據(jù)情況

圖3 研究區(qū)高分二號(hào)1 m分辨率遙感影像

3.2 技術(shù)方法

DEM數(shù)據(jù)包含了大量的地形信息，通過(guò)DEM數(shù)據(jù)可以完成對(duì)地形特征的提取，對(duì)數(shù)字地形分析具有重要作用。通過(guò)GIS高程圖分析，“峰林”是基部低矮，相互獨(dú)立的成群山體，而“峰叢”是基座相連的延續(xù)山峰，當(dāng)盆地被“峰林與峰叢”環(huán)繞時(shí)，就形成了喀斯特盆地。而位于喀斯特中心部位是叢谷，在該中心部位，當(dāng)“峰林與峰叢”狹長(zhǎng)排列，就形成喀斯特谷地；當(dāng)“峰林與峰叢”環(huán)繞封閉時(shí)，就形成喀斯特洼地[5]。通過(guò)對(duì)萬(wàn)峰林內(nèi)地域性特征的分析，將峰頭作為山頂點(diǎn)來(lái)判斷。山頂點(diǎn)是指在特定鄰域分析范圍內(nèi)，該點(diǎn)都比周圍點(diǎn)高的區(qū)域；凹陷點(diǎn)是指特定鄰域分析范圍內(nèi)，該點(diǎn)都比周圍點(diǎn)低的區(qū)域。本研究選用ArcGIS10.4軟件作為數(shù)據(jù)分析平臺(tái)，通過(guò)基于規(guī)則格網(wǎng)DEM典型地形特征山峰頂點(diǎn)的快速提取算法，即通過(guò)等高線、山頂點(diǎn)、凹陷點(diǎn)的計(jì)算和提取，利用等高線提取、鄰域分析和窗口計(jì)算等工具完成柵格數(shù)據(jù)的表面分析，再利用空間分析模塊中的等高線提取功能，完成等高距15 m和75 m的等高線圖繪制，并以此作為山頂點(diǎn)空間分布圖的背景圖，利用鄰域分析和柵格計(jì)算器工具完成山頂點(diǎn)提取，實(shí)現(xiàn)對(duì)萬(wàn)峰林景區(qū)內(nèi)峰頭個(gè)數(shù)的統(tǒng)計(jì)。

利用格網(wǎng)DEM數(shù)據(jù)提取山頂點(diǎn)時(shí)，不同的分析尺度可獲得不同層次的山頂點(diǎn)。首先，對(duì)分析窗口內(nèi)各變量賦值對(duì)原始 DEM 單元矩陣進(jìn)行掃描，設(shè)定一個(gè)鄰域分析的閾值，再對(duì)該分析區(qū)域進(jìn)行重分類，滿足條件的則判定為分析區(qū)域內(nèi)的山頂點(diǎn)，但是在實(shí)際驗(yàn)證過(guò)程中發(fā)現(xiàn)提取的山頂點(diǎn)存在偽噪音點(diǎn)[7]，因此，需要對(duì)山頂點(diǎn)進(jìn)行處理。封閉等高線可以作為輔助判別的山頂點(diǎn)的有效工具，因此采用一種基于形態(tài)分析的山頂點(diǎn)提取方法，可以快速?gòu)腄EM數(shù)據(jù)中提取山頂點(diǎn)，從而對(duì)這些數(shù)據(jù)的可靠性進(jìn)行客觀的評(píng)定，提高其準(zhǔn)確性。可將等高線的間距作為判定山頂點(diǎn)的高差閾值。從DEM數(shù)據(jù)中提取等高線，并將等高線圖轉(zhuǎn)化為面狀圖，將提取到的面狀區(qū)域與之前提取到的山頂點(diǎn)進(jìn)行疊加取交集，得到去除偽噪音點(diǎn)的山頂點(diǎn)圖(圖4)。

圖4 研究路線

3.3 提取步驟及結(jié)果

山頂點(diǎn)提取按照焦點(diǎn)統(tǒng)計(jì)→柵格計(jì)算→重分類→矢量轉(zhuǎn)換→提取山頂點(diǎn)高程的步驟完成。首先打開ArcGIS10.4軟件并加載DEM數(shù)據(jù)，可以將DEM數(shù)據(jù)進(jìn)行裁剪，加快計(jì)算速度([ArcToolbox]→[Spatial Analysis Tools]→[Extraction]→[Extract by Mask])；焦點(diǎn)統(tǒng)計(jì)是指提取分析窗口中柵格的最大值，并將其賦給焦點(diǎn)(待計(jì)算柵格中心)，選擇[Spatial Analyst Tools]→[Neighborhood]→[Focal Statistics]，提取大小為50×50的分析窗口中柵格的最大值；進(jìn)行柵格計(jì)算，提取山頂點(diǎn)，選擇[Spatial Analyst Tools]→[Map Algebra]→[Raster Calculator]，輸入公式，提取山頂點(diǎn)；對(duì)計(jì)算結(jié)果進(jìn)行重分類，選擇[Spatial Analyst Tools]→[Reclass]→[Reclassify]；將重分類結(jié)果轉(zhuǎn)換為矢量點(diǎn)，選擇[Conversion Tools]→[From Raster]→[Raster to Point]；求出所有山頂點(diǎn)的高程，選擇[Spatial Analyst Tools]→[Extraction]→[extract values to points]，輸入?yún)?shù)文件得到結(jié)果。

通過(guò)鄰域分析結(jié)合等高線面狀圖提取得到山頂點(diǎn)的分布圖，為了將偽噪音點(diǎn)去除，提取高差閾值為15 m和75 m的等高線面狀分布圖進(jìn)行疊置分析，最后得到實(shí)際峰頭分布圖(圖5、6)。

圖5 不同高差閾值的等高線面狀分布

基于以上原則，將研究區(qū)內(nèi)萬(wàn)峰林峰頭個(gè)數(shù)進(jìn)行計(jì)算、統(tǒng)計(jì)、分析，得到峰頭個(gè)數(shù)為8972個(gè)，其中景區(qū)峰頭個(gè)數(shù)為3309個(gè)(圖7)。

圖6 山頂點(diǎn)分布

圖7 研究區(qū)峰頭個(gè)數(shù)分析、統(tǒng)計(jì)與評(píng)估

4 結(jié)論與討論

基于規(guī)則格網(wǎng)的DEM在地貌分析中比較方便而且容易與遙感影像數(shù)據(jù)結(jié)合分析。本文依據(jù)山頂點(diǎn)的地貌學(xué)定義及其在DEM中的形態(tài)特征，采用ArcGIS山頂點(diǎn)提取算法，運(yùn)用山頂點(diǎn)高差閾值提取的等高線轉(zhuǎn)換為面狀圖進(jìn)行了疊置分析，從而快速提取了山頂點(diǎn)個(gè)數(shù)。通過(guò)分析與評(píng)估提取萬(wàn)峰林峰頭個(gè)數(shù)為8972個(gè)，其中景區(qū)峰頭個(gè)數(shù)為3309個(gè)。該研究結(jié)果豐富了DEM的分析內(nèi)容，對(duì)貴州喀斯特地域的研究有參考作用。但本研究還有待進(jìn)一步加深，如：可深入研究分辨率及分析窗口大小對(duì)不同地貌的山頂點(diǎn)提取的影響。