发明名称 基于运动检测指导的红外焦平面非均匀性校正方法
摘要 基于运动检测指导的红外焦平面非均匀性校正方法,属于红外焦平面探测器领域,目的在于能够自适应,有效抑制目标退化,解决伪像(鬼影)问题,并且易于硬件实现,本发明包括(1)预处理步骤,(2)校正步骤,(3)迭代步长调整步骤,(4)增益校正系数更新步骤。本发明有效的解决和抑制了场景类算法中产生的目标退化和伪象现象,去除了固定图案噪声,具有自适应性,计算操作简单,利于硬件实现,达到实时性。
申请公布号 CN100433793C 申请公布日期 2008.11.12
申请号 CN200710051920.6 申请日期 2007.04.19
申请人 华中科技大学 发明人 张天序;桑红石;钟胜;李洁珺;袁雅婧;施长城;周泱;刘慧娜
分类号 H04N5/14(2006.01);H04N7/26(2006.01) 主分类号 H04N5/14(2006.01)
代理机构 华中科技大学专利中心 代理人 方放
主权项 1.一种基于运动检测指导的红外焦平面非均匀性校正方法,顺序包括:(1)预处理步骤,用红外焦平面探测器采集M帧均匀辐照的图像,然后对它们进行时域平均得到背景帧图像B,其各象元点灰度值Bi,j;对背景帧图像各象元点利用基于场景的无效象元检测算法进行无效象元检测,检测结果存入坏元模板,根据检测结果计算背景帧图像中有效象元的灰度平均值Bave;设定待校正原始图像序列总帧数为N,N为自然数,初始化待校正原始图像各个象元的增益校正系数<math><mrow><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>1</mn></msubsup><mo>=</mo><mn>1</mn><mo>,</mo></mrow>初始化待校正原始图像序列号n=1,M为10~100;(2)校正步骤,输入第n帧原始图像Vn,其各象元点灰度值为Vi,j n,将Vi,j n减去背景帧对应位置的灰度值Bi,j,得到第n帧待校正图像Xn,其各象元点灰度值为Xi,j n:<math><mrow><msubsup><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo></mrow>其中下标i,j分别为象元在图像中的横坐标和纵坐标;读取坏元模板,对待校正图像Xn对应坏元模板中有效象元位置的象元进行校正,校正公式为:<math><mrow><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>&times;</mo><msubsup><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msub><mi>B</mi><mi>ave</mi></msub><mo>,</mo></mrow>其中,Gi,j n为第n帧增益校正系数;同时对待校正图像Xn对应无效象元位置的像素值用左邻域象元像素和上领域象元像素之和的平均值替代,最后得到校正后图像Yn,其各象元灰度值为Yi,j n;(3)迭代步长调整步骤,判断是否n>1,否则转步骤(4);是则计算前后两帧校正后图像对应象元运动方差σT,ij 2,判断是否<math><mrow><msubsup><mi>&sigma;</mi><mrow><mi>T</mi><mo>,</mo><mi>ij</mi></mrow><mn>2</mn></msubsup><mo>></mo><msubsup><mi>&sigma;</mi><mrow><mi>T</mi><mi>max</mi></mrow><mn>2</mn></msubsup><mo>,</mo></mrow>否则调整迭代步长μi,j n,转步骤(4),是则再次初始化<math><mrow><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mn>1</mn><mo>,</mo></mrow>待校正原始图像序列号n加1,判断是否n>N,是则结束,否则转步骤(2),其中运动方差阈值σTmax 2为(0.5~1)Bave;(4)增益校正系数更新步骤,利用神经网络校正算法对图像的增益非均匀性进行补偿,即更新增益校正系数,公式为:<math><mrow><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><mn>2</mn><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><msubsup><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow><math><mrow><msubsup><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mo>/</mo><mn>4</mn><mo>,</mo></mrow>其中fi,j n为校正后图像Yn中各象元的四邻域象元灰度值均值,增益校正系数更新后,待校正原始图像序列号n加1,判断是否n>N,是则结束,否则转步骤(2)。
地址 430074湖北省武汉市洪山区珞喻路1037号
您可能感兴趣的专利