1. 参考文献
2.Fattal模型实现
2.1 Fattal模型Demo
clc
clear all;
input = imread('6.png');
% input = imread('1.bmp');
% input = imread('plane.jpg');
[h,w,s]=size(input);
figure, imshow(input), title('Input image') ;
imfog=double(input);
Im=double(input)/255;
I=zeros(h*w,3);
I(:,1)=reshape(Im(:,:,1),h*w,1);
I(:,2)=reshape(Im(:,:,2),h*w,1);
I(:,3)=reshape(Im(:,:,3),h*w,1);
B = [0.8 0.8 0.9]';
[est_t est_l est_eta] = estimate(B, I);
A=[175 172 172];
t = reshape(est_t,h,w);
min1=min(min(t));
max1=max(max(t));
t=(t-min1)/(max1-min1); %线性拉伸
t=max(t,0.1);
for j=1:3
im_defog(:,:,j) = uint8(1.0*((imfog(:,:,j) - (1-t)*double(A(j)))./t));
end
figure, imshow(im_defog), title('Output image');