数字图像处理实验报告-3

《数字图像处理》作业报告一．实验内容：⑴图像的锐化：使用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;fori=2:a-1forj=2:b-1r(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));endend调用指令：[cm,map]=imread('C:\DocumentsandSettings\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);fori=2:a-1forj=2:b-1r1(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);endEnd调用指令：I=imread('C:\DocumentsandSettings\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);fori=3:a-2forj=3:b-2r(i,j)=0;fors=-2:2fort=-2:2r(i,j)=r(i,j)+x(i+s,j+t)/16;endendr(i,j)=r(i,j)*16/25;endend调用指令：I=imread('C:\DocumentsandSettings\admistator.AUTO-7FAB4782EA\桌面\数字图像处理\homework2(skeleton_orig).bmp');I1=I-uint8(sobel(I));I2=average(I1);//显示图(d)经5*5均值滤波器平滑滤波之后图像(e)imshow(uint(I2))（4）/*由(c)和(e)相乘形成掩蔽图像(f)，由(a)和(f)求和得到锐化图像(g)，最后对锐化图像（g）应用幂率变换（c=1，r=0.5）*/调用指令：I=imread('C:\DocumentsandSettings\admistator.AUTO-7FAB4782EA\桌面\数字图像处理\homework2(skeleton_orig).bmp');I1=uint8(I-lap(I));I2=I-uint8(sobel(I));I3=uint8(average(I2));I1=double(I1);I3=double(I3);I4=I1.*I3/128;imshow(uint8(I4))//显示掩蔽图像(f)I5=I+uint8(I4);imshow(I5)//显示锐化图像(g)I5=double(I5);I6=sqrt(I5/256);imshow(uint8(I6*256))//显示最后结果图(h)4.2实验结果(图(a)~图(h))（b）图(a)的拉普拉斯变换（b）图(a)的拉普拉斯变换（a）全身骨骼扫描图像（d）图(a)的sobel处理（d）图(a)的sobel处理（c）图(a)和