基于超声rf背散射信号时频分析的弹性参数成像方法
【专利摘要】本发明公开了一种基于超声RF背散射信号时频分析的弹性参数成像方法，包括以下步骤：1、使用宽频超声线阵探头扫描活体肝组织，获取超声回波RF背散射信号；2、解调超声回波RF背散射信号并显示B型图；3、选择适当大小的ROI感兴趣区，大小设置为MXW；4、计算每一条扫描线的三层小波包分解，获取各个子频带能量谱序列；5、采用基于最小二乘法的四阶多项式拟合子频带能量谱序列，获取每条扫描线的弹性参量；6、计算ROI内所有扫描线的弹性参量的均值；7、将获得的均值映射为成像灰度，叠加在原B型图相应ROI位置上，实现ROI的弹性成像。具有能更好地反映肝组织微结构的特征等优点。
【专利说明】
基于超声RF背散射信号时频分析的弹性参数成像方法
技术领域
[0001] 本发明涉及一种超声医学技术，特别涉及一种基于超声RF背散射信号时频分析的 弹性参数成像方法。
【背景技术】
[0002] 肝组织的弹性或硬度改变与肝纤维化和肝硬化密切相关。而超声弹性成像技术是 观察生物组织软硬度的低成本、无损和简便的影像检查手段。
[0003] 目前超声弹性成像方法主要分为两大类，静态弹性成像技术和动态弹性成像技 术，主要采用组织应变率和弹性模量作为成像参数，而且已经应用于医学临床。但是，它们 仍需要改进。主要原因是：（1)静态弹性成像采用人工方式去模拟匀速加压，不可避免地会 受到人为因素的影响，产生的应变和位移可以随着压力大小，频率快慢的不同而变化，从而 导致生成的图像也会随之变化，并且该方法只适合于浅表器官的检查，即使是利用呼吸引 起位移的准静态弹性成像，也会因为假定了边界是无限的理想条件而带来测量误差，因此 静态弹性成像是一种定性影像检查手段，临床应用受限。（2)动态弹性成像首先使用了较高 功率的超声波激励组织，使之在超声辐射力作用下产生剪切波，通过采用超高速平面波捕 捉在剪切波传播中组织的位移，估计剪切波速度，在此计算中是假定了成像组织满足密度 一致性并且各向同性的条件，但实际上生物的组织密度可能不一致，而且大多是各向异性， 从而计算出的弹性模量会存在不同程度的误差；同时，高功率的激励超声波有可能对人体 造成一定程度的损害，目前在医学临床中颇受争议；另一方面该技术使用了超高速成像技 术，因此对超声探头和系统的要求较高，从而增加了探头和相关处理部分的制造成本。 [0004]近年来，有国外学者提出利用超声RF背散射信号进行时频分析，利用声谱提取出 与组织弹性密切相关的特征参量，并针对不同硬度的体模（杨氏模量分别为9.18kPa， 39.90kPa，129.27kPa)进行了实验，实验结果证明，采用超声RF背散射信号时频分析获取的 特征参量可以有效区分不同硬度的体模，区分度可以达到81.28%，但是并没有利用这些弹 性参量进行超声成像的研究。
【发明内容】
[0005] 本发明的目的在于克服现有技术的缺点与不足，提供一种基于超声RF背散射信号 时频分析的弹性参数成像方法，该弹性参数成像方法采用了常见的B超检查中获取的背散 射回波RF信号，利用基于小波包分解的时频分析方法，提取与肝组织软硬度相关的组织特 征，以此特征作为成像参数映射为成像灰度，对组织弹性进行成像。
[0006] 本发明的目的通过下述技术方案实现:一种基于超声RF背散射信号时频分析的弹 性参数成像方法，包括以下步骤：
[0007] (1)使用宽频超声线阵探头扫描活体肝组织，获取超声回波RF背散射信号；
[0008] (2)解调超声回波RF背散射信号并显示B型图；
[0009] (3)选择适当大小的R0I感兴趣区，大小可以设置为MXW，即Μ条扫描线，每条扫描线 有W个采样点；
[0010] (4)取R0I内的第m(m = 0,初始状态)条声束进行基于Mallat算法的三层小波包分 解，小波基选取Daubechiesie，提取八个子频带的能量谱序列。小波包分解中，分解层数的 增加可以让信号的高频成分和低频成分的分解都达到很精细的程度，经过N层小波包分解 之后，原信号的能量被分解到2 N个正交频带上，在这里N=3,所以产生八个子频带分解系数 序列，子频带分解系数序列的集合记作{SN,j(k)}。其中，k = 0,1，2,…(W/2n)_1，j = 0,1，2， 3,S^(k)为扫描线序列经过小波包分解到第N层的第j个子频带中的第k个分解系 数，每个子频带共有f个分解系数，j为小波包分解到第N层的第j个子频带；
[0011] (5)上一步获取的SN,j(k)所对应的子频带能量谱序列集合记作{E N,j(k)}，EN,j(k) 的计算公式如下：
[0012]
[0013]其中，En, j(k)表示分解到第N层的第j个子频带中的第k个分解系数的能量值。取出 如下序列
|对上述序列点进行基于最小二乘法的四阶 多项式拟合，公式如下：
[0014] P(x)=A+Bx+Cx2+Dx3+Ex4，
[0015] 取出弹性参量A;
[0016] (6)从k = 0开始，令k = k+l，重复步骤(5)，得到弹性参量A的序列
[0017]
直到k=(W/2N)_l终止。将弹性参量序列A求取一阶均值，于 是得到这条扫描线的弹性参量:? ,公式如下：
[0018]
[0019] (7)第!11条扫描线的弹性参量，记作1，其中，111 = 0，1，2广，-1。重复步骤(4)(5) (6)，求得每一条扫描线的弹性参量，因此Μ条扫描线的弹性参量组成弹性参量序列 ，计算R0I内所有扫描线的弹性参量均值ΕΤ，计算公式如下：
[0020]
[0021] ET即为所求的整个R0I区域的弹性参量均值；
[0022] (8)将整个R0I区域的弹性参量均值ET，按照一定比例进行调整，映射为成像灰度， 计算公式如下：
[0023]
[0024] 将瓦作为该R0I的象素值，叠加在原B型图相应位置上，实现R0I的弹性成像。其中 α为比例参数，根据具体反的值进行调整，这里并没有统一的标准，在实施例中根据具体情 况而定。
[0025] 本发明首先使用宽频超声线阵探头扫描活体人体肝组织，用同一探头接收超声回 波射频RF背散射信号；解调并且显示其Β型图，选取感兴趣区R0I;对R0I内的每一条扫描线 做三层小波包分解，提取各个子频带能量谱序列;采用四阶多项式拟合能量谱序列，获取每 条扫描线的弹性参量;将R0I内所有扫描线的弹性参量平均，均值映射为成像灰度，象素值 和参数值成比例关系，叠加在原Β型图的相应位置，实现R0I的弹性新参数成像。本发明首次 提出基于超声RF背散射信号时频分析的肝组织弹性新参数成像方法，能更好地反映肝组织 微结构的特征。
[0026] 本发明相对于现有技术具有如下的优点及效果：
[0027] 1、本发明首次提出基于超声RF背散射信号时频分析的肝组织弹性新参数成像方 法，为临床肝硬化诊断以及监测提供了参考。
[0028] 2、本发明首次提出采用超声RF背散射信号的时频特性，提取反映肝组织弹性的参 量，并实现反映弹性的新参数成像，无须增加额外成本，无损，简便，低成本。
【附图说明】
[0029]图1是本发明方法的流程图。
[0030] 图2是肝组织Β型图。
[0031]图3是小波包分解的示意图。
[0032]图4是弹性参量获取原理示意图。
[0033]图5a是正常肝组织感兴趣区R0I的超声弹性新参数成像图。
[0034]图5b是肝硬化肝组织感兴趣区R0I的超声弹性新参数成像图。
【具体实施方式】
[0035]下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限 于此。
[0036] 实施例
[0037] 如图1所示，本实施例基于超声RF背散射信号时频分析的肝组织弹性参量成像方 法详述如下。
[0038] 本实施例使用了加拿大Ultrasonix公司生产的Sonix TOUCH及中心频率为6·6ΜΗζ 的宽频线阵超声探头扫描肝脏包膜下的肝组织区域，获取超声回波RF背散射信号。
[0039]我们对获取的超声回波RF背散射信号进行希尔伯特解调并显示其超声B型图，如 图2所示。
[0040] 在超声B型图上选取大小为20X256的R0I，得到了R0I区域内的20条扫描线，每条扫 描线的采样点为256,我们对每一条扫描线进行基于Mallat算法的三层小波包分解，如图3 所示，小波基选取Daubechiesl6，得到八个子频带的分解系数序列，从而计算出子频带能量 谱序列。
[0041] 利用基于最小二乘法的四阶多项式拟合我们获取的子频带能量谱序列，如图4所 示，提取拟合系数A，计算由A组成的序列的均值，得到该扫描线的弹性参量整个R0I内所 有扫描线的弹性参量组成弹性参量序列不...?，. . ，其中m表示为第m条扫描线，整 个ROI共有20条扫描线。
[0042] 计算R0I内弹性参量序列的均值，计算如下：
[0043]
[0044]将整个R0I区域的弹性参量均值ET，按照一定比例进行调整，映射为成像灰度。在 我们的实验研究中，正常肝样本49例，肝硬化样本41例，正常肝样本ET的变化范围在 450.7755~559.7487，肝硬化样本ET的变化范围在589.6099~680.3542。我们选择比例参 数α = 1200计算R0I图像灰度五f，可以实现在原B超图上叠加弹性影像，如图5a所示，可见新 参数成像灰阶反映了组织弹性，可以区分正常和肝硬化样本，如图5b所示，R0I图像灰度元 的计算如下：
[0045]
[0046] 上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的 限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化， 均应为等效的置换方式，都包含在本发明的保护范围之内。
【主权项】
1. 基于超声RF背散射信号时频分析的弹性参数成像方法，其特征在于，包括以下步骤： (1) 使用宽频超声线阵探头扫描活体肝组织，获取超声回波RF背散射信号； (2) 解调超声回波RF背散射信号并显示B型图； (3) 选择适当大小的ROI感兴趣区，大小设置为MXW，即M为扫描线数，W为每条扫描线上 的采样点数； (4) 计算每一条扫描线的三层小波包分解，获取各个子频带能量谱序列； (5) 采用基于最小二乘法的四阶多项式拟合子频带能量谱序列，获取每条扫描线的弹 性参量； (6) 计算ROI内所有扫描线的弹性参量的均值； (7) 将获得的均值映射为成像灰度，叠加在原B型图相应ROI位置上，实现ROI的弹性成 像。2. 根据权利要求1所述的基于超声RF背散射信号时频分析的弹性参数成像方法，其特 征在于，所述步骤(4)中，获取子频带能量谱序列的方法如下:设感兴趣区ROI大小为MXW， 拥有M条扫描线，对每一条扫描线进行基于Mai lat算法的小波包分解，小波基选择 Daubechiesl6； 小波包分解中，分解层数的增加让信号的高频成分和低频成分的分解都达到很精细的 程度;扫描线序列经过N层小波包分解之后，N表示小波包分解的层数，可以得到信号在任一 频段的特征信息，即原信号的能量被分解到2 N个正交频带上，其中，N = 3,所以产生八个子 频带分解系数序列；子频带分解系数序列的集合记作{S^(k)}，对应的子频带能量谱序列 集合记作{EN,j(k)};能量序列的连续计算公式如下：式中：S^(t)为信号经过小波包分解至第N层的第j个子频带中时刻t的幅值，j为小波 包分解到第N层的第j个子频带;ENjt)表示信号分解到第N层的第j个子频带中时亥Ijt的能 量值； 离散形式表现为：式中，S^(k)为扫描线序列经过小波包分解到第N层的第j个子频带中的第k个分解系 数，每个子频带共有个分解系数，j为小波包分解到第N层的第j个子频带;EN」(k)表示分 解到第N层的第j个子频带中的第k个分解系数的能量值;其中，k = 0，1，2，…(W/2n)-1，j = 0，1，2,3，."2n-1。3. 根据权利要求1所述的基于超声RF背散射信号时频分析的弹性参数成像方法，其特 征在于，所述步骤(5)中，获取每条扫描线的弹性参量的方法如下：获得了一条扫描线的子 频带能量谱序列集合{ENu(k)}后，取出如下序列[心⑷,...心:、，(幻1,对上 述序列点进行基于最小二乘法的四阶多项式拟合，公式如下： P(x) =A+Bx+Cx2+Dx3+Ex4? 取出弹性参量A，从k = 0开始，令k = k+l，重复上述步骤，得到弹性参量A的序列 |為，4,4^1^^1|，直到1^=(1/24-1终止;将4序列求取一阶均值，于是得到这条扫描线 的弹性参量d,公式如下：4. 根据权利要求1所述的基于超声RF背散射信号时频分析的弹性参数成像方法，其特 征在于，所述步骤(6)中，计算ROI内所有扫描线的弹性参量均值的方法如下：获得了第m条 扫描线的弹性参量，记作了，其中，m = 0,1，2, 一M-l ;M条扫描线组成弹性参量序列 m ，计算ROI内所有扫描线的弹性参量均值ET，计算公式如下：ET即为所求的整个R0I区域的弹性参量均值。5. 根据权利要求1所述的基于超声RF背散射信号时频分析的弹性参数成像方法，其特 征在于，所述步骤(7)中，实现超声弹性新参数成像的方法如下:根据权利要求1中所述的超 声RF背散射信号时频分析的肝组织新参数成像方法，步骤(6)获取了整个R0I区域的弹性参 量均值ET，按照一定比例进行调整，映射为成像灰度，计算公式如下：将巧作为该R0I的象素值，叠加在原B型图相应位置上，实现R0I的弹性成像，其中，a根 据具体有效的弹性参数进行调整。
【文档编号】A61B8/08GK106037799SQ201610471908
【公开日】2016年10月26日
【申请日】2016年6月22日
【发明人】林春漪, 潘浩
【申请人】华南理工大学