CN101594825A - 锥束ct动态成像的方法和设备 - Google Patents
锥束ct动态成像的方法和设备 Download PDFInfo
- Publication number
- CN101594825A CN101594825A CNA200780013865XA CN200780013865A CN101594825A CN 101594825 A CN101594825 A CN 101594825A CN A200780013865X A CNA200780013865X A CN A200780013865XA CN 200780013865 A CN200780013865 A CN 200780013865A CN 101594825 A CN101594825 A CN 101594825A
- Authority
- CN
- China
- Prior art keywords
- prime
- scan
- cone
- data
- projection
- 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
- 238000000034 method Methods 0.000 title claims description 30
- 238000003384 imaging method Methods 0.000 title description 12
- 229910052704 radon Inorganic materials 0.000 claims description 18
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 16
- 230000005855 radiation Effects 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 7
- 238000007408 cone-beam computed tomography Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims 2
- 230000000007 visual effect Effects 0.000 claims 1
- 208000030984 MIRAGE syndrome Diseases 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 12
- TVLSRXXIMLFWEO-UHFFFAOYSA-N prochloraz Chemical compound C1=CN=CN1C(=O)N(CCC)CCOC1=C(Cl)C=C(Cl)C=C1Cl TVLSRXXIMLFWEO-UHFFFAOYSA-N 0.000 description 12
- 238000010408 sweeping Methods 0.000 description 9
- 230000009467 reduction Effects 0.000 description 6
- 230000002123 temporal effect Effects 0.000 description 4
- 238000005094 computer simulation Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000021615 conjugation Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000001915 proofreading effect Effects 0.000 description 2
- 150000003257 radon Chemical class 0.000 description 2
- 230000008707 rearrangement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229910052709 silver Inorganic materials 0.000 description 2
- 239000004332 silver Substances 0.000 description 2
- 238000002583 angiography Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000001990 intravenous administration Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010603 microCT Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000004899 motility Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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/04—Positioning of patients; Tiltable beds or the like
- A61B6/0407—Supports, e.g. tables or beds, for the body or parts of the body
- A61B6/0435—Supports, e.g. tables or beds, for the body or parts of the body with means for imaging suspended breasts
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
-
- 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Theoretical Computer Science (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Pulmonology (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种锥束扫描***沿着半圆形扫描。该重建使用了一个对于距离该扫描平面更远的行减小的加权函数以考虑投影数据的冗余度。另一实施例使用圆形加稀疏螺旋扫描几何。图像数据是实时采集的。
Description
相关申请的引用
本申请要求2006年2月27日提交的美国临时专利申请第60/776666号的权益,这里通过引用将其公开内容结合到本说明书中。
政府利益声明
研发本发明的工作部分由NIH Grants 8 R01 EB002775、R01 9HL078181和4 R33 CA94300支持。政府拥有本发明中的某些权利。
技术领域
本发明涉及成像,更特别地涉及使用加权函数增加时间分辨率的半扫描成像技术。
背景技术
近几年来,半扫描方法在锥束CT中的使用成为一个热门课题,这是由于时间分辨率的提高。目前有多种不同类型的锥束半扫描方案,例如基于FDK的,基于锥束滤波反向投影的(CBFBP),和基于Grangeat的。每个方案使用平面扫描轨道(圆形的或非圆形的)来实施该半扫描方案。理论上,只要该重建物体在某个尺寸大小内,圆形半扫描可以基于Grangeat公式的导数,根据一阶导数辐射数据可以近似获得与圆形完全扫描相同的Radon域信息。即使在圆形半扫描范围内,仍然存在冗余。Grangeat型半扫描(GHS)将空间投影数据映射到一阶导数辐射数据并且在Radon域对它们加权。在通过线性内插/外插圆形扫描不能进入的该Radon域的阴影区而添加丢失数据之后,使用3DRadon逆公式来得到该重建图像。
当前用于锥束CT的FDK型半扫描(FDKHSFW)方案使用Parker或者其他基于扇形束几何的加权系数,其中对所有检测行应用相同的加权系数。CBFBP算法使用该Radon域中的冗余投影数据,在滤波反向投影(FBP)的结构中完成半扫描重建并且获得几乎与FDKHSFW相同的性能。当该阴影区被填充线性内插数据时,Grangeat型半扫描方案在校正未扫描(off-scanning)平面衰减系数降低上的性能超过了FDK型半扫描方案。然而,来自GHS的重建图像的空间分辨率比FDKHSFW差,因为在FDK中包含的数据内***GHS少。此外,GHS不能处理纵向方向上的被截(truncated)数据。CBFBP相关的半扫描和FDKHSFW在距离Z=0最远的重建图像位置显示了各种衰减系数降低伪像(artifacts),其中Z是旋转轴。该伪像是实践中不希望出现的。
发明内容
本发明的一个目的是在一定程度上校正该降低问题。
本发明的另一个目的是在这样做的同时保持空间分辨率。
为了实现上述和其他目的,本发明提出了一种具有配合锥束几何(FDKHSCW)的新加权函数的FDK半扫描方案,其中该加权函数是与锥束几何相关的。该具有新锥束加权函数的FDK半扫描算法被描述。实施了一种计算机模拟,并且评估该FDKHSCW以与FDKFS和FDKHSFW相比较。
对基于FDK的圆形半扫描重建(FDKHSCW)启发式提出了一种新的锥束加权方案,用于在一定程度上沿着该旋转轴校正该密度降低伪像,接着使用用于较大锥角的原始FDK算法。当该锥角较大时,根据衰减系数降低使用FDKHSCW来代替FDKFS和FDKHSFW时,对有噪声和无噪声的Shepp-Logan幻象进行的计算机模拟显示产生了改进,同时维持相同的虚拟图象质量。FDKHSCW需要在滤波前进行额外的锥束加权,并且仅使用[β,180+β+2Δ]的扫描范围,其中β是X射线的起始投影角,Δ是完全扇形角的一半,二者都被限定在该扫描平面中。该起始角一旦确定,就可以处理每个投影图像(对半扫描进行锥束加权,由FDK继承的像素加权,滤波)。所以,与完全扫描方案相比较,重建物体需要更少的时间,这在实践中是一个非常希望具备的特征。这里,该半扫描被定义为180°加上扫描几何的锥角。此外,该半扫描方案提供了选择任何起始点用于重建的灵活性,只要该扫描范围被确保即可,这是锥束CT动态成像的另一优选特征。如果执行完全圆形扫描,为了增加该***和锥束CT重建的时间分辨率,可以通过选择不同的起始点来根据该完全圆形扫描执行多个半扫描重建。因此,本发明的另一目的是,使用半扫描算法根据物体的一个或多个圆形扫描投影图像重建多个半扫描重建,以增加锥束CT动态成像扫描的时间分辨率。
基于Silver提出的思想(M.D.Silver,“A method for includingredundant data in computed tomography”,Med.Phys.27,773-774(2000)),我们甚至能够通过使得扫描范围大于180+2Δ来实施扩展的半扫描方案,应用该新锥束加权函数以获得更好的噪声特性。
该圆形锥束半扫描加权方案对于低对比度物体工作得更好。我们从Shepp-Logan幻象的模拟中可以看到,最大补偿位于0.03的衰减系数内。我们预期,FDKHSCW能够显示高对比度幻象例如Defrise盘幻象中的强度降低的改进。但是它对于低对比度幻象更有前途。
其他提出的改进的FDK方法称为T-FDK和FDK-SLANT,它们也在一定程度上以更大的锥角校正FDK中继承的沿着旋转轴衰减系数降低。在这些方法和FDKHSCW之间存在不同。虽然这些方法的结果显示了与FDKHSCW相似的校正,但是FDK-SLANT和T-FDK需要根据锥束数据平行重排(rebinned)。这意味着,直到整组数据采集和平行重排序过程完成才开始该滤波部分,然后进行反向投影以进行图像重建。FDKHSCW的优点在于,一旦获得2D投影数据,该滤波部分就能开始并且立即进行反向投影。只要托架速度和读出速率足够高,这个方案就能在实施连续动态成像时实现几乎实时监测。Wang(G.Wang,“X-ray micro-CT witha displaced detector array”,Med.Phys.29,1634-1636(2002))研发了一种在移位的(displaced)检测器阵列上进行锥束完全圆形扫描重建而不需要重排该投影数据以重建的加权方案。对于该冗余区域,可以通过调节扫描范围中的加权条件来将我们的方案应用到这个算法中。
近来,提出了一种新的圆形3D加权重建算法,基于对直接射线与其共轭射线之间的数据不一致的研究来减少锥束伪像。其基本思想是在反向投影期间将滤波的投影数据乘以与锥束几何相关的校正系数。但是它校正的伪像不是这里FDKHSCW试图校正的,即衰减系数降低。然而,这两种方案可以相结合。
总之,通过结合新的锥束加权方案,基于一个平板检测器提出了一种新的基于FDK的探索性的圆形轨道半扫描近似算法,数值模拟证实了它的可行性。
由于可以以更少的时间来完成该数据的扫描和处理,所以可以更快地完成成像。因而,四维扫描(三维空间和一维时间)就变得可能。
在以下美国专利中公开了相关的***和方法:6987831号,“Apparatus and method for cone beam volume computed tomographybreast imaging”;6618466号,“Apparatus and method for x-rayscatter reduction and correction for fan beam CT and cone beamvolume CT”;6504892号,“Apparatus and method for cone beam volumecomputed tomography using circle-plus-multiple-arc orbit”;6480565号,“Apparatus and method for cone beam volume computedtomography breast imaging”;6477221号,“System and method forfast parallel cone-beam reconstruction using one or moremicroprocessors”;6298110号,“Cone beam volume CT angiographyimaging system and method”;6075836号,“Method of and systemfor intravenous volume tomographic digital angiography imaging”;和5999587号,“Method of and system for cone-beam tomographyreconstruction”。这里通过引用将它们的全部内容结合到本说明书中。这些专利中公开的技术可以用于与这里公开的技术相结合。
附图说明
下面参照附图具体给出本发明的优选实施例,其中:
图1示出了第一优选实施例的扫描几何;
图2示出了图1的扫描几何中的冗余区域;
图3A-3D示出了加权系数;
图4A-4C示出了幻象的重建矢状图像;
图4D示出了该幻象;
图4E和4F沿着图4D中所示的线对图4A-4C中的图像的线剖面(profile)进行比较;
图5A和5B示出了具有不同噪声水平的结果;
图5C和5D示出了对于图5A和5B的结果的剖面比较;
图6示出了第二优选实施例中的Radon域;和
图7A和7B示出了第二优选实施例中的扫描几何。
具体实施方式
下面将参照附图具体给出本发明的优选实施例,其中相似的参考数字表示相似的部件或步骤。
通过对所有倾斜扇形波束对物体的作用求和而在扇形波束算法的基础上扩展FDK算法。该重建是基于滤波和在该锥形内反向投影单个扇形波束的。基于图1中的锥束几何,其中显示了锥束发射器或其他X射线源102和平板检测器104的平面的相对位置,FDK的公式为:
s=-xsinβ+ycosβ;
*号表示卷积;so表示从X射线源到起点(origin)的距离;n,m表示整数值,其中n=0和m=0对应于经过该起点的中心射线;β表示限定在该扫描平面中的投影角;p表示沿着t轴的虚拟检测器采样间隔;ξ表示沿着Z轴的虚拟检测器采样间隔;Rβ(np,mξ)表示实际的离散2D投影数据;h(np)表示沿着t轴的离散一维斜坡滤波器冲击响应。
预加权项可以分解为两个余弦项,即 这意味着FDK将未扫描平面投影数据投影到该扫描平面中,然后执行该2D扇形波束重建算法。在公式(1)中,积分前面的因子被用于抵消当实施完全圆形扫描时的投影冗余。这就表示该未扫描平面投影数据具有与扫描平面中的投影数据相同的冗余度。
该锥束半扫描方案也是与FDK相结合的扇形波束半扫描的扩展,其中将从扫描平面几何计算的加权系数应用到所有投影行,如下:
s=-xsinβ+ycosβ;
这是FDKHSFW方案,其中Δ是沿着t轴的中心扫描平面的完全扇形角的一半。该未扫描平面投影数据仍然被作为它们具有相同冗余度进行处理。ω(β,np)是基于该扫描平面几何计算的离散加权系数,可以表示为Parker加权函数或其他任何加权函数,只要它能够在双重和单独采样区域之间产生平滑过渡以避免这些区域边界处的不连续。毫无疑问,FDKHSFW保持了FDK完全扫描方案所具有的所有特性。
然而,对于该扫描平面之外的锥束投影数据,即使在360°采样投影,单次圆形轨道采集也不能获得完全双重采样的投影。换句话说,当投影行距离该扫描平面更远时,该投影冗余度就变得越来越少。如果FDK算法被直接应用到未加权的半扫描投影数据上,那么重建图像将不可避免地具有伪像。一种处理远离扫描平面的更少冗余投影行数据的加权方式如下所示:
β′是锥形加权角度。β′依赖于Z方向(旋转轴)上的行投影数据的位置。Δ′是G.T.Gullberg与G.L.Zeng,“A cone-beam filterbackprojection reconstruction algorithm for cardiac singlephoton emission computed tomography”,IEEE Trans.Med.Imaging11,91-101(1992)中所采用的倾斜扇形角的一半。要注意当m为零时,该加权函数实际上是用于扇形波束的Parker加权函数。
通过将该锥束加权函数与FDK相结合,可以获得FDKHSCW如下:
s=-xsinβ+ycosβ;
请注意,该投影数据在滤波前必须被加权。因为FDKHSFW是半扫描重建的公知的方案,所以对于FDKHSCW的需要是它不应产生比FDKHSFW更多的伪像。
现在将说明对于半扫描锥束加权的进一步研究。在圆形扇形波束半扫描中,按照投影角度β,在扫描平面中有两个冗余区域。图2显示在区域I中采集的投影射线数据具有区域II中的共轭射线数据。在这两个区域中,投影射线数据是整体或部分冗余的。如果完全扇形角的一半是Δ度,根据限定在扫描平面内的投影角而定的半扫描范围是从0°到180°+2Δ。第一和第二冗余区域分别是从0°到4Δ和从180°-2Δ到180°+2Δ。在传统的FDK锥束半扫描方案中,所有行投影数据被限定在扫描平面内的同一组系数加权,因为希望该扫描平面外的行投影数据具有与扫描平面内的行投影数据相同的冗余度。
该圆形锥束半扫描加权方案的提出是基于这样的思想,对于不同行的投影数据,加权系数应该不同,并且对于距离扫描平面最远的行投影数据,应该加权最少。至今为止,我们还没有看到任何讨论这个问题的文献。我们发现,如果使用
作为不同行投影数据的加权角,那么在扫描平面外的第一冗余区域中的加权系数将不会与该扫描平面内计算的加权系数有很大不同;如果Δ=15°且半锥角也是15°,则最大差别低于百分之0.2。另一方面,当β′被用作第二冗余区域中的加权角时,扫描平面外的加权系数将显著不同于该扫描平面内的加权系数并且在不同行互相不同,从而在重建图像中对于该扫描平面外的位置的密度降低产生了补偿。该加权角β′具有两个特性:第一,它具有由mξ反应的行位置依赖性,间接与锥角信息相关联;第二,当β在第一冗余区域时它与β的差别比β在第二冗余区域时更小。从而,在第二冗余区域实施该锥角相关的加权系数对于实现我们的方案是有利的。
为了使计算机模拟更接近于实际的CBCT配置,将几何的参数根据其物理长度(毫米)设置而不是根据标准化的(normalized)单位。从X射线源到该重建的等角点(iso-center)和到检测器的距离分别是780mm和1109mm。该完全扇形和锥形角是30度。检测器面积是595×595mm2,具有512乘512的矩阵大小。体元(voxel)大小是0.816mm3。笛卡尔坐标(X,Y,Z)被用于描述该物体,其中Z是旋转轴。投影角的采样速率是0.8°,对于完全扫描投影图像总数是450,对于半扫描是262。低对比度Shepp-Logan幻象被使用(几何参数见参考文献7),其所有几何参数被乘以200以模拟该幻象的物理长度(毫米)。
基于上述定义的扫描几何参数,通过选取在冗余区域I中β=46°和在冗余区域II中β=192°来比较与FDKHSFW和FDKHSCW相关的加权系数分布。图3A-3D分别显示了FDKHSFW(β=46°),FDKHSCW(β=46°),FDKHSFW(β=192°)和FDKHSCW(β=192°)。
图4A-4C以显示窗口[1.005 1.05]示出了不同FDK方案(分别为FDKFS、FDKHSFW、FDKHSCW)在X=0mm时的重建矢状图像。图4D示出了幻象。图4E和4F示出了沿着图4D中的幻象图像中的白色垂直和水平实线的剖面比较。在反向投影之前对无噪声的加权投影数据使用斜坡滤波器。
为了测试这种新方案对于实际CBCT数据采集中常见的量子噪声的性能,我们生成了被量子噪声干扰的数据。100kVp的X射线被选择,其对应于2.9972*107光子/cm2·mR的有效光子通量。每个投影的曝光水平被设定为4mR,FDKFS和FDKHSCW的总曝光水平分别是1800mR和1048mR。
图5A-5D示出了不同噪声水平下的重建结果和剖面比较。在滤波期间使用哈明窗来抑制噪声。更特别地,图5A和5B分别示出了总曝光水平为1800mR的FDKFS和总曝光水平为1048mR的FDKHSCW,而图5C和5D示出了沿着图4D中的垂直和水平实线在FDKFS、FDKHSCW和幻象之间的剖面比较。
下面将给出第一优选实施例的一个变体。如图6所示,在Radon域中,该圆形(区域602)中的点表示在圆形扫描中采集的该物体的radon值,而该圆形(区域604)中的阴影区域表示通过圆形扫描不能采集到的丢失的Radon点组。如Hu所主张的,在圆形扫描中,FDK仅使用该物体圆形内和点线弧边界上的冗余点,而不使用该物体圆形内的实线弧边界上的非冗余点。因而,使用以下校正式:
s=-xsinβ+ycosβ
t=xcosβ+ysinβ
如果把上述重建表示为f1(t,s,z),那么最终重建是:
f(t,s,z)=f1(t,s,z)+f2(t,s,z)
现在将说明特别应用于胸部成像的第二优选实施例。
基于当前CBBCT的几何参数,我们设计了一种新的CBBCT扫描方案。如图7A和7B所示,在圆形扫描期间,X射线源102的位置在z=0cm。在圆形扫描后,X射线源102和检测器104同时降低,同时它们仍然在旋转。当该X射线源到达z=10的点时,它开始发射并且在间歇Δl之间保持发射直到它在该螺旋扫描中完成八次发射。对于该螺旋扫描期间的每次发射,该X射线源保持与它在圆形扫描中相同的准直。这种方案称为稀疏螺旋扫描。实际上,为了高效实施稀疏螺旋扫描,该***需要包括位于托架上的滑环,用于提供电源到托架上部件的连接以及二维检测器与计算机***之间的通信,和用于沿着旋转轴上下移动该托架或该物体所在的物体支架的电机。
图7B示出了与螺旋线扫描相关的投影角。在螺旋扫描期间有八个发射点,覆盖了4π的角度范围,Z方向上的移动是从48mm到128mm,递增间隔为10mm,基于该模拟胸部幻象的尺寸。
从这种额外的扫描轨道获得的Radon数据点中的一些仍然可以通过圆形扫描获得,这在Radon域中被称为冗余采样点,并且可以被冗余窗函数(RWF)有效消除,并且如在HL扫描期间的准直器的几何设置能够避开该螺旋线上的那些发射点之间的冗余采样Radon点。由于该螺旋线扫描期间的准直不可避免地导致纵向截断,所以将使用与几何相关的截断窗函数(TWF)来处理这种情况。
其中
用于螺旋线扫描的重建项可以设计成一种滤波反向投影(FBP):
是用于对通过螺旋线扫描采集的但是已经被圆形扫描获得的Radon数据进行分块(block out)的窗函数。基于该扫描设计,螺旋线扫描在该扫描期间将不可避免地遭遇纵向截断。它所采集的一些radon点不能反映实际的Radon数据,应该在反向投影期间被消除。窗口用于实现这一目的。
虽然上面已经公开了优选实施例及其变体,但是阅读了本说明书的本领域技术人员将会很容易认识到,在本发明的范围内也可以实现其他实施例。例如,数值是示例性的而非限制性的。而且,本发明可以在任意适当的扫描设备上实现,包括波束发射器、平板或其他二维检测器或其他适当检测器以及用于该二者相对移动的托架的任意适当的组合,以及用于处理该图像数据以产生图像的计算机和适当输出(例如显示器或打印机)或用于该图像的存储介质。用于执行本发明的软件可以在任意介质上实现为任意适当的格式,例如,物理介质如CD-ROM或通过因特网或内联网的连接。因此,本发明应当被解释为仅由所附的权利要求限定。
Claims (33)
1.一种用于形成物体的图像的方法,该方法包括:
(a)通过使用锥束发射器和检测器执行半圆形扫描几何以获得多行的投影数据,从而在扫描平面内执行该物体的锥束扫描;
(b)利用加权因子对该投影数据加权以产生加权投影数据,其中随着该行远离该扫描平面,该加权因子减小;以及
(c)根据该加权投影数据形成该图像。
2.如权利要求1所述的方法,其中步骤(c)包括滤波该加权投影数据。
3.如权利要求1所述的方法,其中该加权因子是加权角的函数,该加权角是该行和扫描角的函数。
4.如权利要求3所述的方法,其中该加权角是
其中β是扫描角,m是行数,ξ是该检测器的虚拟检测采样间隔,so是从该发射器到该扫描几何的起点的距离。
5.如权利要求4所述的方法,其中该加权因子是
MO是该检测器的一半宽度。
6.如权利要求1所述的方法,还包括在步骤(b)和(c)之间,向该加权投影数据中添加表示通过圆形扫描不能采集到的Radon空间中的信息的校正项。
7.如权利要求1所述的方法,其中步骤(a)、(b)和(c)实时执行。
8.一种用于形成物体的图像的方法,该方法包括:
(a)通过使用锥束发射器和检测器执行圆形或部分圆形加螺旋扫描几何以获得投影数据,从而在扫描平面内执行该物体的锥束扫描,其中该扫描几何包括圆形或部分圆形扫描几何和螺旋扫描几何,其中该投影数据包括仅在沿着该螺旋扫描几何的离散点获取的投影数据;和
(b)根据该投影数据形成该图像。
9.如权利要求8所述的方法,其中该物体包括胸部。
10.如权利要求9所述的方法,其中该圆形或部分圆形扫描几何是围绕胸部的,并且该螺旋扫描几何是沿着胸部的。
11.如权利要求10所述的方法,其中步骤(b)包括添加来自该圆形或部分圆形扫描几何和螺旋扫描几何的项。
12.如权利要求11所述的方法,其中步骤(b)还包括添加表示通过圆形扫描不能采集到的Radon空间中的信息的校正项。
13.一种用于形成物体图像的***,该***包括:
锥束发射器;
检测器;
托架,用于引起该发射器和检测器的相对移动以通过使用锥束发射器和检测器执行半圆形扫描几何以获得多行的投影数据,从而在扫描平面内执行该物体的锥束扫描;和
处理器,接收该投影数据,利用加权因子对该投影数据加权以产生加权投影数据,其中随着该行远离该扫描平面,该加权因子减小,以及根据该加权投影数据形成该图像。
14.如权利要求13所述的***,其中该处理器滤波该加权投影数据。
15.如权利要求13所述的***,其中该加权因子是加权角的函数,该加权角是该行和扫描角的函数。
16.如权利要求15所述的***,其中该加权角是
其中β是扫描角,m是行数,ξ是该检测器的虚拟检测采样间隔,so是从该发射器到该扫描几何的起点的距离。
17.如权利要求16所述的***,其中该加权因子是
MO是该检测器的一半宽度。
18.如权利要求13所述的***,其中该处理器向该加权投影数据中添加表示通过圆形扫描不能采集到的Radon空间中的信息的校正项。
19.如权利要求13所述的***,其中该发射器、检测器、托架和处理器实时运行。
20.一种用于形成物体图像的***,该***包括:
锥束发射器;
检测器;
托架,用于引起该发射器和检测器的相对移动以通过使用锥束发射器和检测器执行圆形或部分圆形加螺旋扫描几何以获得投影数据,从而在扫描平面内执行该物体的锥束扫描,其中该扫描几何包括圆形或部分圆形扫描几何和螺旋扫描几何,其中该投影数据包括仅在沿着该螺旋扫描几何的离散点获取的投影数据;和
处理器,接收该投影数据,用于根据该投影数据形成该图像。
21.如权利要求20所述的***,其中该处理器添加来自该圆形或部分圆形扫描几何和螺旋扫描几何的项。
22.如权利要求21所述的***,其中该处理器还添加表示通过圆形扫描不能采集到的Radon空间中的信息的校正项。
23.如权利要求20所述的***,其中该托架围绕旋转轴旋转,以使得该圆形或部分圆形扫描几何以及二者都围绕该旋转轴旋转并且沿着该旋转轴移动以实现螺旋扫描几何。
24.一种用于产生物体的锥束CT图像的设备,该设备包括:
托架框架;
至少一个用于移动该托架框架以形成数据采集几何的电机;
附着到该托架框架以随着该托架框架移动的锥束辐射源;
附着到该托架框架以随着该托架框架移动的二维检测器,该二维检测器被设置在该辐射路径上;
支架,当获取该物体的投影图像时该物体支撑在其上,该支架支撑该物体以使得该物体被设置在该锥束辐射源和二维检测器之间;和
至少一个计算机***,用于控制物体扫描、执行该物体重建和图像的图像分析;
其特征在于:
该辐射源是锥束辐射源;
该获取的图像是该物体的锥束CT投影图像;
该至少一个电机移动该托架框架以便能够通过使该源和检测器围绕经过该物体的轴同步旋转,以形成用于锥束体积计算机断层扫描的数据采集几何,从而实现该物体的体积扫描,该体积扫描产生图像信号;
该至少一个电机上下移动该托架框架或支架,同时该源和检测器围绕经过该物体的轴同步移动,以形成用于锥束体积计算机断层扫描的数据采集几何,该体积扫描产生图像信号;
通过对该图像信号执行锥束计算机断层扫描重建以产生该物体的三维衰减系数分布,从而根据该图像信号形成该三维锥束CT图像;和
该至少一个电机包括用于移动该源和二维检测器以限定数据采集几何的电机。
25.如权利要求24所述的设备,其特征还在于该数据采集几何是定义为180°加上锥角圆形扫描的半扫描。
26.如权利要求24所述的设备,其特征还在于该数据采集几何是圆形扫描几何。
27.如权利要求24所述的设备,其特征还在于该数据采集几何是螺旋几何。
28.如权利要求24所述的设备,其特征还在于该数据采集几何是圆形加线形几何。
29.如权利要求28所述的设备,其特征还在于该圆形加线形几何包括单条线。
30.如权利要求28所述的设备,其特征还在于圆形加线形几何包括多条线。
31.如权利要求24所述的设备,其特征还在于该数据采集几何是360°扫描,从而可以通过选择不同的起始点而生成用于多个半扫描重建的多组半扫描投影图像。
32.如权利要求24所述的设备,其特征还在于该数据采集几何是N×360°的扫描,其中N是正整数。
33.一种用于产生物体的锥束CT图像的设备,该设备包括:
托架框架;
至少一个用于移动该托架框架以形成数据采集几何的电机;
附着到该托架框架以随着该托架框架移动的辐射源;
附着到该托架框架以随着该托架框架移动的二维检测器,该二维检测器被设置在该辐射路径上;
支架,当获取该物体的投影图像时该物体支撑在其上,该支架支撑该物体以使得该物体被设置在该辐射源和二维检测器之间;和
至少一个计算机***,用于控制胸部扫描、执行该胸部重建和图象的图像分析;
其特征在于:
该辐射源是锥束辐射源;
该获取的图像是锥束CT投影图像;
该至少一个电机移动该托架框架以便能够通过使该源和检测器围绕经过该物体的轴同步旋转,以形成用于锥束体积计算机断层扫描的数据采集几何,从而实现该物体的体积扫描,该体积扫描产生图像信号;
该至少一个电机上下移动该托架框架或支架,同时该源和检测器围绕经过该物体的轴同步移动,以形成用于锥束体积计算机断层扫描的数据采集几何,该体积扫描产生图像信号;
通过对该图像信号执行锥束体积计算机断层扫描重建以产生该胸部的三维衰减系数分布,从而根据该图像信号形成该三维锥束体积CT图像;
该至少一个电机移动该托架框架以使得该二维检测器对该物体进行体积扫描;
该设备还包括该托架框架上的滑环,用于提供电源到托架上部件的连接和该二维检测器与计算机***之间的通信;和
该至少一个电机移动该托架框架以使得该二维检测器对该物体进行体积扫描。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410048385.9A CN103908282B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
CN201510962389.2A CN105596022B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US77666606P | 2006-02-27 | 2006-02-27 | |
US60/776,666 | 2006-02-27 | ||
PCT/US2007/005080 WO2007098284A2 (en) | 2006-02-27 | 2007-02-27 | Method and apparatus for cone beam ct dynamic imaging |
Related Child Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510962389.2A Division CN105596022B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
CN201410048385.9A Division CN103908282B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101594825A true CN101594825A (zh) | 2009-12-02 |
CN101594825B CN101594825B (zh) | 2014-07-23 |
Family
ID=38438027
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410048385.9A Expired - Fee Related CN103908282B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
CN200780013865.XA Active CN101594825B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
CN201510962389.2A Active CN105596022B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410048385.9A Expired - Fee Related CN103908282B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510962389.2A Active CN105596022B (zh) | 2006-02-27 | 2007-02-27 | 锥束ct动态成像的方法和设备 |
Country Status (6)
Country | Link |
---|---|
US (2) | US7643605B2 (zh) |
EP (1) | EP1988826A4 (zh) |
CN (3) | CN103908282B (zh) |
AU (1) | AU2007217369A1 (zh) |
CA (1) | CA2643931A1 (zh) |
WO (1) | WO2007098284A2 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887591A (zh) * | 2010-06-24 | 2010-11-17 | 西北工业大学 | 基于矩形包围盒的锥束ct快速重建方法 |
CN102232844A (zh) * | 2010-04-23 | 2011-11-09 | 东芝It·控制***株式会社 | Ct装置和ct装置的摄影方法 |
CN103285518A (zh) * | 2012-01-31 | 2013-09-11 | 西门子公司 | 放射治疗中用于记录锥束断层照片的***和方法 |
WO2013134949A1 (zh) * | 2012-03-16 | 2013-09-19 | 西安电子科技大学 | 内窥式x射线发光断层成像装置及方法 |
CN104809750A (zh) * | 2015-05-04 | 2015-07-29 | 重庆大学 | 一种直线扫描ct***及图像重建方法 |
CN109620273A (zh) * | 2018-12-10 | 2019-04-16 | 合肥中科离子医学技术装备有限公司 | 一种实时计算短扫描权重的快速cbct重建算法 |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6987831B2 (en) * | 1999-11-18 | 2006-01-17 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
CA2643931A1 (en) * | 2006-02-27 | 2007-08-30 | University Of Rochester | Method and apparatus for cone beam ct dynamic imaging |
EP2005394B1 (en) * | 2006-04-06 | 2018-01-17 | Philips Intellectual Property & Standards GmbH | Method for reconstruction images and reconstruction system for reconstructing images |
EP1986551A1 (en) * | 2006-04-19 | 2008-11-05 | Philips Intellectual Property & Standards GmbH | Cone beam computed tomography with multiple partial scanning trajectories |
US8086010B2 (en) * | 2006-06-30 | 2011-12-27 | Kabushiki Kaisha Toshiba | Medical image diagnosis apparatus and the control method thereof |
GB2465726A (en) | 2007-08-23 | 2010-06-02 | Fischer Medical Technologies Inc | Improved computed tomography breast imaging and biopsy system |
SE0702258L (sv) * | 2007-10-09 | 2009-03-31 | Xcounter Ab | Anordning och metod för att upptaga strålningsbilddata av ett objekt |
US8023767B1 (en) | 2008-03-10 | 2011-09-20 | University Of Rochester | Method and apparatus for 3D metal and high-density artifact correction for cone-beam and fan-beam CT imaging |
US8401144B2 (en) * | 2008-08-07 | 2013-03-19 | Koninklijke Philips Electronics N.V. | Method and apparatus for correcting artifacts in circular CT scans |
US8284892B2 (en) * | 2008-12-22 | 2012-10-09 | General Electric Company | System and method for image reconstruction |
CN102456227B (zh) * | 2010-10-28 | 2015-05-27 | 清华大学 | Ct图像重建方法及装置 |
US8989469B2 (en) * | 2010-12-20 | 2015-03-24 | The Board Of Trustees Of The Leland Stanford Junior University | Systems and methods for simultaneous acquisition of scatter and image projection data in computed tomography |
US8805037B2 (en) | 2011-05-31 | 2014-08-12 | General Electric Company | Method and system for reconstruction of tomographic images |
US8615121B2 (en) | 2011-05-31 | 2013-12-24 | General Electric Company | Reconstruction of projection data to generate tomographic images having improved frequency characteristics |
US8818065B2 (en) * | 2011-07-01 | 2014-08-26 | Carestream Health, Inc. | Methods and apparatus for scatter correction for CBCT system and cone-beam image reconstruction |
US8861829B2 (en) | 2011-09-30 | 2014-10-14 | General Electric Company | Method and system for reconstruction of tomographic images |
US8731269B2 (en) * | 2011-10-19 | 2014-05-20 | Kabushiki Kaisha Toshiba | Method and system for substantially reducing artifacts in circular cone beam computer tomography (CT) |
DE102013200329B4 (de) * | 2013-01-11 | 2014-11-13 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Dejustagekorrektur für Bildgebungsverfahren |
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 |
US8948337B2 (en) | 2013-03-11 | 2015-02-03 | General Electric Company | Computed tomography image reconstruction |
WO2017102607A1 (en) * | 2015-12-18 | 2017-06-22 | Koninklijke Philips N.V. | Tomographic imaging device and method for sparse angular sampling |
EP3254623B1 (en) * | 2016-06-09 | 2020-09-16 | Agfa Nv | Method and system for correcting geometric misalignment during image reconstruction in chest tomosynthesis |
US11219425B2 (en) * | 2016-12-21 | 2022-01-11 | Koninklijke Philips N.V. | Redundancy-weighted image reconstruction for short-scan X-ray tomography with an off-center X-ray detector |
CN109146800B (zh) * | 2018-07-23 | 2019-08-13 | 广州华端科技有限公司 | 锥束计算机断层成像图像校正方法和*** |
CN109360252B (zh) * | 2018-09-13 | 2020-08-14 | 北京航空航天大学 | 基于投影变换的锥束cl投影数据等效转换方法 |
CN116758183B (zh) * | 2023-08-18 | 2023-11-07 | 有方(合肥)医疗科技有限公司 | 一种cbct图像重建方法 |
CN117838169B (zh) * | 2024-03-08 | 2024-05-24 | 江苏一影医疗设备有限公司 | 一种基于站立位cbct的成像方法、***和设备 |
Family Cites Families (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3165630A (en) * | 1961-06-13 | 1965-01-12 | Polaroid Corp | Table for holding and positioning a female subject and film during breast x-ray exposures |
US3973126A (en) * | 1975-07-31 | 1976-08-03 | General Electric Company | Mammography |
US4015836A (en) * | 1975-07-31 | 1977-04-05 | General Electric Company | Mammography table |
US5270923A (en) * | 1989-11-02 | 1993-12-14 | General Electric Company | Computed tomographic image reconstruction method for helical scanning using interpolation of partial scans for image construction |
US5409497A (en) | 1991-03-11 | 1995-04-25 | Fischer Imaging Corporation | Orbital aiming device for mammo biopsy |
US5170439A (en) * | 1991-06-11 | 1992-12-08 | Picker International, Inc. | Cone beam reconstruction using combined circle and line orbits |
US5289520A (en) * | 1991-11-27 | 1994-02-22 | Lorad Corporation | Stereotactic mammography imaging system with prone position examination table and CCD camera |
US5375156A (en) * | 1992-03-31 | 1994-12-20 | Siemens Medical Systems, Inc. | Method and apparatus for 3-D computer tomography |
US5564438A (en) * | 1993-08-09 | 1996-10-15 | Merchant; Thomas E. | Method and apparatus for prone position radiation therapy of the breast |
US5802137A (en) * | 1993-08-16 | 1998-09-01 | Commonwealth Scientific And Industrial Research | X-ray optics, especially for phase contrast imaging |
US5463666A (en) * | 1993-11-12 | 1995-10-31 | General Electric Company | Helical and circle scan region of interest computerized tomography |
US5526394A (en) * | 1993-11-26 | 1996-06-11 | Fischer Imaging Corporation | Digital scan mammography apparatus |
US5448607A (en) * | 1994-02-08 | 1995-09-05 | Analogic Corporation | X-ray tomography system with gantry pivot and translation control |
JPH07303633A (ja) * | 1994-05-11 | 1995-11-21 | Mitsubishi Electric Corp | X線***撮影装置 |
US5459769A (en) * | 1994-11-09 | 1995-10-17 | General Electric Company | Procedure for monitoring contrast agent application in a CT imaging system |
AUPN201295A0 (en) * | 1995-03-28 | 1995-04-27 | Commonwealth Scientific And Industrial Research Organisation | Simplified conditions and configurations for phase-contrast imaging with hard x-rays |
US6002738A (en) * | 1995-07-07 | 1999-12-14 | Silicon Graphics, Inc. | System and method of performing tomographic reconstruction and volume rendering using texture mapping |
US5602891A (en) * | 1995-11-13 | 1997-02-11 | Beth Israel | Imaging apparatus and method with compensation for object motion |
US5784481A (en) * | 1996-06-25 | 1998-07-21 | General Electric Company | CT cone beam image reconstruction with circle and line scan path |
CA2227531C (en) * | 1997-01-20 | 2003-03-18 | Hitachi, Ltd. | Graphics processing unit and graphics processing system |
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 |
US6014419A (en) * | 1997-11-07 | 2000-01-11 | Hu; Hui | CT cone beam scanner with fast and complete data acquistion and accurate and efficient regional reconstruction |
US6047042A (en) * | 1998-03-25 | 2000-04-04 | Continental X-Ray Corporation | Automatic exposure and brightness control for fluoroscopic and radio-graphic imaging |
DE19843812A1 (de) * | 1998-09-24 | 2000-03-30 | Philips Corp Intellectual Pty | Computertomographie-Verfahren mit kegelförmigem Strahlenbündel |
US6987831B2 (en) * | 1999-11-18 | 2006-01-17 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
US6480565B1 (en) * | 1999-11-18 | 2002-11-12 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
US6324247B1 (en) * | 1999-12-30 | 2001-11-27 | Ge Medical Systems Global Technology Company, Llc | Partial scan weighting for multislice CT imaging with arbitrary pitch |
US6327325B1 (en) * | 2000-02-16 | 2001-12-04 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for adaptive interpolation reduced view CT scan |
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 |
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 |
US6778630B2 (en) * | 2001-03-23 | 2004-08-17 | Kabushiki Kaisha Toshiba | Method and system for reconstructing computed tomography images using redundant data |
US6678346B2 (en) * | 2001-10-11 | 2004-01-13 | Ge Medical Systems Global Technology Company Llc | Cone-beam CT scanner with image reconstruction using multiple sub-images |
US6618466B1 (en) * | 2002-02-21 | 2003-09-09 | University Of Rochester | Apparatus and method for x-ray scatter reduction and correction for fan beam CT and cone beam volume CT |
US7245755B1 (en) * | 2002-07-10 | 2007-07-17 | Xiaochuan Pan | Algorithm for image reconstruction and image noise analysis in computed tomography |
FR2847798B1 (fr) * | 2002-11-28 | 2006-02-10 | Ge Med Sys Global Tech Co Llc | Procede pour determiner des parametres fonctionnels dans un dispositif de fluoroscopie |
WO2004058070A1 (ja) * | 2002-12-26 | 2004-07-15 | Atsushi Momose | X線撮像装置および撮像方法 |
US6917663B2 (en) * | 2003-06-16 | 2005-07-12 | Kabushiki Kaisha Toshiba | Cone-beam reconstruction apparatus and computed tomography apparatus |
EP1774475A1 (en) * | 2004-06-28 | 2007-04-18 | Koninklijke Philips Electronics N.V. | Examination apparatus for perfusion studies |
US7412026B2 (en) * | 2004-07-02 | 2008-08-12 | The Board Of Regents Of The University Of Oklahoma | Phase-contrast x-ray imaging systems and methods |
WO2006085253A2 (en) * | 2005-02-11 | 2006-08-17 | Philips Intellectual Property & Standards Gmbh | Computer tomography apparatus, method of examining an object of interest with a computer tomography apparatus, computer-readable medium and program element |
CA2605836A1 (en) * | 2005-04-25 | 2006-11-02 | University Of Rochester | Method and apparatus of global de-noising for ct imaging |
CN100457040C (zh) | 2005-11-17 | 2009-02-04 | 中国科学院高能物理研究所 | 同步辐射x射线相位衬度ct成像装置及实验方法 |
DE102006046034A1 (de) * | 2006-02-01 | 2007-08-16 | Siemens Ag | Röntgen-CT-System zur Erzeugung projektiver und tomographischer Phasenkontrastaufnahmen |
CN101011257B (zh) | 2006-02-01 | 2011-07-06 | 西门子公司 | 产生投影或断层造影相位对比图像的焦点-检测器装置 |
CA2643931A1 (en) * | 2006-02-27 | 2007-08-30 | University Of Rochester | Method and apparatus for cone beam ct dynamic imaging |
-
2007
- 2007-02-27 CA CA002643931A patent/CA2643931A1/en not_active Abandoned
- 2007-02-27 CN CN201410048385.9A patent/CN103908282B/zh not_active Expired - Fee Related
- 2007-02-27 CN CN200780013865.XA patent/CN101594825B/zh active Active
- 2007-02-27 EP EP07751813A patent/EP1988826A4/en not_active Withdrawn
- 2007-02-27 CN CN201510962389.2A patent/CN105596022B/zh active Active
- 2007-02-27 US US11/711,155 patent/US7643605B2/en active Active
- 2007-02-27 AU AU2007217369A patent/AU2007217369A1/en not_active Abandoned
- 2007-02-27 WO PCT/US2007/005080 patent/WO2007098284A2/en active Application Filing
-
2009
- 2009-11-13 US US12/618,357 patent/US8976925B2/en active Active
Non-Patent Citations (2)
Title |
---|
DONG YANG,ETC: "Implementation & Evaluation of the Half-Scan Scheme Based on", 《MEDICAL IMAGING 2004》 * |
DONG YANG,ETC: "Modified FDK half-scan (MFDKHS) scheme on flat panel detectorbased", 《MEDICAL IMAGING 2005》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102232844A (zh) * | 2010-04-23 | 2011-11-09 | 东芝It·控制***株式会社 | Ct装置和ct装置的摄影方法 |
CN102232844B (zh) * | 2010-04-23 | 2013-07-03 | 东芝It·控制***株式会社 | Ct装置和ct装置的摄影方法 |
CN101887591A (zh) * | 2010-06-24 | 2010-11-17 | 西北工业大学 | 基于矩形包围盒的锥束ct快速重建方法 |
CN103285518A (zh) * | 2012-01-31 | 2013-09-11 | 西门子公司 | 放射治疗中用于记录锥束断层照片的***和方法 |
CN103285518B (zh) * | 2012-01-31 | 2016-07-06 | 西门子公司 | 放射治疗中用于记录锥束断层照片的***和方法 |
WO2013134949A1 (zh) * | 2012-03-16 | 2013-09-19 | 西安电子科技大学 | 内窥式x射线发光断层成像装置及方法 |
CN104168830A (zh) * | 2012-03-16 | 2014-11-26 | 西安电子科技大学 | 内窥式x射线发光断层成像装置及方法 |
CN104168830B (zh) * | 2012-03-16 | 2016-06-22 | 西安电子科技大学 | 内窥式x射线发光断层成像装置及方法 |
US10034647B2 (en) | 2012-03-16 | 2018-07-31 | Xidian University | System and method for endoscopic X-ray luminescence computed tomographic imaging |
CN104809750A (zh) * | 2015-05-04 | 2015-07-29 | 重庆大学 | 一种直线扫描ct***及图像重建方法 |
CN104809750B (zh) * | 2015-05-04 | 2017-10-24 | 重庆大学 | 一种直线扫描ct***及图像重建方法 |
CN109620273A (zh) * | 2018-12-10 | 2019-04-16 | 合肥中科离子医学技术装备有限公司 | 一种实时计算短扫描权重的快速cbct重建算法 |
Also Published As
Publication number | Publication date |
---|---|
US8976925B2 (en) | 2015-03-10 |
US20100158335A1 (en) | 2010-06-24 |
WO2007098284B1 (en) | 2008-11-27 |
EP1988826A4 (en) | 2010-05-19 |
WO2007098284A3 (en) | 2008-10-16 |
EP1988826A2 (en) | 2008-11-12 |
US20070253528A1 (en) | 2007-11-01 |
CN105596022A (zh) | 2016-05-25 |
AU2007217369A1 (en) | 2007-08-30 |
CN103908282B (zh) | 2016-05-18 |
CN105596022B (zh) | 2020-05-22 |
CN101594825B (zh) | 2014-07-23 |
CA2643931A1 (en) | 2007-08-30 |
CN103908282A (zh) | 2014-07-09 |
WO2007098284A2 (en) | 2007-08-30 |
US7643605B2 (en) | 2010-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101594825B (zh) | 锥束ct动态成像的方法和设备 | |
Noo et al. | Single-slice rebinning method for helical cone-beam CT | |
CN102596040B (zh) | 计算机断层摄影设备 | |
US7333588B2 (en) | Virtual spherical anode computed tomography | |
CN103384498B (zh) | 探测装置 | |
CN102144928B (zh) | 利用多重x射线源的ct测量 | |
EP2512345B1 (en) | Computed tomography apparatus | |
US9025838B2 (en) | Apparatus and method for hybrid reconstruction of an object from projection data | |
JP2009529394A (ja) | ラドンデータから画像関数を再構成する方法 | |
CN102293660B (zh) | 在心脏ct中改善时间分辨率 | |
CN101036164A (zh) | 采用与体素相关插值的图像重建 | |
US7529335B2 (en) | Voxel-driven spiral reconstruction for cone-beam computer tomography | |
JP2001057976A (ja) | 立体画像再構成方法及び装置並びにctスキャナー | |
CN104644198A (zh) | 用于对对象进行成像的成像*** | |
Yu et al. | Local ROI reconstruction via generalized FBP and BPF algorithms along more flexible curves | |
US20090116721A1 (en) | Method and system for multi-detector-row computed tomography utilizing projection data up-sampling with shift | |
Hsieh | Tomographic reconstruction for tilted helical multislice CT | |
Shechter et al. | The frequency split method for helical cone‐beam reconstruction | |
Kim et al. | Evaluation of algebraic iterative image reconstruction methods for tetrahedron beam computed tomography systems | |
Zeng et al. | Cone-beam iterative reconstruction of a segment of a long object | |
JP3950101B2 (ja) | X線コンピュータ断層撮影装置 | |
Baek et al. | Initial results with a multisource inverse-geometry CT system | |
CN106462983A (zh) | 计算机化层析图像曝光和重构方法 | |
Ashfaq | Data-Driven Iterative Image Reconstruction For Motion Correction In Head Computed Tomography | |
Sun | Rigid motion correction for head CT imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Ning Ruola Inventor after: D *yang Inventor before: R *ning Inventor before: D *yang |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: R. NING D. YANG TO: NING RUOLA D. YANG |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |