[發明專利]一種淺海大型復雜沙波區地形重構方法

有效

• 申請號: 201510071764.4
• 申請日: 2015-02-11
• 公開（公告）號: CN104613945B
• 公開（公告）日: 2017-01-18
• 發明（設計）人: 張華國;傅斌;何謝鍇;厲冬玲;史愛琴;樓琇林
• 申請（專利權）人: 國家海洋局第二海洋研究所
• 主分類號: G01C13/00
• 分類號: G01C13/00
• 代理公司: 浙江英普律師事務所33238
• 代理人: 陳小良

技術領域

本發明屬于遙感技術應用和海底地形測繪領域，具體是一種能夠利用遙感圖像中所含的沙波沙脊特征信息，如所希望那樣對以間隔型測線方式獲取的淺海大型復雜沙波區多波束測深數據進行有效地形重構的方法。

背景技術

海底地形測繪是探索和研究海洋的第一步，是開發和保護海洋的首要條件，也是現今海洋學研究的迫切需求。在回聲測深儀發明之前，主要靠測深桿和測深錘來測量水深，測量精度較差。20世紀20年代回聲測深儀出現之后，現代意義上海圖的繪制才得以實現。但早期的測深儀為單波束發射，一次發射只能得到測量船正下方的水深，因而只能實現點、線測量，無法反映測線之間的地形地貌。1970年代中期出現的多波束測深技術實現了帶狀測量(最大寬度可達水深的7倍)，顯著提高了測深效率。盡管如此，限于測量周期長、人力消耗大和資金需求高等方面的劣勢，多波束測深仍然通常采用間隔型測線方式進行。獲得的水深數據呈條帶狀分布，可用于大范圍淺海地形地貌特征統計分析。而要構建全覆蓋的淺海地形數字高程模型(DEM)須通過反距離權重法、Kriging插值法和樣條插值法等空間數據插值方法等空間插值方法實現。但由于條帶狀多波束測深數據的不連續性，僅通過上述空間插值方法無法構建連續的淺海地在地形，因此僅依賴多波束測深方法仍然很難滿足海底地形測繪需求。

已公開的發明專利(ZL201310043013.2)提出了一種基于遙感圖像特征的淺海沙波區多波束測深地形重構方法，利用遙感技術方法的大范圍覆蓋和高頻率重訪的能力，基于遙感圖像中所包含的因沙波地形對水動力調制而引起的灰度差異特征信息(亮暗條紋信息)，結合間隔型測線方式測量獲取的多波束測深數據，解決了淺海沙波區的地形重構難題。但該發明方法僅對沙波形態簡單且分布規律的小區域沙波區的地形重構能取得較好效果，而對復雜沙波區進行整體地形重構時，由于沙波形態各異、空間布局復雜，仍然無法滿足高精度地形重構的要求。

本發明專門針對淺海大型復雜沙波區地形重構的需要，通過基于沙波信息的分區組合和重疊連續賦值處理，解決了淺海大型復雜沙波區地形重構的技術難題。

發明內容

本發明的目的是提供了一種新的淺海大型復雜沙波區地形重構方法。

本發明通過以下技術方案實現：

一種淺海大型復雜沙波區地形重構方法，其特征在于其包括以下步驟：

(1)獲得多波束測深數據覆蓋區域含有沙波信息的遙感圖像，采用幾何校正方法進行遙感圖像定位配準，實現與多波束測深數據空間基本配準；

(2)根據遙感圖像上沙波紋理分布特征，提取沙波波脊；

(3)根據提取的沙波波脊線的分布數值，確定工作區范圍；

(4)評估多波束測深數據點的分布密度，在沒有多波束測深數據分布的區域均勻生成一定數量的輔助點，生成工作區輔助點；

(5)在每一條沙波波脊線兩端點分別繪制垂直于該沙波波脊線的輔助線，并延長至最近的其他沙波波脊線或工作區邊界線；

(6)用全部沙波波脊線和輔助線將工作區分割成若干分區；

(7)對每一個分區進行子區劃分，所有分區的子區劃分大小保持一致；子區劃分方法為：將每兩條相鄰的沙波波脊線之間的區域劃分為若干個子區；

(8)輔助點水深估算賦值；

(9)將輔助點的水深值和多波束測深數據點水深值一起進行空間插值，獲得工作區沙波地形數字高程模型，完成淺海沙波地形重構。

作為優選，上述一種淺海大型復雜沙波區地形重構方法的步驟(3)中所述的工作區范圍設置為矩形，矩形的四條邊離沙波波脊線的最小距離在100m-200m之間，同時使得矩形的兩條邊盡可能平行于沙波波脊線走向，另兩條邊盡可能垂直于沙波波脊線走向。

作為優選，上述一種淺海大型復雜沙波區地形重構方法的步驟(5)中所述的輔助線垂直于該沙波波脊線，然后分別向兩側延伸到最近的其他沙波波脊線或工作區邊界線。

作為優選，上述一種淺海大型復雜沙波區地形重構方法的步驟(8)中所述的輔助點水深估算賦值方法是：添加的每個輔助點為圓心、以多波束測線間距L為半徑設置搜索范圍；確定該輔助點所在的步驟(7)所劃分的子區，將與該子區及全部相鄰子區作為目標區域，搜索目標區域內的多波束測深數據點；采用反距離權重法確定該輔助點的水深值，即把搜索到的多波束測深數據點的水深值以距離平方的倒數為權重值疊加給該輔助點，疊加公式為：