%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%++++++++++++++++++++++自71王潘改進Reinhard算法+++++++++++++++++++++++++++++
%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%=========================打開參考文件======================================
Ref=double（imread（‘pic1-4.jpg‘））;
RefR=Ref（::1）;
RefG=Ref（::2）;
RefB=Ref（::3）;
%=========================打開局部參考文件==================================
Refj=double（imread（‘pic1-41.jpg‘））;
RefRj=Ref（::1）;
RefGj=Ref（::2）;
RefBj=Ref（::3）;
%=========================打開源文件========================================
Sou=double（imread（‘pic1-5.jpg‘））;
SouR=Sou（::1）;
SouG=Sou（::2）;
SouB=Sou（::3）;
%=========================打開局部文件======================================
Souj=double（imread（‘pic1-51.jpg‘））;
SouRj=Sou（::1）;
SouGj=Sou（::2）;
SouBj=Sou（::3）;
%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%++++++++++++++++++++++++++轉換到lab空間+++++++++++++++++++++++++++++++++++
%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%===============================定義轉換函數===============================
con1=[0.38110.57830.0402;0.19670.72440.0782;0.02410.12880.8444];
con2=[1?1?1;1?1?-2;1?-1?0];
con3=[1/sqrt（3）?0?0;0?1/sqrt（6）?0;0?0?1/sqrt（2）];
%============================參考圖像轉到lab空間============================
[widthheightdepth]=size（Ref）;
l1=zeros（widthheight）;
alpha1=zeros（widthheight）;
beta1=zeros（widthheight）;
%================中轉到LMS==============%
L1=zeros（widthheight）;
M1=zeros（widthheight）;
S1=zeros（widthheight）;
for?i=1:width
????for?j=1:height
????????m=con1*[RefR（ij）;RefG（ij）;RefB（ij）];
????????L1（ij）=m（1）;
????????M1（ij）=m（2）;
????????S1（ij）=m（3）;
????end
end
%================取LMS對數============%
L1=log10（L1）;
M1=log10（M1）;
S1=log10（S1）;
%==============轉至lab空間============%
for?i=1:width
????for?j=1:height
????????m=con3*con2*[L1（ij）;M1（ij）;S1（ij）];
????????l1（ij）=m（1）;
????????alpha1（ij）=m（2）;
????????beta1（ij）=m（3）;
????end
end
%============================局部參考圖像轉到lab空間============================
[widthjheightjdepthj]=size（Refj）;
l1j=zeros（widthjheightj）;
alpha1j=zeros（widthjheightj）;
beta1j=zeros（widthjheightj）;
%================中轉到LMS==============%
L1j=zeros（widthjheightj）;
M1j=zeros（widthjheightj）;
S1j=zeros（widthjheightj）;
for?i=1:widthj
????for?j=1:heightj
????????mj=con1*[RefRj（ij）;RefGj（ij）;RefBj（ij）];
????????L1j（ij）=mj（1）;
????????M1j（ij）=mj（2）;
????????S1j（ij）=mj（3）;
????end
end
%================取LMS對數============%
L1j=log10（L1j）;
M1j=log10（M1j）;
S1j=log10（S1j）;
%==============轉至lab空間============%
for?i=1:widthj
????for?j=1:heightj
????????m=con3*con2*[L1j（ij）;M1j（ij）;S1j（ij）];
????????l1j（ij）=mj（1）;
????????alpha1j（ij）=mj（2）;
????????beta1j（ij）=mj（3）;
????end
end
%=============================源圖像轉到lab空間============================
[widheidep]=size（Sou）;
l2=zer

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

?????文件????????393??2011-01-10?21:08??2007011565_王潘_大作業\Readme.txt

?????文件????1000783??2010-12-28?14:38??2007011565_王潘_大作業\圖像處理大作業實驗報告.pdf

?????文件??????21294??2010-12-28?18:57??2007011565_王潘_大作業\處理結果\c1=0.2c2=0.8.JPG

?????文件??????21568??2010-12-28?18:59??2007011565_王潘_大作業\處理結果\c1=0.5c2=0.5.JPG

?????文件??????21465??2010-12-28?18:59??2007011565_王潘_大作業\處理結果\c1=0.6c2=0.4.JPG

?????文件??????21441??2010-12-28?19:00??2007011565_王潘_大作業\處理結果\c1=0.8c2=0.2.JPG

?????文件??????21472??2010-12-28?18:55??2007011565_王潘_大作業\處理結果\Reinhard.jpg

?????文件??????69271??2010-12-24?21:37??2007011565_王潘_大作業\程序\pic1-4.JPG

?????文件???????9261??2010-12-26?13:39??2007011565_王潘_大作業\程序\pic1-41.JPG

?????文件??????52580??2010-12-24?21:40??2007011565_王潘_大作業\程序\pic1-5.JPG

?????文件???????8533??2010-12-26?14:20??2007011565_王潘_大作業\程序\pic1-51.JPG

?????文件???????9468??2011-01-10?18:41??2007011565_王潘_大作業\程序\Reinhard?improved.m

?????文件???????5876??2011-01-10?18:36??2007011565_王潘_大作業\程序\Reinhard.m

?????目錄??????????0??2011-03-22?19:46??2007011565_王潘_大作業\處理結果

?????目錄??????????0??2011-03-22?19:46??2007011565_王潘_大作業\程序

?????目錄??????????0??2011-03-22?19:46??2007011565_王潘_大作業

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

??????????????1263405????????????????????16