[發明專利]一種數字減影血管造影方法

有效

• 申請號: 201410449767.2
• 申請日: 2014-09-04
• 公開（公告）號: CN104156931B
• 公開（公告）日: 2017-03-29
• 發明（設計）人: 曲建明
• 申請（專利權）人: 成都金盤電子科大多媒體技術有限公司
• 主分類號: G06T5/50
• 分類號: G06T5/50
• 代理公司: 成都金英專利代理事務所(普通合伙)51218
• 代理人: 袁英

1.一種新的數字減影血管造影方法，其特征在于：它包括以下步驟：

S1：圖像序列對數變換：對減影過程中的基象和造影原像分別進行對數變換處理，能使減影圖像的血管更清晰，對第i幀圖像序列I_i中的每個像素進行對數變換，得到第i幀輸出圖像P_i，其計算公式為：

Pi=Ln(IiBW)]]>

式中，BW－圖像幀I_i的像素最大值；

S2：計算背景圖像模板：為了抑制噪聲和造影劑給血管圖像帶來的影響，將圖像序列P_i進行加權求和，而第i幀圖像序列P_i像素值的權重大小與它和上一幀圖像序列P_i-1同位置像素值的差成反比，背景圖像模板mask_i的計算公式如下：

dfi(m,n)=α·|Pi(m,n)-Pi-1(m,n)|mask0(m,n)=P0(m,n)maski(m,n)=(1-dfi(m,n))·maski-1(m,n)+Pi(m,n)·dfi(m,n)]]>

式中，α－常數，且α∈(0,1)；

m、n－圖像幀中像素的位置；

P_i(m,n)－第i幀圖像P_i中位于m，n坐標的像素值；

S3：計算初步血管圖像：對第i幀圖像序列P_i和其背景圖像模板mask_i的差進行指數變換，得到第i幀的初步血管圖像V_i，其計算公式如下：

Vi=e-|Pi-maski|;]]>

S4：計算血管增強模板：對第i幀圖像序列P_i和上一幀圖像序列P_i-1的差進行指數變換，得到第i幀的血管增強模板D_i′，其計算公式如下：

Di=e-|Pi-Pi-1|Di′(m,n)=Min(Di(m,n),Di+1(m,n));]]>

S5：中值濾波處理：對第i幀血管增強模板D_i′進行3×3的中值濾波處理，濾除增強模像中的噪聲，得到第i幀的血管三維增強模板D_i″，其計算公式如下：

D_i″＝medfilt(D_i′,[3?3])；

S6：血管減影圖像增強：用第i幀的血管三維增強模板D_i″對第i幀的初步血管圖像V_i進行圖像增強，得到第i幀的血管減影圖像V_i’(m,n)，其計算公式為：

Vi′(m,n)=Vi(m,n)·(1+Di′′(m,n)Max(Di′′))]]>

式中，Max(D_i″)－第i幀的血管三維增強模板D_i″的最大像素值；

S7：血管圖像合成：對每一幀血管減影圖像V_i’(m,n)進行累加，得到血管合成影像V_t，其計算公式為：

Vt＝∑V_i′；

S8：圖像位寬變換處理：首先計算圖像序列的最大值和最小值，然后根據圖像序列的最大值和最小值計算圖像像素值的最大值和最小值，再根據圖像像素值的最大值和最小值對圖像序列進行限幅處理，最后對限幅后的圖像序列進行平移拉伸，即得到最終輸出顯示的圖像序列；

所述的圖像位寬變換處理包括血管減影圖像的圖像位寬變換處理和血管合成圖像的圖像位寬變換處理；

S9：圖像數據輸出：通過計算機顯示器輸出顯示血管減影圖像和血管合成圖像。