点击模型表面，标记却偏在一旁；普通屏幕正常，换成高分屏后选点位置成倍偏移；二维切片点中了，三维视图里的十字线却不对应。这些现象经常来自同一个根源：多个坐标空间之间缺少明确契约。本文不把拾取简化成一个 API 调用，而是沿着鼠标事件到物理点、再到体素索引的路径逐步解释。
**阅读范围：**面向三维可视化、医学影像与桌面软件开发者。本文使用合成数据与工程示例，资料核对日期为 2026 年 10 月 16 日。所述软件功能不等于实际设备精度；本文不涉及患者数据、诊疗建议或临床合格阈值。
图 1：鼠标点到哪里才算选中？Qt 与 VTK 的拾取坐标、DPI 与体素定位的核心关系（原创工程示意）
一、先列出交互链路中的坐标空间
Qt 控件事件通常提供控件相关的位置，渲染器使用自己的显示坐标，场景对象可能具有模型坐标与世界变换，影像还具有物理空间和连续索引。相同的三个数字在不同空间中含义完全不同，变量应写明来源与单位，避免统一叫 pos 后在函数间隐式转换。
建议接口明确使用 logical_position、display_position、world_point 和 continuous_index 等名称。只有在边界处转换，内部尽量保持同一空间。出现偏移时，记录各层数值比反复调整一个"修正系数"更容易发现根因，也能避免为一台显示器加补丁后破坏另一台。
二、逻辑像素与设备像素不能重复缩放
高 DPI 环境中，逻辑像素和实际绘制像素可能通过 devicePixelRatio 关联。Qt 与 VTK 的适配组件也可能已经完成部分转换。如果应用层再乘一次比例，就会产生双倍缩放。正确策略是确认所用适配接口的输入约定，而不是看到高分屏就统一乘 DPR。
本文给出的 display_y 公式只是明确约定下的示例：假定控件高度与鼠标纵坐标以逻辑像素表示，并需要转换到另一种显示坐标。它没有测试真实 Qt 事件、高 DPI 屏幕和适配组件，因此不能直接作为所有 QVTK 控件的通用转换代码。
三、纵轴原点和像素中心需要一致
某些界面坐标从左上角开始，渲染显示坐标通常有另一种纵向约定。转换时是否使用 height-1，还涉及表示像素中心还是连续边界。只在中间点击测试可能看不出差一像素，应该覆盖四角、边界和非整数 DPR。
控件边框、工具栏和布局边距也会影响事件位置。如果拿窗口坐标直接传给子控件的渲染器，即使 DPI 为一也可能出现固定偏移。先通过 Qt 的坐标映射获得目标控件局部位置，再按适配组件契约传递，避免手工减去一个看似稳定的工具栏高度。
四、一个二维屏幕点不能唯一确定三维点
透视场景中，同一屏幕位置对应一条射线。需要深度、与几何的交点或与指定平面的交点，才能确定三维位置。直接把显示坐标的 z 写成零，通常只得到某个特定深度平面上的点，并不意味着选中了模型表面。
因此拾取前要先定义语义：选择最前方可见表面、最近顶点、某个单元，还是影像平面中的位置。不同拾取器和算法返回的信息不同。界面上标记看起来接近并不能证明选中的是同一种对象，应检查返回的对象、单元或点标识。
五、点拾取与单元拾取回答不同问题
点拾取通常关注离射线足够近的顶点，结果可能随着网格采样密度变化；单元拾取则可以返回几何单元上的交点。用户希望在三角形内部任意选点时，直接吸附最近顶点可能产生明显跳动。容差的单位与作用范围也应按对应 API 理解。
对多个对象共存的场景，应限制可拾取对象集合，避免隐藏辅助线或透明外壳抢先被选中。没有命中时要清除或保持旧选择，并在交互规则中明确。不能继续使用上一次成功点，却让用户以为当前点击已经生效。
六、对象变换决定返回点怎样解释
actor 可能具有平移、旋转和缩放，源网格仍保留模型坐标。拾取得到的世界点与原始点数组不一定在同一空间。若需要访问模型属性或保存标注，应明确是否要使用对象变换的逆，将世界点恢复到对应局部坐标。
非均匀缩放还会影响法向和距离解释。把世界点直接当局部点存入文件，重新加载到另一变换状态下就会偏移。标注数据最好同时记录参考对象和坐标空间，避免"这个点在哪个模型上"只能依靠当前界面状态推断。
七、物理点到体素索引要考虑方向
vtkImageData 的原点、间距与方向矩阵共同决定索引到物理点的映射。仅计算 (point-origin)/spacing 只适用于特定方向条件。本文加入绕 Z 轴的方向变化，再通过 VTK 接口往返转换，检查与显式矩阵计算一致。
示例索引为 (4,5,6)，间距为 (0.7,0.8,2)，原点为 (10,20,30)，方向交换并翻转面内轴，得到物理点 (6,22.8,42)。这个非平凡结果有助于发现遗漏方向的实现；若输入方向为单位阵，错误公式可能在测试中一直通过。
八、连续索引和整数体素要分开
拾取位置通常落在体素之间，连续索引包含小数。显示灰度时可以插值，选择离散标签时可能采用最近邻；两者不能通过随意取整统一处理。取整后还要检查有效范围，尤其是边界附近与负值位置，不能让 Python 负索引悄悄指向数组末尾。
整数索引的排列也要与数组布局一致。VTK 接口使用 XYZ，而 NumPy 三维数组常按 ZYX 组织。完成物理逆变换后，应在单一位置明确交换索引，避免在鼠标层、读取层与显示层各交换一次，最终只在部分切片上碰巧正确。
九、二维切片拾取可利用已知平面
在 MPR 视图中，目标通常位于当前切片平面，可以通过显示点与平面的关系恢复物理位置。平面原点、基向量和输出间距必须与实际重切一致。若界面只保存层号，却允许斜切或缩放，后续十字线联动很难保持正确。
一个可靠状态是保存公共物理点，再让各视图投影显示。鼠标只是修改这个状态的输入方式。这样三维标记、二维十字线和数值面板可以共享同一位置，减少各窗口维护独立索引而逐渐偏离的问题。
十、拾取测试应该覆盖交互边界
至少包括没有命中、被遮挡对象、边缘点击、多个视口、控件缩放与显示器 DPI 变化。若窗口跨屏移动，比例可能变化，不能只在程序启动时读取一次。视口偏移也需要纳入转换，不能把整窗坐标直接用于右侧子视图。
本篇实际运行的代码只验证影像物理点与连续索引的往返，不包含屏幕拾取。把测试范围写清楚有助于下一步补齐验证：先保证几何运算，再通过真实窗口中的已知屏幕投影点检查事件链路，不能把其中一项成功当成全部成功。
十一、测量前先确认选点语义
两点距离应该在同一个具有明确单位的空间中计算。屏幕像素距离随缩放改变，模型坐标距离可能受 actor 缩放影响，物理空间距离才具有约定的单位含义。界面显示更多小数位不会让选点更准确，也不能证明真实设备精度。
对于保存的测量，应记录两个点、空间、参考数据版本和计算规则。若后来模型发生变换或数据被重采样，应用应明确重新映射还是保持原参考。这样工程人员能够解释一个数值从哪里来，而不是只保留最终标签文本。
完整示例与执行方式
本系列使用独立 Python 3.11 环境，实际安装 VTK 9.5.2、PySide6 6.10.1，以及示例需要的 NumPy、SciPy、NiBabel、pydicom。各篇只需安装实际导入的依赖。将下方代码保存为 demo.py 执行；代码涉及的尺寸按正文约定解释，不应直接作为真实系统的参数。
import numpy as np
from vtkmodules.vtkCommonDataModel import vtkImageData
from vtkmodules.vtkCommonMath import vtkMatrix3x3
image=vtkImageData();image.SetDimensions(30,40,20)
image.SetSpacing(.7,.8,2.);image.SetOrigin(10,20,30)
D=np.array([[0.,-1.,0.],[1.,0.,0.],[0.,0.,1.]])
m=vtkMatrix3x3()
for i in range(3):
for j in range(3):m.SetElement(i,j,D[i,j])
image.SetDirectionMatrix(m)
index=[4.,5.,6.];point=[0.,0.,0.]
image.TransformContinuousIndexToPhysicalPoint(index,point)
back=[0.,0.,0.];image.TransformPhysicalPointToContinuousIndex(point,back)
np.testing.assert_allclose(back,index)
expected=np.array([10,20,30])+D@(np.array([.7,.8,2])*index)
np.testing.assert_allclose(point,expected)
widget_height=700;logical_y=120;dpr=1.5
vtk_y=(widget_height-1-logical_y)*dpr
assert vtk_y>=0
print('world_point_mm',point,'recovered_index',back,'illustrative_display_y',vtk_y)
print('PASS: physical/index inverse; display formula is an explicit example convention')补充实验：独立数值检查
下面的补充脚本已独立运行通过。它把正文中的关键关系转成小型断言，便于区分数据处理和显示表现；保存为 extra.py 后运行。断言覆盖范围以代码为准，不代表完整项目已经验证。
import numpy as np
origin = np.array([10.,20.,30.])
spacing = np.array([.7,.8,2.])
direction = np.array([[0.,-1,0],[1,0,0],[0,0,1.]])
shape_xyz = np.array([30,40,20])
def world_to_index(world):
return np.linalg.solve(direction @ np.diag(spacing), np.asarray(world)-origin)
def nearest_valid_voxel(continuous):
if not np.isfinite(continuous).all():
raise ValueError('Nonfinite index')
# Example policy: reject outside center bounds, then round to nearest center.
if np.any(continuous < 0) or np.any(continuous > shape_xyz-1):
raise ValueError('Outside image center bounds')
return np.floor(continuous+0.5).astype(int)
for index in [np.array([0.,0.,0.]),np.array([4.,5.,6.]),shape_xyz-1]:
world = origin + direction @ (spacing * index)
recovered = world_to_index(world)
np.testing.assert_allclose(recovered,index,atol=1e-12)
np.testing.assert_array_equal(nearest_valid_voxel(recovered),index)
for bad in [np.array([-1.,2,3]),np.array([30.,0,0]),np.array([np.nan,0,0])]:
try:
nearest_valid_voxel(bad)
except ValueError:
pass
else:
raise AssertionError('Invalid index accepted')
print('PASS: world/index checks and explicit out-of-bounds rejection')本机执行结果与验证口径
world_point_mm [6.0, 22.8, 42.0] recovered_index [4.0000000000000036, 5.0, 6.0] illustrative_display_y 868.5
PASS: physical/index inverse; display formula is an explicit example convention以上为本任务中的真实执行输出。断言检查的是代码明确表达的性质；不能由一次通过推断其他平台、其他输入或完整业务流程已经验证。示例若说明未渲染、未编译或采用合成文件，应按该限制解释结果。
图 2：常见现象、可能原因与检查行动（原创工程示意）
将示例接入项目：三条执行路径
先建立影像坐标的可逆基准
运行本文脚本，比较 VTK 接口和显式矩阵结果。更换原点、非等方间距及合法方向，检查已知点往返。再加入体外位置与边界取整，明确何时拒绝、何时夹取，避免在真实鼠标事件中才发现数组越界。
再检验真实控件输入约定
在目标 QVTK 适配方式下记录事件局部坐标、DPR、渲染窗口尺寸与视口。对已知世界点投影后的屏幕位置执行交互检查，覆盖四角和多屏移动。确认适配组件完成哪些转换，再决定应用层是否需要额外处理。
最后对齐产品交互语义
明确表面交点、最近顶点、影像平面或标签体素等不同选择模式。限制可拾取对象并处理未命中状态，保存标注时附带坐标空间与参考对象。把几何正确、界面响应和业务选择规则分别验收。
排错记录应如何写？
高分屏点击成倍偏移
可能 DPR 被重复应用或使用了错误尺寸单位。
记录适配接口输入约定与窗口尺寸，检查转换仅发生一次。
所有点击固定偏移
可能拿窗口坐标代替控件局部坐标。
使用明确的控件坐标映射，检查边距和视口原点。
世界点正确但体素错误
可能忽略方向矩阵或混用 XYZ 与 ZYX。
用物理点往返及已知非对称像素验证索引链路。
点击空白仍显示旧结果
可能未命中时继续使用上一条有效拾取。
明确未命中状态，更新界面提示和标记规则。
适用边界与独立验证
本篇实际验证带方向矩阵的 vtkImageData 连续索引与物理点往返；显示纵坐标仅为说明约定的计算示例。没有运行真实鼠标拾取、高 DPI、多视口和 actor 变换场景，这些是接入 GUI 后必须独立检查的部分。
建议保留原始输入、数据单位、参数、软件版本和失败日志。检查数据如果被反复用于选择实现或调参，就不应继续作为完全独立的最终验证集。对外报告应区分功能通过、几何正确、性能测量与环境兼容性，避免用一个成功截图代替全部证据。
参考资料
写在最后
"千里之行，始于足下。"
------ 老子《道德经》第六十四章
把每一步的输入、状态和证据写清楚，比凭经验反复尝试更容易积累可靠的工程判断。