VTK坐标系转换与vtkTransform代数本质

三维可视化与交互系统（如 CAD、三维人脸扫描、医学图像分割）的核心是几何变换。在 VTK 中，开发者经常遇到的困惑包括：

为什么在屏幕上点击一个点，拿到的坐标不能直接用于三维裁剪？
为什么 vtkTransform 先平移后旋转，与先旋转后平移的结果完全不同？
PreMultiply() 和 PostMultiply() 到底改变了什么？

本文将系统梳理 VTK 的坐标系转换链，并从线性代数角度推导 vtkTransform 的变换本质。

一、VTK 五大坐标系体系

VTK 内部维护了五套层层递进的坐标体系：

1.1 坐标系定义与范围

坐标系名称	维度与单位	坐标原点与范围	作用说明
Model 坐标系	3D (毫米/米)	网格或算法创建时的局部原点。	描述几何体自身形状。
World 坐标系	3D (物理单位)	三维场景统一的绝对原点。	多个 Actor 组合时的空间绝对位置。
View 坐标系	3D (无单位)	位于 Camera 焦点，范围。	归一化视锥体（NDC 坐标）。
Display 坐标系	2D (像素)	视口左下角，右上角。	屏幕鼠标点击与 UI 事件坐标。

二、坐标系转换实战：vtkCoordinate 工具类

在屏幕交互（如鼠标点击选点）中，我们需要将 Display 坐标转换为 World 坐标。

cpp

#include <vtkSmartPointer.h>
#include <vtkCoordinate.h>
#include <vtkRenderer.h>
#include <iostream>

// 将屏幕鼠标点击坐标 (x, y) 转换为世界坐标
void ConvertDisplayToWorld(vtkRenderer* renderer, double displayX, double displayY) {
    auto coordinate = vtkSmartPointer<vtkCoordinate>::New();
    coordinate->SetCoordinateSystemToDisplay();
    coordinate->SetValue(displayX, displayY, 0.0);

    // 绑定当前的 Renderer 进行相机逆投影
    double* worldPos = coordinate->GetComputedWorldValue(renderer);

    std::cout << "World Position: (" 
              << worldPos[0] << ", " 
              << worldPos[1] << ", " 
              << worldPos[2] << ")" << std::endl;
}

三、vtkTransform 的代数本质：PreMultiply vs PostMultiply

在 VTK 中，vtkTransform 用于管理齐次变换矩阵。

假设我们有一个列向量点坐标，变换后的坐标为：

当顺次应用多个变换（平移、旋转、缩放）时，矩阵乘法顺序至关重要。

3.1 PreMultiply (前乘 / 右乘模式 - 默认)

在 PreMultiply() 模式下，新增加的变换矩阵会被乘在已有变换矩阵的右侧：

如果依次调用 RotateZ(a) -> Translate(x, y, z)：

作用于点时：

物理几何意义：变换是基于**局部坐标系（Local Reference Frame）**进行的。每步变换都会随几何体一同旋转/平移轴向。

3.2 PostMultiply (后乘 / 左乘模式)

在 PostMultiply() 模式下，新增加的变换矩阵会被乘在已有变换矩阵的左侧：

如果依次调用 RotateZ(a) -> Translate(x, y, z)：

作用于点时：

物理几何意义：变换是基于**固定世界坐标系（Fixed World Frame）**进行的。坐标轴方向始终固定不变。

四、对比总结表格

模式	矩阵相乘顺序	几何参考框架	典型应用场景
`PreMultiply()`		局部移动坐标系 (Local Frame)	机械臂运动学、机械关节层级旋转
`PostMultiply()`		固定世界坐标系 (World Frame)	三维操纵杆（Manipulator）、视角无关平移

五、三维操纵杆变换代码示范

cpp

#include <vtkSmartPointer.h>
#include <vtkTransform.h>
#include <vtkMatrix4x4.h>

// 构建围绕指定质心 (cx, cy, cz) 旋转的变换矩阵
vtkSmartPointer<vtkTransform> CreateCenterRotation(double angleDegree, double cx, double cy, double cz) {
    auto transform = vtkSmartPointer<vtkTransform>::New();

    // 切换为 PostMultiply 模式，按固定世界坐标系思考
    transform->PostMultiply();

    // 1. 平移质心到原点
    transform->Translate(-cx, -cy, -cz);
    // 2. 绕原点旋转
    transform->RotateZ(angleDegree);
    // 3. 平移回原质心位置
    transform->Translate(cx, cy, cz);

    return transform;
}

六、总结

理解 vtkTransform 的代数本质，关键在于抓住矩阵相乘的左右顺序与几何参考框架的对应关系。在开发三维交互控件时，推荐使用 PostMultiply()，这能使基于世界坐标系的变换逻辑推导更加符合直觉。

VTK坐标系转换与vtkTransform代数本质 ​

一、VTK 五大坐标系体系 ​

1.1 坐标系定义与范围 ​

二、坐标系转换实战：vtkCoordinate 工具类 ​

三、vtkTransform 的代数本质：PreMultiply vs PostMultiply ​

3.1 PreMultiply (前乘 / 右乘模式 - 默认) ​

3.2 PostMultiply (后乘 / 左乘模式) ​

四、对比总结表格 ​

五、三维操纵杆变换代码示范 ​

六、总结 ​