_______ _ _____ _ _ _ _____ _
|__ __| | / ____| | (_) | | / ____| | |
| | | |__ ___ | | ___ ___ _ __ __| |_ _ __ __ _| |_ ___ | (___ _ _ ___| |_ ___ _ __ ___
| | | '_ \ / _ \| | / _ \ / _ \| '__/ _` | | '_ \ / _` | __/ _ \ \___ \| | | / __| __/ _ \ '_ ` _ \
| | | | | | __/| |___| (_) | (_) | | | (_| | | | | | (_| | || __/ ____) | |_| \__ \ || __/ | | | | |
|_| |_| |_|\___| \_____\___/ \___/|_| \__,_|_|_| |_|\__,_|\__\___||_____/ \__, |___/\__\___|_| |_| |_|
__/ |
|___/
🌍 语言选择 | Language
坐标系(Coord)框架为三维空间中的坐标变换提供直观且强大的解决方案。其核心思想是将坐标系视为一等代数对象，支持乘法、除法等算术运算，从而简化复杂的坐标操作。
该框架还引入了革命性的内禀梯度算子用于高级微分几何计算，能够在不使用传统Christoffel符号的情况下直接计算曲率。
- 🔄 直观的坐标变换：通过简单操作符实现世界坐标系↔局部坐标系转换
- 🏗️ 层次化系统：为机器人学和图形学提供高效的多节点变换
- 🎯 代数运算：坐标系作为数学对象支持 +、-、*、/ 运算
- ⚡ 高性能：为实时应用优化
- 🧮 高级数学：可选的微分几何功能用于研究
- 🌐 多平台：C++库及Python绑定
坐标系的概念可以追溯到勒内·笛卡尔，他寻求用几何学来描述天体运动。然而，他的方法缺乏精确计算所需的精度。在笛卡尔之前很久，早期文明就已经有了类似坐标的参考概念——特别是"世界中心"的概念。
在希腊化时期，托勒密宇宙学将地球置于宇宙的中心，而哥白尼后来将这个中心参考点转移到太阳。这些模型之间的关键区别不仅仅是原点的选择，而是它们所启用的数学框架。通过将中心重新定位到太阳，科学家们认识到需要一个动态的、基于运动的数学-物理系统。
因此，坐标系的选择深刻影响世界观和计算范式。例如，爱因斯坦的相对论理论可以被视为将坐标从平坦欧几里得空间扩展到弯曲流形的结果。
如果我们研究微分几何，在微分元素的层面，坐标系可以被线性化。这样，就可以形成某种通用坐标系对象的概念，用作标尺，或者可以理解为空间的维度。
坐标系（或标架），在这里称为Coord对象，是一个根植于群论的数学构造，它定义了三维空间中的坐标系。在物理学中，这样的结构通常被称为参考系，而在微分几何中，它通常被称为标架场或运动标架。
从群论角度来看，坐标系可以被视为能够参与群运算的代数对象。Coord对象统一了这些概念，允许坐标系和单个坐标都作为代数运算中的元素，如乘法和除法。
通过为坐标系扩展算术运算，Coord对象实现直接的坐标变换，无需复杂的矩阵操作。这种方法提供直观的几何解释，简化实际应用中的坐标系操作。
┌─────────────────┐
│ coord3 │ ← 完整坐标系（位置+旋转+缩放）
│ (完整型) │
├─────────────────┤
│ vcoord3 │ ← 向量坐标系（旋转+缩放）
│ (缩放型) │
├─────────────────┤
│ ucoord3 │ ← 单位坐标系（仅旋转）
│ (旋转型) │
└─────────────────┘在C++中，三维空间中的坐标系由原点、三个方向轴和三个缩放分量定义如下：
struct coord {
vec3 ux, uy, uz; // 三个基向量
vec3 s; // 缩放
vec3 o; // 原点
};可以使用三个轴或欧拉角构造坐标系：
coord C1(vec3 o); // 仅位置
coord C2(vec3 ux, vec3 uy, vec3 uz); // 从基向量构造
coord C3(vec3 o, vec3 s, quat q); // 完整形式：位置、缩放、旋转- 组合 (*):
C₃ = C₂ ∘ C₁- 顺序变换 - 相对 (/):
R = C₁ · C₂⁻¹- 相对变换 - 逆向 (%):
R = C₁⁻¹ · C₂- 逆相对变换
坐标变换的基本操作：
从局部坐标系变换到父坐标系：
V0 = V1 * C1 // V1在局部C1 → V0在父系统中从父坐标系投影到局部坐标系：
V1 = V0 / C1 // V0在父系统 → V1在局部C1系统中1. 世界坐标系 ↔ 局部坐标系变换
在世界坐标系和局部坐标系之间转换向量：
VL = Vw / C // 世界坐标到局部坐标
Vw = VL * C // 局部坐标到世界坐标2. 多节点层次化系统
这在机器人学、计算机图形学和复杂机械系统中至关重要：
// 正向变换链
C = C3 * C2 * C1
Vw = VL * C
// 反向变换链
VL = Vw / C
V4 = V1 / C2 / C3 / C43. 平行坐标系转换
在共享同一父系统的坐标系之间转换：
C0 { C1, C2 } // C1和C2都是C0的子系统
V2 = V1 * C1 / C2 // 从C1空间转换到C2空间4. 高级操作
标量乘法：
C * k = {C.o, C.s * k, C.u}
其中: C.u = {C.ux, C.uy, C.uz}四元数操作：
C0 = C1 * q1 // 应用四元数旋转
C1 = C0 / q1 // 移除四元数旋转
q0 = q1 * C1 // 从坐标系提取四元数
q1 = q0 / C1 // 相对四元数向量加法：
C2 = C1 + o // 平移坐标系
其中 C2 = {C1.o + o, C1.v}, C1.v = {C1.ux*C1.s.x, C1.uy*C1.s.y, C1.uz*C1.s.z}- 3D场景图：高效的父子变换
- 角色动画：骨骼层次和骨骼动画
- 相机系统：视图和投影矩阵管理
- 对象定位：直观的位置和方向设置
- 正/逆运动学：关节链计算
- 多臂协调：机器人臂之间的相对定位
- SLAM系统：地图坐标变换
- 路径规划：坐标空间导航
- 装配系统：组件定位和约束
- 机械设计：零件关系和公差
- 制造业：刀具路径坐标变换
- 仿真：物理系统坐标管理
- 玩家移动：角色控制器变换
- 物理系统：碰撞检测坐标空间
- UI系统：屏幕到世界坐标转换
- 网络：同步坐标系统
#include "coord.hpp"
int main() {
// 创建坐标系
coord3 world; // 世界坐标系
coord3 robot(0, 0, 1); // 机器人位于位置(0,0,1)
coord3 arm = coord3::from_eulers(45, 0, 0); // 机械臂绕X轴旋转45°
// 创建变换链
coord3 arm_in_world = robot * arm;
// 将点从机械臂空间变换到世界空间
vec3 point_in_arm(1, 0, 0); // 机械臂坐标系中的点
vec3 point_in_world = point_in_arm * arm_in_world;
// 从世界空间转换回机械臂空间
vec3 converted_back = point_in_world / arm_in_world;
return 0;
}// 具有多个关节的机器人臂
coord3 base;
coord3 shoulder(0, 0, 0.5); // 肩关节
coord3 elbow(0, 0, 0.3); // 肘关节
coord3 wrist(0, 0, 0.2); // 腕关节
// 创建变换链
coord3 full_transform = base * shoulder * elbow * wrist;
// 末端执行器上的点
vec3 end_point(0, 0, 0.1);
vec3 world_position = end_point * full_transform;
// 逆向：给定世界位置找到局部坐标
vec3 local_coords = world_position / full_transform;编译坐标系库的C++版本：
# 基本编译
g++ -std=c++17 -I src -o your_program your_program.cpp
# 示例编译和测试
g++ -std=c++17 -I src -o test_coord test_coord.cpp
./test_coord
# 生产环境优化编译
g++ -std=c++17 -O3 -I src -o optimized_program your_program.cpp要在Python(3.13)中使用coordinate_system，可以通过pip轻松安装：
pip install coordinate_systemfrom coordinate_system import vec3, quat, coord3
# 创建坐标系
a = coord3(0, 0, 1, 0, 45, 0) # 位置和旋转
b = coord3(1, 0, 0, 45, 0, 0)
# 组合变换
a *= b
print(a)以下部分描述用于研究和专业应用的高级数学功能。
该框架为微分几何引入了革命性的内禀梯度算子：
G_μ ≈ (c(u+h) - c(u)) / h / c(u)
其中：
c(u,v)：曲面上的内禀标架（正交归一化的切向量和法向量）h：有限差分步长G_μ：内禀梯度算子，方向μ ∈ {u,v}
几何意义：
G_μ.ux, G_μ.uy：切空间的旋转和形变G_μ.uz：法向量的变化率，即 ∂n/∂μG_μ.s：度量缩放因子的变化
使用内禀梯度算子，第二基本形式可以直接计算：
L = -⟨G_u·n, r_u⟩
M = -½(⟨G_u·n, r_v⟩ + ⟨G_v·n, r_u⟩)
N = -⟨G_v·n, r_v⟩
高斯曲率：
K = (L·N - M²) / det(g)
黎曼曲率张量对应的标架可以通过内禀梯度算子的非对易性计算：
Ω_uv = [G_u, G_v] - G_[u,v]
其中：
[G_u, G_v] = G_u ∘ G_v - G_v ∘ G_u（李括号/交换子）G_[u,v]：李导数项（在坐标基下为零）
黎曼曲率张量分量：
R_{ijkl} = ⟨Ω(eₖ, eₗ) eⱼ, eᵢ⟩
对于曲面，关键分量为：
R_{1212} = -√det(g)·⟨Ω_uv r_v, r_u⟩
从黎曼曲率张量中提取高斯曲率：
K = R_{1212} / det(g) = -⟨Ω_uv r_v, r_u⟩ / √det(g)
其中：
r_u = ∂r/∂u,r_v = ∂r/∂v：坐标基向量det(g) = ‖r_u‖²‖r_v‖² - ⟨r_u, r_v⟩²：度量张量行列式√det(g)：面积元缩放因子
球面验证（理论曲率 K=1）：
- 计算得：
⟨Ω_uv r_v, r_u⟩ = -sinθ √det(g) = sinθK = -(-sinθ)/sinθ = 1✓
圆柱面验证（理论曲率 K=0）：
- 计算得：
⟨Ω_uv r_v, r_u⟩ = 0 K = 0/√det(g) = 0✓
| 曲率类型 | 传统方法 | 内禀梯度方法 | 加速比 |
|---|---|---|---|
| 高斯曲率 | O(n⁴) | O(n³) | ~8× |
| 黎曼张量 | O(n⁶) | O(n³) | ~27× |
| 完整曲率分析 | O(n⁶) | O(n³) | ~27× |
该方法在测试曲面上达到**< 0.001%误差**，消除了对复杂Christoffel符号的需求，同时提供完整的黎曼曲率信息。
验证结果：
- ✅ 球面：K = 1/R²（0.0004%误差）
- ✅ 环面：正/负曲率区域正确识别
- ✅ 比安基恒等式：在数值精度内满足
- ✅ 坐标无关性：结果与参数化无关
圆锥面测试案例 - 使用解析解验证：
- 半顶角：θ = 60°
- 联络系数误差 < 2.54×10⁻⁵
- 曲率分量误差 < 2.54×10⁻⁵
- 完美零标量曲率（符合预期）
框架自然嵌入李理论：
- 群乘法：SE(3)运算
- 李代数：切空间运算
- 括号运算：
[C1, C2] = C1*C2 - C2*C1 - 指数映射：代数到群的变换
本项目采用MIT许可证 - 详情请参见LICENSE文件。
https://zenodo.org/records/14435613
- 历史几何学先驱：笛卡尔、高斯、黎曼
- 现代微分几何社区
- 开源计算几何项目
- 所有贡献者和合作者
Coord框架为以下提供了统一的计算语言：
- 日常坐标变换（主要用例）
- 层次化系统管理（机器人学、图形学）
- 高级微分几何（研究应用）
- 物理参考系（物理仿真）
通过将坐标系视为代数对象，它创建了一种镜像数学推理的直观语法，同时保持计算效率。这种方法在抽象数学和实际实现之间架起了桥梁，使复杂的坐标操作对多个领域的开发者都可以轻松使用。
🌟 如果您觉得有用，请为这个仓库点星！