[發明專利]沖洪積扇信息的自動提取方法

有效

• 申請號: 201110191094.1
• 申請日: 2011-07-08
• 公開（公告）號: CN102231190A
• 公開（公告）日: 2011-11-02
• 發明（設計）人: 馮光勝;楊樹文;高山;顧湘生;李小和;趙新益;曹柏樹
• 申請（專利權）人: 中鐵第四勘察設計院集團有限公司
• 主分類號: G06K9/46
• 分類號: G06K9/46;G06K9/36;G06T7/00
• 代理公司: 武漢開元知識產權代理有限公司 42104
• 代理人: 黃行軍

1.一種沖洪積扇信息的自動提取方法，包括如下步驟：

(1)計算遙感影像中地物的遙感提取指數、植被指數和濕度指數，并分別生成遙感提取指數灰度圖像、植被指數灰度圖像和濕度指數灰度圖像；

(2)將所述遙感提取指數灰度圖像、所述植被指數灰度圖像和所述濕度指數灰度圖像合成彩色合成圖像；

(3)對所述彩色合成圖像進行主成分變換并提取第一主成分生成第一主成分灰度圖像；

(4)對所述第一主成分灰度圖像進行閾值分割，得到沖洪積扇的二值圖像，從而提取出沖洪積扇信息。

2.如權利要求1所述的沖洪積扇信息的自動提取方法，其特征在于，所述步驟(1)之前還包括：

對通過衛星獲得的遙感影像的五個波段進行大氣校正預處理。

3.如權利要求2所述的沖洪積扇信息的自動提取方法，其特征在于，所述衛星為TM衛星，遙感影像的所述五個波段為近紅外波段、紅光波段、綠光波段、藍光波段和中紅外波段，遙感影像為沖積扇的遙感影像或者洪積扇的遙感影像。

4.如權利要求1所述的沖洪積扇信息的自動提取方法，其特征在于，所述步驟(1)具體為：

(11)獲取通過衛星檢測的地物的紅光波長、藍光波長、近紅外波長和綠光波長，并根據公式計算遙感影像中地物的遙感提取指數并生成遙感提取指數灰度圖像，

其中，I_AFI為遙感提取指數，λ_RED為遙感影像中地物的紅光波長，λ_NIR為遙感影像中地物的近紅外波長，λ_BLUE為遙感影像中地物的藍光波長，λ_CREEN為遙感影像中地物的綠光波長；

(12)獲取通過衛星檢測的地物的近紅外波段反射率和紅光波段反射率，并根據公式計算遙感影像中地物的歸一化植被指數并生成植被指數灰度圖像，

其中，I_NDVI為歸一化植被指數，ρ_NIR為遙感影像中地物的近紅外波段反射率，ρ_RED為遙感影像中地物的紅光波段反射率；

(13)獲取通過衛星檢測的地物的中紅外波段反射率和綠光波段反射率，并根據公式計算遙感影像中地物的濕度指數并生成濕度指數灰度圖像，

其中，I_NDMI為土壤亮度指數，ρ_GREEN為遙感影像中地物的綠光波段反射率，ρ_MIR為遙感影像中地物的中紅外波段反射率。

5.如權利要求1所述的沖洪積扇信息的自動提取方法，其特征在于，所述步驟(4)具體為：

(41)對所述第一主成分灰度圖像進行灰度值拉伸；

(42)提取拉伸后的所述第一主成分灰度圖像的灰度直方圖，并通過所述灰度直方圖計算分割閾值；

(43)根據所述分割閾值，將所述第一主成分灰度圖像進行二值化，得到沖洪積扇的二值圖像，從而提取出沖洪積扇信息。

6.如權利要求5所述的沖洪積扇信息的自動提取方法，其特征在于，所述步驟(41)具體為：

將所述第一主成分灰度圖像中的灰度值范圍拉伸到0-255。

7.如權利要求5所述的沖洪積扇信息的自動提取方法，其特征在于，所述步驟(42)具體為：

提取拉伸后的所述第一主成分灰度圖像的灰度直方圖，通過迭代法計算所述灰度直方圖的最大灰度值和最小灰度值，并選擇兩個最大灰度值之間的最小灰度值作為分割閾值。

8.如權利要求1所述的沖洪積扇信息的自動提取方法，其特征在于，所述步驟(4)之后還包括：

對所述沖洪積扇的二值圖像進行形態學膨脹濾波、腐蝕濾波以及去噪處理。