4记录和整理实验报告。图像降质的数学模型
图像复原处理的关键问题在于建立退化模型。输入图像f(x, y)经过某个退化系统后输出的是一幅退化的图像。为了讨论方便， 把噪声引起的退化即噪声对图像的影响一般作为加性噪声 考虑， 这也与许多实际应用情况一致，如图像数字化时的量化 噪声、 随机噪声等就可以作为加性噪声，即使不是加性噪声而 是乘性噪声， 也可以用对数方式将其转化为相加形式。
原始图像f(x, y) 经过一个退化算子或退化系统H(x, y) 的作
用， 再和噪声n(x,y)进行叠加，形成退化后的图像g(x, y)。图2-1表示退化过程的输入和输出的关系，其中H(x, y)概括了退化系统的物理过程，就是所要寻找的退化数学模型。
图2-1 图像的退化模型
数字图像的图像恢复问题可看作是： 根据退化图像g(x , y)和退化算子H(x , y)的形式，沿着反向过程去求解原始图像f(x , y)， 或者说是逆向地寻找原始
n (x, y)f (x, y)H(x, y)g (x, y)图像的最佳近似估计。图像退化的过程可以用数学表达式写成如下的形式：
g(x, y)=H［f(x, y)］+n(x, y) (2-1) 在这里，n(x, y)是一种统计性质的信息。在实际应用中， 往往假设噪声是白噪声，即它的频谱密度为常数，并且与图像不相关。
在图像复原处理中， 尽管非线性、 时变和空间变化的系统模型更具有普遍性和准确性，更与复杂的退化环境相接近，但它给实际处理工作带来了巨大的困难， 常常找不到解或者很难用计算机来处理。因此，在图像复原处理中， 往往用线性系统和空间不变系统模型来加以近似。这种近似的优点使得线性系统中的许多理论可直接用于解决图像复原问题，同时又不失可用性。
2.2匀速直线运动模糊的退化模型
在所有的运动模糊中，由匀速直线运动造成图象模糊的复原问题更具有一般性和普遍意义。因为变速的、非直线运动在某些条件下可以被分解为分段匀速直线运动。本节只讨论由水平匀速直线运动而产生的运动模糊。
假设图象
f?x,y?有一个平面运动，令x?t?和y0?t?分别为在x和y方向上
0运动的变化分量，T表示运动的时间。记录介质的总曝光量是在快门打开后到关闭这段时间的积分。则模糊后的图象为：
Tg(x,y)??0f?x?x0(t),y?y0(t)?dt （2-2）
式中g(x,y)为模糊后的图象。以上就是由于目标与摄像机相对运动造成的图象模糊的连续函数模型。
如果模糊图象是由景物在x方向上作匀速直线运动造成的，则模糊后图象任意点的值为：
T?g?x,y???0fx?x0(t),y?dt （2-3）
式中x0?t?是景物在x方向上的运动分量，若图象总的位移量为a，总的时间为T，则运动的速率为x0?t?=at/T。则上式变为：
Tg(x,y)??0?at （2-4） fx?,y?dt?T?????以上讨论的是连续图象，对于离散图象来说，对上式进行离散化得：
?at （2-5） g(x,y)??fx?,y??t??T?i?0L?1?????其中L为照片上景物移动的像素个数的整数近似值。是每个像素对模糊
产生影响的时间因子。由此可知，运动模糊图象的像素值是原图象相应像素值与其时间的乘积的累加。
从物理现象上看，运动模糊图象实际上就是同一景物图象经过一系列的距离延迟后再叠加，最终形成的图象。如果要由一幅清晰图象模拟出水平匀速运动模糊图象，可按下式进行：
1L?1g(x,y)??f(x,y) （2-6）
Li?0这样可以理解此运动模糊与时间无关，而只与运动模糊的距离有关，在这种条件下，使实验得到简化。因为对一幅实际的运动模糊图象，由于摄像机不同，很难知道其曝光时间和景物运动速度。
我们也可用卷积的方法模拟出水平方向匀速运动模糊。其过程可表示为：
g(x,y)?f(x,y)?h(x,y) (2-7)
其中
?1h(x,y)???L??00?x?L?1其它 (2-8)
h(x,y)称为模糊算子或点扩散函数，“*”表示卷积，f(x,y)表示原始(清晰)图象，g(x,y)表示观察到的退化图象。
如果考虑噪声的影响，运动模糊图象的退化模型可以描述为一个退化函数和一个加性噪声项n(x,y)，处理一幅输入图象f(x,y)产生一幅退化图象g(x,y)。
g(x,y)?f(x,y)?h(x,y)?n(x,y) (2-9)
由于空间域的卷积等同于频率域的乘积，所以式(2-9)的频率域描述为：
G(u,v)?H(u,v)F(u,v)?N(u,v) (2-10)
式(2-9)中的大写字母项是式(2-10)中相应项的傅里叶变换。
2.3点扩散函数的确定
不同的点扩散函数（PSF）会产生不同的模糊图象。明确的知道退化函数是很有用的，有关它的知识越精确，则复原结果就越好。首先讨论几个典型的点扩散函数。
2.3.1典型的点扩散函数
运动模糊的点扩散函数：假设图象是通过一个具有机械快门的摄像机获得的。摄像机和拍摄物体在快门打开期间T的相对运动引起物体在图象中的平滑。
假设V是沿x轴方向的衡常速度，时间T内PSF的傅里叶变换H(u,v)由下式给出：
H(u,v)?sin(?vtu) (2-11)
?Vu离焦模糊的点扩散函数：由于焦距不当导致的图象模糊可以用如下函数表示：
H(u,v)?J2221(u,v)ar (2-12)
其中J1是一阶Bessel函数，r?u?v，a是位移。该模型不具有空间不变性。
大气扰动的点扩散函数：大气的扰动造成的图象模糊在遥感和天文中是需要复原的。它是由大气的不均匀性使穿过的光线偏离引起的，以下给出了数学模型，其表达式为：
22H?u,v??e?c(u?v)56 (2-13)
其中c是一个依赖扰动类型的变量，通常通过实验来确定。幂
5/6有时用1代替。
当我们得到一幅退化图象的时候，首先要判断其退化类型然后通过已知的先验知识进行恢复。以下的讨论主要针对运动模糊PSF进行。
2.3.2运动模糊点扩散函数的离散化
对于运动模糊而言，根据相机与目标的相对运动速度，相机的焦距以及相机相对目标的距离等就可以计算出PSF。例如通过计算得到一幅模糊图片的模糊方向是x=6,y=4，连续的PSF为见图2.2(a)。
根据式(2-14)和式(2-15)，可以得到离散化以后的PSF如图2.2(b)。可以看到由于离散化的原因PSF并非是直线。在本文中的PSF均由此方法获得。
?m??y??x???n??ny????m?y???x?0,1,??,n?1 (2-14)
y?0,1,??,n?1 (2-15)
n：y方向的模糊点数；m：x方向的模糊点数。当n>=m的时候用式(2-14)计算当n 3、 运动模糊图象的复原方法及原理 为了抑制退化而利用有关退化性质知识的预处理方法为图象复原。多数图象复原方法是基于整幅图象上的全局性卷积法。图象的退化可能有多种原因：光学透镜的残次、光电传感器的非线性、胶片材料的颗粒度、物体与摄像机间的相对运动、不当的焦距、遥感或天文中大气的扰动、照片的扫描等等。图象复原的目标是从退化图象中重构出原始图象。 运动模糊图象的恢复是图象复原的主要分支之一，它的恢复算法有很多种。有些算法虽然有很好的恢复效果，但算法复杂，恢复时间比较长（如最大熵法）。有些算法虽然计算速度较快，但恢复效果不尽人意（如空间域逆向恢复）。 下面介绍逆滤波、维纳滤波和有约束最小二乘滤波三种恢复方法的原理。 3.1逆滤波复原原理 在六十年代中期，逆滤波（去卷积）开始被广泛地应用于数字图象复原。Nathan用二维去卷积方法来处理由漫游者、探索者等外星探索发射得到的图象。由于和噪声相比，信号的频谱随着频率升高下降较快，因此高频部分主要是噪声。Nathan采用的是限定逆滤波传递函数最大值的方法。 在同一时期，Harris采用PSF的解析模型对望远镜图象总由于大气扰动造成的模糊进行了逆滤波处理，Mcglamery则采用由实验确定的PSF来对大气扰动图
图像退化 图像复原
