X射线ct转台单侧两次螺旋扫描单层重排重建方法
【专利摘要】本发明涉及一种X射线CT转台单侧两次螺旋扫描单层重排重建方法，首先在转台的同一侧进行两次螺旋扫描，将两组螺旋锥形束投影重排成多层平行束投影，然后利用平行束投影的对称性质去除数据的横向截断，得到一组完全覆盖物体横截面的平行束投影数据，再通过滤波反投影方式进行图像重建。本发明通过在转台的同一侧进行两次螺旋扫描成像，能够更大程度上扩展螺旋锥束CT的横向成像视野，通过实验验证，能够扩展至2.56倍的成像视野，在更大程度上扩展成像视野的同时成像质量没有明显下降，具有与传统全覆盖算法相当的成像质量。
【专利说明】
X射线CT转台单侧两次螺旋扫描单层重排重建方法
技术领域
[0001]本发明属于CT扫描成像技术领域，特别涉及一种X射线CT转台单侧两次螺旋扫描 单层重排重建方法。
【背景技术】
[0002] X射线计算机断层成像(Computed Tomography，CT)是指利用X射线穿透物体的投 影数据进行重建来获得物体内部结构信息的技术，其成像过程涉及物理学、数学、计算机 学、图形图像学和机械学等多学科技术领域。随着面阵探测器在采集效率、耐辐射和抗干扰 等技术上的不断进步，以锥形束射线源和面阵探测器为基础的锥束CT在诸多领域获得了越 来越广泛的应用。与CT硬件发展对应的是扫描方式的不断进步。基于面阵探测器和第四代 扫描方式的螺旋锥束CT具有成像时间短、射线利用率高和分辨率各向同性的优势。此外，与 圆轨迹锥束CT相比，螺旋锥束CT能够获得更加充分的投影数据，解决大锥角时重建图像质 量退化的问题，并且能够解决轴向截断投影数据的重建问题，更加适合于对长物体的快速、 连续成像。目前，螺旋锥束CT已经越来越广泛的应用于医学临床诊断和工业无损检测、逆向 工程和材料组织分析等领域。
[0003] 由于成像系统中成像视野与成像分辨率通常相互制约，使得常规CT系统中很难同 时实现大视野和高分辨成像。螺旋锥束CT能够扩展锥束CT轴向的成像视野，但是受面阵探 测器等硬件尺寸的限制，其横向成像视野仍然是受限的。目前，基于视野扩展的成像方法在 很多方面得到了应用，比如对发动机、大型涡轮叶片以及印刷电路板的三维成像中。对于螺 旋锥束CT的横向视野扩展方法，目前的研究较少，主要有:半覆盖螺旋扫描和多次螺旋扫 描。半覆盖螺旋锥束CT成像能够扩展接近2倍的横向视野，但是，随着螺距的增大，这种扫描 方式缺失的投影数据越来越多，会使重建图像质量严重下降，由于半覆盖扫描几何下每个 投影角度的数据存在严重的横向截断，而该算法采用全局的斜坡滤波器，重建结果存在较 为严重的截断伪影。针对半覆盖螺旋roK算法存在的问题，Guo等人提出了改进的半覆盖螺 旋FDK算法，该算法借鉴2维CT中局部滤波的思想，重建过程中采用一种具有局部特性的滤 波器，新滤波器能够在一定程度上抑制由于数据截断和全局滤波带来的截断伪影;进一步 出现的工业双螺旋锥束CT，其扫描方式为:螺旋扫描前先控制载物台沿探测器行方向平移 一段距离，使射束能够覆盖物体横截面的一半，进行第一次螺旋轨迹扫描;然后控制载物台 沿相反的方向平移，使射束能够覆盖物体横截面的另一半，进行第二次螺旋扫描，利用两次 螺旋扫描得到的投影进行图像重建，但是，由于只采集了相对于成像视野两倍的投影数据， 因此其横向视野扩展率的极限只能达到2倍，这在很多情况下仍然不能满足实际应用的需 求。
【发明内容】
[0004] 为克服现有技术中的不足，本发明提供一种X射线CT转台单侧两次螺旋扫描单层 重排重建方法，通过在转台的同一侧进行两次螺旋扫描，将两组螺旋锥形束投影重排成多 层平行束投影，利用平行束投影的对称性质去除数据的横向截断，得到一组完全覆盖物体 横截面的平行束投影数据，通过滤波反投影方式进行图像重建，能够更加有效的扩展锥束 CT的成像视野。
[0005] 按照本发明所提供的设计方案，一种X射线CT转台单侧两次螺旋扫描单层重排重 建方法，包含如下步骤：
[0006] 步骤1、对被扫描物体在转台的同一侧进行两次螺旋扫描，两次螺旋扫描射线对被 扫描物体横截面的覆盖区域之间相互有交叠，两次螺旋扫描总的覆盖区域大于等于被扫描 物体横截面的一半；
[0007] 步骤2、通过螺旋扫描得到相应的螺旋锥形束投影，在成像几何下将螺旋锥形束投 影分别重排成相应的多层扇形束投影；
[0008] 步骤3、将多层扇形束投影分别重排成相应的多层平行束投影，并在重排过程中对 平行束几何进行统一，同时利用平行束投影的对称性质去除数据的横向截断，得到一组完 全覆盖被扫描物体横截面的平行束投影数据；
[0009] 步骤4、利用滤波反投影方法完成图像重建。
[0010] 上述的，步骤2中具体包含如下内容：
[0011] 步骤2.1、构建成像几何的等价几何模型，承载被扫描物体的转台做顺时针螺旋轨 迹运动等价于光源和探测器同步做逆时针螺旋轨迹运动，转台在水平方向的偏移等同于光 源和探测器同步沿相反方向的平移，以螺旋扫描的面阵探测器中心为原点定义旋转的笛卡 尔直角坐标系(1^^%)，坐标轴的方向向量为匕,(/:)、&(幻、〔'(乂），其中，1等于1或2，表示 第:次螺旋扫描，w轴与面阵探测器的法向量平行，u轴和v轴分别沿面阵探测器的行方向和 列方向；在固定坐标系(x，y，z)情况下，坐标轴方向向量表示为：
[0012]
[0013]
[0014] (6-,(/.) = (0,0,1);
[0015] 步骤2.2、根据等价几何模型，沿探测器水平方向同步偏移光源和探测器进行螺旋 扫描，得到螺旋锥形束投为投影角度，（u， v)表示面探测器上的点；
[0016] 步骤2.3、根据重排公式：
，得到相应的重排后的 多层扇形束投影表示扇形投影角度，u表示线阵探测器探元位置，z表示被扫描物 体第ζ层横截面，对于任意ζ层，p、（0，u)包含360角度范围内的投影数据，其中，
.和Θ满足θ = λ%231，％表示求余运算，s表示光源到探测器的距离，Δ Z表示 锥形束源点与虚拟扇形束源点的距离，D表示螺旋扫描中心射束位置到转台中心轴的距离。 [0017]上述的，所述步骤3具体包含如下内容：
[0018]步骤3.1、相应的多层平行束投影^ (β，t)重排公式表示为：
[0019]
[0020] ，β表示平行束投影角度，且^^[?。…，〖表示旋转轴到平行束的距离^:表示第一 次螺旋扫描的旋转半径，R 2表示第二次螺旋扫描的旋转半径，L表示第一次螺旋扫描旋转轴 到第二次螺旋扫描旋转轴的距离，tlmin、t lmax分别表示旋转轴到中射束的最小、最 大距离;t2min、t2max分别表示旋转轴到Ρ 2ζ(λ，11)中射束最小、最大距离；
[0021] 步骤3.2、根据平行束投影的对称性质82(巾3)=82(31+(})，-〇，其中，0和(})的关系 表示为：
，完全覆盖被扫描物体横截面的平行束投影数据重排公式表 \、ψ" ' 八·> 卜* |_ν *· ? ^-y ·· f 示为U
$Ζ(Φ，t)包含180角度范围内完全覆盖被扫描物体横截 面的投影数据，其中，Φ表示平行束的投影角度，且满足φε[0，π)及te[-RQ，RQ]，R Q表示为 转台横截面的半径。
[0022] 上述的，步骤4中利用滤波反投影方法完成图像重建包含:根据平行束投影数据gz (Φ，t)，采用平行束FBP方法完成图像重建，公式如下
，其 中，f(x，y，z)表示三维图像重建结果，ΚΖ(ΦΛ)表示斜坡滤波后的数据，/(扒X,)·)表示点(X， y，ζ)最终投影到平行束几何中位置。
[0023]本发明的有益效果：
[0024] 本发明解决现有技术中受面阵探测器尺寸等硬件条件的限制，锥束CT成像视野有 限，难以满足应用中对大尺寸物体成像的需求的问题，通过在转台的同一侧进行两次螺旋 扫描，将两组螺旋锥形束投影重排成多层平行束投影，利用平行束投影的对称性质去除数 据的横向截断，得到一组完全覆盖物体横截面的平行束投影数据，通过滤波反投影方式进 行图像重建;相比现有成像视野扩展方法，该扫描方式通过在转台的同一侧进行两次螺旋 扫描成像，能够更大程度上扩展螺旋锥束CT的横向成像视野，通过实验验证，能够扩展至 2.56倍的成像视野，在更大程度上扩展成像视野的同时成像质量没有明显下降，具有与传 统全覆盖算法相当的成像质量。
【附图说明】：
[0025] 图1为本发明的流程示意图；
[0026] 图2为螺旋锥束CT转台单侧两次扫描成像几何示意图；
[0027] 图3为螺旋锥束CT转台单侧两次扫描等价成像几何示意图；
[0028] 图4为多层扇形束投影成像几何示意图；
[0029] 图5为扇形束投影到平行束投影重排几何示意图；
[0030] 图6为仿真实验参数示意图；
[0031 ]图7为转台单侧两侧螺旋扫描和标准螺旋扫描第360度投影示意图；
[0032]图8为利用本发明和螺旋FDK方法进行图像重建的结果对比图；
[0033] 图9为图8所示重建结果中的部分垂直剖线示意图；
[0034] 图10为重建图像误差均方根对比示意图。
【具体实施方式】：
[0035] 下面结合附图和技术方案对本发明作进一步详细的说明，并通过优选的实施例详 细说明本发明的实施方式，但本发明的实施方式并不限于此。
[0036] 实施例一，参见图1所示，一种X射线CT转台单侧两次螺旋扫描单层重排重建方法， 包含如下步骤：
[0037] 步骤1、对被扫描物体在转台的同一侧进行两次螺旋扫描，两次螺旋扫描射线对被 扫描物体横截面的覆盖区域之间相互有交叠，两次螺旋扫描总的覆盖区域大于等于被扫描 物体横截面的一半；
[0038] 步骤2、通过螺旋扫描得到相应的螺旋锥形束投影，在成像几何下将螺旋锥形束投 影分别重排成相应的多层扇形束投影；
[0039] 步骤3、将多层扇形束投影分别重排成相应的多层平行束投影，并在重排过程中对 平行束几何进行统一，同时利用平行束投影的对称性质去除数据的横向截断，得到一组完 全覆盖被扫描物体横截面的平行束投影数据；
[0040] 步骤4、利用滤波反投影方法完成图像重建。
[0041] 本发明通过在转台的同一侧进行两次螺旋扫描，将两组螺旋锥形束投影重排成多 层平行束投影，利用平行束投影的对称性质去除数据的横向截断，得到一组完全覆盖物体 横截面的平行束投影数据，通过滤波反投影方式进行图像重建;相比现有成像视野扩展方 法，该扫描方式通过在转台的同一侧进行两次螺旋扫描成像，能够更大程度上扩展螺旋锥 束CT的横向成像视野。
[0042] 实施例二，参见图1~10所示，一种X射线CT转台单侧两次螺旋扫描单层重排重建 方法，包含如下内容：
[0043] 对被扫描物体在转台的同一侧进行两次螺旋扫描，两次螺旋扫描射线对被扫描物 体横截面的覆盖区域之间相互有交叠，两次螺旋扫描总的覆盖区域大于等于被扫描物体横 截面的一半;通过螺旋扫描得到相应的螺旋锥形束投影，在成像几何下将螺旋锥形束投影 分别重排成相应的多层扇形束投影;将多层扇形束投影分别重排成相应的多层平行束投 影，并在重排过程中对平行束几何进行统一，同时利用平行束投影的对称性质去除数据的 横向截断，得到一组完全覆盖被扫描物体横截面的平行束投影数据;利用滤波反投影方法 完成图像重建。具体内容如下：
[0044]螺旋锥束CT的转台位于光源和面阵探测器之间，扫描时，光源和探测器固定，转台 承载着被扫描物体做螺旋轨迹运动，具体成像几何参见图2所示，2(a)为俯视图、2(b)为侧 视图，两次螺旋扫描的轨迹分别对应2(a)中的虚线圆和实线圆，以及2(b)中的虚线螺旋线 和实线螺旋线，其旋转轴分别位于对应的01和02，两次螺旋扫描转台旋转轴所处的水平位置 不同，且都位于光源到探测垂线的一侧。转台做顺时针螺旋轨迹运动等价于光源和探测器 同步做逆时针螺旋轨迹运动;转台在水平方向上的偏移等同于光源和探测器同步沿相反方 向的平移，等价几何模型如图3所示，3(a)表示俯视图、3(b)表示侧视图，以转台中心〇2为原 点建立固定的直角坐标系(x，y，z)，定义转台横截面的半径为R〇,而物体的最大支撑为半径 为R〇的圆柱体，在图3(a)中物体的支撑与转台重合;定义光源垂直探测器的射束为中心射 束，等价几何模型中两次螺旋扫描是通过沿探测器水平方向同步偏移光源和探测器实现 的，其中，第一次中心射束位于〇1处，实现对物体外侧的扫描;第二次中心射束位于转台旋 转轴02处，实现对物体中部的扫描，两次扫描都以转台中心轴02为旋转轴，光源到探测器的 距离表示为S，图中所示， 81到〇1的距离和82到〇2的距离相同且等于0, 〇1到〇2的距离为L，则两 次螺旋扫描的旋转半径分别为/?1=>/^Τ? Τ和R2 = D，如图3(b)所示，两次扫描的起始高度、 螺距是相同的。
[0045] 在固定坐标系(x，y，z)下两次扫描光源的轨迹分别表示为：
[0046]
[0047]
[0048]，其中，λ表示旋转角度，h表示螺距，λ〇为无符号角度且A〇 = arctan(L/D)，以螺旋扫 描的面阵探测器中心为原点定义旋转的笛卡尔直角坐标系(m，Vl，Wl)，坐标轴的方向向量 为匕μ)、& μ)、(1),其中，w轴与面阵探测器的法向量平行，u轴和V轴分别沿面阵探测器 的行方向和列方向；在固定坐标系(X，y，Ζ)情况下，坐标轴方向向量表示为：
[0049]
[0050]
[0051 ] ^ α) = (〇,ο,?)/
[0052] 图3(a)中，定义转台中心到射束方向的距离t，从光源向转台中心，如果射线位于 转台的左侧则t为负的，反之，则t为正的；定义转台中心到两次螺旋扫描中射束的最小和最 足巨 1? '力U 、tlmax矛口 t2min、t2max，?两iS·:
[0053]
[0054] 根据等价几何模型，沿探测器水平方向同步偏移光源和探测器进行螺旋扫描，得 到螺旋锥形束投影为投影角度，（u， v)表示面探测器上的点，两组投影数据的 重排相互独立;根据重排公式：
，得到相应的重排后的多层 扇形束投影表示扇形投影角度，u表示线阵探测器探元位置，z表示被扫描物体第 ζ层横截面;多层扇形束投影各层之间的数据相互独立，从扇形束投影到平行束投影的重排 只涉及到相同ζ层的数据，对于任意ζ层，包含360角度范围内的投影数据，其中，
λ和Θ满足θ = λ%231，％表示求余运算，s表示光源到探测器的距离，Δ Ζ表示 锥形束源点与虚拟扇形束源点的距离，D表示螺旋扫描中心射束位置到转台中心轴的距离。
[0055] 对于任意z层，定义旋转轴到(λ, u)和ρ2ζ (λ, u)中射束的最小和最大距离分别为 匕^、^_和^*山_，重排后9办，〇内所有到旋转轴的距离小于^^的射线全部来源于 ?^，11)，而所有到旋转轴的距离大于等于^?1"的射线全部来源于?22〇， 11)，通过这种方式， 便可以将两组独立的扇形束投影融合成一组完整的平行束投影，并且能够去除两组投影之 间的冗余数据。具体的，重排过程中的几何关系如图5所示，其中，光源到探测器的距离为S， 81到〇1的距离和82到〇2的距离相同且等于D，而 〇1到〇2的距离为L，两次螺旋扫描的旋转半径 分别为
PR2 = D，由图5所示的成像几何可知，
γ表示带符号扇角，Μ表示转台中心到〇1(〇2)的距离。对于第一组扇形束投影，M = L，
:而对于第二组扇形束投影
由于，当t< t2min时射线全部来源于(λ，U);当t彡t2min时射线全部来源于Ρ2ζ(λ，U)，相应的多层平行束 投影qz(i3，t)重排公式表示为：
[0056]
[0057] ，β表示平行束投影角度，且^^[?。…，丨表示旋转轴到平行束的距离^:表示第一 次螺旋扫描的旋转半径，R 2表示第二次螺旋扫描的旋转半径，L表示第一次螺旋扫描旋转轴 到第二次螺旋扫描旋转轴的距离;根据平行束投影的对称性质&(Φ，t)=g z(3I+(})，-t)，其 中，β和Φ的关系表示为
完全覆盖被扫描物体横截面的平行束投影 数据重排公式表示为：
[0058]
gz( Φ，t)包含180角度范围内完全覆盖被扫描物体横 截面的投影数据，其中，Φ表示平行束的投影角度，且满足Φε[0，π)及te[-RQ，RQ]，R Q表示 为转台横截面的半径。
[0059] 利用滤波反投影方法完成图像重建，根据平行束投影数据gz((i)，t)，采用平行束 FBP方法完成图像重建，公式如下：
其中，f(x，y，z)表示三维 图像重建结果，ΚΖ(Φ，t)表示斜坡滤波后的数据，?(扒U)表示点（x，y，z)最终投影到平行 束几何中位置;为减少重排误差在反投影重建中对重建图像质量的影响，反投影重建方法 也可以通过最初的成像几何计算投影位置来完成图像重建。
[0060] 下面通过具体的数字仿真实验进一步验证本发明的有效性，并将实验结果与经典 螺旋FDK方法进行比较：
[0061 ]构建如图3所示的螺旋锥束CT转台单侧两次螺旋扫描成像几何，并对三维Shepp- Logan体模进行数据采集;同时，在标准螺旋成像几何下对Shepp-Logan体模进行数据采集。 两次采得的投影数据分别使用本发明和经典螺旋FDK方法进行重建，具体的成像几何参数 如图6所示，转台单侧使用的探测器的尺寸为200X100，而Shepp-Logan体模尺寸为256 X 256X256,即物体在横向上超出成像视野，投影数据存在数据截断;标准螺旋扫描仿真实验 中探测器尺寸为200 X 256，而物体的尺寸为256 X 256 X 256，即投影数据不存在横向数据截 断。转台单侧两次螺旋扫描的第360度投影数据分别如图7中7(a)和7(b)所示，图7(a)和7 (b) 分别表示第一次螺旋扫描和第二次螺旋扫描;而标准螺旋扫描第360度投影数据如图7 (c) 所示;本发明和经典螺旋FDK方法的重建结果如图8所示，图8中：（&)、（13)、（(：)分别表示 经典螺旋FDK方法重建结果的x-y、y_z和x-z中心切片，（d)、（e)、（f)分别表示本发明重建结 果的x-y、y-z和χ-ζ中心切片;重建结果剖线如图9所示，图9中（a)、（b)、（c)分别表示图8中 所示的x-y、y-z和χ-ζ中心切片的垂直剖线。
[0062]为进一步对重建结果进行评价，通过对重建三维图像的误差均方根进行比较，如 图10，误差均方根可根据式
·算得到，其中，f jPfr分别表示 真值和重建结果，N表示体素个数。
[0063] 由图8的重建结果可以看出，本发明能够有效对螺旋锥束CT转台单侧两次扫描的 投影数据实施重建，重建结果中没有明显的截断伪影和图像不均匀现象;进一步从图9中的 剖线图数值比较可以看出，本发明重建结果在数值上与真值相差很小，表明本发明重建结 果具有较高的数值准确性；图9和图10中与经典螺旋FDK方法使用标准螺旋锥束投影重建结 果的对比可以说明，本发明重建结果没有引入新的误差，重建图像质量与标准螺旋Π )Κ算法 的重建质量相当。
[0064] 以上说明本发明能够有效的处理投影数据的融合和数据截断问题，在扩大CT成像 视野的同时，成像不会引入额外误差，重建质量较好，与标准螺旋roK算法重建质量相当；同 时，与目前螺旋锥束CT扩展视野算法的横向视野扩展率不能大于2倍的现状相比，本发明能 够更有效地扩展螺旋锥束CT的横向成像视野，如图6所示，本次实验的横向视野扩展率为 2.56〇
[0065] 本发明不局限于上述【具体实施方式】，本领域技术人员还可据此做出多种变化，但 任何与本发明等同或者类似的变化都应涵盖在本发明权利要求的范围内。
【主权项】
1. 一种X射线CT转台单侧两次螺旋扫描单层重排重建方法，其特征在于：包含如下步 骤： 步骤1、对被扫描物体在转台的同一侧进行两次螺旋扫描，两次螺旋扫描射线对被扫描 物体横截面的覆盖区域之间相互有交叠，两次螺旋扫描总的覆盖区域大于等于被扫描物体 横截面的一半； 步骤2、通过螺旋扫描得到相应的螺旋锥形束投影，在成像几何下将螺旋锥形束投影分 别重排成相应的多层扇形束投影； 步骤3、将多层扇形束投影分别重排成相应的多层平行束投影，并在重排过程中对平行 束几何进行统一，同时利用平行束投影的对称性质去除数据的横向截断，得到一组完全覆 盖被扫描物体横截面的平行束投影数据； 步骤4、利用滤波反投影方法完成图像重建。2. 根据权利要求1所述的X射线CT转台单侧两次螺旋扫描单层重排重建方法，其特征在 于:步骤2中具体包含如下内容： 步骤2.1、构建成像几何的等价几何模型，承载被扫描物体的转台做顺时针螺旋轨迹运 动等价于光源和探测器同步做逆时针螺旋轨迹运动，转台在水平方向的偏移等同于光源和 探测器同步沿相反方向的平移，W螺旋扫描的面阵探测器中屯、为原点定义旋转的笛卡尔直 角坐标系(心￥：1，'\￥〇，坐标轴的方向向量为^，，（2)、^^桃、e?,.(义).，其中，i等于1或2,表示第i 次螺旋扫描，W轴与面阵探测器的法向量平行，U轴和V轴分别沿面阵探测器的行方向和列方 向；在固定坐标系(x，y，z)情况下，坐标轴方向向量表示为：步骤2.2、根据等价几何模型，沿探测器水平方向同步偏移光源和探测器进行螺旋扫 描，得到螺旋锥形束投影6ι(λ，υ，ν)，λ为投影角度，（u，v)表示面探测器上的点； 步骤2.3、根据重排公式，得到相应的重排后的多层 扇形束投影ρ?ζ(θ，ιι)，θ表示扇形投影角度，U表示线阵探测器探元位置，Z表示被扫描物体第 Ζ层横截面，对于任意Ζ层，piz(0，u)包含360角度范围内的投影数据，其中和Θ满足θ = λ%2π，％表示求余运算，S表示光源到探测器的距离，Δζ表示锥形束源点与虚 拟扇形束源点的距离，D表示螺旋扫描中屯、射束位置到转台中屯、轴的距离。3. 根据权利要求2所述的X射线CT转台单侧两次螺旋扫描单层重排重建方法，其特征在 于:所述步骤3具体包含如下内容： 步骤3.1、相应的多层平行束投影qz (β，t)重排公式表示为：e表示平行束投影角度，且0￡[〇,231)，*表示旋转轴到平行束的距离，町表示第一次螺旋 扫描的旋转半径，R2表示第二次螺旋扫描的旋转半径，L表示第一次螺旋扫描旋转轴到第二 次螺旋扫描旋转轴的距离，分别表示旋转轴到ρ1ζ(λ，11)中射束的最小、最大距离； t2min、t2max分别表示旋转轴到Ρ2ζ(λ，11)中射束最小、最大距离； 步骤3.2、根据平行束投影的对称性质gz(Φ，t)=gz(3?+Φ，-t)，其中，β和Φ的关系表示 为：完全覆盖被扫描物体横截面的平行束投影数据重排公式表示 为：gz( Φ，t)包含180角度范围内完全覆盖被扫描物体横截面 的投影数据，其中，Φ表示平行束的投影角度，且满足φε[〇,3?)及化[-柄，帖]，柏表示为转 台横截面的半径。4.根据权利要求3所述的X射线CT转台单侧两次螺旋扫描单层重排重建方法，其特征在 于:步骤4中利用滤波反投影方法完成图像重建包含:根据平行束投影数据gz( Φ，t)，采用 平行束FBP方法完成图像重建，公式如下：'其中，f(x，y，z)表示Ξ维图像重建结果，Κζ(Φ，t)表示 斜坡滤波后的数据，攻巧X，的表示点(X，y，Z)最终投影到平行束几何中位置。
【文档编号】A61B6/03GK105973917SQ201610512467
【公开日】2016年9月28日
【申请日】2016年6月29日
【发明人】韩玉, 李磊, 程根阳, 闫镔, 席晓琦, 王林元, 王彪
【申请人】中国人民解放军信息工程大学