摄影测量旋转矩阵计算代码QWK三个参数计算九个方向余弦