《Results in Engineering》：Image-based kinematic calibration of a dual-robot x-ray computed tomography system via reprojection error optimization
编辑推荐：
摘要：本文提出一种基于图像的双机械臂X射线计算机断层扫描（Computed Tomography，CT）系统运动学标定方法，该系统配有附加旋转工作台（turntable），无需外部测量仪器。利用已知几何形状的球形标记物（spherical-marker obj
摘要：本文提出一种基于图像的双机械臂X射线计算机断层扫描（Computed Tomography，CT）系统运动学标定方法，该系统配有附加旋转工作台（turntable），无需外部测量仪器。利用已知几何形状的球形标记物（spherical-marker object）的放射投影图像，在平滑有界约束下通过最小化多视角重投影误差（multi-view reprojection error），联合估计两条机械臂运动学链及转台的修正参数。所得修正后的正运动学模型（forward-kinematic model）可生成单一、全局一致的投影几何（projection geometry），无需后续逐投影细化。该方法利用重建数据的调制传递函数（Modulation Transfer Function，MTF）测量进行定量验证，并与基于直接线性变换（Direct Linear Transformation，DLT）的逐投影几何校正进行比较。结果表明，所提方法对经典步进—停止（Stop-and-Go）和连续FlyBy CT采集均适用。
论文解读——《Image-based kinematic calibration of a dual-robot x-ray computed tomography system via reprojection error optimization》
该文由Anton Weiss完成，发表于《Results in Engineering》。工业CT传统上采用固定几何结构，而新兴的双机械臂（twin-robotic）X射线CT系统利用两台六自由度（6-DOF）工业机器人分别搭载X射线源与平板探测器（flat-panel detector），并可集成旋转工作台实现灵活视角的三维CT采集。然而工业机械臂的绝对定位精度受关节间隙、连杆加工与装配误差、线性导轨非线性等因素限制，其固有位姿偏差会直接导致投影几何不一致，进而在重建体中引入伪影（artifact）与分辨率下降。现有校正手段存在明显局限：基于外部测量设备（激光跟踪仪laser tracker、摄影测量photogrammetry）成本高且有视线遮挡问题；基于直接线性变换（DLT）的逐投影校正需重复采集同一位姿且无法适用于连续FlyBy模式，插补DLT矩阵还会引入非物理偏移；基于标志物（fiducial marker）的方法在试样高衰减或遮挡时失效。因此亟需一种能利用系统自身成像数据、不依赖外部仪器、适用于任意轨迹（含FlyBy）的原位（in-situ）图像标定方法。本文针对此需求提出了一种基于已知几何球基准体重投影误差最小化的双机械臂+转台联合运动学标定框架，并通过MTF定量评估其有效性。
主要关键技术方法
研究人员在双7-DOF KUKA KR120 R2900 extra机械臂+线性导轨（KUKA KL4000）+WEISS CR1000C转台系统上，使用Comet VarioFocus微焦点X射线源（225 kV）与iRay NDT1717M3平板探测器（3072×3072，像素间距139 μm）。设计螺旋排列钢球基准体（约φ150 mm×200 mm高）作为标定物，采集正弦调制圆轨迹240幅投影。初始通过多视角射线相交加权融合法重建球标记三维坐标，再经正交Procrustes刚体配准对齐至已知CAD模型。将各关节连杆间的微小平移与旋转偏差建模为轴角（axis-angle）形式修正量，构建含平滑有界惩罚项（smooth bounded limiter）的非线性最小二乘目标函数，以多视图球标记检测椭圆中心与正运动学模型合成投影间的重投影误差（reprojection error）为残差，采用JAX自动微分与Adam优化器全批量联合优化两机械臂及转台修正参数，获得修正后正运动学模型用于任意轨迹投影几何生成。重建采用SART算法（各向同性体素100 μm），并以玻璃珠模体的调制传递函数（MTF）量化空间分辨率，与DLT逐投影校正对比。
研究结果
5.1. Calibration-Trajectory Results（标定轨迹结果）
研究人员在标定采集轨迹上计算优化前后重投影残差。结果显示重投影误差RMSE由优化前3.64 mm降至0.38 mm，残差均值趋近零、标准差显著缩小，证明优化有效提升了该轨迹的几何一致性。
5.2. Experimental Setup（实验装置）
介绍了双机械臂CT硬件配置、光源与探测器参数及SART重建设置，并说明激光跟踪仪安装时测得的TCP位置误差标准差约为1.24–1.39 mm，确认系统存在需校正的几何偏差。
5.3. CT Trajectory Design and Acquisition（CT轨迹设计与采集）
设计了四种采集模式——步进—停止圆轨迹（SNG-C）、连续FlyBy圆轨迹（FB-C）、带方波垂直调制的FlyBy振荡轨迹（FB-O）、带正弦垂直调制的FlyBy正弦轨迹（FB-S），各采集3000幅投影，用于后续MTF对比评估。
5.4. Runtime and Convergence（运行时间与收敛性）
在NVIDIA RTX A4500上全批量优化约2分钟收敛（JIT编译68 s），损失曲线在早期经历短暂震荡后稳定下降至停止阈值10-6以下，表明算法具有良好的数值稳定性与实用性。
5.5. Error Limit Compliance（误差限合规性）
优化所得回转关节平移修正≤±100 μm、转角≤±2 mrad，移动副平移≤±1 mm，均在预设物理合理边界内，说明平滑有界约束有效避免了非物理解。
5.6. Modulation Transfer Function（调制传递函数）
以玻璃珠模体边缘扩散函数（ESF）→线扩散函数（LSF）→MTF流程评估。结果显示经运动学标定（Kinematic Calibration，KC）校正后，除个别FlyBy轨迹外KC重建MTF较未标定大幅改善；与DLT逐投影校正相比，KC在FB-C及SNG-C上接近DLT水平，在FB-O/FB-S略低——这与KC建立单一全局运动学模型而无法像DLT那样吸收逐视角特异偏差有关，但不影响KC对FlyBy轨迹的适用性优势。
5.7. Observability Results（可观性分析）
通过对残差Jacobian的Gauss–Newton信息矩阵做特征值分解，发现118维待估参数中仅47维位于良好可观子空间，其余弱可观方向主要关联转动轴扰动与部分平移分量，符合所用标定轨迹对特定连杆激励不足的预期。
5.8. Laser Tracker Validation（激光跟踪仪验证）
将图像标定所得运动学模型预测路径与激光跟踪仪实测TCP轨迹对齐后发现仅有常值全局平移偏移（x=0.47 mm, y=0.83 mm, z=1.06 mm），属图像自标定缺乏绝对世界坐标参照的正常现象。源机械臂TCP位置误差标准差由名义模型1.654 mm降至0.462 mm，探测器机械臂由1.205 mm降至0.360 mm，证实标定显著提高了实际轨迹复现精度。
5.9. Qualitative Visual Comparison（定性目视比较）
重建切片显示未标定几何产生明显模糊与细节丢失，而KC与DLT校正视觉效果相当，细微MTF差异需靠定量指标判别。
讨论与结论（翻译浓缩）
所提基于图像的标定方法将重投影误差显著降低，获得满足有界约束且具物理意义的关节与工具修正量，生成的单一运动学参数集可产出内部一致的投影几何，克服了逐投影DLT法须重复位姿故不适用于FlyBy采集的缺陷。重建结果表明运动学标定（KC）在各评估轨迹上均可产生几何一致体积；与DLT相比，KC因建立全局而非逐视角投影模型，在部分含强调制FlyBy轨迹中MTF略低，但这属模型假设差异所致。可观性分析指出标定轨迹仅充分激励部分参数子空间，未来可通过设计更优激励轨迹与扩展修正模型改善。该方法单GPU数分钟内收敛，适合日常再标定。综上，研究人员得出结论：所提无外部仪器的双机械臂X射线CT系统图像基运动学标定框架，通过球基准体重投影误差联合优化可有效修正机械臂及转台几何偏差，在Stop-and-Go与连续FlyBy模式下均能生成全局一致投影几何并显著提升重建质量，为机器人CT系统的原位自标定提供了可行方案。
