"I would like to create a function that creates a transformation matrix for a region from the original position to the new one."
Essentially is this what you have and you want to solve for T?
T_original * T = T_new
If so,
T = T_new * (inverse T_orignal)