点扩展函数psf的估计.docxVIP

点扩展函数的估计 一般的点扩展函数估计是图像恢复中的一个非常困难的问题，一些常用的方法是“运用先验知识的方法，运用后验知识的方法以及误差—参数曲线分析法。 运用先验知识的方法 一般来说，大气湍流、光学系统散焦以及照相机与景物之间的相对运动造成的模糊是图像处理中经常遇到的情况，这类退化的点扩展函数可以根据导致模糊的物理过程（先验知识）来确定。研究表明，对于长时间曝光下大气湍流造成的转移函数（它的傅里叶逆变换即为点扩展函数）可近似为高斯型，即表达成 H 其中，C为与湍流性质有关的常数。光学系统散焦造成的转移函数式熟知的“Bes-Sinc”函数，可写成 H 其中，p=u2+v21 下面以摄影中照相机与景物之间的相对移动造成模糊的情况，作为先验知识来确定转移函数的具体例子。 假定照相机不动，图像fx,y在图像面上移动并且图像fx,y除移动外不随时间变化。令x0(t)和y0(t)分别代表位移的 g 对其两边取傅里叶变换，得到 G = 根据傅里叶变换的空间位置平移性质可得 G =f(u,v) 定义 H 那么上式可以表示成 G(u,v)=H(u,v)f(u,v) 可见，H(u,v)的表达式就是移动模糊的转移函数。如果移动只代表为沿着x方向以速度V作匀速运动，那么有 x 将上式代入H(u,v)表达式，可得 H = =T 运用后验判断的方法 如果事先并不知道退化的物理过程，或者这种物理过程过于复杂，难以用来确定hx,y，那么可能的办法只有从退化图像本身来估计h 如果原始景物含有明显的直线，则有时可以从这些线条的退化图像来确定hx,y，为了说明这一方法，可假定原始景物中有一条平行于x轴的理想线源，记做δ(y)，此处δ(y)被看作是二维函数，但不依赖于x。该理想 h 利用δ函数的筛选性质，此式变成 h 对上式做变量置换x- h 这说明线扩展函数在y方向的分布与位置x无关，即在任何一条与x轴平行的线上，h1y的值是一个常数，而h1y在y方向上任一点的数值是点扩展函数在该点沿 若h1y的傅里叶变换为H H 但我们知道 H 如果把u=0代入这一方程并使用上两式可以得到 H 这表示平行于x轴的线扩展函数的傅里叶变换是转移函数Hu,v在频谱平面上验u=0直线所取的值。同理可以证明，与x轴成θ的线扩展函数，其傅里叶变换则是Hu,v在频谱平面上沿斜率为θn+90°的直线所取的值。因此，如果能断定原始景物含有各种取向 如果能肯定点扩展函数是圆对称的，则Hu,v也是圆对称的。因此盒子要知道沿一条辐射线的Hu,v的值，就知道它各处的值。如果没有这种先验知识，一般必须求得沿着紧挨在一起的许多辐射线上的Hu,v值。倘若频谱平面能被足够密集的这种线上的Hu,v值所覆盖，就能构成 假使原始景物不含有点或者线的内容，然而它可能含有明显的界线。现在将要证明界线的退化图像的导数等于平行于该界线的线源的退化函数。 一条沿x轴的理想界线在数学上可用S(y)表示，这里的S(y) 单位阶跃函数，即 S 设hs h = 由于平行于x轴的界线的退化图像与x无关，所以上式中的hsx,y可用 h 上式中两边对y取偏导数并在右边互换积分和微分算符的次序，得到 ? = = 比较h1 ?h(y)?y 上式表明，一条线的退化图像就是平行于此线的界线退化图像的导数，因此，若原始景物中含有各种取向的界线，则可用前面讨论过的方法由这些界线退化图像的导数来确定点扩展函数，事实上，由于偏微分算子??y是线性位移不变 h = 式中Hs(v)是hs H H 还有一个方法可从退化图像本身估计转移函数Hu,v，把退化图像分成n个大小相等的子图像g g 对等式两边取傅里叶变换，有 G 对于i取乘积变换，得 i=1 或者 H 若原始景物各子图像内灰度起伏足够大，且各子图像之间灰度也有相当大的差异，则可期望上式右边的分母接近于一个常数（也即与u、v无关）。另一方面，Gi(u,v)也可由退化子图像gi(x,y)求得。于是根据上式即可估计 误差-参数曲线分析法 对于点扩展函数可用由某一参数来表征的退化方式（如线性移动和散焦等），这里介绍误差-参数曲线分析方法来估计点扩展函数。 如果点扩展函数可用某一个参数表征，点扩展函数的估计就变成了对应参数的估计。比如散焦和移动分别用散焦半径r和移动距离d来表征。现在把点扩展函数写成h(α 选定一个参数搜索范围，设α0为该范围的初始位置，?α为搜索步长 对i=1,2,? α= 由参数α得到点扩展函数h； 由h和退化图像y，实施恢复算法得到恢复图像x； 计算恢复误差E=y-x*h 判断E≤ε(ε为预先设定的阈值)。是，执行(3)；否，则回到( 做出E-α曲线，由此判定实际退化参数的α值，并获得实际的点扩展函数。 在上述算法中，原则上任何一种有效的恢复算法都可以使用，同样，恢复的误差可以E