CN104132953A - 一种双能x射线相位衬度成像装置及其实现方法 - Google Patents
一种双能x射线相位衬度成像装置及其实现方法 Download PDFInfo
- Publication number
- CN104132953A CN104132953A CN201410384492.9A CN201410384492A CN104132953A CN 104132953 A CN104132953 A CN 104132953A CN 201410384492 A CN201410384492 A CN 201410384492A CN 104132953 A CN104132953 A CN 104132953A
- Authority
- CN
- China
- Prior art keywords
- grating
- ray
- energy
- image
- detector
- 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 92
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000003287 optical effect Effects 0.000 claims abstract description 61
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 238000010521 absorption reaction Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 16
- 230000000007 visual effect Effects 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 230000001915 proofreading effect Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 13
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 9
- 239000004926 polymethyl methacrylate Substances 0.000 description 9
- 210000000988 bone and bone Anatomy 0.000 description 5
- 230000010363 phase shift Effects 0.000 description 5
- 239000011358 absorbing material Substances 0.000 description 4
- 238000012113 quantitative test Methods 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 239000003814 drug Substances 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 239000011521 glass Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000037182 bone density Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000000326 densiometry Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- PVBALTLWZVEAIO-UHFFFAOYSA-N diodone Chemical compound OC(=O)CN1C=C(I)C(=O)C(I)=C1 PVBALTLWZVEAIO-UHFFFAOYSA-N 0.000 description 1
- 229960001845 diodone Drugs 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000011707 mineral Substances 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
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种双能X射线相位衬度成像装置及其实现方法,该双能X射线相位衬度成像装置沿光路依次包括X光机、源光栅、分束光栅、样品室、分析光栅和X射线探测器,其中:X光机用于发出X射线;源光栅用于将大焦点的X光源分成为若干个不相干的小焦点光源;分束光栅用于将小焦点光源分成若干束,照射到样品室中的样品上,并在在分析光栅上形成几何投影;样品室用于放置并固定样品,同时带动样品进行旋转;分析光栅用于与分束光栅一起在X射线探测器上形成莫尔条纹;X射线探测器用于获取并记录该莫尔条纹。利用本发明进行双能X射线相位衬度成像时,选取的两个能量值V高和V低可以根据实际情况随意调整,从而扩大双能X射线相位衬度成像的使用范围。
Description
技术领域
本发明涉及X射线成像技术领域,尤其是一种双能X射线相位衬度成像装置及其实现方法。
背景技术
基于传统吸收衬度的双能X射线成像技术,在临床医学中获得了广泛的应用。双能X射线骨密度测量仪,是目前骨质疏松诊断最准确和最有价值的方法[1,2]。双能X射线CT成像技术,能够把骨组织或软组织从X射线图像中除去,得到仅含软组织或骨组织的图像,从而可以便捷医学诊断。双能减影X射线血管造影技术(利用碘造影剂)可以去除骨骼、肌肉和其他组织,只留下单纯血管的清晰图像,因而可以提供更加准确的诊断信息,同时可以保证更为安全的介入手术[3,4]。但是对于由C、H、N、O等轻元素组成的弱吸收物质,由于他们在不同能量下,吸收系数的差异性很小,所以基于传统吸收衬度的双能X射线成像技术,在弱吸收生物组织间的分辨和识别方面,发挥作用有限[5]。
与传统吸收衬度成像相比,对于弱吸收物体,硬X射线相位衬度成像技术能够提供大大增强的图像衬度和明显改善的成像质量[6-9],这是因为在硬X射线波段,弱吸收物质的相移因子是吸收因子的1000多倍[10]。同时,弱吸收物质的相移因子在不同能量下的差异性明显高于吸收因子。因此双能X射线相位衬度成像技术在弱吸收物质间的分辨识别方面有着巨大的应用价值[5]。
在目前众多的X射线相位衬度成像技术中[11-16],基于三块光栅的成像技术[16],是唯一可以利用普通X光源进行成像试验的,因此该技术最有可能在临床医学、社会安全检查和工业无损检测等领域获得推广应用。
目前的双能X射线相位衬度成像装置和试验方法,正是基于三块光栅的相位衬度成像技术[5]。在该双能成像技术中,相位信息提取是基于衍射光栅的Talbot效应,选取的两个X光机管电压V高和V低必须满足以下两个条件:(1)X光机在高电压V高下发出射线的等效能量E高通过相位光栅时产生的相移为(2)X光机在低电压V低下发出射线的等效能量E低通过相位光栅时产生的相移为由于入射射线通过相位光栅时产生的相移Δφ和射线能量E的一一对应性因此,对于一套构建完毕的成像***,进行双能成像时选取的两个管电压就固定不可改变,在目前的机器配置中,选用的两个管电压必须为V高=70KV和V低=40KV,对应射线的平均能量为E高=40KeV和E低=20KeV,不能选取其他的管电压进行双能X射线相位衬度成像[5]。
然而,在实际的双能X射线相位衬度成像中,由于试验对象各种各样,为了得到比较理想的成像效果,需要根据实际情况,选用合理的X射线能量[17,18]。因此,目前双能X射线相位衬度成像技术的使用范围受到了限制。
参考文献:
[1]R.B.Mazess,H.S.Barden,J.P.Bisek,and J.Hanson,″Dual-energy x-ray absorptiometry for total-body and regional bone-mineraland soft-tissue composition,″The American journal of clinical nutrition,vol.51,pp.1106-1112,1990.
[2]H.A.Kotaniemi,P.Vainio,and E.Alhava,″Bonedensitometry of the spine and femur in children by dual-energy x-rayabsorptiometry,″Bone and mineral,vol.17,pp.75-85,1992.
[3]G.Di Chiro,R.A.Brooks,R.M.Kessler,G.S.Johnston,A.E.Jones,J.R.Herdt,et al.,″Tissue Signatures with Dual-Energy ComputedTomography 1,″Radiology,vol.131,pp.521-523,1979.
[4]T.R.Johnson,B.Krauss,M.Sedlmair,M.Grasruck,H.Bruder,D.Morhard,et al.,″Material differentiation by dual energy CT:initialexperience,″European radiology,vol.17,pp.1510-1517,2007.
[5]C.Kottler,V.Revol,R.Kaufmann,and C.Urban,″Dual energyphase contrast x-ray imaging with Talbot-Lau interferometer,″Journal ofApplied Physics,vol.108,p.114906,2010.
[6]G.Schulz,T.Weitkamp,I.Zanette,F.Pfeiffer,F.Beckmann,C.David,et al.,″High-resolution tomographic imaging of a human cerebellum:comparison of absorption and grating-based phase contrast,″Journal of TheRoyal Society Interface,vol.7,pp.1665-1676,2010.
[7]M.Bech,A.Tapfer,A.Velroyen,A.Yaroshenko,B.Pauwels,J.Hostens,et al.,″In-Vivo dark-field and phase-contrast x-ray imaging,″Scientific reports,vol.3,2013.
[8]A.Momose,W.Yashiro,K.Kido,J.Kiyohara,C.Makifuchi,T.Ito,et al.,″X-ray phase imaging:from synchrotron to hospital,″PhilosophicalTransactions of the Royal SocietyA:Mathematical,Physical and EngineeringSciences,vol.372,p.20130023,2014.
[9]A.Tapfer,M.Bech,I.Zabette,P.Symvoulidis,S.Stangl,G.Multhoff,et al.,″Three-dimensional imaging of whole mouse models:comparing nondestructive X-ray phase-contrast micro-CT withcryotome-based planar epi-illumination imaging,″Journalofm icroscopy,vol.253,pp.24-30,2014.
[10]A.Momose,″Recent advances in X-ray phase imaging,″JapaneseJournal of Applied Physics Part l-Regular Papers Brief Communications&Review Papers,vol.44,pp.6355-6367,Sep 2005.
[11]A.Momose,T.Takeda,Y.Itai,and K.Hirano,″Phase-contrastX-ray computed tomography for observing biological soft tissues,″Naturemedicine,vol.2,pp.473-475,1996.
[12]S.Wilkins,T.Gureyev,D.Gao,A.Pogany,and A.Stevenson,″Phase-contrast imaging using polychromatic hard X-rays,″Nature,vol.384,pp.335-338,1996.
[13]T.Davis,D.Gao,T.Gureyev,A.Stevenson,and S.Wilkins,″Phase-contrast imaging of weakly absorbing materials using hard X-rays,″Nature,vol.373,pp.595-598,1995.
[14]C.David,B.H.H.Solak,and E.Ziegler,″Differentialx-ray phase contrast imaging using a shearing interferometer,″Applied physicsletters,vol.81,pp.3287-3289,2002.
[15]A.Momose,S.Kawamoto,I.Koyama,Y.Hamaishi,K.Takai,andY.Suzuki,″Demonstration of X-ray Talbot interferometry,″Japanese journalofapplied physics,vol.42,p.L866,2003.
[16]F.Pfeiffer,T.Weitkamp,O.Bunk,and C.David,″Phase retrievaland differential phase-contrast imaging with low-brilliance X-ray sources,″Nature Physics,vol.2,pp.258-261,Apr 2006.
[17]L.Grodzins,″Optimum energies for x-ray transmission tomographyof small samples:Applications of synchrotron radiation to computerizedtomography I,″Nuclear Instruments and Methods in Physics Research,vol.206,pp.541-545,1983.
[18]L.Grodzins,″Critical absorption tomography of small samples:proposed applications of synchrotron radiation to computerized tomographyII,″Nuclear Instruments and Methods in Physics Research,vol.206,pp.547-552,1983.
发明内容
(一)要解决的技术问题
有鉴于此,本发明的主要目的在于提供一种双能X射线相位衬度成像装置及其实现方法,以构建一套灵活的双能X射线相位衬度成像装置,在利用该装置进行双能X射线相位衬度成像时,选取的两个能量值V高和V低可以根据实际情况随意调整,从而扩大双能X射线相位衬度成像的使用范围。
(二)技术方案
为达到上述目的,本发明提供了一种双能X射线相位衬度成像装置,该装置沿光路依次包括X光机1、源光栅2、分束光栅3、样品室4、分析光栅5和X射线探测器6,其中:X光机1,用于发出X射线;源光栅2,用于将大焦点的X光源分成为若干个不相干的小焦点光源;分束光栅3,用于将小焦点光源分成若干束,照射到样品室4中的样品上,并在在分析光栅5上形成几何投影;样品室4,用于放置并固定样品,同时带动样品进行旋转;分析光栅5,用于与分束光栅3一起在X射线探测器上形成莫尔条纹;X射线探测器6,用于获取并记录该莫尔条纹。
为达到上述目的,本发明还提供了一种实现双能X射线相位衬度成像的方法,该方法包括:
对准X光机的出光点和X射探测器接收平面的中心点;
对X射线探测器进行Offset矫正和Gain矫正;
精确对准源光栅、分束光栅与分析光栅,使这三块光栅的刻线相互平行、三块光栅所在的平面和X射线探测器的平面相互平行,同时使三块光栅的中心点、X光机的出光点和X射线探测器中心点在一条直线上;
调节X光机的管电压到高能V高处,进行相位步进扫描,同时旋转样品室,采集图像并重建出高能V高处样品的三维相位信息;
调节X光机的管电压到低能V低处,进行相位步进扫描,同时旋转样品室,采集图像并重建出低能V低处样品的三维相位信息。
(三)有益效果
从上述技术方案可以看出,利用本发明提供的双能X射线相位衬度成像装置及其实现方法,进行双能X射线相位衬度成像时,主要具有以下几个优点:
(1)由于该双能X射线相位衬度成像装置的信息提取是利用光栅几何投影产生的莫尔条纹,因此在进行双能X射线相位衬度成像时,选取的两个能量V高和V低可以根据实际情况随意调整。与目前的技术相比(选用的两个能量V高和V低不可调节),大大拓展了双能X射线相位衬度成像的应用范围。
(2)由于该双能X射线相位衬度成像装置采用的三块光栅均为几十微米(40至100微米)的大周期吸收光栅,与目前的技术相比(采用的光栅周期为2微米左右),本成像***光栅的研制加工更加容易、成本更低,并且更有可能研制出大周期(40×40厘米2)、高吸收厚度(300微米)和高性能的光栅,从而可以使双能X射线相位衬度成像技术应用在临床医学成像领域。
(3)该双能X射线相位衬度成像装置采用的三块光栅均为几十微米的大周期光栅,与目前采用2微米左右小周期的双能X射线相位衬度成像技术相比,根据本发明构建的双能X射线相位衬度成像装置对机械稳定性要求大大降低,机架重力、机械制造误差、机械装配误差、机械振动和热膨胀等外界因素对***的成像性能影响会大大降低。
附图说明
图1是本发明提供的双能X射线相位衬度成像装置的光路示意图。
图2是本发明提供的实现双能X射线相位衬度成像的方法流程图。
图3是依照本发明实施例的双能X射线相位衬度成像装置的结构示意图。
图4是依照本发明实施例的采用四根有机玻璃圆柱棒的双能X射线相位衬度成像实验结果的示意图。
图5是依照本发明实施例的X光机在40KV管电压下成像结果的定量分析的示意图。
图6是依照本发明实施例的X光机在50KV管电压下成像结果的定量分析的示意图。
图7是依照本发明实施例的X光机在40KV和50KV管电压下的出射能谱的示意图。
图8是依照本发明实施例的X光机在40KV管电压和50KV管电压下的试验结果对比的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明采用的技术思路是:利用40至100微米的大周期光栅构建双能X射线相位衬度成像装置,借助光栅几何投影产生的莫尔条纹提取样品的相位衬度图像,进而实施双能X射线相位衬度成像。
图1是本发明提供的双能X射线相位衬度成像装置的光路示意图,该装置沿光路依次包括X光机1、源光栅2、分束光栅3、样品室4、分析光栅5和X射线探测器6。其中,X光机1用于发出X射线;源光栅2用于将大焦点的X光源分成为若干个不相干的小焦点光源;分束光栅3用于将小焦点光源分成若干束,照射到样品室4中的样品上,并在在分析光栅5上形成几何投影;样品室4用于放置并固定样品,同时带动样品进行旋转;分析光栅5用于与分束光栅3一起在X射线探测器上形成莫尔条纹;X射线探测器6用于获取并记录该莫尔条纹。
其中,所述源光栅2与所述分析光栅5的周期满足其中p0为源光栅2的周期,p2是分析光栅5的周期,l是源光栅2与分束光栅3之间的距离,d为分束光栅3与分析光栅5之间的距离,以使所述分束光栅3在每个小焦点光源下产生的投影像互相错开一个周期并叠加在一起,从而使得条纹对比度得到增强。
同时,所述分束光栅3与所述分析光栅5的周期满足其中p1为分束光栅3的周期,p2是分析光栅5的周期,l是源光栅2与分束光栅3之间的距离,d为分束光栅3与分析光栅5之间的距离,以使所述分束光栅3投影的周期和所述分析光栅5周期相等,从而能够通过相位步进技术完成信息提取。
另外,所述源光栅2、分束光栅3与分析光栅5均为周期为40微米至100微米的大周期吸收光栅,X射线探测器6为X射线平板探测器。
基于图1所示的双能X射线相位衬度成像装置的光路示意图,图2示出了本发明提供的实现双能X射线相位衬度成像的方法流程图,该方法包括以下步骤:
步骤1:对准X光机的出光点和X射探测器接收平面的中心点;
步骤2:对X射线探测器进行Offset矫正和Gain矫正;
步骤3:精确对准源光栅、分束光栅与分析光栅,使这三块光栅的刻线相互平行、三块光栅所在的平面和X射线探测器的平面相互平行,同时使三块光栅的中心点、X光机的出光点和X射线探测器中心点在一条直线上;
步骤4:调节X光机的管电压到高能V高处,进行相位步进扫描,同时旋转样品室,采集图像并重建出高能V高处样品的三维相位信息;
步骤5:调节X光机的管电压到低能V低处,进行相位步进扫描,同时旋转样品室,采集图像并重建出低能V低处样品的三维相位信息。
基于图1所示的双能X射线相位衬度成像装置和图2所示的实现双能X射线相位衬度成像的方法,图3示出了依照本发明实施例的双能X射线相位衬度成像装置的结构示意图,其中,X光机、源光栅2、分束光栅3、样品室和分析光栅5均通过组装的三维光学精密位移台固定在光学平台上,在分析光栅5横向运动方向,配置有超精密的压电陶瓷电机,定位精度可达到10纳米,用于高精密的相位步进扫描,X射线平板探测器通过机械部件固定在桌子上。
在本实施例中,整个双能X射线相位衬度成像装置的电机控制、图像采集、图像处理和装置维护等工作可以由工业计算机上的LabVIEW程序完成。
在本实施例中,X光机的焦点大小为1毫米,管电压可调范围为30-160KV。源光栅周期为100微米,金厚度为200微米,占空比为1∶1,面积为1×1厘米2。分束光栅周期为50微米,金厚度为200微米,占空比为1∶1,面积为10×10厘米2。分析光栅周期为100微米,金厚度为200微米,占空比为1∶1,面积为10×10厘米2。X射线探测器的像素尺寸为0.2×0.2厘米2,X射线探测器的活动区域面积为20.48×20.48厘米2。源光栅距离X光机出光点的距离为1厘米,源光栅和分束光栅的距离为60厘米,样品室紧贴着分束光栅,分析光栅和分束光栅的距离为60厘米,X射线探测器采用X射线平板探测器,其紧贴着分析光栅放置。
基于该双能X射线相位衬度成像装置,其实现双能X射线相位衬度成像的步骤如下:
(1)利用激光器分别在水平方向和竖直方向,对准X光机出射点和X射线探测器活动区域的中心,使X光机的出光点与X射线探测器接收平面的中心点对准。
(2)根据实验对象,选用合适的高能V高和低能V低进行双能X射线相位衬度成像。如实验对象为小白鼠,高能V高一般选为60KV,低能V低一般选为40KV。
(3)把源光栅2、分束光栅3和分析光栅5移出成像视场,在X光机关闭时利用X射线探测器采集50张图像,平均后保存为一张Offset图像(Ioffset);打开X光机,把X光机的管电压调到低能V低处,稳定10分钟后,利用X射线探测器采集50张图像,平均后保存为一张低能Gain图像把X光机的管电压调到高能V高处,稳定10分钟后,利用X射线探测器采集50张图像,平均后保存为一张高能Gain图像
(4)把X光机的管电压维持在高能V高处,调用Ioffset和利用公式对采集到的Offset图像(Ioffset)和高能Gain图像进行offset校正和gain校正。其中I(m,n)是校正完毕后图像中像素(m,n)的灰度值,Iacquire(m,n)是探测器采集的图像中像素(m,n)的灰度值,Ioffset(m,n)是Offset图像中像素(m,n)的灰度值,是Gain图像中像素(m,n)的灰度值。M是探测器水平方向的像素数量,N是探测器垂直方向的像素数量。
(5)把源光栅2移动至视场的正中央,在Z方向轻微旋转源光栅2,观测X射线探测器采集到的图像,使源光栅2的刻线处于竖直方向。
(6)把分束光栅3和分析光栅5移动进入成像视场内,此时在探测器上会看到倾斜的大周期莫尔条纹。
(7)在Z方向轻微旋转分束光栅3,每次旋转0.1°,同时观察探测器接收到的图像,反复调整,使莫尔条纹处于垂直方向。
(8)在Y方向轻微旋转分束光栅3,每次旋转0.1°,同时观察探测器接收到的图像,反复调整,使视场中的莫尔条纹的周期左右相等。
(9)在X方向轻微旋转分析光栅5,每次旋转0.1°,同时观察探测器接收到的图像,反复调整,使视场中的莫尔条纹的周期上下相等。
(10)在Z方向移动分束光栅3,每次移动0.1毫米,使莫尔条纹的周期逐渐增大,反复调整,直到在整个视场区域,光强保持均匀,没有明显的明暗之分。
(11)在X方向高精度移动分析光栅5,每次移动10微米,在移动过程中,观测探测器采集到图像的明暗变化是否均匀。若不均匀,重新在X方向轻微移动分束光栅3,直到在一个周期内连续移动分析光栅5的过程中,成像视场内各个位置的光强从亮到暗的变化是步调一致的。
(12)把X光机的管电压维持在高能V高处,把样品室移动到成像视场中。
(13)在X方向高精度移动分析光栅5,进行相位步进扫描,同时利用X射线探测器采集并保存图像,整个过程由LabVIEW程序自动完成。相位步进扫描完成后,利用公式计算样品在高能V高处的折射信息。公式中θ(m,n)是像素(m,n)处的折射角,p2是分析光栅5的周期,d是样品室和分析光栅5之间的距离,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,k是累加求和过程中的变量,N是相位步进的步数。
(14)在Y方向把样品室旋转1°,重复步骤13。
(15)重复步骤14,直到样品室完成360°的旋转为止。
(16)把样品室移出视场,在X方向高精度移动分析光栅5,进行相位步进扫描,同时利用探测器采集并保存图像,整个过程由LabVIEW程序自动完成。相位步进扫描完成后,利用公式计算出高能V高处成像***背景的折射信息。公式中θ(m,n)是像素(m,n)处的折射角,p2是分析光栅5的周期,d是分束光栅3和分析光栅5之间的距离,N是相位步进的步数,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,k是累加求和过程中的变量。
(17)利用公式扣除***背景,计算出样品室旋转过程中,每个投影方向纯样品的折射信息。
(18)利用优化的滤波反投影算法重建出试验样品在高能V高处的三维相位信息。其中δ(x,y,z)高能为样品位置(x,y,z)处体积单元在高能V高下折射率的实部,φ是样品围绕旋转轴转动的角度变量,xr是随样品转动的坐标变量,x和y是垂直于光轴平面内的两个坐标变量,z是在垂直于xy平面方向的坐标变量,θr(xr,φ,z)是样品在φ角度位置处投影图像中像素(xr,z)处的折射角,F-1代表傅里叶反变换,ρ是频域区间的坐标变量,j是复数δ(xcosφ+ysinφ-xr)是狄拉克δ函数。
(19)把X光机的管电压调到低能V低处,调用Ioffset和利用公式对采集到的图像进行offset校正和gain校正。其中I(m,n)是校正完毕后图像中像素(m,n)的灰度值,Iacquire(m,n)是探测器采集的图像中像素(m,n)的灰度值,Ioffset(m,n)是Offset图像中像素(m,n)的灰度值,是Gain图像中像素(m,n)的灰度值。M是探测器水平方向的像素数量,N是探测器垂直方向的像素数量。
(20)在X方向高精度移动分析光栅5,进行相位步进扫描,同时利用探测器采集并保存图像,整个过程由LabVIEW程序自动完成。相位步进扫描完成后,利用公式计算低能V低处样品的折射信息。公式中θ(m,n)是像素(m,n)处的折射角,N是相位步进的步数,k是累加求和过程中的变量,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,p2是分析光栅5的周期,d是样品室和分析光栅5之间的距离。
(21)在Y方向把样品室旋转1°,重复步骤20。
(22)重复步骤21,直到样品室完成360°的转动为止。
(23)把样品室移出视场,在X方向高精度移动分析光栅5,进行相位步进扫描,同时利用探测器采集并保存图像,整个过程由LabVIEW程序自动完成。相位步进扫描完成后,利用公式计算出高能V低处成像***背景的折射信息。公式中θ(m,n)是像素(m,n)处的折射角,N是相位步进的步数,k是累加求和过程中的变量,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,p2是分析光栅5的周期,d是分束光栅3和分析光栅5之间的距离。
(24)利用公式扣除***背景,计算出样品室旋转过程中,每个投影方向纯样品的折射信息。
(25)利用优化的滤波反投影算法重建出试验样品在低能V低的三维相位信息。其中δ(x,y,z)高能为样品位置(x,y,z)处体积单元在高能V高下折射率的实部,φ是样品围绕旋转轴转动的角度变量,xr是随样品转动的坐标变量,x和y是垂直于光轴平面内的两个坐标变量,z是在垂直于xy平面方向的坐标变量,θr(xr,φ,z)是样品在φ角度位置处投影图像中像素(xr,z)处的折射角,F-1代表傅里叶反变换,ρ是频域区间的坐标变量,j是复数δ(xcosφ+ysinφ-xr)是狄拉克δ函数。
实验结果
图4是依照本发明实施例的采用四根有机玻璃圆柱棒的双能X射线相位衬度成像实验结果的示意图,四根有机玻璃圆柱棒从左到右依次是:直径为10mm的PMMA,直径为5mm的POM,直径为5mm的PMMA和直径为10mm的POM;X光机在50KV管电压下的成像结果为:(A)吸收图像,(B)折射图像和(C)散射图像;X光机在40KV管电压下的成像结果为:(D)吸收图像,(E)折射图像和(F)散射图像。
图5是依照本发明实施例的X光机在40KV管电压下成像结果的定量分析的示意图,其中,空心正方形符号表示的曲线是实验值,黑色实线为拟合值,(A)是直径为10mm PMMA圆柱棒吸收图像的Profile曲线及拟合曲线,(B)是直径为10mm POM圆柱棒吸收图像的Profile曲线及拟合曲线,(C)是直径为10mm PMMA圆柱棒折射图像的Profile曲线及拟合曲线,(D)是直径为10mm POM圆柱棒折射图像的Profile曲线及拟合曲线。
图6是依照本发明实施例的X光机在50KV管电压下成像结果的定量分析的示意图,其中,空心正方形符号表示的曲线是实验值,黑色实线为拟合值,(A)是直径为10mm PMMA圆柱棒吸收图像的Profile曲线及拟合曲线,(B)是直径为10mm POM圆柱棒吸收图像的Profile曲线及拟合曲线,(C)是直径为10mm PMMA圆柱棒折射图像的Profile曲线及拟合曲线,(D)是直径为10mm POM圆柱棒折射图像的Profile曲线及拟合曲线。
图7是依照本发明实施例的X光机在40KV(实心圆圈符号和黑色线条表示的曲线)和50KV(实心正方形符号和黑色线条表示的曲线)管电压下的出射能谱的示意图。
图8是依照本发明实施例的X光机在40KV管电压(黑色实线表示的曲线)和50KV管电压(黑色虚线表示的曲线)下的试验结果对比的示意图,其中,(A)是直径为10mm PMMA圆柱棒的吸收曲线,(B)是直径为10mmPOM圆柱棒的吸收曲线,(C)是直径为10mm PMMA圆柱棒的折射曲线,(D)是直径为10mm POM圆柱棒的折射曲线。
表1是PMMA和POM圆柱棒折射率实部和虚部的实验拟合值和理论计算值。
表1
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (11)
1.一种双能X射线相位衬度成像装置,其特征在于,该装置沿光路依次包括X光机(1)、源光栅(2)、分束光栅(3)、样品室(4)、分析光栅(5)和X射线探测器(6),其中:
X光机(1),用于发出X射线;
源光栅(2),用于将大焦点的X光源分成为若干个不相干的小焦点光源;
分束光栅(3),用于将小焦点光源分成若干束,照射到样品室(4)中的样品上,并在在分析光栅(5)上形成几何投影;
样品室(4),用于放置并固定样品,同时带动样品进行旋转;
分析光栅(5),用于与分束光栅一起在X射线探测器上形成莫尔条纹;
X射线探测器(6),用于获取并记录该莫尔条纹。
2.根据权利要求1所述的双能X射线相位衬度成像装置,其特征在于,所述源光栅(2)与所述分析光栅(5)的周期满足其中p0为源光栅(2)的周期,p2是分析光栅(5)的周期,l是源光栅(2)与分束光栅(3)之间的距离,d为分束光栅(3)与分析光栅(5)之间的距离,以使所述分束光栅(3)在每个小焦点光源下产生的投影像互相错开一个周期并叠加在一起,从而使得条纹对比度得到增强。
3.根据权利要求1所述的双能X射线相位衬度成像装置,其特征在于,所述分束光栅(3)与所述分析光栅(5)的周期满足其中p1为分束光栅(3)的周期,p2是分析光栅(5)的周期,l是源光栅(2)与分束光栅(3)之间的距离,d为分束光栅(3)与分析光栅(5)之间的距离,以使所述分束光栅(3)投影的周期和所述分析光栅(5)周期相等,从而能够通过相位步进技术完成信息提取。
4.根据权利要求1所述的双能X射线相位衬度成像装置,其特征在于,所述源光栅(2)、分束光栅(3)与分析光栅(5)均为周期为40微米至100微米的大周期吸收光栅。
5.根据权利要求1所述的双能X射线相位衬度成像装置,其特征在于,所述X射线探测器(6)为X射线平板探测器。
6.一种实现双能X射线相位衬度成像的方法,应用于权利要求1至5中任一项所述的双能X射线相位衬度成像装置,其特征在于,该方法包括:
对准X光机的出光点和X射探测器接收平面的中心点;
对X射线探测器进行Offset矫正和Gain矫正;
精确对准源光栅、分束光栅与分析光栅,使这三块光栅的刻线相互平行、三块光栅所在的平面和X射线探测器的平面相互平行,同时使三块光栅的中心点、X光机的出光点和X射线探测器中心点在一条直线上;
调节X光机的管电压到高能V高处,进行相位步进扫描,同时旋转样品室,采集图像并重建出高能V高处样品的三维相位信息;
调节X光机的管电压到低能V高处,进行相位步进扫描,同时旋转样品室,采集图像并重建出低能V低处样品的三维相位信息。
7.根据权利要求6所述的实现双能X射线相位衬度成像的方法,其特征在于,所述对准X光机的出光点和X射探测器接收平面的中心点,是利用激光器分别在水平方向和竖直方向,对准X光机出射点和X射线探测器活动区域的中心,使X光机的出光点与X射线探测器接收平面的中心点对准。
8.根据权利要求6所述的实现双能X射线相位衬度成像的方法,其特征在于,所述对X射线探测器进行Offset矫正和Gain矫正,包括:
将源光栅、分束光栅和分析光栅移出成像视场,在X光机关闭时利用X射线探测器采集50张图像,平均后保存为一张Offset图像(Ioffset);打开X光机,把X光机的管电压调到低能V低处,稳定10分钟后,利用X射线探测器采集50张图像,平均后保存为一张低能Gain图像把X光机的管电压调到高能V高处,稳定10分钟后,利用X射线探测器采集50张图像,平均后保存为一张高能Gain图像
将X光机的管电压维持在高能V高处,调用Offset图像(Ioffset)和高能Gain图像利用公式对采集到的Offset图像(Ioffset)和高能Gain图像进行offset校正和gain校正;其中I(m,n)是校正完毕后图像中像素(m,n)的灰度值,Iacquire(m,n)是探测器采集的图像中像素(m,n)的灰度值,Ioffset(m,n)是Offset图像中像素(m,n)的灰度值,是Gain图像中像素(m,n)的灰度值;M是探测器水平方向的像素数量,N是探测器垂直方向的像素数量。
9.根据权利要求6所述的实现双能X射线相位衬度成像的方法,其特征在于,所述精确对准源光栅、分束光栅与分析光栅,使这三块光栅的刻线相互平行、三块光栅所在的平面和X射线探测器的平面相互平行,同时使三块光栅的中心点、X光机的出光点和X射线探测器中心点在一条直线上,包括:
把源光栅移动至视场的正中央,在Z方向轻微旋转源光栅,观测X射线探测器采集到的图像,使源光栅的刻线处于竖直方向;
把分束光栅和分析光栅移动进入成像视场内,此时在探测器上会看到倾斜的大周期莫尔条纹;
在Z方向轻微旋转分束光栅,每次旋转0.1°,同时观察探测器接收到的图像,反复调整,使莫尔条纹处于垂直方向;
在Y方向轻微旋转分束光栅,每次旋转0.1°,同时观察探测器接收到的图像,反复调整,使视场中的莫尔条纹的周期左右相等;
在X方向轻微旋转分析光栅,每次旋转0.1°,同时观察探测器接收到的图像,反复调整,使视场中的莫尔条纹的周期上下相等;
在Z方向移动分束光栅,每次移动0.1毫米,使莫尔条纹的周期逐渐增大,反复调整,直到在整个视场区域,光强保持均匀,没有明显的明暗之分;
在X方向高精度移动分析光栅,每次移动10微米,在移动过程中,观测探测器采集到图像的明暗变化是否均匀;若不均匀,重新在X方向轻微移动分束光栅,直到在一个周期内连续移动分析光栅的过程中,成像视场内各个位置的光强从亮到暗的变化是步调一致的。
10.根据权利要求6所述的实现双能X射线相位衬度成像的方法,其特征在于,所述调节X光机的管电压到高能V高处,进行相位步进扫描,同时旋转样品室,采集图像并重建出高能V高处样品的三维相位信息,包括:
步骤12:把X光机的管电压维持在高能V高处,把样品室移动到成像视场中;
步骤13:在X方向高精度移动分析光栅,进行相位步进扫描,同时利用X射线探测器采集并保存图像,相位步进扫描完成后,利用公式计算样品在高能V高处的折射信息;公式中θ(m,n)是像素(m,n)处的折射角,p2是分析光栅的周期,d是样品室和分析光栅之间的距离,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,k是累加求和过程中的变量,N是相位步进的步数;
步骤14:在Y方向把样品室旋转1°,重复步骤13;
步骤15:重复步骤14,直到样品室完成360°的旋转为止;
步骤16:把样品室移出视场,在X方向高精度移动分析光栅,进行相位步进扫描,同时利用探测器采集并保存图像,相位步进扫描完成后,利用公式计算出高能V高处成像***背景的折射信息;公式中θ(m,n)是像素(m,n)处的折射角,p2是分析光栅的周期,d是分束光栅与分析光栅之间的距离,N是相位步进的步数,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,k是累加求和过程中的变量;
步骤17:利用公式扣除***背景,计算出样品室旋转过程中,每个投影方向纯样品的折射信息;
步骤18:利用优化的滤波反投影算法重建出试验样品在高能V高处的三维相位信息;其中δ(x,y,z)高能为样品位置(x,y,z)处体积单元在高能V高下折射率的实部,φ是样品围绕旋转轴转动的角度变量,xr是随样品转动的坐标变量,x和y是垂直于光轴平面内的两个坐标变量,z是在垂直于xy平面方向的坐标变量,θr(xr,φ,z)是样品在φ角度位置处投影图像中像素(xr,z)处的折射角,F-1代表傅里叶反变换,ρ是频域区间的坐标变量,j是复数5(xcosφ+ysinφ-xr)是狄拉克δ函数。
11.根据权利要求6所述的实现双能X射线相位衬度成像的方法,其特征在于,所述调节X光机的管电压到低能V低处,进行相位步进扫描,同时旋转样品室,采集图像并重建出低能V低处样品的三维相位信息,包括:
步骤19:把X光机的管电压调到低能V低处,调用Ioffset和利用公式对采集到的图像进行offset校正和gain校正;其中I(m,n)是校正完毕后图像中像素(m,n)的灰度值,Iacquire(m,n)是探测器采集的图像中像素(m,n)的灰度值,Ioffset(m,n)是Offset图像中像素(m,n)的灰度值,是Gain图像中像素(m,n)的灰度值;M是探测器水平方向的像素数量,N是探测器垂直方向的像素数量;
步骤20:在X方向高精度移动分析光栅,进行相位步进扫描,同时利用探测器采集并保存图像,相位步进扫描完成后,利用公式计算低能V低处样品的折射信息;公式中θ(m,n)是像素(m,n)处的折射角,N是相位步进的步数,k是累加求和过程中的变量,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,p2是分析光栅的周期,d是样品室和分析光栅之间的距离;
步骤21:在Y方向把样品室旋转1°,重复步骤20;
步骤22:重复步骤21,直到样品室完成360°的转动为止;
步骤23:把样品室移出视场,在X方向高精度移动分析光栅,进行相位步进扫描,同时利用探测器采集并保存图像,相位步进扫描完成后,利用公式计算出高能V低处成像***背景的折射信息;公式中θ(m,n)是像素(m,n)处的折射角,N是相位步进的步数,k是累加求和过程中的变量,Ik(m,n)是第k步采集到的图像中像素(m,n)的灰度值,p2是分析光栅的周期,d是分束光栅与分析光栅之间的距离;
步骤24:利用公式扣除***背景,计算出样品室旋转过程中,每个投影方向纯样品的折射信息;
步骤25:利用优化的滤波反投影算法重建出试验样品在低能V低的三维相位信息;其中δ(x,y,z)高能为样品位置(x,y,z)处体积单元在高能V高下折射率的实部,φ是样品围绕旋转轴转动的角度变量,xr是随样品转动的坐标变量,x和y是垂直于光轴平面内的两个坐标变量,z是在垂直于xy平面方向的坐标变量,θr(xr,φ,z)是样品在φ角度位置处投影图像中像素(xr,z)处的折射角,F-1代表傅里叶反变换,ρ是频域区间的坐标变量,j是复数δ(xcosφ+ysinφ-xr)是狄拉克δ函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410384492.9A CN104132953B (zh) | 2014-08-01 | 2014-08-06 | 一种双能x射线相位衬度成像装置及其实现方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2014103771294 | 2014-08-01 | ||
CN201410377129 | 2014-08-01 | ||
CN201410377129.4 | 2014-08-01 | ||
CN201410384492.9A CN104132953B (zh) | 2014-08-01 | 2014-08-06 | 一种双能x射线相位衬度成像装置及其实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104132953A true CN104132953A (zh) | 2014-11-05 |
CN104132953B CN104132953B (zh) | 2017-03-29 |
Family
ID=51805717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410384492.9A Expired - Fee Related CN104132953B (zh) | 2014-08-01 | 2014-08-06 | 一种双能x射线相位衬度成像装置及其实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104132953B (zh) |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104458777A (zh) * | 2014-12-30 | 2015-03-25 | 中国科学技术大学 | 用于x射线相位衬度成像的积分水桶相位测量方法 |
CN104535595A (zh) * | 2015-01-09 | 2015-04-22 | 中国科学技术大学 | 一种用于x射线光栅相位衬度成像的背景扣除方法 |
CN104869905A (zh) * | 2012-12-21 | 2015-08-26 | 卡尔斯特里姆保健公司 | 基于微分相衬成像的医疗放射照相光栅 |
CN105388170A (zh) * | 2015-12-17 | 2016-03-09 | 中国科学技术大学 | 一种骨密度的测量方法及装置 |
CN105486693A (zh) * | 2016-01-19 | 2016-04-13 | 中国科学技术大学 | 一种无损检测高精度元件缺陷的方法 |
CN105651793A (zh) * | 2016-01-05 | 2016-06-08 | 合肥泰禾光电科技股份有限公司 | 一种克服物体厚度影响的x光检测方法 |
CN105852895A (zh) * | 2016-04-29 | 2016-08-17 | 合肥工业大学 | 单次曝光的硬x射线光栅干涉仪的信息提取新方法 |
CN105935297A (zh) * | 2016-06-23 | 2016-09-14 | 中国科学院深圳先进技术研究院 | X射线光栅相衬成像ct*** |
CN106023107A (zh) * | 2016-05-17 | 2016-10-12 | 中国科学技术大学 | 一种用于x射线光栅相衬成像装置的探测器图像校正方法 |
CN105997127A (zh) * | 2016-06-21 | 2016-10-12 | 深圳先进技术研究院 | 一种静态乳腺双能ct成像***及成像方法 |
CN106404809A (zh) * | 2016-07-27 | 2017-02-15 | 中国科学技术大学 | 一种用于x射线光栅相衬成像装置的图像校正方法 |
CN107024490A (zh) * | 2016-01-29 | 2017-08-08 | 中国科学院高能物理研究所 | 一次曝光光栅剪切成像装置及数据采集与信息提取方法 |
CN107085000A (zh) * | 2017-04-17 | 2017-08-22 | 深圳先进技术研究院 | X射线光栅相衬成像自动校准***及方法 |
CN107144581A (zh) * | 2017-05-05 | 2017-09-08 | 北京航空航天大学 | 基于横向错位吸收光栅的x射线光栅差分相位衬度成像方法及装置 |
CN107748341A (zh) * | 2017-10-23 | 2018-03-02 | 中国科学院苏州生物医学工程技术研究所 | 高衬度低剂量相位衬度ct成像装置 |
CN107807139A (zh) * | 2016-09-05 | 2018-03-16 | 天津工业大学 | 一种无步进装置的双能x射线相衬成像***及其实现方法 |
CN109283200A (zh) * | 2017-07-20 | 2019-01-29 | 株式会社岛津制作所 | X射线相位成像装置和含有纤维的材料的缺陷检测方法 |
CN109425624A (zh) * | 2017-09-05 | 2019-03-05 | 株式会社岛津制作所 | X射线成像装置 |
CN109580667A (zh) * | 2018-12-05 | 2019-04-05 | 中国科学技术大学 | 单光栅相衬成像方法及*** |
CN110068585A (zh) * | 2019-03-15 | 2019-07-30 | 中国工程物理研究院流体物理研究所 | 双能x射线光栅干涉成像***及方法 |
US10578563B2 (en) | 2012-12-21 | 2020-03-03 | Carestream Health, Inc. | Phase contrast imaging computed tomography scanner |
CN111505034A (zh) * | 2020-04-30 | 2020-08-07 | 合肥工业大学 | 一种基于迭代算法的x射线衍射增强成像方法 |
CN111721786A (zh) * | 2019-03-22 | 2020-09-29 | 中国科学院深圳先进技术研究院 | 一种x射线干涉仪及成像*** |
CN112258628A (zh) * | 2020-10-13 | 2021-01-22 | 南京大学 | 一种基于光栅相位衬度的成像***及方法 |
CN112568923A (zh) * | 2020-12-10 | 2021-03-30 | 中国科学院深圳先进技术研究院 | X射线相位衬度图像提取方法、装置、终端及存储介质 |
CN113237901A (zh) * | 2021-05-07 | 2021-08-10 | 中国科学院上海应用物理研究所 | 一种生物特征识别***及生物特征识别方法 |
CN115227275A (zh) * | 2022-09-23 | 2022-10-25 | 中国科学技术大学 | 一种旋转阻挡器和基于该旋转阻挡器的散射伪影校正方法 |
CN116878420A (zh) * | 2023-05-09 | 2023-10-13 | 中国科学院高能物理研究所 | 一种光栅周期失配设计的光栅干涉波前检测方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101467889A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院高能物理研究所 | 光栅剪切相位衬度ct成像数据采集和重建方法 |
WO2010146503A1 (en) * | 2009-06-16 | 2010-12-23 | Koninklijke Philips Electronics N. V. | Correction method for differential phase contrast imaging |
CN102221565A (zh) * | 2010-04-19 | 2011-10-19 | 清华大学 | X射线源光栅步进成像***与成像方法 |
CN102740775A (zh) * | 2010-02-04 | 2012-10-17 | 富士胶片株式会社 | 放射成像*** |
CN103364839A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 基于光栅剪切成像的安检设备及方法 |
-
2014
- 2014-08-06 CN CN201410384492.9A patent/CN104132953B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101467889A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院高能物理研究所 | 光栅剪切相位衬度ct成像数据采集和重建方法 |
WO2010146503A1 (en) * | 2009-06-16 | 2010-12-23 | Koninklijke Philips Electronics N. V. | Correction method for differential phase contrast imaging |
CN102740775A (zh) * | 2010-02-04 | 2012-10-17 | 富士胶片株式会社 | 放射成像*** |
CN102221565A (zh) * | 2010-04-19 | 2011-10-19 | 清华大学 | X射线源光栅步进成像***与成像方法 |
CN103364839A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院高能物理研究所 | 基于光栅剪切成像的安检设备及方法 |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104869905A (zh) * | 2012-12-21 | 2015-08-26 | 卡尔斯特里姆保健公司 | 基于微分相衬成像的医疗放射照相光栅 |
CN104869905B (zh) * | 2012-12-21 | 2019-08-06 | 卡尔斯特里姆保健公司 | 基于微分相衬成像的医疗放射照相光栅 |
US10578563B2 (en) | 2012-12-21 | 2020-03-03 | Carestream Health, Inc. | Phase contrast imaging computed tomography scanner |
CN104458777A (zh) * | 2014-12-30 | 2015-03-25 | 中国科学技术大学 | 用于x射线相位衬度成像的积分水桶相位测量方法 |
CN104535595B (zh) * | 2015-01-09 | 2017-05-31 | 中国科学技术大学 | 一种用于x射线光栅相位衬度成像的背景扣除方法 |
CN104535595A (zh) * | 2015-01-09 | 2015-04-22 | 中国科学技术大学 | 一种用于x射线光栅相位衬度成像的背景扣除方法 |
CN105388170A (zh) * | 2015-12-17 | 2016-03-09 | 中国科学技术大学 | 一种骨密度的测量方法及装置 |
CN105388170B (zh) * | 2015-12-17 | 2018-07-03 | 中国科学技术大学 | 一种骨密度的测量方法及装置 |
CN105651793A (zh) * | 2016-01-05 | 2016-06-08 | 合肥泰禾光电科技股份有限公司 | 一种克服物体厚度影响的x光检测方法 |
CN105651793B (zh) * | 2016-01-05 | 2019-04-02 | 合肥泰禾光电科技股份有限公司 | 一种克服物体厚度影响的x光检测方法 |
CN105486693A (zh) * | 2016-01-19 | 2016-04-13 | 中国科学技术大学 | 一种无损检测高精度元件缺陷的方法 |
CN107024490A (zh) * | 2016-01-29 | 2017-08-08 | 中国科学院高能物理研究所 | 一次曝光光栅剪切成像装置及数据采集与信息提取方法 |
CN107024490B (zh) * | 2016-01-29 | 2019-07-05 | 中国科学院高能物理研究所 | 一次曝光光栅剪切成像装置及数据采集与信息提取方法 |
CN105852895A (zh) * | 2016-04-29 | 2016-08-17 | 合肥工业大学 | 单次曝光的硬x射线光栅干涉仪的信息提取新方法 |
CN106023107B (zh) * | 2016-05-17 | 2020-03-27 | 中国科学技术大学 | 一种用于x射线光栅相衬成像装置的探测器图像校正方法 |
CN106023107A (zh) * | 2016-05-17 | 2016-10-12 | 中国科学技术大学 | 一种用于x射线光栅相衬成像装置的探测器图像校正方法 |
CN105997127A (zh) * | 2016-06-21 | 2016-10-12 | 深圳先进技术研究院 | 一种静态乳腺双能ct成像***及成像方法 |
CN105935297A (zh) * | 2016-06-23 | 2016-09-14 | 中国科学院深圳先进技术研究院 | X射线光栅相衬成像ct*** |
CN106404809A (zh) * | 2016-07-27 | 2017-02-15 | 中国科学技术大学 | 一种用于x射线光栅相衬成像装置的图像校正方法 |
CN107807139A (zh) * | 2016-09-05 | 2018-03-16 | 天津工业大学 | 一种无步进装置的双能x射线相衬成像***及其实现方法 |
CN107085000A (zh) * | 2017-04-17 | 2017-08-22 | 深圳先进技术研究院 | X射线光栅相衬成像自动校准***及方法 |
CN107085000B (zh) * | 2017-04-17 | 2020-06-16 | 深圳先进技术研究院 | X射线光栅相衬成像自动校准***及方法 |
CN107144581B (zh) * | 2017-05-05 | 2019-09-27 | 北京航空航天大学 | 基于横向错位吸收光栅的x射线光栅差分相位衬度成像方法及装置 |
CN107144581A (zh) * | 2017-05-05 | 2017-09-08 | 北京航空航天大学 | 基于横向错位吸收光栅的x射线光栅差分相位衬度成像方法及装置 |
CN109283200A (zh) * | 2017-07-20 | 2019-01-29 | 株式会社岛津制作所 | X射线相位成像装置和含有纤维的材料的缺陷检测方法 |
CN109283200B (zh) * | 2017-07-20 | 2022-09-16 | 株式会社岛津制作所 | X射线相位成像装置和含有纤维的材料的缺陷检测方法 |
CN109425624A (zh) * | 2017-09-05 | 2019-03-05 | 株式会社岛津制作所 | X射线成像装置 |
CN109425624B (zh) * | 2017-09-05 | 2022-06-03 | 株式会社岛津制作所 | X射线成像装置 |
CN107748341A (zh) * | 2017-10-23 | 2018-03-02 | 中国科学院苏州生物医学工程技术研究所 | 高衬度低剂量相位衬度ct成像装置 |
CN109580667A (zh) * | 2018-12-05 | 2019-04-05 | 中国科学技术大学 | 单光栅相衬成像方法及*** |
CN110068585B (zh) * | 2019-03-15 | 2022-03-22 | 中国工程物理研究院流体物理研究所 | 双能x射线光栅干涉成像***及方法 |
CN110068585A (zh) * | 2019-03-15 | 2019-07-30 | 中国工程物理研究院流体物理研究所 | 双能x射线光栅干涉成像***及方法 |
CN111721786A (zh) * | 2019-03-22 | 2020-09-29 | 中国科学院深圳先进技术研究院 | 一种x射线干涉仪及成像*** |
CN111505034A (zh) * | 2020-04-30 | 2020-08-07 | 合肥工业大学 | 一种基于迭代算法的x射线衍射增强成像方法 |
CN112258628A (zh) * | 2020-10-13 | 2021-01-22 | 南京大学 | 一种基于光栅相位衬度的成像***及方法 |
CN112568923A (zh) * | 2020-12-10 | 2021-03-30 | 中国科学院深圳先进技术研究院 | X射线相位衬度图像提取方法、装置、终端及存储介质 |
CN113237901A (zh) * | 2021-05-07 | 2021-08-10 | 中国科学院上海应用物理研究所 | 一种生物特征识别***及生物特征识别方法 |
CN115227275A (zh) * | 2022-09-23 | 2022-10-25 | 中国科学技术大学 | 一种旋转阻挡器和基于该旋转阻挡器的散射伪影校正方法 |
CN115227275B (zh) * | 2022-09-23 | 2022-12-30 | 中国科学技术大学 | 一种旋转阻挡器和基于该旋转阻挡器的散射伪影校正方法 |
CN116878420A (zh) * | 2023-05-09 | 2023-10-13 | 中国科学院高能物理研究所 | 一种光栅周期失配设计的光栅干涉波前检测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104132953B (zh) | 2017-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104132953A (zh) | 一种双能x射线相位衬度成像装置及其实现方法 | |
CN105606633B (zh) | X射线相衬成像***与成像方法 | |
Pfeiffer et al. | Grating-based X-ray phase contrast for biomedical imaging applications | |
CN103913779A (zh) | 多能ct成像***以及成像方法 | |
CN105852895B (zh) | 单次曝光的硬x射线光栅干涉仪的信息提取方法 | |
CN101467889B (zh) | 光栅剪切相位衬度ct成像数据采集和重建方法 | |
RU2695311C2 (ru) | Устройство рентгеновской визуализации | |
CN105935297A (zh) | X射线光栅相衬成像ct*** | |
Pfeiffer | Milestones and basic principles of grating-based x-ray and neutron phase-contrast imaging | |
JP2015213681A (ja) | 等価ファントム、および等価ファントムを用いたx線タルボ撮影装置の品質評価方法 | |
CN115363617A (zh) | 一种微分相位衬度ct的采集和重建方法 | |
Pfeiffer et al. | Region-of-interest tomography for grating-based x-ray differential phase-contrast imaging | |
CN104535595B (zh) | 一种用于x射线光栅相位衬度成像的背景扣除方法 | |
WO2003040712A1 (en) | Phase-contrast enhanced computed tomography | |
CN106618623A (zh) | 一次曝光的硬x射线光栅干涉仪的成像方法 | |
CN104458777A (zh) | 用于x射线相位衬度成像的积分水桶相位测量方法 | |
Fu et al. | In-line phase contrast micro-CT reconstruction for biomedical specimens | |
CN205924051U (zh) | X射线光栅相衬成像ct*** | |
Lin et al. | Microcomputed tomography | |
Pfeiffer et al. | High-sensitivity phase-contrast tomography of rat brain in phosphate buffered saline | |
Schulz et al. | Imaging the human body: Micro-and nanostructure of human tissues | |
Seifert et al. | Measurement and simulative proof concerning the visibility loss in x-ray Talbot-Lau Moiré imaging | |
Momose et al. | Four-dimensional x-ray phase tomography with Talbot interferometer and white synchrotron light | |
Fu et al. | Single-slice reconstruction method for helical cone-beam differential phase-contrast CT | |
Zhang et al. | Correcting Intensity Drift in X-ray Grating-based Phase Contrast Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
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: 20170329 |