基于可变滤波器长度的超声弹性成像方法技术

本发明专利技术提供了一种基于可变滤波器长度的超声弹性成像方法，首先根据变形前后的两帧RF信号，利用数值优化算法，计算出位移；然后通过传统固定滤波器长度的方法计算出初应变；再对该应变场进一步差分运算求出应变梯度，并对此线性变换得到可变的SG差分滤波器长度分布；最后，结合前面计算出的位移值与滤波器长度分布计算得到最终应变场。本发明专利技术获取的应变图中的运动边界更清晰，计算域内CNR值更高，具有更优的目标检测能力。

A method of ultrasonic elastography based on variable filter length

The present invention provides an ultrasound elastography method based on variable length filter, according to the first two frames RF signals before and after deformation, numerical optimization, calculation of displacement; and then through the traditional method of fixed length of filter to calculate the initial strain; the strain field further differential operation to extract the strain gradient, and this linear transform variable SG differential filter length distribution; finally, the displacement with the combination of the calculated value and the filter length distribution calculated ultimate strain field. The motion boundary obtained in the strain graph of the invention is clearer, and the CNR value in the calculation domain is higher, and the target detection capability is better.

【技术实现步骤摘要】
基于可变滤波器长度的超声弹性成像方法
本专利技术涉及超声弹性成像
，具体涉及一种基于可变滤波器长度的超声弹性成像方法。
技术介绍
生物组织的弹性特性，与组织的病理学状态紧密相关，具有重要的临床应用。超声弹性成像能无创定性定量生物组织的弹性，弥补了传统医学成像模态的不足，具有广阔的应用前景。应变计算是超声弹性成像技术的关键步骤，计算准确度将直接影响成像质量。从数学理论的角度，位移到应变的计算相当于一个数值微分的过程，即应变是位移的梯度。数值微分操作被当作一个数字差分滤波器。然而，数值微分操作是一个不稳定和危险的操作，因为它对信号的噪声或误差非常敏感，具有明显的放大作用，即很小的噪声或误差就会引起数值微分结果的很大误差。Savitzky-Golay差分滤波器(SG差分滤波器)的基本思想是移动、滑动、局部、分段、逐点多项式拟合，具有计算过程清晰简洁，滤波器系数容易获取，滤波器长度可以设置为任意值等诸多优点，能很好的应用于弹性成像中的应变计算。其中，SG差分滤波器长度是影响应变计算结果的一个关键参数。如果滤波器长度设置过短，则会增加应变计算错误率。反之，则会降低应变弹性图的分辨率，即会降低目标检测能力。因此，设置一个最优的滤波器长度是在超声弹性成像中准确计算应变的一个关键因素。传统利用SG差分滤波器计算应变，对滤波器长度都是设置为固定值，这在力学线性环境下，能取得较为满意的结果。但实际上生物组织是力学非线性，特别是对于含有异物(如肿瘤)的生物组织，周围组织与异物之间弹性差值更明显，射频数据采集环境复杂，获取的被测区域射频数据在各处的噪声含量都不一样。此时，采用固定的滤波器长度，成像质量会不佳。
技术实现思路
本申请通过提供一种基于可变滤波器长度的超声弹性成像方法，以解决采用固定滤波器长度而造成的成像质量不佳的技术问题。为解决上述技术问题，本申请采用以下技术方案予以实现：一种基于可变滤波器长度的超声弹性成像方法，包括如下步骤：S1：根据变形前后两帧超声射频信号，利用数值优化算法计算出位移场；S2：利用滤波器长度为固定值的SG差分滤波器对位移场进行差分运算，得出初应变场；S3：对初应变场进行差分运算求出含有异物边界信息的应变梯度场；S4：通过映射函数对应变梯度场的绝对值进行线性变换，得出与异物边界信息相关的滤波器长度分布场；S5：按照滤波器长度分布场，对步骤S1中求出的位移场进行SG差分滤波器差分运算，得出最终的应变场。进一步地，步骤S1中的数值优化算法具体为：假设图像窗口在运动前后满足仿射变换关系，运动前窗口中的任意一点P(x,y)对应于运动后窗口中点Q(x′,y′)，两者的坐标满足：式中，向量为窗口的移动和变形，u，v为刚体位移，为窗口的应变，即窗口的一阶变形，将目标函数Ob(Vrs)最小化即可同时得到向量Vrs中的6个参数值，其中，目标函数进一步地，步骤S2中采用1阶或2阶多项式拟合的SG差分滤波器进行应变计算，假设SG差分滤波器的长度为2M+1，则1阶或2阶多项式拟合的SG差分滤波器的输入输出关系为：进一步地，步骤S4中的映射函数为：FLvar＝round{ζ[1-(1-σ)Φ]}，FLvar为滤波器长度，round()为四舍五入函数，式中，ζ为一常数，σ为调整因子，Φ为应变梯度归一化函数，式中，MB为应变梯度场，MBmax为应变梯度场的最大值，MBmin为应变梯度场的最小值。可选地，调整因子σ取值为0.05～0.95。与现有技术相比，本申请提供的技术方案，具有的技术效果或优点是：本专利技术获取的应变图中的运动边界更清晰，计算域内CNR值更高，具有更优的目标检测能力。附图说明图1(a)为弹性图轴线上的位移曲线图；图1(b)为不同滤波器长度下的SG差分滤波器对应图1(a)的应变曲线图；图2(a)为物理模型中A-line的位置示意图；图2(b)为对应A-line的位移曲线图；图2(c)为对应A-line的理想应变曲线图；图2(d)为对应A-line的理想应变梯度曲线图；图3为本专利技术的流程图；图4(a)为单异物体模的侧向轴线上的应变曲线图；图4(b)为单异物体模的轴向轴线上的应变曲线图；图4(c)为单异物体模的CNR图；图5(a)为离体组织的侧向轴线上的应变曲线图；图5(b)为离体组织的轴向轴线上的应变曲线图；图5(c)为离体组织的CNR图。具体实施方式本申请实施例通过提供一种基于可变滤波器长度的超声弹性成像方法，以解决采用固定滤波器长度而造成的成像质量不佳的技术问题。为了更好的理解上述技术方案，下面将结合说明书附图以及具体的实施方式，对上述技术方案进行详细的说明。实施例噪声或误差一般表现为高频成分，理论分析与实验研究表明，采用低通数字差分器SG差分滤波器进行应变计算，能够取得较好的效果。SG差分滤波器的属性与多项式的阶次以及滤波器长度相关。滤波器长度越长，多项式的阶次越低，滤波器的噪声放大因子(滤波器系数的平方和)越小。本实施例采用1或2阶多项式拟合的SG差分滤波器进行应变计算研究。假设滤波器的长度为2M+1，1或2阶多项式拟合的SG差分滤波器的输入输出关系可以表示为：采用滤波器长度较长的SG差分滤波器能有效的去除噪声，但同时降低了目标检测能力。假设一弹性图轴线上的位移曲线如图1(a)所示，利用不同滤波器长度下的SG差分滤波器对该曲线进行应变计算，结果如图1(b)所示。从图1(b)可看出，随着滤波器长度的增加，应变曲线的运动边界(即背景组织与异物之间的边界)就越平滑，即目标检测能力就越低。从以上分析可知，SG差分滤波器长度是影响弹性成像质量的一个关键因素。不同于以往通过设置固定滤波器长度的方法计算应变，本专利技术通过获取的应变梯度作为反馈参数，并利用一个映射函数得到一个可变的滤波器长度分布场，再据此进行应变计算。所谓应变梯度，即是对位移的二次差分结果。定义超声数据图像轴向上的一条轴线为“A-line”，则A-line对应的位移、应变、应变梯度的演化过程如图2所示。位移、应变、应变梯度都包含有异物的位置信息，但应变梯度更能量化和定位异物的边界。基于上述理论，本专利技术提出了一种基于可变滤波器长度的超声弹性成像方法，如图3所示，包括如下步骤：S1：根据变形前后两帧超声射频信号，利用数值优化算法计算出位移场；数值优化算法具体为：假设图像窗口在运动前后满足仿射变换关系，运动前窗口中的任意一点P(x,y)对应于运动后窗口中点Q(x′,y′)，两者的坐标满足：式中，向量为窗口的移动和变形，u，v为刚体位移，为窗口的应变，即窗口的一阶变形，将目标函数Ob(Vrs)最小化即可同时得到向量Vrs中的6个参数值，其中，目标函数S2：利用滤波器长度为固定值的SG差分滤波器对位移场进行差分运算，得出初应变场；即采用1阶或2阶多项式拟合的SG差分滤波器进行应变计算，假设SG差分滤波器的长度为2M+1，则1阶或2阶多项式拟合的SG差分滤波器的输入输出关系为：S3：对初应变场进行差分运算求出含有异物边界信息的应变梯度场；S4：通过映射函数对应变梯度场的绝对值进行线性变换，得出与异物边界信息相关的滤波器长度分布场；映射函数为：FLvar＝round{ζ[1-(1-σ)Φ]}，FLvar为滤波器长度，round()为四舍五入函数，式中，ζ为一本文档来自技高网...

【技术保护点】
一种基于可变滤波器长度的超声弹性成像方法，其特征在于，包括如下步骤：S1：根据变形前后两帧超声射频信号，利用数值优化算法计算出位移场；S2：利用滤波器长度为固定值的SG差分滤波器对位移场进行差分运算，得出初应变场；S3：对初应变场进行差分运算求出含有异物边界信息的应变梯度场；S4：通过映射函数对应变梯度场的绝对值进行线性变换，得出与异物边界信息相关的滤波器长度分布场；S5：按照滤波器长度分布场，对步骤S1中求出的位移场进行SG差分滤波器差分运算，得出最终的应变场。

【技术特征摘要】
1.一种基于可变滤波器长度的超声弹性成像方法，其特征在于，包括如下步骤：S1：根据变形前后两帧超声射频信号，利用数值优化算法计算出位移场；S2：利用滤波器长度为固定值的SG差分滤波器对位移场进行差分运算，得出初应变场；S3：对初应变场进行差分运算求出含有异物边界信息的应变梯度场；S4：通过映射函数对应变梯度场的绝对值进行线性变换，得出与异物边界信息相关的滤波器长度分布场；S5：按照滤波器长度分布场，对步骤S1中求出的位移场进行SG差分滤波器差分运算，得出最终的应变场。2.根据权利要求1所述的基于可变滤波器长度的超声弹性成像方法，其特征在于，步骤S1中的数值优化算法具体为：假设图像窗口在运动前后满足仿射变换关系，运动前窗口中的任意一点P(x,y)对应于运动后窗口中点Q(x′,y′)，两者的坐标满足：式中，向量为窗口的移动和变形，u，...

【专利技术属性】
技术研发人员：朱新建，吴若愚，吴宝明，何庆华，
申请(专利权)人：中国人民解放军第三军医大学第三附属医院，
类型：发明
国别省市：重庆,50