[發(fā)明專(zhuān)利]一種X射線造影方法及基于X射線造影的動(dòng)態(tài)閱片方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201810203889.1
• 申請(qǐng)日: 2018-03-13
• 公開(kāi)（公告）號(hào): CN108852385B
• 公開(kāi)（公告）日: 2022-03-04
• 發(fā)明（設(shè)計(jì)）人: 肖體喬;王飛翔;陳榮昌;王玉丹;杜國(guó)浩;鄧彪;謝紅蘭
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中國(guó)科學(xué)院上海應(yīng)用物理研究所
• 主分類(lèi)號(hào): A61B6/00
• 分類(lèi)號(hào): A61B6/00;G02B27/02
• 代理公司: 上海智信專(zhuān)利代理有限公司 31002
• 代理人: 鄧琪

1.一種X射線造影方法，包括：

步驟S1：以采樣頻率Fs獲取N幀的原始的造影圖像序列，作為當(dāng)前的造影圖像序列；

步驟S2：在原始的造影圖像序列的其中一幀上選取血管特征點(diǎn)和背景特征點(diǎn)，提取這兩個(gè)特征點(diǎn)的時(shí)序信號(hào)并將其時(shí)序信號(hào)傅里葉變換為頻域信號(hào)；

步驟S3：對(duì)血管特征點(diǎn)和背景特征點(diǎn)的頻域信號(hào)進(jìn)行特征分析，依次選取這兩個(gè)特征點(diǎn)在相同頻率下的最大的至少兩個(gè)模值差，第i大的模值差所對(duì)應(yīng)的頻率作為特征頻率xi*Fs/N Hz；

步驟S4：在當(dāng)前的造影圖像序列中遍歷造影區(qū)域的所有像素點(diǎn)，以0Hz和其中一個(gè)特征頻率xi*Fs/N Hz為基準(zhǔn)對(duì)每個(gè)像素點(diǎn)的時(shí)序信號(hào)做傅里葉帶通濾波處理，并通過(guò)得到的信號(hào)計(jì)算每個(gè)像素點(diǎn)的調(diào)制深度MD_i；

步驟S5：輸出以調(diào)制深度MD_i為造影成像參量的造影圖像；

其中xi為正整數(shù)；

所述每個(gè)像素點(diǎn)的調(diào)制深度MD_i通過(guò)下式計(jì)算：

其中，S₀(x,y,t)、S_xi(x,y,t)分別是對(duì)每個(gè)像素點(diǎn)的時(shí)序信號(hào)做傅里葉帶通濾波處理后獲得的直流信號(hào)和對(duì)應(yīng)xi*Fs/N Hz頻率的交流信號(hào)；Mean(abs(S_xi(x,y,t)))為對(duì)應(yīng)xi*Fs/N Hz頻率的交流信號(hào)的絕對(duì)值的平均值；Mean(abs(S₀(x,y,t)))為直流信號(hào)的絕對(duì)值的平均值。

2.根據(jù)權(quán)利要求1所述的X射線造影方法，其特征在于，所述步驟S3中所選取的模值差為2-4個(gè)，且所述特征頻率的個(gè)數(shù)等于模值差的個(gè)數(shù)。

3.根據(jù)權(quán)利要求1所述的X射線造影方法，其特征在于，所述步驟S4還包括：在遍歷造影區(qū)域的所有像素點(diǎn)之前，先對(duì)其每一幀圖像做去噪處理，該去噪處理包括中值去噪、高斯去噪或維納去噪。

4.一種基于X射線造影的動(dòng)態(tài)閱片方法，包括：

步驟S1’：采用權(quán)利要求1所述X射線造影方法，得到特征頻率xi*Fs/NHz，并得到以調(diào)制深度MD_i為造影成像參量的造影圖像，作為原始的造影圖像序列的前N幀所對(duì)應(yīng)的造影圖像，取j＝1；

步驟S2’：選取原始的造影圖像序列的前N-j幀作為當(dāng)前的造影圖像序列，并以步驟S1’得到的特征頻率計(jì)算所述當(dāng)前的造影圖像序列所對(duì)應(yīng)的特征頻率xi*Fs/(N-j)Hz；

步驟S3’：在當(dāng)前的造影圖像序列中遍歷造影區(qū)域的所有像素點(diǎn)，以0Hz和步驟S2’得到的特征頻率xi*Fs/(N-j)Hz為基準(zhǔn)對(duì)每個(gè)像素點(diǎn)的時(shí)序信號(hào)做傅里葉帶通濾波處理，并通過(guò)得到的信號(hào)計(jì)算每個(gè)像素點(diǎn)的調(diào)制深度MD_i；得到以調(diào)制深度MD_i為造影成像參量的造影圖像，作為原始的造影圖像序列的前N-j幀所對(duì)應(yīng)的造影圖像；

步驟S4’：當(dāng)N-j為n時(shí)，依次播放所述前n、n+1…、N-1、N-2、N幀所對(duì)應(yīng)的造影圖像；否則將j+1更新為j并回到步驟S2’。

5.根據(jù)權(quán)利要求4所述的基于X射線造影的動(dòng)態(tài)閱片方法，其特征在于，當(dāng)i為1時(shí)，所述以調(diào)制深度MD₁為造影成像參量的造影圖像為血管造影圖像。