clc;
clear ;
im=imread('1.jpg');
% im(x1:xn,y1:ym,1);
R=im(:,:,1);
G=im(:,:,2);
B=im(:,:,3);
% 高斯滤波
HSIZE= min(size(im,1),size(im,2));%高斯卷积核尺寸
q=sqrt(2);
SIGMA1=15;%论文里面的c
SIGMA2=80;
SIGMA3=250;
F1 = fspecial('gaussian',18,100);
gaus1= imfilter(R, F1, 'replicate');
% gaus2= imfilter(R, F2, 'replicate');
% gaus3= imfilter(R, F3, 'replicate');
% gausa=(gaus1+gaus2+gaus3)/3;
gaus4= imfilter(G, F1, 'replicate');
% gaus5= imfilter(G, F2, 'replicate');
% gaus6= imfilter(G, F3, 'replicate');
% gausb=(gaus4+gaus5+gaus6)/3;
gaus7= imfilter(B, F1, 'replicate');
% gaus8= imfilter(B, F2, 'replicate');
% gaus9= imfilter(B, F3, 'replicate');
% gausc=(gaus7+gaus8+gaus9)/3;
c=cat(3,gaus1,gaus4,gaus7);
figure;
subplot(231),imshow(im),title('原图');
subplot(232),imshow(gaus1),title('R');
subplot(233),imshow(gaus4),title('G');
subplot(234),imshow(gaus7),title('B');
subplot(235),imshow(c),title('综合');
评论1