CN107957429A - 用于生成x射线图像数据的方法 - Google Patents
用于生成x射线图像数据的方法 Download PDFInfo
- Publication number
- CN107957429A CN107957429A CN201710960677.3A CN201710960677A CN107957429A CN 107957429 A CN107957429 A CN 107957429A CN 201710960677 A CN201710960677 A CN 201710960677A CN 107957429 A CN107957429 A CN 107957429A
- Authority
- CN
- China
- Prior art keywords
- ray
- raw data
- data set
- generate
- correction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000005855 radiation Effects 0.000 claims abstract description 65
- 238000012360 testing method Methods 0.000 claims abstract description 20
- 238000012937 correction Methods 0.000 claims description 43
- 230000006870 function Effects 0.000 claims description 34
- 238000005259 measurement Methods 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 17
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 230000007547 defect Effects 0.000 description 6
- 238000007689 inspection Methods 0.000 description 6
- 239000000047 product Substances 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000006096 absorbing agent Substances 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 229920003023 plastic Polymers 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 238000000333 X-ray scattering Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000013067 intermediate product Substances 0.000 description 1
- 238000012432 intermediate storage Methods 0.000 description 1
- 238000009607 mammography Methods 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
- 238000012805 post-processing Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000010408 sweeping 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/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
-
- 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]
-
- 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]
-
- 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/585—Calibration of detector units
-
- 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]
-
- 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/03—Investigating materials by wave or particle radiation by transmission
-
- 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/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- 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
-
- 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/419—Imaging computed tomograph
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pulmonology (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本公开涉及用于生成X射线图像数据的方法。X射线图像数据通过包括X射线源和X射线检测器的X射线系统生成，X射线检测器具有多个有规律排列的像素，其中：指定具有多个参数的数学散射辐射模型；通过X射线系统扫描测试对象，生成第一和第二原始数据集；在生成第二原始数据集时，X射线掩模被插入在X射线源与测试对象之间，其具有对X射线透明的区域和对X射线不透明的区域；基于第一和第二原始数据集确定用于校准散射辐射模型的参数值；通过X射线系统扫描检查对象，以生成第三原始数据集；通过经校准的散射辐射模型处理第三原始数据集，以生成经校正的第三原始数据集；并基于经校正的第三原始数据集，从检查对象生成一组X射线图像数据。
Description
技术领域
本发明涉及一种用于通过一个X射线系统来生成X射线图像数据的方法，X射线系统包括一个X射线源和一个X射线检测器，该X射线检测器具有多个有规律排列的像素。本发明还涉及一种被配置为执行对应方法的一种X射线系统和一种数据处理单元。
背景技术
通常，用于生成X射线图像数据的方法(即使用X射线的成像方法)利用以下情形：在穿过一个检查对象时，X射线的强度通过吸收而被衰减，这种吸收一方面由取决于检查对象的材料结构，另一方面取决于通过检查对象的路径长度。
利用这种方法，在与材料相互作用时还会发生X射线的附加散射，这种附加散射通常也导致生成的这些图像的质量显著降低。为此，期望减少这些散射的X射线(也称为散射辐射)对这种成像方法的影响。
一种已知的解决方案是使用所谓的准直仪或散射辐射栅格，借助于这些准直仪或散射辐射栅格，大部分散射辐射在穿过检查对象之后、并在任何可能的技术测量之前被吸收，并且因此从辐射场中被过滤掉。由此带来的缺点在于：这些准直仪或散射辐射栅格一方面有时占据大量空间，另一方面还会吸收一些不再可用于技术测量的非散射X射线。
另一种方法规定：通过距检查对象足够远处(“气隙”)的一个检测器来执行X射线的技术测量，并且由于散射X射线的主要部分在与非散射X射线不同的传播方向上被传播，从而确保散射X射线的主要部分能够绕开检测器。然而，在许多情况下，检测器与检查对象之间的距离无法按照所期望的那样来选择。
上述方法的进一步描述例如可以在以下出版物中找到：
[1]ROBERT G.GOULD硕士和JOHN HALE博士在American Journal ofRoentgenology(1974年9月，第122卷，第1期)上发表的“CONTROL OF SCATTERED RADIATIONBY AIR GAP TECHNIQUES:APPLICATIONS TO CHEST RADIOGRAPHY”。
[2]Krol A1,Bassano SINCE,Chamberlain CC,Prasad SC在Med Phys.(1996年7月；23(7):1263-70)上发表的“Scatter reduction in mammography with air gap”。
[3]弗里德里希-亚历山大大学(FAU)的医学物理学学院的Matthias Baer等人在iCT 2012上发表的“Scatter Correction Methods in Dimensional CT”。
公布文本DE 10 2011 006 400 A1和DE 100 55 739 B4也涉及散射辐射(即X射线的散射)的主题。在这两个案子中，使用了确定散射辐射分量的方法。
发明内容
基于上文所述，本发明的目的是公开一种用于生成X射线图像数据的有利方法，以及被配置为执行对应方法的一种X射线系统和/或一种数据处理单元。
根据本发明，该方法被实现为：具有权利要求1的多个特征的一种方法、具有权利要求16的多个特征的一种X射线系统以及具有权利要求17的多个特征的一种数据处理单元。多个优选的变形被包含在多个从属权利要求中。关于该方法所描述的多个优点和多个优选实施例也可以被类似地转移到X射线系统和/或数据处理单元，反之亦然。
在本上下文中，一种相应的方法被用于生成X射线图像数据，并且优选地与多种工业应用相适应。因此，该方法特别地被用于检验产品或中间产品的质量缺陷。在本上下文中，该方法有利地被实施以便按照一种尽可能简单的方式、并在相对短的时间内，通过X射线来检查多个非常相似的对象。
在本上下文中，X射线图像数据通过一个X射线系统被生成，例如通过一个计算机断层扫描仪被生成，X射线系统包括一个X射线源和一个X射线检测器，并且通常包括一个数据处理单元。在此，优选地，X射线检测器被实施为一个所谓的间接变换X射线检测器，并且优选地进一步包括多个有规律排列的像素，这些像素通常被排列在多行和多列中。
该方法还被用于生成相对高质量的X射线图像数据，在X射线图像数据中，来自所谓的散射辐射(即来自在背景技术中所述的散射X射线)的不利并且因此不希望的影响得以减少。为此，该方法公开了一个具有多个参数的数学散射辐射模型，其中所述模型基于X射线图像数据的生成，并且为此被保持驻留在数据处理单元中(至少当数据处理单元存在时)。
在本上下文中，数学散射辐射模型通常被指定一次，例如，在设置X射线系统和/或上述数据处理单元时被指定，因此散射辐射模型例如是上述X射线系统的操作系统的一部分，优选地，散射辐射模型仅在维护工作的上下文中被调整或更新。另一方面，在正常操作中，通常仅散射辐射模型参数的参数值被调整，其中对于要通过X射线系统来检查的对象而言，优选地，这种调整并非针对对象的每次检查来执行，而是通常在校准过程的上下文中被执行，在校准过程之后，通常会检查多个相似的对象，而不对这些参数值进行任何调整。
然而，在各种情况下，根据此处所描述的方法，首先指定一个数学散射辐射模型。然后，在一个校准过程的上下文中，通过X射线系统来扫描一个测试对象从而生成第一原始数据集，并且测试对象被又一次扫描从而生成第二原始数据集，其中在X射线系统中，为了生成第二原始数据集，一个X射线掩模被插入在X射线源与测试对象之间，X射线掩模具有：对X射线透明的至少一个区域，以及对X射线不透明的至少一个区域。
在此，关于是首先执行没有X射线掩模的扫描、还是首先执行具有X射线掩模的扫描，在原理上，对于该方法本身而言都是可行的，因此根据一个备选变形，该方法可以逆序执行。
在本上下文中，X射线掩模通常具有一个板状的形状，并且优选地，以与X射线检测器平行的方式被对准。在本上下文中，X射线掩模被用于阻挡X射线在某些方向上的传播，并且以这种方式来遮蔽X射线检测器上的某些区域，由于从X射线源到X射线检测器上的这些受遮蔽区域的直接路径受到阻挡，因此可以假设：仅散射X射线能够碰撞X射线检测器的这些受遮蔽区域。通过使用这种X射线掩模，能够使用已知原理来生成关于散射辐射的统计数据，其中所述统计数据是测试对象的特征，并且这种特征散射辐射将由数学散射辐射模型来描述。
针对散射辐射模型的参数，一方面，基于第一原始数据集、另一方面基于第二原始数据集，来进一步确定这些参数的参数值，其中这些参数值被用于校准散射辐射模型，并且因此适应于作为测试对象特征的统计数据。在此，根据一个变形，对这些参数值的确定借助于数据处理单元的评估工具而被半自动地执行，或者被数据处理单元完全自动地执行。在下文中，经校准的散射辐射模型被用于检查至少一个检查对象(即，待检查对象)，尽管检查对象与测试对象通常不同，但是检查对象通常与测试对象相似，其中通过X射线系统来扫描一个对应的检查对象，以生成一个第三个原始数据集。在本上下文中，为了生成第三原始数据集，X射线掩模从光束路径中被再次移除。备选地，检查对象也被用作测试对象。
然后，第三原始数据集通过经校准的散射辐射模型而被处理，以生成经校正的第三原始数据集，最后，基于经校正的第三原始数据集，一组X射线图像数据从检查对象而被生成。
如上文所述，有利地，以此方式校准的散射辐射模型不仅基于针对一个检查对象的检查，而且还基于针对不同检查对象(这些检查对象优选为彼此相似)的多次检查或扫描。因此，例如，通过X射线系统，可以渐进地检查一个批次中的多个已完成部件或组件的质量缺陷，其中出于以下原因而将一组相同参数值用于这些检查中的每个检查：对应部件或组件彼此非常相似，以使得可以假设在每种情况下均出现非常相似的散射辐射分布，并且可以假设散射辐射分布被经校准的散射辐射模型非常有效地描述，其中经校准的散射辐射模型具有所确定的一组参数值。在本上下文中，对于每个其他检查对象，另外的原始数据集(即第四原始数据集、第五原始数据集等)通过X射线系统而被生成，并且每个另外的原始数据集通过先前校准的散射辐射模型而被处理，以生成一个相关联的校正原始数据集。最后，基于每个校正原始数据集，来从对应的检查对象生成一组相关联的X射线图像数据。特别地，在这种情况下，优选地，多个检查对象中的一个检查对象(即，例如来自该批次的多个已完成部件中的一个部件，或者多个已完成组件中的一个组件)也被用作一个测试对象。
进一步优选地，数学散射辐射模型由校正函数表征，利用校正函数，来校正为获得一个原始数据集、而以X射线检测器的具有坐标u和v的一个像素测量的强度I(u,v)，其中这样校正的强度IK(u,v)被用于生成对应的校正原始数据集。在本上下文中，对应的校正原始数据集并非必须立即被用于生成X射线图像数据，而是例如可以首先被存储在一个永久性存储器中。然后，根据需要(有时在很晚的时间点)来生成相关联的一组X射线图像数据。
根据一个有利的变形，这一校正函数包括一个校正项KS(u,v)，该校正项KS(u,v)描述经测量的强度I(u,v)在校正期间要被去除的一个分量，该分量由来自X射线检测器外部的、碰撞具有坐标u和v的像素的散射辐射引起。
在一个有利变形中，校正项KS(u,v)被指定为一个函数f(P1(GP2*I)(u，v))，其中该函数的自变量被指定为：2D高斯函数GP2(u,v)与测量的强度I(u,v)的卷积、再与第一参数P1的乘积，并且其中第二参数P2对应于2D高斯函数的标准差。
进一步优选地，这里，该函数f被指定为一个简单的线性函数，即具有参数P3和P4的这种类型的函数f(p1(GP2*I)(u，v))＝P3(p1(GP2*I)(u，V))+P4。
如下这样的校正函数也是有利的，该校正函数包括校正项KU(u，v)，该校正项KU(u，v)描述测量的强度I(u，v)在校正期间要被去除的一个分量，其中上述分量由X射线检测器中的环境(特别是由多个相邻像素)对具有坐标u和v的像素的影响而引起。因此，该校正项KU(u，v)描述所谓的串扰，这种串扰特别地出现在具有闪烁体的间接变换X射线检测器中。
在本上下文中，优选地，对应校正项KU(u，v)＝P1(GP2*I)(u，v)还从第一参数P1与一个卷积的乘积而被获得，该卷积是一个归一化的2D高斯函数GP2(u，v)与测量的强度I(u，v)的卷积，其中第二参数P2对应于该2D高斯函数的标准差。
此外，优选地，数学散射辐射模型以及校正函数被保持为尽可能简单，并且因此根据一个优选实施例，校正函数恰好包括上述两个校正项KS(u，v)和KU(u，v)。然后，由此得到：
IK(u，v)＝I(u，v)-KS(u，v)-KU(u，v)＝I(u，v)-P3(p1(GP2*I)(u，v))+P4-P1(GP2*I)(u，v)。
如上文所述，在对数学散射辐射模型校准的上下文中，优选地，针对散射辐射模型的多个参数的多个参数值(即，每次校准有一组参数值)被确定。在本上下文中，优选地，将一组参数值中的、要被确定的多个参数值的数目保持为尽可能少，即，例如小于10，并且优选地小于6。如果散射辐射模型通过上述具有两个校正项KS(u，v)和KU(u，v)的校正函数而获得，则例如需要确定四个参数值，即P1、P2、P3和P4。
为了确定P1，则例如相对于强度IoRM而绘制强度ImRM，其中强度ImRM在利用X射线掩模来生成第二原始数据集的期间被测量，而强度IoRM在未利用X射线掩模的、第一原始数据集的生成期间而被测量。在这种上下文中，仅考虑所有被遮蔽的像素中的如下像素：即，当X射线掩模被施加时，能够被来自X射线检测器外部的散射X射线碰撞的所有像素。这导致一种类型的点云，其中每个点一方面针对在一次扫描期间用于获得X射线图像的每个方向，另一方面针对所考虑的每个像素。然后，P1被确定为：最大强度区域中的ImRM/IoRM的平均最小值，或者对ImRM/IoRM的鲁棒估计。
通过在不同点扩展函数上进行平均，例如可以进一步确定针对参数P2的参数值，其中不同点扩展函数可以被用于描述：一个被遮蔽区域与一个非遮蔽区域之间的过渡区域在X射线检测器上的强度分布。进一步优选地，还根据以下原理来执行上述确定：首先，将ImRM中的被遮蔽区域内的投影图像IoRM设置为0，从而创建一个图像I2。然后，I2与点扩散函数进行卷积，从而将强度转移到先前被设置为0的区域。最后，寻找点扩散函数，对这个点扩散函数而言，经卷积的I2与ImRM之间的最可能一致性出现在这些被遮蔽区域内。这种搜索采取一种简单坐标搜索的形式。
对参数P3和P4的确定也例如基于一个示图来执行，在该示图中，相对于适配强度IA而绘制在具有施加X射线掩模的情况下得到的测量强度ImRM，其中适配强度IA对应于：在没有施加X射线掩模的情况下得到的测量强度IoRM减去校正项KU。这再次产生了一个点云，其中为了确定参数P3和P4，优选地，一条直线被插入，该直线与点云的多个外端相交，以作为多个插值点。然后，对应直线的斜率以及该直线与轴IA的交点产生了针对参数P3或P4的值。
此外，取决于变形和要求规范，校准数学散射辐射模型所需的X射线掩模通常包括：对X射线区域不透明的多个区域，特别是在一个表面上有规律排列、优选地有规律分布的多个区域。此处，通常，就对X射线不透明的对应区域的数目和/或对应区域的布置密度而言，这种数目和/或密度会随着检查对象(即，待检查的对象)的几何形状的复杂性而增加，这些检查对象要通过X射线系统并基于在此所描述的方法被检查。
在本上下文中，X射线掩模本身例如由对X射线透明的材料制成，并且为了形成对X射线不透明的区域，X射线掩模包括多个凹部，由对X射线不透明的材料制成的多个元件被插入在这些凹部中。
在本上下文中，为了本申请的目的，术语“不透明”和“透明”应被理解为意味着一个透明区域仅以较低程度来吸收X射线，也即，例如由塑料或具有相对低密度的另一材料制成，并且非透明区域以一个更大的程度来吸收X射线，即例如由金属或具有较高密度的另一材料制成。
然而，在本上下文中，也可以根据基本原理以相反方式实现X射线掩模，由此X射线掩模可以由对X射线不透明的材料制成、并且包括对X射线透明的多个区域(特别是有规律排列的多个区域)，或者X射线掩模可以被简单地实施为一种穿孔类型的掩模。然后，数学散射辐射模型可以被相应地适配。
此外，优选地，多个有规律排列的区域被实施为条纹状的、矩形的或椭圆形的(例如，也可以是圆形的)，在矩形实施例的情况下，进一步优选地，这些区域以棋盘图案的形式来排列。
X射线掩模的以下实施例也是有利的，其中X射线掩模包括数目少于200(优选地少于100，特别地少于50)的有规律排列的区域，优选地，这些有规律排列的区域还被布置为均匀分布在一个面上。
用以生成原始数据集的方式通常还要与相应的应用相适应，特别是与多个待检查对象的几何形状的复杂性相适应。在大多数情况下，在本上下文中，为了生成每个原始数据集，即在扫描过程期间，像通常的3D成像一样，多个X射线图像从不同方向被获得。
优选地，在本上下文中，为了生成每个原始数据集，从不同方向获得的X射线图像少于100个(进一步优选地少于60个，并且特别地少于40个)。然而，在特别简单的情况下，例如当由于质量缺陷而要检查轴向对称体(即，投影在不同方向上基本不发生改变的实体)或板状体时，来自一个方向的单个记录也是足够的。
附图说明
下面参考示意图来更详细地解释本发明的示例性实施例。在附图中：
图1示出了一个X射线系统的框图，
图2示出了一种方法的多个方法步骤的框图，
图3示出了一个第一强度分布的示图，以及
图4示出了一个第二强度分布的示图。
所有附图中的对应部分被给予相同的附图标记。
具体实施方式
通过下文的示例并在图1中描绘了一个X射线系统2，X射线系统2包括一个X射线源4、一个具有多个像素8的X射线检测器6和一个数据处理单元10，并且X射线系统2被实施为用于多种工业应用的一个计算机断层扫描仪。在此，X射线检测器6被实施为一个间接变换X射线检测器6，间接变换X射线检测器6的多个像素8以覆盖较大面积的方式被排列成多个行和多个列。X射线系统2还包括一个X射线掩模12，X射线掩模12可以根据需要被移动到X射线源4的锥形光束中，或从锥形光束中移出。
在本上下文中，如果被定位于X射线源4的锥形光束中，则X射线掩模12一方面位于X射线检测器6的上游，并且另一方面位于针对多个对象14(测试对象和检查对象两者)的容器的上游，并且被用于遮蔽X射线检测器6上的多个区域。为此，X射线掩模12具有板状形状，并由对X射线透明的基本材料16(例如，塑料)制成，并且X射线掩模12包括有规律排列的多个容器，多个吸收元件18被插入在这些容器中，这些吸收元件18由一种对X射线不透明的材料(例如金属)制成。在本上下文中，多个吸收元件18在基本材料16中有规律地排列，这种排列形成一种类型的棋盘图案。
在该实施例中，X射线系统2可以被用于生成X射线图像数据，并且因此生成多个待检查对象(即多个检查对象)的多个X射线图像，以便例如检验它们的质量缺陷。在本上下文中，X射线系统2允许生成X射线图像数据，利用X射线图像数据，所谓的散射辐射(即散射的X射线)对图像质量的不利和不希望的影响得以减少。通过使用一个数学散射辐射模型，使得能够减少这种不希望的影响，数学散射辐射模型被保持驻留在数据处理单元10中，并且在从多个待检查对象生成X射线图像数据期间被考虑。
优选地，数学散射辐射模型是相对简单的模型，这种模型在在以下情况下被指定为大致是一次性的，这些情况例如是：在数据处理单元10的配置期间，或者在对X射线系统2和/或数据处理单元10维护的上下文中；并且这种模型随后在X射线系统2的操作期间被使用。在本上下文中，对应的数学散射辐射模型包括被用于校准数学散射辐射模型的多个参数，并且因此这些参数可以容易地与不同的待检查对象相适应。
在本上下文中，在工业应用的情况下，由于此时通常存在大量非常相似的多个对象要检查，所以校准通常较不频繁地被执行。在这种情况下，一次性校准被执行，并且随后多个对象14通过X射线系统2被检查，其中相同校准被用作每次检查的基础，即：在每次从一个对应的待检查对象生成X射线图像数据时，相同校准被用作每次检查的基础。例如，如果要通过X射线系统2来检验一个批次的产品的质量缺陷，则针对这一批次来执行一次性校准，并且确定针对数学散射辐射模型的多个参数的一组参数值，并且随后该组参数值和利用该组参数值被校准的数学散射辐射模型被用作对该批次中产品的所有检查的基础。
在本上下文中，为了校准数学散射辐射模型，首先，一个测试对象被定位在X射线系统2中，其中测试对象是多个待检查对象中的一个待检查对象，或者是与待检查对象非常相似的对象。随后，执行两次扫描，并且因此生成两个原始数据集，其中在这两次扫描中的一次扫描中，即在生成一个原始数据集期间，X射线掩模12被使用，并且为此被移动到X射线源4的锥形光束中。
然后，这两个原始数据集被评估，并且这两个原始数据集被用作确定一组参数值的基础。然后，这些参数值被用于校准散射辐射模型，此后X射线图像数据从多个待检查对象被生成，以便检查这些对象的质量缺陷。在本上下文中，这些待检查对象一个接一个地被定位在X射线系统2中，并且通过X射线系统2被扫描，其中针对每个待检查对象14，生成一组原始数据(即原始数据集)。然后，来自这些待检查对象的对应原始数据集借助于先前经校准的数学散射辐射模型被各自处理，并且因此来自这些待检查对象的每个原始数据集被用作生成校正原始数据集的基础。最后，校正原始数据集被用作生成X射线图像数据的基础，以此方式，散射辐射的影响得以减少。
在本上下文中，取决于应用、并且特别地取决于待检查对象的几何形状的复杂性，对应扫描包含来自不同方向的一个或多个X射线图像。
此外，图2以框图的形式描绘了通过X射线系统2实现的方法。该方法包括其基本组成部分A、B、C和D，即：指定一个数学散射辐射模型A，借助于一个测试对象来校准数学散射辐射模型B，从至少一个检查对象生成X射线图像数据C，以及可选的图像后处理或图像评估D。
在此，校准B至少包括：从一个测试对象生成第一原始数据集B1；从测试对象生成第二原始数据集B2，其中为了生成第二原始数据集B2，X射线掩模12被插入在X射线源与测试对象之间；通过第一原始数据集和第二原始数据集的联合评估，来校准数学散射辐射模型B3；以及对经校准的数学散射辐射模型或至少对多个参数值进行中间存储B4，该多个参数值在对该数学散射辐射模型的多个参数进行校准的上下文中被确定。
图像数据的生成C依次至少包括：从一个检查对象生成一个第三原始数据集C1；借助于经校准的散射辐射模型来处理第三原始数据集，并且生成经校正的第三原始数据集C2；以及基于经校正的第三原始数据集，从检查对象生成一组X射线图像数据C3。
如上文所提到的，在校准数学散射辐射模型的上下文中，针对散射辐射模型的多个参数而确定多个参数值，即对于每个校准具有一组参数值。在示例性实施例中，数学散射辐射模型通过一个校正函数被获得，该校正函数包括四个参数P1、P2、P3和P4。
然后，为了确定P1，在具有X射线掩模12的情况下的测量强度ImRM相对于在没有X射线掩模12的情况下的测量强度IoRM被绘制。这考虑了所有被遮蔽的像素8，即：在具有施加X射线掩模12的情况下、仅能够由来自X射线检测器6外部的散射X射线碰撞的所有像素8。这产生一种类型的点云，其中一个点一方面针对已从其获得X射线图像的每个方向，另一方面针对被考虑的每个像素8。相关联的示图被示出在图3中。然后，P1被确定为ImRM/IoRM的一个平均最小值。在图3中，虚线具有斜率P1，并且在没有一个施加X射线掩模12的情况下，在最大强度区域中与点云相交。
通过在不同头部扩展函数上进行平均，例如可以进一步确定针对参数P2的参数值，其中不同头部扩展函数可以被用于描述：一个被遮蔽区域与一个非遮蔽区域之间的过渡区域中在X射线检测器6上的强度分布。
参数P3和P4的确定也例如基于图4所示的示图被执行。此处，在具有施加X射线掩模12的情况下测量强度ImRM相对于一个适配强度IA被绘制，其中适配强度IA对应于在没有施加X射线掩模12的情况下的测量强度IoRM减去校正项KU(u，v，)＝P1(GP2*I)(u，v)，，其中u和v描述多个像素8的位置。这再次产生了一个点云，其中在该示例性实施例中，为了确定参数P3和P4，一条直线被插入，该直线与点云的多个外端相交。对应直线在图4中被示出为虚线。然后，对应直线的斜率及其与轴IA的交点产生参数P3和P4的值。
本发明不限于上述示例性实施例。相反，本领域技术人员也可以在不脱离本发明的主题的情况下从其中得出其他变形。此外，特别地结合该示例性实施例描述的所有单独特征也可以在不脱离本发明的主题的情况下以另一种方式相互组合。
Claims (15)
1.一种用于通过一个X射线系统(2)来生成X射线图像数据的方法，所述X射线系统(2)包括一个X射线源(4)和一个X射线检测器(6)，所述X射线检测器(6)具有多个有规律排列的像素(8)，其中
-一个具有多个参数的数学散射辐射模型被指定(A)，
-一个测试对象(14)通过所述X射线系统被扫描，以生成一个第一原始数据集(B1)，
-所述测试对象(14)通过所述X射线系统被扫描，以生成一个第二原始数据集，其中为了生成所述第二原始数据集，一个X射线掩模(12)被插入在所述X射线源(4)与所述测试对象(14)之间，所述X射线掩模(12)具有对X射线透明的至少一个区域(16)，并且具有对X射线不透明的至少一个区域(18)(B2)，
-对于所述散射辐射模型的所述多个参数的多个参数值的确定一方面是基于所述第一原始数据集，而另一方面是基于所述第二原始数据集，所述散射辐射模型利用所述多个参数值而被校准(B3)，
-一个检查对象(14)通过所述X射线系统(2)被扫描，以生成一个第三原始数据集(C1)，
-所述第三原始数据集通过经校准的所述散射辐射模型被处理，以生成一个经校正的第三原始数据集(C2)，以及
-基于所述经校正的第三原始数据集，一组X射线图像数据从所述检查对象(14)被生成(C3)。
2.根据权利要求1所述的方法，其中
-至少一个另外的检查对象通过所述X射线系统(2)被扫描，以生成一个另外的原始数据集(C1)，
-每个另外的原始数据集通过经校准的散射辐射模型被处理，以生成一个相关联的经校正的原始数据集(C2)，以及
-基于每个经校正的原始数据集，相关联的一组X射线图像数据从对应的所述检查对象(14)被生成(C3)。
3.根据权利要求1或2所述的方法，
其中所述数学散射辐射模型由一个校正函数表征，利用所述校正函数，来校正为获得一个原始数据集、而利用所述X射线检测器(6)的具有坐标u和v的一个像素(8)测量的强度I(u,v)，这样校正的强度IK(u,v)被用于生成对应的校正原始数据集。
4.根据权利要求3所述的方法，
其中所述校正函数包括一个校正项KS(u,v)，所述校正项KS(u,v)描述经测量的所述强度I(u,v)在所述校正期间要被去除的一个分量，所述分量由来自所述X射线检测器(6)外部的、碰撞具有所述坐标u和v的所述像素(8)的散射辐射引起。
5.根据权利要求4所述的方法，
其中所述校正项Ks(u，v)＝f(P1(GP2*I)(u，v))被指定为一个函数f，其中所述函数的自变量被指定为：2D高斯函数GP2(u,v)与经测量的所述强度I(u,v)的卷积再与第一参数P1的乘积，并且其中第二参数P2对应于所述2D高斯函数的标准差。
6.根据权利要求5所述的方法，
其中所述函数f是一个线性函数。
7.根据权利要求3至6中的任一项所述的方法，
其中所述校正函数包括一个校正项KU(u,v)，所述校正项KU(u,v)描述经测量的所述强度I(u,v)的、在所述校正期间要被去除的一个分量，其中所述分量由所述X射线检测器(6)中的环境对具有所述坐标u和v的所述像素(8)的影响引起。
8.根据权利要求7所述的方法，
其中获取2D高斯函数GP2(u,v)与经测量的所述强度I(u,v)的卷积与第一参数P1相乘的乘积作为所述校正项KU(u，v)＝P1(GP2*I)(u，v)，其中第二参数P2对应于所述2D高斯函数的所述标准差。
9.根据权利要求1至8中的任一项所述的方法，
其中所述X射线掩模(12)包括对X射线不透明的多个有规律排列的区域(18)。
10.根据权利要求9所述的方法，
其中所述多个有规律排列的区域(18)被实施为条纹状的、矩形的或椭圆形的。
11.根据权利要求9或10所述的方法，
其中所述X射线掩模(12)包括数目少于200、特别是数目少于100的有规律排列的区域(18)。
12.根据权利要求1至11中的任一项所述的方法，
其中为了生成每个原始数据集，多个X射线图像从不同方向被获取。
13.根据权利要求12所述的方法，
其中为了生成每个原始数据集，数目少于100、特别是少于60的X射线图像从不同方向被获取。
14.一种X射线系统(2)，包括一个X射线源(4)、一个X射线检测器(6)以及一个数据处理单元(10)，所述X射线检测器(6)具有多个有规律排列的像素(8)，所述数据处理单元(10)被配置为执行如前述权利要求中的任一项所述的方法。
15.一种数据处理单元(10)，被配置为执行如前述权利要求中的任一项所述的方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE102016220096.8A DE102016220096B3 (de) | 2016-10-14 | 2016-10-14 | Verfahren zur Generierung von Röntgenbilddaten |
| DE102016220096.8 | 2016-10-14 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107957429A true CN107957429A (zh) | 2018-04-24 |
| CN107957429B CN107957429B (zh) | 2019-12-13 |
Family
ID=60996745
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710960677.3A Active CN107957429B (zh) | 2016-10-14 | 2017-10-16 | 用于生成x射线图像数据的方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US10012602B2 (zh) |
| CN (1) | CN107957429B (zh) |
| DE (1) | DE102016220096B3 (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112788995A (zh) * | 2018-10-04 | 2021-05-11 | 皇家飞利浦有限公司 | 自适应防散射设备 |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6145889B2 (ja) * | 2014-03-24 | 2017-06-14 | 富士フイルム株式会社 | 放射線画像処理装置および方法並びにプログラム |
| DE102018201411A1 (de) * | 2018-01-30 | 2019-08-01 | Robert Bosch Gmbh | Verfahren zum Ermitteln eines zeitlichen Verlaufs einer Messgröße, Prognosesystem, Aktorsteuerungssystem, Verfahren zum Trainieren des Aktorsteuerungssystems,Trainingssystem, Computerprogramm und maschinenlesbares Speichermedium |
| CN108606805B (zh) * | 2018-03-15 | 2022-02-11 | 东软医疗系统股份有限公司 | 一种散射校正方法、装置及设备 |
| US20230281888A1 (en) * | 2022-03-02 | 2023-09-07 | General Electric Company | Computed tomography scatter and crosstalk correction |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1050790A (zh) * | 1989-10-02 | 1991-04-17 | 老代尔夫特光学工业有限公司 | 制造图象增强管校正掩模的方法 |
| US5157700A (en) * | 1988-09-02 | 1992-10-20 | Canon Kabushiki Kaisha | Exposure apparatus for controlling intensity of exposure radiation |
| JPH11219891A (ja) * | 1998-02-04 | 1999-08-10 | Hitachi Ltd | マスクの検査方法および装置 |
| CN101810488A (zh) * | 2009-01-29 | 2010-08-25 | 皇家飞利浦电子股份有限公司 | 探测值校正设备 |
| CN102727234A (zh) * | 2011-03-31 | 2012-10-17 | 西门子公司 | 产生检查对象的图像数据的方法、处理装置和x射线系统 |
| CN102985008A (zh) * | 2010-05-26 | 2013-03-20 | 株式会社岛津制作所 | X射线摄影装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10047720A1 (de) * | 2000-09-27 | 2002-04-11 | Philips Corp Intellectual Pty | Vorrichtung und Verfahren zur Erzeugung eines Röntgen-Computertomogramms mit einer Streustrahlungskorrektur |
| DE10055739B4 (de) * | 2000-11-10 | 2006-04-27 | Siemens Ag | Streustrahlungskorrekturverfahren für eine Röntgen-Computertomographieeinrichtung |
| US7092482B2 (en) * | 2003-04-11 | 2006-08-15 | Fischer Imaging Corporation | Signal profiling for medical imaging systems |
| DE102011006400B4 (de) * | 2011-03-30 | 2020-12-10 | Siemens Healthcare Gmbh | Verfahren zum Ermitteln des Anteils an Streustrahlung in 2D-Röntgenbildern |
-
2016
- 2016-10-14 DE DE102016220096.8A patent/DE102016220096B3/de active Active
-
2017
- 2017-10-16 CN CN201710960677.3A patent/CN107957429B/zh active Active
- 2017-10-16 US US15/784,308 patent/US10012602B2/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5157700A (en) * | 1988-09-02 | 1992-10-20 | Canon Kabushiki Kaisha | Exposure apparatus for controlling intensity of exposure radiation |
| CN1050790A (zh) * | 1989-10-02 | 1991-04-17 | 老代尔夫特光学工业有限公司 | 制造图象增强管校正掩模的方法 |
| JPH11219891A (ja) * | 1998-02-04 | 1999-08-10 | Hitachi Ltd | マスクの検査方法および装置 |
| CN101810488A (zh) * | 2009-01-29 | 2010-08-25 | 皇家飞利浦电子股份有限公司 | 探测值校正设备 |
| CN102985008A (zh) * | 2010-05-26 | 2013-03-20 | 株式会社岛津制作所 | X射线摄影装置 |
| CN102727234A (zh) * | 2011-03-31 | 2012-10-17 | 西门子公司 | 产生检查对象的图像数据的方法、处理装置和x射线系统 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112788995A (zh) * | 2018-10-04 | 2021-05-11 | 皇家飞利浦有限公司 | 自适应防散射设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20180106734A1 (en) | 2018-04-19 |
| US10012602B2 (en) | 2018-07-03 |
| DE102016220096B3 (de) | 2018-02-08 |
| CN107957429B (zh) | 2019-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kruth et al. | Computed tomography for dimensional metrology | |
| CN107957429A (zh) | 用于生成x射线图像数据的方法 | |
| US8150131B2 (en) | Computer tomography method for determining an object region with the aid of projection images | |
| EP2263075B1 (de) | Verfahren zum tomographischen vermessen von mechanischen werkstücken | |
| US10282872B2 (en) | Noise reduction in tomograms | |
| US9375192B2 (en) | Reconstruction of a cone beam scanned object | |
| DE102007016370A1 (de) | Verfahren und eine Messanordnung zum Erzeugen von dreidimensionalen Bildern von Messobjekten mittels invasiver Strahlung | |
| JP6680470B2 (ja) | 画像取得装置及び画像取得方法並びに画像補正プログラム | |
| KR20190060054A (ko) | 치과용 콘-빔 ｃｔ의 산란선 보정방법 및 보정장치 | |
| DE102018105709A1 (de) | Verfahren zur computertomografischen Messungen von Werkstücken | |
| JPWO2019083014A1 (ja) | 光子計数型のｘ線検出データを処理する方法及び装置、並びに、ｘ線装置 | |
| US8824759B2 (en) | Correcting axial tilt based on object positions in axial slices of three dimensional image | |
| Ingacheva et al. | Polychromatic CT Data Improvement with One‐Parameter Power Correction | |
| JP2008541110A (ja) | 被検対象物の材料界面を決定するための方法及び装置 | |
| WO2003062856A1 (de) | Vorrichtung und verfahren zum erstellen einer korrekturkennlinie für eine reduktion von artefakten bei einer tomographie | |
| Reiter et al. | Case study of empirical beam hardening correction methods for dimensional X-ray computed tomography using a dedicated multi-material reference standard | |
| CN113167746A (zh) | 用于无损分析测试对象的动态辐射准直 | |
| EP3258251B1 (en) | Image reconstruction method for x-ray measurement device, method for manufacturing structure, image reconstruction program for x-ray measurement device, and x-ray measurement device | |
| CN108065950B (zh) | 一种放射成像方法及其系统 | |
| CN112204607B (zh) | 用于x射线成像的散射校正 | |
| Flay et al. | Effects of off-focal radiation on dimensional measurements in industrial cone-beam micro-focus X-ray computed tomography systems | |
| DE102022103888A1 (de) | Verfahren und Vorrichtung zur Computertomografiemessung | |
| DE102014008671A1 (de) | Verfahren zur optimalen Anordnung eines Objekts in einer Vorrichtung und Vorrichtung zur Darstellung einer inneren räumlichen Struktur des Objekts | |
| TWI613998B (zh) | 斷層合成影像邊緣假影抑制方法 | |
| JP6281640B2 (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 |