局部与全局坐标系互转方法及实例求解（含CRa参数）
问题背景
需要实现局部坐标系（local coordinate system）与全局坐标系（global coordinate system）的双向转换，输入参数为位置坐标+旋转角组合：(X, Y, Z, RX, RY, RZ)。
已了解旋转矩阵与平移矩阵的概念，但不清楚如何结合两者得到最终的X,Y,Z及旋转角，思路混乱。希望用通俗语言解释转换逻辑，补充实例说明，并获取相关研究建议。
更新：问题场景说明
存在两个等比例坐标系：
A：全局坐标系（world system）A'：局部坐标系（local system），坐标轴做截断处理以提升可读性
品红色球体位于A'中，其在A'的坐标与旋转角（简称CRa）为：(1.2, 0.7, 0, 0, 0, 0)；A'在A中的CRa为：(0, 3, 5, 90, 0, 0)。
核心需求
- 正向转换：已知球体在
A'的CRa、A'在A的CRa，求球体在A的CRa； - 反向转换：已知球体在
A的CRa、A'在A的CRa，求球体在A'的CRa。
通俗转换逻辑拆解
坐标系转换本质是「先转方向，再挪位置」（正向），反向则是「先挪回位置，再转反方向」。
1. 位置坐标转换
正向（局部→全局）
把局部坐标系里的点，先跟着局部坐标系的旋转角度转成全局方向的向量，再加上局部坐标系原点在全局的位置，就是点的全局坐标。
公式（矩阵形式）：全局位置 = 局部坐标系的旋转矩阵 × 局部位置向量 + 局部坐标系原点的全局位置
反向（全局→局部）
先把全局位置减去局部坐标系原点的全局位置，得到相对局部原点的向量，再用局部坐标系旋转矩阵的逆矩阵（旋转矩阵是正交矩阵，逆矩阵等于转置矩阵）转成局部方向的向量，就是局部坐标。
公式：局部位置 = 局部坐标系旋转矩阵的逆矩阵 × (全局位置 - 局部坐标系原点的全局位置)
2. 旋转角转换
旋转角的转换本质是旋转的叠加，用旋转矩阵计算更可靠（避免欧拉角的万向锁问题）。
正向（局部→全局）
物体在局部的旋转，叠加局部坐标系相对于全局的旋转，就是物体在全局的旋转。矩阵形式：全局旋转矩阵 = 局部坐标系的全局旋转矩阵 × 物体的局部旋转矩阵
之后再将组合旋转矩阵转换成对应的欧拉角(RX, RY, RZ)。
反向（全局→局部）
物体在全局的旋转，抵消局部坐标系相对于全局的旋转（即乘以局部坐标系旋转矩阵的逆矩阵），得到物体在局部的旋转。公式：物体的局部旋转矩阵 = 局部坐标系旋转矩阵的逆矩阵 × 物体的全局旋转矩阵
再转成欧拉角即可。
实例计算（基于场景条件）
已知条件：
A'在A的CRa：(0,3,5,90,0,0)→ 原点全局位置(0,3,5)，绕X轴旋转90°（RX=90°）- 球体在
A'的CRa：(1.2,0.7,0,0,0,0)→ 局部位置(1.2,0.7,0)，无旋转（旋转矩阵为单位矩阵）
1. 正向转换：求球体在A的CRa
位置计算
绕X轴旋转90°的旋转矩阵（右手坐标系）：
R_X(90°) = [ [1, 0, 0], [0, 0, -1], [0, 1, 0] ]
球体局部位置向量[1.2, 0.7, 0]^T乘以旋转矩阵，得到全局方向向量：[1.2, 0, 0.7]^T
加上A'原点的全局位置，最终全局位置：(1.2+0, 0+3, 0.7+5) = (1.2, 3, 5.7)
旋转角计算
球体局部旋转矩阵为单位矩阵，因此全局旋转矩阵 = R_X(90°) × 单位矩阵 = R_X(90°)，对应欧拉角(90,0,0)。
最终球体在A的CRa：(1.2, 3, 5.7, 90, 0, 0)
2. 反向转换：已知球体在A的CRa为(1.2,3,5.7,90,0,0)，求在A'的CRa
位置计算
全局位置减去A'原点的全局位置：(1.2-0, 3-3,5.7-5) = (1.2,0,0.7)
乘以R_X(90°)的逆矩阵（即转置矩阵，与原矩阵一致），得到局部位置：[1.2,0.7,0]
旋转角计算
物体全局旋转矩阵为R_X(90°)，乘以R_X(90°)的逆矩阵得到单位矩阵，对应欧拉角(0,0,0)。
最终球体在A'的CRa：(1.2,0.7,0,0,0,0)，与初始条件一致。
研究建议
- 优先用旋转矩阵/四元数处理旋转，避免欧拉角的万向锁问题，按需再转成欧拉角；
- 掌握正交矩阵的性质（旋转矩阵都是正交矩阵，逆=转置），简化计算；
- 参考3D图形学入门资料（如《Real-Time Rendering》基础章节），获取直观的转换实例；
- 用特殊角度（如0°、90°）的测试案例验证逻辑，快速排查错误；
- 编程实现时优先使用成熟线性代数库（如Eigen、GLM），避免手动写矩阵计算出错。
内容的提问来源于stack exchange，提问作者Mika
