專利名稱：位移估計方法、位移估計裝置的制作方法
技術領域：
本發(fā)明導入使用超聲波信號來估計位移的方法。無論是最終結果還是用于進一步處理的中間步驟，在需要位移估計的應用例中都能夠使用，能夠用于醫(yī)療用和工業(yè)用超聲波設備。
背景技術：
超聲波裝置通過對掃描對象介質(zhì)收發(fā)高頻的機械聲波，對介質(zhì)進行非侵襲地診斷。這種裝置的變換器將超聲波發(fā)送給掃描對象介質(zhì)。超聲波通過散射和反射，與該介質(zhì)的下層構造相互作用。另外，下層構造例如是指比介質(zhì)的表面更靠近介質(zhì)的內(nèi)部側的構造等。下層構造例如是人身體上的血管等。被散射和反射的超聲波包含下層構造的有益信息，由變換器接收。然后由超聲波裝置進行處理，例如作為圖像信息提示給用戶。圖8是表示超聲波RF信號的圖。作為超聲波裝置從反射波獲得的最基礎的數(shù)據(jù)的種類之一，列舉了 RF (Radiofrequency Signal 射頻信號)信號(圖 8)。這是將接收波從模擬直接變換成數(shù)字的情況。在不同的應用例中，能夠從超聲波RF信號中導出亮度模式(B模式)圖像、多普勒圖像等其他種類的數(shù)據(jù)。多個應用例中的一個用于分析所掃描的介質(zhì)的下層構造的動作。多普勒效應被用作估計構造的動作方向和強度的簡單方法。但是，多普勒效應的精度非常有限。多普勒超聲波適用于醫(yī)療用超聲波中像血流的監(jiān)測那樣精度要求不嚴密的應用例。但是，在構造小并且動作極小的血管等是對象物、要求高水準的精度的情況下，需要更高靈敏度的技術。作為最新的應用例，具有為了估計構造的彈性而使用構造位移的超聲波彈性率測量法。為了準確估計彈性，需要從所接收的超聲波中準確地估計構造的位移。使用了超聲波的位移估計，與分辨率(即，估計極小位移的能力)不同，也需要高水準的精度。在多個現(xiàn)有技術中，從超聲波B模式圖像中估計構造的位移。但是，所估計的位移的質(zhì)量大大依賴于B模式圖像的畫質(zhì)和分辨率。在大部分的超聲波裝置中，很難用B模式圖像的分辨率來估計微米單位的位移。并且，也有將如下情況作為焦點的技術從所接收到的RF信號中直接估計位移。
最一般的技術之一是1這樣的互相關的技術。然而，在互相關技術中，計算量增大，只能估計與多個采樣點對應的位移。微米單位的位移通常只不過相當于采樣間隔的一個間隔的極小的一部分。因此，無法使用互相關來進行估計。在這種狀況下，存在根據(jù)信號插值來估計位移的技術。但是，該方法使處理時間增大，并且估計質(zhì)量依賴于插值方法。自相關依賴于所接收到的RF信號的正交解調(diào)信號(也作為基帶信號被公知)的相位信息。該方法具有2所記載的能夠估計與RF信號的子采樣對應的位移這樣的優(yōu)點。但是，在該方法中，容易產(chǎn)生噪聲，并且受調(diào)幅效應影響。S卩，所估計的位移在信號功率高的區(qū)域而產(chǎn)生了偏差。為了克服噪聲的影響，可以選擇更大集合的采樣來執(zhí)行自相關，但在該方法中，估計更詳細的位移的能力下降。從上述技術中開發(fā)了一些克服所估計的位移的不準確性的方法。如3那樣，“coarse to fine”方法通過在不同的估計階段利用不同的窗區(qū)域， 在第1階段粗略估計位移，在第2階段更細致地估計，以此來提高精度。但是，無法提供評價估計質(zhì)量的方法。如1那樣，為了提高第1階段的結果的精度，也具有如下的技術通過執(zhí)行第2 估計階段并使其與第1階段的估計結果進行組合，使用基于第1階段的估計結果的信號卷繞(warping)。如2那樣，也能利用其他的位移校正方法。這些方法由于產(chǎn)生發(fā)散的風險高，所以必須限制階段數(shù)。另外，沒有用于克服因引導收斂并且RF信號功率的不均勻分布所引起的噪聲和偏差的影響的方法。為了提供能夠示出位移估計質(zhì)量、不需要限制估計階段數(shù)、引導反復估計的收斂、 針對極小位移精度也很高的更綜合性的反復估計方法，需要對這些方法進行改善。這樣，作為現(xiàn)有例，公知有專利文件1 6的技術?，F(xiàn)有技術文獻專利文獻專利文獻1美國專利第6277074號說明書專利文獻2美國專利第6506158號說明書(美國專利申請公開第2001/0034485 號說明書)專利文獻3美國專利申請公開第2008/0019609號說明書專利文獻4美國專利申請公開第2009/0221916號說明書專利文獻5美國專利第6270459號說明書專利文獻6國際公開第2008/038615號

發(fā)明內(nèi)容
發(fā)明要解決的問題
關于單一階段的位移估計方法，現(xiàn)在的方法在精度和位移分辨率的方面具有課題?；诨ハ嚓P的方法，在沒有信號插值的情況下，無法估計與子采樣間隔對應的小的位移。然而，在該方法中，噪聲的影響較小。因此，在基于互相關的方法中，雖然精度高，但是位移分辨率低。以提高互相關方法的位移分辨率為目的，可以應用信號插值。但是，這在估計對象的位移是微米單位的情況下，會大幅增加處理功率。并且，該方法的精度依賴于插值算法的質(zhì)量。在基于自相關的方法中，能夠估計較小的位移，但非常容易產(chǎn)生噪聲。因此，在該方法中，雖然位移分辨率高，但是精度低。因此，為了使精度提高，可以利用2階段的估計方法。但是，這些方法在綜合性的反復方法的欠缺、評價位移質(zhì)量的方法的欠缺、限制位移精度的要素即消除噪聲和調(diào)幅效應的方法的欠缺這3個方面具有課題。需要在綜合性的反復方法中組入引導反復估計的收斂的技術。本發(fā)明的目的在于提供一種能高精度地估計組織位移的位移估計裝置(位移估計方法)，進而提供一種根據(jù)所估計的位移來提供適于判別惡性腫瘤和良性腫瘤或正常組織的信息的位移估計裝置。另外，本發(fā)明的另一目的在于提供一種能夠用較少的反復次數(shù)使反復的位移估計處理高精度地收斂的位移估計裝置。用于解決課題的手段為了解決上述課題，本發(fā)明的位移估計方法使用超聲波信號反復估計位移，其中， 該位移估計方法包含如下步驟掃描至少1個超聲波信號，并發(fā)送給介質(zhì)；接收所掃描的從所述介質(zhì)反射的所述超聲波信號；計算窗尺寸；根據(jù)所計算的所述窗尺寸，計算窗的邊界； 使用基于所計算的所述邊界的所述窗，估計所述超聲波信號的各個深度的位移；根據(jù)所估計的所述位移，對所述超聲波信號進行卷繞；以及使用所卷繞的所述超聲波信號，引導該位移估計方法的收斂，使得所述超聲波信號的相關值變大。在本發(fā)明中，示出了發(fā)明者使用利用了反復估計方法的超聲波來估計位移的方法。根據(jù)發(fā)明者的實驗，在估計窗尺寸相同或者沒有適當構成所有的反復估計循環(huán) (round)的情況下，所估計的位移不收斂。因此，導入對該收斂有用的窗計算方法。窗計算方法包含如下的兩個部分為了引導反復估計的收斂而計算窗尺寸的部分、以及為了克服因信號功率的不均勻分布而引起的噪聲和偏差的影響而計算窗邊界的部分。通過窗尺寸的計算，按照每個反復估計循環(huán)決定不同的窗尺寸。這具有如下的兩個目的為了在后續(xù)的循環(huán)中獲得更詳細的估計結果、以及為了使估計可靠地收斂。S卩，能夠避免進行不合適的窗尺寸的處理而得不到收斂從而得不到準確的位移量，能夠進行合適的窗尺寸的處理而可靠地得到收斂，從而能夠得到更準確的位移量。在窗邊界的計算中，利用信號功率決定各估計位置的估計窗的邊界。在以往的估計方法中，窗通常以該位置為中心對稱。這是因不均勻的信號功率的分布而引起偏差的一個原因。通過使用信號功率來決定窗邊界，能夠克服這種影響。為了引導進一步的收斂，導入所估計的位移的質(zhì)量的評價，反復估計方法保證始終收斂到最佳的結果。本發(fā)明包含引導收斂的方法、克服限制精度的要素的方法以及評價結果的質(zhì)量的方法，使用綜合性的反復方法來保證位移估計的精度。發(fā)明效果圖10是示出在仿真中應用了回溯算法后的位移估計的改善結果的圖。圖11是表示體模實驗中的相同結果的圖。圖10是示出估計的第1循環(huán)中的所估計的位移(位移1001)脫離所仿真的位移輪廓(輪廓1000)的情況的圖。然而，在使用了引導收斂的窗方式的回溯算法的最終循環(huán)中，所估計的位移(位移100 與所仿真的輪廓非常吻合。圖11是示出生成整體位移的體模實驗的估計結果的圖。為了例示，選擇單線(參照輪廓1100)。由于生成整體位移，所以期待恒定的位移輪廓(輪廓1101)。在估計的第1循環(huán)中，結果(參照數(shù)據(jù)110 不太適合期待的位移。然而，最終循環(huán)的估計結果(參照數(shù)據(jù)1103)比起第1循環(huán)的估計結果，更好地表示了整體的動作。并且，能夠提供一種可進行合適的窗尺寸處理的位移估計裝置。本發(fā)明的目的在于提供一種能高精度地估計組織位移的位移估計裝置(位移估計方法)，進而提供一種根據(jù)所估計的位移來提供適于判別惡性腫瘤和良性腫瘤或正常組織的信息的位移估計裝置。另外，本發(fā)明的另一目的在于提供一種能夠用較少的反復次數(shù)使反復的位移估計處理高精度地收斂的位移估計裝置。


圖1是表示本發(fā)明的位移估計方法的圖。
圖2是表示本發(fā)明的窗計算方法的圖。
圖3是表示窗尺寸計算的例子的圖。
圖4是表示對稱窗和能量均勻窗的例子的圖。
圖5是表示本發(fā)明的RF信號卷繞方法的圖。
圖6是表示RF信號卷繞和RF采樣延遲的圖。
圖7是表示本發(fā)明的引導收斂的方法的圖。
圖8是表示超聲波RF信號的圖。
圖9是表示使用了自相關的位移估計的圖。
圖10是表示在仿真中使用了本發(fā)明的位移估計結果的例子的圖。圖11是表示在體模實驗中使用了本發(fā)明的位移估計結果的例子的圖。圖12是表示下層構造的位移和對RF信號的影響的圖。圖13是本裝置的框圖。
具體實施方式
以下，參照附圖進行詳細說明。實施方式的位移估計方法是使用超聲波信號(圖13的超聲波信號1201s)反復估計位移(參照圖12的位移1207)的位移估計方法(參照圖13的位移估計裝置IX)。S卩，例如通過進行多次位移大小的確定，進行多次確定處理，能夠確定精度比較好的位移大小。并且，掃描至少1個超聲波信號并發(fā)送(參照圖13的超聲波處理部1 中的發(fā)送部1X1)給介質(zhì)(參照介質(zhì)1200M(圖12)以及介質(zhì)1200M所具有的下層構造(位移測定的對象)1203x(圖 12、13))。即，例如向線方向(圖8的方向81L)的多個位置發(fā)送超聲波信號，進行這些多個位置的掃描。并且，接收(參照接收波1X2)所掃描的從所述介質(zhì)反射的所述超聲波信號(超聲波信號1201s)。另外，超聲波處理部1 例如是探頭(探針)等。并且，計算(參照圖13的信息處理部l)(b所包含的尺寸計算部1X3)窗尺寸(圖 13的窗尺寸1X3M 例如圖3的(a)欄、(b)欄的窗尺寸(winSize)等)。S卩，決定各個循環(huán)中的窗尺寸。另外，尺寸計算部1X3例如可以包含窗計算部102(圖1)，也可以是窗計算部102 的至少一部分等。S卩，例如根據(jù)第n-1循環(huán)中的相關值(圖13的相關值1X7M)來決定第η循環(huán)中的窗尺寸(η彡2)。并且，根據(jù)所計算出的所述窗尺寸，計算(參照邊界計算部1X4)窗的邊界(參照邊界的信息1Χ4Μ 圖4的(a)欄、(b)欄所示的窗的始端(winStart)以及終端(winEnd)
寸乂 O另外，邊界計算部1X4例如可以是窗計算部102(圖1)的至少一部分。S卩，確定始端和終端，這些始端和終端間的寬度是所計算的窗尺寸。并且，使用基于所計算的所述邊界的所述窗，估計(參照估計部1X5)所述超聲波信號的各個深度的位移(位移UbM)。另外，估計部1)(5例如可以是位移估計部101(圖1)的至少一部分等。S卩，例如確定與本裝置的距離相互不同的多個位置上的位移。并且，根據(jù)所估計的所述位移，卷繞所述超聲波信號(參照卷繞部1X6)。另外，卷繞部1)(6例如可以是RF信號卷繞部100(圖1)的至少一部分等。并且，使用所卷繞的所述超聲波信號來引導(參照收斂控制部1X7)該位移估計方法的收斂，使得所述超聲波信號的相關值(相關值1X7M)增大。
8
另外，收斂控制部1X7例如可以是收斂引導部103的至少一部分等。S卩，將在根據(jù)所卷繞的超聲波信號計算出的相關值最大的循環(huán)中估計的位移大小確定為精度最高的大小。并且，使所計算的所述窗尺寸變化，使得用于連續(xù)多個循環(huán)的所述窗尺寸逐漸減S卩，例如在各個循環(huán)中計算出的窗尺寸是比在該循環(huán)之前的循環(huán)中計算出的窗尺寸小的尺寸，所計算的窗尺寸向小尺寸變更。并且，計算所述窗尺寸，使得所有(深度)的窗的信號能量相等。并且，超聲波信號的相關值(第n-1循環(huán)中的相關值)越大，計算越小的所述窗尺寸(第η循環(huán)中的窗尺寸)。這里，相關值等各循環(huán)中的收斂程度的評價值和窗尺寸的關系可以參照基于活體上的實驗值的預定關系式或表示兩者的對應關系的表等來決定。并且，將與第n-1循環(huán)中的相關值對應了所決定的對應關系的窗尺寸確定為第η 循環(huán)中的窗尺寸。S卩，例如在本裝置中進行下面的動作。即，可以向作為測定對象物(例如伴隨著新生血管的惡性腫瘤等)的下層構造 1203χ(圖12)發(fā)送超聲波信號1201s，并接收所發(fā)送的超聲波信號1201s。所確定的延遲時間1208是第1脈沖位置(參照圖12的(b)欄的RF信號1204) 和第2脈沖位置(參照RF信號1205)間的延遲時間。這里，第1脈沖是所接收的超聲波信號1201s中下層構造1203在進行位移1207 之前的脈沖。并且，第2脈沖是測定對象物進行了位移1207后的脈沖。并且，根據(jù)所確定的延遲時間1208來決定位移1207。并且，可以判定所確定的位移1207的大小的下層構造1203x的動作是否是預先設定的動作。并且，可以在判定為是預先設定的動作的情況下，判定為該下層構造1203x是惡性腫瘤(癌)，在判定為不是該動作的情況下，判定為是良性腫瘤或正常組織。另外，例如可以通過收斂控制部1χ7(圖1 來進行這種是否是癌的判定。在這種本裝置中，更具體來講，例如可以如在后面詳細說明那樣，在第n-1循環(huán)中確定的相關值越大，在第η循環(huán)中，越能進行更大窗尺寸的處理。另外，這里，基于窗尺寸的處理例如可以是從所接收的超聲波信號1201s中的該窗尺寸的窗部分的數(shù)據(jù)中計算出第η循環(huán)中的相關值的處理。根據(jù)本處理，由于在各循環(huán)中能夠設定與超聲波的接收信號的信號電平對應的適當窗，所以能夠提高收斂的精度，并且以較少的反復次數(shù)獲得收斂，并且能夠以較少的處理量實現(xiàn)高精度的位移估計。以下的實施方式僅說明各發(fā)明步驟的原理。這里說明的具體例子的各種變形，對于本領域技術人員來說是顯而易見的。因此，本發(fā)明的技術范圍并不限定于本說明書記載的具體例示的內(nèi)容，只限定于權利要求書的范圍。圖1是表示位移估計方法的圖。
9
圖1示出本發(fā)明的主要實施方式。RFSig (d，1，f)表示從超聲波部獲得的RF信號，d表示深度方向(圖8的方向81D)， 1表示線方向(方向81L)，f表示幀方向(方向81F)。圖8進一步示出與這些相關的內(nèi)容。這種RF信號的1幀內(nèi)的1條線(參照圖示的801)通過由超聲波變換器800在該線位置向掃描方向(方向802)發(fā)送脈沖的機構生成。這里，該脈沖一邊反射和散射，一邊沿著進路與下層構造相互作用，該超聲波變換器800接收反射和散射的信號。并且，該信號被超聲波部變換為對應的線。通過相同的機構，在各個位置生成多條線，并在各種時間實例(time instance)生成多個幀。圖12的(a)欄表示下層構造的位移，(b)欄表示對RF信號的影響。圖12示出所掃描的介質(zhì)的下層構造的位移被作為RF信號的延遲而反射的機構。以例示為目的，假設發(fā)送1條線的1個脈沖，下層構造只有1個反射邊界。在第1時間實例(參照(a)欄的左側部分)，超聲波變換器1200向所掃描的介質(zhì)發(fā)送脈沖(參照脈沖1201)。該脈沖通過介質(zhì)傳播到到達下層構造1203(1203χ)的反射邊界為止。通過該邊界，脈沖向變換器反射(參照所反射的脈沖1202)。變換器接收該脈沖，變換為RF信號1204(參照圖12的(b)欄的表的第1行)。RF信號1204的整體的脈沖位置表示被發(fā)送的脈沖1201向下層構造1203前進并被反射回來所需的時間。在第2時間實例(圖12的(a)欄的右側)，RF信號1205 (圖12的(b)欄的表的第3行)由同一機構生成，但下層構造位移到新的位置(參照位移后的下層構造1206)。由于位移x(位移1207)，所發(fā)送的脈沖(脈沖1201)被傳播并被下層構造的邊界反射所需的時間變長，其結果，引起RF信號的延遲(參照延遲時間1208)。另外，在該例子中，位移χ的朝向是圖12所示的方向1203d的朝向。如果知道了該延遲時間1208的延遲值，就能夠估計出位移(位移1207)。并且，在具有多個下層構造的介質(zhì)中，能夠生成多個反射脈沖。這些多個下層構造的位移的量相互不同的原因在于RF信號中的反射脈沖的延遲不同。因為這些脈沖重復并相互作用，所以很難估計各下層構造的準確位移。本發(fā)明提示出通過反復進行估計來準確估計下層構造的位移的方法。根據(jù)圖1所示的主要實施方式，本發(fā)明具備以下的主要模塊。即，在本裝置中，具備RF信號卷繞部100 (參照圖5)、位移估計部101 (參照圖9)、 窗計算部102(參照圖2)、收斂引導部103(參照圖7)、輸出開關104。RF信號卷繞部100根據(jù)上次的估計循環(huán)(第n-1循環(huán))的位移估計結果 dispR0imd(d，l，round)，進行對象估計循環(huán)(第η循環(huán))中的卷繞。在該卷繞中，對于對象估計循環(huán)，卷繞RF信號的所選擇的幀內(nèi)的所選擇的線。另外，其在第1循環(huán)中被設定為0。
該模塊的輸出是被用于估計循環(huán)的被卷繞的RF信號RFSigWarp(d，1，f，round)。另外，在第1估計循環(huán)(即，round = 1)中，通過將dispRound(d, 1,0)初始化為 0，使得RFSigffarp (d, 1，f，1)與從超聲波部獲得的RF信號RFSig (d, 1，f)(前述)實質(zhì)上是相同的。由于估計位移，所以可以將預先定義的幀的集合選擇為針對該模塊的輸入。以例示為目的，選擇2個巾貞，并用f 1和f2來表示它們。RF信號卷繞的目的是根據(jù)在先行的循環(huán)中所估計的位移將1個RF信號校正為與其他RF信號相適合。并且，之后通過從所校正的RF信號中估計殘留位移，并追加到先行循環(huán)的位移上，生成新的位移。一些數(shù)量(預先設定的數(shù)量)的估計循環(huán)之后，所估計的位移與所校正的RF信號幾乎一致，將殘留位移收斂成0。位移估計部101通過根據(jù)RF信號的時間延遲進行導出從而根據(jù)RF信號 RFSigffarp (d, 1，f, round)中估計位移。優(yōu)選執(zhí)行該任務的方法是通過winStart (d，1)和WinEnd(d，1)來記述各深度的估計窗的自相關的方法，但并不限于此。窗計算部102計算各深度的用于所選擇的線的估計窗參數(shù)wir^tart(d，1)和 winEnd(d,1)。窗計算部102獲取RF信號RFSigWarp (d，1，f，round)，并且，獲取對象循環(huán)數(shù) curRound和RF信號差分residue(d，1，round)作為來自收斂引導部103的輸入。收斂引導部103進行本技術的引導收斂的處理。即，為了決定所估計的位移的質(zhì)量，計算所卷繞的RF信號間的差分。將該差分最小(包含對象循環(huán))的循環(huán)判定為表示最準確的位移估計結果。該模塊的輸出是對象估計循環(huán)數(shù)curRound、RF信號差分最小的循環(huán)數(shù) minResidueRound(d, 1)、作為所估計的位移質(zhì)量的測定值的RF信號差分residue (d，1， round)。residue (d，1，round)越小，表示質(zhì)量越好，在估計循環(huán)整體中，如果residue (d， 1, round)減少，則表示收斂。另外，作為計算所卷繞的RF信號間的差分的方法，考慮了“平方差分和”和“絕對差分”的方法，但不限于此。輸出開關104 獲取 curRound 和 minResidueRound (d，1)作為輸入。在對象循環(huán)數(shù)到達最大值的情況下，選擇最準確的估計結果作為最終輸出。否則，為了用于后續(xù)的估計循環(huán)，將該最準確的估計結果選擇為RF信號卷繞部 100的輸入，將residue (d，1，round)和curRound作為輸入發(fā)送給位移估計部101。下面的段落示出了一些基于圖1所示的主要實施方式的具體實施方式
。圖5是表示RF信號卷繞方法的圖。作為本發(fā)明的具體實施方式
，主要的實施方式的RF信號卷繞部100通過圖5所示的方法來實現(xiàn)。位移的結果是RF信號被移動并且被解壓縮。
因此，如下所示，在延遲值計算部500中根據(jù)先行循環(huán)的最準確的位移估計結果 dispRound(d, 1, minResidueRound(d, 1))計算使 RFSig(d，1，fl)適合于 RFSig(d，1，f2) 的采樣中的延遲值(參照圖12的延遲時間1208)。算式1這里，fs是RF信號的采樣頻率，c是所掃描的介質(zhì)中聲速。作為結果，將所掃描的介質(zhì)中的位移反映成所獲得的RF信號中的延遲。由于各深度的延遲值是已知的，所以能夠按照使RFSig(d，1，fl)或RFSig(d，1， f2)的任意一個中的各線的(與所有深度對應的)所有采樣適合于其他采樣，來進行延遲。這在RF采樣延遲部501中執(zhí)行。在1個實施方式中，使用信號插值來執(zhí)行RF采樣延遲部501中的延遲。在其他實施方式中，使用分數(shù)延遲濾波器來執(zhí)行RF采樣延遲部501中的延遲。圖6是表示RF信號卷繞和RF采樣延遲的圖。圖6示出RF信號卷繞處理。上述的圖12示出了由于位移而引起了 RF信號的延遲。RF信號卷繞的目的是通過使1個RF信號向與所估計的位移相反的方向延遲而使該效果無效。作為該無效化的結果，期待所卷繞的RF信號相互適合。延遲值計算部500所計算出的delayValueS(d，l)是表示1個RF信號內(nèi)的各采樣為了達成該適合而延遲的值。使用2個RF信號600和601 (圖6)來進行該例示。數(shù)據(jù)602表示在某個估計循環(huán)之后由延遲值計算部500計算出的delayValues (d， 1)?？芍猟elayValues (d，1)的特定深度的值表示用于使RF信號600中的相同深度的采樣適合于延遲值計算部601的延遲所需的量。然而，delayValues (d, 1)中的延遲值并不始終是數(shù)。因此，由于按照delayValues (d，l)能夠延遲RF信號600中的各采樣，所以舉出信號插值或分數(shù)延遲濾波器作為兩個選擇候選。 通過將位移估計部101應用于RF信號，獲得用于第1循環(huán)的位移和用于后續(xù)循環(huán)的殘留位移，為了獲得最終的位移，對它們進行組合。圖9示出優(yōu)選的實施方式，但并不限于此。圖9表示使用了自相關的位移估計。IQ 解調(diào)部 900 將 RF 信號 RFSigffarp (d, 1，f, round)變換為基帶信號 IQSig (d，1， f)。自相關計算部901根據(jù)以下的式子計算autocorHd，1)。算式2
wmEnd ( d J )autocorr(d,l)= J] IQSigfd,Ii.,/Jx conj(IQSig(d,I；,/,))
Ii=WinSTarTl ci J)
這里，conj ()表示復共軛。最后，根據(jù)以下的式子，由位移計算部902將autocorr (d，1)變換為位移。算式3^.-,/； ir^v nr/” u J ”dispRound(d，1，round) = dispRound(d，1，round-1)+residueDisp這里，arg()是計算復數(shù)的輻角的函數(shù)。根據(jù)從1開始的循環(huán)(表示第1循環(huán))，將dispRouncKd，1,0)初始化為0。圖2是表示窗計算方法的圖。通過窗計算部102針對各估計窗中的各線的各深度，決定不同的估計窗。這在圖2中示出。圖3是表示窗尺寸計算的具體例的圖。窗尺寸計算部200(圖2)決定針對各估計循環(huán)的窗尺寸。該模塊的輸出是winSize (d，1)。在優(yōu)選實施方式中，針對各循環(huán)的窗尺寸隨著循環(huán)的函數(shù)的減少而變化。這在圖3中示出。在具有窗尺寸計算部200(圖2)的實施方式中，如圖3的(a)欄所示，窗尺寸由于減少函數(shù)winSizeRound (round)而直接關系到估計循環(huán)(round)數(shù)。關于由curRoimd所確定的對象估計循環(huán)，根據(jù)以下，確定針對所有深度的窗尺寸。關于所有的 d 禾口 1，winSize (d，1) = winSizeRound (curRound)?！肮潭ù胺桨浮笔窃谟糜谖灰乒烙嫷乃猩疃冗x擇成為估計窗尺寸的winSiZe(d， 1)的相同值的處理?！肮潭ù胺桨浮钡膬?yōu)點是適合于所需的運算量不多(比較少)并且RF信號能量沿著深度方向被規(guī)則準確地分配的情況。在窗尺寸計算部200的其他實施方式中，各深度的窗尺寸根據(jù)作為上限的 maxWin (round)、作為下限的 minWin (round)、RF 信號功率 sigPow(d，1)來計算。maxffin (round)禾口 minWin (round)是循環(huán)的減少函數(shù)。這作為“瞬間功率窗方案”，在圖3的(b)欄中示出。應用以下的條件。算式4winSize(d,1) e [min Win(curRound), max Win(curRound)]并且，如圖3 的(b)欄所示，winSize (d，l)也可以與 sigPow(d，1)、即與 sigPow = 0對應的最大窗尺寸和與最大sigPow對應的最小窗尺寸成反比?！八查g功率窗方案”與“固定窗方案”相比，所需的運算量多，但是，特別在RF信號能量沿著深度方向相當不均勻的情況下，精度更高。為了克服信號對噪聲比的較低程度，對具有較低信號功率的區(qū)域指定更大的估計窗，為了更詳細地分布位移，對具有較大信號功率的區(qū)域指定更小的估計窗。然而，關于規(guī)則準確地分配了信號能量的RF信號，“固定窗方案”和“瞬間功率窗
13方案”沒有大的差異。在窗尺寸計算部200的其他實施方式中，也可以計算所有深度的窗尺寸，使得所有的窗內(nèi)的信號能量保持同量，這被稱為“恒定能量窗方案”。信號功率計算部201(圖2)計算沿著幀方向的RF信號的平均功率。信號功率計算部201獲取所卷繞的RF信號RFSigWarp (d，1，f，round)和對象循環(huán)數(shù)curRound作為輸入。在RFSigWarp (d，1，f，curRound)內(nèi)具有多個幀，按照如下方式計算功率。算式5
權利要求
1.一種位移估計方法，使用超聲波信號反復估計位移，其中， 該位移估計方法包含如下步驟掃描至少1個超聲波信號，并發(fā)送給介質(zhì)； 接收所掃描的從所述介質(zhì)反射的所述超聲波信號； 計算窗尺寸；根據(jù)所計算的所述窗尺寸，計算窗的邊界；使用基于所計算的所述邊界的所述窗，估計所述超聲波信號的各個深度的位移； 根據(jù)所估計的所述位移，對所述超聲波信號進行卷繞；以及使用所卷繞的所述超聲波信號，引導該位移估計方法的收斂，使得所述超聲波信號的相關值變大。
2.根據(jù)權利要求1所述的位移估計方法，其中，改變所計算的所述窗尺寸，使得連續(xù)的多個循環(huán)中使用的所述窗尺寸逐漸減小。
3.根據(jù)權利要求1所述的位移估計方法，其中， 該位移估計方法還包含如下步驟使連續(xù)的多個所述循環(huán)的所述窗尺寸的上限和下限分別按照逐漸減小的方式變化；以及根據(jù)各深度的信號功率，計算該深度的上限和下限間的窗尺寸，作為各個深度的所述窗尺寸。
4.根據(jù)權利要求1所述的位移估計方法，其中， 計算所述窗尺寸，使得所有窗中的信號能量相等。
5.根據(jù)權利要求1所述的位移估計方法，其中， 超聲波信號的相關值越大，計算出越小的所述窗尺寸。
6.根據(jù)權利要求1所述的位移估計方法，其中，所述邊界向該深度的兩側延伸，使得基于該邊界的所述窗成為以對應的深度為中心的對稱的窗。
7.根據(jù)權利要求1所述的位移估計方法，其中，所述邊界向?qū)乃錾疃鹊膬蓚妊由?，使得對應深度的基于該邊界的所述窗的兩側的信號能量相互相等?br> 8.根據(jù)權利要求1所述的位移估計方法，其中， 該位移估計方法還包含如下步驟根據(jù)所計算的延遲值使所述超聲波信號內(nèi)的各采樣延遲，由此卷繞所述超聲波信號。
9.根據(jù)權利要求1所述的位移估計方法，其中， 該位移估計方法還包含如下步驟在各循環(huán)之后，計算所卷繞的RF信號間的信號差分；以及決定所計算的所述信號差分為最小的循環(huán)。
10.一種位移估計裝置，使用超聲波信號反復估計位移，其中， 該位移估計裝置包含發(fā)送部，掃描至少1個超聲波信號，并發(fā)送給介質(zhì)； 接收部，接收所掃描的從所述介質(zhì)反射的所述超聲波信號；尺寸計算部，計算窗尺寸；邊界計算部，根據(jù)所計算的所述窗尺寸，計算窗的邊界；估計部，使用基于所計算的所述邊界的所述窗，估計所述超聲波信號的各個深度的位移；卷繞部，根據(jù)所估計的所述位移，對所述超聲波信號進行卷繞；以及收斂控制部，使用所卷繞的所述超聲波信號，引導該位移估計裝置執(zhí)行的位移估計方法的收斂，使得所述超聲波信號的相關值變大。
全文摘要
本發(fā)明公開一種位移估計方法，使用超聲波信號反復估計位移，其中，包含如下步驟掃描至少1個超聲波信號，并發(fā)送給介質(zhì)；接收所掃描的從所述介質(zhì)反射的所述超聲波信號；計算窗尺寸；根據(jù)所計算的所述窗尺寸，計算窗的邊界；使用所計算的所述窗，估計所述超聲波信號的各個深度的位移；根據(jù)所估計的所述位移，對所述超聲波信號進行卷繞；以及使用所卷繞的所述超聲波信號，引導該位移估計方法的收斂，使得所述超聲波信號的相關值變大。
文檔編號G01N29/06GK102361597SQ201180001460
公開日2012年2月22日 申請日期2011年1月19日 優(yōu)先權日2010年1月20日
發(fā)明者張國成, 范樹峰, 陳英俊 申請人:松下電器產(chǎn)業(yè)株式會社