弹性成像方法及系统.pdf

《弹性成像方法及系统.pdf》由会员分享，可在线阅读，更多相关《弹性成像方法及系统.pdf（11页完整版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利 (10)授权公告号 (45)授权公告日 (21)申请号 201410366344.4 (22)申请日 2014.07.29 (65)同一申请的已公布的文献号 申请公布号 CN 105326529 A (43)申请公布日 2016.02.17 (73)专利权人 深圳迈瑞生物医疗电子股份有限 公司 地址 518000 广东省深圳市南山区高新技 术产业园区科技南十二路迈瑞大厦 (72)发明人 袁宇辰 樊睿 (74)专利代理机构 广州三环专利商标代理有限 公司 44202 代理人 郝传鑫 熊永强 (51)Int.Cl. A61B 8/08(200。

2、6.01) (56)对比文件 CN 102824193 A,2012.12.19, CN 101569543 A,2009.11.04, CN 103211625 A,2013.07.24, US 2010305438 A1,2010.12.02, CN 1586408 A,2005.03.02, US 2010106018 A1,2010.04.29, 审查员 沈研研 (54)发明名称 弹性成像方法及系统 (57)摘要 本发明提供了一种弹性成像方法， 包括以下 步骤， 获取压缩前数据与压缩后数据， 所述压缩 前数据与压缩后数据具有排布为数据矩阵的多 个数据点； 将压缩前数据与压缩后数据进行网。

3、格 划分； 以网格为单位计算压缩前数据与压缩后数 据相互匹配的网格之间的块匹配位移， 获取粗位 移结果； 将粗位移结果拟合为以基准数据的数据 点为单位的拟合位移结果； 将拟合位移结果叠加 于基准数据中的数据点， 取叠加后获得的数据点 为校正拟合数据点， 并获取校正拟合数据点的校 正拟合数据； 获取应变场并生成弹性图像。 本发 明通过局部形变获取校正拟合数据， 从而减小两 帧数据之间的整体相对位移， 提高两帧数据相关 性。 本发明还提供一种弹性成像系统。 权利要求书2页 说明书6页 附图2页 CN 105326529 B 2017.09.26 CN 105326529 B 1.一种弹性成像方法，。

4、 其特征在于， 包括以下步骤， 获取压缩前数据与压缩后数据， 所述压缩前数据与压缩后数据具有排布为数据矩阵的 多个数据点； 将压缩前数据与压缩后数据进行网格划分； 以网格为单位计算压缩前数据与压缩后数据相互匹配的网格之间的块匹配位移， 获取 粗位移结果； 设压缩前数据与压缩后数据二者中的一个为基准数据， 将粗位移结果拟合为以所述基 准数据的数据点为单位的拟合位移结果； 将拟合位移结果叠加于基准数据中的数据点， 取叠加后获得的数据点为校正拟合数据 点， 并获取所述校正拟合数据点的校正拟合数据； 获取所述校正拟合数据与所述压缩前数据、 所述校正拟合数据与所述压缩后数据之间 的位移场， 并根据所述校。

5、正拟合数据与所述压缩前数据、 所述校正拟合数据与所述压缩后 数据之间的位移场获取所述压缩前数据与所述压缩后数据之间的位移场与应变场； 根据所述应变场生成弹性图像。 2.如权利要求1所述的弹性成像方法， 其特征在于， 所述压缩前数据与压缩后数据为I/ Q基带信号帧数据。 3.如权利要求2所述的弹性成像方法， 其特征在于， 当将拟合位移结果叠加于基准数据 中的数据点， 取叠加后获得的数据点为校正拟合数据点， 并获取校正拟合数据点的校正拟 合数据后， 所述数据处理方法还包括以下步骤， 对校正拟合数据进行相位偏差校正。 4.如权利要求1至3中任一项所述的弹性成像方法， 其特征在于， 当以网格为单位计算。

6、 压缩前数据与压缩后数据相互匹配的网格之间的块匹配位移， 获取粗位移结果时， 所述块 匹配位移的匹配标准采用最小化方差和、 均方误差、 绝对平均误差或最大化互相关函数以 及弹性极限差分分类中的任意一种。 5.如权利要求1至3中任一项所述的弹性成像方法， 其特征在于， 当将拟合位移结果叠 加于基准数据中的数据点， 取叠加后获得的数据点为校正拟合数据点， 并获取校正拟合数 据点的校正拟合数据时， 进一步包括以下步骤， 将拟合位移结果叠加于基准数据中对应的数据点， 获取校正拟合数据点； 根据基准数据获取所述校正拟合数据点的数值为整数的校正拟合数据。 6.如权利要求5所述的弹性成像方法， 其特征在于，。

7、 当根据基准数据获取校正拟合数据 点的数值为整数的校正拟合数据时， 利用插值方法获取数值为整数的校正拟合数据。 7.如权利要求6所述的弹性成像方法， 其特征在于， 所述插值方法采用线性插值、 多项 式插值、 三次样条插值中的任意一种。 8.一种弹性成像系统， 其特征在于， 包括 采样模块， 用于获取压缩前数据与压缩后数据； 网格划分模块， 用于将压缩前数据与压缩后数据进行网格划分； 计算模块， 用于以网格为单位计算压缩前数据与压缩后数据相互匹配的网格之间的块 匹配位移， 获取粗位移结果； 拟合模块， 用于将粗位移结果拟合为以所述压缩前数据或压缩后数据的数据点为单位 权 利 要 求 书 1/2 。

8、页 2 CN 105326529 B 2 的拟合位移结果； 叠加模块， 用于将拟合位移结果叠加于压缩前数据或压缩后数据的数据点并获取校正 拟合数据； 成像模块， 根据所述校正拟合数据、 压缩前数据及压缩后数据获取压缩前数据与压缩 后数据之间的位移场及应变场， 生成弹性图像。 9.如权利要求8所述的弹性成像系统， 其特征在于， 所述弹性成像系统还设有校正模 块， 用于对校正拟合数据进行相位偏差校正。 10.如权利要求8所述的弹性成像系统， 其特征在于， 所述拟合模块设有： 拟合单元， 用于将拟合位移结果叠加于压缩前数据中对应的数据点， 获取校正拟合数 据点； 插值单元， 用于利用插值方法获取校正。

9、拟合数据。 权 利 要 求 书 2/2 页 3 CN 105326529 B 3 弹性成像方法及系统 技术领域 0001 本发明涉及一种弹性成像方法及系统。 背景技术 0002 医用超声弹性成像主要指以显示组织弹性差异为目的的一系列成像和信号处理 技术。 如图1所示， 其通常过程为： 通过手持超声探头对目标组织施加压力(pressure)， 探头 以系统预先设定好的扫描规则进行超声波发射并接收回波信息， 经过波束合成环节后输出 射频(radio frequency， RF)信号， 通过计算出压缩前后所采集的RF(Radio Frequency， 射 频)数据中前后两帧数据对应位置发生的位移(d。

10、isplacement)， 进而得到该组织区域各点 的应变值(strain)， 最后显示输出弹性图像。 所采集的RF数据中前后两帧数据的相对位移 越小， 则两帧的相关性越高， 实际应用中， 由于噪声、 所施加压力的不均匀性和所施加压力 幅度变化等因素的影响， 前后两帧数据之间的相关性难以保证， 导致压缩前后数据间的应 变计算出现误差。 这种误差会导致弹性应变图像局部甚至整帧错误， 影响诊断效果。 发明内容 0003 提供一种可提高数据相关性的弹性成像方法。 0004 一种弹性成像方法， 包括以下步骤， 0005 获取压缩前数据与压缩后数据， 所述压缩前数据与压缩后数据具有排布为数据矩 阵的多个。

11、数据点； 0006 将压缩前数据与压缩后数据进行网格划分； 0007 以网格为单位计算压缩前数据与压缩后数据相互匹配的网格之间的块匹配位移， 获取粗位移结果； 0008 设压缩前数据与压缩后数据二者中的一个为基准数据， 将粗位移结果拟合为以所 述基准数据的数据点为单位的拟合位移结果； 0009 将拟合位移结果叠加于基准数据中的数据点， 取叠加后获得的数据点为校正拟合 数据点， 并获取校正拟合数据点的校正拟合数据； 0010 获取校正拟合数据与压缩前数据、 校正拟合数据与压缩后数据之间的位移场， 并 根据校正拟合数据与压缩前数据、 校正拟合数据与压缩后数据之间的位移场获取压缩前数 据与压缩后数据。

12、之间的位移场与应变场； 0011 根据所述应变场生成弹性图像。 0012 进一步的， 所述压缩前数据与压缩后数据为I/Q基带信号帧数据。 0013 进一步的， 当将拟合位移结果叠加于基准数据中的数据点， 取叠加后获得的数据 点为校正拟合数据点， 并获取校正拟合数据点的校正拟合数据后， 所述数据处理方法还包 括以下步骤， 0014 对校正拟合数据进行相位偏差校正。 0015 进一步的， 所述以网格为单位计算压缩前数据与压缩后数据相互匹配的网格之间 说 明 书 1/6 页 4 CN 105326529 B 4 的块匹配位移， 获取粗位移结果时， 所述块匹配位移的匹配标准采用最小化方差和、 均方误 。

13、差、 绝对平均误差或最大化互相关函数以及弹性极限差分分类中的任意一种。 0016 进一步的， 当将拟合位移结果叠加于基准数据中的数据点， 取叠加后获得的数据 点为校正拟合数据点， 并获取校正拟合数据点的校正拟合数据时， 进一步包括以下步骤， 0017 将拟合位移结果叠加于基准数据中对应的数据点， 获取校正拟合数据点； 0018 根据基准数据获取校正拟合数据点的数值为整数的校正拟合数据。 0019 进一步的， 当根据基准数据获取校正拟合数据点的数值为整数的校正拟合数据 时， 利用插值方法获取数值为整数的校正拟合数据。 0020 进一步的， 所述插值方法采用线性插值、 多项式插值、 三次样条插值中。

14、的任意一 种。 0021 本发明还提供一种弹性成像系统， 包括 0022 采样模块， 用于获取压缩前数据与压缩后数据； 0023 网格划分模块， 用于将压缩前数据与压缩后数据进行网格划分； 0024 计算模块， 用于以网格为单位计算压缩前数据与压缩后数据相互匹配的网格之间 的块匹配位移， 获取粗位移结果； 0025 拟合模块， 用于将粗位移结果拟合为以所述压缩前数据的数据点为单位的拟合位 移结果； 0026 叠加模块， 用于将拟合位移结果叠加于压缩前数据或压缩后数据的数据点并获取 校正拟合数据； 0027 成像模块， 根据所述校正拟合数据、 压缩前数据及压缩后数据获取压缩前数据与 压缩后数据之。

15、间的位移场及应变场， 生成弹性图像。 0028 进一步的， 所述用于弹性成像的数据处理系统还设有校正模块， 用于对校正拟合 数据进行相位偏差校正。 0029 进一步的， 所述拟合模块设有： 0030 拟合单元， 用于将拟合位移结果叠加于压缩前数据中对应的数据点， 获取校正拟 合数据点； 0031 插值单元， 用于利用插值方法获取校正拟合数据。 0032 本发明的弹性成像方法， 通过局部形变获取校正拟合数据， 通过校正拟合数据减 小两帧数据之间的整体相对位移， 提高两帧数据相关性。 由于校正拟合数据通过拟合进行 了局部变形， 其与压缩后数据之间相似性较高， 降低因不相关造成的应变计算误差， 便于。

16、弹 性图像的计算与获取。 附图说明 0033 为了更清楚地说明本发明实施例或现有技术中的技术方案， 下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍， 显而易见地， 下面描述中的附图仅仅是本 发明的一些实施例， 对于本领域普通技术人员来讲， 在不付出创造性劳动的前提下， 还可以 根据这些附图获得其他的附图。 0034 图1是现有技术中超声弹性成像的实施流程示意图； 0035 图2是本发明的较佳实施例的弹性成像方法的流程示意图； 说 明 书 2/6 页 5 CN 105326529 B 5 0036 图3是本发明的较佳实施例的弹性成像系统的构成示意图。 具体实施方式 0037 下面将。

17、结合本发明实施例中的附图， 对本发明实施例中的技术方案进行清楚、 完 整地描述， 显然， 所描述的实施例仅仅是本发明一部分实施例， 而不是全部的实施例。 基于 本发明中的实施例， 本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他 实施例， 都属于本发明保护的范围。 0038 请参阅图2， 本发明较佳实施方式提供一种弹性成像方法， 所述方法包括以下步 骤： 0039 步骤S101， 获取数据， 具体地， 获取一对I/Q基带信号帧数据， 在本实施例中， 在通 过超声波探头对待测组织施加压力并收发信号过程中获取该对I/Q基带信号帧数据， 该对 I/Q基带信号帧数据记为压缩前数据(U数据)。

18、与压缩后数据(C数据)。 在本实施实施例中， 所 述I/Q基带信号帧数据为超声探头输出的射频(radio frequency， RF)信号经过正交解调后 生成。 0040 每一帧I/Q基带信号帧数据， 都分为I、 Q两路数据。 所述压缩前数据与压缩后数据 可以是时间上连续两帧数据， 可以是有一定帧间隔的两帧数据， 帧间隔数量由系统预先设 定。 使用一定间隔的两帧数据， 可以有效调整用于计算的两帧数据间的位移量大小， 以便获 取较佳的应变图像。 0041 在本实施例中， 所述压缩前数据与压缩后数据包括多条数据线， 每条数据线包括 多个数据点。 所述压缩前数据与压缩后数据的多条数据线上及多条数据线。

19、上的多个数据 点。 所述压缩前数据的多个数据线上的多个数据点构成具有多行、 多列数据点的压缩前数 据矩阵； 所述压缩后数据的多个数据线上的多个数据点构成具有多行、 多列数据点的压缩 后数据矩阵。 所述多条数据线的延伸方向可作为纵坐标。 0042 步骤S102， 将压缩前数据与压缩后数据进行网格划分， 从而I/Q基带信号帧数据划 分为多个网格， 多个网格排布构成多个行及多个列， 且单个网格中具有若干个数据点。 0043 网格划分的目的是将压缩前数据与压缩后数据划分为便于操作的小块， 以便之后 的块匹配计算(Block Matching Algorithm， BMA)。 划分时的主要参数为每个网格。

20、的尺寸， 需要综合考虑所需的图像分辨率和计算速度。 网格越小， 分辨率越高， 但计算速度也越慢。 此外网格不能太小， 否则会降低块匹配效果， 造成匹配错误。 其具体划分方法与现有技术一 致， 在此不再赘述。 本实施例中， 步骤S102采用单个网格作为之后偏移量搜索的单位。 0044 步骤S103， 在压缩前数据中， 以网格为单位计算压缩前数据与压缩后数据相互匹 配的网格之间的块匹配位移， 从而获取粗位移结果。 在本步骤中， 可根据数据处理需要采用 各类匹配标准， 如最小化方差和、 均方误差、 绝对平均误差或最大化互相关函数以及弹性极 限差分分类等。 可以理解的是， 在步骤S103中， 在获取压。

21、缩前数据的粗位移结果时， 压缩前 数据除第一行的网格外的其他行网格需要叠加上一行的粗位移结果。 可以理解的是， 压缩 前数据与压缩后数据低位等价， 因此可以采用压缩后数据作为参考对压缩前数据进行粗位 移叠加， 也可采用压缩前数据为参考对压缩后数据进行粗位移叠加。 0045 可以理解的是， 在步骤S103中， 既可以对整帧数据进行粗位移计算， 也可以只取其 中若干列进行计算， 并将计算结果取平均值， 赋予该行各点， 作为该行内各点的粗位移值； 说 明 书 3/6 页 6 CN 105326529 B 6 0046 步骤S104， 设压缩前数据与压缩后数据二者中的一个为基准数据， 获取拟合位移 结。

22、果， 将粗位移结果拟合为以所述基准数据的I/Q基带信号帧数据的数据点为单位的拟合 位移结果。 在本实施例中， 获取压缩前数据以网格为单位的粗位移结果， 对压缩前数据的各 个数据点根据粗位移结果进行拟合， 从而获得以压缩前数据的各数据点为单位的拟合位移 结果。 0047 步骤S105， 将所有数据点的拟合位移结果叠加于基准数据对应的原始数据点， 取 叠加后获得的数据点为校正拟合数据点， 并获取校正拟合数据点的校正拟合数据。 本步骤 中， 通过将原始数据点的位置(或坐标)累加上拟合位移结果， 从而获取校正拟合数据点。 在 本实施例中， 校正拟合数据点由压缩前数据中的原始数据点的位置累加拟合位移结果。

23、而 成。 0048 所述步骤S105进一步包括： 0049 步骤S1051， 将所有数据点的拟合位移结果叠加于基准数据中对应的原始数据点， 获取校正拟合数据点。 在此步骤中， 由于一般来说， 粗位移结果的数值不是整数， 因此叠加 后校正拟合数据点也不是整数。 0050 步骤S1052， 根据基准数据获取数值为整数的校正拟合数据。 在此步骤中， 因此利 用插值的方法将叠加后的纵坐标各整数点的幅值重新算出， 并将各整数点的幅值写入校正 拟合数据点相应的位置， 从而获得完整的校正拟合数据。 所述插值方法可采用线性插值、 多 项式插值、 三次样条插值等各类适用的插值方法。 所述校正拟合数据经过粗位移估。

24、计并进 行局部形变， 可以减少与压缩后数据之间的整体相对位移， 提高两帧数据相关性， 降低因不 相关造成的应变计算误差。 0051 步骤S106， 获取校正拟合数据与压缩前数据、 校正拟合数据与压缩后数据之间的 位移场， 并根据校正拟合数据与压缩前数据、 校正拟合数据与压缩后数据之间的位移场获 取压缩前数据与压缩后数据之间的位移场与应变场。 0052 步骤S107， 根据所述应变场生成弹性图像。 0053 为获取校正拟合数据， 需要根据给定校正拟合数据的位移场对压缩前数据进行变 形操作。 采用射频(radio frequency， RF)信号直接进行信号的局部形变， 信号采样率较高， 计算量较。

25、大。 而本方法中的I/Q基带信号帧数据由于可以用带宽采样， 可以在较低的采样率 下进行， 因此能够节省计算量。 但是对I/Q基带信号帧数据直接进行信号变形会引入相位误 差， 需要进行相位偏差校正。 0054 相位误差的形成原因分析如下： 0055 假设原始信号点为RF0， 信号中心频率为0， 对应深度位置为t， 则正交解调后进行 位移叠加时， 结果为: 0056 0057 上式中D0为解调获得的I/Q基带信号帧数据， I0为I/Q基带信号帧数据的I路数据， Q0为I/Q基带信号帧数据的Q路数据； 表示希尔伯特(Hilbert)变换， RF0为射频信号。 0058 由于RF0为射频信号， 因此其。

26、频谱关于原点对称， 只保留RF0在正半轴的频谱； 而之 后乘以的相当于对希尔伯特变换后的频谱卷积 (-0)， 将频谱移到基带。 这里的I0、 Q0也是直接对原始信号点的I/Q基带信号帧数据进行位移叠加后的值， 即位移叠加只移动 位置， 不改变I/Q基带信号帧数据的幅值。 假设该点位置由t移动至t+t， 这里t即相当于 说 明 书 4/6 页 7 CN 105326529 B 7 对应点的位移值。 0059 相较于本发明中先对原始信号做正交解调生成I/Q基带信号帧数据D0， 而后进行 位移叠加； 下式中的D1,I1,Q1则对应的是先对原始信号做位移叠加， 再进行正交解调。 0060 D1为射频信。

27、号位移叠加后正交解调获取的I/Q基带信号帧数据， I1为I/Q基带信号 帧数据D1的I路数据， Q1为I/Q基带信号帧数据D1的Q路数据； D0与D1之间具有相位偏差， D1较 为准确。 0061 由原始的射频信号直接解调的t+t点的正交解调结果为 0062 0063 则有如下关系 0064 D1I1+iQ1D0ei t (1) 0065 或 0066 I1I0cos(t)-Q0sin(t) 0067 Q1Q0cos(t)+I0sin(t) 0068 上式可以理解为， 直接对I/Q基带信号帧数据进行位移叠加并进行变形处理， 会引 入和变形位移场相关的相位偏差。 因此需要如式(1)中所述， 对I。

28、/Q基带信号帧数据进行位 移叠加多获取的校正拟合数据进行相位偏差校正， 即对校正拟合数据乘上矫正相位因子ei t。 可以理解的是， 可根据实际情况选用适用的方法进行相位偏差校正。 0069 本发明还提供一种用于弹性成像的弹性成像系统10， 包括 0070 采样模块11， 用于获取压缩前数据与压缩后数据。 0071 网格划分模块12， 用于将压缩前数据与压缩后数据进行网格划分。 0072 计算模块13， 用于以网格为单位计算压缩前数据与压缩后数据相互匹配的网格之 间的块匹配位移， 获取粗位移结果。 0073 拟合模块15， 用于将粗位移结果拟合为以所述压缩前数据或压缩后数据的数据点 为单位的拟合。

29、位移结果； 进一步的， 所述拟合模块15设有拟合单元151及插值单元153。 拟合 单元151， 用于将拟合位移结果叠加于压缩前数据或压缩后数据中对应的数据点， 获取拟合 数据点； 插值单元153， 用于利用插值方法获取数值为整数的校正拟合数据。 0074 叠加模块17， 用于将拟合位移结果叠加于压缩前数据或压缩后数据中的数据点， 获取校正拟合数据。 0075 成像模块18， 根据所述校正拟合数据、 压缩前数据及压缩后数据获取压缩前数据 与压缩后数据之间的位移场及应变场， 生成弹性图像。 0076 进一步的， 所述用于弹性成像的数据处理系统10还设有校正模块19， 用于对校正 拟合数据进行相位。

30、偏差校正。 校正模块19可连接于所述叠加模块17与所述成像模块18之 间。 0077 在后续弹性图像的计算、 获取过程中， 压缩前数据与压缩后数据之间的位移场等 于压缩前数据至校正拟合数据的位移场加上校正拟合数据至压缩后数据的位移场。 因此本 发明的弹性成像方法， 通过局部形变获取校正拟合数据， 通过校正拟合数据减小两帧数据 之间的整体相对位移， 提高两帧数据相关性。 由于校正拟合数据通过拟合进行了局部变形， 其与压缩后数据之间相似性较高， 降低因不相关造成的应变计算误差， 便于后续弹性图像 的计算与获取。 同时， 该方法放宽了使用者进行超声波探头使对组织的压力大小的限制， 允 说 明 书 5。

31、/6 页 8 CN 105326529 B 8 许更大的施加压力， 降低了对使用者操作技巧的要求。 且其通过对I/Q基带信号帧数据进行 局部形变， 可通过软件实现， 成本较低， 整体流程设置较容易。 0078 以上所揭露的仅为本发明一种较佳实施例而已， 当然不能以此来限定本发明之权 利范围， 本领域普通技术人员可以理解实现上述实施例的全部或部分流程， 并依本发明权 利要求所作的等同变化， 仍属于发明所涵盖的范围。 说 明 书 6/6 页 9 CN 105326529 B 9 图1 图2 说 明 书 附 图 1/2 页 10 CN 105326529 B 10 图3 说 明 书 附 图 2/2 页 11 CN 105326529 B 11 。