ｒｅｔｉｎｅｘ的ｍａｔｌａｂ程序

clear;
%打开任意一幅图片
[filepath]=uigetfile（{‘*.jpg‘‘（*.jpg） jpg图片‘;‘*.bmp‘‘（*.bmp） BMP图片‘;‘*.jpeg‘‘（*.jpeg） jpeg图片‘;‘*.*‘‘（*.*） 所有文件‘}[‘请选择图片！（^_^）‘]）;   
I=imread（[pathfile]）;
%subplot（121）;
figureimshow（I）;             %figure4待处理图像
%[lcd]=size（I）;
x=I（::1）;
xx=I（::2）;
xxx=I（::3）;
I1=rgb2hsi（I）;   %将RGB中的原图转换到HSV空间
f=I1（::1）;
ff=I1（::2）;
fff=I1（::3）;
%%%%%%%%%%%%%构造高斯核%%%%%%%%%%%%%%%%%%%%%%%
k1=8;
k2=8;
%r=161;
r=100
alf=1600;
nn=floor（（r+1）/2）;
for i=1:r
    for j=1:r
        b（ij） =exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）; 
         %b（ij） =exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））; %高斯函数
   end
end
%%%%%%%%%%%对R分量的处理%%%%%%%%%%%%%
%Img = double（f）;
%K=imfilter（Imgb）;  
%G=log（Img+1）-log（K+1）;
%G=Img.*K;
%mi=min（min（G））;
%ma=max（max（G））;
%       L=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%%%对G分量的处理%%%%%%%%%%%%%
%Img = double（ff）;
%K=imfilter（Imgb）;  
%G=log（Img+1）-log（K+1）;
%G=Img.*K;
%mi=min（min（G））;
%ma=max（max（G））;
       %LL=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%%随B分量的处理%%%%%%%%%