购买

¥15.0

加入VIP
  • 专属下载券
  • 上传内容扩展
  • 资料优先审核
  • 免费资料无限下载

上传资料

关闭

关闭

关闭

封号提示

内容

首页 matlab静态图像分割与边缘检测及图像压缩与编码

matlab静态图像分割与边缘检测及图像压缩与编码.doc

matlab静态图像分割与边缘检测及图像压缩与编码

豆浆
2019-05-18 0人阅读 举报 0 0 0 暂无简介

简介:本文档为《matlab静态图像分割与边缘检测及图像压缩与编码doc》,可适用于高等教育领域

学号光电图像处理实验报告实验三:静态图像分割与边缘检测作者肖剑洪专业电子科学与技术学院物理与电子学院指导老师王晓明完成时间实验三静态图像分割与边缘检测一、实验目的.学习常用的图像分割与边缘检测方法并通过实验使学生体会一些主要的分割算子对图像处理的效果以及各种因素对分割效果的影响.观察图像分割的结果产生对所讲述理论知识的直观认识加深对图像分割与边缘检测相关理论知识的理解。.掌握常用图象分割及边缘检测方法的算法设计及编程实现.学会使用MATLAB软件中关于图像分割与边缘检测的函数二、实验设备联想图像处理工作站三、实验内容及要求.自己编写Mfunction实现图像阈值分割算法要求该程序能对级灰度图像进行处理显示处理前、后图像.自己编写Mfunction实现利用Sobel算子进行图像边缘检测的算法并对图像进行检测显示原图像、处理后的图像。.调用Matlab自带的图像处理函数用不同的算子对图像进行分割、边缘检测比较结果。.结合以上实验内容使用ICETECKDMIDKM实验系统进行相应的动态视频图像分割及边缘检测观察结果。四、实验原理.图像分割图像分割是将图像划分成若干个互不相交的小区域的过程小区域是某种意义下具有共同属性的像素的连通集合。图像分割有三种不同的途径:区域法、边界法、边缘法。最常用的是灰度阈值化处理进行的图像分割:域值T的选取直接影响分割的效果!()直方图双峰域值选择()迭代域值选择迭代思想:选择一个初始估计值通过某种策略不断改进(调制)新的估计值直到满足给定的准则。迭代步骤:⑴选择一个初值估计值T一般为最大灰度值和最小灰度值的中间值⑵使用域值T分割图像得到两组像素G(>=T)和G(<T)⑶计算两组像素的灰度均值:μ1和μ⑷计算新域值T=(μ1和μ)⑸重复步骤~,直到新域值的改变量小于预先定义的参数e。()最大类间方差法(Otsu)思想:⑴将图像中目标和背景看作分属不同类别的像素组成⑵判别分析法的目标是确定域值使分属不同类别的像素的类间方差最大。.边缘检测灰度或结构等信息的突变出称为边缘在空间域借助微分算子通过卷积来完成而空域的微分在离散数字图像可以采用差分来近似。基于一阶导数的边缘检测算子有Robert算子、Sobel算子、Prewitt算子等基于二阶导数的边缘检测算子有Laplace算子LOG算子是一种改进的方式。Sobel算子为一对模板:五、实验步骤.根据实验内容的要求在MATLAB软件中编写相应程序)灰度阈值分割实验(使用迭代阈值选择方法))使用Sobel算子的边缘检测实验)调用Matlab已有的图像处理函数对图像进行分割和边缘检测a.调用边缘检测函数edge并与自己编写的函数运行结果进行比较b.用其他边缘检测算子和分割方法对图像进行处理并对结果进行比较。.调试运行程序并记录结果。.结合以上实验内容使用ICETECKDMIDKM实验系统进行相应的动态视频图像增强处理观察结果。.完成实验报告。六、实验程序清单和实验结果及分析.图像分割clearall,closeallH=imread('f:jpg')I=rgbgray(H)figure(),imshow(I)figure()imhist(I)T=Ibw=imbw(I,T)选择阈值T=对图像二值化figure()subplot(,,),imshow(Ibw)T=graythresh(I)采用Otsu方法计算最优阈值T对图像二值化L=uint(T*)Ibw=imbw(I,T)subplot(,,),imshow(Ibw)()迭代域值选择clcclearallH=imread('f:jpg')I=rgbgray(H)I=double(I)k=(max(max(I))min(min(I)))rowscols=size(I)count=count=fori=:rowsforj=:colsifI(i,j)<kcount=countG(count)I=I(i,j)elsecount=countG(count)I=I(i,j)endendendk=(mean(mean(GI))mean(mean(GI)))while(abs(kk)>())k=kcount=count=fori=:rowsforj=:colsifI(i,j)<kcount=countG(count)I=I(i,j)elsecount=countG(count)I=I(i,j)endendendk=(mean(mean(GI))mean(mean(GI)))endfigure()imshow(I)figure()II=imbw(I,k)imshow(II)边缘检测clearall,closeallH=imread('f:jpg')I=rgbgray(H)BW=edge(I,'sobel')BW=edge(I,'canny')BW=edge(I,'prewitt')BW=edge(I,'roberts')BW=edge(I,'log')figure(),imshow(I),title('OriginalImage')figure(),imshow(BW),title('sobel')figure(),imshow(BW),title('canny')figure(),imshow(BW),title('prewitt')figure(),imshow(BW),title('roberts')figure(),imshow(BW),title('log')八、实验思考题.小结一下本实验所用的边缘检测方法。答:边缘检测的基本算子有:一阶::RobertsCross算子Prewitt算子Sobel算子Kirsch算子罗盘算子二阶:MarrHildreth在梯度方向的二阶导数过零点Canny算子Laplacian算子Canny算子(或者这个算子的变体)是最常用的边缘检测方法。在Canny创造性的工作中他研究了设计一个用于边缘检测最优预平滑滤波器中的问题后来他说明这个滤波器能够很好地被一阶高斯导数核优化。另外Canny引入了非最大抑制概念它是说边缘定义为在梯度方向具有最大梯度值的点。在一个离散矩阵中非最大抑制阶梯能够通过一种方法来实现首先预测一阶导数方向、然后把它近似到度的倍数、最后在预测的梯度方向比较梯度幅度。.基于微分方法的边缘检测算法的依据是什么?答:一般来说不同区域的交界处构成边缘同一区域内像素的灰度差较小而不同区域像素灰度差较大因此用微分可以很好地检测满足这种假设的区域边界。学号光电图像处理实验报告实验四:图像压缩与编码作者肖剑洪专业电子科学与技术学院物理与电子学院指导老师王晓明完成时间实验四图像压缩与编码一、实验目的.了解图像压缩编码原理.掌握常用的图像压缩算法及编程实现.学会使用MATLAB软件中关于图像压缩编码的函数二、实验设备联想图像处理工作站三、实验内容及要求.自己编写Mfunction实现图像无损压缩算法要求该程序能对级灰度图像进行压缩数据计算压缩算法的性能.自己编写Mfunction实现图像有损压缩算法要求该程序能对级灰度图像进行压缩数据计算压缩算法的性能.调用Matlab自带的图像处理函数用不同的算法对图像进行压缩编码比较结果。四、实验原理图像压缩与编码是在满足一定保真度的要求下对图像数据的进行变换、编码和压缩去除冗余数据减少表示数字图像时需要的数据量以便于图像的存储和传输即以较少的数据量有损或无损地表示原来的像素矩阵的技术。图像压缩编码可分为两类:一类压缩是可逆的即从压缩后的数据可以完全恢复原来的图像信息没有损失称为无损压缩编码另一类压缩是不可逆的即从压缩后的数据无法完全恢复原来的图像信息有一定损失称为有损压缩编码。传统数据压缩方法的分类:无损压缩包括统计编码(Huffman编码Shannon编码游程编码算术编码等)和轮廓编码?有损压缩包括预测编码(脉冲编码调制PCM,Differential?PCM,AdaptiveDPCM等)、变换编码(DFTDCTKLTWHT小波变换等)和混合编码。??统计编码是根据信源的概率分布特性分配具有惟一可译性的可变长码字降低平均码字长度以提高信息的传输速度节省存储空间。其基本原理是在信号概率分布情况已知的基础上概率大的信号对应的码字短概率小的信号对应的码字长这样就降低了平均码字长度。其中Huffman编码具体的编码方法为?:①把输入元素按其出现概率的大小顺序排列起来然后把最末两个具有最小概率的元素之概率加起来?②把该概率之和同其余概率由大到小排队然后再把两个最小概率加起来再重新排队③重复②直到最后只剩下两个概率为止?。?变换编码的基本原理是通过数学变换可以改变信号能量的分布从而压缩信息量。以傅里叶变换的概念说明合理的变换可以改变信号能量分布的基本原理。在变换编码中有以下二个问题值得注意:图像变换方法的选取子图像大小的选取。传统的DFTDCTKLT和经典小波变换等变换编码?在图像变换后会产生浮点数因而必须对变换后的数据进行量化处理这样就产生不同程度的失真。而新一代的整数小波变换(第二代小波变换)采用提升方法能够实现整数变换因而能够实现图像的无损压缩。新的静态图像压缩标准JPEG中采用了基于提升方法的整数小波变换。?五、实验步骤.选择一组合适的原始图像作为测试图像.根据设计要求选择合适的图像压缩与编码方法.设计压缩与编码算法显示原始图像和压缩重构图像计算压缩的性能指标:压缩率和保真度.完成实验报告。六、实验结果及分析、利用余弦变换实现图像压缩:?DCT先将整体图像分成N×N像素块(一般N=?即个像素块)再对N×N块像素逐一进行DCT变换。由于大多数图像高频分量较小相应于图像高频成分的失真不太敏感可以用更粗的量化在保证所要求的图质下舍弃某些次要信息。实验代码:a=imread('f:jpg')I=rgbgray(a)imshow(I)title('原始图像')disp('原始图像大小:')whos('I')I=imdouble(I)图像类型存储转换将图像矩阵转换成双精度类型T=dctmtx()离散余弦变换矩阵B=blkproc(I,,'P*x*P',T,T')mask=B=blkproc(B,,'P*x',mask)I=blkproc(B,,'P*x*P',T',T)figure,imshow(I)title('压缩后的图像')disp('压缩图像的大小:')whos('I')原始图像大小:NameSizeBytesClassAttributesIxuint压缩图像的大小:NameSizeBytesClassAttributesIxdouble、利用小波变换实现图像压缩clearalla=imread('f:jpg')I=rgbgray(a)imshow(I)title('原始图像')disp('原始图像I的大小:')whos('I')I=imdouble(I)c,s=wavedec(I,,'bior')ca=appcoef(c,s,'bior',)ch=detcoef('h',c,s,)cv=detcoef('v',c,s,)cd=detcoef('d',c,s,)ca=appcoef(c,s,'bior',)ca=wcodemat(ca,,'mat',)ca=*cafigure,image(ca)title('第一次压缩后的图像')disp('第一次压缩图像的大小为:')whos('ca')ca=appcoef(c,s,'bior',)ca=wcodemat(ca,,'mat',)ca=*cafigure,image(ca)title('第二次压缩后的图像')disp('第二次压缩图像的大小为:')whos('ca')原始图像I的大小:NameSizeBytesClassAttributesIxuint第一次压缩图像的大小为:NameSizeBytesClassAttributescaxdouble第二次压缩图像的大小为:NameSizeBytesClassAttributescaxdouble八、实验思考题.DFT与DCT的性能比较答:DCT是离散傅里叶变换DFT是离散余弦变换。共同点:都将空域的图像数据信息转换到频域中即分离出图像的低频到高频成分。区别:、DCT与DFT转换后的域仅包含频域成分就叫频域DWT转换后的域不仅有频域成分还具有空域成分因此叫小波域。、DCT的频域的低频成分在DCT系数图的中间高频成分在四周离系数图中心越远频率越高DFT与DWT的频域的低频成分在系数图左上方越往右下方频率越unknown

用户评价(0)

关闭

新课改视野下建构高中语文教学实验成果报告(32KB)

抱歉,积分不足下载失败,请稍后再试!

提示

试读已结束,如需要继续阅读或者下载,敬请购买!

评分:

/14

VIP

在线
客服

免费
邮箱

爱问共享资料服务号

扫描关注领取更多福利

2017夜夜干精品 2017天天夜夜干精品 夜夜啪狠狠爱2017精品在线观看