具体实施方式

下面结合附图对本发明做进一步的详细说明，但并不是对本发明保护范围的限制。

在CT校直过程中，操作者通过机械调整球管焦点在x-轴的位置，使得iso通道的位置达到期望的理想值。此过程需要扫描一个偏心放置的截面为圆形的细长金属pin(直径约5mm的圆柱体），得到扫描投影数据，该数据由n个view的投影数据组成。根据投影数据，结合探测器的每一个通道具体的光束的方向的角度信息，通过迭代的方法，可以精确计算出iso通道的位置以及球管焦点在x-轴的位置坐标。

本发明通过仿真计算来比较iso通道的计算结果，如图5所示，针对一种按照特定几何布局的CT系统，按照传统方法，iso通道是330.8081；按照本发明的计算方法，iso通道是330.8297；实际的iso通道是330.8299。本发明的计算方法的结果与实际结果的误差为-0.0002个通道，本发明的计算方法的结果比传统方法更为精确。

如图1所示，CT系统的校直涉及到球管和探测器，此二者围绕旋转中心iso进行旋转，对偏心放置的金属pin进行扫描，收集到的扫描投影数据被送到计算机或者特定处理器进行处理和精确计算，得到实际的iso通道以及球管焦点的实际位置，通过与期望的iso通道进行对比进一步得到球管焦点位置的调整量，然后对球管焦点位置进行进行机械调整，继续扫描金属pin,如此反复，使得iso通道达到期望位置，最终实现对CT系统的精确校直。

本发明的特点在于，根据扫描得到的投影数据，结合探测器的每一个通道具体的光束的方向的角度信息，通过迭代的方法，精确计算出iso通道的位置以及球管焦点在x-轴的位置坐标。

如图6和图1所示，具体的扫描和计算过程为：

1）如图1所示，扫描细长金属pin，得到如图2所示的扫描的投影数据。

2）假定球管焦点位于期望位置。

3）根据CT的几何布局，计算对应的iso通道位置。

4）根据CT的几何布局，计算出探测器的每一个通道的入射光束的方向的角度信息，即该光束与y-轴的夹角，如图3所示。

5）利用扫描的投影数据，计算每一个view的所有通道的金属投影值相对于各个通道的光束角度的重心值，该值表示投影重心所在的光束方向的角度。

其中表示view i,通道j的投影值, 表示view i,通道j的方向的角度，m表示总的通道数，表示view i的投影重心所在的光束方向的角度。

6）平均所有view下的投影重心值，得到iso通道所对应的光束方向的角度。

7） 计算球管焦点的x-轴的新位置f。

其中D表示球管焦点到iso的距离。

8）重复3）到7），直到达到规定的迭代次数(例如30次)， 或者相邻两次计算得到的球管焦点的x-坐标的误差小于某一规定值（例如1.0e-4 mm）。 或者相邻两次计算得到的iso通道的值的误差小于某一规定值（例如1.0e-3 通道），如图4所示。

9）根据计算得到的球管焦点的x-坐标，计算焦点位置的调整量d，对球管焦点位置进行机械调整。

其中d为调整量，为期望的焦点位置，代表期望的iso通道的光束角度所对应的焦点位置。

10）重复1）到9），直到扫描和计算得到的iso通道的值与设计的期望值之间的误差小于某一规定值（例如1.0e-3通道）,完成CT的精确校直。

为了加速迭代计算过程，可以采样多种方法设定球管焦点的初始位置。如图1和图7所示，扫描金属pin后，可以根据传统方法计算iso通道以及对应的焦点位置，将其设定为初始的焦点位置。传统方法根据金属pin在每一个通道的投影值以及其所在的通道数来计算iso通道，这一计算过程忽略了具体通道方向的角度信息，其结果并不精确，但它也是对真实结果的一种逼近或者近似，因此可以用来作为焦点位置的初始值来加速迭代过程。

设定球管焦点的初始位置的第2种方法，如图1和图8所示，可以把机械调整前的球管焦点位置设定为初始焦点位置，来加速迭代计算过程。

设定球管焦点的初始位置的第3种方法，如图1和图9所示，可以把机械调整前的球管焦点位置与期望位置的平均值设定为初始焦点位置，来加速迭代计算过程。