实验三图像分割
- 格式:doc
- 大小:280.00 KB
- 文档页数:8
实验报告课程名称数字图像处理导论专业班级_______________姓名 _______________学号_______________电气与信息学院与谐勤奋求就是创新一.实验目得1.理解图像分割得基本概念;2.理解图像边缘提取得基本概念;3.掌握进行边缘提取得基本方法;4.掌握用阈值法进行图像分割得基本方法.二。
实验内容1.分别用Roberts,Sobel与拉普拉斯高斯算子对图像进行边缘检测。
比较三种算子处理得不同之处;2.设计一个检测图1中边缘得程序,要求结果类似图2,并附原理说明。
3.任选一种阈值法进行图像分割、图1 图2三.实验具体实现1.分别用Roberts,Sobel与拉普拉斯高斯算子对图像进行边缘检测。
比较三种算子处理得不同之处;I=imread(’mri、tif');imshow(I)BW1=edge(I,’roberts’);figure ,imshow(BW1),title(’用Roberts算子’)BW2=edge(I,’sobel’);figure,imshow(BW2),title(’用Sobel算子 ')BW3=edge(I,’log’);figure,imshow(BW3),title(’用拉普拉斯高斯算子’)比较提取边缘得效果可以瞧出,sober算子就是一种微分算子,对边缘得定位较精确,但就是会漏去一些边缘细节.而Laplacian—Gaussian算子就是一种二阶边缘检测方法,它通过寻找图象灰度值中二阶过零点来检测边缘并将边缘提取出来,边缘得细节比较丰富。
通过比较可以瞧出Laplacian-Gaussian算子比sober算子边缘更完整,效果更好。
2.设计一个检测图1中边缘得程序,要求结果类似图2,并附原理说明.i=imread('m83、tif’);subplot(1,2,1);imhist(i);title('原始图像直方图');thread=130/255;subplot(1,2,2);i3=im2bw(i,thread);imshow(i3);title('分割结果’);3.任选一种阈值法进行图像分割、i=imread('trees、tif’);subplot(1,2,1);imhist(i);title('原始图像直方图’);thread=100/255;subplot(1,2,2);i3=im2bw(i,thread);imshow(i3);title('分割结果’)1、分别用Roberts,Sobel与拉普拉斯高斯算子对图像进行边缘检测。
图像分割处理实验报告1. 引言图像分割是计算机视觉中的重要任务之一,其目标是将图像划分成具有相似特征的子区域。
图像分割在很多应用领域中都有着广泛的应用,比如医学影像分析、目标检测和图像编辑等。
本实验旨在探索不同的图像分割算法,并比较它们在不同场景下的效果和性能。
2. 实验方法2.1 实验数据本实验选取了一组包含不同场景的图像作为实验数据集,包括自然景观、人物肖像和城市街景等。
每张图像的分辨率为500x500像素。
2.2 实验算法本实验使用了两种经典的图像分割算法进行比较,分别是基于阈值的分割和基于边缘的分割。
2.2.1 基于阈值的分割基于阈值的分割算法是一种简单而直观的方法,其原理是根据像素值的亮度信息将图像分割成不同的区域。
在本实验中,我们将图像的灰度值与一个事先设定的阈值进行比较,如果大于阈值则设为白色,否则设为黑色,从而得到分割后的图像。
2.2.2 基于边缘的分割基于边缘的分割算法利用图像中的边缘信息进行分割,其原理是检测图像中的边缘并将其作为分割的依据。
在本实验中,我们使用了Canny边缘检测算法来提取图像中的边缘信息,然后根据边缘的位置进行分割。
2.3 实验流程本实验的流程如下:1. 加载图像数据集;2. 对每张图像分别应用基于阈值的分割算法和基于边缘的分割算法;3. 计算分割结果和原始图像之间的相似度,使用结构相似性指标(SSIM)进行评估;4. 分析并比较两种算法在不同场景下的分割效果和性能。
3. 实验结果3.1 分割效果实验结果表明,基于阈值的分割算法在处理简单场景的图像时效果较好,可以比较准确地将图像分割为目标区域和背景。
然而,当图像的复杂度增加时,基于阈值的分割算法的效果明显下降,往往会产生较多的误分割。
相比之下,基于边缘的分割算法在处理复杂场景的图像时表现良好。
通过提取图像的边缘信息,该算法能够较准确地分割出图像中的目标区域,相比于基于阈值的分割算法,其产生的误分割较少。
3.2 性能评估通过计算分割结果和原始图像之间的SSIM指标,我们可以得到两种算法在不同场景下的性能评估。
图像分割实验报告《图像分割实验报告》摘要:图像分割是计算机视觉领域的重要研究方向,它在许多领域都有着重要的应用价值。
本实验旨在探究图像分割算法在不同场景下的表现,并对比不同算法的优缺点,为图像分割技术的进一步发展提供参考。
一、实验背景图像分割是指将图像划分成若干个具有独立语义的区域的过程。
图像分割技术在医学影像分析、自动驾驶、图像识别等领域都有着广泛的应用。
因此,对图像分割算法的研究和优化具有重要意义。
二、实验目的本实验旨在通过对比不同图像分割算法在不同场景下的表现,探究其优劣,并为图像分割技术的进一步发展提供参考。
三、实验内容1. 数据准备:收集不同场景下的图像数据,包括自然景观、医学影像、交通场景等。
2. 算法选择:选择常用的图像分割算法,如基于阈值的分割、边缘检测、区域生长等。
3. 实验设计:将不同算法应用于不同场景的图像数据上,对比它们的分割效果和计算速度。
4. 结果分析:对比不同算法的优缺点,并分析其适用场景和改进空间。
四、实验结果通过实验我们发现,在自然景观图像中,基于阈值的分割算法表现较好,能够有效地将图像分割成不同的颜色区域;而在医学影像中,边缘检测算法表现更为出色,能够准确地识别出器官的边缘;在交通场景中,区域生长算法表现较好,能够有效地区分不同的交通标志和车辆。
五、结论不同的图像分割算法在不同场景下有着不同的表现,没有一种算法能够适用于所有场景。
因此,我们需要根据具体的应用场景选择合适的图像分割算法,或者结合多种算法进行优化,以达到更好的分割效果。
六、展望未来,我们将继续探究图像分割算法的优化和改进,以适应不同场景下的需求。
同时,我们还将研究图像分割算法在深度学习和人工智能领域的应用,为图像分割技术的发展贡献力量。
通过本次实验,我们对图像分割算法有了更深入的了解,也为其在实际应用中的选择提供了一定的指导。
希望我们的研究能够为图像分割技术的发展做出一定的贡献。
图像分割实验报告医学图像处理实验报告实验名称:图像分割设计实验姓名:gaojunqiang学号:20105405班级:生医 1001指导教师:……2013年6月5日一、实验目的1、编程实现下列功能:读出存储的黑白灰度图象并显示,用拉普拉斯算子对图象进行边缘检测,显示处理后图象,存储处理后图象。
2、编程实现下列功能:读出存储的黑白灰度图象并显示,用鼠标点击图象上某一点,以灰度相近原则进行区域生长,显示处理后图象,存储处理后图象。
二、实验原理1、拉普拉斯边缘检测二维函数f(x,y)的拉普拉斯是一个二阶的微分定义为:,2f = [,2f / ,x2 ,,2f / ,y2]一般情况下可以用一个拉普拉斯模板算子。
模板算子分为4邻域和8邻域,如下,2f = 4z5 – (z2 + z4 + z6 + z8)0 1 01 -4 10 1 04邻域模板算子,2f = 8z5 – (z1 + z2 + z3 + z4+z5 + z6 + z7+ z8)1 1 11 -8 11 1 18邻域模板算子2、区域增长区域增长是通过一个起始点作为种子点对他周围的点进行选择。
它采取的是一种迭代的思想。
区域增长也分为四邻域和八邻域两种方式。
通过像素的集合进行区域增长的算法如下:1)根据图像的不同应用选择一个或一组种子,它或者是最亮或最暗的点,或者是位于点簇中心的点。
2)选择一个描述符(条件)3)从该种子开始向外扩张,首先把种子像素加入集合,然后不断将与集合中各个像素连通、且满足描述符的像素加入集合4)上一过程进行到不再有满足条件的新结点加入集合为止。
三、实验代码及结果1、拉普拉斯边缘检测代码如下:%主函数如下:clc;closeall;clearall;Imag = imread('lena.tiff');ImagGray = rgb2gray(Imag); %将彩色图片转换成灰度图片 figure()imshow(ImagGray),title('灰度图像');% T = graythresh(ImagGray); %用大津法自动确定阈值 %I=edge(ImagGray,'log',0.004);% figure(),imshow(I), title('laplace边缘图像'); ImagGray =double(ImagGray); T = 60; %设置阈值LapModType = 8; %设置laplace模板方式ImagLapEdge = LaplaceEdge(ImagGray,LapModType,T); %laplace边缘检测ImagLapEdge = uint8(ImagLapEdge);figure()imshow(ImagLapEdge),title('laplace边缘图像');%拉普拉斯边缘检测函数如下:functionImagLapEdge = LaplaceEdge(ImagGray,MoldType,Thresh)%-----------------参数介绍-------------------- %输入参数:% ImagGray: 输入的灰度图像% MoldType: 摸板类型,包括四邻域和八邻域 % Thresh: 边缘检测阈值%输出参数:% ImagEdge: 边缘像素点,存储的是二值化图像[r,c] = size(ImagGray); ImagLapEdge = zeros(r,c);%四邻域拉普拉斯边缘检测算子if 4 == MoldTypefori = 2:r-1for j = 2:c-1Temp =-4*ImagGray(i,j)+ImagGray(i-1,j)+ImagGray(i+1,j)+ImagGray(i,j-1)+ImagGray(i,j+1);if Temp > ThreshImagLapEdge(i,j) = 255; elseImagLapEdge(i,j) = 0; endendendend%八邻域拉普拉斯边缘检测算子if 8 == MoldTypefori = 2:r-1for j = 2:c-1Temp =-8*ImagGray(i,j)+ImagGray(i-1,j)+ImagGray(i+1,j)+ImagGray(i,j-1)+ImagGray(i,j+1)+ImagGray(i-1,j-1)+ImagGray(i+1,j+1)+ImagGray(i+1,j-1)+ImagGray(i-1,j+1);if Temp > ThreshImagLapEdge(i,j) = 255; elseImagLapEdge(i,j) = 0;endendendend拉普拉斯边缘检测实验结果如下:图1 原灰度图像2、区域增长实验代码:主函数如下:clc;closeall;clearall;Imag = imread('lena.jpg');figure()imshow(Imag), title('原灰度图片');n = 4; %设置区域增长的种子点数 [x,y] = ginput(n); %在图像上取点 V = [y,x];V = uint16(V);thresh = 33; %区域增长的阈值growingtype = 8; %增长方式[regionGet,imout,regionsize]=regiongrowing(V,Imag,thresh,growingtype ); %区域增长函数figure()imshow(imout), title('区域分割后的图片'); holdonplot(x,y,'+');由于区域增长代码老师已给这里就不在写出。
医学图像处理实验报告 ----图像分割医学图像处理实验报告----图像分割一.实验目的:掌握基本的图像分割方法,观察图像分割的结果,加深对边缘检测、模板匹配、区域生长的理解。
二.实验内容:边缘检测、模板匹配、区域生长。
三.实验方法:1.边缘检测:图象Blood边缘检测方法Sobel打开Toolboxes\Image Processing项选Edge Detection并运行选图象Blood边缘检测方法Sobel如图1所示按Apply键观察检测到的边界从上面四幅图像的对比来看,阈值逐渐变大,而满足要求的像素点也逐渐变少,使得图像的边缘提取的效果也越来越差,图像轮廓变得不清楚了。
以下为采用Prewitt方法的边缘提取效果:以下为Roberts方法边缘提取的效果:以下为Laplacian of Gaussian方法边缘提取的效果:以上的各种方法的理论算法有所不同,但总体效果基本一致。
以下是选其他图像重做上面的实验(适当简化)2.模板匹配:在Photoshop中打开一黑白灰度图象文件在滤镜菜单其他子菜单中选自定项在自定界面中输入点模板按好键观察处理后图象。
原始图像:点模板滤镜后的图像:0 0 00 1 00 0 0点模板: -1 -1 -1 -1 8 -1-1 -1 -1线模板: -1 -1 -1 2 2 2-1 -1 -1线模板: -1 2 -1 -1 2 -1-1 2 -1线模板: 2 -1 -1 -1 2 -1-1 -1 2线模板: -1 -1 2 -1 2 -12 -1 -1从上面的四种线模板得比较中可以发现:第一种对检测横向图像更为有效,第二种为竖向,后两种为135和45度。
这是与模板的构成有关的。
方向模板:-1 1 1-1 -2 1-1 1 1可以看出这个方向模板较多地体现出东方向的像素。
方向模板:1 1 -11 -2 -11 1 -1可以看出这个模板较多地体现出西方向的情况。
方向模板:-1 -1 -11 -2 11 1 1这个模板较多地体现了南向的情况。
实验三数字图像分割处理一、实验目的1.熟悉并掌握图像分割处理的方法以及编程实现方法。
2.熟悉图像分割处理的过程。
二、实验内容图像分割是一种重要的图像分析技术。
图像分割是从图像预处理到图像识别和分析理解的关键步骤,在图像处理中展重要位置,一方面它是目标表达的基础,对特征测量有重要的影响,另一方面,图像分割以及基于分割的目标表达、特征提取和参数测量等将原始图像转化为更为抽象更为紧凑的形式,使得到更高层的图像识别、分析和理解成为可能。
典型而传统的分割方法可分为基于阈值的方法、基于边缘的方法和基于区域的分割方法等。
三、实验原理1、灰度阈值法是把图像的灰度分成不同等级,然后用设置灰度阈值的方法确定有意义的区域或欲分割物体的边界。
阈值T一般可写成:T=T[x,y,p(x,y),q(x,y)] (1)其中p(x,y)代表点(x,y)处的灰度值,q(x,y)代表该点邻域的某种局部特性.根据(1)式可将阈值分为全局阈值、局部阈值和动态阈值.如果T的选取只与p(x,y)有关,则是全局阈值,全局阈值是利用全局信息(例如整个图像的灰度直方图)得到的阈值,它仅与全图各像素的本身性质有关,对全图使用;如果T的选取与p(x,y),q(x,y)有关,则是局部阈值,它与图像局部区域性质有关,全局阈值和局部阈值也称为固定阈值;如果T的选取不仅与p(x,y),q(x,y)有关,还与该点的坐标(x,y)有关,则是动态阈值.动态阈值的选取是将原始图像分成若干个子图,然后利用某种固定阈值选取法对每一个子图确定一个阈值,再通过对这些子图所得阈值的插值得到对图像中每个像素进行分割所需的阈值。
2、边缘检测法Robert算子是一种利用局部差分算子寻找边缘的算子:g(x,y)={[f(x,y)-f(x+1,y+1)]^2+[f(x+1,y)-f(x,y+1)]^2}^1/2Laplacian算子是常用的边缘检测算子,它是各向同性的二阶导数。
Canny算子是由Gaussian卷积得到的边缘图像用抑制非极大值点的算法对梯度图像作后续处理。
图像分割实验报告一、实验目的1、掌握图像分割的基本思想,了解其分割技术及其计算策略;2、学会从图像处理到分析的关键步骤,掌握图像分割过程;3、了解图像分割的意义,进一步加深对图像分析与理解;4、掌握基本分割方法:迭代分割与OTSU图像分割,并编程实现。
二、实验原理(一)迭代阈值分割选取的基本思路就是:首先根据图像中物体的灰度分布情况,选取一个近似阈值作为初始阈值,一个较好的方法就就是将图像的灰度均值作为初始阈值,然后通过分割图像与修改阈值的迭代过程获得认可的最佳阈值。
迭代式阈值选取过程可描述如下:1、计算初始化阈值g0=(g max+g min);22、根据g0,将图像分为两部分,分别计算灰度值期望,取其平均值为g1;3、如此反复迭代,当|g n-g n−1|足够小时,停止迭代,取T=g n即为最终阈值。
(二)OTSU图像分割(最大类间方差法)就是一种自适应的阈值确定的方法,就是按图像的灰度特性,将图像分成背景与目标两部分。
背景与目标之间的类间方差越大,说明构成图像的两部分的差别越大, 当部分目标错分为背景或部分背景错分为目标都会导致两部分差别变小。
因此,使类间方差最大的分割意味着错分概率最小。
以最佳门限将图像灰度直方图分割成两部分,使两部分类间方差取最大值,即分离性最大。
OTSU阈值选取过程可描述如下:1、记T为目标与背景的分割阈值,目标点数占图像比例为w1, 平均灰度为u1;背景点数占图像比例为w2,平均灰度为u1;2、图像的总平均灰度为:u=w1*u1+w2*u2;3、目标与背景图象的方差:g=w1*(u1-u)*(u1-u)+w1*(u2-u)*(u2-u)=w1*w2*(u1-u2)*(u1-u2);4、当方差g最大时,可以认为此时前景与背景差异最大,此时的灰度T就是最佳阈值。
二、实验内容1、利用C++编程实现迭代阈值图像分割算法;2、利用C++编程实现OTSU动态阈值图像分割算法。
三、实验框图四、实验结果与分析图1原图图2迭代图像分割图3OTSU图像分割实验分析:由图2可知:迭代所得的阈值分割的图象效果良好,能区分出图像的目标与背景的主要区域,但在图像的细微处还没有很好的区分度。
实验三图像分割
一、实验目的
1、了解图像分割的基本概念;
2、掌握阈值分割、边缘检测的基本分割方法;
3、对检测的目标图像分析其目标特征
二、实验内容
1、实验原理
阈值分割利用图像中要提取的目标物与其背景在灰度特性上的差异,把图像视为具有不同灰度级的两类区域的组合,选择一个合适的阈值,以确定图像中每个像素点应该属于目标区域还是背景区域,从而产生二值图像。
边缘检测是利用边缘灰度变化的一阶或二阶导数的特点,可以将边缘点检测出来。
常用梯度、roberts、sobel、prewitt等算子进行检测。
图1 两种边缘点附近灰度方向导数变化规律
2、MATLAB实现
(1)在处理图像直方图的工具箱中,核心函数为imhist,其语法为:
imhist(f,n)——直接显示
f为输入图像,h为其直方图,n是形成直方图的灰度级个数(默认256)。
(2)阈值分割
BW=im2bw( I,level) ——将灰度图像、RGB图像转换为二值图像
Level为阈值(0~1),当输入图像的亮度小于level时,输出0,大于时输出1。
或不用函数。
直接编程实现。
(3)边缘检测
函数edge提供了几个导数估计器。
该函数基本语法为:
[g,t]=edge(f,‘method’,parameters)
f——输入图像,g——输出图像,t——阈值。
‘method’是具体用到的检测方法(sobel、prewitt、roberts、log、zerocross、canny),parameters对应不同检测方法的参数。
sobel边缘检测器
[g,t] = edge(f,'sobel',T,'dir')
T:指定阈值,dir:检测边缘首选方向(horizontal、vertical、both)
g:检测到的逻辑图像,边缘位置为1,其余位置为0。
t可选,输出参数edge函数所用阈值T。
(一般t和T参数可以不用,dir默认为both)
prewitt边缘检测器
[g,t] = edge(f,'prewitt',T,'dir')
该函数参数与sobel相同。
prewitt检测器比sobel检测器在计算上简单一些,但比较容易产生一些噪声。
Robert边缘检测器
[g,t] = edge(f,'roberts',T,'dir')
该函数参数与sobel相同。
(3)二值图像特征量
a)bwperim(I,CONN)——提取对象边界像素
CONN为邻域类型,4代表四连接,8代表八连接
BW2=bwperim(BW1,CONN),返回二值图像BW1的边界图像BW2。
b)S=bwarea(BW)——返回图像面积
c)Eul=bweuler(BW,n)——返回图像的欧拉数
n为连接类型,4为四连接,8为八连接
三、实验步骤
1、读取图像rice.tif,分别用Roberts, Sobel、prewitt对图像进行边缘检测,显示
检测结果。
比较三种算子处理的不同之处。
clear;
close all;
I=imread('rice.tif');
subplot(2,2,1);
imshow(I);title('原图');
[g,t] = edge(I,'roberts');
subplot(2,2,2);
imshow(g);title('roberts对图像边缘检测结果');
[g,t] = edge(I,'sobel');
subplot(2,2,3);
imshow(g);title('sobel对图像边缘检测结果');
[g,t] = edge(I,'prewitt');
subplot(2,2,4);
imshow(g);title('prewitt对图像边缘检测结果');
Roberts检测水平和垂直边缘的效果好于斜向边缘,定位精度高,对噪声敏感,无法抑制噪声的影响。
Sobel对噪声具有平滑作用,提供较为精确的边缘方向信息,但边缘定位精度不够高。
Prewitt对噪声具有平滑作用,定位精度不够高。
prewitt检测器比sobel 检测器在计算上简单一些,但相对比较sobel容易产生一些噪声。
2、读取图像rice.tif,绘制直方图,双峰法阈值分割。
clear;
close all;
I=imread('rice.tif');
imshow(I);title('原图');
figure(2);
imhist(I,256);
title('直方图');
figure(3);
BW=im2bw(I,0.5);
imshow(BW);title('阈值分割二值图');
3、读取细胞图像‘cancer02.bmp’,根据相关阈值识别癌细胞。
具体过程:
读取图像(imread)——转换二值图像
提取边缘(bwperim)
计算周长L(边界像素的点数)
计算细胞高度H、宽度W、面积A(bwarea)
计算圆度C=4πA / L2、
矩形度R=A/(H*W)
伸长度E=min(H,W)/max(H,W)
设定阈值Ath1=10000 Ath2=50000
Cth=0.5 Rth=0.5 Eth=0.8
进行识别:若面积大于50000,则为可疑小细胞癌细胞。
若面积小于10000,则为正常细胞。
若面积在10000-50000之间,且C大于阈值0.5,R大于0.5,E大于阈值0.8,则为癌细胞。
相关函数:
Sum(BW)
若BW为图像矩阵,每一列的所有像素求和(或所有行求和),生成一维向量。
若BW为一维向量,则对所有元素求和。
在求周长、高度、宽度时可以用。
程序:
clear
close all
f=imread('cancer02.bmp');
subplot(2,1,1)
imshow(f);title('原图像');
BW=bwperim(f,4);
subplot(2,1,2)
imshow(BW);title('图像边缘');
L=sum(sum(BW))
H=max(sum(f)/255)
f1=rot90(f);
W=max(sum(f1)/255)
A=bwarea(f)
C=4*pi*A / L^2
R=A/(H*W)
E=min(H,W)/max(H,W)
ATH1=10000;ATH2=50000;Cth=0.5;Rth=0.5;Eth=0.8;
if A>ATH2
'可疑小细胞癌细胞'
elseif A<ATH1
'正常细胞'
else A>ATH1 & A<A TH2 & Cth>0.5 & Rth>0.5 & Eth>0.8 '癌细胞'
end
end
end
运行结果:
L =595
H =204
W =203
A =3.3241e+004
C =1.1799 R =0.8027 E =0.9951 ans =癌细胞。