CN1062610A - 将锥形束x射线投影数据再现物体三维ct图象用的方法和装置 - Google Patents
将锥形束x射线投影数据再现物体三维ct图象用的方法和装置 Download PDFInfo
- Publication number
- CN1062610A CN1062610A CN91111691A CN91111691A CN1062610A CN 1062610 A CN1062610 A CN 1062610A CN 91111691 A CN91111691 A CN 91111691A CN 91111691 A CN91111691 A CN 91111691A CN 1062610 A CN1062610 A CN 1062610A
- Authority
- CN
- China
- Prior art keywords
- ijk
- plane
- straight line
- line
- circle
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/416—Exact reconstruction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Nuclear Medicine (AREA)
- Image Generation (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了经逆Radon变换再现三维CT图
象用的、将X射线锥形束数据(通过物体的线积分)
变换成Radon数据(平分积分)的方法。各平面积分
的径向导数,是经沿归一化检测器平面上的一对直线
中的每一条进行积分以求出加权的线积分来求的,这
对直线被定义为共用一条旋转轴且彼此相对转动δβ
角的一对相应积分平面与归一化检测器平面的交线,
再将此加权线积分间的差值除以转角δβ。此法可将
锥形束数据变换为表示Radon空间中在任何一组平
面上的平面积分值,例如一组共轴垂直平面。
Description
在此公开及请求保护的本发明,涉及的是以下共同指定的专利申请主题,因而这些专利申请全部的公开,在此合并引作参考:
并此同时由Kwok C.Tam提交的申请号为-,发明名称为“由锥形束投影数据或由平面积分再现物体三维CT图象用的平行处理方法和装置”的专利申请〔RD-19564〕,以及
并此同时由Kwok C.Tam提出的申请号为-,发明名称为“由不完全锥形束投影数据再现物体三维CT图象用的方法和装置”〔RD-19695〕的专利申请。
本发明总地涉及三维计算机X射线断层摄影术(CT),更确切地说,涉及将X射线锥形束数据变换为平面积分,通过逆氡变换(inverse Radon transformation)再现三维图象用的方法和装置。
在传统的既为医学且为工业应用的CT中,采用的是扇形束的X射线和线性列阵的检测器。所获得的是二维成象。当数据组为完整的,而且象质也相当高时,只有物体的单一切片(Slice)能被同时成象。在需要三维图象时,采用的是“切片堆叠”的方法。同时取得三维数据组的二维切片,本来就是缓慢而费时的。此外,在医学应用中,由于相邻的切片无法同时成象,所以会产生移动的赝象。而且,由于切片之间的距离通常都比X射线准直器的孔径小,所以能利用的辐射剂量要比最佳剂量低,其结果是使身体的许多部位要被双重照射。
比较更新一些的基于所谓的锥形束几何结构基础上的方法,所采用的是二维列阵的检测器来替代线性列阵的检测器,以及锥形束的X射线源来替代扇形束的X射线源。整个物体随时都被锥形束的X射线源照射,因而锥形束扫描要比使用扇形束或平行束进行一个切片一个切片地扫描要快得多。而且,由于物体上的每个“点”是靠X射线以三维而不是二维的方式进行观察的,故所达到的对比度要比使用传统的二维X射线CT可能达到的高得多。为了获得锥形束投影数据,最好对物体在360°的角度范围内进行扫描,或者通过让X射线源在适当的扫描轨道(例如围绕该物体的圆形轨道上)上运动,同时保持此二维列阵检测器参照此X射线源固定不动;或者通过转动该物体,同时让X射线源和检测器保持不动。在两种情况下,都是靠X射线源和物体之间的相对运动来实现扫描的。
大多数X射线CT中的图象再现过程,都是基于氡的逆变换(inversion)过程基础上的,其中的物体图象是由物体的氡变换总额再现的。二维物体的氡变换,是由贯穿该物体的各条线上物体密度的积分构成的。三维物体的氡变换,是由平面积分构成的。然而锥形束数据是不能和通过逆氡变换再现图象直接兼容的,它需要使用物体的平面积分作为输入。所以,通过逆变换由锥形束扫描数据再现出图象,一般要包括两个步骤,即(1)将锥形束数据变换为平面积分,以及(2)在平面积分上完成逆氡变换,以获得图象。本发明主要致力于用来将X射线锥形束数据变换成对于氡空间内的一组适宜平面的平面积分或代表平面积分值的有效方法和装置。上述联案申请-〔RD-19564〕公开了用于完成逆氡变换的两步法,首先从对于氡空间中的一组共轴的垂直平面的平面积分开始。因而,在此公开的本发明,可被利用来将X射线锥形束数据变换成能代表对于氡空间中的一组共轴的垂直平面的平面积分的值,而且申请号为-〔RD-19564〕的发明,可被利用来完成三维图象再现的逆氡变换部分。
一种用来将锥形束数据变换成平面积分的方法,已被公开在Cerald N.Minerbo的“由锥形束投影数据进行卷积再现”一文中,参见IEEE Trans.Nucl.Sci.,NS-26卷,No.2,第2682~2684页(1979.4)。然而遗憾的是,如在L.A.Feldkamp,L.C.Davis以及J.W.Kress的“实际锥形束算法”一文中(参见J.Opt.Soc.Am.A.,1卷,No.6,第612-619页,1984,6)讨论过的那样,Minerbo的求导包含有不容易矫正的误差,而且提供的结果不能成立。
在Bruce D.Smith的“由锥形束投影进行图象再现:必要和充分的条件以及再现方法”一文中(参见IEEE Trans.Med.Imag.,MI-44卷,第1425页,1985,3),公开了一种以锥角形(Horn)影响函数由锥形束数据变换成平面积分的一维卷积的方法。由于该卷积是把所有平面上的平面积分混合在一起,所以对卷积结果中一点的计算,对应在一个视角度下,需要有检测器上的所有数据。因此,计算方面的任务是非常紧张的。
在P.Grangeat的“通过从锥形几何结构的X射线摄影中再现三维成象系统的分析”(“Analyse d′un System D-Imagerie 3D par Reconstruction a partir de Radiographiex X en Ceometrie conique”)一文中,参见其在国际电讯学院的博士论文(I-Ecole Natinale Superieure des Telecommunications),法国(1987),公开了用来从锥形束数据计算平面积分的导数方面的技术。然而计算出来的数据点存在于氡空间内球壳的一组大圆上。这些大圆通常并不落在氡空间中任何随意的一组平面上,而且也不落在氡空间中一组共轴的垂直平面上。因而它们不适于作为逆氡变换的输入。为取得有关被用在逆氡变换中的垂直平面的数据,还要求在三维内插方面作广泛的工作,而且内插会把误差引入数据中。
因此,本发明的目的在于提供一种方法和装置,用来将X射线锥形束数据变换为代表对于氡空间中任何随意一组平面的平面积分值,以便通过逆氡变换再现出三维图象。
本发明更加明确的目的,是提供一种方法和装置,用来将X射线锥形束数据变换为能代表对于氡空间中一组共轴的平面的平面积分值,以便通过逆氡变换再现出三维图象。
本发明的另一个目的是提供一种方法和装置,对于每一个视角,与已有技术中的方法要求检测器的全部数据相比,在检测器上仅有两条小间隔的数据线用来计算代表平面积分的值。
本发明的再一个目的是提供一种方法和装置,它是精确的而且并不需要内插,能够将X射线锥形束数据变换成能代表氡空间中对于一组共轴平面的平面积分值,或者变换成能代表氡空间中对于任何随意一组平面的平面积分值,以便通过逆氡变换再现出三维图象。
本发明最后一个目的是提供一种方法和装置,在将X射线锥形束数据变换成能代表氡空间中对于一组共轴的平面或者任何随意的一组平面的平面积分值,以通过逆氡变换再现出三维图象时,它能将所需要的计算量减至最少。
根据本发明提供的由锥形束投影数据再现出物体三维图象的方法,其中的锥形束投影数据,是以穿过物体的直线积分形式,对于多个X射线源位置Si中的每一个位置编制成归一化检测器平面上的二维数据组,此归一化检测器平面则包含有原点并垂直于由每一特定的源位置Si至该原点的直线。此方法包括两个总体步骤,即先确定能代表氡空间中对于一组平面φj的平面积分的值,然后对这些能代表对这组平面φj的平面积分的值完成逆氡变换,以再现出物体的图象。在更具体的实施例中,此平面φj包括含有与原点相交的参考轴的一组共轴的平面。
本发明的一个重要方面是要求出氡空间中在特定点处的平面积分值或者氡数据(实际上是氡数据的径向导数),即通过沿归一化检测器平面上的一对直线中的每一条进行积分以求出加权的线积分,这对直线被定义为共用一根旋转轴并且彼此之间相对转动δβ角的一对相应的积分平面与此归一化检测器平面的交线,然后将此加权的线积分之间的差值除以该转角δβ。
更确切地说,确定代表对氡空间中一组平面φj上的平面积分值的步骤,包括对于每一种X射线源位置Si的如下嵌套(nested)步骤:
在氡空间中确定相应的球壳,球壳上的氡数据可被求出,确定平面φj与对应于特定X射线源位置Si的球壳的交线,并在该球壳上确定一组圆Dij,而且,对于这组圆Dij中的每个圆来说:
将旋转轴定义为,经过此特定的X射线源位置Si而与此特定圆Dij相交,且与此特定圆Dij的平面正交的一条直线;
确定一组共轴的积分平面Qijk,此积分平面Qijk中的每一个都包含此特定的旋转轴,并且与此特定的圆Dij相交以确定氡数据点Rijk的位置,而且此积分平面Qijk在相应的线Lijk上与归一化检测器平面相交,而且,对于此归一化检测器平面上的每条线Lijk来说:
将相应的积分平面Qijk转动一个小角度δβ,确定在相应的直线Lijk′与此归一化检测器平面相交的平面Qijk′;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,并且
将此加权的线积分Jijk和Jijk′之差除以旋转角δβ,以得出在特定点Rijk处氡数据的径向导数。
与此类似,根据本发明的用于从锥形束投影数据再现物体三维图象的装置,包括如可编程序的计算机装置,用来利用上面概括出来的方法,求出能代表在氡空间中一组平面φj的平面积分的数值的装置;以及如在申请号为-〔RD-19564〕中公开的处理装置,用来在表示对于一组平面φj的平面积分的数值上完成逆氡变换,以便再现物体的图象。
在更具体的实施例中,平面φj包括含有与原点相交的参考轴的一组共轴的平面,确定代表对于这组平面φj的平面积分值的步骤,最好包括对于不在参考轴上的每一个源位置Si的嵌套步骤:
在氡空间中对应的球壳上,在包含此特定的源位置Si并垂直于平面φj的平面内,确定出相应的圆Gi,确定此平面φj和圆Dij与该特定圆Gi的交点,并在此圆Gi上确定与圆Dij对应的多个点Pij;
将相应的圆Gi,由特定的源位置Si在为一化的检测器平面上投影为一条直线Mi;将Pij点投影至直线Mi上的对应点Cij,并且,对于归一化检测器平面内的每个投影点Cij:
在此归一化检测平面内经过此被投影的点以多种取向作出一些线Lijk,这些线Lijk,也就是沿着一条通过此特定的X射线源位置Si、特定的点Pij以及特定的被投影点Cij的直线的相应积分平面Qijk(其中每个积分平面Qijk都包含旋转轴)在归一化检测平面上的交线;
在此归一化检测平面内,将这些线Lijk中的每一条绕着该投影点Cij转动一个小角度δθ,以确定线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测平面的交线,并且可根据几何关系从角度δθ求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,并且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
在更加具体的实施例中,平面φj包括含有与原点相交的参考轴的一组共轴的平面,确定代表对于这组平面φj的平面积分值的步骤,包括对于参考轴上的每一种X射线源位置Si的嵌套步骤:
对于每一个同对应于此特定X射线源位置Si的球壳相交的平面φj,确定一个特定的圆Dij;
将此特定的圆Dij从特定的X射线源位置Si投影到归一化检测器平面内的线Lij*上;
在归一化检测平面内作一些与此直线Lij*垂直的平行线Lijk,这些平行线Lijk,也就是沿着一条通过此特定的X射线源位置Si并正交于此特定圆Dij平面的直线的相应积分平面Qijk(其中每个积分平面Qijk都包含旋转轴)在归一化检测器平面上的交线;
将上述平行线Lijk中的每一条平移一小段距离,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并且可根据几何关系从线Lijk与Lijk′之间的距离求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应加权的线积分Jijk和Jijk′,并且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij相交处氡数据的径向导数。
本发明的新特点在所附的权利要求书中具体提出,本发明的构成方式和内容,从以下结合附图所作的详细描述中将会更好地熟悉和理解,其中
图1描绘用于三维CT的锥形束扫描几何图形与实施本发明的再现装置连接的情形;
图2a、2b、2c、2d、2e和2f为描绘氡变换法进行三维CT成象的图解;
图3为在给定点处物体的三维氡变换表示;
图4表示氡空间内的一组共轴的平面φj,其中每个平面都包含一条垂直轴或参考轴,在这些轴上可求出氡数据(平面积分);
图5A为物体的参照系和坐标系;
图5B为归一化检测器参照系和坐标系;
图5C表示图5A和5B的坐标系原点相重合的方式;
图6表示积分平面、积分参照系和对应的坐标系;
图7表示与经过物体的积分平面对应的锥形束数据;
图8与其类似,表示与经过物体的一对相隔很小的相邻积分平面对应的锥形束数据;
图9表示积分平面上的几何关系;
图10表示根据本发明的由X射线锥形束数据计算氡数据的径向导数的方法;
图11表示可从一个X射线源位置上计算出来的平面上的氡数据;
图12描绘一个能代表从一个源位置上计算出来的所有氡数据的球壳或氡球壳;
图13表示如情况之1实例称作绕a轴旋转的操作;
图14表示通过绕a轴旋转产生的氡数据;
图15表示在投影至归一化检测器平面的一条直线上的每一点,通过完成旋转操作产生氡数据的情况;
图16表示如情况之2实例的绕b轴的旋转操作;
图17A表示通过绕b轴旋转产生氡数据的情况;
图17B表示由绕b轴旋转而在归一化检测器平面内产生数据点及数据线的情况;
图18A用线标示图4,表示氡空间内一组共轴的垂直平面;
图18B表示对于特定的锥形束X射线源位置产生的氡球壳;
图18C表示图18A的共轴垂直平面与图18B的氡球壳相交;
图19详细表示如图18C中垂直平面与氡球壳相交处的许多圆中之一,并且说明在垂直平面内由检测数据产生氡数据的情况之1的方法;
图20表示当源位置S处在垂直平面的转轴上时,垂直平面与氡球壳相交处的一些圆,并且说明由检测器数据产生氡数据的情况之2的方法。
图21表示两垂直平面间的相交。
首先参见图1,这是一种典型的采用锥形束几何构形的扫描及获得数据结构,与实施本发明的再现装置连接在一起。物体20位于锥形束X射线点状辐射源22和二维检测器列阵24之间的视场中,能够提供出锥形束的投影数据。旋转轴26则穿过该视场和物体20。中平面28可被定义为包含此X射线点源22并垂直于旋转轴26的平面。按照惯例,旋转轴26被称之为Z轴,旋转轴26与中平面28的交点取作坐标原点。如图所示,X和Y轴位于中平面28内。为了从多种X射线源位置Si上对物体20进行扫描,X射线源22应沿着适当的扫描轨道30相对此物体20及视场运动,同时让检测器24相对此X射线源22保持固定。在图1中,为解释方便起见,扫描轨道30被表示为位于中平面28内的圆形扫描轨道30,然而其它的扫描轨道也可以采用,而且正如在下文中概括讨论的,事实上是更可取的。
检测器列阵24连在数据获得系统(DAS)32上。在工作时,穿透物体的X射线光子被X射线检测器列阵24检测,并被数据获得系统(DAS)32记录。通过光子计数,在由空气信号归一化并变换为负对数(negative of the logarithms)之后,就能表示经过物体20的线积分。因此,这些数据是靠X射线源22和检测器24沿扫描轨道30进行扫描(或者等效于让该物体20旋转,而X射线源22和检测器24保持固定),在绕物体20的许多X射线源位置Si上获得的。
然而应当指出的是,按照这种单一扫描采集到的数据组(集)是不完整的,而且可能容许也可能不能容许的赝象会因而导入再现中,这取决于具体的应用。Smith在上述1985年的文章中曾经指出,如果在穿过研究中的物体的每一平面上都存在来自该X射线源扫描轨道的点(假定检测器相对于X射线源锁定,并且大到足以跨越被检查的物体),那么锥形束数据组就是完整的。由Minerbo在上述1979年的文章中以及由Heang K.Tuy在“用于锥形束再现的反演公式”(参见SIAM J.Math.,43卷,No.3,第546-552页,1983,6)的文章中提出的构形(Smith指出能满足其数据完整性条件的),采用的是彼此垂直的两个圆形的X射线源扫描轨道。另外一种能获得数据完整性的扫描构形,被公开在共同指定的美国专利申请No.07/572,651中,是1990.8.27.由Eberhard等人提交的,名称为“三锥CT中为数据完整性用的方波锥形束扫描轨道”。另一个办法是,替代这种获得完整的锥形束X射线数据组,可以采用上面联案申请的申请号为-〔RD-19695〕的发明,采用光学方式获取物体的边界信息,以对逆氡变换过程中失去的数据进行迭代校正。
数据获得系统32是接在具有代表性的处理机34上的,后者被用来再现物体20的三维图象,即根据本发明从通过该物体20的线积分中,计算出对于一组平面的平面积分,并在此平面积分上完成逆氡变换,以便再现出物体20的三维图象。适合的图象显示器36接在该具有代表性的处理机34上面,用作输出装置。
现在参见图2a至2f以及图3,概括表示的是氡变换进行三维成象的方法。
具体说来,物体本身是以其X射线衰减系数f(x,y,z)定义的(图2a)。测得的锥形束投影数据,于是就对应于该函数在辐射方向上的线积分X(θ)=∫f(r,θ,z0)dr(图2b)。检测器数据的线性积分(又称为检测器积分),是由∫X(θ)dθ=∫∫f(r,θ,z0)dr dθ给出来的(图2c)。在平行束的情况下,这些检测器积分简单地等于物体的氡变换。然而在锥形束的情况下,氡变换代之以下式∫∫f(r,θ,z0)r dr dθ给出(图2d)。氡变换积分中的附加因子r,是由于从笛卡尔坐标变换为极坐标的坐标变换雅可比行列式产生的。如在图2e及2f中表示的那样,逆氡变换过程能从检测器积分中再现出三维CT图象。由于直接的逆氡变换要求以物体的平面积分作为输入,所以将线积分(锥形束检测器积分)变换成平面积分(氡数据)的预备步骤便成为必须,本发明正是指向于此。
如在图3中表示的那样,在点X0,Y0,Z0处物体的三维氡变换,是由X射线衰减系数在通过该点X0,Y0,Z0的平面(即垂直于原点至该点X0,Y0,Z0的直线的平面)上的面积分给出的,并可表示为
R(X0,Y0,Z0)=∫∫f(x,y,z)da (1)
平面
对于二维氡变换来说,除了积分是对于线而不是对于面之外情况是类似的。
此平面积分还可表示为
靠逆氡变换,可从其平面积分R再现出三维物体f(r);逆氡变换可被表示为
如在上述联案申请的申请号为-〔RD-19564〕的申请中详细描述的那样,方程式(3)中表示的逆氡变换,可以通过两步过程来实现。步骤1,包括在氡空间中每一个都包含Z轴的许多垂直平面φj内的二维CT图象再现。步骤2,包括在许多水平平面内的二维CT图象再现。
因此,如在图4中表示的那样,所谓需要作为逆氡变换输入的,是在氡空间中包含参考轴的许多平面φj上求出和构成的平面积分,例如在包含垂直参考轴或Z轴48的垂直平面40,42,44和46上的平面积分。
在此详细描述的方法和装置,用于将X射线锥形束数据,变换为在一组共轴的垂直平面φj上或在任何随意的一组平面上的平面积分,以作为逆氡变换输入。
参见图5A、5B及5C,我们首先规定在此进行分析中使用的参照系及其相关联的坐标系。(应当指出的是,在详细讨论的下文中,参照系及坐标系不同于以图1的典型扫描构形表示的广义坐标系。)具体说来,图5A表示物体的参照系,是相对该物体20固定的。在此参照系中随意一点的空间位置,可被表示为(X,Y,Z)三个一组。方程式(2)和(3)中变量,就是物体参照系中的变量。图5B表示归一化检测器的参照系,是相对于此归一化检测器50固定的。在归一化检测器的参照系中的空间位置,是由(U,V,W)三个一组表示的。(U,V,W)坐标系的原点,位于此归一化检测器50的中心,而且U轴和V轴位于此归一化检测器平面内。X射线源的位置Si总是在W轴上,但其距该归一化检测器50中心的距离,可从一个X射线源位置Si变化到另一个X射线源位置Si。
如在图5C中表示的那样,我们假定(X,Y,Z)坐标系的原点总是和(U,V,W)坐标系的原点重合的。实际上,这等于仅仅是把实际的检测器24读数换算到平面上,该平面通过(X,Y,Z)坐标系的原点,并且与连接X射线源位置Si和原点的直线正交。
为了简化平面积分的变换,我们现在引入第三个参照系。参见图6,Q是随意一个包含X射线源S的平面。令平面Q在线L上与此归一化检测器平面相交。将坐标系(a,b,c)定义为以X射线源S作为原点,以使 为在平面Q中与线L正交的单位向量, 为在平面Q中与L平行的单位向量,而且 = × 。我们将把此坐标系(a,b,c)称为积分参照系中的坐标系。为了进一步简化平面Q上的积分,我们指出,因其坐标C总为0,故平面Q上的每个点均以一对(a,b)表征。通过进行坐标变换:(a,b,c)→(r,θ,c),可以把一对(a,b)变换成与 , 座标系相关的极坐标(r,θ),其中
θ=tan-1 (b)/(a)
图7表示锥形束扫描中的典型情况。假定一个如图7所示的物体的平面Q或者切片。此锥形束X射线源S以直线L的形式将此平面Q投影到归一化检测器平面上。在此积分参照系中,令a轴在C点与直线L相交(图6)。通过作图,SC应与线L正交。令|SC|代表X射线源S和交点C之间的距离。线L上的数据X(t),其中t表示点C沿L的位移,可由下式给出:
X(t)=∫∫f(r,θ,O)δ〔θ-tan-1( (t)/(|SC|) )〕drdθ
=∫f(r,tan-1(t)/(|SC|) ,O)dr
换言之,数据X(t),表示沿Q平面上与a轴成θ=tan-1(t/|SC|)角直线的物体密度的线积分。并且指明变量t与tanθ成正比,通过对直线L上的变量t进行X(t)积分(以某种适合的加权),人们有希望得到f(r,θ,O)对于平面Q上的变量r和θ的积分值,即∫∫f(r,θ,O)drdθ。为此目的,让我们用变量t来表示量I=∫∫f(r,θ,O)drdθ。于是
t=|SC|tanθ
dt=|SC|sec2θdθ
因而积分I可由下式给出:
因此,通过对直线L上的锥形束数据X(t)进行加权积分,便可以得到量I=∫∫f(r,θ,O)drdθ。与此对比,作为该平面的氡数据R(s, ),在物体参照系中,其中的S为平面Q至原点的距离,而且 为法线,则可由下式给出:
由于I≠R(s, ),所以氡数据不能靠沿着归一化检测器平面上的直线对锥形束数据积分来获得。(然而,如果检测器上的数据是由X射线平行束产生的,则沿此归一化检测器平面上的直线对数据进行积分,会能产生氡数据。)
氡数据R(s, )和积分I之间仅有的差别是,积分I中没有r因子。由于在每一r值处的密度值f(r,θ,O)是不知道的(换句话说,则不需要在第一位置进行锥形束扫描),所以这种差别不能够通过以r对此数据加权来补偿。
注意到一种将r因子引进被积分函数中的方法,让平面Q围绕经过原点的平面上的任一条直线旋转,则此平面上的每个点都被平移一个与rsinγ成正比的量,其中r为该点的径向坐标,γ为该点相对此转轴的角度坐标。这种观察促进了参照图8及9作出的下列进展。
现在参见图8,让我们来假定物体中另一个与图7的平面Q非常接近的平面Q′。平面Q′,是靠将平面Q本身,围绕经过X射线源位置S的平面Q上的转轴a′转动一小角度δβ得到的。对于实际转角δβ大小的选择,需要综合考虑每一具体系统中的精度和信噪比之间的关系。比较小的转角δβ会得到较高的精度,但以增加数据中的噪声为代价,且反之亦然。平面Q′在归一化检测平面内投影出另一条线L′。线L和L′相交于C′,此处就是转轴a′与法向检测器平面的相交处。
令δθ表示L和L′间的夹角。(显然,角度δθ和δβ是以简明的几何关系彼此相关的,而且δβ可从δθ求出。下文对于两种具体的情况将提供该计算的实例。)平面Q′上的每个点(r,θ,O),在极坐标中是和积分参照系相关的,可被认为是在平面Q中由相应点平移一个δr量得出的。用α表示转轴a′和a之间的角度,则在点r=(r,θ,O)处的平移量δr,可由下式给出:
改变δr的空间坐标r,将导致密度值f(r)的相应改变,随后再引起积分值I的δI改变,并可由下式给出:
所以
其中
T(c)=∫∫〔f(r,θ,c)sin(θ-α)〕rdrdθ
T(c)这个量除了附加因子sin(θ-α)之外几乎与方程式(5)的氡数据完全一样。然而,现在附加因子可被补偿,因为它只包含角度变量,其值可通过测量来改变。这可以通过规定一个新的积分J来实现,该积分包括在I的被积分函数内的加权因子,以消掉此附加因子,
J=∫∫ (f(r,θ,O)drd0)/(sin(θ-α)) (6)
平面Q内的几何关系表示在图9中。由图9中我们得出
sinα= (△C)/(|SC′|)
cosα= (|SC|)/(|SC′|)
其中△C表示C′和C的位移。因此我们有
sin(θ-α)=sinθcosα-cosθsinα
(7)
于是,我们可以将所要求的积分J表示为包含有加权因子1/sin(θ-α)的变量t。将方程式(4)和(7)代入方程式(6),我们得到
(8)
通过如前面I积分同样的数学处理,我们得到
(dJ)/(dβ) = ∫∫f(r,θ,c)C=0rdrdθ (9)
其中,在物体参照系中,S为平面Q和原点之间的距离,并且
就是函数f在平面Q上的平面积分。
利用方程式(9),原则上我们能从锥形束数据中计算出氡数据的径向导数,而且氡数据本身也可以通过对在径向尺度(dimension)的结果进行积分来获得。此过程表示在图10中。为了求出在点P=s 处氡数据的径向导数值,在物体参照系中,此处s=|OP|,而且 =OP/|OP|,我们进行如下:
1.确定经过点P并垂直于直线OP的平面Q;
2.确定平面Q与归一化检测器平面相交处的直线L;
3.在直线L上找出一点C,使直线SC垂直于L;
4.在直线L上任取一点C′,将旋转轴a′规定为从S到C′的直线。
等效地,让平面Q绕此旋转轴a′转过一个小角度δβ,于是得到平面Q′;让直线L在此检测器平面内绕点C′转过一个小角度δθ,于是得到直线L′,平面Q′则在直线L′处与此归一化检测器平面相交;
5.利用方程式(8),分别计算出在直线L和L′上的量J和J′;
6.根据几何关系由δθ计算出角度δβ;
7.利用下列方程式,从量J,J′和δβ得到点P处氡数据的径向导数:
= (J′-J)/(δβ)
采用上述程序,我们能够得到对于穿过为锥形束X射源辐照的物体的所有平面的氡数据。顺便要说的是,这和前面Smith 1989年文章中叙述的条件是一致的,在锥形束扫描中,为了具有完整的数据,穿过物体的每个平面都应与X射线源的位置相交。
现在参见图11,12及13,在一种X射线源位置上能够产生的氡数据(即平面积分)的范围,可以用这种方法定量地估算出来。令图11的平面为包含锥形束X射线源S和原点O的任意平面,称之为W平面。假定有任一个平面Q与平面W正交且包含此X射线源S。令U为平面Q与平面W相交的直线,也就是说平面Q经过直线U且与平面W即图11的图面正交。让V为平面W上经过原点且与直线U正交的直线,并让直线U和V在点P相交。让 为沿直线V的单位向量。如在附录A中参考图21表示的那样,向量|OP| 垂直于平面Q,因而在Q上的平面积分就是在物体参照系中的氡数据R(|OP|, ),即在P点的氡数据。由于角OPS为直角,所以点P位于平面W内以OS为直径的圆上。通过把同样的操作加在所有垂直于平面W且经过X射线源S的平面上,于是就在整个圆上产生出氡数据,如在图11中那样。
如在图12中表示的那样,通过重复W平面内在包含线段OS的所有其它平面上进行的整个操作,氡数据便在包含OS作为直径的所有圆上产生出来。换而言之,氡数据是在以OS作为直径的球壳上产生出来的。这种球壳又可称之为氡球壳。
现在将仔细考虑两种具体的情况,其区别在于图6中的积分平面Q围绕其转动的旋转轴的取向不同。如在上文中描述的那样,许多旋转轴中的每一条都通过此X射线源位置S,通过此X射线源位置S及原点O的直线SO则与归一化检测器平面正交,而且此归一化检测器平面中包含有原点。
情况之1适用于,当旋转轴不与b轴重合并且包括此旋转轴与a轴重合的特定情况。因而情况之1可被描述为绕a轴的连续旋转。在此描述的具体实施例中,氡数据是在一组共轴的平面φj上产生的,其中每个平面都包含如图4中表示的垂直轴或参考轴,情况之1则适用于不在此垂直轴或参考轴上的所有X射线源位置Si。
情况之2适用于当旋转轴与图6中b轴重合的情况。在这种情况下,旋转轴是和归一化检测器平面平行的。在此描述的具体实施例中,氡数据是在一组共轴的平面φj上产生的,其中每个平面都包含如图4中表示的垂直轴或参考轴,情况之2则适用于在此垂直轴或参考轴上的所有X射线源位置Si。
对于情况之1的旋转轴与a轴重合的特定情况,将不予以考虑。在这种情况下图10中的α=0,而且两条线L和L′相交在a轴与检测器平面相交的C点。如果在归一化检测器平面上过C点以所有取向作出一些直线(图13),那么这些线就是包含线SC作为共用轴线的一些积分平面在此归一化检测器平面上以不同取向的投影。把这些积分平面称为Q1,Q2,Q3……等等。那么根据方程式(9),就能由量dJ/dβ(顾及检测器平面上对于每一对相邻的具有小间隔的直线计算出来的J中的加权函数sinθ)给出投影在这一对相邻直线上的平面Qi上的平面积分的导数。
这种情况表示在图14中,它描绘出在SC方向上的视图,即由X射线源朝向检测器上面的交点C。图14的平面包含有原点O,而且线SC与图14垂直。P点为图14的平面和线SC之间的交点。由于线SC垂直于图14的平面,因此包含线SC的所有平面Qi,都作为构成以P点为中心的极坐标的直线出现。图14中的这些称之为A1,A2,A3,……等等的直线,分别和平面Q1,Q2,Q3,……等等对应。由原点O对这些直线中的每一条线作垂线,并让每一对垂线如图所示在位置B1,B2,B3,……处相交。而且,如在附录A中表明的那样,由原点垂直于线A1,A2,A3,……的直线,同样也垂直于平面Q1,Q2,Q3,……等等。因此在平面Q1,Q2,Q3,……等上面的平面积分,包括在点B1,B2,B3,……等处的氡数据。而且,因为每一个Bi都是一条来自原点O一条来自P点的这样两条垂线的交点,因此所有的Bi点都落在以线段OP为直径的圆周上。
由于Bi点全都位于图14的平面内,所以这些点全都落在其中的圆的平面,与直线SC垂直。此外,由于线段OP垂直于线段SP,所以P点位于以OS为直径的球的表面上。因而P是线段SC与氡球壳的相交点。
图13中表示的操作,可被概括为:
(1)在法向检测器平面上过C点以所有取向作出一些直线;
(2)以加权函数sinθ计算出这些线中每一条线上的J的量;
(3)计算导数dJ/dβ。
其结果是,氡数据将在包含原点O并垂直于直线SC的平面内的圆周上(图14)产生,而且以线段OP作为该圆的直径,其中P为直线SC与氡球壳相交的点。这整个操作在此被称为检测器平面内在C点处的旋转操作。
图15表示对于特定X射线源位置的归一化检测器平面。旋转轴与此归一化检测器平面在Ci点相交,而且直线62,64,66,68,70及72,为代表各积分平面Q与此归一化检测器平面的交线。为了说明图15上下文中的旋转操作,在归一化检测器平面内任意取一条直线M,并在此直线上的每一个Cj点处进行旋转操作。对于每一个Cj点都能产生一个氡数据的圆Dj,其中该圆的平面垂直于直线SCj,而且该圆的直径就是线段OPj,此处Pj是线SCj与氡球壳相交的点。由于所有的Cj点在直线M上排成一线,因此所有的点Pj,都位于包含X射线源S和线M的平面与氡球壳相交处的圆周G上。而且,因为每一个圆Dj的平面都垂直于位于圆G平面内的相应直线SPj,所以圆Dj的平面与圆G的平面垂直。
概括说来,如果旋转操作是在归一化检测器内直线M上的所有点进行的,那么氡数据就在一系列圆Dj上产生。圆Dj的平面,与包含有X射线源S及直线M的平面同氡球壳相交处的圆G的平面垂直,而且圆Dj的直径就是线段OPj,其中Pj为圆G上的点。如果在直线M上点Cj足够细地被采样,那么所产生的这组氡圆Dj就会足够密集地覆盖整个氡球壳。换而言之,整个氡球壳上的数据,可以通过在检测器平面内的一条线上进行旋转操作来产生。
计算函数J时使用的加权函数sinθ,在这种情况下是奇异的,θ=0。这种奇异性可以通过加权滑动平均处理予以消除,这种加权滑动平均处理已成功地应用在其影响函数也包含奇异性的传统X射线CT中滤波后的回填投影(filtered backprojection)的滤波部分中。
现在将考虑情况之2,其中的旋转轴与图6中的b轴重合。在这种情况下,图10中的α=π/2,而且两条直线L和L′彼此平行,因为转轴平行于归一化检测器平面。
参见图16,如果全部定位,都经与L正交的线OC(参见附录B),作出一些平行于L和L′的线,则这些线就是包含有b轴作为公共轴的一些不同取向的平面在检测器平面上的投影。标称这些平面为Q1,Q2,Q3,……等等。那么根据方程式(9),对于在归一化检测器上每一对具有小间隔的相邻直线上以加权函数cosθ对J的计算中,可由dJ/dβ得到投影到这一对相邻直线的平面Qj上的平面积分的导数。
这种情况表示在图17A和图17B中,其中图17A表示在b轴方向上的视图,图17B则表示作为情况之2的归一化检测器平面上的数据点和数据线。图17A的平面是包含有X射线源的c-a平面。在附录B中,它表明原点O也位于该平面内。由于b轴垂直于图17A的面,因此包含b轴作为公共轴的所有平面Qj,均作为构成以X射线源S为中心的极坐标中的线出现。在图17A中标称为A1,A2,A3,……等等的这些线,分别与平面Q1,Q2,Q3,……等等对应。从原点向这些线中的每一条作垂线,并让每一对垂线在B1,B2,B3,……定位上相交,如图所示。而且,如在附录A中表明的那样,在极坐标中由原点垂直于线A1,A2,A3,……的直线,也与平面Q1,Q2,Q3,……等等垂直。因此,对于Q1,Q2,Q3,……的平面积分,便在点B1,B2,B3,……等构成了氡数据。而且,由于每一个Bj都是一条来自原点O一条来自X射线源S的这样两条垂直线的相交处,所以这些点都落在以线段OS为直径的圆上。
鉴于Bj点全部落在图17A的面上,故氡数据圆的平面与垂直该平面的b轴垂直。于是如图17B所示,b轴平行于检测器平面上的一组平行线,其中包括L和L′,而且圆的平面与这组平行线垂直。
图16中表示的操作可概括为:
(1)于检测器平面上在所有覆盖整个平面的定位上,平行于θ方向作出一些线;
(2)在这些线上以加权函数cosθ计算出J的量,并且
(3)计算出导数dJ/dβ。
其结果是,氡数据在包含有原点O和X射线源S且垂直于检测器平面内的一组平行线的平面上的圆(图17A)上产生,而且线段OS为该圆的直径。这整个操作在此被称为以θ角平移的操作。
为了在整个氡球壳上产生氡数据,平移操作要以所有的角度在归一化检测器平面上完成。对于每一种角度θj都能产生一个氡数据圆Dj,此处圆的平面与检测器上成θj角的一些直线垂直,并包含原点O和X射线源S,而且圆的直径就是线段OS。如果此角能得到足够细的采样时,那么所产生的这组氡圆就会足够密集地覆盖整个氡球壳。
最后,我们可以着手处理我们提出要解决的任务:在图4的这组共轴的垂直平面θj上由锥形束数据产生出氡数据。
一般说来,此过程包括下文中由脚标i,j和k指定的嵌套步骤。脚标i对应于各种X射线源位置Si。因此,对于每一种X射线源位置Si,后续的步骤将予以重复。脚标j对应于这组平面θj,它可以是氡空间中指望在其上面产生平面积分的随意一组平面。对于每一种X射线源位置Si来说,在其上面能够得到氡数据的对应的氡球壳被确定,而且这些平面θj与此氡球壳的交线,在氡球壳上确定出一组圆Dij。对于每一特定的X射线源位置Si而言,脚标i保持固定,而脚标j变化以指定这组圆Dij中的每一个圆。
然后对于这些圆Dij中的每一个圆,重复更后续步骤中的每一步骤。具体地说,旋转轴被定义为一条经过此特定的源位置Si与此特定的圆Dij相交且与该圆Dij的平面垂直(并且垂直于相应的平面θj)的直线。在这些旋转轴中的每一条旋转轴上,可以确定一组积分平面Qijk,这些积分平面Qijk是以此特定的旋转轴共轴的,并与此特定的圆Dij相交,以确定对此特定积分平面Qijk的氡数据点Rijk的定位。对于每一特定的源位置Si和圆Dij来说,脚标i和j保持固定,而脚标k变化以指定每一个积分平面Qijk。应当使嵌套步骤的整个过程包括多阶的单个积分平面Qijk,其中每一单个的积分平面Qijk都对应一个单个的氡数据点Rijk。
继续此总的过程,每个多阶积分平面Qijk都在相应的直线Lijk上与归一化检测器平面相交(其取向与特定的源位置Si对应)。然后对于归一化检测器平面上这些线Lijk中的每一条线,将相应的积分平面Qijk绕转轴转动一个小的转角,以确定在相应的直线Lijk′上与此归一化检测器平面相交的平面Qijk′。等同地,让直线Lijk绕该直线上的一点旋转(情况之1),或者平移为平行线(情况之2),这要取决于具体的源位置Si,以便确定直线Lijk′和对应的积分平面Qijk′。
最后,为了求出在特定点Rijk处氡数据的径向导数,可通过分别沿着直线Lijk及Lijk′进行积分求出加权的线积分Jijk及Jijk′,并且将此加权的线积分之间的差值除以转角δβ。
现在考虑上述联案申请-〔RD-19564〕,需要在一组共轴的垂直平面上产生平面积分,以作为其氡的逆变换过程输入的特定情况,图18A象图4一样,表示这组包含Z轴作为公共轴的共轴的垂直平面φj,此处的垂直轴取作物体参照系中的Z轴。图18B表示一般锥形束扫描的情况,具有源S、原点O和在源位置上产生的氡球壳。鉴于源位置S不在Z轴上,故图18B是上文定义的情况之1的例子。图18C说明某种几何结构,表示图18A的平面φj与特定源位置Si下的图18B氡球壳相交的情形。由此可以看出,这些交线于单个平面φj内一些点的氡球壳上确立了一些圆。
图19具体表示如图18C的情况之1的情况下,如何在垂直平面φj上产生氡数据,这里的锥形束数据是在许多源位置Si中的每一位置上取得的,每一种源位置Si都会产生相应的球壳或氡球壳,如图18B、18C及19的氡球壳。而且,对于每一种源位置Si来说,相应的圆Gi,是在包含此源位置Si且垂直于平面φj的平面(即与垂直轴垂直的水平平面)内相应的氡球壳上确立的。
可以看出,对于每一种源位置Si来说,垂直平面φj中的每一个平面,都与对应此特定源位置Si的氡球壳以圆相交,此圆可称之为圆Dij。该圆Dij则在对应此圆Dij的Pij点处与圆Gi相交(换一种说法,这些平面φj在与特定的源位置Si及平面φj对应的Pij点与圆Gi相交)。此圆Dij则通过原点O,并以从O到Pij的直线作为直径。具体说来
Gi-其水平面通过源位置Si并与氡球壳相交的圆;
Dij-其包含Z轴的垂直平面φj与氡球壳相交的圆;
Hi-经过氡球壳中心的水平大圆平面。
由于原点O既在垂直平面φj内又在氡球壳上,所以它在交圆Dij上。根据对称性,经过氡球壳中心的水平大圆平面Hi,必然将垂直平面φj上的圆Dij平分。因此,Dij的中心应位于水平大圆平面Hi内。令OPij为圆Dij的直径。点Pij和水平大圆平面Hi之间的距离,等于原点O和Hi之间的距离。由此可见,点Pij落在圆Gi上,因为包含Gi的水平平面和平面Hi间的距离,等于原点O和平面Hi间的距离。
由于需要在每一个垂直平面φj上产生平面积分,而且对于每一种源位置Si来说,由于这些平面积分只能在对应的氡球壳上产生,所以,对于每种源位置Si下作为此特定源位置Si的每个圆Dij上的平面积分来说,Si将作为一条对应的直线Mi把Gi投影到归一化检测器平面上。交点Pij则在直线Mi上投影到对应的点Cij。
然后于直线Mi上的这些点Cij中的每一个点上进行旋转操作。在每一个Cij点上进行的旋转操作将在与此特定点Cij对应的整个圆Dij上产生出氡数据。因此,图19中的每个圆Dij,可以用图14中具有直径OP的圆表示,而且图14中的P点代表图19中的每个交点Pij。对于每一个Cij点来说,旋转轴就是经过此特定源位置Si、对应点Pij和Cij的一条直线。当对于相应于特定源位置Si的特定直线Mi上的所有点Cij完成旋转操作之后,作为此特定源位置Si的整个氡球壳上的氡数据便产生出来,并在平面φj上以圆Dij出现。当对于每一种源位置Si重复上述操作之后,氡数据便被填在在所有需要的垂直平面上。
总结这种旋转操作,对于图19中的每一个投影点Cij来说,在归一化检测器平面内过此投影点Cij以多种取向作出一些直线Lijk。这些直线Lijk就是图13及14中表示的相应积分平面Qijk在归一化检测器平面上的交线。这些积分平面中的每一个,都包含一个沿着通过该特定源位置Si、特定点Pij及特定投影点Cij的一条直线的旋转轴。如图14中的多阶图形可以绘出,每一图形的平面都与此特定旋转轴垂直。
在归一化检测器平面上,将每一条Lijk直线绕着投影点旋转一个小角度δθ,以便确立线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线。平面Qijk和Qijk′之间的转角δβ,可根据几何关系由角δθ求出来。然后按照上文描述的方法,通过沿直线Lijk及Lijk′积分,分别求出加权的线积分Jijk和Jijk′。最后,将此加权的线积分Jijk和Jijk′之间的差值除以转角δβ,以得出平面Qijk与圆Dij相交处圆Dij上点的氡数据的径向导数。
给定两条检测线之间的角度δθ,确定两积分平面之间的转角δβ,可通过一些几何处置来完成。下面是对于情况之1的情况的几何推导公式:
δβ=cosφjcosη (1+tan2θ)/(1+cos2φjcos2ηtan2θ) δθ
其中,φj-平面φj相对SO的方位角;
η-SO的极坐标角;
θ-直线Lijk和平面φj与归一化检测器平面相交处的参考线之间的夹角。
当源位置Si是在Z轴上,不管是直接在原点O之上还是直接在其之下时(属于上文中规定的情况之2的情形),上述情况之1的几何处置不能应用。
具体说来,参见图20中表示的源Si直接在原点O之下的情况。这种情况也表示在图16,17A及17B中,其中每条转轴都经过源位置Si而与归一化检测器平面平行。在图20中,每个以Z轴作为公共轴或参考轴的垂直平面φj,以OSi作为直径的氡球壳上的大圆Di相交。每个圆Dij都有作为积分平面的对应转轴,与圆Dij及其相应的平面φj垂直。
对于情况之2来说,每个特定的圆Dij,在归一化检测器平面上从源位置Si投影成直线Lij*。由于Dij圆包含O和Si,所以直线Lij*也就是平面φj在归一化检测器平面上的投影。
为了在圆Dij上产生出氡数据,在归一化检测器平面上垂直于直线Lij*作一些平行线Lijk。这些线Lijk由在图16及17B中归一化检测器平面上的一些平行线代表,并且是对应积分平面Qijk在检测器平面上的交线,其中每个积分平面Qijk,都包含一条沿着经过此特定的源位置Si且与特定圆Dij的平面垂直的直线的旋转轴。
将平行线Lijk中的每一条都平移一小段距离,以确定线Lijk′,此线就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线。这两个平面Qijk和Qijk′之间的转角δβ,可根据几何关系由直线Lijk和Lijk′之间的距离求出。然后按上文描述的方法,通过沿直线Lijk及Lijk′进行积分,分别求出加权的线积分Jijk和Jijk′。最后,将此加权线积分Jijk和Jijk′之间的差值除以转角δβ,以便得到平面Qijk与圆Dij相交处圆Dij上点的氡数据的轴向导数。
对于情况之2的情况来说,确定两积分平面之间的转角δβ,可根据以下几何推导公式来完成:
δβ= (︱so︱δy)/(︱so︳2+y2)
其中,|SO|-S和O之间的距离;
y-线Lijk距O的距离;
δy-平移距离。
考虑在本发明的实施例中需要作出计算的次数,在每一种视角下,在氡空间的一点上产生数据,要求在检测器平面内的一条数据线上的计算数次,约包含N个数据点。因此,为了在氡球壳上的圆内产生数据,需要约有N×N=N2次计算。为了在覆盖氡球壳的N个圆上产生数据,其所要求计算的次数因而就约等为N×N2=N3。最后,在所有N种视角下,总的计算量为N3×N=N4。
作出同样估算的快速方法如下:为在氡空间中的一点上产生数据,需要N次计算。由于在氡空间中约存在N3个点,所以总的计算量约等于N3×N=N4。
附录A
参见图21,它表明,图11的向量|OP| 垂直于在图11中的直线U上与平面W相交的平面Q。定义出两个平面Q1和Q2,以使
1-垂直于Q1的单位量
|PP′|2=λ′ 1 2+λ2 2+(λ3-λ3′)
对于一个固定的P来说,|PP′|2的最小值出现在λ1′=0及λ′3=λ3的情况下,即
P′=λ3n3
为直线L上的点,于是
PP′=P-P′
因此,非常明显,PP′垂直于直线L,并且垂直于平面Q2。
附录B
在情况之2的情况下,对图20的几何关系的分析:
所以SC垂直于L,
L在(U,V)平面上
所以OS垂直于L
由于SC垂直于L,而且OS垂直于L,故我们得出结论,包含直线SC及OS的平面垂直于直线L。由于L是沿着b轴,所以这个平面就是经过S的(c,a)平面。此平面包含点O,S和C。
Claims (8)
1、一种用来由锥形束投影数据再现物体三维图象的方法,此锥形束投影数据是以通过物体线积分的形式,对于多个X射线源位置Si中的每一个位置作为归一化检测器平面上的二维数据组而组成的,此归一化检测器平面包含原点,并垂直于由每一特定的源位置Si至该原点的直线,上述方法包括:
确定代表对于氡空间中一组平面φj上的平面积分的值,对于每一种源位置Si:
在氡空间中确定相应的球壳,在该球壳上的氡数据可被求出,确定平面φj与对应于特定源位置Si的球壳的交线,并在该球壳上确定一组圆Dij,而且,对于这组圆Dij中的每个圆来说:
将旋转轴定义为,经过此特定的源位置Si而与此特定圆Dij相交,且与此特定圆Dij的平面正交的一条直线;
确定一组共轴的积分平面Qijk,此积分平面Qijk中的每一个都包含此特定的旋转轴,并且与此特定的圆Dij相交以确定氡数据点Rijk的位置,而且此积分平面Qijk在相应的直线Lijk上与归一化检测器平面相交,而且,对于此归一化检测器平面上的每条直线Lijk来说:
将相应的积分平面Qijk转动一个小角度δβ,以确定以相应的直线Lijk,与此归一化检测器平面相交的平面Qijk′;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,而且
将此加权的线积分Jijk和Jijk′之差除以旋转角δβ,以得出在特定点Rijk处氡数据的径向导数;而且
在表示对于这组平面φj的平面积分值上进行逆氡变换,以便再现该物体的图象。
2、根据权利要求1的方法,其特征在于:平面φj包括含有与原点相交的参考轴的一组共轴的平面;而且,确定代表对于这组平面φj的平面积分值的步骤,包括对于不在参考轴上的每一个源位置Si:
在氡空间中对应的球壳上,在包含此特定的源位置Si并垂直于平面φj的平面内,确定出相应的圆Ci,确定此平面φj和圆Dij与该特定圆Ci的交点,并在此圆Ci上确定与圆Dij对应的多个点Pij;
将相应的圆Ci,由特定的源位置Si在归一化的检测器平面上投影为一条直线Mi;将Pij点投影至直线Mi上的对应点Cij,并且,对于归一化检测器平面内的每个投影点Cij:
在此归一化检测器平面内经过此被投影的点以多种取向作一些直线Lijk,这些线Lijk就是相应积分平面Qijk在归一化检测器平面上的交线,其中每个积分平面Qijk都包含一条沿着通过此特定的源位置Si、特定的点Pij以及特定的被投影点Cij的直线的旋转轴;
在此归一化检测器平面内,将这些线Lijk中的每一条绕着该投影点Cij转动一个小角度δθ,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并且可根据几何关系从角度δθ求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,而且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
3、根据权利要求1的方法,其特征在于:平面φj包括含有与原点相交的参考轴的一组共轴的平面;而且,确定代表对于这组平面φj的平面积分值的步骤,包括对于参考轴上的每一种源位置Si:
对于每一个同对应于此特定源位置Si的球壳相交的平面φj,确定一个特定的圆Dij;
将此特定的圆Dij从特定的源位置Si投影到归一化检测器平面内的线Lij*上;
在归一化检测器平面内作一些与此直线Lij*垂直的平行线Lijk,这些平行线Lijk就是相应积分平面Qijk在归一化检测器平面上的交线,其中每个积分平面Qijk都包含一条沿着通过此特定的X射线源位置Si并垂直于此特定圆Dij平面的直线的旋转轴;
将上述平行线Lijk中的每一条平移一小段距离,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并可根据几何关系从直线Lijk与Lijk′之间的距离求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,并且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
4、根据权利要求2的方法,其特征在于:确定代表对于这组平面φj的平面积分值的步骤,包括对于参考轴上的每一种X射线源位置:
对于每一个同对应于此特定源位置Si的球壳相交的平面φj,确定一个特定的圆Dij;
将此特定的圆Dij从特定的源位置Si投影到归一化检测器平面内的线Lij*上;
在归一化检测器平面内作一些与此直线Lij*垂直的平行线Lijk,这些平行线Lijk就是相应积分平面Qijk在归一化检测器平面上的交线,其中每个积分平面Qijk都包含一条沿着通过此特定的源位置Si并垂直于此特定圆Dij平面的直线的旋转轴;
将上述平行线Lijk中的每一条平移一小段距离,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并可根据几何关系从线Lijk与Lijk′之间的距离求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,并且
将此加权线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
5、一种用来由锥形束投影数据再现物体三维图象的装置,此锥形束投影数据是以通过物体线积分的形式,对于多个X射线源位置Si中的每一个位置作为归一化检测器平面上的二维数据而组成的,此归一化检测器平面包含原点,并垂直于由每一特定的源位置Si至该原点的直线,上述装置包括:
用来确定代表对于氡空间中一组平面φj上的平面积分值的装置,对于每一种源位置Si:
在氡空间中确定相应的球壳,在该球壳上的氡数据可被求出,确定平面φj与对应于特定源位置Si的球壳的交线,并在该球壳上确定一组圆Dij,而且,对于这组圆Dij中的每个圆来说:
将旋转轴定义为,经过此特定的源位置Si而与此特定圆Dij相交,且与此特定圆Dij的平面正交的一条直线;
确定一组共轴的积分平面Qijk,此积分平面Qijk中的每一个都包含此特定的旋转轴,并且与此特定的圆Dij相交以确定氡数据点Rijk的位置,而且此积分平面Qijk在相应的直线Lijk上与归一化检测器平面相交,而且,对于此归一化检测器平面上的每条直线Lijk来说:
将相应的积分平面Qijk转动一个小角度δβ,以确定以相应的直线Lijk′与此归一化检测器平面相交的平面Qijk′;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,而且
将此加权的线积分Jijk和Jijk′之差除以旋转角δβ,以得出在特定点Rijk处氡数据的径向导数;以及
用来在表示对于这组平面φj的平面积分值上进行逆氡变换,以便再现该物体图象的装置。
6、根据权利要求5的装置,其特征在于:平面φj包括含有与原点相交的参考轴的一组共轴的平面;而且,上述用来确定代表对于这组平面φj的平面积分值的装置,是可以操作的,对于不在参考轴上的每一种源位置Si:
在氡空间中对应的球壳上,在包含此特定的源位置Si并垂直于平面φj的平面内,确定出相应的圆Ci,确定此平面φj和圆Dij与该特定圆Ci的交点,并在此圆Ci上确定与圆Dij对应的多个点Pij;
将相应的圆Ci,由特定的源位置Si投影为在归一化的检测器平面上的一条直线Mi;将Pij点投影至直线Mi上的对应点Cij,并且,对于归一化检测器平面内的每个投影点Cij:
在此归一化检测器平面内经过此被投影的点以多种取向作一些直线Lijk,这些线Lijk就是相应积分平面Qijk在归一化检测器平面上的交线,其中每个积分平面Qijk都包含一条沿着通过此特定的源位置Si、特定的点Pij以及特定的被投影点Cij的直线的旋转轴;
在此归一化检测器平面内,将这些线Lijk中的每一条绕着该投影点Cij转动一个小角度δθ,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并且可根据几何关系从角度δθ中求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,而且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
7、根据权利要求5的装置,其特征在于:平面φj包括含有与原点相交的参考轴的一组共轴的平面;而且,上述用来确定代表对于这组平面φj的平面积分值的装置,是可以操作的,对于不在参考轴上的每一种源位置Si:
对于每一个同对应于此特定源位置Si的球壳相交的平面φj,确定一个特定的圆Dij;
将此特定的圆Dij从特定的源位置Si投影到归一化检测器平面内的线Lij*上;
在归一化检测器平面内作一些与此直线Lij*垂直的平行线Lijk,这些平行线Lijk就是相应积分平面Qijk在归一化检测器平面上的交线,其中每个积分平面Qijk都包含一条沿着通过此特定的源位置Si并垂直于此特定圆Dij平面的直线的旋转轴;
将上述平行线Lijk中的每一条平移一小段距离,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并可根据几何关系从线Lijk与Lijk′之间的距离求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,并且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
8、根据权利要求6的装置,其特征在于:所述用来确定代表对于这组平面φj的平面积分值的装置,是可以操作的,对于参考轴上的每一种源位置Si:
对于每一个同对应于此特定源位置Si的球壳相交的平面φj,确定一个特定的圆Dij;
将此特定的圆Dij从特定的源位置Si投影到归一化检测器平面内的线Lij*上;
在归一化检测器平面内作一些与此直线Lij*垂直的平行线Lijk,这些平行线Lijk就是相应积分平面Qijk在归一化检测器平面上的交线,其中每个积分平面Qijk都包含一条沿着通过此特定的源位置Si并垂直于此特定圆Dij平面的直线的旋转轴;
将上述平行线Lijk中的每一条平移一小段距离,以确定直线Lijk′,此线Lijk′就是包含此特定旋转轴的平面Qijk′与归一化检测器平面的交线,并可根据几何关系从线Lijk与Lijk′之间的距离求出平面Qijk和Qijk′之间的转角δβ;
沿着直线Lijk和Lijk′进行积分,以求出相应的加权的线积分Jijk和Jijk′,并且
将此加权的线积分Jijk和Jijk′之差除以转角δβ,以便在圆Dij上得出平面Qijk与该圆Dij交点处氡数据的径向导数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US631,815 | 1990-12-21 | ||
US07/631,815 US5257183A (en) | 1990-12-21 | 1990-12-21 | Method and apparatus for converting cone beam X-ray projection data to planar integral and reconstructing a three-dimensional computerized tomography (CT) image of an object |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011243589A Division CN1254676C (zh) | 1990-12-21 | 2001-07-26 | 由锥形束投影数据再现物体三维ct图象的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1062610A true CN1062610A (zh) | 1992-07-08 |
Family
ID=24532866
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN91111691A Pending CN1062610A (zh) | 1990-12-21 | 1991-12-20 | 将锥形束x射线投影数据再现物体三维ct图象用的方法和装置 |
CNB011243589A Expired - Fee Related CN1254676C (zh) | 1990-12-21 | 2001-07-26 | 由锥形束投影数据再现物体三维ct图象的方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011243589A Expired - Fee Related CN1254676C (zh) | 1990-12-21 | 2001-07-26 | 由锥形束投影数据再现物体三维ct图象的方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US5257183A (zh) |
EP (1) | EP0492895B1 (zh) |
JP (1) | JPH07120411B2 (zh) |
CN (2) | CN1062610A (zh) |
CA (1) | CA2056482A1 (zh) |
DE (1) | DE69129100T2 (zh) |
IL (1) | IL100313A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102988074A (zh) * | 2006-04-14 | 2013-03-27 | 威廉博蒙特医院 | 扫描狭槽锥形束计算机断层摄影以及扫描聚焦光斑锥形束计算机断层摄影 |
CN111190226A (zh) * | 2020-02-18 | 2020-05-22 | 中国石油大学(华东) | 一种三维地震数据面波噪声压制方法 |
Families Citing this family (69)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5390226A (en) * | 1992-07-02 | 1995-02-14 | General Electric Company | Method and apparatus for pre-processing cone beam projection data for exact three dimensional computer tomographic image reconstruction of a portion of an object |
US6219441B1 (en) * | 1993-06-22 | 2001-04-17 | General Electric Company | Reconstruction of images from three-dimensional cone beam data |
US5446776A (en) * | 1993-08-02 | 1995-08-29 | General Electric Company | Tomography with generation of radon data on polar grid points |
US5420788A (en) * | 1993-08-09 | 1995-05-30 | General Electric Company | Method for processing data in 2D and 3D computed X-ray tomography |
US5461650A (en) * | 1993-10-18 | 1995-10-24 | General Electric Company | Method and system for pre-processing cone beam data for reconstructing free of interpolation-induced artifacts a three dimensional computerized tomography image |
US5390111A (en) * | 1993-11-12 | 1995-02-14 | General Electric Company | Method and system for processing cone beam data for reconstructing free of boundary-induced artifacts a three dimensional computerized tomography image |
US5463666A (en) * | 1993-11-12 | 1995-10-31 | General Electric Company | Helical and circle scan region of interest computerized tomography |
US5461651A (en) * | 1993-11-17 | 1995-10-24 | General Electric Company | Reconstruction of images in cone beam scanning with rectangular detector elements |
US5400255A (en) * | 1994-02-14 | 1995-03-21 | General Electric Company | Reconstruction of images from cone beam data |
US5465283A (en) * | 1994-06-06 | 1995-11-07 | General Electric Company | Method and apparatus for generating radon derivative data on a near uniformly space grid in a cone beam computerized tomography implementation |
CA2156141A1 (en) * | 1994-09-28 | 1996-03-29 | Kaveh Azar | Interactive scanning device or system |
US5500883A (en) * | 1994-12-05 | 1996-03-19 | General Electric Company | Method and system for efficiently generating radon derivative data in a cone beam computerized tomography implementation |
US5504792A (en) * | 1994-12-27 | 1996-04-02 | General Electric Company | Method and system for masking cone beam projection data generated from either a region of interest helical scan or a helical scan |
US5521954A (en) * | 1995-02-21 | 1996-05-28 | General Electric Company | Method and apparatus for improved accuracy in computed tomography cone-beam imaging |
US5594766A (en) * | 1995-03-31 | 1997-01-14 | General Electric Company | Method and apparatus for improving efficiency in computing radon derivatives for CT imaging |
US5559846A (en) * | 1995-09-22 | 1996-09-24 | General Electric Company | Technique for improving processor efficiency in CT imaging |
US5805659A (en) * | 1996-09-30 | 1998-09-08 | Siemens Corporate Research, Inc. | Method and apparatus for spiral scan region of interest imaging |
US5748697A (en) * | 1996-12-20 | 1998-05-05 | Siemens Corporate Research, Inc. | Method and apparatus for eliminating boundary errors in cone beam imaging |
US5999587A (en) * | 1997-07-03 | 1999-12-07 | University Of Rochester | Method of and system for cone-beam tomography reconstruction |
US6075836A (en) | 1997-07-03 | 2000-06-13 | University Of Rochester | Method of and system for intravenous volume tomographic digital angiography imaging |
US5901195A (en) * | 1997-09-30 | 1999-05-04 | Siemens Corporate Research, Inc. | Two-step radon inversion processing for φ-planes having local radon origins |
US5862198A (en) * | 1997-09-30 | 1999-01-19 | Siemens Corporate Research, Inc. | Pre-calculated hitlist for reducing run-time processing of an exact cone beam reconstruction algorithm |
US5901196A (en) * | 1997-09-30 | 1999-05-04 | Siemens Corporate Research, Inc. | Reduction of hitlist size in spiral cone beam CT by use of local radon origins |
JP4146907B2 (ja) * | 1998-01-13 | 2008-09-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 診断区域の螺旋走査工程を含むコンピュータトモグラフィ法 |
US5881123A (en) * | 1998-03-31 | 1999-03-09 | Siemens Corporate Research, Inc. | Simplified cone beam image reconstruction using 3D backprojection |
US6009142A (en) * | 1998-03-31 | 1999-12-28 | Siemens Corporate Research, Inc. | Practical cone beam image reconstruction using local regions-of-interest |
US6078638A (en) * | 1998-09-30 | 2000-06-20 | Siemens Corporate Research, Inc. | Pixel grouping for filtering cone beam detector data during 3D image reconstruction |
US6169779B1 (en) * | 1998-10-20 | 2001-01-02 | Analogic Corporation | Data acquisition sequence for computed tomography scanners |
JP2000197627A (ja) * | 1999-01-05 | 2000-07-18 | Hitachi Medical Corp | X線ct装置 |
US6687393B1 (en) | 1999-07-16 | 2004-02-03 | General Electric Company | Efficient methods and apparatus for resampling three dimensional datasets |
US6411670B1 (en) | 1999-11-17 | 2002-06-25 | General Electric Company | Data rebinning to increase resolution in CT image reconstruction |
US6480565B1 (en) | 1999-11-18 | 2002-11-12 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
US6519239B1 (en) | 1999-11-19 | 2003-02-11 | Motorola, Inc. | Method and apparatus for providing dispatch service in a CDMA communication system |
US6671515B1 (en) | 2000-06-06 | 2003-12-30 | Motorola, Inc. | Method and apparatus for selecting communication cells in a wireless communication system |
US6504892B1 (en) | 2000-10-13 | 2003-01-07 | University Of Rochester | System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit |
US6862273B2 (en) * | 2001-01-10 | 2005-03-01 | Motorola, Inc. | Method and apparatus for providing dispatch scan in a CDMA communication system |
US6477221B1 (en) | 2001-02-16 | 2002-11-05 | University Of Rochester | System and method for fast parallel cone-beam reconstruction using one or more microprocessors |
US6928063B2 (en) * | 2001-03-16 | 2005-08-09 | Motorola, Inc. | Method and apparatus for providing a dispatch patch service in a CDMA communication system |
US6965108B2 (en) * | 2001-07-30 | 2005-11-15 | Euro-Celtique, S.A. | Method and apparatus for three dimensional imaging using infrared radiation |
US6574297B2 (en) | 2001-10-30 | 2003-06-03 | Siemens Corporate Research, Inc. | System and method for image reconstruction in a cone beam imaging system |
US6961572B2 (en) | 2001-12-20 | 2005-11-01 | Motorola, Inc. | Method and apparatus for base-initiated, CDMA-dispatch soft handoff |
US6788941B2 (en) | 2001-12-20 | 2004-09-07 | Motorola, Inc. | Method and apparatus for mobile-initiated, CDMA-dispatch soft handoff |
US7099290B2 (en) * | 2001-12-20 | 2006-08-29 | Motorola, Inc. | Method and apparatus for CDMA-dispatch soft handoff |
US6801783B2 (en) | 2001-12-20 | 2004-10-05 | Motorola, Inc. | Base site and method for quickly establishing a CDMA dispatch call |
US6925309B2 (en) | 2001-12-20 | 2005-08-02 | Motorola, Inc. | Method and apparatus for quickly establishing a CDMA dispatch call |
US6763081B2 (en) * | 2002-10-02 | 2004-07-13 | Siemens Corporate Research, Inc. | Cone beam computed tomography imaging system and method providing efficient utilization of detector area |
US6983034B2 (en) | 2003-02-14 | 2006-01-03 | University Of Iowa Research Foundation | Methods and devices for CT reconstruction using a grangeat approach |
DE10313110A1 (de) * | 2003-03-24 | 2004-10-21 | Sirona Dental Systems Gmbh | Röntgeneinrichtung und röntgenstrahlenempfindliche Kamera |
US8804899B2 (en) | 2003-04-25 | 2014-08-12 | Rapiscan Systems, Inc. | Imaging, data acquisition, data transmission, and data distribution methods and systems for high data rate tomographic X-ray scanners |
US8204173B2 (en) * | 2003-04-25 | 2012-06-19 | Rapiscan Systems, Inc. | System and method for image reconstruction by using multi-sheet surface rebinning |
US7209535B2 (en) * | 2003-06-20 | 2007-04-24 | Wisconsin Alumni Research Foundation | Fourier space tomographic image reconstruction method |
US6990167B2 (en) * | 2003-08-29 | 2006-01-24 | Wisconsin Alumni Research Foundation | Image reconstruction method for divergent beam scanner |
DE102004018498A1 (de) * | 2004-04-14 | 2005-11-17 | Siemens Ag | Betriebsverfahren für eine Röntgenanlage, rechnergestütztes Ermittlungsverfahren für mindestens eine 3D-Rekonstruktion eines Objekts und hiermit korrespondierende Einrichtungen |
US7394887B2 (en) * | 2004-10-15 | 2008-07-01 | General Electric Company | Method and apparatus for reconstruction of tilted cone beam data |
EP1828985A1 (en) * | 2004-11-24 | 2007-09-05 | Wisconsin Alumni Research Foundation | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data |
CN101231254B (zh) * | 2007-01-24 | 2010-08-11 | 清华大学 | 双源三维立体成像系统和成像方法 |
JP4991487B2 (ja) * | 2007-11-01 | 2012-08-01 | 株式会社日立製作所 | 情報記憶装置及び記憶媒体 |
US9858716B2 (en) * | 2008-02-28 | 2018-01-02 | International Business Machines Corporation | Fast three-dimensional visualization of object volumes without image reconstruction by direct display of acquired sensor data |
JP5405187B2 (ja) * | 2009-05-08 | 2014-02-05 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
US8160200B2 (en) | 2010-03-30 | 2012-04-17 | General Electric Company | Method and system for image data acquisition |
CN101980302A (zh) * | 2010-10-22 | 2011-02-23 | 南方医科大学 | 投影数据恢复导引的非局部平均低剂量ct重建方法 |
US9364191B2 (en) | 2013-02-11 | 2016-06-14 | University Of Rochester | Method and apparatus of spectral differential phase-contrast cone-beam CT and hybrid cone-beam CT |
US10169909B2 (en) * | 2014-08-07 | 2019-01-01 | Pixar | Generating a volumetric projection for an object |
CN104899827B (zh) * | 2015-05-26 | 2018-04-10 | 大连理工大学 | 基于固定分辨率条件下的离散Radon投影和Mojette投影转换方法 |
CN109521456B (zh) * | 2018-09-25 | 2022-10-21 | 中国辐射防护研究院 | 一种基于正则化最小二乘法的γ辐射源项反演方法及系统 |
GB2586256B (en) * | 2019-08-14 | 2023-02-22 | Adaptix Ltd | An X-ray tomography system and method |
CN110702707A (zh) * | 2019-10-16 | 2020-01-17 | 四川轻化工大学 | 一种获取核废物桶层析γ扫描图像的方法 |
CN110866959B (zh) * | 2019-11-12 | 2023-09-15 | 上海联影医疗科技股份有限公司 | 图像重建方法、系统、装置及存储介质 |
CN114295653B (zh) * | 2021-11-30 | 2023-12-22 | 太仓市林源电线电缆有限公司 | 一种电工铜线的多段式在线连续检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5073910A (en) * | 1990-08-27 | 1991-12-17 | General Electric Company | Square wave cone beam scanning trajectory for data completeness in three-dimensional computerized tomography |
US5068882A (en) * | 1990-08-27 | 1991-11-26 | General Electric Company | Dual parallel cone beam circular scanning trajectories for reduced data incompleteness in three-dimensional computerized tomography |
-
1990
- 1990-12-21 US US07/631,815 patent/US5257183A/en not_active Expired - Lifetime
-
1991
- 1991-11-28 CA CA002056482A patent/CA2056482A1/en not_active Abandoned
- 1991-12-06 JP JP3348602A patent/JPH07120411B2/ja not_active Expired - Fee Related
- 1991-12-10 IL IL10031391A patent/IL100313A/en not_active IP Right Cessation
- 1991-12-11 DE DE69129100T patent/DE69129100T2/de not_active Expired - Fee Related
- 1991-12-11 EP EP91311514A patent/EP0492895B1/en not_active Expired - Lifetime
- 1991-12-20 CN CN91111691A patent/CN1062610A/zh active Pending
-
2001
- 2001-07-26 CN CNB011243589A patent/CN1254676C/zh not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102988074A (zh) * | 2006-04-14 | 2013-03-27 | 威廉博蒙特医院 | 扫描狭槽锥形束计算机断层摄影以及扫描聚焦光斑锥形束计算机断层摄影 |
CN111190226A (zh) * | 2020-02-18 | 2020-05-22 | 中国石油大学(华东) | 一种三维地震数据面波噪声压制方法 |
CN111190226B (zh) * | 2020-02-18 | 2022-03-25 | 中国石油大学(华东) | 一种三维地震数据面波噪声压制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1380542A (zh) | 2002-11-20 |
US5257183A (en) | 1993-10-26 |
CN1254676C (zh) | 2006-05-03 |
IL100313A0 (en) | 1992-09-06 |
EP0492895A3 (en) | 1992-11-25 |
EP0492895B1 (en) | 1998-03-18 |
DE69129100D1 (de) | 1998-04-23 |
JPH07120411B2 (ja) | 1995-12-20 |
JPH0628454A (ja) | 1994-02-04 |
DE69129100T2 (de) | 1998-10-08 |
EP0492895A2 (en) | 1992-07-01 |
IL100313A (en) | 1996-01-19 |
CA2056482A1 (en) | 1992-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1254676C (zh) | 由锥形束投影数据再现物体三维ct图象的方法 | |
CN1658795A (zh) | 用于快速发散波束断层摄影术的方法和装置 | |
CN1788684A (zh) | 图像重建方法和x射线计算机断层摄影设备 | |
CN101048802A (zh) | 从氡数据重建(n+1)维图像函数的方法和设备 | |
CN1289034C (zh) | X射线ct装置、三维图象再次构成方法及权重设定方法 | |
CN101032409A (zh) | 图像显示设备和x射线ct设备 | |
CN1261077C (zh) | 计算机x射线断层造影设备 | |
CN1290513A (zh) | 螺旋锥形束扫描计算机层析x射线摄影系统中立体图象的再现装置和方法 | |
CN1688254A (zh) | X线断层摄影装置 | |
CN1268035A (zh) | 利用螺旋扫描ct再现立体图像 | |
CN1650809A (zh) | 图像重构方法和x线计算机体层摄影设备 | |
CN1063171A (zh) | 再现三维计算机断层扫描图象的方法和装置 | |
CN1518956A (zh) | 三维反向投影方法和x射线ct装置 | |
EP0492896A2 (en) | Reconstructing a 3-D image | |
CN1236729C (zh) | 图像重构方法和射线照相装置 | |
CN1284651A (zh) | 计算机层析x射线摄影系统中的图像再现装置和方法 | |
CN101052352A (zh) | 图像诊断处理设备和图像诊断处理程序 | |
US20030081715A1 (en) | System and method for image reconstruction in a cone beam imaging system | |
CN85108635A (zh) | 图象数据处理方法及其系统 | |
JP4342164B2 (ja) | コンピュータ断層撮影装置 | |
WO2012055147A1 (zh) | 计算机断层成像(ct)图像重建方法及装置 | |
CN101034479A (zh) | 图像重建方法和x射线ct设备 | |
Li et al. | A curve-filtered FDK (C-FDK) reconstruction algorithm for circular cone-beam CT | |
Lesaint | Data consistency conditions in X-ray transmission imaging and their application to the self-calibration problem. | |
JP6243219B2 (ja) | 画像生成装置および放射線断層撮影装置並びにプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |