数字图像处理实验报告_3.doc

《数字图像处理》作业报告 一．实验内容： ⑴ 图像的锐化：使用Sobel，Laplacian 算子分别对图像进行运算，观察并体会运算结果。 ⑵ 综合练习：对需要进行处理的图像分析，正确运用所学的知识，采用正确的步骤，对图像进行各类处理，以得到令人满意的图像效果。 二．实验目的： 深入理解空间域图像增强各类方法的原理，根据原理自行用Matlab软件编程（不调用原函数）实现对数字图像进行处理，感受各种不同的图像处理方法对最终图像效果的影响，最后能达到选择适合方法对图像处理或者使用混合方法对图像进行处理的要求。 三．实验步骤： 1.深入理解图像增强各个方法的原理，仔细阅读Matlab帮助文件中有关以上函数的使用说明，能充分理解其使用方法并能运用它们完成实验内容。 2.将homework2(skeleton_orig).bmp图像文件读入Matlab，对原图像使用拉普拉斯算子变换得到图(b)，再将图(b)与原图相加得到图(c)。(注：图像详见实验结果) 3.对原图像经过sobel算子锐化得到图(d)，将图(d)经过5*5均值滤波器平滑滤波之后得到图(e)。 4.将分别经过拉普拉斯算子锐化和sobel算子锐化后的图(c)与图(e)相乘得到掩蔽图像(f)，将原图(a)和掩蔽图像(f)相加得到(g)，最后对图(g)应用幂律变换得到最终图(h)。 四．实验结果： 4.1代码实现： （1）/*实现对图像的拉普拉斯算子的锐化*/ M文件1： 代码实现： （1）/*实现对图像的拉普拉斯算子的锐化*/ M文件1： function [r] = lap(x) [a,b]=size(x); r=x; for i=2:a-1 for j=2:b-1 r(i,j)=8*x(i,j)-(x(i+1,j)+x(i-1,j)+x(i,j+1)+x(i,j-1)+x(i-1,j-1)+x(i-1,j+1)+x(i+1,j-1)+x(i+1,j+1)); end end 调用指令： [cm,map]=imread(C:\Documents and Settings\admistator.AUTO-7FAB4782EA\桌面\数字图像处理\homework2(skeleton_orig).bmp); //读入给定图像 y=lap(cm); //对读入图像进行拉普拉斯变换 imshow(cm); //显示原图像(a) imshow(y); //显示原图像的拉普拉斯后图像(b) Imshow(cm+y); //显示对原图进行拉普拉斯算子锐化的图像(c) （2）/* 对原图像sobel处理并显示*/ M文件2： function [r] = sobel(x) [a,b]=size(x); r=double(x); r1=double(x); r2=double(x); for i=2:a-1 for j=2:b-1 r1(i,j)=abs(2*x(i+1,j)+x(i+1,j+1)+x(i+1,j-1)-x(i-1,j-1)-2*x(i-1,j)-x(i-1,j+1)); r2(i,j)=abs(2*x(i,j+1)+x(i+1,j+1)+x(i-1,j+1)-x(i-1,j-1)-2*x(i,j-1)-x(i+1,j-1)); r(i,j)=r1(i,j)+r2(i,j); end End 调用指令： I=imread(C:\Documents and Settings\admistator.AUTO-7FAB4782EA\桌面\数字图像处理\homework2(skeleton_orig).bmp); I1=I-uint8(sobel(I)); imshow(I1); //显示经sobel算子锐化的图像(d) （3）/* 对sobel图像进行5*5均值滤波平滑*/ M文件3： function [r] = average(x) [a,b]=size(x); r=x; r=double(r); for i=3:a-2 for j=3:b-2 r(i,j)=0; for s=-2:2