一种采用圆轨道扇形束x射线ct扫描机快速重建断层图像的方法
【专利摘要】本发明涉及一种采用圆轨道扇形束X射线CT扫描机快速重建断层图像的方法，该方法包括以下步骤：(1)获取经对数化处理的投影数据，并定义为数据集A，然后对数据集A进行一维滤波，把将所得到的滤波后投影数据定义为数据集B；(2)构建一幅离散的待重建图像，并将其中所有像素点的像素值初始化为0；(3)将待重建图像中每个像素点进行反投影，即，首先搜索数据集A，寻找该像素点对应的投影值，一旦投影值为零便直接将该像素点的像素值赋为零，然后寻找下一像素点的投影值，而当该像素点的投影值不为零时，则搜索数据集B，寻找该像素点的投影值；如此历遍全部待重建图像，分别将未赋值的各像素点对应于数据集B中的投影值进行累加并赋值，得到重建后的断层图像。
【专利说明】一种采用圆轨道扇形束X射线CT扫描机快速重建断层图像 的方法
【技术领域】
[0001] 本发明涉及一般的图像数据处理，具体涉及2D图像的生成方法，该方法适用于扇 形束X射线CT图像的滤波反投影重建。
【背景技术】
[0002] X射线CT (computed tomography,计算机断层摄影术）广泛应用工业无损检测和 医学影像诊断和治疗等领域。CT图像质量好坏取决于所使用的重建算法。其中，扇束X射 线CT滤波反投影重建算法属于解析算法，相比于迭代算法，它具有计算量明显小的优势， 同时，CT图像质量和空间分辨率也较高，是一种非常实用的成像算法，因此，滤波反投影重 建算法仍是目前商用X线CT系统所采用的主流算法。
[0003] 申请公布号为CN103489206A专利申请公开了"一种基于扇束X光CT滤波反投影 重建的混合滤波方法"，该方法采用以下方法进行滤波反投影重建：（1)获取扇束X光CT投 影数据；（2)对投影数据进行混合滤波处理；（3)沿着X射线方法进行反投影重建得到断层 图像。
[0004] 此外，上述专利申请还公开了对扇束X光CT投影数据离散化的方法以及对离散化 后的投影数据进行混合滤波处理方法和混合滤波所使用的滤波器，其中，
[0005] 所述对扇束X光CT投影数据离散化的方法即是，设置采样点S = nd以间隔d = 1，角度β = 1?2π对离散化前的投影数据进行采样，得到离散后的投影数据p(s，β);
[0006] 所述对离散化后的投影数据进行混合滤波处理方法为：设定h(S)为混合滤波器， 根据P(S，β) =p(s，i3)**h(S)对投影数据p(s，β)进行卷积混合滤波，得到滤波后的投 影数据冲，/
[0007] 所述混合滤波所使用的滤波器为：
[0008] h (nd) = λ ^sl (nd) + λ 2hL (nd)；
[0009] 加权系数：λ 1彡0, λ 2彡0, λ彳λ 2 = 1，用户根据实际需要确定；
[0010] SL滤波器：
【权利要求】
1. 一种采用圆轨道扇形束X射线CT扫描机快速重建断层图像的方法，该方法包括以下 步骤： (1) 获取采用圆轨道扇形束X射线CT扫描机在每一预设角度下经对数化处理的投影数 据，并定义为数据集A，然后对数据集A中每一预设角度下经对数化处理后的投影数据进行 一维滤波，把将所得到的滤波后投影数据定义为数据集B ; (2) 构建一幅离散的待重建图像，并将其中所有像素点的像素值初始化为0 ; (3) 将待重建图像中每个像素点沿X射线方向在0?360°范围内依次按所述预设角 度进行反投影，该反投影步骤为：首先搜索数据集A，寻找该像素点在当前角度下所对应的 投影值，一旦投影值为零便直接将该像素点在当前角度下的像素值赋为零，然后寻找下一 像素点在当前角度下所对应的投影值，而当该像素点在当前角度下所对应的投影值不为零 时，则搜索数据集B，寻找该像素点在当前角度下所对应的投影值；如此历遍全部待重建图 像，分别将未赋值的各像素点在每一角度下对应于数据集B中的投影值进行累加，并将所 得到的累加值作为对应像素点的像素值，得到重建后的断层图像。
2. 根据权利要求1所述的一种采用圆轨道扇形束X射线CT扫描机快速重建断层图像 的方法，其特征在于，所述一维滤波的滤波器为SL滤波器、RL滤波器或Lewitt滤波器。
【文档编号】G06T11/00GK104361615SQ201410568553
【公开日】2015年2月18日 申请日期:2014年10月22日 优先权日:2014年10月22日
【发明者】周凌宏, 齐宏亮, 卢文婷, 李翰威, 李斌, 郭静钰, 陈梓嘉 申请人:南方医科大学