MSR MATLAB

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=