一个用于从心脏瓣膜三维图像进行非侵入性患者特异性瓣叶应变估计的几何特征追踪方法

《Annals of Biomedical Engineering》：A Geometric Feature Tracking Approach for Noninvasive Patient-Specific Estimation of Leaflet Strain from 3D Images of Heart Valves

　　

摘要
目的：心脏瓣膜病（valvular heart disease）普遍存在，是心力衰竭的主要诱因之一。瓣叶应变（leaflet strain）是评估瓣膜病理起始与进展背后力学机制的有前景指标。然而，从临床获取的患者图像中非侵入性量化瓣膜应变的通用方法仍然有限。本研究旨在开发一个稳健的特征追踪（feature tracking）框架，能够实现跨多变瓣膜形态的精确形状匹配，并可从三维超声心动图（three-dimensional echocardiographic, 3DE）图像中量化儿童和成人患者的在体房室瓣叶应变。
方法：研究人员开发了一个几何特征追踪框架，用于从3DE图像量化在体瓣叶应变，并评估整个心动周期内的解剖变形。该方法整合了一个基于队列的几何参考图谱（cohort-derived geometric reference atlas）以建立几何对应关系，并在高斯混合模型（Gaussian mixture model, GMM）框架内引入了一种新颖的距离加权一致性点漂移（distance-weighted coherent point drift, DW-CPD）算法用于非刚体配准。研究人员将性能与一个有限元（finite element, FE）基准模型进行了对比，并将该方法与传统的基于点的追踪方法进行了比较。该框架被应用于儿童和成人患者数据集（N = 31），以评估在多变瓣膜形态下的稳健性。
结果：所提出的方法在量化解剖对齐和瓣叶应变方面比传统基于点的方法展现出更高的准确性。与有限元基准的验证确认了应变估计的改进。该框架实现了对儿童和成人患者中不同形态瓣膜变形的可靠相间追踪。分析发现了一个与瓣叶膨出（脱垂）相关的第一主应变（1^st principal strain）的一致分布模式。
结论：这一特征追踪框架提供了一种从临床3DE图像中非侵入性量化房室瓣瓣叶应变的通用方法。生物力学应变模式的特征描述可能改善预后评估，并支持对心脏瓣膜病的纵向评价。对心脏瓣膜病生物力学特征的进一步研究有可能增强对心脏瓣膜病的预后评估和纵向评价。

以下是根据论文主体内容撰写的论文解读文章：

**研究背景与问题**

心脏瓣膜病（valvular heart disease）是一种普遍存在的心脏疾病，显著促进各年龄段心衰的发生与发展。超声心动图是临床评估瓣膜形态异常（如瓣叶膨出、脱垂和帐篷样改变）的金标准。三维超声心动图（3DE）能够可视化动态的三维瓣膜几何结构，并量化反流容积和有效 orifice 面积等功能指标。然而，仅凭全局形态特征不足以刻画驱动疾病起始、进展或手术修复长期耐久性的生物力学环境。相比之下，应变分析能提供反映组织硬度和边界条件的局部组织变形详细图谱，可作为识别退行性变和潜在组织损伤高风险区域的有效指标。目前，从3DE图像（包括经食管超声心动图和经胸超声心动图）中可靠地非侵入性评估瓣叶应变，对于推进瓣膜病的预后评价和提高手术修复耐久性预测具有巨大潜力。现有的特征追踪研究主要集中在二尖瓣（mitral valve, MV）上，且大多算法为专有或局限于商业软件。此外，由于瓣叶在开放与闭合时间帧之间巨大且快速的几何变化，先前研究主要追踪从收缩中期到收缩末期的瓣膜，未能捕捉几乎无应变的舒张期开放状态。目前缺乏一个通用框架来跨所有三维模态分析瓣叶应变。因此，本研究旨在开发一个稳健的开源特征追踪应变量化框架，用于评估包括先天性心脏病患者在内的多种瓣膜形态。

**研究内容与结论**

研究人员开发了一个几何特征追踪框架，用于从3DE图像量化在体瓣叶应变。该方法整合了基于队列的几何参考图谱以建立几何对应关系，并引入了一种新颖的距离加权一致性点漂移（distance-weighted coherent point drift, DW-CPD）算法进行非刚体配准。通过有限元（finite element, FE）基准模型验证了该框架的性能，并应用于三个患者队列（N=31，包括儿童和成人）以评估其稳健性。结果显示，与现有基于点的方法相比，所提出的框架在形状匹配和应变估计方面表现出更高的精度，并且成功应用于不同形态的房室瓣（包括二尖瓣和三尖瓣）。研究人员发现，第一主应变分布的中位数和四分位距（interquartile range, IQR）与瓣叶膨出（脱垂）相关，可作为潜在的生物力学标志。该研究为从临床3DE图像非侵入性量化瓣叶应变提供了通用方法，具有改善瓣膜病预后评估的价值。论文发表在《Annals of Biomedical Engineering》。

**主要关键技术方法**

本研究主要采用以下关键技术方法：1）**距离加权一致性点漂移（DW-CPD）算法**：在经典一致性点漂移（coherent point drift, CPD）框架基础上，引入基于指数距离衰减函数的局部性权重，增强空间邻近点对的对应概率，抑制远距离虚假匹配，从而实现对大幅变形瓣叶的稳健非刚体配准。2）**基于队列的几何参考图谱构建**：在每队列内部随机选择一个舒张期瓣叶作为参考，通过DW-CPD将所有同解剖标签的瓣叶配准至该参考，采用广义普鲁克分析（generalized Procrustes analysis）计算平均形状（参考图谱），为不同患者和心相间建立点对点对应关系。3）**表面应变计算**：通过计算参考与变形构型之间的变形梯度，推导格林-拉格朗日应变张量（Green–Lagrange strain tensor），提取第一主应变；同时通过三角形单元面积变化计算面应变（areal strain）。4）**有限元（FE）基准验证**：利用图像衍生的FE二尖瓣模型（Lee-Sacks超弹性材料，施加10 kPa跨瓣压）生成地面真实应变，评估DW-CPD与传统方法（BCPD、经典CPD）的配准精度（均值对称距离MSD、95% Hausdorff距离HD₉₅）和应变估计误差（平均绝对误差MAE）。样本队列来源：数据集来自费城儿童医院（含左心发育不全综合征HLHS患儿的三尖瓣TV和二尖瓣MV）与宾夕法尼亚大学（成人正常二尖瓣MV）的现存数据库。

**研究结果**

**有限元分析验证**
通过FE基准模型对比，所有特征追踪方法均能大致捕获瓣叶变形。DW-CPD在形状相似性（MSD和HD₉₅）上优于BCPD和经典CPD，并取得最低的面应变MAE（0.138）和第一主应变MAE（0.055）。与FE地面真实值的散点图显示强相关性：面应变Pearson r=0.853，第一主应变Pearson r=0.710（均p<0.001）。误差主要集中于高曲率区域（如连合褶皱附近）。

**节段性瓣叶心房表面评估**
三个队列中各瓣叶舒张期与收缩期表面面积变化幅度差异较大。HLHS三尖瓣（TV）和儿童二尖瓣（MV）组显示正负百分比变化并存，而成人MV组多数瓣叶面积为负变化，提示可能显著的环状收缩。

**特征追踪精度定量评估**
在所有队列中，DW-CPD配准后的瓣叶几何与分割地面真实之间一致性良好，MSD均低于0.4 mm，HD₉₅低于1 mm。前叶在HLHS TV组中误差最大（MSD平均0.357±0.068 mm）。与经典CPD相比，DW-CPD消除了叶缘自穿透等伪影。

**病理三尖瓣应变评估（N=11）**
在11例HLHS患儿中，瓣叶膨出（脱垂）特征表现为较宽且平坦的第一主应变分布，中位数和IQR升高。轻度以上反流组前叶的平均面应变和第一主应变高于轻度以下组，但后叶反之。队列水平主成分分析（PCA）显示，与平均形状相比，最拴系（tethered）和最膨出（billowing）的瓣膜均表现出更高的平均应变和应变变异度（标准差SD）。

**正常与病变儿童二尖瓣应变评估（N=10）**
10例儿童患者中，前后叶应变分布差异显著。病变组（如膨出、风湿性、马凡综合征）中，病例10（后叶膨出）表现出中位数和IQR均>0.5的第一主应变分布。正常组与病变组平均应变差异较小（每叶<0.05）。PCA显示，前后叶间距最小（-2 SD）和最大（+2 SD）的瓣膜均比平均形状具有更高的第一主应变。

**成人二尖瓣应变评估（N=10）**
10例成人患者（均无二尖瓣膨出）显示狭窄的应变分布。前叶和后叶平均面应变均为负值（-0.215和-0.107），后叶第一主应变平均值为前叶的2.8倍。PCA显示，随着前后叶间距增大，平均面应变和第一主应变均递增，且应变SD增加。

**讨论总结**

研究人员开发了一个开源的特征追踪框架，用于量化房室瓣瓣叶应变，并通过FEA基准验证了其准确性和稳健性。在三个儿童和成人队列（共31例）中成功应用，范围涵盖正常瓣膜到严重膨出形态。该框架实现了跨队列一致的几何配准，无需调整参数，展现了通用性。第一主应变分布的中位数和IQR被识别为瓣叶生物力学改变的潜在标志，在膨出瓣叶中显著升高。与现有文献一致，后叶应变普遍高于前叶。但成人MV中的负的面应变可能反映了成像和分割的局限性（如闭合区分辨率不足）或真实的环状收缩效应。未来研究需要更大样本验证应变标志物与病理的关联，并改进图像分辨率和分割精度以减少偏差。研究结论部分翻译如下：

**结论**
研究人员开发了一种新颖的基于点的特征追踪方法，用于追踪舒张期与收缩期帧之间的瓣叶运动。该方法提供了与FEA相当的应变估计。我们成功地将该方法应用于三组房室瓣（共31例儿童和成人患者），涵盖了从健康二尖瓣到伴有瓣叶受限和显著膨出的广泛瓣膜解剖谱。尽管未发现能清晰区分健康与病变瓣膜的统一应变模式，但观察到一个有希望的趋势：应变分布的中位数和IQR升高似乎是瓣叶膨出的强标志。经过进一步验证，该指标可作为一种新的定量工具，增强对心脏瓣膜病的预后判断。