一种数字减影血管造影运动伪影消除方法及其系统

摘要

一种数字减影血管造影运动伪影消除方法，步骤如下：步骤1、读取；步骤2、选点；步骤3、构建DSA时空体；步骤4、时空切片；步骤5、连轨迹；步骤6，分析DSA像素的时空运动特性；步骤7、三角剖分；步骤8、仿射变换；步骤9、时空卷绕；步骤10、优化；步骤11、配准；步骤12、灰度校正；步骤13，对数减影。本发明属于图像处理技术。采用时空分析的方法，使DSA图像配准更加精确，有效去除运动伪影，得到清晰的血管造影图像。进而提高医生的诊断准确性，提高工作效率。

主权项

一种数字减影血管造影运动伪影消除方法，其特征在于步骤如下：步骤1、读取；将数字减影血管造影仪器与计算机连接，在计算机中安装操作系统和图像处理系统，用图像处理软件读入一组DSA图像，并将DSA图像进行预处理；步骤2、选点；采用边缘检测算法找到每幅DSA图像中血管影像的边缘和边缘点，计算出边缘点的梯度值，根据阈值来判断该点是否作为备选控制点，并选取至少10个备选控制点；步骤3、构建DSA时空体；将该组DSA二维图像以时间T为轴按拍摄顺序间隔排列成平行六面体空间；步骤4、时空切片；用切片指令对DSA时空体做垂直方向的时空切片；步骤5、连轨迹；用图像处理软件将各时空切片对应的各控制点连线，得到各个控制点的运动轨迹；步骤6，分析DSA像素的时空运动特性；选择运动轨迹变化最大的10个反映DSA图像背景结构运动规律的控制点；步骤7、三角剖分；用三角剖分的方法将所有的控制点作为三角形的顶点，把图像剖分成不均匀的三角形网格；步骤8、仿射变换；采用仿射变换法进行初步几何校正；步骤9、时空卷绕；用Warping非线性变换法参考前一幅图像中控制点的变形位置，逐幅确定图像非控制点的位移；步骤10、优化；采用Powell优化算法来优化相似度准则函数求极值的过程；步骤11、配准；选择三次样条插值算法进行亚像素级的配准；步骤12、灰度校正；逐幅对图像的灰度失真进行回朔校正；步骤13，对数减影；用对数减影法进行减影，获得消除伪影的减影图像。