发明名称 扫描型红外成像系统的非均匀校正方法
摘要 本发明公开了一种扫描型红外成像系统的非均匀校正方法,包括以下步骤:从红外探测器获取红外视频流,并设置L=1;对红外视频流的第L帧图像进行垂直于扫描方向上的图像分割,以形成k个区域,并保证每个区域中沿扫描方向上的像素数大于或等于3;设置m=1;判断第L帧图像的第m个区域是否为第一帧图像的第1个区域;若第L帧图像的第m个区域是第一帧图像的第1个区域,则设置n=0,并读取线列探测器每个像元的增益系数初值G0(i)和偏移系数初值O0(i),i表示不同的像元;根据增益系数Gn(i)和偏移系数On(i)对第m个区域中每个像素进行非均匀校正,以得到非均匀校正结果;对第m个区域中每个像素的非均匀校正结果执行邻域中值处理,以得到非均匀校正期望结果。
申请公布号 CN102385701B 申请公布日期 2013.09.18
申请号 CN201110312154.0 申请日期 2011.10.14
申请人 华中科技大学 发明人 马泳;樊凡;周波;梁琨
分类号 G06K9/32(2006.01)I 主分类号 G06K9/32(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 朱仁玲
主权项 1.一种扫描型红外成像系统的非均匀校正方法,其特征在于,包括以下步骤:分别获取高温T<sub>H</sub>和低温T<sub>L</sub>下的均匀辐射场景图像的像素灰度值x<sub>i`j`</sub>(T<sub>H</sub>)和x<sub>i`j`</sub>(T<sub>L</sub>),其中i`为所述均匀辐射场景图像的像素的列坐标,j`为所述像素的行坐标;沿扫描方向逐行计算所述均匀辐射场景图像的灰度均值,以分别得线列扫描型探测器中第i个像元在高温T<sub>H</sub>下的平均灰度响应X<sub>i</sub>(T<sub>H</sub>)和低温T<sub>L</sub>下的平均灰度响应X<sub>i</sub>(T<sub>L</sub>),计算公式如下:<maths num="0001"><![CDATA[<math><mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><msup><mi>j</mi><mo>`</mo></msup></munder><msub><mi>x</mi><mrow><msup><mi>i</mi><mo>&prime;</mo></msup><msup><mi>j</mi><mo>,</mo></msup></mrow></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mrow></mrow><msup><mi>j</mi><mo>`</mo></msup></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>H</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><msup><mi>j</mi><mo>`</mo></msup></munder><msub><mi>x</mi><mrow><msup><mi>i</mi><mo>`</mo></msup><msup><mi>j</mi><mo>`</mo></msup></mrow></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>H</mi></msub><mo>)</mo></mrow></mrow><msup><mi>j</mi><mo>`</mo></msup></mfrac><mo>,</mo></mrow></math>]]></maths>其中i`=i;对所述高温T<sub>H</sub>下的平均灰度值X<sub>i</sub>(T<sub>H</sub>)和所述低温T<sub>L</sub>下的平均灰度值X<sub>i</sub>(T<sub>L</sub>)进行响应规格化处理,以得到增益系数初值G<sub>0</sub>(i)和偏移系数初值O<sub>0</sub>(i),具体公式为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>G</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>V</mi><msub><mi>T</mi><mi>H</mi></msub></msub><mo>-</mo><msub><mi>V</mi><msub><mi>T</mi><mi>L</mi></msub></msub></mrow><mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>H</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>O</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>V</mi><msub><mi>T</mi><mi>H</mi></msub></msub><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><msub><mi>T</mi><mi>L</mi></msub></msub><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>H</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>H</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中规格化的<img file="FDA00003062870200015.GIF" wi="61" he="73" />和<img file="FDA00003062870200016.GIF" wi="60" he="76" />为均匀辐射场景图像的焦平面线列响应的均值:<maths num="0005"><![CDATA[<math><mrow><msub><mi>V</mi><msub><mi>T</mi><mi>H</mi></msub></msub><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>H</mi></msub><mo>)</mo></mrow></mrow><mi>i</mi></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>V</mi><msub><mi>T</mi><mi>L</mi></msub></msub><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>L</mi></msub><mo>)</mo></mrow></mrow><mi>i</mi></mfrac><mo>;</mo></mrow></math>]]></maths>从红外探测器获取红外视频流,并设置L=1;对所述红外视频流的第L帧图像进行垂直于扫描方向上的图像分割,以形成k个区域,并保证每个所述区域中沿所述扫描方向上的像素数大于或等于3;设置m=1;判断所述第L帧图像的第m个区域是否为第一帧图像的第1个区域;若所述第L帧图像的第m个区域是所述第一帧图像的第1个区域,则设置n=0,并读取红外探测器每个像元的增益系数初值G<sub>0</sub>(i)和偏移系数初值O<sub>0</sub>(i),i表示不同的像元;根据所述增益系数G<sub>n</sub>(i)和所述偏移系数O<sub>n</sub>(i)对所述第m个区域中每个像素进行非均匀校正,以得到非均匀校正结果;对所述第m个区域中每个像素的非均匀校正结果执行邻域中值处理,以得到非均匀校正期望结果,是采用以下公式:f(h,j)=median{y(h+p,j+q)},其中f(h,j)为非均匀矫正期望结果,y(h+p,j+q)代表非均匀校正的结果,p∈{-1,0,1},q∈{-1,0,1},median表示对集合{y(h+p,j+q)}中的所有数取中值作为非均匀校正期望结果;根据所述第m个区域中每个像素的非均匀校正结果和非均匀校正期望结果,使用最陡下降法对所述增益系数G<sub>n</sub>(i)和所述偏移系数O<sub>n</sub>(i)进行处理,以得到G<sub>n+1</sub>(i)和O<sub>n+1</sub>(i),是采用以下公式:G<sub>n+1</sub>(i)=G<sub>n</sub>(i),<maths num="0007"><![CDATA[<math><mrow><msub><mi>O</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mn>2</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><msub><mi>O</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>&alpha;</mi><mrow><mo>(</mo><mi>y</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中h=i,α为迭代步长,N为所述第m个区域在扫描方向上的像素的数量;设置n=n+1;判断是否已经处理完所有k个区域;若已经处理完所有k个区域,则将k个区域的所有像素的非均匀校正结果进行合并,以得到第L个帧的红外校正图像;设置L=L+1,并重复所述对所述红外视频流的第L帧图像进行垂直于扫描方向上的图像分割,以形成k个区域,并保证每个所述区域中沿所述扫描方向上的像素数大于或等于3的步骤及其随后的所有步骤,直到红外视频流结束为止。
地址 430074 湖北省武汉市洪山区珞喻路1037号