融合SAR图像信息的光学图像去云方法技术

本发明专利技术提供了一种融合SAR图像信息的光学图像去云方法，包括以下步骤：通过SAR图像及光学图像计算优先级函数，并根据优先级按先后顺序选取待修复样本块；对于选取的待修复样本块，利用对应的SAR图像中的样本块和光学图像样本块中的已知像素点从光学图像已知区域中选取最相似的匹配块；找到最相似的匹配块后，用其中对应的像素点对样本块中缺失像素点进行替换，得到修复后的样本块。本发明专利技术从优先级计算及最佳匹配块的搜索方面对基于块相似性的算法进行改进。联合SAR影像辅助光学影像进行云下缺失信息的修复，并引入更多相似性度量准则，从而得到更准确而不仅仅是看起来合理的修复结果。

【技术实现步骤摘要】
融合SAR图像信息的光学图像去云方法
本专利技术涉及遥感图像处理
，尤其涉及一种信息融合与光学图像修复的方法。
技术介绍
近年来，随着现代遥感技术的飞速发展，人们可获取到的遥感影像数据量迅速增长。然而光学图像中云层的干扰会导致图像信息的丢失，影响其后续解译及应用。目前根据辅助数据源的不同，常见的图像修复算法主要分为三类：基于多光谱的方法、基于空间的方法和基于多时相的方法。基于多光谱的方法利用多光谱含云波段及部分云穿透能力较强的冗余波段间的相关性，建立相应的经验或物理公式，来进行云层检测和信息重建；基于多时相的方法通过对于同一场景不同时刻的遥感图像，通过无云区域间的组合来实现云层修复，依据云层多存在于低频的特性，对多幅配准图像进行多尺度变换，在低频基于归一化低频指数设置阈值进行云层抑制，高频提取拥有更高能量的地物信息，最后进行反变换即可重构出无云影像；基于空间的方法假设缺失数据和剩余数据共享相同的统计或几何结构，于是使用数据中的已知部分填充缺失的数据区域来进行云层去除。基于块相似性的算法是一种经典的基于空间的图像修复算法，其结合了纹理合成和扩散填充的优点，可对较大面积的缺失区域进行修复并得到较为自然的结果。其首先在已知区域与缺失区域边界处通过优先级选取像素点，并以此像素点为中心取待修复样本块；然后借助待修复样本块中已知像素点来寻找源区域中匹配的目标块；最后将样本块中缺失像素点用匹配目标块中像素点进行替换，直到整个区域被修复。合成孔径雷达(SyntheticApertureRadar，SAR)作为一种主动式微波遥感手段，具有全天时、全天候的观测能力，不受云层、雨雪、光照条件等环境因素影响，于是可借助SAR图像来辅助光学图像进行修复。如现有技术中在传统基于块相似性算法的基础上提出了一种以SAR图像为目标区域内真实结构的光学图像修复方法，将图像分为海陆两部分，利用SAR图像来进行优先级的计算并在对应的光学图像区域选取待修复样本块，然后通过平方差之和在光学图像中寻找对应匹配块进行云下信息的修复。基于空间的方法对于填补小规模的数据缺失效果较好，但对于较大尺寸的信息缺失，由于缺乏辅助信息，只能得到在视觉上尽可能合理而不是准确的结果。第一，在优先级计算阶段，基于块相似性的算法中优先权函数由鼓励纹理传播的置信度项和鼓励几何结构的传播的数据项相乘得到。但其同时传播几何和纹理的特点有时会导致图像几何结构的破坏，且置信度项在修复过程中迅速下降到零的缺点也会导致后期修复顺序的错乱。同时，如前一步出现错误的填充，仅利用待修复光学图像来进行优先权函数计算的策略也会导致误差的积累，从而使得修复结果不准确。现有技术的算法虽然利用了SAR图像信息，但由于SAR与光学图像辐射特性存在较大差异，仅利用SAR图像信息的方法对于较复杂的图像修复效果不好。第二，在最优匹配块的搜索阶段，基于块相似性的算法选取SSD最小作为匹配标准，利用样本块中的已知像素点从源区域选取最相似的匹配块。但此标准较为单一，仅对光谱相似度进行了度量，在要修复区域纹理较为简单时效果较好，但在图像较为复杂时匹配错误的概率明显增大。且仅利用部分已知像素点进行匹配的方法也大大增加了出错率，经过不断的传播累积后仅能得到看似合理而不是准确的修复结果。
技术实现思路
有鉴于此，本专利技术的主要目的在于提供一种融合SAR图像信息的光学图像去云方法，以期部分地解决上述技术问题中的至少之一。为了实现上述目的，作为本专利技术的一方面，提供了一种融合SAR图像信息的光学图像去云方法，包括以下步骤：通过SAR图像及光学图像计算优先级函数，并根据优先级按先后顺序选取待修复样本块；对于选取的待修复样本块，利用对应的SAR图像中的样本块和光学图像样本块中的已知像素点从光学图像已知区域中选取最相似的匹配块；找到最相似的匹配块后，用其中对应的像素点对样本块中缺失像素点进行替换，得到修复后的样本块。其中，所述通过SAR图像及光学图像计算优先级函数，并根据优先级按先后顺序选取待修复样本块的步骤包括：由数据项将图像几何结构传播到目标区域；由混合项合成纹理。其中，所述优先级函数的计算公式如下：其中P(p)为优先级函数，D(p)为数据项，Pcd(p)为混合项，step为当前运行步数，为步数阈值。其中，在当前运行步数step小于等于时，优先级由数据项D(p)确定。其中，所述数据项D(p)为改进的数据项，由SAR图像数据项与光学图像数据项加权求和得到，所述数据项为边界在p处的单位法向量与等照线强度的乘积。其中，在当前运行步数step大于时，优先级由混合项Pcd(p)确定。其中，所述混合项由改进的数据项D(p)与正则化后的置信度项加权求和得到。其中，所述对于选取的待修复样本块，利用对应的SAR图像中的样本块和光学图像样本块中的已知像素点从光学图像已知区域中选取最相似的匹配块的步骤，包括：SAR图像中最佳匹配块通过平方差之和和结构相似性搜索，光学图像中最佳匹配块通过平方差之和搜索；对上述两部分进行加权求和；加入距离因子来使距离缺失块更近的样本块有更高的权重，最终综合得到最相似的匹配块。其中，先通过联合平方差之和筛选出前n个匹配块，再选出最优匹配快，以此来优化修复时间。其中，所述光学图像去云方法还包括重复全部步骤直到所有缺失区域修复完成。基于上述技术方案可知，本专利技术的光学图像去云方法相对于现有技术至少具有如下有益效果的一部分：本专利技术从优先级计算及最佳匹配块的搜索方面对基于块相似性的算法进行改进。联合SAR影像辅助光学影像进行云下缺失信息的修复，并引入更多相似性度量准则，从而得到更准确而不仅仅是看起来合理的修复结果。1.对于当前基于块相似性的修复算法仅利用光学图像来计算优先权函数从而导致后期修复顺序的错乱的问题，选择加入具有全天时、全天候，不受云层、雨雪等环境因素影响的SAR图像来作为辅助，其在光学图像被云层遮挡的部分仍可以正常成像，所以在计算数据项时加入部分SAR图像信息可以避免前期填充误差导致的后期修复顺序的错乱。2.基于块相似性的修复算法由鼓励纹理传播的置信度项和鼓励几何结构的传播的数据项相乘得到，但是同时传播几何和纹理的特点有时会导致图像几何结构的破坏，且置信度项在修复过程中迅速下降到零的缺点也会导致后期修复顺序的错乱。本专利技术选择分离的优先级计算策略，先传播几何结构，再修复纹理部分。在传播几何结构时，因为SAR图像的结构信息更加完整和准确，所以在数据项中使得SAR图像占有更高的权重；在传播纹理时，因为需要使得修复后图像的纹理信息连续性更高，可视性更好，所以在数据项加入更多的光学图像信息。对于置信度项在修复过程中会迅速下降到零的问题，对其进行正则化处理，并将混合项由相乘改为相加，从而使得置信度项快速下降的影响变小。3.在最优匹配块的搜索阶段，对于仅利用光学图像样本块中的部分已知像素点来进行匹配块选择，导致错误率较本文档来自技高网...

【技术保护点】
1.一种融合SAR图像信息的光学图像去云方法，其特征在于，包括以下步骤：/n通过SAR图像及光学图像计算优先级函数，并根据优先级按先后顺序选取待修复样本块；/n对于选取的待修复样本块，利用对应的SAR图像中的样本块和光学图像样本块中的已知像素点从光学图像已知区域中选取最相似的匹配块；/n找到最相似的匹配块后，用其中对应的像素点对样本块中缺失像素点进行替换，得到修复后的样本块。/n

【技术特征摘要】
1.一种融合SAR图像信息的光学图像去云方法，其特征在于，包括以下步骤：
通过SAR图像及光学图像计算优先级函数，并根据优先级按先后顺序选取待修复样本块；
对于选取的待修复样本块，利用对应的SAR图像中的样本块和光学图像样本块中的已知像素点从光学图像已知区域中选取最相似的匹配块；
找到最相似的匹配块后，用其中对应的像素点对样本块中缺失像素点进行替换，得到修复后的样本块。


2.根据权利要求1所述的光学图像去云方法，其特征在于，所述通过SAR图像及光学图像计算优先级函数，并根据优先级按先后顺序选取待修复样本块的步骤包括：
由数据项将图像几何结构传播到目标区域；
由混合项合成纹理。


3.根据权利要求1所述的光学图像去云方法，其特征在于，所述优先级函数的计算公式如下：



其中P(p)为优先级函数，D(p)为数据项，Pcd(p)为混合项，step为当前运行步数，为步数阈值。


4.根据权利要求3所述的光学图像去云方法，其特征在于，在当前运行步数step小于等于TΩs时，优先级由数据项D(p)确定。


5.根据权利要求4所述的光学图像去云方法，其特征在于，所述数据项D(p)为改进的数据...

【专利技术属性】
技术研发人员：陈子涵，许宁，尤红建，
申请(专利权)人：中国科学院空天信息创新研究院，
类型：发明
国别省市：北京;11