def METHOD_NAME(src, dst):
    assert np.array_equal(src.shape, dst.shape) and len(src.shape) == 2, \
        '2d input arrays are expected, got {}'.format(src.shape)
    src_col_mean, src_col_std = FaceIdentifier.normalize(src, axis=0)
    dst_col_mean, dst_col_std = FaceIdentifier.normalize(dst, axis=0)
    u, _, vt = np.linalg.svd(np.matmul(src.T, dst))
    r = np.matmul(u, vt).T
    transform = np.empty((2, 3))
    transform[:, 0:2] = r * (dst_col_std / src_col_std)
    transform[:, 2] = dst_col_mean.T - np.matmul(transform[:, 0:2], src_col_mean.T)
    return transform