see https://math.stackexchange.com/questions/78137/decomposition-of-a-nonsquare-affine-matrix
def getComponents(normalised_homography):
'''((translationx, translationy), rotation, (scalex, scaley), shear)'''
a = normalised_homography[0,0]
b = normalised_homography[0,1]
c = normalised_homography[0,2]
d = normalised_homography[1,0]
e = normalised_homography[1,1]
f = normalised_homography[1,2]
p = math.sqrt(a*a + b*b)
r = (a*e - b*d)/(p)
q = (a*d+b*e)/(a*e - b*d)
translation = (c,f)
scale = (p,r)
shear = q
theta = math.atan2(b,a)
return (translation, theta, scale, shear)
与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…