CN107233105A - 一种用于ct图像重建的修正方法及修正系统 - Google Patents
一种用于ct图像重建的修正方法及修正系统 Download PDFInfo
- Publication number
- CN107233105A CN107233105A CN201710373329.6A CN201710373329A CN107233105A CN 107233105 A CN107233105 A CN 107233105A CN 201710373329 A CN201710373329 A CN 201710373329A CN 107233105 A CN107233105 A CN 107233105A
- Authority
- CN
- China
- Prior art keywords
- mrow
- detector
- mtr
- mtd
- coordinate system
- 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
Links
- 238000002715 modification method Methods 0.000 title claims abstract description 16
- 239000013598 vector Substances 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000001914 filtration Methods 0.000 claims description 9
- 238000009795 derivation Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 239000007787 solid Substances 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 4
- 238000013170 computed tomography imaging Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 1
- 201000006549 dyspepsia Diseases 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种用于CT图像重建的修正方法及修正系统。提出了一种用于CT图像重建的修正方法,其特征在于根据探测器偏置角度φ,得到探测器上的任意点在所述固定坐标系XYZ中的新坐标。XYZ坐标系的空间坐标原点O与CT系统中的探测器上任意点所构成的矢量用表示为:在探测器偏置角度φ后,空间坐标原点O与CT系统中的探测器上任意点所构成的矢量用表示为:从而得到偏置后的探测器上任意点在XYZ坐标系中的新坐标。所述方法可用于对CT系统几何校准方法以及CT重建算法进行修正。本发明提出的修正方法和修正系统可以显著改善CT成像质量,可用于医学和工业CT成像。
Description
技术领域
本发明涉及医学和工业CT成像,更具体地涉及一种用于CT图像重建的修正方法及修正系统。
背景技术
计算机断层成像(Computed Tomography,简称CT)是用X射线对人体的特定部位按一定厚度的层面进行扫描,由于不同的人体组织对X射线的吸收能力不同,可以用计算机重建出断层面的影像。CT是通过无损方式获取物体内部结构信息的一种重要成像手段,拥有高分辨率、高灵敏度以及多层次等众多优点,被广泛应用于各个医疗临床检查领域。
CT图像重建技术是CT成像系统的核心技术之一,涉及到物理学、数学、计算机图形学等多个学科领域。通过投影重建图像的理论与算法是CT的核心理论与基础算法。随着CT成像技术的应用越来越广泛,人们对CT重建图像质量的要求也越来越高,这就需要开发更先进的重建算法,以适应不同情况下的投影数据,从而重建出更高质量的图像。
CT图像重建算法在定义坐标系时,需要知道各个几何参数的准确数值。CT重建时需要如下7个几何参数:光源在探测器平面的投影中心点坐标(横坐标和纵坐标)、探测器的扭转角、探测器的倾斜角、探测器的旋转角、光源—旋转轴Z的距离(R)、光源—探测器的距离(D),如图1所示。如果几何参数存在误差,会严重影响CT重建图像质量。而几何校准可以获得CT系统几何参数信息。因此,几何校准是CT系统调试的重要组成部分,也是获得精确CT重建图像的前提条件。然而,在实际CT系统中,探测器经常存在偏置角度,而几何校准方法要么没有考虑偏置角度,要么对偏置角度的校准误差较大。
发明内容
本发明的解决方案的目的是解决前面强调的问题。
因此,本专利提出了一种探测器偏置后探测器上的任意点在XYZ坐标系中的新坐标的求解方法以及系统,可用于对CT系统几何校准方法以及CT重建算法进行修正。
本发明的一方面提供了一种用于CT图像重建的修正方法,CT系统包括光源和探测器。在所述方法中定义固定坐标系XYZ,其中Z为旋转轴,R为光源至旋转轴Z的距离,D为光源至探测器的距离,U和V为探测器上像素的行和列方向,分别与坐标轴Y和Z平行。在固定坐标系XYZ中,旋转坐标系中的三个单位矢量可以表示为:
并且所述固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的矢量表示为:所述修正方法的特征在于:根据探测器偏置角度φ,得到探测器上的任意点在所述固定坐标系XYZ中的新坐标。
在一些实施方式中,所述探测器偏置角度φ为探测器围绕U轴旋转角度φ。
在一些实施方式中,在探测器围绕U轴旋转角度φ后,旋转坐标系的三个单位矢量变更为:
在一些实施方式中,在探测器围绕U轴旋转角度φ后的用表示为:
在一些实施方式中,在探测器围绕U轴旋转角度φ后,探测器上的任意点在固定坐标系XYZ中的新坐标为:
在一些实施方式中,角度φ在-90度到90度之间。
在一些实施方式中,角度φ在-5度到5度之间。
本发明的另一方面提供了一种CT重建算法,其特征在于包括根据权利要求1所述的用于CT图像重建的修正方法。
在一些实施方式中,所述CT重建算法包括反投影滤波(BPF)算法,所述反投影滤波(BPF)算法包括以下步骤:对探测器上获得的物体图像(序列集f)进行反log操作:以得到处理后的投影图像;对每个角度下的锥束投影数据进行求导;由求导后的投影数据,对螺旋轨道内每条PI线上的物体图像点进行加权反投影;对加权反投影得到的PI线上的数据,沿着该条PI线求Hilbert变换,得到物体在该线上的重建结果;以及对螺旋轨道内所有PI线上的物体重建结果进行重采样,得到物体的三维重建图像。
在一些实施方式中,所述CT重建算法包括滤波反投影(FBP)算法,所述滤波反投影(FBP)算法包括以下步骤:对探测器上获得的物体图(序列集f)像进行反log操作:以得到处理后的投影图像;以及对所述投影图像进行加权、卷积滤波以及加权反投影,以获得CT图像。
本发明的另一方面提供了一种CT系统几何校准方法,其特征在于包括根据权利要求1所述的用于CT图像重建的修正方法。
在一些实施方式中,所述CT系统几何校准方法的特征在于选自:基于小球体模投影椭圆方程的CT系统几何校准方法、以及基于投影矩阵的CT系统几何校准方法。
本发明的另一方面提供了一种用于CT图像重建的修正系统,其可以包括:固定坐标系定义模块,其用于定义固定坐标系XYZ,并获取旋转坐标系中的单位矢量以及固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的矢量的表示;修正模块,其用于获取修正后的旋转坐标系中的单位矢量以及修正后的坐标原点O与探测器上的任意点所构成的矢量的表示;新坐标获取模块,其用于根据修正后的所述矢量的表示得到探测器上的任意点在所述固定坐标系XYZ中的新坐标。
在一些实施方式中,所述固定坐标系定义模块定义固定坐标系XYZ,其中Z为旋转轴,U轴和V轴为探测器上像素的行和列方向,分别与坐标轴Y和Z平行,在固定坐标系XYZ中,旋转坐标系中的所述单位矢量表示为:
并且
所述固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的所述矢量表示为:
其中,R为光源至旋转轴Z的距离,D为光源至探测器的距离
在一些实施方式中,所述修正模块根据探测器偏置角度φ获得修正后的单位矢量以及矢量探测器偏置角度φ为探测器围绕U轴旋转角度φ,在探测器围绕U轴旋转角度φ后,旋转坐标系的所述单位矢量变更为:
所述矢量表示为:
在一些实施方式中,所述新坐标获取模块得到探测器上的任意点在所述固定坐标系XYZ中的新坐标为:
本发明提出的修正方法、CT重建算法、CT系统几何校准方法以及修正系统可以显著改善CT成像质量,可用于医学和工业CT成像。
本领域技术人员在阅读整个说明书和权利要求书时将理解本发明的这些优点和其它优点。
附图说明
为了更好地理解本发明,现在仅通过非限制性示例的方式参照所附附图描述其优选实施方式,其中:
图1示出了CT系统几何结构示意图。
图2示出了探测器倾斜后的CT系统几何结构示意图。
图3示出了探测器倾斜后的CT系统几何结构的俯视图。
图4示出了探测器倾斜后的CT系统几何结构的侧视图。
具体实施方式
下面结合附图对本发明的具体实施例进行说明。在下文所描述的本发明的具体实施例中,为了能更好地理解本发明而描述了一些很具体的技术特征,但显而易见的是,对于本领域的技术人员来说,并不是所有的这些技术特征都是实现本发明的必要技术特征。下文所描述的本发明的一些具体实施例只是本发明的一些示例性的具体实施例,其不应被视为对本发明的限制。另外,为了避免使本发明变得难以理解,对于一些公知的技术没有进行描述。
在CT成像系统中,由CT系统空间几何关系可以得到公式(1):
a[u,v,1]T=A[x,y,z,1]T (1)
其中A为系统矩阵,a为系数;空间任意点坐标(x,y,z)在探测器上的投影点坐标用(u,v)表示。本发明旨在推导当探测器倾斜φ角度时,探测器上任意点的更新坐标。
CT成像系统几何示意图如图1-2所示。在图1中,其中S为光源,光源为X射线发生装置,探测器为X射线接收装置。定义固定坐标系XYZ,Z为旋转轴,(u0,v0)表示光源在探测器平面的投影中心点坐标。R为光源至旋转轴Z的距离,D为光源至探测器的距离,U和V为探测器上像素的排列方向,分别与坐标轴Y和Z平行,定义探测器的左下角为坐标原点。因此,在固定坐标系XYZ中,旋转坐标系中的三个单位矢量可以表示为:
固定坐标系XYZ的坐标原点O与探测器上任意点所构成的矢量用表示,
现在参考图2-4,图2表示探测器围绕U轴(探测器行方向)旋转后的系统几何结构;图3是探测器倾斜后的CT系统几何结构的俯视图,θ表示ew与X轴的夹角;图4是探测器倾斜后的CT系统几何结构的侧视图。其中,eu、ev和ew分表表示探测器行方向、列方向和垂直方向的三个方向角。ev′表示探测器旋转后的列方向角。ev与ev′的夹角为φ,即探测器旋转角度为φ,如图4所示。
如图2所示,当探测器围绕U轴旋转φ角度后,三个单位矢量变更为:
其中,eu和ew保持不变,变化后的ev用ev′来表示。在一些实施方式中,角度φ在-90度到90度之间。在一些实施方式中,角度φ在-5度到5度之间。
空间坐标原点O与探测器上任意点所构成的矢量用表示,当探测器旋转φ角度后,变化后的用表示:
即,
因此,我们可以得到旋转后的探测器上的任意点在XYZ坐标系中的新坐标:
在重建算法中用新的x、y和z值来进行重建,从而修正由于探测器偏置所带来的误差。
本发明提供的用于CT图像重建的修正方法可以用于CT重建算法。在一些实施方式中,所述CT重建算法包括反投影滤波(BPF)算法,所述反投影滤波(BPF)算法包括以下步骤:对探测器上获得的物体图像(序列集f)进行反log操作:以得到处理后的投影图像;对每个角度下的锥束投影数据进行求导;由求导后的投影数据,对螺旋轨道内每条PI线上的物体图像点进行加权反投影;对加权反投影得到的PI线上的数据,沿着该条PI线求Hilbert变换,得到物体在该线上的重建结果;以及对螺旋轨道内所有PI线上的物体重建结果进行重采样,得到物体的三维重建图像。在一些实施方式中,所述CT重建算法包括滤波反投影(FBP)算法,所述滤波反投影(FBP)算法包括以下步骤:对探测器上获得的物体图(序列集f)像进行反log操作:以得到处理后的投影图像;以及对所述投影图像进行加权、卷积滤波以及加权反投影,以获得CT图像。
本发明提供的用于CT图像重建的修正方法还可以用于CT系统几何校准方法。在一些实施方式中,所述几何校准方法是基于小球体模投影椭圆方程的CT系统几何校准方法。在一些实施方式中,所述几何校准方法是基于投影矩阵的CT系统几何校准方法。
本发明还提供了一种用于CT图像重建的修正系统:其可以包括:固定坐标系定义模块,其用于定义固定坐标系XYZ,并获取旋转坐标系中的单位矢量以及固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的矢量的表示;修正模块,其用于获取修正后的旋转坐标系中的单位矢量以及修正后的坐标原点O与探测器上的任意点所构成的矢量的表示;新坐标获取模块,其用于根据修正后的所述矢量的表示得到探测器上的任意点在所述固定坐标系XYZ中的新坐标。
尽管已经根据优选的实施方案对本发明进行了说明,但是存在落入本发明范围之内的改动、置换以及各种替代等同方案。还应当注意的是,存在多种实现本发明的工艺的可选方式。因此,意在将随附的权利要求书解释为包含落在本发明的主旨和范围之内的所有这些改动、置换以及各种替代等同方案。
Claims (16)
1.一种用于CT图像重建的修正方法,CT系统包括光源和探测器;定义固定坐标系XYZ,其中Z为旋转轴,U轴和V轴为探测器上像素的行和列方向,分别与坐标轴Y和Z平行,在固定坐标系XYZ中,旋转坐标系中的三个单位矢量表示为:
并且
所述固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的矢量表示为:
其中,R为光源至旋转轴Z的距离,D为光源至探测器的距离;
所述修正方法的特征在于:根据探测器偏置角度φ,得到探测器上的任意点在所述固定坐标系XYZ中的新坐标。
2.根据权利要求1所述的方法,其中探测器偏置角度φ为探测器围绕U轴旋转角度φ。
3.根据权利要求2所述的方法,在探测器围绕U轴旋转角度φ后,旋转坐标系的三个单位矢量变更为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>e</mi>
<mi>v</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>e</mi>
<mi>v</mi>
</msub>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>e</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
4.根据权利要求3所述的方法,其中在探测器围绕U轴旋转角度φ后的用表示为:
5.根据权利要求4所述的方法,其中在探测器围绕U轴旋转角度φ后,探测器上的任意点在所述固定坐标系XYZ中的新坐标为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mo>-</mo>
<mi>u</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mi>R</mi>
<mo>+</mo>
<mi>v</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>u</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mi>R</mi>
<mo>+</mo>
<mi>v</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>z</mi>
<mo>=</mo>
<mi>v</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&phi;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
6.根据权利要求1所述的方法,其中角度φ在-90度到90度之间。
7.根据权利要求1所述的方法,其中角度φ在-5度到5度之间。
8.一种CT重建算法,其特征在于包括根据权利要求1所述的用于CT图像重建的修正方法。
9.根据权利要求8所述的CT重建算法,其特征在于:所述CT重建算法包括反投影滤波(BPF)算法,所述反投影滤波(BPF)算法包括以下步骤:
对探测器上获得的物体图像进行反log操作:以得到处理后的投影图像;
对每个角度下的锥束投影数据进行求导;
由求导后的投影数据,对螺旋轨道内每条PI线上的物体图像点进行加权反投影;
对加权反投影得到的PI线上的数据,沿着该条PI线求Hilbert变换,得到物体在该线上的重建结果;以及
对螺旋轨道内所有PI线上的物体重建结果进行重采样,得到物体的三维重建图像。
10.根据权利要求8所述的CT重建算法,其特征在于:所述CT重建算法包括滤波反投影(FBP)算法,所述滤波反投影(FBP)算法包括以下步骤:
对探测器上获得的物体图像进行反log操作:以得到处理后的投影图像;以及
对所述投影图像进行加权、卷积滤波以及加权反投影,以获得CT图像。
11.一种CT系统几何校准方法,其特征在于包括根据权利要求1所述的用于CT图像重建的修正方法。
12.根据权利要求11所述的几何校准方法,其特征在于所述几何校准方法选自:基于小球体模投影椭圆方程的CT系统几何校准方法、以及基于投影矩阵的CT系统几何校准方法。
13.一种用于CT图像重建的修正系统,其特征在于至少包括:
固定坐标系定义模块,其用于定义固定坐标系XYZ,并获取旋转坐标系中的单位矢量以及固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的矢量的表示;
修正模块,其用于获取修正后的旋转坐标系中的单位矢量以及修正后的坐标原点O与探测器上的任意点所构成的矢量的表示;
新坐标获取模块,其用于根据修正后的所述矢量的表示得到探测器上的任意点在所述固定坐标系XYZ中的新坐标。
14.根据权利要求13所述的系统,其特征在于,所述固定坐标系定义模块定义固定坐标系XYZ,其中Z为旋转轴,U轴和V轴为探测器上像素的行和列方向,分别与坐标轴Y和Z平行,在固定坐标系XYZ中,旋转坐标系中的所述单位矢量表示为:
并且
所述固定坐标系XYZ的坐标原点O与探测器上的任意点所构成的所述矢量表示为:
其中,R为光源至旋转轴Z的距离,D为光源至探测器的距离。
15.根据权利要求13所述的系统,其特征在于:所述修正模块根据探测器偏置角度φ获得修正后的单位矢量以及矢量探测器偏置角度φ为探测器围绕U轴旋转角度φ,在探测器围绕U轴旋转角度φ后,旋转坐标系的所述单位矢量变更为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>e</mi>
<mi>v</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>e</mi>
<mi>v</mi>
</msub>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>e</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
2
所述矢量表示为:
16.根据权利要求15所述的系统,所述新坐标获取模块得到探测器上的任意点在所述固定坐标系XYZ中的新坐标为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mo>-</mo>
<mi>u</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mi>R</mi>
<mo>+</mo>
<mi>v</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>u</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<mi>R</mi>
<mo>+</mo>
<mi>v</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>z</mi>
<mo>=</mo>
<mi>v</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&phi;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
3
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710373329.6A CN107233105B (zh) | 2017-05-24 | 2017-05-24 | 一种用于ct图像重建的修正方法及修正系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710373329.6A CN107233105B (zh) | 2017-05-24 | 2017-05-24 | 一种用于ct图像重建的修正方法及修正系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107233105A true CN107233105A (zh) | 2017-10-10 |
CN107233105B CN107233105B (zh) | 2020-12-11 |
Family
ID=59985195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710373329.6A Active CN107233105B (zh) | 2017-05-24 | 2017-05-24 | 一种用于ct图像重建的修正方法及修正系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107233105B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109671128A (zh) * | 2018-12-07 | 2019-04-23 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
CN110517330A (zh) * | 2019-08-07 | 2019-11-29 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN114708344A (zh) * | 2022-03-14 | 2022-07-05 | 北京理工大学 | 一种基于场论的ct图像局部重建的系统及方法 |
CN114758079A (zh) * | 2022-05-19 | 2022-07-15 | 国网河南省电力公司电力科学研究院 | 电气设备内部结构局部成像三维重建方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120051626A1 (en) * | 2010-08-30 | 2012-03-01 | Yong Long | Three-dimensional forward and back projection methods |
CN102973291A (zh) * | 2012-12-20 | 2013-03-20 | 电子科技大学 | C型臂半精确滤波反投影断层成像方法 |
CN103006251A (zh) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
CN103445803A (zh) * | 2013-09-09 | 2013-12-18 | 深圳先进技术研究院 | 基于正弦图的ct系统射束硬化消除方法及其系统 |
CN105849537A (zh) * | 2013-11-28 | 2016-08-10 | 尼康计量公众有限公司 | 校准设备和计算机断层扫描方法 |
CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
-
2017
- 2017-05-24 CN CN201710373329.6A patent/CN107233105B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120051626A1 (en) * | 2010-08-30 | 2012-03-01 | Yong Long | Three-dimensional forward and back projection methods |
CN103006251A (zh) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
CN102973291A (zh) * | 2012-12-20 | 2013-03-20 | 电子科技大学 | C型臂半精确滤波反投影断层成像方法 |
CN103445803A (zh) * | 2013-09-09 | 2013-12-18 | 深圳先进技术研究院 | 基于正弦图的ct系统射束硬化消除方法及其系统 |
CN105849537A (zh) * | 2013-11-28 | 2016-08-10 | 尼康计量公众有限公司 | 校准设备和计算机断层扫描方法 |
CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109671128A (zh) * | 2018-12-07 | 2019-04-23 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
CN110517330A (zh) * | 2019-08-07 | 2019-11-29 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN110517330B (zh) * | 2019-08-07 | 2021-05-28 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN114708344A (zh) * | 2022-03-14 | 2022-07-05 | 北京理工大学 | 一种基于场论的ct图像局部重建的系统及方法 |
CN114708344B (zh) * | 2022-03-14 | 2024-06-25 | 北京理工大学 | 一种基于场论的ct图像局部重建的系统及方法 |
CN114758079A (zh) * | 2022-05-19 | 2022-07-15 | 国网河南省电力公司电力科学研究院 | 电气设备内部结构局部成像三维重建方法 |
CN114758079B (zh) * | 2022-05-19 | 2024-02-20 | 国网河南省电力公司电力科学研究院 | 电气设备内部结构局部成像三维重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107233105B (zh) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107233105A (zh) | 一种用于ct图像重建的修正方法及修正系统 | |
CN108122203A (zh) | 一种几何参数的校正方法、装置、设备及系统 | |
CN101133962B (zh) | 再现三维立体图像的方法和x射线设备 | |
JP4415762B2 (ja) | 断層撮影装置 | |
JP2008520257A5 (zh) | ||
JP4360817B2 (ja) | 放射線断層撮影装置 | |
CN107592935A (zh) | X射线ct的自动系统校准方法 | |
CN105769233A (zh) | 一种几何校正方法 | |
US7978810B2 (en) | Imaging method for variable pitch spiral CT and a CT machine for carrying out the method | |
JP7154611B2 (ja) | インテリアct画像生成方法 | |
CN104665862B (zh) | 在cbct中消除几何伪影的方法以及使用该方法的cbct系统 | |
US7737972B2 (en) | Systems and methods for digital volumetric laminar tomography | |
JP2019158534A (ja) | 計測用x線ct装置、及び、その断層画像生成方法 | |
CN103717137B (zh) | 图像生成装置 | |
Herbst et al. | Dynamic detector offsets for field of view extension in C‐arm computed tomography with application to weight‐bearing imaging | |
US8861829B2 (en) | Method and system for reconstruction of tomographic images | |
CN103620393B (zh) | 成像装置 | |
JP4498023B2 (ja) | X線ct装置 | |
JP2007198866A (ja) | 広義サドルコーンビームct装置および3次元再構成法 | |
CN106097411B (zh) | Ct机图像重建方法及高分辨ct扫描机 | |
US8577113B2 (en) | Image reconstruction for truncated detector array | |
Sun et al. | Estimation of local data-insufficiency in motion-corrected helical CT | |
JPWO2007096936A1 (ja) | 断層撮影装置および演算処理プログラム | |
JP2010184086A (ja) | 放射線ct装置及びその制御方法 | |
Xia et al. | Image reconstruction in reduced circular sinusoidal cone-beam CT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |