CN106023107B - 一种用于x射线光栅相衬成像装置的探测器图像校正方法 - Google Patents

一种用于x射线光栅相衬成像装置的探测器图像校正方法 Download PDF

Info

Publication number
CN106023107B
CN106023107B CN201610331299.8A CN201610331299A CN106023107B CN 106023107 B CN106023107 B CN 106023107B CN 201610331299 A CN201610331299 A CN 201610331299A CN 106023107 B CN106023107 B CN 106023107B
Authority
CN
China
Prior art keywords
image
grating
detector
gain
correction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610331299.8A
Other languages
English (en)
Other versions
CN106023107A (zh
Inventor
刘刚
胡仁芳
昝贵彬
韩华杰
高昆
陆亚林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201610331299.8A priority Critical patent/CN106023107B/zh
Publication of CN106023107A publication Critical patent/CN106023107A/zh
Application granted granted Critical
Publication of CN106023107B publication Critical patent/CN106023107B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G06T5/77
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/0007Image acquisition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image

Abstract

本发明提供了用于X射线光栅相衬成像装置的探测器图像校正方法,方法首先将源光栅、分束光栅和分析光栅从光路中移除,在X射线管关闭的状态下,采集探测器输出的图像,作为暗电流校正图像Ioffset,然后将源光栅放置于X射线管和探测器之间,打开X射线管,采集探测器输出的图像,作为采集增益校正图像Igain,接着将源光栅、分束光栅和分析光栅放回光路中,并将样品放置于源光栅和分束光栅或者分束光栅和分析光栅之间的光路上,采集探测器输出的图像,作为数据图像Iacquire,最后根据暗电流校正图像Ioffset和增益校正图像Igain对需校正数据图像Iacquire进行暗电流校正和增益校正,得到校正后的数据图像I。

Description

一种用于X射线光栅相衬成像装置的探测器图像校正方法
技术领域
本发明属于辐射成像领域,尤其涉及一种用于X射线光栅相衬成像装置的探测器图像校正方法。
背景技术
自从1895年伦琴教授发现X射线以来,X射线成像技术作为一个最重要的检测手段,已经在生物医学、材料科学、工业检测以及国土安全等众多领域得到了令人瞩目的应用。传统的X射线成像方法是基于X射线对不同物质的吸收差异而实现的,对于如金属、骨骼等由重元素或高密度材料组成的物质,传统的X射线吸收衬度成像技术具有很好的效果。但是对于由碳、氢、氧等轻元素或低密度材料组成的物质,比如医学成像中的乳腺、血管、脂肪以及软骨等人体软组织,传统的X射线成像技术得到的图像衬度极低,分辨能力有限。在用于X射线成像的能力范围内(10-100keV),对于碳、氢、氧等轻元素,X射线的相移系数是各自吸收系数的1000倍以上。因此,利用X射线穿透那些由轻元素组成的物体产生的相位变化(即相移)来进行的X射线相位成像技术能够得到远高于传统吸收方法的图像衬度和测量灵敏度。
X射线相位衬度成像方法,因其对于人体软组织等低原子序数物质相对于传统成像方法有明显优势,自提出以来引起了各方面的广泛关注。从上世纪90年代开始,X射线相衬成像技术主要发展了晶体干涉仪方法、衍射增强方法、同轴方法以及光栅干涉仪方法。由于X射线相移探测要求X射线光源有比较高的相干性,所以刚开始X射线相衬成像技术都是在同步辐射上或者微焦点光源上完成的。在2006年,Pfeiffer等人从可见光的相位测量方法中得到启发,在原有的基于两块光栅的Talbot干涉仪基础上增加了一块源光栅提出了可以利用普通光源的Talbot-Lau干涉仪。由于该方法摆脱了庞大而昂贵的同步辐射光源以及低功率的微焦点光源,真正使得X射线相衬成像应用于医学成像、无损检测领域等成为了可能。Talbot-Lau干涉仪是利用相干光照明产生的Talbot自成像和Lau效应,而同样的在可见光领域有一套利用非相干光照明下光栅投影产生的莫尔条纹进行相位测量的方法。清华大学工程物理系的王振天等人基于上述的几何投影方法,也提出了一套基于非相干光照明的三块光栅组成的相衬成像装置。与Talbot-Lau干涉仪相比,该装置进一步降低了对光源相干性以及光栅的要求,成为了另外一种应用潜力巨大的相衬成像方法。光栅相衬成像方法,其最大特点就是可以同时获取物体的吸收,折射以及暗场图像,三种信息可以反应物质的不同信息,相互补充。
使用常规光源的光栅相衬成像装置和普通的X射线吸收衬度成像装置相比,除了增加了三块X射线透射光栅以外,其也是由普通的X光源、平板探测器以及样品台等组成。对于平板探测器,由于X射线源的不同、接收器内电子线路的不一致性及其正常变化,都会引起不同像素在同样X射线剂量辐射的情况下具有不同的输出信号。这方面的原因主要包括随机噪声、偏置误差、像素响应不一致以及瑕疵像素。因此,为了得到正确和准确的图像,光栅相衬成像得到的数据也需要进行探测器图像校正。
发明内容
(一)要解决的技术问题
本发明提供一种用于X射线光栅相衬成像装置的探测器图像校正方法,能有效地进行随机噪声校正、暗电流校正、增益校正和瑕疵像素校正。
(二)技术方案
本发明提供一种用于X射线光栅相衬成像装置的探测器图像校正方法,X射线光栅相衬成像装置包括X射线管、源光栅、分束光栅、分析光栅和探测器,X射线管产生X射线依次经过源光栅、分束光栅和分析光栅后,在探测器上形成探测器图像,探测器图像校正方法包括:
S1,采集暗电流校正图像Ioffset:将源光栅、分束光栅和分析光栅从光路中移除,在X射线管关闭的状态下,采集探测器输出的图像,作为暗电流校正图像Ioffset
S2,采集增益校正图像Igain:将源光栅放置于X射线管和探测器之间,打开X射线管,采集探测器输出的图像,作为采集增益校正图像Igain
S3,采集数据图像Iacquire:将源光栅、分束光栅和分析光栅放回光路中,并将样品放置于源光栅和分束光栅之间或者分束光栅和分析光栅之间的光路上,采集探测器输出的图像,作为数据图像Iacquire
S4,校正数据图像Iacquire:根据暗电流校正图像Ioffset和增益校正图像Igain对校正数据图像Iacquire进行暗电流校正和增益校正,得到校正后的数据图像I。
进一步,在步骤S4中,采用如下公式进行暗电流校正和增益校正:
Figure BDA0000992123480000031
其中I是校正完毕后的图像,mean表示对图像内所有像素值求平均值。
进一步,方法还包括:S5,校正瑕疵像素:根据亮场图像ILF设定一个灰度区间,若一像素的像素值与一阈值的差值在所述灰度区间外,则判定该像素为瑕疵像素,并将该瑕疵像素的像素值取其邻域中正常像素的灰度平均值,其中,亮场图像ILF=Igain-Ioffset
进一步,采集暗电流校正图像Ioffset时,设定探测器的曝光时间,连续采集多张探测器输出的图像,并对该多张图像进行平均处理后作为暗电流校正图像Ioffset
进一步,采集增益校正图像Igain时,设定探测器的曝光时间、X射线管的电压、电流及焦点,使得探测器输出像素值在探测器最大动态范围的70%,连续采集多张探测器输出的图像,并对该多张图像进行平均处理后作为增益校正图像Igain
进一步,采集数据图像Iacquire时,设定探测器的曝光时间、X射线管的电压、电流及焦点,连续采集多张探测器输出的图像,并对该多张图像进行平均处理后作为数据图像Iacquire
(三)有益效果
本发明提供的用于X射线光栅相衬成像装置的探测器图像校正方法,由于将源光栅放入光源和探测器直接参与了增益校正图像Igain的采集,取得了简化实验操作、消除源光栅对成像光场的影响、提高探测器动态的利用范围的效果。
附图说明
图1为本发明实施例提供的X射线光栅相衬成像装置的示意图;
图2为本发明实施例提供的用于X射线光栅相衬成像装置的探测器图像校正方法的流程图;
图3为本发明实施例提供的探测器图像校正方法中增益校正图像采集的示意图;
图4为本发明实施例提供的亮场图;
图5为本发明实施例提供的实验样品实物图;
图6为本发明实施例提供的本发明得到的样品吸收、折射和暗场图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
图1为本发明实施例提供的X射线光栅相衬成像装置的示意图,如图1所示,装置包括X射线管,源光栅G0,光阑,分束光栅G1,分析光栅G2,探测器、光学精密位移台、光学平台以及控制计算机等零部件。三块光栅均通过光学精密位移台固定于光学平台之上,其中分析光栅与光学精密位移台之间安装有压电陶瓷精密位移电机,用以完成分析光栅的步进运动。X射线源参数设置为:焦点s=1.0mm,电压60KV,电流22.5mA,曝光时间2s。三块光栅参数:源光栅周期p0=120um,占空比1∶2;分束光栅周期p1=60um,占空比1∶2;分析光栅周期p2=120um,占空比1∶1。探测器像素大小200um,总共1024×1024像素,动态范围16bit。源光栅距离光源铍窗z0=0.5cm,与分束光栅距离z1=55cm,分束光栅和分析光栅距离z2=55cm,样品在分束光栅后面10cm,分析光栅紧贴探测器(实际距离闪烁体10cm左右,前面有保护平板探测器的玻璃板)。
图2为本发明实施例提供的用于X射线光栅相衬成像装置的探测器图像校正方法的流程图,如图2所示,方法包括:
S100,采集暗电流校正图像Ioffset,打开探测器预热10分钟之后在X光源关闭时,打开数据采集程序设置单张图像曝光时间为2s,连续采集50张图像进行平均后保存为一张图像作为暗电流校正图像。
S200,采集增益校正图像Igain:如图3所示,将源光栅G0放置在X光源铍窗后0.5cm处,打开X光源设置电压60kV、电流22.5mA、焦点1.0mm,通过对输出图像的观察不断调整电流在12mA,使得探测器输出不出现65556的过曝区域,此时探测器的平均输出为14000左右,然后调整电流至11.2mA使得探测器平均输出在12000左右,然后调整光阑使得探测器上的有效视场为12×12cm左右以能够包含最后的分析光栅,稳定10分钟之后连续采集50张图像进行平均后保存作为增益校正图像。
S300,采集数据图像Iacquire:关闭X射线管,将光学元件移回光路,校准光路之后,即可以进行数据图像Iacquire采集,将样品竖直固定于样品台上,关闭实验铅房防护门,打开X射线管,管电压设置为60KV,初始化并设置探测器参数后,移动样品台至视野中央,控制分析光栅进行相位步进,垂直于分析光栅栅线方向分10步均匀移动分析光栅一个周期,由之前的参数描述可知,即每步步进12um,对于每一步,均采集20张图做平均后保存,记为
Figure BDA0000992123480000051
k=1,2,...,10。样品图像采集完成后,控制样品台电机,将样品移出视野,并将分析光栅移回原位置,按照同样的操作流程,获取背景图像,记为
Figure BDA0000992123480000052
k=1,2,...,10。
S400,校正数据图像Iacquire:根据暗电流校正图像Ioffset和增益校正图像Igain对校正数据图像Iacquire进行暗电流校正和增益校正,得到校正后的数据图像I:
Figure BDA0000992123480000053
其中I是校正完毕后的图像,mean表示对图像内所有像素值求平均值。
S500,校正瑕疵像素:根据亮场图像设定一个灰度区间,若一像素的像素值与一阈值的差值在所述灰度区间外,则判定该像素为瑕疵像素,并将该瑕疵像素的像素值取其领域中正常像素的灰度平均值,其中,增益校正图像Igain减去暗电流校正图像Ioffset得到的就是亮场图像ILF,如图4所示,亮场图像ILF是个均匀的光场,根据亮场图像可以确定瑕疵像素的位置。
S600,计算样品吸收图像、折射图像和暗场图像:按照如下公式获取样品吸收图像:
Figure BDA0000992123480000061
按如下公式获取折射图像:
Figure BDA0000992123480000062
按如下公式获取暗场图像:
Figure BDA0000992123480000063
其中,m,n表示探测器像素的行列位置,In表示求以自然常数e为底的对数,arg表示求复数的幅角,abs表示求复数的膜。
本发明实施例采用的样品如图5所示,由三个有机玻璃的圆柱体组成,从左到右分别是直径2cm的PMMA圆柱棒、直径1cm的POM圆柱棒和直径5mm的PMMA圆柱棒。图6A和图6B是本发明实施例得到的样品吸收、折射和暗场图像,其中图6A是校正前的探测器图像,图6B是校正后的探测器图像,对比图6A和图6B可以看出,在同样灰度范围内,经过校正的吸收、折射和暗场图像比未校正的结果显示出更多细节,正确的反映样品的特征和信息。特别是暗场图像,未校正图像完全看不出样品的特征,不能得到样品正确的信息。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种用于X射线光栅相衬成像装置的探测器图像校正方法,所述X射线光栅相衬成像装置包括X射线管、源光栅、分束光栅、分析光栅和探测器,X射线管产生X射线依次经过源光栅、分束光栅、和分析光栅后,在探测器上形成探测器图像,其特征在于,探测器图像校正方法包括:
S1,采集暗电流校正图像Ioffset:将源光栅、分束光栅和分析光栅从光路中移除,在X射线管关闭的状态下,采集探测器输出的图像,作为暗电流校正图像Ioffset
S2,采集增益校正图像Igain:将源光栅放置于X射线管和探测器之间,打开X射线管,采集探测器输出的图像,作为采集增益校正图像Igain
S3,采集数据图像Iacquire:将源光栅、分束光栅和分析光栅放回光路中,并将样品放置于源光栅和分束光栅之间或者分束光栅和分析光栅之间的光路上,采集探测器输出的图像,作为数据图像Iacquire
S4,校正数据图像Iacquire:根据暗电流校正图像Ioffset和增益校正图像Igain对数据图像Iacquire进行暗电流校正和增益校正,得到校正后的数据图像I;
S5,校正瑕疵像素:根据亮场图像ILF设定一个灰度区间,若一像素的像素值与一阈值的差值在所述灰度区间外,则判定该像素为瑕疵像素,并将该瑕疵像素的像素值取其邻域中正常像素的灰度平均值,其中,亮场图像ILF=Igain-Ioffset
S6,计算样品吸收图像、折射图像和/或暗场图像:按照如下公式获取样品吸收图像:
Figure FDA0002262932440000011
按如下公式获取折射图像:
Figure FDA0002262932440000012
按如下公式获取暗场图像:
Figure FDA0002262932440000021
其中,m,n表示探测器像素的行列位置,In表示求以自然常数e为底的对数,arg表示求复数的幅角,abs表示求复数的膜,N为相位步进总步数,k表示步数,
Figure FDA0002262932440000022
表示相位步进的第k步样品图像强度,
Figure FDA0002262932440000023
表示第k步无样品图像强度,p2表示分析光栅的周期,z2表示分束光栅和分析光栅距离。
2.根据权利要求1所述的用于X射线光栅相衬成像装置的探测器图像校正方法,其特征在于,在所述S4中,采用如下公式进行暗电流校正和增益校正:
Figure FDA0002262932440000024
其中I是校正完毕后的图像,mean表示对图像内所有像素值求平均值。
3.根据权利要求1所述的用于X射线光栅相衬成像装置的探测器图像校正方法,其特征在于,采集暗电流校正图像Ioffset时,设定探测器的曝光时间,连续采集多张探测器输出的图像,并对该多张图像进行平均处理后作为暗电流校正图像Ioffset
4.根据权利要求3所述的用于X射线光栅相衬成像装置的探测器图像校正方法,其特征在于,采集增益校正图像Igain时,设定探测器的曝光时间、X射线管的电压、电流及焦点,使得探测器输出像素值在探测器最大动态范围的70%,连续采集多张探测器输出的图像,并对该多张图像进行平均处理后作为增益校正图像Igain,采集数据图像Iacquire时,设定探测器的曝光时间、X射线管的电压、电流及焦点,连续采集多张探测器输出的图像,并对该多张图像进行平均处理后作为数据图像Iacquire
CN201610331299.8A 2016-05-17 2016-05-17 一种用于x射线光栅相衬成像装置的探测器图像校正方法 Active CN106023107B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610331299.8A CN106023107B (zh) 2016-05-17 2016-05-17 一种用于x射线光栅相衬成像装置的探测器图像校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610331299.8A CN106023107B (zh) 2016-05-17 2016-05-17 一种用于x射线光栅相衬成像装置的探测器图像校正方法

Publications (2)

Publication Number Publication Date
CN106023107A CN106023107A (zh) 2016-10-12
CN106023107B true CN106023107B (zh) 2020-03-27

Family

ID=57098819

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610331299.8A Active CN106023107B (zh) 2016-05-17 2016-05-17 一种用于x射线光栅相衬成像装置的探测器图像校正方法

Country Status (1)

Country Link
CN (1) CN106023107B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109996006B (zh) * 2017-12-29 2021-09-17 江苏康众数字医疗科技股份有限公司 光信号长积分模式成像方法及成像系统
US20240046629A1 (en) * 2018-08-08 2024-02-08 Rensselaer Polytechnic Institute Enhancing contrast sensitivity and resolution in a grating interferometer by machine learning
CN110133012B (zh) * 2019-07-02 2022-01-18 合肥工业大学 基于三探测器光栅干涉仪的单次曝光多模式x射线成像方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104132953A (zh) * 2014-08-01 2014-11-05 中国科学技术大学 一种双能x射线相位衬度成像装置及其实现方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104132953A (zh) * 2014-08-01 2014-11-05 中国科学技术大学 一种双能x射线相位衬度成像装置及其实现方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
X射线平板探测器数字成像及其图像校准;周正干 等;《北京航空航天大学学报》;20040831;第30卷(第8期);第3节 *
基于变能量成像的能量参数预测;毕䶮;《中国优秀硕士学位论文全文数据库信息科技辑》;20150715(第07期);第16页 *
平板探测器成像系统的成像条件优化;安宏庆;《陕西省第十三届无损检测年会》;20121231;237-244 *

Also Published As

Publication number Publication date
CN106023107A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
US8523434B2 (en) Image processing apparatus and image processing method for removing residual image
JP2525162B2 (ja) X線撮影装置
JP6133231B2 (ja) X線エネルギースペクトル測定方法およびx線エネルギースペクトル測定装置およびx線ct装置
JP7042806B2 (ja) スペクトルコンピュータ断層撮影法(ct)のスペクトル較正
JP3560374B2 (ja) X線撮像における散乱照射線補正のために部分的に透過性である物体で遮蔽する方法
US10176601B2 (en) Methods and systems for improving image quality of projection image data acquired using flat panel detectors
CN106023107B (zh) 一种用于x射线光栅相衬成像装置的探测器图像校正方法
WO2009142166A1 (ja) X線診断装置
JP2012120652A (ja) 放射線撮影システム
US10779776B2 (en) Apparatus for X-ray imaging an object
Liu et al. Spectral response model for a multibin photon-counting spectral computed tomography detector and its applications
WO2012164901A1 (ja) 放射線被曝量取得方法および装置並びに放射線画像撮影システム
JP2012115577A (ja) 放射線撮影システム
JPH10234724A (ja) X線ct装置
Schmid et al. Modeling vibrations of a tiled talbot-lau interferometer on a clinical CT
JPH0614911A (ja) X線診断方法及びその装置
US8023620B2 (en) X-ray arrangement with a converter and associated X-ray method
CN106404809A (zh) 一种用于x射线光栅相衬成像装置的图像校正方法
US7949174B2 (en) System and method for calibrating an X-ray detector
JP2020514776A (ja) 放射線医学システムによって投与された放射線量を推定する方法
Meng et al. Development of a dual-detector X-ray imaging system for phase retrieval study
Park et al. Effects of image processing on the detective quantum efficiency
JP2012120650A (ja) 放射線撮影システム及び放射線位相コントラスト画像生成方法
Hauke et al. Construction and evaluation of a high-energy grating-based x-ray phase-contrast imaging setup
Schmid et al. Modeling vibrations of a tiled Talbot-Lau interferometer on a clinical CT system

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