matlab彩色图像锐化

锐化的目的是突出图像的细节。本文使用Laplacian算子进行锐化处理，其他锐化算子的处理类似。

锐化过程如下：

读取RGB彩色图像

分别提取R、G、B通道的分量

设置锐化模板

对图像三个分量分别锐化滤波

将滤波后的三分量组合

% 彩色图像锐化


clc;
clear all;
close all;

rgb1= imread('football.jpg');
rgb=im2double(rgb1);
rgb_R=rgb(:,:,1);
rgb_G= rgb(:,:,2);
rgb_B= rgb(:,:,3);
lapMatrix=[1 1 1;
            1 -8 1;
            1 1 1];
        %模板1 1*8
Matrix2=[1 0 1;
    0 -4 0;
    1 0 1];
%模板21*4

%模板1*8 3个不同的第三参数的锐化滤波
%replicate

f_R=imfilter(rgb_R,lapMatrix,'replicate');
f_G=imfilter(rgb_G,lapMatrix,'replicate');
f_B=imfilter(rgb_B,lapMatrix,'replicate');

%symmetric
f_R2=imfilter(rgb_R,lapMatrix,'symmetric');
f_G2=imfilter(rgb_G,lapMatrix,'symmetric');
f_B2=imfilter(rgb_B,lapMatrix,'symmetric');

%circular
f_R3=imfilter(rgb_R,lapMatrix,'circular');
f_G3=imfilter(rgb_G,lapMatrix,'circular');
f_B3=imfilter(rgb_B,lapMatrix,'circular');

%模板 1*4 3个不同的第三参数的锐化滤波
%replicate
f2_R=imfilter(rgb_R,Matrix2,'replicate');
f2_G=imfilter(rgb_G,Matrix2,'replicate');
f2_B=imfilter(rgb_B,Matrix2,'replicate');

%symmetric
f2_R2=imfilter(rgb_R,Matrix2,'symmetric');
f2_G2=imfilter(rgb_G,Matrix2,'symmetric');
f2_B2=imfilter(rgb_B,Matrix2,'symmetric');

%circular
f2_R3=imfilter(rgb_R,Matrix2,'circular');
f2_G3=imfilter(rgb_G,Matrix2,'circular');
f2_B3=imfilter(rgb_B,Matrix2,'circular');

%模板1 锐化后的3分量
rgb_tmp=cat(3,f_R,f_G,f_B);
rgb_tmp2=cat(3,f_R2,f_G2,f_B2);
rgb_tmp3=cat(3,f_R3,f_G3,f_B3);

%模板2 锐化后的3分量
rgb2_tmp=cat(3,f2_R,f2_G,f2_B);
rgb2_tmp2=cat(3,f2_R2,f2_G2,f2_B2);
rgb2_tmp3=cat(3,f2_R3,f2_G3,f2_B3);

%模板1 3种参数的处理结果与原图相减
rgb_sharped =imsubtract(rgb,rgb_tmp);
rgb_sharped2=imsubtract(rgb,rgb_tmp2);
rgb_sharped3=imsubtract(rgb,rgb_tmp3);

%模板2 3种参数的处理结果与原图相减
rgb2_sharped =imsubtract(rgb,rgb2_tmp);
rgb2_sharped2=imsubtract(rgb,rgb2_tmp2);
rgb2_sharped3=imsubtract(rgb,rgb2_tmp3);

figure;
subplot(341);imshow(rgb);title('原图');
subplot(342);imshow(rgb_R);title('R');
subplot(343);imshow(rgb_G);title('G');
subplot(344);imshow(rgb_B);title('B');

subplot(345);imshow(rgb_tmp);title('cat_r_g_b ');
subplot(346);imshow(rgb_tmp2);title('cat_r_g_b2');
subplot(347);imshow(rgb_tmp3);title('cat_r_g_b3');

subplot(348);imshow(rgb_sharped);title('1*8 replicate');
subplot(349);imshow(rgb_sharped2);title('1*8 symmetric');
subplot(3,4,10);imshow(rgb_sharped3);title('1*8 3 circular');

rgb_diff=imsubtract(rgb_sharped,rgb_sharped2);%replicate sub symmetric

subplot(3,4,11);imshow(rgb_diff);title('replicateSubSymmetric');

figure;
subplot(341);imshow(rgb);title('原图');
subplot(342);imshow(rgb_R);title('R');
subplot(343);imshow(rgb_G);title('G');
subplot(344);imshow(rgb_B);title('B');

subplot(345);imshow(rgb2_tmp);title('cat_r_g_b ');
subplot(346);imshow(rgb2_tmp2);title('cat_r_g_b2');
subplot(347);imshow(rgb2_tmp3);title('cat_r_g_b3');

subplot(348);imshow(rgb2_sharped);title('1*4 replicate');
subplot(349);imshow(rgb2_sharped2);title('1*4 symmetric');
subplot(3,4,10);imshow(rgb2_sharped3);title('1*4 3 circular');

rgb2_diff=imsubtract(rgb2_sharped,rgb2_sharped2);%replicate sub symmetric

subplot(3,4,11);imshow(rgb2_diff);title('replicateSubSymmetric2');

rgb3_diff=imsubtract(rgb2_sharped2,rgb_sharped2);%1*8Sub1*4
figure;
imshow(rgb3_diff);title('1*8Sub1*4');

图1

图2


图3