主权项 |
1.一种X射线冠脉造影图像血管增强方法,其特征在于:包括以下步骤:A、对X射线冠脉造影图像I进行二维经验模态分解A1、利用输入的X射线冠脉造影图像I作为分解初始值,设定当前层为i=1,第i层本征模态函数的估计为r<sub>i</sub>=I;A2、令第i层本征模态函数的第j次估计量为h<sub>j</sub>=r<sub>i</sub>,j=1;A3、求取h<sub>j</sub>的局部极大值和极小值;A4、对所获得的极大和极小值分别用三次样条插值得到上包络env<sub>max</sub>和下包络env<sub>min</sub>;A5、计算上包络env<sub>max</sub>和下包络env<sub>min</sub>的均值,<img file="FDA0000159059030000011.GIF" wi="409" he="98" />A6、从第i层本征模态函数的第j次估计量h<sub>j</sub>中减去m,h<sub>j+1</sub>=h<sub>j</sub>-m;A7、令第i层本征模态函数的估计次数j=j+1;A8、判断:如果||h<sub>j</sub>‖<ε或j>=4不成立,则转到步骤A3;如果||h<sub>j</sub>‖<ε或j>=4成立,则获得X射线冠脉造影图像I的第i层本征模态函数F<sub>i</sub>=h<sub>j</sub>;A9、得到第i+1层本征模态函数的初始估计<img file="FDA0000159059030000012.GIF" wi="370" he="73" />A10、令本征模态函数的层序数i=i+1;A11、如果r<sub>i</sub>为单调函数,则R=r<sub>i</sub>为残余分量,分解过程结束;否则转到步骤A2;B、分析各层经验模态分解的本征模态函数进行去噪,并利用能够代表冠脉造影图像的特征层来加权重构去除背景噪声的图像B1、分析X射线冠脉造影图像I的各层本征模态函数F<sub>i</sub>,计算强度直方图得到强度分布点数最大的强度值hist<sub>i</sub>;B2、计算各层本征模态函数F<sub>i</sub>的强度方差σ<sub>i</sub>;B3、将方差σ<sub>i</sub>或强度值hist<sub>i</sub>的最大者作为第i层本征模态函数的阈值Th<sub>i</sub>;B4、对每一层本征模态函数利用阈值Th<sub>i</sub>将小于阈值的点置0;B5、设定进行重构的本征模态函数层数为M<sub>ite</sub>=3,迭代初始值ite=1;B6、设定过程变量temp为0;B7、将去噪后的第ite层本征模态函数F<sub>i</sub>取绝对值abs(F<sub>ite</sub>);B8、过程变量<img file="FDA0000159059030000021.GIF" wi="657" he="107" />其中abs(·)为取绝对值操作;B9、令ite=ite+1;B10、如果ite<(M<sub>ite</sub>+1),则转到B8;否则重构过程结束,输出重构冠脉造影图像I<sub>re</sub>=temp;C、对重构的图像进行Hessian矩阵分解,利用血管测度函数或血管特征拟合函数达到增强血管结构的目的C1、构造多尺度高斯滤波器集<img file="FDA0000159059030000022.GIF" wi="88" he="44" />σ<sub>n</sub>为高斯滤波器的方差,σ<sub>n</sub>∈[σ<sub>1</sub>,...,σ<sub>N</sub>];C2、对多尺度高斯滤波器集<img file="FDA0000159059030000023.GIF" wi="58" he="45" />求二阶导数,用<img file="FDA0000159059030000024.GIF" wi="58" he="45" />的二阶导数与重构后的冠脉造影图像I<sub>re</sub>进行卷积得到<img file="FDA0000159059030000025.GIF" wi="376" he="63" />和<img file="FDA0000159059030000026.GIF" wi="176" he="63" />C3、利用图像I<sub>re</sub>的二阶导数计算Hessian矩阵的特征值;C4、若Hessian矩阵的两个特征值|λ<sub>1</sub>>|λ<sub>2</sub>|,利用血管测度函数或血管特征拟合函数计算血管特征图<img file="FDA0000159059030000027.GIF" wi="805" he="223" />C5、得到血管特征图<img file="FDA0000159059030000028.GIF" wi="484" he="100" />C6、输出血管结构增强结果图<img file="FDA0000159059030000029.GIF" wi="242" he="76" /> |