CN105852895B - 单次曝光的硬x射线光栅干涉仪的信息提取方法 - Google Patents
单次曝光的硬x射线光栅干涉仪的信息提取方法 Download PDFInfo
- Publication number
- CN105852895B CN105852895B CN201610291248.7A CN201610291248A CN105852895B CN 105852895 B CN105852895 B CN 105852895B CN 201610291248 A CN201610291248 A CN 201610291248A CN 105852895 B CN105852895 B CN 105852895B
- Authority
- CN
- China
- Prior art keywords
- grating
- ray
- hard
- imaged
- formula
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000010521 absorption reaction Methods 0.000 claims abstract description 32
- 230000010363 phase shift Effects 0.000 claims abstract description 27
- 230000011664 signaling Effects 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 5
- 230000003287 optical effect Effects 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000000605 extraction Methods 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 7
- 239000003574 free electron Substances 0.000 claims description 4
- 241000216843 Ursus arctos horribilis Species 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 32
- 230000005855 radiation Effects 0.000 abstract description 11
- 239000003814 drug Substances 0.000 abstract description 8
- 239000000284 extract Substances 0.000 abstract description 7
- 210000004872 soft tissue Anatomy 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000002591 computed tomography Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229920005439 Perspex® Polymers 0.000 description 2
- 239000004411 aluminium Substances 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 2
- 238000009659 non-destructive testing Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 239000004926 polymethyl methacrylate Substances 0.000 description 2
- 208000019155 Radiation injury Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000000333 X-ray scattering Methods 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 239000002861 polymer material Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 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
-
- 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]
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- 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/44—Constructional features of apparatus for radiation diagnosis
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
-
- 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/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- 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/54—Control of apparatus or devices for radiation diagnosis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种单次曝光的硬X射线光栅干涉仪的信息提取方法，其特征是按如下步骤进行：1移动任一光栅，使硬X射线光栅干涉仪工作在光强变化曲线的左半腰或右半腰位置处；2分别获取背景投影图像和被成像物体的投影图像；3获取归一化的被成像物体的投影图像，对其进行取对数处理；4根据对数结果，构建表达式进行傅立叶变换；5根据傅里叶变换结果，利用逆傅里叶变换，分别提取被成像物体的吸收和相移信号。本发明摒弃了繁琐的光栅扫描流程，简化了硬X射线光栅干涉仪的图像获取流程，能够实现快速、低辐射剂量的硬X射线相位衬度成像，从而提高了成像效率，为低辐射剂量临床医学成像提供新途径。
Description
技术领域
本发明涉及X射线成像物理领域，具体的说是一种单次曝光的硬X射线光栅干涉仪的信息提取方法。
背景技术
传统X射线成像技术，包括投影成像术(Radiography)和计算机断层成像技术(Computed Tomography，简称CT)，是现代科学技术对人类健康的巨大贡献。一百年以来，该项技术在投影数据采集和重建理论等方面不断取得新的进展，逐渐发展成为生命医学、材料科学、考古科学、工业无损检测等领域不可或缺的重要工具。
传统X射线成像技术的物理基础是物体对X射线吸收性质的差异。因此，在对金属、骨骼等重元素组成的物体进行成像时，能够获得很高的图像衬度，得到清晰的图像。但是，临床医学成像中的人体软组织、无损检测中的聚合物材料等物体，主要由碳、氢、氧等轻元素组成，它们对X射线的吸收非常弱，以至于利用现有X成像技术无法有效观察它们。因此，有必要发展新的X射线成像技术和方法来实现对弱吸收物体的成像。
科学家们发现，在硬X射线波段(10-100keV)，对碳、氢、氧等软组织的主要组成元素，其折射率相移项是其吸收项的1000多倍。因此，对人体软组织等弱吸收物质来说，测量X射线穿过物体时的相移信息要比检测衰减信息有效得多。X射线相位衬度成像技术正是通过记录X射线穿过物体后相位的改变量(即相移)来形成图像衬度的。与传统的吸收衬度成像技术相比，X射线相位衬度成像技术能够获得更高的图像衬度、更高的测量灵敏度和更低的辐射剂量。20世纪90年代中期以来，基于X射线相位衬度，陆续发展了四种投影成像方法，包括晶体干涉仪、光栅干涉仪、衍射增强成像、相位传播成像。这些方法在对生物软组织等轻元素样品成像时，都获得了比传统X射线成像技术高得多的图像衬度。
其中，硬X射线光栅干涉仪是21世纪初才发展起来的，目前被认为是最有发展潜力的X射线相位衬度成像方法之一。与其他成像方法相比，硬X射线光栅干涉仪是目前唯一能够有效利用常规X射线源获得相位衬度的方法，因而在临床医学诊疗、材料科学等领域具有广阔的应用前景。硬X射线光栅干涉仪具有很高的测量灵敏度，能够从一组投影数据中同时获得物体的吸收、相移和散射图像。三种不同的物体图像互为补充，从不同角度反映了物体的内部结构信息。
在硬X射线光栅干涉仪中，一方面，探测器记录的投影图像中同时包含物体了的吸收、相移和散射信号。另一方面，定量分析和研究物体内部结构、计算机三维断层重建等均要求独立、纯粹的吸收、相移和散射信号。因此，物体信息提取是不可或缺的中间步骤。必须发展相应信息提取方法，能够准确、定量地从实验投影图像中提取物体的吸收、相位和散射信息。
目前，硬X射线光栅干涉仪普遍采用相位步进法进行实验数据采集和物体信息提取。这种方法要求在横向扫描光栅的同时，采集多幅投影图像(实验上不少于4张)。这样不仅导致很长的图像数据采集时间，更严重的是，多次曝光增加了物体的辐射剂量和辐射损伤的几率。这与X射线成像设备快速、低辐射剂量的发展趋势是矛盾的，阻碍了硬X射线光栅干涉仪方法在临床医学诊断等领域的推广应用。因此，发展新的信息提取方法，克服相位步进法光栅扫描、多次物体曝光的局限性，已经发展成为硬X射线光栅干涉仪实用化进程中必须解决的瓶颈问题之一。
发明内容
本发明为避免现有技术存在的不足之处，提供一种单次曝光的硬X射线光栅干涉仪的信息提取方法，摒弃了繁琐的光栅扫描流程，简化了硬X射线光栅干涉仪的图像获取流程，能够实现快速、低辐射剂量的硬X射线相位衬度成像，从而提高成像效率，为低辐射剂量临床医学成像提供新途径。
为解决上述技术问题，本发明采用的技术方案是：
本发明种单次曝光的硬X射线光栅干涉仪的信息提取方法，所述硬X射线光栅干涉仪包括：X射线源、源光栅G0、分束光栅G1、分析光栅G2、探测器；其特点是，所述信息提取方法按如下步骤进行：
步骤1、固定所述源光栅G0、分束光栅G1和分析光栅G2中的任意两个光栅，并沿着移动方向将另一个光栅移动自身四分之一的光栅周期，使得所述硬X射线光栅干涉仪工作在光强变化曲线的左半腰或右半腰位置处；所述移动方向为同时垂直于光轴和光栅栅条的方向；
步骤2、启动所述X射线源和所述探测器，设置曝光时间为t；
利用所述探测器按照所述曝光时间t获取初始背景的投影图像I0后，关闭所述X射线源；
步骤3、将被成像物体放置到所述分束光栅G1的视场中央后，启动所述X射线源，并利用所述探测器按照所述曝光时间t获取所述被成像物体的投影图像ID后，关闭所述X射线源；
步骤4、对所述被成像物体的投影图像ID进行背景归一化处理，得到归一化投影图像I′D，其中I′D＝ID/I0；
步骤5、对归一化投影图像I′D进行取对数处理，得到处理结果ln(I′D)；
步骤6、对表达式K1-ln(I′D)作一维傅里叶变换，得到变换结果F；K1为常数，且满足K1＝lnS(p2/4)或K1＝lnS(-p2/4)，其中，p2是所述分析光栅G2的光栅周期；S(p2/4)为所述硬X射线光栅干涉仪的光强变化曲线在p2/4处的数值；
步骤7、根据变换结果F，利用式(1)得到所述被成像物体的电荷密度投影的一维傅里叶变换结果F1：
F1＝F/(σKN+iK2λ2reu/M) (1)
式(1)中，σKN表示单个自由电子的康普顿散射截面，i表示虚数单位，i2＝-1，K2为常数，且满足或其中S′(p2/4)表示所述硬X射线光栅干涉仪光强变化曲线的一阶空间导数在p2/4处的数值，dT是所述分束光栅G1到所述分析光栅G2之间的轴向距离，λ表示所述X射线源的等效波长，re是经典电子半径，u表示空间频率，M是所述硬X射线光栅干涉仪的几何放大率；
步骤8、利用式(2)提取所述被成像物体的吸收信号T：
式(2)中，表示一维逆傅里叶变换；
步骤9、利用式(3)提取所述被成像物体的相移信号Φ：
以所述被成像物体的吸收信号T和相移信号Φ作为信息提取结果。
与已有技术相比，本发明的有益效果是：
1、本发明基于硬X射线波段，人体乳腺等软组织的相移-衰减二象性(Phase-Attenuation Duality)的物理基础，提出了单次曝光物体信息提取新方法，解决了单次曝光下物体吸收、相移信息的定量提取问题，克服了现有相位步进法的光栅步进扫描、多次物体曝光、高辐射剂量的局限性，实现了快速、低辐射剂量的硬X射线相位衬度成像。
2、与现有的相位步进法相比，本发明在获取物体图像图像时，通过将硬X射线光栅干涉仪工作在左半腰或右半腰位置处，避免了繁琐的光栅步进扫描，极大地简化了实验步骤，降低了硬X射线光栅干涉仪对机械稳定性的要求；
3、与现有的相位步进法相比，本发明只要求对物体进行一次曝光，避免了对物体的多次(实验上至少4次)曝光，减少了辐射剂量，降低了辐射损伤风险，更加安全可靠；
4、与现有的相位步进法相比，本发明只需要获取1张投影图像，就能够同时提取物体的吸收和相移信号，避免了获取多张(实验上至少4张)投影图像，减少了实验时间，提高了实验效率；
5、本发明摒弃了繁琐的光栅步进扫描，每个投影角度下只需要获取1张物体的投影图像，使得能够与计算机断层扫描技术直接结合，获取物体的线性吸收系数、相移系数的三维断层重建。
6、与传统X射线成像方法相比，本发明通过获取物体的相移信号，能够更加清晰地观察人体软组织等的精细结构，更加准确区别正常与病变组织。
7、本发明提出的物体信息提取新方法同样适用于二维硬X射线光栅干涉仪、中子光栅干涉仪。
附图说明
图1为现有技术中以大焦点源与源光栅G0的组合为硬X射线源的一维光栅干涉仪示意图；
图2为现有技术中硬X射线光栅干涉仪的光强变化曲线图；
图3为现有技术中有机玻璃的硬X射线散射截面图；
图4为本发明物体吸收信号的提取结果图；
图5为本发明物体相移信号的提取结果图。
具体实施方式
本实施例中，如图1所示，硬X射线光栅干涉仪包括：X射线源、源光栅G0、分束光栅G1、分析光栅G2、探测器；一种单次曝光的硬X射线光栅干涉仪的信息提取方法是按如下步骤进行：
步骤1、固定源光栅G0、分束光栅G1和分析光栅G2中任意两个光栅，并沿着移动方向将另一个光栅移动自身四分之一的光栅周期，使得硬X射线光栅干涉仪工作在光强变化曲线的左半腰或右半腰位置处(图2中箭头指示位置)；移动方向为同时垂直于光轴和光栅栅条的方向；
硬X射线光栅干涉仪工作在光强变化曲线的左半腰或右半腰位置处，使得被成像物体的相移信号的测量灵敏度达到最大，能够更加准确地提取被成像物体的相移信号。
步骤2、启动X射线源和探测器，设置曝光时间为t；
利用探测器按照曝光时间t获取初始背景的投影图像I0后，用于被成像物体的投影图像的归一化处理，关闭X射线源；
步骤3、将被成像物体放置到分束光栅G1的视场中央后，启动X射线源，并利用探测器按照曝光时间t获取被成像物体的投影图像ID后，关闭X射线源；
以硬X射线光栅干涉仪工作在右半腰p2/4为例。探测器获取的被成像物体的投影图像ID满足：
式(3.1)中，T是待提取的被成像物体的吸收信号；S()表示硬X射线光栅干涉仪的光强变化曲线；p2是分析光栅G2的光栅周期；λ表示X射线源的等效波长；dT是分束光栅G1到分析光栅G2之间的轴向距离；M是硬X射线光栅干涉仪的几何放大率；表示一阶空间微分运算；Φ是待提取的被成像物体的相移信号。
步骤4、对被成像物体的投影图像ID进行背景归一化处理，得到归一化投影图像I′D，其中I′D＝ID/I0；
步骤5、对归一化投影图像I′D进行取对数处理，得到处理结果ln(I′D)；
以硬X射线光栅干涉仪工作在右半腰p2/4为例。根据式(3.1)，归一化投影图像I′D的对数处理结果ln(I′D)满足：
在临床医学成像等应用领域，人体软组织等被成像物体对高能硬X射线的折射信号(正比于被成像物体的相移信号的一阶空间导数)是非常小的，远小于p2/4dT；步骤1将硬X射线光栅干涉仪工作在光强变化曲线的左半腰或右半腰位置处，对于小折射信号，能够合理地对光强变化曲线作一阶线性近似，使得非线性成像方程简化为线性成像方程，
式(5.2)代入式(5.1)可得到，
式(5.3)中K1为常数，且满足K1＝lnS(p2/4)，其中p2是分析光栅G2的光栅周期；S(p2/4)为硬X射线光栅干涉仪的光强变化曲线在p2/4处的数值；K2为常数，且满足
其中S′(p2/4)表示硬X射线光栅干涉仪光强变化曲线的一阶空间导数在p2/4处的数值，dT是分束光栅G1到分析光栅G2之间的轴向距离，λ表示X射线源的等效波长，表示一阶空间微分运算。
T是待提取的被成像物体的吸收信号。在临床医学成像等领域感兴趣的硬X射线能量范围，被成像物体(如乳腺等人体软组织)的吸收信号主要来源于光电吸收截面和康普顿散射截面。对于高能硬X射线(光子能量不低于60keV)，光电吸收截面对被成像物体的吸收信号的贡献非常小，完全可以被忽略。以人体软组织的等效替代物—有机玻璃为被成像物体。如图3所示，随着X射线能量的提高，光电吸收截面迅速衰减。对于60keV的硬X射线，康普顿散射截面对被成像物体的吸收信号的贡献超过了95％。因此，被成像物体的吸收信号T能够合理地近似表达为：
式(5.4)中σKN表示单个自由电子的康普顿散射截面，可利用Klein-Nishina公式得到，
式(5.5)中η＝E/511keV，其中E表示硬X射线的等效光子能量，re是经典电子半径，re＝2.8*10^-15米；ρe是被成像物体的电荷密度。
Φ是待提取的被成像物体的相移信号。当X射线能量远离被成像物体的组成元素的吸收边时(硬X射线成像时的普遍情形)，相移信号来自于相干散射效应，可表达为，
Φ(x,y)＝λre∫ρe(x,y,z)dz (5.6)
式(5.4)、(5.6)代入式(5.3)得到，
步骤6、对表达式K1-ln(I′D)作一维傅里叶变换，得到变换结果F；K1为常数，且满足K1＝lnS(p2/4)或K1＝lnS(-p2/4)，其中，p2是分析光栅G2的光栅周期；S(p2/4)为X射线光栅干涉仪的光强变化曲线在p2/4处的数值；
对式(5.7)作一维傅里叶变换，得到变换结果F：
式(6.1)中，表示一维傅里叶变换，F1是被成像物体的电荷密度投影的一维傅里叶变换结果，i表示虚数单位，i2＝-1，u表示空间频率，决定于探测器的像素尺寸。
步骤7、根据变换结果F，利用式(1)得到被成像物体的电荷密度投影的一维傅里叶变换结果F1：
F1＝F/(σKN+iK2λ2reu/M) (1)
式(1)中，σKN表示单个自由电子的康普顿散射截面，i表示虚数单位，i2＝-1，K2为常数，且满足或其中S′(p2/4)表示硬X射线光栅干涉仪光强变化曲线的一阶空间导数在p2/4处的数值，dT是分束光栅G1到分析光栅G2之间的轴向距离，λ表示X射线源的等效波长，re是经典电子半径，u表示空间频率，M是硬X射线光栅干涉仪的几何放大率；
步骤8、利用式(2)提取被成像物体的吸收信号T：
式(2)中，表示一维逆傅里叶变换；
根据式(5.4)，被成像物体的吸收信号的一维傅里叶变换结果FT为：
将式(1)代入式(8.1)得到，
对式(8.2)作一维逆傅里叶变换，提取被成像物体的吸收信号T：
图4为物体吸收信号的提取结果图，被成像物体是直径2毫米的铝棒(左侧)和直径4毫米的有机玻璃棒(右侧)。入射X射线等效能量为60keV。根据图4，利用式(2)提取的吸收信号与理论值符合的很好，证实了本发明能够有效提取被成像物体的吸收信号。
步骤9、利用式(3)提取被成像物体的相移信号Φ：
根据式(5.6)，被成像物体的相移信号Φ的一维傅里叶变换结果FΦ满足：
将式(1)代入式(9.1)得到：
对式(9.2)作一维逆傅里叶变换，提取被成像物体的相移信号Φ：
图5为物体相移信号的提取结果图，被成像物体是直径2毫米的铝棒(左侧)和直径4毫米的有机玻璃棒(右侧)。入射X射线等效能量为60keV。根据图5，利用式(3)提取的相移信号与理论值符合的很好，证实了本发明能够有效提取被成像物体的相移信号。尽管在被成像物体的边界位置，提取结果值相比理论值有一些偏差。其原因是式(3)中分母上的低通滤波器造成被成像物体边界处高频信息的丢失。这些偏差可通过进一步减小探测器像素尺寸等途径来减小。这并不影响本发明提出的信息提取新方法的可行性。
以被成像物体的吸收信号T和相移信号Φ作为信息提取结果。
Claims (1)
1.一种单次曝光的硬X射线光栅干涉仪的信息提取方法，所述硬X射线光栅干涉仪包括：X射线源、源光栅G0、分束光栅G1、分析光栅G2、探测器；其特征是，所述信息提取方法按如下步骤进行：
步骤1、固定所述源光栅G0、分束光栅G1和分析光栅G2中的任意两个光栅，并沿着移动方向将另一个光栅移动自身四分之一的光栅周期，使得所述硬X射线光栅干涉仪工作在光强变化曲线的左半腰或右半腰位置处；所述移动方向为同时垂直于光轴和光栅栅条的方向；
步骤2、启动所述X射线源和所述探测器，设置曝光时间为t；
利用所述探测器按照所述曝光时间t获取初始背景的投影图像I0后，关闭所述X射线源；
步骤3、将被成像物体放置到所述分束光栅G1的视场中央后，启动所述X射线源，并利用所述探测器按照所述曝光时间t获取所述被成像物体的投影图像ID后，关闭所述X射线源；
步骤4、对所述被成像物体的投影图像ID进行背景归一化处理，得到归一化投影图像I′D，其中I′D＝ID/I0；
步骤5、对归一化投影图像I′D进行取对数处理，得到处理结果ln(I′D)；
步骤6、对表达式K1-ln(I′D)作一维傅里叶变换，得到变换结果F；K1为常数，且满足K1＝lnS(p2/4)或K1＝lnS(-p2/4)，其中，p2是所述分析光栅G2的光栅周期；S(p2/4)为所述硬X射线光栅干涉仪的光强变化曲线在p2/4处的数值；
步骤7、根据变换结果F，利用式(1)得到所述被成像物体的电荷密度投影的一维傅里叶变换结果F1：
F1＝F/(σKN+iK2λ2reu/M) (1)
式(1)中，σKN表示单个自由电子的康普顿散射截面，i表示虚数单位，i2＝-1，K2为常数，且满足或其中S′(p2/4)表示所述硬X射线光栅干涉仪光强变化曲线的一阶空间导数在p2/4处的数值，dT是所述分束光栅G1到所述分析光栅G2之间的轴向距离，λ表示所述X射线源的等效波长，re是经典电子半径，u表示空间频率，M是所述硬X射线光栅干涉仪的几何放大率；
步骤8、利用式(2)提取所述被成像物体的吸收信号T：
式(2)中，表示一维逆傅里叶变换；
步骤9、利用式(3)提取所述被成像物体的相移信号Φ：
以所述被成像物体的吸收信号T和相移信号Φ作为信息提取结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610291248.7A CN105852895B (zh) | 2016-04-29 | 2016-04-29 | 单次曝光的硬x射线光栅干涉仪的信息提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610291248.7A CN105852895B (zh) | 2016-04-29 | 2016-04-29 | 单次曝光的硬x射线光栅干涉仪的信息提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105852895A CN105852895A (zh) | 2016-08-17 |
| CN105852895B true CN105852895B (zh) | 2018-07-31 |
Family
ID=56631118
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610291248.7A Active CN105852895B (zh) | 2016-04-29 | 2016-04-29 | 单次曝光的硬x射线光栅干涉仪的信息提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105852895B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106618623B (zh) * | 2017-01-11 | 2019-08-30 | 合肥工业大学 | 一次曝光的硬x射线光栅干涉仪的成像方法 |
| CN107290360B (zh) * | 2017-06-22 | 2019-07-19 | 合肥工业大学 | 双波长中子光栅干涉仪的成像方法 |
| CN110133012B (zh) * | 2019-07-02 | 2022-01-18 | 合肥工业大学 | 基于三探测器光栅干涉仪的单次曝光多模式x射线成像方法 |
| CN111595877B (zh) * | 2020-05-27 | 2022-03-29 | 合肥工业大学 | 一种x射线衍射增强成像的多衬度图像提取方法 |
| CN112415030B (zh) * | 2020-11-18 | 2022-02-15 | 首都师范大学 | 一种x射线微分相移ct的感兴趣区域重建方法 |
| CN113063809B (zh) * | 2021-03-24 | 2022-05-10 | 合肥工业大学 | 一种基于霍夫变换法的x射线光栅干涉仪成像方法 |
| CN114886445B (zh) * | 2022-07-15 | 2022-12-13 | 康达洲际医疗器械有限公司 | 一种基于多叶光栅动态可调的双c臂三维成像方法与*** |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101467889B (zh) * | 2007-12-26 | 2010-08-25 | 中国科学院高能物理研究所 | 光栅剪切相位衬度ct成像数据采集和重建方法 |
| CN101943668B (zh) * | 2009-07-07 | 2013-03-27 | 清华大学 | X射线暗场成像***和方法 |
| CN103188996B (zh) * | 2010-10-29 | 2015-06-24 | 富士胶片株式会社 | 放射线照相相衬成像设备 |
| US20130259194A1 (en) * | 2012-03-30 | 2013-10-03 | Kwok L. Yip | Hybrid slot-scanning grating-based differential phase contrast imaging system for medical radiographic imaging |
| CN104869905B (zh) * | 2012-12-21 | 2019-08-06 | 卡尔斯特里姆保健公司 | 基于微分相衬成像的医疗放射照相光栅 |
| EP2978377B1 (en) * | 2013-03-26 | 2021-05-05 | Institute of Experimental and Applied Physics | Method of phase gradient radiography and arrangement of an imaging system for application of the method |
| WO2015168473A1 (en) * | 2014-05-01 | 2015-11-05 | Sigray, Inc. | X-ray interferometric imaging system |
| CN104132953B (zh) * | 2014-08-01 | 2017-03-29 | 中国科学技术大学 | 一种双能x射线相位衬度成像装置及其实现方法 |
-
2016
- 2016-04-29 CN CN201610291248.7A patent/CN105852895B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105852895A (zh) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105852895B (zh) | 单次曝光的硬x射线光栅干涉仪的信息提取方法 | |
| US8121249B2 (en) | Multi-parameter X-ray computed tomography | |
| JP6529968B2 (ja) | Ｘ線格子位相コントラストイメージング装置及び方法 | |
| US9494534B2 (en) | Material differentiation with phase contrast imaging | |
| Schulz et al. | High-resolution tomographic imaging of a human cerebellum: comparison of absorption and grating-based phase contrast | |
| WO2016070771A1 (zh) | X射线相衬成像***与成像方法 | |
| EP3548876B1 (de) | Dunkelfeldtensortomographieverfahren | |
| CN104569002B (zh) | 基于光子计数的x射线相衬成像***、方法及其设备 | |
| CN106575533B (zh) | 用于x射线成像设备的光栅设备 | |
| Lang et al. | Experimental comparison of grating-and propagation-based hard X-ray phase tomography of soft tissue | |
| CN106618623B (zh) | 一次曝光的硬x射线光栅干涉仪的成像方法 | |
| Ando et al. | X-ray dark-field phase-contrast imaging: Origins of the concept to practical implementation and applications | |
| Qiu et al. | Internal elemental imaging by scanning X-ray fluorescence microtomography at the hard X-ray microprobe beamline of the SSRF: Preliminary experimental results | |
| Wagermaier et al. | Understanding Hierarchy and Functions of Bone Using Scanning X‐ray Scattering Methods | |
| Castro et al. | Coherent scattering X-ray imaging at the Brazilian National Synchrotron Laboratory: Preliminary breast images | |
| Rawlik et al. | Refraction beats attenuation in breast CT | |
| Huang et al. | Strategy of extraction methods and reconstruction algorithms in computed tomography of diffraction enhanced imaging | |
| Zhu et al. | Three‐dimensional reciprocal space x‐ray coherent scattering tomography of two‐dimensional object | |
| Xue et al. | Microscopic identification of Chinese medicinal materials based on X-ray phase contrast imaging: from qualitative to quantitative | |
| Liu et al. | Phase-retrieval-based synchrotron X-ray micro-tomography for 3D structural characterization and quantitative analysis of agalloch eaglewood | |
| Li et al. | Quantitative characterization of ex vivo breast tissue via x-ray phase-contrast tomography | |
| Liu et al. | Quantitative characterization and diagnosis via hard X‐ray phase‐contrast microtomography | |
| Endrizzi et al. | X-ray phase-contrast tomography of breast tissue specimen with a multi-aperture analyser synchrotron set-up | |
| Dreossi et al. | Synchrotron radiation micro-tomography: a non-invasive tool for the characterization of archaeological wood | |
| Daskalaki et al. | Image quality evaluation of phase contrast mammographic techniques |
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 |