数字图像处理实验之拉普拉斯算子锐化图像

代码如下：
主函数：

% 输入图像
img_1 = imread('Fig0338(a)(blurry_moon).tif');

% 两种拉普拉斯模板
model_1 = [0, 1, 0; 1, -4, 1; 0, 1, 0];
model_2 = [1, 1, 1; 1, -8, 1; 1, 1, 1];

% 调用函数
[img_2, img_3] = m_13_Sharpening_Laplace(img_1, model_1);
[img_6, ~] = m_13_Sharpening_Laplace(img_1, model_2);

% 使用第一种拉普拉斯模板锐化
img_4 = img_1 - img_2;
% 使用第二种拉普拉斯模板锐化
img_5 = img_1 - img_6;

% 显示结果
subplot(151), imshow(img_1);
subplot(152), imshow(img_2);
subplot(153), imshow(img_3);
subplot(154), imshow(img_4);
subplot(155), imshow(img_5);

m_13_Sharpening_Laplace.m 代码如下：

function [img_2, img_3] = m_13_Sharpening_Laplace(img_1, model)
% 该函数返回使用指定拉普拉斯模板滤波后的结果
% img_2: 未标定的结果，即负值直接定为0
% img_3: 标定后的结果
% 解释一下标定。例如结果中灰度级的范围为-100到100
% 如果不标定，那么直接将所有负数直接指定为0，然后对比度拉伸，即从0到100拉伸到0到255。
% 如果标定，那么就是将-100到100拉伸到0到255，-100对应0，100对应255。

[r, c] = size(img_1);
[~, m] = size(model);
% 不解释，看我之前博客
img_1s = m_10_Smooth_create_img(img_1, r, c, m - 1);

% 创建两个新图像
img_2 = zeros(r, c);
img_3 = zeros(r, c);

% 逐个像素处理
for x = 1 : r
    for y = 1 : c
        % 滤波器范围内所有像素点之和
        sum = 0;
        for i = 1:m
            for j = 1:m
                sum = sum + model(i, j) * img_1s(x + i - 1, y + j - 1);
            end
        end
        % 不标定，负数直接指定为0
        if sum <= 0
            img_2(x, y) = 0;
        else
            img_2(x, y) = sum;
        end
        % 标定，不处理
        img_3(x, y) = sum;
    end
end

% 使用下面两个函数进行标定，看我前面博客
img_2 = im2uint8(mat2gray(img_2));
img_3 = im2uint8(mat2gray(img_3));
end

结果如下：