图像处理课件

1、2022年5月3日1模拟图像和数字图像模拟图像和数字图像 1)模拟图像模拟图像 模拟图像可用连续函数来描述。模拟图像可用连续函数来描述。 其特点：光照位置和光照强度均为连续变化的。其特点：光照位置和光照强度均为连续变化的。 2)数字图像数字图像 可用矩阵或数组来描述可用矩阵或数组来描述 像素或像元的属性：空间位置和灰度。像素或像元的属性：空间位置和灰度。),(yxFI 1, 11 , 10 , 11, 11 , 10 , 11, 01 , 00 , 0,NMMMNNiiiiiiiiinmII2022年5月3日3数字图像处理空域法 把图像看作是平面中各个像素组成的集合，把图像看作是平面中各个像素

2、组成的集合，直接对这一二维函数进行相应的处理。直接对这一二维函数进行相应的处理。 分类：分类： 邻域处理法：梯度运算邻域处理法：梯度运算 (gradient (gradient Algorithm)Algorithm)、拉普拉斯算子运算、拉普拉斯算子运算 (Laplacian (Laplacian operator)operator)、平滑算子运算、平滑算子运算 (Smoothing (Smoothing operator)operator)卷积运算卷积运算 (Convolution algorithm)(Convolution algorithm) 点处理法：点处理法：灰度处理（灰度处理（g

3、rey processing），面），面积、周长、体积、重心运算等等。积、周长、体积、重心运算等等。2022年5月3日4数字图像处理变换域法 首先对图像进行正交变换，得到变换域系对图像进行正交变换，得到变换域系数阵列，然后再施行各种处理，处理后再数阵列，然后再施行各种处理，处理后再反变换到空间域，得到处理结果。反变换到空间域，得到处理结果。 包括：滤波、数据压缩、特征提取等处理。包括：滤波、数据压缩、特征提取等处理。2.1.1 图像的表示 数学表示 二维离散亮度函数f(x,y) x,y说明图像像素的坐标 函数值f代表(x,y)处像素的亮度值 二维矩阵Am,n m,n为图像的宽和高 矩阵元素a(

4、I,j)的值，表示图像在第i行，第j列的像素的灰度值。灰度级：表示像素明暗程度的整数量。层次：表示灰度级的数量。图像数据的实际层次越多视觉效果就越好。2.1.2 图像的质量：层次2.1.2 图像的质量：对比度 对比度：指一幅图像中灰度反差的大小。 对比度=最大亮度/最小亮度。图像的颜色：RGB模型 扬-赫姆霍尔兹视觉三基色假说 C=R+G+B 视网膜椎体细胞感红、感绿、感蓝色素 相加混色：相加混色： 红绿黄红绿黄 红蓝紫红蓝紫 蓝绿青蓝绿青 红绿红绿蓝白蓝白图像的颜色： CMYK CMYK减色法 青（Cyan）、洋红（Magenta）、黄（Yellow）、黑（Black） 2.1.4 图像的像

5、素：邻域象素的邻域象素的邻域4-邻域N4(p)： 对角邻域ND(p)：8-邻域N8(p)：prrsssrsrprrrrpssss2.1.4 图像的像素：连通像素的连通性 4-连通像素的连通性像素的连通性像素的连通性为空象素间的距离距离量度函数距离量度函数 (1) 欧氏（Euclidean）距离 (2) 城区（city-block）距离 (3) 棋盘（chessboard）距离2/1 22E)()(),(tysxqpD ),(4tysxqpD) , ( max),(8tysxqpD2.2.1图像的数字化：统一的采样和量化 图像采样的形式化定义 设Z表示整数集合 采样处理：将xy平面分配到一个网格

6、上，且每一个网格中心的坐标是一个笛卡儿乘积ZxZ的元素对，即所有有序元素对(a,b)的集合，其中a和b属于整数集合Zxy平面(a,b)2.2.1图像的数字化：统一的采样和量化 图像的量化 取值的数字化被称为图像灰度级量化 量化处理：将f 映射到Z的处理 Z的最大取值，确定像素的灰度级数G = 2m， 如256f19第三章 数字图像处理基础 3.1 图像运算算术运算逻辑运算 3.2 空域变换几何变换非几何变换 3.3 频域变换傅立叶变换导言傅立叶变换的特性快速傅立叶变换第三章 数字图象处理基础203.1.1 图象运算：算术运算 去除“叠加性”噪音 对于原图象f(x,y),有一个噪音图象集 gi(

7、x,y) i =1,2,.M其中：gi(x,y) = f(x,y) + h(x,y)iM个图象的均值定义为：g(x,y) = 1/M (g0(x,y)+g1(x,y)+ gM(x,y)当：噪音h(x,y)i为互不相关，且均值为0时，上述图象均值将降低噪音的影响。第三章数字图象处理基础 第一节图象运算213.2.2 非几何变换 3.2.2 非几何变换非几何变换的定义模板运算灰度级变换直方图223.2.2 非几何变换：非几何变换的定义 非几何变换的定义对于原图象f(x,y)，灰度值变换函数T(f(x,y)唯一确定了非几何变换： g(x,y) = T(f(x,y) （2）g(x,y)是目标图象没有几

8、何位置的改变。2022年5月3日23/36 目的 改善图像的视觉效果，提高图像的清晰度 将图像转换成一种更适合于人或机器分析处理的形式 在图像增强的过程中，没有新信息的增加，只是通过压制一部分信息，突出另一部分信息。有选择地突出其中感兴趣的某些特有选择地突出其中感兴趣的某些特征，衰减其中不需要的特征征，衰减其中不需要的特征第四章第四章 图像增强图像增强2022年5月3日24/36 方法 空域处理（空间滤波）空域处理（空间滤波） 点运算：直接对原图像各像素的灰度值进行逐点运算；代数运算（算术运算）：一般整个图像范围内运算。 模板运算 (局部运算)：在像点的邻域范围内运算； 频域处理（频域滤波）借

9、助傅氏变换，增强感兴趣的频域处理（频域滤波）借助傅氏变换，增强感兴趣的频率分量频率分量 具体操作 去除噪声、边缘增强、提高对比度、增加亮度改善颜去除噪声、边缘增强、提高对比度、增加亮度改善颜色效果、改善细微层次等色效果、改善细微层次等 算法针对性强：理论上缺乏统一性，图像质量评算法针对性强：理论上缺乏统一性，图像质量评价：主观判断价：主观判断引言引言2022年5月3日25/36 点运算（点运算（point operation）定义）定义 对于一幅输入图像，将产生一幅输出图像，输出对于一幅输入图像，将产生一幅输出图像，输出图像的每个像素点的灰度值由输入像素点决定。图像的每个像素点的灰度值由输入像

10、素点决定。点运算由灰度变换函数（点运算由灰度变换函数（gray-scale transformation, GST）确定。）确定。备注：备注： （1）与局部（邻域）运算的差别，输入像素）与局部（邻域）运算的差别，输入像素-输出像素一一对应；输出像素一一对应； （1）与几何运算的差别，不改变图像的空间关系；）与几何运算的差别，不改变图像的空间关系； （2）又称为对比度增强，对比度拉伸或灰度变换。）又称为对比度增强，对比度拉伸或灰度变换。( , ),B x yfA x y点运算点运算2022年5月3日26/362022年5月3日27/362022年5月3日28/362022年5月3日294. 直方

11、图均衡的实现直方图均衡的实现 直方图均衡的步骤：直方图均衡的步骤： （1 1）计算原图像的归一化灰度级别及其分布概）计算原图像的归一化灰度级别及其分布概率率p pr r(r(rk k)=n)=nk k/n/n。 （2 2）根据直方图均衡化公式（）根据直方图均衡化公式（4.134.13）求变换函）求变换函数的各灰度等级值数的各灰度等级值s sk k。2022年5月3日304. 直方图均衡的实现直方图均衡的实现( (续续1)1) 直方图均衡的步骤：直方图均衡的步骤： （3 3）将所得的变换函数的各灰度等级值转化成标准）将所得的变换函数的各灰度等级值转化成标准的灰度级别值。的灰度级别值。即把第（2）

12、步求得的各sk值，按靠近原则近似到与原图像灰度级别相同的标准灰度级别中。此时获得的即是均衡化后的新图像中存在的灰度级别值，其对应的像素个数不为零；对于那些在变换过程中“被丢失了的”灰度级别值，将其像素个数设为零。2022年5月3日314. 直方图均衡的实现直方图均衡的实现（续（续2 2） （4 4）求新图像的各灰度级别值）求新图像的各灰度级别值s sl l(l=0,1,(l=0,1,L-1),L-1)的像数数目。的像数数目。在前一步的计算结果中，如果不存在灰度级别值sl，则该灰度级别的像素数目为零；如果存在灰度级别值sl，则根据其与之相关的sk=T(rk)和sk的对应关系，确定该灰度级别sl的

13、像数数目。 （5 5）用）用s sk k代替代替s sl l(k,l=0,1,(k,l=0,1,L-1),L-1)，并进而求新，并进而求新图像中各灰度级别的分布概率图像中各灰度级别的分布概率p ps s(s(sk k)=m)=mk k/n/n。 （6 6）画出经均衡化后的新图像的直方图，）画出经均衡化后的新图像的直方图，2022年5月3日32例例4.2.14.2.1 已知有一幅大小为已知有一幅大小为64646464的图像，灰度级的图像，灰度级为为8 8。图像中各灰度级的像素数目如表。图像中各灰度级的像素数目如表4-14-1所示。要所示。要求：求： (a) (a) 画出原图像的直方图；画出原图像

14、的直方图；表表4.1 6464图像灰度分布图像灰度分布 (b) (b) 利用直方图均衡方法利用直方图均衡方法求出均衡化后新图像的直方图。求出均衡化后新图像的直方图。 knk0 07907901 1102310232 28508503 36566564 43293295 52452456 61221227 781812022年5月3日33解：解：（1 1）画原图像的直方图）画原图像的直方图 归一化灰度级，即求归一化灰度级，即求r rk k=k/(L-1)=k/7=k/(L-1)=k/7，结果如，结果如表表4.24.2所示所示。 表表4.2 归一化灰度分布及概率归一化灰度分布及概率 计算第计算第k

15、 k个灰度级出现个灰度级出现的概率的概率p pr r(r(rk k)=n)=nk k/n=n/n=nk k/4096/4096，结果如表结果如表4.24.2所示。所示。 所画的原图像的直方所画的原图像的直方图如图图如图4.94.9所示。所示。 例例4.2.1 4.2.1 ( (续续1)1)k kr rk kPr(rPr(rk k)=n)=nk k/n/n0 0=0=00.190.191 1=1/7=1/70.250.252 2=2/7=2/70.210.213 3=3/7=3/70.160.164 4=4/7=4/70.080.085 5=5/7=5/70.060.066 6=6/7=6/70

16、.030.037 7=1=10.020.022022年5月3日34例例4.2.1 4.2.1 (续续2) 0.250.25 0.20 0.20 0.15 0.15 0.10 0.10 0.05 0.05 0 0)(krrpkr7172737475761图图4.92022年5月3日35解：解：（2 2）利用直方图均衡化方法求出均衡化后的新）利用直方图均衡化方法求出均衡化后的新图像的直方图。图像的直方图。 根据直方图均衡化公式（根据直方图均衡化公式（4.134.13）求变换函数的）求变换函数的各灰度等级值。各灰度等级值。 例例4.2.1 4.2.1 (续续3)19. 04096790)(0000j

17、jnnrTs44. 0409610234096790)(1011jjnnrTs65. 04096850409610234096790)(2022jjnnrTs同理有：同理有： 0 . 1;98. 0;95. 0;89. 0;81. 076543sssss2022年5月3日36解：解：对应的变换函数如图对应的变换函数如图4.104.10所示。所示。 例例4.2.1 4.2.1 (续续4) 1.001.00 0.80 0.80 0.60 0.60 0.40 0.40 0.20 0.20 0 0kskr71727374757612022年5月3日37解：解： 将所得的变换函数的各灰度等级值转化成标将

18、所得的变换函数的各灰度等级值转化成标准的灰度级别值准的灰度级别值 先将先将s sk k值按靠近原则对应到原灰度级别中值按靠近原则对应到原灰度级别中: :分数值：分数值： 0 1/7 2/7 3/7 4/7 5/7 6/7 1十进制值十进制值： 0 0.143 0.286 0.429 0.571 0.714 0.857 1 比较可得：比较可得：例例4.2.1 4.2.1 (续续5)1;76;75;73;7176543210ssssssss2022年5月3日38解：解： 求新图像的各灰度级别值求新图像的各灰度级别值s sl l（l=0,1,l=0,1,7,7）的像数数目的像数数目 用用s sk k

19、代替代替s sl l（k,l=0,1,k,l=0,1,7,7），并求新图像中），并求新图像中各灰度级别的概率各灰度级别的概率p ps s(s(sk k)=m)=mk k/n=m/n=mk k/4096/4096，结果如表所，结果如表所示。示。 例例4.2.1 4.2.1 (续续6)2022年5月3日39例例4.2.1 4.2.1 (续续10)kkskmnmspkks)(0 00 00 00 01 11/71/77907900.190.192 22/72/70 00 03 33/73/7102310230.250.254 44/74/70 00 05 55/75/78508500.210.216

20、 66/76/79859850.240.247 71 14484480.110.11解：解：2022年5月3日40例例4.2.1 4.2.1 (续续11)解：解： 画出经均衡化后的新图像的直方图，如图画出经均衡化后的新图像的直方图，如图4.114.11所示。所示。 0.250.250.200.200.150.150.100.100.050.05 0 07172737475761)(ksspks图图4.11 4.11 新图像的直方图新图像的直方图 2022年5月3日411. 问题的提出问题的提出背景：直方图均衡化的缺陷：不能用于交互方式的图像增强应用，因为直方图均衡化只能产生唯一一个结果，恒定值

21、直方图近似希望通过一个指定的函数（如高斯函数）或用交互图形产生一个特定的直方图。根据这个直方图确定一个灰度级变换T(r),使由T产生的新图像的直方图符合指定的直方图2022年5月3日42/121空域图像增强:空域过滤器 点运算增强 直方图增强空域过滤器1) 空域过滤器的基本概念 2) 钝化过滤器 3) 锐化过滤器2022年5月3日43/1211) 空域过滤处理的基本概念空域过滤及过滤器的定义 使用空域模板进行的图像处理，被称为空域过滤。模板本身被称为空域过滤器2022年5月3日44/121 模板运算(邻域操作) 去除或衰减图像中噪声和假轮廓 邻域平均法 中值滤波 突出图像的边缘与特征 一阶微分

22、 二阶微分空间域增强空间域增强2022年5月3日45/121 通常用一种称为“模板卷积”的运算来实现线性邻域增强。 其步骤为： 将模板在图中漫游，并将模板中心与图中某个象素位置重合； 将模板上系数与模板下对应象素相乘； 将所有乘积相加； 将和赋给图中对应模板中心位置的象素。模板卷积模板卷积2022年5月3日46/121 主观上使图像变得柔和 过于突出的图像边缘细节及轮廓使得图像缺乏柔和性。 消除图像噪声 图像噪声由于具有随机性，一般体现为图像的高频成分。平滑滤波器可以抑制图像噪声。邻域平滑邻域平滑2022年5月3日47/121邻域平滑邻域平滑2022年5月3日48/121邻域平滑邻域平滑202

23、2年5月3日49/121 滤波模板的选择滤波模板的选择 模板尺寸越大，滤波程度越高；噪声消除越彻底，而同时，图像细节的损失也越大。 通常模板的所有系数之和为1，这是为了保证模板卷积完成后目标图像的总体亮度不发生改变。 模板系数的值可正可负，可大可小。系数值越大，代表该系数对应象素的权重越高。负值则可用于锐化滤波器。邻域增强邻域增强2022年5月3日50/121 非线性滤波器 最大值滤波 主要用途：寻找最亮点 计算公式： g(x,y)= max (Nf(x,y) 最小值滤波 主要用途：寻找最暗点 计算公式： g(x,y)= min (Nf(x,y) 中值滤波器 基本思想是用象素点邻域灰度值的中值

24、来代替该象素点的灰度值， 主要用途：钝化图像、去除噪音。 计算公式： g(x,y)= mid(Nf(x,y)邻域增强邻域增强2022年5月3日51/121 中值滤波算法的实现 计算以点i,j为中心的函数窗象素中值步骤如下： 将模板区域内的象素排序， 选择排序象素集的中间值作为点i,j的新值。2022年5月3日52/121因为噪声（如椒盐噪声）的出现，使该点像素比因为噪声（如椒盐噪声）的出现，使该点像素比周围的像素亮（暗）许多。周围的像素亮（暗）许多。如果在某个模板中，对像素进行由小到大排列的如果在某个模板中，对像素进行由小到大排列的重新排列，那么最亮的或者是最暗的点一定被排重新排列，那么最亮的

25、或者是最暗的点一定被排在两侧。在两侧。取模板中排在中间位置上的像素的灰度值替代待取模板中排在中间位置上的像素的灰度值替代待处理像素的值，就可以达到滤除噪声的目的。处理像素的值，就可以达到滤除噪声的目的。2022年5月3日53/121锐化锐化微分能产生与钝化相反的效果，即锐化的效果。2022年5月3日54/121 图像锐化图像锐化 1 ）微分法2022年5月3日55/121 在图像处理中应用微分最常用的方法是计算梯度。f = f / x , f / y 函数f(x,y)在(x,y)处的梯度为一个向量：梯度向量的计算：f = mag(f ) = (f / x)2 +(f / y)21/2 图像锐化

26、的方法图像锐化的方法梯度法梯度法2022年5月3日56/121 计算这个向量的大小为： f = mag(f ) = (f / x)2 +(f / y)21/2 考虑一个3x3的图像区域，z代表灰度级，上式在点z5的f值可用数字方式近似。 (f / x) 用（z5 z6）近似 (f / y) 用（z5 z8）近似，组合为： f (z5 - z6)2 + (z5 - z8)21/2z2z8z5z3z9z6z1z7z42022年5月3日57/121用绝对值替换平方和平方根 有： f |z5 - z6| + |z5 - z8|另外一种计算方法是使用交叉差： f (z5 - z9)2 + (z6 - z

27、8)21/2 f |z5 - z9| + |z6 - z8|z2z8z5z3z9z6z1z7z42022年5月3日58/121水平垂直差分方法和罗伯特差分方法的图示形式为：水平垂直差分方法和罗伯特差分方法的图示形式为： ), ( jif) 1, (jif), 1(jif ) 1, 1( jif ), ( jif) 1, (jif), 1(jif ) 1, 1( jif(a) (a) 水平垂直差分水平垂直差分 (b) (b) 罗伯特差分罗伯特差分 2022年5月3日59/121微分过滤器模板 Roberts交叉梯度算子 Prewitt梯度算子 Sobel梯度算子2022年5月3日60/121 梯

28、度计算由两个模板组成，第一个求 得梯度的第一项，第二个求得梯度的 第二项，然后求和，得到梯度。z2z8z5z3z9z6z1z7z401-10-1001 Roberts 交叉梯度算子 f |z5-z9| + |z6-z8|两个模板称为Roberts 交叉梯度算子。2022年5月3日61/121 f |(z7 +z8 + z9) - (z1 + z2 + z3) | + |(z3 +z6 + z9) - (z1 + z4 + z7) |z2z8z5z3z9z6z1z7z4-110-110-110000-1-1-1111 Prewitt梯度算子3x3的梯度模板2022年5月3日62/121 f |(

29、z7 +2z8 + z9) - (z1 + 2z2 + z3) | + |(z3 +2z6 + z9) - (z1 + 2z4 + z7) |z2z8z5z3z9z6z1z7z4-220-110-110000-1-1-2112 Sobel梯度算子3x3的梯度模板 2022年5月3日63/121锐化高通滤波过滤器效果的分析 常数或变化平缓的区域，结果为0或很小，图像很暗，亮度被降低了 在暗的背景上边缘被增强了 图像的整体对比度降低了 计算时会出现负值，归0处理为常见2022年5月3日64/121锐化高通滤波基本高通空域滤波的缺点和问题 高通滤波在增强了边的同时，丢失了图像的层次和亮度2022年5

30、月3日65/121锐化高通滤波(3) 高频补偿过滤高频补偿过滤的原理过滤器扩大因子及模板系数的设计高频补偿过滤模板尺寸的选定高频补偿过滤器效果的分析2022年5月3日66/121高频补偿过滤的原理 弥补高通滤波的缺陷，在增强边和细节的同时，不丢失原图像的低频成分。高通滤波可看作为：高通 = 原图 低通在上式原图上乘一个扩大因子A，有高频补偿过滤：高频补偿 = A原图 低通2022年5月3日67/121高频补偿过滤的原理高频补偿 = A原图 低通 = (A 1)原图 + （原图 低通） = (A 1)原图 + 高通 当A = 1时，高频补偿就是高通过滤， 当A 1 时，原图像的一部分被加到高通中

31、。2022年5月3日68/121高频补偿过滤的原理高频补偿 = (A 1) * 原图 + 高通2022年5月3日69/121过滤器扩大因子及模板系数设计 对于 3x3的模板，设 w = 9A 1；（高通时 w = 8）A的值决定了过滤器的特性 当 A = 1.1时，意味着把 0.1个原图像加到基本高通上。-1-1w-1-1-1-1-1-11/9 *高频补偿模板70/72频域图像增强引言 频域增强的原理 频率平面与图像空域特性的关系 图像变化平缓的部分靠近频率平面的圆心，这个区域为低频区域 图像中的边、噪音、变化陡峻的部分，以放射方向离开频率平面的圆心，这个区域为高频区域2022年5月3日71/

32、1214.54.5频率域图像增强频率域图像增强 4.5.1 4.5.1 基本实现思想和实现方法基本实现思想和实现方法 2022年5月3日72/1214.5.14.5.1基本实现思想和实现方法基本实现思想和实现方法 设设f(x,y)f(x,y)为输入图像，为输入图像，F(u,v)F(u,v)为输入图像的傅为输入图像的傅立叶变换，立叶变换，H(u,v)H(u,v)为转移函数为转移函数( (也称为滤波函数也称为滤波函数) )， G(u,v)G(u,v)为对为对F(u,v)F(u,v)进行频率域滤波后的输出，进行频率域滤波后的输出， g(x,y)g(x,y)为经频率域滤波后的输出图像，则有：为经频率域

33、滤波后的输出图像，则有： ),(),),(vuHvuFvuG（),(),(1vuGFyxg（4.56） （4.57） 2022年5月3日73/1214.5.14.5.1基本实现思想和实现方法基本实现思想和实现方法 g(x,y)g(x,y)增强后的图像增强后的图像f(x,y)f(x,y)输入图像输入图像F(u,v)H(u,v)F(u,v)傅立叶傅立叶变变 换换频率滤波频率滤波H(u,v)傅立叶傅立叶反变换反变换预处理预处理后处理后处理图图4.27 频率域图像增强步骤频率域图像增强步骤 第五章第五章 图像复原图像复原 概述概述 图像退化图像退化/ /复原过程的模型复原过程的模型 噪声模型噪声模型

34、空间域滤波复原（唯一退化是噪声）空间域滤波复原（唯一退化是噪声） 频率域滤波复原（削减周期噪声）频率域滤波复原（削减周期噪声） 逆滤波逆滤波 维纳滤波维纳滤波 几何变换几何变换 退化退化 成像过程中的成像过程中的”退化退化”，是指由于成像系，是指由于成像系统各种因素的影响，使得图像质量降低。统各种因素的影响，使得图像质量降低。2. 图像复原：图像复原：利用图像本身的有关知识去除或减轻图像品质利用图像本身的有关知识去除或减轻图像品质下降（退化）的处理方法。下降（退化）的处理方法。5.1 概述概述 改善图像质量的两种方法： 图像的增强图像的增强 不考虑图像降质的原因，只将图像中感兴趣的特征有选择的

35、突出，而衰减次要信息； 特点：能提高图像的可读性，但改善后的图像不一定逼近原始图像， 如衰减各种噪声、突出目标的轮廓等。 图像的复原图像的复原 针对图像降质的原因，设法去补偿降质因素，从而使改善后的图像尽可能的逼近原始图像。5.1 概述概述图像增强：主观过程，而图像复原：客观过程图像增强：主观过程，而图像复原：客观过程。 图像增强被认为是一种对比度拉伸等，提供给图像增强被认为是一种对比度拉伸等，提供给用户喜欢接收的图像；而图像复原技术追求恢用户喜欢接收的图像；而图像复原技术追求恢复原始图像的最优估值。复原始图像的最优估值。图像复原技术可以使用空间域或频率域滤波器图像复原技术可以使用空间域或频率

36、域滤波器实现。实现。5.1 5.1 概述概述5.2 图像退化图像退化/复原过程的模型复原过程的模型 f(x,y)表示一幅输入图像表示一幅输入图像 g(x,y)是是f(x,y)产生的一幅退化图像产生的一幅退化图像 H表示退化函数表示退化函数 (x,y )表示外加噪声表示外加噪声给定给定g(x,y)，关于退化函数，关于退化函数H的一些知识和外加噪声项的一些知识和外加噪声项(x,y ), 怎样获得关于原始图像的近似估计怎样获得关于原始图像的近似估计 ？),(yxf5.4 空间域滤波复原（唯一退化是噪空间域滤波复原（唯一退化是噪声）声）5.4 空间域滤波复原（唯一退化是噪空间域滤波复原（唯一退化是噪声

37、）声） 图像复原的空间滤波器图像复原的空间滤波器 均值滤波器均值滤波器 算术均值滤波器、几何均值滤波器、谐波均值滤算术均值滤波器、几何均值滤波器、谐波均值滤波器波器 顺序统计滤波器顺序统计滤波器 中值滤波器、最大值滤波器、最小值滤波器、中中值滤波器、最大值滤波器、最小值滤波器、中点滤波器、修正后的阿尔法均值滤波器点滤波器、修正后的阿尔法均值滤波器 自适应滤波器自适应滤波器 自适应局部噪声消除滤波器、自适应中值滤波器自适应局部噪声消除滤波器、自适应中值滤波器5.4.1 均值滤波器均值滤波器1. 算术均值滤波器算术均值滤波器5.4.1 均值滤波器均值滤波器2. 几何均值滤波器几何均值滤波器5.4.

38、2 顺序统计滤波器顺序统计滤波器1. 中值滤波器中值滤波器5.4.2 顺序统计滤波器顺序统计滤波器5. 修正后的阿尔法均值滤波器修正后的阿尔法均值滤波器mn-1，5.4.3 自适应滤波器自适应滤波器2. 自适应中值滤波器自适应中值滤波器定义下列符号：定义下列符号：5.4.3 自适应滤波器自适应滤波器算法算法：5.4.3 自适应滤波器自适应滤波器5.4.3 自适应滤波器自适应滤波器逆滤波没有说明怎样处理噪声逆滤波没有说明怎样处理噪声. 维纳滤波综合考虑退化函数和噪声统计特征维纳滤波综合考虑退化函数和噪声统计特征. ,f认为图像和噪声是随机函数目标是找未污染图像 的估计值f,使它们之间的均方误差最

39、小.22() eEff(5.8.1)221|( , )|( , )( , )( , ) |( , )|( , )/( , )fH u vF u vG u vH u vH u vS u vSu v(5.8.2)式式(5.8.1)中误差函数的最小值在频率中用下式表达中误差函数的最小值在频率中用下式表达:维纳滤波维纳滤波, 括号中的项组成的滤波器通常称为最小均方误差滤波器括号中的项组成的滤波器通常称为最小均方误差滤波器,或最小二乘方或最小二乘方误差滤波器误差滤波器.处理白噪声处理白噪声(噪声的傅立叶谱为常量噪声的傅立叶谱为常量)时时,谱谱|N(u,v)|2是一个常是一个常数数,问题可以简化问题可以简

40、化,但但|F(u,v)|2未知未知.*2*22( , )( , )( , )|( , )|( , )( , )( , ) |( , )|( , ) |( , )|fH u vHu vH u vH u vHu v H u vS u vN u vSu vF u v退化函数的复共轭噪声的功率谱未退化图像的功率谱第5章图像复原 5.8 维纳滤波维纳滤波221|( , )|( , )( , )( , ) |( , )|H u vF u vG u vH u vH u vKK为特殊常数为特殊常数.经常用下式近似经常用下式近似:(5.8.2)的维纳滤波要求的维纳滤波要求: 未退化图像和噪声的功率必须是已知的未

41、退化图像和噪声的功率必须是已知的.虽然用虽然用(5.8.3)近似的方法能得到好的结果近似的方法能得到好的结果,但功率谱比的常数但功率谱比的常数K的估的估计一般没有合适的解计一般没有合适的解.(5.8.3)第5章图像复原 5.8 维纳滤波维纳滤波第六章 图像分割 6.1 图像分割引言 6.2 边界分割法 6.3 边缘连接分割法 6.4 阈值分割法 6.5 面向区域的分割 6.6 数学形态学图像处理6.2 图像分割:边界分割法 点的检测 用空域的高通滤波器来检测孤立点例： R = (-1 * 8 * 8 + 128 * 8) / 9 = (120 * 8) / 9 = 960 / 9 = 106设

42、 ：阈值：T = 64 R T8888128 8888图像-1-1-1-18-1-1-1-1模板6.2 图像分割:边界分割法 点的检测算法描述 设定阈值 T，如T = 32、64、128等,并计算高通滤波值R 如果R值等于0，说明当前检测点的灰度值与周围点的相同 当R的值足够大时，说明该点的值与周围的点非常不同，是孤立点。通过阈值T来判断 |R| T 检测到一个孤立点6.2 图像分割:边界分割法 线的检测 通过比较典型模板的计算值，确定一个点是否在某个方向的线上-1-1-1222-1-1-1水平模板-1-12-12-12-1-145度模板-12-1-12-1-12-1垂直模板2-1-1-12-1-1-12135度模板6.2 图像分割:边界分割法 线的检测算法描述 依次计算4个方向的典型检测模板，得到Ri i=1,2,3,4 如 |Ri| |Rj| 对于所有的j = i，那么这个点被称为在方向上更接近模板i 所