一种图像条带噪声及坏线消除方法
【技术领域】
[0001] 本发明属于遥感图像处理领域，具体涉及一种图像条带噪声及坏线消除方法。
【背景技术】
[0002] 目前多数航空或航天光学遥感图像均以线阵推扫方式获取，当传感器线阵上某一 传感器出现故障时，就会造成采集的图像中有坏线出现。一部分坏线完全没有数值，无法使 用；另一部分则数值显著高于其它像素形成亮线，或数据显著低于其它像素形成暗线，这一 部分数据则可以通过处理使其成为正常像素继续使用。
[0003] 目前的主要的环线去除算法为均值替换法和矩滤波方法。均值替换法是利用与坏 线相临的两例数据求平均值代替坏线数据，这种方法完全舍弃了亮线或暗线中有用的数 据，造成了信息的损失。矩滤波方法则是利用相临图像的均值和方差为参考，使亮线或暗线 的均值与方差与之相等，但实际上坏线的均值和方差不可能与相临像素相同，因此处理后 可能会有明显的误差。
[0004]数据传感器的光电子元件，在接受的光信号之后，通过放大器放大，将光信号转化 成电信号，测量的地物辐射亮度值为：
[0005] 〇Ν=α*Ι+β (1)
[0006] 其中α代表放大器的增益，I代表传感器接收到的光信号强度，β代表为消除传感器 中的暗电流所加入的偏移量。α值随传感器的不同有所变化，β同α-样除随传感器发生变化 外，还受到传感器温度、时间的变化影响。线阵推扫成像方式的传感器，因此沿扫描带方向 的同列数据为同一传感器接收的数据，具有相同的α值和Μ直，因此坏线产生的原因主要由α 和邱勺误差造成。所以只要确定合适的α、β即可以实现对坏线的去除。一般情况下，传感器在 出厂以前α值已经经过准确的校正，造成的误差往往较小。β值则不停的变化，需要传感器在 成像前对暗电流进行测量确定，因此产生的误差可能性较大。因此合理调整m直可以实现对 坏线的消除。
【发明内容】
[0007] 本发明的目的是克服现有技术的缺陷，提供一种应用于线阵推扫方式成像的传感 器的图像条带噪声及坏线消除方法，通过对图像进行计算，消除影像中的非零值或饱和值 坏线，以及图像中的条带噪声。
[0008] 本发明是这样实现的：
[0009] -种图像条带噪声及坏线消除方法，具体包括如下步骤：
[0010] 步骤1:设图像为Μ行N列，计算图像的均值和各列的像素均值；
[0011 ]步骤2 :选取列像素均值与图像均值最接近的一列作为参考列，假设为第J列，令jL =jR = J〇
[0012] 步骤3:对第jL_l列像素，求数值仇，使得
[0013]
[0015] β表示像素值应加入的偏移量，见公式（l);p(i，j)为为图像中第i行的第j列像素 的值;将图像以第J列为中心，分为两部分，左侧的计算用九表示，右侧计算用j R表示；
[0016] 步骤4:计算p(jL-l，i)=p(jL-l，i)+βl，p(jR+l，i)=p(jR+l，i)+β2;i = l，2…M· [0017]步骤5:将步骤4得到的jL_l和jR+l作为新参考列，即分别令九=九-1，加=加+1。重 复步骤3和步骤4;若九=2,则结束九对应算式的计算;j_r = N-1贝lj结束j_r对应算式的计算；
[0018] 如上所述的步骤30中βι和02计算可利用如下步骤进行：
[0019] (1)图像中的相邻两列像素，对于i = l，2,…Μ，计算Δ pi = p( j，i)-p( j+1，i)，则Δ P={ ΔΡ1，Δρ2…ΔρΜ};其中M为图像的行数，p(i, j)为图像中第i行的第j列像素的值，APi 是同一行相邻的两个像素的差值。
[0020] (2)对ΔΡ中的元素按大小进行排序，开$成新序列ΔΡ={ Δρ^，Δρ、··· Δρ、}。
[0021 ] (3)贝旭=Δ pV〗，即Δ ρ向量各元素的中位数。
[0022] 本发明的有益效果是：
[0023] 本发明可以有效的消除影像中的垂直或水平条带噪声，经过处理后的图像不存在 明显的条带噪声现象。为了进一步验证方法的准确性，对质量较高的无条带噪声数据进行 验证计算，目的是为了解算法是否会对无条带数据进行错误的调整，对Aster数据进行的计 算表明，各列的m直大部分为〇，仅有极个别为1或2，说明本方法确定的β值精度较高。
【附图说明】
[0024] 图1是本发明的一种图像条带噪声及坏线消除方法的流程图。
【具体实施方式】
[0025]下面结合附图和实施例对本发明进行进一步描述。
[0026]如图1所示，一种图像条带噪声及坏线消除方法，具体包括如下步骤：
[0027]步骤1:设图像为Μ行Ν列，计算图像的均值和各列的像素均值；
[0028] 步骤2:选取列像素均值与图像均值最接近的一列作为参考列，假设为第J列，令九 =jR = J〇
[0029] 步骤3:对第jL_l列像素，求数值仇，使得
[0032] β表示像素值应加入的偏移量，见公式（1);1)(1，」）为为图像中第1行的第」列像素 的值;将图像以第J列为中心，分为两部分，左侧的计算用jL表示，右侧计算用j R表示；
[0033] 步骤 4:计算 p( jL_l，i)=p( jL_l，i)+0i，p( j_R+l，i)=p( jR+l，i)+&;i = l，2......M.
[0034] 步骤5:将步骤4得到的jL_l和jR+l作为新参考列，即分别令jL = jL-l，jR = jR+l<^ 复步骤3和步骤4;若j_l = 2,则结束j_l对应算式的计算;j_r = N-1贝lj结束j_r对应算式的计算； [0035] 其中步骤30中&和说计算可利用如下步骤进行：
[0036] (1)图像中的相邻两列像素，对于i = l，2,......M，计算Δ pi = p( j，i)-p( j+1，i)，则 AP={ ΔΡ1，Δρ2…ΔρΜ};其中M为图像的行数，p(i, j)为图像中第i行的第j列像素的值，Δ Pi是同一行相邻的两个像素的差值。
[0037] (2)对ΔΡ中的元素按大小进行排序，形成新序列ΔΡΗΔ?Λ,Δρ'^Δρ'Μ}。
[0038] (3)贝旭=Δρ'Μ/2,即ΔΡ向量各元素的中位数。
[0039]试验表明，该方法可以有效的消除影像中的垂直或水平条带噪声，经过处理后的 图像不存在明显的条带噪声现象。为了进一步验证方法的准确性，对质量较高的无条带噪 声数据进行验证计算，目的是为了解算法是否会对无条带数据进行错误的调整，对Aster数 据进行的计算表明，各列的m直大部分为〇，仅有极个别为1或2，说明本方法确定的β值精度 较高。
【主权项】
1. 一种图像条带噪声及坏线消除方法，具体包括如下步骤： 步骤1:设图像为Μ行N列，计算图像的均值和各列的像素均值； 步骤2:选取列像素均值与图像均值最接近的一列作为参考列，假设为第J列，令jL = jR =J〇 步骤3:对第jL-1列像素，求数值βι，使得为最小值;对第jR+1列像素，求数值防，使得为最小值。 β表示像素值应加入的偏移量，见公式（1);P(i，j)为为图像中第i行的第j列像素的值； 将图像W第J列为中屯、，分为两部分，左侧的计算用jL表示，右侧计算用jR表示； 步骤 4:计算 p( jL-l，i)=p(jL-l，i)+0i，p( jR+l，i)=p(jR+l，i)+02;i = l，2...M. 步骤5:将步骤4得到的jL-1和jR+1作为新参考列，即分别令北=北-1〇'1^叫1^+1。重复步 骤3和步骤4;若北=2,则结束北对应算式的计算;jR = N-l则结束jR对应算式的计算。2. 根据权利要求1所述的图像条带噪声及坏线消除方法，其特征在于:所述的步骤30中 &和防计算可利用如下步骤进行： (1) 图像中的相邻两列像素，对于i = l，2,…M，计算Δ pi = p( j，i)-p( j+1，i)，则Δ P = { Δρι，Δρ2. . . Αρμ};其中Μ为图像的行数，p(i, j)为图像中第i行的第j列像素的值，Api是 同一行相邻的两个像素的差值。 (2) 对ΔP中的元素按大小进行排序，形成新序列ΔP={Δp'l，Δp'2...ΔpΔM}。 (3) 贝化=Δρ'Μ/2,即ΔΡ向量各元素的中位数。
【专利摘要】本发明属于遥感图像处理领域，具体涉及一种图像条带噪声及坏线消除方法，的目的是克服现有技术的缺陷，提供一种应用于线阵推扫方式成像的传感器的图像条带噪声及坏线消除方法，通过对图像进行计算，消除影像中的非零值或饱和值坏线，以及图像中的条带噪声。本发明可以有效的消除影像中的垂直或水平条带噪声，经过处理后的图像不存在明显的条带噪声现象。为了进一步验证方法的准确性，对质量较高的无条带噪声数据进行验证计算，目的是为了解算法是否会对无条带数据进行错误的调整，对Aster数据进行的计算表明，各列的β值大部分为0，仅有极个别为1或2，说明本方法确定的β值精度较高。
【IPC分类】G06T5/00
【公开号】CN105551003
【申请号】CN201510961793
【发明人】陆冬华, 赵英俊, 张东辉, 周觅
【申请人】核工业北京地质研究院
【公开日】2016年5月4日
【申请日】2015年12月21日