CN108135557A - 计算机断层摄影x射线成像 - Google Patents
计算机断层摄影x射线成像 Download PDFInfo
- Publication number
- CN108135557A CN108135557A CN201780003448.0A CN201780003448A CN108135557A CN 108135557 A CN108135557 A CN 108135557A CN 201780003448 A CN201780003448 A CN 201780003448A CN 108135557 A CN108135557 A CN 108135557A
- Authority
- CN
- China
- Prior art keywords
- ray
- detector
- projection
- data
- datas
- 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 107
- 238000003384 imaging method Methods 0.000 title claims abstract description 40
- 241001269238 Data Species 0.000 claims abstract description 71
- 230000005855 radiation Effects 0.000 claims abstract description 63
- 238000010606 normalization Methods 0.000 claims abstract description 31
- 238000012937 correction Methods 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 27
- 238000005259 measurement Methods 0.000 claims description 34
- 238000000034 method Methods 0.000 claims description 29
- 238000004590 computer program Methods 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims 1
- 230000008569 process Effects 0.000 description 5
- 230000002238 attenuated effect Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007430 reference method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 239000012780 transparent material Substances 0.000 description 1
- 238000012800 visualization Methods 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- 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
-
- 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
- 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
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
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)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- General Physics & Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及计算机断层摄影X射线成像。为了提供用于重建的进一步改进的数据,提供了一种用于计算机断层摄影X射线成像的系统(10)。所述系统包括数据接口(12)和处理单元(14)。所述数据接口被配置为提供针对至少第一X射线能量范围和第二X射线能量范围的至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据,所述第一X射线能量范围与所述第二X射线能量范围彼此不同。所述处理单元被配置为确定针对切片归一化的校正,并且对所述第一CT X射线投影数据和所述第二CT X射线投影数据应用均等切片归一化,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据。为了校正,均等切片归一化基于外面的探测器元件的测得的数据。此外,数据接口被配置为提供经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。在范例中,所述系统还包括计算机断层摄影X射线成像采集装置(20),所述计算机断层摄影X射线成像采集装置具有被配置为生成X射线束的X射线源(22)和被配置为能量分辨X射线探测器的X射线探测器(26)以同时提供分别针对至少两个不同的X射线能量范围的X射线辐射投影数据。所述计算机断层摄影X射线成像采集装置被配置为采集针对所述至少第一X射线能量范围和第二X射线能量范围的所述对象的感兴趣区域的至少第一CT X射线投影数据和第二CT X射线投影数据。
Description
技术领域
本发明涉及用于计算机断层摄影X射线成像的系统和用于计算机断层摄影X射线成像的方法,以及用于控制装置的计算机程序单元,以及计算机可读存储介质。
背景技术
在谱计算机断层摄影(CT)中,针对几个能量水平或范围探测X射线辐射,以提供增强的图像投影数据用于重建图像。在谱CT成像中,当组合若干幅投影数据图像时,作为预处理流程,可以提供调整。例如,US 2016/0225169描述了一种局部调整正则化参数的方法,用于完全3D迭代的CT重建中的图像质量优化。这可能基于图像测量并具有补偿差异的目的。但是,已经表明这可能会受到噪声的影响。
发明内容
可能需要针对重建提供进一步改进的数据。
本发明的目的通过独立权利要求的主题来解决;在从属权利要求中并入了另外的实施例。应该注意,本发明的以下描述的方面也适用于:用于计算机断层摄影X射线成像的系统和用于计算机断层摄影X射线成像的方法以及用于控制设备的计算机程序单元和计算机可读介质。
根据本发明,提供了一种用于计算机断层摄影X射线成像的系统。所述系统包括数据接口和处理单元。所述数据接口被配置为提供对象的针对至少第一X射线能量范围和第二X射线能量范围的至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据,所述第一X射线能量范围与所述第二X射线能量范围彼此不同。所述对象的所述至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据包括来自不同角度的X射线测量结果的多个切片,以产生所述对象的截面图像。所述处理单元被配置为确定针对X射线测量结果的多个切片的归一化的校正,以关于所述第一CT X射线投影数据和所述第二CT X射线投影数据的像素强度值来对它们进行改变。所述处理单元还被配置为应用针对所述第一CT X射线投影数据和所述第二CT X射线投影数据的均等切片归一化,其中,所述切片归一化中的X射线测量结果中的所述多个切片的所述像素强度值是能改变的,使得提供对预定值的参考,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据。为了校正,所述均等切片归一化基于来自探测如下的X射线辐射的探测器区域的探测器测量的测得的数据,所述X射线辐射在成像期间未受被布置在X射线源与X射线探测器之间的辐射路径中的所述对象影响。此外,所述数据接口被配置为提供经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。
术语“未受影响”是指有意免于影响,即免于衰减。因此,X射线辐射不被对象(例如患者)衰减。X射线辐射也不被其他对象衰减,例如针对患者或设备的保持和支持装置。
换句话说,为了校正,均等切片归一化基于在成像期间不受布置在辐射路径中的对象影响的探测器测量结果的测得的数据。
这意味着,为了校正,这些探测器测量结果的测得的数据被提供作为均等切片归一化的基础,这些测量结果是在CT X射线投影数据采集期间从具有到X射线源的直接视线的探测器区域得到的,即在X射线图像数据采集期间。
归一化所基于的测得的数据是在实际成像扫描期间采集的,即与对象的图像数据同时采集,但是在探测到对象的图像数据的探测器区域外被探测到。
简而言之,为了校正,用作均等切片归一化的基础的探测器测量结果的测得的数据来自具有未被衰减的X射线辐射的区域。术语“未被衰减”是指基本上未被衰减的X射线辐射,即辐射不通过诸如患者的任何相关对象。当然,辐射可以穿过基本上X射线透明材料,例如用作壳体或覆盖物的,衣服和支撑表面,并且当然也可以穿过空气。
术语“均等”是指针对两个能量范围(即,针对第一X射线能量范围和第二X射线能量范围)应用校正形式的相同的归一化。
在范例中,提供了一种用于计算机断层摄影X射线成像的系统。所述系统包括数据接口和处理单元。所述数据接口被配置为提供至少第一X射线能量范围和第二X射线能量范围的至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据,所述第一X射线能量范围与所述第二X射线能量范围彼此不同。处理单元被配置为确定用于切片归一化的校正。所述处理单元还被配置为对第一CT X射线投影数据和第二CT X射线投影数据应用均等切片归一化,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据。为了校正,均等切片归一化基于不受对象影响的探测器测量结果的测得的数据,例如来自外面的探测器元件的数据。此外,数据接口被配置为提供经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。
由此,执行了切片归一化,例如以补偿管通量波动。有利地,针对不同能量范围的切片归一化处理导致测量之间的较小差异,并且导致在代表校正或归一化投影数据重建图像时改进的结果。其还改进以减少由例如外面的探测器元件受到患者或系统内任何其他对象影响而引起的伪影。由于提供了相同的归一化,因此需要检查是否探测到对象影响,这可能导致需要关闭切片归一化过程,这将导致经准备的数据突然变化,基于本发明的过程较不容易出错并且提供的数据是更一致的数据并且伪影减少。
术语“切片”是指多个采集的数据,其作为CT中的相邻切片被采集以理想地覆盖感兴趣区域的整个体积。结果是3D数据集,3D数据集然后用于投影以实现可视化。
术语“归一化”是指图像处理中的像素强度值的范围改变时的处理。例如,范围被改变,使得像素强度值与针对预定辐射衰减的预定标准像素强度值对齐。在另一个范例中,范围被改变以改善进一步的图像处理过程,例如,可以通过归一化来增强对比度。
术语“切片归一化”是指CT中测量的切片的归一化,即在重建以进行可视化之前的测得的数据。“预定”值可以被提供为标准化值。因此,归一化使得能够以优化的方式改变测得的像素强度值的范围,以允许信息内容的改进的使用,即用于CT中的重建过程。“预定”值可以是针对特定系统设置或针对特定检查或介入设置(类型)或针对CT中的特定类型的重建而定义的值。“预定”值也可以是通常使用的标准化的值。
根据范例,用作均等切片归一化的基础的由探测器测量结果提供的测得的数据来自对象投影外部的探测器的区域,使得X射线辐射直接到达探测器。
根据范例,测得的数据由外面的探测器元件提供。
术语“外面的探测器元件”是指布置在探测器的主要(或中心)部分之外的探测器元件,所述中心部分被提供用于探测透过患者/对象的辐射。术语“外面的”是相对于被提供用于对象的投影的探测器部分而言的。根据探测器的大小,“外面的”探测器元件可能位于外围设备上(当几乎所有的探测器表面都用于图像数据采集时)。如果只使用探测器的较小部分来对对象进行投影,则“外面的”探测器元件可以是探测器表面区域内的探测器区域。
根据范例,所述系统还包括计算机断层摄影X射线成像采集装置,所述计算机断层摄影X射线成像采集装置具有被配置为生成X射线束的X射线源和被配置为能量分辨X射线探测器的X射线探测器以同时提供分别针对至少两个不同的X射线能量范围的X射线辐射投影数据。所述计算机断层摄影X射线成像采集装置被配置为采集针对所述至少第一X射线能量范围和第二X射线能量范围的所述对象的感兴趣区域的至少第一CT X射线投影数据和第二CT X射线投影数据。
根据范例,所述能量分辨X射线探测器是具有第一探测器层和第二探测器层的双层X射线探测器,其中,所述第一探测器层和所述第二探测器层相对于所述X射线源被彼此前后地布置,并且所述第一探测器层和所述第二探测器层被配置为同时提供针对所述至少两个不同的X射线能量范围的所述X射线投影数据。
根据范例,所述能量分辨X射线探测器是针对不同能量对光子进行计数的光子计数探测器。
根据本发明,还提供了一种用于计算机断层摄影X射线成像的方法。所述方法包括以下步骤:
a)采集具有彼此不同的至少第一X射线能量范围和第二X射线能量范围的对象的感兴趣区域的CT X射线投影数据;其中,所述对象的所述至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据包括来自不同角度的X射线测量结果的多个切片,以产生所述对象的截面图像。
b)确定针对X射线测量结果的所述多个切片的切片归一化的校正因子,以关于所述第一CT X射线投影数据和所述第二CT X射线投影数据的像素强度值来对它们进行改变;
c)对所述第一CT X射线投影数据和所述第二CT X射线投影数据应用均等切片归一化,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据;并且
d)提供所述经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。
其中,对于步骤b),均等切片归一化基于来自探测如下的X射线辐射的探测器区域的探测器测量结果的测得的数据,所述X射线辐射在成像期间未受被布置在X射线源与X射线探测器之间的辐射路径中的所述对象影响。
根据一个范例,提供了一种用于计算机断层摄影X射线成像的方法。所述方法包括以下步骤:
a)采集具有彼此不同的至少第一X射线能量范围和第二X射线能量范围的对象的感兴趣区域的CT X射线投影数据;
b)确定针对切片归一化的校正因子;
c)对所述第一CT X射线投影数据和所述第二CT X射线投影数据应用均等切片归一化,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据;并且
d)提供所述经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。
其中,对于步骤b),均等切片归一化基于不受对象影响的探测器测量结果的测得的数据,例如外面的探测器元件。
根据范例,用作均等切片归一化的基础的由探测器测量结果提供的测得的数据来自对象投影外部的探测器的区域,使得X射线辐射直接到达探测器。
根据范例,所述测得的数据由外面的探测器元件提供。
根据范例,在步骤a)之前,提供a1)利用能量分辨X射线探测器采集具有至少第一X射线能量范围和第二X射线能量范围的对象的感兴趣区域的CT X射线投影数据的步骤。
根据范例,对于步骤b),针对每个能量水平范围提供外面的探测器元件的测得的数据,并且针对切片归一化的校正确定平均值。
根据另一个范例,对于步骤b),针对每个能量水平范围提供外面的探测器元件的测得的数据,并且针对切片标准化的校正确定第一能量水平范围或第二能量水平范围的测量值。
根据一个方面,针对投影图像数据的两个(或更多个)能量范围提供均等归一化。实现两个探测器层的一致处理的一个选择是使用上部和下部探测器层的外面的探测器元件的测得的数据的平均来确定针对两个层的切片归一化-重新标度。另一种选择是仅使用两个层中的一个的测得的数据来确定针对两个层的切片归一化-重新标度。由于X射线源的不同X射线强度输出,形成归一化基础的测得的数据可能会有所不同。由于数据是在辐射通过患者的辐射部分之外的区域中测量的,因而可以采集测量的辐射并将其与标准化强度值相关。得到的归一化因此考虑了来自X射线管的辐射的强度的可能变化,例如,由X射线管的通量强度变化引起的,并且得到用于重建的改进的数据。
参考下文所描述的实施例,本发明的这些和其他方面将变得显而易见并将得以阐述。
附图说明
下面将参考附图来描述本发明的示范性实施例:
图1示出了用于计算机断层摄影X射线成像的系统的示意性概览。
图2示出了具有计算机断层摄影X射线成像采集装置的系统的另一范例。
图3示出了用于计算机断层摄影X射线成像的方法的范例的步骤。
具体实施方式
图1示出了用于计算机断层摄影X射线成像的系统10。所述系统10包括数据接口12和处理单元14。数据接口12被配置为提供对象的针对至少第一X射线能量范围和第二X射线能量范围的至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据。所述第一X射线能量范围与所述第二X射线能量范围彼此不同。所述对象的所述至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据包括来自不同角度的X射线测量结果的多个切片,以产生所述对象的截面图像。所述处理单元14被配置为确定针对X射线测量结果的多个切片归一化的校正,以关于所述第一CT X射线投影数据和所述第二CT X射线投影数据的像素强度值来对它们进行改变。所述处理单元14还被配置为应用针对所述第一CT X射线投影数据和所述第二CT X射线投影数据的均等切片归一化,其中,所述切片归一化中的X射线测量结果中的所述多个切片的所述像素强度值是能改变的,使得提供对预定值的参考,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据。为了校正,所述均等切片归一化基于来自探测如下的X射线辐射的探测器区域的探测器测量的测得的数据,所述X射线辐射在成像期间未受被布置在X射线源与X射线探测器之间的辐射路径中的所述对象影响。此外,数据接口12被配置为提供经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。
缩写CT是指计算机断层摄影X射线成像。
术语“彼此不同”是指彼此独立的范围或值,或者是相邻的或者部分重叠的范围,或者范围是完全重叠的,但是对范围的不同区域具有不同的贡献。
进一步的处理(例如图像处理)可以是用于生成二维或三维CT X射线图像的图像重建步骤。图像因此根据校正的投影数据而被重建。
术语“校正”是指应用到投影数据的校正。所述校正可以被提供为校正值或校正因子或甚至校正项。校正可以应用于对数强度或强度。
在范例中,用作均等切片归一化的基础的由探测器测量结果提供的测得的数据来自对象投影外部的探测器的区域,使得X射线辐射直接到达探测器。
在一个选项中,测得的数据由外面的探测器元件提供。
在范例中,为了校正,提供了与测得的强度相乘的校正因子。
在范例中,为了校正,提供了被加到对数强度或被从对数强度中减去的校正值。
在范例中,为了校正,针对每个“视图”单独地确定切片归一化校正,其中,视图是对于针对源和探测器的相同位置和取向所采集的至少两个不同能量范围的至少两个投影。
图2示出了还包括计算机断层摄影X射线成像采集装置20的范例,所述计算机断层摄影X射线成像采集装置具有被配置为生成X射线束24的X射线源22和被配置为能量分辨X射线探测器的X射线探测器26以同时提供分别针对至少两个不同的X射线能量范围的X射线辐射投影数据。所述计算机断层摄影X射线成像采集装置20被配置为采集至少针对第一X射线能量范围和第二X射线能量范围的对象的感兴趣区域的至少第一CT X射线投影数据和第二CT X射线投影数据。
X射线源22和X射线探测器26被可旋转地安装到机架28。在相对于彼此的固定关系中,它们可以围绕机架旋转,其由箭头20所指示。在范例中,X射线源提供锥形X射线束或扇形X射线束。取决于X射线束,提供相应的探测器元件的阵列或线探测器。
在扇形束的情况下,术语“切片”是指以切片形式(由于扇形束辐射对象而产生的)形式采集的数据。
在锥形束的情况下,代替切片,提供X射线衰减图像数据,其然后被用于重建。因此,术语“切片”还涉及用锥形X射线辐射采集的原始图像。
数据接口12和处理单元14数据连接到由线连接表示的X射线成像采集装置20。当然,这种连接可以是有线连接,或者也可以是无线连接。没有显示进一步的供应,例如用于控制系统的能量供应或数据供应。值得注意的是,数据接口12和处理单元14当然也连接到X射线探测器26。
术语“能量分辨”是指能够探测具有不同能量的X射线辐射并同时提供该数据的探测器。所述对象因此利用包含不同范围能量的X射线辐射被辐射。
能量分辨X射线探测器可以是具有第一探测器层和第二探测器层的双层X射线探测器。所述第一探测器层和所述第二探测器层相对于所述X射线源被彼此前后地布置,并且所述第一探测器层和所述第二探测器层被配置为同时提供针对所述至少两个不同的X射线能量范围的所述X射线投影数据。
能量分辨X射线探测器可以作为例如多于两层的多层探测器提供,例如作为三层探测器。
在另一个范例中,所述能量分辨X射线探测器是针对不同能量对光子进行计数的光子计数探测器。
图3示出了用于计算机断层摄影X射线成像的方法100的步骤的范例。在也被称为步骤a)的第一步骤102中,采集具有彼此不同的至少第一X射线能量范围和第二X射线能量范围的对象的感兴趣区域的CT X射线投影数据。所述对象的所述至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据包括来自不同角度的X射线测量结果的多个切片,以产生所述对象的截面图像。在也被称为步骤b)的第二步骤104中,确定用于X射线测量结果的多个切片的切片归一化的校正因子以根据其像素强度值来改变第一CT X射线投影数据和第二CT X射线投影数据。在也被称为步骤c)的第三步骤106中,对第一CT X射线投影数据和第二CT X射线投影数据应用均等切片归一化,并由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据。在也被称为步骤d)的第四步骤108中,提供经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据用于进一步处理。对于步骤b),均等切片归一化基于来自探测如下的X射线辐射的探测器区域的探测器测量结果的测得的数据,所述X射线辐射在成像期间未受布置在X射线源与X射线探测器之间的辐射路径中的对象影响。
在范例中,用作均等切片归一化的基础的由探测器测量结果提供的测得的数据来自对象投影外部的探测器的区域,使得X射线辐射直接到达探测器;
在一个选项中,测得的数据由外面的探测器元件提供。
作为选项,由散列线指示,在步骤a)之前,提供a1)初始步骤,利用能量分辨X射线探测器采集具有至少第一X射线能量范围和第二X射线能量范围的对象的感兴趣区域的CTX射线投影数据。
在范例中,在步骤a)中,通过双层X射线探测器采集X射线投影数据,所述双层X射线探测器在第一层中提供具有第一X射线能量范围的第一X射线投影数据,并且在双层X射线探测器的第二层中提供具有第二X射线能量范围的X射线投影数据。
在范例中,在步骤d)之后,提供图像重建的步骤e)以提供对象的感兴趣区域的三维X射线图像。
根据范例(未详细示出),对于步骤b),针对每个能量水平范围提供外面的探测器元件的测得的数据,并且针对切片归一化的校正确定平均值。
在双层的范例中,测量第一层和第二层中的外面的探测器元件的数据。
根据另一个范例(也没有进一步详细示出),对于步骤b),针对每个能量水平范围提供外面的探测器元件的测得的数据,并且针对切片标准化的校正确定第一能量水平范围或第二能量水平范围的测量值。
根据范例,对于步骤c),提供了与测得的强度相乘的校正因子。
根据范例,对于步骤c),提供了被添加到对数强度或从对数强度中减去的校正值。
在范例中,对于步骤c),针对每个“视图”单独地确定切片归一化校正,其中,视图是对于针对源和探测器的相同位置和取向所采集的至少两个不同能量范围的至少两个投影。
在本发明的另一示范性实施例中,提供了一种计算机程序或计算机程序单元,其特征在于,其适于在合适的系统上执行根据前述实施例中的一个的方法的方法步骤。
计算机程序单元因此可以存储在计算单元上,其也可以是本发明的实施例的部分。该计算单元可以适于执行上述方法的步骤或引起上述方法的步骤的执行。此外,其可以适于操作上述装置的部件。所述计算单元可以适于自动地操作和/或执行用户的命令。计算机程序可被加载到数据处理器的工作存储器中。数据处理器因此可以被装备为实施本发明的方法。
本发明的该范例性实施例覆盖了从最开始使用本发明的计算机程序和借助于更新将现有程序转变为使用本发明的程序的计算机程序。
进一步,计算机程序单元可以能够提供用以实现如上所述方法的范例性实施例的过程的所有必要步骤。
根据本发明的另一示范性实施例,提出了一种计算机可读介质,诸如CD-ROM,其中,所述计算机可读介质具有存储在其上的计算机程序单元,所述计算机程序元件由前一部分所描述。
计算机程序可以存储和/或分布在适合的介质上,例如与其他硬件一起被提供或作为其他硬件的部分被提供的光学存储介质或固态介质,但是计算机程序也可以以其他形式分布,例如经由因特网或其他的有线或无线的电信系统分布。
然而,计算机程序也可以通过如万维网的网络来提供并且可以被从这样的网络下载到数据处理器的工作存储器中。根据本发明的另外的示范性实施例,提供了一种用于使得计算机程序单元可供下载的介质,所述计算机程序单元被布置为执行本发明的先前描述的实施例中的一个。
必须指出,本发明的实施例参考不同主题进行描述。具体地,参考方法类型权利要求描述了一些实施例,而参考装置类型权利要求描述了其他实施例。然而,本领域技术人员以上和以下描述可以得出,除非另行指出,除了属于同一类型的主题的特任的任何组合之外,涉及不同主题的特征之间的任何组合也被认为由本申请公开。然而,所有特征能够被组合,提供超过所述特征的简单加和的协同效应。
尽管已经在附图和前面的描述中详细图示和描述了本发明,但是这样的图示和描述应当被认为是说明性或示范性的,而非限制性的。本发明不限于所公开的实施例。本领域技术人员通过研究附图、公开内容以及从属权利要求,在实践请求保护的本发明时能够理解并且实现对所公开的实施例的其他变型。
在权利要求书中,词语“包括”不排除其他元件或步骤,并且词语“一”或“一个”不排除多个。单个处理器或其他单元可以完成权利要求书中所记载的若干个项目的功能。尽管在互不相同的从属权利要求中记载了特定措施,但是这并不指示不能有利地使用这些措施的组合。权利要求书中的任何附图标记均不应被解释为对范围的限制。
Claims (13)
1.一种用于计算机断层摄影X射线成像的系统(10),所述系统包括:
-数据接口(12);以及
-处理单元(14);
其中,所述数据接口被配置为提供对象的针对至少第一X射线能量范围和第二X射线能量范围的至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据,所述第一X射线能量范围与所述第二X射线能量范围彼此不同;
其中,所述对象的所述至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据包括来自不同角度的X射线测量结果的多个切片,以产生所述对象的截面图像;
其中,所述处理单元被配置为确定针对X射线测量结果的所述多个切片的归一化的校正,以关于所述第一CT X射线投影数据和所述第二CT X射线投影数据的像素强度值来对它们进行改变;并且针对所述第一CT X射线投影数据和所述第二CT X射线投影数据应用均等切片归一化,其中,所述切片归一化中的X射线测量结果中的所述多个切片的所述像素强度值是能改变的,使得提供对预定值的参考,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据;其中,为了所述校正,所述均等切片归一化基于来自探测如下的X射线辐射的探测器区域的探测器测量结果的测得的数据:所述X射线辐射在成像期间未受被布置在X射线源与X射线探测器之间的辐射路径中的所述对象影响;并且
其中,所述数据接口被配置为提供所述经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理。
2.根据权利要求1所述的系统,其中,用作所述均等切片归一化的基础的由所述探测器测量结果提供的所述测得的数据得自所述对象的投影外部的探测器的区,使得所述X射线辐射直接到达所述探测器;并且
其中,优选地,所述测得的数据由外面的探测器元件提供。
3.根据权利要求1或2所述的系统,其中,所述系统还包括计算机断层摄影X射线成像采集装置(20),所述计算机断层摄影X射线成像采集装置具有:
-X射线源(22),其被配置为生成X射线束,以及
-X射线探测器(26),其被配置为能量分辨X射线探测器,以同时提供分别针对至少两个不同X射线能量范围的X射线辐射投影数据;
其中,所述计算机断层摄影X射线成像采集装置被配置为采集所述对象的感兴趣区域的针对至少所述第一X射线能量范围和所述第二X射线能量范围的至少所述第一CT X射线投影数据和所述第二CT X射线投影数据。
4.根据权利要求3所述的系统,其中,所述能量分辨X射线探测器是具有第一探测器层和第二探测器层的双层X射线探测器,其中,所述第一探测器层和所述第二探测器层相对于所述X射线源被彼此前后地布置,并且所述第一探测器层和所述第二探测器层被配置为同时提供针对至少两个不同的X射线能量范围的所述X射线投影数据。
5.根据权利要求3所述的系统,其中,所述能量分辨X射线探测器是针对不同能量对光子进行计数的光子计数探测器。
6.一种用于计算机断层摄影X射线成像的方法(100),所述方法包括以下步骤:
a)采集(102)对象的感兴趣区域的具有彼此不同的至少第一X射线能量范围和第二X射线能量范围的CT X射线投影数据;其中,所述对象的所述至少第一CT X射线辐射投影数据和第二CT X射线辐射投影数据包括来自不同角度的X射线测量结果的多个切片,以产生所述对象的截面图像;
b)确定(104)针对X射线测量结果的所述多个切片的切片归一化的校正因子,以关于所述第一CT X射线投影数据和所述第二CT X射线投影数据的像素强度值来对它们进行改变;
c)对所述第一CT X射线投影数据和所述第二CT X射线投影数据应用(106)均等切片归一化,并且由此生成经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据;并且
d)提供(108)所述经准备的第一CT X射线投影数据和经准备的第二CT X射线投影数据以供进一步处理;
其中,对于步骤b),均等切片归一化基于来自探测如下的X射线辐射的探测器区域的探测器测量结果的测得的数据,所述X射线辐射在成像期间未受被布置在X射线源与X射线探测器之间的辐射路径中的所述对象影响。
7.根据权利要求6所述的方法,其中,用作所述均等切片归一化的基础的由所述探测器测量结果提供的所述测得的数据得自所述对象的投影外部的探测器的区,使得所述X射线辐射直接到达所述探测器;并且
其中,优选地,所述测得的数据由外面的探测器元件提供。
8.根据权利要求6或7所述的方法,其中,在步骤a)之前提供以下步骤:
-利用能量分辨X射线探测器采集(110)所述对象的感兴趣区域的具有至少所述第一X射线能量范围和所述第二X射线能量范围的CT X射线投影数据。
9.根据权利要求7或8所述的方法,其中,对于步骤b),针对每个能量水平范围提供所述外面的探测器元件的测得的数据,并且针对所述切片归一化的所述校正来确定平均值。
10.根据权利要求7、8或9所述的方法,其中,对于步骤b),针对每个能量水平范围提供所述外面的探测器元件的测得的数据,并且针对所述切片标准化的所述校正来确定所述第一能量水平范围或所述第二能量水平范围的测得的值。
11.根据权利要求6、7、8、9或10所述的方法,其中,对于步骤c):
i)提供与测得的强度相乘的校正因子;并且/或者
ii)提供被加到对数强度或被从对数强度减去的校正值。
12.一种用于控制根据权利要求1至5中的一项所述的装置的计算机程序单元,所述计算机程序单元当由处理单元运行时适于执行根据权利要求6至11所述的方法的步骤。
13.一种存储有根据权利要求12所述的程序单元的计算机可读介质。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP16187966 | 2016-09-09 | ||
EP16187966.3 | 2016-09-09 | ||
PCT/EP2017/072039 WO2018046414A1 (en) | 2016-09-09 | 2017-09-04 | Computer tomography x-ray imaging |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108135557A true CN108135557A (zh) | 2018-06-08 |
CN108135557B CN108135557B (zh) | 2019-04-30 |
Family
ID=56896398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780003448.0A Active CN108135557B (zh) | 2016-09-09 | 2017-09-04 | 计算机断层摄影x射线成像 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10485501B2 (zh) |
EP (1) | EP3346922B1 (zh) |
JP (1) | JP6641065B2 (zh) |
CN (1) | CN108135557B (zh) |
WO (1) | WO2018046414A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112969411A (zh) * | 2018-10-04 | 2021-06-15 | 皇家飞利浦有限公司 | 用于提供谱图像的系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1284659A (zh) * | 1999-08-17 | 2001-02-21 | 高级光学技术股份有限公司 | 用于双能量x射线成像的设备和方法 |
CN1748648A (zh) * | 2004-09-14 | 2006-03-22 | 株式会社东芝 | X射线计算机断层摄影装置 |
US20090304142A1 (en) * | 2008-06-09 | 2009-12-10 | Arineta Ltd. | Ct scanner with scatter radiation correction and method of using same |
CN101810488A (zh) * | 2009-01-29 | 2010-08-25 | 皇家飞利浦电子股份有限公司 | 探测值校正设备 |
WO2015011587A1 (en) * | 2013-07-23 | 2015-01-29 | Koninklijke Philips N.V. | Hybrid (spectral/non-spectral) imaging detector array and corresponding processing electronics |
US20150178958A1 (en) * | 2013-12-20 | 2015-06-25 | Toshiba Medical Systems Corporation | X ray computer tomography apparatus and image processing apparatus |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5761257A (en) * | 1996-12-02 | 1998-06-02 | General Electric Company | Normalizing projection data in a computed tomography system |
JP4163279B2 (ja) * | 1998-02-27 | 2008-10-08 | 株式会社東芝 | データ収集・処理装置およびこの装置を用いたシステム |
JP2003126078A (ja) * | 2001-10-23 | 2003-05-07 | Hitachi Medical Corp | X線ct装置 |
JP5329204B2 (ja) * | 2008-12-19 | 2013-10-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US8699812B2 (en) * | 2011-06-09 | 2014-04-15 | General Electric Company | System and method for quality improvement in CT imaging |
US8926177B2 (en) * | 2012-05-21 | 2015-01-06 | General Electric Company | Source side monitoring device for an imaging system |
CN105142526A (zh) | 2013-04-16 | 2015-12-09 | 株式会社东芝 | X射线ct装置 |
WO2014209972A1 (en) | 2013-06-24 | 2014-12-31 | Massachusetts Institute Of Technology | Normalization correction for multiple-detection enhanced emission tomography |
WO2015017851A1 (en) | 2013-08-02 | 2015-02-05 | Matthews Dawn C | Method for analyzing and correcting measurement variability in pet images |
US9478049B2 (en) | 2013-09-30 | 2016-10-25 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3D iterative CT reconstruction |
US20160206255A1 (en) * | 2015-01-16 | 2016-07-21 | Kabushiki Kaisha Toshiba | Hybrid passive/active multi-layer energy discriminating photon-counting detector |
US9482630B2 (en) | 2015-01-21 | 2016-11-01 | Toshiba Medical Systems Corporation | Multiple-layered energy-integrating detector in a hybrid computed tomography scanner |
-
2017
- 2017-09-04 JP JP2019512865A patent/JP6641065B2/ja active Active
- 2017-09-04 EP EP17765381.3A patent/EP3346922B1/en active Active
- 2017-09-04 US US15/762,114 patent/US10485501B2/en active Active
- 2017-09-04 WO PCT/EP2017/072039 patent/WO2018046414A1/en active Application Filing
- 2017-09-04 CN CN201780003448.0A patent/CN108135557B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1284659A (zh) * | 1999-08-17 | 2001-02-21 | 高级光学技术股份有限公司 | 用于双能量x射线成像的设备和方法 |
CN1748648A (zh) * | 2004-09-14 | 2006-03-22 | 株式会社东芝 | X射线计算机断层摄影装置 |
US20090304142A1 (en) * | 2008-06-09 | 2009-12-10 | Arineta Ltd. | Ct scanner with scatter radiation correction and method of using same |
CN101810488A (zh) * | 2009-01-29 | 2010-08-25 | 皇家飞利浦电子股份有限公司 | 探测值校正设备 |
WO2015011587A1 (en) * | 2013-07-23 | 2015-01-29 | Koninklijke Philips N.V. | Hybrid (spectral/non-spectral) imaging detector array and corresponding processing electronics |
US20150178958A1 (en) * | 2013-12-20 | 2015-06-25 | Toshiba Medical Systems Corporation | X ray computer tomography apparatus and image processing apparatus |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112969411A (zh) * | 2018-10-04 | 2021-06-15 | 皇家飞利浦有限公司 | 用于提供谱图像的系统 |
Also Published As
Publication number | Publication date |
---|---|
US20190183443A1 (en) | 2019-06-20 |
CN108135557B (zh) | 2019-04-30 |
WO2018046414A1 (en) | 2018-03-15 |
JP2019528861A (ja) | 2019-10-17 |
JP6641065B2 (ja) | 2020-02-05 |
EP3346922A1 (en) | 2018-07-18 |
US10485501B2 (en) | 2019-11-26 |
EP3346922B1 (en) | 2019-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2606561C2 (ru) | Компенсация усечения для итерационной реконструкции в компьютерной томографии (кт) с коническим пучком в комбинированных системах офэкт/кт | |
CN105452852B (zh) | 通过使用相衬ct的单色衰减对比图像生成 | |
US9615812B2 (en) | Calibration of imagers with dynamic beam shapers | |
CN104856708B (zh) | 选择辐射形状过滤器的方法和x射线成像系统 | |
JPWO2010061810A1 (ja) | 放射線撮像装置 | |
US10605933B2 (en) | X-ray spectral calibration technique for cone-beam CT | |
US20110060566A1 (en) | Method and apparatus for scatter correction | |
US8290116B2 (en) | Imaging apparatus including correction unit for scattered radiation | |
EP3170148B1 (en) | Iterative reconstruction method for spectral, phase-contrast imaging | |
CN110072459A (zh) | 用于自校准的自校准ct检测器、系统和方法 | |
CN105916283B (zh) | 自动调整用于x射线检测器的信号分析参数 | |
KR20190085740A (ko) | 단층 촬영 장치, 그 제어 방법, 및 컴퓨터 프로그램 제품 | |
CN106880902A (zh) | 基于单能量图像记录确定材料属性值的空间分布 | |
JP5696778B2 (ja) | ポジトロンct装置および再構成方法 | |
JP2017127638A (ja) | X線ct装置、情報処理装置、および情報処理方法 | |
JP6339684B2 (ja) | フォトンカウンティングct装置および推定被ばく量算出方法 | |
CN106028938A (zh) | X射线ct装置以及拍摄方法 | |
KR102605456B1 (ko) | 치아 이미지 생성 시스템으로부터의 이미지 데이터를 개선하는 방법 | |
CN108135557B (zh) | 计算机断层摄影x射线成像 | |
CN114650774B (zh) | 用于处理由暗场和/或相衬x射线成像系统采集的数据的装置 | |
US10755450B2 (en) | Display of depth location of computed tomography slice images relative to an object to be imaged | |
Ghani et al. | An empirical method for geometric calibration of a photon counting detector-based cone beam CT system | |
US20070133750A1 (en) | Method for operation of a scattered-ray grid | |
NL1033544C2 (nl) | Werkwijze en systeem voor radiografische afbeelding met orgaan-gebaseerd stralingsprofielvoorschrift. | |
JP5268248B2 (ja) | ディジタル及びコンピューティッド・ラジオグラフィ画像についての適応画像処理及び表示法 |
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 |