matlab自编sobel算子,prewitt算子,log算子锐化图像.doc

课程名称 实验项目名称 图象的锐化 学生姓名 专业班级 学号 实验成绩 指导老师（签名 ） 日期 一. 实验目的和要求 二. 　I=imread(coins.png);%读入图像　　subplot(2,2,1),imshow(I);title(原图像);%显示原图像　　H=fspecial(sobel); %应用sobel算子锐化图像　　I2=filter2(H,I); %sobel算子滤波锐化　　subplot(2,2,2);imshow(I2); %显示sobel算子锐化图像　　title(sobel算子锐化图像);　　H=fspecial(prewitt);%应用prewitt算子锐化图像　　I3=filter2(H,I);%prewitt算子滤波锐化　　subplot(2,2,3);imshow(I3); %显示prewitt算子锐化图像　　title(prewitt算子锐化图像);　　H=fspecial(log); %应用log算子锐化图像　　I4=filter2(H,I); %log算子滤波锐化　　subplot(2,2,4);imshow(I4);%显示log算子锐化图像　　title(log算子锐化图像);A=imread(baboon.bmp); A=rgb2gray(A); subplot(2,2,1),imshow(A),title(原始图象); image=double(A); [m,n]=size(A); u=zeros(1,9); k=zeros(1,9); for i=2:m-1, for j=2:n-1, u(1)=0*image(i,j);u(2)=2*image(i,j+1);u(3)=1*image(i-1,j+1); u(4)=0*image(i-1,j);u(5)=-1*image(i-1,j-1);u(6)=-2*image(i,j-1); u(7)=-1*image(i+1,j-1);u(8)=0*image(i+1,j);u(9)=1*image(i+1,j+1); rimage1(i,j)=abs(sum(u)); k(1)=0*image(i,j);k(2)=0*image(i,j+1);k(3)=1*image(i-1,j+1); k(4)=2*image(i-1,j);k(5)=1*image(i-1,j-1);k(6)=0*image(i,j-1); k(7)=-1*image(i+1,j-1);k(8)=-2*image(i+1,j);k(9)=-1*image(i+1,j+1); rimage2(i,j)=abs(sum(k)); xuhua(i,j)=rimage1(i,j)+rimage2(i,j); end, end, subplot(2,2,2), imshow(xuhua,map), title(Sobel算子锐化图象); for i=2:m-1, for j=2:n-1, u(1)=0*image(i,j);u(2)=1*image(i,j+1);u(3)=1*image(i-1,j+1); u(4)=0*image(i-1,j);u(5)=-1*image(i-1,j-1);u(6)=-1*image(i,j-1); u(7)=-1*image(i+1,j-1);u(8)=0*image(i+1,j);u(9)=1*image(i+1,j+1); rimage1(i,j)=abs(sum(u)); k(1)=0*image(i,j);k(2)=0*image(i,j+1);k(3)=1*image(i-1,j+1); k(4)=1*image(i-1,j);k(5)=1*image(i-1,j-1);k(6)=0*image(i,j-1); k(7)=-1*image(i+1,j-1);k(8)=-1*image(i+1,j);k(9)=-1*image(i+1,j+1); rimag