Question

这段基于区域方差的小波图像融合的matlab程序帮看下错，谢了。。

编译的时候第5行就开始有错，无语了，帮下，万分感谢，图像是lena的那个灰度图像
clear all;
close all;
x=imread('1.bmp');
y=imread('2.bmp');
a=rgb2gray(x);
b=rgb2gray(y);
a=double(a);
b=double(b);
[aA,aH,aV,aD]=dwt2(a,'bior2.4');
[bA,bH,bV,bD]=dwt2(b,'bior2.4');
newA=zeros(size(aA));
newH=zeros(size(aH));
newV=zeros(size(aV));
newD=zeros(size(aD));
[m,n]=size(aA);
for i=1:m
for j=1:n
newA(i,j)=(aA(i,j)+bA(i,j))/2;
end;
end;
fun=inline('var(x(:))');
Var_aH=nlfilter(aH,[3 3],fun);
Var_aV=nlfilter(aV,[3 3],fun);
Var_aD=nlfilter(aD,[3 3],fun);
Var_bH=nlfilter(bH,[3 3],fun);
Var_bV=nlfilter(bV,[3 3],fun);
Var_bD=nlfilter(bD,[3 3],fun);
for i=1:m
for j=1:n
if Var_aH(i,j)>=Var_bH(i,j);
newH(i,j)=aH(i,j);
else
newH(i,j)=bH(i,j);
end
end
end
for i=1:m
for j=1:n
if Var_aV(i,j)>=Var_bV(i,j);
newV(i,j)=aV(i,j);
else
newV(i,j)=bV(i,j);
end
end
end
for i=1:m
for j=1:n
if Var_aD(i,j)>=Var_bD(i,j);
newD(i,j)=aD(i,j);
else
newD(i,j)=bD(i,j);
end
end
end
new=idwt2(newA,newH,newV,newD,'bior2.4');
new=uint8(new);
figure;
imshow(new);
title('区域方差取大')

Answer 1

小波分析后，你也没处理，组合起来依然是原信号啊。方差不会变的。

你用连续小波变换，得到了1：43尺度下的小波基系数组z（i）。

如果你指的是在一定分辨率下，该信号在某个尺度上投射的分量的方差。

就用z（i）乘以这个小波φ（τ，σ）基。在用STD求方差呗。
另外,站长团上有产品团购,便宜有保证