请问有没有人有关于 图像拼接 mosaic 的资料？？

grrrgrrr

请问有没有人有关于 图像拼接 mosaic 的资料？？
最好是源程序之类的，谢谢！！

happy

1. % clear all;
   
2. % tic
   
3. % I=imread('finger.jpg');
   
4. % I=I(:,:,1);
   
5. % % I1=im2double(I(1:128,1:128));
   
6. % % I2=im2double(I(31:158,31:158));
   
7. % I1=I(1:128,1:128);
   
8. % I2=I(51:178,51:178);
   
9. % F1=fft2(I1);
   
10. % F2=fft2(I2);
    
11. % F2=conj(F2);
    
12. % F=ifft2(F1.*F2/norm(F1.*F2));
    
13. % F=abs(F);
    
14. % F_max=max(max(F));
    
15. % [m,n]=find(F==F_max)
    
16. % figure,mesh(1:128,1:128,abs(F))
    
17. % toc
    
18. 
    
19. clear all;
    
20. tic
    
21. I1=imread('fingerprint1.bmp');
    
22. I2=imread('fingerprint2.bmp');
    
23. I3=I1(131:258,151:278);
    
24. I4=I2(134:261,201:328);
    
25. F1=fft2(I3);
    
26. F2=fft2(I4);
    
27. F2=conj(F2);
    
28. F=ifft2(F1.*F2/norm(F1.*F2));
    
29. F=abs(F);
    
30. F_max=max(max(F));
    
31. [m,n]=find(F==F_max)
    
32. % figure,mesh(abs(F))
    
33. 
    
34. [Height,Width]=size(I3);
    
35. mos_Height=Height+m-1;
    
36. mos_Width=Width+n-1;
    
37. 
    
38. A=uint8(255*ones(mos_Height,mos_Width));
    
39. 
    
40. for i=1:mos_Height
    
41.     if i<m
    
42.         for j=1:Width
    
43.             A(i,j)=I3(i,j);
    
44.         end
    
45.     elseif i>m-1&i<Height+1
    
46.         for j=1:mos_Width
    
47.             if j<n
    
48.                 A(i,j)=I3(i,j);
    
49.             else
    
50.                 A(i,j)=I4(i-m+1,j-n+1);
    
51.             end
    
52.         end
    
53.     else
    
54.         for j=n:mos_Width
    
55.             A(i,j)=I4(i-m+1,j-n+1);
    
56.         end
    
57.     end
    
58. end
    
59. imshow(A)
    
60.             
    
61. toc

复制代码