CN102253061A - 立式锥束计算机断层成像校准***及应用该***的方法 - Google Patents

立式锥束计算机断层成像校准***及应用该***的方法 Download PDF

Info

Publication number
CN102253061A
CN102253061A CN2011100979583A CN201110097958A CN102253061A CN 102253061 A CN102253061 A CN 102253061A CN 2011100979583 A CN2011100979583 A CN 2011100979583A CN 201110097958 A CN201110097958 A CN 201110097958A CN 102253061 A CN102253061 A CN 102253061A
Authority
CN
China
Prior art keywords
turning axle
vertical
metal ball
rotation
detector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2011100979583A
Other languages
English (en)
Other versions
CN102253061B (zh
Inventor
罗守华
张奎
林晓珍
顾宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN 201110097958 priority Critical patent/CN102253061B/zh
Publication of CN102253061A publication Critical patent/CN102253061A/zh
Application granted granted Critical
Publication of CN102253061B publication Critical patent/CN102253061B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明涉及一种立式锥束CT校准***及应用该***的方法，该校准***包括X射线源（100）、与X射线源（100）相距一定距离的探测器（106）和设置在X射线源（100）与探测器（106）之间的载物台（102）、用于固定载物台（102）的旋转轴（315）、固定在载物台（102）上的金属球校准模型（104）。该方法包括以下步骤：测量旋转轴在旋转过程中的旋转中心；测量旋转轴的旋转一致性；测量旋转轴在垂直移动过程中的一致性；测量旋转轴在水平移动过程中的一致性。通过本发明使提高***的重建效率和重建质量。

Description

立式锥束计算机断层成像校准***及应用该***的方法
技术领域
本发明涉及一种立式锥束CT（Computer Tomography，计算机断层成像技术）***校准的方法，尤其是指射线源与探测器固定，通过使样本与载物台一起随旋转轴旋转、上下前后移动获取扫描数据的锥束CT***的校准方法，更具体的涉及了一种用FDK算法进行重建的锥束CT***校准的方法。FDK算法是Feldkamp，Davis，Kress在1984年提出的基于锥束投影的滤波反投影重建算法。
背景技术
锥束CT***具有空间分辨率高，采集时间短和射线利用率高等显著特点。锥束CT***中一般常用的重建算法是FDK算法。FDK算法是Feldkamp，Davis，Kress在1984年提出的基于锥束投影的滤波反投影重建算法，其特点是理论结构简单，重建速度快，易于实际应用，尤其是在锥角较小的情况下可获得较好的重建结构。然而，FDK算法要求应用该算法的***满足一定的成像几何结构：锥束CT***中，与探测器垂直且经过射线源的平面称为中平面，旋转轴与中平面垂直，同时与探测器像元阵列的列方向平行、行方向垂直；过锥束射线源焦点做垂直探测器的直线，该直线经过旋转轴；旋转轴的旋转保持稳定，在旋转过程中始终保持与中平面垂直；射线源到探测器的垂线的垂足是用于重建的投影图像的中心。
然而由于零部件机械加工精度有限以及安装误差，锥束CT***安装完成之后，***成像几何结构无法精确满足FDK算法的要求，即使经过算法几何校正修正，重建图像的质量也会有所下降。本发明提出的一种立式锥束CT***校准方法评估当前***的成像几何结构，并给出***机械微调的指导意见。锥束CT***经调整后在机械上达到一个比较精确地成像几何结构，减少重建算法中几何校正的负担。锥束CT***的成像几何结构越接近FDK算法要求的成像几何结构，***的重建效率越高，重建图像的质量也越高。所以，锥束CT***校准可以有效提高***的重建效率和重建图像质量。
发明内容
技术问题：本发明要解决的技术问题是提了供一种立式锥束CT***校准***及应用该***的立式锥束CT***校准方法，以提高***的重建效率和重建质量
技术方案：为解决上述技术问题，本发明提供了一种立式锥束CT***校准***，该***包括X射线源、与X射线源相距一定距离的探测器和设置在X射线源与探测器之间的载物台、用于固定载物台的旋转轴、固定在载物台上的金属球校准模型；其中金属球校准模型在扫描过程中与载物台一起随旋转轴旋转；X射线源与探测器的位置固定不动，通过载物台旋转获取投影数据；金属球校准模型上固定一金属球；金属球校准模型、载物台和旋转轴构成样本***。
优选的，金属球校准模型是有机玻璃制成的圆柱形杯子，杯子的上边缘固定金属球。
本发明还提供了一种立式锥束CT校准方法，该方法包括以下步骤：测量旋转轴在旋转过程中的旋转中心；测量旋转轴的旋转一致性；所述旋转一致性是指使样本***固定在某个位置，使金属球校准模型位于合适的旋转采样点，金属球校准模型旋转360°M次，M为自然数，分别在每次旋转后计算金属小球的位置，金属小球位置的一致性就是旋转轴的旋转一致性；测量旋转轴在垂直移动过程中的一致性；所述旋转轴在垂直移动过程中的一致性是指在旋转轴上下移动的过程中旋转中心的一致性；测量旋转轴在水平移动过程中的一致性，所述旋转轴在水平移动过程中的一致性是指在旋转轴前后移动的过程中旋转中心的一致性；
根据上述数实现立体锥束CT校准***的校准。
优选的，所述旋转轴的旋转中心的测量方法是：使样本***固定在某个位置，金属球校准模型随旋转轴旋转，均匀地取N个旋转采样点并取其在探测器上的投影图像；根据金属球校准模型旋转一周的投影数据计算旋转轴的旋转中心，N为偶数。
优选的，所述测量旋转轴的旋转一致性的测量方法是：使样本***固定在某个位置，金属球校准模型随旋转轴旋转，均匀地取N个旋转采样点并取其在探测器上的投影图像，N为偶数；比较M次旋转中，相同位置金属小球的投影位置和旋转中心投影位置的偏移，根据偏移量即可知道旋转轴的旋转一致性，M为自然数。
优选的，旋转轴在垂直移动过程中的一致性的测量方法为：在垂直方向，样本***上下移动使金属球校准模型在旋转过程中金属小球的投影轨迹完整落在探测器上，在这个垂直范围内选取V个采样点，V为偶数，其中V/2个垂直采样点位于中平面上方，V/2个垂直采样点位于中平面下方；在每个垂直采样点，金属球校准模型随旋转轴旋转，在N个均匀间隔的位置取其投影，N为偶数，计算每个旋转采样点金属小球在探测器上的投影位置，并计算旋转轴中心点在探测器上的投影位置，比较V个垂直采样点计算所得的旋转中心点的偏移，即可知道旋转轴在垂直移动过程中的一致性。
优选的，旋转轴在水平移动过程中的一致性的测量方法为：在水平方向，样本***水平移动范围内取H个水平采样点，H为自然数；在每个水平采样点，分别测量该水平位置旋转轴垂直移动的一致性，根据该水平位置V个垂直采样点计算的旋转轴中心投影，V为偶数，计算该水平位置射线源焦点在探测器上的投影点；比较H个水平采样点位置计算得到的射线源焦点的投影的偏移，即可知道旋转轴在水平移动过程中的一致性。
优选的，所述的立式锥束CT***采用FDK算法进行重建； FDK算法对立式锥束CT***的成像几何结构的要求是：旋转轴与中平面垂直，同时与探测器像元阵列的列方向平行、行方向垂直；过X射线源焦点做垂直探测器的直线，该直线经过旋转轴；旋转轴的旋转保持稳定，在旋转过程中始终保持与中平面垂直；X射线源到探测器的垂线的垂足是用于重建的投影图像的中心。
优选的，金属球校准模型是一个有机玻璃制的圆柱形杯子，杯子的上边缘固定一个金属球。
有益效果：  
本发明提出的一种立式锥束CT***校***及方法，在样本***水平运动范围内取若干个样本点，每个样本位置测量载物台垂直方向移动的一致性，获取每个位置旋转中心曲线。在同一平面中显示所有水平采样位置的旋转中心曲线。需要注意的是，样本***从靠近探测器位置向靠近射线源位置移动的过程中，***成像的放大倍数从小变大，探测器上一个像素所表示的实际长度从大变小。所以在比较各个水平位置旋转中心偏移前，需要将旋转中心位置转换成相同的数量单位。当***几何结构存在微小的偏差的时候，在最靠近探测器位置估算的射线源焦点在探测器上的投影位置是最接近真实值的，而在最靠近射线源位置估算的焦点投影值偏离真实值最远，比较两个极限位置的偏差就可以确定射线源偏移的方向和大小。
锥束CT***中水平位置旋转中心的偏移表明***中射线源、探测器与载物台的水平导轨的相对几何关系不满足FDK算法要求的几何结构。根据***成像几何结构的测量，对射线源或CCD进行相应的微调。调整后重新测量***的成像几何结构，验证之前的调整方向是否正确。多次测量与调整至无法通过机械方法精确***的成像几何结构。
附图说明
图1：立式锥束CT***的扫描结构；
图2a：立式锥束CT***校准所用的金属球校准模型的侧面；
图2b：立式锥束CT***校准所用的金属球校准模型的俯视图；
图3：金属球校准模型的投影示意图；
图4：立式锥束CT***校准过程的流程图；
具体实施方式
下面结合附图对本发明做进一步说明。
为了更详尽地表述上述发明的技术方案，以下列举出具体的实施例来说明技术效果。需要强调的是，这些实施例是用于说明本发明而不限制本发明的范围。
本发明提出了一种立式锥束CT校准***及应用于该***的方法，尤其是射线源与探测器固定不动、载物台随旋转轴旋转，载物台与旋转轴构成的整体可以上下左右移动的锥束CT***的成像几何结构的校准方法。
特别指出的是该锥束CT***的重建算法采用FDK算法，该算法是由Feldkamp，Davis和Kress在1984年提出的基于锥束投影的滤波反投影重建算法，参考他们在《美国光学会志，A辑：光学、图像科学与视觉》1984年6月第1卷第6期612至619页发表的文章《实用的锥束重建算法》。
FDK重建算法对***的成像几何结构有一定的要求：旋转轴与锥束射线源的中平面垂直，同时与探测器像元阵列的列方向平行、行方向垂直；过锥束射线源焦点做垂直探测器的直线，该直线经过旋转轴；旋转轴的旋转保持稳定，在旋转过程中始终保持与中平面垂直；射线源到探测器的垂线的垂足是用于重建的投影图像的中心。
本发明提出了一种立式锥束CT***校准方法，尤其是射线源与探测器固定不动、载物台随旋转轴旋转上下左右移动的锥束CT***的成像几何结构的校准方法，特别指出的是该锥束CT***的重建算法采用FDK算法。FDK重建算法对***的成像几何结构有一定的要求：旋转轴与锥束射线源的中平面垂直，同时与探测器像元阵列的列方向平行、行方向垂直；过锥束射线源焦点做垂直探测器的直线，该直线经过旋转轴；旋转轴的旋转保持稳定，在旋转过程中始终保持与中平面垂直；射线源到探测器的垂线的垂足是用于重建的投影图像的中心。
参见图1－图4，本发明提供的立式锥束CT***校准系包括X射线源100、与X射线源100相距一定距离的探测器106和设置在X射线源100与探测器106之间的载物台102、用于固定载物台102的旋转轴315、固定在载物台102上的金属球校准模型104。
本发明中所说的中平面是指：立式锥束CT校准***中与探测器106垂直且经过X射线源100的平面称为中平面。
其中金属球校准模型104在扫描过程中与载物台102一起随旋转轴315旋转；X射线源100与探测器106的位置固定不动，通过载物台102旋转获取投影数据；金属球校准模型104上固定一金属球105；金属球校准模型104、载物台102和旋转轴315构成样本***201。
金属球校准模型104是有机玻璃制成的圆柱形杯子，杯子的上边缘固定金属球105。
本发明提供的立式锥束CT校准方法包括以下步骤：测量旋转轴315在旋转过程中的旋转中心；测量旋转轴315的旋转一致性；所述旋转一致性是指使样本***201固定在某个位置，使金属球校准模型104位于合适的旋转采样点，金属球校准模型104旋转360°M次，M为自然数，分别在每次旋转后计算金属小球105的位置，金属小球105位置的一致性就是旋转轴的旋转一致性；测量315旋转轴在垂直移动过程中的一致性；所述旋转轴315在垂直移动过程中的一致性是指在旋转轴上下移动的过程中旋转中心的一致性；测量旋转轴315在水平移动过程中的一致性，所述旋转轴315在水平移动过程中的一致性是指在旋转轴前后移动的过程中旋转中心的一致性；根据上述数实现立体锥束CT校准***的校准。
所述旋转轴的旋转中心的测量方法是：使样本***201固定在某个位置，金属球校准模型104随旋转轴315旋转，均匀地取N个旋转采样点并取其在探测器106上的投影图像；根据金属球校准模型104旋转一周的投影数据计算旋转轴的旋转中心，N为偶数。
所述测量旋转轴的旋转一致性的测量方法是：使样本***201固定在某个位置，金属球校准模型104随旋转轴315旋转，均匀地取N个旋转采样点并取其在探测器106上的投影图像，N为偶数；比较M次旋转中，相同位置金属小球105的投影位置和旋转中心投影位置的偏移，根据偏移量即可知道旋转轴的旋转一致性，M为自然数。
旋转轴在垂直移动过程中的一致性的测量方法为：在垂直方向，样本***201上下移动使金属球校准模型104在旋转过程中金属小球105的投影轨迹完整落在探测器106上，在这个垂直范围内选取V个采样点，V为偶数，其中V/2个垂直采样点位于中平面上方，V/2个垂直采样点位于中平面下方；在每个垂直采样点，金属球校准模型104随旋转轴315旋转，在N个均匀间隔的位置取其投影，N为偶数，计算每个旋转采样点金属小球105在探测器106上的投影位置，并计算旋转轴中心点在探测器上的投影位置，比较V个垂直采样点计算所得的旋转中心点的偏移，即可知道旋转轴在垂直移动过程中的一致性。
旋转轴在水平移动过程中的一致性的测量方法为：在水平方向，样本***201水平移动范围内取H个水平采样点，H为自然数；在每个水平采样点，分别测量该水平位置旋转轴垂直移动的一致性，根据该水平位置V个垂直采样点计算的旋转轴中心投影，V为偶数，计算该水平位置射线源焦点在探测器上的投影点；比较H个水平采样点位置计算得到的射线源焦点的投影的偏移，即可知道旋转轴在水平移动过程中的一致性。
所述的立式锥束CT***采用FDK算法进行重建； FDK算法对立式锥束CT***的成像几何结构的要求是：旋转轴315与中平面垂直，同时与探测器106像元阵列的列方向平行、行方向垂直；过X射线源100焦点做垂直探测器106的直线，该直线经过旋转轴315；旋转轴315的旋转保持稳定，在旋转过程中始终保持与中平面垂直；X射线源到探测器的垂线的垂足是用于重建的投影图像的中心。
金属球校准模型104是一个有机玻璃制的圆柱形杯子，杯子的上边缘固定一个金属球105。
详细如下：图1是立式锥束CT***的扫描结构，立式锥束CT***的组成部分包括微焦斑X射线源100，载物台102，金属球校准模型104和探测器106。金属球校准模型104固定在载物台102上，载物台102固定在旋转轴315上。该***中金属球校准模型104和载物台102可随旋转轴315旋转。金属球校准模型104、载物台102和旋转轴315构成样本***201。样本***201可沿着水平导轨前后移动，沿垂直导轨上下移动。X射线源100与探测器106的位置在***安装完成后固定不变，但可以根据需要进行微调。
图2a是立式锥束CT***校准所用的金属球校准模型104的侧面；图2b是立式锥束CT***校准所用的金属球校准模型104的俯视图。金属球校准模型104是一个有机玻璃制的圆柱形杯子，杯子的上边缘有一个金属小球105。X射线源100、样本***201和探测器106之间的相对位置不同，***的放大倍数不一样，即当样本***201靠近探测器106时***放大倍数较小，而当样本***201靠近X射线源100时***放大倍数较大。因此，为了让金属校准模型104的投影图像尽量填满探测器106平面，样本***201位于不同水平位置时不同直径的金属校准模型104，越靠近X射线源100的位置，金属球校准模型104杯体的直径就越小。
图3是金属球校准模型在旋转过程的小球投影示意图。金属小球105相对于金属校准模型104尺寸很小，可以视为质点。金属校准模型104与载物台102一起随旋转轴315旋转过程中，小球质点312的运动轨迹是一个圆。在锥束投影中，如果质点在圆形轨迹上旋转，其投影314的轨迹就是一个椭圆。圆形轨迹中心的投影一定在旋转轴的投影线上。在小球质点312旋转的过程中，相隔180°的两个位置的小球质点连线一定经过旋转中心。由于锥束投影中一条直线在探测器上的投影仍然是一条直线，所以相隔180°的两个位置的小球质点的投影在探测器上的连线必定经过圆轨迹中心在探测器上的投影。小球质点经过360°的旋转，在N（N为偶数）个均匀间隔的位置取其投影。相隔180°的两个位置的小球质点的投影连线经过旋转轴中心的投影，通过最小二乘方法可以根据N/2条直线方程求得旋转轴的旋转中心。
图4是立式锥束CT***的校准过程的流程图。
立式锥束CT***的校准过程主要包括旋转轴的旋转中心的测量；旋转轴的旋转一致性；旋转轴在垂直移动过程中的一致性；旋转轴在水平移动过程中的一致性。
旋转轴的旋转中心的测量方法如下：金属球校准模型104与载物台102随旋转轴315旋转，每隔360/N°采集一幅投影图，计算金属小球105在探测器上的投影点。
计算小球在探测器上投影点的方法：在投影图中，由于金属球校准模型104的有机玻璃杯和背景的灰度较高，金属小球105的灰度较低，可以使用阈值分割法将金属球分割出来。然后将分割出来的金属球范围的数据保存，并将灰度值归一化，即使图像中所有像素的灰度值分布在0-255之间。设分割出的金属球区域的点坐标为                                                
Figure 13771DEST_PATH_IMAGE001
，对应像素的灰度值为
Figure 936727DEST_PATH_IMAGE002
，其中
Figure 246486DEST_PATH_IMAGE003
，P是分割出的金属球区域的像素个数。令金属小球的投影点坐标为
Figure 753428DEST_PATH_IMAGE004
，则该坐标可以通过如下公式计算得到：
       
Figure 806835DEST_PATH_IMAGE005
，
Figure 584298DEST_PATH_IMAGE006
                               （1）
计算旋转中心的计算方法：金属球校准模型104旋转360°后获得N个均匀间隔的小球投影。用表示小球质点在第i个位置的投影点坐标，相隔180°的两个位置的小球质点的投影点为
Figure 295082DEST_PATH_IMAGE001
和
Figure 152180DEST_PATH_IMAGE007
，其中
Figure 784149DEST_PATH_IMAGE008
，通过这两个点的直线可表示为：
Figure 435711DEST_PATH_IMAGE009
    
Figure 917245DEST_PATH_IMAGE010
        （2）
其中表示旋转中心的投影点坐标。
由于N是偶数，这样的方程共有N/2个，求解这些方程的最小二乘解就可以得到旋转中心的投影点坐标
Figure 798931DEST_PATH_IMAGE011
。具体做法为：将等式（2）改写为如下形式：
Figure 824655DEST_PATH_IMAGE012
                                                       （3）
其中，
       ，
Figure 228272DEST_PATH_IMAGE014
，
Figure 631571DEST_PATH_IMAGE010
          （4）
则方程组可以写为：
                                                            (5)
方程组（5）称为超定方程组，解此方程的最小二乘解即可得旋转中心的投影点坐标
Figure 98640DEST_PATH_IMAGE011
。
旋转轴的旋转一致性的测量方法如下：金属球校准模型104经过360°的旋转，在N（N为偶数）个均匀间隔的位置取其投影。样本***201固定在某个位置，金属球校准模型104旋转360°M（M为自然数）次，每次旋转中分别保存N个位置的小球质点的投影位置并计算该位置旋转中心的投影位置。比较M次旋转中，相同位置小球质点投影位置和旋转中心投影位置的偏移，根据分析所得的偏移量即可知道旋转轴的旋转一致性。
旋转轴在垂直移动过程中的一致性的测量方法如下：在垂直方向，上下移动样本***201使小球质点在旋转过程的投影轨迹完整落在探测器上，在这个垂直范围内选取V（V为偶数）个采样点，其中V/2个垂直采样点位于射线源中平面上面，V/2个垂直采样点位于射线源中平面下面。在每个垂直采样点，金属球校准模型104经过360°的旋转，在N（N为偶数）个均匀间隔的位置取其投影。记录每个旋转采样点小球质点在探测器上的投影位置，并计算旋转轴中心点在探测器上的投影位置。理想情况下，V个垂直采样点旋转轴中心的投影点的连线是一条平行探测器平面像元阵列的列方向平行、行方向垂直的直线。实际测量的结果，V个垂直采样点旋转轴中心的投影点的连线与探测器平面像元阵列的列方向有一个夹角，即说明旋转轴在样本***201上下移动过程中的不稳定性，其旋转中心沿着某个方向偏移。
旋转轴在水平移动过程中的一致性的测量方法如下：在水平方向，样本***201前后可移动范围内取H（H为自然数）个水平采样点，。在每个水平采样点，分别测量该水平位置旋转轴垂直移动的一致性，根据该水平位置V（V为偶数）个旋转轴中心投影，计算该水平位置射线源焦点在探测器上的投影点。
计算射线源焦点在探测器上的投影点的方法：样本***在某个水平位置测量计算的旋转中心点投影坐标表示为
Figure 773335DEST_PATH_IMAGE016
，其中
Figure 296720DEST_PATH_IMAGE017
。该位置射线源焦点在探测器上的投影点坐标是
Figure 664248DEST_PATH_IMAGE018
，则该坐标可以通过如下公式计算得到：
Figure 640294DEST_PATH_IMAGE019
，
Figure 118680DEST_PATH_IMAGE020
                                                                                     （6）
其中，
       
Figure 762151DEST_PATH_IMAGE021
，
Figure 533536DEST_PATH_IMAGE022
                               （7）
       其中
Figure 262457DEST_PATH_IMAGE023
表示第i个垂直采样位置小球质点椭圆轨迹的短轴。小球质点在旋转过程中其投影的椭圆形轨迹的短轴随着样本***的位置变化而变化：样本***越靠近中平面，小球质点的椭圆轨迹的短轴越短；样本***越远离中平面，小球质点的椭圆轨迹的短轴越长。所以，沿探测器的列方向，第i个旋转中心的投影点与射线源焦点的投影点的距离与第i个小球质点椭圆轨迹的短轴成正比，如等式（7）所述。由于V是偶数且V/2个垂直采样点分别位于中平面的两边，根据等式（7）可以计算得到V/2个射线源焦点投影纵坐标
Figure 278955DEST_PATH_IMAGE024
，再取平均就得到更加准确的射线源焦点投影纵坐标
Figure 776932DEST_PATH_IMAGE025
。
理论上，所有H个水平采样位置的射线源焦点投影点的连线与锥束X射线的中心线平行、与探测器平面垂直。实际测量的结果，所有H个水平采样位置的射线源焦点投影点的连线与探测器平面不垂直而是有一定的夹角。根据这个夹角和载物台相对射线源和探测器的相对位置，可以估算射线源的偏移位置。

Claims (9)

1.一种立式锥束计算机断层成像校准***，其特征在于：该***包括X射线源（100）、与X射线源（100）相对设置且不接触的探测器（106）和设置在X射线源（100）与探测器（106）之间的载物台（102）、用于固定载物台（102）的旋转轴（315）、固定在载物台（102）上的金属球校准模型（104）；其中金属球校准模型（104）在扫描过程中与载物台（102）一起随旋转轴（315）旋转；X射线源（100）与探测器（106）的位置固定不动，通过载物台(102)旋转获取投影数据；金属球校准模型（104）上固定一金属球（105）；金属球校准模型（104）、载物台（102）和旋转轴（315）构成样本***（201）。
2.根据权利要求1所述的立式锥束计算机断层成像校准***，其特征在于：金属球校准模型（104）是有机玻璃制成的圆柱形杯子，杯子的上边缘固定金属球（105）。
3.一种用于权利要求1所述的***的立式锥束计算机断层成像校准方法，其特征在于：该方法包括以下步骤： 测量旋转轴（315）在旋转过程中的旋转中心；测量旋转轴（315）的旋转一致性；所述旋转一致性是指使样本***（201）固定在某个位置，使金属球校准模型（104）位于合适的旋转采样点，金属球校准模型（104）旋转360°M次，M为自然数，分别在每次旋转后计算金属小球（105）的位置，金属小球（105）位置的一致性就是旋转轴的旋转一致性；测量（315）旋转轴在垂直移动过程中的一致性；所述旋转轴（315）在垂直移动过程中的一致性是指在旋转轴上下移动的过程中旋转中心的一致性；测量旋转轴（315）在水平移动过程中的一致性，所述旋转轴（315）在水平移动过程中的一致性是指在旋转轴前后移动的过程中旋转中心的一致性；根据上述数实现立体锥束CT校准***的校准。
4.根据权利要求3所述的立式锥束计算机断层成像***校准方法，其特征在于：所述旋转轴的旋转中心的测量方法是：使样本***（201）固定在某个位置，金属球校准模型（104）随旋转轴（315）旋转，均匀地取N个旋转采样点并取其在探测器（106）上的投影图像；根据金属球校准模型（104）旋转一周的投影数据计算旋转轴的旋转中心，N为偶数。
5.根据权利要求3所述的立式锥束计算机断层成像***校准方法，其特征在于：所述测量旋转轴的旋转一致性的测量方法是：使样本***（201）固定在某个位置，金属球校准模型（104）随旋转轴（315）旋转，均匀地取N个旋转采样点并取其在探测器（106）上的投影图像，N为偶数；比较M次旋转中，相同位置金属小球（105）的投影位置和旋转中心投影位置的偏移，根据偏移量即可知道旋转轴的旋转一致性，M为自然数。
6.根据权利要求3所述的立式锥束计算机断层成像***校准方法，其特征在于：旋转轴在垂直移动过程中的一致性的测量方法为：在垂直方向，样本***（201）上下移动使金属球校准模型（104）在旋转过程中金属小球（105）的投影轨迹完整落在探测器（106）上，在这个垂直范围内选取V个采样点，V为偶数，其中V/2个垂直采样点位于中平面上方，V/2个垂直采样点位于中平面下方；在每个垂直采样点，金属球校准模型（104）随旋转轴（315）旋转，在N个均匀间隔的位置取其投影，N为偶数，计算每个旋转采样点金属小球（105）在探测器（106）上的投影位置，并计算旋转轴中心点在探测器上的投影位置，比较V个垂直采样点计算所得的旋转中心点的偏移，即可知道旋转轴在垂直移动过程中的一致性。
7.根据权利要求3所述的立式锥束计算机断层成像***校准方法，其特征在于：旋转轴在水平移动过程中的一致性的测量方法为：在水平方向，样本***（201）水平移动范围内取H个水平采样点，H为自然数；在每个水平采样点，分别测量该水平位置旋转轴垂直移动的一致性，根据该水平位置V个垂直采样点计算的旋转轴中心投影，V为偶数，计算该水平位置射线源焦点在探测器上的投影点；比较H个水平采样点位置计算得到的射线源焦点的投影的偏移，即可知道旋转轴在水平移动过程中的一致性。
8.根据权利要求3所述的立式锥束计算机断层成像***校准方法，其特征在于：所述的立式锥束CT***采用FDK算法进行重建； FDK算法对立式锥束CT***的成像几何结构的要求是：旋转轴（315）与中平面垂直，同时与探测器（106）像元阵列的列方向平行、行方向垂直；过X射线源（100）焦点做垂直探测器（106）的直线，该直线经过旋转轴（315）；旋转轴（315）的旋转保持稳定，在旋转过程中始终保持与中平面垂直；X射线源到探测器的垂线的垂足是用于重建的投影图像的中心。
9.根据权利要求3到8任一权利要求所述的立式锥束计算机断层成像***校准方法，其特征在于：金属球校准模型（104）是一个有机玻璃制的圆柱形杯子，杯子的上边缘固定一个金属球（105）。
CN 201110097958 2011-04-19 2011-04-19 立式锥束计算机断层成像校准***及应用该***的方法 Expired - Fee Related CN102253061B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110097958 CN102253061B (zh) 2011-04-19 2011-04-19 立式锥束计算机断层成像校准***及应用该***的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110097958 CN102253061B (zh) 2011-04-19 2011-04-19 立式锥束计算机断层成像校准***及应用该***的方法

Publications (2)

Publication Number Publication Date
CN102253061A true CN102253061A (zh) 2011-11-23
CN102253061B CN102253061B (zh) 2013-07-10

Family

ID=44980450

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110097958 Expired - Fee Related CN102253061B (zh) 2011-04-19 2011-04-19 立式锥束计算机断层成像校准***及应用该***的方法

Country Status (1)

Country Link
CN (1) CN102253061B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103226113A (zh) * 2013-03-29 2013-07-31 中国计量科学研究院 锥束3d-ct扫描***重建体素尺寸的自动标定方法
CN106706675A (zh) * 2015-07-16 2017-05-24 中国科学院高能物理研究所 一种基于计算机分层扫描成像cl***的校正方法
CN107507257A (zh) * 2017-08-11 2017-12-22 重庆真测科技股份有限公司 一种大视场分度偏移重建方法及其***
WO2018126335A1 (zh) * 2017-01-03 2018-07-12 苏州海斯菲德信息科技有限公司 基于小球模体的锥束ct***几何参数评价及校正方法
WO2019095838A1 (zh) * 2017-11-14 2019-05-23 深圳先进技术研究院 一种静态锥束ct成像***几何校准装置及方法
CN110786875A (zh) * 2018-08-02 2020-02-14 株式会社岛津制作所 放射线图像处理装置及放射线图像处理方法
CN113945174A (zh) * 2021-10-21 2022-01-18 中国工程物理研究院激光聚变研究中心 一种x射线投影测量图像尺寸校准方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1096426A1 (en) * 1999-10-27 2001-05-02 GE Medical Systems Global Technology Company LLC Methods and apparatus for cone beam multi-slice CT correction
JP2005241575A (ja) * 2004-02-27 2005-09-08 Toshiba Corp Ｘ線断層撮影装置および立体透視画像構成装置
CN101832954A (zh) * 2010-03-15 2010-09-15 中国工程物理研究院应用电子学研究所 锥束xct***用移动组件以及用其进行图像重建坐标系原点标定的方法
CN101936720A (zh) * 2010-07-30 2011-01-05 北京航空航天大学 一种适用于锥束xct***的探测器扭转角的标定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1096426A1 (en) * 1999-10-27 2001-05-02 GE Medical Systems Global Technology Company LLC Methods and apparatus for cone beam multi-slice CT correction
JP2005241575A (ja) * 2004-02-27 2005-09-08 Toshiba Corp Ｘ線断層撮影装置および立体透視画像構成装置
CN101832954A (zh) * 2010-03-15 2010-09-15 中国工程物理研究院应用电子学研究所 锥束xct***用移动组件以及用其进行图像重建坐标系原点标定的方法
CN101936720A (zh) * 2010-07-30 2011-01-05 北京航空航天大学 一种适用于锥束xct***的探测器扭转角的标定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
胡栋材等: "基于平板探测器的锥束CT散射校正方法", 《CT理论与应用研究》 *
黄魁东等: "基于配准模型仿真的锥束CT射束硬化校正方法", 《***仿真学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103226113A (zh) * 2013-03-29 2013-07-31 中国计量科学研究院 锥束3d-ct扫描***重建体素尺寸的自动标定方法
CN103226113B (zh) * 2013-03-29 2015-10-21 中国计量科学研究院 锥束3d-ct扫描***重建体素尺寸的自动标定方法
CN106706675A (zh) * 2015-07-16 2017-05-24 中国科学院高能物理研究所 一种基于计算机分层扫描成像cl***的校正方法
WO2018126335A1 (zh) * 2017-01-03 2018-07-12 苏州海斯菲德信息科技有限公司 基于小球模体的锥束ct***几何参数评价及校正方法
CN107507257A (zh) * 2017-08-11 2017-12-22 重庆真测科技股份有限公司 一种大视场分度偏移重建方法及其***
CN107507257B (zh) * 2017-08-11 2021-01-29 重庆真测科技股份有限公司 一种大视场分度偏移重建方法及其***
WO2019095838A1 (zh) * 2017-11-14 2019-05-23 深圳先进技术研究院 一种静态锥束ct成像***几何校准装置及方法
CN110786875A (zh) * 2018-08-02 2020-02-14 株式会社岛津制作所 放射线图像处理装置及放射线图像处理方法
CN110786875B (zh) * 2018-08-02 2024-03-12 株式会社岛津制作所 放射线图像处理装置及放射线图像处理方法
CN113945174A (zh) * 2021-10-21 2022-01-18 中国工程物理研究院激光聚变研究中心 一种x射线投影测量图像尺寸校准方法
CN113945174B (zh) * 2021-10-21 2023-10-17 中国工程物理研究院激光聚变研究中心 一种x射线投影测量图像尺寸校准方法

Also Published As

Publication number Publication date
CN102253061B (zh) 2013-07-10

Similar Documents

Publication Publication Date Title
CN102253061B (zh) 立式锥束计算机断层成像校准***及应用该***的方法
CN104783824B (zh) X射线成像***的校正方法
CN103784160B (zh) 一种锥束ct***几何位置的校正装置及其校正方法
CN111248934B (zh) 一种用于cbct***机械校正的方法及***
CN103735282B (zh) 一种锥束ct***探测器几何校正装置及其校正方法
CN102652674B (zh) 一种消除ct图像中的几何伪影的方法和***
CN103800032B (zh) 用于锥束ct***几何位置校正的校正***及其校正方法
CN101936720B (zh) 一种适用于锥束xct***的探测器扭转角的标定方法
CN100526829C (zh) 一种杂光测试仪
CN203776924U (zh) 一种锥束ct***几何位置的校正装置
CN103759679B (zh) 一种锥束ct***角度偏差测量方法
CN103226113B (zh) 锥束3d-ct扫描***重建体素尺寸的自动标定方法
CN209032406U (zh) 一种锥束ct***几何校准装置
CN102488528B (zh) 一种层析成像几何参数的校准方法
CN202104929U (zh) 一种用于体ct几何校正的装置
CN105844594A (zh) 一种工业ct扫描***旋转中心校正方法
CN106706675B (zh) 一种基于计算机分层扫描成像cl***的校正方法
CN106821405A (zh) 一种x光机的参数标定方法、装置及***
CN105318844A (zh) 球形状测定方法和装置
JP5537226B2 (ja) 放射線撮像装置
CN107016655A (zh) 锥束cl几何全参数迭代校正方法
CN113812971A (zh) 一种多自由度四维双能锥束ct成像***及方法
CN114295650A (zh) X射线ct检测装置及检测方法
WO2015093248A1 (ja) Ｃ型アームを用いたｘ線断層撮影装置及びその制御方法
AU2019201125B2 (en) Geometric alignment, sample motion correction, and intensity normalization of computed tomography projections using pi-line optimization

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130710

Termination date: 20160419

CF01 Termination of patent right due to non-payment of annual fee