图像锐化MATLAB实现（包含laplacian算子、Prewitt算子、Robert算子、sobel算子）

image=imread('tea.png'); I=rgb2gray(image); [m,n]=size(I); H1=I;H2=I;H3=I;H4=I; model1=[0,-1,0; -1,5,-1; 0,-1,0]; for i=2:m-1 for j=2:n-1 H1(i,j)=-I(i+1,j)-I(i-1,j)-I(i,j+1)-I(i,j-1)+4*I(i,j); end end model2=[-1,-1,-1; -1,9,-1; -1,-1,-1]; for i=2:m-1 for j=2:n-1 H2(i,j)=-I(i+1,j)-I(i-1,j)-I(i,j+1)-I(i,j-1)-I(i-1,j-1)-I(i-1,j+1)-I(i+1,j-1)-I(i+1,j+1)+8*I(i,j); end end model3=[1,-2,1; -2,5,-2; 1,-2,1]; for i=2:m-1 for j=2:n-1 H3(i,j)=-2*I(i+1,j)-2*I(i-1,j)-2*I(i,j+1)-2*I(i,j-1)+I(i-1,j-1)+I(i-1,j+1)+I(i+1,j-1)+I(i+1,j+1)+4*I(i,j); end end model4=[0,-1,0; -1,5,-1; 0,-1,0]; for i=2:m-1 for j=2:n-1 H4(i,j)=-I(i+1,j)-I(i-1,j)-I(i,j+1)-I(i,j-1)+5*I(i,j); end end subplot(3,2,1),imshow(image),title('原图'); subplot(3,2,2),imshow(I),title('原图灰度图'); subplot(3,2,3),imshow(H1),title('拉普拉斯H1变换图'); subplot(3,2,4),imshow(H2),title('拉普拉斯H2变换图'); subplot(3,2,5),imshow(H3),title('拉普拉斯H3变换图'); subplot(3,2,6),imshow(H4),title('拉普拉斯H4变换图');