【技术实现步骤摘要】
本专利技术属于计算机
，涉及ー种旋转X射线造影图像迭代重建方法。
技术介绍
旋转X射线冠状动脉造影成像技术是继双平面血管造影成像以来又ー项得到广泛关注的冠状动脉成像技木。通过该技木，医师能够从多个角度提供完整，准确地观察血管的形状以及运动方式，对脑血管肿瘤，冠心病的诊疗具有重要意义。但仅凭旋转图像序列，难以对冠状动脉的具体结构有一个直观准确的认识，因此血管造影的三维重建成为目前学术领域的研究热点，也是医疗器械厂商亟待解决的重要问题。旋转X射线造影技术的成像几何同三维重建同锥束CT类似。由于旋转造影成像过程中伴随着心脏运动，为了获得准确的血管三维结构，需要采用迭代重建算法。但是由于该成像系统的空间分辨率高于一般的 锥束CT，该重建问题的计算量远远大于锥束CT图象重建。迭代算法的运算主要集中在投影、反投影两个环节，投影、反投影的快速计算问题成为目前旋转X射线血管造影成像三维重建问题的主要瓶颈，投影矩阵的使用可以大大加速投影，反投影操作的运算速度，但受到硬件条件限制，目前的计算机内存存储能力不足以存储旋转造影系统投影矩阵。
技术实现思路
技术问题本专利技术提供了一种...
【技术保护点】
一种旋转X射线造影图像迭代重建方法，其特征在于，包括以下步骤：1）从旋转X射线造影设备读取扫描数据文件，保存投影序列图像并记录如下参数：射线源到探测板距离SDD，射线源到C臂旋转中心距离SOD，各投影采样的旋转角度，投影图像像素边长h，投影图像长度U个像素单位和宽度V个像素单位，所述投影图像尺寸与二维投影空间尺寸一致；确定二维投影空间坐标轴u、v，所述坐标轴u、v分别平行于二维投影空间的长、宽方向，根据用户的精度需求，分别设定三维图像空间长度为X个体素单位、宽度为Y个体素单位、高度为Z个体素单位，三维图像空间体素边长l，确定三维图像空间坐标轴x、y、z，所述坐标轴x、y、z...
【技术特征摘要】
1.一种旋转X射线造影图像迭代重建方法，其特征在于，包括以下步骤 1)从旋转X射线造影设备读取扫描数据文件，保存投影序列图像并记录如下参数射线源到探测板距离SDD，射线源到C臂旋转中心距离SOD，各投影采样的旋转角度，投影图像像素边长h，投影图像长度U个像素单位和宽度V个像素单位，所述投影图像尺寸与二维投影空间尺寸一致；确定二维投影空间坐标轴U、V，所述坐标轴U、V分别平行于二维投影空间的长、宽方向，根据用户的精度需求，分别设定三维图像空间长度为X个体素单位、宽度为Y个体素单位、高度为Z个体素单位，三维图像空间体素边长1，确定三维图像空间坐标轴χ、y、z，所述坐标轴x、y、z分别平行于三维图像空间的长、宽、高方向，且均通过三维图像空间的中心位置； 2)根据步骤I)中设定的三维图像空间长宽高对三维图像空间进行降采样操作，具体方法为选定降采样倍数dif，所述降采样倍数dif为能被三维图像空间长度X、宽度Y、高度Z分 别整除的整数，将三维图像空间中由dif*dif*dif个体素组成的、边长为dif*l的立方体内的体素，归并为降采样体素，对整个三维图像空间完成归并处理后,按照归并顺序排列所述降采样体素，得到降采样三维图像空间，降采样体素边长为三维图像空间体素边长的dif倍； 根据步骤I)中记录的投影图像长宽，对二维投影空间进行降采样操作，具体方法为选定降采样倍数Clpf，所述降采样倍数为能被二维投影空间长度U和宽度V分别整除的整数，将二维投影空间中由dpf*dpf个像素组成的、边长为dpf*h的正方形内的像素，归并为投影空间降采样像素，对整个二维投影空间完成归并处理后，按照归并顺序排列所述投影空间降采样像素，得到降采样二维投影空间，投影空间降采样像素边长为二维投影空间像素边长的Clpf倍； 根据步骤I)中记录的投影图像长宽，对步骤I)获得的投影序列图像进行降采样操作，得到降采样序列投影图像，具体方法为对投影序列中每个投影图像，将其中由dpf*dpf个像素组成的、边长为dpf*h的正方形内的像素，归并为投影图像降采样像素，并且对所述正方形内的像素值进行累加求和，作为投影图像降采样像素值，对整个投影图像完成归并处理后，按照归并顺序排列所述投影图像降采样像素，得到降采样投影图像，投影图像降采样像素边长为投影图像像素边长的Clpf倍； 3)根据步骤I)中所记录的探测板距离SDD，射线源到C臂旋转中心距离SOD，初始采样旋转角度，针对步骤2)得到的降采样三维图像空间和降采样二维投影空间，利用距离驱动算法，构造初始扫描方向低分辨率体素索引投影矩阵，并将得到的初始扫描方向低分辨率体素索引投影矩阵以三元组存放方式保存于内存设备中； 4)针对步骤2)提供的降采样三维图像空间，构造各旋转角度下降采样三维图像空间旋转矩阵R，具体方法为首先根据图像旋转角度位置关系，利用线性插值算法构造并存储降采样图像空间切面层旋转矩阵R°，随后根据降采样三维图像空间旋转矩阵R的分块对称性，利用公式Λ =到降采样图像空间旋转矩阵R ； Z 5)对步骤2)得到的降采样投影序列图像采用顶帽滤波方法进行滤波处理，对滤波处理结果二值化得到第一分割结果，对步骤2)得到的降采样投影序列图像采用Frangi血管滤波方法进行滤波处理，对滤波处理结果二值化得到第二分割结果，对第一分割结果和第二分割结果求并集，得到低分辨率投影序列图像血管分割结果；6)利用步骤3)所得到的三元组格式存储的初始扫描方向低分辨率体素索引投影矩阵和步骤4)所得到的各旋转角度降采样图像空间旋转矩阵，对步骤5)所得到的低分辨率投影序列图像血管分割结果进行反投影操作具体方法为对低分辨率投影图像血管分割结果，根据像素编号排列拉伸低分辨率投影图像血管分割结果，得到低分辨率投影图像血管分割结果向量，将所得向量同初始扫描方向低分辨率体素索引投影矩阵的转置矩阵相乘得到中间低分辨率反投影结果；将该中间低分辨率反投影结果同对应旋转角度降采样图像空间旋转矩阵的转置矩阵相乘，得到该旋转角度下的低分辨率反投影结果； 7)利用步骤6)所得到的低分辨率反投影结果确定低分辨率三维图像空间血管掩模，具体方法为当投影数目小于等于5时，对各旋转角度下低分辨率反投影结果求交集，作为低分辨率三维图像空间血管掩模；当投影数目大于5时，利用对各旋转角度下低分辨率反投影结果求和并进行阈值划分，作为低分辨率三维图像空间血管掩模； 8)对步骤7)得到的低分辨率三维图像空间血管掩模，根据步骤2)所设定的降采样倍数dif进行升采样操作，得到三维图像空间血管掩模，具体方法为将低分辨率三维图像空 间血管掩模的每一个低分辨率体素在X，1，Z方向等分为dif份，拆分成为dif*dif*dif个高分辨率体素，将该低分辨率体素的值作为其拆分成的高分辨率体素的值，并按照拆...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。