一种多光谱遥感反射率影像吸收特征提取方法
【技术领域】
[0001] 本发明涉及遥感图像处理技术领域，尤其涉及一种基于虚拟波段的多光谱遥感反 射率影像吸收特征提取方法。
【背景技术】
[0002] 遥感信息提取主要依赖于地物的光谱和空间特征的差异。由于技术条件的限制， 遥感传感器同时达到很高的空间分辨率和光谱分辨率有较大难度。目前，具有较高空间分 辨率的多光谱遥感传感器仍然是遥感发展的主流和热点，如World View3，Landsat8，ASTER 等等。但是，如何充分发掘多光谱遥感数据的光谱维信息，仍是亟待解决的难题。通常，多光 谱遥感光谱维信息采用主成分分析、波段比值等方法获取，但这些方法需要较多的背景知 识，并且无法逐波段单独提取吸收特征信息。
[0003] 包络线去除法是一种常用的光谱分析方法，它可以将光谱特征归一到一致的光谱 背景上，有利于和其他光谱曲线进行吸收特征参量数值的比较，从而为定量分析奠定基础。 1984年，包络线去除法由Clark首次提出，此后被广泛用于高光谱遥感影像处理以及点光谱 数据处理当中，成为光谱吸收特征提取最重要的方法。然而，传统的包络线去除法在处理多 光谱数据时，由于其方法需要选定首尾波段及中间凸集节点，并且在这些波段无法提取有 效吸收特征信息，从而导致多光谱遥感有限的光谱信息被极大浪费，特征提取效果不佳。
【发明内容】
[0004] 本发明是为了解决上述不足，提供了一种多光谱遥感反射率影像吸收特征提取方 法。
[0005] 本发明的上述目的通过以下的技术方案来实现:一种多光谱遥感反射率影像吸收 特征提取方法，其特征在于:包括以下步骤:根据原始多光谱遥感反射率影像波段设置生成 虚拟波段反射率影像数据；
[0006] 对虚拟波段反射率影像数据逐像素获取包络线；
[0007] 将虚拟波段反射率影像数据获取的包络线按原始遥感影像波段设置进行重采样；
[0008] 采用上述步骤获得的包络线对原始多光谱反射率影像数据进行包络线去除。
[0009] 优选地，所述多光谱遥感反射率影像为经过大气校正等预处理后得到的地表反射 率影像。
[0010] 优选地，所述根据原始多光谱遥感反射率影像波段设置生成虚拟波段反射率影像 具体包括：
[0011]在原始影像首波段之前生成虚拟首波段；
[0012] 在原始影像末波段之后生成虚拟末波段；
[0013] 在原始影像相邻波段之间生成虚拟中间波段；
[0014] 将虚拟首波段、虚拟中间波段以及虚拟末波段按照波长递增顺序组合在一起，生 成虚拟波段反射率影像。
[0015] 优选地，所述在原始影像首波段之前生成虚拟首波段具体包括：
[0016] 计算原始影像首波段与第二波段之间的中心波长间隔；
[0017] 在原始影像首波段短波方向，距其一半上述波长间隔位置，生成虚拟首波段，其各 像素取值与真实影像首波段完全相同。
[0018] 优选地，所述在原始影像末波段之后生成虚拟末波段具体包括：
[0019] 计算原始影像末波段与倒数第二波段之间的中心波长间隔；
[0020] 在原始影像末波段长波方向，距其一半上述波长间隔位置，生成虚拟末波段影像， 其各像素取值与真实影像末波段完全相同。
[0021] 优选地，所述在原始影像两两波段之间生成虚拟中间波段具体包括：
[0022] 取原始影像相邻波段中心波长平均值，作为虚拟中间波段中心波长；
[0023]取原始影像各像素相邻波段较大值，作为虚拟中间波段各像素取值。
[0024] 优选地，所述包络线为对虚拟波段光谱采用外壳系数法处理获得的背景光谱。
[0025] 优选地，所述重采样指按照目标波段设置对原始光谱进行线性插值。
[0026] 优选地，所述包络线去除指将原始光谱与包络线逐波段相除，获得包络线去除光 谱。
[0027] 本发明与现有技术相比的优点是:通过采用本发明公开的一种多光谱遥感反射率 影像吸收特征提取方法，克服了传统包络线去除法在多光谱反射率数据波段吸收特征提取 当中的局限性，提出根据原始多光谱遥感影像波段设置生成虚拟波段数据，采用虚拟波段 数据逐像素生成包络线，并重采样到原始波段设置完成包络线去除，从而实现了多光谱遥 感反射率数据各波段吸收特征定量提取，可以进而实现对矿物、植被、土壤、水体等地面目 标成分定量化研究。
【附图说明】
[0028]图1是本发明的流程图。
[0029] 图2是5种矿物的平均反射率光谱。
[0030] 图3是经过处理后五种矿物的平均光谱。
【具体实施方式】
[0031] 下面结合实施例对本发明进一步详述。
[0032] -种多光谱遥感反射率影像吸收特征提取方法，如图1所示，包括以下步骤：
[0033] 步骤101，根据原始多光谱遥感反射率影像波段设置生成虚拟波段反射率影像数 据；
[0034] 步骤102,对虚拟波段反射率影像数据逐像素获取包络线；
[0035] 步骤103,将虚拟波段反射率影像数据获取的包络线按原始遥感影像波段设置进 行重采样；
[0036] 步骤104,采用上述步骤获得的包络线对原始多光谱反射率影像数据进行包络线 去除。
[0037] 本实施例中Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)采集到的遥感影像数据经大气校正等预处理被转换成地表反射率影像 后，首先基于其波段设置生成虚拟波段反射率影像。然后，对虚拟波段反射率影像数据逐像 素获取包络线，并将其按原始波段设置进行重采样，得到实际使用的包络线。最后，对地表 反射率影像逐像素进行包络线去除处理，从而提取出吸收特征信息。
[0038]本实施例中所用的遥感影像数据为Terra卫星搭载的ASTER传感器（Advanced Spaceborne Thermal Emission and Reflection Radiometer)米集的数据，使用的波段为 可见近红外的3个波段(空间分辨率15m)，以及短波红外的6个波段(空间分辨率30米），共9 个波段数据。将该数据短波红外波段重采样到15m分辨率，与可见近红外波段组合在一起进 行大气校正，得到反射率影像。该反射率影像各波段中心波长如表1所示。
[0039]表1所用ASTER遥感反射率影像各波段中心波长
[0040]
[00411本发明实施例中选择的ASTER遥感影像获取于2005年4月30日，实验区域是位于智 利北部的丘基卡玛塔铜矿。王桂珍等人曾使用该地区数据开展矿物蚀变信息提取方法研 究，为与实际验证结果进行比较，本方法也使用该区域数据进行吸收特征提取。
[0042]将选取的ASTER L1B产品转换成地表反射率影像。
[0043] 在影像上根据之前研究结果，选取白云母、伊利石、高岭石、绿泥石以及黄钾铁矾 等五种矿物的感兴趣区，作为已知样本，用于验证吸收特征提取效果。图2展示了这5种矿物 的平均反射率光谱。
[0044] 本发明实例中，基于ASTER多光谱反射率影像，生成虚拟首波段、虚拟末波段以及8 个虚拟中间波段，并按波长顺序组合在一起，生成虚拟波段反射率影像，共计