遥感融合定量评价matlab程序代码
(2006-05-27 21:11:33)
function[DIF,Std,D3,r,c,G]=digitalanalyse(f,l,h)
Dfusion=imread(f);
Dlow=imread(l);
Dhigh=imread(h);
[rh,ch]=size(Dhigh);
Dlowh=imresize(Dlow,[rh,ch],'bicubic');
mean=mean2(Dlow);
meanf=mean2(Dfusion);
DIF=abs(meanf-mean);%求灰度均值差,差值越小,光谱保真性越好
Std=std2(Dfusion);%求标准偏差,标准差越大灰度级分布越分散,目视效果越好
Std2=std2(Dlowh);
Ds=Std^2;%求方差,标准差的平方
Dl=Std2^2;
D3=abs(Ds-Dl);%求差方差,即融合图像相比于原图像信息丢失或增加的多少
p=imhist(Dfusion(:),8)/numel(Dfusion(:));%求直方图
r=entropy(Dfusion(:));
%求熵,熵越大,说明融合图像携带的信息量越大,信息越丰富
c=corr2(Dfusion(:),Dlowh(:));
%求相关系数,反映融合影像同原多光谱影像特征相似程度,即光谱保持性能
h1=diff(Dfusion);%求影像差分,反映图像清晰度
h=mean2(h1);
g1=diff(Dlow);
g=mean2(g1);
G=h-g;%求融合后影像与原影像差分差值,表示融合后影像清晰度提高或降低情况
分享:
喜欢
0
赠金笔
加载中,请稍候......
评论加载中,请稍候...
发评论
登录名: 密码: 找回密码 注册记住登录状态
昵 称:
评论并转载此博文
发评论
以上网友发言只代表其个人观点,不代表新浪网的观点或立场。