泊松融合

Halcom

泊松融合--视频分享
百度网盘链接：链接：http://pan.baidu.com/s/1skOTRul

具体链接在halcom.cn论坛，联系人QQ：3283892722
欢迎大家录制视频，你可在论坛进行打赏分享。
视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1


泊松融合参考链接:【1】http://cs.brown.edu/courses/csci1950-g/
【2】http://www.cnblogs.com/CVFans/p/4475885.html

1. clc,clear,close all
   
2. warning off
   
3. lena = double(imread('lena.png'));
   
4. girl = double(imread('girl.png'));
   
5. [Lh,Lv] = imgrad(lena);  % 水平Lh和垂直梯度Lv
   
6. [Gh,Gv] = imgrad(girl);  % 水平Gh和垂直梯度Gv
   
7. 
   
8. X = lena;
   
9. Fh = Lh; % 水平梯度Lh
   
10. Fv = Lv; % 垂直梯度Lv
    
11. 
    
12. % 融合的坐标位置
    
13. w = 57;
    
14. h = 16;
    
15. LX = 123;
    
16. LY = 125;
    
17. GX = 89;
    
18. GY = 101;
    
19. 
    
20. % 贴图融合
    
21. X(LY:LY+h,LX:LX+w,:) = girl(GY:GY+h,GX:GX+w,:);
    
22. Fh(LY:LY+h,LX:LX+w,:) = Gh(GY:GY+h,GX:GX+w,:);
    
23. Fv(LY:LY+h,LX:LX+w,:) = Gv(GY:GY+h,GX:GX+w,:);
    
24. 
    
25. figure(1)
    
26. subplot(231),imshow(uint8(lena),[]);xlabel('lena');
    
27. subplot(232),imshow(uint8(girl),[]);xlabel('girl');
    
28. subplot(233),imshow(uint8(X),[]);xlabel('X');
    
29. subplot(234),imshow(im2uint8(Fh),[]);xlabel('Fh');
    
30. subplot(235),imshow(im2uint8(Fv),[]);xlabel('Fv');
    
31. % mask
    
32. msk = zeros(size(X));
    
33. msk(LY:LY+h,LX:LX+w,:) = 1;
    
34. subplot(236),imshow(msk,[]);xlabel('msk');
    
35. %%
    
36. figure(2)
    
37. tic;
    
38. Y = PoissonJacobi( X, Fh, Fv, msk );
    
39. subplot(121),imshow(uint8(Y),[]);xlabel('PoissonJacobi');
    
40. toc
    
41. 
    
42. tic;
    
43. Y = PoissonGaussSeidel( X, Fh, Fv, msk );
    
44. toc
    
45. subplot(122),imshow(uint8(Y),[]);xlabel('PoissonGaussSeidel');

复制代码

梯度计算：

1. function [Fh Fv] = imgrad(X)
   
2. 
   
3. Kh = [ 0,-1, 1 ];
   
4. Kv = [ 0;-1; 1 ];
   
5. 
   
6. Fh = imfilter(X,Kh,'replicate');
   
7. Fv = imfilter(X,Kv,'replicate');
   

复制代码

Poisson Image Reconstruction by Jacobi Algorithm：

1. function dst = PoissonJacobi(src, Fh, Fv, msk, itr, th, verbose)
   
2. 
   
3. if( nargin < 5 )
   
4.     itr = 1024;
   
5. end
   
6. 
   
7. if( nargin < 6 )
   
8.     th = 1E-3;
   
9. end
   
10. 
    
11. if( nargin < 7 )
    
12.     verbose = false;
    
13. end
    
14. 
    
15. K=[0,1,0;1,0,1;0,1,0];
    
16. p = ( msk > 0 );
    
17. lap = grad2lap(Fh,Fv);
    
18. 
    
19. df0 = 1E32;
    
20. dst = src;
    
21. dst0 = dst;
    
22. for i = 1:itr
    
23.     lpf = imfilter(dst,K,'replicate');
    
24.     dst(p) = (lap(p) + lpf(p))/4;
    
25.    
    
26.     dif = abs(dst-dst0);
    
27.     df = max(dif(:));
    
28.    
    
29.     if( verbose )
    
30.         fprintf('%d %g %g\n',i, df, (df0 - df)/df0);
    
31.     end
    
32.    
    
33.     if( (df0 - df)/df0 < th )
    
34.         break;
    
35.     end
    
36.     dst0 = dst;
    
37.     df0 = df;
    
38. end
    
39. 
    
40. function lap = grad2lap(Fh, Fv)
    
41. lap = circshift(Fh,[0,1]) + circshift(Fv,[1,0]) - Fh - Fv;

复制代码

Halcom

>> A = [1 1 0 0; 1 1 0 0; 0 0 0 0; 0 0 0 0]+1

A =

     2     2     1     1
     2     2     1     1
     1     1     1     1
     1     1     1     1

>> Y1 = circshift(A,[1 1])

Y1 =

     1     1     1     1
     1     2     2     1
     1     2     2     1
     1     1     1     1

>> Y2 = circshift(A,[0 1])

Y2 =

     1     2     2     1
     1     2     2     1
     1     1     1     1
     1     1     1     1

>> Y3 = circshift(A,[1 0])

Y3 =

     1     1     1     1
     2     2     1     1
     2     2     1     1
     1     1     1     1