CN107228865B - 计算机断层摄影x-射线显微镜系统中谱表征的方法和系统 - Google Patents

计算机断层摄影x-射线显微镜系统中谱表征的方法和系统 Download PDF

Info

Publication number
CN107228865B
CN107228865B CN201710178999.2A CN201710178999A CN107228865B CN 107228865 B CN107228865 B CN 107228865B CN 201710178999 A CN201710178999 A CN 201710178999A CN 107228865 B CN107228865 B CN 107228865B
Authority
CN
China
Prior art keywords
ray
spectrum
filter
transmission
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
Application number
CN201710178999.2A
Other languages
English (en)
Other versions
CN107228865A (zh
Inventor
黄志峰
托马斯·A·凯斯
洛伦斯·B·斯蒂格
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Carl Zeiss X Ray Microscopy Inc
Original Assignee
Carl Zeiss X Ray Microscopy Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Carl Zeiss X Ray Microscopy Inc filed Critical Carl Zeiss X Ray Microscopy Inc
Publication of CN107228865A publication Critical patent/CN107228865A/zh
Application granted granted Critical
Publication of CN107228865B publication Critical patent/CN107228865B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/06—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 measuring the absorption
    • G01N23/083—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 measuring the absorption the radiation being X-rays
    • G01N23/087—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 measuring the absorption the radiation being X-rays using polyenergetic X-rays
    • 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
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/30—Accessories, mechanical or electrical features
    • G01N2223/303—Accessories, mechanical or electrical features calibrating, standardising
    • G01N2223/3032—Accessories, mechanical or electrical features calibrating, standardising periodic calibration, e.g. with filter wheel
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/30—Accessories, mechanical or electrical features
    • G01N2223/309—Accessories, mechanical or electrical features support of sample holder
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/30—Accessories, mechanical or electrical features
    • G01N2223/313—Accessories, mechanical or electrical features filters, rotating filter disc
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/401—Imaging image processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/418—Imaging electron microscope
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/419—Imaging computed tomograph

Landscapes

  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Pulmonology (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Toxicology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

一种用于断层重建、束硬化校正、双能量CT和系统诊断等等的谱测量和估计方法，包括确定源加速电压、预滤器和/或检测器的组合的谱，并且在测量多个预滤器的透射值以后计算所述源加速电压、预滤器和/或检测器组合的校正谱。

Description

计算机断层摄影X-射线显微镜系统中谱表征的方法和系统
背景技术
X-射线计算机断层摄影(CT)是一种用于检查和分析样品的内部结构的非破坏性技术。一般来说，当X-射线穿过样品时，X-射线被样品吸收或散射。没有被吸收或散射离开的X-射线透射通过样品，并且然后被检测器系统检测。在检测器系统形成的图像被称为X-射线投影。经由标准的CT重建算法，比如滤过反向投影算法(FBP)，从这一系列的不同角度的投影得出层析成像体积数据集。
某些X-射线CT系统使用多色X-射线源来生成所述X-射线投影。多色X-射线源包括X-射线管(实验室光源)、白色同步加速器束或基于加速器的光源。总体而言，从这些光源发射的多色X-射线束包括具有许多不同能量的X-射线。这些X-射线的根据能量的分布规律一般叫做该射束的谱。多色X-射线束与只包括一种能量或很窄的能量范围的X-射线的单色射束是不一样的。
白色同步加速器束即为一种多色X-射线源。随着电子通过所述同步加速器加速，所述电子以很窄的与环相切的角度释放出强烈的“白色”的多色辐射。该射束包括从柔和的紫外线到强烈的X-射线的辐射。可以通过应用滤波器，比如传递/反射选定能量的辐射的晶体单色器和/或X-射线镜，就可以从所述白色多色射束中创建单色X-射线源。
实验室光源，比如X-射线管也可以产生多色X-射线，一般使用固定或旋转阳极。在真空管中，灯丝(比如钨丝)作为阴极而金属靶作为阳极。在所述阴极上施加高电压(加速)，在所述阴极和阳极间创建高电势。这使得电子在真空中从阴极向阳极流动并加速。电子与所述阳极材料碰撞并加速所述阳极材料中的其他电子、离子以及核子。这一过程生成X-射线。X-射线管产生的X-射线的谱是依据所述阳极靶材料和所述加速电压的结果。所述谱的特征在于连续地制动辐射谱(“刹车辐射”)X-射线，以及因为阳极材料芯电子去激发和离子化产生的特定频率的次级X-射线辐射，也叫X-射线荧光(XRF)。所述次级X-射线辐射或X-射线线为阳极材料中所用的金属所特有。所述制动辐射是主要的X-射线，可以用作多色X-射线源，而X-射线滤波器和/或X-射线镜可以应用在所述特征X-射线辐射和/或所述制动辐射以产生单色X-射线，如实施例所述。
X-射线源发射的X-射线谱在使用所述X-射线源的期间会随着时间改变。例如，在X-射线管中，由于所述电子撞击所述阳极靶材料，靶材料的某些部分会烧蚀或剥落。这造成靶子随着时间的过去发生质量损失(例如变得单薄)。这一现象也叫做“光源靶烧毁”。随着所述靶子变得单薄，所述靶子传递的低X-射线能量的X-射线越来越多，由此造成所述X-射线的X-射线谱随着时间过去发生改变。
在X-射线CT系统中使用多色X-射线束是有优势的。使用多色X-射线光的主要优势在于，对于给定的光源，多色X-射线束一般比单色X-射线更强。这是因为避免了有损耗的能量滤波器。
但是使用多色X-射线束也有劣势。与单色光不同，多色X-射线的所述X-射线吸收通常与样品材料的厚度不成比例。这是因为当射束穿过样品的时候，所述多色X-射线束的较低X-射线能量相对较低X-射线能量被样品吸收更多。因此，在使用多色X-射线束生成X-射线投影的时候，使用一种叫做束硬化(BH)的处理。束硬化与X-射线谱中随着X-射线穿过样品朝着更高的X-射线能量的改变有关。
在从多色X-射线到断层摄影重建的期间，束硬化常常产生伪影。束硬化生成的典型伪影包括杯突伪影以及条纹伪影。一般而言，样品中具有较高原子序数(Z)的元素，比如金属，在断层摄影重建影像中比原子序数低的元素会产生更多的BH伪影。
为了减少或防止多色X-射线束创建的断层摄影重建中的伪影，对系统的X-射线源发射的X-射线的精确能量谱有先验知识，以及对不同X-射线能量上的检测系统的一个或多个检测器的灵敏度有先验知识是非常重要的。因此，在使用多色X-射线源的X-射线CT系统创建的样品的断层摄影重建的伪影减少中，X-射线源谱测量或估计通常是关键问题。
对于X-射线成像/CT系统中X-射线源给定的加速电压(kVp)、预滤器和检测器的组合，典型的谱测量和估计方法如下：1)使用已知的标准体模(比如步进式光楔体模和锥体模)获取X-射线透射测量结果；以及2)通过某些迭代算法(例如期望最大化<EM>算法)，基于所述透射测量结果估计该组合的实际谱。一般而言，用蒙特卡洛模拟谱作为初始输入。参见L.Zhang,G.Zhang,Z.Chen,Y.Xing,J.Cheng and Y.Xiao,"X-ray spectrum estimationfrom transmission measurements using the expectation maximization method,"2007IEEE Nuclear Science Symposium Conference Record,M13-293,3089-3993；andE.Sidky,L.Yu,X.Pan,Y.Zou and M.Vannier,"A robust method of x-ray sourcespectrum estimation from transmission measurements:Demonstrated on computersimulated,scatter-free transmission data,"Journal of Applied Physics 97,124701(2005)。这两篇文章都在此引用。
发明内容
当前X-射线CT系统中的X-射线源谱测量和估计有不足之处。因为当前大多数灵活的X-射线CT系统支持多个加速电压、预滤器和/或检测器，所述透射率测量程序必须针对所述加速电压、预滤器和/或检测器的每个组合重复。因此，对所述不同组合的X-射线源谱测量和估计是非常消耗时间的。
而且，基于体模的X-射线源谱测量/校准不适于长期断层摄影。因为X-射线源产生的所述X-射线的能量谱会随着时间改变(甚至在扫描期间就改变)，长期扫描一般要求至少有一次使用标准体模测量X-射线源发射的X-射线的能量谱的改变的扫描中再校准。为了这一目的，操作人员一般会停止扫描，从X-射线CT系统中移除样品(一般也包括样品座)，并将一个标准体模置于光路中以执行X-射线源谱测量。完成所述一次或多次校准之后，操作人员移除所述体模，更换包含样品的样品座并恢复扫描。这种X-射线源谱校准不方便、耗时久，并且可能在样品更换并恢复扫描之后生成的断层摄影重建中造成误差。如果这些样品没有被精准地如谱再校准之前一样地放置在光路中，则可能产生这些误差。
X-射线源发射的X-射线的谱是X-射线CT系统的主要特征。因而，X-射线源谱校准以及测量可以被用作这些系统的问题或维护方面的诊断。
本发明是简化的X-射线谱测量估计的系统和方法。本发明利用一些预滤器，所述预滤器作为X-射线源谱测量和/或校准工具具有不同的材料和/或不同厚度，而非使用传统的步进式光楔体模或锥体模。在一个实施例中，所述X-射线源谱的测量和校准使用X-射线CT系统的标准滤光轮组件。所述滤光轮一般包括一打或以上的位于X-射线源输出点附近的预滤器。所述预滤器用于在X-射线束入射到样品之前将X-射线束中的具有不同能量的X-射线移除。本发明的实施例的另一方面是关于在X-射线CT系统的计算机系统中执行的软件程序。
因为X-射线CT系统不在真空中操作，而是使用空气作为X-射线束的透射媒介，在使用所述滤光轮的具有不同厚度/密度的预滤器的时候，所述X-射线CT系统必须考虑由空气分子造成的少量的X-射线吸收。为了这一目的，优选地，所述预滤器中的一个是所述滤光轮中的一个空斑，也叫做空气预滤器。根据本发明的原理，先获取X-射线源透射测量。出于这个目的，所述X-射线CT系统中没有放置样品或者所述样品被移出所述X-射线源的光路，并且获取对滤波器、X-射线能量以及至少一个检测器的不同组合的透射测量结果。一般先做所述空气预滤器的所述测量，然后进行对不同材料和/或非零厚度的预滤器的透射测量。每个预滤器的透射曲线是根据所述透射测量绘制，每个预滤器的衰减曲线是根据所述相关的透射曲线绘制的，并且随后根据所述衰减曲线创建一个基线谱集。所述基线谱包括谱，创建所述谱是用于每个预滤器的衰减曲线并与所述每个预滤器有关，包括所述空气预滤器。
然后，在所述X-射线CT系统的操作期间(例如在扫描样品的期间)，每个预滤器的透射值，例如空气预滤器的透射值，会被测量并跟同一预滤器在所述基线谱集中的相关透射值进行比较。这一比较被用来决定所述X-射线源的X-射线谱是否已经随着时间发生改变，即自计算了所述基线谱以来是否已经随着时间发生改变。如果所述比较显示所述X-射线谱已经发生改变，本发明可以用来通过估计所述X-射线源发射的X-射线的谱发生的改变来进行校准。
为了这一目的，可以通过仅对一些(例如2或3个)选定的预滤器进行新的透射测量来估计所述X-射线源谱发生的改变。为了这个估计，所述新的透射测量结果随即被“拟合”相应的透射曲线(例如衰减曲线)，从中创建所述基线谱。然后这些经计算的谱被用在所述样品的投影所决定的所述断层摄影重建的计算中。在别的应用中，这些经计算的谱被用来诊断所述X-射线源的问题，或者评估X-射线源的目标(例如实验室X-射线管的阳极)的剩余寿命。
采用本发明的原理可以取得数个优势。首先，所述新的谱测量和校准方法不需要标准的比如步进式光楔体模和锥体模，而所述步进式光楔体模和锥体模对于现有的谱测量方法是必须的。其次，不再需要像现有的谱测量方法那样去暴力测量X-射线源加速电压、预滤器和检测器的每个不同组合的衰减曲线。取而代之的是源加速电压和预滤器和检测器的不同组合的所有的衰减曲线可以基于更少的、现有测量方法对至少一个源加速电压和预滤器和至少一个检测器不同组合的测量来计算。操作者可以将数个或者全部的加速电压、预滤器和检测器的组合的谱存档，但是一次只使用一个用于断层摄影。第三，因为不再需要用体模置换出样品(或者说样品和样品座)然后更换所述样品并恢复扫描，因为所述样品的不完全的对准给所述断层摄影重建造成误差的可能性被大大降低。
本发明提出的谱测量和估计方法可以用来排除和减少所述标准步进式光楔体模或锥体模的使用。取而代之地，现有X-射线CT系统的滤光轮可以被用来完成同样的目标，而可能仅需要针对必要的预滤器组选择和/或改变所述滤波器。这简化了谱测量。
本发明提出的所述系统和方法要求进行的谱测量较之现有系统的方法也更少。相反，使用CT系统X-射线的预滤器、加速电压和检测器的不同选择，可以基于少至2次或3次扫描来估计加速电压和预滤器和检测器的多个不同组合的谱。这极大地节约了时间。
最后，所述提出的系统和方法排除了现有X-射线源谱校准方法因为要向X-射线CT系统中装入体模并更换样品回X-射线CT系统而引起的潜在计算误差。例如，当前方法通常要求从系统中移除所述样品座以及样品并将所述体模置于所述系统以执行所述X-射线源校准，以及随后再将所述样品座和样品放回所述系统以完成对该样品的扫描。这一过程是添乱并耗时的。相反，有了本发明提出的方法，操作人员一般仅需将所述样品转出所述X-射线CT系统的光路以执行所述校准步骤，然后再以之前存储的样品对准设置将所述样品转回到所述光路中以执行所述扫描。所述样品无需从所述X-射线CT系统中移除。
总体而言，根据一个方面，本发明的特征是X-射线CT系统中的X-射线谱测量和估计方法。本方法确定至少一个X-射线加速电压、滤波器和至少一个检测器的多个不同组合的从X-射线CT系统的X-射线源发射的X-射线的基线谱。在所述X-射线CT系统操作期间，监测所述X-射线以确定所述X-射线源谱的改变，并且响应于确定所述X-射线源谱已经发生改变而基于一个或多个滤波器，比如至少2或3个滤波器的测量透射值来计算新的基线谱。
优选地，所述方法用于断层摄影重建和束硬化校正。该方法还能用于多能量计算机断层摄影以及用于诊断所述X-射线CT系统。
在一个实施例中，所述方法监测所述X-射线以通过使用X-射线加速电压、一个滤波器和检测器的组合来获取透射值，再将所获取的透射值与以所述X-射线加速电压、所述一个滤波器以及所述检测器的相同组合的基线谱内的透射值进行比较。
一般而言，所述方法确定所述基线谱是通过使用X-射线加速电压、滤波器和检测器的组合来获取透射测量结果，并从透射曲线计算每个所述滤波器的衰减曲线。
所述方法确定所述基线谱此外或另外还可以通过基于X-射线加速电压、滤波器和检测器的组合获取的透射测量结果计算每个所述滤波器的衰减曲线，再将期望最大化(EM)算法用于所述衰减曲线。
所述方法计算新的基线谱是通过识别用于测量至少2或3个滤波器的透射值的加速电压以及检测器的组合，以所识别的加速电压和检测器组合来来测量空气滤波器的透射值，再使用所述空气滤波器的所述透射值以及所述至少2或3个滤波器的所述透射值来创建透射曲线并从所述透射曲线中创建衰减曲线。
所述方法计算新的基线谱是此外或另外通过从所述基线谱以通常的加速电压和检测器组合来选择空气滤波器谱，针对所选择的空气滤波器谱来拟合所述衰减曲线，再使用EM算法，在所拟合的衰减曲线内以所述通常的加速电压和检测器组合根据所述空气滤波器谱来估计所述空气滤波器的校正谱。
所述方法计算新的基线谱是此外或另外通过使用所述测量的至少2或3个滤波器的透射值来创建空气滤波器的估计的校正谱，再从所述空气滤波器的所述估计的校正谱来计算所述新的基线谱。
所述方法优选在没有幻影体模存在的情况下完成至少基于测量2或3个滤波器的测量透射值来计算新的基线谱模。
所述方法此外还可以在扫描所述样品的期间使用所述基线谱以及所述新的基线谱以校正所述样品的断层摄影重建中的束硬化。
一般而言，根据另一方面，本发明的特征是X-射线显微镜系统中的X-射线谱测量和估计方法，针对诸如以下的应用：断层摄影重建、束硬化校正、多(双)能量计算机断层摄影(CT)以及X-射线显微镜系统诊断。所述方法确定X-射线源加速电压、预滤器和/或检测器的组合的谱。在操作所述X-射线显微镜之前，通过测量多个预滤器的透射值而从所测量的透射值，估计或计算所述源加速电压、预滤器和/或检测器的组合的谱。
一般而言，根据另一方面，本发明的特征是X-射线CT系统，所述X-射线CT系统包括X-射线成像系统以及计算机系统，所述计算机系统具有控制器，配置用来控制所述X-射线CT系统的组件。所述X-射线成像系统包括X-射线源、滤波器以及检测器系统。所述计算机系统确定以X-射线源的X射线加速电压、滤波器以及所述检测器系统的检测器的不同组合的从所述X-射线源发射的X-射线的基线谱；监测所述X-射线以确定所述X-射线源谱的改变；并且响应于确定所述X-射线源谱已经改变而基于测量一个或多个滤波器，例如至少2或3个滤波器的透射值来计算新的基线谱。
优选地，所述新的基线谱用于断层摄影重建和束硬化校正。所述新的基线谱还可以用于多能量计算机断层摄影和所述X-射线CT系统诊断。
在一个实施例中，所述计算机系统确定X-射线源谱的改变是通过使用X-射线加速电压、一个滤波器以及检测器的组合来获取透射值，再将所获取的透射值与以所述X-射线加速电压、所述一个滤波器以及所述检测器的组合的在所述基线谱内的透射值进行比较。
所述计算机系统确定所述基线谱是通过使用X-射线加速电压、滤波器和检测器的组合来获取透射测量结果，使用所述透射测量结果计算每个所述滤波器的透射曲线，再从所述透射曲线计算每个所述滤波器的衰减曲线。
所述计算机系统确定所述基线谱此外或另外还可以通过基于使用X-射线加速电压、滤波器和检测器的组合获取的透射测量结果来计算每个所述滤波器的衰减曲线，再将期望最大化(EM)算法用于所述衰减曲线。
所述计算机系统计算新的基线谱是此外或另外通过识别用于测量所述至少2或3个滤波器的透射值的加速电压和检测器的组合，以所识别的加速电压和检测器组合来测量空气滤波器的透射值，再使用所述空气滤波器的所述透射值以及所述2或3个滤波器的透射值创建透射曲线并从所述透射曲线来创建衰减曲线。
参考附图和权利要求，以下将更具体地说明本发明的上述和其他特征，包括部件的结构和组合的各种新颖细节以及其他优势。应当理解体现本发明的具体方法和设备只是用作说明而非限制本发明。在不脱离本发明的范围的情况下，本发明的原理和特征可以被用于各种和无数的实施例中。
附图说明
在附图中，附图标记通过不同的视图指示相同的部分。附图不一定按照比例绘制；重点在于示出本发明的原理。图中：
图1是本发明的方法可适用的X-射线CT系统的示意图；
图2是具有滤波器的X-射线CT系统的示例性滤波轮(filter wheel)的示意图；
图3A是根据本发明的原理，用于获取从X-射线CT系统的X-射线源发射的X-射线的基线谱集的方法的流程图；
图3B-1是用于X-射线CT系统的谱校准方法的流程图，其中在执行样品的单次短持续时间断层摄影扫描期间，该方法应用来自图3A的方法的基线谱集以生成样品的谱校正投影；
图3B-2是用于样品的长持续时间扫描的谱校准方法的流程图；
图4是示出提供用于从图3A、3B-1和3B-2的方法中的所选择的加速电压、预滤器和检测器的组合计算谱的更多细节的方法的流程图；
图5示出包括作为用于AIR(即，空气预滤器)测量的、HE1预滤器测量的和HE1预滤器计算的滤波器指数的函数的衰减值的不同衰减曲线的图，；
图6是作为X-射线加速能量(kVp)的函数的两个不同的估计的校正X-射线能谱的图；
图7是从使用在70kVp，80kVp，90kVp，…，150kVp(这里仅示出70kVp，90kVp，110kVp，130kVp和150kVp)下的滤波器集的实际衰减曲线测量结果，用于更精细加速电压的衰减曲线(即，衰减曲线是透射曲线的负对数)的内插的示例；以及
图8示出了根据本发明的原理创建的X-射线CT系统的X-射线源谱的衰减曲线的曲线图，其中已经确定X-射线源谱随时间改变，并且其中可以利用仅使用两个预滤器透射测量结果(LE1和HE5)创建的拟合的衰减曲线的图以估计跨所有预滤器设置、加速电压和检测器组合的改变的X-射线能谱的校正。
具体实施方式
现在将参照附图更全面地描述本发明，附图中展示了本发明的说明性实施例。然而，本发明可能以许多不同形式实施，并且不应被解释为限于本文所阐述的实施例；相反，提供这些实施例使得本公开将是详细和完全的，并且将向本领域的技术人员充分地传达本发明的范围。
如本文所使用的，术语“和/或”包括一个或多个相关所列项目的任何和所有组合。此外，除非另有明确说明，本文中冠词“一”，“一个”和“该”的单数形式也意为包括复数形式。还应当理解的是，当在本说明书中使用术语“包括”，“包含”，和/或包括，或指定所述特征的存在，整数，步骤，操作，元件和/或组件的存在，但不排除存在或添加一个或多个其它特征，整数，步骤，操作，元件，和/或组件。进一步地，应当理解的是，当包括组件或子系统的元件被称为和/或示出为连接或耦合到另一元件时，它可以直接连接或耦合到另一元件，或者可以存在中间元件。
对于单色X-射线，比尔定律描述了当X-射线通过纯净物体时X-射线将被衰减：
I＝I0 e-μt
其中I0是入射在物体上的X-射线强度，t是物体厚度，I是透射物体的X-射线强度，μ是物体的衰减系数(包括光电吸收和康普顿散射的效应，以及若X-射线能量>1.02MeV时电子对效应)，其取决于样品密度ρ，原子序数Z，原子质量A和X-射线能量E。
上述公式可以用积分形式改写：
I＝I0e-∫μ(l)dl，和
∫μ(l)dl＝-ln(I/I0)
其中l是穿过物体的束路径。这完全满足了CT重建算法的线积分需求。其结果是，对于单色X-射线不引入束硬化(BH)伪影。
另一方面，多色X-射线源在谱D(E)上生成X-射线。为简单起见，D(E)也包括检测器灵敏度的影响。物体的下游强度I由下式给出：
I＝I0∫D(E)e-∫μ(E,I)dldE
μ(E)通常是E的非线性函数，由下式给出：
Figure GDA0003662177930000101
其中
Figure GDA0003662177930000102
n≈4(光电吸收部分)，以及
Figure GDA0003662177930000103
(康普顿散射部分)，其中fKN是克莱-仁科公式。
由于μ(E)通常是能量E的非线性函数，来自多色X-射线源的X-射线的透射不能满足CT重建算法的线积分需求。μhigh-energy<μlow-energy，意味着该材料相对于高-能量X-射线吸收“更多的”低-能量X-射线。通过样品透射的衰减的X-射线的输出谱看起来难以“移动”，被称为束硬化(BH)。高-Z金属材料具有更严重的束硬化问题，导致在CT重建断层摄影图像中典型的金属伪影。然而，如果从源发射的X-射线的谱是已知的或可以导出，可以在断层摄影重建中减少或去除这些伪影。
图1是本发明的方法可适用的X-射线CT系统100的示意图，并且其通常易受束硬化问题影响。
一般来说，X-射线CT系统100包括生成多色X-射线束104的X-射线源102，以及具有样品保持器112的旋转台110，用于将样品114保持在来自X-射线源系统102的X-射线束104中。图像或X-射线投影通过检测器系统118捕捉。X-射线CT系统102、旋转台110和检测器系统安装到X-射线CT系统100的基座108。计算机系统200通常接收和处理这些图像，并且提供系统100的一般控制。计算机系统200或专用图形处理器通常将使用X-射线投影执行断层摄影重建。
在一个示例中，X-射线源102是多色X-射线源。多色X-射线源优选是实验室X-射线源，因为它的无处不在和相对低的成本。尽管如此，同步加速器源或基于加速器的源是另一种选择。
常见的实验室X-射线源包括X-射线管，其中电子在真空中通过电场加速并射入金属靶片，随着电子在金属中减速而发射X-射线。通常，取决于所使用的金属靶的类型，这样的源产生与在某些能量的强度尖锐峰组合的背景X-射线(即轫致辐射)的连续谱，所述某些能量从靶的特征线得到。另外，X-射线束是分散的，并缺乏空间和时间相干性。
在一个示例中，X-射线源102是具有钨靶的旋转阳极型或微聚焦源。也可以采用包括钼、金、铂、银或铜靶。优选地，使用X-射线源102的透射结构，其中电子束从薄靶103的背侧撞击薄靶103。从靶103的另一侧发射的X-射线用作束104。
在另一个更特定的示例中，源102是结构化的阳极X-射线源，如在2008年10月28日授权给Yun等人的美国专利号7443953中描述的，其内容通过引用整体并入本文。在Yun中，该源具有由期望的靶材制成的薄顶层和由具有最佳热性质的低原子序数和低密度材料制成的厚底层。阳极可以包括，例如，沉积在铍或金刚石基底层上的具有最佳厚度的铜层。
还可以采用产生具有适合于本文所述的断层摄影应用的能量的辐射的X-射线激光器。
在又一个示例中，X-射线源102是金属射流X-射线源，例如可以从瑞典的ExcillumAB，Kista处获得。这种类型的源使用微聚焦X-射线管，其中阳极是液体-金属射流。因此，阳极是连续再生的并已经熔融。
由源102生成的X-射线束104具有通常通过源的操作参数控制的能谱。在实验室源的情况下，支配参数包括靶材和加速电压(kVp)。能谱也由抑制不需要的能量或辐射波长的任何调节滤波器支配。例如，使用例如能量滤波器(设计用于选择期望的X-射线波长范围/带宽)来消除或衰减存在于束中的不期望的波长。
除了X-射线源102之外，本发明依赖于滤波器的可用性，该滤波器在与样品114相互作用之前过滤X-射线束104(预滤器)。
在一个优选地实施例中，当前X-射线CT系统的标准滤波轮150内的滤波器被用作校准工具的预滤器。滤波轮150的滤波器被放置在X-射线束103的路径中，并且被插入X-射线源102和样品114之间。当前X-射线CT系统的操作者使用滤波器以修改和/或补偿对于X-射线输出的束硬化属性。滤波轮150的滤波器用作测量工具的预滤器。
更详细地，滤波轮150由计算机系统200的控制器210控制。滤波轮150包括框架155，其在控制器210的控制下经由其控制界面130在轴154上旋转。经由滤波轮150的控制界面130，操作者可以选择滤波轮150，以使轮150的滤波器中的一个与X-射线源系统102的靶103上的出口孔相邻。以这种方式，所选择的滤波器(在一个示例中，空气预滤器)被对准以在与样品114相互作用之前过滤束104。
图2示出滤波轮150的一个实施例。滤波轮150包括安装在框架155中的16个单独的预滤器。这些预滤器包括六个增加厚度的材料#1滤波器(LE1，LE2，…，LE6)和六个增加厚度的材料#2滤波器(HE1，HE2，…，HE6)。为了测量最高能量，向预滤器集中添加3个增加厚度的材料#3滤波器(BE1，BE2，BE3)和其它预滤器。这里，材料#1、#2和#3具有不同的密度。通常，材料#1具有用于过滤低能量X-射线的轻密度，材料#2具有用于过滤低能量和中能量X-射线的中密度，并且材料#3具有用于过滤全带X-射线的重密度。最后，在滤波轮150中存在一个空的点156(即，空气点或空气预滤器)。
返回到图1，当样品114暴露于X-射线束104时，透射穿过样品的X-射线光子形成由检测器系统118接收的衰减的X-射线束106。在一些其它的示例中，诸如波带板透镜的物镜用于在X-射线成像系统100的检测器系统118上形成图像。在替代的实施例中，检测器系统是平板检测器。
在检测器系统118的最常见结构中，使用几何放大，样品114的放大投影图像以等于源到样品距离与源到检测器距离之间的反比的放大率形成在检测器系统118上。一般来说，由X-射线台提供的几何放大率在2x和100x之间，或者更大。在这种情况下，X-射线图像的分辨率受到X-射线源系统102的焦点尺寸或虚拟尺寸限制。
为了实现高分辨率，X-射线CT系统100的一个实施例进一步利用检测器系统118的非常高分辨率的检测器124-1，结合将样品114定位在靠近X-射线源系统102。在高分辨率的检测器124-1的一个实施方式中，闪烁体与显微镜物镜结合使用以提供在2x和100x之间，或者更大的范围中的附加放大率。
可以包括其它可能的检测器，作为在X-射线CT系统100中所示的检测器系统118的一部分。例如，检测器系统118可以包括一个较低分辨率的检测器124-2，如图1的所示的实施例中所示。在示例中，这可以是平板检测器或具有较低放大率的显微镜物镜的检测器。检测器系统118的一个、两个或甚至多个检测器124的结构是可能的。
优选地，两个或更多个检测器124-1、124-2安装在检测器系统118的转台122上，使得它们可以交替地旋转到来自样品114的衰减的束106的路径中。
通常基于操作者限定的参数，计算机系统200的控制器210经由控制界面130指示旋转台110以在X-射线源系统102校准期间将样品114移出束路径。在完成校准部分之后，控制器210将样品114移回到束路径中并相对于束104旋转样品114，以执行样品114的CT扫描。
在一个示例中，检测器系统118创建来自与检测器系统118的检测器124-1、124-2中的闪烁体交互的衰减的X-射线束106的X-射线光子的像素中的图像表示。在检测器系统118处形成的该图像被称为X-射线投影或X-射线投影图像。
在一个示例中，计算机系统200包括图像处理器220，其分析X-射线投影并且可能执行从X-射线投影创建的断层摄影重建所需的计算。连接到计算机系统200的显示装置240显示来自X-射线CT系统100的信息。诸如触摸屏、键盘和/或计算机鼠标的输入装置250实现在操作者、计算机系统200和显示装置240之间的交互。
计算机系统200从连接到计算机系统200的数据库260来加载信息并将信息存储到数据库260。
在一个示例中，使用在计算机系统200上执行的在显示装置240上显示其界面的用户界面应用，操作者限定/选择CT扫描或校准参数。这些包括X-射线加速电压设置，和用于限定X-射线源系统102上的扫描和曝光时间的X-射线能谱的设置。操作者通常还选择其它设置，例如为样品114创建的X-射线投影图像的数目和用于旋转X-射线束104中的X-射线CT扫描的样品114的旋转台110的角度。
计算机系统200，借助于其图像处理器220，接受来自检测器系统118的与样品114的每个旋转角度相关联的图像或投影信息。图像处理器220为样品114的每个旋转角度创建单独的投影图像，并且使用CT重建算法结合投影图像以为样品创建三维断层摄影体积信息。
为了补偿束硬化，在执行样品114的扫描之前必须知道X-射线源系统102的谱，以减少由X-射线CT系统100用于样品114生成的投影和/或断层摄影重建中的伪影。
图3A示出了用于估计X-射线源系统102的谱的校准方法。然后，这些估计的谱由图像处理器220使用，以生成具有减少的伪影的样品114的断层摄影重建。在该方法期间，样品114旋转或移出束104路径。
在步骤308中，计算机系统200从操作者接收指令，以选择空气滤波器、N个附加滤波器以及X-射线源系统102的不同加速电压和检测器的不同组合/对M，并且将指令传递到控制器210。需要用于一个加速电压和一个检测器的至少一对M。
在步骤310中，控制器210选择X-射线源系统102的检测器124-1/124-2和能量设置(例如加速电压(kVp))的组合/下一对M。然后，在步骤312中，控制器210经由滤波器150的控制器界面130发送信号以旋转滤波轮150，直到AIR预滤器156进入束104的路径中。
在步骤314中，获取用于当前所选择的滤波器(这里是AIR预滤器156)透射值，并且将该透射值储存到用于AIR预滤器的缓冲器。
然后，在步骤316中，滤波轮150旋转到下一个所选择的预滤器，并且该方法转换回到步骤314，以使用在步骤310中所选择的(加速电压、检测器)设置的相同对M来获取用于下一个预滤器的透射值。对于除空气滤波器156之外的N个附加滤波器执行步骤314和316N次。这由附图标记317-N表示。在步骤314的每次通过结束时，然后将在指定(加速电压、检测器)设置下的(当前)预滤器的透射值存储到用于当前预滤器的缓冲器。在迭代处理步骤314/316在(加速电压、检测器)设置的相同对M上重复N次之后，方法转换回到步骤310，以获取(加速电压、检测器)设置的下一对M。这由附图标记319-M表示。然后执行步骤310、312和314/316，并且计算对于空气预滤器156和下一个预滤器的透射值，并将其存储到它们各自的缓冲器。重复步骤310、312和314/316直到M次，对应于(加速电压、检测器)所选择的M对的数量。
在步骤320中，从每个预滤器的透射值的缓冲器创建对于每个预滤器的透射曲线，并且从每个透射曲线计算衰减曲线。每个透射/衰减曲线与不同的预滤器(LE1，LE2，...，LE6，HE1，HE2，...，HE6，BE1，BE2，BE3和AIR 156和可能的其它滤波器)相关联。每个预滤器的衰减曲线是其相关联的透射曲线的负对数。然后，将每个预滤器的透射值、透射曲线和衰减曲线的缓冲器作为吸收参考数据230存储到计算机系统200和/或数据库260。
根据步骤321，然后计算机系统200通过将期望最大化(EM)算法应用于衰减曲线来计算X-射线源谱的估计，也称为基线谱集354。为此目的，在一个示例中，利用具有特定吸收边缘的蒙特卡罗模拟谱作为EM算法的初始输入，使得算法可以更快地收敛。在图4的方法中提供了使用EM算法的基线谱集354的计算的更多细节，其描述包括在下文中。
然后在样品114的后续X-射线扫描期间，可以应用基线谱集354以修复束硬化伪影，否则该束硬化伪影可能出现在从样品114的扫描创建的断层摄影重建中。一般来说，一旦获取基线谱集354，除非更换或更新X-射线管102(或其靶103)或检测器124，否则包括重复步骤310-321的重新校准在每天或每次断层摄影扫描中不是必需的。
注意，所选择的加速电压可以是以10kVp的间隔离散的，例如40kVp，50kVp，60kVp，...，160kVp。如果在实践中需要某些指定的加速电压，例如66kVp或115kVp，则可以通过使用基线谱354内插计算这些加速电压的X-射线的谱。
图3B-1和3B-2示出了用于创建结合X-射线源谱监测的样品114的断层摄影重建的简化的工作流程。从通过扫描样品114生成的X-射线图像投影创建断层摄影重建。图3B-1的方法适用于样品114的单次短持续时间扫描，而图3B-2的方法适用于在较长的扫描持续时间相同样品114的可能的多次扫描。
谱监测方面或功能确定在样品114的扫描期间X-射线源102的谱是否已经充分改变，并且可以使用在图3A的方法中创建的基线谱集354重新校准X-射线谱。图3B-1和图3B-2中的附图标记360表示与谱监测功能相关联的步骤，更简单的称为谱监测器360。谱监测器360包括步骤323和324。
图3B-1和3B-2的方法包括类似的步骤，并且大体上以类似的方式执行。该方法通常由计算机系统200执行和控制，其中任何计算密集型活动可能由图像处理器220或其它协处理系统执行。首先整体描述图3B-1的方法，然后描述图3B-2。
在图3B-1中，根据步骤322，选择指定的加速电压、预滤器和检测器作为样品114的扫描的一部分。然后，在步骤323中获取参考值。对于谱监测器360的设置，测量某个预滤器的透射值，以确定谱是否已经从由基线谱354指定的谱改变。在一个示例中，为了确定谱是否已经改变，将测量的透射值与吸收数据230中所选择的预滤器、加速电压和检测器的相同组合的存储的透射值进行比较。在另一个示例中，为了确定谱是否已经改变，将测量的透射值与在X-射线加速电压、预滤器和检测器的相同组合下的基线谱354内的透射值进行比较。
如果在步骤324中确定X-射线源102的谱没有改变，则在步骤350中将样品114旋转到束路径中，并且然后在步骤352-1中获取投影。
然而，一般来说，使用时，X-射线源102的谱将改变。例如，如果扫描需要几个小时或甚至几天，将经常发生源靶燃烧导致的X-射线谱改变。因此一般来说，在较长持续时间的扫描中或不同扫描之间将需要多参考过程，以及由谱监测器360提供的谱监测。
如果在图3B-1的步骤324中确定谱改变是显著的，例如由于误差或噪声而与预期不同，则需要重新校准X-射线源102的谱，并且该方法转换到步骤328以开始重新校准过程。在当前系统和方法中，必须重复使用幻影的典型谱测量结果以重新校准X-射线源102。相比之下，在X-射线CT系统100的优选实施例中，一旦经由图3A的方法获取基线谱集354，则采用以下谱重新校准方法。
根据步骤328，该方法选择加速电压和检测器的组合。然后通过旋转滤波轮150将AIR(空)滤波器156定位到束104路径中。在步骤328结束之后的各个步骤中使用加速电压和检测器的相同组合或通常的设置。
在步骤330中，使用通常的设置获取AIR滤波器透射测量结果。然后，在步骤332中，选择两个或三个或可能更多的非AIR预滤器用于定位到束104路径中以测量以通常的设置的它们的透射值。在一个示例中，所选择的预滤器之一与测量高透射值相关联，并且所选择的预滤器之一与测量最低透射值(<5％)相关联。
在步骤333中，从在步骤330中获取的空气滤波器156的测量的透射值和从在步骤332中获取的2或3个其它预滤器的测量的透射值创建透射曲线。然后从透射曲线创建衰减曲线。
根据步骤334，然后从基线谱集354选择空气预滤器谱，其中所选择的空气滤波器谱与加速能量和检测器的组合相同，其作为先前步骤328、330和332中的“通常的设置”。所选择的空气预滤器谱由附图标记“354-空气”表示。
然后相对于该所选择的空气预滤器谱，354-空气“拟合”衰减曲线。一般来说，为此目的可使用多项式拟合方法。
在步骤336中，该方法通过使用EM算法从拟合的衰减曲线估计跨加速电压和检测器的所有组合的空气预滤器的校正谱。EM算法根据在来自步骤334的拟合的衰减曲线内的以“通常的设置”的空气预滤器谱，创建空气预滤器156的估计的校正谱。对于空气预滤器的估计的校正谱的创建的更多信息包括在下面的图4的描述中。
根据步骤340，该方法然后基于步骤336中创建的空气预滤器的估计的校正谱，跨(加速电压、预滤器、检测器)的所有组合计算新的基线谱集354’。新的基线谱集354’校正/考虑在步骤324中检测的X-射线源102的能谱的改变。
在完成步骤340时，该方法转换到步骤350。
在步骤350中，样品114旋转或移回到束104路径中，并且获取投影，直到在步骤352-1中的断层摄影扫描的结束或已经获取预定数量的投影。
图3B-2示出了在样品114的长扫描中间如何重新获取参考值。如图3B-1的方法中，步骤322以(能量、预滤器、检测器)设置的所选择的元组执行样品114的扫描，并且应用谱监测器360的步骤323、324以确定X-射线源102的能谱是否已经改变。如果谱已经改变，该方法转换到步骤328以重新校准X-射线源102。否则，该方法转换到步骤350并且将样品114旋转或移到束104路径中，并且该方法转换到步骤352-2以获取样品114的X-射线投影。
然而，与图3B-1的方法中的步骤351-1不同，由于较长的扫描持续时间，步骤352-2获取更多的样品114的X-射线图像投影。例如，在投影数量(预定数量)或预定扫描时间的阈值之后，样品144旋转或移出束路径，并且该方法转换回到步骤323以重新应用谱监测器360。在步骤323中，获取新的参考值并且测量某个预滤器的透射值。这里再次基于测量的透射值与将从基线谱集354预测的透射值进行比较，在步骤324中确定X-射线源谱是否已经改变。
如果X-射线源谱没有改变，在步骤350中样品在束路径中再次旋转或偏移，以在步骤352-2中继续收集预定数量的X-射线图像投影。
图4示出了用于计算图3A的步骤321中的基线谱354的所选择的加速电压、预滤器和检测器的组合的谱，用于计算步骤336中的空气预滤器的估计的校正谱和用于计算在上述图3B-1和图3B-2的步骤340中的新的基线谱354’的方法。
该方法可以诸如由计算机系统200以及可能由图像处理器220实现，作为软件程序以生成具有减少的伪影的断层摄影重建。在步骤340中，该计算将产生新的基线谱集354’，其将替换在图3A的方法中计算的原始基线谱354。然后，例如，在优选实施例中，新的基线谱354’可以通过图像处理器220用于样品114的断层摄影重建，其可以消除或显著减少由于束硬化引起的断层摄影重建中的伪影。
在步骤410中，获取AIR透射曲线测量结果。现在，假设所选择的加速电压+AIR+检测器的组合的谱(即AIR谱)是S0(E)，其中E表示能量。其结果是，在步骤412中，可以通过EM算法基于由在步骤414中的几个预滤器测量的衰减曲线估计该空气谱S0(E)。
例如，滤波轮中有N个预滤器。它们由不同的材料(或一些由相同的材料但具有不同的厚度制成)制成。它们的线性衰减系数分别为μi(E)，i＝1,2,3，...，N，并且它们的厚度分别为Ti，i 1,2,3，...，N。
假设入射X-射线的强度(即，每单位时间的光子计数)是I0，其可以通过使用滤波轮中的空(空气)点测量。通过旋转滤波轮，穿过某个滤波轮的下游X-射线的强度为I1，i＝1,2,3，…，N，，然后测量所选择的加速电压+AIR+检测器的组合的衰减曲线：
(I0/I0，I1/I0，I2/I0，…，Ii/I0，…，IN/I0)
S0(E)可以根据上述衰减曲线估计。
对于每个预滤器，提供S0(E)、I0、Ii、Ti和μi(E)之间的关系：
Ii＝I0∫S0(E)e-μi(E)TidE(1)
可以使用S0(E)来计算所选择的加速电压+预滤器+检测器(即预滤器谱)的所有组合的衰减曲线。
注意谱改变监测和重新校准的工作流程也可用于X-射线CT系统诊断。例如，如果谱监测器360确定X-射线源102的X-射线谱已经改变，这可以表示X-射线管和/或检测器没有有效地操作或具有问题。然后可以使用根据图4的方法创建的新的基线谱集354’以诊断X-射线CT系统100。
有两种不同的方式以获取如图4所示的预滤器谱，由标记为“间接方式”和“直接方式”的路径表示。
1)下面首先描述间接方式(更准确)。
根据步骤416，获取预滤器透射曲线计算。为此目的，以第k个预滤器为例：假设原始滤波轮中的第k个预滤器总是在束路径中，并且包括或可能由相同的预滤器集组成的另一个“虚拟”滤波轮作为谱测量工具，然后穿过原始滤波轮中第k个预滤器和“虚拟”滤波轮中的某个预滤器(由j表示)的下游X-射线的强度是：
Ikj,j＝0,1,2,3,…,N，即，
Figure GDA0003662177930000204
当j＝0时，其意味着“虚拟”滤波轮中的空(空气)点在束路径中并且μ0(E)＝0且T0＝0，然后Ik,j＝0意味着原始滤波轮中第k个预滤器的参考强度。
如果(2)除以(1)，即，
Figure GDA0003662177930000201
然后，
Figure GDA0003662177930000202
根据(3)，我们获取与所选择的加速电压+第k个预滤器+检测器的组合相关的衰减曲线：
(Ik,j＝0/Ik,j＝0,Ik,j＝1/Ik,j＝0,Ik,j＝2/Ik,j＝0,…,Ik,j/Ik,j＝0,…,Ik,j＝N/Ik,j＝0)
结果，在步骤418中，根据上述衰减曲线经由EM算法可以估计所选择的加速电压+第k个预滤器+检测器的组合的谱。
2)现在在下文中描述直接方式(近似的)。
步骤420描述了预滤器谱计算Sk(E)的计算。Sk(E)可以通过预滤器衰减的Sk(E)直接计算如下：
Figure GDA0003662177930000203
直接方式提供具有短得多的计算时间的近似(例如校正的谱354’)。
图5示出了衰减值的三个不同衰减曲线的图，其中测量值作为AIR(即，空气预滤器)测量的、HE1预滤器测量的和HE1预滤器计算的滤波器指数的函数。
X-射线源102使用的加速电压是70kVp。在滤波轮150中有16个滤波器。还通过另一滤波器集(例如用于HE1测量曲线)测量70kVp和HE1预滤器的组合的衰减曲线。
从70kVp和AIR(AIR测量曲线)的组合的衰减曲线计算70kVp和HE1(例如，HE1计算曲线)的组合的衰减曲线。显然，计算结果与测量结果非常一致。
图6是作为能量(kVp)的函数的不同估计的校正X-射线源谱601和602的图。估计的校正X-射线源谱601/602包括在部分地使用图5中示出的衰减曲线创建的新的基线谱集354’中。每个估计的校正X-射线源谱601/602绘制以不同能量/加速电压处以kVp为单位的X-射线光子的归一化计数。
估计的校正谱601基于70kVp和空气预滤器的组合的衰减值，并且其基于图5的“AIR测量的”衰减曲线估计。以类似的方式，基于图5中的“HE1计算曲线”估计70kVp和HE1预滤器602的组合的估计的校正谱。
图7示出了从其创建基线谱354的衰减曲线的示例。假设通过使用具有不同厚度的材料#2滤波器来测量以10kVp的间隔离散的加速电压，诸如70kVp，80kVp，...，150kVp的一些衰减曲线。得到的数据在作为以厘米为单位的滤波器厚度的函数的衰减的左侧图中示出。然后，对于更精细的离散的加速电压，诸如75kVp，89kVp，...，114kVp等，可以平滑地内插衰减曲线的图。这由作为加速电压和滤波器厚度的函数的衰减的右侧图示出。
从X-射线源102发射的X-射线的谱可以随着获取的X-射线图像投影而随时间改变。例如，X-射线源靶燃烧将改变谱。
图8示出了对于X-射线源102的衰减曲线的不同图，其中从X-射线源102发射的X-射线的束的谱确定为随着时间已经改变。衰减曲线包括“先前全测量结果”衰减曲线、“当前全测量结果”衰减曲线和拟合的衰减曲线。尽管衰减曲线在所有X-射线能量/加速电压上创建，但是仅示出了所选择的加速电压(这里为60kVp)的衰减曲线。
在X-射线CT系统102的初始启动期间，使用跨所有预滤器设置，加速电压和检测器组合获取的透射测量结果，例如经由图3A的方法，创建“先前全测量结果”衰减曲线。然后从“先前全测量结果”衰减曲线创建先前的基线谱354。
在确定X-射线源谱已经改变时，还使用跨所有预滤器设置、加速电压和检测器组合获取的透射测量结果创建“当前全测量结果”衰减曲线，并且从“当前全测量结果”衰减曲线创建当前的基线谱354。
重要的是注意，“当前全测量结果”衰减曲线和相关联的当前的基线谱354的创建不是必需的，并且这里仅为了参考和比较的目的而创建。具体地，当前的系统和方法通常必须执行该步骤以校正X-射线源102的X-射线源谱的改变。相反，根据图3B-1/3B-2的方法，操作者可以基于以通常的加速电压和检测器组合的少至2或3个预滤器获取的透射测量结果，估计跨所有预滤器设置，加速电压和检测器124组合的改变的X-射线谱的校正。
假设从先前的基线谱354提取60kVp+AIR+检测器组合的衰减曲线。它在图中被命名为“先前全测量结果”。在几天之后，从X-射线CT系统100的X-射线源102发射的束的谱确定已经改变。
然而，在图3B-1/3B-2的步骤332中仅需要2个，例如预滤器(LE1和HE5)的透射测量结果，以便估计跨预滤器设置、加速电压和检测器124的所有组合的X-射线束103中的X-射线的改变的谱的校正。根据先前的衰减曲线和两个预滤器测量值拟合60kVp+AIR+检测器的估计的衰减曲线，如图中曲线“具有2个滤波器测量结果的拟合曲线”所示。
重要的是注意“拟合曲线”衰减曲线图对应于“当前全测量结果”衰减曲线图是多么接近。拟合的曲线衰减图与对于相同所选择的加速电压的“当前全测量结果”衰减曲线几乎完全重叠。
尽管已经参照的优选实施例具体示出和描述了本发明，但是本领域技术人员将理解，在不脱离所附权利要求所涵盖的的本发明的范围的情况下，可以在形式和细节上进行各种改变。

Claims (20)

1.一种X-射线CT系统中的X-射线源谱测量和估计方法，所述方法包括：
使用滤光轮来测量并校准X-射线源谱，其中所述滤光轮包括滤波器，所述滤波器从发射的X-射线束移除不同能量的X-射线；
确定以至少一个X-射线加速电压、所述滤光轮的多于一个滤波器和至少一个检测器的多个不同组合的从X-射线CT系统的X-射线源发射的X-射线的基线谱；
在所述X-射线CT系统的操作期间，监测所述X-射线以确定所述X-射线源谱的改变；并且
响应于确定所述X-射线源谱已经改变而基于测量至少两个或更多个滤波器的测量透射值来计算新的基线谱。
2.根据权利要求1所述的方法，其特征在于，所述方法用于断层摄影重建和束硬化校正。
3.根据权利要求1所述的方法，其特征在于，所述方法用于多能量计算机断层摄影。
4.根据权利要求1所述的方法，其特征在于，所述方法用于X-射线CT系统诊断。
5.根据权利要求1所述的方法，其特征在于，通过以下来执行所述监测所述X-射线以确定所述X-射线源谱的改变：
使用X-射线加速电压、一个滤波器和检测器的组合来获取透射值；并且
将所获取的透射值与以所述X-射线加速电压、所述一个滤波器以及所述检测器的相同组合的所述基线谱内的透射值进行比较。
6.根据权利要求1所述的方法，其特征在于，确定所述基线谱包括：
使用X-射线加速电压、滤波器和检测器的组合来获取透射测量结果；
使用所述透射测量结果计算每个所述滤波器的透射曲线；并且
从所述透射曲线计算每个所述滤波器的衰减曲线。
7.根据权利要求1所述的方法，其特征在于，确定所述基线谱包括：
基于使用X-射线加速电压、滤波器和检测器的组合获取的透射测量结果计算每个所述滤波器的衰减曲线；以及
将期望最大化(EM)算法应用于所述衰减曲线。
8.根据权利要求1所述的方法，其特征在于，计算新的基线谱包括：
识别用于测量至少2或3个滤波器的透射值的加速电压和检测器的组合；
以所识别的加速电压和检测器组合来测量空气滤波器的透射值；以及
使用所述空气滤波器的所述透射值以及所述至少2或3个滤波器的所述透射值来创建透射曲线并从所述透射曲线来创建衰减曲线。
9.根据权利要求1所述的方法，其特征在于，计算新的基线谱包括：
从所述基线谱以通常的加速电压和检测器组合来选择空气滤波器谱；
针对所选择的空气滤波器谱来拟合衰减曲线；以及
使用EM算法，在所拟合的衰减曲线内以所述通常的加速电压和检测器组合根据所述空气滤波器谱来估计所述空气滤波器的校正谱。
10.根据权利要求1所述的方法，其特征在于，计算新的基线谱包括：
使用所测量的至少2或3个滤波器的透射值来创建空气滤波器的估计的校正谱；以及
从所述空气滤波器的所述估计的校正谱来计算所述新的基线谱。
11.根据权利要求1所述的方法，其特征在于，在没有体模存在的情况下完成基于测量至少2或3个滤波器的透射值来计算新的基线谱。
12.根据权利要求1所述的方法，其特征在于，还包括在样品的扫描期间使用所述基线谱以及所述新的基线谱以校正所述样品的断层摄影重建中的束硬化。
13.一种X-射线CT系统，所述系统包括：
滤光轮，包括用于从发射的X-射线束移除不同能量的X-射线的滤波器；
X-射线成像系统，包括X-射线源、滤波器以及检测器系统；以及
计算机系统，用于确定以X-射线源的至少一个X-射线加速电压、所述滤光轮的多于一个滤波器以及至少一个所述检测器系统的检测器的不同组合的从所述X-射线源发射的X-射线的基线谱；监测所述X-射线以确定所述X-射线源谱的改变；并且响应于确定所述X-射线源谱已经改变而基于测量至少两个或更多个滤波器的透射值来计算新的基线谱。
14.根据权利要求13所述的系统，其特征在于，所述新的基线谱用于断层摄影重建和束硬化校正。
15.根据权利要求13所述的系统，其特征在于，所述新的基线谱用于多能量计算机断层摄影。
16.根据权利要求13所述的系统，其特征在于，所述新的基线谱用于X-射线CT系统诊断。
17.根据权利要求13所述的系统，其特征在于，所述计算机系统通过以下来确定所述X-射线源谱的改变：
使用X-射线加速电压、一个滤波器以及检测器的组合来获取透射值；以及
将所获取的透射值与以所述X-射线加速电压、所述一个滤波器以及所述检测器的组合的在所述基线谱内的透射值进行比较。
18.根据权利要求13所述的系统，其特征在于，所述计算机系统通过以下来确定所述基线谱：
使用X-射线加速电压、滤波器和检测器的组合来获取透射测量结果；
使用所述透射测量结果计算每个所述滤波器的透射曲线；以及
从所述透射曲线计算每个所述滤波器的衰减曲线。
19.根据权利要求13所述的系统，其特征在于，所述计算机系统通过以下来确定所述基线谱：
基于使用X-射线加速电压、滤波器和检测器的组合获取的透射测量结果来计算每个所述滤波器的衰减曲线；以及
将期望最大化(EM)算法应用于所述衰减曲线。
20.根据权利要求13所述的系统，其特征在于，所述计算机系统通过以下确定所述基线谱：
识别用于测量至少2或3个滤波器的透射值的加速电压和检测器的组合；
以所识别的加速电压和检测器组合来测量空气滤波器的透射值；以及
使用所述空气滤波器的所述透射值以及所述2或3个滤波器的所述透射值来创建透射曲线并从所述透射曲线来创建衰减曲线。
CN201710178999.2A 2016-03-23 2017-03-23 计算机断层摄影x-射线显微镜系统中谱表征的方法和系统 Active CN107228865B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201662312042P 2016-03-23 2016-03-23
US62/312,042 2016-03-23

Publications (2)

Publication Number Publication Date
CN107228865A CN107228865A (zh) 2017-10-03
CN107228865B true CN107228865B (zh) 2022-08-26

Family

ID=58410160

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710178999.2A Active CN107228865B (zh) 2016-03-23 2017-03-23 计算机断层摄影x-射线显微镜系统中谱表征的方法和系统

Country Status (3)

Country Link
US (1) US10859515B2 (zh)
EP (1) EP3223002B1 (zh)
CN (1) CN107228865B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3223002B1 (en) 2016-03-23 2021-04-28 Carl Zeiss X-Ray Microscopy, Inc. Method, tool and computer program product for measurement and estimation of an x-ray source spectrum in computed tomography, and x-ray ct system
EP3465180A4 (en) * 2016-05-30 2020-03-04 Southern Innovation International Pty Ltd MATERIAL CHARACTERIZATION SYSTEM AND METHOD
US10685759B2 (en) * 2017-01-16 2020-06-16 Fei Company Statistical analysis in X-ray imaging
EP3428629B1 (en) * 2017-07-14 2022-12-07 Malvern Panalytical B.V. Analysis of x-ray spectra using curve fitting
JP7143567B2 (ja) * 2018-09-14 2022-09-29 株式会社島津テクノリサーチ 材料試験機および放射線ｃｔ装置
DE102018132025A1 (de) * 2018-12-13 2020-06-18 Intom GmbH Industrielles Röntgen-Werkstückvermessungssystem und Verfahren zum Betreiben eines solchen
US10942137B2 (en) * 2019-03-26 2021-03-09 The Government of the United States of America, as represented by the Secretary of Homeland Security Reconstruction of X-ray spectra using integrating detectors (RoXS-ID)
CN113391341B (zh) * 2021-05-25 2023-08-11 首都师范大学 一种考虑散射光子影响的x射线能谱估计方法
CN116831604A (zh) * 2021-12-06 2023-10-03 武汉联影生命科学仪器有限公司 数据校正方法
US20240044811A1 (en) * 2022-08-03 2024-02-08 Baker Hughes Holdings Llc Rapid High-Resolution Computerized Tomography

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104737005A (zh) * 2012-10-18 2015-06-24 卡尔蔡司X射线显微镜公司 具有晶粒取向映射功能的实验室x射线微断层扫描系统
CN104995690A (zh) * 2013-02-15 2015-10-21 卡尔蔡司X射线显微镜公司 多能量x射线显微镜数据采集及图像重建系统及方法
CN105073010A (zh) * 2013-04-04 2015-11-18 株式会社东芝 X射线计算机断层摄影装置

Family Cites Families (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002102217A (ja) * 2000-09-28 2002-04-09 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔシステム、ガントリ装置、コンソール端末及びその制御方法及び記憶媒体
US6748043B1 (en) * 2000-10-19 2004-06-08 Analogic Corporation Method and apparatus for stabilizing the measurement of CT numbers
US7050529B2 (en) 2002-07-23 2006-05-23 Ge Medical Systems Global Technolgy Company, Llc Methods and apparatus for performing a computed tomography scan
US6909771B2 (en) * 2002-11-22 2005-06-21 Board Of Regents, The University Of Texas System Three component x-ray bone densitometry
WO2005009206A2 (en) * 2003-06-25 2005-02-03 Besson Guy M Dynamic multi-spectral imaging system
US7086780B2 (en) * 2004-05-20 2006-08-08 General Electric Company Methods for spectrally calibrating CT imaging apparatus detectors
US7224763B2 (en) * 2004-07-27 2007-05-29 Analogic Corporation Method of and system for X-ray spectral correction in multi-energy computed tomography
US7583779B2 (en) * 2004-11-24 2009-09-01 General Electric Company System and method for acquisition and reconstruction of contrast-enhanced, artifact-reduced CT images
DE102005027436B4 (de) * 2005-06-14 2008-09-04 Siemens Ag Verfahren zur Berechnung von absorberspezifischen Gewichtungskoeffizienten und Verfahren zur Verbesserung eines von einem Absorber abhängigen Kontrast-zu-Rausch-Verhältnisses in einem von einer Röntgeneinrichtung erzeugten Röntgenbild eines zu untersuchenden Objektes
WO2007057841A2 (en) * 2005-11-18 2007-05-24 Koninklijke Philips Electronics N.V. Systems and methods using x-ray tube spectra for computed tomography applications
US7443953B1 (en) 2005-12-09 2008-10-28 Xradia, Inc. Structured anode X-ray source for X-ray microscopy
DE102007020065A1 (de) * 2007-04-27 2008-10-30 Siemens Ag Verfahren für die Erstellung von Massenbelegungsbildern anhand von in unterschiedlichen Energiebereichen aufgenommenen Schwächungsbildern
US8135205B2 (en) * 2008-12-20 2012-03-13 Varian Medical Systems, Inc. Rapid matching of noisy spectra or other signal patterns
KR20110055870A (ko) * 2009-11-20 2011-05-26 삼성전자주식회사 엑스선 디텍터의 에너지 응답 특성을 추정하는 방법 및 장치
FR2961904B1 (fr) * 2010-06-29 2012-08-17 Commissariat Energie Atomique Procede d'identification de materiaux a partir de radiographies x multi energies
US9693742B2 (en) * 2011-01-18 2017-07-04 Siemens Aktiengesellschaft Method for generating a contrast medium-assisted X-ray image and X-ray system
US20130156163A1 (en) * 2011-12-19 2013-06-20 General Electric Company Method and apparatus for reconstructing an image of an object
JP6462262B2 (ja) * 2013-08-30 2019-01-30 キヤノンメディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置及びフォトンカウンティングｃｔ装置
JP6289108B2 (ja) * 2014-01-14 2018-03-07 キヤノンメディカルシステムズ株式会社 フォトンカウンティングｃｔ装置及びフォトンカウンティングｃｔデータ処理方法
US9256938B2 (en) * 2014-02-26 2016-02-09 Kabushiki Kaisha Toshiba Characteristic X-ray escape correction in photon-counting detectors
FR3023001A1 (fr) * 2014-06-30 2016-01-01 Commissariat Energie Atomique Procede d'analyse d'un objet en deux temps utilisant un rayonnement en transmission puis un spectre en diffusion.
US10121267B2 (en) * 2014-07-03 2018-11-06 Duke University Spectral estimation and poly-energetic reconstruction methods and X-ray systems
US10159450B2 (en) * 2014-10-01 2018-12-25 Toshiba Medical Systems Corporation X-ray CT apparatus including a photon-counting detector, and an image processing apparatus and an image processing method for correcting detection signals detected by the photon-counting detector
US9775579B2 (en) * 2015-02-18 2017-10-03 Guy M. Besson Multi-source CT system and imaging method
US10646186B2 (en) * 2016-01-18 2020-05-12 Canon Medical Systems Corporation X-ray CT apparatus, information processing device and information processing method
EP3223002B1 (en) 2016-03-23 2021-04-28 Carl Zeiss X-Ray Microscopy, Inc. Method, tool and computer program product for measurement and estimation of an x-ray source spectrum in computed tomography, and x-ray ct system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104737005A (zh) * 2012-10-18 2015-06-24 卡尔蔡司X射线显微镜公司 具有晶粒取向映射功能的实验室x射线微断层扫描系统
CN104995690A (zh) * 2013-02-15 2015-10-21 卡尔蔡司X射线显微镜公司 多能量x射线显微镜数据采集及图像重建系统及方法
CN105073010A (zh) * 2013-04-04 2015-11-18 株式会社东芝 X射线计算机断层摄影装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A robust method of x-ray source spectrum estimation from transmission measurements: Demonstrated on computer simulated ,scatter-free transmission data;Emil Y. Sidky et.al;《Journal of Applied Physics》;AIP;20050616;第97卷;12 pages *

Also Published As

Publication number Publication date
US10859515B2 (en) 2020-12-08
EP3223002B1 (en) 2021-04-28
CN107228865A (zh) 2017-10-03
EP3223002A1 (en) 2017-09-27
US20170276620A1 (en) 2017-09-28

Similar Documents

Publication Publication Date Title
CN107228865B (zh) 计算机断层摄影x-射线显微镜系统中谱表征的方法和系统
US10297048B2 (en) Segmentation and spectrum based metal artifact reduction method and system
EP3076873B1 (en) System for generating spectral computed tomography projection data
US9383324B2 (en) Laboratory X-ray micro-tomography system with crystallographic grain orientation mapping capabilities
US8971488B2 (en) Systems and methods for detecting an image of an object using multi-beam imaging from an X-ray beam having a polychromatic distribution
US7409043B2 (en) Method and apparatus to control radiation tube focal spot size
EP1933170A1 (en) Method and system for computed tomography using transmission and fluorescence measurements
KR20080100200A (ko) 다색 분포를 가진 엑스선 빔의 사용에 의해 물체의 이미지를 검출하기 위한 시스템 및 방법
US8885910B2 (en) Systems and methods for X-ray imaging
Baldelli et al. A prototype of a quasi-monochromatic system for mammography applications
US7319733B2 (en) System and method for imaging using monoenergetic X-ray sources
US8100584B2 (en) Voltage measurement in an imaging system
CN107212898B (zh) 图像重建方法
CN108022272B (zh) 用于基于分割和光谱的金属伪影降低的计算机程序和ct系统
Boone et al. Full-field spectroscopic measurement of the X-ray beam from a multilayer monochromator using a hyperspectral X-ray camera
JP3208393U (ja) Ｘ線ｃｔシステム
Zhang et al. Improvements to conventional X-ray tube-based cone-beam computed tomography system
CA2575066A1 (en) Systems and methods for detecting an image of an object by use of an x-ray beam having a polychromatic distribution
JP7207856B2 (ja) Ｃｔ撮影装置
Krumm et al. Referenceless beam hardening correction in 3d computed tomography images of multi-material objects
JP4062232B2 (ja) Ｘ線ｃｔ装置及びｘ線ｃｔ装置による撮像方法
EP4177597A1 (en) Method and system for analytical x-ray calibration, reconstruction and indexing using simulation
CZ30426U1 (cs) Rentgenový CT systém
WO2024206092A1 (en) Scintillator defect correction in x-ray microscopy
JP2022128100A (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