一种基于分层知识的不透水面信息提取方法
【技术领域】
[0001] 本发明涉及透水面丰度提取技术领域，尤其涉及一种利用分层知识的提取中等分 辨率遥感影像不透水面丰度方法。
【背景技术】
[0002] 城市化进程正在全球范围内以前所未有的速度进行着，这一过程导致大量自然地 表逐渐向人工不透水面转变。不透水面是城市化程度的重要标志，同时也是环境科学、生态 学、地理学等领域关注的热点问题。
[0003] 目前，利用遥感影像，特别是中低分辨率影像的不透水面丰度提取方法普遍存在 用单一像元或高反射率(high albedo，H)和低反射率(low albedo，L)两类端元表示全部不 透水面类型，对不透水面类内差异考虑不足。此外，部分不透水面和部分土壤等地物存在相 似的光谱特征。较小的类间差异和较大的类内差异导致不透水面丰度提取存在在高丰度地 区低估、低丰度地区高估的问题。譬如，以往线性光谱混合分析(Linear Spectral Mixture Analysis，LSMA)及基于LSMA的不透水面丰度提取方法并没有根据城市生物物理组成的不 同而采取不同的端元模型，而是在整个影响范围内提取一个端元集应用于整个影像。这种 方法的缺点是没有考虑地物的类内差异，同时不同地物的光谱相似性也会导致丰度提取结 果的高估或低估。
【发明内容】
[0004] 本发明要解决的技术问题是提供一种基于分层知识的不透水面丰度提取方法，能 够改善由于相似地物的类间差异和不透水面的类内差异导致的不透水面信息在高丰度区 域低估、低丰度区域高估的现象，提高不透水面丰度信息提取精度。
[0005] 为解决上述问题，本发明采用如下的技术方案：
[0006] -种基于分层知识的不透水面信息提取方法包括以下步骤:步骤S1、基于城市生 物物理组成的影像分层，具体包括以下步骤：
[0007] 步骤S1.1、对城市遥感影像进行预处理，所述预处理包含辐射校正与去除水体掩 膜；
[0008] 步骤S1.2、计算预处理后的城市遥感影像生物物理组成的特征分量，对典型地物 在特征分量的特性进行分析；
[0009] 步骤S1.3、增强所述特征分量；
[0010]步骤S1.4、根据所述增强后特征分量，构建分层规则，将城市遥感影像分为不同层 次；
[0011] 步骤S2、基于LSMA的不透水面丰度提取，具体过程包括以下步骤：
[0012] 步骤S2.1、对各层分别进行MNF变换，分离数据中的噪声；
[0013] 步骤S2.2、在每一层的MNF散点图中选取端元，根据每一层的生物物理组成采用不 同类型的端元集；
[0014] 步骤S2.3、根据所述端元集信息在每一层分别进行LSMA分解，高反射率Η和低反射 率L丰度相加得到每一层不透水面丰度信息，通过对得到的每一层不透水面信息进行合并， 提取整个区域内不透水面丰度信息。
[0015] 2、作为优选，步骤S1.2具体为:在预处理后的图像中，采用IHS变换亮度分量I及归 一化植被指数NDVI作为表征城市地物构成的生物物理特征分量;其中，I表达同一地物的光 谱差异，其定义为：
[0016]
[0017]其中，R，G，B分别是影像中3个不同波段，以不透水面为高亮区域并和其他地物有 最大区分度为选取准则。
[0018]
[0019]其中，NIR为近红外波段，RED为红光波段。
[0020]作为优选，步骤S1.3具体为:对I、NDVI特征分量增强，其数学表达式为：
[0021]
[0022] 其中，ienh为增强后的特征分量，iη。!·为归一化的原始特征分量，λ为灵敏度因子，Θ 为归一化特征分量要增强的目标地物的平均值。
[0023] 作为优选，步骤S1.4具体为:在对特征分量进行增强的基础上，分析其直方图，确 定增强后亮度I分量阈值Ρ1、Ρ3及增强后NDVI分量阈值Ρ2、Ρ4，构建图像分成规则，其定义 为：
[0024] Area 1: Ienh>Pl，NDVIenh〈P2，
[0025] Area 3:Ienh〈P3，NDVIenh>P4，
[0026] Area 2为剩余部分。
[0027]作为优选，步骤S2.2具体为:Area 1选取高反射率-低反射率-土壤端元集，Area 2 选取高反射率-低反射率-土壤-植被端元集，Area 3选取低反射率-土壤-植被端元集。
[0028] 本发明的技术方案以城市遥感影像为主要数据源，根据区域的生物物理组成将研 究区分层，对各层独立采用不同的混合像元端元模型，利用LSMA提取区域不透水面信息，即 通过提取生物物理组成特征分量，构建分层规则，将研究区分层，再针对不同层的特征分别 独立选取端元集进行混合像元分解，得到不透水面丰度信息。
[0029] 本发明与现有技术相比具有如下特点：充分利用城市生物物理组成对遥感影像进 行分层，并且充分利用层级信息在影像每一层分别独立选取不同类型的端元进行混合像元 分解，从而增强不透水面类内差异、相似地物类间差异，提高不透水面丰度提取精度。
【附图说明】
[0030] 图1本发明基于分层知识的不透水面丰度提取方法流程图；
[0031 ]图2为遥感影像分层流程示意图；
[0032]图3为特征分量增强示意图，其中，A:亮度分量直方图;B:增强后亮度分量直方图； C:归一化植被指数直方图;D:增强后归一化植被指数直方图；
[0033] 图4为不透水面丰度计算流程示意图；
[0034] 图5为不透水面丰度信息提取结果示意图，其中，图5a:本发明的方法得到的不透 水面丰度;图5b:简单四端元LSMA得到的不透水面丰度。
【具体实施方式】
[0035]本发明提供一种基于分层知识的不透水面丰度提取方法，根据表征生物物理组成 的特征分量将图像分层，再根据每层地物构成特点独立选取端元集，提高了不透水面丰度 提取的精度，如图1所示，包括以下步骤，
[0036]步骤S1、基于城市生物物理组成的影像分层
[0037]城市不同生物物理成分由于构成地物不同具有明显的光谱特性差异，这种差异是 对图像进行分层的基础;本发明采用增强后的城市生物物理特征分量构建图像分层规则， 将图像分为不同生物物理层次，如图2所示，具体包括以下步骤：
[0038]步骤S1.1、对城市遥感影像进行预处理，所述预处理包含辐射校正与去除水体掩 膜。
[0039]对于城市遥感影像(主要针对Landsat卫星的中等分辨率影像)进行辐射定标，即 将遥感影像DN值转换为地面反射率，同时，水体虽然是城市重要地物类型，但根据V-I-S (Vegetation-Impervious Surface-Soil，植被-不透水面-土壤)模型，城市影像所有像元 均可视为植被、不透水面、土壤这三类的线性组合。因此，为消除水体影响，利用非监督分类 对图像进行去水体掩膜。本实施例中，利用ENVI软件辐射定标工具进行定标，获取地表反射 率；利用ISODATA(Iterative Self-Organizing Data Analysis Technique，迭代自组织数 据分析技术)非监督分类进行去水体掩膜。
[0040]步骤S1.2、计算预处理后的城市遥感影像生物物理组成的特征分量，对地物在特 征分量的特性进行分析。
[0041 ] 在预处理后的图像中，采用IHS( Intensity，Hue，Saturation)变换亮度分量I (Intensity)及归一化植被指数NDVI(Normalized Difference Vegetation Index)作为表 征城市地物构成的生物物理特征分量;其中，I表达同一地物的光谱差异，其定义为：
[0042]
[0043]其中，R，G，B分别是影像中3个不同波段，以不透水面为高亮区域并和其他地物有 最大区分度为选取准则，本实