function?y=MSR（fabc）
f=double（f）;
R=f（::1）;
G=f（::2）;
B=f（::3）;
I=（0.299.*R+0.587.*G+0.114.*B）./255;
[mn]=size（I）;

%?for?i=1:m
%?????for?j=1:n
%?????????N1（ij）=exp（-（i^2+j^2）/a^2）;
%?????????N2（ij）=exp（-（i^2+j^2）/b^2）;
%?????????N3（ij）=exp（-（i^2+j^2）/c^2）;
%?????end
%?end
%?K1=1/sum（sum（N11）2）;
%?K2=1/sum（sum（N21）2）;
%?K3=1/sum（sum（N31）2）;
%?
%?N1=K1*N1;
%?N2=K2*N2;
%?N3=K3*N3;

%?F1=imfilter（IN1‘replicate‘）;
%?F2=imfilter（IN2‘replicate‘）;
%?F3=imfilter（IN3‘replicate‘）;
%%%%%%%%%%設定高斯參數%%%%%%?
%?sigma_1=15;???%三個高斯核?
%?sigma_2=80;?
%?sigma_3=250;?
[x?y]=meshgrid（（-（size（I1）-1）/2）:（size（I1）/2）（-（size（I2）-1）/2）:（size（I2）/2））;???
gauss_1=exp（-（x.^2+y.^2）/（2*a*a））;??%計算高斯函數?
Gauss_1=gauss_1/sum（gauss_1（:））;??%歸一化處理?
gauss_2=exp（-（x.^2+y.^2）/（2*b*b））;?
Gauss_2=