[發(fā)明專利]一種超聲彈性成像性能增強(qiáng)方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210486012.0
• 申請(qǐng)日: 2012-11-26
• 公開（公告）號(hào): CN102920480A
• 公開（公告）日: 2013-02-13
• 發(fā)明（設(shè)計(jì)）人: 崔少國(guó);彭彩碧;黃賢英
• 申請(qǐng)（專利權(quán)）人: 重慶理工大學(xué)
• 主分類號(hào): A61B8/00
• 分類號(hào): A61B8/00
• 代理公司: 重慶博凱知識(shí)產(chǎn)權(quán)代理有限公司 50212
• 代理人: 穆祥維

1.一種超聲彈性成像性能增強(qiáng)方法，其特征在于，該方法包括如下步驟：

（1）、在準(zhǔn)靜態(tài)壓縮組織過(guò)程中，超聲探頭發(fā)射超聲波，取壓縮前后兩幀RF回波，轉(zhuǎn)換成基帶IQ信號(hào)；

（2）、對(duì)壓縮前信號(hào)，以位移估計(jì)點(diǎn)為中心取二維運(yùn)動(dòng)追蹤窗；

（3）、第一行追蹤窗初始迭代位移均設(shè)為0，即d’_m,1,0=?0，其中下標(biāo)m表示列號(hào)，1表示第一行，0表示初始迭代；

（4）、使用二維相位根搜索位移估計(jì)算法，迭代產(chǎn)生第一行軸向位移估計(jì)值，方法如下：

?④

上式中：d_m,n,l表示第m列第n行追蹤窗進(jìn)行l次迭代后的位移估計(jì)值；為側(cè)向窗寬；為軸向窗長(zhǎng)；和分別是相鄰估計(jì)窗之間側(cè)向和軸向的間隔；m，n表示估計(jì)窗的位置下標(biāo)；S_1b和S_2b分別表示壓縮前和壓縮后的基帶信號(hào)；_0?表示中心角頻率；*代表共軛；j是復(fù)信號(hào)；arg是求反正切操作；

（5）、計(jì)算該位移對(duì)應(yīng)的相關(guān)系，為下一行初始迭代位移選擇做準(zhǔn)備；使用復(fù)相關(guān)歸一化相關(guān)系數(shù)法；

⑤

上式中：N是常數(shù)，是設(shè)定的位移估計(jì)迭代次數(shù)；其它參數(shù)涵義與步驟（4）中公式④的參數(shù)涵義相同；

（6）、對(duì)于第二行位移估計(jì)，先從它相鄰上一行中以它為軸左右各取L列的2L+1個(gè)已估計(jì)位移中選擇相關(guān)系數(shù)最高對(duì)應(yīng)的位移作為它的初始位移，0≤L≤2，即：?

????????????????⑥

其中：b滿足，其中0<=b<=L，m+b是上一行2L+1范圍內(nèi)已估計(jì)位移對(duì)應(yīng)相關(guān)系數(shù)最高的列號(hào)，表示求最大值函數(shù)；

（7）、使用步驟（4）中的二維相位根搜索位移估計(jì)算法迭代產(chǎn)生第二行軸向位移估計(jì)值；

（8）、同第二行軸向位移估計(jì)方法一樣，估計(jì)出第三行到最后一行的軸向位移；

（9）、對(duì)于步驟（3）～步驟（8）產(chǎn)生的位移值，使用二維窗位移位置估計(jì)改進(jìn)算法估計(jì)其對(duì)應(yīng)的軸向位置；該二維窗位移位置估計(jì)改進(jìn)算法如下：?

??????⑦

式中，表示第m列第n行追蹤窗對(duì)應(yīng)位移的軸向位置坐標(biāo)，和分別是側(cè)向和軸向相鄰追蹤窗之間的間隔，和分別是追蹤窗的側(cè)向窗寬和軸向窗長(zhǎng)，是追蹤窗的軸向位移，是權(quán)值，，其中和分別是壓縮前和壓縮后信號(hào)的包絡(luò)；

（10）、根據(jù)位移位置估計(jì)，使用線性插值法將位移估計(jì)值修正到追蹤窗中心點(diǎn)，方法如下：

當(dāng)<時(shí)，?????⑧

當(dāng)>時(shí)，????⑨

上式中，表示第m列第n行追蹤窗修正后的位移，表示該窗由步驟（4）得到修正前的位移，是位移對(duì)應(yīng)的縱向坐標(biāo)，是追蹤窗中心點(diǎn)的縱坐標(biāo)；

（11）、使用梯度法進(jìn)行縱向應(yīng)變估計(jì)；

（12）、對(duì)應(yīng)變估計(jì)矩陣進(jìn)行下采樣和灰度映射，掃描轉(zhuǎn)換生成彈性圖像。