im1 = imread('image1.bmp');
info = imfinfo('image1.bmp');
registered = imtransform(im1,t_concord,...
'XData',[1 info.Width], 'YData',[1 info.Height]);
You can do Step 1-4 once. After that, t_concord can be used
repeatably.
The original source of this procedure is written by help of Image Processing Toolbox in MATLAB. Please look and check in detail as you like.