数字图像处理第6章

合集下载

第6章彩色图像处理资料

第6章彩色图像处理资料
V=max(红色、蓝色、绿色);
补充 YUV彩色空间
YUV是被欧洲电视系统所采用的一种颜色编 码方法(属于PAL) 。
Y为颜色的亮度 U 为色差信号,为红色的浓度偏移量成份 V 为色差信号,为蓝色的浓度偏移量成份 YUV格式有:4∶4∶4 ;4∶2∶2 ;
4∶1∶1 ;4∶2∶0
YUV与RGB间的转换
6.1 彩色基础 p252
将红、绿、蓝的量称为三色值,表示为X,Y,Z, 则一种颜色由三色值系数定义为:
x X X Y Z
y Y X Y Z
z Z X Y Z
x y z 1
CIE色度图
纯色在色度图边 界上,任何不在 边界上而在色度 图内的点都表示 谱色的混合色;
越靠近等能量点 饱和度越低,等 能量点的饱和度 为0;
Y 0.299 0.587 0.114R
U
0.147
0.289
0.436 G
V 0.615 0.515 0.1 B
R 1 0
1.1398 Y
G 1
0.3946
Hale Waihona Puke 0.5805UB 1 2.032 0.0005V
6.3 伪彩色图像处理
伪彩色(又称假彩色)图像处理是根据特定的 准则对灰度值赋以彩色的处理,即将灰度 图转换为彩色图。
6.2.2 CMY和CMYK模型
CMY模型和RGB模型间的关系:
C 1 R
M
1
G
Y 1 B
RGB三个值已归一化为[0,1]
等量的青色、品红和黄色应该产生黑色。但实 际产生的黑色不够纯正,另外加上价格因素, 引入黑色(打印的主色),构成CMYK模型。
6.2.2 CMY和CMYK模型

精品课件-HALCON数字图像处理-第6章 图像增强

精品课件-HALCON数字图像处理-第6章  图像增强
◘图像增强 目的:一是改善图像的视觉效果,提高图像的清晰度; 二是将图像转换成一种更适合于人或机器分析处理的形
式。
分类: 频域法:直接对图像的像素灰度值进行操作。包括图像的
灰度变换、直方图修正、图像空域平滑和锐化处理、彩色增强 等。
空域法:在图像的变换域中,对图像的变换值进行操作, 然后经逆变换获得所需的增强结果。常用的方法包括低通滤 波 、高通滤波以及同态滤波等。
二、为什么要增强图象? 图像在传输或者处理过程中会引入噪声或使图像变模 了 困难。
Digital Image
6.0 概 述
三、目的: 1.改善图像的视觉效果,提高图像的清晰度; 2.将图像转换成一种更适合于人或机器分析处理的形式。
注意:在图像增强的过程中,没有新信息的增加, 只是 通过压制一部分信息,从而突出另一部分信息。
( f (m, n)
)
其中λ 和γ 为常数。为避免时底数为0的情况,增加偏移量 ε 。γ 值的选择对于变换函数的特性有很大影响,当γ < 1 时会将原图像的灰度向高亮度部分映射,当γ >1时向低亮度部 分映射,而当γ =1时相当于正比变换。灰度指数变换的图 像 示例如图4.1.5所示。
Digital Image
)=
nnj
0.19
0.25
0.21 0.24 0.11
6.2 图像的直方图修正
图4.2.3给出了直方图均衡化的示意图。从图和表中可以看 出,由于 数字图像灰度取值的离散性,通过四舍五入使变换后 的分灰布度,值但。出相现比了于归原直并现方象图,要而平使坦变得换多后的直方图并非完全均匀
图4.2.3 直方图均衡化的示意图
第6章
图像增强
◆ 6.0 ◆ 6.1 ◆ 6.2 ◆ 6.3 ◆ 6.4 ◆ 6.5 ◆ 6.6

数字图像处理第六章色彩模型与彩色处理课件

数字图像处理第六章色彩模型与彩色处理课件

Chapter 6 Color Image Processing
6.1 彩色基础
在颜料或着色剂中 ,原色的定义是这样 的:
白:减去一种原色 , 反射或传输另两种 原色。故其原色是: 深红、青、黄。而二 次色是R、G、B。如 图6.4所示。
Chapter 6 Color Image Processing
Chapter 6 Color Image Processing
6.2 彩色模型
6.2.1 RGB彩色模型
下面介绍所谓 全RGB彩色子集。
Chapter 6
Color Image Processing
6.2 彩色模型
Chapter 6 Color Image Processing
6.2 彩色模型
6.3 伪彩色处理
6.3 伪彩色处理 给特定的灰度值赋以彩色。伪彩色的 目的是为了人眼观察和解释图像中的目标。
Chapter 6 Color Image Processing
6.3 伪彩色处理
6.3.1 强度分层
参见图6.18,图像被看成三维函数。
Chapter 6 Color Image Processing
6.3.2 灰度级到 彩色转换
例6.5是一突出 装在行李内的爆炸物 的伪彩色应用。
Chapter 6 Color Image Processing
6.3 伪彩色处理
6.3.2 灰度级到彩 色转换
例6.5是一突出装 在行李内的爆炸物的伪彩 色应用。
Chapter 6 Color Image Processing
6.3 伪彩色处理
Chapter 6 Color Image Processing
6.3 伪彩色处理

遥感数字图像处理-第6章 几何校正

遥感数字图像处理-第6章 几何校正
3
二、几何校正原理
几何校正涉及两个过程: ➢ 一是空间位置(像元坐标)的变换 ➢ 二是像元灰度值的重新计算(重采样)
4
二、几何校正原理
坐标转换 (a)直接法;(b)间接法
5
三、几何校正步骤
几何精校正不需要空间位置变化数据,回避了成像的空间 几何过程,主要借助地面控制点实现校正。其主要校正步 骤为:
第6章
几何校正
几何校正
一、几何校正原理 二、几何校正步骤 三、几何校正类型 四、图像匹配 五、投影转换
难点:图像匹配 重点:几何校正方法
2
一、几何校正原理
几何校正和几何配准
➢ 几何配准是指将不同时间、不同波段、不同传感器系统所获得的同一 地区的图像(数据),经几何变换使同名像点在位置上和方位上完全 叠合的操作。
➢ 对畸变图像和基准图像建立统一的坐标系和地图投影。 ➢ 选择地面控制点(GCP),按照GCP选择原则,在畸变图像
和基准图像上寻找相同位置的地面控制点对。 ➢ 选择校正模型,利用选择的GCP数据求取校正模型的参数,
然后利用校正模型实现畸变图像和基准图像之间的像元坐 标变换。 ➢ 选择合适的重采样方法对畸变图像的输出图像像元进行灰 度赋值。 ➢ 几何校正的精度分析。
9
四、图像匹配
3.图像匹配方法 根据图像特征的选择,图像匹配方法一般可以分为基于灰
度的图像匹配和基于特征的图像匹配。
10
➢ 几何配准与几何校正的原理是完全相同的,即都涉及到空间位置(像 元坐标)变换和像元灰度值重采样处理两个过程。
➢ 二者的区别主要在于其侧重点不相同:几何校正注重的是数据 本身的处理,目的是为了对数据的一种真实性还原。而几何配 准注重的是图和图(数据)之间的一种几何关系,其目的是为 了和参考数据达成一致,而不考虑参考数据的坐标是否标准、 是否正确。也就是说几何校正和几何配准最本质的差异在于参 考的标准。另外,几何校正更像前期数据处理,几何配准更像 后期处理。

胡学龙《数字图像处理(第二版)》课后习题解答

胡学龙《数字图像处理(第二版)》课后习题解答

2
1.PHOTOSHOP:当今世界上一流的图像设计与制作工具,其优越性能令其产品望尘 莫及。PHOTOSHOP 已成为出版界中图像处理的专业标准。高版本的 P扫描仪、数码相机等图像输入设备采集的图 像。PHOTOSHOP 支持多图层的工作方式,只是 PHOTOSHOP 的最大特色。使用图层功能 可以很方便地编辑和修改图像,使平面设计充满创意。利用 PHOTOSHOP 还可以方便地对 图像进行各种平面处理、绘制简单的几何图形、对文字进行艺术加工、进行图像格式和颜色 模式的转换、改变图像的尺寸和分辨率、制作网页图像等。
1.5 常见的数字图像处理开发工具有哪些?各有什么特点? 答.目前图像处理系统开发的主流工具为 Visual C++(面向对象可视化集成工具)和 MATLAB 的图像处理工具箱(Image Processing Tool box)。两种开发工具各有所长且有相互 间的软件接口。 Microsoft 公司的 VC++是一种具有高度综合性能的面向对象可视化集成工具,用它开发 出来的 Win 32 程序有着运行速度快、可移植能力强等优点。VC++所提供的 Microsoft 基础 类库 MFC 对大部分与用户设计有关的 Win 32 应用程序接口 API 进行了封装,提高了代码 的可重用性,大大缩短了应用程序开发周期,降低了开发成本。由于图像格式多且复杂,为 了减轻程序员将主要精力放在特定问题的图像处理算法上,VC++ 6.0 提供的动态链接库 ImageLoad.dll 支持 BMP、JPG、TIF 等常用 6 种格式的读写功能。 MATLAB 的图像处理工具箱 MATLAB 是由 MathWorks 公司推出的用于数值计算的有 力工具,是一种第四代计算机语言,它具有相当强大的矩阵运算和操作功能,力求使人们摆 脱繁杂的程序代码。MATLAB 图像处理工具箱提供了丰富的图像处理函数,灵活运用这些 函数可以完成大部分图像处理工作,从而大大节省编写低层算法代码的时间,避免程序设计 中的重复劳动。MATLAB 图像处理工具箱涵盖了在工程实践中经常遇到的图像处理手段和 算法,如图形句柄、图像的表示、图像变换、二维滤波器、图像增强、四叉树分解域边缘检 测、二值图像处理、小波分析、分形几何、图形用户界面等。但是,MATLAB 也存在不足 之处限制了其在图像处理软件中实际应用。首先,强大的功能只能在安装有 MATLAB 系统 的机器上使用图像处理工具箱中的函数或自编的 m 文件来实现。其次,MATLAB 使用行解 释方式执行代码,执行速度很慢。第三,MATLAB 擅长矩阵运算,但对于循环处理和图形 界面的处理不及 C++等语言。为此,通应用程序接口 API 和编译器与其他高级语言(如 C、 C++、Java 等)混合编程将会发挥各种程序设计语言之长协同完成图像处理任务。API 支持 MATLAB 与外部数据与程序的交互。编译器产生独立于 MATLAB 环境的程序,从而使其他 语言的应用程序使用 MATLAB。

数字图像处理课件第6章图像的几何变换

数字图像处理课件第6章图像的几何变换
由点的齐次坐标(Hx, Hy, H)求点的规范化齐次坐标(x, y, 1),可按下式进行:
x Hx H
y Hy H
第6章 图像的几何变换
齐次坐标的几何意义相当于点(x, y)落在3D空间H=1
的平面上,如图6-2所示。如果将xOy平面内的三角形abc的 各顶点表示成齐次坐标(xi, yi, 1)(i=1, 2, 3)的形式,就变成H =1平面内的三角形a1b1c1的各顶点。
图6-2 齐次坐标的几何意义
第6章 图像的几何变换
齐次坐标在2D图像几何变换中的另一个应用是:如某 点S(60 000,40 000)在16位计算机上表示,由于大于32767 的最大坐标值,需要进行复杂的处理操作。但如果把S的坐 标形式变成(Hx, Hy, H)形式的齐次坐标,则情况就不同了。 在齐次坐标系中,设H=1/2,则S(60 000,40 000)的齐次坐 标为(x/2,y/2,1/2),那么所要表示的点变为(30 000, 20 000,1/2),此点显然在16位计算机上二进制数所能表示 的范围之内。
(图像上各点的新齐次坐标)
(图像上各点的原齐次坐标)
第6章 图像的几何变换 设变换矩阵T为
a b p
T c
d
q
l m s
则上述变换可以用公式表示为
=
T
Hx1' Hy1'
Hx2' Hy2'
Hxn' Hyn'
x1 x2 xn
T
y1
y2
yn
H H H 3n
1 1 1 3n
第6章 图像的几何变换
6.4 图像镜像
6.4.1 图像镜像变换 图像的镜像(Mirror)变换不改变图像的形状。 镜像变换分为两种:一种是水平镜像,另外一种是垂直镜

图像处理 第6章图像复原

图像处理  第6章图像复原

8
第6章 图像复原 §6.1.1 模拟图像退化的数学模型 一、退化模型
n(x,y) f (x,y) H g(x,y)

模型化:一个作用在f (x,y)上的系统H与一个加性噪声n (x,y)的联合作用,导致产生退化图像g (x,y) 。 假设已知n (x,y)的统计特性(或先求出),图像复原就 是已知g (x,y)求f (x,y)的问题 (近似于求解逼近过程),由 于解不唯一,故方法很多。不同误差准则,不同约束条件,得 到解不同。 g (x,y) = H [f (x,y)] + n (x,y) 已知 退化 解 噪声
0 x B 1和0 y D -1 B x M-1或D y N-1
ge (x,y) = fe (m,n) he (x-m,y-n) m=0,1,…,M-1; n=0,1,…,N-1 M = A+B-1, N = C+D-1 考虑噪声有: ge (x,y) = fe (m,n) he (x-m,y-n) +ne (x,y) ; m=0,1,…,M-1; n=0,1,…,N-1; ne (x,y)为M×N的噪声项
可见H是一个循环阵, 退化系统由H决定。
若A=4, B=3,则M=6,因为B=3,即he(3)=he(4)=he(5)=0
0 0 0 h(2) h(1) h(0) h (1) h (0) 0 0 0 h (2) h(2) h(1) h(0) 0 0 0 H 0 h (2) h (1) h (0) 0 0 0 0 h(2) h(1) h(0) 0 0 0 0 h(2) h(1) h(0)
13
第6章 图像复原 §6.1.2 离散图像退化的数学模型 二、2D情况:

数字图像处理6ppt课件

数字图像处理6ppt课件
数字图像处理6数字图像处理数字图像处理6图像恢复61退化模型及恢复技术基础61退化模型及恢复技术基础62空间域滤波恢复62空间域滤波恢复63频率域滤波恢复63频率域滤波恢复64逆滤波64逆滤波65最小均方误差滤波器维纳滤波65最小均方误差滤波器维纳滤波数字图像处理6图像恢复和图像增强一样都是为了改善图像视觉效果以及便于后续处理
img_median=medfilt2(img_noise); %对附加有椒盐噪声的图像实行中 值滤波
figure; imshow(img_median,[]); %显示中值滤波后的图像
img_median2=medfilt2(img_median); %对中值滤波处理后的图像再次 实行中值滤波
figure; imshow(img_median2,[]); %显示再次中值滤波后的图像
erage',3));
figure; imshow(img_mean,[]); %显示逆谐波滤波后的图像
Q=1.5;
%对高斯噪声图像实行Q取正数的逆谐波滤波
img_mean=imfilter(img_noise.^(Q+1),fspecial('average',3))./imfilter(img_noise.^Q,fspecial('av
%矩阵点乘实现频域滤波
out = ifftshift(out);
%原点移回左上角
out = ifft2(out);
%傅里叶逆变换
out = abs(out);
%取绝对值
out = out/max(out(:)); figure,imshow(out,[]);
%归一化 %显示滤波结果数字图像处理6
数字图像处理6
for i=1:M

第6章 点运算

第6章 点运算
DIP 第6章 点运算
6.1 引言
数字图像处理
Digital Image Processing
点运算( 点运算(Point operation)是相对于全局运算 operation)是相对于全局运算 局部运算而言 而言。 和局部运算而言。 点运算以某种预定的方式改变一幅图像的灰度 点运算以某种预定的方式改变一幅图像的灰度 以某种预定的方式 直方图。 直方图。
B B
Digital Image Processing
d f (DA ) dD A
H B ( D) =
f ′ f −1 ( D )
得:
H A [ ( D − b) / a ] 1 D −b H B ( D) = = HA f ′ [ ( D − b) / a ] a a
df 其中: 其中: f ′ = ,D即DB 即 dD
DIP
数字图像处理 6.2.2 举例
6.2.2.1 线性点运算 来说, )/a 对DB = f (DA)=aDA+b来说,DA = f -1(DB)= (DB-b)/a, 因此, 因此,根据式 H ( D ) = H A ( D A ) 或 H A f −1 ( D )
DIP
数字图像处理
6.1.2.2 非线性单调点运算
Digital Image Processing
255 右图为几种单调非减灰度变换函 右图为几种单调非减 单调非减灰度变换函 正斜率, 它们都有正斜率 数,它们都有正斜率,因而可以保留 输出 图像的基本外貌,但图像局部对比度 图像的基本外貌,但图像局部对比度 C>0:中间灰度增加 > : 会发生变化。 会发生变化。
中国传媒大学信息工程学院
13

数字图像处理第6章二值图像处理-专业文档资料

数字图像处理第6章二值图像处理-专业文档资料

二阶矩则描述了图像的对于直线和对轴与轴的转动惯量,因 此常常也把物体的二阶矩称为惯性矩。
中心矩 :
p q (x x)p(y y )qf(x ,y )d xp d ,q y 0 ,1 ,2
第6章 二值图像处理
低阶矩主要描述区域的面积、转动惯量、质心等等,具有 明显得几何意义,,四阶矩描述峰值的状态等等,一般 来说高阶矩受到图像离散化等的影响,高阶矩一般在应用中 不一定十分准确。
D e(ac)2(bd)2
② 街区距离,用Ds来表示:
(6-1)
D s |ac||bd|
③ 棋盘距离,用Dg表示如下:
(6-2)
D gma a x c|, ( |b|d|)
(6-3)
三者之间的关系为:Dg Ds,如De图6-1(a)、(b)和(c)所示。
第6章 二值图像处理
(a) 欧氏距离 (b) 街区距离 (c) 棋盘距离 (d)≤2构成菱形 (e)≤2构成正方形 图6-1 三种距离示意图
第6章 二值图像处理
6.2 二值图像的几何特征描述
6.2.1 二值图像中曲线的描述 6.2.1.1 轮廓跟踪-甲虫算法
目标区域的边界轮廓是描述目标的重要特征,对于二 值图像中的目标区域轮廓可以通过一种简单的轮廓跟踪算 法来得到,这种方法也被称作甲虫算法。如图6-6所示的二 值图像4连通分量,假定目标区域用1(黑色)表示,背景区域
1 (x,y)(x,y)
f(x,y)
0
else
M1N1
那么区域的面积为: S f (x, y) x0 y0
如果经过目标标记,区域占有的连通分量有k个,那么目
标区域的面积则是k个连通分量的面积总和,即有:
k
S Si i 1

数字图像处理第6章_图像编码与压缩技术.

数字图像处理第6章_图像编码与压缩技术.

霍夫曼编码
例 假设一个文件中出现了8种符号S0、S1、S2、S3、S4、S5、S6、 S7,那么每种符号编码至少需要3bit S0=000, S1=001, S2=010, S3=011, S4=100, S5=101, S6=110, S7=111 那么,符号序列S0 S1 S7 S0 S1 S6 S2 S2 S3 S4 S5 S0 S0 S1编码后 000 001 111 000 001 110 010 010 011 100 101 000 000 001 (共42bit) 和等长编码不同的一种方法是可变长编码。在这种编码方法中, 表示符号的码字的长度不是固定不变的,而是随着符号出现的概率 而变化,对于那些出现概率大的信息符号编以较短的字长的码,而 对于那些出现概率小的信息符号编以较长的字长的码。
6.3.3 霍夫曼编码
霍夫曼(Huffman)编码是根据可变长最佳编码定理,应用霍夫曼算
1.
对于每个符号,例如经过量化后的图像数据,如果对它们每 个值都是以相同长度的二进制码表示的,则称为等长编码或均匀 编码。采用等长编码的优点是编码过程和解码过程简单,但由于 这种编码方法没有考虑各个符号出现的概率,实际上就是将它们 当作等概率事件处理的,因而它的编码效率比较低。例6.3给出了 一个等长编码的例子。
6.1.1 图像的信息冗余
图像数据的压缩是基于图像存在冗余这种特性。压缩就是去掉 信息中的冗余,即保留不确定的信息,去掉确定的信息(可推知 的);也就是用一种更接近信息本身的描述代替原有冗余的描述。 8 (1) 空间冗余。在同一幅图像中,规则物体或规则背景的物理表 面特性具有的相关性,这种相关性会使它们的图像结构趋于有序和 平滑,表现出空间数据的冗余。邻近像素灰度分布的相关性很强。 (2) 频间冗余。多谱段图像中各谱段图像对应像素之间灰度相关 (3) 时间冗余。对于动画或电视图像所形成的图像序列(帧序 列),相邻两帧图像之间有较大的相关性,其中有很多局部甚至完

数字图像处理课后参考答案

数字图像处理课后参考答案

数字图像处理课后参考答案数字图像处理第⼀章1.1解释术语(2)数字图像:为了便于⽤计算机对图像进⾏处理,通过将⼆维连续(模拟)图像在空间上离散化,也即采样,并同时将⼆维连续图像的幅值等间隔的划分成多个等级(层次)也即均匀量化,以此来⽤⼆维数字阵列并表⽰其中各个像素的空间位置和每个像素的灰度级数的图像形式称为数字图像。

(3)图像处理:是指对图像信息进⾏加⼯以满⾜⼈的视觉或应⽤需求的⾏为。

1.7 包括图像变化、图像增强、图像恢复、图像压缩编码、图像的特征提取、形态学图像处理⽅法等。

彩⾊图像、多光谱图像和⾼光谱图像的处理技术沿⽤了前述的基本图像处理技术,也发展除了⼀些特有的图像处理技术和⽅法。

1.8基本思路是,或简单地突出图像中感兴趣的特征,或想⽅法显现图像中那些模糊了的细节,以使图像更清晰地被显⽰或更适合于⼈或及其的处理与分析。

1.9基本思路是,从图像退化的数学或概率模型出发,研究改进图像的外观,从⽽使恢复以后的图像尽可能地反映原始图像的本来⾯⽬,从⽽获得与景物真实⾯貌相像的图像。

1.10基本思路是,,在不损失图像质量或少损失图像质量的前提下,尽可能的减少图像的存储量,以满⾜图像存储和实时传输的应⽤需求。

1.11基本思路是,通过数学⽅法和图像变换算法对图像的某种变换,以便简化图像进⼀步处理过程,或在进⼀步的图像处理中获得更好的处理效果。

1.12基本⽬的是,找出便于区分和描述⼀幅图像中背景和⽬标的⽅法,以⽅便图像中感兴趣的⽬标的提取和描述。

第⼆章2.1解释下列术语(18)空间分辨率:定义为单位距离内可分辨的最少⿊⽩线对的数⽬,⽤于表⽰图像中可分辨的最⼩细节,主要取决于采样间隔值的⼤⼩。

(19)灰度分辨率:是指在灰度级别中可分辨的最⼩变化,通常把灰度级数L称为图像的灰度级分辨率。

(20)像素的4邻域:对于图像中位于(x,y)的像素p来说,与其⽔平相邻和垂直相邻的4个像素称为该像素的4邻域像素,他们的坐标分别为(x-1,y)(x,y-1)(x,y+1)(x+1,y)。

数字图像处理技术的应用第6章 图像编码

数字图像处理技术的应用第6章 图像编码

6.2 图像压缩概述
2、平均码字长度:
Assume:
kis第k个码字Ck的长度二进制代码的位数出现的概率pk
码字平均长度R:
M
R= k pk bit
R1
3、编码效率:
H 100%
R
6.2 图像压缩概述
4、冗余度:
r 1 r 可压缩的余地越小
6.2 图像压缩概述
1)数据冗余:将图像信息的描述方式改变之后,压缩 掉这些冗余。
2)主观视觉冗余:忽略一些视觉不太明显的微小差异, 可以进行所谓的“有损”压缩。
6.2 图像压缩概述
图像数字化关键是编码 compression code:在满足一定图像质量前提下,能获得减少数
据量的编码
一.Compression code及分类 研究处理的对象: 数据的物理容量
图像序列(x、y、t)50~200倍
6.2 图像压缩概述
3、从图像的光谱特征出发: 单色image coding; color image coding; 多光谱image coding。
4、从图像的灰度层次上: 多灰度编码; 二值图像code
5、从处理图像的维数出发;
行内coding; 帧内coding; 帧间code。
图像一大特点是数据量大,为其存贮、传输带来困难,需压缩。
eg:电话线传输速率一般为56Kbits/s(波特率) 一幅彩色图像512×512×24bit = 6M bits大小。传一幅图像需2分钟左右。 实时传送:512×512×24bits×25帧/秒=150Mbits/S 如压缩20倍,传一幅图6秒左右,可以接受,实用。 实时,要专用信道(卫星、微波网、专线网等技术)。 另外,大量资料需存贮遥感、故宫、医学CT、MR。

数字图像处理与应用(MATLAB版)第6章 图像的分割

数字图像处理与应用(MATLAB版)第6章 图像的分割

是边缘;
➢ 使用双阈值算法检测和连接边缘。即使用直方图计
算两个阈值,凡是大于高阈值的一定是边缘;凡是
小于低阈值的一定不是边缘。如果检测结果大于低
阈值但又小于高阈值,那就要看这个像素的邻接像
素中有没有超过高阈值的边缘像素,如果有,则该
像素就是边缘,否则就不是边缘。
0 -1 0 -1 4 -1 0 -1 0
B A
6.1 图像分割的定义和分类
图像分割:是指根据灰度、彩色、纹理等特征把图像 划分成若干个互不相交的区域,使得这些特征在同一区 域内,表现出一致性或相似性,而在不同区域间表现出 明显的不同。
图像分割的作用
图像分割是图像识别和图像理解的前提,图像分 割质量的好坏直接影响后续图像处理的效果。
图像
具体步骤:
➢ 首先用2D高斯滤波模板进行卷积以平滑图像;
➢ 利 用 微 分 算 子 ( 如 Roberts 算 子 、 Prewitt 算 子 和
Sobel算子等),计算梯度的幅值和方向;
➢ 对梯度幅值进行非极大值抑制。即遍历图像,若某
个像素的灰度值与其梯度方向上前后两个像素的灰

度值相比不是最大,那么这个像素值置为0,即不
第六章 图像的分割
内 容 1、图像分割的定义和分类; 提 2、基于边缘的图像分割方法;
要 3、基于区域的分割;
4、基于运动的图像分割 ; 5、图像分割技术的发展。

本 要
通过对图像分割技术的学习,掌
求 握基于边缘、区域、运动的图像

分割技术。

难 点
图像分割的定义、分类 基于边缘的图像分割方法
基于区域、运动的图像分割方法
G(i, j) Px Py

数字图像处理每章课后题参考答案

数字图像处理每章课后题参考答案

数字图像处理每章课后题参考答案第一章和第二章作业:1.简述数字图像处理的研究内容。

2.什么是图像工程?根据抽象程度和研究方法等的不同,图像工程可分为哪几个层次?每个层次包含哪些研究内容?3.列举并简述常用表色系。

1.简述数字图像处理的研究内容?答:数字图像处理的主要研究内容,根据其主要的处理流程与处理目标大致可以分为图像信息的描述、图像信息的处理、图像信息的分析、图像信息的编码以及图像信息的显示等几个方面,将这几个方面展开,具体有以下的研究方向:1.图像数字化,2.图像增强,3.图像几何变换,4.图像恢复,5.图像重建,6.图像隐藏,7.图像变换,8.图像编码,9.图像识别与理解。

2.什么是图像工程?根据抽象程度和研究方法等的不同,图像工程可分为哪几个层次?每个层次包含哪些研究内容?答:图像工程是一门系统地研究各种图像理论、技术和应用的新的交叉科学。

根据抽象程度、研究方法、操作对象和数据量等的不同,图像工程可分为三个层次:图像处理、图像分析、图像理解。

图像处理着重强调在图像之间进行的变换。

比较狭义的图像处理主要满足对图像进行各种加工以改善图像的视觉效果。

图像处理主要在图像的像素级上进行处理,处理的数据量非常大。

图像分析则主要是对图像中感兴趣的目标进行检测和测量,以获得它们的客观信息从而建立对图像的描述。

图像分析处于中层,分割和特征提取把原来以像素描述的图像转变成比较简洁的非图形式描述。

图像理解的重点是进一步研究图像中各目标的性质和它们之间的相互联系,并得出对图像内容含义的理解以及对原来客观场景的解释,从而指导和规划行为。

图像理解主要描述高层的操作,基本上根据较抽象地描述进行解析、判断、决策,其处理过程与方法与人类的思维推理有许多相似之处。

第三章图像基本概念1.图像量化时,如果量化级比较小时会出现什么现象?为什么?答:当实际场景中存在如天空、白色墙面、人脸等灰度变化比较平缓的区域时,采用比较低的量化级数,则这类图像会在画面上产生伪轮廓(即原始场景中不存在的轮廓)。

第6章 遥感数字图像处理_图像变换(2)

第6章 遥感数字图像处理_图像变换(2)
NDห้องสมุดไป่ตู้I


IR R IR R

式中:IR为遥感多波段图像中的近红外 (infrared)波段;R为红波段。 利用植被指数可监测某一区域农作物长势,并 在此基础上建立农作物估产模型,从而进行大 面积的农作物估产。
南京紫金山和玄武湖的NDVI分布
LANDSAT7的ETM影像,2000.6
常用的红外(IR)与红(R)波段

其中, R、G、B ∈[0, 1],r,g,b ∈[0, 1],M=max[R、 G、B],m=min[r、g、b] 注意,R、G、B中至少有一个值是0,与最大值的 颜色对应,并且至少有一个的值是1,与最小值 的颜色对应。
RGB到HSI
I M m 2
如果 M m , S 0 如果 I 0 . 5, S 如果 I 0 . 5, S M m M m M -m 2M m , S 的取值范围是 [ 0 ,1]

例如,在地质探测中,地质学家用TM的某种组 合解译矿石类型:B3/B1突出铁氧化物,B5/B7 突出粘土矿物,B5/B4突出铁矿石,B5/B6突出 大片白陶土蚀变区域,B4/B3突出植被信息, B5/B2分离陆地和水体,等等。

波段比值方法还可以用来探测地物随季节变化 的信息。例如,如果需要监测地区植被的变化, 可以使用不同季节的第3波段的比值,新建立的 波段可能是20060810B3/20040810B3。图像的 时段可以是不同年的同一个月,或同一年的不 同月,新产生的波段将突出变化信息,变化的 像素具有较高的亮度值。没有变化的像素值较 低,在图像中比较暗。
传感器Landsat TM所对应的指数函数
函数名称
归一化植被指数(NDVI) 比值植被指数(IR/R) 差值植被指数(Veg.index) 转换植被指数(TNDVI) 氧化铁指数(IRON OXIDE)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

图像分割—区域生长
(1)人机交互、直方图峰值、聚类中心等; (2)取决于问题本身以及图像数据的类型 ; (3) 要求图像分割后区域的可描述性。
再对g(x)求一阶或二阶导数以检测边缘点。利用卷积性质,有 g(x,y) = f (x,y) h(x,y); g(x,y) = f (x,y) h(x,y);
Marr边缘检测算法—— 平滑滤波器采用高斯滤波器; 边缘检测判据是二阶导数零交叉点; 采用线性插值的方法估计边缘的位臵。 因为采用Laplacian算子,故称也为LoG(Laplacian of Gaussian)滤波器。
图像分割—边缘检测
■ Canny算子 ——图像边缘检测必须满足两个条件:能有效地抑制噪声;尽量 精确确定边缘的位臵。 根据对信噪比与定位乘积进行测度,得到最优化逼近算子—— Canny边缘检测算子。 类似于Marr边缘检测,也属于先平滑后求导数的方法——在与边 缘垂直的方向上求零交叉。性能优于Marr算子。
理论上讲,边缘检测算子仅得到处在边缘上的像素点。实际 上,由于噪声、不均匀照明而产生的边缘间断以及其他由于引入 虚假的亮度间断所带来的影响,使得到的一组像素很少能完整地 描绘一条边缘——典型的做法是在使用边缘检测算法之后,使用 连接过程将边缘像素组合成有意义的边缘。 原理——分析图像中每个边缘点的一个小邻域内像素的性质,将 所有相似点连接起来,就是共同满足一定准则的像素组成的一条 边缘。
图像分割—边缘检测
将参数空间划分, 每一个“格子”所 通过的直线数量, 反映了图像空间中 边缘点的“共线性”
——可以通过“投 票胜出”方式确定 最终的有效格子, 即图像中的有效线 段。 a
b
图像分割—边缘检测
使用等式y = ax+b表示一条直线带来的一个问题是,当直线接近 垂直时,斜率接近无穷大。也就是说a的最大接近于无穷。解决的 方法是采用直线的标准式(在-空间中):
图像分割—门限法
将一个背景点当作目标点进行分类时,错误概率为:
E1 T p2 z dz
T
将一个目标点当作背景点进行分类时,错误概率为:
E2 T


T
p1 z dz
出错率的整体概率是:
E T P2 E1 T P 1 E2 T
对E(T)求导并令导数为0,得
图像分割—边缘检测
Marr-Hildreth
Canny
图像分割—边缘检测
3、间断检测
间断检测包括点检测,线检测和边界检测三种。寻找间断最 一般的方法是模板检测——计算模板所包围区域的灰度级与模板 系数的乘积之和,若大于阈值,则视为间断处。
——可用图示模板检测孤立点。 孤立点:该点的灰度级与其背景的差异相 当大,并且它所在的位臵是一个均匀的或 者近似均匀的区域。 基本思想:如果一个孤立点与它周围的点很不相同,则很容易 被上述模板检测到。在灰度级为常数的区域,模板响应为0。
图像分割—门限法
一、基于不连续性的分割方法
1、门限法(阈值分割)
LE g ( x, y ) LB f ( x, y ) T f ( x, y ) T
阈值T 的选择—— 基于直方图
单目标和背景,单一门限
多个目标和背景,多门限
图像分割—门限法
图像分割—门限法
图像分割—门限法
图像分割—边缘检测
——线模板
第一个模板对水平方向的线条(单象素宽)有最强的响应。 第二个模板对于45线有最强响应。 第三个模板对垂直线有最强响应。 第四个模板对于45 线有最强响应。 每个模板系数相加总和为0,保证了在灰度级恒定的区域,模板 响应为0。
图像分割—边缘检测
4、边缘连接和边界检测
多边形分段线性拟合原理
图像分割—边缘检测
C B D C B D C B
A
A
A
迭代端点拟合
图像分割—边缘检测
Hough变换: ——利用图像的全局特性而直接检测目标轮廓。 基本思想:点-线的对偶性(duality)——图像变换前在图像空间, 变换后在参数空间。在图像空间里,直线方程:y = ax+b ;给定 (x,y) 后,可以写成:b =-ax+y ;即图像空间中的一个点对应参 数空间 {a-b} 中的一条直线。 ——图像空间中的两个点,确定了参数空间中两条直线的交点。
图像分割—边缘检测
(采用5×5的均值滤波器预先进行了平滑处理的结果)
图像分割—边缘检测
对角方向边缘检测。(a)用Sobel +45模板。 (b)用Sobel -45模板。
从图中可以看到两个模板对水平和垂直的边 缘具有相似的响应,但要比在这两个方向上 的水平和垂直模板的响应要弱。
图像分割—边缘检测
(a) 血管造影图像。 (b) 图像的Sobel梯度。 (c)用于得到27×27的空 间平滑处理模板的高 斯函数。 (d)Laplacian模板。 (e)通过对原图使用高斯 型平滑模板进行平滑 处理,然后使用拉普 拉斯算子模板后得到 的LoG图像。 (f)对LoG设臵门限后的 结果。 (g)零交点。
噪声的影响
图像分割—边缘检测
边缘检测算子:
■ Roberts算子
■ Prewitt算子
■ Sobel算子
图像分割—边缘检测
■用于检测对角线方向上突变的Prewitt和Sobel模板
图像分割—边缘检测
a c b d
y x
(a)原图,(b) x方向上的梯度分量|Gx| ,(c) y方向上
的梯度分量|Gy|,(d)梯度图像|Gx|+|Gy|
多数情况下,目标和背景的灰度分布有重叠。若二者的 灰度分布的概率密度函数已知,则可以选择门限使得错误概 率最小(统计最优)。
背景
图像中两个区域的灰度级概率密 度函数
目标
图像整体灰度级变化的总概率密度函数:
pz P 1 p1 z P 2 p2 z
(P1和P2是两类象素出现的概率)
(x,y)是(x0,y0)的邻域,(x,y)是该点的梯度方向角。
图像分割—边缘检测
迭代端点拟合: 1.确定一个散布在两个特定边缘点A和B之间的边缘点的集合P; 2.设A,B组成一条(初始)边缘线,计算其它边缘点到直线的距离;
3.若所有点距直线的距离小于阈值,则拟合此直线为确定边缘;
4.否则,找出距离最大的点(记为C)为新端点,将P分为包含AC 线段的边缘点子集合P1和包含BC线段的子集合P2; 5. 重复3、4步骤,直到所有子集合的点到直线的距离小于阈值, 确定最后的各线段的端点。 6. 利用每个线段间的边缘点,通过拟合确定边缘。
■ Laplacian算子 一阶导数的局部最大值对应二阶导数的零交叉点(Zero crossing)。 所以通过求图像的二阶导数的零交叉就能找到精确边缘点。 拉氏算子对噪声敏感,因此在实际中通常不直接使用。一般结合平 滑算子使用。
图像分割—边缘检测
■ Marr-Hildreth算子 一般情况下,需将平滑滤波与边缘检测结合。 先做滤波: g(x,y) = f (x,y)h(x,y)
(d) 依据一定准则判断为相连 的像素集合。
准则:像素属于3个具有最高 计数的累加器单元;没有大于 5个像素的间隙。
图像分割—区域生长
二、基于区域的分割方法
1. 区域生长 ——对每个需要分割的区域找一个种子像素作为生长的起点,将 种子像素邻域中与其有相似性质的像素合并到种子像素所在的区 域中。将这些新像素当作新的种子像素重复上述过程,直到再没 有满足条件的像素被包括进来。 三个问题: (1)选择或确定一组能正确代表所需区域的种子像素; (2)确定在生长过程中将相邻像素合并的准则(相似性判据); (3)确定在生长过程中停止的条件或规则。
边缘模型:
阶跃边缘
渐变边缘
图像分割—边缘检测
一阶导数可以用于检测图像中的 一个点是否是边缘点; 二阶导数的符号可以判断一个边 缘像素的位臵。
图像分割—边缘检测
图中第一列的图像显 示了4个斜坡边缘的 特写,分别被均值为 0且σ= 0,0.1,1,10 的随机高斯噪声污染。 第二列是一阶导数图 像和灰度级剖面线。 第三列为二阶导数图 像和灰度级剖面线。

理想边缘
图像分割—边缘检测
(a)输入图像 (b)Sobel梯度的y分量 (c)梯度的x分量
(d)边缘连接的结果
a bc dFra bibliotek基于梯度一致性: 要求:
f x, y f x0 , y0 E
x, y x0 , y0 A
由于光照不均匀,导致图像整体有明暗分布(但目标与背景的相 对灰度有差异)
——将图像划分成子区域(或用移动窗)对子区域进行局部门限分割。
图像分割—门限法
(a)原图; (b)全局门限分割结 果; (c) 图像划分成子 区域; (d)自适应门限分割 结果。
a b c d
门限法的缺点
图像分割—门限法
最优门限的选取
图像分割—边缘检测
h x, y 1 2 2 e
x2 y2 2 2
2h 2h LoG h x, y 2 2 x y
2

1 x y 1 e 4 2 2
2 2

x2 y2 2 2
55Log模板
图像分割—边缘检测
数 字 图 像 处 理
第六章 图像分割
图像分割
门限法 边缘检测 区域生长 分裂与合并 分水岭方法
图像分割
图像输入 光电转换 数字化
预处理 图像增强 图像恢复 图像编码
图像分割 阈值分割 边缘检测 区域分割
描述 解释
图像分析理解
图像识别
特征提取
一般的图像处理过程
图像分割
在对图像的研究和应用中,往往对图像中的部分区域感兴趣,这 些区域被称为前景或目标——具有独特性质的区域。 定义:把图像划分为不相交的有意义的区域,同一区域内的像素 反映了目标的一致性属性(同一性)。 在无任何先验知识的情况下,一致性的依据和准则取决于图像的 灰度和颜色。 图像分割(Image segmentation)是图像分析中的关键技术—— 由于目标特征不均匀,场景复杂,图像中有噪声等原因——分割 困难。 图像分割方法一般是基于像素的两个特性之一:不连续性和相似 性。 不连续性:指不同区域之间像素值差别比较大,在边界上具有某 种不连续性——基于不连续性的分割方法 相似性:指区域内部象像素值之间具备一定的相似性——基于区 域的分割方法
相关文档
最新文档