求基于雾气遮罩去雾算法在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('原始图像');
[h,w,s]=size(I);
min_I=zeros(h,w);

%下面取得暗影通道图像
for i=1:h
for j=1:w
dark_I(i,j)=min(I(i,j,:));
end
end

subplot(223)
imshow(dark_I);
title('dark channnel的图形');

Max_dark_channel=double(max(max(dark_I))) %天空亮度
dark_channel=double(dark_I);
t=1-w0*(dark_channel/Max_dark_channel); %取得透谢分布率图

subplot(224)
T=uint8(t*255);