Matlab registration

%*************************************************************************% %****************** Author: Vlachos Marios ******************% %****************** Designed by: Vlachos Marios ******************% %****************** Electrical and Computer Engineer ******************% %****************** PhD Student ******************% %****************** University of Patras, Greece ******************% %****************** 8 November 2005 ******************% %*************** email: mvlachos@george.wcl2.ee.upatras.gr ***************% %*************************************************************************% clear all; close all; % Start time tic; % Read base image and find its size %************************************************************************** I=imread('P10001.bmp'); I=I(:,:,1); [M N]=size(I); %************************************************************************** % Read transform image %************************************************************************** I_transform=imread('P11111.bmp'); I_transform=I_transform(:,:,1); %************************************************************************** % Initial conditions (scale parameters) %************************************************************************** initial_scale=[pi/12 1 1]; %************************************************************************** % Find the scale parameters that minimizes the function rescale and display % the min value of the function %************************************************************************** [scale,Fval]=fminsearch('rescale2',initial_scale,[],I,I_transform); disp(Fval); %************************************************************************** % Scale parameters %************************************************************************** u=scale(1); dx=scale(2); dy=scale(3); %************************************************************************** % Construct T matrix for rotation and translation by multiplying the two % above matrixs %************************************************************************** Trotation=[cos(u) sin(u) 0; -sin(u) cos(u) 0; 0 0 1;]; Ttranslation=[1 0 0; 0 1 0; dx dy 1;]; T=Trotation*Ttranslation; %************************************************************************** % Align the transform image %************************************************************************** Tform=maketform('affine',T); I_aligned=imtransform(I_transform,Tform,'Xdata',[1 N],'Ydata',[1 M]); %************************************************************************** % t1 is the time of all process t1=toc; % View images on the same figure %************************************************************************** subplot(1,3,1); imshow(I); title('Original Image'); subplot(1,3,2); imshow(I_transform); title('Transformed Image'); subplot(1,3,3); imshow(I_aligned); title('Aligned Image'); %************************************************************************** % View three different figures %************************************************************************** figure, imshow(I); title('Original Image'); figure, imshow(I_transform); title('Transformed Image'); figure, imshow(I_aligned); title('Aligned Image'); %**************************************************************************