一种基于迭代外推的4f光学系统图像复原方法.doc

一种基于迭代外推的4f光学系统图像复原方法李鹏，田逢春，陈建军，李显利（重庆大学通信工程学院，重庆，400030)摘要：大多数图像复原算法需经多次迭代运算才能获得满意的结果,其耗时较长。针对该问题，本文引入4f迭代函数系统来实现图像迭代运算，提出了一种理论模型，并结合迭代外推算法进行频谱的延拓，对图像进行降噪处理。实验结果证明，该方法能有效改善图像质量、提高图像峰值信噪比。关键词：图像；复原；噪声；迭代外推；4f系统中图分类号：O438.2 文献标识码：AAn Improved Image Restoration Method for 4f Optical System Based on Iterative ExtrapolationLI Peng, TIAN Feng-chun,CHEN Jian-jun,LI Xian-li (College of Communication Engineering, Chongqing University, Chongqing 400030, CHN)ABSTRACT：Most image restoration algorithms are subject to many iteration to get a satisfactory result, time-consuming long. To solve the problem, we introduces the 4f iterative system to achieve the image iterative calculations, proposed a theoretical model, combined with iterative spectrum extrapolation method extension, the image noise reduction processing.Experimental results show that for low-noise image that the method can effectively improve the image quality, improve the image of peak signal to noise ratio.KEYWORDS：image; restoration; noise; iterative extrapolation; 4f system 1引言：与传统信息手段相比，光信息处理有速度快、并行性和大容量的优点。但由于光信息处理在图像获取过程中受到许多因素的影响会导致图像质量的下降,如光学系统的衍射、像差、离焦、噪声干扰、光学器件口径限制等,它们会造成图像的模糊和失真。目前研究者根据图像退化模型建立了一些约束条件和最优准则,提出了一系列复原方法,如逆滤波法、约束最小平方滤波法等线性方法和凸集投影法、最大似然法等非线性方法1-2,针对光学系统的带宽有限这一特点，本文采用波场延拓法3对图像频谱尽量外推，使图像质量得到提高。同时对4f迭代系统在理想和有噪声两种情况进行了分析，最后通过4f系统实拍光学图片和经过消除光源不均匀性4预处理的4f系统图片验证该方法的有效性，图像的主观质量和峰值信噪比（PSNR）都得到了较大提升。2 迭代外推算法52.1 迭代外推算法基本原理迭代外推算法是在空域和频域之间来回迭代，在空域中和频域中进行修正以加强预先就知道的知识或实际测得的数据。图1是该算法的实现框图。收稿日期：2000-00-00； 收到修改稿日期：2000-00-00基金项目：中国博士后科学基金（20080430096）资助项目作者简介：李 鹏（1982-)，男，江西新余，重庆大学信号与信息处理专业硕士生，主要从事光学信息处理方面的研究。E-mail:cqmobile_。田逢春（1963-），男，重庆人，教授，博士生导师，主要研究领域为信号处理、信息光学、生物信息学、通信技术。E-mail：FengchunT（通信联系人）。陈建军（1978-），男，浙江台州，博士，主要研究领域为图像处理。图1 迭代外推算法的框图Fig.1 Block diagram of iterative extrapolation method在频域中，成像系统的通带之内的物谱是默认已知的，因为这个数据可以实际测得。首先对图像进行第一次傅里叶变换，可以发现处于成像系统通带之内的那一部分物谱，掌握了这两组数据，就可以实行迭代。由于空间的截断改变新像的频谱，引入了成像系统通带外的频谱分量，此外通带内的频谱分量也会变化。重复空间限界的过程，再次做变换，重复对已知频谱分量的加强。这种方法的优点是引入了超出通带外的频谱分量，并使它们得到逐步改善，在没有噪声和一定噪声时，这个算法是收敛的6。2.2 一种频域波场延拓的方法3设光学系统中的介质为各向同性的，取三维标量光波方程作为延拓的基本方程： （1）式中为三维光波场，为介质的折射率，为光速，为时间。为了简化问题的复杂度，我们这里把光学系统传播介质都看成了均匀介质，忽略了光学系统中透镜和空气介质折射率的差异。将（1）式对、分别进行傅里叶变换，考虑到对于某个坐标轴的逆傅氏变换得到的结果。 （2） （3）将（2），（3）两式带入（1）可得： （4）式中为关于、的三维傅里叶变换，令 （5）由前面的假设，在方向上是常数，可得 （6） 当时，即有相位移延拓公式： （7） 将相位移延拓公式降维至二维可得 （8） （9）通过式（8）和式（9），我们可以发现通过对变化域上的频谱的进行适当的相位改变，可以对光学系统的频谱有一定的延拓作用，可以接收到图像更多的频谱信息，使得图像的画质和PSNR值都得到一定的提高。3 4f系统及其噪声模型31经典的4f系统6及其噪声模型信息光学中典型的f系统如图2所示， f为透镜焦距，S为相干点光源，Lo为准直透镜，P1为输入平面，其上放置空间光调制器（SLM）来加载输入图像。该输入经傅里叶透镜L1做傅里叶变换， 在P2面得到其频谱。在P2 面放置另一个SLM，其上加载所需频域滤波器，经滤波后再通过第二个傅里叶透镜L2做傅里叶逆变换，在P3面上得到输出图像，用CCD采集该图像并转换成电信号进行后处理。图2 典型的4f光学系统 Fig.2 Typical 4f optical system实际f系统的噪声情况和系统传递函数非常复杂，包含着系统器件的物理特性。为简化系统，将f系统看成一个整体，仅将加性噪声分离出来。系统的简化模型如图所示。图中，传递函数定义为，理想情况下它应是冲激函数，但是由于透镜孔径等因素的影响，实际上是一个具有低通特性的函数。基于该简化模型，系统的输出可写为 （0）式中为卷积符号。图3 4f系统简化模型Fig.3 Simplified model of 4f system3.2 4f迭代系统8及其噪声模型通过对频谱相位的调节，进行多次迭代最终实现系统通带的延拓，显然典型的光学4f系统无法满足这个需求。因此引入二反射镜型反馈系统来实验迭代功能，二反射镜型反馈系统是在常规的相干光处理系统的基础上构成的。如图4所示，图中实线表示正向通道，虚线表示反馈通道，输入面P1和输出面P3上各置一块反射率较高、透射率较低的分光镜M1和M2，它们相对与光轴对称的倾斜一个较小的角。L1和L2分别为正向光路和反馈回路的傅里叶透镜。由于我们采用的迭代算法是对反馈的频谱面进行相位的调节，故假设在P2面上针对反馈回路上放置相位检测和相位调节装置，对于反馈回路频谱进行调节后进入正向通道，和原来的输入图像迭加成为一个新的输入，再进入反馈回路，经过若干次迭加后系统得到稳定的输出，从而实现对于图像质量的提高。图4 4f迭代系统Fig.4 The 4f iterative system我们采用的迭代外推算法是针对于反馈变化域相位调节的一种算法，所以噪声模型是在频域中讨论的。现将二反射镜型反馈系统看成一个理想系统，仅将加性噪声分离出来。系统的简化模型如图4所示。图4中，传递函数定义为，理想情况下它应是冲激函数，但是由于透镜孔径等因素的影响，实际上是一个具有低通特性的函数。基于该简化模型，系统初始输出为 （11）式中为系统的初始频谱输入，为系统的初始频谱输出，为系统噪声的频谱。第一次迭代输出为 （12）图5 4f迭代系统简化模型 Fig.5 Simplified model of 4f iterative system4 实验的结果与分析4.1仿真图像实验我们通过仿真得到的4f系统的输出作为迭代系统的初始输入，进行多次迭代，来验证迭代外推算法对于图像PSNR和图像质量是否有提高。第一组：实验中分别选用 Boat、Lena和Baboon作为无噪模拟4f系统的输入，此时模拟4f系统的输出。同时把作为迭代系统的初始输入，通过多次迭代，最终得到稳定图像。仿真结果如表1所示： 表1 基于无噪模拟4f系统迭代外推前、后图像的PSNR及其差值 Tab.1 Before and after iterative extrapolation image of PSNR based onnon-noise simulation 4f system and their difference 信噪比图像仿真获得图像迭代获得图像PSNR增加值Boat61620.897dB29.012dB8.1155dBLena61622.678dB30.521dB7.843dBBaboon61620.95dB28.893dB7.9433dB通过表1仿真结果可以看出，无噪的模拟4f系统输出经过多次迭代，图像的PSNR值得到了较大的提升。 (a) (b) 图6 仿真图像和迭代图像对比 Fig.6 Simulation images and iterative image contrast其中图(a)基于无噪模拟4f系统的输出图像，图（b）则是使用图（a）作为迭代系统初始输入，经过多次迭代得到的图像。由于在传输过程中没有引入噪声，所以两图的图像质量都很好，在细节的表现上图（b）略好于图（a）。第二组：实验中分别选用 Boat、Lena和Baboon作为含噪模拟4f系统的输入，此时模拟4f系统的输出。均值为0，归一化方差为0.001的高斯白噪声。同时把作为迭代系统的初始输入，通过多次迭代，最终得到稳定图像。仿真结果如表2所示。表2 基于含噪模拟4f系统迭代外推前、后图像的信噪比Tab.2 Before and after iterative extrapolation image of PSNR based onsimulation 4f system and their difference 信噪比图像仿真获得图像迭代获得图像PSNR增加值Boat61618.259dB20.103dB1.8437dBLena61619.227dB20.384dB1.1571dBBaboon61618.271dB20.034dB1.7628dB通过表2仿真结果可以看出，含噪的模拟4f系统输出经过多次迭代，图像的PSNR值得到了一定的提升。 图7仿真图像和迭代图像对比Fig.7 Simulation images and iterative image contrast其中图(a)基于含噪模拟4f系统的输出图像，图（b）则是使用图（a）作为迭代系统初始输入，经过多次迭代得到的图像。由于在传输过程中引入噪声，所以两图的图像质量比图（6）中的无噪图片都要差些，总体而言两图细节差异不大。通过两组仿真实验结果的对比，我们可以看出在无噪或者噪声较低的系统中，通过频谱延拓的方法，进行多次迭代，最终得到的稳定图像相对于未迭代图像的PSNR值确实有较大的提高，同时若系统存在一定的噪声影响，通过该方法也可以提高图像的PSNR值，但提高的数值远比无噪声系统时要低。 (a) (b)图8 无噪模拟4f系统输出图像频谱Fig.8 Spectrum of the non-noise simulation 4f system output image其中(a)图为无噪模拟4f系统未进行迭代的图像频谱，(b)图为经过迭代后输出图像频谱。从迭代前后的频谱图上可以看出,经过多次迭代后，对于通带内的频谱分量得到了加强，对于通带外的频谱有一定的延拓作用，从而达到提高图像PSNR的作用。4.2 光学4f系统实测图像实验实验平台如图7，采用DH-HNKW250P型He-Ne激光器作为光源，波长；采用Sony-LCX038空间光调制器，分辨率：1024*768，对角线尺寸为：0.7inch；GCO-0203M型消色差傅里叶变换透镜，焦距400；以及采集设备为Canon EOS350D等设备构成f系统实验平台。输入面加载的图像分辨率均为616 Pixel616 Pixel大小。 图9试验使用的4f光学系统 Fig.9 the 4f system of our experiment为了避免实验数据的特殊性，选取实验效果中等的输出作为实验用的数据。减小在算法的验证时与实际中差异较大的可能性。我们设计了两组实验来验证。第一组：在实际的4f系统中，由于光源的不均匀性和系统噪声等因素的影响，造成系统的实测图不能看成经过理想4f系统得到的输出。因此我们把图像先进行预处理来去除光源的不均匀性，把处理后的图像近似看成为无噪4f系统的输出图像，进行多次迭代。迭代结果如表3所示： 表3 4f系统实测图像、预处理后图像和迭代获得图像Tab.3 Output images of 4f system、pre-processing and iterative extrapolation 信噪比图像原始图像预处理后图像预处理图像为初始输入，迭代获得图像Boat61621.665dB23.538dB26.275dBLena61623.521dB23.692dB25.346dBBaboon61619.541dB20.896dB21.645dB从表3可以看出，由于预处理过程无法做到完全地去除光源的不均匀性等因素的影响，所以使用预处理后的图像作为迭代初始输入，迭代多次后获得图像的PSNR与理想状态时还是存在一定的差距，但PSNR仍有明显提高。第二组：直接用4f系统实际获得图像进行迭代，迭代结果如表4所示： 表4 4f系统实测图像和迭代获得图像Tab.4 Output images of 4f system and iterative extrapolation 信噪比图像原始图像迭代获得图像PSNR增加值Boat61621.665dB21.745dB0.08dBLena61623.521dB23.722dB0.201dBBaboon61619.541dB19.816dB0.275dB由于本组4f系统实测图像受到光源不均匀性和噪声等因素的影响较大，采用这组数据作为迭代系统的初始输入对于图像PSNR值的提升比较有限，同时说明在实验效果一般时，迭代算法对于图像PSNR值提升有限；在实验效果较好时，迭代算法才可对图像PSNR值有较大提升。5 结论针对光学系统的低通特性导致图像的模糊和失真，根据迭代外推算法的基本原理，采用频率域波场外延的方法，通过对频谱进行适当的相位改变，经过多次迭代使得系统获得通带外的部分频谱分量，使得