欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

数字图像处理实验之梯度增强边缘

程序员文章站 2022-07-14 21:42:24
...

只提供代码,不提供思路,思路看课本。
主函数:

% 读入图像
img_1 = imread('Fig0342(a)(contact_lens_original).tif');
% Sobel模板
model_1 = [-1, -2, -1; 0, 0, 0; 1, 2, 1];
model_2 = [-1, 0, 1; -2, 0, 2; -1, 0, 1];

% 调用函数处理
img_2 = m_15_Gradient_Enhancement_process(img_1, model_1, model_2);

% 显示结果
subplot(121), imshow(img_1);
subplot(122), imshow(img_2);

m_15_Gradient_Enhancement_process.m 代码如下:

function img_2 = m_15_Gradient_Enhancement_process(img_1, model_1, model_2)

% 信息获取
[r, c] = size(img_1);
[~, m] = size(model_1);
% 获得填充边后的图像
img_1s = m_10_Smooth_create_img(img_1, r, c, m - 1);
% 创建新图像
img_2 = zeros(r, c);

for x = 1 : r
    for y = 1 : c
        % 求 g_x 和 g_y
        sum_1 = 0;
        sum_2 = 0;
        for i = 1 : m
            for j = 1 : m
                sum_1 = sum_1 + model_1(i, j) * img_1s(x + i - 1, y + j - 1);
                sum_2 = sum_2 + model_2(i, j) * img_1s(x + i - 1, y + j - 1);
            end
            % 绝对值之和
            img_2(x, y) = abs(sum_1) + abs(sum_2);
        end
    end
end

% 处理图像
img_2 = im2uint8(mat2gray(img_2));
end

实验结果如下:
数字图像处理实验之梯度增强边缘