图像锐化（三）

图像锐化的目的是加强图像中景物的边缘和轮廓，突出图像中的细节或增强被模糊了的细节。

方法三：高斯滤波

% 高斯曲线及其一阶二阶导数
clear;
a=0;sigma=2;
x=-10:0.0001:10;
figure(1)
y=(1/((sqrt(2*pi))*sigma))*exp(-((x-a).^2)/(2*sigma.^2));
plot(x,y,'b','LineWidth',1.5);
z = diff(y);
figure(2)
plot(z,'b','LineWidth',1.5);
z1 = diff(z);
figure(3)
plot(z1,'b','LineWidth',1.5);

从上图可以看出，总权值的95%包含在2sigma 的中间范围内，这个特性使得高斯函数常被用来作为权值；

 

Log 算子

即H为高斯函数的二阶导

clc;
clear;
src=imread('test.jpg');
src=rgb2gray(src);
src=im2double(src);
figure,subplot(231),imshow(src),title('orgin');

%----- log 算子 -----%
H1=fspecial('log', 7, 1);
R1=imfilter(src,H1);
edge1=abs(R1);
subplot(232),imshow(R1),title('noneedge1');
subplot(233),imshow(edge1),title('absedge1');
afterLog1=src+edge1;
subplot(234),imshow(afterLog1),title('afterLog1');

x=[-3 -2 -1 0 1 2 3];
y=[-3 -2 -1 0 1 2 3];
sigma=1;
for i=1:7
    for j=1:7
        H2(i,j)=1/(pi*sigma.^4)*((x(i)*x(i)+y(j)*y(j))/(2*sigma.^2)-1)*exp(1).^(-(x(i)*x(i)+y(j)*y(j))/(2*sigma.^2));
    end
end
R2=imfilter(src,H2);
edge2=abs(R2);
subplot(235),imshow(edge2),title('absedge2');
afterLog2=src+edge2;
subplot(236),imshow(afterLog2),title('afterLog2');