CN103365067B - 可实现三维动态观测的光栅剪切成像装置和方法 - Google Patents
可实现三维动态观测的光栅剪切成像装置和方法 Download PDFInfo
- Publication number
- CN103365067B CN103365067B CN201310110905.XA CN201310110905A CN103365067B CN 103365067 B CN103365067 B CN 103365067B CN 201310110905 A CN201310110905 A CN 201310110905A CN 103365067 B CN103365067 B CN 103365067B
- Authority
- CN
- China
- Prior art keywords
- grating
- array
- projection
- sample
- data
- 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 - Fee Related
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种可实现三维动态观测的光栅剪切成像装置和方法,装置包括:样品台;光源阵列和光源光栅阵列,均排列在以样品台为圆心的半个第一圆周上、且光源阵列贴近光源光栅阵列设置;探测器阵列和分析光栅阵列,均排列在以样品台为圆心的另外半个第一圆周上、且探测器阵列贴近分析光栅阵列设置;分束光栅阵列排列在以样品台为圆心靠近光源光栅阵列或靠近分析光栅阵列的半个第二圆周上,第一圆周的半径大于第二圆周;在每条连接一光源和一探测器的直线上依次设置的各部件均构成一套成像子系统。本发明提高了图像采集速度,降低样品辐射剂量,有利于对样品内部结构的局部区域实现三维动态观测,可满足医学成像检测等应用需求。
Description
技术领域
本发明涉及成像技术领域,特别是涉及一种可实现三维动态观测的光栅剪切成像装置和方法。
背景技术
传统的X射线投影成像术(Radiography)和由此发展起来的计算机断层成像术(Computed Tomography,简称CT)是现代科学技术对人类健康的巨大贡献。一百年以来,该项技术在投影数据采集和重建理论等方面不断地取得新的进展,成为医学、生物、考古、材料科学、安全检查和工业无损检测等领域不可或缺的重要工具。
传统的X射线投影成像术的物理基础是物质对X射线吸收的差别,基于吸收机制的X射线投影成像技术仅对人体骨骼观察的比较清楚,而对人体软组织成像模糊,特别是难以观察到人体内部毫米量级的早期恶性肿瘤(如癌症等)。这其中的原因在于人体软组织主要是由轻元素组成的,轻元素对于硬X射线几乎没有吸收,就象可见光透过水中的一个小玻璃球,几乎没有留下可以察觉的痕迹。因而需要发展一种新方法,能清楚地分辨癌组织和正常组织,这就好象要想办法看清楚水中的一个小玻璃球。科学家发现,X射线经过轻元素物质时,相位变化是吸收变化的一千倍到十万倍,因此X射线相位衬度成像可以提供比吸收衬度成像敏感得多的信号,特别适合观察轻元素构成的物体。
近二十年来,基于X射线相位衬度,发展了四种投影成像方法,分别是晶体干涉仪成像方法、相位传播成像方法、衍射增强成像方法和光栅剪切成像方法。其中,光栅剪切成像方法虽然是本世纪才发展起来,却最有发展前途。它不但能获得样品的折射角图像,重建出样品折射率实部衰减率和折射率实部衰减率导数的三维分布,而且具有和常规X射线光源结合、向医学临床应用发展的前景。除此之外,光栅剪切成像方法中还能利用样品中的散射信号,对样品中的粉末、泡沫、多孔或多晶结构进行成像,并重建出这些结构线性消光系数的三维分布。
目前利用光栅扫描提取相位信息和散射信息是国际上发展的主流,然而,光栅扫描的方法不符合医学成像简便快速的要求。
发明内容
本发明提供一种可实现三维动态观测的光栅剪切成像装置和方法。
一方面,本发明提供了一种可实现三维动态观测的光栅剪切成像装置,包括:
样品台、包括多个光源的光源阵列、包括多个光源光栅的光源光栅阵列、包括多个分束光栅的分束光栅阵列、包括多个分析光栅的分析光栅阵列以及包括多个探测器的探测器阵列;
所述光源阵列和所述光源光栅阵列均排列在以所述样品台为圆心的半个第一圆周上、且所述光源阵列贴近所述光源光栅阵列设置;所述探测器阵列和所述分析光栅阵列均排列在以所述样品台为圆心的另外半个第一圆周上、且所述探测器阵列贴近所述分析光栅阵列设置;所述分束光栅阵列排列在以所述样品台为圆心靠近所述光源光栅阵列的半个第二圆周上,或所述分束光栅阵列排列在以所述样品台为圆心靠近所述分析光栅阵列的半个第二圆周上,所述第一圆周的半径大于所述第二圆周;
在每条连接一光源和一探测器的直线上依次设置的光源、光源光栅、分束光栅、分析光栅和探测器,均构成一套成像子系统。
另一方面,本发明还提供了一种基于上述可实现三维动态观测的光栅剪切成像装置的光栅剪切成像方法,在各组成像子系统中分别执行以下步骤:
调整光源,使光源发射出的光束对准沿连接该光源到与其对称设置的探测器的直线上设置的光源光栅、分束光栅、样品台、分析光栅和探测器,并沿第一圆周或第二圆周的径向垂直方向和径向平行方向分别调整光源光栅、分束光栅和分析光栅位置,使分束光栅的一维多缝产生一维光束阵列,使分析光栅的一维多缝对准所述一维光束阵列;
测量位移曲线,沿第一圆周或第二圆周的径向垂直方向和径向平行方向分别调整光源光栅、分束光栅和分析光栅位置,以调整所述分析光栅和所述一维光束阵列之间发生的剪切位移,通过探测器的各探测单元均采集光强随剪切位移变化的位移曲线,并用三角函数曲线拟合光强随剪切位移变化的位移曲线;
通过探测器探测光强背景和空间位置变化,沿第一圆周或第二圆周的切向调整光源光栅、分束光栅或分析光栅位置,直至把分析光栅和相应的所述一维光束阵列之间的剪切位移调整在光强背景满足成像要求的采集位置;
通过所述探测器采集固定在样品台上的样品对应所述光强背景下的投影像,并将探测器阵列同时采集的多个角度的投影像构成一套投影数据。
本发明提供的技术方案,无须转动样品,即可同时对样品各个视角进行成像,图像采集的速度快,基于光栅剪切原理实现相位衬度成像,基于计算机断层成像原理实现三维成像,可以对样品的内部结构实现三维动态观测,可满足医学检测等方面的应用需求。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1(a)和图1(b)为本发明提供的可实现三维动态观测的光栅剪切成像装置的可选结构示意图;在图1(a)中的光栅剪切成像装置中,分束光栅阵列排列在样品台和光源光栅阵列之间的半个第二圆周上,在图1(b)中的光栅剪切成像装置中,分束光栅阵列排列在样品台和分析光栅阵列之间的半个第二圆周上;
图2为本发明一实施例中样品对X射线光束吸收衰减作用的示意图,其中I0为入射光,I为出射光;
图3为本发明一实施例中样品对X射线光束产生折射作用的示意图;
图4为本发明一实施例中样品对X射线光束产生散射作用的示意图;
图5为本发明一实施例中光强随分析光栅(四条黑色)相对分束光栅产生的一维光束阵列(条纹填充)沿X轴剪切位移变化的位移曲线,(从左到右)分析光栅(四条黑色)和分束光栅一维光束阵列(条纹填充)之间剪切位移分别固定在暗场位置、左半亮场位置、亮场位置、右半亮场位置、暗场位置;其中为通光面积一半时的剪切位移;
图中标记示意为:1-光源;2-分束光栅;3-样品;4-分析光栅;5-探测器;6-光源光栅。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。应当注意,为了清楚的目的,附图和说明中省略了与本发明无关的、本领域普通技术人员已知的部件和处理的表示和描述。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1(a)和图1(b)为本发明提供的可实现三维动态观测的光栅剪切成像装置的可选结构示意图;在图1(a)中的光栅剪切成像装置中,分束光栅阵列排列在样品台和光源光栅阵列之间的半个第二圆周上,在图1(b)中的光栅剪切成像装置中,分束光栅阵列排列在样品台和分析光栅阵列之间的半个第二圆周上。如图所示,本发明实施例提供的可实现三维动态观测的光栅剪切成像装置包括:用于为样品3提供承载的样品台、包括多个光源1的光源阵列、包括多个光源光栅6的光源光栅阵列、包括多个分束光栅2的分束光栅阵列、包括多个分析光栅4的分析光栅阵列以及包括多个探测器5的探测器阵列。
所述光源阵列和所述光源光栅阵列均排列在以所述样品台为圆心的半个第一圆周上、且所述光源阵列贴近所述光源光栅阵列设置;所述探测器阵列和所述分析光栅阵列均排列在以所述样品台为圆心的另外半个第一圆周上、且所述探测器阵列贴近所述分析光栅阵列设置;所述分束光栅阵列排列在以所述样品台为圆心靠近所述光源光栅阵列的半个第二圆周上,或所述分束光栅阵列排列在以所述样品台为圆心靠近所述分析光栅阵列的半个第二圆周上,所述第一圆周的半径大于所述第二圆周。在每条连接一光源和一探测器的直线上依次设置的光源、光源光栅、分束光栅、分析光栅和探测器,均构成一套成像子系统。
样品可通过样品台装载、调整和固定。光源阵列用于产生多视角同时照明样品的光束阵列。所述光源阵列中每一个光源均能发射一光束,在每一光束上,都有自己的光源光栅、分束光栅、分析光栅和探测器,形成在独立视角获取样品信息的成像子系统;所述光源光栅用于产生多缝光源,其中每一缝光源通过分束光栅均能产生和分析光栅高度对准的一维光束阵列,相邻两缝光源通过分束光栅产生的两套一维光束阵列,以错位一个周期的方式在分析光栅处重合;所述光源数目、光源光栅数目、分束光栅数目、分析光栅数目和探测器数目均相等,且与独立获取样品信息的成像子系统数目相等。
上述光栅剪切成像装置无须转动样品,即可同时对样品各个视角进行成像,图像采集的速度快,基于光栅剪切原理实现相位衬度成像,基于计算机断层成像原理实现三维成像,可以对样品的内部结构实现三维动态观测,可满足医学检测等方面的应用需求。
光栅剪切成像装置中,每套成像子系统相对独立。可选的,在每套成像子系统中:
所述光源为碳纳米管X射线光源,用于在一个视角产生照明样品的X射线光束。
所述光源光栅为一维周期结构的多缝吸收光栅,用于在所述光源产生的光束的照射下产生一维多缝光源;其中每一缝光源通过分束光栅均能产生和分析光栅高度对准的一维光束阵列。
所述分束光栅为一维周期结构的多缝吸收光栅或相位光栅,用于在所述一维多缝光源的照射下,将所述一维多缝光源中来自每一条缝光源的光束分束为一维光束阵列,相邻两缝光源各自产生的一维光束阵列之间错位一个周期。
所述分析光栅为一维周期结构的多缝吸收光栅,用于对所述一维光束阵列进行滤波,产生满足成像要求的光强背景;分析光栅既可以是折射角信号加强或者抑制的调制器,又可以是散射角信号加强或者抑制的调制器。
所述探测器包括多个探测单元构成的一维线阵或二维面阵,用于探测所述光源发出的、经所述光源光栅、分束光栅、样品台承载的样品和所述分析光栅的光束的光强分布。
该方案用于动态观测样品内部的三维结构。
光源光栅、分束光栅和分析光栅是提取样品折射角和散射角信息的关键部件。可选的,所述光源光栅、所述分束光栅和所述分析光栅均为吸收光栅,或所述分束光栅为相位光栅,所述光源光栅和所述分析光栅为吸收光栅。各光栅的厚度可根据实际吸光程度要求设计。一种可选的实施方式例如,所述光源光栅、所述分束光栅或所述分析光栅为吸收光栅时,其栅条厚度为至少使透过光强衰减到入射光强的10%所需的厚度;吸收光栅可采用金或其他重金属制成,以提高其吸光性能;或所述分束光栅为相位光栅时,其栅条厚度需能够使透过光束获得π或π/2的相移。
可选的,所述光源光栅、分束光栅和分析光栅的光栅栅条均竖直设置。
为了提高探测器的角探测灵敏度,所述分束光栅和分析光栅之间的距离为0.1~5米。可选的,第一圆周的半径为1-5米;和/或,所述第二圆周的半径为0.1-4.5米。
可选的,任意两个相邻的光源之间的圆心角均为π/10-π/100。
可选的,所述光源光栅的栅条宽和缝宽的比值为1:1-3:1。
可选的,所述光源光栅的周期与所述分析光栅的周期形成针孔成像关系,针孔是分束光栅上的任意一条缝。
可选的,所述分束光栅的周期为1~100微米。
可选的,所述分束光栅和/或分析光栅的栅条宽和缝宽的比值为1:1的一维多缝光栅,即所述分束光栅的栅条宽和缝宽相等,和/或,所述分析光栅的栅条宽和缝宽相等。所述分析光栅周期由分束光栅周期决定,例如,所述分析光栅的周期等于所述光源中心光线对所述分束光栅周期的几何投影或几何投影的二分之一。
本发明基于光栅剪切成像装置的成像方法中,各组成像子系统中分别执行以下步骤:
步骤1:调整光源,使光源发射出的光束对准沿连接该光源到与其对称设置的探测器的直线上设置的光源光栅、分束光栅、样品台、分析光栅和探测器,并沿第一圆周或第二圆周的径向垂直方向和径向平行方向分别调整光源光栅、分束光栅和分析光栅位置,使分束光栅的一维多缝产生一维光束阵列,使分析光栅的一维多缝对准所述一维光束阵列。
步骤2:测量位移曲线,沿第一圆周或第二圆周的径向垂直方向和径向平行方向分别调整光源光栅、分束光栅和分析光栅位置,以调整所述分析光栅和所述一维光束阵列之间发生的剪切位移,通过探测器的各探测单元均采集光强随剪切位移变化的位移曲线,并用三角函数曲线拟合光强随剪切位移变化的位移曲线;
步骤3:通过探测器探测光强背景和空间位置变化,沿第一圆周或第二圆周的切向调整光源光栅、分束光栅或分析光栅位置,直至把分析光栅和相应的所述一维光束阵列之间的剪切位移调整在光强背景满足成像要求的采集位置。
步骤4:通过所述探测器采集固定在样品台上的样品对应所述光强背景下的投影像,并将探测器阵列同时采集的多个角度的投影像构成一套投影数据。
基于上述光栅剪切成像装置在成像过程中,无须转动样品,即可同时对样品各个视角进行成像,图像采集的速度快,基于光栅剪切原理实现相位衬度成像,基于计算机断层成像原理实现三维成像,可以对样品的内部结构实现三维动态观测,可满足医学检测等方面的应用需求。
为了实现对样品内部局部区域的实时动态观测,可选的,所述光栅剪切三维成像方法还可包括:
步骤5:根据物函数和位移曲线的卷积运算,建立光栅剪切成像方程,根据所述光栅剪切成像方程求出探测器采集的投影数据的数学表达式,根据得到的数学表达式确定至少一种预定的物理参量投影数据的半定量表达式和/或定量表达式。
步骤6:利用CT三维成像原理,根据所述至少一种物理参量投影数据的半定量表达式重建样品的相应物理参量的三维分布,并选择一种物理参量重建后的三维分布对样品内部局部区域进行三维动态观测;或者,利用CT三维成像原理,根据所述至少一种物理参量投影数据的定量表达式重建样品的相应物理参量的三维分布,并选择一种物理参量重建后的三维分布对样品内部局部区域进行三维动态观测。
可选的,所述光强背景包括:亮场背景、左半亮场背景、右半亮场背景和/或暗场背景;所述投影像包括:所述样品在所述亮场背景下的亮场投影数据、在所述左半亮场背景下的左半亮场投影数据、在所述右半亮场背景下的右半亮场投影数据、和/或在所述暗场背景下的暗场投影数据;所述至少一种物理参量投影数据包括:线性吸收系数投影数据、线性散射系数投影数据或线性消光系数投影数据、折射率实部衰减率导数投影数据;所述的根据上述投影数据重建的物理参量包括:线性吸收系数、线性散射系数或线性消光系数、折射率实部衰减率和折射率实部衰减率导数。
换而言之,上述方案基于光栅剪切成像装置成像方法可分为两个层次,第一层次是三维动态半定量成像方法,只须一次曝光,就可以获得三维动态观测的半定量投影数据;第二层次是三维快速精确成像方法,只须三次曝光,就可以获得三维快速精确重建的定量投影数据。
例如,按照光束传播顺序,每个成像子系统的光栅剪切成像依次由三个简单的物理过程组成:(1)分束光栅对入射光束进行分束,产生照射样品的一维光束阵列,(2)样品对一维光束阵列产生吸收、折射和散射作用,(3)分析光栅选择合适位置,产生满足成像要求的光强背景,并对经过样品作用的一维光束阵列进行滤波,探测器采集样品吸收、折射和散射产生的光强信号。成像过程一般可以写成物函数和成像系统脉冲响应函数的卷积。因此,上述三个简单的物理过程可以重组为无样品过程和有样品过程。无样品过程反映的是成像系统的特性,通过它可以获得描述成像系统特性的脉冲响应函数;有样品过程描述样品经过成像系统作用的成像过程,可以用物函数和成像系统脉冲响应函数的卷积来表达。要对光栅剪切成像进行数学描述,首先要建立物函数的数学模型,其次求出描述成像系统特性的脉冲响应函数,然后计算物函数和脉冲响应函数的卷积,完成对整个成像过程的数学描述。下面的推导步骤为:第一步,建立样品对X射线作用的数学模型,写出物函数的数学表达,第二步,推导分析光栅对分束光栅一维光束阵列滤波的数学表达,求出描述成像系统特性的脉冲响应函数(即位移曲线),第三步,计算物函数和脉冲响应函数(即位移曲线)的卷积,推导出光栅剪切成像方程。
第一步:在建立数学模型之前,先对样品中一点进行定义。在二维成像中,样品所在物面上一点(x,y),不是一个二维几何点,而是一个以(x,y)为中心的物面积元ΔxΔy,Δx和Δy的大小由探测器分辨率决定;在三维成像中,样品体内一点(x,y,z),不是一个三维几何点,而是一个以(x,y,z)为中心的物体积元ΔxΔyΔz,Δx、Δy和Δz的大小由探测器分辨率决定。在此特别说明,下面提到的物点,其含义均为物面积元或物体积元。
样品对X射线产生吸收、折射和散射三种作用。
吸收(包括非弹性散射)是一个X射线能量在样品中转化为热能的耗散过程,如图2所示,描述了样品对X射线光束吸收衰减作用图像。根据图2,样品中一点(x,y)对通过该点光线的吸收可以表达为:
上式左边表示入射光束,右边表示出射光束,其中表示光束角度矢量,
其中μ(x,y,z)为样品的线性吸收系数。(1)式的物理意义为,吸收引起光线强度衰减,但不改变光线方向。(1)式还可以表示为:
折射是一个能量守恒的过程,如图3所示,描述了样品对X射线光束折射作用的图像。根据图3,样品中一点(x,y)对通过该点光线的折射可以表达为:
上式左边表示入射光束,右边表示出射光束,其中表示光束的角度矢量,为折射角矢量,其数学表达为:
其中δ(x,y,z)为样品折射率实部衰减率。(4)式的物理意义为,折射改变光线方向,但不改变光线强度。(4)式还可以写为分量表达式:
散射由面积元内部很多小颗粒的多重折射引起,也是一个能量守恒的过程,如图4所示,描述了样品对入射X射线散射作用的图像。散射和折射的不同之处在于,折射把样品物面上一个面积元作为一个整体来研究,即把样品物面上一个面积元作为一个微小棱镜,而散射则研究这个面积元内部的不均匀性质,即相当于研究微小棱镜内部的气泡、颗粒、微孔、微晶和杂质等不均匀。因此,对于每个面积元,只有一条折射光线和一个折射角,却有多条散射光线和多个散射角。换言之,散射是一个光束分散的过程。因为样品有一定厚度,在面积元内部沿着光束传播方向,各小颗粒分布是随机的,前后两个小颗粒产生的折射是相互独立的,小颗粒每次折射使入射光线偏离入射方向的角度是随机的,所以根据中心极限定理,散射角是以入射角(或者折射角)为中心服从二维正态统计分布,可以用方差来描述散射角分布范围。根据图4,一光线射入样品时,由于散射引起分散,出射光线分为两部分,散射光线和未散射光线,未散射光线仍然沿着入射方向传播,而散射光线偏离入射方向传播。随着光线在样品中穿行,散射事件的不断发生,散射光线不断产生,散射能量逐渐增强,而未散射光线逐渐消弱,未散射能量逐渐减弱,称为消光衰减。需要特别说明的是,每条光线都可能遇到面积元内部多个小颗粒的折射,需要把第一次被小颗粒折射产生第一次偏离的散射光线和该散射光线继续被后续小颗粒折射产生进一步偏离的散射光线分别考虑,这是因为散射能量是由一次折射产生的一次散射决定的,而以后的多次折射产生的多次散射仅仅使散射能量分布范围更大,增加散射角方差,而对增加或减少散射能量几乎不起作用。简言之,小颗粒一次折射决定散射能量和未散射能量的比例,小颗粒多次折射决定散射角方差。因此,消光衰减和吸收衰减一样遵循比尔定律。设入射光线能量为1,未散射光线继续沿着入射方向传播,其所携带能量,即消光衰减为exp(-Γ(x,y)),散射光线能量为1-exp(-Γ(x,y))。在散射中心对称的条件下,样品一点(x,y)对通过该点光线的散射可以表达为:
上式左边表示入射光束,右边表示出射光束,其中表示光束角度矢量。(7)式还可以写为分量表达式:
在(7)式或(8)式右边第一项中
其中γ(x,y,z)为样品的线性消光系数,第二项中σ2(x,y)为(x,y)点处样品整体厚度产生的散射角方差。因为样品整体厚度的散射角方差σ2(x,y),是光线传播路径上一系列厚度为Δzi薄片的微分散射角方差Δσ2(x,y,z)之和,所以样品整体厚度的散射角方差可以表示为
其中ω(x,y,z)为线性散射系数。为了建立线性散射系数和线性消光系数之间的关系,把(10)式与(9)式进行比较,可得:
ω(x,y,z)=ε(x,y,z)γ(x,y,z), (11)
其中ε(x,y,z)为扩散因子。若样品是由散射性质相同的材料构成时,扩散因子ε(x,y,z)就是常量,则下式成立:
σ(x,y)=ε·Γ(x,y)。 (12)
此时就可以从一种信号获得另一种信号。换言之,若样品是由散射性质相同的材料构成时,则两个几何意义不同的散射信号可以归结为一个信号。
综合考虑上述三种作用,在样品散射中心对称的条件下,样品中一点(x,y)对通过该点光束的作用可以用物函数表达,
其分量表达式为:
根据(14)式,可知仅在x方向起作用的物函数为:
仅在y方向起作用的物函数为:
根据(13)式或(14)式,可知出射X射线携带了以下四种样品信号:
(一)吸收衰减exp(-Μ(x,y)),其中Μ(x,y)为线性吸收系数μ(x,y,z)的投影路径积分
(二)折射角 ,其中为X方向的单位矢量,为Y方向的单位矢量,θx(x,y)为样品折射率实部衰减率δ(x,y,z)X方向的偏导数的投影路径积分θy(x,y)为样品折射率实部衰减率δ(x,y,z)Y方向的偏导数的投影路径积分
(三)消光衰减exp(-Γ(x,y)),其中Γ(x,y)为线性消光系数γ(x,y,z)的投影路径积分
(四)散射角方差σ2(x,y),是线性散射系数的投影路径积分:
线性散射系数和线性消光系数之间的关系为:
ω(x,y,z)=ε(x,y,z)γ(x,y,z),
其中ε(x,y,z)为扩散因子。若样品是由散射性质相同的材料构成,扩散因子ε就为常数,则线性散射系数和线性消光系数之间的关系为:
ω(x,y,z)=ε·γ(x,y,z),
消光衰减和散射角方差之间的关系为:
σ2(x,y)=ε·Γ(x,y)。
第二步:推导分析光栅对分束光栅一维光束阵列滤波作用的数学表达。因为分析光栅和分束光栅一维光束阵列都是一维周期函数,具有几乎相同的函数分布,分析光栅相对分束光栅一维光束阵列的剪切位移可以调整两者之间的相关性,所以分析光栅对分束光栅一维光束阵列的滤波作用在数学上是互相关运算。当沿x轴方向移动光源光栅或者分束光栅或者分析光栅时,都会引起一维光束阵列和分析光栅之间发生剪切位移,探测器每个像素都可以测到光强随剪切位移而变化的位移曲线,图5为一维光束阵列和分析光栅之间沿x轴方向进行剪切位移时的位移曲线;因为位移曲线形状类似余弦曲线,为了简化提取折射和散射信息的数学表达,所以用余弦曲线拟合位移曲线,其表达式为:
其中为分析光栅相对分束光栅产生的一维光束阵列沿X轴方向的剪切角位移,xg为分析光栅相对分束光栅产生的一维光束阵列沿X轴方向的剪切位移,D为分束光栅和分析光栅之间在光束传播方向上的距离,p为位移曲线周期,也是分析光栅的周期,为无样品时位移曲线平均值,Smax和Smin分别为位移曲线的最大值和最小值,为无样品时位移曲线的可见度。图5中位移曲线上a、b、c、d四点分别对应亮场、右半亮场、左半亮场、暗场。亮场代表分束光栅产生的一维光束阵列几乎全部通过分析光栅,暗场代表分束光栅产生的一维光束阵列几乎全被分析光栅阻挡,半亮场代表分束光栅产生的一维光束阵列中,一半被分析光栅阻挡,一半通过分析光栅。
位移曲线(17)式是描述光栅剪切成像系统的特性函数,它就是光栅剪切成像系统的脉冲响应函数。虽然位移曲线是一条无界的周期性曲线,不满足脉冲响应函数是局域函数的一般要求,但是实际上位移曲线的作用范围满足局域函数的要求。为了说明这一点,来考察物函数的作用范围。样品的折射角和散射角为微弧度量级,样品到探测器的距离为米量极,因而样品折射和散射引起的光线位移为微米量极。由此可知,虽然光栅剪切成像系统不存在一个有界的特性函数,但是样品每一点的物函数本身是局域函数,它决定了位移曲线在成像过程中,只能在物函数限制的局域范围内发挥作用,因而光栅剪切成像系统的脉冲响应函数是位移曲线上一有限范围的分段函数。因为位移曲线上下起伏,不同位置的分段函数具有不同的脉冲响应,所以本发明为光栅剪切成像系统提出四种脉冲响应函数:
(1)位于位移曲线峰位的亮场脉冲响应函数,这是一段偶函数,参见图5中位移曲线上a点附近函数曲线,其主要作用是抑制样品的散射信号;
(2)位于位移曲线谷位的暗场脉冲响应函数,这是一段偶函数,参见图5中位移曲线上d点附近函数曲线,其主要作用是增强样品的散射信号;
(3)位于位移曲线腰位的右半亮场脉冲响应函数,这是一段奇函数,参见图5中位移曲线上c点附近函数曲线,其主要作用是提取样品向右的折射信号;
(4)位于位移曲线腰位的左半亮场脉冲响应函数,这是一段奇函数,参见图5中位移曲线上b点附近函数曲线,其主要作用是提取样品向左的折射信号。
上述四种脉冲响应函数,其对折射和散射信号的探测灵敏度取决于分束光栅和分析光栅之间的距离和它们的周期,距离越长、周期越小,探测灵敏度越高。
第三步:推导光栅剪切成像方程。前面已经讨论清楚,当样品放在样品台上后,探测器在分析光栅后面测得的光强分布是物函数和位移曲线的卷积。光栅剪切成像方程可从x轴方向作用的物函数Ox(x,y,θg)和位移曲线S(θg)的卷积推导而出:
其中,I0为无样品时分束板前的光强,V(x,y)为放入样品后位移曲线的可见度,又称为样品的可见度像,且有
为无样品时位移曲线的可见度,exp(-Γ(x,y))为样品的消光衰减像,σ2(x,y)为样品的散射角方差像。
本发明根据光栅剪切成像方程,提出了一套采集投影数据的方法。先考虑一个成像子系统的成像,然后考虑整个成像装置多个成像子系统同时成像多视角投影数据。
在一个成像子系统中,位移曲线的峰位和谷位、右腰位和左腰位具有特殊的物理意义。在位移曲线峰位,分束光栅产生的一维光束阵列几乎全部通过分析光栅,成像背景是亮场;在位移曲线谷位,分束光栅产生的一维光束阵列几乎全部被分析光栅阻挡,成像背景是暗场;在位移曲线右腰位,分束光栅产生的一维光束阵列中,每一微光束的左半边被分析光栅阻挡,而其右半边通过分析光栅照亮探测器,形成右半亮场成像背景,在位移曲线左腰位,分束光栅产生的一维光束阵列中,每一微光束的右半边被分析光栅阻挡,而其左半边通过分析光栅照亮探测器,形成左半亮场成像背景;当把样品分别放入亮场、暗场、右半亮场和左半亮场,探测器可以分别拍摄到亮场投影像、暗场投影像、右半亮场投影像和左半亮场投影像;亮场投影像的特点是:经过样品方向不变的信号强,而方向改变的信号弱;暗场投影像的特点是:经过样品方向不变的信号弱,而方向改变的信号、特别是杂乱无章的散射信号强;右半亮场投影像的特点是:向右偏转的光强信号得到加强,而向左偏转的光强信号得到抑制;左半亮场投影像的特点是:向左偏转的光强信号得到加强,而向右偏转的光强信号得到抑制。
上述成像子系统的成像机制和成像过程适用于成像装置内的其他成像子系统,只须为描述一个成像子系统的二维成像数学表达式中增加一个角度坐标,就可描述多个成像子系统同时采集多视角投影数据的成像过程。下面推导出包括亮场投影数据、暗场投影数据、右半亮场投影数据和左半亮场投影数据等投影像的数学表达式。
通过所述探测器采集固定在样品台上的样品在亮场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在亮场背景位置,即剪切位移xg=0,剪切角位移各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到亮场投影数据其数学表达式为:
通过所述探测器采集固定在样品台上的样品在暗场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在暗场背景位置,即剪切位移剪切角位移各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到暗场投影数据其数学表达式为:
通过所述探测器采集固定在样品台上的样品在右半亮场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在右半亮场背景位置,即剪切位移剪切角位移各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到右半亮场投影数据其数学表达式为:
通过所述探测器采集固定在样品台上的样品在左半亮场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在左半亮场背景位置,即剪切位移剪切角位移各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到左半亮场投影数据其数学表达式为:
本发明根据样品对X射线的吸收、折射和散射三种相互作用,在光栅剪切成像中定义了三种物理参量投影数据,分别为线性吸收系数投影数据
折射率实部衰减率导数投影数据
线性散射系数投影数据
或线性消光系数投影数据
下面描述获得三种物理参量投影数据半定量表示式、定量表达式的获取方法。
线性吸收系数投影数据半定量表示式的获取方法:在忽略样品折射和散射的条件下,有
根据亮场投影数据(20)式或暗场投影数据(21)式,可得
或者,
折射率实部衰减率导数投影数据半定量表示式的获取方法:在忽略样品吸收和散射的条件下,有
根据右半亮场投影数据(22)式或左半亮场投影数据(23)式,可得
或者,
在忽略样品吸收和折射的条件下,有:
根据亮场投影数据(20)式,获得可见度的半定量表达式为:
或者,根据暗场投影数据(21)式,获得可见度的半定量表达式为:
在弱散射条件下,
有
将(31)式代入(19)式,可得可见度和散射角方差的关系为:
线性散射系数投影数据半定量表示式的获取方法:把可见度像(28)式代入(32)式,可得:
或者,把可见度像(29)式代入(32)式,可得:
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性消光系数投影数据的半定量表达式为:
或者,
在强散射条件下,
有:
有:
把(36)式代入(19)式,得:
V(x,y)≈V0exp(-Γ(x,y)), (40)
线性散射系数投影数据半定量表示式的获取方法:把可见度像(28)式代入(40)式,可得:
或者,把可见度像(29)式代入(40)式,可得:
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性散射系数投影数据的半定量表达式为:
或者,
线性吸收系数投影数据定量表示式的获取方法:根据亮场投影数据(20)式和暗场投影数据(21)式,可得:
根据右半亮场投影数据(22)式和左半亮场投影数据(23)式,得
折射率实部衰减率导数投影数据定量表示式的获取方法:根据右半亮场投影数据(22)式、左半亮场投影数据(23)式、亮场投影数据(20)式和暗场投影数据(21)式,可从下列方程组获得折射率实部衰减率导数投影数据定量表示式:
在弱散射条件下,线性散射系数投影数据定量表达式的获取方法:根据右半亮场投影数据(22)式、左半亮场投影数据(23)式、亮场投影数据(20)式和暗场投影数据(21)式,可从下列方程组获得线性散射系数投影数据定量表示式:
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性消光系数投影数据的定量表达式为:
在强散射条件下,线性消光系数投影数据定量表达式的获取方法:根据右半亮场投影数据(22)式、左半亮场投影数据(23)式、亮场投影数据(20)式和暗场投影数据(21)式,可从下列方程组获得线性消光系数投影数据定量表示式:
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,确定所述线性散射系数投影数据的定量表达式的方程组为:
通过光源阵列、光源光栅阵列、分束光栅阵列、分析光栅阵列和探测器阵列,针对样品内的局部区域采集到的亮场投影数据、半亮场投影数据和暗场投影数据,以及由这些投影数据表达的半定量投影数据和定量投影数据,对样品整体而言为局部照明(或称为照明截断)的不完备投影数据,对样品内局部区域而言为多个分离角度照明(或称为角度截断)的不完备投影数据,利用传统CT重建理论不可能精确重建样品内的局部区域。然而,根据内部CT局部重建理论和压缩传感理论,特别像人体内部器官,在待重建的器官内存在部分已知区域的条件下,可以精确重建人体内部器官。因此,利用一次曝光获得的半定量投影数据可以对人体内部器官实现三维动态观测和成像;利用三次曝光获得的定量投影数据可以对人体内部器官实现快速精确三维重建。
根据亮场投影数据、暗场投影数据、右半亮场投影数据和左半亮场投影数据各自的数学表达式,获得的半定量投影数据,再根据CT三维成像原理,获得重建线性吸收系数、线性散射系数或线性消光系数、折射率实部衰减率和折射率实部衰减率导数的半定量表达式的方法为:
设(x,y,z)为X射线光束坐标系,(x′,y′,z′)为样品坐标系,两套坐标系之间的关系为
在忽略样品折射和散射的条件下,线性吸收系数的重建算法公式为:
或者,
其中F-1代表一维逆傅立叶变换,或者为得耳塔函数,ρ为与空间坐标x对应的空间频谱坐标;
在忽略样品吸收和折射、并在弱散射条件下,线性散射系数的重建算法公式为:
或者,
在忽略样品吸收和折射、并在强散射条件下,线性消光系数的重建算法公式为:
或者,
在忽略样品吸收和散射的条件下,折射率实部衰减率的重建算法公式为:
或者,
其中i为虚数;
在忽略样品吸收和散射的条件下,折射率实部衰减率导数的重建算法公式为:
或者,
其中-N/2<k≤N/2,改变k,可以重建不同方向的折射率实部衰减率导数。
根据亮场投影数据、暗场投影数据、右半亮场投影数据和左半亮场投影数据之间的定量关系,获得的定量投影数据,再根据CT三维成像原理,获得重建线性吸收系数、线性散射系数或线性消光系数、折射率实部衰减率和折射率实部衰减率导数的定量表达式的方法为:
设(x,y,z)为X射线光束坐标系,(x′,y′,z′)为样品坐标系,两套坐标系之间的关系为
线性吸收系数的重建算法公式为:
其中线性吸收系数投影数据可选以下两个定量表达式之一:
或者,
在弱散射条件下,线性散射系数的重建算法公式为:
其中所述线性散射系数投影数据的定量表达式可从下列方程组获得:
在强散射条件下,线性消光系数的重建算法公式为:
其中所述线性消光系数投影数据的定量表达式可从下列方程组获得:
折射率实部衰减率的重建算法公式为:
折射率实部衰减率导数的重建算法公式为:
其中折射率实部衰减率导数投影数据的定量表达式可从下列方程组获得:
其中,i为虚数,-N/2<k≤N/2,改变k,可以重建不同方向的折射率实部衰减率导数。
通过上述分析可知,本发明提供的剪切三维成像装置和方法无须转动样品,即可同时对样品各个视角进行成像,图像采集的速度快,基于光栅剪切原理实现相位衬度成像,基于计算机断层成像原理实现三维成像,密度分辨率高,密度不均匀性分辨率高,例如,一次曝光就可采集到进行三维动态观测的半定量投影数据;三次曝光就可以采集到进行快速精确重建的定量投影数据,可实现对样品的内部结构实现三维动态观测,可满足医学检测等方面的应用需求。
此外,本发明还提供了一种医学检测设备,该医学检测设备包括有如上述任一种光栅剪切成像装置。有关光栅剪切成像装置和方法的技术方案,可参见上文描述,在此不再赘述。
在本发明上述各实施例中,实施例的序号仅仅便于描述,不代表实施例的优劣。对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
在本发明的装置和方法等实施例中,显然,各部件或各步骤是可以分解、组合和/或分解后重新组合的。这些分解和/或重新组合应视为本发明的等效方案。同时,在上面对本发明具体实施例的描述中,针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用,与其它实施方式中的特征相组合,或替代其它实施方式中的特征。
应该强调,术语“包括/包含”在本文使用时指特征、要素、步骤或组件的存在,但并不排除一个或更多个其它特征、要素、步骤或组件的存在或附加。
最后应说明的是:虽然以上已经详细说明了本发明及其优点,但是应当理解在不超出由所附的权利要求所限定的本发明的精神和范围的情况下可以进行各种改变、替代和变换。而且,本发明的范围不仅限于说明书所描述的过程、设备、手段、方法和步骤的具体实施例。本领域内的普通技术人员从本发明的公开内容将容易理解,根据本发明可以使用执行与在此所述的相应实施例基本相同的功能或者获得与其基本相同的结果的、现有和将来要被开发的过程、设备、手段、方法或者步骤。因此,所附的权利要求旨在在它们的范围内包括这样的过程、设备、手段、方法或者步骤。
Claims (15)
1.一种可实现三维动态观测的光栅剪切成像装置,其特征在于,包括:
样品台、包括多个光源的光源阵列、包括多个光源光栅的光源光栅阵列、包括多个分束光栅的分束光栅阵列、包括多个分析光栅的分析光栅阵列以及包括多个探测器的探测器阵列;
所述光源阵列和所述光源光栅阵列均排列在以所述样品台为圆心的半个第一圆周上、且所述光源阵列贴近所述光源光栅阵列设置;所述探测器阵列和所述分析光栅阵列均排列在以所述样品台为圆心的另外半个第一圆周上、且所述探测器阵列贴近所述分析光栅阵列设置;所述分束光栅阵列排列在以所述样品台为圆心靠近所述光源光栅阵列或靠近所述分析光栅阵列的半个第二圆周上,所述第一圆周的半径大于所述第二圆周;
在每条连接一光源和一探测器的直线上依次设置的光源、光源光栅、分束光栅、分析光栅和探测器,均构成一套成像子系统。
2.根据权利要求1所述的光栅剪切成像装置,其特征在于,在每套成像子系统中:
所述光源为碳纳米管X射线光源;
所述光源光栅为一维周期结构的多缝吸收光栅,用于在所述光源产生的光束的照射下产生一维多缝光源;
所述分束光栅为一维周期结构的多缝吸收光栅或相位光栅,用于在所述一维多缝光源的照射下,将所述一维多缝光源中来自每一条缝光源的光束分束为一维光束阵列,相邻两缝光源各自产生的一维光束阵列之间错位一个周期;
所述分析光栅为一维周期结构的多缝吸收光栅,用于产生成像所要求的光强背景,对所述一维光束阵列进行滤波;
所述探测器包括多个探测单元构成的一维线阵或二维面阵,用于探测所述光源发出的、经所述光源光栅、分束光栅、样品台承载的样品和所述分析光栅的光束的光强分布。
3.根据权利要求1或2所述的光栅剪切成像装置,其特征在于,
所述光源光栅、分束光栅和分析光栅的光栅栅条均竖直设置;和/或,
所述第一圆周的半径为1-5米;和/或,
所述第二圆周的半径为0.1-4.5米;和/或,
任意两个相邻的光源之间的圆心角均为π/10-π/100;和/或,
所述光源光栅的栅条宽和缝宽的比值为1:1-3:1;和/或,
所述光源光栅的周期与所述分析光栅的周期形成针孔成像关系,针孔是分束光栅上的任意一条缝;和/或,
所述分束光栅的周期为1~100微米;和/或,
所述分束光栅的栅条宽和缝宽相等;和/或,
所述分析光栅的周期等于所述光源中心光线对所述分束光栅周期的几何投影或几何投影的二分之一;和/或,
所述分析光栅的栅条宽和缝宽相等。
4.根据权利要求3所述的光栅剪切成像装置,其特征在于,所述光源光栅、所述分束光栅或所述分析光栅为吸收光栅时,其栅条厚度为至少使透过光强衰减到入射光强的10%所需的厚度;或所述分束光栅为相位光栅时,其栅条厚度需能够使透过光束获得π或π/2的相移。
5.一种基于如权利要求1-4任一所述的可实现三维动态观测的光栅剪切成像装置的光栅剪切成像方法,其特征在于,在各组成像子系统中分别执行以下步骤:
调整光源,使光源发射出的光束对准沿连接该光源到与其对称设置的探测器的直线上设置的光源光栅、分束光栅、样品台、分析光栅和探测器,并沿第一圆周或第二圆周的径向垂直方向和径向平行方向分别调整光源光栅、分束光栅和分析光栅位置,使分束光栅的一维多缝产生一维光束阵列,使分析光栅的一维多缝对准所述一维光束阵列;
测量位移曲线,沿第一圆周或第二圆周的径向垂直方向和径向平行方向分别调整光源光栅、分束光栅和分析光栅位置,以调整所述分析光栅和所述一维光束阵列之间发生的剪切位移,通过探测器的各探测单元均采集光强随剪切位移变化的位移曲线,并用三角函数曲线拟合光强随剪切位移变化的位移曲线;
通过探测器探测光强背景和空间位置变化,沿第一圆周或第二圆周的切向调整光源光栅、分束光栅或分析光栅位置,直至把分析光栅和相应的 所述一维光束阵列之间的剪切位移调整在光强背景满足成像要求的采集位置;
通过所述探测器采集固定在样品台上的样品对应所述光强背景下的投影像,并将探测器阵列同时采集的多个角度的投影像构成一套投影数据。
6.根据权利要求5所述的光栅剪切成像方法,其特征在于,所述通过所述探测器采集固定在样品台上的样品对应所述光强背景下的投影数据之后,还包括:
根据物函数和位移曲线的卷积运算,建立光栅剪切成像方程,根据所述光栅剪切成像方程求出探测器采集的投影数据的数学表达式,根据得到的数学表达式确定至少一种预定的物理参量投影数据的半定量表达式和/或定量表达式;
利用CT三维成像原理,根据所述至少一种预定的物理参量投影数据的半定量表达式重建样品的相应物理参量的三维分布,并选择一种重建的物理参量对样品内部局部区域进行三维动态观测;或者,利用CT三维成像原理,根据所述至少一种预定的物理参量投影数据的定量表达式重建样品的相应物理参量的三维分布,并选择一种重建的物理参量对样品内部局部区域进行三维动态观测。
7.根据权利要求6所述的光栅剪切成像方法,其特征在于,
所述光强背景包括:亮场背景、左半亮场背景、右半亮场背景和/或暗场背景;
所述投影像包括:所述样品在所述亮场背景下的亮场投影数据、在所述左半亮场背景下的左半亮场投影数据、在所述右半亮场背景下的右半亮场投影数据、和/或在所述暗场背景下的暗场投影数据;
所述至少一种预定的物理参量投影数据包括:线性吸收系数投影数据、线性散射系数投影数据或线性消光系数投影数据、折射率实部衰减率导数投影数据;所述的根据上述投影数据重建的物理参量包括:线性吸收系数、线性散射系数或线性消光系数、折射率实部衰减率和折射率实部衰减率导数。
8.根据权利要求7所述的光栅剪切成像方法,其特征在于,
所述分束光栅阵列产生的一维光束阵列和分析光栅在垂直光束传播方 向的平面内,沿垂直于栅条的方向进行剪切位移时,所述位移曲线用余弦函数近似表达为:
其中为分析光栅相对于一维光束阵列的剪切角位移,xg分析光栅相对于一维光束阵列的剪切位移,D为分束光栅和分析光栅之间的距离,p为位移曲线周期,也是分析光栅的周期,为无样品时位移曲线平均值,为无样品时位移曲线的可见度,Smax和Smin分别为位移曲线的最大值和最小值。
9.根据权利要求8所述的光栅剪切成像方法,其特征在于,在所述样品散射中心对称的条件下,用物函数表达样品中一点(x,y)对通过该点光线的作用,所述物函数表达式为:
其中,指的是:样品中一点对通过该点光线的折射;
或
其中仅在x方向作用的物函数为:
仅在y方向作用的物函数为:
其中,表示光束的角度矢量,和分别为平行和垂直于圆周平面的分量;
在物函数中,吸收衰减像的数学表达为:
exp(-Μ(x,y)),
其中Μ(x,y)为线性吸收系数μ(x,y,z)的投影路径积分:
折射角像的数学表达为:
其中为x方向的单位矢量,为y方向的单位矢量,为样品垂直于光束传播方向的折射角矢量,θx(x,y)为其沿X方向的分量,平行于圆周平面,也是折射率实部衰减率δ(x,y,z)沿X方向偏导数的投影路径积分:
θy(x,y)为其沿Y方向的分量,垂直于圆周平面,
也是样品折射率实部衰减率δ(x,y,z)沿Y方向偏导数的投影路径积分:消光衰减像的数学表达为:
exp(-Γ(x,y)),
其中Γ(x,y)为线性消光系数γ(x,y,z)的投影路径积分:
散射角方差像的数学表达为:
σ2(x,y),
它是各微分薄层散射角方差dσ2(x,y,z)的投影路径积分:
其中ω(x,y,z)为线性散射系数,其与线性消光系数之间的关系为:
ω(x,y,z)=ε(x,y,z)γ(x,y,z),
其中ε(x,y,z)为扩散因子,若样品由散射性质相同的材料构成,
ε(x,y,z)为常量,则下式成立:
σ2(x,y)=ε·Γ(x,y)。
10.根据权利要求9所述的光栅剪切成像方法,其特征在于:当一维光束阵列和分析光栅沿垂直于栅条的方向进行剪切位移时,光栅剪切成像方程由仅在x轴方向作用的物函数和位移曲线的卷积推导得出:
其中Ox(x,y,θx)表示仅在x轴方向作用的物函数,为分析光栅相对于一维光束阵列的剪切角位移,xg为分析光栅相对于一维光束阵列的剪切位移,D为分束光栅和分析光栅之间的距离,p为位移曲线周期,也是分析光栅的周期,I(x,y,θg)为样品内一点(x,y)在剪切角位移为θg时的光强,I0为入射光光强,为无样品时位移曲线平均值,Smax和Smin分别为位移曲线的最大值和最小值,θx(x,y)表示折射角在x轴的分量,也是样品折射率实 部衰减率δ(x,y,z)在x方向的偏导数的投影路径积分,表达式为:
V(x,y)为有样品时位移曲线的可见度,又称样品的可见度像,其表达式为:
为无样品时位移曲线的可见度,exp(-Γ(x,y))为消光衰减像,σ2(x,y)为散射角方差像。
11.根据权利要求10所述的光栅剪切成像方法,其特征在于:
探测器采集固定在样品台上的样品在亮场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在亮场背景位置,即剪切位移xg=0,剪切角位移 各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到亮场投影数据其数学表达式为:
探测器采集固定在样品台上的样品在暗场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光 栅之间的剪切位移固定在暗场背景位置,即剪切位移剪切角位移 各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到暗场投影数据其数学表达式为:
探测器采集固定在样品台上的样品在右半亮场光强背景下的投影数据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在右半亮场背景位置,即剪切位移剪切角位移各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到右半亮场投影数据其数学表达式为:
探测器采集固定在样品台上的样品在左半亮场光强背景下的投影数 据,包括:在每个成像子系统中,把分束光栅产生的一维光束阵列和相应分析光栅之间的剪切位移固定在左半亮场背景位置,即剪切位移剪切角位移各成像子系统的光源,即按圆心角排列的光源阵列,同时发光,各成像子系统的探测器,即与光源阵列对称排列的探测器阵列,同时拍摄样品的投影像,获得探测器阵列在光源阵列一次曝光瞬间采集到左半亮场投影数据其数学表达式为:
12.根据权利要求11所述的光栅剪切成像方法,其特征在于:
确定所述线性吸收系数投影数据的半定量表达式的方法为:在忽略样品折射和散射的条件下,即有:
线性吸收系数投影数据的半定量表达式为:
或者,
确定所述折射率实部衰减率导数投影数据的半定量表达式的方法为:在忽略样品吸收和散射的条件下,即有:
折射率实部衰减率导数投影数据的半定量表达式为:
或者,
在忽略样品吸收和折射的条件下,即有:
可见度的半定量表达式为:
或者
在弱散射条件下,可见度和散射角方差的关系为:
确定线性散射系数投影数据的半定量表达式的方法为:
或者,
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性消光系数投影数据的半定量表达式为:
或者,
在强散射条件下,可见度和散射角方差的关系为:
确定线性消光系数投影数据的半定量表达式为:
或者,
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性散射系数投影数据的半定量表达式为:
或者,
13.根据权利要求11所述的光栅剪切成像方法,其特征在于:
确定所述线性吸收系数投影数据的定量表达式为:
或者,
确定所述折射率实部衰减率导数投影数据的定量表达式的方程组为:
确定所述线性散射系数投影数据和所述线性消光系数投影数据的定量表达式方法为:
在弱散射条件下,可从以下方程组获得所述线性散射系数投影数据的定量表达式:
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,确定所述线性消光系数投影数据的定量表达式的方程组为:
在强散射条件下,可从以下方程组获得所述线性消光系数投影数据的 定量表达式:
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,确定所述线性散射系数投影数据的定量表达式的方程组为:
14.根据权利要求11所述的光栅剪切成像方法,其特征在于,所述利用CT三维成像原理,根据所述至少一种物理参量投影数据的半定量表达式重建样品的相应物理参量的三维分布,包括:
设(x,y,z)为X射线光束坐标系,(x′,y′,z′)为样品坐标系,两套坐标系之间的关系为:
在忽略样品折射和散射的条件下,线性吸收系数的重建算法公式为:
或者,
其中F-1代表一维逆傅立叶变换,或者为得耳塔函数,ρ为与空间坐标x对应的空间频谱坐标;
在忽略样品吸收和折射、且在弱散射条件下,线性散射系数的重建算法公式为:
或者,
在忽略样品吸收和折射、并在强散射条件下,线性消光系数的重建算法公式为:
或者,
在忽略样品吸收和散射的条件下,折射率实部衰减率的重建算法公式为:
或者,
其中i为虚数;
在忽略样品吸收和散射的条件下,折射率实部衰减率导数的重建算法公式为:
或者,
其中-N/2<k≤N/2,改变k,可以重建不同方向的折射率实部衰减率导数。
15.根据权利要求11所述的光栅剪切成像方法,其特征在于,
所述利用CT三维成像原理,根据所述至少一种物理参量投影数据的定量表达式重建样品的相应物理参量的三维分布包括:
设(x,y,z)为X射线光束坐标系,(x′,y′,z′)为样品坐标系,两套坐标系之间的关系为:
线性吸收系数的重建算法公式为:
其中所述线性吸收系数投影数据可选以下两个定量表达式之一:
或者,
在弱散射条件下,线性散射系数的重建算法公式为:
其中所述线性散射系数投影数据的定量表达式可从下列方程组获得:
在强散射条件下,线性消光系数的重建算法公式为:
其中所述线性消光系数投影数据的定量表达式可从下列方程组获 得:
折射率实部衰减率的重建算法公式为:
折射率实部衰减率导数的重建算法公式为:
其中折射率实部衰减率导数投影数据的定量表达式可从下列方程组获得:
其中,i为虚数,-N/2<k≤N/2,改变k,可以重建不同方向的折射率实部衰减率导数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310110905.XA CN103365067B (zh) | 2012-04-01 | 2013-04-01 | 可实现三维动态观测的光栅剪切成像装置和方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210096606 | 2012-04-01 | ||
CN201210096606.0 | 2012-04-01 | ||
CN2012100966060 | 2012-04-01 | ||
CN201310110905.XA CN103365067B (zh) | 2012-04-01 | 2013-04-01 | 可实现三维动态观测的光栅剪切成像装置和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103365067A CN103365067A (zh) | 2013-10-23 |
CN103365067B true CN103365067B (zh) | 2016-12-28 |
Family
ID=49359104
Family Applications (7)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310110905.XA Expired - Fee Related CN103365067B (zh) | 2012-04-01 | 2013-04-01 | 可实现三维动态观测的光栅剪切成像装置和方法 |
CN201310111218.XA Expired - Fee Related CN103356223B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的 ct 成像系统及方法 |
CN201310111203.3A Expired - Fee Related CN103365068B (zh) | 2012-04-01 | 2013-04-01 | 光栅剪切三维成像系统及光栅剪切三维成像方法 |
CN201310111040.9A Expired - Fee Related CN103356207B (zh) | 2012-04-01 | 2013-04-01 | 基于光栅剪切成像的医学检测设备和方法 |
CN201310111197.1A Expired - Fee Related CN103364839B (zh) | 2012-04-01 | 2013-04-01 | 基于光栅剪切成像的安检设备及方法 |
CN201310110920.4A Expired - Fee Related CN103364418B (zh) | 2012-04-01 | 2013-04-01 | 光栅剪切二维成像系统及光栅剪切二维成像方法 |
CN201310111387.3A Expired - Fee Related CN103356208B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的二维成像系统及方法 |
Family Applications After (6)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310111218.XA Expired - Fee Related CN103356223B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的 ct 成像系统及方法 |
CN201310111203.3A Expired - Fee Related CN103365068B (zh) | 2012-04-01 | 2013-04-01 | 光栅剪切三维成像系统及光栅剪切三维成像方法 |
CN201310111040.9A Expired - Fee Related CN103356207B (zh) | 2012-04-01 | 2013-04-01 | 基于光栅剪切成像的医学检测设备和方法 |
CN201310111197.1A Expired - Fee Related CN103364839B (zh) | 2012-04-01 | 2013-04-01 | 基于光栅剪切成像的安检设备及方法 |
CN201310110920.4A Expired - Fee Related CN103364418B (zh) | 2012-04-01 | 2013-04-01 | 光栅剪切二维成像系统及光栅剪切二维成像方法 |
CN201310111387.3A Expired - Fee Related CN103356208B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的二维成像系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (7) | CN103365067B (zh) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103365067B (zh) * | 2012-04-01 | 2016-12-28 | 中国科学院高能物理研究所 | 可实现三维动态观测的光栅剪切成像装置和方法 |
CN104586415B (zh) * | 2013-10-31 | 2019-10-08 | Ge医疗系统环球技术有限公司 | 准直器对准偏差确定方法及计算机化断层成像系统 |
US10393681B2 (en) | 2014-02-14 | 2019-08-27 | Canon Kabushiki Kaisha | X-ray Talbot interferometer and X-ray Talbot interferometer system |
CN104132953B (zh) * | 2014-08-01 | 2017-03-29 | 中国科学技术大学 | 一种双能x射线相位衬度成像装置及其实现方法 |
CN106999137B (zh) * | 2014-11-24 | 2021-03-05 | 皇家飞利浦有限公司 | 用于x射线相衬断层合成成像的探测器和成像系统 |
JP6460226B2 (ja) * | 2015-04-09 | 2019-01-30 | 株式会社島津製作所 | X線撮影装置 |
KR20170087320A (ko) | 2016-01-20 | 2017-07-28 | 삼성전자주식회사 | 단층 영상 생성 장치 및 그에 따른 단층 영상 복원 방법 |
CN107024490B (zh) * | 2016-01-29 | 2019-07-05 | 中国科学院高能物理研究所 | 一次曝光光栅剪切成像装置及数据采集与信息提取方法 |
CN105935297A (zh) * | 2016-06-23 | 2016-09-14 | 中国科学院深圳先进技术研究院 | X射线光栅相衬成像ct系统 |
CN106197553A (zh) * | 2016-08-26 | 2016-12-07 | 中云智慧(北京)科技有限公司 | 核生化检测装置及安检装置 |
KR102260802B1 (ko) * | 2017-01-05 | 2021-06-07 | 제너럴 일렉트릭 캄파니 | 단층촬영 재구성에서 사용하기 위한 데이터의 딥러닝 기반 추정 |
JP6943090B2 (ja) * | 2017-09-05 | 2021-09-29 | 株式会社島津製作所 | X線イメージング装置 |
CN107748341A (zh) * | 2017-10-23 | 2018-03-02 | 中国科学院苏州生物医学工程技术研究所 | 高衬度低剂量相位衬度ct成像装置 |
CN107966113B (zh) * | 2017-12-04 | 2023-08-11 | 济南大学 | 一种模体成像试验台控制方法 |
CN108169257B (zh) * | 2017-12-08 | 2019-08-13 | 中国科学院高能物理研究所 | 具有三种定量成像机制的纳米ct三维成像方法 |
EP3498171A1 (en) * | 2017-12-15 | 2019-06-19 | Koninklijke Philips N.V. | Single shot x-ray phase-contrast and dark field imaging |
CN108042110A (zh) * | 2017-12-22 | 2018-05-18 | 深圳先进技术研究院 | 多模成像系统 |
CN109631765B (zh) * | 2018-12-07 | 2020-12-04 | 李里 | 影像位移传感器及其测量方法 |
US10881371B2 (en) * | 2018-12-27 | 2021-01-05 | Medtronic Navigation, Inc. | System and method for imaging a subject |
CN110133010B (zh) * | 2019-04-04 | 2020-10-27 | 中国科学技术大学 | 一种x射线相位衬度成像方法 |
CN110133012B (zh) * | 2019-07-02 | 2022-01-18 | 合肥工业大学 | 基于三探测器光栅干涉仪的单次曝光多模式x射线成像方法 |
CN110833427B (zh) * | 2019-11-29 | 2021-01-29 | 清华大学 | 光栅成像系统及其扫描方法 |
CN114460805B (zh) * | 2020-10-21 | 2024-05-28 | 中国科学院国家空间科学中心 | 一种基于高通滤波的遮挡物散射成像系统 |
CN113729593B (zh) * | 2021-09-28 | 2022-11-01 | 上海交通大学 | 基于多角度散射随机矩阵的3d内窥镜用血流成像方法 |
CN115171464B (zh) * | 2022-08-26 | 2024-05-03 | 首都医科大学宣武医院 | 一种可视化穿刺模具 |
CN115684222B (zh) * | 2022-12-21 | 2023-04-11 | 济南汉江光电科技有限公司 | 一种快速低剂量的x射线多模态ct系统及成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4841555A (en) * | 1987-08-03 | 1989-06-20 | University Of Chicago | Method and system for removing scatter and veiling glate and other artifacts in digital radiography |
CN101467889A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院高能物理研究所 | 光栅剪切相位衬度ct成像数据采集和重建方法 |
CN101532969A (zh) * | 2007-11-23 | 2009-09-16 | 同方威视技术股份有限公司 | X射线光栅相衬成像系统及方法 |
US20110243302A1 (en) * | 2010-03-30 | 2011-10-06 | Fujifilm Corporation | Radiation imaging system and method |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1447046A1 (en) * | 2003-02-14 | 2004-08-18 | Paul Scherrer Institut | Apparatus and method to obtain phase contrast x-ray images |
EP1731099A1 (en) * | 2005-06-06 | 2006-12-13 | Paul Scherrer Institut | Interferometer for quantitative phase contrast imaging and tomography with an incoherent polychromatic x-ray source |
CN100457040C (zh) * | 2005-11-17 | 2009-02-04 | 中国科学院高能物理研究所 | 同步辐射x射线相位衬度ct成像装置及实验方法 |
CN101011253B (zh) * | 2006-02-01 | 2011-06-15 | 西门子公司 | 产生投影或断层造影的相位对比照片的焦点-检测器装置 |
EP1879020A1 (en) * | 2006-07-12 | 2008-01-16 | Paul Scherrer Institut | X-ray interferometer for phase contrast imaging |
EP2060909B1 (en) * | 2007-11-15 | 2011-09-07 | CSEM Centre Suisse d'Electronique et de Microtechnique SA - Recherche et Développement | Interferometer device and method |
CN201191275Y (zh) * | 2007-11-23 | 2009-02-04 | 同方威视技术股份有限公司 | 一种x射线光栅相衬成像系统 |
WO2009076700A1 (en) * | 2007-12-14 | 2009-06-25 | Commonwealth Scientific And Industrial Research Organisation | Phase-contrast imaging method and apparatus |
JP2010164373A (ja) * | 2009-01-14 | 2010-07-29 | Konica Minolta Medical & Graphic Inc | X線撮影装置、およびx線撮影方法 |
CN101943668B (zh) * | 2009-07-07 | 2013-03-27 | 清华大学 | X射线暗场成像系统和方法 |
JP5586899B2 (ja) * | 2009-08-26 | 2014-09-10 | キヤノン株式会社 | X線用位相格子及びその製造方法 |
JP5438649B2 (ja) * | 2010-03-26 | 2014-03-12 | 富士フイルム株式会社 | 放射線撮影システム及び位置ずれ判定方法 |
KR101180067B1 (ko) * | 2010-05-25 | 2012-09-05 | 단국대학교 산학협력단 | 엑스선 직렬 격자 간섭계 |
CN103365067B (zh) * | 2012-04-01 | 2016-12-28 | 中国科学院高能物理研究所 | 可实现三维动态观测的光栅剪切成像装置和方法 |
CN202794067U (zh) * | 2012-09-19 | 2013-03-13 | 同方威视技术股份有限公司 | 一种行李物品ct安检系统及其探测器装置 |
-
2013
- 2013-04-01 CN CN201310110905.XA patent/CN103365067B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111218.XA patent/CN103356223B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111203.3A patent/CN103365068B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111040.9A patent/CN103356207B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111197.1A patent/CN103364839B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310110920.4A patent/CN103364418B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111387.3A patent/CN103356208B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4841555A (en) * | 1987-08-03 | 1989-06-20 | University Of Chicago | Method and system for removing scatter and veiling glate and other artifacts in digital radiography |
CN101532969A (zh) * | 2007-11-23 | 2009-09-16 | 同方威视技术股份有限公司 | X射线光栅相衬成像系统及方法 |
CN101467889A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院高能物理研究所 | 光栅剪切相位衬度ct成像数据采集和重建方法 |
US20110243302A1 (en) * | 2010-03-30 | 2011-10-06 | Fujifilm Corporation | Radiation imaging system and method |
Non-Patent Citations (1)
Title |
---|
Study of OSEM with different subsets grating-based X-ray differential phase-contrast imaging;Kai Zhang等;《Analytical and Bioanalytical Chemistry》;20110528;第401卷(第3期);第837-843页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103364839A (zh) | 2013-10-23 |
CN103356207B (zh) | 2016-08-03 |
CN103356207A (zh) | 2013-10-23 |
CN103365068B (zh) | 2016-03-09 |
CN103364418A (zh) | 2013-10-23 |
CN103365067A (zh) | 2013-10-23 |
CN103356208B (zh) | 2015-01-14 |
CN103356208A (zh) | 2013-10-23 |
CN103365068A (zh) | 2013-10-23 |
CN103364418B (zh) | 2016-08-03 |
CN103356223A (zh) | 2013-10-23 |
CN103364839B (zh) | 2015-12-09 |
CN103356223B (zh) | 2015-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103365067B (zh) | 可实现三维动态观测的光栅剪切成像装置和方法 | |
EP2453226B1 (en) | X-ray dark-field imaging system and method | |
US8121249B2 (en) | Multi-parameter X-ray computed tomography | |
CN106725565B (zh) | 一种稀疏投影下的锥束xct成像质量评估方法 | |
CN105575455B (zh) | 一种x射线衰减器设计方法和应用以及利用该方法设计的带有衰减器的ct装置 | |
CN104622492A (zh) | 一种x射线光栅相位衬度成像装置和方法 | |
CN107144581A (zh) | 基于横向错位吸收光栅的x射线光栅差分相位衬度成像方法及装置 | |
CN102906783A (zh) | 边缘保护的噪声滤波 | |
Bliznakova et al. | A software platform for phase contrast x-ray breast imaging research | |
US20200380738A1 (en) | Three-Dimensional Reconstruction Method Based on Half-Peak Probability Density Distribution | |
CN105852895A (zh) | 单次曝光的硬x射线光栅干涉仪的信息提取新方法 | |
US11464470B2 (en) | X-ray backscatter systems and methods for performing imaging tomosynthesis | |
CN101785029A (zh) | 根据能量相关投影数据对感兴趣区域进行成像的成像系统 | |
CN106618623B (zh) | 一次曝光的硬x射线光栅干涉仪的成像方法 | |
US20220221413A1 (en) | Detection scheme for x-ray small angle scattering | |
Pietersoone et al. | Combining Wave and Particle Effects in the Simulation of X-ray Phase Contrast—A Review | |
Nguyen et al. | The development of Radon transforms associated to Compton scatter imaging concepts | |
Li et al. | Detection performance study for cone-beam differential phase contrast CT | |
Pelliccia et al. | Comparison of different numerical treatments for x-ray phase tomography of soft tissue from differential phase projections | |
Fasina | Implementation and Validation of a GPU-based X-ray Diffraction Monte Carlo Simulator for in-vivo Breast imaging applications | |
Jadick et al. | Cramér-Rao lower bound in the context of spectral x-ray imaging with propagation-based phase contrast | |
Suhonen | Scattering and refraction as contrast mechanisms in x-ray imaging | |
Ulseth | Parallelized X-Ray Tracing with GPU Ray-Tracing Engine | |
Niknami et al. | Design and simulation of a semiconductor detector-based Compton imaging system with efficiency analysis | |
Agasthya | Low-dose imaging of liver diseases through neutron stimulated emission computed tomography: simulations in GEANT4 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161228 Termination date: 20180401 |