Great code, thank you a lot Mr Kroon.
I have one problem with the function affine_warp.m :
When I choose a different size for the output image of the image_interpolation.m function I get the error
??? Error using ==> reshape
To RESHAPE the number of elements must not change.
Error in ==> affine_warp>image_interpolation at 241
Error in ==> affine_warp at 64
Is it something that can be solved? I want to display the whole transformed image. In case I use the original size, I loose part of the transformed image.