CN109414237A - X射线成像数据处理设备和方法 - Google Patents
X射线成像数据处理设备和方法 Download PDFInfo
- Publication number
- CN109414237A CN109414237A CN201880002656.3A CN201880002656A CN109414237A CN 109414237 A CN109414237 A CN 109414237A CN 201880002656 A CN201880002656 A CN 201880002656A CN 109414237 A CN109414237 A CN 109414237A
- Authority
- CN
- China
- Prior art keywords
- branch mailbox
- energy branch
- scale factor
- phase shift
- ray
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 62
- 238000012545 processing Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 title claims description 14
- 230000010363 phase shift Effects 0.000 claims abstract description 78
- 238000005259 measurement Methods 0.000 claims abstract description 9
- 238000003672 processing method Methods 0.000 claims description 6
- 238000013508 migration Methods 0.000 claims description 5
- 230000001737 promoting effect Effects 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 abstract description 6
- 238000006073 displacement reaction Methods 0.000 abstract description 4
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 230000002349 favourable effect Effects 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000000454 anti-cipatory effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 230000010358 mechanical oscillation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum 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/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
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- 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
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
Abstract
本发明涉及处理X射线成像数据,特别是相衬X射线成像数据,并且通过利用光谱信息解决与例如基于光栅的差分相衬成像中的振动相关联的问题。已经认识到可以利用光谱信息来基于多能量分箱X射线探测器、振动引起的莫尔图案相移的能量独立性、物体引起的莫尔图案相移的能量依赖性以及被设计为允许测量物体根据能量分箱引起的相移的依赖性的体模来确定莫尔图案的振动状态。特别地,发现通过适当的数据处理,能够从物体属性中解脱出振动引起的移位等,因此减少了图像中的振动伪影,同时并不特别需要改变成像设备的基本物理结构。
Description
技术领域
本发明涉及处理X射线成像数据,特别是相衬X射线成像数据。
背景技术
基于光栅的差分相衬成像近年来获得了增长的势头。这种技术能够提供与常规的衰减图像互补的额外信息,即,折射信号和暗场信号。其中,全景数字乳腺摄影(FFDM)是该技术的有前景的应用。最近,第一个临床FFDM设备成功配备了基于飞利浦MicroDose L30的X射线光栅硬件(参见T.Koehler等人的“Slit-scanning differential x-ray phase-contrast mammography:Proof-of-concept experimental studies”(Medical Physics42,第1959页,2015年)。第一个结果表明,狭缝扫描系统满足FOV、扫描时间和剂量方面的相关要求。
然而,差分相衬成像容易受到机械不稳定性和振动的影响。为了提取物体参数(衰减信号、折射信号和暗场信号),通常通过考虑没有物体的参考扫描(“空白扫描”,参见上文引用的T.Koehler等人的文章)来处理记录的数据。因此,所得到的图像质量受到重复扫描期间获得的莫尔图案的可重复性的限制。常规地,例如由空白扫描与物体扫描之间的振动引起的偏差被错误地归因于物体并导致图像伪影。
例如在DE 102014213817 A1中讨论了如下尝试,其旨在借助于对成像装置进行结构修改来减少振动等的发生。然而,这种结构修改会导致更复杂的机械装置(即,更高的成本),而即使复杂的尝试也不能完全消除不希望的影响。
发明内容
本发明的目的是提供一种用于解决与特别是基于光栅的差分相衬成像中的振动等相关联的上述问题的方法。
在本发明的第一方面中,提出了一种用于能量分箱的差分相衬X射线成像数据的处理设备,其包括:比例因子提供单元,其被布置为提供针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体的引起的相移之间的比例因子;以及暗场和相移确定单元,其被布置为使用由所述比例因子提供单元提供的所述比例因子根据被包括在所述成像数据中的、针对所述第一能量分箱的像素的探测强度值以及针对所述第二能量分箱的所述像素的探测强度值来确定暗场信号和物体引起的相移中的至少一项。在要求保护的发明的背景中,术语“暗场和相移确定单元”还可以意指被布置为既确定暗场信号又确定物体引起的相移的单个单元,或者被布置为在单独的子单元中独立地确定暗场信号和物体引起的相移的确定单元。
在本发明的第二方面中,提出了一种比例因子确定设备,其用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,所述比例因子确定设备包括:计算单元,其被布置为:使用所述第一能量分箱、具有预定形状的体模根据通过差分相衬X射线成像获得的相位图像来计算针对所述第一能量分箱的第一相移,并且使用所述第二能量分箱、所述体模根据通过差分相衬X射线成像获得的相位图像来计算针对所述第二能量分箱的第二相移;以及获得单元,其被布置为根据所述第一相移和所述第二相移来获得所述比例因子。
在本发明的第三方面中,提出了一种用于能量分箱的差分相衬X射线成像数据的处理方法,其包括:比例因子提供步骤,其提供针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子;以及确定步骤,其使用通过所述比例因子提供步骤提供的所述比例因子根据被包括在所述成像数据中的、针对所述第一能量分箱的像素的探测强度值以及针对所述第二能量分箱的所述像素的探测强度值来确定暗场信号和物体引起的相移中的至少一项。
在本发明的第四方面中,提出了一种比例因子确定方法,其用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,所述比例因子确定方法包括:计算步骤,其使用所述第一能量分箱、具有预定形状的体模根据通过差分相衬X射线成像获得的相位图像来计算针对所述第一能量分箱的第一相移,并且使用所述第二能量分箱、所述体模根据通过差分相衬X射线成像获得的相位图像来计算针对所述第二能量分箱的第二相移;以及获得步骤,其根据所述第一相移和所述第二相移来获得所述比例因子。
本发明通过利用光谱信息解决了与例如基于光栅的差分相衬成像中的振动相关联的问题。发明人已经认识到可以利用光谱信息来基于多能量分箱X射线探测器、振动引起的莫尔图案相移的能量独立性、物体引起的莫尔图案相移的能量依赖性和被设计为允许测量物体根据能量分箱引起的相移的依赖性的体模来确定莫尔图案的振动状态。
特别地,发现通过适当的数据处理,能够从物体属性中解脱出振动引起的移位等,因此减少了图像中的振动伪影,同时并不特别需要改变成像设备的基本物理结构。
这里应当注意,本发明不限于仅两个能量分箱的情况,因为也预想到多于两个能量分箱。当考虑相应的探测强度值时,可以使用三个或甚至更多个能量分箱来确定因振动引起的相移,同时使用另外的能量分箱来提供允许单独(稍微独立)确定的益处,例如,通过考虑第一能量分箱和第二能量分箱,考虑第一能量分箱和(另外的)第三能量分箱以及分别考虑第二能量分箱和第三能量分箱,如本文中仅针对两个能量分箱明确讨论的情况。而且,可以使用两个能量分箱来识别因振动引起的相移,然后将上述相移用于使用第三能量分箱的同时差分相衬成像中的补偿。
在优选实施例中,所述比例因子提供单元被布置为根据针对所述第一能量分箱的衰减值和针对所述第二能量分箱的衰减值来提供所述比例因子。
在提供分别根据第一能量分箱的衰减值和根据第二能量分箱的衰减值的比例因子时,可以实现提高的准确度,同时(根据情况)可以仅提供因子的近似值,这并不改变不同的衰减(例如在只有很小的影响的情况下)。
在上述实施例的优选修改方案中,所述比例因子提供单元包括查找表。
通过使用衰减值作为查找参数,查找表允许快速访问比例因子。尽管如此,提供因子的其他方式也是可能的,例如提供用于计算的公式或计算与查找表的组合。
在优选实施例中,所述暗场和相移确定单元被布置为根据所述第一能量分箱的测量数据mi,1和所述第二能量分箱的测量数据mi,2来确定针对所述第一能量分箱的衰减值A1、针对所述第二能量分箱的衰减值A2、针对所述第一能量分箱的暗场信号D1、针对所述第二能量分箱的暗场信号D2、物体引起的相移φ1,以及干扰引起的相移ψi,从而使针对所选择的图像值获得所述测量数据的可能性最大化。
这包括例如使惩罚测量数据与预期数据之间的差异的成本函数最小化。
在上述实施例的优选修改方案中,所述暗场和相移确定单元被布置为采用根据下式的最小二乘最小化:
Δ(A1,A2,D1,D2,φ1,ψi)=∑(mi,1-A1Bi,1(1+D1Vi,1cos(αi,1+ψi+φ1)))2+∑(mi,2-A2Bi,2(1+D2Vi,2cos(αi,2+ψi+c·φ1)))2,
其中,Bi,1、Bi,2、Vi,1、Vi,2、αi,1和αi,2分别是在空白扫描期间获得的莫尔图案的强度、可见度和相位。
然而,使用如最小二乘最小化的方法只是一种选择。如在给定情况(处理光子计数,使负泊松对数似然最小化)下,如果空白扫描强度变化,则可能还想要使用加权最小二乘成本函数。
在可以特别考虑作为上述实施例的优选修改方案的另外的优选实施例中,所述暗场和相移确定单元被布置为在所述确定中(例如在所述最小二乘最小化中)包括关于促使相邻像素的干扰引起的相移之间的差异小的另外的限制。特别地,特别是在最小二乘最小化的背景中,这可能对应于关于相邻像素的干扰引起的相移之间的绝对差值小于预定值的限制。
可以认为具有相邻像素将遭受类似的干扰或振动的先验知识,使得这种信息可以被包括在所述确定中。可以注意到,振动可以改变光栅相对于彼此的相对位置,其中,进一步地,干涉仪可以对这种变化非常敏感(发现即使在μm范围内的变化也可以产生可见效果)。
在优选实施例中,提供了一种比例因子确定系统,其用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,并且所述比例因子确定系统包括:根据本发明的比例因子确定设备;以及体模,其中,所述体模包括多个部分,每个部分具有以预定方式在横向于所述X射线的射束方向的方向上改变的高度,所述体模的所述部分的高度是所述射束方向上的部分的延伸,其中,所述体模的不同部分具有不同的平均高度,其中,所述比例因子确定设备被布置为根据针对所述第一能量分箱的衰减值和针对所述第二能量分箱的衰减值来确定所述比例因子。
在优选实施例中,提供了一种差分相衬X射线成像系统,并且所述差分相衬X射线成像系统包括:差分相衬成像单元;X射线源,其被布置为提供第一能量分箱的X射线和第二能量分箱的X射线;根据本发明的处理设备,其被布置为接收来自所述差分相衬成像单元的能量分箱的差分相衬X射线成像数据。
在上述实施例的优选修改方案中,所述成像系统还包括根据上述实施例的比例因子确定系统,其中,所述比例因子确定系统被布置为向所述处理设备的所述比例因子提供单元提供关于所述比例因子的信息。
预想到这样的情况,其中,比例因子可以是恒定的(或者至少足够接近恒定值(或值的集合)),因此不需要特定的装置来确定该值(除了用于存储预定值或值的集合的单元),而在其他情况下,用于专门确定关于比例因子信息的比例因子确定系统是有利的,例如在能量分箱尺寸太大而使得由这些箱测得的有效能量因物体的射束硬化而明显改变的情况下就是如此。然而,对于窄的能量分箱,仅使用公式(5)(参见下文)来计算比例因子可能就足够了。
在本发明的另外的方面中,提出了一种用于能量分箱的差分相衬X射线成像数据的处理方法的软件产品,所述软件产品包括程序代码模块,所述程序代码模块用于当所述软件产品在计算机上运行时使得所述计算机执行根据本发明的方法的步骤,并且提出了一种用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子的软件产品,所述软件产品包括程序代码模块,所述程序代码模块用于当所述软件产品在计算机上运行时使得所述计算机执行根据本发明的方法的步骤。
应当理解,根据权利要求1所述的处理设备、根据权利要求7所述的确定单元、根据权利要求11所述的处理方法、根据权利要求12所述的确定方法以及根据权利要求13和14所述的计算机程序具有相似和/或相同的优选实施例,特别地,如从属权利要求中所定义的实施例。
应当理解,本发明的优选实施例也能够是从属权利要求或上述实施例与相应的独立权利要求的任意组合。
参考下文描述的实施例,本发明的这些方面和其他方面将变得明显并且得到阐明。
附图说明
在以下图中:
图1示出了根据本发明的实施例的差分相衬X射线成像系统,并且
图2示出了图示根据本发明的实施例的用于能量分箱的差分相衬X射线成像数据的处理方法的流程图。
具体实施方式
图1示出了根据本发明的实施例的差分相衬X射线成像系统。
系统10包括X射线源15和对应的差分相衬成像单元20,X射线源15被布置为提供第一能量分箱的X射线和第二能量分箱的X射线,例如包括第一能量范围的光子和第二能量范围的光子。
在图1中,在X射线源15与成像单元20之间提供了体模25,而在用于对物体成像的系统10的操作中,在X射线源15与成像单元20之间提供了物体(未示出)而不是体模25。
体模25包括若干不同高度的支柱25a。这些支柱25a的顶部是三角形,所述三角形在射束方向上具有恒定的投射材料高度斜率。应当注意,(通过射束方向上的材料高度实现的)投影电子密度的恒定梯度产生莫尔图案的恒定相移。
另外,在系统10中,成像单元20被耦合到比例因子确定设备30和处理设备35。
比例因子确定设备30包括计算单元40和获得单元45,其操作在下面讨论。
处理单元35包括比例因子提供单元50以及暗场信号和相移确定单元55,其操作从下面的讨论中也变得清楚。
比例因子提供单元50被耦合到比例因子确定设备30,以从比例因子确定设备30接收关于比例因子的数据或信息。
图2示出了图示根据本发明的实施例的用于能量分箱的差分相衬X射线成像数据的处理方法的流程图。
如下面还要讨论的,该实施例的方法包括计算步骤60,其中,根据相位图像来计算针对第一能量分箱的第一相移,所述相位图像是通过对体模进行差分相衬X射线成像而获得的,如图1所示,其中,相位图像使用或对应于第一能量分箱。计算步骤60还包括根据相位图像来计算针对第二能量分箱的第二相移,所述相位图像是通过对体模进行差分相衬X射线成像而获得的,这里使用第二能量分箱。
基于在计算步骤60中计算出的信息,后续的获得步骤65包括根据第一相移和第二相移来获得比例因子,同时额外地基于体模来考虑衰减的影响。
在比例因子提供步骤70中提供所获得的比例因子,并且所获得的比例因子指示针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例。
最后,提供确定步骤75,其使用通过比例因子提供步骤70提供的比例因子根据被包括在成像数据中的、针对第一能量分箱的像素的探测强度值以及针对第二能量分箱的该像素的探测强度值来确定暗场信号和物体引起的相移中的至少一项。
在下文中,进一步讨论了本发明的实施例的上述概述。
针对普通几何射线的像素i,用于在Talbot-Lau干涉仪后面测量的强度的常用模型是:
Ii=ABi(1+DVicos(αi+φ))(1)
其中,Bi、Vi和αi分别指代在空白扫描期间获得的莫尔图案的强度、可见度和相位。物体参数A、D和φ分别对应于衰减、暗场信号和相移。通过使用测量数据mi的最小二乘最小化来获得物体参数:
Δ(A,D,φ)=∑(mi-ABi(1+DVicos(αi+φ)))2 (2)
然而,如上所述,机械振动产生对应于下式的莫尔图案的额外相移ψi:
Ii=ABi(1+DVicos(αi+ψi+φ)) (3)
利用非光谱数据累积,在最小化期间不能分离φ与ψi。然而,使用具有两个能量分箱(e=[1,2])的X射线探测器,公式(3)为:
Ii,e=AeBi,e(1+DeVi,ecos(αi,e+ψi+φe)) (4)
注意,除了振动移位ψi以外,所有参数都取决于光子能量。
对于单能量X射线光子,相移φ的能量依赖性产生:
对于多色X射线,比例因子c经由下式描述了在两个能量分箱中获得的相移φ1与φ2之间的关系:
φ1=c(s1,s2)·φ2 (6)
并且比例因子c取决于对应于两个能量分箱的有效能谱s1和s2。它能够使用图1中以草图描绘的模型进行测量。
单独对两个能量分箱使用公式(2)产生针对两个箱的相位图像。通过对与相位图像中的各个支柱相对应的区求平均,经由公式(6)获得针对每个支柱的c。由于不同的高度,支柱产生不同程度的射束硬化,其被编码在能够从衰减图像中提取的针对箱A1、A2的衰减信号中。遵循该策略,可以获得针对c(A1,A2)的查找表。
在略微不同的实施例(未图示)中,与单能量X射线光子的情况相似,也针对对应于能量分箱的平均能量的商的多色X射线,比例因子c相等,而相应地也可以获得针对c(A1,A2)的这种查找表。
为了从物体参数中解脱出振动引起的移位,我们为两个能量分箱明确地写出公式(3)并利用公式(6)。
Ii,1=A1Bi,1(1+D1Vi,1cos(αi,1+ψi+φ1)) (7)
Ii,2=A2Bi,2(1+D2Vi,2cos(αi,2+ψi+c(A1,A2)·φ1)) (8)
通过使用来自两个箱的测量数据mi,e的最小二乘最小化,能够从振动分量ψi中解脱出物体相位φ1:
Δ(A1,A2,D1,D2,φ1,ψi)=∑(mi,1-A1Bi,1(1+D1Vi,1cos(αi,1+ψi+φ1)))2+∑(mi,2-A2Bi,2(1+D2Vi,2cos(αi,2+ψi+c·φ1)))2 (9)
注意,上面提供了已经足以在每个像素的基础上分离物体引起的相移与振动引起的相移的信息。然而,通过增加一些另外的先验知识能够使估计更加鲁棒,即,如果数据是与相邻探测器像素同时采集的,则振动引起的相位对于相邻像素应当相同(或几乎相同)。如果用x对垂直于扫描方向的图像中的一行中的像素进行索引,则最小化能够被表示为:
其中,β是正则化参数,它控制允许振动引起的相位的变化程度。在这种情况下,也考虑相邻像素的振动相位。粗体符号用于指代向量。最后一项促使相邻像素的干扰引起的相移之间的差异小,这与先验知识一致,即,这些差异具有振动光栅的共同原因。
发明人提出通过利用基于多能量分箱X射线探测器的光谱信息来避免涉及机械不稳定性和变化的问题。莫尔图案的振动状态能够通过利用振动引起的莫尔图案相移的能量独立性和物体引起的位移的能量依赖性来确定。发明人提出了一种用于测量物体引起的相移的能量依赖性的体模。通过适当的数据处理,能够从物体属性中解脱出振动引起的移位,从而使图像伪影最小化。
虽然已经在附图和前面的描述中详细图示和描述了本发明,但是这样的图示和描述应被认为是图示性或示例性的而非限制性的;本发明不限于所公开的实施例。
特别地,需要指出,上述讨论涉及使用两个能量分箱,而本发明不应被理解为以这种方式受到限制,例如,在额外的能量分箱的情况下,也可以相应地提供额外的比例因子。
本领域技术人员通过研究附图、公开内容以及权利要求,在实践请求保护的发明时能够理解并实现对所公开的实施例的其他变型。
在权利要求中,“包括”一词不排除其他元件或步骤,并且词语“一”或“一个”不排除多个。
单个处理器、设备或其他单元可以实现在权利要求中记载的若干项的功能。尽管某些措施被记载在互不相同的从属权利要求中,但是这并不指示不能有利地使用这些措施的组合。
诸如提供、确定、计算、获得和处理的操作能够被实施为计算机程序的程序代码模块和/或专用硬件。
计算机程序可以被存储和/或被分布在合适的介质上,例如与其他硬件一起或作为其他硬件的部分供应的光学存储介质或固态介质,但是也可以被以其他形式分布,例如经由互联网或其他有线或无线的电信系统。
权利要求中的任何附图标记都不应被解释为对范围的限制。
Claims (14)
1.一种用于能量分箱的差分相衬X射线成像数据的处理设备(35),包括:
比例因子提供单元(50),其被布置为提供针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,以及
暗场和相移确定单元(55),其被布置为使用由所述比例因子提供单元提供的所述比例因子根据被包括在所述成像数据中的、针对所述第一能量分箱的像素的探测强度值以及针对所述第二能量分箱的所述像素的探测强度值来确定暗场信号和物体引起的相移中的至少一项。
2.根据权利要求1所述的处理设备(35),其中,所述比例因子提供单元(50)被布置为根据针对所述第一能量分箱的衰减值和针对所述第二能量分箱的衰减值来提供所述比例因子。
3.根据权利要求2所述的处理设备(35),其中,所述比例因子提供单元(50)包括查找表。
4.根据权利要求1所述的处理设备(35),其中,所述暗场和相移确定单元(55)被布置为根据所述第一能量分箱的测量数据mi,1和所述第二能量分箱的测量数据mi,2来确定针对所述第一能量分箱的衰减值A1、针对所述第二能量分箱的衰减值A2、针对所述第一能量分箱的暗场信号D1、针对所述第二能量分箱的暗场信号D2、物体引起的相移φ1,以及干扰引起的相移ψi,从而使获得所述测量数据的可能性最大化。
5.根据权利要求1所述的处理设备(35),其中,所述暗场和相移确定单元(55)被布置为采用根据下式的最小二乘最小化:
Δ(A1,A2,D1,D2,φ1,ψi)=∑(mi,1-A1Bi,1(1+D1Vi,1cos(αi,1+ψi+φ1)))2+∑(mi,2-A2Bi,2(1+D2Vi,2cos(αi,2+ψi+c·φ1)))2,
其中,Bi,1、Bi,2、Vi,1、Vi,2、αi,1和αi,2分别是在空白扫描期间获得的莫尔图案的强度、可见度和相位。
6.根据权利要求1所述的处理设备(35),其中,所述暗场和相移确定单元(55)被布置为在所述确定中包括关于促使相邻像素的干扰引起的相移之间的差异小的另外的限制。
7.一种比例因子确定设备(30),用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,所述比例因子确定设备包括:
计算单元(40),其被布置为:使用所述第一能量分箱、具有预定形状的体模(25)根据通过差分相衬X射线成像获得的相位图像来计算针对所述第一能量分箱的第一相移,并且使用所述第二能量分箱、所述体模(25)根据通过差分相衬X射线成像获得的相位图像来计算针对所述第二能量分箱的第二相移,以及
获得单元(45),其被布置为根据所述第一相移和所述第二相移来获得所述比例因子。
8.一种比例因子确定系统,用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,所述比例因子确定系统包括:
根据权利要求7所述的比例因子确定设备(30),以及
体模(25),其中,所述体模包括多个部分,每个部分具有以预定方式在横向于所述X射线的射束方向的方向上改变的高度,所述体模的所述部分的高度是所述射束方向上的部分的延伸,
其中,所述体模(25)的不同部分具有不同的平均高度,
其中,所述比例因子确定设备(30)被布置为根据针对所述第一能量分箱的衰减值和针对所述第二能量分箱的衰减值来确定所述比例因子。
9.一种差分相衬X射线成像系统(10),包括:
差分相衬成像单元(20),
X射线源(15),其被布置为提供第一能量分箱的X射线和第二能量分箱的X射线,
根据权利要求1所述的处理设备(35),其被布置为接收来自所述差分相衬成像单元(20)的能量分箱的差分相衬X射线成像数据。
10.根据权利要求9所述的成像系统(10),还包括根据权利要求8所述的比例因子确定系统,其中,所述比例因子确定系统被布置为向所述处理设备的所述比例因子提供单元(50)提供关于所述比例因子的信息。
11.一种用于能量分箱的差分相衬X射线成像数据的处理方法,包括:
比例因子提供步骤(70),其提供针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子,以及
确定步骤(75),其使用通过所述比例因子提供步骤提供的所述比例因子根据被包括在所述成像数据中的、针对所述第一能量分箱的像素的探测强度值以及针对所述第二能量分箱的所述像素的探测强度值来确定暗场信号和物体引起的相移中的至少一项。
12.根据权利要求11所述的方法,其中,所述比例因子提供步骤包括:
计算步骤(60),其使用所述第一能量分箱、具有预定形状的体模根据通过差分相衬X射线成像获得的相位图像来计算针对所述第一能量分箱的第一相移,并且使用所述第二能量分箱、所述体模根据通过差分相衬X射线成像获得的相位图像来计算针对所述第二能量分箱的第二相移,以及
获得步骤(65),其根据所述第一相移和所述第二相移来获得所述比例因子。
13.一种用于能量分箱的差分相衬X射线成像数据的处理方法的软件产品,所述软件产品包括程序代码模块,所述程序代码模块用于当所述软件产品在计算机上运行时使得所述计算机执行根据权利要求11所述的方法的步骤。
14.一种用于在能量分箱的差分相衬X射线成像的背景中确定针对第一X射线能量分箱的物体引起的相移与针对第二X射线能量分箱的物体引起的相移之间的比例因子的软件产品,所述软件产品包括程序代码模块,所述程序代码模块用于当所述软件产品在计算机上运行时使得所述计算机执行根据权利要求12所述的方法的步骤。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP17162475.2A EP3378396A1 (en) | 2017-03-23 | 2017-03-23 | X-ray imaging data processing device and method |
EP17162475.2 | 2017-03-23 | ||
PCT/EP2018/057391 WO2018172501A1 (en) | 2017-03-23 | 2018-03-23 | X-ray imaging data processing device and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109414237A true CN109414237A (zh) | 2019-03-01 |
CN109414237B CN109414237B (zh) | 2020-02-21 |
Family
ID=58412912
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880002656.3A Expired - Fee Related CN109414237B (zh) | 2017-03-23 | 2018-03-23 | X射线成像数据处理设备和方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11529112B2 (zh) |
EP (2) | EP3378396A1 (zh) |
JP (1) | JP7117316B2 (zh) |
CN (1) | CN109414237B (zh) |
WO (1) | WO2018172501A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3922178A1 (en) * | 2020-06-08 | 2021-12-15 | Koninklijke Philips N.V. | Spectral dark-field imaging |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101416073A (zh) * | 2006-03-29 | 2009-04-22 | 皇家飞利浦电子股份有限公司 | 用于重建图像的双能量衰减数据的信噪比的动态优化 |
US20140185746A1 (en) * | 2012-12-28 | 2014-07-03 | Pavlo Baturin | Spectral grating-based differential phase contrast system for medical radiographic imaging |
CN104039227A (zh) * | 2012-01-12 | 2014-09-10 | 皇家飞利浦有限公司 | 在x射线系统中生成衰减图像数据和相位图像数据 |
CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
CN105705097A (zh) * | 2013-11-08 | 2016-06-22 | 皇家飞利浦有限公司 | 针对差分相位衬度ct的经验性射束硬化校正 |
WO2016177903A1 (en) * | 2015-05-07 | 2016-11-10 | Koninklijke Philips N.V. | Beam hardening correction for scanning dark field and phase contrast imaging |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006015356B4 (de) * | 2006-02-01 | 2016-09-22 | Siemens Healthcare Gmbh | Verfahren zur Erzeugung projektiver und tomographischer Phasenkontrastaufnahmen mit einem Röntgen-System |
ATE524056T1 (de) | 2007-11-15 | 2011-09-15 | Suisse Electronique Microtech | Interferometervorrichtung und verfahren |
CN104869905B (zh) | 2012-12-21 | 2019-08-06 | 卡尔斯特里姆保健公司 | 基于微分相衬成像的医疗放射照相光栅 |
US9907524B2 (en) * | 2012-12-21 | 2018-03-06 | Carestream Health, Inc. | Material decomposition technique using x-ray phase contrast imaging system |
DE102014213817A1 (de) | 2014-07-16 | 2015-06-18 | Siemens Aktiengesellschaft | Vorrichtung und Verfahren zur Gewinnung eines Phasenkontrastbildes |
US10314556B2 (en) * | 2015-05-06 | 2019-06-11 | Koninklijke Philips N.V. | Optimal energy weighting of dark field signal in differential phase contrast X-ray imaging |
-
2017
- 2017-03-23 EP EP17162475.2A patent/EP3378396A1/en not_active Withdrawn
-
2018
- 2018-03-23 CN CN201880002656.3A patent/CN109414237B/zh not_active Expired - Fee Related
- 2018-03-23 WO PCT/EP2018/057391 patent/WO2018172501A1/en active Application Filing
- 2018-03-23 US US16/322,509 patent/US11529112B2/en active Active
- 2018-03-23 EP EP18718705.9A patent/EP3442423B1/en active Active
- 2018-03-23 JP JP2019551561A patent/JP7117316B2/ja active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101416073A (zh) * | 2006-03-29 | 2009-04-22 | 皇家飞利浦电子股份有限公司 | 用于重建图像的双能量衰减数据的信噪比的动态优化 |
CN104039227A (zh) * | 2012-01-12 | 2014-09-10 | 皇家飞利浦有限公司 | 在x射线系统中生成衰减图像数据和相位图像数据 |
US20140185746A1 (en) * | 2012-12-28 | 2014-07-03 | Pavlo Baturin | Spectral grating-based differential phase contrast system for medical radiographic imaging |
CN105705097A (zh) * | 2013-11-08 | 2016-06-22 | 皇家飞利浦有限公司 | 针对差分相位衬度ct的经验性射束硬化校正 |
CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
WO2016177903A1 (en) * | 2015-05-07 | 2016-11-10 | Koninklijke Philips N.V. | Beam hardening correction for scanning dark field and phase contrast imaging |
Also Published As
Publication number | Publication date |
---|---|
JP2020512080A (ja) | 2020-04-23 |
EP3378396A1 (en) | 2018-09-26 |
JP7117316B2 (ja) | 2022-08-12 |
US11529112B2 (en) | 2022-12-20 |
EP3442423B1 (en) | 2019-09-04 |
CN109414237B (zh) | 2020-02-21 |
EP3442423A1 (en) | 2019-02-20 |
WO2018172501A1 (en) | 2018-09-27 |
US20210321970A1 (en) | 2021-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NL193506C (nl) | Inrichting voor het automatisch corrigeren van foutieve registratie. | |
CN100496402C (zh) | 图像处理方法、图像处理系统以及x-射线ct系统 | |
CN104508704B (zh) | 人体测量方法 | |
JP4644819B2 (ja) | 微小変位計測法及び装置 | |
US20220130084A1 (en) | Systems and methods for medical image processing using deep neural network | |
US10395380B2 (en) | Image processing apparatus, image processing method, and storage medium | |
EP3915511A1 (en) | Motion blur compensation | |
CN112396560A (zh) | 用于使用深度神经网络对医疗图像进行去模糊的系统和方法 | |
Miettinen et al. | NRStitcher: non-rigid stitching of terapixel-scale volumetric images | |
WO2020146490A1 (en) | Depth reconstruction in additive fabrication | |
US10290102B2 (en) | Image quality driven non-rigid image registration | |
CN106408616A (zh) | 一种ct成像中投影图背景不一致的校正方法 | |
CN108416734A (zh) | 基于边缘驱动的文本图像超分辨率重建方法和装置 | |
JP4102785B2 (ja) | Mtf測定方法および装置 | |
JP5154751B2 (ja) | 医用画像処理装置 | |
CN109414237A (zh) | X射线成像数据处理设备和方法 | |
US10997712B2 (en) | Devices, systems, and methods for anchor-point-enabled multi-scale subfield alignment | |
JPH10320544A (ja) | 診断画像処理方法 | |
US9886749B2 (en) | Apparatus and method for parameterizing a plant | |
CN110211091A (zh) | 一种高分辨图像重构方法、装置及隐裂无损检测系统 | |
JP2016164709A (ja) | 画像処理装置、撮像装置および画像処理プログラム | |
JP2019011984A (ja) | 変位取得装置、変位取得方法、およびプログラム | |
Tabb et al. | Using cameras for precise measurement of two-dimensional plant features: CASS | |
US6683972B1 (en) | Image process, computer to evaluate data and nuclear resonance tomograph equipped with such a computer | |
CN103340644A (zh) | 一种图像位置补偿方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200221 |
|
CF01 | Termination of patent right due to non-payment of annual fee |