图像评价指标PSNRSSIMQNRSAMERGAS

多种常见的图像评价指标峰值信噪比，结构相似性指数，光谱角等PSNR,SSIM,QNR,SAM,ERGAS

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Description: 
%           Quality with No Reference （QNR）. Spectral distortion index. 
% 
% Interface:
%           D_lambda_index = D_lambda（I_FI_MSSp）
%
% Inputs:
%           I_F:                Pansharpened image;
%           I_MS:               MS image resampled to panchromatic scale;
%           S:                  Block size （optional）; Default value: 32;
%           p:                  Exponent value （optional）; Default value: p = 1.
% 
% Outputs:
%           D_lambda_index:     D_lambda index.
% 
% References:
%           [Alparone08]        L. Alparone B. Aiazzi S. Baronti A. Garzelli F. Nencini and M. Selva “Multispectral and panchromatic data fusion assessment without reference“
%                               Photogrammetric Engineering and Remote Sensing vol. 74 no. 2 pp. 193?00 February 2008. 
%           [Vivone14]          G. Vivone L. Alparone J. Chanussot M. Dalla Mura A. Garzelli G. Licciardi R. Restaino and L. Wald 揂 Critical Comparison Among Pansharpening Algorithms? 
%                               IEEE Transaction on Geoscience and Remote Sensing 2014. （Accepted）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function D_lambda_index = D_lambda（I_FI_MSSp）

if （size（I_F） ~= size（I_MS））
    error（‘The two input images must have the same dimensions‘）
end

[NMNb] = size（I_F）;

if （rem（NS） ~= 0）
    error（‘The number of rows must be multiple of the block size‘）
end

if （rem（MS） ~= 0）
    error（‘The number of columns must be multiple of the block size‘）
end

D_lambda_index = 0;
for i = 1:Nb-1
    for j = i+1:Nb 
        band1 = I_MS（::i）;
        band2 = I_MS（::j）;
        fun_uqi = @（bs） uqi（bs.data...
            band2（bs.location（1）:bs.location（1）+S-1...
            bs.location（2）:bs.location（2）+S-1））;
        Qmap_exp = blockproc（band1[S S]fun_uqi）;
        Q_exp = mean2（Qmap_exp）;
        band1 = I_F（::i）;
        band2 = I_F（::j）;
        fun_uqi = @（bs） uqi（bs.data...
            band2（bs.location（1）:bs.location（1）+S-1...
            bs.location（2）:bs.location（2）+S-1））;
        Qmap_fused = blockproc（band1[S S]fun_uqi）;
        Q_fused = mean2（Qmap_fused）;
        D_lambda_index = D_lambda_index + abs（Q_fused-Q_exp）^p;
    end
end
s = （（Nb^2）-Nb）/2;
D_lambda_index = （D_lambda_index/s）^（1/p）;

end

%%%%%%% Q-index on x and y images
function Q = uqi（xy）

x = double（x（:））;
y = double（y（:））;
mx = mean（x）;
my = mean（y）;
C = cov（xy）;

Q = 4 * C（12） * mx * my / （C（11）+C（22）） / （mx^2 + my^2）;  

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2618  2014-10-11 14:57  D_lambda.m
     文件        2801  2014-11-08 16:39  D_s.m
     文件        1455  2014-10-11 14:58  ERGAS.m
     文件        4160  2013-07-17 10:19  img_qi.m
     文件        2623  2015-01-28 08:31  MTF.m
     文件        1987  2014-11-23 11:29  MTF_PAN.m
     文件         147  2014-09-27 15:41  norm_blocco.m
     文件         762  2014-09-27 15:41  odd_extension.m
     文件        1721  2014-09-27 15:45  onions_quality.m
     文件         616  2014-09-27 15:42  onion_mult.m
     文件         719  2014-09-27 15:43  onion_mult2D.m
     文件        1244  2014-10-11 15:05  Q.m
     文件        3047  2015-01-28 08:35  q2n.m
     文件        2030  2014-11-23 11:30  QNR.m
     文件        1606  2014-10-11 15:06  SAM.m
     文件        6377  2013-03-24 15:51  ssim.m