超聲彈性成像位移值的計(jì)算方法
【專利摘要】本發(fā)明提供一種超聲彈性成像位移值的計(jì)算方法,其中,所述方法包括:S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;S3、選取至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值,并將每條數(shù)據(jù)線的先驗(yàn)值作為計(jì)算同一條數(shù)據(jù)線或相鄰數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。本發(fā)明的超聲彈性成像位移值的計(jì)算方法,能夠針對(duì)不同的應(yīng)變大小自適應(yīng)地調(diào)整搜索范圍,大大減少了計(jì)算量,同時(shí)具有很高的魯棒性來(lái)避免錯(cuò)誤擴(kuò)散,保證了彈性成像的實(shí)時(shí)性和準(zhǔn)確性。
【專利說(shuō)明】超聲彈性成像位移值的計(jì)算方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種超聲彈性成像位移值的計(jì)算方法,主要應(yīng)用于醫(yī)療超聲診斷【技術(shù)領(lǐng)域】。
【背景技術(shù)】
[0002]超聲彈性成像是醫(yī)學(xué)超聲成像的新模式,其主要利用超聲對(duì)組織的力學(xué)特性參數(shù)進(jìn)行檢查并成像。
[0003]在人體組織發(fā)生病變時(shí),其力學(xué)特性參數(shù)通常也發(fā)生改變。例如,當(dāng)乳腺癌、前列腺癌等疾病發(fā)生時(shí),其對(duì)應(yīng)的彈性系數(shù)(楊氏模量,表現(xiàn)為硬度)也會(huì)發(fā)生改變。彈性成像為臨床上許多疾病,尤其是腫瘤的檢查與診斷提供了重要幫助。
[0004]彈性成像通常先通過(guò)對(duì)組織施加激勵(lì)(利用探頭壓迫或者利用聲壓進(jìn)行激勵(lì)),然后對(duì)組織發(fā)生形變前后的超聲信號(hào)進(jìn)行分析,計(jì)算出組織的位移分布,然后再利用梯度運(yùn)算計(jì)算出應(yīng)變分布,最后將應(yīng)變分布用偽彩的形式顯示。
[0005]其中,計(jì)算組織的位移分布是彈性成像過(guò)程中的核心環(huán)節(jié),該算法的好壞直接影響了成像的質(zhì)量和成像的速度。一般地,如圖1所示,圖1是現(xiàn)有技術(shù)中超聲彈性成像位移值的計(jì)算方法的流程圖,相應(yīng)的,彈性成像的通常步驟為:
S1、使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào)。
[0006]S2、在兩幀數(shù)字信號(hào)上分段選取位移估計(jì)點(diǎn)。
`[0007]S3、對(duì)其中一幀上的位移估計(jì)點(diǎn)代表的一小段信號(hào)和另一幀上對(duì)應(yīng)位置的信號(hào)做相關(guān)運(yùn)算并搜索相關(guān)系數(shù)的最大值,根據(jù)相關(guān)系數(shù)最大值位置即可獲取該位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0008]S4、對(duì)獲取的絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,獲取應(yīng)變信號(hào)。
[0009]S5、對(duì)應(yīng)變信號(hào)進(jìn)行偽彩處理,作為彈性成像結(jié)果輸出。
[0010]其中,獲取位移估計(jì)點(diǎn)的位移值是彈性成像中最重要的也是最耗時(shí)的步驟,不僅因?yàn)槠溆?jì)算本身耗時(shí)較多,而且由于搜索第一幀某位移估計(jì)點(diǎn)代表的一段信號(hào)在第二幀的位置時(shí),沒(méi)有一個(gè)搜索范圍,需要使用第二幀信號(hào)中一個(gè)很大范圍內(nèi)的所有信號(hào)與該位移估計(jì)點(diǎn)的信號(hào)做匹配計(jì)算,大大增加了計(jì)算時(shí)間,成為實(shí)時(shí)彈性成像的一個(gè)重要阻礙。
【發(fā)明內(nèi)容】
[0011]為了解決上述問(wèn)題,本發(fā)明的目的在于提供一種超聲彈性成像位移值的計(jì)算方法,該方法通過(guò)區(qū)域邊界限制的先驗(yàn)值搜索算法計(jì)算位移估計(jì)點(diǎn)的絕對(duì)位移值,減少了超聲彈性成像的計(jì)算時(shí)間。
[0012]相應(yīng)地,本發(fā)明一實(shí)施方式的超聲彈性成像位移值的計(jì)算方法,所述方法包括以下步驟:
S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;
S3、選取至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值,并將每條數(shù)據(jù)線的先驗(yàn)值作為計(jì)算同一條數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。
[0013]作為本發(fā)明的進(jìn)一步改進(jìn),所述回波信號(hào)劃分的多個(gè)部分中包括第一部分和第二部分,所述S3步驟包括:
選取第一部分,計(jì)算所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第一部分先驗(yàn)值;
根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值作為同一條數(shù)據(jù)線上第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0014]作為本發(fā)明的進(jìn)一步改進(jìn),所述第一部分和所述第二部分之間間隔有至少一個(gè)部分。
[0015]作為本發(fā)明的進(jìn)一步改進(jìn),所述回波信號(hào)劃分的多個(gè)部分還包括第三部分,所述第三部分在軸向上位于所述第一部分和第二部分之間;所述S3步驟包括:
根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第二部分先驗(yàn)值;
根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成同一條數(shù)據(jù)線上第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0016]作為本發(fā)明的進(jìn)一步改進(jìn),“根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值”步驟具體包括:
濾除已計(jì)算出的絕對(duì)位移值中的噪聲;
取濾除噪聲后的絕對(duì)位移值的中值或取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0017]相應(yīng)地,本發(fā)明另一實(shí)施方式的超聲彈性成像位移值的計(jì)算方法,所述方法包括以下步驟:
51、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;
52、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn);
53、通過(guò)搜索算法計(jì)算至少一條數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值;
54、根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出先驗(yàn)值,并根據(jù)所述先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0018]作為本發(fā)明的進(jìn)一步改進(jìn),“根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出先驗(yàn)值”步驟具體包括:
濾除已計(jì)算出的絕對(duì)位移值中的噪聲;
取濾除噪聲后的絕對(duì)位移值的中值或取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0019]作為本發(fā)明的進(jìn)一步改進(jìn),“根據(jù)所述先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值”步驟具體包括:
沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分,根據(jù)相鄰數(shù)據(jù)線所對(duì)應(yīng)的部分在軸向上的位置關(guān)系及其先驗(yàn)值形成搜索區(qū)間,以通過(guò)搜索算法計(jì)算出相鄰數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0020]相應(yīng)地,本發(fā)明再一實(shí)施方式的超聲彈性成像位移值的計(jì)算方法,所述方法包括以下步驟:
所述方法包括以下步驟:
S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;
S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;
S3、通過(guò)搜索算法計(jì)算某一條數(shù)據(jù)線中某一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值;
S4、根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值,并根據(jù)所述部分先驗(yàn)值計(jì)算該條數(shù)據(jù)線的另一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,直至獲得該條數(shù)據(jù)線所有位移估計(jì)點(diǎn)的絕對(duì)位移值;
S5、根據(jù)該條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出整體先驗(yàn)值,并根據(jù)所述整體先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0021]作為本發(fā)明的進(jìn)一步改進(jìn),所述回波信號(hào)劃分的多個(gè)部分中包括第一部分和第二部分,所述S3步驟包括:
選取某一條數(shù)據(jù)線的第一部分,計(jì)算所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出第一部分先驗(yàn)值;
根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將第一部分先驗(yàn)值作為第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0022]作為本發(fā)明的進(jìn)一步改進(jìn),所述第一部分和所述第二部分之間間隔有至少一個(gè)部分。
[0023]作為本發(fā)明的進(jìn)一步改進(jìn),所述回波信號(hào)劃分的多個(gè)部分還包括第三部分,所述第三部分在軸向上位于所述第一部分和第二部分之間,所述S3步驟包括:
根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出第二部分先驗(yàn)值;
根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將第一部分先驗(yàn)值和第二部分先驗(yàn)值形成第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0024]作為本發(fā)明的進(jìn)一步改進(jìn),“根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值”步驟具體包括:
濾除已計(jì)算出的絕對(duì)位移值中的噪聲;
取濾除噪聲后的絕對(duì)位移值的中值或取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0025]與現(xiàn)有技術(shù)相比,本發(fā)明的超聲彈性成像位移值的計(jì)算方法,能夠針對(duì)不同的應(yīng)變大小自適應(yīng)地調(diào)整搜索范圍,大大減少了計(jì)算量,同時(shí)具有很高的魯棒性來(lái)避免錯(cuò)誤擴(kuò)散,保證了彈性成像的實(shí)時(shí)性和準(zhǔn)確性。
【專利附圖】
【附圖說(shuō)明】
[0026]圖1是現(xiàn)有技術(shù)中超聲彈性成像位移值的計(jì)算方法的流程圖;圖2是本發(fā)明第一實(shí)施方式中超聲彈性成像位移值的計(jì)算方法的流程圖;
圖3是本發(fā)明第二實(shí)施方式中超聲彈性成像位移值的計(jì)算方法的流程圖;
圖4是本發(fā)明第三實(shí)施方式中超聲彈性成像位移值的計(jì)算方法的流程圖;
圖5是單幀回波信號(hào)區(qū)域劃分示意圖;
圖6是本發(fā)明超聲彈性成像位移值的計(jì)算裝置的結(jié)構(gòu)示意圖。
【具體實(shí)施方式】
[0027]以下將結(jié)合附圖所示的各實(shí)施方式對(duì)本發(fā)明進(jìn)行詳細(xì)描述。但這些實(shí)施方式并不限制本發(fā)明,本領(lǐng)域的普通技術(shù)人員根據(jù)這些實(shí)施方式所做出的結(jié)構(gòu)、方法、或功能上的變換均包含在本發(fā)明的保護(hù)范圍內(nèi)。
[0028]如圖2所示,在本發(fā)明第一實(shí)施方式中的超聲彈性成像位移值的計(jì)算方法,所述方法包括以下步驟:
S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;
具體的,在準(zhǔn)靜態(tài)壓縮組織的過(guò)程中,使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào),相應(yīng)的,所述數(shù)字信號(hào)為選取壓縮前后的兩幀RF回波信號(hào)或者IQ信號(hào);在以下具體實(shí)施例中,將以RF回波信號(hào)作為具體實(shí)施例來(lái)介紹。
[0029]S2、在每條數(shù)據(jù)線上`選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;
相應(yīng)的,結(jié)合圖5所示,圖5為單幀回波信號(hào)分區(qū)劃分示意圖。
[0030]具體的,沿著單幀回波信號(hào)的軸向方向,即超聲傳播方向,分別在每條數(shù)據(jù)線上按照一定的間隔選取多個(gè)位移估計(jì)點(diǎn)。相應(yīng)的,由于獲取的信號(hào)在軸向上采樣點(diǎn)較密,而位移一般都比較小,故,通常以一定的間隔選取采樣點(diǎn)作為該局部區(qū)域的位移估計(jì)點(diǎn)。
[0031]在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn)后,將每一數(shù)據(jù)線上的位移估計(jì)點(diǎn)沿著軸向方向分為大小相當(dāng)?shù)亩鄠€(gè)部分??梢岳斫獾氖?,在實(shí)際計(jì)算中,可以按照實(shí)際計(jì)算需求,將所述多個(gè)部分的數(shù)值具體設(shè)定。相應(yīng)的,為了方便說(shuō)明,本實(shí)施例中分為5部分為例。
[0032]相應(yīng)的,將單幀回波信號(hào)劃分的5部分,依次為第I部分,第2部分,第3部分,第4部分,第5部分。
[0033]S3、選取至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值,并將每條數(shù)據(jù)線的先驗(yàn)值作為計(jì)算同一條數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。
[0034]相應(yīng)的,若所述回波信號(hào)劃分的多個(gè)部分中僅包括第一部分和第二部分。
[0035]選取第一部分,計(jì)算所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第一部分先驗(yàn)值;根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值作為同一條數(shù)據(jù)線上第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0036]若所述第一部分和所述第二部分之間間隔有至少一個(gè)第三部分。
[0037]則根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第二部分先驗(yàn)值;根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成同一條數(shù)據(jù)線上第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0038]依次類推,計(jì)算出每條數(shù)據(jù)線上各部分的位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0039]具體的,結(jié)合圖5所示,在對(duì)每一條數(shù)據(jù)線計(jì)算其包含的所有絕對(duì)位移值的過(guò)程中,按照先計(jì)算第I部分,然后第5部分,第3部分,第2部分,第4部分的順序進(jìn)行計(jì)算。每當(dāng)一個(gè)部分的絕對(duì)位移值計(jì)算完成后,選取這部分絕對(duì)位移值的中值作為先驗(yàn)值,供相鄰部分使用;也可以選取對(duì)絕對(duì)位移值取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0040]優(yōu)選的,濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲后的絕對(duì)位移值的中值作為先驗(yàn)值。
[0041]或?yàn)V除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲并取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0042]相應(yīng)的,選取第I部分,采用傳統(tǒng)搜索算法計(jì)算所述第I部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第I部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的
第一部分先驗(yàn)值。
[0043]為了盡可能的縮小搜索范圍,在得到第一部分先驗(yàn)值后,以獲取的第一部分先驗(yàn)值為基礎(chǔ),選取當(dāng)前數(shù)據(jù)線上的最后一部分進(jìn)行計(jì)算,在本實(shí)施例中,首先選取第5部分進(jìn)行計(jì)算。當(dāng)然,實(shí)際計(jì)算中,在得到第一部分先驗(yàn)值后,也可以選取當(dāng)前數(shù)據(jù)線上的第2部分,第3部分或第4部分先進(jìn)行計(jì)算,在得到第一部分先驗(yàn)值后,選取其他部分進(jìn)行計(jì)算相較與選取最后一部分進(jìn)行計(jì)算,會(huì)擴(kuò)大其計(jì)算過(guò)程中的搜索范圍,增加工作量。
[0044]相應(yīng)的,以第一部分先驗(yàn)值為基礎(chǔ),計(jì)算出第5部分中包含的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,按照第一部分的計(jì)算方法,計(jì)算出每條數(shù)據(jù)線上的第二部分先驗(yàn)值。
[0045]相應(yīng)的,根據(jù)當(dāng)前數(shù)據(jù)線各部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成當(dāng)前數(shù)據(jù)線上第3部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出每條數(shù)據(jù)線上第3部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0046]在本發(fā)明的【具體實(shí)施方式】中,為了方便描述,以一條數(shù)據(jù)線為例做具體介紹:
假設(shè)以第一個(gè)位移估計(jì)點(diǎn)為參考零點(diǎn),選取第I部分,獲得了各個(gè)位移估計(jì)點(diǎn)的絕對(duì)
位移值為:1.4,13,1.6,2,31,…5.3,9,5.6 ;本領(lǐng)域技術(shù)人員公知:所述絕對(duì)位移值應(yīng)該是單調(diào)的,故判斷出13,31,9為噪聲等原因引起的錯(cuò)誤,濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算出的中值為4.2,則以第I部分計(jì)算出的絕對(duì)位移值4.2作為當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值。
[0047]根據(jù)第I部分、第2部分、第3部分、第4部分,以及第5部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值作為第5部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第5部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0048]計(jì)算第5部 分各位移估計(jì)點(diǎn)的絕對(duì)位移值過(guò)程中,將從第一部分先驗(yàn)值4.2以后開(kāi)始進(jìn)行相關(guān)運(yùn)算;相應(yīng)的,第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值分別為:10.2,10.5….13.5,在濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算得出的中值為12.2;則以第5部分計(jì)算出的絕對(duì)位移值12.2作為當(dāng)前數(shù)據(jù)線的第二部分先驗(yàn)值。
[0049]接著計(jì)算第3部分,在區(qū)間[4.2,12.2]上面計(jì)算各位移估計(jì)點(diǎn)的絕對(duì)位移值,獲得一系列值,假設(shè)其計(jì)算出的中值為8.5。則以第3部分計(jì)算出的絕對(duì)位移值8.5作為當(dāng)前數(shù)據(jù)線的第三部分先驗(yàn)值。
[0050]最后計(jì)算第2部分和第4部分,分別在[4.2,8.5],[8.5,12.2]的區(qū)間上進(jìn)行運(yùn)算,直到獲取全部數(shù)據(jù)線上的所有絕對(duì)位移值。采用該種算法獲取全部絕對(duì)位移值,大大減少了搜索范圍。而且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性。
[0051]S4、依據(jù)獲取的所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0052]在本發(fā)明的第一實(shí)施方式中,僅是一幀回波信號(hào)中的一條數(shù)據(jù)線為例做具體介紹,實(shí)際運(yùn)算過(guò)程中,一幀回波信號(hào)中所包含的每條數(shù)據(jù)線同時(shí)進(jìn)行計(jì)算,最后將每條數(shù)據(jù)線中的全部所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0053]本發(fā)明的超聲彈性成像位移值的計(jì)算方法,能夠針對(duì)不同的應(yīng)變大小自適應(yīng)地調(diào)整搜索范圍,大大減少了計(jì)算量,且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性,具有很高的魯棒性來(lái)避免錯(cuò)誤擴(kuò)散,保證了彈性成像的實(shí)時(shí)性和準(zhǔn)確性。
[0054]如圖3所示,在本發(fā)明二實(shí)施方式中的超聲彈性成像位移值的計(jì)算方法。
[0055]所述第二實(shí)施方式與第一實(shí)施方式的區(qū)別在于,首先采用傳統(tǒng)搜索算法計(jì)算至少一條數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值;之后根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出先驗(yàn)值,并根據(jù)所述 先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0056]具體的,本發(fā)明第二實(shí)施方式中的超聲彈性成像位移值的計(jì)算方法包括以下步驟:
S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;
具體的,在準(zhǔn)靜態(tài)壓縮組織的過(guò)程中,使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào),相應(yīng)的,所述數(shù)字信號(hào)為選取壓縮前后的兩幀RF回波信號(hào)或者IQ信號(hào);在以下具體實(shí)施例中,將以RF回波信號(hào)作為具體實(shí)施例來(lái)介紹。
[0057]S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn);
S3、通過(guò)搜索算法計(jì)算至少一條數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值;
相應(yīng)的,步驟S2、S3中所述求取一條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值的算法,采用傳統(tǒng)的算法求取,在此不做詳細(xì)贅述。
[0058]S4、根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出先驗(yàn)值,并根據(jù)所述先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0059]相應(yīng)的,在傳統(tǒng)搜索算法中,濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲后的絕對(duì)位移值的中值作為先驗(yàn)值。
[0060]或?yàn)V除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲并取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0061]相應(yīng)的,沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分,根據(jù)相鄰數(shù)據(jù)線所對(duì)應(yīng)的部分在軸向上的位置關(guān)系及其先驗(yàn)值形成搜索區(qū)間,以通過(guò)搜索算法計(jì)算出相鄰數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0062]具體的,結(jié)合圖5所示,例如:通過(guò)傳統(tǒng)搜索算法計(jì)算得出,第一數(shù)據(jù)線上的位移估計(jì)點(diǎn)的所有絕對(duì)位移值依次為:1.4,13,1.6,2,31,…5.3,9,5.6 ;之后沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為5個(gè)部分,劃分后,本領(lǐng)域技術(shù)人員公知:所述絕對(duì)位移值應(yīng)該是單調(diào)的,故判斷出13,31,9為噪聲等原因引起的錯(cuò)誤,濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,求取第一數(shù)據(jù)線的各部分中值依次為:第I部分的中值為4.2,第2部分的中值為6.3,第3部分的中值為8.5,第4部分的中值為10.3,第5部分的中值為12.2。
[0063]通過(guò)第一數(shù)據(jù)線的各部分中值,即所述第一數(shù)據(jù)線的各部分的先驗(yàn)值,可以確定相鄰數(shù)據(jù)線各部分位移估計(jì)點(diǎn)所在區(qū)域的搜索范圍。相應(yīng)的,計(jì)算與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線中所有位移估計(jì)點(diǎn)的絕對(duì)位移值時(shí),其使用所述第一數(shù)據(jù)線的先驗(yàn)值的算法有很多種,在本實(shí)施方式中,采用所述第一數(shù)據(jù)線的N+1部分的先驗(yàn)值作為與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線的位移最大值這一算法,以求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0064]具體的,采用此種算法,求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上各部分位移估計(jì)點(diǎn)的絕對(duì)位移值所分布的區(qū)間為:第I部分的區(qū)間為:[0,6.3],第2部分[4.2,8.5],第3部分[6.3,10.2],第4部分的區(qū)間為:[8.5,10.3],第5部分大于10.3。
[0065]參照各部分所處于的相關(guān)區(qū)域,直接求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的各位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0066]S5、依據(jù)獲取的所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0067]步驟S5為本領(lǐng)域技術(shù)人員所公知的技術(shù),在此不做詳細(xì)贅述。
[0068]本發(fā)明第二實(shí)施方式中所描述的超聲彈性成像位移值的計(jì)算方法,相較于現(xiàn)有技術(shù),同樣能夠針對(duì)不同的應(yīng)變大小自適應(yīng)地調(diào)整搜索范圍,大大減少了計(jì)算量,且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性,具有很高的魯棒性來(lái)避免錯(cuò)誤擴(kuò)散,保證了彈性成像的實(shí)時(shí)性和準(zhǔn)確性。
[0069]相應(yīng)的,如圖4所示,圖4為本發(fā)明第三實(shí)施方式的超聲彈性成像位移值的計(jì)算方法。所述第三實(shí)施方式是結(jié)合了第一實(shí)施方式和第二實(shí)施方式。其原理為:通過(guò)第一實(shí)施方式求取其中I條數(shù)據(jù)線上各部分的先驗(yàn)值,之后根據(jù)所述先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0070]相應(yīng)的,所述第三方式的超聲彈性成像位移值的計(jì)算方法包括以下步驟:
S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;
具體的,在準(zhǔn)靜態(tài)壓縮組織的過(guò)程中,使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào),相應(yīng)的,所述數(shù)字信號(hào)為選取壓縮前后的兩幀RF回波信號(hào)或者IQ信號(hào);在以下具體實(shí)施例中,將以RF回波信號(hào)作為具體實(shí)施例來(lái)介紹。
[0071]S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分。
[0072]相應(yīng)的,結(jié)合圖5所示,圖5是單幀回波信號(hào)區(qū)域劃分示意圖。
[0073]具體的,沿著單幀回波信號(hào)的軸向方向,即超聲傳播方向,分別在每條數(shù)據(jù)線上按照一定的間隔選取多個(gè)位移估計(jì)點(diǎn)。相應(yīng)的,由于獲取的信號(hào)在軸向上采樣點(diǎn)較密,而位移一般都比較小,故,通常以一定的間隔選取采樣點(diǎn)作為該局部區(qū)域的位移估計(jì)點(diǎn)。
[0074]在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn)后,將每一數(shù)據(jù)線上的位移估計(jì)點(diǎn)沿著軸向方向分為大小相當(dāng)?shù)亩鄠€(gè)部分??梢岳斫獾氖?,在實(shí)際計(jì)算中,可以按照實(shí)際計(jì)算需求,將所述多個(gè)部分的數(shù)值具體設(shè)定。相應(yīng)的,為了方便說(shuō)明,本實(shí)施例中分為5部分為例。
[0075]相應(yīng)的,將單幀回波信號(hào)劃分的5部分,依次為第I部分,第2部分,第3部分,第4部分,第5部分。
[0076]S3、通過(guò)搜索算法計(jì)算某一條數(shù)據(jù)線中某一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0077]具體的,通過(guò)傳統(tǒng)的搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線的第I部分的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0078]S4、根據(jù)每個(gè)位 移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值,并根據(jù)所述部分先驗(yàn)值計(jì)算該條數(shù)據(jù)線上的另一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,直至獲得該條數(shù)據(jù)線上所有位移估計(jì)點(diǎn)的絕對(duì)位移值;
需要說(shuō)明的是,此時(shí)獲取的部分先驗(yàn)值,只是根據(jù)已計(jì)算出的其中一條數(shù)據(jù)線所求取的。
[0079]相應(yīng)的,選取其中一條數(shù)據(jù)線的至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)當(dāng)前數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的先驗(yàn)值,并將所述先驗(yàn)值作為計(jì)算當(dāng)前數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。
[0080]相應(yīng)的,若所述回波信號(hào)劃分的多個(gè)部分中僅包括第一部分和第二部分。
[0081]選取第一部分,計(jì)算當(dāng)前數(shù)據(jù)線中的所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值;
根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值作為當(dāng)前數(shù)據(jù)線上第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0082]若所述第一部分和所述第二部分之間間隔有至少一個(gè)第三部分。
[0083]則根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的第二部分先驗(yàn)值;之后,根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成當(dāng)前數(shù)據(jù)線上第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0084]具體的,結(jié)合圖5所示,在對(duì)當(dāng)前數(shù)據(jù)線計(jì)算其包含的所有絕對(duì)位移值的過(guò)程中,按照先計(jì)算第I部分,然后第5部分,第3部分,第2部分,第4部分的順序進(jìn)行計(jì)算。每當(dāng)一個(gè)部分的絕對(duì)位移值計(jì)算完成后,選取這部分絕對(duì)位移值的中值作為先驗(yàn)值,供相鄰部分使用;也可以選取對(duì)絕對(duì)位移值取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0085]優(yōu)選的,濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲后的絕對(duì)位移值的中值作為先驗(yàn)值。
[0086]或?yàn)V除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲并取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0087]相應(yīng)的,選取第I部分,采用傳統(tǒng)搜索算法計(jì)算所述第I部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第I部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的
第一部分先驗(yàn)值。
[0088]為了盡可能的縮小搜索范圍,在得到第一部分先驗(yàn)值后,以獲取的第一部分先驗(yàn)值為基礎(chǔ),選取當(dāng)前數(shù)據(jù)線上的最后一部分進(jìn)行計(jì)算,在本實(shí)施例中,首先選取第5部分進(jìn)行計(jì)算。當(dāng)然,實(shí)際計(jì)算中,在得到第一部分先驗(yàn)值后,也可以選取當(dāng)前數(shù)據(jù)線上的第2部分,第3部分或第4部分先進(jìn)行計(jì)算,在得到第一部分先驗(yàn)值后,選取其他部分進(jìn)行計(jì)算相較與選取最后一部分進(jìn)行計(jì)算,會(huì)擴(kuò)大其計(jì)算過(guò)程中的搜索范圍,增加工作量。
[0089]相應(yīng)的,以第一部分先驗(yàn)值為基礎(chǔ),計(jì)算出第5部分中包含的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,按照第一部分的計(jì)算方法,計(jì)算出當(dāng)前數(shù)據(jù)線上的第二部分先驗(yàn)值。
[0090]相應(yīng)的,根據(jù)當(dāng)前數(shù)據(jù)線各部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成當(dāng)前數(shù)據(jù)線上第3部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第3部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0091]具體的,在當(dāng)前數(shù)據(jù)線中,假設(shè)以第一個(gè)位移估計(jì)點(diǎn)為參考零點(diǎn),選取第I部分,獲得了各個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值為:1.4,13,1.6,2,31,…5.3,9,5.6 ;本領(lǐng)域技術(shù)人員公知:所述絕對(duì)位移值應(yīng)該`是單調(diào)的,故判斷出13,31,9為噪聲等原因引起的錯(cuò)誤,濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算出的中值為4.2,則以第I部分計(jì)算出的絕對(duì)位移值4.2作為當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值。
[0092]根據(jù)第I部分、第2部分、第3部分、第4部分,以及第5部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值作為第5部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第5部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0093]計(jì)算第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值過(guò)程中,將從第一部分先驗(yàn)值4.2以后開(kāi)始進(jìn)行相關(guān)運(yùn)算;相應(yīng)的,第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值分別為:10.2,10.5….13.5,在濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算得出的中值為12.2;則以第5部分計(jì)算出的絕對(duì)位移值12.2作為當(dāng)前數(shù)據(jù)線的第二部分先驗(yàn)值。
[0094]接著計(jì)算第3部分,在區(qū)間[4.2,12.2]上面計(jì)算各位移估計(jì)點(diǎn)的絕對(duì)位移值,獲得一系列值,假設(shè)其計(jì)算出的中值為8.5。則以第3部分計(jì)算出的絕對(duì)位移值8.5作為當(dāng)前數(shù)據(jù)線的第三部分先驗(yàn)值。
[0095]最后計(jì)算第2部分和第4部分,分別在[4.2,8.5],[8.5,12.2]的區(qū)間上進(jìn)行運(yùn)算,直到獲取當(dāng)前數(shù)據(jù)線上的所有絕對(duì)位移值。采用該種算法獲取全部絕對(duì)位移值,大大減少了搜索范圍。而且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性。
[0096]S5、根據(jù)該條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出整體先驗(yàn)值,并根據(jù)所述整體先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。[0097]具體的,結(jié)合圖5所示,通過(guò)第一數(shù)據(jù)線的各部分中值,即所述第一數(shù)據(jù)線的各部分的先驗(yàn)值,可以確定與所述第一數(shù)據(jù)線相鄰數(shù)據(jù)線各部分位移估計(jì)點(diǎn)所在區(qū)域的搜索范圍。相應(yīng)的,計(jì)算與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線中所有位移估計(jì)點(diǎn)的絕對(duì)位移值時(shí),其使用所述第一數(shù)據(jù)線的先驗(yàn)值的算法有很多種,在本實(shí)施方式中,采用所述第一數(shù)據(jù)線的N+1部分的先驗(yàn)值作為與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線的位移最大值這一算法,以求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0098]具體的,采用此種算法,求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上各部分位移估計(jì)點(diǎn)的絕對(duì)位移值所分布的區(qū)間為:第I部分的區(qū)間為:[0,6.3],第2部分[4.2,8.5],第3部分[6.3,10.2],第4部分的區(qū)間為:[8.5,10.3],第5部分大于10.3。
[0099]參照各部分所處于的相關(guān)區(qū)域,直接求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的各位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0100]S6、依據(jù)獲取的所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0101]步驟S6為本領(lǐng)域技術(shù)人員所公知的技術(shù),在此不做詳細(xì)贅述。
[0102]本發(fā)明的超聲彈性成像位移值的計(jì)算方法,相較于現(xiàn)有技術(shù),能夠針對(duì)不同的應(yīng)變大小自適應(yīng)地調(diào)整搜索范圍,大大減少了計(jì)算量,且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性。
[0103]相應(yīng)的,如圖6所示,圖6是本發(fā)明實(shí)施方式中超聲彈性成像位移值的計(jì)算裝置的結(jié)構(gòu)示意圖。
[0104]在本發(fā)明的第一實(shí)施方式中,所述裝置包括:傳感模塊100,用于獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;處理模塊200,用于在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;選取至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值,并將每條數(shù)據(jù)線的先驗(yàn)值作為計(jì)算同一條數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。
[0105]具體的,所述傳感模塊100在準(zhǔn)靜態(tài)壓縮組織的過(guò)程中,使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào),相應(yīng)的,所述數(shù)字信號(hào)為選取壓縮前后的兩幀RF回波信號(hào)或者IQ信號(hào);在以下具體實(shí)施例中,將以RF回波信號(hào)作為具體實(shí)施例來(lái)介紹。
[0106]相應(yīng)的,結(jié)合圖5所不,圖5為單幀回波信號(hào)分區(qū)劃分不意圖。
[0107]具體的,所述處理模塊200沿著單幀回波信號(hào)的軸向方向,即超聲傳播方向,分別在每條數(shù)據(jù)線上按照一定的間隔選取多個(gè)位移估計(jì)點(diǎn)。相應(yīng)的,由于獲取的信號(hào)在軸向上采樣點(diǎn)較密,而位移一般都比較小,故,通常以一定的間隔選取采樣點(diǎn)作為該局部區(qū)域的位移估計(jì)點(diǎn)。
[0108]所述處理模塊200在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn)后,將每一數(shù)據(jù)線上的位移估計(jì)點(diǎn)沿著軸向方向分為大小相當(dāng)?shù)闹辽?個(gè)部分。可以理解的是,在實(shí)際計(jì)算中,可以按照實(shí)際計(jì)算需求,將所述多個(gè)部分的數(shù)值具體設(shè)定。相應(yīng)的,為了方便說(shuō)明,本實(shí)施例中分為5部分為例。
[0109]相應(yīng)的,將單幀回波信號(hào)劃分的5部分,依次為第I部分,第2部分,第3部分,第4部分,第5部分。
[0110]相應(yīng)的,若所述回波信號(hào)劃分的多個(gè)部分中僅包括第一部分和第二部分。
[0111]所述處理模塊200選取第一部分,計(jì)算所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第一部分先驗(yàn)值;根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值作為同一條數(shù)據(jù)線上第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0112]若所述第一部分和所述第二部分之間間隔有至少一個(gè)第三部分。
[0113]所述處理模塊200則根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第二部分先驗(yàn)值;根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成同一條數(shù)據(jù)線上第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0114]依次類推,計(jì)算出每條數(shù)據(jù)線上各部分的位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0115]具體的,結(jié)合圖5所示,所述處理模塊200在對(duì)每一條數(shù)據(jù)線計(jì)算其包含的所有絕對(duì)位移值的過(guò)程中,按照先計(jì)算第I部分,然后第5部分,第3部分,第2部分,第4部分的順序進(jìn)行計(jì)算。每當(dāng)一個(gè)部分的絕對(duì)位移值計(jì)算完成后,選取這部分絕對(duì)位移值的中值作為先驗(yàn)值,供相鄰部分使用;也可以選取對(duì)絕對(duì)位移值取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0116]優(yōu)選的,濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲后的絕對(duì)位移值的中值作為先驗(yàn)值。
[0117]或?yàn)V除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲并取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0118]相應(yīng)的,所述處理模塊200選取第I部分,采用傳統(tǒng)搜索算法計(jì)算所述第I部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第I部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第一部分先驗(yàn)值。
[0119]為了盡可能的縮小搜索范圍,所述處理模塊200在得到第一部分先驗(yàn)值后,以獲取的第一部分先驗(yàn)值為基礎(chǔ),選取當(dāng)前數(shù)據(jù)線上的最后一部分進(jìn)行計(jì)算,在本實(shí)施例中,首先選取第5部分進(jìn)行計(jì)算。當(dāng)然,實(shí)際計(jì)算中,在得到第一部分先驗(yàn)值后,也可以選取當(dāng)前數(shù)據(jù)線上的第2部分,第3部分或第4部分先進(jìn)行計(jì)算,在得到第一部分先驗(yàn)值后,選取其他部分進(jìn)行計(jì)算相較與選取最后一部分進(jìn)行計(jì)算,會(huì)擴(kuò)大其計(jì)算過(guò)程中的搜索范圍,增加工作量。
[0120]相應(yīng)的,所述處理模塊200以第一部分先驗(yàn)值為基礎(chǔ),計(jì)算出第5部分中包含的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,按照第一部分的計(jì)算方法,計(jì)算出每條數(shù)據(jù)線上的第二部分先驗(yàn)值。[0121]相應(yīng)的,所述處理模塊200根據(jù)當(dāng)前數(shù)據(jù)線各部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成當(dāng)前數(shù)據(jù)線上第3部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出每條數(shù)據(jù)線上第3部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0122]在本發(fā)明的【具體實(shí)施方式】中,為了方便描述,以一條數(shù)據(jù)線為例做具體介紹:
假設(shè)以第一個(gè)位移估計(jì)點(diǎn)為參考零點(diǎn),所述處理模塊200選取第I部分,獲得了各個(gè)位
移估計(jì)點(diǎn)的絕對(duì)位移值為:1.4,13,1.6,2,31,…5.3,9,5.6 ;本領(lǐng)域技術(shù)人員公知:絕對(duì)位移值應(yīng)該是單調(diào)的,故判斷出13,31,9為噪聲等原因引起的錯(cuò)誤,濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算出的中值為4.2,則以第I部分計(jì)算出的絕對(duì)位移值4.2作為當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值。
[0123]所述處理模塊200根據(jù)第I部分、第2部分、第3部分、第4部分,以及第5部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值作為第5部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第5部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0124]所述處理模塊200計(jì)算第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值過(guò)程中,將從第一部分先驗(yàn)值4.2以后開(kāi)始進(jìn)行相關(guān)運(yùn)算;相應(yīng)的,第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值分別為:10.2,10.5....13.5,在濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算得出的中值為12.2 ;則以第5部分計(jì)算出的絕對(duì)位移值12.2作為當(dāng)前數(shù)據(jù)線的第二部分先驗(yàn)值。 [0125]接著計(jì)算第3部分,所述處理模塊200在區(qū)間[4.2,12.2]上面計(jì)算各位移估計(jì)點(diǎn)的絕對(duì)位移值,獲得一系列值,假設(shè)其計(jì)算出的中值為8.5。則以第3部分計(jì)算出的絕對(duì)位移值8.5作為當(dāng)前數(shù)據(jù)線的第三部分先驗(yàn)值。
[0126]所述處理模塊200最后計(jì)算第2部分和第4部分,所述處理模塊200分別在[4.2,8.5],[8.5,12.2]的區(qū)間上進(jìn)行運(yùn)算,直到獲取全部數(shù)據(jù)線上的所有絕對(duì)位移值。采用該種算法獲取全部絕對(duì)位移值,大大減少了搜索范圍。而且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性。
[0127]相應(yīng)的,所述處理模塊200還用于依據(jù)獲取的所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0128]在本發(fā)明的第一實(shí)施方式中,僅是一幀回波信號(hào)中的一條數(shù)據(jù)線為例做具體介紹,實(shí)際運(yùn)算過(guò)程中,一幀回波信號(hào)中所包含的每條數(shù)據(jù)線同時(shí)進(jìn)行計(jì)算,最后所述處理模塊200將每條數(shù)據(jù)線中的全部所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0129]相應(yīng)的本發(fā)明二實(shí)施方式中超聲彈性成像位移值的計(jì)算裝置包括:傳感模塊100,用于獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;處理模塊200,用于在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;通過(guò)搜索算法計(jì)算某一條數(shù)據(jù)線中某一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值;根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值,計(jì)算該條數(shù)據(jù)線上的另一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,直至獲得該條數(shù)據(jù)線上所有位移估計(jì)點(diǎn)的位移值;根據(jù)該條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出整體先驗(yàn)值,并根據(jù)所述整體先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。[0130]具體的,所述傳感模塊100,在準(zhǔn)靜態(tài)壓縮組織的過(guò)程中,使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào),相應(yīng)的,所述數(shù)字信號(hào)為選取壓縮前后的兩幀RF回波信號(hào)或者IQ信號(hào);在以下具體實(shí)施例中,將以RF回波信號(hào)作為具體實(shí)施例來(lái)介紹。
[0131]相應(yīng)的,所述處理模塊200求取一條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值的算法,采用傳統(tǒng)的算法求取,在此不做詳細(xì)贅述。
[0132]所述處理模塊200還用于在傳統(tǒng)搜索算法中,濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲后的絕對(duì)位移值的中值作為先驗(yàn)值。
[0133]或所述處理模塊200用于濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲并取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0134]相應(yīng)的,所述處理模塊200沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分,根據(jù)相鄰數(shù)據(jù)線所對(duì)應(yīng)的部分在軸向上的位置關(guān)系及其先驗(yàn)值形成搜索區(qū)間,以通過(guò)搜索算法計(jì)算出相鄰數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0135]具體的,結(jié)合圖5所示,例如:所述處理模塊200通過(guò)傳統(tǒng)搜索算法計(jì)算得出,第一數(shù)據(jù)線上的位移估計(jì)點(diǎn)的所有絕對(duì)位移值依次為:1.4,13,1.6,2,31,…5.3,9,5.6 ;之后沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為5個(gè)部分,劃分后,本領(lǐng)域技術(shù)人員公知:所述絕對(duì)位移值應(yīng)該是單調(diào)的,故判斷出13,31,9為噪聲等原因引起的錯(cuò)誤,濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,求取第一數(shù)據(jù)線的各部分中值依次為--第I部分的中值為4.2,第2部分的中值為6.3,第3部分的中值為8.5,第4部分的中值為10.3,第5部分的中值為12.2。
[0136]所述處理模塊200通過(guò)第一數(shù)據(jù)線的各部分中值,即所述第一數(shù)據(jù)線的各部分的先驗(yàn)值,可以確定相鄰數(shù)據(jù)線`各部分位移估計(jì)點(diǎn)所在區(qū)域的搜索范圍。相應(yīng)的,計(jì)算與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線中所有位移估計(jì)點(diǎn)的絕對(duì)位移值時(shí),其使用所述第一數(shù)據(jù)線的先驗(yàn)值的算法有很多種,在本實(shí)施方式中,采用所述第一數(shù)據(jù)線的N+1部分的先驗(yàn)值作為與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線的位移最大值這一算法,以求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,所述處理模塊200繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0137]具體的,所述處理模塊200采用此種算法,求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上各部分位移估計(jì)點(diǎn)的絕對(duì)位移值所分布的區(qū)間為:第I部分的區(qū)間為:[0,6.3],第2部分[4.2,8.5],第3部分[6.3,10.2],第4部分的區(qū)間為:[8.5,10.3],第5部分大于10.3。
[0138]參照各部分所處于的相關(guān)區(qū)域,所述處理模塊200直接求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的各位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,所述處理模塊200繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0139]所述處理模塊200還用于:依據(jù)獲取的所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0140]相應(yīng)的本發(fā)明三實(shí)施方式中超聲彈性成像位移值的計(jì)算裝置包括:傳感模塊100,用于獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線;處理模塊200用于:在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分;通過(guò)搜索算法計(jì)算某一條數(shù)據(jù)線中某一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值;根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值,計(jì)算該條數(shù)據(jù)線上的另一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,直至獲得該條數(shù)據(jù)線上所有位移估計(jì)點(diǎn)的位移值;根據(jù)該條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出整體先驗(yàn)值,并根據(jù)所述整體先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0141]具體的,所述傳感模塊100用于在準(zhǔn)靜態(tài)壓縮組織的過(guò)程中,使用探頭壓迫或者超聲產(chǎn)生的聲壓力等沿著一個(gè)方向壓迫組織,以此獲得組織受壓過(guò)程中的兩幀數(shù)字信號(hào),相應(yīng)的,所述數(shù)字信號(hào)為選取壓縮前后的兩幀RF回波信號(hào)或者IQ信號(hào);在以下具體實(shí)施例中,將以RF回波信號(hào)作為具體實(shí)施例來(lái)介紹。
[0142]具體的,結(jié)合圖5所示,圖5是單幀回波信號(hào)區(qū)域劃分示意圖。
[0143]具體的,所述處理模塊200沿著單幀回波信號(hào)的軸向方向,即超聲傳播方向,分別在每條數(shù)據(jù)線上按照一定的間隔選取多個(gè)位移估計(jì)點(diǎn)。相應(yīng)的,由于獲取的信號(hào)在軸向上采樣點(diǎn)較密,而位移一般都比較小,故,通常以一定的間隔選取采樣點(diǎn)作為該局部區(qū)域的位移估計(jì)點(diǎn)。
[0144]所述處理模塊200在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn)后,將每一數(shù)據(jù)線上的位移估計(jì)點(diǎn)沿著軸向方向分為大小相當(dāng)?shù)亩鄠€(gè)部分??梢岳斫獾氖?,在實(shí)際計(jì)算中,可以按照實(shí)際計(jì)算需求,所述處理模塊200將所述多個(gè)部分的數(shù)值具體設(shè)定。相應(yīng)的,為了方便說(shuō)明,本實(shí)施例中分為5部分為例。
[0145]相應(yīng)的,所述處理模塊200將單幀回波信號(hào)劃分的5部分,依次為第I部分,第2部分,第3部分,第4部分,第5部分。
[0146]具體的,所述處`理模塊200通過(guò)傳統(tǒng)的搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線的第I部分的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0147]所述處理模塊200根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值,并根據(jù)所述部分先驗(yàn)值計(jì)算該條數(shù)據(jù)線上的另一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,直至獲得該條數(shù)據(jù)線上所有位移估計(jì)點(diǎn)的絕對(duì)位移值;需要說(shuō)明的是,此時(shí)獲取的部分先驗(yàn)值,只是根據(jù)已計(jì)算出的其中一條數(shù)據(jù)線所求取的。
[0148]相應(yīng)的,所述處理模塊200選取其中一條數(shù)據(jù)線的至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)當(dāng)前數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的先驗(yàn)值,并將所述先驗(yàn)值作為計(jì)算當(dāng)前數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。
[0149]相應(yīng)的,若所述回波信號(hào)劃分的多個(gè)部分中僅包括第一部分和第二部分。
[0150]所述處理模塊200選取第一部分,計(jì)算當(dāng)前數(shù)據(jù)線中的所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值;
所述處理模塊200根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值作為當(dāng)前數(shù)據(jù)線上第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0151 ] 若所述第一部分和所述第二部分之間間隔有至少一個(gè)第三部分。[0152]所述處理模塊200則根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的第二部分先驗(yàn)值;
之后,所述處理模塊200根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成當(dāng)前數(shù)據(jù)線上第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0153]具體的,結(jié)合圖5所示,所述處理模塊200在對(duì)當(dāng)前數(shù)據(jù)線計(jì)算其包含的所有絕對(duì)位移值的過(guò)程中,按照先計(jì)算第I部分,然后第5部分,第3部分,第2部分,第4部分的順序進(jìn)行計(jì)算。每當(dāng)一個(gè)部分的絕對(duì)位移值計(jì)算完成后,選取這部分絕對(duì)位移值的中值作為先驗(yàn)值,供相鄰部分使用;中值可以采用中值作為先驗(yàn)值的算法,也可以選取對(duì)絕對(duì)位移值取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0154]優(yōu)選的,所述處理模塊200濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲后的絕對(duì)位移值的中值作為先驗(yàn)值。
[0155]或所述處理模塊200濾除已計(jì)算出的絕對(duì)位移值中的噪聲;取濾除噪聲并取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
[0156]相應(yīng)的,所述處理模塊200選取第I部分,采用傳統(tǒng)搜索算法計(jì)算所述第I部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第I部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值。
[0157]為了盡可能的縮小搜索范圍,所述處理模塊200在得到第一部分先驗(yàn)值后,以獲取的第一部分先驗(yàn)值為基礎(chǔ),選取當(dāng)前數(shù)據(jù)線上的最后一部分進(jìn)行計(jì)算,在本實(shí)施例中,首先選取第5部分進(jìn)行計(jì)算。當(dāng)然,實(shí)際計(jì)算中,在得到第一部分先驗(yàn)值后,也可以選取當(dāng)前數(shù)據(jù)線上的第2部分,第3部分或第4部分先進(jìn)行計(jì)算,在得到第一部分先驗(yàn)值后,選取其他部分進(jìn)行計(jì)算相較與選取最后一部分進(jìn)行計(jì)算,會(huì)擴(kuò)大其計(jì)算過(guò)程中的搜索范圍,增加工作量。
[0158]相應(yīng)的,所述處理模塊200以第一部分先驗(yàn)值為基礎(chǔ),計(jì)算出第5部分中包含的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,按照第一部分的計(jì)算方法,計(jì)算出當(dāng)前數(shù)據(jù)線上的第二部分先驗(yàn)值。
[0159]相應(yīng)的,所述處理模塊200根據(jù)當(dāng)前數(shù)據(jù)線各部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成當(dāng)前數(shù)據(jù)線上第3部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第3部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0160]具體的,在當(dāng)前數(shù)據(jù)線中,假設(shè)以第一個(gè)位移估計(jì)點(diǎn)為參考零點(diǎn),所述處理模塊200選取第I部分,獲得了各個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值為:1.4,13,1.6,2,31,…5.3,9,5.6 ;本領(lǐng)域技術(shù)人員公知:絕對(duì)位移值應(yīng)該是單調(diào)的,故判斷出13,31,9為噪聲等原因引起的錯(cuò)誤,濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算出的中值為4.2,則以第I部分計(jì)算出的絕對(duì)位移值4.2作為當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值。
[0161]所述處理模塊200根據(jù)第I部分、第2部分、第3部分、第4部分,以及第5部分在軸向上的位置關(guān)系,將當(dāng)前數(shù)據(jù)線的第一部分先驗(yàn)值作為第5部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出當(dāng)前數(shù)據(jù)線上第5部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0162]所述處理模塊200計(jì)算第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值過(guò)程中,將從第一部分先驗(yàn)值4.2以后開(kāi)始進(jìn)行相關(guān)運(yùn)算;相應(yīng)的,第5部分各位移估計(jì)點(diǎn)的絕對(duì)位移值分別為:10.2,10.5....13.5,在濾除已計(jì)算出的絕對(duì)位移值中的噪聲后,假設(shè)計(jì)算得出的中值為
12.2 ;則以第5部分計(jì)算出的絕對(duì)位移值12.2作為當(dāng)前數(shù)據(jù)線的第二部分先驗(yàn)值。
[0163]接著計(jì)算第3部分,所述處理模塊200在區(qū)間[4.2,12.2]上面計(jì)算各位移估計(jì)點(diǎn)的絕對(duì)位移值,獲得一系列值,假設(shè)其計(jì)算出的中值為8.5。則以第3部分計(jì)算出的絕對(duì)位移值8.5作為當(dāng)前數(shù)據(jù)線的第三部分先驗(yàn)值。
[0164]最后計(jì)算第2部分和第4部分,所述處理模塊200分別在[4.2,8.5], [8.5,12.2]的區(qū)間上進(jìn)行運(yùn)算,直到獲取當(dāng)前數(shù)據(jù)線上的所有絕對(duì)位移值。所述處理模塊200采用該種算法獲取全部絕對(duì)位移值,大大減少了搜索范圍。而且由于采取的是使用相鄰部分的中值作為先驗(yàn)值,只要這部分里面的多數(shù)位移估計(jì)點(diǎn)是正確的,都能保證該先驗(yàn)值的有效性。
[0165]所述處理模塊200根據(jù)該條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出整體先驗(yàn)值,并根據(jù)所述整體先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0166]具體的,結(jié)合圖5所示,所述處理模塊200通過(guò)第一數(shù)據(jù)線的各部分中值,即所述第一數(shù)據(jù)線的各部分的先驗(yàn)值,可以確定相鄰數(shù)據(jù)線各部分位移估計(jì)點(diǎn)所在區(qū)域的搜索范圍。相應(yīng)的,所述處理模塊200計(jì)算與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線中所有位移估計(jì)點(diǎn)的絕對(duì)位移值時(shí),其使用所述第一數(shù)據(jù)線的先驗(yàn)值的算法有很多種,在本實(shí)施方式中,所述處理模塊200采用所述第一數(shù)據(jù)線的N+1部分的先驗(yàn)值作為與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線的位移最大值 這一算法,以求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,所述處理模塊200繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0167]具體的,所述處理模塊200采用此種算法,求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上各部分位移估計(jì)點(diǎn)的絕對(duì)位移值所分布的區(qū)間為:第I部分的區(qū)間為:[0,6.3],第2部分[4.2,8.5],第3部分[6.3,10.2],第4部分的區(qū)間為:[8.5,10.3],第5部分大于
10.3。參照各部分所處于的相關(guān)區(qū)域,所述處理模塊200直接求取與所述第一數(shù)據(jù)線相鄰的當(dāng)前數(shù)據(jù)線上的各位移估計(jì)點(diǎn)的絕對(duì)位移值。并按照此種方法,所述處理模塊200繼續(xù)順序求取與當(dāng)前數(shù)據(jù)線相鄰的其他數(shù)據(jù)線上的所有位移估計(jì)點(diǎn)的絕對(duì)位移值,直到獲取整個(gè)數(shù)據(jù)幀上所有位移估計(jì)點(diǎn)的絕對(duì)位移值。
[0168]所述處理模塊200還用于:依據(jù)獲取的所述絕對(duì)位移值沿軸向方向進(jìn)行梯度運(yùn)算,并獲取相應(yīng)的應(yīng)變數(shù)據(jù),以偽彩的形式顯示應(yīng)變數(shù)據(jù)圖像。
[0169]與現(xiàn)有技術(shù)相比,本發(fā)明的超聲彈性成像位移值的計(jì)算方法及裝置,能夠針對(duì)不同的應(yīng)變大小自適應(yīng)地調(diào)整搜索范圍,大大減少了計(jì)算量,同時(shí)具有很高的魯棒性來(lái)避免錯(cuò)誤擴(kuò)散,保證了彈性成像的實(shí)時(shí)性和準(zhǔn)確性。
[0170]為了描述的方便,描述以上裝置時(shí)以功能分為各種模塊分別描述。當(dāng)然,在實(shí)施本申請(qǐng)時(shí)可以把各模塊的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)。
[0171]通過(guò)以上的實(shí)施方式的描述可知,本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請(qǐng)可借助軟件加必需的通用硬件平臺(tái)的方式來(lái)實(shí)現(xiàn)。基于這樣的理解,本申請(qǐng)的技術(shù)方案本質(zhì)上或者說(shuō)對(duì)現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來(lái),該計(jì)算機(jī)軟件產(chǎn)品可以保存在保存介質(zhì)中,如ROM/RAM、磁碟、光盤等,包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī),信息推送服務(wù)器,或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請(qǐng)各個(gè)實(shí)施方式或者實(shí)施方式的某些部分所述的方法。
[0172]以上所描述的裝置實(shí)施方式僅僅是示意性的,其中所述作為分離部件說(shuō)明的模塊可以是或者也可以不是物理上分開(kāi)的,作為模塊顯示的部件可以是或者也可以不是物理模塊,即可以位于一個(gè)地方,或者也可以分布到多個(gè)網(wǎng)絡(luò)模塊上??梢愿鶕?jù)實(shí)際的需要選擇其中的部分或者全部模塊來(lái)實(shí)現(xiàn)本實(shí)施方式方案的目的。本領(lǐng)域普通技術(shù)人員在不付出創(chuàng)造性勞動(dòng)的情況下,即可以理解并實(shí)施。
[0173]本申請(qǐng)可用于眾多通用或?qū)S玫挠?jì)算系統(tǒng)環(huán)境或配置中。例如:個(gè)人計(jì)算機(jī)、信息推送服務(wù)器計(jì)算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理模塊系統(tǒng)、基于微處理模塊的系統(tǒng)、置頂盒、可編程的消費(fèi)電子設(shè)備、網(wǎng)絡(luò)PC、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計(jì)算環(huán)境等等。
[0174]本申請(qǐng)可以在由計(jì)算機(jī)執(zhí)行的計(jì)算機(jī)可執(zhí)行指令的一般上下文中描述,例如程序模塊。一般地,程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對(duì)象、組件、數(shù)據(jù)結(jié)構(gòu)等等。也可以在分布式計(jì)算環(huán)境中實(shí)踐本申請(qǐng),在這些分布式計(jì)算環(huán)境中,由通過(guò)通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來(lái)執(zhí)行任務(wù)。在分布式計(jì)算環(huán)境中,程序模塊可以位于包括保存設(shè)備在內(nèi)的本地和遠(yuǎn)程計(jì)算機(jī)保存介質(zhì)中。
[0175]應(yīng)當(dāng)理解,雖然本說(shuō)明書(shū)按照實(shí)施方式加以描述,但并非每個(gè)實(shí)施方式僅包含一個(gè)獨(dú)立的技術(shù)方案,說(shuō)明書(shū)的這種敘述方式僅僅是為清楚起見(jiàn),本領(lǐng)域技術(shù)人員應(yīng)當(dāng)將說(shuō)明書(shū)作為一個(gè)整體,各實(shí)施 方式中的技術(shù)方案也可以經(jīng)適當(dāng)組合,形成本領(lǐng)域技術(shù)人員可以理解的其他實(shí)施方式。
[0176]上文所列出的一系列的詳細(xì)說(shuō)明僅僅是針對(duì)本發(fā)明的可行性實(shí)施方式的具體說(shuō)明,它們并非用以限制本發(fā)明的保護(hù)范圍,凡未脫離本發(fā)明技藝精神所作的等效實(shí)施方式或變更均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。
【權(quán)利要求】
1.一種超聲彈性成像位移值的計(jì)算方法,其特征在于,所述方法包括以下步驟: S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線; S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分; S3、選取至少一部分,計(jì)算在該部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值,并將每條數(shù)據(jù)線的先驗(yàn)值作為計(jì)算同一條數(shù)據(jù)線上某一部分或某幾部分的絕對(duì)位移值的計(jì)算基礎(chǔ)。
2.根據(jù)權(quán)利要求1所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,所述回波信號(hào)劃分的多個(gè)部分中包括第一部分和第二部分,所述S3步驟包括: 選取第一部分,計(jì)算所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第一部分先驗(yàn)值; 根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值作為同一條數(shù)據(jù)線上第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
3.根據(jù)權(quán)利要求2所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,所述第一部分和所述第二部分之間間隔有至少一個(gè)部分。
4.根據(jù)權(quán)利要求3所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,所述回波信號(hào)劃分的多個(gè)部分還包括第三部分,所述第三部分在軸向上位于所述第一部分和第二部分之間; 所述S3步驟包括: 根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的第二部分先驗(yàn)值; 根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將每條數(shù)據(jù)線的第一部分先驗(yàn)值和第二部分先驗(yàn)值形成同一條數(shù)據(jù)線上第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出同一條數(shù)據(jù)線上第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
5.根據(jù)權(quán)利要求1所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,“根據(jù)每條數(shù)據(jù)線上的絕對(duì)位移值計(jì)算出每條數(shù)據(jù)線的先驗(yàn)值”步驟具體包括: 濾除已計(jì)算出的絕對(duì)位移值中的噪聲; 取濾除噪聲后的絕對(duì)位移值的中值或取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
6.一種超聲彈性成像位移值的計(jì)算方法,其特征在于,所述方法包括以下步驟: S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線; S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn); S3、通過(guò)搜索算法計(jì)算至少一條數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值; S4、根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出先驗(yàn)值,并根據(jù)所述先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
7.根據(jù)權(quán)利要求6所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,“根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出先驗(yàn)值”步驟具體包括: 濾除已計(jì)算出的絕對(duì)位移值中的噪聲; 取濾除噪聲后的絕對(duì)位移值的中值或取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
8.根據(jù)權(quán)利要求6所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,“根據(jù)所述先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值”步驟具體包括: 沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分,根據(jù)相鄰數(shù)據(jù)線所對(duì)應(yīng)的部分在軸向上的位置關(guān)系及其先驗(yàn)值形成搜索區(qū)間,以通過(guò)搜索算法計(jì)算出相鄰數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
9.一種超聲彈性成像位移值的計(jì)算方法,其特征在于,所述方法包括以下步驟: S1、獲取壓縮組織前后的兩幀回波信號(hào),每幀回波信號(hào)包括多條數(shù)據(jù)線; S2、在每條數(shù)據(jù)線上選取多個(gè)位移估計(jì)點(diǎn),并沿?cái)?shù)據(jù)線的軸向?qū)⑺龌夭ㄐ盘?hào)劃分為多個(gè)部分; S3、通過(guò)搜索算法計(jì)算某一條數(shù)據(jù)線中某一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值; S4、根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值,并根據(jù)所述部分先驗(yàn)值計(jì)算該條數(shù)據(jù)線的另一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,直至獲得該條數(shù)據(jù)線所有位移估計(jì)點(diǎn)的絕對(duì)位移值; S5、根據(jù)該條數(shù)據(jù)線上每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出整體先驗(yàn)值,并根據(jù)所述整體先驗(yàn)值計(jì)算相鄰數(shù)據(jù)線上的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
10.根據(jù)權(quán)利要求9所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,所述回波信號(hào)劃分的多個(gè)部分中包括第一部分和第二部分,所述S3步驟包括: 選取某一條數(shù)據(jù)線的第一部分,計(jì)算所述第一部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值,并根據(jù)所述第一部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出第一部分先驗(yàn)值; 根據(jù)第一部分和第二部分在軸向上的位置關(guān)系,將第一部分先驗(yàn)值作為第二部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間的一端值,以通過(guò)搜索算法計(jì)算出第二部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
11.根據(jù)權(quán)利要求10所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,所述第一部分和所述第二部分之間間隔有至少一個(gè)部分。
12.根據(jù)權(quán)利要求11所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,所述回波信號(hào)劃分的多個(gè)部分還包括第三部分,所述第三部分在軸向上位于所述第一部分和第二部分之間,所述S3步驟包括: 根據(jù)所述第二部分內(nèi)的每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出第二部分先驗(yàn)值; 根據(jù)第一部分、第二部分、第三部分在軸向上的位置關(guān)系,將第一部分先驗(yàn)值和第二部分先驗(yàn)值形成第三部分內(nèi)位移估計(jì)點(diǎn)的搜索區(qū)間,以通過(guò)搜索算法計(jì)算出第三部分內(nèi)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值。
13.根據(jù)權(quán)利要求9所述的超聲彈性成像位移值的計(jì)算方法,其特征在于,“根據(jù)每個(gè)位移估計(jì)點(diǎn)的絕對(duì)位移值計(jì)算出部分先驗(yàn)值”步驟具體包括: 濾除已計(jì)算出的絕對(duì)位移值中的噪聲; 取濾除噪聲后的絕對(duì)位移值的中值或取整后出現(xiàn)次數(shù)最多的絕對(duì)位移值作為先驗(yàn)值。
【文檔編號(hào)】A61B8/08GK103829975SQ201410096941
【公開(kāi)日】2014年6月4日 申請(qǐng)日期:2014年3月17日 優(yōu)先權(quán)日:2014年3月17日
【發(fā)明者】周一蘭, 凌濤, 陳惠人 申請(qǐng)人:飛依諾科技(蘇州)有限公司