CN106456094B - 用于生成对象的图像的成像系统 - Google Patents
用于生成对象的图像的成像系统 Download PDFInfo
- Publication number
- CN106456094B CN106456094B CN201580029374.9A CN201580029374A CN106456094B CN 106456094 B CN106456094 B CN 106456094B CN 201580029374 A CN201580029374 A CN 201580029374A CN 106456094 B CN106456094 B CN 106456094B
- Authority
- CN
- China
- Prior art keywords
- projection data
- spectral
- spectrally
- estimation
- image
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 72
- 230000003595 spectral effect Effects 0.000 claims abstract description 191
- 238000001514 detection method Methods 0.000 claims abstract description 62
- 239000000463 material Substances 0.000 claims abstract description 55
- 238000005259 measurement Methods 0.000 claims abstract description 54
- 238000009826 distribution Methods 0.000 claims abstract description 48
- 230000005855 radiation Effects 0.000 claims description 61
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 12
- 238000002591 computed tomography Methods 0.000 abstract description 15
- 238000013170 computed tomography imaging Methods 0.000 abstract description 3
- 238000000701 chemical imaging Methods 0.000 abstract description 2
- 239000010410 layer Substances 0.000 description 18
- 238000000034 method Methods 0.000 description 10
- 239000002355 dual-layer Substances 0.000 description 8
- 210000000988 bone and bone Anatomy 0.000 description 7
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 6
- 229910052791 calcium Inorganic materials 0.000 description 6
- 239000011575 calcium Substances 0.000 description 6
- 238000005316 response function Methods 0.000 description 5
- 230000001419 dependent effect Effects 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 208000006386 Bone Resorption Diseases 0.000 description 3
- 230000024279 bone resorption Effects 0.000 description 3
- 239000002872 contrast media Substances 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000011896 sensitive detection Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 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/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/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/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/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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明涉及用于生成对象的图像的成像系统(17)，如计算机断层摄影系统。由具有谱探测元件和非谱探测元件的探测器(6)生成谱测量投影数据和非谱测量投影数据，并且通过使用能够已经引起所述非谱测量投影数据的模型材料分布并且通过基于所述模型材料分布模拟谱估计投影数据的测量结果，来估计谱估计投影数据。图像是基于测量和估计的谱投影数据重建的。除所述谱测量投影数据之外，使用所述谱估计投影数据可以得到高质量谱成像，尤其是高质量谱计算机断层摄影成像，其使用不仅具有通常更复杂的谱探测元件而且具有较简单的非谱探测元件的简化探测器。
Description
技术领域
本发明涉及用于生成对象的图像的成像系统、成像方法和成像计算机程序。本发明具体而言涉及用于生成对象的图像的谱计算机断层摄影系统、方法和计算机程序。
背景技术
US 2007/0147574 A1公开了一种用于采集包括能量积分和能量鉴别数据测量结果的图像数据集的系统。该系统适于在采集周期期间获得能量积分测量数据和能量鉴别测量数据并且以选定的方式组合和重建能量积分测量数据和能量鉴别测量数据以生成能量鉴别分量图像。
US 2013/0101156 A1公开了一种包括具有带至少一个谱通道的至少一个模态的无创成像设备的装置。该装置还包括：存储器，其已经在其中存储了用于不同材料和可行性准则的多个模型；和控制电路，其被配置为处理如由无创成像设备所提供的针对对象的成像信息，其中，多个模型被用于识别针对成像信息的各部分的候选材料，并且其中，可行性准则被用于通过避免不可能的材料或材料的组合来减少候选材料。
R.Carmi等人的文章“Material Separation with Dual-Layer CT”IEEE NuclearScience Symposium Conference Record，第4卷、第1876至1878页(2005)公开了一种谱计算机断层摄影系统，其使用用于生成谱投影数据的双层探测器，以及用于基于所生成的谱投影数据重建诸如人的对象的图像的重建单元。双层探测器使用具有在堆叠中的两个敏感探测层的技术复杂的配置以便提供能量分离。
发明内容
本发明的目标是提供一种用于对对象成像的成像系统、方法和计算机程序，其允许基于谱投影数据对对象的图像的生成，其中，可以使用技术上较不复杂的探测器。
在本发明的第一方面中，呈现了一种用于生成对象的图像的成像系统，其中，所述成像系统包括：
-测量投影数据提供单元，其用于提供谱测量投影数据和非谱测量投影数据，其中，已经基于已经穿过包括所述对象的成像区域的辐射由具有谱探测元件和非谱探测元件的探测器分别生成了所述谱测量投影数据和所述非谱测量投影数据，
-估计单元，其用于估计对应于结束于所述非谱探测元件处的辐射路径的谱估计投影数据，其中，为了估计针对结束于非谱探测元件处的辐射路径的谱估计投影数据值，所述估计单元适于：a)基于所述非谱测量投影数据来确定所述成像区域内的模型材料分布，所述模型材料分布能够已经引起对应于所述辐射路径的非谱测量投影数据值；并且b)基于所述模型材料分布来模拟所述谱估计投影数据值的测量结果，
-重建单元，其用于基于所述谱测量投影数据和所述谱估计投影数据来重建所述对象的图像。
由于所述成像系统可以基于所提供的谱测量投影数据和非谱测量投影数据生成所述对象的所述图像，因而被用于生成这些测量投影数据的探测器不需要仅包括在技术上比非谱探测元件更复杂的谱探测元件，从而降低所述探测器的总体技术复杂性。而且，由于所述估计单元可以估计对应于结束于所述非谱探测元件处的辐射路径的谱估计投影数据，其中，这些谱估计投影数据与所提供的谱测量投影数据一起被用于重建所述对象的图像，因此仍然可以生成所述对象的高质量图像。例如，由于所述探测器不仅包括谱探测元件，因此所述谱测量投影数据可能被截断，这可以导致所生成的图像中的图像伪影。该截断和所述对应的图像伪影可以通过不仅使用用于重建所述图像的所述谱测量投影数据而且通过使用所述谱估计投影数据来减轻，尤其是，消除。
所述测量投影数据提供单元可以是投影数据采集系统，其用于通过使用具有谱探测元件和非谱探测元件的探测器来采集所述谱测量投影数据和所述非谱测量投影数据。然而，所述测量投影数据提供单元还可以是仅存储单元，其用于存储所述谱测量投影数据和所述非谱测量投影数据并且用于检索所存储的谱测量投影数据和非谱测量投影数据以用于提供其。所述测量投影数据提供单元还可以是接收单元，其用于从投影数据采集系统接收所述谱测量投影数据和非谱测量投影数据并且用于提供所接收的谱测量投影数据和非谱测量投影数据。
所述对象优选地是如人或动物的活对象。然而，所述对象还可以是技术对象。例如，所述成像系统可以出于安全原因被用于生成机场的行李的图像。
在优选的实施例中，为了确定用于估计谱估计投影数据值的所述模型材料分布，所述估计单元适于基于对应的非谱投影数据值，来确定沿着对应于要被确定的谱估计投影数据值的辐射路径的模型材料分布。具体而言，所述估计单元适于估计水分布作为沿着所述辐射路径的所述模型材料分布。例如，对于非谱测量投影数据值而言，等价水长度可以基于相应的非谱测量投影数据值来确定。所述等价水长度然后可以被用于通过模拟沿着所述等价水长度的辐射的行进来确定所述谱估计投影数据。这允许基于所述非谱测量投影数据和投影数据生成参数相对快速地估计谱估计投影数据。
所述谱估计投影数据值的测量结果的所述模拟优选地是模拟使用还已经被用于测量谱测量投影数据的一种探测器元件的测量结果的模拟。因此，其优选地考虑被用于测量所述谱测量投影数据的所述探测技术。例如，如果被用于测量所述谱测量投影数据的所述探测器元件是具有在堆叠中的两个敏感探测层的双层探测器元件以便提供光子能量分离，则在模拟所述测量过程时，优选地考虑所述层的已知吸收特性。
在另外的实施例中，为了确定用于估计谱估计投影数据值的所述模型材料分布，所述估计单元适于基于所述谱测量投影数据和所述非谱投影数据来重建中间图像，并且基于经重建的中间图像来确定所述材料分布。具体而言，所述估计单元适于将所述谱测量投影数据积分从而生成积分投影数据，并且基于所述积分投影数据和所述非谱测量投影数据来重建所述中间图像。优选地，所述估计单元适于通过分割经重建的中间图像中的不同材料来确定所述材料分布。例如，如果所述对象是人，则可以在经重建的中间图像中分割骨和组织区，并且该分割可以被用于确定沿着相应的辐射路径的钙和水的分布，其中，该钙/水材料分布可以被用于模拟谱估计投影数据值的测量结果，以用于确定其。这允许估计谱估计投影数据的进一步改进的准确度。
还优选的是，所述估计单元适于将所述谱测量投影数据和所述谱估计投影数据在所述谱测量投影数据与所述谱估计投影数据之间的边界处平滑。从所述谱测量投影数据到所述谱估计投影数据的过渡可能不是非常平滑的。通过在该过渡或边界处将这些投影数据平滑，还可以改进所述对象的最终经重建的图像的质量。
在本发明的另一方面中，呈现了一种用于生成对象的图像的成像方法，其中，所述成像方法包括：
-由测量投影数据提供单元提供谱测量投影数据和非谱测量投影数据，其中，已经基于已经穿过包括所述对象的成像区域的辐射由具有谱探测元件和非谱探测元件的探测器分别生成了所述谱测量投影数据和所述非谱测量投影数据，
-由估计单元估计对应于结束于所述非谱探测元件处的辐射路径的谱估计投影数据，其中，为了估计针对结束于非谱探测元件处的辐射路径的谱估计投影数据值，所述估计单元：a)基于所述非谱测量投影数据来确定所述成像区域内的模型材料分布，所述模型材料分布能够已经引起对应于所述辐射路径的非谱测量投影数据值；并且b)基于所述模型材料分布来模拟所述谱估计投影数据值的测量结果，
-由重建单元基于所述谱测量投影数据和所述谱估计投影数据来重建所述对象的图像。
在本发明的由一方面中，呈现了一种用于对感兴趣区域进行成像的计算机程序，其中，所述计算机程序包括程序代码模块，当在控制根据权利要求1所述的成像系统的计算机上运行所述计算机程序时，所述程序代码模块用于令所述成像系统执行根据权利要求9所述的成像方法的步骤。
应该理解，根据权利要求1所述的成像系统、根据权利要求9所述的成像方法和根据权利要求10所述的计算机程序具有相似和/或相同的优选实施例，尤其是如从属权利要求中定义的。
应该理解，本发明的优选实施例也可以是从属权利要求或以上实施例与各自的独立权利要求的任何组合。
本发明的这些和其他方面将根据下文描述的实施例变得显而易见，并且将参考下文描述的实施例得到阐述。
附图说明
在以下附图中：
图1示意性并且示范性地示出了用于生成对象的图像的成像系统的实施例；
图2示意性并且示范性地示出了具有谱探测元件和非谱探测元件的探测器；
图3示意性并且示范性地图示了截断的谱测量投影数据和谱估计投影数据；并且
图4示出了示范性地图示用于生成对象的图像的成像系统的实施例的流程图。
具体实施方式
图1示意性并且示范性地示出了用于生成对象的图像的成像系统的实施例。在该实施例中，成像系统17是用于生成人的计算机断层摄影图像的谱计算机断层摄影系统。计算机断层摄影系统17包括机架1，机架1能够关于平行于z方向延伸的旋转轴R旋转。在该实施例中，作为X射线管的辐射源2被安装在机架1上。辐射源2被提供有准直器3，准直器3在该实施例中根据由辐射源2生成的辐射而形成圆锥形辐射束4。辐射穿过成像区域5内的人(未示出在图1中)，所述成像区域也可以被认为是检查区并且在该实施例中是圆柱形的。在已经穿过成像区域5之后，辐射束4入射在包括二维探测表面的探测器6上。探测器6被安装在机架1上。
计算机断层摄影系统17包括两个电机7、8。由电机7以优选地恒定但是可调节的角速度来驱动机架1。电机8被提供用于将被布置在如平行于旋转轴R或z轴的方向的成像区域5内的台的支持模块上的人位移。例如，这些电机7、8由控制单元9控制，使得辐射源2和人沿着螺旋形轨迹相对于彼此移动。然而，还可能的是，不移动人，而是仅旋转辐射源2，即辐射源2相对于人沿着圆形轨迹移动。而且，在另一实施例中，准直器3可以适于形成另一射束形状(尤其地扇形射束)，并且探测器6可以包括对应于其他射束形状(尤其地扇形射束)整形的探测表面。
如在图2中示意性并且示范性地图示的，在由辐射源2发射的辐射已经穿过人14之后，由探测器6的非谱探测元件15和谱探测元件16来探测辐射。因此，探测器6不仅包括谱探测元件16而且包括非谱探测元件15。在辐射源2和人的相对移动期间，探测器6根据入射在探测器6的探测表面上的辐射，通过使用谱探测元件16来生成谱测量投影数据并且通过使用非谱探测元件15来生成非谱测量投影数据。因此，辐射源2、用于相对于人移动辐射源2的元件(尤其地，电机7、8)和机架1和探测器6形成用于提供谱测量投影数据和非谱测量投影数据的测量投影数据提供单元18，其中，已经根据描述投影数据的生成的投影数据生成参数、基于已经穿过包括人14的成像区域4的辐射4、由具有谱探测元件16和非谱探测元件15的探测器16分别生成了谱测量投影数据和非谱测量投影数据已经基于由。测量投影数据提供单元18还可以被认为是投影数据采集系统。
投影数据生成参数描述投影数据的生成，即其描述例如被用于生成谱测量投影数据和非谱测量投影数据的辐射的采集几何结构和强度。采集几何结构定义关于对象的辐射路径，辐射在由探测器的谱探测元件和非谱探测元件探测到之前沿着所述辐射路径行进。
谱测量投影数据19可以被截断，如在图3中示意性并且示范性地所图示的。图3示范性并且示意性地示出了取决于探测器6上的相应的探测元件d的谱测量投影数据值P。如在图3中可以看到的，在非谱探测元件15与谱探测元件16之间的边界处，谱测量投影数据19被截断。
因此，计算机断层摄影系统17还包括估计单元10，估计单元10用于估计对应于结束于非谱探测元件15处的辐射路径的谱估计投影数据22，其中，为了估计针对结束于非谱探测元件15处的辐射路径的谱估计投影数据值，估计单元10适于：a)基于非谱测量投影数据并且任选地还基于投影数据生成参数，来确定能够已经引起对应于辐射路径的非谱测量投影数据值的成像区域5内的模型材料分布；和b)基于模型材料分布并且任选地还基于投影数据生成参数，来模拟谱估计投影数据值的测量结果。具体而言，估计单元10适于基于对应的非谱投影数据值和投影数据生成参数，来确定沿着对应于要被确定的谱估计投影数据值的辐射路径的模型材料分布。例如，估计单元10适于水分布估计作为沿着辐射路径的模型材料分布，即，等价水长度。备选地，估计单元10可以适于基于谱测量投影数据、非谱投影数据和投影数据生成参数来重建中间图像，并且基于经重建的中间图像来确定材料分布。在这种情况下，估计单元10优选地适于积分谱测量投影数据从而生成积分投影数据，并且适于基于积分投影数据、非谱测量投影数据和投影数据生成参数来重建中间图像。估计单元10可以适于通过分割经重建的中间图像中的不同材料来确定材料分布。例如，可以在经重建的中间图像中分割骨和组织区域以便确定材料分布。
估计单元10还适于将谱探测元件16和非谱探测元件15之间的边界处的谱测量投影数据19和谱估计投影数据22平滑。在图3中，平滑的谱估计投影数据由虚线指示。
计算机断层摄影系统17还包括重建单元11，其用于基于谱测量投影数据、谱估计投影数据和投影数据生成参数来重建人14的图像。在该实施例中，重建单元11适于使用滤波反投影算法，以用于重建人的图像。重建单元11可以适于使用材料分解技术以便基于谱投影数据来重建对应于人内的不同材料的不同的图像。例如，如果造影剂已经注入人中，则可以重建仅示出人内的造影剂的第一图像，并且可以重建示出在没有造影剂的情况下人的第二图像。重建单元11还可以适于基于谱投影数据，来重建对应于诸如康普顿效应和光电效应的不同的物理效应的不同的图像。当然，重建单元11还可以适于使用其他谱重建技术。例如，在E.Roessl和R.Proksa的文章“K-edge imaging in x-ray computed tomographyusing multi-bin photon counting detectors”、Physics in Medicine and Biology、第52卷、第4679到4696页(2007)中公开了可以由重建单元使用的已知重建技术，通过引用将其并入本文。
估计单元10和重建单元11是还可以由控制单元9控制的图像生成设备23的部分。计算机断层摄影系统17还包括如键盘、计算机鼠标和触摸板等的输入单元12和用于显示经重建的图像的显示器13。
在以下中，将参考图4中所示的流程图示范性地描述用于生成对象的图像的成像方法的实施例。
在该实施例中，用于生成对象的图像的成像方法是用于生成的人的图像的谱计算机断层摄影成像方法。在步骤101中，辐射源2围绕旋转轴R旋转并且人14未移动，即辐射源2沿着人14周围的圆形轨迹行进。在另一实施例中，辐射源2可以沿着相对于人的另一轨迹(例如，螺旋轨迹)移动。辐射源2发射穿过人14的辐射，并且已经穿过人14的辐射由探测器6探测到，探测器6生成谱测量投影数据和非谱测量投影数据。步骤101可以被认为是用于提供谱测量投影数据和非谱测量投影数据的投影数据提供步骤。
在步骤102中，通过估计单元10估计对应于结束于非谱探测元件15处的辐射路径的谱估计投影数据22，其中，为了估计针对结束于非谱探测元件15处的辐射路径的谱估计投影数据值，估计单元10：a)基于非谱测量投影数据和投影数据生成参数，来确定能够已经引起对应于所示辐射路径的非谱测量投影数据值的成像区域5内的模型材料分布；并且b)基于模型材料分布和投影数据生成参数，来模拟谱估计投影数据值的测量结果。在步骤103中，人14的图像是基于谱测量投影数据、谱估计投影数据和投影数据生成参数来重建的，其中，在步骤140中，在显示器13上示出经重建的图像。
探测器6的谱探测元件16优选地是双层探测器元件。双层探测器元件具有在堆叠中的两个敏感探测层以便提供光子能量分离。双层探测器元件是公知的并且在例如R.Carmi等人的上述文章中被公开，通过引用将其并入本文。在其他实施例中，探测器可以包括如具有能量鉴别的光子计数探测元件的其他谱探测元件。
谱探测元件的产生与非谱探测元件的产生相比是技术上相对复杂并且因此相对昂贵，所述非谱探测元件还可以被认为是常规探测元件。为了降低复杂度和因此成本，探测器6是混合探测器，其中，探测器的一部分装备有常规的技术上较简单并且较便宜的探测元件，并且另一部分装备有谱探测元件。不同的部分可以非对称地被分布(如在图2中示意性并且示范性地图示的)或者其可以对称地分布，其中，优选地在每种情况下，探测器关于圆周方向包括谱内部分和两个外部非谱部分。
除基于谱测量投影数据和谱估计投影数据来重建图像之外，重建单元11还可以适于重建常规图像，即基于仅非谱投影数据来重建图像，所述常规图像对应于来自常规探测器的图像并且可以通过使用如滤波反投影技术的标准重建技术来重建。具体而言，非谱测量投影数据可以与在谱积分之后还形成非谱投影数据的谱积分的谱测量投影数据一起使用，以便重建常规图像。如果谱探测元件是双层探测元件，来自双层探测元件的两层的读数可以被求和，以便生成非谱投影数据值。
所提供的谱测量投影数据是截断的投影数据，其中，截断可以导致最终经重建的图像中的图像伪影。尤其是在滤波反投影算法的情况下，截断的投影数据能够产生两个问题。高频伪影能够存在于截断界限附近，并且经滤波的投影数据的低频分量可能是不正确的。图像伪影可以导致经重建的图像中的偏置，其中，偏置从临床应用观点看是实质性缺点。基于谱投影数据的谱计算机断层摄影(尤其是基于双能量投影数据的双能量计算机断层摄影)提供生成如量化碘映射的量化临床信息的有价值的机会。经重建的图像中的偏置将使量化退化。因此，上文参考图1和图2所描述的计算机断层摄影系统允许量化误差的校正(尤其是对于混合双能量/常规断层摄影探测器系统)。
估计单元可以适于根据常规探测数据(即，根据非谱测量投影数据)来估计截断区域中的等价水长度，并且根据给定水长度和双层元件的层的已知吸收特性来估计预期的双层信号(即，谱估计投影数据)。所估计和所测量的双层信号(即，所测量和所估计的谱投影数据)可以然后被组合为非截断的投影数据，其可以在截断边界附近平滑以便避免尖锐过渡。然后，组装和平滑的谱投影数据可以在反投影之前被滤波，以便重建图像。
在实施例中，估计单元适于使用以下等式来确定沿着对应于要被确定的谱估计投影数据值的辐射路径的模型材料分布：
其中，I是所测量的探测信号，即，例如对应的非谱投影数据值，I0是在没有吸收对象的情况下的已知探测信号，E是能量，S(E)是整个成像系统的已知谱响应函数，尤其地包括X射线管发射谱、射束滤波、探测器的谱敏感度函数、在能量积分探测元件的情况下的因子E和影响总能量敏感度的所有其他因子，L是相应的射束中的材料的长度(即，沿着相应的辐射路径)，并且μ(E)是材料的已知线性吸收系数。成像系统具有三个不同的谱响应函数，即针对非谱探测元件的第一谱响应函数SN(E)，针对谱探测元件的下层的第二谱响应函数SL(E)和针对谱探测元件的上层的第三谱响应函数SU(E)。
为了确定等价水长度，非谱测量结果IN因此可以由以下来等式建模：
其中，估计单元可以适于针对等价水长度LH2O数字地求解该等式。在等式(2)中，μH2O指代水的线性吸收系数。
谱估计投影数据值是能量相关值，即对于不同的能量或能量分布而言，谱估计投影数据值是不同的。例如，其可以包括：第一值IL，其对应于双层探测器的下层和由下层探测到的对应的能量分布；以及第二值IU，其对应于双层探测器的上层和由上层探测到的对应的能量分布。估计单元还可以适于根据以下等式使用等价水长度LH2O来估计包括IL和IU的谱估计投影数据值(即，模拟包括IL和IU的谱估计投影数据值的测量结果)：
估计单元还可以适于重建常规图像(即，基于非谱投影数据重建的中间图像)，分割经重建的常规图像内的骨区和组织区并且估计钙/水线积分以替换简单的“仅水”假定。所估计的钙/水线积分可以被认为是谱估计投影数据，其可以由重建单元与谱测量投影数据一起用于生成最终图像。
例如，非光谱中间图像可以被分割，以便将图像分裂为骨和非骨(尤其是，软组织)分量。对于每个非谱探测元件而言，针对骨分量的相关材料厚度L骨和针对非骨分量的LST可以通过使用描述通过不同的分割分量的辐射路径(即，描述相应的分割分量的前向投影)的投影数据生成参数来确定。估计单元然后可以根据以下等式来估计谱投影数据(即，模拟谱投影数据的测量结果)。
其中，在等式(5)和(6)中，μBone指代线性骨吸收系数并且μST指代线性非骨吸收系数。由于骨的吸收与钙的吸收基本上类似并且组织的吸收与水的吸收基本上类似，因此针对μBone可以使用钙线性吸收系数并且针对μST可以使用水线性吸收系数。
尽管在上文所描述的实施例中，成像系统是计算机断层摄影系统，但是在其他实施例中，成像系统还可以是如X射线C型臂系统的另一成像系统。
尽管在上文所描述的实施例中，成像系统是包括用于采集投影数据的模块和用于处理投影数据以便生成对象的图像的模块的完整的成像系统，但是在其他实施例中，成像系统还可以是不适于采集投影数据(即，不包括投影数据采集系统)的系统。在这种情况下，成像系统包括用于提供测量投影数据的测量投影数据提供单元，其中，测量投影数据提供单元可以仅是用于存储测量投影数据和用于检索量投影数据以用于提供其的存储单元，或者测量投影数据提供单元可以仅是用于接收测量投影数据和用于提供接收的测量投影数据的接收单元。
尽管在上文所描述的实施例中，估计单元适于估计被用于估计谱投影数据值的水分布(尤其是等价水长度)，但是在其他实施例中，估计单元可以适于估计可以被用于估计谱投影数据值的其他材料分布(尤其是其他等价材料长度)。
本领域技术人员通过研究附图、公开内容以及权利要求书，在实践请求保护的本发明时能够理解并且实现对所公开的实施例的其他变型。
在权利要求书中，词语“包括”不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。
单个单元或其他单元可以履行权利要求书中所记载的若干个项目的功能。尽管在互不相同的从属权利要求中记载了特定措施，但是这并不指示不能有利地使用这些措施的组合。
可以通过任何其他数目的单元或设备执行由一个或若干单元或设备执行的操作，如测量投影数据的提供、谱估计投影数据的估计、投影数据的平滑、基于投影数据的图像的重建等。例如，可以通过单个单元或通过任何其他数目的不同的单元执行步骤102和步骤103。根据成像方法的成像系统的这些操作和和/或控制可以被实施为计算机程序的程序代码模块和/或专用硬件。
计算机程序可以存储/分布在适合的介质上，例如与其他硬件一起被提供或作为其他硬件的部分被提供的光学存储介质或固态介质，但是计算机程序也可以以其他形式分布，例如经由因特网或其他的有线或无线的电信系统分布。
权利要求中的任何附图标记都不应被解释为对范围的限制。
本发明涉及一种如用于生成对象的图像的计算机断层摄影系统的成像系统。谱测量投影数据和非谱测量投影数据是通过具有谱探测元件和非谱探测元件的探测器来生成的，并且谱估计投影数据是通过使用能够已经引起非谱测量投影数据的模型材料分布并且通过基于模型材料分布模拟谱估计投影数据的测量结果来估计的。图像是基于测量和估计的谱投影数据来重建的。除谱测量投影数据之外，使用谱估计投影数据可以实现高质量谱成像，尤其是高质量谱计算机断层摄影成像，其使用不仅具有通常更复杂的谱探测元件而且具有更简单的非谱探测元件的简化探测器。
Claims (10)
1.一种用于生成对象的图像的成像系统，所述成像系统(17)包括：
-测量投影数据提供单元(18)，其用于提供谱测量投影数据和非谱测量投影数据，其中，基于已经穿过包括所述对象(14)的成像区域(5)的辐射(4)由具有谱探测元件(16)和非谱探测元件(15)的探测器(6)分别已经生成了所述谱测量投影数据和所述非谱测量投影数据，其特征在于，
-估计单元(10)，其用于估计对应于结束于所述非谱探测元件(15)处的辐射路径的谱估计投影数据，其中，为了估计针对结束于非谱探测元件(15)处的辐射路径的谱估计投影数据的值，所述估计单元(10)适于：a)基于所述非谱测量投影数据来确定所述成像区域(5)内的模型材料分布，所述模型材料分布应当能够引起对应于所述辐射路径的非谱测量投影数据值；并且b)基于所述模型材料分布来模拟所述谱估计投影数据的值的测量结果，
-重建单元(11)，其用于基于所述谱测量投影数据和所述谱估计投影数据来重建所述对象(14)的图像。
2.根据权利要求1所述的成像系统，其中，为了确定用于估计谱估计投影数据的值的所述模型材料分布，所述估计单元(10)适于基于对应的非谱测量投影数据的值来确定沿着对应于要被确定的所述谱估计投影数据的值的辐射路径的模型材料分布。
3.根据权利要求2所述的成像系统，其中，所述估计单元(10)适于估计水分布作为沿着所述辐射路径的所述模型材料分布。
4.根据权利要求1所述的成像系统，其中，为了确定用于估计谱估计投影数据的值的所述模型材料分布，所述估计单元(10)适于基于所述谱测量投影数据和非谱测量投影数据来重建中间图像，并且适于基于经重建的中间图像来确定所述模型材料分布。
5.根据权利要求4所述的成像系统，其中，所述估计单元(10)适于将所述谱测量投影数据进行积分，从而生成积分投影数据，并且适于基于所述积分投影数据和所述非谱测量投影数据来重建所述中间图像。
6.根据权利要求4所述的成像系统，其中，所述估计单元(10)适于通过分割经重建的中间图像中的不同材料来确定所述模型材料分布。
7.根据权利要求1所述的成像系统，其中，所述估计单元(10)适于使所述谱测量投影数据和所述谱估计投影数据在所述谱测量投影数据与所述谱估计投影数据之间的边界处平滑。
8.根据权利要求1所述的成像系统，其中，所述重建单元(11)适于通过使用滤波反投影算法来重建所述对象(14)的所述图像。
9.一种用于生成对象的图像的成像方法，所述成像方法包括：
-由测量投影数据提供单元(18)提供谱测量投影数据和非谱测量投影数据，其中，基于已经穿过包括所述对象(14)的成像区域(5)的辐射(4)由具有谱探测元件(16)和非谱探测元件(15)的探测器(6)分别已经生成了所述谱测量投影数据和所述非谱测量投影数据，其特征在于，
-由估计单元(10)估计对应于结束于所述非谱探测元件(15)处的辐射路径的谱估计投影数据，其中，为了估计针对结束于非谱探测元件(15)处的辐射路径的谱估计投影数据的值，所述估计单元(10)：a)基于所述非谱测量投影数据来确定所述成像区域(5)内的模型材料分布，所述模型材料分布应当能够引起对应于所述辐射路径的非谱测量投影数据值；并且b)基于所述模型材料分布来模拟所述谱估计投影数据的值的测量结果，
-由重建单元(11)基于所述谱测量投影数据和所述谱估计投影数据来重建所述对象(14)的图像。
10.一种存储用于对感兴趣区域进行成像的计算机程序的存储介质，其特征在于，所述计算机程序包括程序代码模块，所述程序代码模块用于当所述计算机程序在控制根据权利要求1所述的成像系统的计算机上运行时令所述成像系统执行根据权利要求9所述的成像方法的步骤。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP14171093.9 | 2014-06-04 | ||
| EP14171093 | 2014-06-04 | ||
| PCT/EP2015/060679 WO2015185343A1 (en) | 2014-06-04 | 2015-05-13 | Imaging system for generating an image of an object |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106456094A CN106456094A (zh) | 2017-02-22 |
| CN106456094B true CN106456094B (zh) | 2020-03-24 |
Family
ID=50980124
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201580029374.9A Active CN106456094B (zh) | 2014-06-04 | 2015-05-13 | 用于生成对象的图像的成像系统 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US10143434B2 (zh) |
| EP (1) | EP3152596B1 (zh) |
| JP (1) | JP6546201B2 (zh) |
| CN (1) | CN106456094B (zh) |
| WO (1) | WO2015185343A1 (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6479967B2 (ja) * | 2014-09-19 | 2019-03-06 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | スペクトル投影拡張 |
| WO2019059632A1 (ko) * | 2017-09-25 | 2019-03-28 | 한국과학기술원 | 프리즘을 이용한 초분광 영상 재구성 방법 및 시스템 |
| WO2019141527A1 (en) * | 2018-01-16 | 2019-07-25 | Koninklijke Philips N.V. | Spectral imaging with a non-spectral imaging system |
| CN108903961A (zh) * | 2018-07-19 | 2018-11-30 | 深圳市倍康美医疗电子商务有限公司 | 一种cbct成像方法、存储介质及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1284659A (zh) * | 1999-08-17 | 2001-02-21 | 高级光学技术股份有限公司 | 用于双能量x射线成像的设备和方法 |
| CN1575760A (zh) * | 2003-07-03 | 2005-02-09 | Ge医药系统环球科技公司 | 用于截断补偿的方法和设备 |
| CN1985765A (zh) * | 2005-12-22 | 2007-06-27 | 通用电气公司 | 利用混合的计算机断层摄影检测器来执行图像重构的方法 |
| CN101166997A (zh) * | 2005-04-26 | 2008-04-23 | 皇家飞利浦电子股份有限公司 | 光谱ct的检测器阵列 |
| CN101529275A (zh) * | 2006-10-27 | 2009-09-09 | 皇家飞利浦电子股份有限公司 | 用于对对象进行成像的成像系统 |
| JP2010172590A (ja) * | 2009-01-30 | 2010-08-12 | Ge Medical Systems Global Technology Co Llc | Ｘ線ｃｔ装置および画像処理装置 |
| JP2012215235A (ja) * | 2011-03-31 | 2012-11-08 | Piolax Inc | ボルトクリップ |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5755031B2 (ja) * | 2011-05-30 | 2015-07-29 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成装置およびｘ線ｃｔ装置並びにプログラム |
| US9970890B2 (en) * | 2011-10-20 | 2018-05-15 | Varex Imaging Corporation | Method and apparatus pertaining to non-invasive identification of materials |
-
2015
- 2015-05-13 WO PCT/EP2015/060679 patent/WO2015185343A1/en active Application Filing
- 2015-05-13 US US15/312,372 patent/US10143434B2/en active Active
- 2015-05-13 CN CN201580029374.9A patent/CN106456094B/zh active Active
- 2015-05-13 JP JP2016571088A patent/JP6546201B2/ja active Active
- 2015-05-13 EP EP15722214.2A patent/EP3152596B1/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1284659A (zh) * | 1999-08-17 | 2001-02-21 | 高级光学技术股份有限公司 | 用于双能量x射线成像的设备和方法 |
| CN1575760A (zh) * | 2003-07-03 | 2005-02-09 | Ge医药系统环球科技公司 | 用于截断补偿的方法和设备 |
| CN101166997A (zh) * | 2005-04-26 | 2008-04-23 | 皇家飞利浦电子股份有限公司 | 光谱ct的检测器阵列 |
| CN1985765A (zh) * | 2005-12-22 | 2007-06-27 | 通用电气公司 | 利用混合的计算机断层摄影检测器来执行图像重构的方法 |
| CN101529275A (zh) * | 2006-10-27 | 2009-09-09 | 皇家飞利浦电子股份有限公司 | 用于对对象进行成像的成像系统 |
| JP2010172590A (ja) * | 2009-01-30 | 2010-08-12 | Ge Medical Systems Global Technology Co Llc | Ｘ線ｃｔ装置および画像処理装置 |
| JP2012215235A (ja) * | 2011-03-31 | 2012-11-08 | Piolax Inc | ボルトクリップ |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2017525402A (ja) | 2017-09-07 |
| CN106456094A (zh) | 2017-02-22 |
| JP6546201B2 (ja) | 2019-07-17 |
| EP3152596B1 (en) | 2019-07-31 |
| US20170079605A1 (en) | 2017-03-23 |
| US10143434B2 (en) | 2018-12-04 |
| WO2015185343A1 (en) | 2015-12-10 |
| EP3152596A1 (en) | 2017-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10034646B2 (en) | Material decomposition of multi-spectral X-ray projections using neural networks | |
| Sellerer et al. | Dual-energy CT: a phantom comparison of different platforms for abdominal imaging | |
| US9159462B2 (en) | Detection apparatus | |
| EP2091437B1 (en) | Ct imaging system | |
| US10542956B2 (en) | Estimation of an attenuation map based on scattered coincidences in a PET system | |
| EP3215015B1 (en) | Computed tomography system | |
| RU2444764C2 (ru) | Система визуализации для визуализации объекта | |
| CN106880902B (zh) | 基于单能量图像记录确定材料属性值的空间分布 | |
| CN105517492A (zh) | 用于处理心脏数据的处理设备 | |
| CN106456094B (zh) | 用于生成对象的图像的成像系统 | |
| CN103717137B (zh) | 图像生成装置 | |
| CN105793894A (zh) | 根据图像数据来进行骨骼分割 | |
| US9177397B2 (en) | Imaging apparatus | |
| JP6021347B2 (ja) | 医用画像撮像装置及び医用画像撮像方法 | |
| US7924968B2 (en) | Imaging system for imaging a region of interest from energy-dependent projection data | |
| KR102325343B1 (ko) | 의료 영상 장치 및 의료 영상 처리 방법 | |
| JP2015532140A (ja) | スペクトル画像診断のための方法及びシステム | |
| CN112969411A (zh) | 用于提供谱图像的系统 | |
| JP5920770B2 (ja) | Ｘ線コンピュータ断層撮影方法及びｘ線コンピュータ断層撮影システム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |