实验六 利用Matlab进行图象处理

《实验六 利用Matlab进行图象处理》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数字图象处理实验指导书实验六利用Matlab进行图象处理实验六利用Matlab进行图象处理一、实验条件PC机Matlab软件（已安装图象处理等工具箱）大量样图（软件已提供）二、实验目的1、熟悉和理解使用Matlab图象处理工具箱函数进行简单图象处理的方法；2、观察图象处理效果，复习课堂教学内容，加深对图象处理原理和算法的理解；3、分析软件提供的实例（DEMOS），了解工程实际中图象处理的一般步骤。三、实验原理利用Matlab图象处理工具箱函数进行图象处理操作，相关图象处理原理及算法参见教材。四、实验内容1、以灰度图象rice.tif为例，说明利用Matlab图象处理工具箱函数进行图象处

2、理的一般步骤（详细步骤见后面“实验步骤”），此例涉及的图象处理操作有：图象读入，图象存储，图象显示，图象直方图及均衡，图象频谱分析，图象形态学分析，图象代数运算，图象灰度增强，图象二值化及分析等；2、观察并分析Matlab软件提供的一些图象处理实例（DEMOS）；3、对学有余力的同学，可以操作Matlab提供的其他高级应用，如小波变换等。五、实验步骤第1部分内容：利用Matlab对灰度图象rice.tif进行处理，处理的主要目的是消除灰度图象中亮度不一致的背景，并使用阈值将修改后的图象转换为二值图象，使用成员标记方法计算二值图象中的对象个数并分析相关统计特性。具体步骤如下：步骤1：双击

3、桌面上的Matlab图标，进入Matlab集成环境，在命令窗口依次键入下述命令；步骤2：读取及显示灰度图象。清除工作平台的所有变量，关闭已打开的图象窗口，读取及显示灰度图象rice.tif：%清除工作平台的所有变量，关闭已打开的图象窗口clear,closeall;%读取及显示灰度图象I=imread('rice.tif');数字图象处理实验指导书实验六利用Matlab进行图象处理imshow(I),title('原始灰度图象');%检查内存中的图象数据whos;步骤2步骤3：对灰度图象进行频谱分析。%灰度图象的傅立叶频谱图F=fft2(I,256,256);%对图象进行傅立叶变换F2

4、=fftshift(F);%将频谱坐标原点移至图形的中心figure,imshow(log(abs(F2)),[]),title('图象的傅立叶频谱图');%%灰度图象的DCT频谱图D=dct2(I,256,256);%对图象进行DCT变换figure,imshow(log(abs(D)),[]),title('图象的DCT频谱图');步骤3步骤4：灰度图象直方图及均衡。%显示图象的直方图figure,imhist(I),title('灰度图象的直方图');E=histeq(I);%对图象进行直方图均衡数字图象处理实验指导书实验六利用Matlab进行图象处理figure,imhist(

5、E),title('灰度图象均衡后的直方图');步骤4步骤5：图象存储及图象文件信息查询。%将均衡后的图象保存为jpg格式文件imwrite(E,'rice.jpg');%查询rice.jpg图象文件的信息imfinfo('rice.jpg')%注意后面不加“；”步骤6：进行灰度图象形态学（开启）操作，提取图象背景。%估计灰度图象的背景亮度并显示background=imopen(I,strel('disk',15));figure,imshow(background),title('灰度图象的背景亮度估计');步骤7：图象代数运算（图象相减）。%从原始图象中减去背景图象并显示I2=i

6、msubtract(I,background);figure,imshow(I2),title('原始灰度图象减去背景后的图象');步骤6步骤7步骤8：调节对比度，对图象进行灰度增强。%调节图象对比度I3=imadjust(I2,stretchlim(I2),[01]);数字图象处理实验指导书实验六利用Matlab进行图象处理figure,imshow(I3),title('调节图象对比度后的图象');步骤9：对灰度图象进行二值化。%使用阈值操作将灰度图象转化为二值图象level=graythresh(I3)bw=im2bw(I3,level);figure,imshow(bw),ti

7、tle('使用阈值操作将灰度图象转化为二值图象');步骤8步骤9步骤10：对二值图象进行连接成份标记，查看图象中的对象书目。%检查二值图象中的对象个数[labeled,numobjects]=bwlabel(bw,4)%bwlabel为标记函数步骤11：将图象标记矩阵转化成伪彩色图象，便于观察图象中的各个对象。%将标记矩阵转化成伪彩色图象rgb_label=label2rgb(labeled,'spring','c','shuffle');fi