wangwangla 发表于 2007-5-24 13:27

请教:这是一个图像拼接的程序,用相位相关的方法作的。

这是一个图像拼接的程序,用相位相关的方法作的。但是结果出来不是已经哦拼接好的图像,是一个三维的曲线图,不知道这个过程中哪一步错了?哪位高手帮忙看一下指正一下,非常感谢啦。
a3=imread('a.jpg');%读入图像
a4=imread('b.jpg');
a1=a3(1:178,1:178);
a2=a4(51:228,51:228);
subplot(2,2,1),imshow(a1)
subplot(2,2,2),imshow(a2) %显示规定大小的图像
b1=a1(:,:,1); %取R分量
b2=a2(:,:,1);
g1=fft2(b1);%求FFT变换
g2=fft2(b2);
g2=conj(g2);
c=ifft2(g1.*g2/norm(g1.*g2));%求户功率谱,求平移量
c=abs(c);%求绝对值
c_max=max(max(c));%求最大值
=find(c==c_max);%求最大值处的位置
figure,mesh(abs(c))               %创建拼接图,绘二维网线
=size(a1);
mos_height=height+m-1;
mos_width=width+n-1;
a=uint8(255*ones(mos_height,mos_width))%全1数组
k=1;
for i=1:width               %融合
for j=1:heighth
    if a1(i,j)==a2(i+m,j+n)
         if i<k
             k=i;
       a(i,j)=0.5a1(i,j)+a2(i+m,j+n) %显示重叠部分
             l=j;
end
end
for i=1:width                     
for j=1:height                              
    if i>k and j<l
      a(i,j)=a1(i,j) %显示属于a1的部分
end
end
for i=1:width
for j=1:heighth
    if i<(k+m )and j>l
   a(i+m,j+n)=a2(i,j) %显示属于a2的部分
end
end
subplot(2,2,3),imshow(a)
toc

eight 发表于 2007-5-24 14:45

原帖由 wangwangla 于 2007-5-24 13:27 发表 http://www.chinavib.com/forum/images/common/back.gif
这是一个图像拼接的程序,用相位相关的方法作的。但是结果出来不是已经哦拼接好的图像,是一个三维的曲线图,不知道这个过程中哪一步错了?哪位高手帮忙看一下指正一下,非常感谢啦。
a3=imread('a.jpg');%读入 ...

请先看看 置顶贴:聚宝盆

wangwangla 发表于 2007-5-29 16:48

我看过了,这是第二次发了,帮我看一下哈

eight 发表于 2007-5-29 16:50

原帖由 wangwangla 于 2007-5-29 16:48 发表 http://www.chinavib.com/forum/images/common/back.gif
我看过了,这是第二次发了,帮我看一下哈

如果这是你看过的结果,我建议你就以下问题再仔细看看:

1. 你不提供出错信息,如何得知你程序哪里出错了?
2. 你不提供数据文件,别人如何帮你测试?
页: [1]
查看完整版本: 请教:这是一个图像拼接的程序,用相位相关的方法作的。