clear?all
clc
f=imread（‘jiedao.jpg‘）;
figure（1）;
subplot（121）;
imshow（f）;title（‘原圖像‘）;
fr=f（:?:?1）;?fg=f（:?:?2）;?fb=f（:?:?3）;%RGB通道
mr=mat2gray（im2double（fr））;?mg=mat2gray（im2double（fg））;?mb=mat2gray（im2double（fb））;%數(shù)據(jù)類型歸一化
[pq]=size（mr）;[st]=size（mg）;[mn]=size（mb）;
mr=ones（pq）.*（1/255）+mr;mg=ones（st）.*（1/255）+mg;mb=ones（mn）.*（1/255）+mb;
alf1=1458;?%定義標(biāo)準(zhǔn)差alf=a^2/2???a=54
n=161;%定義模板大小?
n1=floor（（n+1）/2）;%計(jì)算中心?
for?i=1:n?
for?j=1:n?
???b（ij）?=exp（-（（i-n1）^2+（j-n1）^2）/（4*alf1））/（pi*alf1）;?%高斯函數(shù)
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;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1947??2011-05-03?09:04??Retinex.m

-----------?---------??----------?-----??----

?????????????????1947????????????????????1