一种卫星遥感图像区域条带噪声抑制装置及其方法.pdf

《一种卫星遥感图像区域条带噪声抑制装置及其方法.pdf》由会员分享，可在线阅读，更多相关《一种卫星遥感图像区域条带噪声抑制装置及其方法.pdf（10页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103761718 A (43)申请公布日 2014.04.30 CN 103761718 A (21)申请号 201410048541.1 (22)申请日 2014.02.12 G06T 5/00(2006.01) G06T 5/10(2006.01) (71)申请人 北京空间机电研究所 地址 100076 北京市丰台区南大红门路 1 号 9201 信箱 5 分箱 (72)发明人 岳春宇 何红艳 王治强 鲍云飞 邢坤 周楠 齐文雯 李岩 李方琦 (74)专利代理机构 中国航天科技专利中心 11009 代理人 安丽 (54) 发明名称 一种卫星遥感图像区域条带噪声抑制。

2、装置及 其方法 (57) 摘要 本发明属于卫星应用领域， 涉及一种卫星遥 感图像区域条带噪声抑制装置及其方法。该装置 包括参数设置模块， 灰度复原模块和图像信息重 建模块 ； 参数设置模块读入待处理图像， 并设置 处理流程的各项参数 ； 灰度复原模块包含滤波单 元， 复原系数计算单元和噪声抑制单元， 抑制待 处理图像中的区域条带噪声， 得到灰度复原图像 ； 图像信息重建模块包含参考图像生成单元和重建 结果输出单元， 在灰度信息复原结果图像基础上， 对待处理图像中无噪声区域灰度重建， 保持待处 理图像无噪声区域辐射信息不变， 生成区域条带 噪声抑制结果图像。该方法具有保持图像辐射信 息， 去除区。

3、域条带噪声效果好的特点。 (51)Int.Cl. 权利要求书 2 页 说明书 6 页 附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书6页 附图1页 (10)申请公布号 CN 103761718 A CN 103761718 A 1/2 页 2 1. 一种卫星遥感图像区域条带噪声抑制装置， 其特征在于 ： 包括参数设置模块， 灰度 复原模块和图像信息重建模块 ； 参数设置模块读入待处理图像 Iori， 并设置流程迭代次数值 it、 第一平滑滤波器 h1、 第 二平滑滤波器 h2、 第三平滑滤波器 h3, 和第一灰度复原系数阈值 kth1、 第二灰。

4、度复原系数阈 值 kth2, 将待处理图像 Iori、 第一平滑滤波器 h1和第一灰度复原系数阈值 kth1、 第二灰度复 原系数阈值 kth2送到灰度复原模块 ； 将待处理图像 Iori和第二平滑滤波器 h2、 第三平滑滤 波器 h3送到图像信息重建模块 ； 灰度复原模块包含滤波单元， 复原系数计算单元和噪声抑制单元 ； 滤波单元根据参数 设置模块提供的第一平滑滤波器 h1对待处理图像 Iori滤波处理， 得到滤波后图像 Ilp， 待处 理图像 Iori与滤波后图像 Ilp相减得到差值图像 Idif， 将滤波后图像 Ilp和差值图像 Idif输 出给复原系数计算单元 ； 复原系数计算单元根据。

5、差值图像 Idif每一个像素灰度值与第一灰 度复原系数阈值 kth1之间的大小关系确定该像素对应位置第一复原系数， 得到复原系数 矩阵 k1； 复原系数计算单元计算待处理图像 Iori和滤波后图像 Ilp列均值之差， 得到数组 Colmeandif， 计算差值图像 Idif每一像素灰度值与该像素所在列在数组 Colmeandif中对应值 之差， 比较该差与第二灰度复原系数阈值 kth2之间的大小关系确定该像素对应位置第二复 原系数， 得到复原系数矩阵 k2， 将滤波后图像 Ilp、 差值图像 Idif和灰度复原系数矩阵 k1、 k2 输出给噪声抑制单元 ； 在噪声抑制单元， 差值图像 Idif。

6、中每一像素灰度值与灰度复原系数 矩阵 k1、 k2中对应位置值相乘， 再与滤波后图像 Ilp中对应位置像素灰度值相加， 得到灰度 复原图像 Idns， 将灰度复原图像 Idns输出给图像信息重建模块 ； 图像信息重建模块包括参考图像生成单元和重建结果输出单元 ； 参考图像生成单元对 灰度复原模块提供的灰度复原图像 Idns和参数设置模块提供的待处理图像 Iori分别采用 参数设置模块提供的第二平滑滤波器 h2、 第三平滑滤波器 h3进行滤波处理， 生成第一参考 图像第二参考图像输出给重建结果输出单元 ； 重建结果输出单元比较参数设置 模块提供的待处理图像 Iori和参考图像生成单元提供的第一参。

7、考图像第二参考图像 对应位置像素灰度值之间的大小关系， 生成图像重建结果图像 Iresult； 在重建结果输出 单元中， 当迭代次数大于或者等于参数设置模块提供的迭代次数值 it 时， 输出图像重建结 果图像 Iresult为区域条带噪声抑制结果图像， 当迭代次数小于参数设置模块提供的迭代次 数值 it 时， 令图像重建结果图像 Iresult为下一次迭代计算的待处理图像， 输出给灰度复原 模块继续计算。 2. 根据权利要求 1 所述的一种卫星遥感图像区域条带噪声抑制装置， 其特征在于 ： 所 述的参数设置模块包括图像获取单元和参数准备单元， 图像获取单元获取待处理图像 Iori 并输出给参数。

8、准备单元 ； 参数准备单元根据图像获取单元提供的待处理图像 Iori的大小及 噪声强度设置流程迭代次数值 it、 第一平滑滤波器 h1、 第二平滑滤波器 h2、 第三平滑滤波器 h3, 和第一灰度复原系数阈值 kth1、 第二灰度复原系数阈值 kth2。 3. 一种卫星遥感图像区域条带噪声抑制方法， 其特征在于 ： 由计算准备阶段、 灰度复 原阶段、 图像信息重建阶段和结果输出阶段配合完成， 所述计算准备阶段实现步骤如下 ： 步骤 （1） ， 读取待处理图像 Iori和 Iori尺寸大小 MN ； 权 利 要 求 书 CN 103761718 A 2 2/2 页 3 步骤 （2） ， 根据步骤。

9、 （1） 读取的 Iori大小 MN 及区域噪声强度， 设定处理流程迭代次数 it， 及每次迭代的第一灰度复原系数阈值 kth1、 第二灰度复原系数阈值 kth2和第一平滑滤 波器 h1、 第二平滑滤波器 h2、 第三平滑滤波器 h3； 所述灰度复原阶段实现步骤如下 ： 步骤 （3） ， 使用步骤 （2）设定的第一平滑滤波器 h1对步骤 （1）得到的待处理图像 Iori低通滤波， 得到滤波后图像 Ilp， 计算待处理图像 Iori与滤波后图像 Ilp的差值图像 Idif=Iori-Ilp； 步骤 （4） ， 计算步骤 （1） 得到的待处理图像 Iori和步骤 （3） 得到的滤波后图像 Ilp每。

10、一 列的均值数组 Colmeanori和 Colmeanlp， Colmeanori和 Colmeanlp对应列相减得到其差数组 Colmeandif Colmeanori-Colmeanlp； 步骤 （5） ， 对步骤 （3）得到的差值图像 Idif的第 i 行 j 列个像素的灰度值 Idif(i,j) 进行判断， 若 abs(Idif(i,j)kth1， 则该像素对应位置第一复原系数 k1(i,j)=1， 否则 k1(i,j)=abs(Idif(i,j)/kth1； 步 骤 （6） ， 步 骤 （4）得 到 的 差 数 组 Colmeandif第 j 列 的 值 为 Colmeandif(。

11、j)， 若 abs(Colmean(j)-Idif(i,j)kth2， 则该像素对应位置第二复原系数 k2(i,j)=1， 否则 k2(i,j )=abs(Colmean(j)-Idif(i,j)/kth2； 步骤 （7） ， 根据步骤 （3） 得到的滤波后图像 Ilp、 差值图像 Idif和步骤 （5） 、 步骤 （6） 分别 得到的复原系数矩阵 k1、 k2， 计算灰度复原图像 Idns=Ilp+k1k2Idif； 所述图像信息重建阶段实现步骤如下 ： 步骤 （8） ， 使用步骤 （2） 设定的第二平滑滤波器 h2对步骤 （1） 得到的待处理图像 Iori和 步骤 （7） 得到的灰度复原图。

12、像 Idns低通滤波， 得到对待处理图像 Iori和灰度复原图像 Idns的 第二平滑滤波器 h2滤波后图像和计算第一参考图像 步骤 （9） ， 使用滤波器 h3对待处理图像 Iori和灰度复原图像 Idns低通滤波， 得到对待 处理图像 Iori和灰度复原图像 Idns的 h3滤波器滤波后图像和计算第二参考图像 步骤 （10） ， 根据步骤 （1） 得到的待处理图像 Iori、 步骤 （8） 得到的第一参考图像 步骤 （9）得到的第二参考图像对图像重建结果图像 Iresult的第 i 行 j 列个像素的 灰 度 值 Iresult(i,j) 进 行 判 断， 若则 否则 所述结果输出阶段实现。

13、步骤如下 ： 步骤 （11） ， 当迭代次数小于步骤 （2） 设定的流程迭代次数值 it 时， 令 Iori=Iresult， 转入 灰度复原阶段， 直到迭代结束， 输出结果 ； 当迭代次数大于等于步骤 （2） 设定的流程迭代次 数值 it 时， 输出结果 Iresult为区域条带抑制结果图像。 权 利 要 求 书 CN 103761718 A 3 1/6 页 4 一种卫星遥感图像区域条带噪声抑制装置及其方法 技术领域 0001 本发明属于卫星应用领域， 涉及一种卫星遥感图像区域条带噪声抑制装置及其方 法。该装置针对辐射校正后仍存在区域性条带噪声的遥感图像， 通过区域条带噪声邻域信 息恢复条带。

14、噪声区域的灰度信息， 抑制区域条带噪声， 同时又对抑噪结果进行灰度信息重 建， 保持图像无噪声区域的灰度信息不变。 该方法具有保持图像辐射信息， 去除区域条带噪 声效果好的特点。 背景技术 0002 条带噪声是由于推扫成像传感器各 CCD 在光谱响应区内的响应函数不一致造成 的， 严重降低了遥感图像的辐射质量和清晰度， 影响图像中地貌信息的解译和后续应用。 0003 目前遥感图像条带噪声消除方法主要是基于星上定标灯的内定标法、 直方图统计 法及基于图像分析的矩匹配法、 直方图匹配法等。星上内定标法和直方图统计法是生成遥 感图像辐射校正产品的常用方法。但是星上内定标方法由于星上环境复杂， 不能完。

15、全去除 条带， 会有残余现象， 而且许多星上并无定标系统。 直方图统计法需要大量统计正常图像样 本， 不同地区地理位置及地物不同需要统计不同的样本， 工作量很大， 当出现新情况时会校 正失败。 直方图匹配和矩匹配方法一般是在遥感图像辐射校正之后， 解决残余的条带噪声。 直方图匹配法和矩匹配法的前提是假设各传感器所探测的地物具有相同均衡的辐射分布。 将各传感器图像的直方图调整到参考直方图中， 或将所有 CCD 探元形成的列图像的均值和 方差调整到一致。该类方法不适用于复杂地物， 同时由于在辐射校正后的一级遥感图像产 品上根据图像直方图的概率分布或参考列图像均值和方差进行校正， 会改变图像灰度的物。

16、 理性质， 对后续定量应用产生不确定的影响。 0004 目前对辐射校正后的遥感图像的条带噪声去除方法均为全局性的， 会改变整幅遥 感图像灰度信息， 且主要处理单列的条带噪声， 对于多列的宽区域条带噪声如 CCD 搭接区 整体与图像辐射不一致的情况没有考虑。 发明内容 0005 本发明的技术解决问题 ： 克服现有技术的不足， 提供一种卫星遥感图像区域条带 噪声抑制装置及其方法， 该装置首先对待处理图像大范围低通滤波， 将区域条带噪声分离 出来， 根据条带噪声区域与邻域灰度信息的连续性， 恢复噪声区域的灰度信息， 抑制区域条 带噪声， 同时对抑噪结果进行灰度重建， 保持待处理图像无噪声区域的灰度信。

17、息不变。 使用 该方法抑制区域条带噪声能够在保持图像辐射信息不变的情况下较好的去除区域条带噪 声。 0006 本发明的技术解决方案 ： 一种卫星遥感图像区域条带噪声抑制装置， 包括参数设 置模块， 灰度复原模块和图像信息重建模块 ； 0007 参数设置模块读入待处理图像 Iori， 并设置流程迭代次数值 it、 第一平滑滤波器 h1、 第二平滑滤波器h2、 第三平滑滤波器h3,和第一灰度复原系数阈值kth1、 第二灰度复原系 说 明 书 CN 103761718 A 4 2/6 页 5 数阈值 kth2, 将待处理图像 Iori、 第一平滑滤波器 h1和第一灰度复原系数阈值 kth1、 第二灰。

18、 度复原系数阈值 kth2送到灰度复原模块 ； 将待处理图像 Iori和第二平滑滤波器 h2、 第三平 滑滤波器 h3送到图像信息重建模块 ； 0008 灰度复原模块包含滤波单元， 复原系数计算单元和噪声抑制单元 ； 滤波单元根据 参数设置模块提供的第一平滑滤波器 h1对待处理图像 Iori滤波处理， 得到滤波后图像 Ilp， 待处理图像 Iori与滤波后图像 Ilp相减得到差值图像 Idif， 将滤波后图像 Ilp和差值图像 Idif 输出给复原系数计算单元 ； 复原系数计算单元根据差值图像 Idif每一个像素灰度值与第一 灰度复原系数阈值 kth1之间的大小关系确定该像素对应位置第一复原系。

19、数， 得到复原系 数矩阵 k1； 复原系数计算单元计算待处理图像 Iori和滤波后图像 Ilp列均值之差， 得到数组 Colmeandif， 计算差值图像 Idif每一像素灰度值与该像素所在列在数组 Colmeandif中对应值 之差， 比较该差与第二灰度复原系数阈值 kth2之间的大小关系确定该像素对应位置第二复 原系数， 得到复原系数矩阵 k2， 将滤波后图像 Ilp、 差值图像 Idif和灰度复原系数矩阵 k1、 k2 输出给噪声抑制单元 ； 在噪声抑制单元， 差值图像 Idif中每一像素灰度值与灰度复原系数 矩阵 k1、 k2中对应位置值相乘， 再与滤波后图像 Ilp中对应位置像素灰度。

20、值相加， 得到灰度 复原图像 Idns， 将灰度复原图像 Idns输出给图像信息重建模块 ； 0009 图像信息重建模块包括参考图像生成单元和重建结果输出单元 ； 参考图像生成单 元对灰度复原模块提供的灰度复原图像 Idns和参数设置模块提供的待处理图像 Iori分别采 用参数设置模块提供的第二平滑滤波器 h2、 第三平滑滤波器 h3进行滤波处理， 生成第一参 考图像第二参考图像输出给重建结果输出单元 ； 重建结果输出单元比较参数设 置模块提供的待处理图像 Iori和参考图像生成单元提供的第一参考图像第二参考图 像对应位置像素灰度值之间的大小关系， 生成图像重建结果图像 Iresult； 在重。

21、建结果输 出单元中， 当迭代次数大于或者等于参数设置模块提供的迭代次数值 it 时， 输出图像重建 结果图像 Iresult为区域条带噪声抑制结果图像， 当迭代次数小于参数设置模块提供的迭代 次数值 it 时， 令图像重建结果图像 Iresult为下一次迭代计算的待处理图像， 输出给灰度复 原模块继续计算。 0010 所述的参数设置模块包括图像获取单元和参数准备单元， 图像获取单元获取待 处理图像 Iori并输出给参数准备单元 ； 参数准备单元根据图像获取单元提供的待处理图像 Iori的大小及噪声强度设置流程迭代次数值 it、 第一平滑滤波器 h1、 第二平滑滤波器 h2、 第 三平滑滤波器 。

22、h3, 和第一灰度复原系数阈值 kth1、 第二灰度复原系数阈值 kth2。 0011 一种卫星遥感图像区域条带噪声抑制方法， 由计算准备阶段、 灰度复原阶段、 图像 信息重建阶段和结果输出阶段配合完成， 所述计算准备阶段实现步骤如下 ： 0012 步骤 （1） ， 读取待处理图像 Iori和 Iori尺寸大小 MN ； 0013 步骤 （2） ， 根据步骤 （1） 读取的 Iori大小 MN 及区域噪声强度， 设定处理流程迭代 次数 it， 及每次迭代的第一灰度复原系数阈值 kth1、 第二灰度复原系数阈值 kth2和第一平 滑滤波器 h1、 第二平滑滤波器 h2、 第三平滑滤波器 h3。 。

23、0014 所述灰度复原阶段实现步骤如下 ： 0015 步骤 （3） ， 使用步骤 （2）设定的第一平滑滤波器 h1对步骤 （1）得到的待处理图 像 Iori低通滤波， 得到滤波后图像 Ilp， 计算待处理图像 Iori与滤波后图像 Ilp的差值图像 说 明 书 CN 103761718 A 5 3/6 页 6 Idif=Iori-Ilp； 0016 步骤 （4） ， 计算步骤 （1） 得到的待处理图像 Iori和步骤 （3） 得到的滤波后图像 Ilp每 一列的均值数组 Colmeanori和 Colmeanlp， Colmeanori和 Colmeanlp对应列相减得到其差数组 Colmean。

24、dif=Colmeanori-Colmeanlp； 0017 步骤 （5） ， 对步骤 （3） 得到的差值图像 Idif的第 i 行 j 列个像素的灰度值 Idif(i,j) 进行判断， 若 abs(Idif(i,j)kth1， 则该像素对应位置第一复原系数 k1(i,j)=1， 否则 k1(i,j)=abs(Idif(i,j)/kth1； 0018 步骤 （6） ， 步骤 （4）得到的差数组 Colmeandif第 j 列的值为 Colmeandif(j)， 若 abs(Colmean(j)-Idif(i,j)kth2， 则该像素对应位置第二复原系数 k2(i,j)=1， 否则 k2(i,j。

25、 )=abs(Colmean(j)-Idif(i,j)/kth2； 0019 步骤 （7） ， 根据步骤 （3） 得到的滤波后图像 Ilp、 差值图像 Idif和步骤 （5） 、 步骤 （6） 分别得到的复原系数矩阵 k1、 k2， 计算灰度复原图像 Idns=Ilp+k1k2Idif； 0020 所述图像信息重建阶段实现步骤如下 ： 0021 步骤 （8） ， 使用步骤 （2）设定的第二平滑滤波器 h2对步骤 （1）得到的待处 理图像 Iori和步骤 （7）得到的灰度复原图像 Idns低通滤波， 得到对待处理图像 Iori和 灰度复原图像 Idns的第二平滑滤波器 h2滤波后图像和计算第一参。

26、考图像 0022 步骤 （9） ， 使用滤波器 h3对待处理图像 Iori和灰度复原图像 Idns低通滤波， 得到对 待处理图像 Iori和灰度复原图像 Idns的 h3滤波器滤波后图像和计算第二参考图 像 0023 步骤 （10） ， 根据步骤 （1）得到的待处理图像 Iori、 步骤 （8）得到的第一参考图像 步骤 （9） 得到的第二参考图像对图像重建结果图像 Iresult的第 i 行 j 列个像素 的灰度值 Iresult(i,j) 进行判断， 若则 否则 0024 所述结果输出阶段实现步骤如下 ： 0025 步骤 （11） ， 当迭代次数小于步骤 （2） 设定的流程迭代次数值 it 。

27、时， 令 Iori=Iresult， 转入灰度复原阶段， 直到迭代结束， 输出结果 ； 当迭代次数大于等于步骤 （2） 设定的流程迭 代次数值 it 时， 输出结果 Iresult为区域条带抑制结果图像。 0026 本发明与现有技术相比具有如下优点 ： 0027 （1） 首先对图像大范围低通滤波， 将区域条带噪声分离出来， 根据条带噪声区域与 邻域灰度信息的连续性， 恢复噪声区域的灰度信息， 抑制区域条带噪声， 同时对抑噪结果进 行灰度重建， 保持待处理图像无噪声区域的灰度信息不变。 0028 （2） 使用该方法抑制区域条带噪声能够在保持图像辐射信息不变的情况下较好的 去除区域条带噪声。 说 。

28、明 书 CN 103761718 A 6 4/6 页 7 附图说明 0029 图 1 为本发明的装置结构示意图 ； 0030 图 2 为本发明的操作流程图 ； 具体实施方式 0031 本发明的基本思路 ： 针对辐射校正后仍存在区域性条带噪声的卫星遥感图像， 通 过区域条带噪声邻域信息恢复条带噪声区域的灰度信息， 抑制区域条带噪声， 同时又对抑 噪结果进行灰度重建， 保持待处理图像无噪声区域的灰度信息不变。该方法具有保持图像 辐射信息， 去除区域条带噪声效果好的特点。 0032 下面结合附图对本发明做进一步详细描述， 如图 1 所示一种卫星遥感图像区域条 带噪声抑制装置， 包括参数设置模块， 灰。

29、度复原模块和图像信息重建模块 ； 0033 参数设置模块读入待处理图像 Iori， 并设置流程迭代次数值 it、 第一平滑滤波器 h1、 第二平滑滤波器h2、 第三平滑滤波器h3,和第一灰度复原系数阈值kth1、 第二灰度复原系 数阈值 kth2, 将待处理图像 Iori、 第一平滑滤波器 h1和第一灰度复原系数阈值 kth1、 第二灰 度复原系数阈值 kth2送到灰度复原模块 ； 将待处理图像 Iori和第二平滑滤波器 h2、 第三平 滑滤波器 h3送到图像信息重建模块 ； 0034 灰度复原模块包含滤波单元， 复原系数计算单元和噪声抑制单元 ； 滤波单元根据 参数设置模块提供的第一平滑滤波。

30、器 h1对待处理图像 Iori滤波处理， 得到滤波后图像 Ilp， 待处理图像 Iori与滤波后图像 Ilp相减得到差值图像 Idif， 将滤波后图像 Ilp和差值图像 Idif 输出给复原系数计算单元 ； 复原系数计算单元根据差值图像 Idif每一个像素灰度值与第一 灰度复原系数阈值 kth1之间的大小关系确定该像素对应位置第一复原系数， 得到复原系 数矩阵 k1； 复原系数计算单元计算待处理图像 Iori和滤波后图像 Ilp列均值之差， 得到数组 Colmeandif， 计算差值图像 Idif每一像素灰度值与该像素所在列在数组 Colmeandif中对应值 之差， 比较该差与第二灰度复原系。

31、数阈值 kth2之间的大小关系确定该像素对应位置第二复 原系数， 得到复原系数矩阵 k2， 将滤波后图像 Ilp、 差值图像 Idif和灰度复原系数矩阵 k1、 k2 输出给噪声抑制单元 ； 在噪声抑制单元， 差值图像 Idif中每一像素灰度值与灰度复原系数 矩阵 k1、 k2中对应位置值相乘， 再与滤波后图像 Ilp中对应位置像素灰度值相加， 得到灰度 复原图像 Idns， 将灰度复原图像 Idns输出给图像信息重建模块 ； 0035 图像信息重建模块包括参考图像生成单元和重建结果输出单元 ； 参考图像生成单 元对灰度复原模块提供的灰度复原图像 Idns和参数设置模块提供的待处理图像 Ior。

32、i分别采 用参数设置模块提供的第二平滑滤波器 h2、 第三平滑滤波器 h3进行滤波处理， 生成第一参 考图像第二参考图像输出给重建结果输出单元 ； 重建结果输出单元比较参数设 置模块提供的待处理图像 Iori和参考图像生成单元提供的第一参考图像第二参考图 像对应位置像素灰度值之间的大小关系， 生成图像重建结果图像 Iresult； 在重建结果输 出单元中， 当迭代次数大于或者等于参数设置模块提供的迭代次数值 it 时， 输出图像重建 结果图像 Iresult为区域条带噪声抑制结果图像， 当迭代次数小于参数设置模块提供的迭代 次数值 it 时， 令图像重建结果图像 Iresult为下一次迭代计算。

33、的待处理图像， 输出给灰度复 原模块继续计算。 说 明 书 CN 103761718 A 7 5/6 页 8 0036 所述的参数设置模块包括图像获取单元和参数准备单元， 图像获取单元获取待 处理图像 Iori并输出给参数准备单元 ； 参数准备单元根据图像获取单元提供的待处理图像 Iori的大小及噪声强度设置流程迭代次数值 it、 第一平滑滤波器 h1、 第二平滑滤波器 h2、 第 三平滑滤波器 h3, 和第一灰度复原系数阈值 kth1、 第二灰度复原系数阈值 kth2。 0037 本发明优选的处理流程如图 2 所示， 一种卫星遥感图像区域条带噪声抑制方法的 实施可以分为四个阶段， 计算准备阶。

34、段、 灰度复原阶段、 图像信息重建阶段和结果输出阶 段。 0038 1. 计算准备阶段 0039 （1） 由参数设置模块图像获取单元读取待处理图像 Iori及 Iori尺寸大小 （2867 像 素 3036 像素） ； 0040 （2） 根据 Iori大小及区域噪声强度， 在参数设置模块参数准备单元设定迭代次数 it=2， 设定第一次迭代的两个灰度恢复阈值 kth1=80、 kth2=20， 第一次迭代的滤波器 h1为 1401 的均值滤波器、 h2为 4011 的均值滤波器、 h3为 1011 的均值滤波器 ； 设定第二次 迭代的两个灰度恢复阈值kth1=200、 kth2=20， 第二次迭。

35、代的滤波器h1为121的均值滤 波器、 h2 为 1011 的均值滤波器、 h3 为 4011 的均值滤波器。 0041 2. 灰度复原阶段 0042 （3） 由灰度复原模块调用滤波单元使用滤波器 h1对 Iori低通滤波， 得到 Ilp， 计算 Iori与 Ilp的差值图像 Idif=Iori-Ilp； 0043 （4）在灰度复原模块复原系数计算单元计算 Iori和 Ilp每一列的均值数组 Colmeanori和 Colmeanlp， Colmeanori和 Colmeanlp对应列相减得到其差数组 Colmeandif=Colm eanori-Colmeanlp； 0044 （5）在灰度复。

36、原模块复原系数计算单元中， 对于 Idif的第 i 行 j 列个像素的灰 度值 Idif(i,j)， 若 abs(Idif(i,j)kth1， 则该像素对应位置第一复原系数 k1(i,j)=1， 否则 k1(i,j)=abs(Idif(i,j)/kth1； 0045 （6） Colmeandif(j) 为 Colmeandif第 j 列的值， 若 abs(Colmean(j)-Idif(i,j)kth2， 则该像素对应位置第二复原系数 k2(i,j)=1， 否则 k2(i,j)=abs(Colmean(j)-Idif(i,j)/ kth2； 0046 （7） 根据滤波后图像 Ilp、 差值图像。

37、 Idif和步骤 （5） 、 步骤 （6） 分别得到的复原系数 矩阵 k1、 k2， 计算灰度复原图像 Idns=Ilp+k1k2Idif。 0047 3. 图像信息重建阶段 0048 （8） 在图像信息重建模块参考图像生成单元， 使用滤波器 h2对 Iori和 Idns低通滤 波， 得到和计算第一参考图像 0049 （9） 在图像信息重建模块调用滤波单元， 使用滤波器 h3对 Iori和 Idns低通滤波， 得 到对待处理图像 Iori和灰度复原图像 Idns的 h2滤波器滤波后图像和计算第二参 考图像 0050 （10） 在图像信息重建模块调用滤波单元， 对于图像重建结果图像 Iresul。

38、t的第 i 行 j 列个像素的灰度值 Iresult(i,j)， 若则 说 明 书 CN 103761718 A 8 6/6 页 9 否则 0051 4. 结果输出阶段 0052 （11） 迭代次数为 1， 小于 it， 令 Iori=Iresult， 转入灰度复原阶段， 直到步骤 （10） ； 0053 （12） 迭代次数为 2， 等于 it， 输出结果。 0054 本实施例中， 由本发明抑制区域条带噪声结果良好。辐射效果保持情况见下表 : 0055 像素数 （个）占全图比例 辐射信息改变3341373.84% 辐射信息不变837007596.16% 0056 本发明未详细阐述部分属于本领域公知技术。 说 明 书 CN 103761718 A 9 1/1 页 10 图 1 图 2 说 明 书 附 图 CN 103761718 A 10 。