Transformation matrix between two dicom image coordinates

This code is used to calculate transformation matrix between two dicom image sets acquired from an MR scanner using dicom info. Keep in mind that these two image sets should be taken with the same localizer and there should be no motion artifact.
The code first finds the transformation matrix between the two image coordinates and patient coordinate separately. With these two transformation matrices, you can match any point in two image sets.
I used it for TurboFlash and EPI DTI image sets.
Here is the link to Matlab codes: download