位移估计方法、位移估计装置.pdf

1、(10)申请公布号 CN 102361597 A (43)申请公布日 2012.02.22 CN 102361597 A *CN102361597A* (21)申请号 201180001460.0 (22)申请日 2011.01.19 2010-010336 2010.01.20 JP A61B 8/08(2006.01) G01N 29/06(2006.01) (71)申请人 松下电器产业株式会社 地址 日本大阪府 (72)发明人 陈英俊 张国成 范树峰 (74)专利代理机构 永新专利商标代理有限公司 72002 代理人 王成坤 胡建新 (54) 发明名称 位移估计方法、 位移估计装置 (5

2、7) 摘要 本发明公开一种位移估计方法， 使用超声波 信号反复估计位移， 其中， 包含如下步骤 ： 扫描至 少 1 个超声波信号， 并发送给介质 ； 接收所扫描 的从所述介质反射的所述超声波信号 ； 计算窗尺 寸 ； 根据所计算的所述窗尺寸， 计算窗的边界 ； 使 用所计算的所述窗， 估计所述超声波信号的各个 深度的位移 ； 根据所估计的所述位移， 对所述超 声波信号进行卷绕 ； 以及使用所卷绕的所述超声 波信号， 引导该位移估计方法的收敛， 使得所述超 声波信号的相关值变大。 (30)优先权数据 (85)PCT申请进入国家阶段日 2011.09.20 (86)PCT申请的申请数据 PCT/J

3、P2011/000259 2011.01.19 (87)PCT申请的公布数据 WO2011/089898 JA 2011.07.28 (51)Int.Cl. (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 15 页 附图 13 页 CN 102361614 A1/2 页 2 1. 一种位移估计方法， 使用超声波信号反复估计位移， 其中， 该位移估计方法包含如下步骤 ： 扫描至少 1 个超声波信号， 并发送给介质 ； 接收所扫描的从所述介质反射的所述超声波信号 ； 计算窗尺寸 ； 根据所计算的所述窗尺寸， 计算窗的边界 ； 使用基于所计算的所述边界的所述窗

4、， 估计所述超声波信号的各个深度的位移 ； 根据所估计的所述位移， 对所述超声波信号进行卷绕 ； 以及 使用所卷绕的所述超声波信号， 引导该位移估计方法的收敛， 使得所述超声波信号的 相关值变大。 2. 根据权利要求 1 所述的位移估计方法， 其中， 改变所计算的所述窗尺寸， 使得连续的多个循环中使用的所述窗尺寸逐渐减小。 3. 根据权利要求 1 所述的位移估计方法， 其中， 该位移估计方法还包含如下步骤 ： 使连续的多个所述循环的所述窗尺寸的上限和下限分别按照逐渐减小的方式变化 ； 以 及 根据各深度的信号功率， 计算该深度的上限和下限间的窗尺寸， 作为各个深度的所述 窗尺寸。 4. 根据权

5、利要求 1 所述的位移估计方法， 其中， 计算所述窗尺寸， 使得所有窗中的信号能量相等。 5. 根据权利要求 1 所述的位移估计方法， 其中， 超声波信号的相关值越大， 计算出越小的所述窗尺寸。 6. 根据权利要求 1 所述的位移估计方法， 其中， 所述边界向该深度的两侧延伸， 使得基于该边界的所述窗成为以对应的深度为中心的 对称的窗。 7. 根据权利要求 1 所述的位移估计方法， 其中， 所述边界向对应的所述深度的两侧延伸， 使得对应深度的基于该边界的所述窗的两侧 的信号能量相互相等。 8. 根据权利要求 1 所述的位移估计方法， 其中， 该位移估计方法还包含如下步骤 ： 根据所计算的延迟值

6、使所述超声波信号内的各采样延迟， 由此卷绕所述超声波信号。 9. 根据权利要求 1 所述的位移估计方法， 其中， 该位移估计方法还包含如下步骤 ： 在各循环之后， 计算所卷绕的 RF 信号间的信号差分 ； 以及 决定所计算的所述信号差分为最小的循环。 10. 一种位移估计装置， 使用超声波信号反复估计位移， 其中， 该位移估计装置包含 ： 发送部， 扫描至少 1 个超声波信号， 并发送给介质 ； 接收部， 接收所扫描的从所述介质反射的所述超声波信号 ； 权 利 要 求 书 CN 102361597 A CN 102361614 A2/2 页 3 尺寸计算部， 计算窗尺寸 ； 边界计算部， 根据

7、所计算的所述窗尺寸， 计算窗的边界 ； 估计部， 使用基于所计算的所述边界的所述窗， 估计所述超声波信号的各个深度的位 移 ； 卷绕部， 根据所估计的所述位移， 对所述超声波信号进行卷绕 ； 以及 收敛控制部， 使用所卷绕的所述超声波信号， 引导该位移估计装置执行的位移估计方 法的收敛， 使得所述超声波信号的相关值变大。 权 利 要 求 书 CN 102361597 A CN 102361614 A1/15 页 4 位移估计方法、 位移估计装置 技术领域 0001 本发明导入使用超声波信号来估计位移的方法。 无论是最终结果还是用于进一步 处理的中间步骤， 在需要位移估计的应用例中都能够使用，

8、能够用于医疗用和工业用超声 波设备。 背景技术 0002 超声波装置通过对扫描对象介质收发高频的机械声波， 对介质进行非侵袭地诊 断。 0003 这种装置的变换器将超声波发送给扫描对象介质。 0004 超声波通过散射和反射， 与该介质的下层构造相互作用。 0005 另外， 下层构造例如是指比介质的表面更靠近介质的内部侧的构造等。下层构造 例如是人身体上的血管等。 0006 被散射和反射的超声波包含下层构造的有益信息， 由变换器接收。然后由超声波 装置进行处理， 例如作为图像信息提示给用户。 0007 图 8 是表示超声波 RF 信号的图。 0008 作 为 超 声 波 装 置 从 反 射 波

9、获 得 的 最 基 础 的 数 据 的 种 类 之 一， 列 举 了 RF(Radiofrequency Signal ： 射频信号 ) 信号 ( 图 8)。 0009 这是将接收波从模拟直接变换成数字的情况。 0010 在不同的应用例中， 能够从超声波 RF 信号中导出亮度模式 (B 模式 ) 图像、 多普勒 图像等其他种类的数据。 0011 多个应用例中的一个用于分析所扫描的介质的下层构造的动作。 0012 多普勒效应被用作估计构造的动作方向和强度的简单方法。 0013 但是， 多普勒效应的精度非常有限。 0014 多普勒超声波适用于医疗用超声波中像血流的监测那样精度要求不严密的应用 例。

10、 0015 但是， 在构造小并且动作极小的血管等是对象物、 要求高水准的精度的情况下， 需 要更高灵敏度的技术。 0016 作为最新的应用例， 具有为了估计构造的弹性而使用构造位移的超声波弹性率测 量法。 0017 为了准确估计弹性， 需要从所接收的超声波中准确地估计构造的位移。 0018 使用了超声波的位移估计， 与分辨率 ( 即， 估计极小位移的能力 ) 不同， 也需要高 水准的精度。 0019 在多个现有技术中， 从超声波 B 模式图像中估计构造的位移。 0020 但是， 所估计的位移的质量大大依赖于 B 模式图像的画质和分辨率。 0021 在大部分的超声波装置中， 很难用 B 模式图像

11、的分辨率来估计微米单位的位移。 0022 并且， 也有将如下情况作为焦点的技术 ： 从所接收到的 RF 信号中直接估计位移。 说 明 书 CN 102361597 A CN 102361614 A2/15 页 5 0023 最一般的技术之一是 【1】 这样的互相关的技术。 0024 然而， 在互相关技术中， 计算量增大， 只能估计与多个采样点对应的位移。 0025 微米单位的位移通常只不过相当于采样间隔的一个间隔的极小的一部分。 0026 因此， 无法使用互相关来进行估计。 0027 在这种状况下， 存在根据信号插值来估计位移的技术。 0028 但是， 该方法使处理时间增大， 并且估计质量依赖

12、于插值方法。 0029 自相关依赖于所接收到的 RF 信号的正交解调信号 ( 也作为基带信号被公知 ) 的 相位信息。 0030 该方法具有 【2】 所记载的能够估计与 RF 信号的子采样对应的位移这样的优点。 0031 但是， 在该方法中， 容易产生噪声， 并且受调幅效应影响。 0032 即， 所估计的位移在信号功率高的区域而产生了偏差。 0033 为了克服噪声的影响， 可以选择更大集合的采样来执行自相关， 但在该方法中， 估 计更详细的位移的能力下降。 0034 从上述技术中开发了一些克服所估计的位移的不准确性的方法。 0035 如 【3】 那样，“coarse to fine” 方法通过

13、在不同的估计阶段利用不同的窗区域， 在第 1 阶段粗略估计位移， 在第 2 阶段更细致地估计， 以此来提高精度。 0036 但是， 无法提供评价估计质量的方法。 0037 如 【1】 那样， 为了提高第 1 阶段的结果的精度， 也具有如下的技术 ： 通过执行第 2 估计阶段并使其与第 1 阶段的估计结果进行组合， 使用基于第 1 阶段的估计结果的信号卷 绕 (warping)。 0038 如 【2】 那样， 也能利用其他的位移校正方法。 0039 这些方法由于产生发散的风险高， 所以必须限制阶段数。 0040 另外， 没有用于克服因引导收敛并且 RF 信号功率的不均匀分布所引起的噪声和 偏差的

14、影响的方法。 0041 为了提供能够示出位移估计质量、 不需要限制估计阶段数、 引导反复估计的收敛、 针对极小位移精度也很高的更综合性的反复估计方法， 需要对这些方法进行改善。 0042 这样， 作为现有例， 公知有专利文件 1 6 的技术。 0043 现有技术文献 0044 专利文献 0045 【专利文献 1】 美国专利第 6277074 号说明书 0046 【专利文献 2】 美国专利第 6506158 号说明书 ( 美国专利申请公开第 2001/0034485 号说明书 ) 0047 【专利文献 3】 美国专利申请公开第 2008/0019609 号说明书 0048 【专利文献 4】 美国

15、专利申请公开第 2009/0221916 号说明书 0049 【专利文献 5】 美国专利第 6270459 号说明书 0050 【专利文献 6】 国际公开第 2008/038615 号 发明内容 0051 发明要解决的问题 说 明 书 CN 102361597 A CN 102361614 A3/15 页 6 0052 关于单一阶段的位移估计方法， 现在的方法在精度和位移分辨率的方面具有课 题。 0053 基于互相关的方法， 在没有信号插值的情况下， 无法估计与子采样间隔对应的小 的位移。 0054 然而， 在该方法中， 噪声的影响较小。 0055 因此， 在基于互相关的方法中， 虽然精度高，

16、 但是位移分辨率低。 0056 以提高互相关方法的位移分辨率为目的， 可以应用信号插值。 0057 但是， 这在估计对象的位移是微米单位的情况下， 会大幅增加处理功率。 0058 并且， 该方法的精度依赖于插值算法的质量。 0059 在基于自相关的方法中， 能够估计较小的位移， 但非常容易产生噪声。 0060 因此， 在该方法中， 虽然位移分辨率高， 但是精度低。 0061 因此， 为了使精度提高， 可以利用 2 阶段的估计方法。 0062 但是， 这些方法在综合性的反复方法的欠缺、 评价位移质量的方法的欠缺、 限制位 移精度的要素即消除噪声和调幅效应的方法的欠缺这 3 个方面具有课题。 00

17、63 需要在综合性的反复方法中组入引导反复估计的收敛的技术。 0064 本发明的目的在于提供一种能高精度地估计组织位移的位移估计装置 ( 位移估 计方法 )， 进而提供一种根据所估计的位移来提供适于判别恶性肿瘤和良性肿瘤或正常组 织的信息的位移估计装置。 0065 另外， 本发明的另一目的在于提供一种能够用较少的反复次数使反复的位移估计 处理高精度地收敛的位移估计装置。 0066 用于解决课题的手段 0067 为了解决上述课题， 本发明的位移估计方法使用超声波信号反复估计位移， 其中， 该位移估计方法包含如下步骤 ： 扫描至少 1 个超声波信号， 并发送给介质 ； 接收所扫描的从 所述介质反射

18、的所述超声波信号 ； 计算窗尺寸 ； 根据所计算的所述窗尺寸， 计算窗的边界 ； 使用基于所计算的所述边界的所述窗， 估计所述超声波信号的各个深度的位移 ； 根据所估 计的所述位移， 对所述超声波信号进行卷绕 ； 以及使用所卷绕的所述超声波信号， 引导该位 移估计方法的收敛， 使得所述超声波信号的相关值变大。 0068 在本发明中， 示出了发明者使用利用了反复估计方法的超声波来估计位移的方 法。 0069 根据发明者的实验， 在估计窗尺寸相同或者没有适当构成所有的反复估计循环 (round) 的情况下， 所估计的位移不收敛。 0070 因此， 导入对该收敛有用的窗计算方法。 0071 窗计算方

19、法包含如下的两个部分 ： 为了引导反复估计的收敛而计算窗尺寸的部 分、 以及为了克服因信号功率的不均匀分布而引起的噪声和偏差的影响而计算窗边界的部 分。 0072 通过窗尺寸的计算， 按照每个反复估计循环决定不同的窗尺寸。 0073 这具有如下的两个目的 ： 为了在后续的循环中获得更详细的估计结果、 以及为了 使估计可靠地收敛。 0074 即， 能够避免进行不合适的窗尺寸的处理而得不到收敛从而得不到准确的位移 说 明 书 CN 102361597 A CN 102361614 A4/15 页 7 量， 能够进行合适的窗尺寸的处理而可靠地得到收敛， 从而能够得到更准确的位移量。 0075 在窗边

20、界的计算中， 利用信号功率决定各估计位置的估计窗的边界。 0076 在以往的估计方法中， 窗通常以该位置为中心对称。 0077 这是因不均匀的信号功率的分布而引起偏差的一个原因。 0078 通过使用信号功率来决定窗边界， 能够克服这种影响。 0079 为了引导进一步的收敛， 导入所估计的位移的质量的评价， 反复估计方法保证始 终收敛到最佳的结果。 0080 本发明包含引导收敛的方法、 克服限制精度的要素的方法以及评价结果的质量的 方法， 使用综合性的反复方法来保证位移估计的精度。 0081 发明效果 0082 图 10 是示出在仿真中应用了回溯算法后的位移估计的改善结果的图。 0083 图 1

21、1 是表示体模实验中的相同结果的图。 0084 图 10 是示出估计的第 1 循环中的所估计的位移 ( 位移 1001) 脱离所仿真的位移 轮廓 ( 轮廓 1000) 的情况的图。 0085 然而， 在使用了引导收敛的窗方式的回溯算法的最终循环中， 所估计的位移 ( 位 移 1002) 与所仿真的轮廓非常吻合。 0086 图 11 是示出生成整体位移的体模实验的估计结果的图。 0087 为了例示， 选择单线 ( 参照轮廓 1100)。 0088 由于生成整体位移， 所以期待恒定的位移轮廓 ( 轮廓 1101)。 0089 在估计的第 1 循环中， 结果 ( 参照数据 1102) 不太适合期待的

22、位移。 0090 然而， 最终循环的估计结果 ( 参照数据 1103) 比起第 1 循环的估计结果， 更好地表 示了整体的动作。 0091 并且， 能够提供一种可进行合适的窗尺寸处理的位移估计装置。 0092 本发明的目的在于提供一种能高精度地估计组织位移的位移估计装置 ( 位移估 计方法 )， 进而提供一种根据所估计的位移来提供适于判别恶性肿瘤和良性肿瘤或正常组 织的信息的位移估计装置。 0093 另外， 本发明的另一目的在于提供一种能够用较少的反复次数使反复的位移估计 处理高精度地收敛的位移估计装置。 附图说明 0094 图 1 是表示本发明的位移估计方法的图。 0095 图 2 是表示本

23、发明的窗计算方法的图。 0096 图 3 是表示窗尺寸计算的例子的图。 0097 图 4 是表示对称窗和能量均匀窗的例子的图。 0098 图 5 是表示本发明的 RF 信号卷绕方法的图。 0099 图 6 是表示 RF 信号卷绕和 RF 采样延迟的图。 0100 图 7 是表示本发明的引导收敛的方法的图。 0101 图 8 是表示超声波 RF 信号的图。 0102 图 9 是表示使用了自相关的位移估计的图。 说 明 书 CN 102361597 A CN 102361614 A5/15 页 8 0103 图 10 是表示在仿真中使用了本发明的位移估计结果的例子的图。 0104 图 11 是表示

24、在体模实验中使用了本发明的位移估计结果的例子的图。 0105 图 12 是表示下层构造的位移和对 RF 信号的影响的图。 0106 图 13 是本装置的框图。 具体实施方式 0107 以下， 参照附图进行详细说明。 0108 实施方式的位移估计方法是使用超声波信号 ( 图 13 的超声波信号 1201s) 反复估 计位移 ( 参照图 12 的位移 1207) 的位移估计方法 ( 参照图 13 的位移估计装置 1X)。 0109 即， 例如通过进行多次位移大小的确定， 进行多次确定处理， 能够确定精度比较好 的位移大小。 0110 并且， 扫描至少 1 个超声波信号并发送 ( 参照图 13 的超

25、声波处理部 1Xa 中的发送 部 1X1) 给介质 ( 参照介质 1200M( 图 12) 以及介质 1200M 所具有的下层构造 ( 位移测定的 对象 )1203x( 图 12、 13)。 0111 即， 例如向线方向 ( 图 8 的方向 81L) 的多个位置发送超声波信号， 进行这些多个 位置的扫描。 0112 并且， 接收 ( 参照接收波 1X2) 所扫描的从所述介质反射的所述超声波信号 ( 超声 波信号 1201s)。 0113 另外， 超声波处理部 1Xa 例如是探头 ( 探针 ) 等。 0114 并且， 计算 ( 参照图 13 的信息处理部 1Xb 所包含的尺寸计算部 1X3) 窗

26、尺寸 ( 图 13 的窗尺寸 1X3M ： 例如图 3 的 (a) 栏、 (b) 栏的窗尺寸 (winSize) 等 )。 0115 即， 决定各个循环中的窗尺寸。 0116 另外， 尺寸计算部 1X3 例如可以包含窗计算部 102( 图 1)， 也可以是窗计算部 102 的至少一部分等。 0117 即， 例如根据第 n-1 循环中的相关值 ( 图 13 的相关值 1X7M) 来决定第 n 循环中的 窗尺寸 (n 2)。 0118 并且， 根据所计算出的所述窗尺寸， 计算 ( 参照边界计算部 1X4) 窗的边界 ( 参照 边界的信息 1X4M ： 图 4 的 (a) 栏、 (b) 栏所示的窗的

27、始端 (winStart) 以及终端 (winEnd) 等 )。 0119 另外， 边界计算部 1X4 例如可以是窗计算部 102( 图 1) 的至少一部分。 0120 即， 确定始端和终端， 这些始端和终端间的宽度是所计算的窗尺寸。 0121 并且， 使用基于所计算的所述边界的所述窗， 估计 ( 参照估计部 1X5) 所述超声波 信号的各个深度的位移 ( 位移 1XbM)。 0122 另外， 估计部 1X5 例如可以是位移估计部 101( 图 1) 的至少一部分等。 0123 即， 例如确定与本装置的距离相互不同的多个位置上的位移。 0124 并且， 根据所估计的所述位移， 卷绕所述超声波信

28、号 ( 参照卷绕部 1X6)。 0125 另外， 卷绕部 1X6 例如可以是 RF 信号卷绕部 100( 图 1) 的至少一部分等。 0126 并且， 使用所卷绕的所述超声波信号来引导(参照收敛控制部1X7)该位移估计方 法的收敛， 使得所述超声波信号的相关值 ( 相关值 1X7M) 增大。 说 明 书 CN 102361597 A CN 102361614 A6/15 页 9 0127 另外， 收敛控制部 1X7 例如可以是收敛引导部 103 的至少一部分等。 0128 即， 将在根据所卷绕的超声波信号计算出的相关值最大的循环中估计的位移大小 确定为精度最高的大小。 0129 并且， 使所计

29、算的所述窗尺寸变化， 使得用于连续多个循环的所述窗尺寸逐渐减 小。 0130 即， 例如在各个循环中计算出的窗尺寸是比在该循环之前的循环中计算出的窗尺 寸小的尺寸， 所计算的窗尺寸向小尺寸变更。 0131 并且， 计算所述窗尺寸， 使得所有 ( 深度 ) 的窗的信号能量相等。 0132 并且， 超声波信号的相关值 ( 第 n-1 循环中的相关值 ) 越大， 计算越小的所述窗尺 寸(第n循环中的窗尺寸)。 这里， 相关值等各循环中的收敛程度的评价值和窗尺寸的关系 可以参照基于活体上的实验值的预定关系式或表示两者的对应关系的表等来决定。 0133 并且， 将与第 n-1 循环中的相关值对应了所决定

30、的对应关系的窗尺寸确定为第 n 循环中的窗尺寸。 0134 即， 例如在本装置中进行下面的动作。 0135 即， 可以向作为测定对象物 ( 例如伴随着新生血管的恶性肿瘤等 ) 的下层构造 1203x( 图 12) 发送超声波信号 1201s， 并接收所发送的超声波信号 1201s。 0136 所确定的延迟时间 1208 是第 1 脉冲位置 ( 参照图 12 的 (b) 栏的 RF 信号 1204) 和第 2 脉冲位置 ( 参照 RF 信号 1205) 间的延迟时间。 0137 这里， 第 1 脉冲是所接收的超声波信号 1201s 中下层构造 1203 在进行位移 1207 之前的脉冲。 013

31、8 并且， 第 2 脉冲是测定对象物进行了位移 1207 后的脉冲。 0139 并且， 根据所确定的延迟时间 1208 来决定位移 1207。 0140 并且， 可以判定所确定的位移1207的大小的下层构造1203x的动作是否是预先设 定的动作。 0141 并且， 可以在判定为是预先设定的动作的情况下， 判定为该下层构造 1203x 是恶 性肿瘤 ( 癌 )， 在判定为不是该动作的情况下， 判定为是良性肿瘤或正常组织。 0142 另外， 例如可以通过收敛控制部 1x7( 图 13) 来进行这种是否是癌的判定。 0143 在这种本装置中， 更具体来讲， 例如可以如在后面详细说明那样， 在第 n-

32、1 循环中 确定的相关值越大， 在第 n 循环中， 越能进行更大窗尺寸的处理。 0144 另外， 这里， 基于窗尺寸的处理例如可以是从所接收的超声波信号 1201s 中的该 窗尺寸的窗部分的数据中计算出第 n 循环中的相关值的处理。 0145 根据本处理， 由于在各循环中能够设定与超声波的接收信号的信号电平对应的适 当窗， 所以能够提高收敛的精度， 并且以较少的反复次数获得收敛， 并且能够以较少的处理 量实现高精度的位移估计。 0146 以下的实施方式仅说明各发明步骤的原理。 0147 这里说明的具体例子的各种变形， 对于本领域技术人员来说是显而易见的。 0148 因此， 本发明的技术范围并不

33、限定于本说明书记载的具体例示的内容， 只限定于 权利要求书的范围。 0149 图 1 是表示位移估计方法的图。 说 明 书 CN 102361597 A CN 102361614 A7/15 页 10 0150 图 1 示出本发明的主要实施方式。 0151 RFSig(d， l， f)表示从超声波部获得的RF信号， d表示深度方向(图8的方向81D)， l 表示线方向 ( 方向 81L)， f 表示帧方向 ( 方向 81F)。 0152 图 8 进一步示出与这些相关的内容。 0153 这种 RF 信号的 1 帧内的 1 条线 ( 参照图示的 801) 通过由超声波变换器 800 在该 线位置向

34、扫描方向 ( 方向 802) 发送脉冲的机构生成。 0154 这里， 该脉冲一边反射和散射， 一边沿着进路与下层构造相互作用， 该超声波变换 器 800 接收反射和散射的信号。 0155 并且， 该信号被超声波部变换为对应的线。 0156 通过相同的机构， 在各个位置生成多条线， 并在各种时间实例(time instance)生 成多个帧。 0157 图 12 的 (a) 栏表示下层构造的位移， (b) 栏表示对 RF 信号的影响。 0158 图 12 示出所扫描的介质的下层构造的位移被作为 RF 信号的延迟而反射的机构。 0159 以例示为目的， 假设发送 1 条线的 1 个脉冲， 下层构造

35、只有 1 个反射边界。 0160 在第 1 时间实例 ( 参照 (a) 栏的左侧部分 )， 超声波变换器 1200 向所扫描的介质 发送脉冲 ( 参照脉冲 1201)。 0161 该脉冲通过介质传播到到达下层构造 1203(1203x) 的反射边界为止。 0162 通过该边界， 脉冲向变换器反射 ( 参照所反射的脉冲 1202)。 0163 变换器接收该脉冲， 变换为 RF 信号 1204( 参照图 12 的 (b) 栏的表的第 1 行 )。 0164 RF 信号 1204 的整体的脉冲位置表示被发送的脉冲 1201 向下层构造 1203 前进并 被反射回来所需的时间。 0165 在第 2 时

36、间实例 ( 图 12 的 (a) 栏的右侧 )， RF 信号 1205( 图 12 的 (b) 栏的表的 第 3 行 ) 由同一机构生成， 但下层构造位移到新的位置 ( 参照位移后的下层构造 1206)。 0166 由于位移 x( 位移 1207)， 所发送的脉冲 ( 脉冲 1201) 被传播并被下层构造的边界 反射所需的时间变长， 其结果， 引起 RF 信号的延迟 ( 参照延迟时间 1208)。 0167 另外， 在该例子中， 位移 x 的朝向是图 12 所示的方向 1203d 的朝向。 0168 如果知道了该延迟时间 1208 的延迟值， 就能够估计出位移 ( 位移 1207)。 0169

37、 并且， 在具有多个下层构造的介质中， 能够生成多个反射脉冲。 0170 这些多个下层构造的位移的量相互不同的原因在于 RF 信号中的反射脉冲的延迟 不同。 0171 因为这些脉冲重复并相互作用， 所以很难估计各下层构造的准确位移。 0172 本发明提示出通过反复进行估计来准确估计下层构造的位移的方法。 0173 根据图 1 所示的主要实施方式， 本发明具备以下的主要模块。 0174 即， 在本装置中， 具备 RF 信号卷绕部 100( 参照图 5)、 位移估计部 101( 参照图 9)、 窗计算部 102( 参照图 2)、 收敛引导部 103( 参照图 7)、 输出开关 104。 0175

38、RF 信号卷绕部 100 根据上次的估计循环 ( 第 n-1 循环 ) 的位移估计结果 dispRound(d， l， round)， 进行对象估计循环 ( 第 n 循环 ) 中的卷绕。在该卷绕中， 对于对象 估计循环， 卷绕 RF 信号的所选择的帧内的所选择的线。 0176 另外， 其在第 1 循环中被设定为 0。 说 明 书 CN 102361597 A CN 102361614 A8/15 页 11 0177 该模块的输出是被用于估计循环的被卷绕的 RF 信号 RFSigWarp(d， l， f， round)。 0178 另外， 在第 1 估计循环 ( 即， round 1) 中， 通

39、过将 dispRound(d， l， 0) 初始化为 0， 使得 RFSigWarp(d， l， f， 1) 与从超声波部获得的 RF 信号 RFSig(d， l， f)( 前述 ) 实质上 是相同的。 0179 由于估计位移， 所以可以将预先定义的帧的集合选择为针对该模块的输入。 0180 以例示为目的， 选择 2 个帧， 并用 f1 和 f2 来表示它们。 0181 RF 信号卷绕的目的是根据在先行的循环中所估计的位移将 1 个 RF 信号校正为与 其他 RF 信号相适合。 0182 并且， 之后通过从所校正的 RF 信号中估计残留位移， 并追加到先行循环的位移 上， 生成新的位移。 01

40、83 一些数量 ( 预先设定的数量 ) 的估计循环之后， 所估计的位移与所校正的 RF 信号 几乎一致， 将残留位移收敛成 0。 0184 位移估计部 101 通过根据 RF 信号的时间延迟进行导出从而根据 RF 信号 RFSigWarp(d， l， f， round) 中估计位移。 0185 优选执行该任务的方法是通过winStart(d， l)和winEnd(d， l)来记述各深度的估 计窗的自相关的方法， 但并不限于此。 0186 窗计算部 102 计算各深度的用于所选择的线的估计窗参数 winStart(d， l) 和 winEnd(d， l)。 0187 窗计算部 102 获取 RF

41、 信号 RFSigWarp(d， l， f， round)， 并且， 获取对象循环数 curRound 和 RF 信号差分 residue(d， l， round) 作为来自收敛引导部 103 的输入。 0188 收敛引导部 103 进行本技术的引导收敛的处理。 0189 即， 为了决定所估计的位移的质量， 计算所卷绕的 RF 信号间的差分。 0190 将该差分最小 ( 包含对象循环 ) 的循环判定为表示最准确的位移估计结果。 0191 该 模块的输出是对象估计循环数 curRound、 RF 信号差分最小的循环 数 minResidueRound(d， l)、 作为所估计的位移质量的测定值的

42、 RF 信号差分 residue(d， l， round)。 0192 residue(d， l， round) 越小， 表示质量越好， 在估计循环整体中， 如果 residue(d， l， round) 减少， 则表示收敛。 0193 另外， 作为计算所卷绕的 RF 信号间的差分的方法， 考虑了 “平方差分和” 和 “绝对 差分” 的方法， 但不限于此。 0194 输出开关 104 获取 curRound 和 minResidueRound(d， l) 作为输入。 0195 在对象循环数到达最大值的情况下， 选择最准确的估计结果作为最终输出。 0196 否则， 为了用于后续的估计循环， 将该

43、最准确的估计结果选择为 RF 信号卷绕部 100 的输入， 将 residue(d， l， round) 和 curRound 作为输入发送给位移估计部 101。 0197 下面的段落示出了一些基于图 1 所示的主要实施方式的具体实施方式。 0198 图 5 是表示 RF 信号卷绕方法的图。 0199 作为本发明的具体实施方式， 主要的实施方式的 RF 信号卷绕部 100 通过图 5 所示 的方法来实现。 0200 位移的结果是 RF 信号被移动并且被解压缩。 说 明 书 CN 102361597 A CN 102361614 A9/15 页 12 0201 因此， 如下所示， 在延迟值计算部

44、 500 中根据先行循环的最准确的位移估计结果 dispRound(d， l， minResidueRound(d， l) 计算使 RFSig(d， l， f1) 适合于 RFSig(d， l， f2) 的采样中的延迟值 ( 参照图 12 的延迟时间 1208)。 0202 【算式 1】 0203 0204 这里， fs 是 RF 信号的采样频率， c 是所扫描的介质中声速。 0205 作为结果， 将所扫描的介质中的位移反映成所获得的 RF 信号中的延迟。 0206 由于各深度的延迟值是已知的， 所以能够按照使 RFSig(d， l， f1) 或 RFSig(d， l， f2) 的任意一个中的

45、各线的 ( 与所有深度对应的 ) 所有采样适合于其他采样， 来进行延迟。 0207 这在 RF 采样延迟部 501 中执行。 0208 在 1 个实施方式中， 使用信号插值来执行 RF 采样延迟部 501 中的延迟。 0209 在其他实施方式中， 使用分数延迟滤波器来执行 RF 采样延迟部 501 中的延迟。 0210 图 6 是表示 RF 信号卷绕和 RF 采样延迟的图。 0211 图 6 示出 RF 信号卷绕处理。 0212 上述的图 12 示出了由于位移而引起了 RF 信号的延迟。 0213 RF 信号卷绕的目的是通过使 1 个 RF 信号向与所估计的位移相反的方向延迟而使 该效果无效。

46、 0214 作为该无效化的结果， 期待所卷绕的 RF 信号相互适合。 0215 延迟值计算部 500 所计算出的 delayValues(d， l) 是表示 1 个 RF 信号内的各采样 为了达成该适合而延迟的值。 0216 使用 2 个 RF 信号 600 和 601( 图 6) 来进行该例示。 0217 数据602表示在某个估计循环之后由延迟值计算部500计算出的delayValues(d， l)。 0218 可知 delayValues(d， l) 的特定深度的值表示用于使 RF 信号 600 中的相同深度的 采样适合于延迟值计算部 601 的延迟所需的量。 0219 然而， delay

47、Values(d， l) 中的延迟值并不始终是数。 0220 因此， 由于按照 delayValues(d， l) 能够延迟 RF 信号 600 中的各采样， 所以举出信 号插值或分数延迟滤波器作为两个选择候选。 0221 通过将位移估计部 101 应用于 RF 信号， 获得用于第 1 循环的位移和用于后续循环 的残留位移， 为了获得最终的位移， 对它们进行组合。 0222 图 9 示出优选的实施方式， 但并不限于此。 0223 图 9 表示使用了自相关的位移估计。 0224 IQ 解调部 900 将 RF 信号 RFSigWarp(d， l， f， round) 变换为基带信号 IQSig(

48、d， l， f)。 0225 自相关计算部 901 根据以下的式子计算 autocorr(d， l)。 0226 【算式 2】 0227 说 明 书 CN 102361597 A CN 102361614 A10/15 页 13 0228 这里， conj() 表示复共轭。 0229 最后， 根据以下的式子， 由位移计算部 902 将 autocorr(d， l) 变换为位移。 0230 【算式 3】 0231 0232 dispRound(d， l， round) dispRound(d， l， round-1)+residueDisp 0233 这里， arg() 是计算复数的辐角的函数。 0234 根据从 1 开始的循环 ( 表示第 1 循环 )， 将 dispRound(d， l， 0) 初始化为 0。 0235 图 2 是表示窗计算方法的图。 0236 通过窗计算部 102 针对各估计窗中的各线的各深度， 决定不同