[發(fā)明專利]實時超聲彈性成像方法和系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201410326830.3
• 申請日: 2014-07-09
• 公開（公告）號: CN105266849A
• 公開（公告）日: 2016-01-27
• 發(fā)明（設(shè)計）人: 孫新;趙明昌
• 申請（專利權(quán)）人: 無錫祥生醫(yī)學(xué)影像有限責(zé)任公司
• 主分類號: A61B8/08
• 分類號: A61B8/08
• 代理公司: 無錫市大為專利商標(biāo)事務(wù)所(普通合伙) 32104
• 代理人: 曹祖良

1.一種實時超聲彈性成像方法，其特征在于，包括下述步驟：

S10.對生物組織進行緩速擠壓時對生物組織目標(biāo)區(qū)域進行超聲波掃查并接收回波信號；

S20.對步驟S10中接收的回波信號進行處理，形成線數(shù)據(jù)；

S30.對步驟S20得到的線數(shù)據(jù)進行處理，形成生物組織受外力緩速擠壓后產(chǎn)生形變的不同狀態(tài)下的B模式圖像；

S40.對生物組織受外力緩速擠壓后產(chǎn)生形變的不同狀態(tài)下的兩幀B模式圖像選定兩個在不同幀圖像的同一位置區(qū)域，即選定兩幀圖像的ROI區(qū)域；

S50.通過第一CPU數(shù)據(jù)處理模塊對該兩幀圖像的ROI區(qū)域使用光流法求取位移場；并在計算過程中將含重復(fù)計算的子步驟的參數(shù)信息傳輸至第一GPU數(shù)據(jù)緩存模塊，進行緩存；

S60.第一GPU數(shù)據(jù)緩存模塊將接收到的子步驟的參數(shù)信息傳輸給各自的GPU數(shù)據(jù)處理模塊進行數(shù)據(jù)處理，S50中的含重復(fù)計算的各子步驟的運算分配給至少一個GPU核工作單元，由至少一個GPU核工作單元組成GPU工作組；

S70.將每個GPU工作組計算需要的參數(shù)信息映射到每個GPU工作組的顯存中，等待所有GPU工作組的參數(shù)信息映射完成，然后各GPU工作組進行步驟S50中的各子步驟的運算處理；

S80.將GPU工作組的運算結(jié)果信息傳輸至第二GPU數(shù)據(jù)緩存模塊；

S90.第二GPU數(shù)據(jù)緩存模塊將運算結(jié)果信息傳輸至第二CPU數(shù)據(jù)處理模塊，第二CPU數(shù)據(jù)處理模塊讀出運算結(jié)果信息并帶入光流法運算過程，最終求出ROI區(qū)域的位移場即光流場；

S100.使用低通濾波器對軸向光流場求取應(yīng)變，得到圖像ROI區(qū)域的軸向應(yīng)變場；

S105.對經(jīng)過濾波得到的圖像ROI區(qū)域的軸向應(yīng)變場信息進行降噪聲處理；

S110.對降噪后的ROI區(qū)域軸向應(yīng)變場信息進行可視化、彩色處理得到生物組織彩色彈性圖像。

2.如權(quán)利要求1所述的實時超聲彈性成像方法，其特征在于：

所述步驟S20中，通過波束合成模塊對回波信號具體進行A/D轉(zhuǎn)換、同相位疊加處理；

所述步驟S30中，通過B模式圖像處理模塊對步驟S20得到的線數(shù)據(jù)具體進行降噪、濾波、提高信噪比處理。

3.如權(quán)利要求1所述的實時超聲彈性成像方法，其特征在于：

所述步驟S50中，含重復(fù)計算的子步驟的參數(shù)信息是指：構(gòu)建圖像金字塔、中心梯度計算、計算目標(biāo)圖像及其導(dǎo)數(shù)的形變、估計雙變量、估計位移這些計算的參數(shù)信息。