一 图像退化
由此,我们给出图像退化的描述(图像退化及其过程描述)如下:
图像的退化是指图像在形成、传输和记录过程中,由于成像系统、传输介质和设备的不完善,使图像的质量下降(变坏)。其典型表现为:模糊、失真、有噪声。产生原因:成像系统像差、传感器拍摄姿态和扫描非线性、成像设备与物体运动的相对运动、大气湍流、成像和处理过程中引入的噪声等。
二 图像修复
针对这些问题,我们需要对退化后的图像进行复原。这是我们本节内容的第二个关键词图像复原,图像复原就是尽可能恢复退化图像的本来面目,它是沿图像退化的逆过程进行处理,也就是如果我们知道图像是经历了什么样的过程导致退化,就可以按其逆过程来复原图像。因此,图像复原过程流程如下:
找退化原因→建立退化模型→反向推演→恢复图像
典型的图像复原是根据图像退化的先验知识,建立退化现象的数学模型,再根据模型进行反向的推演运算,以恢复原来的景物图像。因此,图像复原的关键是知道图像退化的过程,即图像退化模型。并据此采用相反的过程求得原始图像。
针对不同的退化问题,图像复原的方法主要有:代数方法恢复、运动模糊恢复、逆滤波恢复、维纳滤波恢复、功率谱均衡恢复、约束最小平方恢复、最大后验恢复、最大熵恢复、几何失真恢复等。
这里也许同学们会有一个疑问,那就是图像复原和前面讲过的图像增强有什么区别呢?区别如下:
图像增强不考虑图像是如何退化的,而是主观上试图采用各种技术来增强图像的视觉效果。
因此,图像增强可以不顾增强后的图像是否失真,只要达到想要的目视效果就可以。而图像复原就完全不同,需知道图像退化的机制和过程等先验知识,客观上找出一种相应的逆处理方法,从而得到复原的图像。
如果图像已退化,应先作复原处理,再作增强处理。二者的目的都是为了改善图像的质量。
1 图像退化的数学模型
输入图像f(x,y)经过某个退化系统后输出的是一幅退化的图像。为了讨论方便,把噪声引起的退化即噪声对图像的影响一般作为加性噪声考虑。 原始图像f(x,y)经过一个退化算子或退化系统H(x,y)的作用,再和噪声n(x,y)进行叠加,形成退化后的图像g(x,y)。图1 表示退化过程的输入和输出之间的关系,其中H(x,y)概括了退化系统的物理过程,它就是我
们要寻找的退化数学模型。
这里的H 就是退化系统通常,我们假设图像经过的退化系统是线性非时变系统,线性非时变系统具有如下四个基本性质:
H 的四个性质(假设噪声n(x,y) = 0):
<<<原始图像偏移多少,响应的退化图像也偏移多少,说明线性系统在图像任意位置的响应只与在该位置的输入值有关而与位置本身无关根据这些特点,输入信号与其经过线性非时变系统的输出信号之间的关系,以及傅里叶变换的性质,我们得到如下时域及频域的关系表达式:
即:在时域上分析时,原始图像经过退化系统后得到的退化图像g(x,y)等于原始输入图像f(x,y)与系统冲激响应h(x,y)的卷积再加上噪声信号。
在频(率)域上分析时,退化图像的傅里叶变换G(u,v)等于原始图像的傅里叶变换F(u,v)与退化系统的频率响应H(u,v)相乘,再加上噪声信号的傅里叶变换N(u,v)这里我们还要为大家解释一下为什么采用线性位移不变系统模型来描述图像退化过程,主要是如
下三个原因:
1)由于许多种退化都可以用线性位移不变模型来近似,这样线性系统中的许多数学工具如线性代
数,能用于求解图像复原问题,从而使运算方法简捷和快速。
2)当退化不太严重时,一般用线性位移不变系统模型来复原图像,在很多应用中有较好的复原结
果,且计算大为简化。
3)实际上,尽管非线性和位移可变的情况能更加准确而普遍地反映图像复原问题的本质,但在数
学上求解困难。只有在要求很精确的情况下才用位移可变的模型去求解,其求解也常以位移不变
的解法为基础加以修改而成。
换句话说,就是利用了:线性系统的计算简洁快速,同时它也是非线性和位移可变系统的基础,
非线性位移可变系统的求解可以转化为线性移不变求解。
2 图像逆滤波复原
1)逆滤波恢复原理
根据上一节讲解的图像退化与复原的概念以及图像退化模型,如下图所示,原始图像f(x,y)经过退化系统H 以及噪声n(x,y)的影响后得到退化后的图像g(x,y),经过复原系统M 将g(x,y)尽最大可能恢复到近似原图f’(x,y) ,因此,对图像的复原可以看作是图像退化的逆过程。
由图像的退化模型及图像复原的基本过程可见,复原处理的关键在于对系统H 的基本了解。由信号与系统的时域和频域分析与处理理论可知,对退化图像的恢复可以在时域或者频域来完成,时域是通过卷积实现,卷积是积分求和的过程,根据傅里叶变换的时域卷积性质,可以知道,频域处理是通过乘积来实现。乘法计算比卷积计算简单,同时傅里叶变换的快速实现(FFT)使得对图像的频率域处理更为简单。因此,可以在频率域来解决该问题。
设图像退化前的傅里叶变换为F(u,v),退化后的傅里叶变换为G(u,v),系统函数即退化函数的傅里叶变换为H(u,v),从频率域角度看,它使图像退化,因而反映了成像系统的性能。由前面得知的图像退化频域描述模型(如下式)
该过程就是用退化函数除退化图像的傅里叶变换,得到退化前图像的傅里叶变换的估计,因此叫直接逆滤波图像复原,该方法是对经退化函数H 退化的图像进行复原的最简单方法。
由该式可知,即使知道退化函数,也不能准确的复原图像,因为N(u,v)未知,甚至有更糟的情况是如果退化函数是零或是非常小的值时,则噪声与退化函数的壁纸N(u,v) / H(u,v)比较大,很容易支配F(u, v)的估计值,会对逆滤波复原的图像产生很大的影响,有可能使恢复得到的图像与原始f(x,y)相差很大,甚至面目全非。
解决这个问题的一种方法是限制滤波的频率,从频谱图可知,高频分量(通常对应的是噪声信号)的值接近 0,而H(0,0)在频率域中通常是H(u, v)的最高值。因此可缩短滤波半径,使通过的频率接近原点,减少遇到零值的概率。
2)逆滤波复原步骤
逆滤波复原过程可归纳如下:
若噪声存在,为减少噪声对复原信号的影响,则将H-1(u,v)进行调整,
(1) 在H(u,v) =0附近,人为地仔细设置 H-1(u,v)使得复原信号受噪声影响变小。
(2) 使H(u,v)具有低通滤波性质。
通过上面的傅里叶逆变换得到的复原图像。以Lena 图像为例,采用上述逆滤波方法进行图像复原的结果如图所示: