位移估计方法、位移估计装置的制作方法

专利名称：位移估计方法、位移估计装置的制作方法
技术领域：
本发明导入使用超声波信号来估计位移的方法。无论是最终结果还是用于进一步处理的中间步骤，在需要位移估计的应用例中都能够使用，能够用于医疗用和工业用超声波设备。
背景技术：
超声波装置通过对扫描对象介质收发高频的机械声波，对介质进行非侵袭地诊断。这种装置的变换器将超声波发送给扫描对象介质。超声波通过散射和反射，与该介质的下层构造相互作用。另外，下层构造例如是指比介质的表面更靠近介质的内部侧的构造等。下层构造例如是人身体上的血管等。被散射和反射的超声波包含下层构造的有益信息，由变换器接收。然后由超声波装置进行处理，例如作为图像信息提示给用户。图8是表示超声波RF信号的图。作为超声波装置从反射波获得的最基础的数据的种类之一，列举了 RF (Radiofrequency Signal 射频信号)信号(图 8)。这是将接收波从模拟直接变换成数字的情况。在不同的应用例中，能够从超声波RF信号中导出亮度模式(B模式)图像、多普勒图像等其他种类的数据。多个应用例中的一个用于分析所扫描的介质的下层构造的动作。多普勒效应被用作估计构造的动作方向和强度的简单方法。但是，多普勒效应的精度非常有限。多普勒超声波适用于医疗用超声波中像血流的监测那样精度要求不严密的应用例。但是，在构造小并且动作极小的血管等是对象物、要求高水准的精度的情况下，需要更高灵敏度的技术。作为最新的应用例，具有为了估计构造的弹性而使用构造位移的超声波弹性率测量法。为了准确估计弹性，需要从所接收的超声波中准确地估计构造的位移。使用了超声波的位移估计，与分辨率(即，估计极小位移的能力)不同，也需要高水准的精度。在多个现有技术中，从超声波B模式图像中估计构造的位移。但是，所估计的位移的质量大大依赖于B模式图像的画质和分辨率。在大部分的超声波装置中，很难用B模式图像的分辨率来估计微米单位的位移。并且，也有将如下情况作为焦点的技术从所接收到的RF信号中直接估计位移。
最一般的技术之一是1这样的互相关的技术。然而，在互相关技术中，计算量增大，只能估计与多个采样点对应的位移。微米单位的位移通常只不过相当于采样间隔的一个间隔的极小的一部分。因此，无法使用互相关来进行估计。在这种状况下，存在根据信号插值来估计位移的技术。但是，该方法使处理时间增大，并且估计质量依赖于插值方法。自相关依赖于所接收到的RF信号的正交解调信号(也作为基带信号被公知)的相位信息。该方法具有2所记载的能够估计与RF信号的子采样对应的位移这样的优点。但是，在该方法中，容易产生噪声，并且受调幅效应影响。S卩，所估计的位移在信号功率高的区域而产生了偏差。为了克服噪声的影响，可以选择更大集合的采样来执行自相关，但在该方法中，估计更详细的位移的能力下降。从上述技术中开发了一些克服所估计的位移的不准确性的方法。如3那样，“coarse to fine”方法通过在不同的估计阶段利用不同的窗区域， 在第1阶段粗略估计位移，在第2阶段更细致地估计，以此来提高精度。但是，无法提供评价估计质量的方法。如1那样，为了提高第1阶段的结果的精度，也具有如下的技术通过执行第2 估计阶段并使其与第1阶段的估计结果进行组合，使用基于第1阶段的估计结果的信号卷绕(warping)。如2那样，也能利用其他的位移校正方法。这些方法由于产生发散的风险高，所以必须限制阶段数。另外，没有用于克服因引导收敛并且RF信号功率的不均勻分布所引起的噪声和偏差的影响的方法。为了提供能够示出位移估计质量、不需要限制估计阶段数、引导反复估计的收敛、 针对极小位移精度也很高的更综合性的反复估计方法，需要对这些方法进行改善。这样，作为现有例，公知有专利文件1 6的技术。现有技术文献专利文献专利文献1美国专利第6277074号说明书专利文献2美国专利第6506158号说明书(美国专利申请公开第2001/0034485 号说明书)专利文献3美国专利申请公开第2008/0019609号说明书专利文献4美国专利申请公开第2009/0221916号说明书专利文献5美国专利第6270459号说明书专利文献6国际公开第2008/038615号

发明内容
发明要解决的问题
关于单一阶段的位移估计方法，现在的方法在精度和位移分辨率的方面具有课题。基于互相关的方法，在没有信号插值的情况下，无法估计与子采样间隔对应的小的位移。然而，在该方法中，噪声的影响较小。因此，在基于互相关的方法中，虽然精度高，但是位移分辨率低。以提高互相关方法的位移分辨率为目的，可以应用信号插值。但是，这在估计对象的位移是微米单位的情况下，会大幅增加处理功率。并且，该方法的精度依赖于插值算法的质量。在基于自相关的方法中，能够估计较小的位移，但非常容易产生噪声。因此，在该方法中，虽然位移分辨率高，但是精度低。因此，为了使精度提高，可以利用2阶段的估计方法。但是，这些方法在综合性的反复方法的欠缺、评价位移质量的方法的欠缺、限制位移精度的要素即消除噪声和调幅效应的方法的欠缺这3个方面具有课题。需要在综合性的反复方法中组入引导反复估计的收敛的技术。本发明的目的在于提供一种能高精度地估计组织位移的位移估计装置(位移估计方法)，进而提供一种根据所估计的位移来提供适于判别恶性肿瘤和良性肿瘤或正常组织的信息的位移估计装置。另外，本发明的另一目的在于提供一种能够用较少的反复次数使反复的位移估计处理高精度地收敛的位移估计装置。用于解决课题的手段为了解决上述课题，本发明的位移估计方法使用超声波信号反复估计位移，其中， 该位移估计方法包含如下步骤扫描至少1个超声波信号，并发送给介质；接收所扫描的从所述介质反射的所述超声波信号；计算窗尺寸；根据所计算的所述窗尺寸，计算窗的边界； 使用基于所计算的所述边界的所述窗，估计所述超声波信号的各个深度的位移；根据所估计的所述位移，对所述超声波信号进行卷绕；以及使用所卷绕的所述超声波信号，引导该位移估计方法的收敛，使得所述超声波信号的相关值变大。在本发明中，示出了发明者使用利用了反复估计方法的超声波来估计位移的方法。根据发明者的实验，在估计窗尺寸相同或者没有适当构成所有的反复估计循环 (round)的情况下，所估计的位移不收敛。因此，导入对该收敛有用的窗计算方法。窗计算方法包含如下的两个部分为了引导反复估计的收敛而计算窗尺寸的部分、以及为了克服因信号功率的不均勻分布而引起的噪声和偏差的影响而计算窗边界的部分。通过窗尺寸的计算，按照每个反复估计循环决定不同的窗尺寸。这具有如下的两个目的为了在后续的循环中获得更详细的估计结果、以及为了使估计可靠地收敛。S卩，能够避免进行不合适的窗尺寸的处理而得不到收敛从而得不到准确的位移量，能够进行合适的窗尺寸的处理而可靠地得到收敛，从而能够得到更准确的位移量。在窗边界的计算中，利用信号功率决定各估计位置的估计窗的边界。在以往的估计方法中，窗通常以该位置为中心对称。这是因不均勻的信号功率的分布而引起偏差的一个原因。通过使用信号功率来决定窗边界，能够克服这种影响。为了引导进一步的收敛，导入所估计的位移的质量的评价，反复估计方法保证始终收敛到最佳的结果。本发明包含引导收敛的方法、克服限制精度的要素的方法以及评价结果的质量的方法，使用综合性的反复方法来保证位移估计的精度。发明效果图10是示出在仿真中应用了回溯算法后的位移估计的改善结果的图。图11是表示体模实验中的相同结果的图。图10是示出估计的第1循环中的所估计的位移(位移1001)脱离所仿真的位移轮廓(轮廓1000)的情况的图。然而，在使用了引导收敛的窗方式的回溯算法的最终循环中，所估计的位移(位移100 与所仿真的轮廓非常吻合。图11是示出生成整体位移的体模实验的估计结果的图。为了例示，选择单线(参照轮廓1100)。由于生成整体位移，所以期待恒定的位移轮廓(轮廓1101)。在估计的第1循环中，结果(参照数据110 不太适合期待的位移。然而，最终循环的估计结果(参照数据1103)比起第1循环的估计结果，更好地表示了整体的动作。并且，能够提供一种可进行合适的窗尺寸处理的位移估计装置。本发明的目的在于提供一种能高精度地估计组织位移的位移估计装置(位移估计方法)，进而提供一种根据所估计的位移来提供适于判别恶性肿瘤和良性肿瘤或正常组织的信息的位移估计装置。另外，本发明的另一目的在于提供一种能够用较少的反复次数使反复的位移估计处理高精度地收敛的位移估计装置。


图1是表示本发明的位移估计方法的图。
图2是表示本发明的窗计算方法的图。
图3是表示窗尺寸计算的例子的图。
图4是表示对称窗和能量均勻窗的例子的图。
图5是表示本发明的RF信号卷绕方法的图。
图6是表示RF信号卷绕和RF采样延迟的图。
图7是表示本发明的引导收敛的方法的图。
图8是表示超声波RF信号的图。
图9是表示使用了自相关的位移估计的图。
图10是表示在仿真中使用了本发明的位移估计结果的例子的图。图11是表示在体模实验中使用了本发明的位移估计结果的例子的图。图12是表示下层构造的位移和对RF信号的影响的图。图13是本装置的框图。
具体实施方式
以下，参照附图进行详细说明。实施方式的位移估计方法是使用超声波信号(图13的超声波信号1201s)反复估计位移(参照图12的位移1207)的位移估计方法(参照图13的位移估计装置IX)。S卩，例如通过进行多次位移大小的确定，进行多次确定处理，能够确定精度比较好的位移大小。并且，扫描至少1个超声波信号并发送(参照图13的超声波处理部1 中的发送部1X1)给介质(参照介质1200M(图12)以及介质1200M所具有的下层构造(位移测定的对象)1203x(图 12、13))。即，例如向线方向(图8的方向81L)的多个位置发送超声波信号，进行这些多个位置的扫描。并且，接收(参照接收波1X2)所扫描的从所述介质反射的所述超声波信号(超声波信号1201s)。另外，超声波处理部1 例如是探头(探针)等。并且，计算(参照图13的信息处理部l)(b所包含的尺寸计算部1X3)窗尺寸(图 13的窗尺寸1X3M 例如图3的(a)栏、(b)栏的窗尺寸(winSize)等)。S卩，决定各个循环中的窗尺寸。另外，尺寸计算部1X3例如可以包含窗计算部102(图1)，也可以是窗计算部102 的至少一部分等。S卩，例如根据第n-1循环中的相关值(图13的相关值1X7M)来决定第η循环中的窗尺寸(η彡2)。并且，根据所计算出的所述窗尺寸，计算(参照边界计算部1X4)窗的边界(参照边界的信息1Χ4Μ 图4的(a)栏、(b)栏所示的窗的始端(winStart)以及终端(winEnd)
寸乂 O另外，边界计算部1X4例如可以是窗计算部102(图1)的至少一部分。S卩，确定始端和终端，这些始端和终端间的宽度是所计算的窗尺寸。并且，使用基于所计算的所述边界的所述窗，估计(参照估计部1X5)所述超声波信号的各个深度的位移(位移UbM)。另外，估计部1)(5例如可以是位移估计部101(图1)的至少一部分等。S卩，例如确定与本装置的距离相互不同的多个位置上的位移。并且，根据所估计的所述位移，卷绕所述超声波信号(参照卷绕部1X6)。另外，卷绕部1)(6例如可以是RF信号卷绕部100(图1)的至少一部分等。并且，使用所卷绕的所述超声波信号来引导(参照收敛控制部1X7)该位移估计方法的收敛，使得所述超声波信号的相关值(相关值1X7M)增大。
8
另外，收敛控制部1X7例如可以是收敛引导部103的至少一部分等。S卩，将在根据所卷绕的超声波信号计算出的相关值最大的循环中估计的位移大小确定为精度最高的大小。并且，使所计算的所述窗尺寸变化，使得用于连续多个循环的所述窗尺寸逐渐减S卩，例如在各个循环中计算出的窗尺寸是比在该循环之前的循环中计算出的窗尺寸小的尺寸，所计算的窗尺寸向小尺寸变更。并且，计算所述窗尺寸，使得所有(深度)的窗的信号能量相等。并且，超声波信号的相关值(第n-1循环中的相关值)越大，计算越小的所述窗尺寸(第η循环中的窗尺寸)。这里，相关值等各循环中的收敛程度的评价值和窗尺寸的关系可以参照基于活体上的实验值的预定关系式或表示两者的对应关系的表等来决定。并且，将与第n-1循环中的相关值对应了所决定的对应关系的窗尺寸确定为第η 循环中的窗尺寸。S卩，例如在本装置中进行下面的动作。即，可以向作为测定对象物(例如伴随着新生血管的恶性肿瘤等)的下层构造 1203χ(图12)发送超声波信号1201s，并接收所发送的超声波信号1201s。所确定的延迟时间1208是第1脉冲位置(参照图12的(b)栏的RF信号1204) 和第2脉冲位置(参照RF信号1205)间的延迟时间。这里，第1脉冲是所接收的超声波信号1201s中下层构造1203在进行位移1207 之前的脉冲。并且，第2脉冲是测定对象物进行了位移1207后的脉冲。并且，根据所确定的延迟时间1208来决定位移1207。并且，可以判定所确定的位移1207的大小的下层构造1203x的动作是否是预先设定的动作。并且，可以在判定为是预先设定的动作的情况下，判定为该下层构造1203x是恶性肿瘤(癌)，在判定为不是该动作的情况下，判定为是良性肿瘤或正常组织。另外，例如可以通过收敛控制部1χ7(图1 来进行这种是否是癌的判定。在这种本装置中，更具体来讲，例如可以如在后面详细说明那样，在第n-1循环中确定的相关值越大，在第η循环中，越能进行更大窗尺寸的处理。另外，这里，基于窗尺寸的处理例如可以是从所接收的超声波信号1201s中的该窗尺寸的窗部分的数据中计算出第η循环中的相关值的处理。根据本处理，由于在各循环中能够设定与超声波的接收信号的信号电平对应的适当窗，所以能够提高收敛的精度，并且以较少的反复次数获得收敛，并且能够以较少的处理量实现高精度的位移估计。以下的实施方式仅说明各发明步骤的原理。这里说明的具体例子的各种变形，对于本领域技术人员来说是显而易见的。因此，本发明的技术范围并不限定于本说明书记载的具体例示的内容，只限定于权利要求书的范围。图1是表示位移估计方法的图。
9
图1示出本发明的主要实施方式。RFSig (d，1，f)表示从超声波部获得的RF信号，d表示深度方向(图8的方向81D)， 1表示线方向(方向81L)，f表示帧方向(方向81F)。图8进一步示出与这些相关的内容。这种RF信号的1帧内的1条线(参照图示的801)通过由超声波变换器800在该线位置向扫描方向(方向802)发送脉冲的机构生成。这里，该脉冲一边反射和散射，一边沿着进路与下层构造相互作用，该超声波变换器800接收反射和散射的信号。并且，该信号被超声波部变换为对应的线。通过相同的机构，在各个位置生成多条线，并在各种时间实例(time instance)生成多个帧。图12的(a)栏表示下层构造的位移，(b)栏表示对RF信号的影响。图12示出所扫描的介质的下层构造的位移被作为RF信号的延迟而反射的机构。以例示为目的，假设发送1条线的1个脉冲，下层构造只有1个反射边界。在第1时间实例(参照(a)栏的左侧部分)，超声波变换器1200向所扫描的介质发送脉冲(参照脉冲1201)。该脉冲通过介质传播到到达下层构造1203(1203χ)的反射边界为止。通过该边界，脉冲向变换器反射(参照所反射的脉冲1202)。变换器接收该脉冲，变换为RF信号1204(参照图12的(b)栏的表的第1行)。RF信号1204的整体的脉冲位置表示被发送的脉冲1201向下层构造1203前进并被反射回来所需的时间。在第2时间实例(图12的(a)栏的右侧)，RF信号1205 (图12的(b)栏的表的第3行)由同一机构生成，但下层构造位移到新的位置(参照位移后的下层构造1206)。由于位移x(位移1207)，所发送的脉冲(脉冲1201)被传播并被下层构造的边界反射所需的时间变长，其结果，引起RF信号的延迟(参照延迟时间1208)。另外，在该例子中，位移χ的朝向是图12所示的方向1203d的朝向。如果知道了该延迟时间1208的延迟值，就能够估计出位移(位移1207)。并且，在具有多个下层构造的介质中，能够生成多个反射脉冲。这些多个下层构造的位移的量相互不同的原因在于RF信号中的反射脉冲的延迟不同。因为这些脉冲重复并相互作用，所以很难估计各下层构造的准确位移。本发明提示出通过反复进行估计来准确估计下层构造的位移的方法。根据图1所示的主要实施方式，本发明具备以下的主要模块。即，在本装置中，具备RF信号卷绕部100 (参照图5)、位移估计部101 (参照图9)、 窗计算部102(参照图2)、收敛引导部103(参照图7)、输出开关104。RF信号卷绕部100根据上次的估计循环(第n-1循环)的位移估计结果 dispR0imd(d，l，round)，进行对象估计循环(第η循环)中的卷绕。在该卷绕中，对于对象估计循环，卷绕RF信号的所选择的帧内的所选择的线。另外，其在第1循环中被设定为0。
该模块的输出是被用于估计循环的被卷绕的RF信号RFSigWarp(d，1，f，round)。另外，在第1估计循环(即，round = 1)中，通过将dispRound(d, 1,0)初始化为 0，使得RFSigffarp (d, 1，f，1)与从超声波部获得的RF信号RFSig (d, 1，f)(前述)实质上是相同的。由于估计位移，所以可以将预先定义的帧的集合选择为针对该模块的输入。以例示为目的，选择2个巾贞，并用f 1和f2来表示它们。RF信号卷绕的目的是根据在先行的循环中所估计的位移将1个RF信号校正为与其他RF信号相适合。并且，之后通过从所校正的RF信号中估计残留位移，并追加到先行循环的位移上，生成新的位移。一些数量(预先设定的数量)的估计循环之后，所估计的位移与所校正的RF信号几乎一致，将残留位移收敛成0。位移估计部101通过根据RF信号的时间延迟进行导出从而根据RF信号 RFSigffarp (d, 1，f, round)中估计位移。优选执行该任务的方法是通过winStart (d，1)和WinEnd(d，1)来记述各深度的估计窗的自相关的方法，但并不限于此。窗计算部102计算各深度的用于所选择的线的估计窗参数wir^tart(d，1)和 winEnd(d,1)。窗计算部102获取RF信号RFSigWarp (d，1，f，round)，并且，获取对象循环数 curRound和RF信号差分residue(d，1，round)作为来自收敛引导部103的输入。收敛引导部103进行本技术的引导收敛的处理。即，为了决定所估计的位移的质量，计算所卷绕的RF信号间的差分。将该差分最小(包含对象循环)的循环判定为表示最准确的位移估计结果。该模块的输出是对象估计循环数curRound、RF信号差分最小的循环数 minResidueRound(d, 1)、作为所估计的位移质量的测定值的RF信号差分residue (d，1， round)。residue (d，1，round)越小，表示质量越好，在估计循环整体中，如果residue (d， 1, round)减少，则表示收敛。另外，作为计算所卷绕的RF信号间的差分的方法，考虑了“平方差分和”和“绝对差分”的方法，但不限于此。输出开关104 获取 curRound 和 minResidueRound (d，1)作为输入。在对象循环数到达最大值的情况下，选择最准确的估计结果作为最终输出。否则，为了用于后续的估计循环，将该最准确的估计结果选择为RF信号卷绕部 100的输入，将residue (d，1，round)和curRound作为输入发送给位移估计部101。下面的段落示出了一些基于图1所示的主要实施方式的具体实施方式
。图5是表示RF信号卷绕方法的图。作为本发明的具体实施方式
，主要的实施方式的RF信号卷绕部100通过图5所示的方法来实现。位移的结果是RF信号被移动并且被解压缩。
因此，如下所示，在延迟值计算部500中根据先行循环的最准确的位移估计结果 dispRound(d, 1, minResidueRound(d, 1))计算使 RFSig(d，1，fl)适合于 RFSig(d，1，f2) 的采样中的延迟值(参照图12的延迟时间1208)。算式1这里，fs是RF信号的采样频率，c是所扫描的介质中声速。作为结果，将所扫描的介质中的位移反映成所获得的RF信号中的延迟。由于各深度的延迟值是已知的，所以能够按照使RFSig(d，1，fl)或RFSig(d，1， f2)的任意一个中的各线的(与所有深度对应的)所有采样适合于其他采样，来进行延迟。这在RF采样延迟部501中执行。在1个实施方式中，使用信号插值来执行RF采样延迟部501中的延迟。在其他实施方式中，使用分数延迟滤波器来执行RF采样延迟部501中的延迟。图6是表示RF信号卷绕和RF采样延迟的图。图6示出RF信号卷绕处理。上述的图12示出了由于位移而引起了 RF信号的延迟。RF信号卷绕的目的是通过使1个RF信号向与所估计的位移相反的方向延迟而使该效果无效。作为该无效化的结果，期待所卷绕的RF信号相互适合。延迟值计算部500所计算出的delayValueS(d，l)是表示1个RF信号内的各采样为了达成该适合而延迟的值。使用2个RF信号600和601 (图6)来进行该例示。数据602表示在某个估计循环之后由延迟值计算部500计算出的delayValues (d， 1)。可知delayValues (d，1)的特定深度的值表示用于使RF信号600中的相同深度的采样适合于延迟值计算部601的延迟所需的量。然而，delayValues (d, 1)中的延迟值并不始终是数。因此，由于按照delayValues (d，l)能够延迟RF信号600中的各采样，所以举出信号插值或分数延迟滤波器作为两个选择候选。 通过将位移估计部101应用于RF信号，获得用于第1循环的位移和用于后续循环的残留位移，为了获得最终的位移，对它们进行组合。图9示出优选的实施方式，但并不限于此。图9表示使用了自相关的位移估计。IQ 解调部 900 将 RF 信号 RFSigffarp (d, 1，f, round)变换为基带信号 IQSig (d，1， f)。自相关计算部901根据以下的式子计算autocorHd，1)。算式2
wmEnd ( d J )autocorr(d,l)= J] IQSigfd,Ii.,/Jx conj(IQSig(d,I；,/,))
Ii=WinSTarTl ci J)
这里，conj ()表示复共轭。最后，根据以下的式子，由位移计算部902将autocorr (d，1)变换为位移。算式3^.-,/； ir^v nr/” u J ”dispRound(d，1，round) = dispRound(d，1，round-1)+residueDisp这里，arg()是计算复数的辐角的函数。根据从1开始的循环(表示第1循环)，将dispRouncKd，1,0)初始化为0。图2是表示窗计算方法的图。通过窗计算部102针对各估计窗中的各线的各深度，决定不同的估计窗。这在图2中示出。图3是表示窗尺寸计算的具体例的图。窗尺寸计算部200(图2)决定针对各估计循环的窗尺寸。该模块的输出是winSize (d，1)。在优选实施方式中，针对各循环的窗尺寸随着循环的函数的减少而变化。这在图3中示出。在具有窗尺寸计算部200(图2)的实施方式中，如图3的(a)栏所示，窗尺寸由于减少函数winSizeRound (round)而直接关系到估计循环(round)数。关于由curRoimd所确定的对象估计循环，根据以下，确定针对所有深度的窗尺寸。关于所有的 d 禾口 1，winSize (d，1) = winSizeRound (curRound)。“固定窗方案”是在用于位移估计的所有深度选择成为估计窗尺寸的winSiZe(d， 1)的相同值的处理。“固定窗方案”的优点是适合于所需的运算量不多(比较少)并且RF信号能量沿着深度方向被规则准确地分配的情况。在窗尺寸计算部200的其他实施方式中，各深度的窗尺寸根据作为上限的 maxWin (round)、作为下限的 minWin (round)、RF 信号功率 sigPow(d，1)来计算。maxffin (round)禾口 minWin (round)是循环的减少函数。这作为“瞬间功率窗方案”，在图3的(b)栏中示出。应用以下的条件。算式4winSize(d,1) e [min Win(curRound), max Win(curRound)]并且，如图3 的(b)栏所示，winSize (d，l)也可以与 sigPow(d，1)、即与 sigPow = 0对应的最大窗尺寸和与最大sigPow对应的最小窗尺寸成反比。“瞬间功率窗方案”与“固定窗方案”相比，所需的运算量多，但是，特别在RF信号能量沿着深度方向相当不均勻的情况下，精度更高。为了克服信号对噪声比的较低程度，对具有较低信号功率的区域指定更大的估计窗，为了更详细地分布位移，对具有较大信号功率的区域指定更小的估计窗。然而，关于规则准确地分配了信号能量的RF信号，“固定窗方案”和“瞬间功率窗
13方案”没有大的差异。在窗尺寸计算部200的其他实施方式中，也可以计算所有深度的窗尺寸，使得所有的窗内的信号能量保持同量，这被称为“恒定能量窗方案”。信号功率计算部201(图2)计算沿着帧方向的RF信号的平均功率。信号功率计算部201获取所卷绕的RF信号RFSigWarp (d，1，f，round)和对象循环数curRound作为输入。在RFSigWarp (d，1，f，curRound)内具有多个帧，按照如下方式计算功率。算式5
权利要求
1.一种位移估计方法，使用超声波信号反复估计位移，其中， 该位移估计方法包含如下步骤扫描至少1个超声波信号，并发送给介质； 接收所扫描的从所述介质反射的所述超声波信号； 计算窗尺寸；根据所计算的所述窗尺寸，计算窗的边界；使用基于所计算的所述边界的所述窗，估计所述超声波信号的各个深度的位移； 根据所估计的所述位移，对所述超声波信号进行卷绕；以及使用所卷绕的所述超声波信号，引导该位移估计方法的收敛，使得所述超声波信号的相关值变大。
2.根据权利要求1所述的位移估计方法，其中，改变所计算的所述窗尺寸，使得连续的多个循环中使用的所述窗尺寸逐渐减小。
3.根据权利要求1所述的位移估计方法，其中， 该位移估计方法还包含如下步骤使连续的多个所述循环的所述窗尺寸的上限和下限分别按照逐渐减小的方式变化；以及根据各深度的信号功率，计算该深度的上限和下限间的窗尺寸，作为各个深度的所述窗尺寸。
4.根据权利要求1所述的位移估计方法，其中， 计算所述窗尺寸，使得所有窗中的信号能量相等。
5.根据权利要求1所述的位移估计方法，其中， 超声波信号的相关值越大，计算出越小的所述窗尺寸。
6.根据权利要求1所述的位移估计方法，其中，所述边界向该深度的两侧延伸，使得基于该边界的所述窗成为以对应的深度为中心的对称的窗。
7.根据权利要求1所述的位移估计方法，其中，所述边界向对应的所述深度的两侧延伸，使得对应深度的基于该边界的所述窗的两侧的信号能量相互相等。
8.根据权利要求1所述的位移估计方法，其中， 该位移估计方法还包含如下步骤根据所计算的延迟值使所述超声波信号内的各采样延迟，由此卷绕所述超声波信号。
9.根据权利要求1所述的位移估计方法，其中， 该位移估计方法还包含如下步骤在各循环之后，计算所卷绕的RF信号间的信号差分；以及决定所计算的所述信号差分为最小的循环。
10.一种位移估计装置，使用超声波信号反复估计位移，其中， 该位移估计装置包含发送部，扫描至少1个超声波信号，并发送给介质； 接收部，接收所扫描的从所述介质反射的所述超声波信号；尺寸计算部，计算窗尺寸；边界计算部，根据所计算的所述窗尺寸，计算窗的边界；估计部，使用基于所计算的所述边界的所述窗，估计所述超声波信号的各个深度的位移；卷绕部，根据所估计的所述位移，对所述超声波信号进行卷绕；以及收敛控制部，使用所卷绕的所述超声波信号，引导该位移估计装置执行的位移估计方法的收敛，使得所述超声波信号的相关值变大。
全文摘要
本发明公开一种位移估计方法，使用超声波信号反复估计位移，其中，包含如下步骤扫描至少1个超声波信号，并发送给介质；接收所扫描的从所述介质反射的所述超声波信号；计算窗尺寸；根据所计算的所述窗尺寸，计算窗的边界；使用所计算的所述窗，估计所述超声波信号的各个深度的位移；根据所估计的所述位移，对所述超声波信号进行卷绕；以及使用所卷绕的所述超声波信号，引导该位移估计方法的收敛，使得所述超声波信号的相关值变大。
文档编号G01N29/06GK102361597SQ201180001460
公开日2012年2月22日 申请日期2011年1月19日 优先权日2010年1月20日
发明者张国成, 范树峰, 陈英俊 申请人:松下电器产业株式会社