CN104039227B - 在x射线系统中生成衰减图像数据和相位图像数据 - Google Patents
在x射线系统中生成衰减图像数据和相位图像数据 Download PDFInfo
- Publication number
- CN104039227B CN104039227B CN201280066654.3A CN201280066654A CN104039227B CN 104039227 B CN104039227 B CN 104039227B CN 201280066654 A CN201280066654 A CN 201280066654A CN 104039227 B CN104039227 B CN 104039227B
- Authority
- CN
- China
- Prior art keywords
- image data
- pixel
- ray
- measurement
- phase
- 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
- 238000005259 measurement Methods 0.000 claims abstract description 140
- 238000000034 method Methods 0.000 claims description 37
- 238000001228 spectrum Methods 0.000 claims description 32
- 238000012545 processing Methods 0.000 claims description 25
- 230000035945 sensitivity Effects 0.000 claims description 25
- 238000003384 imaging method Methods 0.000 claims description 10
- 241001269238 Data Species 0.000 claims 3
- 230000003595 spectral effect Effects 0.000 description 8
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000001419 dependent effect Effects 0.000 description 5
- 230000010363 phase shift Effects 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 239000006096 absorbing agent Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 101001044908 Cairina moschata Chymotrypsin inhibitor Proteins 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 238000012360 testing method 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/484—Diagnostic techniques involving phase contrast X-ray 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
在生成相位图像数据38中，接收包括第一逐像素测量信号值30a的第一X射线图像数据28a和包括第二逐像素测量信号值30b的第二X射线图像数据28b。以第一测量模式获得所述第一X射线图像数据28a，并且以不同于所述第一测量模式的第二测量模式同时获得所述第二X射线图像数据28b。通过根据以第一测量模式在像素31处获得的第一测量信号值30a和以第二测量模式在所述像素31处获得的第二测量信号值30b确定所述像素31处的相位值34来确定包括来自所述第一X射线图像数据28a和所述第二X射线图像数据28b的逐像素相位值34的相位图像数据38。
Description
技术领域
本发明涉及一种用于生成相位图像数据的方法、计算机程序以及计算机可读介质，并涉及一种X射线系统。
背景技术
X射线照相及断层摄影是用于多种应用的重要方法，例如大批量样品的非破坏性研究、工业产品的质量检验以及患者身体内部的感兴趣解剖结构和组织区域的无创检查。
基于X射线的衰减的X射线成像可以得到良好的结果，其中，诸如骨骼的高吸收解剖结构嵌在相对弱吸收材料的组织中。这是由于以下事实，即硬X射线束的穿透深度可以相当高，这允许记录衰减系数的尖锐突出。
在US2009/0304149Al和US2011/0216878Al中，示出了适于探测基于衰减的图像数据的X射线成像系统。
当对具有相似吸收截面的不同类型的组织进行检查时(例如在乳房摄影或血管造影中)，X射线吸收的对比度可能相对较差。在这种情况下，可以采用相位对比X射线照相及断层摄影，其中，对穿透感兴趣对象的X射线的相位变化进行检查。
存在若干方法来检索相位信息。目前，最流行的是Talbot干涉仪类型的方法，其可以是基于光栅的。但是，也存在替代方案，其中，同轴全息摄影方法是吸引人的，因为它可以不需要任何光学仪器，如需要仔细对准的单色仪或光栅。
在WO 2008/006470 Al中示出了用于相位对比成像的X射线干涉仪。
发明内容
本发明的目的之一是提供生成相位图像数据的另外的可能性。
该目的可以由独立权利要求的主题来实现。进一步的示范性实施例根据从属权利要求和以下的描述将是显而易见的。
本发明的一个方面涉及一种用于生成相位图像数据的方法。
根据本发明的实施例，该方法包括以下步骤：接收包括第一逐像素测量信号值的第一X射线图像数据；接收包括第二逐像素测量信号值的第二X射线图像数据；其中，所述第一X射线图像数据是以第一测量模式获得的，并且所述第二X射线图像数据是以不同于所述第一测量模式的第二测量模式获得的。通常，该方法可以包括以下步骤：接收至少两组X射线图像数据，其中，以不同的测量模式获得每个组。
该方法包括另外的步骤：通过根据以第一测量模式在像素处获得的第一测量信号值和以第二测量模式在该像素处获得的第二测量信号值确定该像素处的相位值来根据第一X射线图像数据和第二X射线图像数据确定包括逐像素相位值的相位图像数据。
X射线图像数据的相位和衰减信息通常混合在由探测器获得或探测到的强度信息或数据中，该强度信息或数据被存储在X射线图像数据中。以下可以被看作是本发明的主旨，即当以不同的测量模式例如以不同的能量或利用不同的探测器类型获得两组图像数据时，该混合信息或数据可以被分离。该分离可以基于以下事实来执行，即存在将像素处的相位值(和任选的衰减值)映射到强度值或测量信号值的函数关系(通常是可以利用数学模型确定的函数)。利用基于两种不同测量模式的两种测量，可以从该函数关系推导出两个独立的方程，这两个方程可以被求解以得到相位值和/或衰减值。
本发明的另一方面涉及一种用于生成相位图像数据的程序单元，其在由处理器执行时适于执行如上面和下面所描述的方法的步骤，并且涉及一种计算机可读介质，在该计算机可读介质中存储了这种计算机程序。计算机可读介质可以是软盘、硬盘、USB(通用串行总线)存储设备、RAM(随机存取存储器)、ROM(只读存储器)和EPROM(可擦除可编程只读存储器)。计算机可读介质还可以是数据通信网络，例如因特网，其允许下载程序代码。
本发明的又一方面涉及一种用于X射线系统的控制器，所述X射线系统适于执行如上面和下面所描述的方法。
本发明的再一个方面涉及一种X射线系统。根据本发明的实施例，该X射线系统包括X射线源、X射线探测器和控制器。该X射线系统适于以不同的测量模式获得X射线图像数据。该控制器适于基于以不同的测量模式获得的X射线图像数据生成相位图像数据。必须理解的是，如上面和下面描述的方法的特征可以是如上面和下面描述的系统或控制器的特征，反之亦然。
本发明的这些和其他方面将从下文描述的实施例变得显而易见，并将参照下文描述的实施例得以阐述。
附图说明
图1示意性示出根据本发明的实施例的X射线系统。
图2示意性示出根据本发明的实施例的处理单元。
图3示出根据本发明的实施例的用于控制X射线系统的流程图。
图4示出根据本发明的实施例的用于控制X射线系统的流程图。
图5示意性示出一种差分相位成像系统。
图6示出根据本发明的实施例的用于控制X射线系统的流程图。
原则上，在附图中为相同的部件提供相同的附图标记。
具体实施方式
图1示出包括探测器装置12和评估单元14的X射线系统10。探测器装置12包括X射线源16和X射线探测器18。对象20被放置在X射线源16与X射线探测器18之间，并且从X射线源16发出的X射线能够穿透对象20并落在探测器18上。在穿透对象20的过程中，X射线被衰减和移相。然而，探测器18可能只适于以逐像素的方式测量X射线的强度。例如，探测器18是线探测器或二维探测器。
可以包括工作站14的评估单元14包括适于控制X射线源16和探测器18的控制器22。例如，控制器22可以控制X射线源16的操作持续时间和/或从X射线源16发出的X射线的能级。此外，控制器22可以在测量过程中接收从X射线探测器18生成的图像数据28。
图像数据28可以在处理单元24中进行处理，并且可以显示在具有显示器的HMI 26上。X射线系统10的操作可以由操作人员经由HMI 26进行控制。
在下面，解释了允许分离图像数据中的吸收和相移的不同方法。通常，这可以通过以不同的测量模式获得图像数据28来实现。例如，图像数据28可以以不同的能量水平或利用不同的测量方法获得。
为了以不同的能量水平获得图象数据28，X射线源16可以适于以不同的能量水平生成X射线，和/或探测器18可以适用于区分不同的能量水平处的X射线。在后者的情况下，探测器18可以是分光探测器。
然而，也可以通过以不同的测量模式(特别是在一次测量过程中同时地)探测X射线来获得图像数据28。这可以通过使用具有不同谱加权的至少两种测量和/或通过使用适用于同时以18种测量模式工作的探测器18来实现。在这些情况下，可以使用非分光探测器18，例如计数与积分X射线(CIX)探测器18或积分和与坎贝尔处理(Campbelling)X射线(ICX)探测器18。后面的这两个范例可以比分光探测器更便宜。此外，这种探测器18可以产生具有相关噪声的图像数据28，这可以有利于吸收和相位对比度的分离。
US2009/0304149 A1示出适用于计数与积分的探测器。
US2011/0216878 Al示出适用于积分与坎贝尔处理的探测器。
谱加权并不局限于如分光探测器可能具有的高尖峰函数，而也可以使用像对ICX探测器的CIX有效的一般加权。因此，可以利用这些探测器的不同信道中的噪声相关性的益处。
图2示出处理单元24中的图像数据28的处理。所探测的图像数据28是由探测器18获得的，图像数据28包括与不同测量模式相关的至少两组图像数据28a、28b。每个图像数据28a、28b包括针对图像数据28a、28b中的每个像素31的测量信号值30a、30b。例如，测量信号值30a、30b可以指示探测器18的特定像素31处的X射线强度。
对于每个像素31，处理单元24基于测量信号值30a、30b确定衰减值32和/或相位值34，如下面将解释的。根据衰减值32和/或相位值34，由处理单元24重建衰减图像数据36和/或相位图像数据38。必须注意，处理单元24能够仅生成衰减图像数据36或相位图像数据38。
处理单元24可以包括具有处理器的计算机。在处理器上运行并存储在处理单元24的存储设备中的软件可以执行如上面和下面解释的方法。同样，图像数据28a、28b以及衰减图像数据36和相位图像数据38可以被存储在处理单元24的存储设备中。
根据本发明的实施例，用于生成衰减图像数据36和/或相位图像数据38的方法包括以下步骤：接收包括处理单元24中的第一逐像素测量信号值30a的第一X射线图像数据28a，以及接收包括第二逐像素测量信号值30b的第二X射线图像数据28b。
根据本发明的实施例，X射线系统10包括X射线源16、X射线探测器18和适于从相位图像数据中分离衰减图像数据的控制器22。X射线系统10可以适于以不同测量模式获得X射线图像数据28a、28b，例如，X射线系统10包括CIX或ICX探测器18。
根据本发明的实施例，X射线探测器18适于执行不同的测量方法，特别是同时执行不同的测量方法。
根据本发明的实施例，X射线探测器18适于以至少两种不同的方式对所探测的光子进行谱加权，例如通过对光子计数、对光子能量积分或对光子进行坎贝尔处理。坎贝尔处理可以通过探测光子流的波动并确定波动的平均值来执行，其为光子能量的平方的度量。
在下面将解释使衰减值32和相位值34与测量信号值30a、30b相关的函数关系所基于的理论背景。
强度I和相位在垂直于光轴A且紧随对象20之后(假设处于z＝0)的平面中的分布能够以对象20的复折射率n(x,y,z,λ)：n＝1–δ+iβ的线积分的形式来表达为
I(x,y,z＝0,λ)＝I0(x,y,z)e-M(x,y,z＝0，λ)
其中
M是依赖于折射率n的虚部β的衰减函数。I0(x,y,λ)是X射线源16的空间依赖主谱。
此外，波前的相位函数根据下式受到干扰：
并且依赖于折射率n的实部δ。应该注意到，对于折射率n，其实部δ以及其虚部β依赖于波长λ并因此依赖于X射线的能量。
相位函数并不是直接可探测的，但是在波前进一步沿光轴A传播的情况下将导致强度变化。
任务是根据在探测器平面z＝R处的测量结果检索M和两者。特定像素的衰减值32基于探测器平面中的像素的位置x,y处的衰减函数M，并且相位值34基于探测器平面中的像素的位置x,y处的相位函数
向探测器平面z＝R处的传播能够通过菲涅耳传播，或者近似通过强度方程的传递来描述：
其能够在强度在近场中的对象后面不剧烈变化的情况下被进一步简化。在以下假设条件下：
这导致
是已知的。在对象中没有吸收K边缘的情况下并且针对具有平均波长λ0的相当小的能量范围，这种依赖性由下式给出：
其中，对于康普顿效应能够被忽略的低能量来说p为3，并且能量越高，p越小。使用σ＝λ/λ0，γ＝Rλ/2π并且为了简洁而忽略空间坐标，则可以得到
该方程需要被再次转换以针对给定波长提供预期强度：
然后测量信号值30a、30b为：
其中，谱灵敏度函数Φ对探测器18的数据采集进行建模。指数i显示出针对不同的测量模式存在不同的函数关系。
例如，不同测量模式的最常见依赖性是：
项能够被省略，这得到
最后的函数是将衰减值32(即平均波长λ0处的M)和相位值34(即平均波长λ0处的)映射到测量信号S的函数关系或函数。该关系可以是逐像素评估的。
该函数关系依赖于谱灵敏度函数Φi(λ)和/或主谱I0(λ)。
谱灵敏度函数Φi(λ)可以由探测器18的校准测量获得。备选地或额外地，可以使用具有已知谱灵敏度函数Φi(λ)的探测器18(参见上述内容)。
同样，可以利用校准测量来确定主谱I0(λ)。
可能只需要S的两个独立测量值30a、30b来通过求解上述函数关系推导出两个量M(λ0)和
从至少两组图像数据28a、28b生成衰减图像数据36和相位图像数据38的不同可能性将参考图3和图4进行解释。
根据本发明的实施例，该方法包括以下步骤：通过根据以第一测量模式在像素31处获得的第一测量信号值30a和以第二测量模式在像素31处获得的第二测量信号值30b确定像素31处的相位值34，来根据第一X射线图像数据28a和第二X射线图像数据28b确定包括逐像素衰减值32的衰减图像数据36和/或包括逐像素相位值34的相位图像数据38。
根据本发明的实施例，图像数据36和/或38是基于将衰减值32和相位值34映射到第一测量信号值30a和第二测量信号值30b的函数关系获得的。该函数关系依赖于第一测量模式和第二测量模式，例如依赖于用于以各自模式进行测量的探测器18的谱灵敏度函数Φi。
根据本发明的实施例，该函数关系基于将衰减值32和/或相位值34映射到测量信号值S的逐像素函数。
参考图3，图像数据36、36是基于对上述函数关系的数值求解而确定的。
在步骤S10中，执行诊断扫描，其得到至少两组图像数据28a、28b。具体而言，第一图像数据28a是以第一测量模式获得的，并且第二图像数据28b是以第二测量模式获得的。例如，第一图像数据28a和第二图像数据28b可以利用计数探测器18、积分探测器或坎贝尔处理探测器18获得，并且第二图像数据28b利用积分探测器获得。
可以在不移动探测器18或X射线系统10的机械部件的情况下获得或探测第一X射线图像数据28a和第二X射线图像数据28b。可以同时探测两组图像数据28a、28b。
同样，可以例如例用分光镜探测器18以不同的能级获得第一图像数据28a和第二图像数据28b。
从临床角度来看，可能进一步期望的是同时采集对患者运动不灵敏的两种图像数据28a、28b。使用分光镜探测器或使用适于在一次测量过程中执行上述测量模式的探测器(如上述CIX或ICX探测器18)以不同的能量同时采集图像数据28a、28b是可能的。
根据本发明的实施例，该方法包括以下步骤：利用适于以不同谱加权探测光子的一个探测器18同时获得第一X射线图像数据28a和第二X射线图像数据28b。
根据本发明的实施例，第一X射线图像数据28a是以第一测量模式获得的，并且第二X射线图像数据28b是以与第一测量模式不同的第二测量模式获得的。
根据本发明的实施例，第一X射线图像数据28a是以光子能量的第一谱加权获得的，并且第二X射线图像数据28b是以与第一谱加权不同的第二谱加权获得的。谱加权可以被视为以不同方式考虑的针对具有不同能量的光子的测量方法。
根据本发明的实施例，第一谱加权和第二谱加权包括对光子进行计数、对光子能量进行积分或者确定光子能量波动。
在步骤S12中，图像数据28a、28b在处理单元24中被接收并且被处理以生成图像数据36、38。
针对每个像素，处理单元24将第一测量信号值30a、针对第一测量模式的谱灵敏度函数Φ1(λ)和针对第一能级的主谱I0(λ)插入到该函数关系中以获得针对两个未知值32、34的第一方程。
类似地，处理单元24将第二测量信号值30a、针对第二测量模式的谱灵敏度函数Φ2(λ)和针对第二能级的主谱I0(λ)插入到该函数关系中以获得针对两个未知值32、34的第二方程。
可以由处理单元24对这两个方程进行数值求解以确定值32、34。
根据本发明的实施例，该函数关系依赖于对光子能量进行加权的(谱)灵敏度函数Φi。第一测量模式和第二测量模式可以包括具有不同谱灵敏度函数Φi的不同谱加权。
根据本发明的实施例，相位值(34)是基于对第一测量模式和第二测量模式进行建模的至少两个(谱)灵敏度函数Φi(λ)确定的。具体而言，(谱)灵敏度函数Φi(λ)可以对在特定测量模式中光子能量的探测器18的加权进行建模。
根据本发明的实施例，衰减值32和相位值34是通过对该函数关系进行数值反转而确定的。
衰减图像数据36是根据所有像素的衰减值32生成的，并且相位图像数据38是根据所有像素的相位值34生成的。在此，拉普拉斯算子▽2可以被反转以得到相位图像数据38。
参考图4，图像数据34、36是基于经验模型确定的。
在步骤S20中，执行校准扫描。在该步骤中具有已知折射率n的对象被放置在X射线源16与探测器18之间并且执行测量。测量结果被评估并且根据这些测量结果确定经验模型。针对每种测量模式，经验模型可以对将第一测量信号值和第二测量信号值映射到衰减值和相位值的函数进行建模。该经验模型可以基于查找表。
可能必需针对在后续数据分析期间应该被评估的相同测量模式执行在校准扫描期间的测量。
应该注意到，校准扫描可以仅执行一次，而后续的诊断扫描可以利用在校准扫描期间生成的相同校准数据执行若干次。
在步骤S22中，类似于图3中的步骤S12执行诊断扫描，这得到以第一测量模式获得第一图像数据28a以及以第二测量模式获得第二图像数据28b。
在步骤S24中，图像数据28a、28b在处理单元24中被接收并且被处理以生成图像数据36、38。针对每个像素，第一测量值30a和第二测量值30b被输入到经验模型中以确定衰减值32和相位值34。
根据本发明的实施例，该函数关系基于使衰减值32和相位值34与基于测量模式的测量信号值30a、30b相互关联的经验模型。
根据本发明的实施例，该经验模型是利用校准测量确定的。
如上面和下面所描述的方法也可以利用图5所示的差分相位成像系统50来执行。该提出的方法的主要益处可以在于其在对象扫描期间不再需要机械相位步进。因此，扫描速度可以被提高一定的数量级。此外，该方法可以用于利用连续移动的源探测器布置的断层摄影。
图5示意性示出具有X射线源16和探测器18的差分相位成像系统10。
X射线源16可以包括诸如X射线管的非相干X射线源52以及用于实现空间波束相干的源栅格54。探测器18可以包括相位栅格56、吸收器栅格58以及适于通过从X射线源16发出穿过对象20的X射线来探测图像数据28a、28b的X射线探测器元件60。
相位栅格56用作相移分束器，其将X射线束的相前(phase front)的变化转换成强度调制，即具有在吸收器栅格58的位置处的相位栅格56的半间距的典型长度标度的干涉图案。吸收器栅格58生成来自离开相位栅格56的X射线的莫尔干涉图案。探测器元件60上的莫尔干涉图案包含关于在穿过对象20和相位栅格56两者后的经偏转和相移的X射线的相移的信息。
不同于移动栅格54、56、58之一，在下面描述了一种生成衰减图像数据36和相位图像数据38的方法，其使用不同的能级来获得X射线图像数据28a、28b。该方法可以被视为用于差分相位对比度CT的校准和处理方案。例如，可以使用能量选择性探测器元件60。
由于相位步进可能需要高机械精度并且可能增加扫描时间，避免机械步进可以得到更块的断层摄影设置。同样可以实现连续旋转的系统10。
在下面将解释关于DPCI系统10的理论背景，使衰减值32和相位值34与测量信号值30a、30b相关的函数关系基于该理论背景。
如前所述，折射率n依赖于波长和能量。
n＝1–δ+iβ
使用投影近似，根据下式由对象衰减X射线束的强度：
波前的相位根据下式被进一步干扰：
其中，所有积分都是沿着X射线的几何路径执行的。假设系统设计能量是E0。作为网格56和58的相对位置ξ的函数的每个像素31中的强度变化(归一化为1)可以被建模为：
I(ξ，E0)＝A0(cos(2πξ+Φ0)+1)
其中，Φ是梯度函数并且ξ是栅格取向。量A0等于平均强度I(E0)。
为了说明所提出的处理方案的基本思想，首先描述具有完美网格和相当低的X射线能量的理想情况(更复杂的通用方案将在后面概述)：首先回想在没有K-边缘的情况下并且针对相当小的能量(其中，康普顿散射能够被忽略)，折射率具有如下能量依赖性：
使用这些关系，能够总结出：
获知这一关系，能够根据不同能量处的观察强度估计感兴趣的量，即A0和Φ0，但是每次都有ξ＝0。注意到在上述方程中，量A0已经被归一化到空气(校准)扫描，以便补偿不同能量处的入射X射线束的不同强度。
在真实系统中，针对作为网格位移函数的强度依赖性的模型可能不会很好地满足。具体而言，X射线束可以是部分非相干的，网格56可能不是施加相移π的完美相位网格，网格58可能不会完美地吸收X射线，并且探测器60的谱响应可能不是完美的。
所有这些效应可以通过在校准扫描期间执行常规网格步进来补偿。更具体地，在该校准扫描期间，测量作为网格位移函数的强度的函数依赖性。最终的函数gE可以依赖于能量区间，能量区间由下标E来反映：
I(ξ，E)＝AgE(2πξ+Φ)
由于具有gE是周期性的的先验知识，将期望将测量的值扩展在有限傅里叶级数中。针对更高的能量可能出现另一个重要的难题，即吸收系数的能量依赖性变得更复杂。通常将总吸收分成光子效应和康普顿效应的贡献。
该方程提供使相位值与强度值30a、30b或测量信号值30a、30b相关的函数关系。该方程可以被视为将参考能量E处的期望相位梯度Φ0与测量数据I连接的前向模型。注意到ξ是网格56和58的相对位置并且可以被设置为0。
图6示出用于确定相位图像数据38的方法的流程图。
在步骤S30中执行校准扫描。在该扫描中，灵敏度函数gE被测量或确定。注意到在校准扫描中不需要移动栅格56、58。具体而言，通过利用不同能量处的强度测量结果拟合函数gE的傅里叶系数来确定有限傅里叶级数，可以确定函数gE。额外的拟合参数是参数B光子和B康普顿。这些参数和有限傅里叶级数可以被视为探测器18的经验模型。
在步骤S32中执行诊断扫描。在诊断扫描期间，测量针对一对能量Ei(例如至少两种能量Ei)的强度I(0,Ei)。针对每种能量Ei生成X射线图像数据28a、28b。注意到在诊断扫描中不需要移动栅格56、58。
根据本发明的实施例，以第一能级E1获得第一X射线图像数据28a，并且以与第一能级E1不同的第二能级E2获得第二X射线图像数据28b。
根据本发明的实施例，该方法包括以下步骤：以第一能级E1利用基于栅格的差分相位成像系统50获得第一X射线图像数据28a；以及以第二能级E2并且在相同的栅格位置处利用该基于栅格的差分相位成像系统50获得第二X射线图像数据28b。
在步骤S34中，诊断扫描期间生成的X射线图像数据28a、28b在处理单元24中被接收并且被评估以生成相位图像数据38。
具有所确定的函数gE的上述函数关系被用于根据以不同的能量Ei在像素31处获得的测量信号值30a、30b来确定像素31处的相位值34或Φ0。例如，最终的方程可以被数值反转。
根据本发明的实施例，相位值34是基于对探测器18的能量行为进行建模的(谱)灵敏度函数gE而确定的。
尽管已经在附图和前面的描述中详细说明和描述了本发明，但这种说明和描述被认为是说明性或示范性的而非限制性的；本发明并不局限于所公开的实施例。通过研究附图、说明书和权利要求书，本领域技术人员在实践所要求保护的本发明时能够理解和实现对所公开实施例的其他变型。在权利要求中，词语“包括”不排除其他元件或步骤，并且量词“一”或“一个”不排除多个。单个处理器或控制器或其他单元可以实现权利要求中记载的若干项的功能。在互不相同的从属权利要求中记载特定措施并不指示不能有利地使用这些措施的组合。权利要求书中的任何附图标记不得被解释为对范围的限制。
Claims (11)
1.一种用于生成相位图像数据(38)的方法，所述方法包括如下步骤：
接收包括第一逐像素测量信号值(30a)的第一X射线图像数据(28a)；
接收包括第二逐像素测量信号值(30b)的第二X射线图像数据(28b)；
其中，所述第一X射线图像数据(28a)是以第一测量模式获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一测量模式的第二测量模式获得的；
通过根据以第一测量模式在像素(31)处获得的第一测量信号值(30a)和以第二测量模式在所述像素(31)处获得的第二测量信号值(30b)确定所述像素(31)处的相位值(34)，来根据所述第一X射线图像数据(28a)和所述第二X射线图像数据(28b)确定包括逐像素相位值(34)的相位图像数据(38)；其中，X射线探测器(18)用于同时执行所述第一测量模式和所述第二测量模式；
其中，所述第一X射线图像数据(28a)是以光子能量的第一谱加权获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一谱加权的第二谱加权获得的，并且
其中，所述相位值(34)是基于对所述第一测量模式的能量行为进行建模的第一灵敏度函数和对所述第二测量模式的能量行为进行建模的第二灵敏度函数而确定的，所述第一灵敏度函数不同于所述第二灵敏度函数。
2.根据权利要求1所述的方法，
其中，所述第一谱加权和所述第二谱加权包括对光子进行计数、对光子能量进行积分或确定光子能量波动。
3.根据前述权利要求之一所述的方法，
其中，所述相位值(34)是通过对方程系进行数值反转而确定的，所述方程系包括针对所述相位值(34)和所述第一测量信号值(30a)的方程以及针对所述相位值(34)和所述第二测量信号值(30b)的方程。
4.根据权利要求1-2之一所述的方法，
其中，所述相位值是基于使所述相位值(34)与基于测量模式的测量信号值(30a、30b)相互关联的经验模型而确定的。
5.根据权利要求1-2之一所述的方法，还包括以下步骤：
根据以第一测量模式在像素(31)处获得的第一测量信号值(30a)和以第二测量模式在所述像素(31)处获得的第二测量信号值(30b)来确定所述像素(31)处的衰减值(32)。
6.根据权利要求1-2之一所述的方法，还包括：
利用适于以不同的谱加权探测光子的一个探测器(18)来同时获得所述第一X射线图像数据(28a)和所述第二X射线图像数据(28b)。
7.根据权利要求1-2之一所述的方法，
其中，所述第一X射线图像数据(28a)是以第一能量水平获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一能量水平的第二能量水平获得的。
8.根据权利要求1-2之一所述的方法，还包括：
利用基于光栅的差分相位成像系统以第一能量水平来获得所述第一X射线图像数据(28a)；
利用所述基于光栅的差分相位成像系统以第二能量水平并在相同的光栅位置处来获得所述第二X射线图像数据(28b)。
9.一种用于生成相位图像数据(38)的装置，所述装置包括：
用于接收包括第一逐像素测量信号值(30a)的第一X射线图像数据(28a)的单元；
用于接收包括第二逐像素测量信号值(30b)的第二X射线图像数据(28b)的单元；
其中，所述第一X射线图像数据(28a)是以第一测量模式获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一测量模式的第二测量模式获得的；
用于通过根据以第一测量模式在像素(31)处获得的第一测量信号值(30a)和以第二测量模式在所述像素(31)处获得的第二测量信号值(30b)确定所述像素(31)处的相位值(34)，来根据所述第一X射线图像数据(28a)和所述第二X射线图像数据(28b)确定包括逐像素相位值(34)的相位图像数据(38)的单元；其中，X射线探测器(18)用于同时执行所述第一测量模式和所述第二测量模式；
其中，所述第一X射线图像数据(28a)是以光子能量的第一谱加权获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一谱加权的第二谱加权获得的，并且
其中，所述相位值(34)是基于对所述第一测量模式的能量行为进行建模的第一灵敏度函数和对所述第二测量模式的能量行为进行建模的第二灵敏度函数而确定的，所述第一灵敏度函数不同于所述第二灵敏度函数。
10.一种X射线系统(10)，包括：
X射线源(16)；
X射线探测器(18)；
处理单元(24)；
其中，所述X射线系统(10)的所述X射线探测器适于(i)获得包括第一逐像素测量信号值(30a)的第一X射线图像数据(28a)；(ii)获得包括第二逐像素测量信号值(30b)的第二X射线图像数据(28b)；其中，所述第一X射线图像数据(28a)是以第一测量模式获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一测量模式的第二测量模式获得的；并且所述X射线系统(10)的所述处理单元适于通过根据以第一测量模式在像素(31)处获得的第一测量信号值(30a)和以第二测量模式在所述像素(31)处获得的第二测量信号值(30b)确定所述像素(31)处的相位值(34)，来根据所述第一X射线图像数据(28a)和所述第二X射线图像数据(28b)确定包括逐像素相位值(34)的相位图像数据(38)；
其中，所述处理单元适于基于以所述第一测量模式和所述第二测量模式获得的所述X射线图像数据(28a、28b)来生成相位图像数据(38)；
并且其中，所述X射线探测器(18)被配置为同时执行所述第一测量模式和所述第二测量模式；
其中，所述第一X射线图像数据(28a)是以光子能量的第一谱加权获得的，并且所述第二X射线图像数据(28b)是以不同于所述第一谱加权的第二谱加权获得的，并且
其中，所述相位值(34)是基于对所述第一测量模式的能量行为进行建模的第一灵敏度函数和对所述第二测量模式的能量行为进行建模的第二灵敏度函数而确定的，所述第一灵敏度函数不同于所述第二灵敏度函数。
11.根据权利要求10所述的X射线系统(10)，
其中，所述X射线探测器(18)适于执行不同的测量模式。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201261585761P | 2012-01-12 | 2012-01-12 | |
| US61/585,761 | 2012-01-12 | ||
| PCT/IB2012/057717 WO2013104966A1 (en) | 2012-01-12 | 2012-12-26 | Generating attenuation image data and phase image data in an x-ray system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104039227A CN104039227A (zh) | 2014-09-10 |
| CN104039227B true CN104039227B (zh) | 2018-02-02 |
Family
ID=47741197
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201280066654.3A Active CN104039227B (zh) | 2012-01-12 | 2012-12-26 | 在x射线系统中生成衰减图像数据和相位图像数据 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US9538970B2 (zh) |
| EP (1) | EP2802270B1 (zh) |
| JP (1) | JP6177800B2 (zh) |
| CN (1) | CN104039227B (zh) |
| RU (1) | RU2014133014A (zh) |
| WO (1) | WO2013104966A1 (zh) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103168228B (zh) * | 2010-10-19 | 2015-11-25 | 皇家飞利浦电子股份有限公司 | 微分相位对比成像 |
| US9001967B2 (en) * | 2012-12-28 | 2015-04-07 | Carestream Health, Inc. | Spectral grating-based differential phase contrast system for medical radiographic imaging |
| US10096098B2 (en) | 2013-12-30 | 2018-10-09 | Carestream Health, Inc. | Phase retrieval from differential phase contrast imaging |
| US10578563B2 (en) | 2012-12-21 | 2020-03-03 | Carestream Health, Inc. | Phase contrast imaging computed tomography scanner |
| US9724063B2 (en) | 2012-12-21 | 2017-08-08 | Carestream Health, Inc. | Surrogate phantom for differential phase contrast imaging |
| US9907524B2 (en) | 2012-12-21 | 2018-03-06 | Carestream Health, Inc. | Material decomposition technique using x-ray phase contrast imaging system |
| US9357975B2 (en) | 2013-12-30 | 2016-06-07 | Carestream Health, Inc. | Large FOV phase contrast imaging based on detuned configuration including acquisition and reconstruction techniques |
| US9700267B2 (en) | 2012-12-21 | 2017-07-11 | Carestream Health, Inc. | Method and apparatus for fabrication and tuning of grating-based differential phase contrast imaging system |
| US9494534B2 (en) | 2012-12-21 | 2016-11-15 | Carestream Health, Inc. | Material differentiation with phase contrast imaging |
| US9155510B2 (en) * | 2013-10-12 | 2015-10-13 | Wisconsin Alumni Research Foundation | Systems and methods for generating x-ray phase contrast images using a conventional x-ray imaging system |
| JP6495945B2 (ja) | 2014-06-02 | 2019-04-03 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 微分位相コントラストイメージングにおけるスペクトル位相アンラッピングのためのバイアスフリー正則化 |
| US10426417B2 (en) * | 2014-06-16 | 2019-10-01 | Koninklijke Philips N.V. | Computed tomography (CT) hybrid data acquisition |
| CN109073613B (zh) * | 2016-03-31 | 2021-04-27 | 株式会社岛津制作所 | 峰检测方法以及数据处理装置 |
| EP3541285B1 (en) * | 2016-11-16 | 2020-06-24 | Koninklijke Philips N.V. | Apparatus for generating multi energy data from phase contrast imaging data |
| EP3378396A1 (en) | 2017-03-23 | 2018-09-26 | Koninklijke Philips N.V. | X-ray imaging data processing device and method |
| JP7069670B2 (ja) * | 2017-12-04 | 2022-05-18 | コニカミノルタ株式会社 | Ｘ線撮影システム |
| US10902650B2 (en) * | 2018-04-19 | 2021-01-26 | Fei Company | X-ray beam-hardening correction in tomographic reconstruction using the Alvarez-Macovski attenuation model |
Family Cites Families (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5802137A (en) * | 1993-08-16 | 1998-09-01 | Commonwealth Scientific And Industrial Research | X-ray optics, especially for phase contrast imaging |
| EP0954951B1 (en) * | 1996-12-24 | 2009-02-04 | XRT Limited | Phase retrieval in phase contrast imaging |
| US5812629A (en) * | 1997-04-30 | 1998-09-22 | Clauser; John F. | Ultrahigh resolution interferometric x-ray imaging |
| JP4118535B2 (ja) * | 2001-07-03 | 2008-07-16 | 株式会社日立メディコ | Ｘ線検査装置 |
| EP1447046A1 (en) * | 2003-02-14 | 2004-08-18 | Paul Scherrer Institut | Apparatus and method to obtain phase contrast x-ray images |
| JP2006334069A (ja) * | 2005-06-01 | 2006-12-14 | Toshiba Corp | Ｘ線画像撮影方法および装置 |
| JP2009502227A (ja) * | 2005-07-22 | 2009-01-29 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 多色スペクトルによるｘ線検出器イメージング |
| DE102006015356B4 (de) * | 2006-02-01 | 2016-09-22 | Siemens Healthcare Gmbh | Verfahren zur Erzeugung projektiver und tomographischer Phasenkontrastaufnahmen mit einem Röntgen-System |
| DE102006063048B3 (de) * | 2006-02-01 | 2018-03-29 | Siemens Healthcare Gmbh | Fokus/Detektor-System einer Röntgenapparatur zur Erzeugung von Phasenkontrastaufnahmen |
| EP1879020A1 (en) | 2006-07-12 | 2008-01-16 | Paul Scherrer Institut | X-ray interferometer for phase contrast imaging |
| CN101413905B (zh) * | 2008-10-10 | 2011-03-16 | 深圳大学 | X射线微分干涉相衬成像系统 |
| EP2370836B1 (en) | 2008-11-25 | 2018-05-30 | Koninklijke Philips N.V. | Spectral imaging |
| JP2010188035A (ja) * | 2009-02-20 | 2010-09-02 | Fujifilm Corp | Ｘ線画像撮影システムおよびｘ線画像撮影装置 |
| DE102009019514A1 (de) * | 2009-04-30 | 2010-11-11 | Siemens Aktiengesellschaft | CT-System und Verfahren zur Phasenkontrast- und Absorptionsbildgebung |
| JP2012143548A (ja) * | 2010-12-21 | 2012-08-02 | Fujifilm Corp | 放射線画像取得方法および放射線画像撮影装置 |
-
2012
- 2012-12-26 RU RU2014133014A patent/RU2014133014A/ru not_active Application Discontinuation
- 2012-12-26 WO PCT/IB2012/057717 patent/WO2013104966A1/en active Application Filing
- 2012-12-26 US US14/370,301 patent/US9538970B2/en active Active
- 2012-12-26 EP EP12824918.2A patent/EP2802270B1/en active Active
- 2012-12-26 CN CN201280066654.3A patent/CN104039227B/zh active Active
- 2012-12-26 JP JP2014551687A patent/JP6177800B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US20140355740A1 (en) | 2014-12-04 |
| JP6177800B2 (ja) | 2017-08-09 |
| WO2013104966A1 (en) | 2013-07-18 |
| EP2802270A1 (en) | 2014-11-19 |
| EP2802270B1 (en) | 2017-10-11 |
| CN104039227A (zh) | 2014-09-10 |
| US9538970B2 (en) | 2017-01-10 |
| RU2014133014A (ru) | 2016-03-10 |
| JP2015503430A (ja) | 2015-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104039227B (zh) | 在x射线系统中生成衰减图像数据和相位图像数据 | |
| US9795350B2 (en) | Material differentiation with phase contrast imaging | |
| CN105452852B (zh) | 通过使用相衬ct的单色衰减对比图像生成 | |
| JP6590381B2 (ja) | Ｘ線装置、データ処理装置及びデータ処理方法 | |
| US8005185B2 (en) | Method to determine phase and/or amplitude between interfering, adjacent x-ray beams in a detector pixel in a talbot interferometer | |
| JP6105586B2 (ja) | エネルギー高感度検出の微分位相コントラストイメージング | |
| JP5127247B2 (ja) | Ｘ線装置の焦点‐検出器装置 | |
| US20160309068A1 (en) | Spatial frequency domain imaging using custom patterns | |
| CN105628718A (zh) | 多能谱x射线光栅成像系统与成像方法 | |
| JP2009525084A (ja) | Ｘ線システムを用いた投影画像および断層撮影画像の作成方法 | |
| CN110770607B (zh) | 处理光子计数型x射线检测数据的方法及装置、以及x射线装置 | |
| EP3170148B1 (en) | Iterative reconstruction method for spectral, phase-contrast imaging | |
| JP6168753B2 (ja) | 演算装置、演算プログラム、ｘ線測定システム、およびｘ線測定方法 | |
| JP2016506267A (ja) | 画像処理装置及び画像をフィルタリングするための方法 | |
| EP3541285B1 (en) | Apparatus for generating multi energy data from phase contrast imaging data | |
| WO2021206076A1 (ja) | 計測信号処理装置、計測信号処理方法およびプログラム | |
| CN114650774B (zh) | 用于处理由暗场和/或相衬x射线成像系统采集的数据的装置 | |
| EP3468472B1 (en) | Phantom device, dark field imaging system and method for acquiring a dark field image | |
| JP2014102087A (ja) | 光学的測定方法 | |
| JP2015040744A (ja) | 光学的測定方法 | |
| CN105916449B (zh) | Mri指标的推算方法以及生物体测定装置 | |
| JP2014092438A (ja) | 光学的測定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |