奇异值分解在红外弱小目标背景抑制中的应用

#奇异值分解在红外弱小目标背景抑制中的

应用2010010208张翠翠（控制科学与工程学院 控制理论与控制工程 2010010208）摘要:复杂背景的抑制是红外弱小目标检测技术的一个难题。为解决这个问题,本文提、引言红外告警系统因其被动探测、 高度隐蔽的特点而受到广泛的重视。 当红外弱小目标距离很远时,其成像面积非常小,且目标与背景的对比度、信噪比较低 ,常表现为淹没在复杂背景 （例如大面积的云层和地面建筑物 ）中的几个像点,即为弱小目标。如果要可靠、稳定、准确地检测并跟踪这类目标 ,则必须对图像进行预处理,而高性能的背景抑制是其中重要而关键的一项预处理。近20多年来,红外图像背景抑制技术得到较大的发展 ,主要有时域滤波、空域滤波、频域滤波、数学形态学滤波和小波域滤波等滤波方法。但是 ,当背景比较复杂时 ,这类滤波算法不能完全平滑边缘 ,从而导致检测概率的降低 ,虚警率增大。在这种情况下 ,为了使有用的目标特征被保留并得到有效增强 ,则必须要对复杂背景实行自适应的抑制。为此 ,本文提出了一种基于奇异值分解的背景抑制算法。作为一种非线性滤波 ,其从矩阵的角度出发 ,对图像矩阵进行奇异值分解 ,并根据定义的偏差指数确定有效奇异值来重构图像 ,从而达到平滑复杂背景 ,增强其突变部分 ,即目标信号的目的。 用真实的红外图像序列进行实验 ,结果验证了本算法能对复杂背景有效地平滑 ,增强其突变部分 ,即抑制了复杂背景 ,增强了目标信号。、基于奇异值分解的弱小目标背景抑制奇异值分解□□□□□□A是m口n的实矩阵，且A的秩为r,其中 r口min（m,n）,则存在m阶酉矩阵U和n阶酉矩阵V使得A的奇异值分解可用下式来表示 ：(1)a=urvh(1)mxn mxmmxnnxn式中r式中rmxn□是r□□□□□，口=diag（a,a，…,a）;U=[U,U,…，U],由m维列向量U=[u,u,…，u[T（i=1,2，…，m）构成;V=[V,V,…，V],由n维列向量V=[V,Vi,/，vk（i=1,2，…，n）构成;r□An谴零奇异值个数 n，其中入i（入1口入i2TOC\o"1-5"\h\z1i2i ni□□入r>0）是AtA并且也是 AAt的非零特征值的全体 ；Ui,Vi分别是 AAt和AtA□□□□□□□ 入i□□□□□□□□ ，式⑴写成如下形式：A=XoUVt （2）mxn iiii=1对式（2）可以理解为图像矩阵 A□□□□□□□□□□□□□□□□□□□□□息,可以想象那些接近零的奇异值也只含有少量矩阵重构信息 ,因此重构图像矩阵时可以将其忽略,只利用携带其信息的非零奇异值进行重构即可。U=£oUVt （3）mxn iiii=1式中，s（s口r）表示重构A时所需的奇异值数目。2.2基于奇异值分解的背景抑制由于红外探测器接受红外辐射只能反映其空间的明暗程度， 红外图像反映在□□□□□□□□□□□, □□□□□□□□ BMP图像□□□□□□□□□□□□□□BMP□□□□□□ ,□□□□□□□□□□□□□□□□□□□□□的情况（而不是彩色图像中的色彩） 。这正反应了红外辐射规律：温度高、辐射强的地方的亮度高， 在图像上像素的反映就更接近白色， 反之就越接近黑色。 BMP图像是用彩色成分 RGB来描述的， 即像素的色彩值包含的红、 绿、蓝三种颜色的多少，而灰度图像的明暗程度，是用灰色（就是红、绿、蓝三种颜色成分完全相等的颜色）来表示的，这一点在图像处理中是需要注意的。灰度也可认为是亮度， 简单的说就是色彩的深浅程度。 实际上在我们的日常生活中，任何颜色都有红、绿、蓝三原色组成，通过三原色色彩深浅的组合，可□□□□□□□□□□□□□□□□, R、G、B值均相等，称为灰度值，每一个像素都对应一个灰度值。对于 8位的灰度图像，其灰度值范围为 0—255。产品能够展现的灰度数量越多， 也就意味着这款产品的色彩表现力更加丰富， 能够实现更强的色彩层次。例如三原色16级灰度，能显示的颜色就是16016016=4096色。空中红外运动弱小目标图像由目标、 背景（主要指天空和云层 ）和红外成像器的电子噪声构成。其中 ,具有云层反射的复杂背景 ,背景在空间上有缓变部分 ,也有剧变部分 ,背景是在时间上体现为缓变 ;噪声按照产生的物理机理的不同 ,可分为热噪声、散粒噪声、1/f噪声三大类,前二者是白噪声 ,主要出现在中、 高频段,幅度与频率无关 ,而1/f噪声是色噪声 ,出现在低频段 ,幅度与频率成反比。 因此,图像中的目标即使在整幅图像中强度不是最强的 ,但在它所处的小区域中与局部背景的差异较明显 ,一般要高于局部背景的辐射强度 ,另一方面 ,背景中强度较高的像素虽然灰度值较大 ,但在它所处的局域中与周围背景无明显差异。

图像的大量细节信息体现在图像矩阵A的两个正交矩阵图像的大量细节信息体现在图像矩阵A的两个正交矩阵U和V中，口含有弱小目标的红外图像的大量弱小目标信息体现在的前几个最大的奇异值所对应的□□□□U=[U1,U2,…，Um]、右奇异向量V=[V,V,…,V]□□□□ 1(a)□□□□□□□□□□□原]□□□□□ ，对其曲式12(3)作奇异值分解。图1(b)~(e)别是采用不同数目的非零奇异值重构的结果 ,□□□□□□□□□□值所对应的左、右奇异向量的变化将引起图像的巨大变化 ,而图像对较小的部分奇异值所对应的奇异向量的变化是不敏感的。 这主要是因为背景杂波的能量比较小,所以它对应的奇异值比较小 ,可以通过去除小奇异值滤掉背景子空间 ,然后在有效的弱小目标子空间上重构原图像 ,就可以实现背景抑制的目的。(al 〃>|仲=KIfiltp'in-Ml1 11filtp'in-Ml1 11rl哨二5(11图i利用svn取不同数目非零奇异值重构图像效果□□□□□□□ ,□□□□□□□□□ ,□□□□□□□□□□□□□□□阵u□□□□□□ ,□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□ ，在式(3)中,□□□□□□□□□□□□□ S可以从1~r中选择 ,它的选择是基于不同的背景复杂度情况 ,为了确定选取多少个奇异值 ,这需要通过实验来比较利用不同奇异值数目重构后图像保存目标信号抑制背景噪声的性能 ,来确定最终选取的重构奇异值数目。 因此,定义了偏差指数来确定在不同应用场合精确重构图像时所需的奇异值数目 ,以达到最佳背景抑制效果。2.3基于偏差指数的非零奇异值数目选择偏差指数定义为重构图像各个像素灰度值与原图像相应像素灰度值差的绝对值与原图像相应像素灰度值比值的平均值 ,其值的大小表示重构图像与原图像平均灰度值的相对差异 ,它用来反映重构图像与原图像在弱小目标信息上的匹配程度和将原高空间分辨率图像的细节 (即目标信息 )传递给重构图像的能力 ,如式TOC\o"1-5"\h\z(4)所示 :£工l[A(i,j)-A(i,j)]/A(i,j)1p=^14^ (4)r MXN其中,pD偏差指数，A为原始含有弱小目标的图像 ，ri使用r□□□□□

值个数,MDN为一帧图像的大小。在图2中,将原始图像分解后采用不同的奇异值重构发现在偏差指数达到 0.997时对应的重构图像弱小目标信息保存最好,通过大量的统计实验,所获得的、实验结果与分析3.1评价准则本文选用了五个参数来衡量算法保存目标信号、去除复杂地面背景的有效性。它们是 :(a)信噪比SNR=|Gt-Gb|/obc□5)(b)信噪比增益 ISNR=SNRout/SNRin (6)(c)对比度,通过大量的统计实验,所获得的、实验结果与分析3.1评价准则本文选用了五个参数来衡量算法保存目标信号、去除复杂地面背景的有效性。它们是 :(a)信噪比SNR=|Gt-Gb|/obc□5)(b)信噪比增益 ISNR=SNRout/SNRin (6)(c)对比度SCR=|Gt-Gb|/|Gt+Gb|D100% (7)(d)对比度增益 ISCR=SCRout/SCRin (8)(e)背景抑制因子 BSF=oout/ain(9)式中，Gt,Gb分别是目标的平均灰度值和□□□□□值;SNRin,SNRoEiiD□□□□□□□□□□□□□□□□□□□□□□□□□□□□ ;abc□□□□□□□□□□域内灰度值的均;SCRin,SCRo分1别代;ain□aout分别为原始图像的均方差和滤波后图像的均方差。3.2结果分析选择了三组有代表性的单目标红外图像序列进行实验验证。实验图像大小为128]128像素，像素灰度8位,□□□□ 2以内,□□□□7%左右□基于CPU内存为512M的PC机，通过matlab7.000的图像□□□□□□□□□实验。具体过程如下：图像读取功能主要包括： 在用户界面的右上输入图像地址， 通过该地址可以读取图像数据，对图像进行处理□ 当选中处理图像后, 通过get函数将地址值赋给strl变量。在 callpush□□,□□imread□□□□□□□□□□□□□□的图像。图 3分别为三个序列中的一幅原图 ,图4为采用二维最小均方误差算法的处理结果 ,图5为本文算法的处理效果。处理前后评价参数数值如表 1所示。□□□□□□□□□□□ l，:lLMS算法)自适应滤波算法，该算法是一种以期望□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□：一个是滤波过程， 另一个是自适应过程。 在滤波过程中，自适应滤波器计算其对输入的响应，并且通过与期望响应比较， 得到估计的误差信号。 在自适应过程中，系统估计误差自动调整滤波器自身的参数。该算法的迭代公式如下：e(n)=d(n)-y(n)y(n)=X(n)TW(n)X(n)=[x(n),x(nT),…x(n—L+1)]TW(n+1)=W(n)+2ue(n)X(n)表示自适应迭代下一时刻权系数矢量由当前时刻权系数矢量加上以误差函数为比例因子的输入矢量W(n)t=[W(n),W(n),W(n)]表示时刻n的自适■波权系数L——滤波器的阶数0 1 L-1d(n)——期望输出值y(n)——实际输出值e(n)——误差u——步长因子，为控制稳定性和收敛速度的参数表1不同背景下不同方法的性能比较原图信息二雄最小均方估计市文方法编号4Atin,tS.\liISCtiESTJSAliISCfitiSi'KSTa1.7210.K?M312.923.XIJ0.934.hl10.977,与0.01?II1.S?4.W2.2<)3.74?.960.96?.0316.7110.060.016c1.5.262.2.6X3.ys0.<)64.7211.?47.010.016从实验结果可以看出 ,在原始图像中,目标仅占一个或几个像元 ,呈现为点状□□□□□□□ ,□□□□□□□□□□□ ，特别是图3(b)和图3(c)中的目标处在云层和地面干扰之中。二维最小均方误差算法处理后 ,目标虽然被增强了,但同时背景也被增强 ,而没有被较好地抑制 ,这样必然会影响后续的检测处理。 而经本文算法处理后,图像信噪比增益提高了 4倍以上,图像对比度增益提高 10倍以上,背景抑制因子达到了 7倍以上,较好地抑制了背景,很好地保存并增强了弱小目标信号。特别是对第 2,3图像序列,目标受到大面积云层和地面建筑物的干扰,与二维最小均方误差算法结果相比较 ,本文算法不但平滑了云层和地面建筑物内部 ,而且将云层和地面建筑边缘也平滑掉 ,使得图像整体对比度得到了很大的改善。 所以,在不同的背景、 不同对比度和信噪比发生变化的条件下 ,本文算法都可以得到较好的处理效果。4、结论矩阵的奇异值分解实际是矩阵在酉等价下的一种标准形式， 鉴于本文中 2.1中的( 1)(2)(3)式，本文提出的基于奇异值分解的红外弱小目标背景抑制算法，□□□□□ □3口式中的与两个酉矩阵所对应的左奇异向量 U=[U,U,…，U]、□□□□V=[V,V,…,V],□□□□□□□□□□□□□□□, mm 1借助于前几 m个较大的非零奇异值在有曲弱小目而空间上重构原曲 ,□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□ 0.02s,□□□□□□□量较小,便于算法软硬件的实时实现。但该算法也存在一个缺点：只能处理一些固定不动的红外弱小目标，而对于运动的红外弱小目标还不能得到较理想的处理，有待于进一步的研究。参考文献：[1]罗军辉，姬红兵，刘靳.一种新的红外小目标检测算法及其应用 [J].半导体光电,2007,28(2):290-293.[2]周铭，许少辉.一种红外小目标的图像检测方法 [J].半导体光电 ,2004,25(3):224—227.[3]苏赋，杨文淑，徐智勇 .红外小目标小波多尺度相关检测方法 [J].半导体光电 ,2007,28(4)：5