%常用图像锐化算法clcclearcloseallima=imread('132.jpg');%读入图像ima=rgb2gray(ima);%转为灰度图像ima=double(ima);bw1=edge(ima,'sobel');%sobel算子锐化figure;subplot(121);imshow(uint8(ima));title('原始图像');%图像显示subplot(122);imshow(bw1);title('sobel算子锐化');bw2=edge(ima,'prewitt');%prewitt算子锐化figure;subplot(121);imshow(uint8(ima));title('原始图像');subplot(122);imshow(bw2);title('prewitt算子锐化');bw3=edge(ima,'roberts');%roberts算子锐化figure;subplot(121);imshow(uint8(ima));title('原始图像');subplot(122);imshow(bw3);title('roberts算子锐化');bw4=edge(ima,'log');%log算子锐化figure;subplot(121);imshow(uint8(ima));title('原始图像');subplot(122);imshow(bw4);title('log算子锐化');bw5=edge(ima,'canny');%nny算子锐化figure;subplot(121);imshow(uint8(ima));title('原始图像');subplot(122);imshow(bw5);title('canny算子锐化');h1=fspecial('gaussian',[99]);%gaussian低通滤波器锐化bw6=imfilter(ima,h1);figure;subplot(121);imshow(uint8(ima));title('原始图像');subplot(122);imshow(uint8(bw6));title('gaussian低通滤波器');h2=fspecial('laplacian');%laplacian算子锐化bw7=imfilter(ima,h1);figure;subplot(121);imshow(uint8(ima));title('原始图像');subplot(122);imshow(uint8(bw7));title('laplacian算子锐化');