% 实验4 空域图像增强 图像锐化
% 用Laplacian锐化算子(分k=1和k=2两种情况)对256*256大小、256级灰度的数字图像lena进行锐化处理,显示处理前、后图像。
clc
clear
close all
I = imread('kg.jpg');
% figure,imshow(I);
[h w] = size(I);
I=double(I);
II=I;
% 用Laplacian锐化算子(分k=1和k=2两种情况)对256*256大小、256级灰度的数字图像lena进行锐化处理,显示处理前、后图像。
for i=2:h-1
for j=2:w-1
II(i,j)=I(i-1,j)+I(i+1,j)+I(i,j-1)+I(i,j+1)-4*I(i,j);
end
end
% 处理得到g11(k=1) g12(k=2)
k=1;
g11=I-k*II;
% figure,imshow(uint8(g11));
k=2;
g12=I-k*II;
% figure,imshow(uint8(g12));
figure
subplot(1,2,1),imshow(uint8(g11));
subplot(1,2,2),imshow(uint8(g12));
% g21(k=1)边缘 g22(k=2)边缘 显示比较
g21=I;
k=1;
for i=2:h-1
for j=2:w-1
g21(i,j)=4*k*I(i,j)-k*(I(i-1,j)+I(i+1,j)+I(i,j-1)+I(i,j+1));
% g21(i,j)=k*(I(i-1,j)+I(i+1,j)+I(i,j-1)+I(i,j+1))-4*k*I(i,j);
end
end
g22=I;
k=2;
for i=2:h-1
for j=2:w-1
g22(i,j)=4*k*I(i,j)-k*(I(i-1,j)+I(i+1,j)+I(i,j-1)+I(i,j+1));
% g22(i,j)=k*(I(i-1,j)+I(i+1,j)+I(i,j-1)+I(i,j+1))-4*k*I(i,j);
end
end
figure
subplot(1,2,1),imshow(uint8(g21));
subplot(1,2,2),imshow(uint8(g22));