Step 2: Making geometric transformation matrix

The second step is making geometric matrix as

t_concord = cp2tform(input_points,base_points,'projective');

cp2tform will find the parameters of the distortion that best fits the stray input_points and base_points you picked.