暗原色原理去雾matlab代码（Matlab）

简要说明:
推荐学习暗原色去雾的原理两篇：
第一篇为pdf格式 900k可以下载：
https://zkxb.jsu.edu.cn/CN/article/downloadArticleFile.do?attachType=PDF&id=2170
第二篇为：
https://blog.csdn.net/yangdashi888/article/details/54571669
本文代码参考上面两篇文章，码出来的。效果不错。

figure,
f=imread("小树林.png");
%f=imresize(f,1/5); 
f=double(f);
[m,n,o]=size(f);
z=zeros(m,n);
for i=1:m
    for j=1:n
        temp=[0,0,0];
        temp(1)=f(i,j,1);
        temp(2)=f(i,j,2);
        temp(3)=f(i,j,3);
        z(i,j)=min(temp);
    end
end
zz=z;
Y=ordfilt2(z,1,ones(3,3),'symmetric');%，相当于3×3的最小值滤波
%subplot(131);
imshow(uint8(f));
title("原始图像");
%subplot(132);
% imshow(uint8(Y));
% title("暗通道图像");

%获取千分之一最大的值 从 z 中
z_max=zeros(m,n);

for t=1:150
    max_value=max(max(z));
    for i=1:m
        for j=1:n
          if max_value == z(i,j)
              z_max(i,j)=1;
              z(i,j)=0;
              max_value=256;
          end
        end
    end
end
num=0;
for i=1:m
    for j=1:n
        if z_max(i,j)==1
            num=num+1;
        end
    end
end

%获取 z_max中点中的最大值
temp=1:1:150;
num=0;
zong_value=0;
for i=1:m
    for j=1:n
        if z_max(i,j)==1
            temp1=[0,0,0];
            temp1(1)=f(i,j,1);
            temp1(2)=f(i,j,2);
            temp1(3)=f(i,j,3);
            num=num+1;
            temp(num)=0;
            temp(num)=max(temp1);
            zong_value=zong_value+temp(num);
        end
    end
end
A=zong_value/150  %此值为A值
%https://blog.csdn.net/yangdashi888/article/details/54571669
%根据 A 求 t
t_x=zeros(m,n);
for i=1:m
    for j=1:n
        temp1=[0,0,0];
        temp1(1)=f(i,j,1)/A;
        temp1(2)=f(i,j,2)/A;
        temp1(3)=f(i,j,3)/A;
        t_x(i,j)=min(temp1);
    end
end
Y=ordfilt2(t_x,1,ones(3,3),'symmetric');%，相当于3×3的最小值滤波
for i=1:m
    for j=1:n
        Y(i,j)=1-0.95*Y(i,j);
    end
end
%在根据 f,A,t求无雾霾图像
t=Y;
J=zeros(m,n,o);
for i=1:m
    for j=1:n
        for k=1:o
            J(i,j,k)=(f(i,j,k)-A)./(max(t(i,j),0.1))+A;
        end
    end
end

figure,
imshow(uint8(J),[]);
title("处理后图像");

结果图：
原始图像：

结果图：