本发明专利技术公开了一种基于混合亚像元解混的遥感多光谱与全色图像融合方法，包括将原始低分辨率多光谱图像MS0采用cubic方式上采样到原始全色图像PAN的分辨率，得到MS；并基于该MS计算归一化植被指数NDVI；利用PAN和NDVI，对MS中的亚像元进行分类，得到植被-非植被混合亚像元CM和其他亚像元CO；利用PAN和NDVI，判断植被-非植被混合亚像元的地物类型类别；对于MS中每个亚像元t，根据其对应的亚像元类别和地物类型类别进行融合。本发明专利技术的有益效果：通过对多光谱图像中的混合亚像元，在融合过程中进行光谱解混处理，从而降低融合图像中混合像元数量，锐化融合图像边界和降低光谱失真。
【技术实现步骤摘要】
本专利技术涉及遥感图像处理
，具体来说，涉及一种基于混合亚像元解混的 遥感多光谱与全色图像融合方法。
技术介绍
遥感数据融合是将不同传感器数据所提供的信息加以综合，以获得更全面、准确 的图像，它能实现优势互补和消除冗余，克服单一传感器在几何、光谱、空间和时间分辨率 等方面的局限，从而有利于影像解译、分类和目标识别。随着航天、传感器和通信等技术的 发展，遥感数据在空间、光谱、时间分辨率上快速提高。由于传感器物理特性和数据传输能 力等的限制，高空间分辨率和高光谱空间分辨率影像难以同时获取。因此，大量目前在轨的 卫星（例如Landsat7ETM+、QuickBird、SP0T_5和WorldView-2/3等)会同时提供高空间分 辨率的全色波段（PAN)和低空间分辨率的多光谱波段（MS)。由于大量的应用需要使用高空 间分辨率的多光谱图像，近年来国内外研宄人员发展了大量遥感融合技术来融合MS和PAN 图像以得到高空间分辨率的MS图像。 现有的MS与PAN融合方法可以分为三类：基于成分替代的方法、基于调制的 方法和基于多尺度分析的方法。基于成分替代融合的典型代表算法有IHS变换、PCA， Gram-Schmidt等，而基于调制技术的融合算法主要有Brovey变换、Pradines'、合成变化 比率（SyntheticVariableRatio)、SmoothingFilter-basedIntensityModulation、 PANSHARP(PS)等。基于成分替代和调制技术的融合算法的特点是快速且易于实现，但其 融合结果会导致一定程度的光谱失真。基于多尺度变换的融合方法，有效地保持了多光谱 图像的光谱信息。然而这类方法的融合图像可能会出现空间畸变，典型的现象有振铃效应、 虚景混淆、边缘以及纹理模糊。为综合成分替代和多分辨率分析等方法分别在空间和光谱 信息保留方面的优势，一些研宄人员提出了成分替代（PCA、HIS、ICA等）或调制技术与多分 辨率分析相结合的融合算法，这类方法能获得优于成分替代和标准多分辨率分析融合的结 果，但增加了计算复杂度。有一些研宄利用正则化优化的方法对融合结果进行迭代。尽管 国内外研宄人员提出了大量新的融合方法来降低融合图像光谱失真，如何进一步降低光谱 失真仍然是遥感融合研宄面临的主要挑战。 由于MS和PAN图像空间分辨率的差异，低分辨率MS图像的一个混合像元在重采 样到PAN分辨率的MS图像中对应多个亚像元。这些亚像元中部分像元对应PAN图像中的 纯像元，另一部分则对应PAN图像中的边界像元。在现有融合方法中，对应PAN纯像元的MS 混合亚像元在融合图像中仍然为混合像元，这使得融合图像中植被和非植被像元的边界模 糊，存在较大的光谱失真。因此，为了改善融合图像质量，需要在融合过程中将MS混合亚像 元解混到对应PAN纯像元的地物类别。然而，现有融合方法鲜有考虑混合亚像元的解混。在 少量考虑了混合亚像元解混的融合方法中，主要使用面对对象分类的方法来获取地物类别 信息。但面向对象分类在参数设置、样本训练等方面的复杂性限制了这些方法的实际应用。 针对相关技术中的问题，目前尚未提出有效的解决方案。
技术实现思路
针对相关技术中的上述技术问题，本专利技术提出一种基于混合亚像元解混的遥感多 光谱与全色图像融合方法，该方法对MS图像中的混合亚像元（即对应PAN纯像元的混合MS 像元的亚像元)，在融合过程中进行光谱解混处理，从而降低融合图像中混合像元数量，锐 化融合图像边界和降低光谱失真。 本专利技术的技术方案是这样实现的： 一种，包括以下步骤： 将原始低分辨率多光谱图像（#&)采用cubic方式上采样到原始全色图像韵分辨 率，得到并基于该#5计算归一化植被指数MWJ; 利用/"崩和柳KJ，对#5中的亚像元进行分类，得到植被-非植被混合亚像元和其他 亚像元 利用产崩和柳KJ，判断植被-非植被混合亚像元的地物类型类别，其中，所述地物类型 类别包括植被类、或非植被类； 对于#5中每个亚像元（，根据其对应的亚像元类别和地物类型类别进行融合。 此外，上述的还包括：将 采用平均法下采样到R倍空间分辨率，再采用立方卷积cubic方式上采用到原始分辨 率，得到产骑;其中，R为与烈崩空间分辨比率。 另外，在利用/"崩和柳KJ，对#5中的亚像元进行分类，得到植被-非植被混合亚像 元0卩其他亚像元印寸，可利用边缘检测Canny算子分别计算和AZU的边界石胃 万胃/;并利用直径为的圆盘结构元素，对万胃/进行形态学膨胀得到宽的边界图像 ，胃/ ;同时将在各《中值为1但在中值为〇的像元设置为〇,得到PAN边界图像 再采用直径为的圆盘结构元素对进行形态学膨胀得到宽的/边界图像，利 用^^^和对#5中的亚像元进行分类，得到植被-非植被混合亚像元0卩其他亚像元 C。。 此外，在利用，和，对MS中的亚像元进行分类，得到植被-非植被混合亚 像元^卩其他亚像元印寸，可利用中的亚像元分为两类，得到亚像元分 类图，并根据亚像元分类图，得到植被-非植被混合亚像元^卩其他亚像元Q其中，亚像 元分类图戽:的分类规则如下：其中，如果亚像元（在#^中对应值为1，则被分为植被-非植被混合亚像元'如果t#中对应值为0,则被分为其它亚像元类G。 另外，在利用烈辦卩M^KJ，判断植被-非植被混合亚像元的地物类型类别时，可根 据PAN边界图像iT_，计算植被边界像元图像&和非植被像元图像并利用直径为 3的圆盘结构元素，对&进行形态学膨胀，得到膨胀后的植被边界像元图像馬，并将馬中 值为1且在^r中值为1的像元设置为〇,并重复该步骤R-1次，使得馬中值为1的像元 均为植被像元；同时利用直径为3的圆盘结构元素，对.^进行形态学膨胀，得到JV，并将iV中值为1且在中值为1的像元设置为〇,并重复该步骤R-1次，使得1中值为1的 像元均为非植被像元；而对于1和&中每个值为1的像元（，分别在以（为中心、大小 为的领域窗口W((，匁内计算所有值为1的像元的值，并根据该值，采用 以下公式，分别得到图像垃和f;并对于mpi中每个属于混合亚像元类的像元（，根据和rf，采用以下公 式判断其属于植被G或非植被Cl类别，得到混合亚像元分类图 此外，在根据PAN边界图像iT_计算植被边界像元图像^和非植被像元图像 4时，可对于每个亚像元（，在其8邻域内查找与其在PAN图像中灰度差异最大的邻居 像元广得到像元对f};其中，如果则」Cr⑵=1，且( f) =0 ;如果⑵<M?KJ(f)，则拉观⑵=〇, 且C( (，）=1。 另外，在对于#5中每个亚像元（，根据其对应的亚像元类别和地物类型类别进行 融合时，可采用暗像元法计算#5和图像的雾气值贫和尽，其中，i等于1，……N，N 为#5图像的波段数量；并判断亚像元（是否属于其它亚像元类&并在判断结果为是的情 况下，利用该像元的原始光谱，采用以下公式进行融合，其中，藤天为#5的第i波段，MSf为融合图像MS/的第i波段；同时判断亚像元缇否属于 植被-非植被混合亚像元心并在判断结果为是的情况下，利用其邻域内的纯植被或纯非植 被像元的光谱，采用以下公式进行融合：其中，邻域内的纯植被或纯非植本文档来自技高网...
【技术保护点】
一种基于混合亚像元解混的遥感多光谱与全色图像融合方法，其特征在于，包括以下步骤：将原始低分辨率多光谱图像MS0采用立方卷积cubic方式上采样到原始全色图像PAN的分辨率，得到MS；并基于该MS计算归一化植被指数NDVI；利用PAN和NDVI，对MS中的亚像元进行分类，得到植被‑非植被混合亚像元CM和其他亚像元CO；利用PAN和NDVI，判断植被‑非植被混合亚像元的地物类型类别，其中，所述地物类型类别包括植被类、或非植被类；对于MS中每个亚像元t，根据其对应的亚像元类别和地物类型类别进行融合。
【技术特征摘要】
【专利技术属性】
技术研发人员：李慧，荆林海，王黎明，
申请(专利权)人：中国科学院遥感与数字地球研究所，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。