CN201191275Y - 一种x射线光栅相衬成像*** - Google Patents
一种x射线光栅相衬成像*** Download PDFInfo
- Publication number
- CN201191275Y CN201191275Y CN 200720190341 CN200720190341U CN201191275Y CN 201191275 Y CN201191275 Y CN 201191275Y CN 200720190341 CN200720190341 CN 200720190341 CN 200720190341 U CN200720190341 U CN 200720190341U CN 201191275 Y CN201191275 Y CN 201191275Y
- Authority
- CN
- China
- Prior art keywords
- ray
- grating
- mrow
- emitting device
- absorption
- 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.)
- Expired - Lifetime
Links
- 238000003384 imaging method Methods 0.000 claims abstract description 75
- 238000010521 absorption reaction Methods 0.000 claims abstract description 51
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000005305 interferometry Methods 0.000 claims description 10
- 239000010931 gold Substances 0.000 claims description 9
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 5
- 229910052737 gold Inorganic materials 0.000 claims description 5
- 229910001385 heavy metal Inorganic materials 0.000 claims description 5
- 239000011358 absorbing material Substances 0.000 claims description 3
- 238000004846 x-ray emission Methods 0.000 claims 1
- 230000005855 radiation Effects 0.000 description 28
- 238000000034 method Methods 0.000 description 24
- 238000005516 engineering process Methods 0.000 description 15
- 230000000694 effects Effects 0.000 description 13
- 238000013170 computed tomography imaging Methods 0.000 description 11
- 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 11
- 230000008859 change Effects 0.000 description 10
- 239000000126 substance Substances 0.000 description 10
- 230000001427 coherent effect Effects 0.000 description 9
- 230000010363 phase shift Effects 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 230000005469 synchrotron radiation Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- PLXMOAALOJOTIY-FPTXNFDTSA-N Aesculin Natural products OC[C@@H]1[C@@H](O)[C@H](O)[C@@H](O)[C@H](O)[C@H]1Oc2cc3C=CC(=O)Oc3cc2O PLXMOAALOJOTIY-FPTXNFDTSA-N 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 241000238631 Hexapoda Species 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 125000004429 atom Chemical group 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002657 fibrous material Substances 0.000 description 1
- 210000000548 hind-foot Anatomy 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000012770 industrial material Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本实用新型提供一种X射线光栅相衬成像***,用于对物体进行检测,包括:X射线发射装置,用于向被检测物体发射X射线;第一和第二吸收光栅,二者平行地依次位于X射线方向上,被检测物体折射的X射线经由该第一和/或第二吸收光栅形成强度变化的X射线信号,其中该X射线信号包含被折射X射线的折射角信息;检测单元,用于接收强度变化的X射线,并将X射线信号转换为电信号;以及数据处理单元用于从检测单元接收电信号,提取该电信号中折射角信息,并根据折射角信息利用预定算法来重构物体内部的折射率分布的断层图像。通过本实用新型,可以摆脱对射线源相干性的依赖,没有Tablot距离的限制,能使用十微米量级以上周期的光栅实现近分米量级视场的相衬成像。
Description
技术领域
本实用新型涉及图形成像领域,具体而言,涉及用X射线对物体进行CT成像的领域。
背景技术
传统的X射线成像及其CT技术是利用物质材料对X射线的衰减特性来非破坏性地检查物体的内部结构。如果物体内部的各部分结构组成的密度差异明显,则X射线成像的效果尤为显著。但是,以轻元素(例如氢、炭、氮和氧)构成的物质,它们对X射线来说是弱吸收物质,所以用传统的X射线成像技术几乎看不到它们内部的具体结构。即使是用其它辅助的手段,例如给生物组织打上造影剂也很难得到清晰的图像。二十世纪九十年代中期,由于第三代同步辐射装置的发展,硬X射线位相衬度成像(phase-contrastimaging,简称相衬成像)技术应运而生。相衬成像技术就是通过捕捉X射线的相移信息来观察物体内部的电子密度变化,从而揭示物体的内部结构。相衬成像技术使X射线成像的空间分辨率由毫米量级推进到微米量级甚至纳米量级,并将X射线成像可检测的物质范围由对X射线高吸收的重元素物质扩展到弱吸收的轻元素物质。
至今十年间,借助同步辐射射线源的高亮度和良好的相干性,相衬成像技术已经发展了至少4种的成像方法:干涉法、类同轴成像、衍射增强成像和光栅相衬成像等,这些技术与传统X射线成像及其CT技术、MRI技术、超声技术等各种现有成像技术相比具有独特的优势,因此,相衬成像技术已经成为X射线成像领域里最前沿的技术之一。然而,相衬成像技术对X射线源的苛刻要求和自身的成像特点,极大限制了它在医学上的临床应用。
首先,从X射线源方面来说,同步辐射装置的造价昂贵、体积庞大,视场小(十毫米量级),限制了它的应用范围。另外一个选择是微焦点X射线源,它发出的X射线具有部分相干性,可以实现相衬成像。但是,微焦点X射线源亮度非常低,因此探测器需要曝光时间相当长(几十秒甚至几分钟),这也是临床应用所不能容忍的。如果能在通用X射线机上实现相衬成像,那么将是非常有实际价值和意义的事情。
其次,从相衬成像方法的自身特点而言,干涉法、类同轴在实际使用中受到一定的限制。由于干涉法、类同轴成像通过观察相干X射线的干涉现象或Fresnel衍射现象来获取位相信息,都需要具有很高空间相干性的X射线源和分辨率达微米量级的探测器。但是这两种装置都很昂贵,而且微米分辨率的探测器的面积一般都很小(几平方厘米左右),这就决定了整个成像***的视场也比较小,从而做不了大样品的成像检查。衍射增强成像方法虽然可以使用通用X射线机和分辨率不那么高的探测器来实现高对比度(密度分辨率)的物体边缘增强的成像,但是光路中单色晶体的单色作用使得视场扁小,单色光亮度降低。因此,这些相衬成像方法在医学临床应用中具有一定的局限性。
瑞士PSI(Paul Scherrer Institut.)的David C等人于2002年在ESRF上用2相位光栅和1个分析晶体形成的光路上首次实现了基于光栅衍射Tabot效应的硬X射线光栅相衬成像方法。日本的Momose A等于2003年后在Spring-8上也展开了基于位相光栅和吸收光栅的光栅相衬成像方法的研究,提出用两个光栅提取一阶位相信息的技术。瑞士PSI的Weitkamp T和Pfeiffer F等于2003在DavidC工作的基础上在SLS和ESRF上同样实现了基于2个光栅的相衬成像技术。然而,上述方法都是基于同步辐射射线源上进行的,视场都很小,极大限制了相衬成像的应用。
Pfeiffer F等人于2006年取得了突破,他们使用三块功能不同的光栅实现了基于通用X射线机的光栅相衬成像方法,克服了衍射增强成像方法视场小的缺点,得到了接近64mm×64mm大小视场,从此为相衬成像技术的实际应用开拓了切实可行的道路;但是,他们所建的光栅相衬成像***的成像时间比较长(曝光时间为40秒),无法进行活体实验;另外这套***并没有实现相衬CT成像。然而,这种方法使用了源光栅(source grating)Talbot将X射线机射线源分成一系列宽度在25微米到50微米之间的、互不相干的线射线源。单个线射线源的发出的X射线是部分相干的,与相位光栅作用产生Talbot效应。在本质上,这种成像方法利用源光栅使通用X射线机发出的非相干光产生部分相干性,利用光栅衍射的Talbot效应来实现相衬成像。
因此,现有的所有的相衬成像方法(包括光栅相衬成像方法)的本质是利用相干或者部分相干的X射线的干涉或衍射现象来增强辐射图像的低对比度分辨率(也就是密度分辨率)。因为通用X射线机发出的X射线是多色、非空间相干的,根据Talbot效应原理,对于现有的光栅相衬成像技术来说,必须使用3个不同功能的光栅来实现,而且,相位光栅和吸收光栅之间的距离(被称为Talbot距离)必须选取恰当,才能满足多色条件下Talbot效应的成像要求。为了实现在硬X射线的光栅衍射,对光栅的制作提出了很高的要求:微米量级周期的精度、大深宽比。这些不利因素将极大限制光栅相衬成像技术在医学和工业上的实际应用。
发明内容
本实用新型提出一种采用相衬成像技术的相衬CT成像***,其使用两个吸收光栅就可以在多色、非空间相干的射线源(通用X射线机)下获得一阶微分相移信息。
本实用新型提供一种X射线光栅相衬成像***,用于对物体进行检测,该***包括:X射线发射装置,用于向被检测物体发射X射线;第一和第二吸收光栅,二者平行地依次位于X射线方向上,被检测物体折射的X射线经由该第一和/或第二吸收光栅形成强度变化的X射线信号,其中该X射线信号包含被折射X射线的折射角信息;检测单元,用于接收所述强度变化的X射线,并将X射线信号转换为电信号;以及数据处理单元用于从所述检测单元接收所述电信号,提取该电信号中折射角信息,并根据所述折射角信息利用预定算法来重构物体内部的折射率分布的断层图像。
本实用新型最大的优点是完全摆脱了对射线源相干性的依赖,没有Tablot距离的限制,而且能使用十微米量级以上的周期的光栅实现近分米量级视场的相衬成像。与传统X射线成像相比,本***能够对弱吸收物质(例如乳腺、血管和肌肉等软组织、纤维材料、昆虫等)进行高对比度的成像。与现有相衬成像相比,无需考虑光源的相干性条件,同时也降低了微米量级周期、大深宽比光栅的制作难度要求,并可很容易地推广到使用高能量(>40keV)X光进行相衬成像。而且,本实用新型的相衬CT成像***与传统CT成像***天然兼容,只需要在普通CT成像***上架设两光栅平动旋转***即可实现。本实用新型将进一步降低相衬成像实际应用的门槛,为相衬成像走向医学、生物学、工业材料等领域应用开拓崭新的思路和途径,具有重大的实际意义和应用价值。
附图说明
图1示出了本实用新型的光栅相衬CT成像***的示意图;
图2示出了在不同焦点尺寸的射线源照射条件下,两个光栅使用相位步进技术时所对应的模拟光强变化曲线;
图3示出了X光与物体相互作用后波阵面变化情况;
图4示出了X射线被物体折射的示意图;
图5示出了两个吸收光栅的对X射线的“光闸”作用;
图6a和图6b分别示出了CT模式的两种方式;
图7示出了塑料模型内部某像素点与背景中某象素点的光强变化曲线比较;
图8a-8d分别示出了折射角图像的三个示例;
图9示出莫尔条纹示意图;
图10示出了在射线源前面加一个多缝准直器的示意图;以及
图11示出了被测物体放置在两吸收光栅之间的示意图。
具体实施方式
见附图1,本实用新型的光栅相衬CT成像***,主要包括以下组成部分:X射线发射装置(在图中表示为射线源),用于向被检测物体发射X射线;第一和第二吸收光栅(在图中分别表示为光栅A和光栅B),二者平行地依次位于X射线方向上,被检测物体折射的X射线经由该第一和/或第二吸收光栅形成强度变化的X射线信号,其中该X射线信号包含被折射X射线的折射角信息;检测单元(在图中表示为探测器),用于接收所述强度变化的X射线,并将X射线信号转换为电信号;以及数据处理单元(未示出),用于从所述检测单元接收所述电信号,提取该电信号中折射角信息,并根据所述折射角信息利用预定算法来重构物体内部的折射率分布的断层图像。
为了更清楚的描述本实用新型,将本实用新型的技术方案与Pfeiffer F等人于2006利用三块光栅来检测物体的技术方案进行对比说明,以更突出本实用新型的特点。
在Pfeiffer F等人于2006提出的技术方案中,能够利用三块光栅而使用通用X射线机来实现光栅相衬成像方法的根本原因是源光栅将不相干的射线源细分为多个25-50微米宽的线射线源,它们发出具有部分相干性的X射线,这些X射线与4微米周期位相光栅发生光栅衍射效应,使得位相光栅在Talbot距离上产生2微米周期的自成像。对该自成像和最后一块1微米周期的吸收光栅采用相位步进(Phase stepping)技术或莫尔干涉测量法(Moireinterferometry)的相位信息提取技术,从而计算出被测物体的、与一阶折射率相关的折射角信息。要使X射线与光栅A产生光栅衍射的空间相干条件为:
其中p1为光栅A的周期,l为射线源到光栅A距离并且λ为X射线波长。假设光栅周期p1为4微米,射线源到光栅A的距离l为2米,使用12.4keV,则要求射线源的焦点尺寸至少为
这正是被源光栅细分的线射线源的宽度要求。
不同于上面基于3个光栅的相衬成像技术,本实用新型的相衬成像技术使用非相干的射线源直接照射,而不需要源光栅。也就是说,本实用新型的相衬成像技术不需要射线源的相干性(包括时间相干性和空间相干性),同样不需要考虑光栅衍射或者Talbot效应,也能得到相衬图像。当不需要考虑X射线的相干性时,成像***可用几何光学近似理论来描述。根据式(1)使得X射线与光栅A不产生衍射的条件是
从方程式(3)中可以看出,当光栅Ade周期即p1越大,则方程式右侧的值越小,因此本实用新型采用较大周期的光栅也能够得到满足要求的相衬图像。假设光栅A的周期为20微米,射线源到光栅距离2米,则
而通用X射线机的焦点尺寸至少几十微米或几百微米以上,所以通用X射线机满足***的要求。20微米周期的光栅与现有光栅相衬成像技术要求2或4微米周期的光栅相比,制作难度将大大降低。
然而,并非射线源焦点尺寸越大就越好,它受制于两个吸收光栅的相位步进技术的精度要求。图2是两吸收光栅的周期分别为20微米和22微米,源到光栅A距离2米,两光栅之间距离0.2米时,在不同焦点尺寸的射线源照射条件下,两光栅使用相位步进技术时所对应的模拟光强变化曲线。横坐标为步进大小(0.5微米/步),纵坐标为光强值。要注意的是,根据式(4)估算的结果,焦点尺寸为1微米的射线源发出的X射线是相干的,但为了方便和较大焦点的射线源进行比较,也同样在几何光学近似下的获得的。从图2发现,焦点尺寸越大,光强变化曲线就越平滑,直至焦点尺寸等于220微米时成为一条直线。也就是说相位步进技术不再起作用。当焦点尺寸大于220微米时,光强变化曲线呈现其它变化趋势。
根据相位步进技术的原理,可以推算出这射线源焦点尺寸的临界值为p0,critical:
根据图2所示,因为当焦点尺寸大于这个临界值时仍有光强变化,仍可根据这一点获得相衬效果的图像,但衬度效果没有当焦点尺寸小于p0,critical的情况要好。根据上述结果,可以得出,在p0,critical范围内,焦点尺寸越小,光强变化曲线越陡,那么提取出来的相衬图像的衬度也越高。优选情况下,射线源焦点尺寸大小p0要求不超过p0,critical的一半。当然,如果射线源焦点尺寸p0大于该临界值的话,也可以采用其他方式来解决此问题。
其中一种方法是在焦点尺寸大于p0的射线源之前放置多缝准直器,这样将大射线源细分为一系列等效的小焦点线射线源,为了提高图像衬度,每个小焦点线射线源的焦点尺寸大小p0,i也要求不超过p0,critical的一半。当然,射线源焦点大小也可以大于p0,critical,这样相衬图像的衬度由当时的光强变化曲线来决定。
自然,本领域技术人员可以理解,在采用多缝准直器的情况下,所述距离l实际上是光栅A与该多缝准直器之间的距离,而不是与实际射线源之间的距离。在射线源之前放置多缝准直器来形成具有期望焦点大小的光源的方法为本领域技术人员所公知,因此这里将不再详述。
上面讨论了对射线源的要求,下面讨论光栅之间以及光栅与射线发射装置之间的关系,需要理解的是,这里所说的射线装置可以仅包括焦点尺寸满足要求的射线源,还可以包括焦点尺寸大于预定值的射线源和多缝准直器。
因为通用X射线机的X射线束流是锥束,因此本实用新型的相衬成像***所采用的两个吸收光栅的周期要求也成几何比例关系:
因此,与原有光栅相衬成像技术相比,本实用新型的相衬成像***中两光栅之间的距离d,与X射线的波长λ无关,因此与Talbot距离无关,可根据实际情况需求由式(5)(6)而定。
下面结合本实用新型以及上述方程式来具体讨论本实用新型所使用的技术方案。
由于射线源可以是焦点尺寸较大的X射线源,因此通常的商用X射线机即可满足要求。区别于传统X射线成像,本实用新型的相衬成像技术要求射线装置的工作电压最好设置在5-40kVp条件下(具体数值由吸收光栅的金高度决定,也就是说金的高度越高,允许使用的工作电压则越高)。通用X射线机的焦点尺寸至少几十微米或几百微米以上,它发出的X射线是多色而且是非空间相干的。
优选地,本***适合使用在小焦点源上,射线源的焦点尺寸优选在50-1000微米之间。但如果使用微焦点设备(焦点小于50微米),则由于射线源的空间相干性增加,吸收光栅将产生Talbot效应。这种情况和Pfeiffer F的光栅相衬方法的原理一致,这样,另外一块吸收光栅只有放置在Talbot距离的位置点上才能有比较好的成像效果。如果使用焦点大于1000微米的射线装置,那么图像边缘模糊会增加,从而不利于图像衬度的增强。对于这种情况,可以在射线源前面放置多缝准直器(每条缝的宽度大小在50-200微米之间)如图10所示,这样将大射线源细分为一系列等效的小焦点尺寸线射线源,虽然图像的空间分辨率仍由射线源焦点尺寸决定,但各线射线源产生的衬度高的图像叠加在一起,仍能得到一个衬度相当高的相衬图像。
根据本实用新型优选实施方式,光栅部分包括两个吸收光栅,即图1中的光栅A和光栅B,用于提取X射线的一阶相位信息。吸收光栅的制作通常有两种:一般利用MEMS技术在硅片衬底上电镀上重金属材料(例如Au)线条作为光栅不透光部分,也可使用在重金属材料上刻划出光栅。光栅A和光栅B平面都平行于探测器平面(XY方向)。光栅的面积大小决定了整个相衬成像***的视场大小。举例说明,若光栅的衬底为4硅片,则***的视场在直径10cm的范围内。
根据本实用新型优选实施方式的本***的两个吸收光栅的周期在2-30微米之间。所述光栅使用重金属作为吸收材料,以金(Au)为例,金的高度由使用的X射线的能量决定,在10-100微米之间。例如,对20keV的X射线来说,金高大于16微米能阻挡90%的X射线,这样能获得衬度比较好的相衬图像。其中,光栅A的周期为p1,光栅B的周期为p2,占空比(duty cycle)一般为1。两个光栅的距离为d,如图1所示。
下面讨论所述探测器部分。
探测器***记录的是X射线强度变化。本***可使用高分辨率的X射线CCD相机或平板探测器。本***与其它相衬成像方法相比,对探测器的空间分辨率(象素大小)的要求不高:不要求十几微米、甚至几微米的高分辨率,可只使用几十微米或几百微米量级的分辨率。当然,如果探测器的空间分辨率越高,则最终相衬图像的信息越清晰。由于相衬成像的优点是能提高图像的密度分辨率(对比度分辨率),所以对探测器的动态范围有比较高的要求(>12bits)。
下面来讨论从接收到的信号中提取折射角信息(或X射线的一阶相位信息)的技术。
众所周知,X射线具有“波粒二象性”(wave-particle duality),也就是说X射线既有粒子性,又有波动性。以波动性来看,X射线同样具备透射、反射、折射、偏振、相干和不相干散射、衍射等光学性质。当X射线穿过某一物体的时候,X射线的振幅因为物体对它的吸收而减弱;同时,X射线与物质原子相干散射使X射线产生了相移(phase shift),宏观上表现为X射线的波阵面产生变形。也就是说,X射线在穿过物体时发生了折射现象,如图3所示。在宏观上,X射线与物质之间的相互作用可以用物质对X射线的复折射率n来描述:
n=1-δ-iβ (7)
其中δ与物质的相移截面p相关,β与物质的线性吸收系数μ相关,它们之间的关系如下表示
其中λ为X射线的波长,δ被称为位相因子,β被称为吸收因子。对X射线来说δ一般为10-5左右,所以n小于1,因而X射线在物质分界面上只发生微小的折射,可用折射角表示(约为微弧度量级)。假设在二维XY平面上,X射线沿笛卡儿坐标系X轴方向传播,如图4所示,被物体折射的X射线的传播偏离原来方向的折射角近似为
其中Φ(y)表示X射线沿传播路径M上的总相移。经过物体内部某几何点的X射线的相移为p(x,y)dx,则沿传播路径积分可求得总相移为
结合式(9)和(10),有
式(11)表明,折射角Δθ是物体内部折射率一阶微分的沿传播路径的积分。只要我们能够测量出折射X射线的折射角Δθ,就可根据式(11)并结合相衬CT重建方法求解出物体内部折射率n的分布。
因此,如何测量折射角Δθ是本实用新型的两光栅相衬成像技术的核心。根据几何光学近似理论,两吸收光栅的作用相当“光闸”一样。当两吸收光栅相对运动的时候,光栅对X射线来说就像闸门一样时开时闭,从而引起探测器的光强分布的强弱变化。如图5所示,从两光栅相衬成像***的二维平面俯视图观察。举例说明,在没有物体的情况下,X射线1被光栅A遮挡即被Au吸收,X射线2不被两光栅遮挡而顺利到达探测器。当在光路上放上物体后,X射线1被物体折射后的折射角为Δθ1,X射线2被物体折射后的折射角为Δθ2。这时情况发生变化,X射线1因折射而顺利避开两光栅的遮挡到达探测器,X射线2则被光栅B遮挡从而被Au吸收。这样一来,探测器所探测到的光强变化反映了X射线被物体折射的情况。通过相位步进(Phase stepping)技术和莫尔干涉测量法(Moireinterferometry)可计算出折射角信息。
为帮助更好地理解本实用新型,下面讨论上述的相位步进技术和莫尔干涉测量法。
1)相位步进(Phase stepping)技术
两光栅相对运动:光栅A固定不动,光栅B沿X方向在光栅周期p2距离范围内平移N(N>1)步(也可以光栅B固定不动,光栅A沿X方向步进)。光栅B每平移一步,***采集一次图像;在平移距离范围内采集N张图像后,可获得每个象素(探测器的探测面上的每个点)在一个光栅周期内光强变化曲线的分布情况。光强变化函数的形状类似于正弦函数:
其中A表示振幅,B与光栅周期相关,表示相位。举例说明,取塑料模型(见附图8a)内部的某个象素与背景部分的某个象素的光强变化曲线进行比较,见图7(为更好地表现差异,两光栅相对平移了一个半周期)。从图7可发现,当X射线被物体折射后,其光强变化曲线相对于背景的光强变化曲线发生了移动。光强变化曲线可用正弦函数来描述,那么光强变化曲线的移动Δf对应着正弦函数的相位变化显然,当两光栅相对运动一个周期p2,则正弦函数相位变化2π,所以有
而Δf是由于X射线被折射引起的,它与折射角的关系是:
Δf=dΔθ (14)
联合式(13)和(14),有
这样就可得到折射角信息。
举例说明,实验条件是:X射线球管焦点尺寸为400微米,光栅A和光栅B周期分别为20和22微米,源与光栅A距离l为1916.4mm,两光栅距离d为191.6mm,探测器分辨率127微米。X射线球管设置高压27kVp,电流9.6mA。两光栅相对运动的步长1微米,共步进22步,每步各采集一张投影数据。被测物体是直径1厘米的塑料模型(如图8a所示)。为尽量消除噪声的影响,每张投影图像共采集400张图像进行叠加。最后根据式(15)计算的折射角图像如图8b所示。从图8b发现,塑料模型内部的三个孔地结构的边缘清晰地表现出来了。图8c和8d显示了在同样***条件下测量出来的小老鼠的后爪和尾巴的折射角图像。
2)莫尔干涉测量法(Moire interferometry)
光栅A和光栅B平行,两者在相对旋转一微小角度ε。当X射线照射到这两块光栅后,在探测器上会观测到的莫尔条纹,其周期为
如图9所示。可见,ε越小,w越大。当X射线经物体折射后,其对应的莫尔条纹也将产生相应的形变,设其偏移量为Δf,则折射角可通过式(17)求出:
综合上述两种信息提取技术,均可以得到折射角信息。
2.CT重建算法
源光栅探测器三者不动,旋转物体;或者物体不动,源光栅探测器三者同时绕物体旋转,均可实现对物体的相衬CT成像。
对每一个投影角度φ,均采用相位步进(Phase stepping)技术或莫尔干涉测量法(Moire interferometry)方法计算当前投影角度的折射角信息Δθφ:
其中l表示X射线传播路径。
而传统CT技术的基于线性衰减系数的投影数据表达公式为
如图6所示,本实用新型的两光栅相衬CT成像***,根据物体旋转中心轴与光栅平行方向(图1中Z方向)相互关系,可将CT模式分为2种方式:物体旋转中心轴与光栅平行方向互相平行的方式(平行方式,如图6a所示);物体旋转中心轴与光栅平行方向互相垂直的方式(垂直方式,如图6b所示)。这两种方式的相衬CT重建公式是不一样的。
1)平行方式
在平行模式里,折射角信息Δθψ可看作一个在XOY平面里的矢量,其矢量方向始终垂直于投影方向。因此,Δθφ的矢量方向会随着投影方向的改变而改变,因此就不能直接使用传统衰减CT的重建算法(μ不随投影方向改变而改变)来计算物体的折射率n(或相位因子δ)断层图像。而因采用针对一阶微分投影数据的CT重建算法,需要使用希耳伯特(Hilbert)滤波器对一阶微分投影数据进行滤波,然后再反投影,重建出折射率n(或相位因子δ)断层图像。
若射线源和物体之间的距离比较远,且射线源张角较小的情况下,可使用一阶微分投影数据的平行束CT重建公式来近似重建折射率n(或相位因子δ)断层图像:
其中 式(20)是在频域极坐标系描述。-jsign(ρ)即是希尔伯特(Hilbert)滤波函数。ΘP′(ρ,β)是折射角Δθφ在频域极坐标上的表示形式。
当然也可直接采用一阶微分投影数据的扇束CT重建公式重建折射率n(或相位因子δ)断层图像,其重建公式在实数域极坐标描述为
其中,Θ(φ,s)是折射角Δθφ在实数域极坐标表示。R是源和物体的距离,D是源和探测器的距离,s是探测器上的像素点距离探测器中心的距离,U为反投影权重因子:
实际上射线源发射出来的是锥束X射线,因此可以直接采用一阶微分投影数据的锥束CT重建公式折射率n(或相位因子δ)断层图像,其重建公式用FDK型表述:
其中,
Θφ(Y,Z)表示折射角,是空间向量所指定的点,是X射线传播方向在x-y平面的投影,是x-y平面中和垂直的方向,R是射线源和物体间距离,D是射线源和探测器间的距离,Y表示探测器单元的横坐标,Z表示探测器单元的纵坐标。
2)垂直方式
需要说明的是,本实用新型的两光栅相衬CT成像***中,被检测物体也可以放置在两个吸收光栅之间,如图11所示,原理如上所述,同样也能获得相衬图像。
上述处理过程可以通过通用计算机或专用计算机来实现,其中可以对获得的投影数据进行处理,包括原始数据处理和显示、折射角信息提取、对折射角投影数据重建出物体图像以及绘制三维立体图像等等,并通过显示器显示出来。计算机可以是高性能的单个PC,也可以是工作站或机群。显示器可以是CRT传统显示器也可以是液晶显示器。
根据本实用新型,CT扫描方式可以分为两种,即源和探测器不动,物体自身旋转;或者物体不动,源和探测器绕物体旋转。为实现自动检测的功能,本实用新型的X射线光栅相衬成像***进一步包括机械控制单元,该机械控制***包括:物体台旋转装置,用于仅使被检测物体发生旋转。此时,仅有被检测的物体发生旋转,而X射线发射装置-光栅-检测单元保持不动。或者,机械控制***包括X射线发射装置-光栅-检测单元旋转装置,用于仅使发射装置-光栅-检测单元发生旋转,而此时物体保持不动,这在被检测物体不方便旋转,例如被检测物体过分庞大时尤其方便。
当然,所述物台旋转装置和X射线发射装置-光栅-检测单元旋转装置可以同时存在于该机械控制***中,或者可以仅采用其中一种装置。
已经详细描述了本实用新型,但是本领域技术人员将理解的是,对于目前所给出的公开内容,在不脱离这里所描述的本实用新型概念的精神范围内可以做出变形。因此,并不意味着本实用新型的范围局限于所示出的和所描述的特定实施例。
Claims (17)
1.一种X射线光栅相衬成像***,用于对物体进行检测,特征在于,该***包括:
X射线发射装置,用于向被检测物体发射X射线;
第一和第二吸收光栅,二者平行地依次位于X射线方向上,被检测物体折射的X射线经由该第一和/或第二吸收光栅形成强度变化的X射线信号,其中该X射线信号包含被折射X射线的折射角信息;
检测单元,用于接收所述强度变化的X射线,并将X射线信号转换为电信号;以及
数据处理单元用于从所述检测单元接收所述电信号,提取该电信号中折射角信息,并根据所述折射角信息利用预定算法来重构物体内部的折射率分布的断层图像。
2.如权利要求1的***,其中特征在于,所述X射线发射装置的焦点尺寸满足方程式 其中λ为X射线的波长,p0为X射线发射装置的焦点尺寸,p2为所述第一吸收光栅的周期,l为所述第一吸收光栅与所述X射线发射装置的距离。
3.如权利要求2的***,其中特征在于,所述X射线发射装置包括X射线源,优选地,该X射线源的焦点尺寸不大于所述X射线发射装置的临界焦点尺寸的一半,其中所述X射线发射装置的临界焦点尺寸p0,critical为 其中p2为所述第二吸收光栅的周期,d为所述第一吸收光栅与所述第二吸收光栅之间的距离。
4.如权利要求3的***,其中特征在于,所述第一和第二吸收光栅之间的距离满足方程式
5.如权利要求4的***,其中特征在于,所述X射线发射装置的焦点尺寸大小为50-1000微米。
6.如权利要求4的***,其中特征在于,所述第一吸收光栅和第二吸收光栅的周期在2微米-30微米之间。
7.如权利要求2的***,其中特征在于,所述X射线发射装置包括X射线源和多缝准直器,其中所述X射线源的焦点尺寸大于所述X射线发射装置的临界焦点尺寸的一半,所述X射线发射装置的临界焦点尺寸为 所述多缝准直器放置在所述X射线发射装置与所述第一吸收光栅之间,并且优选地,所述多缝准直器每条缝的宽度不大于所述X射线发射装置的临界焦点尺寸的一半。
8.如权利要求7的***,其中特征在于,所述多缝准直器每个缝的大小为50-200微米。
9.如权利要求2的***,其中特征在于,所述第一和第二吸收光栅采用重金属作为吸收材料,所述吸收材料的厚度为10微米-100微米。
10.如权利要求9的***,其中特征在于,所述重金属为金。
11.如权利要求1-10中任一项的***,其中特征在于,所述被检测物***于所述X射线发射装置与所述第一和第二吸收光栅之间。
12.如权利要求1-10中任一项的***,其中特征在于,所述第一吸收光栅和第二吸收光栅分别位于被检测物体的两侧。
13.如权利要求1的***,其中特征在于,所述检测单元的动态范围大于12位。
14.如权利要求1的***,其中特征在于,提取所述折射角信息包括使用相位步进法和莫尔干涉测量法。
15.如权利要求1的***,其中特征在于,重构物体内部的折射率分布的断层图像是通过CT重构算法来实现的,包括平行束,扇束和锥束形式的一阶微分投影数据重建算法。
16.如权利要求15的***,其中特征在于,所述CT重构算法包括平行方式和垂直方式。
17.如权利要求1的***,其中特征在于,所述***进一步包括机械控制单元,该机械控制***包括:物体台旋转装置,用于仅使被检测物体发生旋转;和/或X射线发射装置-光栅-检测单元旋转装置,用于仅使发射装置-光栅-检测单元发生旋转。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200720190341 CN201191275Y (zh) | 2007-11-23 | 2007-11-23 | 一种x射线光栅相衬成像*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200720190341 CN201191275Y (zh) | 2007-11-23 | 2007-11-23 | 一种x射线光栅相衬成像*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN201191275Y true CN201191275Y (zh) | 2009-02-04 |
Family
ID=40335438
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200720190341 Expired - Lifetime CN201191275Y (zh) | 2007-11-23 | 2007-11-23 | 一种x射线光栅相衬成像*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN201191275Y (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011003278A1 (zh) * | 2009-07-07 | 2011-01-13 | 清华大学 | X射线暗场成像***和方法 |
WO2011130896A1 (zh) * | 2010-04-19 | 2011-10-27 | 清华大学 | X射线源光栅步进成像***与成像方法 |
CN102365687A (zh) * | 2009-03-27 | 2012-02-29 | 皇家飞利浦电子股份有限公司 | 消色差的相衬成像 |
CN103356208A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 用于人体医学检测的二维成像***及方法 |
CN104095645A (zh) * | 2013-04-03 | 2014-10-15 | 中国科学院高能物理研究所 | 高分辨率和高信噪比的x射线成像***和成像方法 |
CN104095646A (zh) * | 2013-04-03 | 2014-10-15 | 中国科学院高能物理研究所 | 高分辨率和高信噪比的双旋ct成像***和成像方法 |
CN104321805A (zh) * | 2012-05-14 | 2015-01-28 | 皇家飞利浦有限公司 | 暗场计算机断层摄影成像 |
WO2016070771A1 (zh) * | 2014-11-04 | 2016-05-12 | 清华大学 | X射线相衬成像***与成像方法 |
CN105935297A (zh) * | 2016-06-23 | 2016-09-14 | 中国科学院深圳先进技术研究院 | X射线光栅相衬成像ct*** |
CN109580667A (zh) * | 2018-12-05 | 2019-04-05 | 中国科学技术大学 | 单光栅相衬成像方法及*** |
CN111474190A (zh) * | 2020-05-21 | 2020-07-31 | 中国科学院上海光学精密机械研究所 | 二氧化硅微球三维全形貌检测装置及其检测方法 |
CN111721786A (zh) * | 2019-03-22 | 2020-09-29 | 中国科学院深圳先进技术研究院 | 一种x射线干涉仪及成像*** |
WO2020199194A1 (zh) * | 2019-04-04 | 2020-10-08 | 中国科学技术大学 | 一种x射线相位衬度成像方法 |
US11266363B2 (en) * | 2016-02-17 | 2022-03-08 | Rensselaer Polytechnic Institute | Energy-sensitive multi-contrast cost-effective CT system |
CN114264227A (zh) * | 2021-11-26 | 2022-04-01 | 武汉联影生命科学仪器有限公司 | 焦点尺寸和位置的测量装置、方法和测量模组 |
-
2007
- 2007-11-23 CN CN 200720190341 patent/CN201191275Y/zh not_active Expired - Lifetime
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102365687A (zh) * | 2009-03-27 | 2012-02-29 | 皇家飞利浦电子股份有限公司 | 消色差的相衬成像 |
CN102365687B (zh) * | 2009-03-27 | 2015-08-19 | 皇家飞利浦电子股份有限公司 | 消色差的相衬成像 |
WO2011003278A1 (zh) * | 2009-07-07 | 2011-01-13 | 清华大学 | X射线暗场成像***和方法 |
CN101943668B (zh) * | 2009-07-07 | 2013-03-27 | 清华大学 | X射线暗场成像***和方法 |
US8515002B2 (en) | 2009-07-07 | 2013-08-20 | Tsinghua University | X-ray dark-field imaging system and method |
WO2011130896A1 (zh) * | 2010-04-19 | 2011-10-27 | 清华大学 | X射线源光栅步进成像***与成像方法 |
US9134259B2 (en) | 2010-04-19 | 2015-09-15 | Tsinghua University | X-ray source grating stepping imaging system and image method |
CN103356223A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 用于人体医学检测的 ct 成像***及方法 |
CN103356208A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 用于人体医学检测的二维成像***及方法 |
CN103364418A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 光栅剪切二维成像***及光栅剪切二维成像方法 |
CN103364418B (zh) * | 2012-04-01 | 2016-08-03 | 中国科学院高能物理研究所 | 光栅剪切二维成像***及光栅剪切二维成像方法 |
CN103356208B (zh) * | 2012-04-01 | 2015-01-14 | 中国科学院高能物理研究所 | 用于人体医学检测的二维成像***及方法 |
CN103356207B (zh) * | 2012-04-01 | 2016-08-03 | 中国科学院高能物理研究所 | 基于光栅剪切成像的医学检测设备和方法 |
CN103356207A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 基于光栅剪切成像的医学检测设备和方法 |
CN104321805A (zh) * | 2012-05-14 | 2015-01-28 | 皇家飞利浦有限公司 | 暗场计算机断层摄影成像 |
CN104321805B (zh) * | 2012-05-14 | 2017-12-15 | 皇家飞利浦有限公司 | 暗场计算机断层摄影成像 |
CN104095646A (zh) * | 2013-04-03 | 2014-10-15 | 中国科学院高能物理研究所 | 高分辨率和高信噪比的双旋ct成像***和成像方法 |
CN104095645A (zh) * | 2013-04-03 | 2014-10-15 | 中国科学院高能物理研究所 | 高分辨率和高信噪比的x射线成像***和成像方法 |
US10267752B2 (en) | 2014-11-04 | 2019-04-23 | Tsinghua University | X-ray phase-contrast imaging system and imaging method |
WO2016070771A1 (zh) * | 2014-11-04 | 2016-05-12 | 清华大学 | X射线相衬成像***与成像方法 |
US11266363B2 (en) * | 2016-02-17 | 2022-03-08 | Rensselaer Polytechnic Institute | Energy-sensitive multi-contrast cost-effective CT system |
CN105935297A (zh) * | 2016-06-23 | 2016-09-14 | 中国科学院深圳先进技术研究院 | X射线光栅相衬成像ct*** |
CN109580667A (zh) * | 2018-12-05 | 2019-04-05 | 中国科学技术大学 | 单光栅相衬成像方法及*** |
CN111721786A (zh) * | 2019-03-22 | 2020-09-29 | 中国科学院深圳先进技术研究院 | 一种x射线干涉仪及成像*** |
WO2020199194A1 (zh) * | 2019-04-04 | 2020-10-08 | 中国科学技术大学 | 一种x射线相位衬度成像方法 |
CN111474190A (zh) * | 2020-05-21 | 2020-07-31 | 中国科学院上海光学精密机械研究所 | 二氧化硅微球三维全形貌检测装置及其检测方法 |
CN114264227A (zh) * | 2021-11-26 | 2022-04-01 | 武汉联影生命科学仪器有限公司 | 焦点尺寸和位置的测量装置、方法和测量模组 |
CN114264227B (zh) * | 2021-11-26 | 2023-07-25 | 武汉联影生命科学仪器有限公司 | 焦点尺寸和位置的测量装置和方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN201191275Y (zh) | 一种x射线光栅相衬成像*** | |
CN101576515B (zh) | X射线光栅相衬成像***及方法 | |
CN103365067B (zh) | 可实现三维动态观测的光栅剪切成像装置和方法 | |
US8972191B2 (en) | Low dose single step grating based X-ray phase contrast imaging | |
CN102221565B (zh) | X射线源光栅步进成像***与成像方法 | |
US8515002B2 (en) | X-ray dark-field imaging system and method | |
EP2586373B1 (en) | X-ray interferometer | |
JP5601909B2 (ja) | X線撮像装置及びこれを用いるx線撮像方法 | |
US20050286680A1 (en) | X-ray imaging system and imaging method | |
US20120099702A1 (en) | Correction method for differential phase contrast imaging | |
WO2007125833A1 (ja) | X線撮像装置及びx線撮像方法 | |
US20010038680A1 (en) | X-ray phase-contrast medical micro-imaging methods | |
JP5665834B2 (ja) | X線撮像装置 | |
Seifert et al. | Measurement and simulative proof concerning the visibility loss in x-ray Talbot-Lau Moiré imaging | |
Sung et al. | Wave optics simulation of grating‐based X‐ray phase‐contrast imaging using 4D Mouse Whole Body (MOBY) phantom | |
KR101827726B1 (ko) | 멀티플 및 엑스선회절격자 선택형 위상차 엑스선 영상장치 | |
Alloo et al. | Multimodal Intrinsic Speckle-Tracking (MIST) to extract rapidly-varying diffuse X-ray scatter | |
WO2018168621A1 (ja) | 放射線画像生成装置 | |
Pyakurel | Phase and dark field radiography and CT with mesh-based structured illumination and polycapillary optics | |
JP6797762B2 (ja) | 放射線画像生成装置及び放射線画像生成方法 | |
Lider | Talbot and Talbot–Lau X-ray interferometers | |
Włodarczyk et al. | Analytical reconstructions of intensity modulated x-ray phase-contrast imaging of human scale phantoms | |
KR101678596B1 (ko) | 엑스선회절격자 선택형 위상차 엑스선 영상장치 | |
Gullberg et al. | Tensor tomography of dark field scatter using x-ray interferometry with bi-prisms | |
KR101755153B1 (ko) | 3차원 영상의 생성 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CX01 | Expiry of patent term |
Granted publication date: 20090204 |