數字圖像處理——Sobel算子銳化、Prewitt算子銳化
一、Sobel算子銳化
%函數名稱爲Image_Sobel,輸入參數Image,輸出參數IMAGE
function [IMAGE] = Image_Sobel(Image)
%獲取矩陣的行、列、波段數
[m,n,bands] = size(Image);
%定義模板大小,假設模板大小3×3
A = 1;
%定義Sobel算子x,y方向矩陣
Sobelx = [-1 -2 -1;0 0 0;1 2 1];
Sobely = [-1 0 1;-2 0 2;-1 0 1];
%初始化矩陣
Image1 = zeros(m,n,bands);
IMAGE = Image;
%Sobel算子
for k = 1:bands
for i = 1+A:m-A
for j = 1+A:n-A
temp = Image(i-A:i+A,j-A:j+A,k);
Image1(i,j,k) = abs(sum(sum(Sobelx.*temp)))+abs(sum(sum(Sobely.*temp)));
end
end
end
IMAGE = Image + Image1;
%畫圖,左右分別表示原圖和兩幅處理後的圖像
figure(1);
subplot(1,3,1);
imshow(uint8(Image));
title('原始圖像');
subplot(1,3,2);
imshow(uint8(Image1));
title('邊緣提取圖像');
subplot(1,3,3);
imshow(uint8(IMAGE));
title('Sobel算子銳化後的圖像')
end
二、prewitt算子銳化
%函數名稱爲Image_Prewitt,輸入參數Image,輸出參數IMAGE
function [IMAGE] = Image_Prewitt(Image)
%獲取矩陣的行、列、波段數
[m,n,bands] = size(Image);
%定義模板大小,假設模板大小3×3
A = 1;
%定義Prewitt算子x,y方向矩陣
Prewittx = [-1 -2 -1;0 0 0;1 2 1];
Prewitty = [-1 0 1;-2 0 2;-1 0 1];
%初始化矩陣
Image1 = zeros(m,n,bands);
IMAGE = Image;
%Sobel算子
for k = 1:bands
for i = 1+A:m-A
for j = 1+A:n-A
temp = Image(i-A:i+A,j-A:j+A,k);
Image1(i,j,k) = abs(sum(sum(Prewittx.*temp)))+abs(sum(sum(Prewitty.*temp)));
end
end
end
IMAGE = Image + Image1;
%畫圖,左中右分別表示原圖和兩幅處理後的圖像
figure(1);
subplot(1,3,1);
imshow(uint8(Image));
title('原始圖像');
subplot(1,3,2);
imshow(uint8(Image1));
title('邊緣提取圖像');
subplot(1,3,3);
imshow(uint8(IMAGE));
title('Prewitt銳化後的圖像')
end
調用方法
Image = double(imread('當前目錄下的圖片名稱.jpg'));
Sobel = Image_Sobel(Image);
Prewitt = Image_Prewitt(Image);