红外图像的边缘提取

摘 要：提出了一种基于人眼微动机理的边缘提取算法。通过模拟眼球的微动，提取图像的微动边缘，同时为了减少伪边缘的产生，对其微动边缘图像进行均值滤波处理，最后应用非极大值抑制和双阈值检测边缘连接提取图像的二值化边缘。实验结果显示，该算法边缘提取效果较好，达到了预期的效果。

　　传统算法大多都是基于局部窗口的微分梯度算子，对噪声敏感，因此不适合处理受噪声干扰严重的图像。由于人眼微动具有超分辨率的特性，对红外图像处理提供了良好的启示。本文就是根据人眼微动机理研究红外图像边缘的提取，并通过实验加以验证。实验结果表明，基于人眼微动机理的红外图像边缘提取算法不但能够精确提取图像的边缘，同时能够很好地减少图像中的伪边缘，具有良好的边缘提取效果。

　　1 人眼微动成像原理

　　早在1952年，DITCHBURN和GINSBORG等人就注意到人眼在固视状态下具有无意识的微小运动，即人眼微动[5,6],它分为三种模式：高频振颤、飘移运动和闪动。由于前两种幅度不大，同时目前的眼球检测技术难以精确地测量，因此这里提到的人眼微动主要指闪动。国外学者对眼球微动进行了大量的实验研究。实验表明，人眼在停止眼球所有运动的时候，静止的图像将突然变模糊继而消失，眼球微动与视觉的产生存在直接的关系，且双目微动优于单目微动。近年来，CONDE M等[7]学者通过测量图像消失或再现前后微动发生的概率、速率以及振幅的变化，进一步揭示了眼球微动与图像消失有着直接关系，微动与固视图像的清晰度有因果联系。

　　人眼微动的成像原理[8,9]为：人眼微动使感兴趣信息更准确地落入视网膜的中央凹区，信息越精确地落于中央凹区，人眼所感受到的图像锐化程度越强；人眼微动会使近凹区反应增强，它们通过横向连接具有抑制作用的神经网络对中央凹区的信息进行修正，如一次修正不够理想，人眼微动会使信息重新更精确地回到中央凹区，同时随着信息尺度的大小调整微动幅度，如此反复直到精确地辨识出信息为止。信息尺度越小，微动幅度就越小，反之亦然。同时随着微动速率的提高，对应视网膜上感受视野的锐化能力就会越强。本文通过对图像进行平移来近似模拟人眼这种微动机制，并将其应用于红外图像的边缘提取。

　　2 算法基本过程

　　2.1 算法基本原理

　　算法可以理解为有基本运算及规定的运算顺序所构成的完整的解题步骤。或者看成按照要求设计好的有限的确切的计算序列，并且这样的步骤和序列可以解决一类问题。

　　人眼主要依靠微动机制来分辨图像的边缘，眼球微动幅度越小，图像的边缘越细致，随着幅度增大，图像边缘线条变粗，但是对大尺度边缘突出能力强于小幅度的情形。人眼微动理论具有强大的边缘提取能力。结合视网膜节细胞对方向的敏感性，首先选择人眼微动的方向，然后通过微动图像计算微动方向的边缘图像，之后各个微动方向边缘图像进入竞争环节，竞争的结果则为各个微动方向最优的整体灰度边缘图像，最后进行二值化处理生成二值边缘图像。

　　设原始图像为f（x,y），则在某一微动方向上图像的微动边缘图像由下式表示：

　　g（x,y）={f（x,y）-f?兹（x+k1·？驻x,y+k2·？驻y）} （1）

　　式中，？驻x、？驻y分别表示图像向x、y方向移动的距离单元，k1、k2表示移动的大小，？兹表示移动的方向，计算公式为：

　　?兹=arctan（k2/k1） （2）

　　以r （x,y）表示竞争网络的输出，则竞争灰度边缘图像可表示为：

　　r（x,y）=max{gi（x,y）|?坌i} （3）

　　设定阈值t，得到二值化边缘图像：

　　R（x,y）=1 r（x,y）≥t0 r（x,y）t （4）

　　2.2 计算步骤

　　（1）输入原始图像A，通过对原始图像A在微动方向上平移，生成综合微动图像F。F=[Ah,Av,Ad]，其中h、v、d分别代表水平、垂直和倾斜方向。本文分别将图像A向8个方向平移，移动距离为一个像素单位。

　　（2）计算各微动方向的边缘图像H：

　　Ci=Fi-A， i=h,v,d （5）

　　（3）计算竞争灰度边缘图像H：

　　H=max（Ci）， i=h,v,d （6）

　　（4）将竞争灰度边缘图像H重新量化到[0，255]。

　　（5）为了减少伪边缘的产生，对竞争边缘图像H进行均值滤波处理：

　　G=mean（H） （7）

　　（6）对量化滤波后的灰度边缘图像，采用非极大值抑制和双阈值检测边缘连接处理，得到二值边缘图像。

　　2.3 非极大值抑制

直接对经过量化滤波的竞争灰度边缘图像进行二值化操作并不能准确地提取出图像的边缘，因此需要对经过量化滤波的竞争灰度边缘图像的幅值进行非极大值抑制来进一步确定边缘点。若图像G（x,y）上（i,j）像素点的边缘强度G（i,j）小于沿平移线方向上的两个相邻像素点的边缘强度，则认为该像素点为非边缘点，将其灰度值设为