方法文章

基于相位分辨功能肺部磁共振成像的肺通气与灌注(V/Q)评估

3.4K 次观看

DOI:

10.3791/66380

2024年8月9日

本文内容

摘要

本文介绍了一种无需对比剂的质子磁共振技术——相位分辨功能性肺部MRI,用于评估肺通气和灌注动力学。该技术已在不同磁场强度和年龄组中得到验证并具有适用性,未来可能通过辅助疾病定量分析和治疗监测,提升临床决策水平。

摘要

傅里叶分解是一种无需对比剂的¹H MRI 方法,用于评估肺灌注(Q)和通气(V)。在图像配准后,对每个体素的时间序列进行分析,以提取心脏和呼吸频率成分。

采用标准的二维扰相梯度回波序列,时间分辨率为约300 ms,开发了一种图像排序算法,以实现具有更高时间分辨率的相位解析功能肺成像(PREFUL)。因此,可在潮气量呼吸期间评估区域性流量-容积环(FVL),并描绘心动周期中脉搏波的传播过程。该方法可在1.5T或3T磁共振系统上,使用标准MR硬件实施,无需进行序列编程,因为所述方案可在大多数系统上通过默认的SPGRE序列实现。

PREFUL 通气 MRI 已通过129Xe 和19F 气体成像技术进行了验证,区域一致性良好。PREFUL 灌注加权 MRI 已通过单光子发射计算机断层扫描(SPECT)以及动态对比增强(DCE)MRI 进行了验证。PREFUL 已在双中心、双厂商环境下进行了测试,目前正应用于多项正在进行的多中心试验中。此外,该技术在多种磁场强度(0.55T–3T)和不同年龄组(包括新生儿)中均具有可行性。

定量通气/灌注(V/Q)PREFUL MRI 已被用于囊性纤维化、慢性阻塞性肺疾病、慢性血栓栓塞性肺动脉高压以及冠状病毒病-2019患者,以量化疾病程度并监测治疗后的变化。此外,PREFUL V/Q 成像已被证明可预测肺移植术后因慢性肺移植物功能障碍导致的移植功能丧失。综上所述,PREFUL MRI 是一种经过验证的技术,可用于区域性肺部疾病的定量通气及肺脉搏波/灌注成像,在疾病检测、量化评估和治疗监测方面具有潜在的临床应用价值,可为当前临床常规提供额外信息。

引言

呼吸系统因其复杂的机制而容易受到多种疾病的侵袭。其中,慢性阻塞性肺疾病(COPD)、囊性纤维化(CF)和慢性血栓栓塞性肺动脉高压(CTEPH)等慢性呼吸系统疾病显著降低了患者的预期寿命1。因此,早期诊断、病情监测以及治疗反应评估变得至关重要。

肺功能测试(PFT)可获得诸如蒂夫诺-皮内利指数(Tiffeneau-Pinelli index)等整体肺功能参数,该指数定义为1秒用力呼气容积(FEV1)与用力肺活量(FVC)的比值2。此类参数在临床常规中已得到广泛应用,但缺乏区域性信息,且对患者的配合度要求较高。在此方面,影像学技术可提供额外的洞察力,并有望实现更为敏感的参数评估。计算机断层扫描(CT)能够对肺实质形态进行高分辨率成像,而近年来发展的参数响应映射等技术还可获取功能信息3。然而,单光子发射计算机断层扫描(SPECT)仍是目前显示肺通气与灌注(V/Q)的金标准4。上述所有影像学方法的共同局限在于需要暴露于电离辐射,在监测性应用及脆弱人群中的使用需特别审慎。因此,目前正持续推动磁共振成像(MRI)作为替代性影像技术的发展。

由于肺部固有的低质子密度和快速信号衰减,肺部本身就是磁共振成像(MRI)中具有挑战性的器官5。在众多成像方法中,目前应用最广泛的技术包括使用超极化气体(例如129Xe MRI)评估通气功能6,以及静脉注射基于钆的对比剂以显示灌注情况7。这些方法具有较高的信噪比(SNR),被磁共振领域广泛视为金标准。一种较新的方法无需使用任何对比剂,可在自由呼吸状态下利用常规质子磁共振成像实现,每层图像的总采集时间约为1分钟。因此,该方法避免了对比剂可能引起的不良事件以及近年来备受争议的长期影响,同时无需额外的超极化设备或多核硬件,更易于推广应用。此外,由于采用自由呼吸采集方式,避免了寻找合适肺部充气状态这一难题,而该问题可能会影响所获得的通气缺陷值8

这种基于间接MR信号的方法最初由Zapke等人提出,他们利用了质子加权信号S与肺容积V之间的倒数关系:S~1/V。9 该方法基于将自由呼吸状态下采集的图像转换至一个共同的膨胀状态(通常位于呼气末与吸气末之间的中间位置),从而补偿运动影响,并允许对每个体素内的信号时间序列进行分析。随后,可通过使用Klimeš等人提出的公式(110,从这些所谓的配准图像中推导出通气量测量结果:

肺通气量方程;呼吸容积、比值及调节分析;公式展示。     (1)

在吸气(Insp)、呼气(Exp)和注册状态(Reg)下的容积/信号基础上,该方法进一步引入了傅里叶分解,以区分与呼吸频率(通气)和脉搏频率(灌注)相关的信号调制,从而仅通过一次采集即可生成完全空间匹配的通气/血流(V/Q)图11。由于呼吸频率与心率之间通常存在明显差异,原本在时域中重叠的两个信号成分可通过傅里叶分析在频域中有效分离,从而实现上述目标。随着设备从低场强(0.35 T)过渡到1.5 T,并采用优化的平衡稳态自由进动序列(bSSFP)12,该方法逐渐受到更多关注,并在多项后续研究中得到验证13,14,15

由于呼吸和脉搏存在变异性,且在1.5T下商业可用的平衡稳态自由进动(bSSFP,梯度补偿)成像可能产生显著的带状伪影(明显的信号缺失条带),因此提出了一种相关方法,即采用扰相梯度回波序列(SPGRE)结合宽频低通和高通滤波技术16,17这能够捕捉到与真实呼吸和脉搏相关的更复杂的信号调制谱。随后在时域中计算振幅,避免了必须选择某一特定频率峰值的需要。通过将传统的一次性向一个参考状态配准分解为两个独立步骤,进一步实现了优化。该方法利用了自由呼吸过程中,在吸气末与呼气末之间采集到的不同呼吸相位,其向固定状态所需形变程度各不相同这一事实。在选择若干组图像并确定各单幅图像所属的呼吸状态组后,执行以下流程:1)在各自呼吸状态组内进行配准;2)逐组进行组间配准,从一个相邻组依次配准至下一组(例如,1->2, 2->3,…) 分配给代表参考组的组别。该方法进一步通过为每幅图像进行相位估计,以建立更高的表观时间分辨率,从而促进对通气和灌注动态的分析,由此形成了相位解析功能肺(PREFUL)磁共振术语,以将此技术分支与其他相关技术区分开来18后续研究利用了完整的呼吸和心动周期所提供的附加信息,显示出这些参数的潜在敏感性可能有所提高19,20,21.

与金标准SPECT进行验证,结果显示缺陷区域的Dice系数≥67%22,且使用129Xe进行更直接的通气测量显示,在混合的慢性阻塞性肺病/囊性纤维化/健康人群中,通气缺陷百分比的相关性≥62%23 以及在CF多中心、多厂商队列中的84%24,该研究还显示PREFUL与129Xe的肺清除指数具有相似的相关性(r = 0.82 和 r = 0.91)。同一研究的灌注分析表明,在所评估的各中心之间,与DCE在空间重叠性方面无显著差异。25. 一项纳入九个中心的前瞻性子研究还报告了PREFUL结果与DCE的一致性以及各中心间PREFUL结果的符合率26在慢性阻塞性肺疾病患者中的可重复性分析结果显示,所有参数的变异系数均低于15%27目前的研究表明,与仅考虑吸气末和呼气末阶段的“静态”通气参数相比,区域流量-容积环(FVL)参数在检测治疗变化方面具有更高的预测能力和敏感性。在慢性阻塞性肺疾病(COPD)患者中,使用茚达特罗-格隆溴铵(indacaterol-glycopyrronium, IND/GLY)吸入治疗后,通过区域流量-容积环(FVL)测量已证实其对治疗具有响应性。28与之相应,FVL 参数可预测双肺移植患者的移植物丢失,而肺量计检查则不能(P = .02 对比 P = 0.33)。29初步可行性研究表明,使用标准临床 MRI 设备,可在自由呼吸的婴儿和新生儿中实现基于 PREFUL 的功能性肺部成像30,31Glandorf 等人比较了在 1.5T 和 3T 场强下(SPGRE 序列)PREFUL 参数的表现,发现大多数参数无显著差异,且尽管磁场强度不同,这些参数仍具有高度可重复性32这一点可能具有重要意义,因为并非所有机构都能获得1.5T或更低磁场强度的扫描仪。最近,通过使用PREFUL方法分析0.55T下的bSSFP数据,已证明了在感染COVID-19后检测持续症状的可行性和有效性。33.

总之,尽管PREFUL仍是一种相对较新的技术,但已得到广泛研究。研究人员已评估了若干重要标准,包括与更直接且成熟测量方法的验证对比、可重复性、对病理变化的敏感性,以及对治疗和病情进展变化的响应能力。然而,尽管该技术对设备要求较低,目前仍仅有少数专业中心在使用。因此,本研究的目的是以文字和可视化形式总结PREFUL MR的最新操作方法。这些信息可用于在更多机构中推广该技术,从而从长远来看推动该技术走向成熟。

方案

本研究获得了汉诺威医学院伦理委员会的批准,确保在整个研究过程中严格遵守伦理标准。该研究严格遵循《赫尔辛基宣言》中规定的指导原则,强调医学研究的伦理规范。此外,在所有受试者(或其父母或法定监护人)参与 MRI 检查前,均已认真获得知情同意。有关核心实验流程的简化概述(包括数据采集、配准、滤波和排序,以及最终的心脏和呼吸周期合成),请参见图1。以下各节将详细描述所有相关步骤。

心脏成像过程:层面采集、配准、滤波、周期合成;各阶段示意图。
图1:PREFUL核心组件的示意图。 1) 自由呼吸状态下的数据采集,这需要进行2) 基于配准的运动补偿(如轮廓图所示),从而实现对傅里叶分量进行逐体素分析,如步骤3) 滤波与排序所示。经过低通(通气)和高通(灌注)滤波后,4) 利用估计出的相位将图像按更高的表观时间分辨率进行排序,并合成一个完整的心动周期和呼吸周期。请注意,此图为简化示意图,完整细节见原文。本图中省略了其他步骤,包括参数量化和报告生成。缩写:PREFUL = 相位分辨功能肺成像。 请点击此处查看该图的放大版本。

1. 招募

  1. 招募 成人患者或未成年人基于肺功能检查和/或既往影像学检查(如X光或CT扫描)诊断为肺部疾病,纳入PREFUL MRI试验的患者及健康对照需满足以下条件:能够提供知情同意;能够在MRI扫描期间平躺并保持静止;无MRI检查已知禁忌症(如金属植入物、幽闭恐惧症)。
  2. 排除患者 如果受试者处于妊娠状态、在过去6个月内接受过肺部手术、存在严重呼吸窘迫或需要持续吸氧,或既往对钆基对比剂有不良反应(若除PREFUL外还需进行增强MRI检查),则排除。
  3. 招募 健康对照组 如果没有 n无已知的肺部疾病史,依据临床检查和自我报告确认;肺功能测试结果正常;能够提供知情同意;能够平躺并在整个 MRI 扫描过程中保持静止;且无 已知的 MRI 禁忌证。
  4. 排除 个体作为 健康对照组 根据当前或既往吸烟史、已知肺部毒素或职业危害暴露史、遗传性肺部疾病的家族史、任何可能影响肺功能的已知慢性疾病,以及妊娠情况。
  5. 获取 知情同意书。若已从参与者处获得包含研究目的和程序、潜在风险与益处、保密保证、研究持续时间以及 有权随时退出且不产生任何后果。
    1. 当招募对象包括未成年人时,需应用以下附加步骤。
      1. 除未成年人的同意外,还需至少获得其一位家长或监护人的书面知情同意。
      2. 以适合未成年人理解的方式呈现整个过程。
      3. 对于年龄较小的未成年人,应使用视觉辅助工具、故事书或简化的解释方法。
      4. 确保研究与目标年龄组相关,并避免不必要地纳入未成年人。
      5. 允许未成年人在所有与试验相关的操作过程中由一名可信赖的成年人(如父母、监护人)陪同,除非此举会影响试验的完整性。
        注意:上述步骤假设该试验已通过伦理委员会的审查并获得批准。

2. 数据采集

  1. 执行预扫描安全检查程序。
    1. 在进行MRI之前,应对患者进行详细的预筛查,以识别潜在的禁忌症,如既往手术史、植入式器械、纹身或金属碎片接触史。
    2. 到达后,向患者讲解该操作的磁性特性及相关风险。
    3. 指导患者取下所有个人金属物品,包括珠宝、手表和某些衣物部件,必要时为其提供检查服。
    4. 由受过培训的专业人员对患者进行目视检查,以发现被遗漏的金属物体。
    5. 检查所有医疗器械或植入物是否与 MRI 兼容。
    6. 确保制定应急方案,以保障工作人员和患者的安全。
  2. 摆放受试者并进行扫描
    1. 将患者头先进入,仰卧位置于0.55T、1.5T或3T系统上。
    2. 提供听力保护装置、紧急呼叫铃、防护垫和毯子,以确保安全与舒适。
    3. 将多通道柔性线圈置于下颌下方,以确保对所有肺部区域实现最佳线圈灵敏度。
    4. 固定线圈位置以保持稳定,同时不妨碍患者的呼吸。
    5. 指导患者闭上眼睛,然后使用磁共振激光标记肺部中心位置。
    6. 将肺部中心对准等中心,并让患者再次睁开眼睛。
    7. 扫描初始定位像以确定大致方位,随后进行横断面形态学扫描,以识别气管分叉。
    8. 以气管分叉处作为一致的解剖标志,固定第一张冠状切片,以提高可重复性。
    9. 根据扫描方案,选择捕获 三片 间隔一个 层间距 (从边缘到边缘测量)的 100% (切片厚度)或采集 多个切片 贯穿整个肺部,距离为 20%33%.
    10. 完整获取每个切片 分别地 且不交错。
    11. 将重建后的图像上传至影像归档与通信系统(PACS),以便后续访问与分析。或者,为遵守隐私与数据保护标准,可手动将图像导出至指定的网络驱动器或类似的存储解决方案。
      注意:有关序列方案和参数的详细信息,请参阅 表1 表2有关切片定位的可视化示意图,参见 图2.

CT扫描切片定位;正确与不充分缩放的直方图;信号分析示意图。
图2横断面三维梯度回波图像中显示的 PREFUL 实验典型切片定位。 请注意,第一层切片以气管分叉处为可重复的解剖标志进行定位。2nd 和 3rd 切片在前、后方向上以 100% 的切片间距放置。代表性直方图显示了适当与不适当(过低)缩放条件下的数值分布情况。后者会导致动态范围降低和精度损失。同样应避免过高的不适当缩放,因其会引起信号截断(此处未显示)。缩写:PREFUL = 相位解析功能性肺成像。 请点击此处以查看此图的放大版本。

表1:PREFUL方案的典型流程。 在完成定位像扫描后,以横断面方向获取肺部的三维容积图像。该图像用于规划后续的PREFUL扫描(序列参数设置参见表2,层面定位参见图2)。根据研究需要,可增加其他序列。缩写:PREFUL = 相位解析功能性肺成像。 请点击此处下载该表格。

表2:采用扰相梯度回波序列进行PREFUL采集的序列参数汇总。 缩写:PREFUL = 相位解析功能性肺成像;SPGRE = 扰相梯度回波序列 请点击此处下载该表格。

3. 后处理

排序指标图表;数据分析;时间范围;组内配准过程;统计图。
图3:组导向配准的示意图,用于最小化运动补偿所需的形变。 将图像(以圆圈表示)根据排序指标(例如分割的肺区面积)分为10组后,对每组内的图像进行组内配准,使其对齐至一个中间位置(以第1组的虚线表示)。随后,将配准后的图像进行平均,并逐步向中间组进行组间配准的最终步骤。缩写:GOREG = 组导向配准。 请点击此处查看此图的放大版本。

灌注通气图结合信号分析,显示采集数据与周期插值。
图 4:灌注与通气排序算法示意图。 对于灌注(左侧),通过分段拟合(上排)估算相位,并对采集数据重新排序(下排)。对于通气(右侧),首先排除异常值(上排),然后根据基于振幅及其微分的余弦模型(下排)进行排序,以区分呼气与吸气阶段。 请点击此处查看该图的放大版本。

肺通气分析;压力-容积环图;区域性流量-容积评分彩色图谱;呼吸研究。
图 5:一名43岁女性慢性阻塞性肺疾病(COPD)患者的典型流量-容积环及其对应的FVL-相关性度量指标。 请注意,随着流量-容积环(FVL)的变化,FVL-CM值随之降低。缩写:FVL = 流量-容积环;FVL-CM = FVL-相关性度量指标;COPD = 慢性阻塞性肺疾病。 请点击此处查看该图的放大版本。

表3:参数图的阈值应用及缺陷图的合并。 请点击此处下载该表格。

  1. 注册
    注意:以下部分描述了采用基于组的定向(GOREG)方案,将图像非刚性配准至参考体积(优选中等呼吸水平)以补偿呼吸和心脏运动的步骤。
    1. 获取图像。
    2. 首先,对所有未配准的图像使用训练好的U-Net进行肺部分割,并应用截止频率为0.7 Hz的低通滤波器,以估算呼吸时相。
      注意:这将启动全自动处理流程,后续步骤将在后台自动执行。所演示和描述的处理步骤不限于特定应用程序或编程语言,因此可通过多种编程语言在自定义应用程序中实现。
    3. 将大量分割出的体素归类为吸气相,少量分割出的体素归类为呼气相。
    4. 将图像按十分位数进行分组,确保生成的10个组中每组包含相同数量的图像。
    5. 选择 ANTs34 (使用互相关度量的B样条SyN)Forsberg35,36 (弹性与流体正则化的多项式展开)作为 注册算法.
    6. 对每组分别进行组内配准,以该组各自的中间肺部位置为目标。
    7. 对各组结果进行平均处理,以获得每组的一张图像。
    8. 将每组图像沿5方向向下一个相邻图像进行配准th
    9. 根据需要,将形变链应用于原始图像,以达到第5组平均图像所代表的呼吸相位置。例如,对于属于第3组的图像36,需执行以下步骤: 形变场 应用: 36->图像 57(第3组中等肺部位置)-> 步骤 3->4 -> 步骤 4->5.
      注意:有关 GOREG 程序的详细描述,请参见 图3使用并行计算进行配准以减少处理时间。采用GOREG配准方法,旨在最小化每一步配准所需的形变量,从而确保算法的稳定收敛。然而,也可仅通过一次形变步骤直接向中间肺部位置进行配准。尽管所有已发表的2D PREFUL研究均使用ANTs,Forsberg方法的处理速度可提升高达6倍,且结果质量与一项3D PREFUL研究中报告的结果相当。37本报告中的示例结果均采用 Forsberg 配准生成。请参见 图1 以观察膈肌运动的配准效果。
  2. 通用过滤
    1. 使用图像引导滤波对配准后的图像进行去噪38,以时间平均配准图像作为引导图像。应用以下设置: 邻域大小 = [10, 10],平滑程度 = 1.
    2. 用于通气和灌注分析时,采用具有低通或高通滤波的滤波器 截断值 0.7 Hz 以抑制各自的另一组分。除定量灌注计算外,所有后续处理步骤均排除前20幅图像,以确保所包含的时间序列达到稳态。
      注意:如果受试者的呼吸频率高于约40次/分钟,可能需要调整截断值。
  3. 分割
    注意:最终分割使用中间肺部位置的配准图像,通过以下所述的两步法进行。
    1. 对时间平均配准图像使用训练好的U-Net或手动进行肺边界(肺ROI)分割。
    2. 然后,排除较大的中央血管,以优化肺边界分割,并获取肺实质的感兴趣区域(ROI)。
  4. 灌注
    注意:以下步骤用于精确估算所获取序列中每张图像的心动周期时相,该序列的采样频率相对较低(约3-5帧/秒),尤其与心率(通常为40-90次/分钟)相比。通过将数据按照确定的心动周期时相进行重排,可获得完整的心动周期,并实现超越数据采集采样率的时间分辨率18 (参见 图 4 以图示说明分选流程)。进行相位估计时,需要一个具有强灌注加权信号的感兴趣区域(ROI)。采用如下迭代搜索算法:22.
    1. 将肺边界ROI连接起来,以在搜索ROI中包含纵隔。
    2. 通过计算图像序列的标准差生成一幅简单的灌注加权图。
    3. 在搜索感兴趣区域(ROI)内,确定该图谱第98百分位数对应的区域,作为后续步骤的种子感兴趣区域(seed ROI)。
    4. 执行分段拟合,逐步增大种子点的尺寸,直至拟合效果不再提升。
    5. 根据拟合性能对扩展的种子ROI进行排序。
    6. 迭代地将最优扩展的种子感兴趣区域(ROI)与次优、第三优等依次组合,直至组合不再提升评估指标或所有种子ROI均已考虑。将最终得到的ROI视为用于心脏周期估计的血管ROI。
    7. 在优化的相位估计区域内对信号进行空间平均,以生成一个用于相位估计的信号-时间序列。
    8. 通过信号的局部极大值将信号分割成较小的部分,进行分段估计,随后采用分段正弦拟合,考虑振幅、相位偏移和频率等参数(图4).
    9. 对图像进行相位排序,以表示一个心动周期。
    10. 采用 Nadaraya-Watson 核回归方法,结合高斯核(sigma = 0.1),将 15 个相位插值到覆盖单个心动周期的均匀时间网格上。
      注意:请参阅 图1 从舒张期开始,经收缩期,再回到舒张期的健康志愿者部分代表性合成完整心动周期。
  5. 通气
    注意:在灌注分析中,需注意心率保持相对稳定,且振幅变化可忽略不计。相比之下,通气过程往往在潮气量和频率上表现出更多变异,导致不同的呼吸状态,相同呼吸时相的信号振幅并不总是相同。受Fischer等人自门控方法的启发,根据信号振幅对通气进行分类至关重要。
    1. 使用经验法则排除极端异常值(低于第5百分位或高于第97百分位的数据)。
    2. 根据为注册分组创建的信号-时间序列推导出振幅范围 R 和偏移量 C。
    3. 定义一个具有任意选定频率 f 的模型函数 A(t)呼吸作用 (此处为 0.3 Hz):
      呼吸周期方程 A(t) = Rcos(2πfRespirationt) + C;用于分析的数学表达式。
    4. 根据斜率将数据分类为吸气和呼气状态。
    5. 根据模型函数实现更精确的相位测定:
      呼吸周期时相方程,t_吸气 × f_呼吸 = cos⁻¹((A(t)-C)/R) / (2π)。
      呼吸周期方程;公式:t<sub>Expiration</sub> * f<sub>Respiration</sub> 与反余弦函数。
    6. 随后,根据相位对样本进行对齐,并应用Nadaraya-Watson核回归,在呼吸周期内均匀间隔的时间点上计算通气量。
    7. 参照公式1,将吸气相替换为相应相位,计算每个相位的局部通气量(RVent)。
      静力平衡方程,RVent(相位),数学关系示意图,用于教学。
      注意:请参阅 图1 查看一名健康志愿者完整呼吸周期的示例合成子集,该周期始于呼气,过渡至吸气,再返回呼气。有关排序算法的图示,请参见 图4.
  6. 参数计算
    注意:利用合成的完整呼吸和心动周期,可以进一步推导出更多参数。以下将介绍其中一些最重要的参数。
  7. 区域通气(RVent)
    1. 根据吸气相计算 RVent:
      呼吸通气方程显示通气量(RVent)的计算:R_Vent = S_Mid/S_Insp - S_Mid/S_Exp。
  8. 流量-容积环相关性度量(FVL-CM)
    注意:为评估所有呼吸时相,需执行一系列步骤,以生成相当于用力-容积环(FVL)分析的 MRI 图像,类似于肺功能检测。
    1. 使用区域通气量(RVent)的一阶时间导数计算其斜率,作为流量的替代指标。采用步长为 h 的对称差商法:
      数值微分公式,df/dt≈(f(t+h)−f(t−h))/2h,数学概念。
    2. 可选择性地将区域或平均的 RVent 斜率表示为相应 RVent 的函数,从而生成与 FVL 分析等效的 PREFUL 结果。
    3. 通过在肺实质ROI内识别相对通气(RVent)值位于第80至第90百分位数范围内的最大连通区域,确定一个参考ROI。
    4. 对确定的参考ROI内的流量-容积环路进行平均。
    5. 为确定每个肺部FVC与参考信号的相似性,将肺实质ROI内的每个FVC与参考信号进行零滞后互相关分析:
      互相关公式 Rxy(0)=Σxnyn,信号处理,数学方程
    6. 根据以下要求进行标准化:
      静力平衡方程;Rxy,norm(0) = 1/√[Rxx(0)Ryy(0)] Rxy(0);数学表达式
      此处,x 和 y 分别代表参考曲线和相应的 RVent 流量曲线。
      注意:采用零滞后,因此延迟通气会导致相关性降低。
      注意:请参阅 图5 有关FVL计算及其衍生的FVL-CM指标的图示。
  9. 定量灌注
    注意:定量分析根据 Glandorf 等人的方法,使用瞬态过程中获取的第一幅图像进行36.
    1. 使用步骤 3.1.2 中计算得到的肺部体素数量 A,将前四幅已配准的图像归一化至中期吸气水平。这种对由质子密度变化引起的调制效应的校正可表示为:
      信号转换方程;公式:S_Phase→mid=A(phase)/A(mid)×S_Phase;科学分析
    2. 使用以下模型对信号进行指数拟合,以估算与最大磁化相关的信号:
      物理学或化学分析中的指数衰减公式 \(M=S_0 \exp(-bt)\)。
    3. 使用肺实质ROI中显示最强信号的心动周期时相,确定与实质灌注相关的图Q。
    4. 为估算区域血流分数(BF),将 S 进行归一化处理0 通过计算搜索感兴趣区域(全血体素)中高于第99.99百分位数值的平均值来确定该值:
      贝叶斯因子公式 BF=S₀/S₉₉,₉₉% 用于统计分析和概率计算。
    5. 通过考虑稳态(SS)与初始状态之间最大中位信号差异的比值来估算心动周期中的交换分数(EF),该比值由 S 确定0 以及与血流相关的信号差异 Q:
      效应量计算公式:EF = Q/Median(S₀−SS);统计分析概念
    6. 测定心率 f心脏 使用傅里叶分析(对应最大峰值的频率)从血管ROI中提取1/s单位的数据。
    7. 计算最终定量灌注(QQ),单位为 mL∙min-1∙100 mL-1 如下所示:
      心输出量方程;使用射血分数(EF)、血流(BF)、容积(VV)和心率的公式;数学平衡分析。
      注意:此处体素体积(VV)被约去 [mL/mL],并使用换算因子 60 s/min 和约定因子 100/100,以将最终结果表示为 [mL∙min⁻¹∙100 g⁻¹]-1∙100 mL-1]。EF 和 BF 是无量纲比值。
  10. 阈值分割与统计分析
    1. 对上述参数进行统计学描述,包括中层切片及所有切片的平均值(所有值/中层切片值)和标准差(所有值/仅中层切片值)。
    2. 此外,将标准差归一化为变异系数,以获得离散程度的相对度量。选择均值和变异系数作为最终的统计输出结果。
    3. 对参数图应用阈值以生成缺陷图,并计算缺陷百分比值(参见 表3).
    4. 将低于阈值的区域分类为通气或灌注缺损(VD/QD)。进一步合并这些图谱以量化缺损区与正常区的重叠情况(V/Q 分类),包括如四格表中所示的以下组合: 表3:
    5. 计算通气缺陷(VD)、灌注缺陷(QD)和通气/灌注(V/Q)类别的缺陷百分比,即相应类别的体素数量占肺实质总体素数量的比例:
      缺陷百分比 = 缺陷体素数 / 肺实质体素数
    6. 为每张切片及复合冠状切片计算该缺陷百分比。在本研究中,选择联合方法,其中通气缺陷由“或”运算确定:VD = VD(RVent) OR VD(FVL-CM)。
      注意:上述分析使用一款商业软件应用程序完成(参见 材料清单) 使用 Forsberg 配准工具箱。
参数阈值说明
RVent90th 百分位 × 0.4自适应阈值
FVL-CM90%固定阈值
Q90th 百分位 × 0.15自适应阈值

V/Q无 QDQD
无 VD正常VQ 失配(仅 QD)
VDVQ 失配(仅 VD)VQ 缺陷-匹配

结果

图2 的下半部分展示了适当和不适当缩放对动态范围所产生的相应影响。图6 显示了信号分布的不均匀性,分别代表了未进行线圈归一化和已进行线圈归一化的扫描结果。建议避免低动态范围以及未进行线圈归一化的图像。

MRI线圈校正效果;校正前后对比;胸部扫描比较;成像结果分析。
图6:未进行线圈校正(错误)与进行线圈校正(正确)后采集的代表性图像。 注意靠近线圈元件的体表边界处存在人为的信号增强现象。 请点击此处查看该图的放大版本。

图7展示了成功与失败的自动分割结果。请注意,失败的分割未能包含全部肺部体素,这将导致后续分析和统计结果失真。对于存在渗出病变的病例需格外注意,因为此类体素可能因信号强度较高而被人工智能模型误判为血管,或完全未被分割出来。

CT肺部分割过程示意图,比较受试者A和B的切片分割阶段。
图7:自动化分割结果示例,显示失败结果(受试者A:83岁男性,慢性阻塞性肺疾病)与成功结果(受试者B:30岁女性健康对照)。 第一行为输入AI模型的图像。第二行为第一阶段分割结果,即肺边界识别结果。第三行为排除血管后的最终结果。如蓝色箭头所示,算法因高信号肺部变异而受到挑战,导致肺边界检测错误。需注意,图像已通过最大信号进行归一化处理,但由于受试者A的扫描未进行线圈归一化,导致结果存在差异。红色区域显示自动检测出的感兴趣区域(ROI),用于灌注期排序。缩写:COPD = 慢性阻塞性肺疾病;ROIs = 感兴趣区域。 请点击此处查看该图的放大版本。

图8图9 展示了一名健康对照者(30岁,女性)和一名慢性阻塞性肺疾病(COPD)患者(60岁,男性)的代表性参数图。请注意,健康对照者的通气和灌注更为均一,因此缺陷体素较少。表4表5 中列出了相应感兴趣区域(ROI)统计结果的报告。

灌注-通气分析示意图;肺部成像,定量灌注,区域性通气,V/Q
图830岁健康女性对照者的PREFUL参数图。 灌注(1st 行),区域通气(2nd 行),流量-容积环相关度量指标(3rd 行)以及阈值化后的V/Q图(4th 行)。注意实质值的均匀分布及较低的缺损百分比。缩写:PREFUL = 相位解析功能肺成像;V = 通气;Q = 灌注。 请点击此处以查看此图的放大版本。

利用流量-容积评分在彩色编码解剖图像中进行肺灌注与通气显像。
图960岁男性COPD患者的PREFUL参数图。 灌注(1st 行),区域通气(2nd 行),流速-容积环相关度量指标(3rd 行)以及经阈值处理的V/Q图(4th 行)。注意实质值的异质性分布及较高的缺损百分比。缩写:COPD = 慢性阻塞性肺疾病;PREFUL = 相位解析功能肺;V = 通气;Q = 灌注。 请点击此处以查看此图的放大版本。

表4:健康对照者(30岁女性)的PREFUL参数示例报告。 注意变异系数和缺陷值均较低,这与同一受试者在图8中所示的参数图一致。另见表5图9。缩写:PREFUL = 相位解析功能性肺成像。 请点击此处下载该表格。

表5:一名慢性阻塞性肺疾病患者(60岁男性)的PREFUL参数示例报告。 请注意较高的变异系数和缺陷值,这些结果与同一受试者在图9中所示的参数图一致。另见表4图8请点击此处下载该表格。

补充材料:PREFUL 算法的动画说明。请点击此处下载该文件。

讨论

关键步骤
采集过程中最常见的问题之一是信号缩放不当,这会导致在 DICOM 转换过程中因数字数据表示的精度降低而丢失信息,从而可能在后续处理阶段引发问题。另一个更为严重的问题是以交错方式采集多个层面,这会显著降低各个层面的有效时间分辨率。此外,根据层面之间的距离,还可能影响灌注对比度和定量结果,因为流入依赖于无磁化历史的新鲜自旋。在序列设置过程中需格外注意,特别是梯度强度、非对称回波、带宽和平行成像等参数。若其中任意一个参数偏离建议设置,均可能导致回波时间(TE)和时间分辨率不理想。

后处理包含多个步骤,必须按照所述顺序依次进行。例如,在低通滤波之后进行配准是没有意义的。因此,任何一个步骤的失败都会导致后续步骤的失败,这使得配准阶段尤为重要。由于并不存在单一的配准算法,不同实现方式下需要根据具体情况经验性地设置参数。若未对这些参数进行精细调整,错误的配准将导致无法生成任何有意义的结果。后处理中另一个可能耗时且关键的步骤是分割。错误的分割可能导致最终报告中的参数计算完全错误(例如,因包含了非肺部区域)。当使用深度学习算法时,这类错误分割更易发生,因为这些算法通常适应于特定的图像外观,而当应用于来自其他厂商或设备、图像外观略有差异的数据时,性能可能下降。因此,必须对分割的准确性进行视觉质量检查,并在必要时进行人工修正。

故障排除
典型的故障排除流程是逐项执行所有步骤,并检查中间结果的合理性。主要步骤的处理流程如下:首先确认图像是在自由呼吸状态下使用正确的序列和参数采集的。接着,检查信号的动态范围是否合适(肺实质中的信号强度约为 50 AU)。如果原始数据仍然可用,可使用适当的缩放因子重新重建图像,以避免重新采集数据。然后检查图像配准是否完成,且无明显伪影和残余运动。接下来,检查小感兴趣区(ROI)的时间序列是否呈现出符合预期的通气和灌注相关调制。随后,检查所应用的滤波器是否以预期方式改变了图像(例如,低通滤波后的数据中不应存在高频调制)。检查合成的呼吸和心脏周期是否符合生理特征,且无突然的跳跃。评估分割的准确性。请注意,一旦确定问题发生的主要步骤,可能需要在更精细的分辨率层级上进行进一步排查。

局限性
尽管本方案已知可产生可重复且敏感的结果,但在数据采集和后处理过程中涉及的步骤和参数众多,几乎存在无限的优化可能,且这些因素相互关联。因此,应采用自下而上的方法,首先优化序列协议(例如,关于信噪比和功能对比噪声比)。对于后续的后处理优化,采用数字肺模型作为预定义的基准真值可能具有帮助40。如文中所述,该模型可模拟自由呼吸状态下的采集过程,并包含多种类别以模拟通气/血流灌注缺陷。通过引入已知的运动所致形变,还可直接测试配准算法。尽管具有上述优势,但每个模型本质上都受限于将复杂现实映射到有限且简化模型的准确性。

本方案中提供的阈值通过经验分析发现,在健康志愿者及不同患者队列中均表现出合理的结果。然而,如前所述,根据序列、磁场强度和队列的不同,可能需要进行调整。

PREFUL 的一个普遍局限性在于其复杂的后期处理,目前尚无现成的医疗产品可用,尽管西门子医疗(Siemens Healthineers)和 BioVisioneers 已在科研合作或商业环境下为科研目的提供了初步的开发中版本。计算过程通常涉及并行处理,对 CPU 和内存提出了较高要求,可能需要现代工作站或服务器解决方案才能有效处理大量数据。此外,耗时的后期处理步骤目前阻碍了结果的即时呈现,而即时呈现对于临床工作流程而言是十分理想的。

与其他方法的比较
存在多种类似的方法,如 PREFUL,包括其前身傅里叶分解法及其其他衍生方法,例如矩阵束分解法41,以及略有不同的自门控无对比剂功能肺部 MRI(SENCEFUL MRI)42。虽然傅里叶分解法及类似方法在频域中操作,但 PREFUL 采用较宽松的傅里叶滤波,并在时域中后续计算振幅。因此,无需选择对应于通气/灌注的特定峰值。这可能使 PREFUL 对呼吸变异性具有较低的敏感性,而呼吸变异性在人体受试者中是已知存在的现象。

尽管PREFUL通过图像排序实现,但SENCEFUL采用k空间线的排序,从而提供了更高的灵活性。然而,SENCEFUL需要具备自门控能力的序列,而PREFUL则可使用常规的扰相梯度回波序列实现。类似地,在基于傅里叶分解的方法中常用的bSSFP序列以更高的信噪比和血流对比度著称,但在肺部成像时通常需要更多的优化,尤其是在3T场强下43。然而,除此之外,并无其他原因阻止将PREFUL与bSSFP采集相结合44

所有这些基于信号的方法均假设某些不希望的信号影响(包括T1、T2/T2*、扩散、层面内运动以及非正交灌注的体素)可以忽略不计。尽管PREFUL的进一步验证间接表明这些影响确实并不关键,Triphan等人却指出,其结果依赖于有效T1和TE,这可通过血液与实质成分在不同TE下的权重差异来解释45。鉴于此,bSSFP利用T2/T1对比度显示血液的初始优势,与SPGRE更简单的对比机制相比,可能为实现精确量化带来额外挑战。然而,仍需要开展更多直接探讨各种MR变量影响的研究——例如Glandorf等人针对对比剂所开展的工作46,47——因为这类研究可直接量化这些因素对PREFUL的影响。

重要性
作为一种无需屏气、无需对比剂的方法,PREFUL 与前述相关方法具有诸多共同优势:1)无需电离辐射和对比剂使用,2)无需额外的硬件或人员,3)采集过程仅需患者最低限度的配合。这些优势使 PREFUL 成为一种便捷的监测工具,尤其适用于慢性肺病儿童等脆弱人群。尽管 SPGRE 序列的信噪比较低,但该方法具有良好的可及性,且无需额外的序列编程或共享,进一步推动了其广泛应用。

正如引言部分所述,越来越多的研究显示出良好的验证性、可重复性、灵敏度结果以及监测能力,表明该技术及其相应动态参数的重要性正在不断提升,并将通过广泛传播而得到进一步支持。

披露

Filip Klimeš、Andreas Voskrebenzev 和 Jens Vogel-Claussen 是 BioVisioneers GmbH 公司的股东,该公司对肺部磁共振成像方法具有商业利益。

致谢

本工作由德国肺研究中心(DZL)资助。作者谨向所有为PREFUL的进一步开发做出贡献并提供支持的人士致以诚挚的感谢,特别是:Marcel Gutberlet、Till F. Kaireit、Lea Behrendt、Julian Glandorf、Sonja Lüdiger、Tawfik Moher Alsady、Katharina Bünemann、Marius Wernz、Robin Müller、Maximilian Zubke、Gesa Pöhler、Agilo Kern、Cristian Crisosto、Milan Speth、Arnd Obert、Julienne Scheller、Jim Wild、Edwin van Beek、Helen Marshall、Jens Gottlieb、Martha Dohna、Diane Renz、Anna-Maria Dittrich、Tobias Welte、Jens Hohlfeld、Patrick Zardo、Giles Santyr、Franz Wolfgang Hirsch、Robert Grimm、Bastian Bier、Bassem Ismail、André Fischer、Berthold Kiefer、Gregor Thoermer 和 Rebecca Ramb。此外,作者还要感谢放射科技师以及研究参与者。特别感谢汉诺威医学院放射科的 Frank Schröder 和 Sven Thiele 在实施MRI检查过程中提供的卓越技术支持。

材料

本文使用的材料清单
姓名公司目录编号评论
高级归一化工具(Advanced Normalization Tools, ANTs)-图像配准工具箱(https://stnava.github.io/ANTs/; https://github.com/fordanic/image-registration)
Forsberg-图像配准工具箱
MRI西门子医疗公司(Siemens Healthineers AG),德国慕尼黑0.55T / 1.5T / 3T 扫描仪
PREFUL 应用程序(PREFUL App)德国 BioVisioneers 有限公司,德国拉岑PREFUL 分析、图表与报告

参考文献

  1. Vogelmeier, C. F., et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive lung disease 2017 Report: GOLD Executive Summary. Eur Respir J. 49 (3), 1700214(2017).
  2. Ranu, H., Wilde, M., Madden, B. Pulmonary function tests. Ulster Med J. 80 (2), 84-90 (2011).
  3. Galbán, C. J., et al. Computed tomography-based biomarker provides unique signature for diagnosis of COPD phenotypes and disease progression. Nat Med. 18 (11), 1711-1715 (2012).
  4. Roach, P. J., Schembri, G. P., Bailey, D. L. V/Q Scanning using SPECT and SPECT/CT. J Nucl Med. 54 (9), 1588-1596 (2013).
  5. Wild, J. M., et al. MRI of the lung (1/3): methods. Insights Imaging. 3 (4), 345-353 (2012).
  6. Albert, M. S., et al. Biological magnetic resonance imaging using laser-polarized 129Xe. Nature. 370 (6486), 199-201 (1994).
  7. Berthezène, Y., et al. Contrast-enhanced MR imaging of the lung: assessments of ventilation and perfusion. Radiology. 183 (3), 667-672 (1992).
  8. Hughes, P. J. C., et al. Assessment of the influence of lung inflation state on the quantitative parameters derived from hyperpolarized gas lung ventilation MRI in healthy volunteers. J Appl Physiol. 126 (1), 183-192 (2019).
  9. Zapke, M., et al. Magnetic resonance lung function - a breakthrough for lung imaging and functional assessment? A phantom study and clinical trial. Respir Res. 7 (1), 106(2006).
  10. Klimeš, F., et al. Free-breathing quantification of regional ventilation derived by phase-resolved functional lung (PREFUL) MRI. NMR Biomed. 32 (6), e4088(2019).
  11. Deimling, M., Jellus, V., Geiger, B., Chefd'hotel, C. Time resolved lung ventilation imaging by Fourier decomposition. Proc. Intl. Soc. Mag. Reson. Med. 16, 2639(2008).
  12. Bauman, G., et al. Non-contrast-enhanced perfusion and ventilation assessment of the human lung by means of fourier decomposition in proton MRI. Magn Reson Med. 62 (3), 656-664 (2009).
  13. Bauman, G., et al. Lung ventilation- and perfusion-weighted Fourier decomposition magnetic resonance imaging: In vivo validation with hyperpolarized 3He and dynamic contrast-enhanced MRI. Magn Reson Med. 69 (1), 229-237 (2013).
  14. Bauman, G., et al. Validation of Fourier decomposition MRI with dynamic contrast-enhanced MRI using visual and automated scoring of pulmonary perfusion in young cystic fibrosis patients. Eur J Radiol. 82 (12), 2371-2377 (2013).
  15. Sommer, G., et al. Non-contrast-enhanced preoperative assessment of lung perfusion in patients with non-small-cell lung cancer using Fourier decomposition magnetic resonance imaging. Eur J Radiol. 82 (12), e879-e887 (2013).
  16. Voskrebenzev, A., Gutberlet, M., Becker, L., Wacker, F., Vogel-Claussen, J. Reproducibility of fractional ventilation derived by Fourier decomposition after adjusting for tidal volume with and without an MRI compatible spirometer. Magn Reson Med. 76 (5), 1542-1550 (2016).
  17. Voskrebenzev, A., Gutberlet, M., Kaireit, T. F., Wacker, F., Vogel-Claussen, J. Low-pass imaging of dynamic acquisitions (LIDA) with a group-oriented registration (GOREG) for proton MR imaging of lung ventilation. Magn Reson Med. 78 (4), 1496-1505 (2017).
  18. Voskrebenzev, A., et al. Feasibility of quantitative regional ventilation and perfusion mapping with phase-resolved functional lung (PREFUL) MRI in healthy volunteers and COPD, CTEPH, and CF patients. Magn Reson Med. 79 (4), 2306-2314 (2018).
  19. Voskrebenzev, A., et al. Imaging-based spirometry in chronic obstructive pulmonary disease (COPD) patients using phase resolved functional lung imaging (PREFUL). Proc. Intl. Soc. Mag. Reson. Med. 26, 1079(2018).
  20. Alsady, T. M., et al. MRI-derived regional flow-volume loop parameters detect early-stage chronic lung allograft dysfunction. J Magn Reson Imaging. 50 (6), 1873-1882 (2019).
  21. Pöhler, G. H., et al. Chronic thromboembolic pulmonary hypertension perioperative monitoring using phase-resolved functional lung (PREFUL)-MRI. J Magn Reson Imaging JMRI. 52 (2), 610-619 (2020).
  22. Behrendt, L., et al. Validation of automated perfusion-weighted phase-resolved functional lung (PREFUL)-MRI in patients with pulmonary diseases. J Magn Reson Imaging. 52 (1), 103-114 (2020).
  23. Kaireit, T. F., et al. Flow volume loop and regional ventilation assessment using phase-resolved functional lung (PREFUL) MRI: comparison with 129 xenon ventilation MRI and lung function testing. J Magn Reson Imaging JMRI. 53 (4), 1092-1105 (2021).
  24. Marshall, H., et al. 129Xe and free-breathing 1H ventilation MRI in patients with cystic fibrosis: a dual-center study. J Magn Reson Imaging. 57 (6), 1908-1921 (2023).
  25. Behrendt, L., et al. A dual center and dual vendor comparison study of automated perfusion-weighted phase-resolved functional lung magnetic resonance imaging with dynamic contrast-enhanced magnetic resonance imaging in patients with cystic fibrosis. Pulm Circ. 12 (2), e12054(2022).
  26. Moher Alsady, T., et al. Multicenter standardization of phase-resolved functional lung MRI in patients with suspected chronic thromboembolic pulmonary hypertension. J Magn Reson Imaging. , (2023).
  27. Pöhler, G. H., et al. Repeatability of phase-resolved functional lung (PREFUL)-MRI ventilation and perfusion parameters in healthy subjects and COPD patients. J Magn Reson Imaging. 53 (3), 915-927 (2021).
  28. Voskrebenzev, A., et al. PREFUL MRI depicts dual bronchodilator changes in COPD: a retrospective analysis of a randomized controlled trial. Radiol Cardiothorac Imaging. 4 (2), e210147(2022).
  29. Vogel-Claussen, J., et al. Phase-resolved functional lung (PREFUL) MRI-derived ventilation and perfusion parameters predict future lung transplant loss. Radiology. 307 (4), e221958(2023).
  30. Zanette, B., et al. Clinical feasibility of structural and functional MRI in free-breathing neonates and infants. J Magn Reson Imaging. 55 (6), 1696-1707 (2022).
  31. Dyke, J. P., et al. Assessment of lung ventilation of premature infants with bronchopulmonary dysplasia at 1.5 Tesla using phase-resolved functional lung magnetic resonance imaging. Pediatr Radiol. 53 (6), 1076-1084 (2023).
  32. Glandorf, J., et al. Comparison of phase-resolved functional lung (PREFUL) MRI derived perfusion and ventilation parameters at 1.5T and 3T in healthy volunteers. PLOS ONE. 15 (12), e0244638(2020).
  33. Lévy, S., et al. Free-breathing low-field MRI of the lungs detects functional alterations associated with persistent symptoms after COVID-19 infection. Invest Radiol. 57 (11), 742-751 (2022).
  34. Avants, B. B., et al. A reproducible evaluation of ANTs similarity metric performance in brain image registration. NeuroImage. 54 (3), 2033-2044 (2011).
  35. Forsberg, D., Andersson, M., Knutsson, H. Extending image registration using polynomial expansion to diffeomorphic deformations. SSBA Symposium on Image Analysis. , Stockholm, Sweden. (2012).
  36. Forsberg, D. fordanic/image-registration. , Available from: https://github.com/fordanic/image-registration (2022).
  37. Klimeš, F., et al. Evaluation of image registration algorithms for 3D phase-resolved functional lung ventilation magnetic resonance imaging in healthy volunteers and chronic obstructive pulmonary disease patients. NMR Biomed. 36 (3), e4860(2023).
  38. He, K., Sun, J., Tang, X. Guided image filtering. IEEE Trans Pattern Anal Mach Intell. 35 (6), 1397-1409 (2013).
  39. Glandorf, J., et al. Perfusion quantification using voxel-wise proton density and median signal decay in PREFUL MRI. Magn Reson Med. 86 (3), 1482-1493 (2021).
  40. Voskrebenzev, A., Gutberlet, M., Klimeš, F., Wacker, F., Vogel-Claussen, J. Introduction of a Digital Lung Model for Validation and Refinement of Functional Lung Imaging Methods. Proc. Intl. Soc. Mag. Reson. Med. 31, 4820(2023).
  41. Bauman, G., Bieri, O. Matrix pencil decomposition of time-resolved proton MRI for robust and improved assessment of pulmonary ventilation and perfusion. Magn Reson Med. 77 (1), 336-342 (2017).
  42. Fischer, A., et al. SElf-gated Non-Contrast-Enhanced FUnctional Lung imaging (SENCEFUL) using a quasi-random fast low-angle shot (FLASH) sequence and proton MRI. NMR Biomed. 27 (8), 907-917 (2014).
  43. Bauman, G., Pusterla, O., Bieri, O. Functional lung imaging with transient spoiled gradient echo. Magn Reson Med. 81 (3), 1915-1923 (2019).
  44. Rotärmel, A., et al. GRE bSSFP vs. FLASH based Fourier decomposition lung MRI at 1.5T: evaluation of image quality, fractional ventilation and lung perfusion in healthy volunteers. Prog Intl Soc Mag Reson. 26, 2-4 (2018).
  45. Triphan, S. M. F., et al. Echo time-dependent observed lung T1 in patients with chronic obstructive pulmonary disease in correlation with quantitative imaging and clinical indices. J Magn Reson Imaging. 54 (5), 1562-1571 (2021).
  46. Glandorf, J., et al. Effect of intravenously injected gadolinium-based contrast agents on functional lung parameters derived by PREFUL MRI. Magnetic Resonance in Medicine. 83 (3), 1045-1054 (2020).
  47. Glandorf, J., et al. Influence of gadolinium, field-strength and sequence type on quantified perfusion values in phase-resolved functional lung MRI. PloS One. 18 (8), e0288744(2023).

重印与许可

标签

3D PREFUL