一种X射线冠脉造影图像降噪方法 .pdfVIP

(19)中华人民共和国国家知识产权局

(12)发明专利说明书

(10)申请公布号CN102663710A

(43)申请公布日2012.09.12

(21)申请号CN201210132698.3

(22)申请日2012.04.28

(71)申请人大连海事大学

地址116026辽宁省大连市高新园区凌海路1号

(72)发明人尹清波陈燕鲁明羽申丽然陈飞孟祥申

(74)专利代理机构大连东方专利代理有限责任公司

代理人李洪福

(51)Int.CI

G06T5/00

权利要求说明书说明书幅图

(54)发明名称

一种X射线冠脉造影图像降噪方法

(57)摘要

本发明公开了一种X射线冠脉造影

图像降噪方法，包括以下步骤：对X射线

冠脉造影图像I进行二维经验模态分解；

分析经验模态分解得到的各层本征模态函

数，进行去噪；利用能够代表冠脉造影图

像血管结构的本征模态函数进行加权求和

来重构去噪后的冠脉造影图像。本发明不

依赖先验知识进行对冠脉造影图像自适应

分解，根据图像特点自适应产生随机噪声

去除阈值。本发明利用经验模态分解方法

对X射线冠脉造影图像进行自适应分解，

对得到的各IMF分量进行去噪，并选择能

够代表血管结构的IMF高频分量进行造影

图像重构即可达到去除噪声的目的，且在

去噪过程中无须先验知识、能够根据造影

图像特点进行自适应降噪。

法律状态

法律状态公告日法律状态信息法律状态

权利要求说明书

1.一种X射线冠脉造影图像降噪方法，其特征在于：包括以下步骤：

A、对X射线冠脉造影图像I进行二维经验模态分解

A1、利用输入的X射线冠脉造影图像I作为分解初始值，设定当前层为i＝1，

第i层本征模态函数的估计为rsubi/sub＝I；

A2、令第i层本征模态函数的第j次估计量为hsubj/sub＝rsubi/sub，j＝1；

A3、求取hsubj/sub的局部极大值和极小值；

A4、对所获得的极大和极小值分别用三次样条插值得到上包络envsubmax/sub

和下包络envsubmin/sub；

A5、计算上包络envsubmax/sub和下包络envsubmin/sub的均值，

A6、从第i层本征模态函数的第j次估计量hsubj/sub中减去m，

hsubj+1/sub＝hsubj/sub-m；

A7、从第i层本征模态函数的估计次数j＝j+1；

A8、判断：如果||hsubj/sub||＜ε或j＞＝4不成立，则转到步骤A3；如果

||hsubj/sub‖＜ε或j＞＝4成立，则获得X射线冠脉造影图

像I的第i层本征模态函数Fsubi/sub＝hsubj/sub；

A9、得到第i+1层本征模态函数的初始估计

A10、本征模态函数的层序数i＝i+1；

A11、如果rsubi/sub为单调函数，则R＝rsubi/sub为残余分量，分解过程

结束；否则转到步骤A2；

B、分析经验模态分解得到的各层本征模态函数，进行去噪

B1、分析X射线冠脉造影图像I的各层本征模态函数Fsubi/sub，计算强度直

方图得到强度分布点数最大的强度值histsubi/sub；

B2、计算各层本征模态函数Fsubi/sub的强度方差σsubi/sub；

B3、将方差σsubi/sub或强度值histsubi/sub的最大者作为第i层本征模态

函数的阈值Thsubi/sub；

B4、对每一层本征模态函数利用阈值Thsubi/sub将小于阈值的点置0；

C、利用能够代表冠脉造影图像血管结构的本征模