%%%????暗通道去霧函數(shù)，可自動(dòng)檢測(cè)天空，采用不同的算法。
%%%????src--------------輸入彩色圖像
%%%????rad1rad2rad3---最小值濾波、引導(dǎo)濾波、天空區(qū)域形態(tài)學(xué)所用算子尺寸
%%%????Amax-------------大氣光最大值
%%%????thre-------------天空區(qū)域梯度閾值
%%%????tao--------------天空區(qū)域亮度閾值
%%%????ratio------------天空區(qū)域比例閾值
%%%????enhance----------是否增強(qiáng)
function?imgout=DarkchDehazed（srcrad1rad2rad3Amaxthretaoratioenhance）
[HWC]=size（src）;
if?C~=3
????error（‘Not?3?channels!‘）;
end
gray=double（rgb2gray（src））;
%%??是否含天空
[FXFY]?=?gradient（gray）;
Indense=sqrt（FX.^2+FY.^2）;
candid=gray>tao?&?Indense<=thre;???????????%???候選區(qū)域：亮度超過閾值，梯度幅度為0.
%?imshow（candid）
SE1=strel（‘square‘rad3）;
issky=imclose（candidSE1）;????????????%???天空標(biāo)為1，非天空0.
%?figure;imshow（issky）
disp（sum（sum（issky））/（H*W））;
if?sum（sum（issky））/（H*W）>=ratio
????disp（‘sky?exist.‘）;
????flag=true;
else
????disp（‘none‘）;
????flag=false;
end
chr=double（src（::1））;
chg=double（src（::2））;
chb=double（src（::3））;
perc=0.001;
numPt=floor（perc*H*W）;
A=zeros（31）;
omg=0.95;
SE2=strel（‘square‘rad1）;???????????%???最小值濾波形態(tài)學(xué)
%???含天空，以天空部分估計(jì)A，分別計(jì)算透射率。
%???預(yù)估透射率
if?flag
????red_s=double（chr）.*issky;
????green_s=double（chg）.*issky;
????blue_s=double（chb）.*issky;
????sky=cat（3red_sgreen_sblue_s）;
????skygray=rgb2gray（uint8（sky））;
????skygray=reshape（skygray1[]）;
????red_s=reshape（red_s1[]）;
????green_s=reshape（green_s1[]）;
????blue_s=reshape（blue_s1[]）;
????[~ix]=sort（skygray‘descend‘）;
????sortedred_s=red_s（ix）;
????sortedgreen_s=green_s（ix）;
????sortedblue_s=blue_s（ix）;
????A（1）=min（mean（sortedred_s（1:numPt））Amax）;
????A（2）=min（mean（sortedgreen_s（1:numPt））Amax）;
????A（3）=min（mean（sortedblue_s（1:numPt））Amax）;
????nonsky=~issky;
????red_ns=double（chr）.*nonsky/A（1）;
????green_ns=double（chg）.*nonsky/A（2）;
????blue_ns=double（chb）.*nonsky/A（3）;
????red_ns=reshape（red_ns1[]）;
????green_ns=reshape（green_ns1[]）;
????blue_ns=reshape（blue_