CN101960298A - X射线成像装置、x射线成像方法和x射线成像装置的控制方法 - Google Patents
X射线成像装置、x射线成像方法和x射线成像装置的控制方法 Download PDFInfo
- Publication number
- CN101960298A CN101960298A CN2009801081499A CN200980108149A CN101960298A CN 101960298 A CN101960298 A CN 101960298A CN 2009801081499 A CN2009801081499 A CN 2009801081499A CN 200980108149 A CN200980108149 A CN 200980108149A CN 101960298 A CN101960298 A CN 101960298A
- Authority
- CN
- China
- Prior art keywords
- ray
- subject
- imaging device
- ray imaging
- cutting element
- 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 53
- 238000000034 method Methods 0.000 title claims description 27
- 238000001514 detection method Methods 0.000 claims description 11
- 239000013078 crystal Substances 0.000 claims description 10
- 239000000463 material Substances 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 229910052737 gold Inorganic materials 0.000 claims description 3
- 229910052745 lead Inorganic materials 0.000 claims description 3
- 229910052715 tantalum Inorganic materials 0.000 claims description 3
- 229910052697 platinum Inorganic materials 0.000 claims description 2
- 229910052721 tungsten Inorganic materials 0.000 claims description 2
- 238000004846 x-ray emission Methods 0.000 claims 1
- 238000002834 transmittance Methods 0.000 abstract 2
- 238000009826 distribution Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- CNQCVBJFEGMYDW-UHFFFAOYSA-N lawrencium atom Chemical compound [Lr] CNQCVBJFEGMYDW-UHFFFAOYSA-N 0.000 description 10
- 238000002591 computed tomography Methods 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 229910002804 graphite Inorganic materials 0.000 description 1
- 239000010439 graphite Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 229910010272 inorganic material Inorganic materials 0.000 description 1
- 239000011147 inorganic material Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K1/00—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating
- G21K1/10—Scattering devices; Absorbing devices; Ionising radiation filters
-
- 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
-
- 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
- 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
-
- 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/20—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 using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
-
- 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
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K2207/00—Particular details of imaging devices or methods using ionizing electromagnetic radiation such as X-rays or gamma rays
- G21K2207/005—Methods and devices obtaining contrast from non-absorbing interaction of the radiation with matter, e.g. phase contrast
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- High Energy & Nuclear Physics (AREA)
- Pulmonology (AREA)
- Heart & Thoracic Surgery (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Toxicology (AREA)
- Biophysics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一种简化的X射线成像装置即使对于轻元素也能够以小的误差因数计算确定有效原子序数。X射线成像装置具有用于产生X射线的X射线产生单元101(400)和用于检测透过被检体104(403)的X射线的检测器105(405)。计算单元106(406)从由检测器检测的数据计算确定可归因于被检体的X射线相位的量和被检体的X射线透射率。计算单元还从由X射线相位的量确定的ρet和由X射线透射率确定的μt计算确定被检体的有效原子序数。
Description
技术领域
本发明涉及使用X射线的成像装置、成像方法和成像装置的控制方法。
背景技术
使用X射线的无损检查方法被普遍用于从工业领域到医疗领域的各种领域。X射线是具有在约0.01~(10-12~10-8m)的范围内的波长的电磁波。具有相对较短的波长()的X射线被称为硬X射线,而具有相对较长的波长()的X射线被称为软X射线。
例如,利用X射线吸收能力的差异的X射线吸收对比方法在诸如钢材料的内部裂纹的检查和手提包的检查的安全性检查的领域中获得应用以利用X射线的高透射能力。
另一方面,虽然X射线吸收对比方法对于检查低密度物体仅相对较差地执行,但是基于X射线的相位偏移的X射线相位对比成像的技术对于检查低密度物体是有效的。因此,X射线相位对比成像技术在聚合物材料的相位分离结构的成像的领域中获得应用。
同时,当在为粒子束癌症治疗制订计划时确定体内的粒子的穿透深度时,需要电子密度分布。
日本专利申请公开No.2007-082663描述了通过使用X射线计算确定用于制备电子密度分布的图像的数据和有效原子序数的方法。更具体而言,上面引用的专利文献描述了用于通过利用具有两种或更多种能量的单色X射线从吸收对比图像获取用于制备电子密度分布的图像的数据和有效原子序数的方法。
发明内容
但是,在日本专利申请公开No.2007-082663描述的方法中伴随有这样的问题:由于需要不可缺少地使用具有两个或更多个波长的X射线,因此对于该方法需要使用复杂的装置。
另外,由于在上面引用的专利文献描述的方法仅利用吸收对比图像,因此对于轻元素(light element)不能获得足够的对比度,使得它会引起误差。
鉴于上述情况,因此,本发明的一个目的是,提供可简化装置并且以小的误差因数计算确定轻元素的有效原子序数的X射线成像装置、X射线成像方法和这种X射线成像装置的控制方法。
根据本发明,通过提供一种X射线成像装置实现以上的目的,该X射线成像装置包括:用于产生X射线的X射线产生单元;用于检测从X射线产生单元发射的并透过被检体的X射线的检测器;和用于从由检测器检测到的数据计算确定可归因于被检体的X射线相位的量和被检体的X射线透射率,并进一步从由X射线相位的量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数的计算单元。
在本发明的另一方面中,提供一种X射线成像方法,该X射线成像方法包括:检测透过被检体的X射线的步骤;和从在检测步骤中检测到的数据计算确定可归因于被检体的X射线相位的量和被检体的X射线透射率,并进一步从由X射线相位的量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数的步骤。
因此,根据本发明,提供可简化装置并且即使对于轻元素仍以小的误差因数计算确定有效原子序数的X射线成像装置、X射线成像方法和这种X射线成像装置的控制方法。
参照附图阅读示例性实施例的以下说明,本发明的其它特征将变得十分清晰。
附图说明
图1是在第一实施例中描述的装置的示意性框图。
图2是被X射线照射的检测位置的示意图。
图3是在第一实施例中描述的计算单元的处理流程图。
图4是在第二实施例中描述的CT装置的示意图。
图5是在第二实施例中描述的计算单元的处理流程图。
图6是在例子中描述的装置的示意性框图。
图7是在第三实施例中描述的装置的示意性框图。
图8是示出通过旋转分光晶体(analyzer crystal)获得的X射线强度分布的示例性曲线图。
图9是在第三实施例中描述的计算单元的处理流程图。
具体实施方式
现在,参照示出可通过利用X射线吸收和X射线的相位偏移两者获得关于有效原子序数的信息的本发明的X射线成像装置和X射线成像方法的实施例的附图,更详细地描述本发明。有效原子序数是物质的平均原子序数。
第一实施例
这里描述通过利用X射线相位的量和X射线透射率获取有效原子序数分布的成像装置和图像拾取方法。
X射线对于物质的复折射率由下式(1)表达。
n=1-δ-iβ ...(1)
式(1)中的实部是与相位对应的项,虚部是与吸收对应的项。对于被检体的X射线吸收,可通过下式(2)定义线性衰减系数μ和式(1)的复折射率的虚部β:
这里,λ是X射线波长,当使用连续X射线时其是有效波长。符号t表示被检体的厚度。
当通过X射线检测器观察透射的X射线的强度时,并且,如果X射线的强度在不存在任何被检体的情况下为I0并在存在被检体的情况下为I,那么这些强度与被检体的线性衰减系数μ的关系由下式(3)表达。
上式(3)中的线性衰减系数μ可由下式(4)表达(参见Daphne F.Jackson and D.J.Hawkes,PHSICS REPORT 70,No.3(1981)169-233):
这里,ρ是被检体的密度,NA是阿伏加德罗常数,A是原子量,σel是光电吸收截面,σcoh是相干散射截面,σincoh是非相干散射截面。式(4)中的光电吸收截面σel可由下式(5)表达。
这里,ΦKN是Klein-Nishina式,b是参数。如Daphne F.Jackson and D.J.Hawkes,PHSICS REPORT 70,No.3(1981)169-233所示,参数b一般使用0.5的值。
应指出,σcoh(Z′、E′)是标准元素(standard element)的相干散射截面,Z′是标准元素的原子序数,E′是标准元素的基准能量。标准元素的基准能量E′可依据标准元素的原子序数Z′、被检体的有效原子序数Z和X射线能量E由如下示出的式(7)限定。
如Daphne F.Jackson and D.J.Hawkes,PHSICS REPORT 70,No.3(1981)169-233中所示,使用氧作为标准元素。
通过将式(5)~(7)引入式(4)来重新布置式(4),如Y.Ohno,M.Torikoshi,T.Tsunoo,K.Hyodo,Nuclear Instruments & Methods in Physics Research A 548(2005)72~77所示的那样,线性衰减系数μ可由式(8)表达:
μ=ρe(Z4F(Z,E)+G(Z,E)) ...(8)
这里,F是与光电吸收截面有关的项,G是与散射截面有关的项。
如式(8)所示,线性衰减系数可被解释为被检体的有效原子序数Z和电子密度ρe的函数。此外,从上式(8)可见,有效原子序数Z可由下式(9)表达。
另一方面,与相位对应的δ示出与电子密度ρe的相关性,如下式(10)表示:
这里,r0是经典电子半径。
如式(10)所示,可从与X射线的相位有关的信息计算确定电子密度。另外,如式(8)所示,可从与X射线的吸收有关的信息计算确定电子密度和有效原子序数两者。因此,可通过利用测量同时获得与X射线的相位有关的信息和关于X射线的吸收的信息获得μ/ρe的值。然后,可通过使用该μ/ρe执行对式(9)求解的计算过程,获得被检体的有效原子序数和电子密度作为图像。
图1是第一实施例中描述的装置的示意性框图,该框图示出其配置。在从X射线源101或X射线产生单元发射的X射线的光路上,布置X射线分割元件103、被检体104和检测器105。从X射线源101发射的X射线在空间上被X射线分割元件103分割。
从共同地拾取X射线的偏移量和偏移方向的图像的观点看,X射线分割元件103优选为具有二维布置的孔的元件,但它可替代地可具有狭缝状孔径。对于布置在X射线分割元件中的孔,假如它们具有允许X射线透过其中的形式,那么它们可以不贯穿该元件的基板。X射线分割元件103的材料选自显示高的X射线吸收率的Pt、Au、Pb、Ta和W等。
在空间上被X射线分割元件103分割的X射线照射到被检体104上,并被被检体104折射和吸收。通过检测器105检测其相位由于折射而偏移的X射线的位置和强度。通过检测器105获取的关于X射线的信息通过计算装置106被数值处理,并被输出到可典型地为监视器的显示单元107。
被检体104可以是人体,或者,如果它不是人体,那么可以是无机材料或无机/有机复合材料。可另外设置用于移动被检体104的物体移动机构(未示出)。当这种物体移动机构被设置以适当地移动被检体104时,可以获得被检体104的特定部位的图像。
检测器105典型地选自X射线CCD照相机或二维直接转换X射线检测器等。
当使用单色X射线时,可以在X射线源101和X射线分割元件103之间布置诸如单色器的单色设备或多层X射线反射镜。
图2是由检测器105检测的检测图像的示意图。应指出,检测图像的水平方向和垂直方向分别由X轴和Y轴限定。
由空心圆表示的检测位置110是在不存在被检体104的状态下获得的检测图像的检测位置。检测图像表示当计算确定偏移量时用作基准的信息。另一方面,由实心圆表示的检测位置120是在存在被检体104的状态下观察到的检测图像的检测位置。这两个图像中的每一个的检测位置可从X射线中的对应的X射线的强度分布的重心确定。
首先,通过计算装置106比较检测位置110和检测位置120,以获得沿X方向的位置偏移量(ΔX)和沿Y方向的位置偏移量(ΔY)。可通过使用位置偏移量(ΔX,ΔY)、被检体104和检测器105之间的距离L,由式(11)获得各折射的X射线的折射角度Δθx和Δθy。
可通过事先在存储器中存储在不存在任何被检体104的状态下观察到的检测图像并在各观察中参照它,计算确定位置偏移量。可通过使用被检体104的保持器的位置作为基准导出距离L。可利用通过观察其折射角度已知的标准物质并使用下式(11)来计算确定的距离L。
这里,t是X射线透过其中的被检体的厚度。换句话说,可从由X射线分割元件103分割原始X射线产生的各X射线的位置偏移确定透射的X射线的相位分布,并且,可从式(10)和(13)导出下式(14)。
现在,以下将参照图3描述由本实施例采用的有效原子序数分布图像重建方法。
首先,从透过被检体的X射线获取X射线强度数据和位置数据(S100)。
由于被检体和X射线分割元件之间的距离(L)是已知的,因此可从位置数据确定折射角度(Δθ)。另外,可从通过使用折射角度获得的微分相位对比图像确定相位然后,可通过使用式(14)获得被检体的ρet分布图像(S101)。换句话说,可从X射线相位的量确定ρet。
然后,通过使用式(3)从被检体的X射线透射率确定μt分布图像(S102)。
因此,可从ρet分布图像和μt分布图像获得μ/ρe图像。μ/ρe和有效原子序数Z之间的关系可由式(9)表达。也可通过将式(4)修改为下式(15)来表达它。
换句话说,可通过给出临时的有效原子序数Z计算确定μ/ρe的值。可通过重复执行计算以便使μ/ρe的观察值和从临时的有效原子序数Z确定的计算值匹配,来获得有效原子序数Z的收敛值(S103)。可通过使用计算的残数是否满足转换要求的确定基准终止通过计算单元的重复计算。
可通过显示单元107显示作为与有效原子序数分布有关的信息被获得的重建的有效原子序数分布图像(S104)。显示单元107可显示数值信息作为与有效原子序数分布有关的信息。另外,可以显示X射线吸收图像、X射线微分相位对比图像、X射线相位对比图像和μ/ρe图像。另外,可以设置用于控制上述的单元和设备中的每一个的控制装置。
当X射线分割元件103具有狭缝状孔径时,以上的观察的描述也是适用的,并且,可从沿与狭缝的延伸方向垂直的方向的X射线位置数据确定折射角度(Δθ)以确定相位作为替代方案,它们可被如下地计算确定:观察被检体,将X射线分割元件103在元件的平面内旋转90°以获取X方向上的数据和Y方向上的数据。
由于本实施例利用上述的X射线的相位的量,因此可提供对于轻元素仅包含低程度的误差因数的X射线成像装置、X射线成像方法和控制这种成像装置的方法。另外,由于不必使用两个波长的X射线,因此装置可在结构上被简化。
第二实施例
由于被检体的厚度t是未知的,因此,如果如式(14)所示的那样获得相位那么不能通过第一实施例直接获得电子密度分布图像。将在第二实施例中描述基于计算断层照相法(CT)的原理直接获得电子密度分布图像的技术。
图4是在本实施例中描述的CT装置的示意图,该示意图示出该CT装置的配置。X射线源201、X射线分割元件203和二维X射线检测器205被布置为通过移动机构被同步驱动以环绕被检体204移动。被X射线分割元件203空间分割的X射线照射到被检体204上,并且,透射的X射线被二维X射线检测器205检测。通过计算单元206处理通过以这种方式观察被检体204的部分而获得的数据,并且,通过显示单元207显示计算确定的电子密度分布图像。
现在,参照图5描述在本实施例中采用的电子密度分布的重建方法。
首先,从透过被检体的X射线获取X射线强度数据和位置数据(S200)。
由于被检体和X射线分割元件之间的距离(L)是已知的,因此,可从位置数据确定折射角度(Δθ)。另外,可从通过使用折射角度获得的微分相位对比图像确定相位然后,可通过使用式(14)获得被检体的ρet分布图像(S201)。
然后,通过使用式(3)从透过被检体的X射线的强度的变化确定μt分布图像(S202)。
对通过CT装置获得的所有投影数据进行以上的处理操作,以获得各个ρet分布图像和μt分布图像。随后,从这些分布图像,通过可用于计算断层照相法的图像重建方法(例如,滤波反投影法)获得ρe和μ的断层图像(S205)。
因此,虽然通过第一实施例只能确定线性衰减系数μ与电子密度ρe的比,但是,在步骤S205中可获得电子密度ρe的分布图像和线性衰减系数μ的分布图像。
然后,从电子密度ρe和线性衰减系数μ的断层图像计算确定μ/ρe图像,并且,通过与在第一实施例中使用的方法类似的方法获得有效原子序数Z的收敛值(S203)。
可通过使用计算的残数(residual)是否满足转换要求的确定基准来终止通过计算单元的重复计算。
可通过显示单元107显示关于重建的有效原子序数分布的信息(S204)。
不仅可在显示单元107上显示关于有效原子序数分布的信息,而且可在其上面显示X射线吸收图像、X射线微分相位对比图像、X射线相位对比图像、电子密度分布图像、线性衰减系数分布图像和μ/ρe图像。
因此,可以以无损方式获得被检体的电子密度的三维断层图像以及有效原子序数的三维断层图像,作为被检体的物理性质的图像。
第三实施例
在第三实施例中描述与在第一实施例使用的技术不同的获取微分相位对比图像和吸收图像以获得有效原子序数分布的技术。
图7是在本实施例中描述的装置的示意性框图,该框图示出该装置的配置。从X射线源400发射的白色X射线401通过作为单色装置操作的单色器402被单色化,并且入射到被检体403上。透过被检体403的X射线被通过使用诸如Si的单晶材料形成的分光晶体404衍射,并且通过二维检测器405被检测。
图8的曲线图中的实线示出当在不存在被检体403的状态下分光晶体404被驱动以旋转满足X射线衍射条件的角度时获得的检测器的像素的强度分布(摇摆曲线)。在图8中,横轴表示分光晶体404的旋转角度,纵轴表示X射线的强度。
另一方面,图8的曲线图中的虚线表示在存在被检体403的状态下的摇摆曲线。
冲击分光晶体404的X射线的入射角度由于被检体403的折射效果而改变。换句话说,出现X射线强度光谱的从该实线向该虚线的变化,使得可从峰值位置偏移的偏移量确定X射线的折射角度。
该实线和该虚线之间的积分强度的差通过被检体403的X射线吸收效果产生。因此,在存在被检体403的状态下通过驱动分光晶体一次旋转非常小的角度以获得其图像,可从峰值偏移量和积分强度的变化获得被检体403的微分相位对比图像和吸收图像。另外,可通过对微分相位对比图像进行积分以获得相位对比图像。
从获得的相位对比图像计算确定ρet分布图像(S301),并然后从吸收图像计算确定μt图像(S302)。然后,可从ρet分布图像和μt分布图像获得μ/ρe图像。换句话说,分别从X射线相位的量和X射线透射率确定ρet和μt。
另外,可通过给出临时的有效原子序数Z计算确定μ/ρe的值。通过重复执行计算以便使μ/ρe的观察值和从临时的有效原子序数Z确定的计算值匹配,可获得有效原子序数Z的收敛值(S303)。可通过显示单元407显示以上述的方式重建的有效原子序数分布图像(S304)。另外,可以显示X射线吸收图像、X射线微分相位对比图像、X射线相位对比图像和μ/ρe图像。
例子
图6是在例子中描述的装置的示意性框图,该示意性框图示出该装置的配置。
使用指示为X射线源301的具有Mo靶子的旋转阳极型X射线产生器作为X射线产生器。使用石墨单色器或Ni/C多层X射线反射镜作为X射线单色器302以仅提取Mo的特性线。
通过单色器302被单色化的X射线被X射线分割元件303二维分割。通过以150μm的间距通过100μm厚Pt板切割总共40×20个直径为50μm的孔,形成X射线分割元件303,但Pt可被Au、Pb、Ta、W或一些其它的类似的材料替代。
被X射线分割元件分割的X射线照射到被检体304上。通过布置在与被检体分开2m的位置上的X射线检测器305检测透过被检体304的X射线,以获得分割的X射线的位置和强度。X射线检测器305可以是X射线CCD照相机。
基于在不存在被检体的状态下观察的各X射线的位置和强度,对于各X射线的位置,通过计算单元306计算确定由式(11)限定的折射角度和由式(3)限定的μt的值。然后,通过式(12)从折射角度计算确定微分相位的量,并且,通过对微分相位的量进行积分获得相位对比图像。此外,通过式(14)获得ρet分布图像,并且,从ρet分布图像和μt分布图像获得μ/ρe分布图像。可通过给出临时的有效原子序数Z从式(9)和(15)计算确定μ/ρe的值。当获得μ/ρe的观察值和计算值时,可通过重复执行计算以便使它们匹配,来重建有效原子序数分布图像。
在作为显示单元307操作的PC监视器上显示通过计算单元306获得的有效原子序数分布图像。还可在PC监视器上显示X射线吸收图像、X射线微分相位对比图像、X射线相位对比图像和μ/ρe图像。
虽然已参照示例性实施例说明了本发明,但应理解,本发明不限于公开的示例性实施例。以下的权利要求的范围应被赋予最宽的解释以包含所有的变更方式和等同的结构和功能。
本申请要求在2008年3月12日提交的日本专利申请No.2008-062788的权益,在此以引用方式引入其全部内容。
Claims (20)
1.一种X射线成像装置,包括:
X射线产生单元,用于产生X射线;
X射线分割元件,用于空间分割从X射线产生单元发射的X射线;
检测器,用于检测由X射线分割元件分割的并透过被检体的X射线;和
计算单元,用于从由检测器检测的数据计算确定可归因于被检体的X射线偏移量和被检体的X射线透射率,并进一步从由X射线偏移量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数。
2.根据权利要求1的X射线成像装置,其中,
计算单元基于从ρet和μt计算确定的μ/ρe计算确定有效原子序数。
3.根据权利要求1的X射线成像装置,还包括:
移动机构,用于同步移动X射线产生单元、X射线分割元件和检测器。
4.根据权利要求3的X射线成像装置,其中,
计算单元通过从ρet图像重建来计算确定被检体的电子密度。
5.根据权利要求4的X射线成像装置,其中,
图像重建通过滤波反投影法被实现。
6.根据权利要求1~5中的任一项的X射线成像装置,其中,
计算单元通过使用位置偏移量以及被检体和检测器之间的距离来计算确定被检体的微分相位对比图像,所述位置偏移量是通过比较在不存在被检体的状态下观察的X射线的检测位置和透过被检体的X射线的检测位置而获得的。
7.根据权利要求6的X射线成像装置,其中,
计算单元通过对微分相位对比图像进行积分获得被检体的X射线相位对比图像。
8.根据权利要求1~7中的任一项的X射线成像装置,其中,
X射线分割元件是具有二维布置的孔的元件。
9.根据权利要求1~8中的任一项的X射线成像装置,其中,
X射线分割元件是通过使用选自Pt、Au、Pb、Ta和W的材料形成的。
10.根据权利要求1~9中的任一项的X射线成像装置,还包括:
显示单元,用于显示关于由计算单元计算确定的被检体的有效原子序数的信息。
11.根据权利要求1~10中的任一项的X射线成像装置,还包括:
被检体移动机构,用于移动被检体。
12.根据权利要求1~11中的任一项的X射线成像装置,还包括:
单色器,被布置在X射线产生单元和X射线分割元件之间以将X射线单色化。
13.一种X射线成像方法,包括:
通过X射线产生单元产生X射线的步骤;
通过X射线分割元件空间分割X射线的步骤;
通过检测器检测被X射线分割元件分割的并透过被检体的X射线的步骤;
从由检测器检测的数据计算确定可归因于被检体的X射线偏移量和被检体的X射线透射率,并进一步从由X射线偏移量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数的步骤。
14.一种控制X射线成像装置的方法,该方法包括:
控制X射线产生单元以产生X射线的步骤;
控制X射线检测器以检测被X射线分割元件分割并随后透过被检体的X射线的步骤;
控制计算单元以计算确定可归因于被检体的X射线偏移量和被检体的X射线透射率的步骤;和
从由X射线偏移量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数的步骤。
15.一种X射线成像装置,包括:
X射线产生单元,用于产生X射线;
检测器,用于检测从X射线产生单元发射的并透过被检体的X射线;和
计算单元,用于从由检测器检测的数据计算确定可归因于被检体的X射线相位的量和被检体的X射线透射率,并进一步从由X射线相位的量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线相位透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数。
16.根据权利要求15的X射线成像装置,其中,所述X射线成像装置还包含
X射线分割元件,用于空间分割从X射线产生单元发射的X射线,并且
所述X射线成像装置从通过使用X射线分割元件获得的X射线偏移量获得X射线相位的量。
17.根据权利要求15的X射线成像装置,其中,X射线成像装置还包括
单色器,用于使从X射线产生单元发射的X射线单色化,和
分光晶体,被布置在被检体和检测器之间,并且
所述X射线成像装置从通过使用分光晶体获得的X射线的折射角度获得X射线相位的量。
18.一种X射线成像方法,包括:
检测透过被检体的X射线的步骤;
从在检测步骤中检测到的数据计算确定可归因于被检体的X射线相位的量和被检体的X射线透射率的步骤;和
从由X射线相位的量的量确定的ρet(ρe:电子密度,t:被检体的厚度)和由X射线透射率确定的μt(μ:线性衰减系数)计算确定被检体的有效原子序数的步骤。
19.根据权利要求18的X射线成像方法,其中,
从通过使用用于空间分割X射线的X射线分割元件获得的X射线偏移量确定可归因于被检体的X射线相位的量。
20.根据权利要求18的X射线成像方法,其中,
可归因于被检体的X射线相位的量由通过使用分光晶体获得的X射线的折射角度确定。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008062788 | 2008-03-12 | ||
JP2008-062788 | 2008-03-12 | ||
PCT/JP2009/055215 WO2009113713A1 (en) | 2008-03-12 | 2009-03-11 | X-ray imaging apparatus, x-ray imaging method and method of controlling x-ray imaging apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101960298A true CN101960298A (zh) | 2011-01-26 |
CN101960298B CN101960298B (zh) | 2013-01-16 |
Family
ID=40639236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009801081499A Expired - Fee Related CN101960298B (zh) | 2008-03-12 | 2009-03-11 | X射线成像装置、x射线成像方法和x射线成像装置的控制方法 |
Country Status (5)
Country | Link |
---|---|
US (2) | US8214158B2 (zh) |
EP (1) | EP2257792B1 (zh) |
JP (1) | JP4512660B2 (zh) |
CN (1) | CN101960298B (zh) |
WO (1) | WO2009113713A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104427938A (zh) * | 2012-06-28 | 2015-03-18 | 西门子公司 | 用于产生相位对比显示的方法和x射线系统 |
CN105452852A (zh) * | 2013-07-30 | 2016-03-30 | 皇家飞利浦有限公司 | 通过使用相衬ct的单色衰减对比图像生成 |
CN106462983A (zh) * | 2014-05-12 | 2017-02-22 | 爱克发医疗保健公司 | 计算机化层析图像曝光和重构方法 |
TWI613436B (zh) * | 2015-07-27 | 2018-02-01 | 日立全球先端科技股份有限公司 | 缺陷判定方法、及x射線檢查裝置 |
CN108072672A (zh) * | 2017-12-14 | 2018-05-25 | 清华大学 | 一种烧蚀结构形貌及产物的在线监测装置及其监测方法 |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4512660B2 (ja) * | 2008-03-12 | 2010-07-28 | キヤノン株式会社 | X線撮像装置、x線撮像方法、x線撮像装置の制御方法 |
JP4847568B2 (ja) * | 2008-10-24 | 2011-12-28 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
JP5143147B2 (ja) * | 2009-01-15 | 2013-02-13 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
GB2479328B (en) * | 2009-01-15 | 2013-08-14 | Canon Kk | X-ray imaging apparatus and method of X-ray imaging |
JP5675169B2 (ja) | 2009-06-18 | 2015-02-25 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
JP5269041B2 (ja) * | 2009-12-04 | 2013-08-21 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
JP5213923B2 (ja) * | 2010-01-29 | 2013-06-19 | キヤノン株式会社 | X線撮像装置およびx線撮像方法 |
JP5733908B2 (ja) * | 2010-04-27 | 2015-06-10 | キヤノン株式会社 | X線撮像装置 |
JP5646906B2 (ja) * | 2010-08-06 | 2014-12-24 | キヤノン株式会社 | X線装置およびx線測定方法 |
JP2013224878A (ja) * | 2012-04-23 | 2013-10-31 | Canon Inc | X線撮像装置およびx線撮像方法 |
CN104244828B (zh) * | 2012-04-24 | 2017-06-30 | 西门子公司 | X射线设备 |
JP6168753B2 (ja) | 2012-11-12 | 2017-07-26 | キヤノン株式会社 | 演算装置、演算プログラム、x線測定システム、およびx線測定方法 |
JP6292988B2 (ja) | 2014-06-12 | 2018-03-14 | キヤノン株式会社 | 物体情報取得方法及び物体情報取得装置 |
US20180188190A1 (en) * | 2016-12-30 | 2018-07-05 | Nikon Corporation | Single-projection photon-counting x-ray quantitative phase-contrast measurements |
JP6391776B2 (ja) * | 2017-06-29 | 2018-09-19 | キヤノン株式会社 | 演算装置、演算プログラム、x線測定システム、およびx線測定方法 |
WO2020178994A1 (ja) * | 2019-03-05 | 2020-09-10 | N-Emラボラトリーズ株式会社 | 顕微鏡観察方法、透過型顕微鏡装置、および走査透過型顕微鏡装置 |
JP2023009382A (ja) * | 2021-07-07 | 2023-01-20 | 浜松ホトニクス株式会社 | 元素分布測定装置および元素分布測定方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3003175B2 (ja) * | 1990-07-26 | 2000-01-24 | 株式会社島津製作所 | X線分光による定量分析方法 |
US5651047A (en) * | 1993-01-25 | 1997-07-22 | Cardiac Mariners, Incorporated | Maneuverable and locateable catheters |
JP2006084275A (ja) * | 2004-09-15 | 2006-03-30 | Hitachi Ltd | 爆発物等の探知方法および装置 |
JP2006214952A (ja) * | 2005-02-07 | 2006-08-17 | Hitachi Ltd | 爆発物探知装置および方法 |
JP4839050B2 (ja) * | 2005-09-21 | 2011-12-14 | 独立行政法人放射線医学総合研究所 | 多色x線測定装置 |
FR2897444B1 (fr) * | 2006-02-10 | 2008-03-14 | Commissariat Energie Atomique | Procede d'estimation de rayonnement diffuse dans un detecteur bidimensionnel |
US20080080670A1 (en) * | 2006-09-29 | 2008-04-03 | Ge Security, Inc. | Systems and methods for classifying a substance |
JP4512660B2 (ja) * | 2008-03-12 | 2010-07-28 | キヤノン株式会社 | X線撮像装置、x線撮像方法、x線撮像装置の制御方法 |
-
2009
- 2009-02-18 JP JP2009035629A patent/JP4512660B2/ja active Active
- 2009-03-11 CN CN2009801081499A patent/CN101960298B/zh not_active Expired - Fee Related
- 2009-03-11 WO PCT/JP2009/055215 patent/WO2009113713A1/en active Application Filing
- 2009-03-11 US US12/526,418 patent/US8214158B2/en active Active
- 2009-03-11 EP EP09718875.9A patent/EP2257792B1/en not_active Not-in-force
-
2010
- 2010-07-12 US US12/834,355 patent/US8036336B2/en not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104427938A (zh) * | 2012-06-28 | 2015-03-18 | 西门子公司 | 用于产生相位对比显示的方法和x射线系统 |
CN105452852A (zh) * | 2013-07-30 | 2016-03-30 | 皇家飞利浦有限公司 | 通过使用相衬ct的单色衰减对比图像生成 |
CN105452852B (zh) * | 2013-07-30 | 2019-02-01 | 皇家飞利浦有限公司 | 通过使用相衬ct的单色衰减对比图像生成 |
CN106462983A (zh) * | 2014-05-12 | 2017-02-22 | 爱克发医疗保健公司 | 计算机化层析图像曝光和重构方法 |
TWI613436B (zh) * | 2015-07-27 | 2018-02-01 | 日立全球先端科技股份有限公司 | 缺陷判定方法、及x射線檢查裝置 |
CN108072672A (zh) * | 2017-12-14 | 2018-05-25 | 清华大学 | 一种烧蚀结构形貌及产物的在线监测装置及其监测方法 |
Also Published As
Publication number | Publication date |
---|---|
US8036336B2 (en) | 2011-10-11 |
WO2009113713A1 (en) | 2009-09-17 |
EP2257792A1 (en) | 2010-12-08 |
CN101960298B (zh) | 2013-01-16 |
US20100318302A1 (en) | 2010-12-16 |
US8214158B2 (en) | 2012-07-03 |
EP2257792B1 (en) | 2015-03-04 |
US20100278304A1 (en) | 2010-11-04 |
JP2009240764A (ja) | 2009-10-22 |
JP4512660B2 (ja) | 2010-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101960298B (zh) | X射线成像装置、x射线成像方法和x射线成像装置的控制方法 | |
Cantatore et al. | Introduction to computed tomography | |
RU2467355C1 (ru) | Устройство получения рентгеновских изображений и способ получения рентгеновских изображений | |
CN101405596B (zh) | 使用具有多色分布的x射线束检测对象图像的系统和方法 | |
CA2348150C (en) | Non-rotating x-ray system for three-dimensional, three-parameter imaging | |
US8477904B2 (en) | X-ray diffraction and computed tomography | |
US20100040192A1 (en) | Image reconstruction method for high-energy, dual-energy ct system | |
KR101378757B1 (ko) | 물질 원소 정보 획득 및 영상 차원의 선택이 가능한 방사선 영상화 장치 | |
Hanke et al. | X-ray microtomography for materials characterization | |
US10132764B2 (en) | System and method for reconstructing the surface topography of an object embedded within a scattering medium | |
CN1346438A (zh) | 用x射线辐射获取物体内部结构图像的方法及其实施设备 | |
CN102272860B (zh) | X射线成像装置和x射线成像方法 | |
Truong et al. | Recent developments on Compton scatter tomography: theory and numerical simulations | |
WO2009097052A1 (en) | Identification of molecular and crystalline materials by the doppler broadening of nuclear states | |
CN102272861B (zh) | X射线成像装置和x射线成像的方法 | |
JP2007508559A (ja) | 扇ビーム・コヒーレント散乱コンピュータ断層撮影法 | |
Zhu et al. | Three‐dimensional reciprocal space x‐ray coherent scattering tomography of two‐dimensional object | |
JP5030056B2 (ja) | 非破壊検査方法及び装置 | |
JP4062232B2 (ja) | X線ct装置及びx線ct装置による撮像方法 | |
Grubsky et al. | Compton imaging tomography technique for NDE of large nonuniform structures | |
Cruvinel | Using a Silicon Drift Detector to Improve an Agricultural Compton Scattering Tomograph for Soil Compaction Analysis. | |
Di Trapani et al. | Speckle-based imaging (SBI) applications with spectral photon counting detectors at the newly established OPTIMATO (OPTimal IMAging and TOmography) laboratory | |
Liew | An enhanced X-ray imaging contrast using a grating-based talbot interferometry | |
Pham | Quantitative absorption imaging using polychromatic X-ray sources | |
Nüesch | Element-sensitive X-ray imaging and computed tomography using energy-selective photodetection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130116 Termination date: 20170311 |