clc;
clear?all;
img_name=‘jiedao.jpg‘;
%?原始圖像
I=double（imread（img_name））/255;
%?獲取圖像大小
[hwc]=size（I）;
win_size?=?7;
img_size=w*h;
figure?imshow（I）;
win_dark=ones（hw）;
%計算分塊darkchannel
?for?j=1+win_size:w-win_size
????for?i=win_size+1:h-win_size
????????m_pos_min?=?min（I（ij:））;
????????for?n=j-win_size:j+win_size
????????????for?m=i-win_size:i+win_size
????????????????if（win_dark（mn）>m_pos_min）
????????????????????win_dark（mn）=m_pos_min;
????????????????end
????????????end
????????end
????end
?end
?%選定精確dark?value坐標
%?win_b?=?zeros（img_size1）;
?figure?imshow（win_dark）;
?win_t=1-0.95*win_dark;
?win_b=zeros（img_size1）;
for?ci=1:h
????for?cj=1:w
????????if（rem（ci-815）<1）
????????????if（rem（cj-815）<1）
????????????????win_b（ci*w+cj）=win_t（ci*w+cj）;
????????????end
????????end
????end
end
?
%顯示分塊darkchannel
%figure?imshow（win_dark）;
neb_size?=?9;
win_size?=?1;
epsilon?=?0.000001;
%指定矩陣形狀
indsM=reshape（1:img_sizehw）;
%計算矩陣L
tlen?=?

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

?????文件??????14451??2012-10-16?21:42??wu1.JPG

?????文件?????192981??2012-10-24?09:56??9new.PNG

?????文件???????3954??2012-10-24?11:26??dark_channelnew.m

?????文件??????15032??2011-04-21?20:09??jiedao.jpg

?????文件??????93838??2012-10-11?21:43??7-1.PNG

?????文件?????100093??2012-10-11?21:49??81.PNG

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

???????????????420349????????????????????6