用于超声引导经皮手术的3D即时（JIT）建模：基于主成分分析（PCA）压缩与曲率感知约束插值的实时体积组织变形预测

《Simulation Modelling Practice and Theory》：3D just-in-time modelling of real-time volumetric tissue deformation for ultrasound-guided percutaneous surgery

　　

三维即时（3D just-in-time，3D JIT）建模可实时预测体积组织变形。主成分分析（principal component analysis，PCA）压缩降低密集场预测成本。曲率感知加权提高局部插值鲁棒性。该方法支持变形一致的超声切片渲染。在多个体积网格上实现了毫秒级运行时间。

论文围绕超声引导经皮手术（ultrasound-guided percutaneous surgery，usPS）中实时体积组织变形预测问题，提出三维即时（3D just-in-time，3D-JIT）建模框架，并结合主成分分析（principal component analysis，PCA）输出压缩与曲率感知约束插值。研究背景方面，usPS使用细器械经皮肤小穿刺到达组织目标，常用于活检、近距离放疗、肾取石和肿瘤消融等，超声因便携、无辐射、易获取且对脆弱患者安全而成为主要成像方式。但usPS学习曲线陡峭，需要同时成像工具并操控工具，工具引起的组织与目标变形会移动解剖结构。二维超声（2D ultrasound，2D-US）需要反复扫查并心理组合平面图像；三维超声（3D ultrasound，3D-US）可提供感兴趣区全景视图，但帧率低，且未消除工具操控与成像之间空间协调的学习曲线。模拟训练可提供无风险练习，但赛博物理usPS模拟器需要实时更新工具-组织相互作用和相应图像反馈。传统有限元模型（finite element model，FEM）虽准确，却计算昂贵，大变形非线性问题困难，标准硬件实时受限。因此，有必要将变形生成离线化，并在在线阶段以快速代理预测替代重复模拟。已有非FEM和数据驱动方法如位置型动力学、质量-弹簧模型、无网格方法、机器学习代理等，但存在训练需求、可解释性、内存成本和泛化性折衷。JIT建模可利用与查询最匹配的局部训练样本即时建立局部预测模型。然而，3D体积组织变形中每个数据库样本对应密集位移场，若网格为??×??×??节点，则每样本含3MNP位移分量，存储和在线重组成本随空间分辨率立方增长。为此，研究人员提出3D PCA-JIT框架，将PCA压缩与JIT建模结合，在线插值在降维系数空间进行后再重建全体积位移场。

研究人员开展的核心工作是：构建离线历史变形数据库，在线检索相似针插入姿态，计算曲率感知约束权重，并在PCA降维空间预测后重建全3D位移场。该方法不依赖特定在线求解器，数据库可由任意能返回体积位移场的模拟引擎生成；本文采用3D多谐样条（3D polyharmonic spline，PHS）作为参考变形引擎。结论表明，PCA-JIT在测试网格分辨率下与PHS参考变形保持接近，预测达到毫秒级，速度随体积网格分辨率增加而提高，在42³网格达到89.0×加速；预测位移场还可耦合超声体积变形与切片提取。其意义在于为交互式超声模拟提供实时、可扩展的体积组织变形代理，支持变形一致可视化，并降低对在线FEM式求解的依赖。

关键技术方法概括如下，不超过250字：离线数据库由PHS生成，样本覆盖针入口点、插入方向和插入深度；入口面采用10×10采样，方向采用10×10采样，深度采用4个采样值，得到N_h=40,000个历史样本/网格分辨率。PCA保留秩k=100。在线采用6D针姿态，门控阈值γ_in=0.05、γ_en=0.15，选择N_opt=40个邻居；曲率感知权重通过Hessian描述符进行正则，并求解Karush–Kuhn–Tucker（KKT）约束系统以满足无偏性和一阶一致性。评估使用100个held-out查询，比较PHS参考、FULL-JIT消融和曲率-free消融；超声渲染使用Blue Phantom肾活检超声体模的200×200×87体素体积；离体猪组织实验作为PHS验证背景。

研究结果按原文小标题概括如下。“2. 3D Just-in-Time modelling”：通过离线PCA压缩、在线邻居检索、约束权重、降维预测和全场重建，预测查询针姿态对应的密集体积位移场。“2.1 Offline compression and local curvature”：PCA将位移场压缩到k维系数空间，并估计有效曲率矩阵；Hessian描述符由保留PCA系数对6D针姿态的二阶变化计算，经聚合和高斯平滑后存储。“2.2 Online optimal neighbourhood selection”：通过插入尖端门控和入口门控筛选候选，再按6D姿态欧氏距离选择N_opt=40个邻居，以保证几何相似性。“2.3 Curvature-aware regularized weights”：求解正则化KKT系统，权重满足1^Tw=1和D^Tw=0，曲率项降低对局部非线性较强邻居的影响。“3. Historical database construction”：采用PHS模拟针插入位移，按入口点、方向和深度采样生成N_h=40,000样本；数据库特定于组织响应、交互参数、几何和边界条件。“4. Evaluation criteria”：以PHS为参考，计算ε_MAE、ε₉₅、ε_max和速度S=t_PHS/t_JIT；100个held-out查询属于覆盖姿态空间内的插值测试。“5. Simulation results”：代表性插入中PCA-JIT与PHS接近，误差主要局部于最大变形区；PCA-JIT与FULL-JIT接近，说明PCA保留主导变形模式；误差随网格分辨率保持较小；在线运行毫秒级，PHS为数百毫秒，42³速度达89.0×；存储从约35.6 GB降至约0.11 GB；归一化误差10^?3对应40×40 mm视野约0.04 mm，平均绝对误差约0.053–0.070 mm；曲率感知相对曲率-free在所有分辨率下略降低ε_MAE、ε₉₅、ε_max，运行时间基本不变。“6. Ultrasound volume deformation and slice rendering”：将预测位移插值到超声体素网格，采用拉回式重采样生成变形超声体积，并用三线性插值提取切片；以200×200×87体素肾活检超声体模体积定性展示直线和斜插入，切片显示针轴和轨迹。“7. Limitations”：数据库限于直线针插入，由入口和尖端坐标参数化，未含弯曲针、斜面尖端偏转、针弯曲或动态转向；组织为均质归一化块，未显式包括异质、各向异性、断裂、切割或患者特异解剖和材料；PHS作为参考引擎，PCA-JIT独立于引擎但训练模型反映该引擎假设，未独立生物力学或临床验证；未全面基准所有代理模型；超声变形示例为定性渲染，定量图像验证为未来工作。“8. Conclusion”：研究结论翻译为：本文提出用于超声引导经皮介入中实时预测体积组织变形的3D即时建模框架。所提方法使用参考变形场的离线数据库，在线检索几何相似的针插入，计算曲率感知约束插值权重，并在降维PCA系数空间预测变形场，然后重建全体积位移场。在测试网格分辨率上，PCA-JIT与所选PHS参考变形保持接近一致，同时实现毫秒级预测时间。速度随体积网格分辨率增加，在42³达到89.0×，证明降维空间预测对密集体积输出的优势。预测位移场还与超声体积变形和切片提取耦合，说明该框架支持交互式超声模拟中变形一致可视化。未来工作聚焦更广实验验证、扩展数据库覆盖、组织特异和患者特异参数化、与额外代理基线比较以及集成力反馈和触觉模拟。讨论部分进一步强调，当前框架在直线插入、均质组织、PHS参考和定性超声渲染条件下评估，未来需扩大实验验证、扩展数据库覆盖、组织特异与患者特异参数化、与额外代理基线比较，并集成力反馈与触觉模拟。作者为Kian Zalzalah、Carlos Rossa。