CN107427272A - 锥形射束计算机断层摄影投影值提供系统和方法 - Google Patents
锥形射束计算机断层摄影投影值提供系统和方法 Download PDFInfo
- Publication number
- CN107427272A CN107427272A CN201780001069.8A CN201780001069A CN107427272A CN 107427272 A CN107427272 A CN 107427272A CN 201780001069 A CN201780001069 A CN 201780001069A CN 107427272 A CN107427272 A CN 107427272A
- Authority
- CN
- China
- Prior art keywords
- ray
- value
- detecting element
- projection value
- probe
- 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
- 238000002591 computed tomography Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims description 19
- 239000000523 sample Substances 0.000 claims abstract description 95
- 230000005855 radiation Effects 0.000 claims description 38
- 238000001228 spectrum Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 239000004744 fabric Substances 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 description 6
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 4
- 229910052710 silicon Inorganic materials 0.000 description 4
- 239000010703 silicon Substances 0.000 description 4
- 239000002872 contrast media Substances 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种投影值确定设备(14),所述投影值确定设备适于基于由已经被锥形射束(4)的X射线穿过的探测元件的三维布置的探测元件生成的探测值而不是基于由未被所述X射线穿过的探测元件生成的探测值来确定针对所述X射线的投影值。因此,所述投影值确定设备并不基于由未被各自的X射线穿过的探测元件生成的探测值来确定针对各自的X射线的投影值。具体地,同样地,在锥形方向上,仅考虑由已经被所述各自的X射线实际穿过的探测元件已经生成的探测值,以用于生成投影值。这能够引起减少的串扰和具有改善的图像质量的计算机断层摄影图像。
Description
技术领域
本发明涉及锥形射束计算机断层摄影投影值提供系统和方法。本发明还涉及用于提供锥形射束计算机断层摄影投影值的计算机程序。
背景技术
在谱锥形射束计算机断层摄影系统中,用于提供由X射线形成的锥形射束的辐射设备和包括用于生成探测值的探测元件的三维布置的探测器绕与对象对齐的旋转轴旋转,其中,基于所生成的探测值来重建对象的计算机断层摄影图像。
探测器包括探测元件的若干二维正交布置,其中,探测元件的每个二维正交布置形成探测元件的各自的平面,并且其中,探测元件的平面平行于旋转轴并且在X射线的扇形方向上与X射线对齐。
已经在探测器内的不同深度处生成的探测值对应于不同的能量,使得探测值是谱探测值。由于锥形方向上的串扰会降低谱探测值的质量,这可能导致计算机断层摄影图像中的伪影,因此会导致图像质量降低。
发明内容
本发明的目的是要提供一种允许提供锥形射束计算机断层摄影投影值的锥形射束计算机断层摄影投影值提供系统,所述锥形射束计算机断层摄影投影值能够用于重建具有改善的图像质量的计算机断层摄影图像。本发明的另外的目的是要提供对应的锥形射束计算机断层摄影投影值提供系统和用于提供锥形射束计算机断层摄影投影值的计算机程序。
在本发明的第一方面中,提出了一种用于提供对象的锥形射束计算机断层摄影投影值的锥形射束计算机断层摄影投影值提供系统,其中,所述系统包括:
-用于提供探测值的探测值提供设备,所述探测值在旋转设备使得至少辐射设备和所述对象绕旋转轴相对于彼此旋转时通过使用所述辐射设备和探测器已经生成,其中,所述辐射设备适于提供由X射线形成的锥形射束,并且所述探测器包括探测元件的三维布置,其中,所述探测元件被布置用于被所述X射线穿过,并且适于生成取决于已经穿过所述探测元件的所述X射线的探测值,
-用于确定投影值的投影值确定设备,其中,所述投影值确定设备适于基于由被X射线穿过的探测元件生成的探测值而不是基于由未被所述X射线穿过的探测元件生成的探测值来确定针对所述X射线的投影值。
在已知的锥形射束计算机断层摄影系统中,当处理探测值时,假设在探测器的深度方向上被在彼此之上的探测元件生成的探测值属于同一X射线,尽管该假设在锥形方向上并不正确。这可能导致计算机断层摄影图像中的串扰和对应的伪影。因为根据本发明,所述投影值确定设备适于基于由被X射线穿过的探测元件生成的探测值而不是基于由未被X射线穿过的探测元件生成的探测值来确定针对X射线的投影值,为了确定针对X射线的投影值,如果这些探测元件已经被X射线实际穿过,则仅考虑由在深度方向上被布置在彼此之上的探测元件生成的探测值。具体地,同样地,在锥形方向上,仅考虑由已经被X射线实际穿过的探测元件已经生成的探测值,以用于生成针对各自的X射线的投影值。这引起减少的串扰和具有改善的图像质量的计算机断层摄影图像。
所述探测值提供设备能够是存储单元,在所述探测值提供设备中已经存储了探测值,并且能够从中检索探测值以提供所述探测值。所述探测值提供设备也能够是用于从例如探测值采集单元接收探测值并用于提供接收到的探测值的接收单元。所述探测值提供设备也能够是探测值采集单元本身。
所述探测器优选是直接转换半导体探测器。具体地,所述探测器是如X.Liu等人的文章“A Silicon-Strip Detector for Photon-Counting Spectral CT:EnergyResolution From 40keV to 120keV”(IEEE Transactions on Nuclear Science,第61卷,第3期,第1099至1105页,2014年)中所公开的边缘硅条探测器一样的边缘硅条探测器,通过引用将其并入本文。
所述投影值确定设备优选适于仅基于由被X射线穿过的一个或若干探测元件生成的一个或若干探测值来确定针对各自的X射线的投影值,并且不基于由未被X射线穿过的探测元件生成的一个或若干探测值来确定任何投影值。因此,优选地,所有投影值都是基于由被X射线穿过的一个或若干探测元件生成的一个或若干探测值而不是基于由未被X射线穿过的探测元件生成的一个或若干探测值来确定的。
优选地,所述系统还包括重建设备,所述重建设备用于基于所确定的投影值来重建计算机断层摄影图像。
锥形射束具有扇形方向和锥形方向,其中,对于锥形射束的每条X射线,能够分配扇形角和锥形角。各自的X射线的锥形角能够被定义为在锥形方向上的各自的X射线与锥形射束的中心X射线之间的角度。中心X射线优选与旋转轴正交。锥形方向优选平行于旋转轴。各自的X射线的扇形角可以被定义为在扇形方向上的各自的X射线与锥形射束的中心X射线之间的角度,其中,扇形方向优选被定义为与锥形方向正交的方向。探测器内的探测元件的深度优选被定义为探测器的入射表面与各自的探测元件之间的距离。探测器能够包括探测元件的三维布置,其中,探测元件的三维位置可以由锥形方向上的位置、扇形方向上的另外的位置和深度来定义。
优选地,所述探测元件在扇形方向上与所述X射线对齐,并且在锥形方向上与所述X射线不对齐。具体地,所述探测器包括探测元件的若干二维布置,其中,探测元件的每个二维布置形成探测元件的各自的平面,其中,探测元件的所述平面平行于所述旋转轴,并且在扇形方向上与所述X射线对齐,并且在锥形方向上与所述X射线不对齐。由于探测元件在扇形方向上与X射线对齐,因此能够通过探测器的几何结构来减小扇形方向上的串扰。在另外的实施例中,探测元件也能够在锥形方向上与X射线对齐。
在实施例中,探测元件的所述平面在扇形方向上被连续地布置并被交替地移位,使得所述平面被交替地布置在第一位置处和第二位置处,其中,所述第一位置比所述第二位置更靠近所述辐射设备。这允许平面在扇形方向上彼此非常紧密地布置,并且同时提供用于根据期望电连接探测元件的大量空间。所述平面在扇形方向上的这种窄布置允许在生成探测值时进行准确的空间采样,从而进一步提高投影值的质量,并因此进一步提高最终重建的计算机断层摄影图像的质量。
优选地,所述辐射设备适于提供所述锥形射束,使得所述锥形射束从不同的发射位置交替地发射,所述不同的发射位置被布置为沿着与平行于所述旋转轴的方向相对于彼此具有一偏移。这能够增加计算机断层摄影成像过程的视场,其中,投影值的质量以及因此经重建的计算机断层摄影图像的质量仍然相对较高,这是因为同样地在这种情况下为了确定针对X射线的投影值,仅被该X射线实际穿过的探测值被考虑。
优选地,所述探测器和所述投影值确定设备适于使得各自的投影值指示所述X射线的光子数量。具体地,优选地,所述探测器和所述投影值确定设备适于使得所述各自的投影值指示处于各自的能量范围内的各自的X射线的光子数量,以便提供谱投影值。在优选实施例中,所述各自的能量范围取决于已经生成了这样的各自的探测值的各自的探测元件的深度:对所述各自的投影值的所述确定基于所述各自的探测值。为了确定针对X射线且针对某一能量范围的投影值,能够使用由被布置在对应于所述能量范围的一深度处且已经被所述X射线穿过的探测元件已经生成的一个或若干探测值。如果使用若干探测值来确定投影值,则这些探测值可以通过例如内插来组合。所述内插能够包括对探测值的加权,其中,各自的权重能够取决于各自的探测元件的中心到各自的X射线的中心的距离,或者能够取决于另一量度,如各自的X射线和各自的探测元件的重叠程度。在实施例中,使用最近邻内插。
在实施例中,所述各自的能量范围取决于所述各自的X射线从所述探测器的入射表面到已经生成了这样的各自的探测值的所述各自的探测元件的长度:对所述各自的投影值的所述确定基于所述各自的探测值。因此,不仅可以考虑各自的探测元件的深度,而且当确定了针对各自的投影值的能量范围时,可以考虑光子在被探测元件探测到之前的实际透射长度。这能够引起针对各自的光子的能量范围的改进的确定,并因此引起针对各自的投影值的改进的确定。所述各自的能量范围可以取决于所述各自的X射线的锥形角。由于通过探测器的各自的光子在被探测元件探测到之前的透射长度取决于各自的X射线的锥形角,因此能够通过使用锥形角来相对简单地考虑透射长度。
在本发明的另外的方面中,提出了一种用于提供对象的锥形射束计算机断层摄影投影值的锥形射束计算机断层摄影投影值提供方法,其中,所述方法包括:
-由探测值提供设备提供探测值,其中,所提供的探测值在旋转设备使得至少辐射设备和所述对象绕旋转轴相对于彼此旋转时通过使用所述辐射设备和探测器已经生成,其中,所述辐射设备适于提供由X射线形成的锥形射束,并且所述探测器包括探测元件的三维布置,其中,所述探测元件被布置用于被所述X射线穿过,并且适于生成取决于已经穿过探测元件的X射线的探测值,
-由投影值确定设备确定投影值,其中,所述投影值确定设备基于由被X射线穿过的探测元件生成的探测值而不是基于由未被所述X射线穿过的探测元件生成的探测值来确定针对所述X射线的投影值。
在本发明的另一方面中,提出了一种用于提供对象的锥形射束计算机断层摄影投影值的计算机程序,其中,所述计算机程序包括程序代码模块,所述程序代码模块用于当所述计算机程序在根据权利要求1所述的锥形射束计算机断层摄影投影值提供系统上运行时令所述锥形射束计算机断层摄影投影值提供系统执行根据权利要求13所述的锥形射束计算机断层摄影投影值提供方法。
应当理解,权利要求1所述的锥形射束计算机断层摄影投影值提供系统,权利要求13所述的锥形射束计算机断层摄影投影值提供方法和权利要求14所述的计算机程序具有相似和/或相同的优选实施例,具体地,如在从属权利要求中定义的那样。
应当理解,本发明的优选实施例也能够是从属权利要求或上述实施例与各自的独立权利要求的任何组合。
参考下文描述的实施例,本发明的这些方面和其它方面将会变得明显且得以阐明。
附图说明
在下列附图中:
图1示意性且示范性地示出了锥形射束计算机断层摄影投影值提供系统的实施例,
图2示意性且示范性地示出了探测器的探测元件的二维布置和锥形射束的X射线,
图3示意性且示范性地图示了沿着旋转轴看到的扇形方向上的探测器几何形状,
图4示意性且示范性地图示了从辐射设备看到的相同的探测器几何形状,
图5示出了示范性地图示锥形射束计算机断层摄影投影值提供方法的实施例的流程图,
图6示意性且示范性地示出了探测器的探测元件的二维布置和从第一发射位置发射的X射线,并且
图7示意性且示范性地示出了探测器的探测元件的二维布置和从第二发射位置发射的X射线。
具体实施方式
图1示意性且示范性地示出了锥形射束计算机断层摄影投影值提供系统的实施例,在该实施例中,所述锥形射束计算机断层摄影投影值提供系统是谱锥形射束计算机断层摄影成像系统。系统10包括能够围绕平行于z方向延伸的旋转轴R旋转的机架1。包括辐射源2的辐射设备40被安装在机架1上,在该实施例中辐射源2是X射线管和准直器3。准直器3根据由X射线管2生成的辐射形成锥形射束4。锥形射束4穿过检查区5内的对象,在该实施例中检查区5为圆柱形的。在穿过检查区5并因此穿过对象之后,锥形射束4入射到探测器6上,所述探测器6包括二维入射表面,其中,探测器6被安装在机架1上。探测器6优选是直接转换半导体探测器,特别是边缘硅条探测器。
系统10包括两个电动机7、8。机架1由电动机7以优选恒定但可调节的角速度驱动。电动机8被提供用于使被布置在平行于旋转轴R或z轴的方向的检查区5中的如患者台的支撑件上的对象移位。例如,这些电动机7、8由控制单元9控制,使得辐射设备40和检查区5内的对象沿着螺旋轨迹相对于彼此移动。然而,辐射设备40和对象也能够沿着如圆形轨迹的另一轨迹相对于彼此移动。电动机7和机架1用于使辐射设备40和探测器6绕对象旋转,因此可以被认为是旋转设备的部件。
探测器6包括探测元件的三维布置,其中,探测元件被布置用于被锥形射束4的X射线穿过,并且适于在辐射设备40和对象的相对运动期间根据已经穿过探测元件的X射线来生成探测值。
探测器6包括探测元件20的若干二维布置,如图3和图4中示意性且示范性地图示的那样,其中,所述若干二维布置形成三维布置,并且在图2中示出了单个二维布置。探测元件20的每个二维布置23形成探测元件20的各自的平面23,所述探测元件20的所述各自的平面23也能够认为是探测元件20的各自的厚片23,其中,探测元件20的平面23平行于旋转轴R,并且在扇形方向26上与X射线21对齐,并且在锥形方向25上与X射线21不对齐。探测元件20的平面23在扇形方向26上被连续地布置并被交替地移位,使得平面23被交替地布置在第一位置处和第二位置处,其中,所述第一位置比所述第二位置更靠近辐射设备40,如图3中示意性且示范性地图示的那样。图3图示了在平面23上在旋转轴R或z轴的方向上的视图,图4图示了从辐射设备40看到的平面23上的视图。如图2所图示的,在各自的平面23内,探测元件20形成矩形网格,其中,探测元件被布置在锥形方向25上且被布置在不同的深度处。深度增加的方向由箭头27指示。
由于辐射设备40和探测器6在它们绕对象旋转时生成探测值,因此辐射设备40和探测器6能够被认为是探测值提供设备的组件,在该实施例中,探测值提供设备是探测值采集单元12。
探测值被提供给用于确定投影值的投影值确定设备13。投影值确定设备13适于基于由被X射线22穿过的探测元件20生成的探测值而不是基于由未被X射线22穿过的探测元件20生成的探测值来确定针对各自的X射线22的投影值。优选地,投影值确定设备13适于基于由被各自的X射线22穿过的探测元件20生成的探测值来针对每条各自的X射线22确定与不同的能量范围相对应的若干投影值。针对各自的X射线22和各自的能量范围确定的投影值优选指示在各自的能量范围内的各自的X射线22的光子数量。探测元件20被布置在不同的深度处,即,例如在图2中示意性且示范性地图示的五个不同深度处。能够向每个深度分配能量范围,使得针对每条X射线22,能够确定若干投影值,所述若干投影值对应于不同的深度并且因此对应于不同的能量范围。如果X射线22在相同的深度处穿过若干探测元件,则能够组合对应的探测值,以用于通过例如内插来生成针对该深度的投影值。
形成锥形射束4的X射线能够通过在锥形方向25上定义多条X射线并通过相应地划分锥形射束4使得它们具有相同的发散宽度来定义。在图2中示意性地示出的范例中,锥形射束4被划分成具有相同发散宽度的九条X射线,其中,在图2中,箭头21、22指代具有发散宽度的各自的X射线的各自的中心线,并且被X射线22穿过的探测元件20被突出显示。X射线22仅在三个最小深度位置以及最大深度位置处穿过单个探测元件20,但是由于也考虑了X射线22的宽度,因此在第二最大深度位置处,两个探测元件20被X射线22穿过。因此,在该范例中,由第二最大深度位置中的这些被穿过的探测元件生成的探测值被组合以用于确定针对对应的能量范围和X射线22的投影值。
在实施例中,投影值确定设备13适于根据各自的X射线22从探测器6的入射表面到探测元件20的长度来将能量范围分配给探测元件20。例如,在图2中,为了确定针对在探测元件的二维布置内的左下角中的探测元件且针对X射线22的能量范围,X射线22在探测器6内从入射表面22到探测元件20的长度被确定,并且基于该长度来确定能量范围,其中,随着长度的增加,能量范围覆盖较大的能量。在实施例中,通过基于各自的射线22的锥形角定义各自的能量范围来考虑各自的长度,其中,假设采集几何形状,特别是辐射设备40相对于探测元件20的位置是已知的。
系统10还包括重建设备14,所述重建设备14用于基于所确定的谱投影值来重建对象的计算机断层摄影图像,其中,重建设备14以及投影值确定设备13可以由控制单元9或另一控制单元来控制。在该实施例中,重建设备14适于使用滤波反投影算法来重建计算机断层摄影图像。重建设备14能够适于使用材料分解技术,以便基于谱投影值来重建与人体内的不同物质相对应的不同的计算机断层摄影图像。例如,如果将造影剂注射到对象中,则能够重建第一图像,所述第一图像仅示出对象内的造影剂,并且能够在不使用造影剂的情况下重建示出对象的第二图像。重建设备14也能够适于基于谱投影值来重建与如康普顿效应和光电效应的不同物理效应相对应的不同的计算机断层摄影图像。当然,重建设备14也能够适用于使用其他谱重建技术。可以由重建设备14使用的已知的重建技术例如在E.Roessl和R.Proksa的文章“K-edge imaging in X-ray computed tomography using multi-binphoton counting detector”(Physics in Medicine and Biology,第52卷,第4679至4696页,2007年)中被公开。重建图像能够被示出在显示器11上。
在下文中,将参考图5所示的流程图来示例性地描述锥形射束计算机断层摄影投影值提供方法的实施例。
在该实施例中,锥形射束计算机断层摄影投影值提供方法是用于生成计算机断层摄影图像的谱锥形射束计算机断层摄影成像方法。在步骤101中,由探测值采集单元12提供探测值。具体地,辐射设备40提供由X射线21形成的锥形射束4,并且探测器6的探测元件20生成取决于已经穿过探测元件20的X射线21的探测值。
在步骤102中,投影值确定设备13基于所提供的探测值来确定投影值,其中,针对各自的X射线22的投影值是基于由被X射线22穿过的探测元件20生成的探测值而不是基于由未被X射线22穿过的探测元件生成的探测值来确定的。在该实施例中,投影值确定设备13针对每条X射线并针对不同的能量范围确定在各自的能量范围内的表示各自的X射线的光子数量的投影值。在步骤103中,重建设备14基于在步骤102中确定的谱投影值来重建一幅或多幅计算机断层摄影图像。在步骤104中,在显示器11上示出一幅或多幅重建的计算机断层摄影图像。
在实施例中,辐射设备40适于提供锥形射束4,使得所述锥形射束4从不同的发射位置交替地发射,所述不同的发射位置被布置为沿着平行于旋转轴R的方向相对于彼此具有一偏移。这在图6和图7中示意性且示范性地进行图示,其中,在图6中,包括X射线34的锥形射束4从第一发射位置30发射,并且其中,在图7中,具有X射线34的锥形射束4从第二位置31发射,其中,所述第一位置和所述第二位置在与锥形方向25相对应的z方向上相对于彼此具有一偏移。在图6和图7中,突出显示了用于确定针对X射线32的投影值的、被X射线32穿过并因此生成探测值的探测元件20。从图6中能够看出,由于X射线32的锥形角相对较小,因此仅考虑彼此叠置的探测元件20。在图7中,X射线33具有相对较大的锥形角,使得考虑了由不同列的探测元件20生成的探测值来确定针对X射线33的投影值。
在已知的谱锥形射束计算机断层摄影系统中,如果锥形角变大和/或如果使用立体X射线管(即,辐射源从不同的发射位置交替地发射锥形射束,所述不同的发射位置被布置为沿着平行于旋转轴的方向相对于彼此具有一偏移),则上述串扰效应变得更加显著。与用于读出具有探测元件的三维布置的探测器的已知读出样式相反,其中假设在探测器的深度方向上被布置在彼此之上的探测元件已经被相同的X射线穿过,投影值确定设备使用适于实际X射线几何形状的修改的读出样式。
如果X射线源不改变相对于辐射源的发射锥形射束的发射位置,则相对于探测器的X射线几何形状在不同的视图中保持恒定,即,当绕对象旋转辐射设备时辐射源位置保持恒定。因此,投影值确定设备用于基于探测值确定投影值的读出样式也能够是恒定的。然而,如果辐射源是立体X射线管,则当发射位置相对于辐射源改变时,X射线几何形状相对于探测器改变,如图6和图7所图示的。因此,用于确定投影值的投影值确定设备使用的读出样式也相应地改变,以便根据当前的X射线几何形状来考虑用于确定针对X射线的投影值的探测值。这确保针对X射线的投影值是基于由被X射线实际穿过的一个或若干探测元件已经生成的一个或若干探测值而不是基于由未被X射线穿过的探测元件生成的探测值来确定的。因此,优选根据由立体管提供的锥形射束的发射位置的切换来切换读出样式。读出几何形状针对X射线几何形状的这种调整优选仅在锥形角方向上(即,在锥形方向上)执行,这是因为优选地探测元件已经在扇形角方向(即,在扇形方向上)与X射线几何形状物理地对齐。
图3示意性且示范性地图示了沿着旋转轴或z轴看到的扇形方向26上的探测器几何形状。在图4中图示了与从辐射设备看到的相同的探测器几何形状。图2图示了包括探测元件的二维布置的单个平面,其中,该平面可以被认为是单边缘硅条探测器阵列。纯笛卡尔读出几何形状与X射线的几何形状并不匹配。为了补偿当使用纯笛卡尔读出结果时会发生的串扰,投影值确定设备优选适于使用与X射线几何形状一致的读出几何形状。
在上述参照图2所述的实施例中,存在探测元件的五个深度位置,其中,每个深度位置能够对应于各自的能量范围。在另一实施例中,被布置在不同深度位置处的探测元件能够被一起读出,使得并不区分这些不同的深度位置,即,能够组合邻近的深度位置,并且能够一起读出被布置在这些深度位置处的探测元件,例如,能够由相同的专用集成电路(ASIC)读出,以便降低用于以减少的深度以及因此的降低的能量分辨率为代价来读出探测元件的技术努力。对于不同的探测区,在包括被一起读出的探测元件的各自的探测区的深度方向上的物理宽度能够是不同的。具体地,该宽度能够随着深度的增加而增加。因此,积分区间(即,探测值被积分的深度区间)优选随探测器中的深度增加而增加,以便考虑吸收深度的能量依赖性。
本领域技术人员通过研究附图、公开内容以及权利要求,在实践请求保护的发明时能够理解并实现对所公开的实施例的其他变型。
在权利要求中,“包括”一词不排除其他元件或步骤,并且词语“一”或“一个”不排除多个。
单个单元或设备可以实现在权利要求中记载的若干项的功能。尽管某些措施被记载在互不相同的从属权利要求中,但是这并不指示不能有利地使用这些措施的组合。
能够通过任何其他数量的单元或设备来执行由一个或多个单元或设备执行的诸如对投影值的确定或对计算机断层摄影图像的重建的操作。锥形射束计算机断层摄影投影值提供系统根据锥形射束计算机断层摄影投影值提供方法的这些操作和/或控制能够被实施为计算机程序的程序代码模块和/或专用硬件。
计算机程序可以被存储/分布在合适的介质上,例如与其他硬件一起或作为其他硬件的部分供应的光学存储介质或固态介质,但是也可以被以其他形式分布,例如经由互联网或其他有线或无线的电信系统。
权利要求中的任何附图标记都不应被解释为对范围的限制。
本发明涉及一种投影值确定设备,所述投影值确定设备适于基于由已经被锥形射束的X射线穿过的探测元件的三维布置的探测元件生成的探测值而不是基于由未被所述X射线穿过的探测元件生成的探测值来确定针对所述X射线的投影值。因此,所述投影值确定设备并不基于由未被各自的X射线穿过的探测元件生成的探测值来确定针对各自的X射线的投影值。具体地,同样地,在锥形方向上,仅考虑由已经被所述各自的X射线实际穿过的探测元件已经生成的探测值,以用于生成投影值。这能够引起减少的串扰和具有改善的图像质量的计算机断层摄影图像。
Claims (14)
1.一种用于提供对象的锥形射束计算机断层摄影投影值的锥形射束计算机断层摄影投影值提供系统,所述系统包括:
-用于提供探测值的探测值提供设备(12),所述探测值在旋转设备(1、8)使得至少辐射设备(40)和所述对象绕旋转轴(R)相对于彼此旋转时通过使用所述辐射设备(40)和探测器(6)已经生成,其中,所述辐射设备(40)适于提供由X射线(21;34)形成的锥形射束(4),并且所述探测器(6)包括探测元件(20)的三维布置,其中,所述探测元件(20)被布置用于被所述X射线(21;34)穿过,并且适于生成取决于已经穿过所述探测元件(20)的所述X射线(21;34)的探测值,
-用于确定投影值的投影值确定设备(13),其中,所述投影值确定设备(13)适于基于由被X射线(22;32;33)穿过的探测元件(20)生成的探测值而不是基于由未被所述X射线(22;32;33)穿过的探测元件(20)生成的探测值来确定针对所述X射线(22;32;33)的投影值。
2.根据权利要求1所述的系统,其中,所述探测元件(20)在扇形方向上与所述X射线(21;34)对齐,并且在锥形方向上与所述X射线(21;34)不对齐。
3.根据权利要求2所述的系统,其中,所述探测器(6)包括探测元件(20)的若干二维布置,其中,探测元件(20)的每个二维布置形成探测元件(20)的各自的平面(23、24),其中,探测元件(20)的所述平面(23、24)平行于所述旋转轴(R),并且在扇形方向上与所述X射线(21;34)对齐,并且在锥形方向上与所述X射线(21;34)不对齐。
4.根据权利要求3所述的系统,其中,探测元件(20)的所述平面(23、24)在扇形方向上被连续地布置并被交替地移位,使得所述平面(23、24)被交替地布置在第一位置处和第二位置处,其中,所述第一位置比所述第二位置更靠近所述辐射设备(40)。
5.根据权利要求1所述的系统,其中,所述辐射设备(40)适于提供所述锥形射束(4),使得所述锥形射束从不同的发射位置交替地发射,所述不同的发射位置被布置为沿着平行于所述旋转轴(R)的方向相对于彼此具有一偏移。
6.根据权利要求1所述的系统,其中,所述探测器(6)和所述投影值确定设备(13)适于使得各自的投影值指示所述X射线(22;32;33)的光子数量。
7.根据权利要求6所述的系统,其中,所述探测器(6)和所述投影值确定设备(13)适于使得所述各自的投影值指示处于各自的能量范围内的各自的X射线(22;32;33)的光子数量,以便提供谱投影值。
8.根据权利要求7所述的系统,其中,所述各自的能量范围取决于已经生成了这样的各自的探测值的各自的探测元件(20)的深度:对所述各自的投影值的所述确定基于所述各自的探测值。
9.根据权利要求7所述的系统,其中,所述各自的能量范围取决于所述各自的X射线(22;32;33)从所述探测器(6)的入射表面到已经生成了这样的各自的探测值的所述各自的探测元件(20)的长度:对所述各自的投影值的所述确定基于所述各自的探测值。
10.根据权利要求7所述的系统,其中,所述各自的能量范围取决于所述各自的X射线(22;32;33)的锥形角。
11.根据权利要求1所述的系统,其中,所述探测器(6)是直接转换半导体探测器(6)。
12.根据权利要求1所述的系统,其中,所述系统还包括重建设备(14),所述重建设备用于基于所确定的投影值来重建计算机断层摄影图像。
13.一种用于提供对象的锥形射束计算机断层摄影投影值的锥形射束计算机断层摄影投影值提供方法,所述方法包括:
-由探测值提供设备(12)提供探测值,其中,所提供的探测值在旋转设备(1、8)使得至少辐射设备(40)和所述对象绕旋转轴(R)相对于彼此旋转时通过使用所述辐射设备(40)和探测器(6)已经生成,其中,所述辐射设备(40)适于提供由X射线(21;34)形成的锥形射束(4),并且所述探测器(6)包括探测元件(20)的三维布置,其中,所述探测元件(20)被布置用于被所述X射线(21;34)穿过,并且适于生成取决于已经穿过所述探测元件(20)的所述X射线(21;34)的探测值,
-由投影值确定设备(13)确定投影值,其中,所述投影值确定设备(13)基于由被X射线(22;32;33)穿过的探测元件(20)生成的探测值而不是基于由未被所述X射线(22;32;33)穿过的探测元件生成的探测值来确定针对所述X射线(22;32;33)的投影值。
14.一种用于提供对象的锥形射束计算机断层摄影投影值的计算机程序,所述计算机程序包括程序代码模块,所述程序代码模块用于当所述计算机程序在根据权利要求1所述的锥形射束计算机断层摄影投影值提供系统上运行时令所述锥形射束计算机断层摄影投影值提供系统执行根据权利要求13所述的锥形射束计算机断层摄影投影值提供方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP16153341 | 2016-01-29 | ||
EP16153341.9 | 2016-01-29 | ||
PCT/EP2017/051017 WO2017129459A1 (en) | 2016-01-29 | 2017-01-19 | Cone beam computed tomography projection values providing system and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107427272A true CN107427272A (zh) | 2017-12-01 |
CN107427272B CN107427272B (zh) | 2022-06-14 |
Family
ID=55361328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780001069.8A Active CN107427272B (zh) | 2016-01-29 | 2017-01-19 | 锥形射束计算机断层摄影投影值提供系统和方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10485498B2 (zh) |
EP (1) | EP3267893B1 (zh) |
JP (1) | JP6585851B2 (zh) |
CN (1) | CN107427272B (zh) |
WO (1) | WO2017129459A1 (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070140414A1 (en) * | 2005-12-07 | 2007-06-21 | Walter Garms | Apparatus and method for providing a near-parallel projection from helical scan data |
CN101039623A (zh) * | 2004-10-15 | 2007-09-19 | 皇家飞利浦电子股份有限公司 | 计算机层析成像方法 |
US20090016485A1 (en) * | 2007-07-10 | 2009-01-15 | Kabushiki Kaisha Toshiba | X-ray computed tomography apparatus, reconstruction processing apparatus, and image processing apparatus |
CN103065279A (zh) * | 2011-10-19 | 2013-04-24 | 株式会社东芝 | 圆轨道锥束计算机断层摄影(ct)中用于大幅度地减轻伪影的方法以及系统 |
CN103262124A (zh) * | 2010-12-17 | 2013-08-21 | 皇家飞利浦电子股份有限公司 | 用于对感兴趣区域成像的成像系统 |
CN103308535A (zh) * | 2012-03-09 | 2013-09-18 | 同方威视技术股份有限公司 | 用于射线扫描成像的设备和方法 |
CN103717137A (zh) * | 2011-07-28 | 2014-04-09 | 皇家飞利浦有限公司 | 图像生成装置 |
CN105093342A (zh) * | 2014-05-14 | 2015-11-25 | 同方威视技术股份有限公司 | 螺旋ct系统及重建方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6720812B2 (en) | 1995-06-02 | 2004-04-13 | Nova R&D, Inc. | Multi-channel integrated circuit |
US6018561A (en) | 1998-07-27 | 2000-01-25 | Siemens Corporate Research, Inc. | Mask boundary correction in a cone beam imaging system using simplified filtered backprojection image reconstruction |
DE19956585A1 (de) | 1999-11-25 | 2001-05-31 | Philips Corp Intellectual Pty | Computertomographie-Verfahren |
JP2004271333A (ja) | 2003-03-07 | 2004-09-30 | Hamamatsu Photonics Kk | シンチレータパネル、イメージセンサ及びエネルギー弁別器 |
JP3942178B2 (ja) * | 2003-07-29 | 2007-07-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ctシステム |
US7095028B2 (en) | 2003-10-15 | 2006-08-22 | Varian Medical Systems | Multi-slice flat panel computed tomography |
US7333587B2 (en) | 2004-02-27 | 2008-02-19 | General Electric Company | Method and system for imaging using multiple offset X-ray emission points |
JP2007021021A (ja) | 2005-07-20 | 2007-02-01 | Ge Medical Systems Global Technology Co Llc | 画像処理装置およびx線ct装置 |
JP4884765B2 (ja) | 2005-12-27 | 2012-02-29 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US8290116B2 (en) | 2008-09-16 | 2012-10-16 | Koninklijke Philips Electronics, N.V. | Imaging apparatus including correction unit for scattered radiation |
CA2650066A1 (en) | 2009-01-16 | 2010-07-16 | Karim S. Karim | Photon counting and integrating pixel readout architecture with dynamic switching operation |
US8466421B2 (en) | 2010-07-30 | 2013-06-18 | Varian Medical Systems Inc. | Radiation detector with multiple operating schemes |
US9347893B2 (en) | 2011-07-26 | 2016-05-24 | Robert Sigurd Nelson | Enhanced resolution imaging systems for digital radiography |
WO2014180487A1 (en) | 2013-05-07 | 2014-11-13 | Cern - European Organization For Nuclear Research | A detector configuration with semiconductor photomultiplier strips and differential readout |
US9488739B2 (en) * | 2013-12-18 | 2016-11-08 | The Board Of Trustees Of The Leland Stanford Junior University | Spectral imaging system and method |
JP6448930B2 (ja) * | 2014-06-30 | 2019-01-09 | キヤノンメディカルシステムズ株式会社 | フォトンカウンティング型x線ct装置及びフォトンカウンティング型イメージングプログラム |
-
2017
- 2017-01-19 CN CN201780001069.8A patent/CN107427272B/zh active Active
- 2017-01-19 WO PCT/EP2017/051017 patent/WO2017129459A1/en active Application Filing
- 2017-01-19 US US15/554,439 patent/US10485498B2/en active Active
- 2017-01-19 EP EP17702309.0A patent/EP3267893B1/en active Active
- 2017-01-19 JP JP2018537784A patent/JP6585851B2/ja active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101039623A (zh) * | 2004-10-15 | 2007-09-19 | 皇家飞利浦电子股份有限公司 | 计算机层析成像方法 |
US20070140414A1 (en) * | 2005-12-07 | 2007-06-21 | Walter Garms | Apparatus and method for providing a near-parallel projection from helical scan data |
US20090016485A1 (en) * | 2007-07-10 | 2009-01-15 | Kabushiki Kaisha Toshiba | X-ray computed tomography apparatus, reconstruction processing apparatus, and image processing apparatus |
CN103262124A (zh) * | 2010-12-17 | 2013-08-21 | 皇家飞利浦电子股份有限公司 | 用于对感兴趣区域成像的成像系统 |
CN103717137A (zh) * | 2011-07-28 | 2014-04-09 | 皇家飞利浦有限公司 | 图像生成装置 |
CN103065279A (zh) * | 2011-10-19 | 2013-04-24 | 株式会社东芝 | 圆轨道锥束计算机断层摄影(ct)中用于大幅度地减轻伪影的方法以及系统 |
CN103308535A (zh) * | 2012-03-09 | 2013-09-18 | 同方威视技术股份有限公司 | 用于射线扫描成像的设备和方法 |
CN105093342A (zh) * | 2014-05-14 | 2015-11-25 | 同方威视技术股份有限公司 | 螺旋ct系统及重建方法 |
Non-Patent Citations (2)
Title |
---|
MUDERS等: "Projection-Wise Filter Optimization for Limited-Angle Cone-Beam CT Using the Approximate Inverse", 《IEEE TRANSACTIONS ON NUCLEAR SCIENCE》 * |
R SCOTT BROCK等: "Reconstruction of a cone-beam CT image via forward iterative projection matching", 《MEDICAL PHYSICS》 * |
Also Published As
Publication number | Publication date |
---|---|
US20190105003A1 (en) | 2019-04-11 |
CN107427272B (zh) | 2022-06-14 |
JP2019502484A (ja) | 2019-01-31 |
EP3267893B1 (en) | 2018-08-15 |
WO2017129459A1 (en) | 2017-08-03 |
US10485498B2 (en) | 2019-11-26 |
EP3267893A1 (en) | 2018-01-17 |
JP6585851B2 (ja) | 2019-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103384498B (zh) | 探测装置 | |
CN106488744B (zh) | X射线拍摄装置以及图像重建方法 | |
CN100380405C (zh) | 使用三维背景投射的计算机断层扫描装置和方法 | |
CN106030293A (zh) | X射线检查装置以及x射线检查方法 | |
CN103959097B (zh) | 包括用于探测x射线辐射的两个闪烁体的探测装置 | |
CN1794951A (zh) | 扇形射束相干散射计算机断层摄影 | |
US20110200172A1 (en) | Method and apparatus for computed imaging backscatter radiography | |
CN1946342A (zh) | 锥束相干散射计算机断层摄影装置 | |
JP2008501463A (ja) | 干渉性散乱コンピュータ断層撮影装置及び方法 | |
KR20080084798A (ko) | 이중 판독 스캐너를 이용한 영상화 방법 및 장치 | |
US20100119034A1 (en) | Computed tomography device and method using circular-pixel position-adaptive interpolation | |
CN106901767A (zh) | X射线系统和图像重建方法 | |
CN102715914B (zh) | G型臂x光机的三维图像生成方法及装置与g型臂x光机 | |
JP6730424B2 (ja) | X線検出システム、x線装置、並びに、x線検出データを処理する装置及び方法 | |
CN105828719A (zh) | X射线成像装置 | |
CN1352919A (zh) | 准直器控制方法和设备,以及x线ct设备 | |
CN106896122A (zh) | 液体检测方法和系统 | |
US7515675B2 (en) | Apparatus and method for providing a near-parallel projection from helical scan data | |
CN102175699B (zh) | 检查系统、检查方法、ct装置以及探测装置 | |
US20070297575A1 (en) | Systems and methods for determining object position | |
CN101785029A (zh) | 根据能量相关投影数据对感兴趣区域进行成像的成像系统 | |
CN101453951A (zh) | 锥束ct半周期闭合螺旋轨迹 | |
CN105074501B (zh) | 放射线检测器和具备该放射线检测器的x射线ct装置 | |
CN101291627A (zh) | 先进的csct检测器形状 | |
CN107427272A (zh) | 锥形射束计算机断层摄影投影值提供系统和方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |