[發明專利]一種城市不透水面遙感提取方法及系統

在審

• 申請號: 201911416291.1
• 申請日: 2019-12-31
• 公開（公告）號: CN111122449A
• 公開（公告）日: 2020-05-08
• 發明（設計）人: 陳艷玲;宮阿都;陳云浩;雷添杰;曾婷婷;楊雨晴
• 申請（專利權）人: 北京師范大學;中國水利水電科學研究院
• 主分類號: G01N21/17
• 分類號: G01N21/17
• 代理公司: 北京高沃律師事務所 11569
• 代理人: 韓雪梅

1.一種城市不透水面遙感提取方法，其特征在于，包括：

獲取Landsat數據；

對所述Landsat數據進行預處理，得到預處理后的遙感數據；

根據所述遙感數據，分別計算NDUI指數、MNDWI指數和SAVI指數；

將所述NDUI指數、所述MNDWI指數和所述SAVI指數進行拉伸處理，得到拉伸后的NDUI指數、拉伸后的MNDWI指數和拉伸后的SAVI指數；

根據所述拉伸后的NDUI指數、所述拉伸后的MNDWI指數和所述拉伸后的SAVI指數，計算歸一化差值城市綜合指數；

根據所述歸一化差值城市綜合指數采用閾值法，提取不透水面信息。

2.根據權利要求1所述的城市不透水面遙感提取方法，其特征在于，所述對所述Landsat數據進行預處理，得到預處理后的遙感數據，具體包括：

對所述Landsat數據進行輻射定標、大氣校正預處理，得到預處理后的遙感數據，所述遙感數據包括：藍光的反射率、近紅外的反射率、第二短波紅外的反射率、綠光的反射率、紅光的反射率和第一短波紅外的反射率。

3.根據權利要求1所述的城市不透水面遙感提取方法，其特征在于，所述根據所述遙感數據，分別計算NDUI指數、MNDWI指數和SAVI指數，具體包括：

根據所述遙感數據采用公式和分別計算NDUI指數、MNDWI指數和SAVI指數；

其中，BLUE為藍光的反射率，NIR為近紅外的反射率，SWIR2為第二短波紅外的反射率，GREEN為綠光的反射率，RED為紅光的反射率，SWIR1為第一短波紅外的反射率，l為土壤調節系數。

4.根據權利要求1所述的城市不透水面遙感提取方法，其特征在于，所述根據所述拉伸后的NDUI指數、所述拉伸后的MNDWI指數和所述拉伸后的SAVI指數，計算歸一化差值城市綜合指數，具體包括：

根據所述拉伸后的NDUI指數、所述拉伸后的MNDWI指數和所述拉伸后的SAVI指數采用公式計算歸一化差值城市綜合指數；

其中，NDUII為歸一化差值城市綜合指數，NDUI^*為拉伸后的NDUI指數，MNDWI^*為拉伸后的MNDWI指數，SAVI^*為拉伸后的SAVI指數。

5.根據權利要求1所述的城市不透水面遙感提取方法，其特征在于，所述根據所述歸一化差值城市綜合指數采用閾值法，提取不透水面信息，具體包括：

采用目視解譯和人工選取相結合的方式，確定閾值；

根據所述閾值將所述歸一化差值城市綜合指數進行二值化處理，得到不透水面信息。

6.一種城市不透水面遙感提取系統，其特征在于，包括：

Landsat數據獲取模塊，用于獲取Landsat數據；

預處理模塊，用于對所述Landsat數據進行預處理，得到預處理后的遙感數據；

各指數計算模塊，用于根據所述遙感數據，分別計算NDUI指數、MNDWI指數和SAVI指數；

拉伸處理模塊，用于將所述NDUI指數、所述MNDWI指數和所述SAVI指數進行拉伸處理，得到拉伸后的NDUI指數、拉伸后的MNDWI指數和拉伸后的SAVI指數；

歸一化差值城市綜合指數計算模塊，用于根據所述拉伸后的NDUI指數、所述拉伸后的MNDWI指數和所述拉伸后的SAVI指數，計算歸一化差值城市綜合指數；

提取模塊，用于根據所述歸一化差值城市綜合指數采用閾值法，提取不透水面信息。

7.根據權利要求6所述的城市不透水面遙感提取系統，其特征在于，所述預處理模塊，具體包括：

預處理單元，用于對所述Landsat數據進行輻射定標、大氣校正預處理，得到預處理后的遙感數據，所述遙感數據包括：藍光的反射率、近紅外的反射率、第二短波紅外的反射率、綠光的反射率、紅光的反射率和第一短波紅外的反射率。