[發(fā)明專利]聲輻射力脈沖成像估算方法和系統(tǒng)

有效

• 申請?zhí)枺?/td>: 201310404937.0
• 申請日: 2013-09-06
• 公開（公告）號: CN103431874A
• 公開（公告）日: 2013-12-11
• 發(fā)明（設計）人: 鄭海榮;馮歌;王叢知;曾成志;楊戈;曾博
• 申請（專利權(quán)）人: 中國科學院深圳先進技術(shù)研究院
• 主分類號: A61B8/00
• 分類號: A61B8/00;A61B8/08
• 代理公司: 廣州華進聯(lián)合專利商標代理有限公司 44224
• 代理人: 吳平

1.一種聲輻射力脈沖成像估算方法，包括：

數(shù)據(jù)劃分步驟，讀取通過聲輻射力脈沖成像所采集的輸入數(shù)據(jù)，并將所述輸入數(shù)據(jù)分成預定數(shù)量數(shù)據(jù)矩陣；

數(shù)據(jù)濾波變換步驟，分別對每個數(shù)據(jù)矩陣進行前向濾波處理和逆向濾波處理，并對濾波處理后的數(shù)據(jù)矩陣進行希爾伯特變換得到解析數(shù)據(jù)；

位移求取步驟，對所述解析數(shù)據(jù)進行互相關(guān)運算，求取得到預定數(shù)量位移矩陣；

位移合并步驟，對每個位移矩形進行時間累加，取累積和最大的列位移，再將預定數(shù)量位移矩陣抽取的列位移合并，生成新的位移矩陣；

插值步驟，對得到的新的位移矩陣進行時間方向上的三次樣條插值；

擬合步驟，對插值結(jié)果進行拉東變換，得到拉東變換矩陣；

估算步驟，獲取所述聲輻射力脈沖成像的參數(shù)，根據(jù)所述參數(shù)及拉東變換矩陣求取剪切波速率。

2.根據(jù)權(quán)利要求1所述的聲輻射力脈沖成像估算方法，其特征在于，所述數(shù)據(jù)濾波變換步驟中，采用多個流按順序依次對預定數(shù)量數(shù)據(jù)矩陣進行前向濾波處理、逆向濾波處理和希爾伯特變換處理。

3.根據(jù)權(quán)利要求1所述的聲輻射力脈沖成像估算方法，其特征在于，所述數(shù)據(jù)濾波變換步驟中前向濾波處理的公式為：

當i≤iParalen時：

tOutPut(k,i)=Σj=0i(param[j]*input(k,i-j))+Σj=i+1iParalen-2(param[j]*input(k,0))]]>

當i＞iParalen時：

tOutPut(k,i)=Σj=0iParalen-2(param[j]*input(k,i-j))]]>

其中，param是濾波器，iParalen是濾波器長度，input為待濾波矩陣，k、i分別表示待濾波矩陣的第k行和第i列，tOutPut是濾波結(jié)果矩陣；

所述數(shù)據(jù)濾波變換步驟中逆向濾波處理的公式為：

當i≥iParalen時：

output(k,Width-1-i)=Σj=0iParalen-1(Param[j]*input(k,Width-1+j-i))]]>

其中，input為待濾波矩陣，k、i分別表示待濾波矩陣的第k行和第i列，output為逆向濾波的濾波結(jié)果矩陣，Width為濾波器寬度；

所述前向濾波處理的內(nèi)核函數(shù)和逆向濾波處理均采用圖形處理器并行處理，所述圖形處理器調(diào)用的內(nèi)核函數(shù)的blocks數(shù)量和threads數(shù)量相同，所述blocks數(shù)量為待濾波矩陣的行數(shù)，所述threads數(shù)量為待濾波矩陣的列數(shù)。