Matlab数字图像处理.ppt
合集下载
第五篇Matlab图像处理PPT课件
3)immovie 以动画方式显示图像各帧,既将各帧图片转换为matlab动画格式;只能用索引图
像; mov=immovie(X,map)
实例:5-7.m
4)subimage subplot以子图方式在一个图像窗口中显示多幅图像,但整个图像窗口采用一个调
色板,subimage可以在各个子图中采用各自的调色板; 实例:5-8.m
3)pixval 交互显示单个像素的亮度值;显示欧几里得距离;
4)imfinfo
第4页/共40页
5)size(I) 获得图像矩阵大小;
6)保存图像 imwrite(I,filename,fmt) imwrite(I,map,filename,fmt) I—待保存的图像矩阵;filename—图像保存路
第20页/共40页
5)zoom 缩放图像
zoom on,zoom off 打开或关闭缩放功能
zoom out
恢复图像的原始尺寸
zoom reset
以当前图像尺寸作为缩放起点
zoom xon,zoom yon 设置X或Y轴缩放功能
6)warp
纹理映射
imshow显示的图像在二维平面上,纹理映射可以把图像显示在其他类型的表 面,如柱面、球面等等;
为uint8和uint16,分别对应[0 255]、[0 65535]
2)二值图像 像素取值只能为0或1;逻辑矩阵;
B=logical(A) 非零为1;
第6页/共40页
3)索引图像 索引图像包括调色板和图像数据两部分,是把颜色进行排列、编号,图像数据
对应为该点像素的颜色序号而非颜色本身; 调色板为m×3矩阵,每一行代表一种颜色,各元素的值介于[0,1]之间,乘以
2)im2bw
数字图象处理及matlab的实现 ppt 课件 第0608章
x a b x y c d y mod N
• 保证|ad-bc|=1,如置换矩阵系数设为a=b=1,c=2,d=3, 置换次数n=20,则他们被当作密钥key,用于解密。
n
27
function [outImg]=arnold(inImg,key) %key=[Times,a,b,c,d]; figure(1), imshow(inImg); title('Original'); [iH iW]=size(inImg); outImg=uint8(zeros(iH,iW)); tempImg=inImg; if iH ~= iW error('The cover must be a square !'); return; elseif size(key,2) ~=5; error('The key must be 5 numbers !'); return; elseif (key(2)*key(5)-key(3)*key(4)) ~=1; error('The Arnold Matrix is not valid !'); return; end for i=1:key(1) % iTime for u=1:iH for v=1:iW temp=tempImg(u,v); ax=mod((key(2)*(u-1)+key(3)*(v-1)),iW)+1; ay=mod((key(4)*(u-1)+key(5)*(v-1)),iW)+1; outImg(ax,ay)=temp; end end tempImg=outImg; end outImg=tempImg; figure(2), imshow(outImg); title('Permuted'); imwrite(outImg, 'Encrypted.bmp');
Matlab数字图像处理PPT课件
Matlab数字图像处理PPT课件
/2、图像的点运算
三、灰度的对数变换
tclo kg s
c为尺度比例常数,s为源灰度值,t为变换后的目标灰 度值。k为常数。灰度的对数变换可以增强一幅图像 中较暗部分的细节,可用来扩展被压缩的高值图像中 的较暗像素。广泛应用于频谱图像的显示中。
Warning:log函数会对输入图像矩阵s中的每个元素进行
1边界选项symmetricreplicatecircular2尺寸选项samefull3模式选项corrconv三滤波器设计4空间域图像增强hfspecialtypeparameters合法取值功能average平均模板disk圆形领域的平均模板gaussian高斯模板laplacian拉普拉斯模板log高斯拉普拉斯模板prewittprewitt水平边缘检测算子sobelsobel水平边缘检测算子parameters为可选项是和所选定的滤波器类型type相关的配置参数如尺寸和标准差等
I=fft2(x);%快速傅里叶变换 I=fft2(x,m,n);
x为输入图像;m和n分别用于将x的第一和第二维规整到指定的长度。 当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。
I1=abs(I);%计算I的幅度谱 I2=angle(I);%计算I的相位谱
W(1,-1) W(1,0) W(1, 1)
f(x-1,y-1) f(x-1,y) f(x-1,y+1)
f(x,y-1) f(x,y)
f(x,y+1)
f(x+1,y-1) f(x+1,y) f(x+1,y+1)
B=imfilter(f,w,option1,option2,…);
f 为要进行滤波操作的图像。 w为滤波操作使用的模板,为一个二维数组,可自己定义。 option1……是可选项,包括: 1、边界选项(’symmetric’、’replicate’、’circular’) 2、尺寸选项(’same’、’full’) 3、模式选项(’corr’、’conv’)
/2、图像的点运算
三、灰度的对数变换
tclo kg s
c为尺度比例常数,s为源灰度值,t为变换后的目标灰 度值。k为常数。灰度的对数变换可以增强一幅图像 中较暗部分的细节,可用来扩展被压缩的高值图像中 的较暗像素。广泛应用于频谱图像的显示中。
Warning:log函数会对输入图像矩阵s中的每个元素进行
1边界选项symmetricreplicatecircular2尺寸选项samefull3模式选项corrconv三滤波器设计4空间域图像增强hfspecialtypeparameters合法取值功能average平均模板disk圆形领域的平均模板gaussian高斯模板laplacian拉普拉斯模板log高斯拉普拉斯模板prewittprewitt水平边缘检测算子sobelsobel水平边缘检测算子parameters为可选项是和所选定的滤波器类型type相关的配置参数如尺寸和标准差等
I=fft2(x);%快速傅里叶变换 I=fft2(x,m,n);
x为输入图像;m和n分别用于将x的第一和第二维规整到指定的长度。 当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。
I1=abs(I);%计算I的幅度谱 I2=angle(I);%计算I的相位谱
W(1,-1) W(1,0) W(1, 1)
f(x-1,y-1) f(x-1,y) f(x-1,y+1)
f(x,y-1) f(x,y)
f(x,y+1)
f(x+1,y-1) f(x+1,y) f(x+1,y+1)
B=imfilter(f,w,option1,option2,…);
f 为要进行滤波操作的图像。 w为滤波操作使用的模板,为一个二维数组,可自己定义。 option1……是可选项,包括: 1、边界选项(’symmetric’、’replicate’、’circular’) 2、尺寸选项(’same’、’full’) 3、模式选项(’corr’、’conv’)
MATLAB的图像处理介绍.ppt
可在matlab命令窗口输入:help imwrite,来获得提示信息。
2)imwrite获得的文件存放在matlab的work文件夹中。
(2)以图像的形式输出,函数为 image(imshow)
矩阵A
Image/imshow
Matlab输出的图像
image的格式是: image(C) image(x,y,C) image(…,’PropertyName’,PropertyValue,…) image(‘PropertyName’,PropertyValue,…) handle=image(…) 其中,x,y分别表示图像显示位置的左上角坐标,C表
4)RGB图像:图像中每个象素的颜色用三个数据来存储,分别指 定红、绿、蓝三原色在象素颜色中的比例关系,组成一个三维数 组,读入matlab后是一个三维的矩阵。
注意:美术教科书中称红、黄、蓝为三原色,讲的是绘画颜 料的使用。一般电视光色等光色是红、绿、蓝。RGB图像就 是采用红、绿、蓝作为三原色的,其中R为红色,G为绿色, B为蓝色。
3)象素:是图像在计算机显示中的度量单位,可以变 化,可大可小。 4)分辨率:是用于度量图像在显示器中清晰程度的一 个参数,分辨率越高,图像越清晰。分辨率是与象素相 关的,即单位长度上的象素数就是分辨率。由此可知, 分辨率越高,象素的几何尺寸就越小。 5)图像文件的大小:指一幅图像在计算机中保存时所 占用的磁盘空间,其大小与所用的颜色模式有关。灰度 图像中的每一个灰度象素只占用一个字节(8位), RGB图像中红、绿、蓝各占用一个字节。另外,图像文 件的大小也直接与其分辨率有关,原因是当分辨率增加 时,一幅图像所包含的象素量急剧增加。 6)句柄:就是对象的代号或标志,它能使计算机方便 地找到所需要的对象并加以相应的操作。MATLAB中的 句柄图形对象包括轴、文本、菜单、控制框、图像等。
2)imwrite获得的文件存放在matlab的work文件夹中。
(2)以图像的形式输出,函数为 image(imshow)
矩阵A
Image/imshow
Matlab输出的图像
image的格式是: image(C) image(x,y,C) image(…,’PropertyName’,PropertyValue,…) image(‘PropertyName’,PropertyValue,…) handle=image(…) 其中,x,y分别表示图像显示位置的左上角坐标,C表
4)RGB图像:图像中每个象素的颜色用三个数据来存储,分别指 定红、绿、蓝三原色在象素颜色中的比例关系,组成一个三维数 组,读入matlab后是一个三维的矩阵。
注意:美术教科书中称红、黄、蓝为三原色,讲的是绘画颜 料的使用。一般电视光色等光色是红、绿、蓝。RGB图像就 是采用红、绿、蓝作为三原色的,其中R为红色,G为绿色, B为蓝色。
3)象素:是图像在计算机显示中的度量单位,可以变 化,可大可小。 4)分辨率:是用于度量图像在显示器中清晰程度的一 个参数,分辨率越高,图像越清晰。分辨率是与象素相 关的,即单位长度上的象素数就是分辨率。由此可知, 分辨率越高,象素的几何尺寸就越小。 5)图像文件的大小:指一幅图像在计算机中保存时所 占用的磁盘空间,其大小与所用的颜色模式有关。灰度 图像中的每一个灰度象素只占用一个字节(8位), RGB图像中红、绿、蓝各占用一个字节。另外,图像文 件的大小也直接与其分辨率有关,原因是当分辨率增加 时,一幅图像所包含的象素量急剧增加。 6)句柄:就是对象的代号或标志,它能使计算机方便 地找到所需要的对象并加以相应的操作。MATLAB中的 句柄图形对象包括轴、文本、菜单、控制框、图像等。
数字图像处理ch01(MATLAB)-课件
2024/10/12
第一章 绪论
17
2024/10/12
第一章 绪论
18
2024/10/12
第一章 绪论
19
2024/10/12
第一章 绪论
20
<2>几何处理
放大、缩小、旋转,配准,几何校正,面积、周长计算。
请计算台湾的陆地面积
2024/10/12
第一章 绪论
21
<3>图象复原
由图象的退化模型,求出原始图象
图像处理是指按照一定的目标,用一系列的操 作来“改造”图像的方法.
2024/10/12
第一章 绪论
7
➢图象处理技术的分类(从方法上进行分类)[2]
1.模拟图象处理(光学图像处理等)
用光学、电子等方法对模拟信号组成的图像,用光学器 件、电子器件进行光学变换等处理得到所需结果(哈哈 镜、望远镜,放大镜,电视等).
2024/10/12
第一章 绪论
22
<4>图象重建[3]
[3]此图像来自罗立民,脑成像,
2024/10/12
第一章 绪论
23
/zhlshb/ct/lx.htm
2024/10/12
第一章 绪论
图形用户界面,动画,网页制作等
2024/10/12象处理的基本概念,和基 本问题,以及一些典型的应用。
2024/10/12
第一章 绪论
33
提问
摄像头(机),扫描仪,CT成像装置,其他图象成像装置
2)图象的存储
各种图象存储压缩格式(JPEG,MPEG等),海量图象数据库技术
3)图象的传输
内部传输(DirectMemoryAccess),外部传输(主要是网络)
Matlab数字图像处理-02-PPT课件
可选参数’nothinning’,指定时可以通过跳过边缘细化算法来加快算法 运行的速度。默认是’thinning’,即进行边缘细化。 2)、基于高斯-拉普拉斯算子的边缘检测
BW=edge(I,’log’,thresh,sigma)
sigma指定生成高斯滤波器所使用的标准差。默认时,标准差为2。滤 镜大小n*n,n的计算方法为:n=ceil(sigma*3)*2+1。
/5、频率域图像增强 一、傅里叶变换
I=fft2(x);%快速傅里叶变换 I=fft2(x,m,n);
x为输入图像;m和n分别用于将x的第一和第二维规整到指定的长度。 当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。
I1=abs(I);%计算I的幅度谱 I2=angle(I);%计算I的相位谱 Y=fftshift(I);%频谱平移
高斯低通
function out = imgaussflpf(I,sigma) %imgaussflpf函数 构造高斯频域低通滤波器 [M,N] = size(I); out = ones(M,N); for i=1:M forj=1:N out(i,j) = exp(-((i-M/2)^2+(j-N/2)^2)/2/sigma^2); end end
R 1 C G 1 M B 1 Y
cmy=imcomplement(rgb);%rgb转成cmy rgb=imcomplement(cmy);cmy转成rgb
HSI模型
HIS模型是从人的视觉系统出发,直接使用颜色三要素色调 (Hue)、饱和度(Saturation)和亮度(Intensity)来描述 颜色。
桥接有单个像素缝隙分割的前景像素 清楚孤立的前景像素 围绕对角线相连的前景像素进行填充 填充单个像素的孔洞
matlab与数字图像处理 PPT课件
hjian@
数字图像处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
The MATLAB Image Processing 2 Toolbox
• The Image Processing Toolbox is a collection
of MATLAB functions (called M-functions or Mfiles) that extend the capability of the MATLAB environment for the solution of digital image processing problems.
数字图像处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
10
Displaying an image(cont.)
• Spatial domain
数字图像处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
11
Displaying an image(cont.)
• pixval : 功能:cursor on image to show pixel values 用法: imshow(I),pixval
eddins36数字图像处理2008年中山大学信息科学与技术学院计算机系黄剑3个projects?电能表读数识别?人脸检测?印刷线路板缺陷检测37数字图像处理2008年中山大学信息科学与技术学院计算机系黄剑电能表读数识别?对字轮进行精确定位与识别
Matlab with DIP 教学
计算机科学系 黄剑
数字图ቤተ መጻሕፍቲ ባይዱ处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
4
How do I know M-function?
数字图像处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
The MATLAB Image Processing 2 Toolbox
• The Image Processing Toolbox is a collection
of MATLAB functions (called M-functions or Mfiles) that extend the capability of the MATLAB environment for the solution of digital image processing problems.
数字图像处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
10
Displaying an image(cont.)
• Spatial domain
数字图像处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
11
Displaying an image(cont.)
• pixval : 功能:cursor on image to show pixel values 用法: imshow(I),pixval
eddins36数字图像处理2008年中山大学信息科学与技术学院计算机系黄剑3个projects?电能表读数识别?人脸检测?印刷线路板缺陷检测37数字图像处理2008年中山大学信息科学与技术学院计算机系黄剑电能表读数识别?对字轮进行精确定位与识别
Matlab with DIP 教学
计算机科学系 黄剑
数字图ቤተ መጻሕፍቲ ባይዱ处理,2008年
中山大学信息科学与技术学院计算机系,黄剑
4
How do I know M-function?
基于MATLAB的数字图像处理_毕业设计PPT
二值图像分析
1.膨胀 膨胀是使区域从四周向 外扩大,给图像中的对象边 界添加像素。 2.腐蚀 腐蚀则是使区域从四周 同时向内缩小,删除对象边 界某些像素。
二值图像分析
3.综合应用 (1)噪声滤除 开运算:开运算是先对图像进行腐蚀,然后再进行膨胀的结果。 闭运算:闭运算是先对图像进行膨胀,然后再进行腐蚀的结果。 将开启运算和闭合运算结合起来可构成噪声滤除器。 MATLAB中还提供了预定义的形态函数bwmorph。 代码: I3=imnoise(I2,'salt & pepper'); I4=bwmorph(I3,'open'); I5=bwmorph(I4,'close');
图像复原
原始图像
模糊后的图像
加噪后的图像
恢复后的图像
ቤተ መጻሕፍቲ ባይዱ
运行后:
图像增强
空域滤波增强:
1.线性平滑滤波器 线性低通滤波器是最常用的线性平滑滤波器。这种滤波器的所有系数
都是正的。对3*3的模板来说,最简单的操作是取所有系数都为1。为保证输
出图像仍在原来的灰度范围内,在计算R后要将其除以9再进行赋值。这种 方法称为邻域平均法。
代码:I=imread('saturn.tif'); J=imnoise(I,'salt & pepper',0.02); imshow(I) figure,imshow(J) K1=filter2(fspecial('average',3),J)/255; figure,imshow(K1) ` title('3*3的均值滤波器')
MATLAB新特点: 1.数据类型和面向对象编程技术 2.控制流和函数类型 3.大量引入图形用户界面 4.引入了全方位帮助系统 5. M 文件编辑、调试 的集成环境和性能剖析 6.MATLAB 环境可运行文件的多样化
数字图像处理MATLAB图像处理PPT课件
C(x, y) A(x, y) B(x, y) 差值图像提供了图像间的差值信息,能用于指导动态监测、 运动目标的检测和跟踪、图像背景的消除及目标识别等。
主要应用举例: ➢ 差影法(检测同一场景两幅图像之间的变化) ➢ 混合图像的分离
第17页/共41页
(1)检测同一场景两幅图像之间的变化
设:时刻1的图像为 T1(x,y), 时刻2的图像为 T2(x,y) g(x,y) = T2 (x,y) - T1(x,y)
第25页/共41页
除法运算(Division)
4、除法运算
C(x, y) A(x, y) B(x, y)
简单的除法运算可用于改变图像的灰度级, 常用于遥感图像处理中。
在四种算术运算中,减法与加法在图像增强 处理中最为有用。
第26页/共41页
几何变换
1)简单变换 • 问题描述:图像的平移、放缩和旋转。 • 解题思路:从易到难。工具:线性代数中的齐次坐标。
x
象
素
y1
填y
充
映
射
f(x1,y1) (x1,y1)非整型
f(x,y) (x,y)整型
第32页/共41页
图像的缩放
• 两种映射方法的对比 • 对于向前映射:每个输出图像的灰度要经过多次运算; • 对于向后映射:每个输出图像的灰度只要经过一次运算。
实际应用中,更经常采用向后映射法。 其中,根据四个相邻像素灰度值计算某位置的像素灰度
1、加法运算
C(x, y) A(x, y) B(x, y)
主要应用举例: ➢ 去除“叠加性”随机噪音 ➢ 生成图像叠加效果
第12页/共41页
(1)去除“叠加性”噪音
对于原图象f(x,y),有一个噪音图像集 { g i (x ,y) } 其中:g i (x ,y) = f(x,y) + ei(x,y)
主要应用举例: ➢ 差影法(检测同一场景两幅图像之间的变化) ➢ 混合图像的分离
第17页/共41页
(1)检测同一场景两幅图像之间的变化
设:时刻1的图像为 T1(x,y), 时刻2的图像为 T2(x,y) g(x,y) = T2 (x,y) - T1(x,y)
第25页/共41页
除法运算(Division)
4、除法运算
C(x, y) A(x, y) B(x, y)
简单的除法运算可用于改变图像的灰度级, 常用于遥感图像处理中。
在四种算术运算中,减法与加法在图像增强 处理中最为有用。
第26页/共41页
几何变换
1)简单变换 • 问题描述:图像的平移、放缩和旋转。 • 解题思路:从易到难。工具:线性代数中的齐次坐标。
x
象
素
y1
填y
充
映
射
f(x1,y1) (x1,y1)非整型
f(x,y) (x,y)整型
第32页/共41页
图像的缩放
• 两种映射方法的对比 • 对于向前映射:每个输出图像的灰度要经过多次运算; • 对于向后映射:每个输出图像的灰度只要经过一次运算。
实际应用中,更经常采用向后映射法。 其中,根据四个相邻像素灰度值计算某位置的像素灰度
1、加法运算
C(x, y) A(x, y) B(x, y)
主要应用举例: ➢ 去除“叠加性”随机噪音 ➢ 生成图像叠加效果
第12页/共41页
(1)去除“叠加性”噪音
对于原图象f(x,y),有一个噪音图像集 { g i (x ,y) } 其中:g i (x ,y) = f(x,y) + ei(x,y)
数字图像处理及MATLAB实现PPT课件
第12页/共69页
8.3.2梯度图像二值化
• 如果用适中的阈值对一幅梯度图像进行二值化,Kirsch的分割法利用了这种现象。 • 算法步骤 • 用一个中偏低的灰度阈值对梯度图像进行二值化从而检测出物体和背景,物体与背景被处于阈值之上的边
界点带分开。随着阈值逐渐提高,就引起物体和背景的同时增长。当它们接触上而又不至于合并时,可用 接触点来定义边界。这是分水岭算法在梯度图像中的应用。
OTSU算法定义:该算法是在灰度直方图的基础上用最小二乘法原理推导出来的, 具有统计意义上的最佳分割阈值。
第8页/共69页
• OTSU基本原理:以最佳阈值将图像的灰度直方图分割成两部分,是两部分之 间的方差取最大值,即分离性最大。
第9页/共69页
3. 迭代法求阈值
原理:图像中前景与背景之间的灰度分布为相互不重叠,在该前提下,实现对 两类对象的阈值分割方法。
除非图像中的物体有陡峭的边沿,否则灰度阈值的取值对所抽取物体的边界的 定位和整体的尺寸有很大的影响。这意味着后续的尺寸(特别是面积)的测量 对于灰度阈值的选择很敏感。由于这个原因,我们需要一个最佳的,或至少是 具有一致性的方法确定阈值。
第5页/共69页
1.直方图技术
• 含有一个与背景明显对比的物体的图像其有包含双峰的灰度直方图
第21页/共69页
Sobel边缘算子图
第22页/共69页
Prewitt边缘算子
Prewitr边缘算子 第23页/共69页
Kirsch边缘算子
图像中的每个点均与这8个模板进行卷积,每个掩模对某个特定边缘方向作出最大响应。所有8个方向中的 最大值作为边缘幅度图像的输出。最大响应掩模的序号构成了对边缘方向的编码。 Kirsch算子的梯度幅度值
8.3.2梯度图像二值化
• 如果用适中的阈值对一幅梯度图像进行二值化,Kirsch的分割法利用了这种现象。 • 算法步骤 • 用一个中偏低的灰度阈值对梯度图像进行二值化从而检测出物体和背景,物体与背景被处于阈值之上的边
界点带分开。随着阈值逐渐提高,就引起物体和背景的同时增长。当它们接触上而又不至于合并时,可用 接触点来定义边界。这是分水岭算法在梯度图像中的应用。
OTSU算法定义:该算法是在灰度直方图的基础上用最小二乘法原理推导出来的, 具有统计意义上的最佳分割阈值。
第8页/共69页
• OTSU基本原理:以最佳阈值将图像的灰度直方图分割成两部分,是两部分之 间的方差取最大值,即分离性最大。
第9页/共69页
3. 迭代法求阈值
原理:图像中前景与背景之间的灰度分布为相互不重叠,在该前提下,实现对 两类对象的阈值分割方法。
除非图像中的物体有陡峭的边沿,否则灰度阈值的取值对所抽取物体的边界的 定位和整体的尺寸有很大的影响。这意味着后续的尺寸(特别是面积)的测量 对于灰度阈值的选择很敏感。由于这个原因,我们需要一个最佳的,或至少是 具有一致性的方法确定阈值。
第5页/共69页
1.直方图技术
• 含有一个与背景明显对比的物体的图像其有包含双峰的灰度直方图
第21页/共69页
Sobel边缘算子图
第22页/共69页
Prewitt边缘算子
Prewitr边缘算子 第23页/共69页
Kirsch边缘算子
图像中的每个点均与这8个模板进行卷积,每个掩模对某个特定边缘方向作出最大响应。所有8个方向中的 最大值作为边缘幅度图像的输出。最大响应掩模的序号构成了对边缘方向的编码。 Kirsch算子的梯度幅度值
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
在MATLAB软件中,结构元素定义为一个STREL对象。如果nhood为结构 元素定义的领域,则任意大小和维数的结构元素的原点坐标为: origin=floor((size(nhood)+1)/2)。在MATLAB软件中,采用函数strel( )创建 任意大小和形状的STREL对象。函数strel( )支持常用的形状,例如线型 (line)、矩形(rectangle)、方形(square)、球形(ball)、钻石型 (diamond)和自定义的任意型(arbitrary)等。
结构元素的映射为或,定义为:
的补集记作或,定义为:
两个集合和的差集记作,定义为:
对于两幅图像和,如果,则称击中,记作;否则,如果,则称击不中A。
12.1.2 结构元素
结构元素是膨胀和腐蚀的最基本组成部分,用于测试输入图像。二维结构 元素是由数值0和1组成的矩阵。结构元素的原点指定了图像中需要处理的 像素范围,结构元素中数值为1的点决定结构元素的领域像素在进行膨胀 或腐蚀操作时是否参与计算。
数学形态学方法比其他空域或频域图像处理和分析方法具有一些明显的优 势。例如,基于数学形态学的边缘信息提取处理优于基于微分运算的边缘 提取算法,它不像微分算法对噪声那样敏感,提取的边缘比较光滑;利用 数学形态学方法提取的图像骨架也比较连续,断点少等;数学形态学易于 用并行处理方法有效的实现,而且硬件实现容易。
12.2.4 图像的边界测定
对于灰度图像可以通过形态学的膨胀和腐蚀来获取图像的边缘。通过形态 学获取灰度图像的边缘的优点是对边缘的方向性依赖比较小。下面通过例 子程序进行说明。
ቤተ መጻሕፍቲ ባይዱ
AB x : Bx A
12.1.3 膨胀与腐蚀
膨胀是将与物体接触的所有背景点合并到该物体中,使边界向外部扩张的 过程。通过膨胀,可以填充图像中的小孔,及在图像边缘处的小凹陷部分。 结构元素对图像的膨胀,记作,定义为:
腐蚀和膨胀是对偶操作。腐蚀是一种消除边界点,使边界向内部收缩的过 程。利用腐蚀操作,可以消除小且无意义的物体。集合被结构元素腐蚀, 记作,定义为:
12.2.2 图像填充操作
在MATLAB软件中,采用函数imfill( )对二值图像或灰度图像进行填充操 作。函数imfill( )的调用格式为:
BW2=imfill(BW):该函数对二值图像BW进行填充操作,对于二维图像允 许用户通过鼠标选择填充的点。通过键盘上面的Backspace键或Delete键可 以取消当前选择的点,通过键盘上的Return键可以结束交互式的选择。
BW=imregionalmax(I):该函数获取灰度图像I的局部极大值,返回值BW 为和原图像大小相同的二值图像,BW中元素1对应的是极大值,其它元素 值为0。
BW=imregionalmax(I, conn):该函数中参数conn为连通类型。对于二维图 像conn可以取值为4和8,默认值为8。对于三维图像,conn可取值为6、18 和26,默认值为26。
I2=imfill(I):该函数对灰度图像进行填充操作,返回值I2也是灰度图像。
12.2.3 最大值和最小值
对于一幅图像可以有多个局部极大值或极小值,但只有一个最大值或最小 值。在MATLAB软件中,采用函数imregionalmax( )获取图像的所有局部极 大值,采用函数imregionalmin( )获取局部极小值。函数imregionalmax ( )的 调用格式为:
第12章 形态学图像处理
形态学是一种应用于图像处理和模式识别领域的新的方法,是一门建立在 严格的数学理论基础上而又密切联系实际的科学。由于形态学具有完备的 数学基础,这为形态学用于图像分析和处理等奠定了坚实的基础。本章详 细的介绍利用MATLAB软件进行形态学图像处理,主要内容包括基本运算、 组合形态学运算以及二值图像的其他形态学运算等内容。
12.1 基本的形态学运算
数学形态学可以看作是一种特殊的数字图象处理方法和理论,主要以图象 的形态特征为研究对象。它通过设计一整套运算、概念和算法,用以描述 图象的基本特征。这些数学工具不同于常用的频域或空域算法,而是建立 在微分几何以及随机集论的基础之上的。数学形态学作为一种用于数字图 象处理和识别的新理论和新方法,它的理论虽然很复杂,但它的基本思想 却是简单而完美的。
12.1.4 开运算和闭运算
结构元素对的开运算,记作,定义为:
即首先采用结构元素对做腐蚀运算,然后再做膨胀运算,使用相同的结构 元素。
闭运算是开运算的对偶运算,记作,定义为:
即首先用结构元素对做膨胀运算,然后再做腐蚀运算,使用相同的结构元 素。
12.2 组合形态学运算
图像的膨胀和腐蚀是最基本的形态学运算。下面介绍组合形态学运算,主 要内容包括高帽滤波和低帽滤波、边界提取、区域填充等等。
12.2.1 高帽滤波和低帽滤波
图像的形态学高帽滤波(top-hat filtering)定义为:
其中为输入的图像,为采用的结构元素,即从图像中减去形态学开操作后 的图像。通过高帽滤波可以增强图像的对比度。
图像的形态学低帽滤波(bottom-hat filtering)定义为:
其中为输入的图像,为采用的结构元素,即从图像中减去形态学闭操作后 的图像。通过低帽滤波可以获取图像的边缘。
12.1.1 基本概念
集合论是数学形态学的基础,首先对集合论的一些基本概念做一个简单介 绍。腐蚀运算和膨胀运算是数学形态学的两个基本变换。参加运算的对象 有两个:图像A(感兴趣目标)和结构集合B,B称为结构元素。结构元素 通常是个圆盘,但它其实可以是任何形状。
设和是的子集,则把图像沿矢量平移一段距离记作或,其定义为:
[BW2, locations]=imfill(BW):该函数中返回值locations包含了交互式选择 时的点的坐标。
BW2=imfill(BW, locations):该函数中通过参数locations指定了进行填充时 的点的坐标。
BW2=imfill(BW, ‘holes’):该函数通过参数holes可以填充二值图像中的 空洞。
结构元素的映射为或,定义为:
的补集记作或,定义为:
两个集合和的差集记作,定义为:
对于两幅图像和,如果,则称击中,记作;否则,如果,则称击不中A。
12.1.2 结构元素
结构元素是膨胀和腐蚀的最基本组成部分,用于测试输入图像。二维结构 元素是由数值0和1组成的矩阵。结构元素的原点指定了图像中需要处理的 像素范围,结构元素中数值为1的点决定结构元素的领域像素在进行膨胀 或腐蚀操作时是否参与计算。
数学形态学方法比其他空域或频域图像处理和分析方法具有一些明显的优 势。例如,基于数学形态学的边缘信息提取处理优于基于微分运算的边缘 提取算法,它不像微分算法对噪声那样敏感,提取的边缘比较光滑;利用 数学形态学方法提取的图像骨架也比较连续,断点少等;数学形态学易于 用并行处理方法有效的实现,而且硬件实现容易。
12.2.4 图像的边界测定
对于灰度图像可以通过形态学的膨胀和腐蚀来获取图像的边缘。通过形态 学获取灰度图像的边缘的优点是对边缘的方向性依赖比较小。下面通过例 子程序进行说明。
ቤተ መጻሕፍቲ ባይዱ
AB x : Bx A
12.1.3 膨胀与腐蚀
膨胀是将与物体接触的所有背景点合并到该物体中,使边界向外部扩张的 过程。通过膨胀,可以填充图像中的小孔,及在图像边缘处的小凹陷部分。 结构元素对图像的膨胀,记作,定义为:
腐蚀和膨胀是对偶操作。腐蚀是一种消除边界点,使边界向内部收缩的过 程。利用腐蚀操作,可以消除小且无意义的物体。集合被结构元素腐蚀, 记作,定义为:
12.2.2 图像填充操作
在MATLAB软件中,采用函数imfill( )对二值图像或灰度图像进行填充操 作。函数imfill( )的调用格式为:
BW2=imfill(BW):该函数对二值图像BW进行填充操作,对于二维图像允 许用户通过鼠标选择填充的点。通过键盘上面的Backspace键或Delete键可 以取消当前选择的点,通过键盘上的Return键可以结束交互式的选择。
BW=imregionalmax(I):该函数获取灰度图像I的局部极大值,返回值BW 为和原图像大小相同的二值图像,BW中元素1对应的是极大值,其它元素 值为0。
BW=imregionalmax(I, conn):该函数中参数conn为连通类型。对于二维图 像conn可以取值为4和8,默认值为8。对于三维图像,conn可取值为6、18 和26,默认值为26。
I2=imfill(I):该函数对灰度图像进行填充操作,返回值I2也是灰度图像。
12.2.3 最大值和最小值
对于一幅图像可以有多个局部极大值或极小值,但只有一个最大值或最小 值。在MATLAB软件中,采用函数imregionalmax( )获取图像的所有局部极 大值,采用函数imregionalmin( )获取局部极小值。函数imregionalmax ( )的 调用格式为:
第12章 形态学图像处理
形态学是一种应用于图像处理和模式识别领域的新的方法,是一门建立在 严格的数学理论基础上而又密切联系实际的科学。由于形态学具有完备的 数学基础,这为形态学用于图像分析和处理等奠定了坚实的基础。本章详 细的介绍利用MATLAB软件进行形态学图像处理,主要内容包括基本运算、 组合形态学运算以及二值图像的其他形态学运算等内容。
12.1 基本的形态学运算
数学形态学可以看作是一种特殊的数字图象处理方法和理论,主要以图象 的形态特征为研究对象。它通过设计一整套运算、概念和算法,用以描述 图象的基本特征。这些数学工具不同于常用的频域或空域算法,而是建立 在微分几何以及随机集论的基础之上的。数学形态学作为一种用于数字图 象处理和识别的新理论和新方法,它的理论虽然很复杂,但它的基本思想 却是简单而完美的。
12.1.4 开运算和闭运算
结构元素对的开运算,记作,定义为:
即首先采用结构元素对做腐蚀运算,然后再做膨胀运算,使用相同的结构 元素。
闭运算是开运算的对偶运算,记作,定义为:
即首先用结构元素对做膨胀运算,然后再做腐蚀运算,使用相同的结构元 素。
12.2 组合形态学运算
图像的膨胀和腐蚀是最基本的形态学运算。下面介绍组合形态学运算,主 要内容包括高帽滤波和低帽滤波、边界提取、区域填充等等。
12.2.1 高帽滤波和低帽滤波
图像的形态学高帽滤波(top-hat filtering)定义为:
其中为输入的图像,为采用的结构元素,即从图像中减去形态学开操作后 的图像。通过高帽滤波可以增强图像的对比度。
图像的形态学低帽滤波(bottom-hat filtering)定义为:
其中为输入的图像,为采用的结构元素,即从图像中减去形态学闭操作后 的图像。通过低帽滤波可以获取图像的边缘。
12.1.1 基本概念
集合论是数学形态学的基础,首先对集合论的一些基本概念做一个简单介 绍。腐蚀运算和膨胀运算是数学形态学的两个基本变换。参加运算的对象 有两个:图像A(感兴趣目标)和结构集合B,B称为结构元素。结构元素 通常是个圆盘,但它其实可以是任何形状。
设和是的子集,则把图像沿矢量平移一段距离记作或,其定义为:
[BW2, locations]=imfill(BW):该函数中返回值locations包含了交互式选择 时的点的坐标。
BW2=imfill(BW, locations):该函数中通过参数locations指定了进行填充时 的点的坐标。
BW2=imfill(BW, ‘holes’):该函数通过参数holes可以填充二值图像中的 空洞。