隻提供代碼,不提供思路,思路看課本。
主函數:
% 讀入圖像
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
實驗結果如下:
![](https://img.laitimes.com/img/__Qf2AjLwojIjJCLyojI0JCLiAzNfRHLGZkRGZkRfJ3bs92YsYTMfVmepNHL9M3RihGZzgFbOdVW350MMBjVtJWd0ckW65UbM5WOHJWa5kHT20ESjBjUIF2X0hXZ0xCMx81dvRWYoNHLrdEZwZ1Rh5WNXp1bwNjW1ZUba9VZwlHdssmch1mclRXY39CXldWYtlWPzNXZj9mcw1ycz9WL49zZwpmL2YTN4UDO0UTMzIzNwkTMwIzLc52YucWbp5GZzNmLn9Gbi1yZtl2Lc9CX6MHc0RHaiojIsJye.jpg)