clc
clear
f=imread（‘2.jpg‘）;
fr=f（:?:?1）;?fg=f（:?:?2）;?fb=f（:?:?3）;%RGB通道
mr=mat2gray（im2double（fr））;?mg=mat2gray（im2double（fg））;?mb=mat2gray（im2double（fb））;%數據類型歸一化
alf1=1458;?%定義標準差alf=a^2/2??a=54
n=161;%定義模板大小?
n1=floor（（n+1）/2）;%計算中心?
for?i=1:n?
????for?j=1:n?
??????b（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（4*alf1））/（pi*alf1）;?%高斯函數
????end?
end?
nr1?=?imfilter（mrb‘conv‘?‘replicate‘）;ng1?=?imfilter（mgb‘conv‘?‘replicate‘）;nb1?=?imfilter（mbb‘conv‘?‘replicate‘）;%卷積濾波
ur1=log（nr1）;?ug1=log（ng1）;?ub1=log（nb1）;
tr1=log（mr）;tg1=log（mg）;tb1=log（mb）;
yr1=（tr1-ur1）/3;yg1=（tg1-ug1）/3;yb1=（tb1-ub1）/3;
alf2=53.38;?%定義標準差alf=a^2/2????a=10.3325
x=31;%定義模板大小?
x1=floor（（n+1）/2）;%計算中心?
for?i=1:n?
????for?j=1:n?
??????a（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（4*alf2））/（6*pi*alf2）;?%高斯函數
????end?
end?
nr2?=?imfilter（mra‘conv‘?‘replicate‘）;ng2?=?imfilter（mga‘conv‘?‘replicate‘）;nb2?=?imfilter（mba‘conv‘?‘replicate‘）;%卷積濾波
ur2=log（nr2）;?ug2=log（ng2）;?ub2=log（nb2）;
tr2=log（mr）;tg2=log（mg）;tb2=log（mb）;
yr2=（tr2-ur2）/3;yg2=（tg2-ug2）/3;yb2=（tb2-ub2）/3;
alf3=13944.5;?%定義標準差alf=a^2/2??a=167
l=501;%定義模板大小?
l1=floor（（n+1）/2）;%計算中心?
for?i=1:n?
????for?j=1:n?
??????e（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（4*alf3））/（4*pi*alf3）;?%高斯函數
????end?
end?
nr3?=?imfilter（mre‘conv‘?‘replicate‘）;ng3?=?imfilter（mge‘conv‘?‘replicate‘）;nb3?=?imfilter（mbe‘conv‘?‘replicate‘）;%卷積濾波
ur3=log（nr3）;?ug3=log（ng3）;?ub3=log（nb3）;
tr3=log（mr）;tg3=log（mg）;tb3=log（mb）;
yr3=（tr3-ur3）/3;yg3=（tg3-ug3）/3;yb3=（tb3-ub3）/3;
dr=yr1+yr2+yr3;dg=yg1+yg2+yg3;db=yb1+yb2+yb3;
????
cr=im2uint8（dr）;?cg=im2uint8（dg）;?cb=im2uint8（db）;
z=cat（3?cr?cg?cb）;??????????
figure?imshow（z）
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-11-15?09:14??retinex?matlab代碼\
?????文件??????267239??2017-10-23?15:15??retinex?matlab代碼\1.jpg
?????文件???????12564??2017-10-23?15:15??retinex?matlab代碼\10.jpg
?????文件??????305425??2017-10-23?15:15??retinex?matlab代碼\11.jpg
?????文件???????25001??2017-10-23?15:15??retinex?matlab代碼\12.jpg
?????文件???????14086??2017-10-23?15:15??retinex?matlab代碼\13.jpg
?????文件??????178720??2017-10-23?15:15??retinex?matlab代碼\14.jpg
?????文件??????326323??2017-10-23?15:15??retinex?matlab代碼\2.jpg
?????文件??????325289??2017-10-23?15:15??retinex?matlab代碼\3.jpg
?????文件??????544839??2017-10-23?15:15??retinex?matlab代碼\4.JPG
?????文件??????202104??2017-10-23?15:15??retinex?matlab代碼\5.jpg
?????文件??????312216??2017-10-23?15:15??retinex?matlab代碼\6.jpg
?????文件???????98152??2017-10-23?15:15??retinex?matlab代碼\7.jpg
?????文件???????43970??2017-10-23?15:15??retinex?matlab代碼\8.jpg
?????文件????????6635??2017-10-23?15:15??retinex?matlab代碼\9.jpg
?????文件?????????664??2017-10-23?15:15??retinex?matlab代碼\Untitled.m
?????文件????????1415??2017-10-23?15:15??retinex?matlab代碼\Untitled2.m
?????文件????????1774??2017-11-15?09:14??retinex?matlab代碼\retinex.m
?????文件????????2009??2017-10-23?15:15??retinex?matlab代碼\retinex_matlab.txt
?????文件????????2736??2017-10-23?15:15??retinex?matlab代碼\tssrcolor20130517.m