function?darktest（filename）
%暗影去霧算法
%filename------文件名或文件絕對路徑
%用法：darktest（‘7.png‘）
close?all
clc

w0=0.65;???%0.65??乘積因子用來保留一些霧，1時完全去霧????
t0=0.1;

I=imread（filename）;
figure;
set（gcf‘outerposition‘get（0‘screensize‘））;
subplot（221）
imshow（I）;
title（‘原始圖像‘）;
[hws]=size（I）;
min_I=zeros（hw）;???????????

%下面取得暗影通道圖像
for?i=1:h?????????????????
????for?j=1:w
????????dark_I（ij）=min（I（ij:））;
????end
end

subplot（223）
imshow（dark_I）;
title（‘dark?channnel的圖形‘）;

Max_dark_channel=doubl