《数字图像处理》实验指导书2013

《数字图像处理》实验指导书山东建筑大学信息与电气工程学院2013年3月目录实验守则……………………………………………………………………1实验一图像增强与平滑…………………………………………………2实验二图像分割与边缘检测…………………………………………5实验三图像的几何变换………………………………………………7实验四图像形态学处理………………………………………………9实验五频域变换………………………………………………………12实验六数字图像处理应用………………………………………14附录1MATLAB简介……………………………………………………16附录2实验报告范例---LOG滤波器提取边缘………………………171实验守则为了加强实验室管理，提高实验教学质量，培养学生理论联系实际的学风，提高实际操作技能，确保人身和设备的安全，特制定如下守则：1．实验前要认真预习，了解实验目的、实验原理、实验线路、实验步骤及欲测的物理量。熟悉和实验有关的仪器设备的用途、使用方法及注意事项。能回答指导教师的提问。2．每个实验小组要选出组长1名，实验时要认真分工，轮流当记录员和操作者，以便每个人都能得到较全面的训练。实验线路接好后，各组学生要自行检查，然后请指导教师检查，经检查无误后方可接通电源。改接线路必须先断开电源，以免触电及损坏仪表。改接后的线路仍需经教师检查。3．实验时应听从教师指导，严格按实验步骤进行，如实记录实验数据，积极思考和分析实验中发生的现象。实验中如遇事故应立即关断电源、保持现场，报告指导教师及时处理。如有损坏仪器或器件，要填写事故原因及报损单。4．实验完毕应立即关断电源，把实验记录交指导教师检查合格后方可拆除线路。如发现实验记录有不当和错误之外，应立即重作。要把所有仪器、仪表、工具、导线整理好。经指导教师清点后方可离开实验室。5．要爱护公物。实验室所有的仪表、工具、器件等一律不得擅自带出。实验中如需更换仪器、仪表、工具、器件等，必须经指导教师同意。不得擅自更换或乱拿他组物品。未经允许，不得擅自拆卸实验用的仪器、设备。6．实验时要严肃认真，讲究文明、礼貌，要保持安静和室内整洁，不随地吐痰和乱抛杂物。实验完毕要打扫实验室卫生。7．书写实验报告时，语言要力求简练，书写应端正，作图应正规。实验报告一般应包括以下内容：(1)实验目的；所用仪器、仪表、器件的名称及规格型号；实验电路原理图。(2)实验项目(指导书中的详细内容不必抄录)；测试记录表；波形图及现象记录。(3)实验数据整理，计算后的结果、及根据实验数据绘制的曲线。(4)实验分析，应包括与理论计算值的比较、现象分析。8．关于实验考核：学生在实验课中的实际技能水平与实验报告的质量是评定本课程实验成绩的标准。因各种原因缺做实验和缺交实验报告者，实验成绩要受影响。2实验一图像增强与平滑一．实验目的及要求1．了解MATLAB的操作环境和基本功能。2．掌握MATLAB中图像增强与平滑的函数的使用方法。3．加深理解图像增强与平滑的算法原理。二、实验内容（一）研究以下程序，分析程序功能；输入执行各命令行，认真观察命令执行的结果。熟悉程序中所使用函数的调用方法，改变有关参数，观察试验结果。（可将每段程序保存为一个.m文件）1．直方图均衡化clearall;closeall%CleartheMATLABworkspaceofanyvariables%andcloseopenfigurewindows.I=imread('pout.tif');%Readsthesampleimages‘pout.tif’,andstoresitinimshow(I)%anarraynamedI.displaytheimagefigure,imhist(I)%Createahistogramoftheimageanddisplayitin%anewfigurewindow.[I2,T]=histeq(I);%Histogramequalization.figure,imshow(I2)%Displaythenewequalizedimage,I2,inanewfigurewindow.figure,imhist(I2)%CreateahistogramoftheequalizedimageI2.figure,plot((0:255)/255,T);%plotthetransformationcurve.imwrite(I2,'pout2.png');%WritethenewlyadjustedimageI2toadiskfilenamed%‘pout2.png’.imfinfo('pout2.png')%Checkthecontentsofthenewlywrittenfile2．直接灰度变换clearall;closeallI=imread('cameraman.tif');J=imadjust(I,[00.2],[0.51]);imshow(I)figure,imshow(J)[X,map]=imread('forest.tif');figure,imshow(X,map)I2=ind2gray(X,map);注意：imadjust()功能：调整图像灰度值或颜色映像表，也可实现伽马校正。语法：J=imadjust(I,[low_inhigh_in],[low_outhigh_out],gamma)newmap=imadjust(map,[low_inhigh_in],[low_outhigh_out],gamma)RGB2=imadjust(RGB1,...)3J2=imadjust(I2,[],[],0.5);figure,imshow(I2)figure,imshow(J2)J3=imadjust(I2,[],[],1.5);figure,imshow(J3)helpimadjust%Displaytheimadjust()functioninformation.3．空域平滑滤波（模糊、去噪）clearall;closeallI=imread('eight.tif');h1=ones(3,3)/9;h2=ones(5,5)/25;I1=imfilter(I,h1);I2=imfilter(I,h2);figure(1),imshow(I),title('OriginalImage');figure(2),imshow(I1),title('FilteredImageWith3*3')figure(3),imshow(I2),title('FilteredImageWith5*5')%加入Gaussian噪声J1=imnoise(I,'gaussian',0,0.005);%加入椒盐噪声J2=imnoise(I,'salt&pepper',0.02);%对J1、J2进行平均值平滑滤波K1=imfilter(J1,fspecial('average',3));K2=imfilter(J2,fspecial('average',3));figure(4);subplot(2,2,1),imshow(J1),title('gaussian');subplot(2,2,2),imshow(J2),title('salt&pepper');subplot(2,2,3),imshow(K1),title('average');subplot(2,2,4),imshow(K2);%对J1、J2进行中值滤波K3=medfilt2(J1,[33]);K4=medfilt2(J2,[33]);figure(5);subplot(2,2,1),imshow(J1),title('gaussian');subplot(2,2,2),imshow(J2),title('salt&pepper');subplot(2,2,3),imshow(K3),title('Medianfiltering');4subplot(2,2,4),imshow(K4)4．空域锐化滤波clearall;closeallI=imread('moon.tif');w=fspecial('laplacian',0)w8=[1,1,1;1,-8,1;1,1,1]I1=imfilter(I,w,'replicate');figure(1);imshow(I),title('OriginalImage');figure(2),imshow(I1),title('LaplacianImage');f=im2double(I);f1=imfilter(f,w,'replicate');figure(3),imshow(f1,[]),title('LaplacianImage');f2=imfilter(f,w8,'replicate');f4=f-f1;f8=f-f2;figure(4),imshow(f4);figure(5),imshow(f8);（二）采用MATLAB底层函数编程实现1．灰度变换之动态范围扩展假定原图像f(x,y)的灰度范围为［a,b］，希望变换后图像g(x,y)的灰度范围扩展至［c,d］，则线性变换可表示为：cayxfabcdyxg]),([),(用MATLAB底层函数编程实现上述变换函数。观察图像‘pout.tif’的灰度直方图，选择合适的参数［a,b］、［c,d］对图像‘pout.tif’进行灰度变换，以获得满意的视觉效果。2．非锐化掩蔽和高升滤波从原图像中减去其非锐化（平滑过的）图像的过程称为非锐化掩蔽，其基本步骤为：⑴对原图像进行平滑滤波得到模糊图像(,)fxy；⑵从原图像中减去模糊图像，产生的差值图像称为模板(,)maskgxy；⑶将模板加到原图像上，得到锐化后的图像(,)gxy。即，(,)(,)-(,)maskgxyfxyfxy5(,)(,)(,)1maskgxyfxykgxyk；用MATLAB函数编程实现上述功能。三、实验设备1．PIII以上微机；2．MATLAB6.5及以上；四、预习与思考1．预习实验内容，阅读教材熟悉实验原理；2．查阅资料，熟悉MATLAB的操作环境和基本功能。熟悉实验中涉及的有关函数。3．利用课余时间，用MATLAB底层函数编程实现实验内容（二）中的灰度线性变换。4．你能否给出实现样例程序功能的其它方法？五、实验报告要求1.简述试验的目的和试验原理；2.叙述各段程序功能，改变有关函数的参数，分析比较实验结果；3.打印出所编写的实验程序。4.写出本实验的心得体会及意见。6实验二图像分割与边缘检测一．实验目的及要求1．利用MATLAB研究图像分割与边缘检测的常用算法原理；2．掌握MATLAB图像域值分割与边缘检测函数的使用方法；3．了解边缘检测的算法和用途，比较Sobel、Prewitt、Canny等算子边缘检测的差异。二、实验内容（一）研究以下程序，分析程序功能；输入执行各命令行，认真观察命令执行的结果。熟悉程序中所使用函数的调用方法，改变有关参数，观察试验结果。1．图像阈值分割clearall,closeall;I=imread('rice.tif');figure(1),imshow(I)figure(2);imhist(I)T=120/255;Ibw1=im2bw(I,T);%选择阈值T=120/255对图像二值化；figure(3);subplot(1,2,1),imshow(Ibw1);T=graythresh(I);%采用Otsu方法计算最优阈值T对图像二值化；L=uint8(T*255)Ibw2=im2bw(I,T);subplot(1,2,2),imshow(Ibw2);helpim2bw;helpgraythresh;（令Ｔ取不同值，重做上述试验，观察试验结果）2．边缘检测clearall,closeall;I=imread('rice.tif');BW1=edge(I,'sobel');BW2=edge(I,'canny');BW3=edge(I,'prewitt