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

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

【技术实现步骤摘要】
一种卫星遥感图像区域条带噪声抑制装置及其方法
本专利技术属于卫星应用领域，涉及一种卫星遥感图像区域条带噪声抑制装置及其方法。该装置针对辐射校正后仍存在区域性条带噪声的遥感图像，通过区域条带噪声邻域信息恢复条带噪声区域的灰度信息，抑制区域条带噪声，同时又对抑噪结果进行灰度信息重建，保持图像无噪声区域的灰度信息不变。该方法具有保持图像辐射信息，去除区域条带噪声效果好的特点。
技术介绍
条带噪声是由于推扫成像传感器各CCD在光谱响应区内的响应函数不一致造成的，严重降低了遥感图像的辐射质量和清晰度，影响图像中地貌信息的解译和后续应用。目前遥感图像条带噪声消除方法主要是基于星上定标灯的内定标法、直方图统计法及基于图像分析的矩匹配法、直方图匹配法等。星上内定标法和直方图统计法是生成遥感图像辐射校正产品的常用方法。但是星上内定标方法由于星上环境复杂，不能完全去除条带，会有残余现象，而且许多星上并无定标系统。直方图统计法需要大量统计正常图像样本，不同地区地理位置及地物不同需要统计不同的样本，工作量很大，当出现新情况时会校正失败。直方图匹配和矩匹配方法一般是在遥感图像辐射校正之后，解决残余的条带噪声。直方图匹配法和矩匹配法的前提是假设各传感器所探测的地物具有相同均衡的辐射分布。将各传感器图像的直方图调整到参考直方图中，或将所有CCD探元形成的列图像的均值和方差调整到一致。该类方法不适用于复杂地物，同时由于在辐射校正后的一级遥感图像产品上根据图像直方图的概率分布或参考列图像均值和方差进行校正，会改变图像灰度的物理性质，对后续定量应用产生不确定的影响。目前对辐射校正后的遥感图像的条带噪声去除方法均为全局性的，会改变整幅遥感图像灰度信息，且主要处理单列的条带噪声，对于多列的宽区域条带噪声如CCD搭接区整体与图像辐射不一致的情况没有考虑。
技术实现思路
本专利技术的技术解决问题：克服现有技术的不足，提供一种卫星遥感图像区域条带噪声抑制装置及其方法，该装置首先对待处理图像大范围低通滤波，将区域条带噪声分离出来，根据条带噪声区域与邻域灰度信息的连续性，恢复噪声区域的灰度信息，抑制区域条带噪声，同时对抑噪结果进行灰度重建，保持待处理图像无噪声区域的灰度信息不变。使用该方法抑制区域条带噪声能够在保持图像辐射信息不变的情况下较好的去除区域条带噪声。本专利技术的技术解决方案：一种卫星遥感图像区域条带噪声抑制装置，包括参数设置模块，灰度复原模块和图像信息重建模块；参数设置模块读入待处理图像Iori，并设置流程迭代次数值it、第一平滑滤波器h1、第二平滑滤波器h2、第三平滑滤波器h3,和第一灰度复原系数阈值kth1、第二灰度复原系数阈值kth2,将待处理图像Iori、第一平滑滤波器h1和第一灰度复原系数阈值kth1、第二灰度复原系数阈值kth2送到灰度复原模块；将待处理图像Iori和第二平滑滤波器h2、第三平滑滤波器h3送到图像信息重建模块；灰度复原模块包含滤波单元，复原系数计算单元和噪声抑制单元；滤波单元根据参数设置模块提供的第一平滑滤波器h1对待处理图像Iori滤波处理，得到滤波后图像Ilp，待处理图像Iori与滤波后图像Ilp相减得到差值图像Idif，将滤波后图像Ilp和差值图像Idif输出给复原系数计算单元；复原系数计算单元根据差值图像Idif每一个像素灰度值与第一灰度复原系数阈值kth1之间的大小关系确定该像素对应位置第一复原系数，得到复原系数矩阵k1；复原系数计算单元计算待处理图像Iori和滤波后图像Ilp列均值之差，得到数组Colmeandif，计算差值图像Idif每一像素灰度值与该像素所在列在数组Colmeandif中对应值之差，比较该差与第二灰度复原系数阈值kth2之间的大小关系确定该像素对应位置第二复原系数，得到复原系数矩阵k2，将滤波后图像Ilp、差值图像Idif和灰度复原系数矩阵k1、k2输出给噪声抑制单元；在噪声抑制单元，差值图像Idif中每一像素灰度值与灰度复原系数矩阵k1、k2中对应位置值相乘，再与滤波后图像Ilp中对应位置像素灰度值相加，得到灰度复原图像Idns，将灰度复原图像Idns输出给图像信息重建模块；图像信息重建模块包括参考图像生成单元和重建结果输出单元；参考图像生成单元对灰度复原模块提供的灰度复原图像Idns和参数设置模块提供的待处理图像Iori分别采用参数设置模块提供的第二平滑滤波器h2、第三平滑滤波器h3进行滤波处理，生成第一参考图像第二参考图像输出给重建结果输出单元；重建结果输出单元比较参数设置模块提供的待处理图像Iori和参考图像生成单元提供的第一参考图像第二参考图像对应位置像素灰度值之间的大小关系，生成图像重建结果图像Iresult；在重建结果输出单元中，当迭代次数大于或者等于参数设置模块提供的迭代次数值it时，输出图像重建结果图像Iresult为区域条带噪声抑制结果图像，当迭代次数小于参数设置模块提供的迭代次数值it时，令图像重建结果图像Iresult为下一次迭代计算的待处理图像，输出给灰度复原模块继续计算。所述的参数设置模块包括图像获取单元和参数准备单元，图像获取单元获取待处理图像Iori并输出给参数准备单元；参数准备单元根据图像获取单元提供的待处理图像Iori的大小及噪声强度设置流程迭代次数值it、第一平滑滤波器h1、第二平滑滤波器h2、第三平滑滤波器h3,和第一灰度复原系数阈值kth1、第二灰度复原系数阈值kth2。一种卫星遥感图像区域条带噪声抑制方法，由计算准备阶段、灰度复原阶段、图像信息重建阶段和结果输出阶段配合完成，所述计算准备阶段实现步骤如下：步骤(1)，读取待处理图像Iori和Iori尺寸大小M×N；步骤(2)，根据步骤(1)读取的Iori大小M×N及区域噪声强度，设定处理流程迭代次数it，及每次迭代的第一灰度复原系数阈值kth1、第二灰度复原系数阈值kth2和第一平滑滤波器h1、第二平滑滤波器h2、第三平滑滤波器h3。所述灰度复原阶段实现步骤如下：步骤(3)，使用步骤(2)设定的第一平滑滤波器h1对步骤(1)得到的待处理图像Iori低通滤波，得到滤波后图像Ilp，计算待处理图像Iori与滤波后图像Ilp的差值图像Idif＝Iori-Ilp；步骤(4)，计算步骤(1)得到的待处理图像Iori和步骤(3)得到的滤波后图像Ilp每一列的均值数组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(j)，若abs(Colmeandif(j)-Idif(i,j))>kth2，则该像素对应位置第二复原系数k2(i,j)＝1，否则k2(i,j)＝abs(Colmeandif(j)-Idif(i,j))/kth2；步骤本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种卫星遥感图像区域条带噪声抑制装置，其特征在于：包括参数设置模块，灰度复原模块和图像信息重建模块；参数设置模块读入待处理图像Iori，并设置流程迭代次数值it、第一平滑滤波器h1、第二平滑滤波器h2、第三平滑滤波器h3,和第一灰度复原系数阈值kth1、第二灰度复原系数阈值kth2,将待处理图像Iori、第一平滑滤波器h1和第一灰度复原系数阈值kth1、第二灰度复原系数阈值kth2送到灰度复原模块；将待处理图像Iori和第二平滑滤波器h2、第三平滑滤波器h3送到图像信息重建模块；灰度复原模块包含滤波单元，复原系数计算单元和噪声抑制单元；滤波单元根据参数设置模块提供的第一平滑滤波器h1对待处理图像Iori滤波处理，得到滤波后图像Ilp，待处理图像Iori与滤波后图像Ilp相减得到差值图像Idif，将滤波后图像Ilp和差值图像Idif输出给复原系数计算单元；复原系数计算单元根据差值图像Idif每一个像素灰度值与第一灰度复原系数阈值kth1之间的大小关系确定该像素对应位置第一复原系数，得到复原系数矩阵k1；复原系数计算单元计算待处理图像Iori和滤波后图像Ilp列均值之差，得到数组Colmeandif，计算差值图像Idif每一像素灰度值与该像素所在列在数组Colmeandif中对应值之差，比较该差与第二灰度复原系数阈值kth2之间的大小关系确定该像素对应位置第二复原系数，得到复原系数矩阵k2，将滤波后图像Ilp、差值图像Idif和灰度复原系数矩阵k1、k2输出给噪声抑制单元；在噪声抑制单元，差值图像Idif中每一像素灰度值与灰度复原系数矩阵k1、k2中对应位置值相乘，再与滤波后图像Ilp中对应位置像素灰度值相加，得到灰度复原图像Idns，将灰度复原图像Idns输出给图像信息重建模块；图像信息重建模块包括参考图像生成单元和重建结果输出单元；参考图像生成单元对灰度复原模块提供的灰度复原图像Idns和参数设置模块提供的待处理图像Iori分别采用参数设置模块提供的第二平滑滤波器h2、第三平滑滤波器h3进行滤波处理，生成第一参考图像第二参考图像输出给重建结果输出单元；重建结果输出单元比较参数设置模块提供的待处理图像Iori和参考图像生成单元提供的第一参考图像第二参考图像对应位置像素灰度值之间的大小关系，生成图像重建结果图像Iresult；在重建结果输出单元中，当迭代次数大于或者等于参数设置模块提供的迭代次数值it时，输出图像重建结果图像Iresult为区域条带噪声抑制结果图像，当迭代次数小于参数设置模块提供的迭代次数值it时，令图像重建结果图像Iresult为下一次迭代计算的待处理图像，输出给灰度复原模块继续计算。2.根据权利要求1所述的一种卫星遥感图像区域条带噪声抑制装置，其特征在于：所述的参数设置模块包括图像获取单元和参数准备单元，图像获取单元获取待处理图像Iori并输出给参数准备单元；参数准备单元根据图像获取单元提供的待处理图像Iori的大小及噪声强度设置流程迭代次数值it、第一平滑滤波器h1、第二平滑滤波器h2、第三平滑滤波器h3,和第一灰度复原系数阈值kth1、第二...

【专利技术属性】
技术研发人员：岳春宇，何红艳，王治强，鲍云飞，邢坤，周楠，齐文雯，李岩，李方琦，
申请(专利权)人：北京空间机电研究所，
类型：发明
国别省市：