何凯明去雾MATLAB代码

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