%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????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