專利名稱：采用兩種尺度的生物組織位移估計方法
技術領域：
本發(fā)明屬于超聲彈性成像技術領域，特別涉及生物組織位移估計方法。
背景技術：
生物組織彈性模量的變化通常與其病理現(xiàn)象有關。例如，惡性的病理損害，例如乳房硬癌、前列腺癌、甲狀腺癌及肝轉(zhuǎn)移等，通常表現(xiàn)為硬的小結(jié)。乳房硬癌是乳腺癌的最常見形式，大約占乳腺癌總數(shù)的四分之三，由于其基質(zhì)密度增大而表現(xiàn)為致密的硬塊。而其他類型的乳腺癌如導管內(nèi)癌和乳頭狀瘤則表現(xiàn)為柔軟的組織，良性的乳腺纖維囊性病也很少表現(xiàn)為硬塊。
生物組織的彈性模量信息對于疾病的診斷過程具有重要的參考價值。然而，包括X射線成像、超聲成像、計算機斷層成像(CT)和磁共振成像(MRI)等在內(nèi)的傳統(tǒng)醫(yī)學成像模態(tài)都不能直接提供關于彈性模量這一組織的基本力學屬性的信息。1991年，J.Ophir提出超聲彈性成像(ultrasound elastography)的方法，對組織的彈性模量分布進行定量估計、成像。目前，超聲彈性模量已經(jīng)成為醫(yī)學超聲成像的一個研究熱點，廣泛應用于乳房、前列腺、動脈粥樣斑塊、心肌動力學以及高強度聚焦超聲與射頻消融引起的損害(lesion)的檢測與評估。
超聲彈性成像的基本原理為將超聲探頭嵌于一塊擠壓平板中，沿著探頭的縱向壓縮組織，分別采集組織壓縮前、后的射頻信號；組織被壓縮時，組織內(nèi)將會產(chǎn)生一個沿壓縮方向的應變，如果組織內(nèi)部彈性模量分布不均勻，組織內(nèi)的應變分布也會有所差異；彈性模量較大的區(qū)域，引起的應變比較小；反之，彈性模量較小的區(qū)域，相應的應變比較大。通過一些方法估計出組織內(nèi)部不同位置的位移，從而計算出組織內(nèi)部的應變分布情況，用來間接描述組織內(nèi)部的彈性模量分布，從而描述組織的生理、病理狀態(tài)。
對于二維超聲彈性成像，一般采用線陣的B型超聲探頭，采集組織壓縮前、后的探頭每一條掃描線的射頻信號，分別進行上面描述的位移估計，從而計算出每一條掃描線對應組織的一維應變分布。最后把所有掃描線對應的一維應變分布按掃描線順序組成一個二維應變分布，以灰度圖或者偽彩圖的形式表示，用來間接描述組織內(nèi)部的彈性模量分布。
一般的超聲彈性成像方法包括以下步驟1.利用商用B型超聲儀器(一般采用線陣探頭)得到待測生物組織(一般為人體組織，也可以為動物組織，以下簡稱組織)壓縮前的一幅數(shù)字化的二維射頻信號(可以采用模擬射頻信號輸出端接信號放大器，再接高速數(shù)據(jù)采集卡，得到數(shù)字化的二維射頻信號；也可以在數(shù)字化B型超聲儀器上直接得到數(shù)字化的二維射頻信號)；2.手持該B型超聲儀器的探頭或者利用步進電機或者螺旋裝置驅(qū)動該探頭，沿著探頭的縱向?qū)υ摻M織施加一個微小的擠壓(對組織施加的壓縮量一般控制在為1％的數(shù)量級)，得到組織壓縮后的一幅數(shù)字化的二維射頻信號；3.從步驟1和2的得到的組織壓縮前、后的二維射頻信號中分別取出第一條掃描線的數(shù)據(jù)，設為s1(n)和s2(n)，n表示該兩條掃描線上的數(shù)據(jù)序號，1≤n≤nmax，n的最大值nmax由該B型超聲儀器的探查深度、發(fā)射的超聲波在組織中的傳播速度以及射頻信號的采樣頻率決定；4.從該掃描線數(shù)據(jù)s1(n)中取一小段長度為T的數(shù)據(jù)d1，其數(shù)據(jù)個數(shù)為U，U＝round(T×U1)，其中，T的單位為mm，U1代表1mm的組織對應的數(shù)據(jù)個數(shù)，由發(fā)射的超聲波在組織中的傳播速度以及射頻信號的采樣頻率決定，round(·)代表四舍五入的取整操作，該數(shù)據(jù)d1的序號從n1到n1+U-1，n1可在1≤n1≤U的范圍內(nèi)選擇；在τ1到τ2確定的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R(τ)，計算公式如下R(τ)=Σi=n1n1+U-1s1(i)s2(i-τ)Σi=n1n1+U-1s12(i)·Σi=n1n1+U-1s22(i-τ)(τ1≤τ≤τ2)]]>其中i為計算過程中表示數(shù)據(jù)序號的循環(huán)變量，τ1為0，τ2為對組織施加的壓縮量，以采樣數(shù)據(jù)的個數(shù)表示；為了提高位移估計的精度，一般還需要對計算得到互相關函數(shù)進行插值，如拋物線等常規(guī)插值方法；5.確定步驟4得到的互相關函數(shù)R(τ)的最大值對應的位置t1t1就是數(shù)據(jù)d1在組織壓縮后的位移(即s1(n)中的序號從n1到n1+U-1的小段數(shù)據(jù)d1的在組織壓縮后移動到s2(n)中的序號從n1-t1到n1+U-1-t1的位置)；6.依次從掃描線數(shù)據(jù)s1(n)中取一小段長度為T即數(shù)據(jù)個數(shù)為U的數(shù)據(jù)d2、d3、…、dN，每段數(shù)據(jù)的序號依次錯開V個采樣數(shù)據(jù)，1≤V≤U，直到再錯開V個采樣數(shù)據(jù)將超出s1(n)的范圍，按步驟4、5相同的方法依次得到各段數(shù)據(jù)對應的位移t2、t3、…、tN，其中N為小段數(shù)據(jù)的總數(shù)；則位移序列t1、t2、…、tN為第一條掃描線數(shù)據(jù)s1(n)對應的組織的位移估計；7.利用與步驟3-6相同的方法，依次得到第2、3、…、M條掃描線數(shù)據(jù)對應的組織的位移估計，其中M為表示探頭的掃描線總數(shù)，由探頭決定；8.對第一條掃描線數(shù)據(jù)s1(n)對應的組織的位移估計序列t1、t2、…、tN求差分，得到組織第一條掃描線s1(n)對應組織的應變分布，計算公式如下，ϵ1=t2-t1V,ϵ2=t3-t2V,···,ϵN-1=tN-tN-1V]]>其中，ε1、ε2、…、εN-1分別為d1、d2、…、dN-1對應的組織應變；9.利用與步驟8相同的方法，依次得到組織第2、3、…、M條掃描線數(shù)據(jù)對應的組織的應變分布；10.將步驟9得到的M條掃描線數(shù)據(jù)對應的應變分布，按照掃描線的順序組合成一個二維數(shù)據(jù)，并以灰度圖或者偽彩圖的形式表示出來，就得到組織的二維應變分布圖。
在超聲彈性成像中，關鍵的問題在于對組織的位移分布進行估計，也就是上面描述的方法的步驟3-7。互相關函數(shù)的值越大，說明壓縮前、后的小段數(shù)據(jù)吻合得越好，互相關函數(shù)的最大值位置代表了壓縮前的小段數(shù)據(jù)在壓縮后對應的位置，從而可以求出該小段數(shù)據(jù)的位移，也就是該小段數(shù)據(jù)對應的組織的位移。
在超聲彈性成像中，選擇的壓縮前信號中的小段數(shù)據(jù)用來跟蹤對應的小段組織的位移，稱為跟蹤波段。其長度稱為跟蹤波段長度，或者尺度。選取合適的尺度在超聲彈性成像的位移估計中非常重要。當對組織施加的壓縮量比較小(如小于1％)的時候，尺度越大，包含的信息越多，位移估計的精度越高，對隨機噪聲的干擾也越不敏感。當對組織施加的壓縮量比較大(如1-5％)的時候，尺度越大，由于組織被壓縮，組織壓縮前、后信號重合的部分越小，利用互相關函數(shù)最大值位置進行位移估計的精度越低，同時，尺度越大，又越不容易受到隨機噪聲的干擾；而當尺度較小的時候，組織壓縮前、后信號重合的部分較大，位移估計的精度較高，然而此時的位移估計又對隨機噪聲的干擾比較敏感。
為了增加組織壓縮前、后信號重合的程度，一般采用線性插值的方法，把被“壓縮”的信號“拉伸”(stretching)成與原信號同長，從而增加壓縮前后波形的重合度。在波形拉伸中，拉伸系數(shù)與對組織施加的壓縮量或應變(以百分比表示)有關。當拉伸系數(shù)接近組織的應變時，波形的重合度最大，從而互相關函數(shù)最大值最大。但是組織的應變在時延估計前是未知的，拉伸系數(shù)的選擇只能根據(jù)施加的位移量和組織的深度作大致判斷。并且由于組織的不均勻性，組織內(nèi)部應變的分布也不均勻，因此特定的拉伸系數(shù)在組織的某些部分可能拉伸不足，在某些部分又拉伸過量，這反過來又引入誤差。因此，有人提出一種對不同位置的波段進行自適應拉伸系數(shù)的方法，使得拉伸后的波段與壓縮前波段的互相關函數(shù)最大值最大，然后利用拉伸系數(shù)直接估計出該位置組織的應變。但是這種自適應拉伸的方法計算時間較長，在實際應用中可能受到限制。
從實際應用的角度，希望對組織施加的壓縮量大一些好，這樣可以更能“看清”組織內(nèi)部彈性模量差異引起的應變分布差異的細節(jié)。一般采用多次壓縮(multi-compression)的方法，將大位移分解為小位移之和，逐次對組織施壓一個小的位移量，計算出相應的應變分布。最后將多幅應變分布疊加，得到大位移時的應變分布。然而，該方法增加了數(shù)據(jù)采集時間。

發(fā)明內(nèi)容
本發(fā)明為了解決對組織施加的壓縮量比較大的時候大尺度和小尺度的矛盾，提出一種采用兩種尺度的組織位移估計方法。在組織壓縮量比較大時，將選擇大尺度或小尺度進行組織位移估計的兩種方法結(jié)合起來，可對這兩種組織位移估計的優(yōu)缺點實現(xiàn)揚長避短，提高組織位移估計的精度。
本發(fā)明提出的一種采用兩種尺度的組織位移估計方法，包括以下步驟1.從組織壓縮前、后的二維射頻信號中分別取出第一條掃描線的數(shù)據(jù)，設為s1(n)和s2(n)，n表示該兩條掃描線上的數(shù)據(jù)序號，1≤n≤nmax，n的最大值nmax由該B型超聲儀器的探查深度、發(fā)射的超聲波在組織中的傳播速度以及射頻信號的采樣頻率決定；2.從該掃描線數(shù)據(jù)s1(n)中取一小段尺度為Ta的數(shù)據(jù)d1，其數(shù)據(jù)個數(shù)為Ua，Ua＝round(Ta×U1)，其中，Ta的單位為mm，U1代表1mm的組織對應的數(shù)據(jù)個數(shù)，由發(fā)射的超聲波在組織中的傳播速度以及射頻信號的采樣頻率決定，round(·)代表四舍五入的取整操作，該數(shù)據(jù)d1的序號從n1到n1+Ua-1，n1可在1≤n1≤Ua的范圍內(nèi)選擇；在τ1到τ2的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R(τ)，計算公式如下R′(τ)=Σi=n1n1+U-1s1(i)s2(i-τ)Σi=n1n1+U-1s12(i)·Σi=n1n1+U-1s22(i-τ)(τ1≤τ≤τ2)]]>其中i為計算過程中表示數(shù)據(jù)序號的循環(huán)變量，τ1為0，τ2為對組織施加的壓縮量，以采樣數(shù)據(jù)的個數(shù)表示；為了提高位移估計的精度，最好還需要對計算得到互相關函數(shù)進行插值處理，如拋物線插值等常規(guī)插值方法；3.確定步驟2得到的互相關函數(shù)R′(τ)的最大值對應的位置t′1，t′1為利用尺度Ta進行位移估計時數(shù)據(jù)d1在組織壓縮后的位移(即s1(n)中的序號從n1到n1+Ua-1的小段數(shù)據(jù)d1的在組織壓縮后移動到s2(n)中的序號從n1-t′1到n1+Ua-1-t′1的位置)；4.從該掃描線數(shù)據(jù)s1(n)中的數(shù)據(jù)段d1的中央位置再取尺度為Tb、序號從n1+round(Ua-Ub2)]]>到n1+round(Ua-Ub2)+Ub-1]]>的一小段數(shù)據(jù)e1，該數(shù)據(jù)個數(shù)為Ub，Ub＝round(Tb×U1)，其中，Tb的單位為mm，并且滿足Tb＜Ta，最好取116Ta≤Tb≤12Ta;]]>在τ3到τ4確定的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R″(τ)，計算公式如下R′′(τ)=Σi=m1m1+Ub-1s1(i)s2(i-t′1-τ)Σi=m1m1+Ub-1s12(i)·Σi=m1m1+Ub-1s22(i-t′1-τ)(τ3≤τ≤τ4)]]>其中τ3和τ4為組織位移精細調(diào)整的范圍，滿足如下關系τ4＜τ2，τ3＝-τ4，最好取116τ2≤τ4≤12τ2;]]>為了提高位移估計的精度，最好還對計算得到互相關函數(shù)進行插值，如拋物線插值等常規(guī)插值方法；
5.確定步驟4得到的互相關函數(shù)R″(τ)的最大值對應的位置t″1，t″1為利用尺度Tb進行位移估計時數(shù)據(jù)e1的位移相對于t′1的偏移；6.利用尺度Tb估計得到的位移偏移t″1對尺度Ta估計得到的位移t′1進行調(diào)整，得到的小段數(shù)據(jù)d1的精細調(diào)整后的位移t1，即t1＝t′1+t″17.依次從掃描線數(shù)據(jù)s1(n)中取一小段尺度為Ta即數(shù)據(jù)個數(shù)為Ua的數(shù)據(jù)d2、d3、…、dN，每段數(shù)據(jù)的序號依次錯開V個采樣數(shù)據(jù)，1≤V≤Ua，直到再錯開V個采樣數(shù)據(jù)將超出s1(n)的范圍，按步驟2-6相同的方法依次得到各段數(shù)據(jù)對應的位移t2、t3、…、tN，其中N為小段數(shù)據(jù)的總數(shù)；則位移序列t1、t2、…、tN為第一條掃描線數(shù)據(jù)s1(n)對應的組織的位移估計；8.利用與步驟1-7相同的方法，依次得到第2、3、…、M條掃描線數(shù)據(jù)對應的組織的位移估計，其中M為表示探頭的掃描線總數(shù)，由探頭決定。
本發(fā)明的原理本發(fā)明采用兩種尺度的跟蹤波段進行位移估計，先用大尺度的跟蹤波段Ta進行位移的初估計，然后用小尺度的跟蹤波段Tb進行細調(diào)。這種采用兩種尺度的方法的思想模仿了人眼觀察事物的方式先用大尺度粗看，再用小尺度細看，也就是先整體后局部的思想。
本發(fā)明利用大的尺度來消除隨機噪聲的影響，從而確定了組織位移的大致范圍，然后利用小的尺度來對組織位移進行精確估計。因為利用小尺度的跟蹤波段進行位移估計的時候，是在大尺度跟蹤波段的估計結(jié)果的基礎上進行，可以在一個小的范圍內(nèi)計算互相關函數(shù)，所以隨機噪聲的影響比較小。這種方法打破了尺度固定不變的思路，很好地解決了對組織施加的壓縮量較大的時候大尺度和小尺度的矛盾，提高了組織位移估計的精度，又不容易受隨機噪聲的影響。
從誤差分析的角度來看，位移估計的誤差可以分成相位誤差和偽峰誤差。相位誤差表現(xiàn)為壓縮前、后射頻掃描線數(shù)據(jù)之間的互相關函數(shù)的最大值位置在對應組織的實際位移值的附近擾動，誤差較小。而偽峰誤差表現(xiàn)為該互相關函數(shù)的最大值位置已經(jīng)嚴重偏離對應組織的實際位移值，此時的組織位移估計已經(jīng)不包含足夠的信息量。當對組織施加的壓縮量比較大的時候，尺度越大，越不容易受隨機噪聲的干擾，越不容易出現(xiàn)偽峰誤差，但是由于組織壓縮前、后信號重合的部分越小，此時的相位誤差較大；而尺度較小時，相位誤差得到減小，但是此時的偽峰誤差較大。
本發(fā)明提出的采用兩種尺度的方法，就是先利用大尺度來減小偽峰誤差，然后利用小尺度來減小相位誤差，從而綜合了兩種尺度的優(yōu)點，減小了誤差的干擾，提高了組織位移估計的精度。
本發(fā)明的特點1采用兩種不同尺度的跟蹤波段，利用先粗略估計再精細調(diào)整、先整理后局部的思路進行位移估計；2先用大尺度的跟蹤波段進行位移的粗略估計，然后用小尺度的跟蹤波段進行精細調(diào)整；3利用精細調(diào)整的結(jié)果對粗略估計的結(jié)果進行修正，從而得到更為精確的組織位移估計。


圖1為一般的組織位移估計方法得到的組織應變分布的計算機仿真結(jié)果；圖2為本發(fā)明提出的采用兩種尺度的組織位移估計方法的實施例得到的組織應變分布的計算機仿真結(jié)果。
具體實施例方式
本發(fā)明提出的采用兩種尺度的組織位移估計方法結(jié)合具體實施例及附圖詳細說明如下本發(fā)明的實施例利用計算機程序和一般的超聲散射模型仿真得到一塊模擬的組織在壓縮前和壓縮后的一維射頻信號。組織厚度為60mm，其中0-30mm處的彈性模量是30-60mm處的彈性模量的3倍，采用簡單的一維模型，因此0-30mm處的應變是30-60mm處的應變的1/3。對組織施加的壓縮量為3.33％，即1.98mm；探頭中心頻率為3.5MHz，-3dB帶寬為2.0MHz，因為采用簡單的一維模型，所以只有一條掃描線；射頻信號的采樣頻率為20MHz，假設超聲波在組織內(nèi)的傳播速度為1540m/s，因此1mm的組織對應1mm/(1540×103mm/s/(20×106Hz)/2≈26個數(shù)據(jù)，因為組織深度為60mm，所以，每一條掃描線的數(shù)據(jù)為60×26＝1560個，對組織施加的壓縮量以采樣數(shù)據(jù)的個數(shù)表示為60×3.33％×26≈52個采樣數(shù)據(jù)。
本實施例采用的大尺度Ta為4mm，搜索范圍為0到52；采用的小尺度Tb為1mm(也可采用2mm或在1mm至2mm范圍內(nèi)的任何尺度)，搜索范圍為-10到10(也可為-20到-20或滿足10≤τ4≤26和τ3＝-τ4的從τ3到τ4的任何范圍)。
本實施例的具體步驟如下1.設組織組織壓縮前、后的一維射頻信號(即掃描線數(shù)據(jù))分別為s1(n)和s2(n)；n表示該兩條掃描線上的數(shù)據(jù)序號，1≤n≤1560；2.從該掃描線數(shù)據(jù)s1(n)中取一小段尺度為Ta的數(shù)據(jù)d1，Ta＝4mm，其數(shù)據(jù)個數(shù)為U，U＝104，該數(shù)據(jù)的序號從32到135；在0到52的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R(τ)，計算公式如下R′(τ)=Σi=n1n1+U-1s1(i)s2(i-τ)Σi=n1n1+U-1s12(i)·Σi=n1n1+U-1s22(i-τ)(0≤τ≤52)]]>3.確定該互相關函數(shù)R′(τ)的最大值對應的位置t′1，t′1就是利用尺度Ta進行位移估計時數(shù)據(jù)d1在組織壓縮后的位移(即s1(n)中的序號從32到135的小段數(shù)據(jù)d1的在組織壓縮后移動到s2(n)中的序號從32-t′1到135-t′1的位置)；
4.取s1(n)中的序號從71到96的一小段數(shù)據(jù)e1，該數(shù)據(jù)e1位于數(shù)據(jù)d1的中央位置；該數(shù)據(jù)e1的尺度為Tb，Tb＝1mm，數(shù)據(jù)個數(shù)為Ub，Ub＝26；在-10到10的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R″(τ)，計算公式如下R′′(τ)=Σi=m1m1+Ub-1s1(i)s2(i-t′1-τ)Σi=m1m1+Ub-1s12(i)·Σi=m1m1+Ub-1s22(i-t′1-τ)(-10≤τ≤10)]]>5.確定步驟4得到的互相關函數(shù)R″(τ)的最大值對應的位置t″1，t″1就是利用尺度Tb進行位移估計時數(shù)據(jù)e1的位移相對于t′1的偏移；6.利用尺度Tb估計得到的位移偏移t″1對尺度Ta估計得到的位移t′1進行調(diào)整，得到的小段數(shù)據(jù)d1的精細調(diào)整后的位移t1，即t1＝t′1+t″17.依次從掃描線數(shù)據(jù)s1(n)中取一小段尺度為Ta即數(shù)據(jù)個數(shù)為Ua的數(shù)據(jù)d2、d3、…、dN，每段數(shù)據(jù)的序號依次錯開32個采樣數(shù)據(jù)，直到再錯開32個采樣數(shù)據(jù)將超出s1(n)的范圍，按步驟2-6相同的方法依次得到各段數(shù)據(jù)對應的位移t2、t3、…、tN，其中N為小段數(shù)據(jù)的總數(shù)，N＝45；則位移序列t1、t2…、tN為掃描線數(shù)據(jù)s1(n)對應的組織的位移估計。
本實施例與一般方法的位移估計效果比較如下圖1和圖2為本實施例的結(jié)果與一般方法的結(jié)果的比較。圖1和圖2中，橫坐標表示組織的不同位置，縱坐標表示應變的大??；圖1中的虛線11和圖2中的虛線21表示組織的實際應變分布；圖1中的12為采用一般方法估計得到的組織應變分布，圖2中的22為采用本發(fā)明提出的采用兩種尺度的方法估計得到的組織應變分布。
由圖1和圖2可見，本發(fā)明提出的采用兩種尺度的方法采用兩種不同尺度的跟蹤波段進行位移估計，結(jié)果比采用尺度固定不變的跟蹤段進行位移估計的一般方法的結(jié)果要好，位移估計的精度更高，特別是在組織應變比較大的區(qū)域(30-60mm位置)更為明顯。
權利要求
1.一種采用兩種尺度的組織位移估計方法，包括以下步驟1)從組織壓縮前、后的二維射頻信號中分別取出第一條掃描線的數(shù)據(jù)，設為s1(n)和s2(n)，n表示該兩條掃描線上的數(shù)據(jù)序號，1≤n≤nmax，nmax為n的最大值；2)從該掃描線數(shù)據(jù)s1(n)中取一小段尺度為Ta的數(shù)據(jù)d1，其數(shù)據(jù)個數(shù)為Ua，Ua＝round(Ta×U1)，其中，Ta的單位為mm，U1代表1mm的組織對應的數(shù)據(jù)個數(shù)，round(·)代表四舍五入的取整操作，該數(shù)據(jù)d1的序號從n1到n1+Ua-1，1≤n1≤Ua；在τ1到τ2的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R(τ)，計算公式如下R′(τ)=Σi=n1n1+U-1s1(i)s2(i-τ)Σi=n1n1+U-1s12(i)·Σi=n1n1+U-1s22(i-τ)(τ1≤τ≤τ2)]]>其中i為計算過程中表示數(shù)據(jù)序號的循環(huán)變量，τ1為0，τ2為對組織施加的壓縮量，以采樣數(shù)據(jù)的個數(shù)表示；3)確定步驟2)得到的互相關函數(shù)R′(τ)的最大值對應的位置t′1，則t′1為利用尺度Ta進行位移估計時數(shù)據(jù)d1在組織壓縮后的位移(即s1(n)中的序號從n1到n1+Ua-1的小段數(shù)據(jù)d1的在組織壓縮后移動到s2(n)中的序號從n1-t′1到n1+Ua-1-t′1的位置)；4)再從該掃描線數(shù)據(jù)s1(n)中的數(shù)據(jù)段d1的中央位置取尺度為Tb、序號從n1+round(Ua-Ub2)]]>到n1+round(Ua-Ub2)+Ub-1]]>的一小段數(shù)據(jù)e1，該數(shù)據(jù)個數(shù)為Ub，Ub＝round(Tb×U1)，其中，Tb的單位為mm，并且滿足Tb＜Ta；在τ3到τ4的搜索范圍內(nèi)求該小段數(shù)據(jù)與掃描線數(shù)據(jù)s2(n)的互相關函數(shù)R″(τ)，計算公式如下R′′(τ)=Σi=m1m1+Ub-1s1(i)s2(i-t′1-τ)Σi=m1m1+Ub-1s12(i)·Σi=m1m1+Ub-1s22(i-t′1-τ)(τ3≤τ≤τ4)]]>其中τ3和τ4為組織位移精細調(diào)整的范圍，滿足如下關系τ4＜τ2；5)確定步驟4得到的互相關函數(shù)R″(τ)的最大值對應的位置t″1，則t″1為利用尺度Tb進行位移估計時數(shù)據(jù)e1的位移相對于t′1的偏移；6)利用尺度Tb估計得到的位移偏移t″1對尺度Ta估計得到的位移t′1進行調(diào)整，得到的小段數(shù)據(jù)d1的精細調(diào)整后的位移t1，即t1＝t′1+t″17)依次從掃描線數(shù)據(jù)s1(n)中取一小段尺度為Ta即數(shù)據(jù)個數(shù)為Ua的數(shù)據(jù)d2、d3、…、dN，每段數(shù)據(jù)的序號依次錯開V個采樣數(shù)據(jù)，1≤V≤Ua，直到再錯開V個采樣數(shù)據(jù)將超出s1(n)的范圍，按步驟2-6相同的方法依次得到各段數(shù)據(jù)對應的位移t2、t3、…、tN，其中N為小段數(shù)據(jù)的總數(shù)；則位移序列t1、t2、…、tN為第一條掃描線數(shù)據(jù)s1(n)對應的組織的位移估計；8)利用與步驟1-7相同的方法，依次得到第2、3、…、M條掃描線數(shù)據(jù)對應的組織的位移估計，其中M為表示探頭的掃描線總數(shù)。
2.如權利要求1所述的，其特征在于，還進一步包括在所述第2)、4)步驟中對所述計算得到的互相關函數(shù)進行插值處理。
3.如權利要求1所述的，其特征在于，所述第4)步驟中的尺度Tb的取值范圍為116Ta≤Tb≤12Ta;]]>所述搜索范圍的τ4取值為116τ2≤τ4≤12τ2.]]>
全文摘要
本發(fā)明屬于超聲彈性成像技術領域，涉及采用兩種尺度的生物組織位移估計方法，該方法包括從組織壓縮前二維射頻信號中第一條掃描線的數(shù)據(jù)取出尺度為T
文檔編號A61B8/08GK1586409SQ20041005690
公開日2005年3月2日 申請日期2004年8月20日 優(yōu)先權日2004年8月20日
發(fā)明者白凈, 羅建文 申請人:清華大學