图像除雾——Fattal模型

    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