[发明专利]基于归一化的高分辨率遥感影像条带噪声去除方法

摘要

本发明公开了一种基于归一化的高分辨率遥感影像条带噪声去除方法，涉及遥感图像处理领域。该方法首先对影像均一地物作掩膜处理，并剔除影像0值和坏像元；然后根据高分卫星CCD推扫成像方式，通过选用矩匹配算法，将不同CCD器件的增益和偏置参数归一化为同一个值，从而有效地去除影像条带噪声；最后，通过动态范围匹配，将校正后的影像DN值范围进行拉伸至原始影像范围，不仅能够很好地去除国产高分影像的条带噪声，不会产生新的条带噪声，而且校正后的影像在色彩和对比度上与原始影像能够保持一致，不会破坏影像原始辐射信息。

主权项

1.一种基于归一化的高分辨率遥感影像条带噪声去除方法，其特征在于，包括如下步骤：S1，根据原始影像中没有条带噪声的均一地物，建立掩膜矩阵M1；根据原始影像中的0值和坏像元，建立掩膜矩阵M2；S2，将原始影像与掩膜矩阵M1、掩膜矩阵M2作掩膜运算，得到原始影像有效像元矩阵M3；S3，统计原始影像b波段有效像元DN值的动态范围[Tbmin，Tbmax]，其中Tbmin表示b波段最小的DN值，Tbmax表示b波段最大的DN值，b波段表示原始影像的任意波段；S4，统计原始影像b波段有效像元矩阵M3b的平均DN值μb和标准差σb，M3b为M3的子矩阵；S5，统计原始影像b波段有效像元矩阵M3b第c列的平均DN值μc和标准差σc，第c列表示M3b的任意列；S6，对原始影像b波段有效像元矩阵M3b中第c列有效像元逐一进行DN值归一化处理，完成对第c列有效像元的校正；采用如下公式进行DN值归一化处理：Yi＝(Ti‑μc)*σc/σb+μb，其中，Ti表示像元i的原始DN值，Yi表示像元i经过处理后的DN值，μb和σb分别表示原始影像b波段的平均DN值和标准差，μc和σc分别表示原始影像b波段第c列的平均DN值和标准差；S7，循环执行S6，直至完成对原始影像b波段有效像元矩阵M3b的全部有效像元的校正；S8，循环执行S7，直至完成对原始影像全部波段的有效像元的校正，得到校正后影像；S9，统计校正后影像b波段有效像元DN值的动态范围[Ybmin，Ybmax]，其中Ybmin表示校正后b波段最小的DN值，Ybmax表示校正后b波段最大的DN值；S10，将S9中校正后影像b波段DN值的动态范围与S3中原始影像b波段DN值的动态范围进行匹配，将校正后的影像DN值拉伸至原始影像范围，保持影像色彩和对比度一致；S11，循环执行S10，直至完成将校正后影像全部波段DN值拉伸至原始影像范围。