图像拉普拉斯金字塔融合

Halcom

图像拉普拉斯金字塔融合：
【1】http://blog.csdn.net/abcjennifer/article/details/7628655
【2】http://www.cnblogs.com/silence-hust/p/4193208.html
【3】http://blog.csdn.net/yiluoyan/article/details/46386107
【4】http://www.ilovematlab.cn/thread-442912-1-1.html

1. close all
   
2. clear
   
3. imga = im2double(imread('apple1.jpg'));
   
4. imgb = im2double(imread('orange1.jpg')); % size(imga) = size(imgb)
   
5. figure,
   
6. subplot(121),imshow(imga,[])
   
7. subplot(122),imshow(imgb,[])
   
8. 
   
9. imga = imresize(imga,[size(imgb,1) size(imgb,2)]);
   
10. [M,N] = size(imga(:,:,1));
    
11. 
    
12. v = 230;
    
13. level = 5;
    
14. limga = genPyr(imga,'lap',level); % the Laplacian pyramid
    
15. limgb = genPyr(imgb,'lap',level);
    
16. 
    
17. maska = zeros(size(imga));
    
18. maska(:,1:v,:) = 1;
    
19. maskb = 1-maska;
    
20. blurh = fspecial('gauss',30,15); % feather the border
    
21. maska = imfilter(maska,blurh,'replicate');
    
22. maskb = imfilter(maskb,blurh,'replicate');
    
23. 
    
24. limgo = cell(1,level); % the blended pyramid
    
25. for p = 1:level
    
26.         [Mp Np ~] = size(limga{p});
    
27.         maskap = imresize(maska,[Mp Np]);
    
28.         maskbp = imresize(maskb,[Mp Np]);
    
29.         limgo{p} = limga{p}.*maskap + limgb{p}.*maskbp;
    
30. end
    
31. imgo = pyrReconstruct(limgo);
    
32. figure,imshow(imgo) % blend by pyramid
    
33. imgo1 = maska.*imga+maskb.*imgb;
    
34. figure,imshow(imgo1) % blend by feathering

复制代码

pyr_expand( img ):

1. function [ imgout ] = pyr_expand( img )
   
2. %PYR_EXPAND  Image pyramid expansion
   
3. %   B = PYR_EXPAND( A )  If A is M-by-N, then the size of B
   
4. %        is (2*M-1)-by-(2*N-1). Support gray or rgb image.
   
5. %        B will be transformed to double class.
   
6. %        Results the same w/ MATLAB func impyramid.
   
7. 
   
8. kw = 5; % default kernel width
   
9. cw = .375; % kernel centre weight, same as MATLAB func impyramid. 0.6 in the Paper
   
10. ker1d = [.25-cw/2 .25 cw .25 .25-cw/2];
    
11. kernel = kron(ker1d,ker1d')*4;
    
12. 
    
13. % expand [a] to [A00 A01;A10 A11] with 4 kernels
    
14. ker00 = kernel(1:2:kw,1:2:kw); % 3*3
    
15. ker01 = kernel(1:2:kw,2:2:kw); % 3*2
    
16. ker10 = kernel(2:2:kw,1:2:kw); % 2*3
    
17. ker11 = kernel(2:2:kw,2:2:kw); % 2*2
    
18. 
    
19. img = im2double(img);
    
20. sz = size(img(:,:,1));
    
21. osz = sz*2-1;
    
22. imgout = zeros(osz(1),osz(2),size(img,3));
    
23. 
    
24. for p = 1:size(img,3)
    
25.         img1 = img(:,:,p);
    
26.         img1ph = padarray(img1,[0 1],'replicate','both'); % horizontally padded
    
27.         img1pv = padarray(img1,[1 0],'replicate','both'); % horizontally padded
    
28.        
    
29.         img00 = imfilter(img1,ker00,'replicate','same');
    
30.         img01 = conv2(img1pv,ker01,'valid'); % imfilter doesn't support 'valid'
    
31.         img10 = conv2(img1ph,ker10,'valid');
    
32.         img11 = conv2(img1,ker11,'valid');
    
33.        
    
34.         imgout(1:2:osz(1),1:2:osz(2),p) = img00;
    
35.         imgout(2:2:osz(1),1:2:osz(2),p) = img10;
    
36.         imgout(1:2:osz(1),2:2:osz(2),p) = img01;
    
37.         imgout(2:2:osz(1),2:2:osz(2),p) = img11;
    
38. end
    
39. 
    
40. end

复制代码

[ pyr ] = genPyr( img, type, level )

1. function [ pyr ] = genPyr( img, type, level )
   
2. %GENPYR generate Gaussian or Laplacian pyramid
   
3. %   PYR = GENPYR(A,TYPE,LEVEL) A is the input image,
   
4. %        can be gray or rgb, will be forced to double.
   
5. %        TYPE can be 'gauss' or 'laplace'.
   
6. %        PYR is a 1*LEVEL cell array.
   
7. % Yan Ke @ THUEE, xjed09@gmail.com
   
8. 
   
9. pyr = cell(1,level);
   
10. pyr{1} = im2double(img);
    
11. for p = 2:level
    
12.         pyr{p} = pyr_reduce(pyr{p-1});
    
13. end
    
14. if strcmp(type,'gauss'), return; end
    
15. 
    
16. for p = level-1:-1:1 % adjust the image size
    
17.         osz = size(pyr{p+1})*2-1;
    
18.         pyr{p} = pyr{p}(1:osz(1),1:osz(2),:);
    
19. end
    
20. 
    
21. for p = 1:level-1
    
22.         pyr{p} = pyr{p}-pyr_expand(pyr{p+1});
    
23. end
    

复制代码

[ imgout ] = pyr_reduce( img )：

1. function [ imgout ] = pyr_reduce( img )
   
2. %PYR_REDUCE  Image pyramid reduction
   
3. %   B = PYR_REDUCE( A )  If A is M-by-N, then the size of B
   
4. %        is ceil(M/2)-by-ceil(N/2). Support gray or rgb image.
   
5. %        B will be transformed to double class.
   
6. %        Results the same w/ MATLAB func impyramid.
   
7. 
   
8. kernelWidth = 5; % default
   
9. cw = .375; % kernel centre weight, same as MATLAB func impyramid. 0.6 in the Paper
   
10. ker1d = [.25-cw/2 .25 cw .25 .25-cw/2];
    
11. kernel = kron(ker1d,ker1d');
    
12. 
    
13. img = im2double(img);
    
14. sz = size(img);
    
15. imgout = [];
    
16. 
    
17. for p = 1:size(img,3)
    
18.         img1 = img(:,:,p);
    
19.         imgFiltered = imfilter(img1,kernel,'replicate','same');
    
20.         imgout(:,:,p) = imgFiltered(1:2:sz(1),1:2:sz(2));
    
21. end
    
22. 
    
23. end

复制代码

[ img ] = pyrReconstruct( pyr )

1. function [ img ] = pyrReconstruct( pyr )
   
2. %PYRRECONSTRUCT Uses a Laplacian pyramid to reconstruct a image
   
3. %   IMG = PYRRECONSTRUCT(PYR) PYR should be a 1*level cell array containing
   
4. %   the pyramid, SIZE(PYR{i}) = SIZE(PYR{i-1})*2-1
   
5. 
   
6. for p = length(pyr)-1:-1:1
   
7.         pyr{p} = pyr{p}+pyr_expand(pyr{p+1});
   
8. end
   
9. img = pyr{1};
   
10. 
    
11. end

复制代码

linuvzg

感谢楼主，好人一生平安