Another way to approach this is to compute the LDU decomposition of the rotation matrix. L and U are shears while D is the scaling in the horizontal and vertica

Another way to approach this is to compute the LDU decomposition of the rotation matrix. L and U are shears while D is the scaling in the horizontal and vertical directions.