数字图像处理6图象重建

1、主要内容： 1.图像(t xin)退化模型 2.图像复原方法 无约束复原 有约束最小二乘复原 逆滤波恢复法 运动图像复原 维纳滤波复原方法 第五章 图象(t xin)复原与重建* 2022/7/251第五章 图像复原与重建共六十二页3.图象重建(zhn jin)的基本概念4.图象重建方法 解联立方程组方法傅立叶变换法滤波-逆投影法5.典型应用 2022/7/252第五章 图像复原与重建(zhn jin)共六十二页5.1 图像退化(tuhu)模型 一.图像退化1.图像退化 图像在形成、传输和记录过程中，由于成像系统、传输介质和设备的不完善，使图像的质量变坏，这一过程称为图像的退化。2.图像复原

2、图像的复原就是要尽可能恢复(huf)退化图像的本来面目，它是沿图像降质的逆向过程进行。 2022/7/253第五章 图像复原与重建共六十二页 典型的图像复原是根据图像退化的先验知识建立一个退化模型，以此模型为基础，采用各种逆退化处理方法进行恢复(huf)，使图像质量得到改善。可见，图像复原主要取决于对图像退化过程的先验知识所掌握的精确程度。3.图像复原的一般过程：确定退化(tuhu)原因建立退化模型反向推演恢复图像2022/7/254第五章 图像复原与重建共六十二页4.图像复原评价准则: 最小均方准则 加权均方准则 最大熵准则 评价准则是用来规定(gudng)复原后的图像与原图像相比较的质量标

3、准。 2022/7/255第五章 图像复原与重建(zhn jin)共六十二页5.图像复原与图像增强目的都是为了改善图像的质量。图像增强不考虑图像是如何退化的，只通过试 探各种技术来增强图像的视觉效果。图像复原则需知道图像退化的机制和过程的先验知识(zh shi)，据此找出一种相应的逆过程方法，从而得到复原的图像。如果图像已退化，应先作复原处理，再作增强处理。 2022/7/256第五章 图像复原与重建(zhn jin)共六十二页二.图像(t xin)退化的数学模型 g(x,y)=Tf(x,y)成像系统：输入-f(x,y),输出-g(x,y)，成像系统的作用(zuyng)-T退化图像为：线性位移

4、不变成像系统线性：Tk1f1 (x,y)+k2 (x,y)=k1 Tf1 (x,y)+ k2Tf2 (x,y)若:g(x,y)=Tf(x,y),则对任一个f(x,y)和常数a,b,均有：Tf(x-a , y-b )=g(x-a , y-b )即运算结果仅取决于输入值，而与位置无关空间位移不变：2022/7/257第五章 图像复原与重建共六十二页二.图像(t xin)退化的数学模型 1.线性位移(wiy)不变成像系统图像退化模型 g(x,y)=f(x,y)*h(x,y)+n(x,y) g(x,y)退化图像f(x,y)-理想图像h(x,y)-点扩散函数n(x,y)-加性噪声Hf(x,y)g(x,y

5、)n(x,y)降质系统2022/7/258第五章 图像复原与重建共六十二页典型(dinxng)的降质原因(1). 光学(gungxu)衍射降质的H(u,v)其中，l为所考虑光的波长，Do为光曈的直径，Di为出射光曈到影像平面的距离。 这是由于相机聚焦不准确引起的，虽然不聚焦由许多参数决定，如相机的焦距、相机孔的大小、形状、物体和相机之间的距离等，但在研究中为了简单起见，我们用下列函数表示聚焦不准引起的模糊。共六十二页通常模糊算子相当于一个(y )低通滤波器，因此当模糊算子作用于原始图像时，会引起图像中边缘和轮廓的模糊。77均匀二维模糊算子作用于图像Camera的结果如下图所示：原始(yunsh

6、)图像 77均匀二维模糊算子作用后的结果 2022/7/2510第五章 图像复原与重建共六十二页(2). 相机(xingj)与景物相对运动的降质设T为曝光时间，a(t)和b(t)分别为对象曝光期间在 x 和 y 方向位移(wiy)的分量通常在拍摄过程中，相机或物体移动造成的运动模糊可以用一维均匀邻域像素灰度的平均值来表示 。共六十二页(3). 大汽湍流(tunli)造成的图像降质这种模糊经常出现在遥感和航空摄影中，由于曝光(bo gung)时间过长引起的模糊可用高斯点扩散函数来表示： 式中K是一个归一化常数，保证模糊的大小为单位值，2可以决定模糊的程度。2022/7/2512第五章 图像复原与

7、重建共六十二页2.线性位移不变系统离散化的退化模型(mxng) 若对图像f(x,y)和点扩散函数h(x,y)均匀采样就可以得到离散的退化模型。假设数字图像f(x,y)和点扩散函数h(x,y)的大小分别为AB 、CD，可先对它们作大小为MN的周期延拓图像，其方法是添加零。即：2022/7/2513第五章 图像复原与重建(zhn jin)共六十二页把周期延拓的fe(x,y)和he(x,y)作为(zuwi)二维周期函数来处理，即在x和y方向上，周期分别为M和N，则由此得到离散的退化模型为两函数的卷积：加上一个(y )延拓为MN的离散噪声项，则：式中x=0,1,2,M-1，y=0,1,2,N-1。函数

8、ge(x,y)为周期函数，其周期与fe(x,y)和he(x,y)的周期一样。2022/7/2514第五章 图像复原与重建共六十二页为便于讨论(toln)和求解，用堆叠方式将二维信号表为一维向量，即：将f、g和n皆相应各行顺时针旋转90度堆叠成MN维向量，则退化模型可被表示为向量矩阵形式：g=Hf+n f =fe(0,0)fe(0,1) fe(0,N-1) fe(1,N-1)fe(1,0)fe(1,1)fe(M-1,N-1)fe(i,N-1)g =ge(0,0)ge(0,1) ge(0,N-1) ge(1,N-1)ge(1,0)ge(1,1)ge(M-1,N-1)ge(i,N-1)MN行MN行从

9、而(cng r)使 g=Hf共六十二页则 H 为 MNMN 阶的分块矩阵(j zhn)，这一矩阵为MM分块循环阵。H 的任意元素Hj是由h(x,y)第j行循环构成，且Hj为NN 循环阵。MM个分块其子块：NN 维共六十二页因此， H 为分块循环矩阵，分块 Hj 中元素(yun s)的第二个下标也是按循环方式变化的。于是(ysh)，二维离散的降质模型化为：解决降质的问题化为求解线性方程组的问题。可以用线性代数和数值分析的方法进行处理。由于h为循环矩阵，有可能进行对角化处理，简化求解。共六十二页例：设 若忽略噪声，求退化(tuhu)图像。 解：周期(zhuq)延拓M=4,N=5 2022/7/25

10、18第五章 图像复原与重建共六十二页用堆叠(dudi)行列式计算(j sun)g得到退化图像。2022/7/2519第五章 图像复原与重建共六十二页5.2 图像复原方法(fngf)一.代数复原方法图像复原的目的是在假设具备有关g、h和n的某些知识的情况下，寻求估计原图像f的某些方法。本部分讨论(toln)在均方误差最小意义下，原图像f的最佳估计， 2022/7/2520第五章 图像复原与重建共六十二页1. 无约束复原由g=Hf+n式可得退化模型(mxng)中的噪声项为 n未知，有意义的准则函数是寻找一个 使得在最小二乘意义上近似于f,即要使噪声(zoshng)项的范数尽可能小,也就是使 把这一

11、问题等效地看作为求准则函数 关于 最小的问题 2022/7/2521第五章 图像复原与重建共六十二页根据(gnj)范数定义：对向量(xingling)求导的两个性质：设a和b为两个列向量，A为对称矩阵，J 为一标量则：(1) 若 J aTb = bTa， 则：(2) 若 J aTA a， 则：因此可将图像复原问题看作是对 求下式的最小值：为此通过J 对 求偏导数，并将结果设为零而达到。将上述性质用于 对 求偏导，有下面结果：共六十二页于是(ysh)求出：因此(ync)：向量 性质 1 性质 1 性质 2共六十二页若设 MN， 即 H 为方阵(fn zhn)，并设H1存在，则：-逆滤波(lb)恢

12、复法 若H为奇异的，无论H-1或（HTH）-1都不存在。这会导致恢复问题的病态性或奇异性。 共六十二页2. 有约束最小二乘复原 为了(wi le)克服恢复问题的病态性质，常需要在恢复过程中施加某种约束。这就是有约束复原。令Q为f的线性算子，有约束最小二乘法复原问题是使形式为 服从约束条件 的最小化问题。这可以归结为寻找一个 ，使下面准则函数最小。 其中，Q是选来对 进行某种线性运算的矩阵，为一常数(chngsh),叫作拉格朗日系数。2022/7/2525第五章 图像复原与重建共六十二页令： 对 的导数(do sh)为零，有 其中=1/。这是求有约束最小二乘复原图像(t xin)的通用方程式。

13、解得：2022/7/2526第五章 图像复原与重建共六十二页通过指定不同的Q,可以达到不同的复原图像。下面便利用通用方程式给出几种具体恢复方法(fngf)。(1)能量约束恢复若取线性运算 Q =I 则得此解的物理意义是在约束条件为 时，复原图像能量应保持不变。 2022/7/2527第五章 图像复原与重建(zhn jin)共六十二页(2)平滑约束恢复(huf)把 考虑成x,y的二维函数，即同时考虑在x和y两个方向上的差分，那么光滑约束取为原图像的拉普拉斯运算。即用卷积形式表示(biosh)如下 2022/7/2528第五章 图像复原与重建共六十二页约束条件为使下式为最小。 令Q=C,复原(f

14、yun)就是在约束条件 下使 为最小。最佳复原(f yun)解为 2022/7/2529第五章 图像复原与重建(zhn jin)共六十二页(3)均方误差最小滤波（维纳滤波）将f和n视为随机变量，并选择Q为噪声与信号之比其中 和 ，分别(fnbi)为信号和噪声的协方差矩阵。则解为一般把1时为含参维纳滤波，=1时为标准维纳滤波。2022/7/2530第五章 图像复原与重建(zhn jin)共六十二页二. 频率域恢复方法 1.逆滤波(lb)恢复法基本原理对于线性移不变系统而言上式两边进行(jnxng)傅立叶变换得 式中G(u,v),F(u,v),H(u,v)和N(u,v)分别是g(x,y), f(x

15、,y), h(x,y) 和n(x,y)的二维傅立叶变换。H(u,v)称为系统的传递函数。从频率域角度看，它使图像退化，因而反映了成像系统的性能。 2022/7/2531第五章 图像复原与重建共六十二页在无噪声(zoshng)的理想情况下，上式可简化为 则：1/H(u,v)称为逆滤波器。对上式再进行傅立叶反变换可得到f(x,y)。但实际上碰到的问题都是有噪声(zoshng)，因而只能求F(u,v)的估计值然后再作傅立叶逆变换得 2022/7/2532第五章 图像复原与重建共六十二页复原步骤： (a)对退化图像g(x，y)作二维离散傅立叶变换，得到G(u,v)； (b)计算系统点扩散函数h(x，y

16、)的二维傅立叶变换，得到H(u,v)；（通常h（x,y）的尺寸小于g(x,y)的尺寸。为了消除(xioch)混叠效应引起的误差，需要把h(x,y)的尺寸延拓。）（c）按 计算（d）计算 的逆傅立叶变换，求得2022/7/2533第五章 图像复原与重建(zhn jin)共六十二页 问题：若噪声为零，则采用(ciyng)逆滤波恢复法能完全再现原图像。若噪声存在，而且H（u,v）很小或为零时，则噪声被放大。这意味着退化图像中小噪声的干扰在H(u,v)较小时，会对逆滤波恢复的图像产生很大的影响，有可能使恢复的图像和f(x,y)相差很大，甚至面目全非。改进的方法：(a)在H（u,v）=0及其附近，人为地

17、仔细设置H-1(u,v)的值，使N(u,v)*H-1(u,v)不会对产生太大影响。2022/7/2534第五章 图像复原与重建(zhn jin)共六十二页下图给出了H(u,v)、H-1(u,v)同改进的滤波特性(txng)HI(u,v)的一维波形，从中可看出与正常的滤波的差别。 (a)图像退化(tuhu)响应 (b)逆滤波器响应 (c)改进的逆滤波器响应 (b)使H(u,v)具有低通滤波性质，即使 2022/7/2535第五章 图像复原与重建共六十二页2.去除由匀速运动引起的模糊原理 在获取图像过程中，由于景物和摄像机之间的相对运动，往往造成图像的模糊。其中由匀速直线运动所造成的模糊图像的恢复

18、问题更具有一般性和普遍意义。因为变速的、非直线的运动在某些条件下可以看成是均匀的、直线运动的合成结果。 设图像f(x,y)有一个平面运动，令x0(t)和y0(t)分别为在x和y方向上运动的变化分量。t表示运动的时间。记录(jl)介质的总曝光量是在快门打开到关闭这段时间的积分。则模糊后的图像g(x,y)为:2022/7/2536第五章 图像复原与重建(zhn jin)共六十二页令G(u,v)为模糊图像(t xin)g(x,y)的傅立叶变换，对上式两边傅立叶变换得 改变上式的积分(jfn)次序，则有： 2022/7/2537第五章 图像复原与重建共六十二页由傅立叶变换(binhun)的位移性质，可

19、得 ：令：则：G(u,v)=H(u,v)F(u,v)-已知退化(tuhu)模型的傅立叶变换式若x(t)、y(t)的性质已知，传递函数可直接由式 求出，因此，f(x,y)可以恢复出来。2022/7/2538第五章 图像复原与重建共六十二页沿水平方向(fngxing)和垂直方向(fngxing)匀速直线运动造成的图像模糊的模型及其恢复的近似表达式：（a)由水平方向均匀直线运动造成的图像模糊的模型及其恢复用以下两式表示： 式中a为总位移量，T为总运动时间(shjin)，m是 的整数部分，L=ka(k为整数)是x的取值范围， A= 2022/7/2539第五章 图像复原与重建共六十二页离散(lsn)式

20、如下： 2022/7/2540第五章 图像复原与重建(zhn jin)共六十二页(b)由垂直方向均匀直线运动造成(zo chn)的图像模糊模型及恢复用以下两式表示： 水平匀速运动(yns yndng)模糊图像的恢复 2022/7/2541第五章 图像复原与重建共六十二页3.维纳滤波复原方法逆滤波(lb)复原方法数学表达式简单，物理意义明确。然而存在着上面讲到的缺点，且难以克服。因此，在逆滤波(lb)理论基础上，不少人从统计学观点出发，设计一类滤波(lb)器用于图像复原，以改善复原图像质量。复原思想：Wienner滤波恢复的思想是在假设图像信号可近似看作平稳随机过程的前提下，按照使恢复的图像与原

21、图像f(x,y)的均方差最小原则来恢复图像。即2022/7/2542第五章 图像复原与重建(zhn jin)共六十二页当采用线性滤波来恢复时，恢复问题(wnt)就归结为找合适的点扩散函数hw(x,y)，使上式成立。由Andrews和Hunt推导满足这一要求的传递函数为： 则有： 这里，H*(u,v)是成像系统传递函数的复共轭；Hw(u,v)就是维纳滤波器的传递函数。Pn（u,v）是噪声功率谱；Pf（u,v）是输入(shr)图像的功率谱。 2022/7/2543第五章 图像复原与重建共六十二页复原过程步骤：(a)计算图像g（x,y）的二维离散傅立叶变换得到 G(u,v)；(b)计算点扩散函数hw

22、(x,y)的二维离散傅立叶变换。同逆滤波一样，为了(wi le)避免混叠效应引起的误差，应将尺寸延拓；(c)估算图像的功率谱密度Pf和噪声的谱密度Pn；(d)由公式 计算图像的估计值 (e)计算 的逆付氏变换，得到恢复后的图像 2022/7/2544第五章 图像复原与重建(zhn jin)共六十二页方法特点：（a）当H（u,v）0或幅值很小时(xiosh)，分母不为零，不会造成严重的运算误差；（b）当Pn 0时，维纳滤波复原方法就是前述的逆虑波复原方法； (c)当Pf 0时， 0，这表示图像无有用信息存在，因而不能从完全是噪音的信号中来“复原”有用信息。对于噪声功率谱Pn(u,v)，可在图像上

23、找一块恒定灰度的区域，然后测定区域灰度图像的功率谱作为Pn(u,v)。 2022/7/2545第五章 图像复原与重建(zhn jin)共六十二页一.图象(t xin)重建的基本概念由物体(wt)截面投影数据，重建该截面的图象图象重建计算机断层扫描CT(Computerized Tomography)技术是一种根据人体截面的投影数据，经过计算机处理来重建该器官截面的图象的技术图象重建的历史英国科学家Housfield在六十年代发明了断层扫描重建图象技术。5.3 图象重建2022/7/2546第五章 图像复原与重建共六十二页图象投影(tuyng)重建基本原理xytsP(t1)P(t2)t1t2f(

24、x,y)设f(x,y)是“截面图象”函数，xOy是建立在物体上的直角坐标系，tOs是建立在扫描装置上的直角坐标系，它们的原点重合，t轴与x轴的夹角为 。两者的坐标变换关系见左边公式。在(t,s)坐标系中，s方向(fngxing)称为射线方向(fngxing)，沿射线方向(fngxing)将图象f(x,y)进行积分P(t)=s f(x,y)ds换句话说， P(t)为f(x,y)在角位置的平行投影。上式也称为雷顿变换(Radon,1917年)2022/7/2547第五章 图像复原与重建共六十二页为了获得(hud)P(t)平行线投影数据，设投影射线为x射线，当x射线穿过物体时，由于物体内部组织的吸收

25、、散射、反射等作用，会使射线强度衰减，射线的衰减规律可表达为硬件(yn jin)实现方法S为射线方向，u(x,y)为射线穿过物体时的衰减系数，它正好反映了物体内部的组织性质，这样，u(x,y)表征了物体断面的图象f(x,y)的组织信息，Nin是射线入射剂量，Nout是穿过物体后射线的剂量，他们都可以通过物理测量而得到。上式可变换为2022/7/2548第五章 图像复原与重建共六十二页如此，投影重建(zhn jin)问题就转化成了如下两个简单问题：1)如何(rh)取Nin,Nout2)如何解方程二.图象重建方法从历史演变来看，通常图象重建的方法有如下5种1）解联立方程组方法2）反投影法3）Fou

26、rier变换法4）卷积-逆滤波法5）逐次逼近法2022/7/2549第五章 图像复原与重建共六十二页解联立方程(lin l fn chn)组法设的步进长度为，t的步进长度为t，令P(i,j)=P(i ,jt)，则经过扫描投影(tuyng)后，可得矩阵例对图象1.jpg进行Radon变换，然后再反变换，观察中间数据。2022/7/2550第五章 图像复原与重建共六十二页2022/7/2551第五章 图像复原与重建(zhn jin)共六十二页ABCDEF6151291290o60o120o612915912解联立方程(lin l fn chn)组得2022/7/2552第五章 图像复原与重建(zh

27、n jin)共六十二页把以上方法推广到多像素多扫描线的情况，便可以(ky)得到重建图象的一般解联立方程组方法。傅立叶变换(binhun)法投影切片定理：设f(x,y)的傅立叶变换为F(u,v)，令u=cos,v=sin, F(,)=F(cos,sin) ，则对以上公式做进一步的推导，可得2022/7/2553第五章 图像复原与重建共六十二页令则滤波(lb)-逆投影法2022/7/2554第五章 图像复原与重建(zhn jin)共六十二页其中(qzhng)称为图象重建滤波器，d是一个足够(zgu)大的实数。三.典型应用CT, MRI, B超等2022/7/2555第五章 图像复原与重建共六十二页

28、四. 三维形状的复原 以上所述的计算机断层，是从多个一维投影数据重建二维断面的方法。为了测出三维物体的形状，可以一面一点点地移动位置，一面求出多个垂直于通过物体中心线的断面，然后把它们依次连接起来，也就是根据一系列二维图像的位置变化构成三维图像。一旦这样的物体三维信息被恢复，就可以求出关于具有任意倾斜度平面的断面，或者可以由三维的任意方向来看物体，从而使对物体形状的判读变得非常容易。使用像计算机X线断层那样的医用图像，由各种各样的视点来观看人体器官的三维图像，对疾病(jbng)的诊断特别重要，为了这一目的的专用装置也正在开发中。2022/7/2556第五章 图像复原与重建(zhn jin)共六十二页从多个断面恢复三维形状(xngzhun)的方法有：2022/7/2557第五章 图像复原与重建(zhn jin)共六十二页1. Voxel 法（体素法） 如果把断面的间隔加密，让断面内的抽样间隔和断面间隔相等，断面内的各像素就可以看成三维空间的小的立方体（voxel）。因此，在多个断面图像(t xin)中，给予相当于这个立方体高度的厚度，把它们三维地堆积起来就可以表现物体的三维图像(t xin)。最近，