[发明专利]一种数字减影血管造影运动伪影消除方法及其系统

权利要求书

1.一种数字减影血管造影运动伪影消除方法，其特征在于步骤如下：

步骤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、对数减影：用对数减影法进行减影，获得消除伪影的减影图像。

2.根据权利要求1所述的一种数字减影血管造影运动伪影消除方法，其特征在于：所述边缘检测算法采用Canny算子或高斯-拉普拉斯算子。