点扩展函数psf的估计

点扩展函数的估计一般的点扩展函数估计是图像恢复中的一个非常困难的问题，一些常用的方法是“运用先验知识的方法，运用后验知识的方法以及误差参数曲线分析法运用先验知识的方法一般来说，大气湍流、光学系统散焦以及照相机与景物之间的相对运动造成的模糊是图像处理中经常遇到的情况，这类退化的点扩展函数可以根据导致模糊的物理过程(先验知识)来确定。研究表明，对于长时间曝光下大气湍流造成的转移函数(它的傅里叶逆变换即为点扩展函数)可近似为高斯型，即表达成H(u,v)UexpC(U2+v2)5/6其中，C为与湍流性质有关的常数。光学系统散焦造成的转移函数式熟知的“Bes-Sinc”函数，可写成H(u,v)=J1(ndp)/(ndp)其中，p=(u2+V2)2,d为光学系统的散焦点扩展函数(在线性移不变系统中是圆函数)的直径，J,)是第一类一阶贝塞尔函数。下面以摄影中照相机与景物之间的相对移动造成模糊的情况，作为先验知识来确定转移函数的具体例子。假定照相机不动，图像f(x,y)在图像面上移动并且图像f(x,y)除移动外不随时间变化。令x0(t)和y0(t)分别代表位移的x分量和y分量，那么在快门开启的时间T内，胶片上某点的总曝光量是图像在移动过程中一系列相应像素的亮度对该点作用之总和。如果快门开启时间与关闭时间可以忽略不计，且光学系统假设是完善的，且有下列关系存在：Tg(x,y)=Ifxx0(t),yy0(t)dt对其两边取傅里叶变换，得到gs00lglg0T|+g|gfx一X0(t),yy0(t)expj2n(ux+vy)dxdydt根据傅里叶变换的空间位置平移性质可得TG(u,v)=IF(u,v)expj2nux0(t)+vyo(t)dt0T=f(u,v)Iexpj2nux0(t)+vyo(t)dt000定义TH(u,v)=Iexpj2nuxo(t)+vyo(t)dt000那么上式可以表示成G(u,v)=H(u,v)f(u,v)可见，H(u,v)的表达式就是移动模糊的转移函数。如果移动只代表为沿着x方向以速度V作匀速运动，那么有xo=Vt，o(t)=Cly将上式代入H(u,v)表达式，可得TH(u,v)=Iexp(j2nuVt)dtnuV=Texp(-jnuVt)sine(nuVt)运用后验判断的方法如果事先并不知道退化的物理过程，或者这种物理过程过于复杂，难以用来确定h(x,y)，那么可能的办法只有从退化图像本身来估计h(x,y)。例如，若有把握断定原始景物某部位有一个清晰的点，于是那个点在退化图像上的模糊影像就是点扩展函数，天文图片会有这种情况，图片上某颗细小星体的退化图像可用来估计点扩展函数。如果原始景物含有明显的直线，则有时可以从这些线条的退化图像来确定h(x,y)，为了说明这一方法，可假定原始景物中有一条平行于x轴的理想线源记做8(y)，此处6(y)被看作是二维函数，但不依赖于X。该理想线源的退化图像则称为线扩展函数，记做h/y)，可表示成J8(P)h(x-a,y-b)dadp-g-g利用6函数的筛选性质，此式变成h1gh(x一a,y)dagggg对上式做变量置换x-a=x，可得h1(yg)=Ih(x,y)dxg这说明线扩展函数在y方向的分布与位置x无关，即在任何一条与x轴平行的线上，h1(y)的值是一个常数，而hi(y)在y方向上任一点的数值是点扩展函数在该点沿x方向的积分。显然，如果点扩展函数式圆对称函数，则线扩展函数与线源的取向无关；否则，就与线源的取向有关。若h1(y)的傅里叶变换为HJV),则gH1(V)=Ih1(y)exp(j2nvy)dyg但我们知道H(u,v)=h(x,y)exp-j2n(ux+vy)dxdygg如果把u=0代入这一方程并使用上两式可以得到H(0,v)=gh(x,y)dxexpj2nvydy=H这表示平行于x轴的线扩展函数的傅里叶变换是转移函数H(u,v)在频谱平面上验u=0直线所取的值。同理可以证明，与x轴成G的线扩展函数，其傅里叶变换则是H(u,v)在频谱平面上沿斜率为On+90的直线所取的值。因此，如果能断定原始景物含有各种取向o1，e2,，o的线，就能从这种集购物退化的12n图像上推到出H(u,v)沿着过原点具有斜率0+90。，o+90。，o+90的那些12n辐射形直线上的值。如果能肯定点扩展函数是圆对称的，则H(u,v)也是圆对称的。因此盒子要知道沿一条辐射线的H(u,v)的值，就知道它各处的值。如果没有这种先验知识,一般必须求得沿着紧挨在一起的许多辐射线上的H(u,v)值。倘若频谱平面能被足够密集的这种线上的H(u,v)值所覆盖，就能构成H(u,v)的精密近似值。并通过内插法求得频谱面坐标网络交点上的值，然后通过傅里叶变换即可求得h(x,y)。假使原始景物不含有点或者线的内容，然而它可能含有明显的界线。现在将要证明界线的退化图像的导数等于平行于该界线的线源的退化函数。一条沿x轴的理想界线在数学上可用S(y)表示，这里的S(y)单位阶跃函数,即1,S(y)=01,y0y0设h(y)是该界线的退化图像，那么shs(x,y)=h(xa,yP)S(P)dadpggh(a,P)S(y卩)dadp由于平行于h(y)代替，变成sx轴的界线的退化图像与x无关，所以上式中的hs(x,y)可用rgrgjh(a,卩)S(y-卩)dadpgg上式中两边对y取偏导数并在右边互换积分和微分算符的次序，得到dh(y)fgfgd=jjh(a，卩)dyS(y一卩)dad卩ggrgrgjh(a,P)8(y一卩)dadp-g-gh(a,y)da-g比较(y)=jgh(x,y)dx和上式，可以看出1丿一gdh(y)/dy=h(y)上式表明，一条线的退化图像就是平行于此线的界线退化图像的导数，因此，若原始景物中含有各种取向的界线，则可用前面讨论过的方法由这些界线退化图像的导数来确定点扩展函数，事实上，由于偏微分算子d/dy是线性位移不变算子，所以根据上式可得dh(y)dfgh1(y)=dy=dyjHs(v)exp(j2nvy)dvLgg=IHs(v)(j2nv)exp(j2nvy)dvg式中Hs(v)是hs(y)得傅里叶变换，由傅里叶变换定义可以看出斗(v)=j2nvH(v)H(0,v)=j2nvHs(v)还有一个方法可从退化图像本身估计转移函数H(u,v),把退化图像分成n个大小相等的子图像g.(x,y),i=12,n，假设点扩展函数取值范围与上述子图i像的尺寸相比足够小，那么对于每个子图像，可以得到gi(x,y)=Jfi(a,卩)h(x-a,y-卩)dad,-g-g对等式两边取傅里叶变换，有Gi(u,v)=H(u,v)Fi(u,v),=1,2,n对于i取乘积变换，得nn口9(叩)=口Fi(u,v)Hn(u,v)i=1i=1或者1/ni=11/ni=1若原始景物各子图像内灰度起伏足够大，且各子图像之间灰度也有相当大的差异，则可期望上式右边的分母接近于一个常数(也即与u、v无关)。另一方面，Gi(u,v)也可由退化子图像gi(x,y)求得。于是根据上式即可估计H(u,v)，从而求得点扩展函数h(x,y)。误差-参数曲线分析法对于点扩展函数可用由某一参数来表征的退化方式(如线性移动和散焦等)，这里介绍误差-参数曲线分析方法来估计点扩展函数。如果点扩展函数可用某一个参数表征，点扩展函数的估计就变成了对应参数的估计。比如散焦和移动分别用散焦半径r和移动距离d来表征。现在把点扩展函数写成h(a)的形式。误差-参数曲线分析法是用下面的方法产生一个误差-参数曲线，通过对曲线的分析来决定点扩展函数的参数，其基本步骤如下。(1) 选定一个参数搜索范围，设a0为该范围的初始位置，Aa为搜索步长，K为搜索步数。(2) 对i=1,2,-k，有a=a0+(i-1)Aa由参数a得到点扩展函数h；由h和退化图像y，实施恢复算法得到恢复图像x；计算恢复误差E=|yx*h|2或其他误差度量。判断Ee(e为预先设定的阈值)。是，执行(3)；否，则回到(2)。(3) 做出E-a曲线，由此判定实际退化参数的a值，并获得实际的点扩展函数。在上述算法中，原则上任何一种有效的恢复算法都可以使用，同样，恢复的误差可以E=|y-x*h|2，也可以是其他任何更为合适的度量误差。在步骤(3)中，观察误差-参数曲线，从曲线右端(最大的参数值)开始，曲线变化率开始显著减小的位置对应的即为真实参数的近似。在一定的信噪比下，误差-参数曲线分析法对于散焦和移动的情形可以进行良好的点扩展函数的估计，但对于其他情形(比如高斯退化)，并不能简单地加以推广，还必须辅之以其他的方法，另外，需要人为参与，这也是该方法的缺点。