[發(fā)明專利]多源異構遙感影像控制點自動采集方法

有效

• 申請?zhí)枺?/td>: 201310143369.3
• 申請日: 2013-04-23
• 公開（公告）號: CN103218787A
• 公開（公告）日: 2013-07-24
• 發(fā)明（設計）人: 唐新明;王光輝;王伶俐;王華斌;李參海;樊文峰
• 申請（專利權）人: 國家測繪地理信息局衛(wèi)星測繪應用中心
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 北京市柳沈律師事務所 11105
• 代理人: 劉貝

1.一種多源異構遙感影像控制點自動采集方法，其特征在于，包括：

從經過幾何精糾正處理的多個數據源獲取遙感影像；

對所述獲取的遙感影像進行分析歸納，統(tǒng)計并分析其屬性信息，該屬性信息包括影像格式、影像分辨率、坐標系統(tǒng)、成圖比例尺、影像時相；

對所述獲取的遙感影像進行優(yōu)化格網設計，根據影像分辨率、幅寬將影像分為不同尺寸的格網；

利用Wallis變換對每一個格網內的影像進行預處理；

對經Wallis變換處理后的影像進行增強，綜合運用Moravec算子、Harris算子、Forstner算子、SUSAN算子、尺度不變特征及最穩(wěn)定極值區(qū)域檢測算法提取特征點；

以控制點影像片采集窗口來遍歷統(tǒng)計所述特征點的個數，選擇特征點數量最大的窗口為采集區(qū)域；

針對所述選擇的采集區(qū)域進行遙感控制點影像片的裁切保存，同時獲取其對應區(qū)域的數字高程模型數據。

2.根據權利要求1所述的影像控制點自動采集方法，其特征在于，所述數據源包括航空影像、資源系列衛(wèi)星影像、SPOT系列衛(wèi)星影像、雷達衛(wèi)星影像中的一個或多個。

3.根據權利要求1所述的影像控制點自動采集方法，其特征在于，所述Wallis變換包括：

將數字圖像分為互不重疊矩形區(qū)域，每區(qū)域的尺度對應于要增強的紋理模式的尺度；

計算每一塊區(qū)域的灰度均值和方差；

將灰度均值和方差的目標值分別設定為127和40-70之間的數值；

計算出每一塊所述區(qū)域的Wallis變換的乘性常數和加性常數；

由雙線性內插計算所述數字圖像的任一像素的系數，并計算變換后的灰度值。

4.根據權利要求1所述的影像控制點自動采集方法，其特征在于，所述對經Wallis變換處理后的影像進行增強包括：先利用Moravec算子、Harris算子、Forstner算子、SUSAN算子提取特征點，然后利用尺度不變特征變換及最穩(wěn)定極值區(qū)域算法來檢測所提取到的特征點，記錄特征點個數，并且利用描述算法對提取出的特征點進行具體描述。

5.根據權利要求4所述的影像控制點自動采集方法，其特征在于，所述Moravec算子的計算包括：

計算各格網區(qū)域影像的興趣值；

給定一個閾值，將所述計算的興趣值大于該閾值的點作為候選點；

選取所述候選點中的極值點作為特征點。

6.根據權利要求4所述的影像控制點自動采集方法，其特征在于，所述Harris算子的計算包括：

求出影像上所有像素點的梯度，即對每個像素點的灰度值進行一階差分運算；

確定一個n×n大小的窗口，生成n×n大小的高斯卷積模板，其中，高斯模板的σ取0.3-0.9之間的值；

利用生成的高斯模板對所述梯度值進行高斯濾波，并計算強度值；

選取所述強度值的局部極值點為特征點，如果所述強度值大于閾值0.6，則將該點取為特征點，并且進行排序保存。

7.根據權利要求4所述的影像控制點自動采集方法，其特征在于，所述Forstner算子的計算包括：

計算各像素的Robert梯度；

根據所計算的Robert梯度值來計算l×l窗口中灰度的協方差矩陣；

利用所計算的協方差矩陣來計算興趣值；

確定特征點，其中，如果興趣值大于給定的閾值，則該像素為特征點。

8.根據權利要求4所述的影像控制點自動采集方法，其特征在于，所述SUSAN算子的計算包括：

用SUSAN模板在圖像上滑動，將模板內部每個圖像像素點的灰度值分別與模板中心像素的灰度值作比較，計算模板中各點的相似度；

基于所計算的相似度來計算圖像中每一點的USAN區(qū)域大??；

基于所計算的USAN區(qū)域大小，用角點響應函數來計算角點初始響應；

通過所述USAN區(qū)域的質心與所述模板中心的距離去除偽角點；

在去除偽角點后，在通過非極大值抑制來尋找角點。