数字图像处理——Sobel算子锐化、Prewitt算子锐化

一、Sobel算子锐化

%函数名称为Image_Sobel,输入参数Image,输出参数IMAGEfunction [IMAGE] = Image_Sobel(Image)%获取矩阵的行、列、波段数[m,n,bands] = size(Image);%定义模板大小，假设模板大小3×3A = 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:bandsfor i = 1+A:m-Afor j = 1+A:n-Atemp = Image(i-A:i+A,j-A:j+A,k);Image1(i,j,k) = abs(sum(sum(Sobelx.*temp)))+abs(sum(sum(Sobely.*temp)));endendendIMAGE = 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,输出参数IMAGEfunction [IMAGE] = Image_Prewitt(Image)%获取矩阵的行、列、波段数[m,n,bands] = size(Image);%定义模板大小，假设模板大小3×3A = 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:bandsfor i = 1+A:m-Afor j = 1+A:n-Atemp = Image(i-A:i+A,j-A:j+A,k);Image1(i,j,k) = abs(sum(sum(Prewittx.*temp)))+abs(sum(sum(Prewitty.*temp)));endendendIMAGE = 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);