%?S=imread（strcat（‘2.bmp‘））;
S=imread（‘2.jpg‘）;
%?S=im2gray（filename）;
S=rgb2gray（S）
[mn]=size（S）;
S=double（S）;

size=m*n;
mean=sum（sum（S）2）/size;
R=zeros（mn）;
R=repmat（meanmn）;

%做高斯函數
c=80;
for?i=1:m
????for?j=1:n
????????N（ij）=exp（-（i^2+j^2）/c^2）;
????end
end
K=1/sum（N（:））;
N1=K*N;
G=conv2（SN1‘same‘）;??%尺度一

c=120;
for?i=1:m
????for?j=1:n
????????N（ij）=exp（-（i^2+j^2）/c^2）;
????end
end
K=1/sum（N（:））;
N1=K*N;
G1=conv2（SN1‘same‘）;??%尺度二

c=250;
for?i=1:m
????for?j=1:n
????????N（ij）=exp（-（i^2+j^2）/c^2）;
????end
end
K=1/sum（N（:））;
N1=K*N;
G2=conv2（SN1‘same‘）;??%尺度三

%計算MSR

SS=zeros（mn）;
for?i=1:m
???for?j=1:n
????SS（ij）=log（S（ij）+0.0005）;
????GG（ij）=log（G（ij）+0.0005）;
????GG1（ij）=log（G1（ij）+0.0005）;
????GG2（ij）=log（G2（ij）+0.0005）;
????V（ij）=?SS（ij）-GG（ij）;
????V1（ij）=?SS（ij）-GG1（ij）;
????V2（ij）=?SS（ij）-GG2（ij）;
???end
end

Min=min（min（V））;
Max=max（max（V））;
Min=min（min（V1））;
Max=max（max（V1））;
Min=min（min（V2））;
Max=max（max（V2））;
for?i=1:m
???for?j=1:n
??????VV（ij）=（V（ij）-Min）*255/（Max-Min）;
??????VV1（ij）=（V1（ij）-Min）*255/（Max-Min）;
??????VV2（ij）=（V2（ij）-Min）*255/（Max-Min）;
???end
end
SumVV=VV/3+VV1/3+VV2/3;
VV=uint8（VV）;
VV1=uint8（VV1）;
VV2=uint8（VV2）;
SumVV=uint8（SumVV）;
S=uint8（S）;
figuresubplot（151）imshow（S）;
subplot（152）imshow（VV）;
subplot（153）imshow（VV1）;
subplot（154）imshow（VV2）;
subplot（155）imshow（SumVV）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????740??2010-04-15?20:14??Retinex\SSR.m

?????文件???????1306??2010-04-15?20:17??Retinex\MySSR.m

?????文件???????1492??2010-04-15?20:19??Retinex\MSR.m

?????目錄??????????0??2010-04-20?09:23??Retinex

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

?????????????????3538????????????????????4