mosaic.m

图像对齐（inage registration）工具箱

function [rep1, rep2]=mosaic(im1,im2,sz);
% creates the mosaic from the images im1 and im2, each block sized sz x zs
% for the registration quality evaluation
%
% [rep, rep2]=mosaic(im1,im2,sz);

vel=size(im1);
st=floor(vel/sz);

pom1=ones(sz,sz);
mask1=~pom1;

for i=1:(st(2)-1),
   mask1=[mask1 pom1];
   pom1=~pom1;
end;

mask=~mask1;

for i=1:(st(1)-1),
   mask=[mask;mask1];
   mask1=~mask1;
end;

mask(vel(1),vel(2))=0;

rep1=im1.*mask+im2.*~mask;
rep2=im1.*~mask+im2.*mask;