基于可见边的去雾图像质量评价

该代码是基于可见边的图像质量评价程序，matlab实现，里面附着相应的文献及测试图像。

% 23/02/2011
% IFSTTAR copyright
%
% The approach is described in details in 
%
% “Blind Contrast Restoration Assessment by Gradient Ratioing at Visible Edges“
% by N. Hautiere J.-P. Tarel D. Aubert and E. Dumont
% in proceedings of International Congress for Stereology （ICS‘07） 
% Saint Etienne France August 30-September 7 2007.
% http://perso.lcpc.fr/tarel.jean-philippe/publis/ics07.html
%

%%%% Cleaning
clc
clear all
close all

%%%% Images reading: the input 2 images must be grayscale

NameOri=‘Original.pgm‘;
%NameResto=‘Restored.pgm‘;
NameResto=‘Restored2.pgm‘;
%NameResto=‘Restored3.pgm‘;
%NameResto=‘Restored4.pgm‘;
%NameResto=‘Restored5.pgm‘;

I1=Imread（NameOri）;
I1=double（I1）;
% if the input image is a color image use following line
% I1=double（rgb2gray（uint8（I1）））;

[nlnc]=size（I1）;

R1=Imread（NameResto）;
R1=double（R1）;
% if the input image is a color image use following line
% R1=double（rgb2gray（uint8（R1）））;

%%%% Figure 1
figure（1）
colormap gray
subplot（121）
imagesc（I1）
axis image
title（‘Original‘）
subplot（122）
imagesc（R1）
axis image
title（‘Restored‘）

%%%% Sobel Gradient
Sy = double（[1 2 1;0 0 0;-1 -2 -1]）;
GIy = imfilter（I1Sy‘symmetric‘）;
GRy1 = imfilter（R1Sy‘symmetric‘）;

Sx = double（[1 0 -1; 2 0 -2; 1 0 -1]）;
GIx = imfilter（I1Sx‘symmetric‘）;
GRx1 = imfilter（R1Sx‘symmetric‘）;

GI=sqrt（（GIx.^2）+（GIy.^2））;
GR1=sqrt（（GRx1.^2）+（GRy1.^2））;

minGI=min（GI（:））;
maxGI=max（GI（:））;

%%%% Figure 2
figure（2）
colormap gray
subplot（121）
imagesc（GI[minGI maxGI]）;
title（[‘Gradients of the original image‘]）; 
axis image
colorbar
subplot（122）
imagesc（GR1[minGI maxGI]）;
title（[‘Gradients of the restored image‘]）; 
axis image
colorbar

%%%% Contrast Computation at 5%
tic
[C1 Crr1]=functionContrastAt5PerCent（R1）;
[Ci Crri]=functionContrastAt5PerCent（I1）;
toc

minCrri=min（Crri（:））;
maxCrri=max（Crri（:））;

%%%% Figure 3
figure（3）
colormap gray
subplot（121）
imagesc（Crri[minCrri  maxCrri]）;
axis image
title（[‘Visible edge in the original image‘]）;
colorbar
subplot（122）
imagesc（Crr1[minCrri  maxCrri]）;
axis image
title（[‘Visible edge in the restored image‘]）;
colorbar

%%%% Visible Gradients Ratio
Cratio1=zeros（nlnc）;
Cratio1（Crr1>0）=GR1（Crr1>0）./GI（Crr1>0）;

rmin=1;
rmax=10;

%%%% Figure 4
figure（4）
imagesc（Cratio1[rmin rmax]）;
axis image
title（[‘Visible gradients ratio between ‘num2str（rmin） ‘ and ‘num2str（rmax）]）; 
colormap jet
colorbar

%%%% Descriptor computation

% Rate of new visible edges
whitePixels1=sum（C1（:））;
whitePixelsi=sum（Ci（:））;
e1=（whitePixels1-whitePixelsi）/whitePixelsi;

% Number of saturated pixels after restoration
ns1=sum（R1（:）==255 |R1（:）==0）;
ns1=ns1/（nl*nc）;

% Restoration quality （geometric mean ratios of visibility level）
XX=log（Cratio1）;
r1=exp（（1/（whitePixels1））*nansum（XX（isfinite（XX））））;

%%%% Figure 5: Final result with the visible edges at 5% and descriptors
f

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     733002  2014-04-02 11:33  EvaluationDescriptor\EvaluationDescriptor\Blind Contrast Restoration Assessment by Gradient Ratioing at Visible Edges.pdf

     文件     406643  2014-04-03 16:33  EvaluationDescriptor\EvaluationDescriptor\CALIS20140331530100.pdf

     文件       3441  2014-04-09 16:22  EvaluationDescriptor\EvaluationDescriptor\EvaluationDescriptorCalculation.m

     文件       5356  2014-04-09 20:24  EvaluationDescriptor\EvaluationDescriptor\functionContrastAt5PerCent.asv

     文件       5356  2014-04-09 20:28  EvaluationDescriptor\EvaluationDescriptor\functionContrastAt5PerCent.m

     文件      41884  2014-04-09 14:50  EvaluationDescriptor\EvaluationDescriptor\Logarithmic image processing.pdf

     文件     437775  2010-02-11 15:58  EvaluationDescriptor\EvaluationDescriptor\Original.pgm

     文件     437775  2010-03-07 23:29  EvaluationDescriptor\EvaluationDescriptor\Restored.pgm

     文件     437775  2010-03-04 17:43  EvaluationDescriptor\EvaluationDescriptor\Restored2.pgm

     文件     437775  2010-03-04 17:44  EvaluationDescriptor\EvaluationDescriptor\Restored3.pgm

     文件     437775  2010-03-04 17:44  EvaluationDescriptor\EvaluationDescriptor\Restored4.pgm

     文件     437775  2010-03-04 17:44  EvaluationDescriptor\EvaluationDescriptor\Restored5.pgm

     文件     488938  2014-04-09 15:15  EvaluationDescriptor\EvaluationDescriptor\The logarithmic image processing model and its applications .pdf

     目录          0  2014-11-06 22:15  EvaluationDescriptor\EvaluationDescriptor

     目录          0  2014-11-06 22:15  EvaluationDescriptor

----------- ---------  ---------- -----  ----

              4311270                    15