CN105093342B - 螺旋ct系统及重建方法 - Google Patents
螺旋ct系统及重建方法 Download PDFInfo
- Publication number
- CN105093342B CN105093342B CN201410202941.3A CN201410202941A CN105093342B CN 105093342 B CN105093342 B CN 105093342B CN 201410202941 A CN201410202941 A CN 201410202941A CN 105093342 B CN105093342 B CN 105093342B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- data
- projection
- detector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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
-
- 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 for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4266—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a plurality of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
-
- 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]
Abstract
公开了一种螺旋CT系统及其重建方法。在一些实施例中提出了用螺旋CT系统获得的投影数据的互补投影数据的加权来补偿因为大螺距造成的数据缺失。在将数据补全后,将投影数据重排为带锥角的平行束数据,进而进行锥角余弦加权和一维滤波,最后进行平行束反投影来得到重建的图像。在一些实施例中,利用上述的方法,可以在现有探测器面积和滑环速度不变的情况下,将皮带速度提高一倍以上,从而提高行李通过率,并且保持重建的图像质量不变。
Description
技术领域
本发明的实施例一般涉及辐射成像,具体涉及大螺距CT系统及其重建方法。
背景技术
CT型行李安检系统已经成为一种重要的爆炸物检测手段,广泛应用于机场、车站等公共场所。CT检测系统的原理是通过扫描获取行李的投影数据,利用重建算法获得其断层图像,根据识别算法对爆炸物进行识别和报警。这种扫描的方式通常采用安装了X光源和探测器的滑环旋转同时物体通过皮带前进的方式,因此扫描轨道属于螺旋轨道。为了满足安检场所对于行李通过率的要求,皮带前进速度必须达到某一限值,例如0.3m/s~0.5m/s。为了采集重建图像所需要的完备的投影数据,需要提高滑环转速或者增加探测器排数。对于滑环,考虑到机械强度和稳定性,其转速不能无限提高;另一方面,由于X光源锥角的限值和探测器硬件成本考虑,也不能无限增加探测器排数。综合以上各方面因素,只有增加扫描螺距是提高行李通过率的最有效的方式。
近年来,在螺旋CT重建算法领域出现了一系列重要成果。例如:Katsevich算法、PI算法、CB-FBP算法等。但是,这些算法都需要满足一定的螺距条件,当螺距和锥角增大时会导致投影数据缺失,重建结果的误差也会增大,甚至引入伪影。为了保证图像质量满足要求,通常螺距因子不大于1.5。
在安检CT系统中,重建算法的实时性也是一项重要的指标,在重建算法中,非一维移不变形式的滤波、具有距离权重因子的锥束反投影、求解大量非线性方程、使用较大的反投影角度范围都会导致重建算法效率降低,应尽量避免使用这一类的算法。因此,上述的一些重建算法都具有影响算法效率的环节,并且也主要用于投影数据完备或者冗余的情况。
综上所述,目前已有技术均未考虑大螺距条件下的数据缺失问题,因此无法直接用于大螺距CT系统。
发明内容
鉴于现有技术中的一个或多个问题,提出了一种螺旋CT系统及其重建方法,能够满足大螺距情况下的图像重建要求。
在本发明的一个方面,提出了一种锥束螺旋CT的重建方法,包括步骤:根据锥束螺旋CT系统的螺距和多排探测器的排间距计算覆盖Tam窗所需要的最少探测器排数;在所述锥束螺旋CT系统的探测器排数小于所述最少探测器排数的情况下,通过对互补的投影数据进行加权来弥补缺失的投影数据;将补全后的投影数据重排为带锥角的平行束数据;对重排后的平行束数据进行锥角余弦加权,然后沿着重排为平行束数据时定义的虚拟探测器的排方向进行一维滤波;以及对滤波后的数据进行无加权的带锥角的平行束反投影,得到重建图像。
在本发明的另一方面,提出了一种锥束螺旋CT系统,包括:根据锥束螺旋CT系统的螺距和多排探测器的排间距计算覆盖Tam窗所需要的最少探测器排数的装置;在所述锥束螺旋CT系统的探测器排数小于所述最少探测器排数的情况下,通过对互补的投影数据进行加权来弥补缺失的投影数据的装置;将补全后的投影数据重排为带锥角的平行束数据的装置;对重排后的平行束数据进行锥角余弦加权,然后沿着重排为平行束数据时定义的虚拟探测器的排方向进行一维滤波的装置;以及对滤波后的数据进行无加权的带锥角的平行束反投影,得到重建图像的装置。
在一些实施例中,利用上述的方案,可以在现有探测器面积和滑环速度不变的情况下,将皮带速度提高一倍以上,从而提高行李通过率,并且保持重建的图像质量不变。
附图说明
为了更好地理解本发明,将根据以下附图对本发明进行详细描述:
图1示出了根据本发明实施例的CT系统中螺旋轨道扫描示意图,其中探测器上的灰色区域表示Tam窗;
图2示出了根据本发明实施例的扇束扫描中的共轭射线对定义;
图3是根据本发明实施方式的CT设备的结构示意图;
图4示出了如图3所示的计算机数据处理器的结构框图;
图5示出了根据本发明实施方式的控制器的结构框图;
图6示出了根据本发明实施例的重建方法的流程图;
图7示出了根据本发明实施例的螺旋轨道扫描下的共轭投影插值示意图
图8示出了重排成带锥角的平行束及虚拟探测器的示意图;
图9示出了在本发明的一个具体实施例中,螺距为3.9cm时的重建结果;以及
图10示出了在本本发明的一个具体实施例中,螺距为6cm时的重建结果。
具体实施方式
下面将详细描述本发明的具体实施例,应当注意,这里描述的实施例只用于举例说明,并不用于限制本发明。在以下描述中,为了提供对本发明的透彻理解,阐述了大量特定细节。然而,对于本领域普通技术人员显而易见的是:不必采用这些特定细节来实行本发明。在其他实例中,为了避免混淆本发明,未具体描述公知的结构、材料或方法。
在整个说明书中,对“一个实施例”、“实施例”、“一个示例”或“示例”的提及意味着:结合该实施例或示例描述的特定特征、结构或特性被包含在本发明至少一个实施例中。因此,在整个说明书的各个地方出现的短语“在一个实施例中”、“在实施例中”、“一个示例”或“示例”不一定都指同一实施例或示例。此外,可以以任何适当的组合和/或子组合将特定的特征、结构或特性组合在一个或多个实施例或示例中。此外,本领域普通技术人员应当理解,这里使用的术语“和/或”包括一个或多个相关列出的项目的任何和所有组合。
针对现有技术无法满足大螺距情况下的重建要求,尤其是螺距因子大于1.5时,本发明的实施例提出了用螺旋CT系统获得的投影数据的互补投影数据的加权来补偿因为大螺距造成的数据缺失。在将数据补全后,将投影数据重排为带锥角的平行束数据,进而进行锥角余弦加权和一维滤波,最后进行平行束反投影来得到重建的图像。在一些实施例中,利用上述的方法,可以在现有探测器面积和滑环速度不变的情况下,将皮带速度提高一倍以上,从而提高行李通过率,并且保持重建的图像质量不变。
图1示出了根据本发明实施例的CT系统中螺旋轨道扫描示意图。如图1所示,定义X光源110的旋转半径为R,滑环旋转1圈皮带前进的距离(称为螺距)为h,柱面等角探测器120的扇束张角为2αm,X光源110到探测器120的距离为D。探测器120的排数为Nrow,排间距为srow。定义螺距因子为:
滑环的圆周旋转和皮带的平移形成了相对运动的螺旋轨道131,在以扫描物体的坐标系中,X光源的运动轨迹可以表示为:
在螺旋扫描轨道下,精确重建所需要的投影数据为相距最近的两段螺旋线在探测器上的投影即Tam-Danielsson窗(以下简称Tam窗)121所覆盖的投影数据,如图1所示。由此,可以计算出在排间距srow和螺距h等参数确定的条件下,覆盖Tam窗所需要的最少探测器排数为:
根据该公式(3)也可以反推出对于一个固定了探测器排数、排间距和扇角的CT系统,其允许的最大螺距和最大螺距因子为:
例如,当扇角2αm=π/3时,最大螺距因子为1.3。
图2示出了根据本发明实施例的扇束扫描中的共轭射线对定义。如图2所示,定义X光源210的旋转半径为R,滑环旋转1圈皮带前进的距离(称为螺距)为h,柱面等角探测器220的扇束张角为2αm,X光源210到探测器220的距离为D。在二维扇束扫描下,穿过同一条直线的两条射线称为共轭射线对L-L',即从从射线源210到位置210'的射线和从射线源位于位置210'时到位置210的射线,如图2所示。在锥束螺旋扫描下,不存在严格意义的共轭射线对,但可以定义具有类似关系的互补投影:
这样,螺旋重建算法在反投影时,可以利用互补的投影数据进行冗余加权或者补偿缺失的数据。
图3是根据本发明实施方式的CT设备的结构示意图。如图3所示,根据本实施方式的CT设备包括:机架20、承载机构40、控制器50、计算机数据处理器60等。机架20包括发出检查用X射线的射线源10,诸如X光机,以及探测和采集装置30。承载机构40承载被检查行李70穿过机架20的射线源10与探测和采集装置30之间的扫描区域,同时机架20围绕被检查行李70的前进方向转动,从而由射线源10发出的锥束射线能够透过被检查行李70,对被检查行李70进行CT扫描。
探测和采集装置30例如是具有整体模块结构的探测器及数据采集器,例如多排探测器,用于探测透射被检物品的射线,获得模拟信号,并且将模拟信号转换成数字信号,从而输出被检查行李70针对X射线的投影数据。控制器50用于控制整个系统的各个部分同步工作。计算机数据处理器60用来处理由数据采集器采集的数据,对数据进行处理并重建,输出结果。
如图3所示,射线源10置于可放置被检物体的一侧,探测和采集装置30置于被检查行李70的另一侧,包括探测器和数据采集器,用于获取被检查行李70的透射数据和/或多角度投影数据。数据采集器中包括数据放大成形电路,它可工作于(电流)积分方式或脉冲(计数)方式。探测和采集装置30的数据输出电缆与控制器50和计算机数据处理器60连接,根据触发命令将采集的数据存储在计算机数据处理器60中。
图4示出了如图3所示的计算机数据处理器60的结构框图。如图4所示,数据采集器所采集的数据通过接口单元68和总线64存储在存储器61中。只读存储器(ROM)62中存储有计算机数据处理器的配置信息以及程序。随机存取存储器(RAM)63用于在处理器66工作过程中暂存各种数据。另外,存储器61中还存储有用于进行数据处理的计算机程序。内部总线64连接上述的存储器61、只读存储器62、随机存取存储器63、输入装置65、处理器66、显示装置67和接口单元68。
在用户通过诸如键盘和鼠标之类的输入装置65输入的操作命令后,计算机程序的指令代码命令处理器66执行预定的数据重建算法,在得到数据处理结果之后,将其显示在诸如LCD显示器之类的显示装置67上,或者直接以诸如打印之类硬拷贝的形式输出处理结果。
图5示出了根据本发明实施方式的控制器的结构框图。如图5所示,控制器50包括:控制单元51,根据来自计算机60的指令,来控制射线源10、承载机构40和探测和采集装置30;触发信号产生单元52,用于在控制单元的控制下产生用来触发射线源10、探测和采集装置30以及承载机构40的动作的触发命令;第一驱动设备53,它在根据触发信号产生单元52在控制单元51的控制下产生的触发命令驱动承载机构40传送被检查行李70;第二驱动设备54,它根据触发信号产生单元52在控制单元51的控制下产生的触发命令机架20旋转。探测和采集装置30获得的投影数据存储在计算机60中进行CT断层图像重建,从而获得被检查行李70的断层图像数据。根据其他实施例,上述的CT成像系统也可以是双能CT系统,也就是机架20的X射线源10能够发出高能和低能两种射线,探测和采集装置30探测到不同能量水平下的投影数据后,由计算机数据处理器60进行双能CT重建,得到被检查行李70的各个断层的等效原子序数和等效电子密度数据。
图6示出了根据本发明实施例的重建方法的流程图。如图6所示,在步骤S61,根据锥束螺旋CT系统的螺距和多排探测器的排间距计算覆盖Tam窗所需要的最少探测器排数。例如,根据螺距、排间距等参数计算覆盖Tam窗所需要的最少探测器排数
在步骤S62,在所述锥束螺旋CT系统的探测器排数小于所述最少探测器排数的情况下,通过对互补的投影数据进行加权来弥补缺失的投影数据。例如,如果实际探测器排数小于则利用互补投影补偿缺失的数据。如图7所示,分两种情况:
1)当s>smax时,
P(λ,α,s)=(1-w1)·P(λ,α,smax)
+w1·P(λ+π-2α,-α,smin) (6)
2)类似地,当s<smin时,
P(λ,α,s)=(1-w2)·P(λ,α,smin)
+w2·P(λ-π-2α,-α,smax) (7)
其中,s表示探测器排(锥角)方向坐标、α表示探测器列(扇角)方向坐标、λ表示投影角、smin表示探测器排(锥角)方向坐标最小值、smax表示探测器排(锥角)方向坐标最大值、R表示X光源旋转半径、h表示滑环旋转一圈皮带前进的距离即螺距、D表示X光源到柱面探测器的距离。
在步骤S63,将补全后的投影数据P(λ,α,s)重排为带锥角的平行束数据Pp(λ,t,s)。重排的虚拟探测器定义为过旋转中心平行于z轴的矩形区域,其宽度等于真实柱面探测器的扇角所覆盖的视野直径2Rsinαm,如图8所示。
在步骤S64,对重排后的数据首先进行锥角余弦加权,然后进行一维滤波,滤波方向沿着虚拟探测器的排方向。由于滤波过程为一维移不变形式,因此可以使用快速傅里叶变换实施滤波。滤波核选择为平行束滤波反投影算法中的R-L卷积核:
Pf(λ,t,s)=(Pp(λ,t,s)cosγ)*g(t) (8)
在步骤S65,进行无加权的带锥角的平行束反投影,得到重建图像:
根据本实施例的方法,可以在现有探测器面积、滑环转速等硬件条件不改变的情况下,将皮带速度提高一倍左右,进而提高行李通过率,而重建图像质量保持不变。另一方面,利用本专利所述方法,可以在特定使用场合下针对现有CT系统设计进行更改,适当降低探测器的排数,从而降低系统成本。
例如,X光源旋转半径为50cm,源到探测器距离为80cm,探测器扇角为60°,探测器排数为32,排间距为1.5mm。若需要满足采集完备投影数据的要求,则最大螺距因子为1.3,对应的最大螺距为3.9cm。若使用本专利所述方法,则最大螺距可以增加至6cm,螺距因子为2。
图9示出了在本发明的一个具体实施例中,螺距为3.9cm时的重建结果。图10示出了在本本发明的一个具体实施例中,螺距为6cm时的重建结果。从图中可见,本发明在提高了扫描速度的同时,保持了重建的图像质量。
根据一些实施例,可以在现有探测器、滑环等硬件条件下进一步提高皮带速度,进而提高通过率,可以用于爆炸物检测、毒品检测等安检领域。此外,所提出的使用互补投影补偿缺失数据的方法可以适用于投影数据不完备的情况,在螺距因子大于1.5时仍然保证重建图像的质量。同时该实施例所提出的重建算法具有移不变滤波形式,使用不带距离权重因子的平行束反投影,使用最少的180度反投影角度范围,不需要求解非线性方程组,相对于现有技术具有较低的复杂性,能够满足安检的实时性要求。
以上的详细描述通过使用示意图、流程图和/或示例,已经阐述了重建方法和螺旋CT系统的众多实施例。在这种示意图、流程图和/或示例包含一个或多个功能和/或操作的情况下,本领域技术人员应理解,这种示意图、流程图或示例中的每一功能和/或操作可以通过各种结构、硬件、软件、固件或实质上它们的任意组合来单独和/或共同实现。在一个实施例中,本发明的实施例所述主题的若干部分可以通过专用集成电路(ASIC)、现场可编程门阵列(FPGA)、数字信号处理器(DSP)、或其他集成格式来实现。然而,本领域技术人员应认识到,这里所公开的实施例的一些方面在整体上或部分地可以等同地实现在集成电路中,实现为在一台或多台计算机上运行的一个或多个计算机程序(例如,实现为在一台或多台计算机系统上运行的一个或多个程序),实现为在一个或多个处理器上运行的一个或多个程序(例如,实现为在一个或多个微处理器上运行的一个或多个程序),实现为固件,或者实质上实现为上述方式的任意组合,并且本领域技术人员根据本公开,将具备设计电路和/或写入软件和/或固件代码的能力。此外,本领域技术人员将认识到,本公开所述主题的机制能够作为多种形式的程序产品进行分发,并且无论实际用来执行分发的信号承载介质的具体类型如何,本公开所述主题的示例性实施例均适用。信号承载介质的示例包括但不限于:可记录型介质,如软盘、硬盘驱动器、紧致盘(CD)、数字通用盘(DVD)、数字磁带、计算机存储器等;以及传输型介质,如数字和/或模拟通信介质(例如,光纤光缆、波导、有线通信链路、无线通信链路等)。
虽然已参照几个典型实施例描述了本发明,但应当理解,所用的术语是说明和示例性、而非限制性的术语。由于本发明能够以多种形式具体实施而不脱离发明的精神或实质,所以应当理解,上述实施例不限于任何前述的细节,而应在随附权利要求所限定的精神和范围内广泛地解释,因此落入权利要求或其等效范围内的全部变化和改型都应为随附权利要求所涵盖。
Claims (5)
1.一种锥束螺旋CT的重建方法,包括步骤:
根据锥束螺旋CT系统的螺距和多排探测器的排间距计算覆盖Tam窗所需要的最少探测器排数;
在所述锥束螺旋CT系统的探测器排数小于所述最少探测器排数的情况下,通过对互补的投影数据进行加权来弥补缺失的投影数据;
将补全后的投影数据重排为带锥角的平行束数据;
对重排后的平行束数据进行锥角余弦加权,然后沿着重排为平行束数据时定义的虚拟探测器的排方向进行一维滤波;以及
对滤波后的数据进行无加权的带锥角的平行束反投影,得到重建图像。
2.如权利要求1所述的重建方法,其中通过对互补的投影数据进行加权来弥补缺失的投影数据的步骤包括:
当s>smax时,
P(λ,α,s)=(1-w1)·P(λ,α,smax)
+w1·P(λ+π-2α,-α,smin)
<mrow>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>h</mi>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mo>-</mo>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>R</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>s</mi>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
当s<smin时,
P(λ,α,s)=(1-w2)·P(λ,α,smin)
+w2·P(λ-π-2α,-α,smax)
<mrow>
<msub>
<mi>w</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>s</mi>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>h</mi>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mo>+</mo>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>R</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,P表示投影数据、s表示探测器排(锥角)方向坐标、α表示探测器列(扇角)方向坐标、λ表示投影角、smin表示探测器排(锥角)方向坐标最小值、smax表示探测器排(锥角)方向坐标最大值、R表示X光源旋转半径、h表示滑环旋转一圈皮带前进的距离即螺距、D表示X光源到柱面探测器的距离。
3.如权利要求1所述的重建方法,其中所述一维滤波使用R-L卷积核。
4.一种锥束螺旋CT系统,包括:
根据锥束螺旋CT系统的螺距和多排探测器的排间距计算覆盖Tam窗所需要的最少探测器排数的装置;
在所述锥束螺旋CT系统的探测器排数小于所述最少探测器排数的情况下,通过对互补的投影数据进行加权来弥补缺失的投影数据的装置;
将补全后的投影数据重排为带锥角的平行束数据的装置;
对重排后的平行束数据进行锥角余弦加权,然后沿着重排为平行束数据时定义的虚拟探测器的排方向进行一维滤波的装置;以及
对滤波后的数据进行无加权的带锥角的平行束反投影,得到重建图像的装置。
5.如权利要求1所述的系统,其中通过对互补的投影数据进行加权来弥补缺失的投影数据的装置包括通过以下方式来弥补缺失的投影数据的装置:
当s>smax时,
P(λ,α,s)=(1-w1)·P(λ,α,smax)
+w1·P(λ+π-2α,-α,smin)
<mrow>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>h</mi>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mo>-</mo>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>R</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>s</mi>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
当s<smin时,
P(λ,α,s)=(1-w2)·P(λ,α,smin)
+w2·P(λ-π-2α,-α,smax)
<mrow>
<msub>
<mi>w</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>s</mi>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>h</mi>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mo>+</mo>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>R</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&alpha;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,P表示投影数据、s表示探测器排(锥角)方向坐标、α表示探测器列(扇角)方向坐标、λ表示投影角、smin表示探测器排(锥角)方向坐标最小值、smax表示探测器排(锥角)方向坐标最大值、R表示X光源旋转半径、h表示滑环旋转一圈皮带前进的距离即螺距、D表示X光源到柱面探测器的距离。
Priority Applications (15)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410202941.3A CN105093342B (zh) | 2014-05-14 | 2014-05-14 | 螺旋ct系统及重建方法 |
US14/708,866 US9466137B2 (en) | 2014-05-14 | 2015-05-11 | Spiral CT systems and reconstruction methods |
EP15167283.9A EP2945125B1 (en) | 2014-05-14 | 2015-05-12 | Spiral ct systems and reconstruction methods |
ES15167283T ES2741577T3 (es) | 2014-05-14 | 2015-05-12 | Sistemas de CT de espiral y métodos de reconstrucción |
BR102015010870-2A BR102015010870B1 (pt) | 2014-05-14 | 2015-05-12 | Método de reconstrução de um sistema de tomografia computadorizada, e, sistema de tomografia computadorizada |
RU2015117682/28A RU2595312C1 (ru) | 2014-05-14 | 2015-05-12 | Системы спиральной кт и способы реконструкции |
AU2015202554A AU2015202554B2 (en) | 2014-05-14 | 2015-05-12 | Spiral CT systems and reconstruction methods |
JP2015098129A JP6177831B2 (ja) | 2014-05-14 | 2015-05-13 | 螺旋ctシステム及びその再構成方法 |
KR1020150066885A KR101704390B1 (ko) | 2014-05-14 | 2015-05-13 | 나선형 ct 시스템 및 재구성 방법 |
MX2016012501A MX358143B (es) | 2014-05-14 | 2015-05-14 | Sistema de tc en espiral y metodo de reconstruccion. |
MYPI2016703561A MY177410A (en) | 2014-05-14 | 2015-05-14 | Spiral ct systems and reconstruction methods |
SG11201608163TA SG11201608163TA (en) | 2014-05-14 | 2015-05-14 | Spiral ct system and reconstruction method |
PCT/CN2015/078956 WO2015172727A1 (zh) | 2014-05-14 | 2015-05-14 | 螺旋ct系统及重建方法 |
HK16105821.6A HK1217769A1 (zh) | 2014-05-14 | 2016-05-23 | 螺旋 系統及重建方法 |
IL247916A IL247916B (en) | 2014-05-14 | 2016-09-19 | Spiral computed tomography systems and reconstruction methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410202941.3A CN105093342B (zh) | 2014-05-14 | 2014-05-14 | 螺旋ct系统及重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093342A CN105093342A (zh) | 2015-11-25 |
CN105093342B true CN105093342B (zh) | 2017-11-17 |
Family
ID=53276717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410202941.3A Active CN105093342B (zh) | 2014-05-14 | 2014-05-14 | 螺旋ct系统及重建方法 |
Country Status (15)
Country | Link |
---|---|
US (1) | US9466137B2 (zh) |
EP (1) | EP2945125B1 (zh) |
JP (1) | JP6177831B2 (zh) |
KR (1) | KR101704390B1 (zh) |
CN (1) | CN105093342B (zh) |
AU (1) | AU2015202554B2 (zh) |
BR (1) | BR102015010870B1 (zh) |
ES (1) | ES2741577T3 (zh) |
HK (1) | HK1217769A1 (zh) |
IL (1) | IL247916B (zh) |
MX (1) | MX358143B (zh) |
MY (1) | MY177410A (zh) |
RU (1) | RU2595312C1 (zh) |
SG (1) | SG11201608163TA (zh) |
WO (1) | WO2015172727A1 (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10485498B2 (en) * | 2016-01-29 | 2019-11-26 | Koninklijke Philips N.V. | Cone beam computed tomography projection values providing system and method |
CN106097411B (zh) * | 2016-06-29 | 2019-07-05 | 沈阳开普医疗影像技术有限公司 | Ct机图像重建方法及高分辨ct扫描机 |
CN106526686B (zh) * | 2016-12-07 | 2019-05-07 | 同方威视技术股份有限公司 | 螺旋ct设备和三维图像重建方法 |
WO2019033390A1 (en) * | 2017-08-18 | 2019-02-21 | Shenzhen United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR IMAGE RECONSTRUCTION |
CN112085829A (zh) * | 2019-05-27 | 2020-12-15 | 清华大学 | 基于神经网络的螺旋ct图像重建方法和设备及存储介质 |
CN111476860B (zh) * | 2020-04-22 | 2023-10-24 | 上海联影医疗科技股份有限公司 | 图像重建方法、装置、计算机设备和存储介质 |
KR20220026235A (ko) | 2020-08-25 | 2022-03-04 | 모니터코퍼레이션 주식회사 | 의료영상 간 분석데이터 매칭시스템 및 방법, 화면을 갖는 전자 장치의 그래픽 사용자 인터페이스 |
CN113298903A (zh) * | 2021-05-25 | 2021-08-24 | 苗闯 | 一种大螺距螺旋ct的重建方法、装置、设备和介质 |
CN113397521B (zh) * | 2021-06-15 | 2022-05-27 | 中国科学院自动化研究所 | 螺旋扫描磁粒子投影断层成像方法、系统、设备 |
CN115097535A (zh) * | 2021-07-07 | 2022-09-23 | 清华大学 | 检查系统和方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1172999A (zh) * | 1996-04-10 | 1998-02-11 | 西门子公司 | 计算机层析x射线摄影机的图像再现方法 |
US6252926B1 (en) * | 1998-07-17 | 2001-06-26 | Siemens Aktiengesellschaft | Method for the reconstruction of images from measured values acquired with a CT apparatus by spiral scan and CT apparatus for the implementation of the method |
US7856078B2 (en) * | 2007-05-25 | 2010-12-21 | Siemens Aktiengesellschaft | Method and X-ray CT system for generating computed tomography displays |
CN102096905A (zh) * | 2009-12-10 | 2011-06-15 | 西门子公司 | 双源计算机断层成像拍摄中的降噪 |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5233518A (en) * | 1989-11-13 | 1993-08-03 | General Electric Company | Extrapolative reconstruction method for helical scanning |
WO1998030980A1 (en) * | 1997-01-14 | 1998-07-16 | Edholm, Paul | Technique and arrangement for tomographic imaging |
US5841141A (en) * | 1997-06-03 | 1998-11-24 | The University Of Utah | Image reconstruction from V-projections acquired by Compton camera |
US6201849B1 (en) * | 1999-08-16 | 2001-03-13 | Analogic Corporation | Apparatus and method for reconstruction of volumetric images in a helical scanning cone-beam computed tomography system |
US6301325B1 (en) * | 1999-12-22 | 2001-10-09 | Ge Medical Systems Global Technology Company Llc | Half-scan algorithm for use with a high speed multi-row fan beam helical detector |
US6775346B2 (en) * | 2002-10-21 | 2004-08-10 | Koninklijke Philips Electronics N.V. | Conebeam computed tomography imaging |
DE10252662A1 (de) * | 2002-11-11 | 2004-05-27 | Philips Intellectual Property & Standards Gmbh | Computertomographie-Verfahren mit kohärenten Streustrahlen und Computertomograph |
ATE465471T1 (de) | 2003-02-14 | 2010-05-15 | Koninkl Philips Electronics Nv | System und verfahren für spiralförmige 3d- computertomographie mit exakter rekonstruktion |
ATE341046T1 (de) * | 2003-02-19 | 2006-10-15 | Koninkl Philips Electronics Nv | Verfahren und einrichtung zur volumetrischen bildrekonstruktion |
EP1663004A2 (en) * | 2003-09-05 | 2006-06-07 | Philips Intellectual Property & Standards GmbH | Computer tomography method using a cone-shaped bundle of rays |
EP1719082A1 (en) * | 2004-02-13 | 2006-11-08 | Philips Intellectual Property & Standards GmbH | Motion artifact compensation |
WO2005121836A1 (en) * | 2004-06-09 | 2005-12-22 | Philips Intellectual Property & Standards Gmbh | Computerized tomography method with helical relative movement and conical beam |
US7359478B2 (en) * | 2004-11-18 | 2008-04-15 | Toshiba Medical Systems Corporation | Method for restoring truncated helical cone-beam computed tomography data |
US7778387B2 (en) * | 2005-05-13 | 2010-08-17 | Koninklijke Philips Electronics N. V. | Reconstruction method for helical cone-beam CT |
CN101071110B (zh) * | 2006-05-08 | 2011-05-11 | 清华大学 | 一种基于螺旋扫描立体成像的货物安全检查方法 |
CN101071111B (zh) * | 2006-05-08 | 2011-05-11 | 清华大学 | 一种多视角航空集装箱安全检查系统及方法 |
WO2008041945A1 (en) * | 2006-10-06 | 2008-04-10 | Sidec Technologies Ab | Parallel beam local tomography reconstruction method |
US8184891B2 (en) * | 2006-10-18 | 2012-05-22 | Koninklijke Philips Electronics, N.V. | Imaging system, imaging method and computer program for imaging a region of interest |
EP2078216B1 (en) * | 2006-10-27 | 2015-07-08 | Koninklijke Philips N.V. | Imaging system for imaging an object |
DE102007021023A1 (de) * | 2007-05-04 | 2008-11-13 | Siemens Ag | Verfahren zur Bilderstellung für die Spiral-CT mit veränderlichem Pitch und CT-Gerät zur Durchführung des Verfahrens |
US7724866B2 (en) * | 2007-06-27 | 2010-05-25 | Analogic Corporation | Method of and system for variable pitch computed tomography scanning for baggage screening |
US8050480B2 (en) * | 2008-06-20 | 2011-11-01 | Siemens Aktiengesellschaft | Method and image reconstruction device for generating computed tomography image data |
US8306304B2 (en) * | 2009-08-31 | 2012-11-06 | Siemens Aktiengesellschaft | Precise image reconstruction of spiral CT images at optional pitch values |
US8824760B2 (en) * | 2009-10-20 | 2014-09-02 | University Of Utah Research Foundation | Modification and elimination of back projection weight during the CT image reconstruction |
-
2014
- 2014-05-14 CN CN201410202941.3A patent/CN105093342B/zh active Active
-
2015
- 2015-05-11 US US14/708,866 patent/US9466137B2/en active Active
- 2015-05-12 RU RU2015117682/28A patent/RU2595312C1/ru active
- 2015-05-12 ES ES15167283T patent/ES2741577T3/es active Active
- 2015-05-12 BR BR102015010870-2A patent/BR102015010870B1/pt active IP Right Grant
- 2015-05-12 AU AU2015202554A patent/AU2015202554B2/en active Active
- 2015-05-12 EP EP15167283.9A patent/EP2945125B1/en active Active
- 2015-05-13 KR KR1020150066885A patent/KR101704390B1/ko active IP Right Grant
- 2015-05-13 JP JP2015098129A patent/JP6177831B2/ja active Active
- 2015-05-14 MY MYPI2016703561A patent/MY177410A/en unknown
- 2015-05-14 MX MX2016012501A patent/MX358143B/es active IP Right Grant
- 2015-05-14 SG SG11201608163TA patent/SG11201608163TA/en unknown
- 2015-05-14 WO PCT/CN2015/078956 patent/WO2015172727A1/zh active Application Filing
-
2016
- 2016-05-23 HK HK16105821.6A patent/HK1217769A1/zh unknown
- 2016-09-19 IL IL247916A patent/IL247916B/en active IP Right Grant
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1172999A (zh) * | 1996-04-10 | 1998-02-11 | 西门子公司 | 计算机层析x射线摄影机的图像再现方法 |
US6252926B1 (en) * | 1998-07-17 | 2001-06-26 | Siemens Aktiengesellschaft | Method for the reconstruction of images from measured values acquired with a CT apparatus by spiral scan and CT apparatus for the implementation of the method |
US7856078B2 (en) * | 2007-05-25 | 2010-12-21 | Siemens Aktiengesellschaft | Method and X-ray CT system for generating computed tomography displays |
CN102096905A (zh) * | 2009-12-10 | 2011-06-15 | 西门子公司 | 双源计算机断层成像拍摄中的降噪 |
Also Published As
Publication number | Publication date |
---|---|
RU2595312C1 (ru) | 2016-08-27 |
JP6177831B2 (ja) | 2017-08-09 |
EP2945125A2 (en) | 2015-11-18 |
EP2945125A3 (en) | 2015-11-25 |
SG11201608163TA (en) | 2016-11-29 |
KR20150130940A (ko) | 2015-11-24 |
IL247916A0 (en) | 2016-11-30 |
AU2015202554A1 (en) | 2015-12-03 |
IL247916B (en) | 2021-02-28 |
JP2015219239A (ja) | 2015-12-07 |
US9466137B2 (en) | 2016-10-11 |
US20150332486A1 (en) | 2015-11-19 |
BR102015010870A2 (pt) | 2015-12-08 |
EP2945125B1 (en) | 2019-07-10 |
AU2015202554B2 (en) | 2016-09-29 |
MX358143B (es) | 2018-08-06 |
HK1217769A1 (zh) | 2017-01-20 |
WO2015172727A1 (zh) | 2015-11-19 |
ES2741577T3 (es) | 2020-02-11 |
KR101704390B1 (ko) | 2017-02-08 |
MX2016012501A (es) | 2017-05-09 |
MY177410A (en) | 2020-09-14 |
BR102015010870B1 (pt) | 2022-06-21 |
CN105093342A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093342B (zh) | 螺旋ct系统及重建方法 | |
JP5588849B2 (ja) | X線画像の時間分解能を高めるシステム | |
US20070253523A1 (en) | Method, apparatus, and computer program product for sinogram completion | |
JPH03188832A (ja) | 断層撮影像作成方法および装置 | |
CN100512758C (zh) | X射线计算机断层摄影装置、图象处理装置 | |
US20060109952A1 (en) | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data | |
US9665953B2 (en) | Efficient quasi-exact 3D image reconstruction algorithm for CT scanners | |
EP1489559B1 (en) | Cone-beam reconstruction apparatus and computed tomography apparatus | |
CN1198537C (zh) | 对对象进行计算机断层成像的方法 | |
US7529335B2 (en) | Voxel-driven spiral reconstruction for cone-beam computer tomography | |
CN102456227A (zh) | Ct图像重建方法及装置 | |
CN104198506B (zh) | 小角度自摆式大型多层螺旋ct设备和检查方法 | |
JP6176828B2 (ja) | 画像再構成装置、画像再構成方法およびx線コンピュータ断層撮影装置 | |
US9704223B2 (en) | Method and system for substantially reducing cone beam artifacts based upon adaptive scaling factor in circular computer tomography (CT) | |
US20050175144A1 (en) | Methods and apparatus for artifact reduction in cone beam CT image reconstruction | |
US6937689B2 (en) | Methods and apparatus for image reconstruction in distributed x-ray source CT systems | |
US20080240342A1 (en) | Advanced Csct Detector Shapes | |
US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object | |
Yin et al. | 3D analytic cone-beam reconstruction for multiaxial CT acquisitions | |
US20160048983A1 (en) | Method and system for substantially reducing cone beam artifacts based upon image domain differentiation in circular computer tomography (ct) | |
CN102258383A (zh) | X射线计算机断层摄影装置以及图像处理方法 | |
Kothari et al. | 3D Image Reconstruction Using X-Rays for CT Scan | |
Li et al. | A new supershort-scan algorithm for fan beam and cone beam reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1217769 Country of ref document: HK |
|
GR01 | Patent grant | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1217769 Country of ref document: HK |