多元特征線遙感識別的灘涂DEM構建

楊立君，王得玉，姜 杰，李文梅，施佳鳳

(1. 南京郵電大學，江蘇 南京 210023； 2. 南京師范大學，江蘇 南京 210023)

由于受觀測環境、調査成本等條件的制約，灘涂地形測量研究長期受到國內外學者關注。1997年Mason等[1-3]假設水邊線是一條等高線，通過疊加不同潮汐條件下的水邊線來生成灘涂數字高程模型，并分析了高程與植被、泥沙粒徑之間的關系；2007年鄭宗生等[4-5]依據多時相遙感影像，利用水邊線構建了崇明東灘潮間帶DEM，并用兩條實測剖面線進行了DEM精度檢查。2014年Li Zhen等[6]利用水邊線生成了年度潮間帶地形圖，實現了德國瓦登海北部灘涂的發展情況定量估算和監測。2014年丁賢榮等[7-8]]通過潮位測站基線計算水邊線高程，利用水邊線構建了江蘇沿海輻射沙脊群DEM。2016年Xu Zhen等[9]借助水邊線方法研究了季節變化對GOMSO灣灘涂地形的影響。上述灘涂地形的遙感測量方法稱之為水邊線復合技術，主要是利用遙感技術提取水邊線，并且將其作為等高線結合高程數據反演灘涂地形。該技術水邊線數量越多、分布越均勻、影像時間跨度越短，地形反演精度越高，反之越差。水邊線復合技術是反演灘涂地形的常用方法，但傳統的水邊線復合技術有效范圍具有局限性且需要大量遙感影像。遙感影像數量越多經濟成本越高，水邊線時間跨度一般也會越長。灘涂地形變化頻繁長時間序列的水邊線反演灘涂地形精度降低，獲取的灘涂地形也丟失了現勢性。挖掘遙感影像上其他高程信息是解決上述問題的重要途徑。

灘涂地區的水邊線、沖刷痕跡線、植被發展下界及植被間的分界線在遙感影像上均清晰可見。這些線狀目標分布在灘涂不同的特征空間，具有顯著的空間分布規律，利用該分布規律可獲取多元線狀目標高程值。本文結合專家經驗知識提取遙感影像線狀目標，提出一種基于線狀地物空間分布規律反演灘涂地形的遙感技術方法，旨在提高水邊線復合技術反演精度與作業范圍。

1 數據與方法

1.1 研究區與數據

1.1.1 研究區

選擇資料條件較為完善的上海崇明東灘為研究區(圖1中箭頭指向的矩形區域)，在中等尺度下開展灘涂地形反演方法研究。崇明東灘是由長江北支和南支兩條水道的徑流，以及漲、落潮流所挾帶的泥沙在崇明島東部逐漸沉積而形成，灘內潮溝密布，高、中、低潮灘分帶十分明顯，屬于典型的淤泥質潮灘[10]。崇明東灘位于長江入?？冢蔽徊钶^大，不同潮時的遙感影像，水邊線高程差異明顯，便于構建垂直差異明顯的灘涂高程模型。

圖1 研究區位置

1.1.2 數據及預處理

研究數據主要有影像數據、光譜數據、點位數據及上海市灘涂專題地形圖。影像數據由7張衛星影像組成，詳細情況見表1。光譜數據是2011年8月應用FieldSpec?3 Hi-Res光譜儀采集不同植被、水體、裸灘等典型地物的光譜數據。點位數據則是應用Trimble5700 GPS測量了研究區3條高程控制剖面線。專題地形圖成圖時間為2007年，成圖比例尺為1∶5000，水上部分采用陸地測量方法，水下部分采用GPS、聲吶探測儀配合驗潮站進行水位改正的測量方法。點位數據的平面坐標系為1954北京坐標系，參考高程基準為吳淞高程。

表1 影像數據

數據預處理主要包括大氣校正、幾何糾正及影像增強等遙感影像的預處理工作。使用FLAASH模塊進行大氣校正，選用NLAPS處理產品的對應定標參數完成輻射定標，大氣模型選擇中緯度各季節，氣溶膠模型選擇海洋模型。在堤壩拐角及農田邊界選取11個特征點作為幾何糾正點，幾何糾正精度均小于半個像元。

1.2 研究方法

灘涂潮間帶高程對植被和沖刷痕跡線的空間分布具有明顯的控制作用，在遙感影像上水邊線、植被發展下邊界、植被間分界線及沖刷痕跡線相互平行，具有高程指示性，均可作為灘涂DEM反演信息源。從實測光譜數據可以看出，不同地物之間(不同植被、植被/陸地、陸地/水體)具有較為明顯的光譜特征差異，利用遙感影像能夠有效地提取各類邊界信息。沖刷痕跡線尺度小，在中等分辨率的衛星影像上并不足以顯示，但在高空間分辨率(如IKONOS)的衛星影像上顯示十分明顯(圖2箭頭所指的位置)。研究的主要思路是利用遙感技術從遙感影像上提取上述特征線，依據地物的空間分布規律空間賦值方法獲取特征線高程，應用GIS工具空間插值法構建灘涂DEM，最后，利用實測的點位數據和地形圖空間插值生成的DEM對反演結果進行精度評價。研究方法如圖3所示。

圖2 特征線遙感識別結果

圖3 研究方法

2 結果與分析

2.1 邊線特征遙感識別

2.1.1 水邊線提取

利用面向對象技術的水邊線提取方法提取水邊線[11-13]，該方法是利用TM6波段確定分類樣本，基于面向對象的思想分別利用灰度、幾何形狀及上下文信息實現水陸分離。由于利用ENVI二次開發語言編程實現水邊線提取，水邊線提取速度快、精度高、人工參與少，結果如圖2所示。

2.1.2 植被邊線提取

研究區植被按照高程由高到低依次為蘆葦、互花米草與海三棱藨草，灘涂濕地植被分類一直是遙感圖像分類的難點，植被之間的分界線提取較為復雜。利用密度分割法和神經網絡綜合分類方法[14]分別提取研究區植被發展下界和植被間邊界，保留分類植被下界(植被下界高程指示精度高)并矢量化，編輯后結果如圖2所示。

2.1.3 沖刷痕跡線提取

Canny算子邊緣檢測能力強[15]，結合局部自相關影像增強、輔助紋理信息邊緣檢測法實現沖刷痕跡線提取，提取結果如圖2所示。

2.2 特征線水平性分析

特征線理論上高程應一致為一條水平線，實際上復雜的地理環境致使特征線并非水平，特征線的水平性誤差直接影響DEM的構建質量。利用上海市灘涂專題地形圖制作與遙感影像相同格網尺寸的DEM數據，將特征線與DEM數據進行疊加，利用GIS工具提取特征線位置對應的DEM高程數據。計算每條特征線的最大高程(Hmax)、最小高程(Hmin)、高程差極大值(Hsub)、平均高程值(Hever)及高程標準差(SD)等信息，以此評價特征線水平性，結果見表2。從表中可以看出，沖刷痕跡線水平性最好，水邊線與植被邊線水平性沒有顯著性差別，高程精度滿足中小比例尺制圖要求。

表2 邊線水平誤差統計

2.3 特征線高程賦值

潮位校正法是水邊線高程賦值的常用方法，但是不能有效推算植被邊線、沖刷痕跡線高程。臨近插值法及擬合法均可以對特征線進行高程賦值。但臨近插值法精度較低，而擬合法應用高程剖面的測量點擬合高程趨勢線，通過擬合曲線獲得高程，能夠有效消除地形隨機波動的影響，精度較高[16]。因此，本文基于RTK技術實測研究區3條高程剖面線，用擬合法對特征線高程賦值。當有多條擬合曲線穿過特征線時，取高程平均值，表3為特征線高程計算結果。

表3 特征線高程值

注：“—”表示高程剖面未通過特征線，“WL”指水邊線。

2.4 灘涂DEM構建

將具有高程屬性的特征線作為數據源，應用ANUDEM空間插值方法[17]生成灘涂DEM。該方法集成在ArcGIS軟件中，利用Topo to Raster工具實現，圖4為灘涂DEM遙感反演結果。

圖4 灘涂DEM反演結果

2.5 DEM質量分析

2.5.1 質量評價

利用專題地形圖制作相同格網尺寸的DEM數據，作為參考數據。將反演結果與參考DEM作差值運算，獲得地形誤差分布圖，如圖5所示。其中，誤差平均值為-0.15 m，中誤差為±0.26 m。

利用式(1)評價反演DEM數據的質量。式中，Zi為反演DEM數據高程，Hi為對應的參考DEM數據高程。選取46個樣本點，計算樣本點數據中誤差為RMSE=±0.12 m。

(1)

圖5 誤差統計分布

2.5.2 誤差分析

影像幾何糾正誤差(平面誤差)，特征線起算高程誤差、特征線水平性誤差及潮溝等影響是構建灘涂DEM誤差的主要來源。多元線性特征遙感識別的DEM構建方法并未考慮潮溝的影響，用一個高程常量來表示特征線高程存在一定的誤差。對于大范圍的灘涂地形反演，特征線水平性誤差及潮溝的影響是主要的誤差來源。利用已知地形圖空間插值法構建參考DEM含有空間插值誤差，但該評價方法能夠從全局了解誤差分布情況，掌握誤差分布特征。利用式(1)評價反演DEM，樣本數據有意避開潮溝，不包含DEM空間插值帶來的誤差僅包含測量誤差，評價質量可信度高。從誤差統計分布圖中可以看出，構建DEM誤差主要集中在潮溝豐富的低潮灘區域、植被與裸灘交匯地帶。前者是由于本文構建灘涂DEM并未考慮潮溝的影響，后者是由于參考DEM誤差較大的原因造成的。由于專題地形圖植被區域采用陸地測量方法，裸灘區域采用水下測量方法，兩種測量方法的交會地帶成圖精度較低，因此，植被與裸灘交會地帶誤差大，并不可信，利用式(1)評價DEM精度較高這一點即可證明。

3 結 論

以水邊線、植被邊線及沖淤痕跡線等異質空間特征線為基礎，結合實測數據擬合線法高程賦值，借助遙感與地理信息新技術構建灘涂DEM，在理論上可行，精度可信，結果符合灘涂地形基本規律，解決了灘涂地區實地測量困難問題。該方法實現了特征線從二維到三維的轉換，除大比例尺工程測量的精度要求外，基本滿足潮灘演變、沖淤計算等基礎研究的精度要求，可作為主要數據資料使用。基于多元線性特征遙感識別的DEM構建方法，本質上是水邊線復合技術的拓展，通過引入植被邊界線及沖淤痕跡線解決了植被潮灘地形遙感反演問題，并且降低了生產成本，但還有以下問題值得說明：① 各種對地觀測衛星捕捉大潮高、低潮位線困難，反演范圍受限，需進一步挖掘單幅遙感影像高程信息；②對較大研究區域，特征線水平性誤差較大，應考慮特征線水平性誤差消除方法，提高反演精度;③不同空間分辨率的遙感影像，空間地理信息量不同，高程反演精度也不一樣，遙感影像空間分辨率對DEM的精度影響本文未進行研究，以后將作為進一步的研究方向。