论文技术总结：双心室修复手术中半自动生成的室间隔补片几何结构及其血流动力学评估

1. 研究背景与问题 (Problem)

先天性心脏病（CHD）与双心室修复的挑战：
先天性心脏病是常见的出生缺陷。对于部分复杂病例（如右室双出口，DORV），传统的单心室姑息治疗伴随较高的长期发病率和死亡率。当解剖条件允许时，**双心室修复（Biventricular Repair）**是更优策略，旨在恢复双心室生理功能。

核心难点：
DORV 患者通常伴有室间隔缺损（VSD），且主动脉和肺动脉均主要起源于右心室。修复手术需要构建一个心室内补片（Intraventricular Baffle），将左心室（LV）的血液通过 VSD 引导至主动脉，同时避免血流梗阻。

• 当前临床局限： 补片的几何形状通常在手术中（心脏停跳状态下）由外科医生实时手工塑形。
• 缺乏术前规划： 这种手动方式限制了术前对替代方案的评估，难以生成解剖贴合、水密（watertight）且适合定量血流动力学评估的计算模型。
• 解剖异质性： DORV 患者的解剖结构差异巨大，导致补片设计高度依赖个体解剖，缺乏标准化的计算设计流程。

2. 方法论 (Methodology)

本文提出了一种半自动计算框架，用于生成和评估患者特异性的心室内补片。该流程包含三个主要步骤：

2.1 患者特异性 3D 模型重建

• 数据源： 使用术前收缩期心脏 CT 血管造影（CTA）图像。
• 分割： 采用基于 nnU-Net 的深度学习框架进行半自动分割，将心脏结构分为 7 类（左/右心室、左/右心房、心肌、主动脉、肺动脉）。
• 人工修正： 由专家对分割结果进行手动修正，特别是 VSD 边缘和主动脉瓣环区域。
• 域构建： 将分割结果合并为连续的血流域，通过截断和封顶（Capping）操作，生成包含 VSD 入口和主动脉出口的水密计算域（$M_{cappedRV}$）。

2.2 补片几何结构生成 (核心算法)

这是该框架的创新核心，旨在生成符合生理要求的补片表面：

1. 边界定义： 外科医生在 3D Slicer 软件中交互式选择缝合线点（Suture line），系统将其转换为密集的样条曲线作为补片边界。
2. 中心线构建 (Centerline)：
   • 摒弃了标准的 VMTK 中心线（适用于管状血管），构建了自定义中心线。
   • 使用 PCHIP (分段三次 Hermite 插值多项式) 连接 VSD 入口和主动脉出口，确保切线方向与入口/出口法向对齐，避免过冲，保证生理路径的平滑性。
3. 截面面积约束与形状调整：
   • 目标面积 ($A_{target}$)： 基于患者体表面积（BSA）和主动脉瓣环直径的 Z 评分，计算符合生理规范的目标横截面积。
   • 截面重塑： 沿中心线提取截面，固定右心室壁部分，仅调整补片部分。通过抛物线位移模型，求解满足目标面积约束的最小变形幅度。
   • 凸性平滑 (Convexity Smoothing)： 为避免局部截面调整导致的几何凹陷和曲率突变（这不符合心包补片的力学特性），对变形幅度序列应用**上凸包络（Upper Convex Envelope）**平滑处理，生成全局平滑的隧道轮廓。
4. 表面生成： 利用 薄板样条 (Thin-Plate Spline, TPS) 插值技术，将边界曲线和平滑后的截面轮廓融合，生成 $C^2$ 连续、最小弯曲能量的最终补片表面。

2.3 血流动力学评估 (CFD)

• 模拟设置： 将生成的补片整合到计算域中，使用 SimVascular 进行稳态计算流体动力学（CFD）模拟。
• 边界条件： 基于术前超声心动图数据设定入口流量，出口设定阻力边界条件。
• 评估指标： 计算跨补片路径的压力梯度（$\Delta P$），并与术后临床测量值进行对比。

3. 关键贡献 (Key Contributions)

1. 半自动计算流程： 首次提出了一套从医学影像到 CFD 就绪模型的完整半自动管道，解决了传统手动建模效率低、难以量化评估的问题。
2. 生理约束的几何设计：
   • 引入了基于 BSA 的目标横截面积约束，确保重建的流出道符合生理规范。
   • 提出了凸性平滑算法，解决了截面独立优化导致的几何不连续问题，生成的补片形状更符合手术中使用的自体或牛心包补片的力学行为（张力主导的薄膜）。
3. 鲁棒的中心线参数化： 针对心室非管状、入口出口大小差异大的特点，设计了基于 PCHIP 的自定义中心线构建方法，优于传统血管中心线提取算法。
4. 水密与拓扑一致性： 确保了生成的几何模型是水密的且拓扑一致，可直接用于高精度的 CFD 模拟。

4. 研究结果 (Results)

研究团队回顾性应用该框架于 4 例 曾接受双心室修复的 DORV 患者：

• 几何生成： 成功为所有 4 例患者生成了光滑、无自交、与原生解剖无缝集成的补片模型。
• 血流动力学表现：
  • 模拟显示血流从 VSD 到主动脉的传输连贯，无大规模回流。
  • 压力梯度： 计算得出的跨补片压力梯度范围为 2 - 5 mmHg。
• 临床对比：
  • 将模拟预测的压力梯度 ($\Delta P_{sim}$) 与术后超声测量的实际梯度 ($\Delta P_{echo}$) 进行对比。
  • 结果显示，所有 4 例患者的模拟梯度均显著低于实际术后测量的梯度（差异在 3.04 至 9.88 mmHg 之间）。
  • 例如，患者 1 的模拟梯度为 2.79 mmHg，而术后超声为 12.67 mmHg。
• 结论： 该计算设计方法能够生成比实际手术中手工塑形更优（阻力更小）的几何结构，且结果在临床可接受范围内。

5. 意义与展望 (Significance & Future Work)

• 术前规划工具： 该框架为复杂先天性心脏病的双心室修复提供了定量的术前决策支持工具，允许外科医生在手术前评估和优化补片几何形状。
• 个性化医疗： 能够处理高度异质性的患者解剖结构，为“边缘性”单心室患者转化为双心室生理提供数据驱动的评估平台。
• 局限性： 当前模型基于刚性壁假设，未考虑心脏收缩运动（FSI）和补片材料的力学特性（如缝合张力、材料刚度）。
• 未来方向：
  • 引入流固耦合（FSI）和电 - 机械 - 流体多物理场模拟，以捕捉心脏运动对血流的影响。
  • 开发全自动边界选择算法，减少对外科医生手动输入的依赖。
  • 结合补片材料力学模型，更准确地预测术后植入形态。

总结： 这项工作通过结合生理约束的几何设计与 CFD 评估，迈出了实现复杂先天性心脏病（特别是 DORV）计算辅助手术规划的重要一步，有望通过优化流出道重建来改善患者的长期血流动力学预后。