I?=?imread（‘111.bmp‘）;
subplot（121）;
imshow（I）;


I=double（I）;
f=I（::1）;
ff=I（::2）;
fff=I（::3）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k1=4;
k2=5;
r=161;
alf=1458;
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）;????????%高斯函數(shù)1
???end
end

k1=8;??????????????????????????????????????????????????????????????????%高斯函數(shù)2
k2=8;
r=161;
alf=1458;
nn=floor（（r+1）/2）;
for?i=1:r
????for?j=1:r
????????bb（ij）?=exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）;?????
???end
end

k1=0.5;
k2=0.5;
r=161;
alf=1458;
nn=floor（（r+1）/2）;
for?i=1:r
????for?j=1:r
????????bbb（ij）?=exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）;????????%高斯函數(shù)23
end
%%%%%%%%%%%對(duì)R分量的處理%%%%%%%%%%%%%
Img?=?double（f）;
[mn]=size（f）;

aa=125;
for?i=1:m
????for?j=1:n?
????????C（ij）=log（1+aa*（Img（ij）/I（ij）））;
????end
end

K=imfilter（Imgb）;
KK=imfilter（Imgbb）;
KKK=imfilter（Imgbbb）;

for?i=1:m
????for?j=1:n??????
???????G（ij）=1/3*（log（Img（ij）+1）-log（K（ij）+1））;
????????G（ij）=1/3*（log（Img（ij）+1）-log（KK（ij）+1））+G（ij）;
?????????G（ij）=C（ij）*（1/3*（log（Img（ij）+1）-log（KKK（ij）+1））+G（ij））;
????end
end

mi=min（min（G））;
ma=max（max（G））;
???????L=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%%%對(duì)G分量的處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Img?=?double（ff）;
[mn]=size（ff）;

aa=125;
for?i=1:m
????for?j=1:n?
????????CC（ij）=log（1+aa*（Img（ij）/I（ij）））;
????end
end

K=imfilter（Imgb）;
KK=imfilter（Imgbb）;
KKK=imfilter（Imgbbb）;
for?i=1:m
????for?j=1:n??????
???????G（ij）=1/3*（log（Img（ij）+1）-log（K（ij）+1））;
????????G（ij）=1/3*（log（Img（ij）+1）-log（KK（ij）+1））+G（ij）;
?????????G（ij）=CC（ij）*（1/3*（log（Img（ij）+1）-log（KKK（ij）+1））+G（ij））;
????end
end

mi=min（min（G））;
ma=max（max（G））;
???????LL=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%%隨B分量的處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Img?=?double（fff）;
[mn]=size（fff）;

aa=125;
for?i=1:m
????for?j=1:n?
????????CCC（ij）=log（1+aa*（Img（ij）/I（ij）））;
????end
end

K=imfilter（Imgb）;
KK=imfilter（Imgbb）;
KKK=imfilter（Imgbbb）;

for?i=1:m
????for?j=1:n??????
???????G（ij）=1/3*（log（Img（ij）+1）-log（K（ij）+1））;
????????G（ij）=1/3*（log（Img（ij）+1）-log（KK（ij）+1））+G（ij）;
?????????G（ij）=CCC（ij）*（1/3*（log（Img（ij）+1）-log（KKK（ij）+1））+G（ij））;
????end
end

mi=min（min（G））;
ma=max（max（G））;

???????LLL=（G-mi）*255/（ma-mi）;
%%%%%%%%%%%%對(duì)彩色圖像的綜合處理?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
msrcr=cat（3LLLLLL）;
subplot（122）;
imshow（uint8（msrcr））;
imwrite（uint8（msrcr）‘msrcrsretnx.bmp‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2612??2009-09-06?20:55??retinex\MSRCR.m
?????文件????????3183??2009-09-07?10:07??retinex\MSR_new.m
?????文件??????????24??2018-04-21?10:52??retinex\Readme.txt
?????文件????????2693??2010-04-23?21:00??retinex\retinex.m
?????文件????????2775??2009-03-08?18:33??retinex\retinex_frankle_mccann.m
?????文件????????5392??2009-03-08?18:33??retinex\retinex_mccann99.m
?????目錄???????????0??2018-04-21?10:51??retinex\