图像锐化matlab算法

%常用图像锐化算法

clc

clear

close all

ima=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',[9 9]);%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算子锐化');