a3=imread('a.jpg');%讀入圖像<br />
a4=imread('b.jpg'); <br />
a1=a3(1:178,1:178);<br />
a2=a4(51:228,51:228);<br />
subplot(2,2,1),imshow(a1)<br />
subplot(2,2,2),imshow(a2) %顯示規定大小的圖像<br />
b1=a1(:,:,1); %取R分量<br />
b2=a2(:,:,1); <br />
g1=fft2(b1);%求FFT變換<br />
g2=fft2(b2);<br />
g2=conj(g2);<br />
c=ifft2(g1.*g2/norm(g1.*g2));%求戶功率譜,求平移量<br />
c=abs(c);&nbsp;&nbsp;%求絕對值<br />
c_max=max(max(c));&nbsp;&nbsp;%求最大值<br />
[m,n]=find(c==c_max);&nbsp;&nbsp;%求最大值處的位置<br />
figure,mesh(abs(c))&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;%創建拼接圖,繪二維網線<br />
[heighth,width]=size(a1);<br />
mos_height=height+m-1;<br />
mos_width=width+n-1;<br />
a=uint8(255*ones(mos_height，mos_width))%全1數組<br />
k=1;<br />
for i=1:width&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;%融合<br />
&nbsp;&nbsp;for j=1:heighth<br />
&nbsp; &nbsp; if a1(i,j)==a2(i+m,j+n)<br />
&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;if i&lt;k<br />
&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp; k=i;<br />
&nbsp; &nbsp;&nbsp; &nbsp; a(i,j)=0.5a1(i,j)+a2(i+m,j+n) %顯示重疊部分<br />
&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp; l=j;<br />
&nbsp;&nbsp;end<br />
end<br />
for i=1:width&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp; <br />
&nbsp;&nbsp;for j=1:height&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;<br />
&nbsp; &nbsp; if i&gt;k and j&lt;l<br />
&nbsp; &nbsp;&nbsp; &nbsp;a(i,j)=a1(i,j) %顯示屬于a1的部分<br />
&nbsp;&nbsp;end<br />
end<br />
for i=1:width<br />
&nbsp;&nbsp;for j=1:heighth<br />
&nbsp; &nbsp; if i&lt;(k+m )and j&gt;l&nbsp;&nbsp;<br />
&nbsp; &nbsp;&nbsp;&nbsp;a(i+m,j+n)=a2(i,j) %顯示屬于a2的部分<br />
&nbsp;&nbsp;end<br />
end<br />
subplot(2,2,3),imshow(a)<br />
toc<br />
<br />