专利名称: |
一种城市不透水面遥感提取方法及系统 |
摘要: |
本发明涉及一种城市不透水面遥感提取方法及系统。所述方法包括:获取Landsat数据;对所述Landsat数据进行预处理,得到预处理后的遥感数据;根据所述遥感数据,分别计算NDUI指数、MNDWI指数和SAVI指数;将所述NDUI指数、所述MNDWI指数和所述SAVI指数进行拉伸处理,得到拉伸后的NDUI指数、拉伸后的MNDWI指数和拉伸后的SAVI指数;根据所述拉伸后的NDUI指数、所述拉伸后的MNDWI指数和所述拉伸后的SAVI指数,计算归一化差值城市综合指数;根据所述归一化差值城市综合指数采用阈值法,提取不透水面信息。本发明能够提高不透水面的提取精度。 |
专利类型: |
发明专利 |
申请人: |
北京师范大学 |
发明人: |
陈艳玲;宫阿都;陈云浩;雷添杰;曾婷婷;杨雨晴 |
专利状态: |
有效 |
申请日期: |
1900-01-20T00:00:00+0805 |
发布日期: |
1900-01-20T08:00:00+0805 |
申请号: |
CN201911416291.1 |
公开号: |
CN111122449A |
代理机构: |
北京高沃律师事务所 |
代理人: |
韩雪梅 |
分类号: |
G01N21/17;G;G01;G01N;G01N21;G01N21/17 |
申请人地址: |
100875 北京市海淀区新街口外大街19号 |
主权项: |
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数据进行辐射定标、大气校正预处理,得到预处理后的遥感数据,所述遥感数据包括:蓝光的反射率、近红外的反射率、第二短波红外的反射率、绿光的反射率、红光的反射率和第一短波红外的反射率。 8.根据权利要求6所述的城市不透水面遥感提取系统,其特征在于,所述各指数计算模块,具体包括: 各指数计算单元,用于根据所述遥感数据采用公式和分别计算NDUI指数、MNDWI指数和SAVI指数; 其中,BLUE为蓝光的反射率,NIR为近红外的反射率,SWIR2为第二短波红外的反射率,GREEN为绿光的反射率,RED为红光的反射率,SWIR1为第一短波红外的反射率,l为土壤调节系数。 9.根据权利要求6所述的城市不透水面遥感提取系统,其特征在于,所述归一化差值城市综合指数计算模块,具体包括: 归一化差值城市综合指数计算单元,用于根据所述拉伸后的NDUI指数、所述拉伸后的MNDWI指数和所述拉伸后的SAVI指数采用公式计算归一化差值城市综合指数; 其中,NDUII为归一化差值城市综合指数,NDUI*为拉伸后的NDUI指数,MNDWI*为拉伸后的MNDWI指数,SAVI*为拉伸后的SAVI指数。 10.根据权利要求6所述的城市不透水面遥感提取系统,其特征在于,所述提取模块,具体包括: 阈值确定单元,用于采用目视解译和人工选取相结合的方式,确定阈值; 提取单元,用于根据所述阈值将所述归一化差值城市综合指数进行二值化处理,得到不透水面信息。 |
所属类别: |
发明专利 |