[發(fā)明專利]一種CT冠脈圖像與XA造影圖像的圖像融合方法及系統(tǒng)

有效

• 申請(qǐng)?zhí)枺?/td>: 201210592283.4
• 申請(qǐng)日: 2012-12-28
• 公開（公告）號(hào): CN103914814A
• 公開（公告）日: 2014-07-09
• 發(fā)明（設(shè)計(jì)）人: 王民漢;周宇;李曉光
• 申請(qǐng)（專利權(quán)）人: 北京思創(chuàng)貫宇科技開發(fā)有限公司
• 主分類號(hào): G06T5/50
• 分類號(hào): G06T5/50
• 代理公司: 北京律誠(chéng)同業(yè)知識(shí)產(chǎn)權(quán)代理有限公司 11006
• 代理人: 祁建國(guó);梁揮

1.一種CT冠脈圖像與XA造影圖像的圖像融合方法，其特征在于，包括：

步驟1，提取XA圖像的血管圖像，對(duì)所述血管圖像進(jìn)行擴(kuò)散處理，得到XA擴(kuò)散圖像；

步驟2，提取CT冠脈圖像，對(duì)所述CT冠脈圖像進(jìn)行二維投影，得到二維CT圖像，然后對(duì)所述二維CT圖像進(jìn)行膨脹處理，得到CT膨脹圖像；

步驟3，以XA擴(kuò)散圖像作為參考圖像，以CT膨脹圖像做浮動(dòng)圖像，根據(jù)相似性測(cè)度函數(shù)對(duì)CT冠脈圖像進(jìn)行處理，獲得CT模擬XA圖像；

步驟4，將配準(zhǔn)后的CT模擬XA圖像加載到XA圖像上進(jìn)行融合顯示。

2.如權(quán)利要求1所述的圖像融合方法，其特征在于，所述步驟1包括：

步驟11，從XA圖像中選擇一幀圖像后，對(duì)所述一幀圖像進(jìn)行增強(qiáng)處理；

步驟12，對(duì)增強(qiáng)處理后的一幀圖像采用分水嶺算法進(jìn)行血管分割，獲得所述血管圖像；

步驟13，對(duì)所述血管圖像進(jìn)行距離變換，公式為：

XA_Diffusion(i)＝max(0,Thread-DT(i))，其中DT(i)為所述血管圖像距離變換后圖像上每一點(diǎn)的數(shù)值，XADiffusion是XA擴(kuò)散圖像，Thread是XADiffusion的閾值；

步驟14，利用閾值對(duì)所述血管圖像進(jìn)行距離變換獲得所需的擴(kuò)散圖像。

3.如權(quán)利要求2所述的圖像融合方法，其特征在于，所述步驟11中采用多尺度Hessian卷積對(duì)所述血管圖像進(jìn)行增強(qiáng)，其公式為

V(s)=0ifλ2>0,exp(-RB22β2)(1-exp(-S22c2))]]>

其中|λ₁|≤|λ₂|，R_B為λ₁和λ₂的比值，S為λ₁和λ₂的F-范數(shù)，為了來抑制圖像背景中的噪聲，參數(shù)β和c是調(diào)節(jié)R_B和S作用的加權(quán)因子，V(s)為單一尺度下圖像增強(qiáng)的數(shù)值。

4.如權(quán)利要求1所述的圖像融合方法，其特征在于，所述步驟2包括：

步驟21，采用區(qū)域增長(zhǎng)法對(duì)心臟CT圖像進(jìn)行冠脈提取，獲得冠脈血管數(shù)據(jù)；

步驟22，根據(jù)XA圖像中的主角度與副角度信息以及發(fā)射源和接收板的位置信息,確定CT冠脈圖像進(jìn)行投影的相關(guān)參數(shù)；

步驟23，所述冠脈血管數(shù)據(jù)和所述相關(guān)參數(shù)在剛性配準(zhǔn)矩陣的作用下，得到二維CT圖像。