研究超声背散射零差 K 模型参数成像(简称零差 K 成像)监测微波消融凝固区的可行性,对比 2 种最新的零差 K 模型参数估算方法,即 RSK 法(基于信号包络振幅的信噪比、偏度和峰度的方法)和 XU 法(基于信号强度的一阶矩、X 统计量和 U 统计量的方法)。首先,将离体猪肝微波消融实验中采集的超声背散射信号经噪声辅助互相关算法、包络检测、滑动窗口法、数字扫描变换和颜色映射等操作得到零差 K 成像。然后利用 20 例离体猪肝微波消融实验来评估零差 K 成像监测凝固区的效果。结果显示 RSK 法和 XU 法的受试者工作特征曲线下面积平均值 ± 标准差分别为 0.77 ± 0.06、0.83 ± 0.08。RSK 法和 XU 法监测凝固区的精度平均值 ± 标准差分别为(86 ± 10)% 和(90 ± 8)%。XU 法估算凝固区面积值与实际猪肝组织凝固区面积的 Bland-Altman 一致性强于 RSK 法。XU 法参数估算并成像所耗时间少于 RSK 法。研究显示超声背散射零差 K 成像可用于监测微波消融过程中凝固区的变化,XU 法优于 RSK 法。
引用本文: 宋爽, 张英华, 周著黄, 吴水才. 超声背散射零差 K 成像监测微波消融. 生物医学工程学杂志, 2021, 38(3): 520-527. doi: 10.7507/1001-5515.202003032 复制
引言
肝肿瘤微波消融已应用于临床[1-3]。微波消融治疗是在医学影像的引导下将消融针经皮刺入肿瘤内,通过消融针发射的电磁微波搅动周围组织中的水分子,产生摩擦和热量,使得肿瘤组织凝固性坏死,从而达到原位治疗肿瘤的目的[4]。
为了最大范围地杀灭肿瘤同时保护正常肝组织,需要通过医学影像对微波消融治疗进行术中实时监测。超声成像因实时、无创、廉价、易普及等特点,已用于肝肿瘤微波消融治疗的影像引导和术中监测,其中 B 模式超声(B 超)在临床应用最广。B 超利用接收到的背散射回波信号的幅度信息进行成像,由于消融过程中组织受热产生的微气泡会导致伪影,且病灶残存肿瘤组织与凝固性坏死组织即凝固区的回声相近,故仅用传统 B 超成像难以确定肿瘤是否消融完全[5]。
声学上,生物组织可建模为一系列散射声波的微小粒子即散射子的集合,超声背散射信号中蕴含着散射子数目、尺寸、排列等信息。超声背散射信号的概率分布模式与组织的微结构(散射子)有关[6-7],因此超声背散射统计成像可用于定量评估组织微结构的变化[8-10]。零差 K 分布是最具物理意义的超声背散射信号概率统计模型[11-12],其参数可反映超声分辨单元内的有效散射子个数以及相干散射和弥漫散射成分。但由于零差 K 模型参数估算比较复杂,目前有关超声背散射零差 K 模型参数成像(简称零差 K 成像)的研究相对较少,主要包括乳腺肿瘤良恶性分类[13]、脂肪肝评估[12,14]、皮肤组织中溃疡愈合的评估[15]、颈动脉斑块成分的评估[16]等。
微波消融过程中,组织受热发生凝固并产生微气泡,导致组织微结构发生变化,引起凝固区组织的散射子状态发生改变,而零差 K 模型参数可定征生物组织的散射子个数及相干/弥漫散射特性,因此可以利用超声背散射零差 K 成像,在微波消融过程中对凝固区进行监测。本文基于离体猪肝微波消融实验采集到的超声背散射信号,使用两种最新方法[17-18]估算超声背散射零差 K 模型参数,并通过受试者工作特征(receiver operating characteristic,ROC)曲线、监测精度、Bland-Altman 一致性检验和算法耗时等性能指标,对比验证超声背散射零差 K 成像监测微波消融的效果。
1 研究方法
1.1 实验数据的采集
1.1.1 微波消融实验平台
离体猪肝微波消融实验平台如图 1 所示。实验设备包括:微波消融仪(型号:KY-2000,南京康友公司),超声仪(型号:Terason T3000,美国 Terason 公司),线阵换能器(型号:12L5A,中心频率为 7.5 MHz,通道数为 256,脉冲长度约为 0.7 mm[19]),水冷式微波消融针(型号:EN868-5),用于固定超声探头的支架,用于盛放猪肝组织的亚克力盒(尺寸:6 cm × 6 cm × 6 cm)等。
1.1.2 超声数据的获取
实验前,将新鲜猪肝组织置于 0.9% 的氯化钠溶液中静置 2 min,以去除猪肝组织内气体。为了降低亚克力盒底部的强回声,在亚克力盒的底部放置了厚度为 1 cm 的纸板和 1 cm 的胶皮垫。
实验开始,将去气的猪肝组织放到亚克力盒中。为防止实验过程中氯化钠溶液从盒子侧壁预留用于插针的圆孔流出,使用双面胶对圆孔进行封堵,然后将消融针穿过双面胶和圆孔后水平插入猪肝组织中,再将 0.9% 的氯化钠溶液注入到亚克力盒中。将线阵换能器沿垂直于消融针方向,固定于猪肝组织上方并保证换能器与氯化钠溶液充分接触(见图 1),此时利用超声仪自带软件可以对猪肝组织进行 B 超成像,利用消融针在 B 超中的高回声来定位消融针位置,如图 2 所示,图中红色箭头指向消融针。将微波消融仪设定好加热参数(消融功率 P = 80 W,消融时间 t = 60 s)开始微波消融实验。实验过程中,利用我们自行编写的软件[9](在超声仪中运行)以 2 帧/秒的速率采集超声背散射信号并保存在超声仪硬盘中,所得超声背散射信号矩阵为 256(扫描线数)× 1 558(采样点数)。
1.2 超声背散射零差 K 成像
1.2.1 理论分析
零差 K 模型由 Jakeman[20]首次提出。估算超声背散射零差 K 模型参数的主要方法有两种,即 RSK 法[17](基于信号包络振幅的信噪比、偏度和峰度的方法)和 XU 法[18](基于信号强度的一阶矩、X 统计量和 U 统计量的方法),但目前还无文献报道 RSK 法和 XU 法零差 K 成像在监测微波消融方面的效果对比。为此,本文对比研究这两种方法监测微波消融的效果。因微波消融过程中加热产生的微气泡导致散射子数目发生改变[10,21],因此使用零差 K 模型的有效散射子数目参数,即 μ(RSK 法参数)或 α(XU 法参数),来监测消融凝固区的变化。
(1)RSK 法:RSK 法的具体流程可参阅文献[17]。它借助超声背散射信号的信噪比、偏度和峰度来估算参数(μ,s2,),其中,μ 表示超声分辨单元内的有效散射子数目,s2表示相干散射信号的能量,表示弥漫散射信号的能量。信噪比 R、偏度 S 和峰度 K 与信号包络振幅 A 的任意阶矩 v 的关系见式(1)~(3),式中 E[·]表示数学期望。
(2)XU 法:XU 法的具体流程可参阅文献[18]。它借助超声背散射信号的强度、U 统计和 X 统计来估算参数(α,,),其中 α 是散射子聚集参数,可以反映超声分辨率单元内有效散射子的个数,表示相干信号能量,表示弥漫散射信号的能量。U 统计和 X 统计的定义如式 (4)~(5),式中,I 表示背散射信号的强度,I = A2。
1.2.2 超声背散射零差 K 成像
超声参数成像可直观显示组织定征结果[12]。为了将超声背散射零差 K 成像用于猪肝组织微波消融消融监测,我们通过 Matlab 2016a 软件实现超声背散射信号零差 K 分布参数值估算及参数成像。
超声背散射零差 K 成像流程如图 3 所示:(1)对采集到的原始超声背散射信号进行希尔伯特变换,取模后得到所监测组织的包络信号。包络信号经对数压缩和数字扫描变换后可重建出 B 超成像。(2)将采集到的原始超声背散射信号进行噪声辅助互相关算法(noise-assisted correlation algorithm,NCA)[19,22]预处理后,再进行包络检测。然后利用滑动窗口法[12,19]估算超声背散射零差 K 模型散射子数目参数(μ 和 α)。最后,对数字扫描变换后得到的估算值矩阵使用颜色映射,得到 μ 或 α 参数图像,即为超声零差 K 成像。
由于在无回声区域存在的噪声信号会导致零差 K 成像伪影,为避免噪声的影响,我们引入 NCA[19,22]。NCA 是将随机人工白噪声分 2 次分别加入每帧超声背散射信号,然后计算这 2 帧包含额外噪声的超声背散射信号的相关系数,最后将相关系数小于阈值所对应的背散射信号的振幅设为零,即用零值代替无回声区域的噪声信号。本文通过 Matlab 软件中的 awgn 函数向原始超声射频信号中添加白噪声,awgn 函数的输入参数 snr 设置为 10。将阈值设置为 80%,将幅值为零的超声背散射信号估算出的零差 K 模型参数也设置为零,就能得到经 NCA 处理后的超声零差 K 成像。
滑动窗口是用一个方形窗口在包络信号中滑动,通过窗口内的超声背散射数据估计 μ 或 α 值,这些估值被赋给窗口中心的像素。窗口以一定的重叠率滑过整个包络信号后,得到 μ 或 α 的估算值矩阵。由于使用滑动窗口算法得到的估算值矩阵的尺寸会小于原始包络信号矩阵的尺寸,因此还需要对 μ 和 α 的估算值矩阵进行数字扫描变换。滑动窗口的大小会影响超声零差 K 成像的分辨率和参数估值的稳定性[19],较大的窗口可获得更稳定的参数估算,但会降低成像分辨率;较小的窗口可获得更高的成像分辨率,但会降低参数估算的稳定性。综合考虑成像的分辨率和参数估算的稳定性,我们参阅文献[12,19]并反复实验,最终选择窗口边长等于入射超声脉冲长度的 3 倍,窗口的重叠率为 50%。
2 实验结果
2.1 B 超成像监测凝固区
图 4 是利用离体猪肝微波消融实验消融过程中(0~60 s)和消融停止后(80~120 s)的超声背散射信号重建出的 B 超图像。消融开始后,组织因受热产生的微气泡在消融针附近形成高回声区,使得图像亮度增强,同时造成消融针的后方声影。消融停止后,由于微气泡逐渐消散,B 超图像亮度有所减弱。通过图 4 的 B 超图像,我们难以确定凝固区的大小及边界。
2.2 RSK 法零差 K 成像监测凝固区
图 5 为对应图 4 的 RSK 法 μ 参数成像。μ 参数图像右侧的色阶条中,颜色从蓝到黄到红,代表 μ 参数的值逐渐增大,对应物理意义是超声波分辨单元内的有效散射子个数从 0 逐渐增加(因其表示个数,因此是无单位的物理量),为了提高图像的对比度,我们将色阶条的动态范围设为 0~4;XU 法的 α 参数成像同理。通过图 5 可看出,从消融开始到消融停止(0~60 s),消融针附近 μ 参数值逐渐增大,μ 参数成像可以反映凝固区的变化。消融停止后(80~120 s),μ 参数值有所减小。
2.3 XU 法零差 K 成像监测凝固区
图 6 为对应图 4 的 XU 法 α 参数成像。与 RSK 法相似,通过图 6 可以看出,从消融开始到消融停止(0~60 s),消融针附近 α 参数值逐渐增大,消融停止后(80~120 s),α 参数值有所减小。
2.4 监测效果评估
采集 20 例离体猪肝微波消融实验的超声背散射信号,通过对比超声零差 K 成像的 ROC 曲线,估算超声零差 K 成像监测凝固区的精度,比较超声零差 K 成像估算凝固区与实际猪肝组织凝固区的 Bland-Altman 一致性来评估超声零差 K 成像监测微波消融的效果。
离体猪肝微波消融实验结束后,将猪肝组织沿超声扫描平面切开的大体病理图像(图 7a)作为金标准,根据大体病理图像,测量并记录组织消融凝固区的长轴和短轴。通过 ROC 评估 RSK 法和 XU 法零差 K 成像监测凝固区的性能,ROC 曲线下面积(area under the ROC curve,AUROC)用作性能指标。以消融针所在位置为中心,在大体病理图像、μ 参数矩阵和 α 参数矩阵中选取一个 30 mm × 30 mm 的感兴趣区(region of interest,ROI),如图 7 中红色虚线框分别为 ROI(g)、ROI(μ)、ROI(α)。
将 ROI(g)二值化得到 ROI-1(g),对 ROI(μ) 和 ROI(α) 内矩阵插值或重采样得到与组织剖面像素点相同的参数矩阵,ROI-1(μ) 和 ROI-1(α)。将 ROI-1(g)与 ROI-1(μ) 和 ROI-1(α) 逐像素预测得到 ROC 曲线,如图 8 所示,图中 RSK 法 μ 参数和 XU 法 α 参数的 AUROC 分别为 0.72、0.75。20 例猪肝组织微波消融实验 RSK 法和 XU 法的 AUROC 平均值 ± 标准差分别为 0.77 ± 0.06、0.83 ± 0.08。
为了在零差 K 成像中识别出消融凝固区的范围,估算凝固区面积值,我们引入多项式拟合(polynomial approximation,PAX)技术[19,23]。用 Ix,z 代表参数成像中的像素点,其中 x、z 分别是成像的横向宽度和轴向深度。PAX 流程是首先假设参数成像每个方向上的函数都是一个 p 阶多项式,通过最小二乘法确定最优多项式,然后用由最优多项式计算得到的值替代 Ix,z。沿轴向和横向的每条线拟合后重建就得到 PAX 成像,将 PAX 与等高线方法相结合,估算出消融凝固区的范围,如图 9 所示。图 9a 和 9b 为与消融停止时刻 μ 和 α 参数成像(图 7b 和 7c)相对应的 PAX 成像,其中黑色实线包围的区域即为估算的消融凝固区。实验发现 PAX 阶数设置为 5 阶、等高线设置为-6 dB 时监测效果最佳。
利用式 (6) 计算出实际猪肝组织凝固区面积,μ 参数 PAX 成像估算出凝固区面积 SRSK,α 参数 PAX 成像估算出凝固区面积 SXU。
式中,S 为凝固区面积,l1 和 l2 为凝固区的长轴和短轴。根据式 (7) 可得到超声零差 K 参数 PAX 成像监测凝固区的精度。
式中,Accuracy 为监测精度;SHK 为 PAX 成像估算出的凝固区面积值,即 SRSK 或 SXU;Sgold 为猪肝组织的实际凝固区面积值。RSK 法和 XU 法对 20 例猪肝组织凝固区监测的精度平均值 ± 标准差分别为(86 ± 10)% 和(90 ± 8)%。
利用 Bland-Altman 回归分析,以 Sgold 为金标准,可以比较超声零差 K 成像估算凝固区的一致性。针对现有的 20 例离体猪肝数据,以面积值 Sgold 为横轴,SRSK 或 SXU 与 Sgold 的差值为纵轴绘制 Bland-Altman 图,如图10 所示。RSK 法和 XU 法的 Bland-Altman 图中只有 5%(1/20)差值点位于 95% 一致性界限区间([mean − 1.96*STD,mean + 1.96*STD],mean 和 STD 分别代表均值和标准差)的外侧,说明 SRSK 与 Sgold、SXU 与 Sgold 都具有一致性,但 XU 法中均值线比 RSK 法更加接近参考线 reference line(代表差值平均数为 0 的线),且 XU 法的 95% 一致性界限的区间范围更小。
通过实验还发现,对 1 例猪肝组织消融停止时刻(t=60 s 处)的 1 帧超声背散射信号,利用 RSK 法估算得到参数成像和 PAX 成像所耗时间为 81 s,利用 XU 法估算得到参数成像和 PAX 成像所耗时间为 17 s。
3 讨论
本文将超声背散射零差 K 成像用于监测微波消融,并对比了 2 种最新的零差 K 模型参数估算方法,即 RSK 法和 XU 法。通过离体猪肝微波消融实验验证了超声背散射零差 K 成像监测微波消融的可行性。
目前用于热消融定量监测的超声技术主要分为超声弹性成像和定量超声两大类[10,21,23-29]。超声弹性成像是对组织施加一个外部或内部的激励,然后利用超声波对组织产生的响应进行跟踪,估算出组织产生的应变、剪切波传播速度或弹性模量[24-26]。定量超声技术从超声背散射信号中估算出组织的声学参数,主要包括超声背散射统计成像[9-10,19]、散射子平均间距成像[27]、超声衰减系数成像[28]、超声回波去相关成像[23]以及基于卷积神经网络结构的超声成像等[29]。超声背散射统计成像中的超声 Nakagami 成像已用于消融术中监测,但 Nakagami 模型 m 参数在散射子数密度较高时存在饱和现象[8]。零差 K 模型能克服这一不足,但因其算法复杂,目前尚未应用于微波消融术中监测。
为避免无回声区的噪声在超声零差 K 成像中形成的伪影,本文引入 NCA 对超声背散射信号进行预处理。NCA 已应用于超声 Nakagami 成像,可有效消除成像无回声区域的伪影,同时不牺牲 Nakagami 成像的帧率[22]。NCA 的原理[21]是添加的白噪声幅值很小,对超声背散射信号没有明显的影响,因此对同一个背散射信号添加 2 次白噪声,得到的 2 条含额外噪声的超声背散射信号,其对应散射介质区域的波形几乎相同,它们之间的相关系数接近 1。相比之下,由于无回声区自身的噪声信号幅度较小,而每次添加的白噪声都具有随机性且其对无回声区信号作用较为明显,因此在无回声区内,2 条含额外噪声的超声背散射信号的相关系数会变小,即小于 1。通过设置相关阈值,NCA 可以自动区分无回声区域的白噪声和散射介质中的背散射信号。
将图 5、图 6 与图 4 对比发现,在监测组织微波消融时,相比于传统 B 超成像,超声零差 K 成像可以更好地监测消融过程中组织凝固区的变化,但当消融停止后,超声零差 K 成像不能用于消融结束的凝固区检测评估。这是由于微波消融过程中产生的微气泡可视为有效散射子,当消融开始后,随着微气泡的产生,组织消融区内有效散射子个数增加,超声零差 K 成像消融区亮度增强,当消融停止后,随着微气泡的消散,组织消融区的有效散射子个数减少,超声零差 K 成像消融区亮度逐渐减弱。
本文通过 20 例猪肝组织微波消融的超声背散射信号,对比验证 RSK 法和 XU 法超声零差 K 成像监测凝固区的效果。结果显示 RSK 法的平均 AUROC 和平均监测精度均小于 XU 法。XU 法估算的凝固区面积值与实际猪肝组织凝固区 Bland-Altman 一致性更强。XU 法估算参数并成像所耗时间更短。由此我们可以得出结论,零差 K 成像可用于监测猪肝组织微波消融过程中凝固区的变化,XU 法优于 RSK 法。
本研究尚有不足之处。首先,由于算法是利用 Matlab 实现,所以比较耗时,未来需考虑 C/C++语言或其他加速算法,以利实时监测。其次,本研究只对离体组织的微波消融进行监测,且实验数据只有 20 例,未来还需加大实验数据并考虑通过在体实验验证算法的有效性。另外,本研究中的零差 K 成像仅对消融结束时刻的监测精度和凝固区面积的 Bland-Altman 一致性进行分析,原因是其他时刻的参考标准不易获得,未来可在不同的消融时刻停止消融并剖开猪肝进行定量分析,以评估本文方法在不同时刻定量监测凝固区的效果。
4 结论
超声背散射零差 K 成像可用于监测离体猪肝微波消融过程中凝固区的变化。RSK 法和 XU 法的 AUROC 平均值 ± 标准差分别为 0.77 ± 0.06、0.83 ± 0.08。基于 RSK 法和 XU 法的零差 K 成像监测凝固区的精度平均值 ± 标准差分别为(86 ± 10)% 和(90 ± 8)%。XU 法估算凝固区面积值与实际猪肝组织凝固区面积的 Bland-Altman 一致性强于 RSK 法。另外,XU 法估算参数并成像所耗时间也少于 RSK 法。因此,在监测肝组织微波消融凝固区时,XU 法优于 RSK 法。
利益冲突声明:本文全体作者均声明不存在利益冲突。
引言
肝肿瘤微波消融已应用于临床[1-3]。微波消融治疗是在医学影像的引导下将消融针经皮刺入肿瘤内,通过消融针发射的电磁微波搅动周围组织中的水分子,产生摩擦和热量,使得肿瘤组织凝固性坏死,从而达到原位治疗肿瘤的目的[4]。
为了最大范围地杀灭肿瘤同时保护正常肝组织,需要通过医学影像对微波消融治疗进行术中实时监测。超声成像因实时、无创、廉价、易普及等特点,已用于肝肿瘤微波消融治疗的影像引导和术中监测,其中 B 模式超声(B 超)在临床应用最广。B 超利用接收到的背散射回波信号的幅度信息进行成像,由于消融过程中组织受热产生的微气泡会导致伪影,且病灶残存肿瘤组织与凝固性坏死组织即凝固区的回声相近,故仅用传统 B 超成像难以确定肿瘤是否消融完全[5]。
声学上,生物组织可建模为一系列散射声波的微小粒子即散射子的集合,超声背散射信号中蕴含着散射子数目、尺寸、排列等信息。超声背散射信号的概率分布模式与组织的微结构(散射子)有关[6-7],因此超声背散射统计成像可用于定量评估组织微结构的变化[8-10]。零差 K 分布是最具物理意义的超声背散射信号概率统计模型[11-12],其参数可反映超声分辨单元内的有效散射子个数以及相干散射和弥漫散射成分。但由于零差 K 模型参数估算比较复杂,目前有关超声背散射零差 K 模型参数成像(简称零差 K 成像)的研究相对较少,主要包括乳腺肿瘤良恶性分类[13]、脂肪肝评估[12,14]、皮肤组织中溃疡愈合的评估[15]、颈动脉斑块成分的评估[16]等。
微波消融过程中,组织受热发生凝固并产生微气泡,导致组织微结构发生变化,引起凝固区组织的散射子状态发生改变,而零差 K 模型参数可定征生物组织的散射子个数及相干/弥漫散射特性,因此可以利用超声背散射零差 K 成像,在微波消融过程中对凝固区进行监测。本文基于离体猪肝微波消融实验采集到的超声背散射信号,使用两种最新方法[17-18]估算超声背散射零差 K 模型参数,并通过受试者工作特征(receiver operating characteristic,ROC)曲线、监测精度、Bland-Altman 一致性检验和算法耗时等性能指标,对比验证超声背散射零差 K 成像监测微波消融的效果。
1 研究方法
1.1 实验数据的采集
1.1.1 微波消融实验平台
离体猪肝微波消融实验平台如图 1 所示。实验设备包括:微波消融仪(型号:KY-2000,南京康友公司),超声仪(型号:Terason T3000,美国 Terason 公司),线阵换能器(型号:12L5A,中心频率为 7.5 MHz,通道数为 256,脉冲长度约为 0.7 mm[19]),水冷式微波消融针(型号:EN868-5),用于固定超声探头的支架,用于盛放猪肝组织的亚克力盒(尺寸:6 cm × 6 cm × 6 cm)等。
1.1.2 超声数据的获取
实验前,将新鲜猪肝组织置于 0.9% 的氯化钠溶液中静置 2 min,以去除猪肝组织内气体。为了降低亚克力盒底部的强回声,在亚克力盒的底部放置了厚度为 1 cm 的纸板和 1 cm 的胶皮垫。
实验开始,将去气的猪肝组织放到亚克力盒中。为防止实验过程中氯化钠溶液从盒子侧壁预留用于插针的圆孔流出,使用双面胶对圆孔进行封堵,然后将消融针穿过双面胶和圆孔后水平插入猪肝组织中,再将 0.9% 的氯化钠溶液注入到亚克力盒中。将线阵换能器沿垂直于消融针方向,固定于猪肝组织上方并保证换能器与氯化钠溶液充分接触(见图 1),此时利用超声仪自带软件可以对猪肝组织进行 B 超成像,利用消融针在 B 超中的高回声来定位消融针位置,如图 2 所示,图中红色箭头指向消融针。将微波消融仪设定好加热参数(消融功率 P = 80 W,消融时间 t = 60 s)开始微波消融实验。实验过程中,利用我们自行编写的软件[9](在超声仪中运行)以 2 帧/秒的速率采集超声背散射信号并保存在超声仪硬盘中,所得超声背散射信号矩阵为 256(扫描线数)× 1 558(采样点数)。
1.2 超声背散射零差 K 成像
1.2.1 理论分析
零差 K 模型由 Jakeman[20]首次提出。估算超声背散射零差 K 模型参数的主要方法有两种,即 RSK 法[17](基于信号包络振幅的信噪比、偏度和峰度的方法)和 XU 法[18](基于信号强度的一阶矩、X 统计量和 U 统计量的方法),但目前还无文献报道 RSK 法和 XU 法零差 K 成像在监测微波消融方面的效果对比。为此,本文对比研究这两种方法监测微波消融的效果。因微波消融过程中加热产生的微气泡导致散射子数目发生改变[10,21],因此使用零差 K 模型的有效散射子数目参数,即 μ(RSK 法参数)或 α(XU 法参数),来监测消融凝固区的变化。
(1)RSK 法:RSK 法的具体流程可参阅文献[17]。它借助超声背散射信号的信噪比、偏度和峰度来估算参数(μ,s2,),其中,μ 表示超声分辨单元内的有效散射子数目,s2表示相干散射信号的能量,表示弥漫散射信号的能量。信噪比 R、偏度 S 和峰度 K 与信号包络振幅 A 的任意阶矩 v 的关系见式(1)~(3),式中 E[·]表示数学期望。
(2)XU 法:XU 法的具体流程可参阅文献[18]。它借助超声背散射信号的强度、U 统计和 X 统计来估算参数(α,,),其中 α 是散射子聚集参数,可以反映超声分辨率单元内有效散射子的个数,表示相干信号能量,表示弥漫散射信号的能量。U 统计和 X 统计的定义如式 (4)~(5),式中,I 表示背散射信号的强度,I = A2。
1.2.2 超声背散射零差 K 成像
超声参数成像可直观显示组织定征结果[12]。为了将超声背散射零差 K 成像用于猪肝组织微波消融消融监测,我们通过 Matlab 2016a 软件实现超声背散射信号零差 K 分布参数值估算及参数成像。
超声背散射零差 K 成像流程如图 3 所示:(1)对采集到的原始超声背散射信号进行希尔伯特变换,取模后得到所监测组织的包络信号。包络信号经对数压缩和数字扫描变换后可重建出 B 超成像。(2)将采集到的原始超声背散射信号进行噪声辅助互相关算法(noise-assisted correlation algorithm,NCA)[19,22]预处理后,再进行包络检测。然后利用滑动窗口法[12,19]估算超声背散射零差 K 模型散射子数目参数(μ 和 α)。最后,对数字扫描变换后得到的估算值矩阵使用颜色映射,得到 μ 或 α 参数图像,即为超声零差 K 成像。
由于在无回声区域存在的噪声信号会导致零差 K 成像伪影,为避免噪声的影响,我们引入 NCA[19,22]。NCA 是将随机人工白噪声分 2 次分别加入每帧超声背散射信号,然后计算这 2 帧包含额外噪声的超声背散射信号的相关系数,最后将相关系数小于阈值所对应的背散射信号的振幅设为零,即用零值代替无回声区域的噪声信号。本文通过 Matlab 软件中的 awgn 函数向原始超声射频信号中添加白噪声,awgn 函数的输入参数 snr 设置为 10。将阈值设置为 80%,将幅值为零的超声背散射信号估算出的零差 K 模型参数也设置为零,就能得到经 NCA 处理后的超声零差 K 成像。
滑动窗口是用一个方形窗口在包络信号中滑动,通过窗口内的超声背散射数据估计 μ 或 α 值,这些估值被赋给窗口中心的像素。窗口以一定的重叠率滑过整个包络信号后,得到 μ 或 α 的估算值矩阵。由于使用滑动窗口算法得到的估算值矩阵的尺寸会小于原始包络信号矩阵的尺寸,因此还需要对 μ 和 α 的估算值矩阵进行数字扫描变换。滑动窗口的大小会影响超声零差 K 成像的分辨率和参数估值的稳定性[19],较大的窗口可获得更稳定的参数估算,但会降低成像分辨率;较小的窗口可获得更高的成像分辨率,但会降低参数估算的稳定性。综合考虑成像的分辨率和参数估算的稳定性,我们参阅文献[12,19]并反复实验,最终选择窗口边长等于入射超声脉冲长度的 3 倍,窗口的重叠率为 50%。
2 实验结果
2.1 B 超成像监测凝固区
图 4 是利用离体猪肝微波消融实验消融过程中(0~60 s)和消融停止后(80~120 s)的超声背散射信号重建出的 B 超图像。消融开始后,组织因受热产生的微气泡在消融针附近形成高回声区,使得图像亮度增强,同时造成消融针的后方声影。消融停止后,由于微气泡逐渐消散,B 超图像亮度有所减弱。通过图 4 的 B 超图像,我们难以确定凝固区的大小及边界。
2.2 RSK 法零差 K 成像监测凝固区
图 5 为对应图 4 的 RSK 法 μ 参数成像。μ 参数图像右侧的色阶条中,颜色从蓝到黄到红,代表 μ 参数的值逐渐增大,对应物理意义是超声波分辨单元内的有效散射子个数从 0 逐渐增加(因其表示个数,因此是无单位的物理量),为了提高图像的对比度,我们将色阶条的动态范围设为 0~4;XU 法的 α 参数成像同理。通过图 5 可看出,从消融开始到消融停止(0~60 s),消融针附近 μ 参数值逐渐增大,μ 参数成像可以反映凝固区的变化。消融停止后(80~120 s),μ 参数值有所减小。
2.3 XU 法零差 K 成像监测凝固区
图 6 为对应图 4 的 XU 法 α 参数成像。与 RSK 法相似,通过图 6 可以看出,从消融开始到消融停止(0~60 s),消融针附近 α 参数值逐渐增大,消融停止后(80~120 s),α 参数值有所减小。
2.4 监测效果评估
采集 20 例离体猪肝微波消融实验的超声背散射信号,通过对比超声零差 K 成像的 ROC 曲线,估算超声零差 K 成像监测凝固区的精度,比较超声零差 K 成像估算凝固区与实际猪肝组织凝固区的 Bland-Altman 一致性来评估超声零差 K 成像监测微波消融的效果。
离体猪肝微波消融实验结束后,将猪肝组织沿超声扫描平面切开的大体病理图像(图 7a)作为金标准,根据大体病理图像,测量并记录组织消融凝固区的长轴和短轴。通过 ROC 评估 RSK 法和 XU 法零差 K 成像监测凝固区的性能,ROC 曲线下面积(area under the ROC curve,AUROC)用作性能指标。以消融针所在位置为中心,在大体病理图像、μ 参数矩阵和 α 参数矩阵中选取一个 30 mm × 30 mm 的感兴趣区(region of interest,ROI),如图 7 中红色虚线框分别为 ROI(g)、ROI(μ)、ROI(α)。
将 ROI(g)二值化得到 ROI-1(g),对 ROI(μ) 和 ROI(α) 内矩阵插值或重采样得到与组织剖面像素点相同的参数矩阵,ROI-1(μ) 和 ROI-1(α)。将 ROI-1(g)与 ROI-1(μ) 和 ROI-1(α) 逐像素预测得到 ROC 曲线,如图 8 所示,图中 RSK 法 μ 参数和 XU 法 α 参数的 AUROC 分别为 0.72、0.75。20 例猪肝组织微波消融实验 RSK 法和 XU 法的 AUROC 平均值 ± 标准差分别为 0.77 ± 0.06、0.83 ± 0.08。
为了在零差 K 成像中识别出消融凝固区的范围,估算凝固区面积值,我们引入多项式拟合(polynomial approximation,PAX)技术[19,23]。用 Ix,z 代表参数成像中的像素点,其中 x、z 分别是成像的横向宽度和轴向深度。PAX 流程是首先假设参数成像每个方向上的函数都是一个 p 阶多项式,通过最小二乘法确定最优多项式,然后用由最优多项式计算得到的值替代 Ix,z。沿轴向和横向的每条线拟合后重建就得到 PAX 成像,将 PAX 与等高线方法相结合,估算出消融凝固区的范围,如图 9 所示。图 9a 和 9b 为与消融停止时刻 μ 和 α 参数成像(图 7b 和 7c)相对应的 PAX 成像,其中黑色实线包围的区域即为估算的消融凝固区。实验发现 PAX 阶数设置为 5 阶、等高线设置为-6 dB 时监测效果最佳。
利用式 (6) 计算出实际猪肝组织凝固区面积,μ 参数 PAX 成像估算出凝固区面积 SRSK,α 参数 PAX 成像估算出凝固区面积 SXU。
式中,S 为凝固区面积,l1 和 l2 为凝固区的长轴和短轴。根据式 (7) 可得到超声零差 K 参数 PAX 成像监测凝固区的精度。
式中,Accuracy 为监测精度;SHK 为 PAX 成像估算出的凝固区面积值,即 SRSK 或 SXU;Sgold 为猪肝组织的实际凝固区面积值。RSK 法和 XU 法对 20 例猪肝组织凝固区监测的精度平均值 ± 标准差分别为(86 ± 10)% 和(90 ± 8)%。
利用 Bland-Altman 回归分析,以 Sgold 为金标准,可以比较超声零差 K 成像估算凝固区的一致性。针对现有的 20 例离体猪肝数据,以面积值 Sgold 为横轴,SRSK 或 SXU 与 Sgold 的差值为纵轴绘制 Bland-Altman 图,如图10 所示。RSK 法和 XU 法的 Bland-Altman 图中只有 5%(1/20)差值点位于 95% 一致性界限区间([mean − 1.96*STD,mean + 1.96*STD],mean 和 STD 分别代表均值和标准差)的外侧,说明 SRSK 与 Sgold、SXU 与 Sgold 都具有一致性,但 XU 法中均值线比 RSK 法更加接近参考线 reference line(代表差值平均数为 0 的线),且 XU 法的 95% 一致性界限的区间范围更小。
通过实验还发现,对 1 例猪肝组织消融停止时刻(t=60 s 处)的 1 帧超声背散射信号,利用 RSK 法估算得到参数成像和 PAX 成像所耗时间为 81 s,利用 XU 法估算得到参数成像和 PAX 成像所耗时间为 17 s。
3 讨论
本文将超声背散射零差 K 成像用于监测微波消融,并对比了 2 种最新的零差 K 模型参数估算方法,即 RSK 法和 XU 法。通过离体猪肝微波消融实验验证了超声背散射零差 K 成像监测微波消融的可行性。
目前用于热消融定量监测的超声技术主要分为超声弹性成像和定量超声两大类[10,21,23-29]。超声弹性成像是对组织施加一个外部或内部的激励,然后利用超声波对组织产生的响应进行跟踪,估算出组织产生的应变、剪切波传播速度或弹性模量[24-26]。定量超声技术从超声背散射信号中估算出组织的声学参数,主要包括超声背散射统计成像[9-10,19]、散射子平均间距成像[27]、超声衰减系数成像[28]、超声回波去相关成像[23]以及基于卷积神经网络结构的超声成像等[29]。超声背散射统计成像中的超声 Nakagami 成像已用于消融术中监测,但 Nakagami 模型 m 参数在散射子数密度较高时存在饱和现象[8]。零差 K 模型能克服这一不足,但因其算法复杂,目前尚未应用于微波消融术中监测。
为避免无回声区的噪声在超声零差 K 成像中形成的伪影,本文引入 NCA 对超声背散射信号进行预处理。NCA 已应用于超声 Nakagami 成像,可有效消除成像无回声区域的伪影,同时不牺牲 Nakagami 成像的帧率[22]。NCA 的原理[21]是添加的白噪声幅值很小,对超声背散射信号没有明显的影响,因此对同一个背散射信号添加 2 次白噪声,得到的 2 条含额外噪声的超声背散射信号,其对应散射介质区域的波形几乎相同,它们之间的相关系数接近 1。相比之下,由于无回声区自身的噪声信号幅度较小,而每次添加的白噪声都具有随机性且其对无回声区信号作用较为明显,因此在无回声区内,2 条含额外噪声的超声背散射信号的相关系数会变小,即小于 1。通过设置相关阈值,NCA 可以自动区分无回声区域的白噪声和散射介质中的背散射信号。
将图 5、图 6 与图 4 对比发现,在监测组织微波消融时,相比于传统 B 超成像,超声零差 K 成像可以更好地监测消融过程中组织凝固区的变化,但当消融停止后,超声零差 K 成像不能用于消融结束的凝固区检测评估。这是由于微波消融过程中产生的微气泡可视为有效散射子,当消融开始后,随着微气泡的产生,组织消融区内有效散射子个数增加,超声零差 K 成像消融区亮度增强,当消融停止后,随着微气泡的消散,组织消融区的有效散射子个数减少,超声零差 K 成像消融区亮度逐渐减弱。
本文通过 20 例猪肝组织微波消融的超声背散射信号,对比验证 RSK 法和 XU 法超声零差 K 成像监测凝固区的效果。结果显示 RSK 法的平均 AUROC 和平均监测精度均小于 XU 法。XU 法估算的凝固区面积值与实际猪肝组织凝固区 Bland-Altman 一致性更强。XU 法估算参数并成像所耗时间更短。由此我们可以得出结论,零差 K 成像可用于监测猪肝组织微波消融过程中凝固区的变化,XU 法优于 RSK 法。
本研究尚有不足之处。首先,由于算法是利用 Matlab 实现,所以比较耗时,未来需考虑 C/C++语言或其他加速算法,以利实时监测。其次,本研究只对离体组织的微波消融进行监测,且实验数据只有 20 例,未来还需加大实验数据并考虑通过在体实验验证算法的有效性。另外,本研究中的零差 K 成像仅对消融结束时刻的监测精度和凝固区面积的 Bland-Altman 一致性进行分析,原因是其他时刻的参考标准不易获得,未来可在不同的消融时刻停止消融并剖开猪肝进行定量分析,以评估本文方法在不同时刻定量监测凝固区的效果。
4 结论
超声背散射零差 K 成像可用于监测离体猪肝微波消融过程中凝固区的变化。RSK 法和 XU 法的 AUROC 平均值 ± 标准差分别为 0.77 ± 0.06、0.83 ± 0.08。基于 RSK 法和 XU 法的零差 K 成像监测凝固区的精度平均值 ± 标准差分别为(86 ± 10)% 和(90 ± 8)%。XU 法估算凝固区面积值与实际猪肝组织凝固区面积的 Bland-Altman 一致性强于 RSK 法。另外,XU 法估算参数并成像所耗时间也少于 RSK 法。因此,在监测肝组织微波消融凝固区时,XU 法优于 RSK 法。
利益冲突声明:本文全体作者均声明不存在利益冲突。