CN103356223B - 用于人体医学检测的 ct 成像系统及方法 - Google Patents
用于人体医学检测的 ct 成像系统及方法 Download PDFInfo
- Publication number
- CN103356223B CN103356223B CN201310111218.XA CN201310111218A CN103356223B CN 103356223 B CN103356223 B CN 103356223B CN 201310111218 A CN201310111218 A CN 201310111218A CN 103356223 B CN103356223 B CN 103356223B
- Authority
- CN
- China
- Prior art keywords
- grating
- projection
- data
- sample
- splitting optical
- 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
Abstract
本发明公开了一种用于人体医学检测的CT成像系统及方法,包括:光源装置;具有狭缝的扇形光束产生装置;分束光栅;样品台,用于承载人体;分析光栅,用于产生成像所要求的背景光强,增强或抑制样品的折射信号或散射信号;探测器,用于探测不同的背景光强变化,采集所述样品在不同背景光强下的投影数据;旋转装置,用于承载并旋转光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器。上述的光栅剪切CT成像系统及方法能够快速采集图像,并且密度分辨率高,密度不均匀性分辨率高,满足人体医学检测的应用需求;而且样品既可以在分束光栅前,也可以在分束光栅后,样品所受到的辐射剂量较低。
Description
技术领域
本发明涉及医学成像技术领域,特别是涉及一种用于人体医学检测的CT成像系统及方法。
背景技术
伦琴在1895年发现了X射线,并于1901年12月10日荣获第一届诺贝尔物理学奖。广为流传的伦琴夫人手的X射线照片揭示了X射线具有强大的穿透力,表明X射线直接成像就可以看到样品的内部结构。这种基于物质对X射线吸收差异的成像机制在上世纪五十年代广泛用于人体医学成像,并于上世纪八十年代开始用于人体三维成像。虽然这种基于X射线吸收机制的成像技术,在观察重元素构成物品时,可以获得衬度(即对比度)足够高的图像,但是在观察轻元素构成物品时,仅能获得模糊的图像。其主要原因在于轻元素原子所含电子数少,轻元素构成物品密度差别小,变动范围在1%—5%之间,不但对X射线吸收弱,而且对X射线吸收差别小,不能形成足够高的衬度。因而基于传统吸收衰减的成像机制在检查人体骨骼时,可以获得衬度足够高的影像,然而在检查由轻元素构成的肿瘤时,不能获得衬度足够高的影像。
X射线相位衬度成像研究始于上世纪九十年代,到目前已经有二十余年。X射线相位衬度成像在探测轻元素构成物质方面,其探测灵敏度比X射线吸收成像高得多,在医学成像方面具有广阔的发展前景。已经发展了四种X射线相位衬度成像方法,其中利用光栅提取样品相位信息的光栅剪切成像方法最具有实际应用的可能性,其最大的优势在于可以和常规X射线光源结合。在X射线光栅剪切成像研究中,研究人员还发现了散射成像机制,它是由样品中很多微小颗粒对X射线的多重折射形成,这种成像机制 对样品中微孔、微泡、微粒、微晶和粉末等结构比较敏感。
目前利用光栅扫描提取相位信息和散射信息是国际上发展的主流,然而,光栅扫描的方法不符合医学成像简便快速的要求。
发明内容
本发明的一个目的是提供一种用于人体医学检测的CT成像系统,以便实现简便快速成像,满足医学检测的应用需求。
本发明的另一个目的是提供一种用于人体医学检测的CT成像方法,以便实现简便快速成像,满足人体医学检测应用需求。
本发明的用于人体医学检测的CT成像系统,包括:
光源装置,用于产生一维多缝X射线光源;
具有狭缝的扇形光束产生装置,用于在X射线光源产生的光束的照射下产生扇形光束,并照射分束光栅;
分束光栅,用于将所述扇形光束分束为一维周期性光束阵列;
样品台,用于承载人体,并沿着垂直于扇形光束扇面的方向移动人体;
分析光栅,用于产生成像所要求的背景光强,增强或抑制样品的折射信号或散射信号;
探测器,用于探测不同的背景光强变化,采集所述样品在不同背景光强下的投影数据;
旋转装置,用于承载光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器,使光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器在保持相对位置不变的情况下围绕人体在与扇面平行的平面内旋转。
本发明的用于人体医学检测的CT成像方法,包括:
调整光源装置,使所述光源装置产生一维多缝光源;
调整具有狭缝的扇形光束产生装置:使所述狭缝产生扇形光束照射分 束光栅,并根据成像的不同需要,使狭缝限制的扇形光束或仅照射分析光栅三部分中的一部分,或照射分析光栅三部分中的两部分,或照射整个分析光栅;
调整分束光栅,使分束光栅平面垂直于所述光束中心传播方向,使分束光栅将所述光束分束为一维周期性光束阵列;
调整分析光栅,使所述分析光栅对准所述分束光栅产生的一维光束阵列;
探测器采集样品的投影数据:把分析光栅和所述分束光栅产生的一维光束阵列之间的剪切位移调整在背景光强满足成像要求的采集位置,放入样品,并使扇形光束围绕样品旋转,通过探测器采集预定旋转角度范围内所述样品在不同背景光强下的投影数据。
本发明的用于人体医学检测的CT成像系统及方法只需采集一套投影数据,就能实现三维半定量成像;只需采集三套背景光强不同的投影数据,就能实现三维定量成像;与国内外目前盛行的方法比较,不仅密度分辨率高,密度不均匀性分辨率高,而且方法简便,所需采集投影数据少,样品所需的辐射剂量低,能够快速采集投影数据,快速重建样品三维结构,能够满足人体医学X射线CT检测应用需求。
附图说明
图1(a)、图1(b)为本发明用于人体医学检测的CT成像系统的结构示意图,其中,在图1(a),当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行,三条探测器能够并行采集三种不同背景光强的投影数据;在图1(b)中,具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条垂直;
图2为本发明样品对X射线光束吸收衰减作用的示意图,其中I0为入射光强,I为出射光强;
图3为本发明样品对X射线光束产生折射作用的示意图;
图4为本发明样品对X射线光束产生散射作用的示意图;
图5(a)和图5(b)为光强随分析光栅相对分束光栅产生的一维光束阵列剪切位移变化的位移曲线,图5(a)为光强随分析光栅(四条黑色)相对分束光栅产生的一维光束阵列(条纹填充)沿X轴剪切位移变化的位移曲线,(从左到右)分析光栅(四条黑色)和分束光栅一维光束阵列(条纹填充)之间剪切位移分别固定在暗场位置、左半亮场位置、亮场位置、右半亮场位置、暗场位置;图5(b)为光强随分析光栅(四条黑色)相对分束光栅产生的一维光束阵列(条纹填充)沿Y轴剪切位移变化的位移曲线,(从下到上)分析光栅(四条黑色)和分束光栅一维光束阵列(条纹填充)之间剪切位移分别固定在暗场位置、下半亮场位置、亮场位置、上半亮场位置、暗场位置;
图6为平行束投影成像几何结构,图中O为样品旋转中心;
图7为扇形光束投影成像几何结构,图中O为样品旋转中心,探测器为一维线阵像素阵列或二维平面或弧面像素阵列;
图8为一条扇形光束射线用平行束的几何参数来描述的示意图。
图9为扇形光束共轭光线之间的关系示意图。
图中标记示意为:1-扩展光源;2-光源光栅;3-具有狭缝的扇形光束产生装置;4-分束光栅;5-样品台;6-分析光栅;7-探测器;8-平行束射线;9-扇形光束射线;10-扇形光束光源旋转轨道。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。应当注意,为了清楚的目的,附图和说明中省略了与本发明无关的、本领域普通技术人员已知的部件和处理的表示和描述。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本 发明保护的范围。
参见图1(a),图1(b),本实施例中的用于人体医学检测的CT成像系统,包括:
光源装置,用于产生一维多缝X射线光源;
具有狭缝的扇形光束产生装置3,用于在X射线光源产生的光束的照射下产生扇形光束,并照射分束光栅4;
分束光栅4,用于将所述扇形光束分束为一维周期性光束阵列;
样品台5,用于承载人体,并沿着垂直于扇形光束扇面的方向移动人体;
分析光栅6,用于产生成像所要求的背景光强,增强或抑制样品的折射信号或散射信号;
探测器7,用于探测不同的背景光强变化,采集所述样品在不同背景光强下的投影数据;
旋转装置(图中未示出),用于承载光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器,使光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器在保持相对位置不变的情况下,围绕人体在与扇面平行的平面内旋转(图中箭头方向为旋转装置的旋转方向)。
所述X射线光源产生的X射线能量最好在40keV-100keV之间,以便更好的对人体的信息进行采集。
上述用于人体医学检测的CT成像系统只需采集一套投影数据,就能实现三维半定量成像;只需采集三套背景光强不同的投影数据,就能实现三维定量成像;与国内外目前盛行的方法比较,不仅密度分辨率高,密度不均匀性分辨率高,而且方法简便,所需采集投影数据少,样品所需的辐射剂量低,能够快速采集投影数据,快速重建样品三维结构,能够满足人体医学X射线CT检测应用需求。
可选的,所述光源装置包括点光源或缝光源;或所述光源装置包括扩展光源1和光源光栅2,或所述光源装置为具有光源光栅互补结构的栅条 靶;所述光源光栅用于将所述扩展光源分割成一维多缝光源,或所述栅条靶可以直接产生一维多缝光源。
可选的,所述样品台5设置于分束光栅4与光源装置之间并且紧邻分束光栅4设置;或所述样品台5设置于分束光栅4与分析光栅6之间并且紧邻分束光栅4设置。
可选的,所述光源光栅2贴近所述光源放置;和/或,所述光源光栅2的栅条宽度大于或等于缝宽,或所述栅条靶的栅条宽度小于或等于缝宽;
具有狭缝的扇形光束产生装置3的狭缝宽度可调;和/或所述扇形光束产生装置为具有狭缝的板结构,所述板结构的厚度为至少使透过光强衰减到入射光强的1%所需的厚度;
所述光源光栅2、所述分束光栅4和所述分析光栅6均为吸收光栅,或所述分束光栅4为相位光栅,所述光源光栅2和所述分析光栅6为吸收光栅;和/或,
所述光源光栅2或栅条靶的周期与所述分析光栅6的周期形成针孔成像关系,针孔是分束光栅4上的任意一条缝;和/或,
所述分束光栅4和所述分析光栅6之间的距离为0.1~5米;和/或,
所述分束光栅4的周期为1~100微米;和/或,
所述分束光栅4的栅条宽度和缝宽相等;和/或,
当具有狭缝的扇形光束产生装置3上的狭缝与分束光栅栅条垂直时,所述分束光栅为圆柱面光栅,圆柱的中心轴位于光源光栅或栅条靶中心;当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行时,所述分束光栅为平面光栅;
当所述分束光栅4为圆柱面光栅时,所述分析光栅6为与分束光栅平行的同心轴圆柱面光栅,产生一种背景光强;当所述分束光栅4为平面光栅时,所述分析光栅6为平行于分束光栅的平面光栅,并且分为上、中、下三部分,分别产生不同的背景光强;
所述分析光栅6的周期等于所述光源中心光线对所述分束光栅周期 的几何投影或几何投影的二分之一;和/或,
所述分析光栅6的栅条宽度和缝宽相等;和/或,
所述探测器7,当分析光栅6产生一种背景光强时,为贴近所述分析光栅的一维线阵或二维面阵探测器,或所述探测器,当分析光栅6产生三种不同的背景光强时,为三条独立的线阵探测器,分别贴近所述分析光栅6上、中、下部分放置;和/或,
所述探测器7包括多个探测单元构成的多个独立并行探测的一维线阵或二维面阵,
所述栅条靶为将靶光源与光源光栅集成为一体设置的结构。
可选的,在所述光源光栅2或分束光栅4或所述分析光栅6为吸收光栅时,其栅条厚度为至少使透过光强衰减到入射光强的10%所需的厚度;在所述分束光栅4为相位光栅时,其栅条厚度需能够使透过光束获得π或π/2的相移。
本实施例的用于人体医学检测的CT成像方法,包括:
调整光源装置,使所述光源装置产生一维多缝光源;
调整具有狭缝的扇形光束产生装置:使所述狭缝产生扇形光束照射分束光栅,并根据成像的不同需要,使狭缝限制的扇形光束或仅照射分析光栅三部分中的一部分,或照射分析光栅三部分中的两部分,或照射整个分析光栅;
调整分束光栅,使分束光栅平面垂直于所述光束中心传播方向,使分束光栅将所述光束分束为一维周期性光束阵列;
调整分析光栅,使所述分析光栅对准所述分束光栅产生的一维光束阵列;
探测器采集样品的投影数据:把分析光栅和所述分束光栅产生的一维光束阵列之间的剪切位移调整在背景光强满足成像要求的采集位置,放入样品,并使扇形光束围绕样品旋转,通过探测器采集预定旋转角度范围内所述样品在所述背景光强下的投影数据。
上述的用于人体医学检测的CT成像方法只需采集一套投影数据,就 能实现三维半定量成像;只需采集三套背景光强不同的投影数据,就能实现三维定量成像;与国内外目前盛行的方法比较,不仅密度分辨率高,密度不均匀性分辨率高,而且方法简便,所需采集投影数据少,样品所需的辐射剂量低,能够快速采集投影数据,快速重建样品三维结构,能够满足人体医学X射线CT检测应用需求。
可选的,所述预定旋转角度范围为0度旋转到度的角度范围, 为扇形光束角宽度,或所述预定旋转角度范围为0度旋转到360度的角度范围。
可选的,所述光源装置包括扩展光源和光源光栅,或所述光源装置为具有光源光栅互补结构的栅条靶,所述“调整光源装置,使所述光源装置产生一维多缝光源”具体为“调整扩展光源和光源光栅,光源光栅将所述扩展光源分割成一维多缝光源,或调整栅条靶产生的一维多缝光源,使每条缝光源都能产生光束照射分束光栅”。
可选的,所述背景光强包括:亮场背景、暗场背景和/或半亮场背景;所述半亮场背景包括右半亮场背景和/或左半亮场背景;或包括上半亮场背景和/或下半亮场背景;
所述采集投影数据包括:所述样品在所述亮场背景下旋转的亮场投影数据、在所述暗场背景下旋转的暗场投影数据、和/或在所述半亮场背景下旋转的半亮场投影数据;所述半亮场投影数据包括右半亮场投影数据和/或左半亮场投影数据,或者包括上半亮场投影数据和/或下半亮场投影数据。
可选的,所述调整分析光栅之后,还包括:在样品台上无样品的条件下,通过探测器阵列探测背景光强的变化,在垂直于扇形光束中心传播方向的平面内沿着垂直于光栅栅条的方向移动所述光源光栅或栅条靶或分束光栅或分析光栅,调整分析光栅和分束光栅产生的一维光束阵列之间的剪切位移,通过探测器的各探测单元采集随剪切位移变化的光强,形成光强随剪切位移变化的位移曲线,并用余弦函数曲线拟合位移曲线;
在探测器采集样品的投影数据之后还包括根据拟合的位移曲线,从所述采集的投影数据中提取所述样品的半定量或定量描述信息,然后进行三 维重建的步骤。
可选的,从所述采集的投影数据中提取所述样品的半定量或定量描述信息,具体包括:
(h1)建立光栅剪切成像方程:用余弦函数曲线拟合测得的位移曲线,建立物函数数学模型、根据物函数和拟合位移曲线的卷积运算,建立光栅剪切成像方程;
(h2)求得所述采集投影数据的数学表达式:根据光栅剪切成像方程分别求得所述亮场投影数据、暗场投影数据和半亮场投影数据的数学表达式;
(h3)提取样品的半定量投影数据:分别把所述亮场投影数据、暗场投影数据和半亮场投影数据的数学表达式进行变形,求得线性吸收系数投影数据、折射率实部衰减率导数投影数据、线性散射系数投影数据或线性消光系数投影数据的半定量表达式;
或,
(h4)提取所述样品的定量投影数据:根据所述亮场投影数据、暗场投影数据和半亮场投影数据的数学表达式之间的定量关系,获得所述样品的线性吸收系数投影数据、折射率实部衰减率导数投影数据、线性散射系数投影数据或线性消光系数投影数据的定量表达式;
(h5)建立重建算法公式:根据傅立叶中心切片定理,建立线性吸收系数、折射率实部衰减率和折射率实部衰减率导数、线性散射系数或线性消光系数的重建算法公式;
(h6)三维重建:根据重建算法公式,重建样品的线性吸收系数、折射率实部衰减率和折射率实部衰减率导数、线性散射系数或线性消光系数的三维分布。
本发明实施例中,扇形光束中心传播方向为Z方向,在垂直扇形光束中心传播方向的平面内,平行于扇形光束扇面的方向为X方向,垂直于扇形光束扇面的方向为Y方向,光栅栅条或平行于扇形光束扇面或垂直于扇 形光束扇面。
例如,光栅剪切成像方法还可包括半定量成像方法和/或定量成像方法。本发明提出的半定量成像方法中,采集一套投影数据,便可获得或与线性吸收系数投影数据、或与折射率实部衰减率导数投影数据、或与线性散射系数投影数据或与线性消光系数投影数据明显相关的半定量投影数据;定量成像方法中,至多采集三套不同背景光强的投影数据,便可从中提取样品的线性吸收系数投影数据、折射率实部衰减率导数投影数据、线性散射系数投影数据或线性消光系数等定量投影数据。
步骤(h1)中,建立光栅剪切成像方程的过程为:
第一步,数学描述所述分析光栅对所述分束光栅产生的一维光束阵列的滤波作用,求出描述成像系统特性的脉冲响应函数;第二步,建立样品对X射线作用的数学模型,写出物函数的数学表达;第三步,计算物函数和脉冲响应函数的卷积,建立光栅剪切成像方程。
数学描述所述分析光栅对所述分束光栅产生的一维光束阵列的滤波作用。因为分束光栅产生的一维光束阵列和分析光栅都是一维周期函数,具有相同的周期,分析光栅相对分束光栅产生的一维光束阵列的剪切位移可以调整两者之间的相关性,所以分析光栅对分束光栅产生的一维光束阵列的滤波作用在数学上是互相关运算。
在图1(a)或图1(b)描绘的光栅剪切成像系统中,各光栅栅条与X或Y轴平行,当沿着垂直于栅条的方向移动光源光栅或栅条靶或分束光栅或分析光栅时,就会引起分束光栅产生的一维光束阵列和分析光栅之间发生剪切位移,探测器每个像素(或称为探测单元)都可以测到背景光强随剪切位移而变化的位移曲线。因为每个像素测得的位移曲线都相同,在成像面满足平移不变性,所以位移曲线就是成像系统的脉冲响应函数。图5(a)或图5(b)为光栅栅条平行于Y或X轴方向时,分束光栅产生的一维光束阵列和分析光栅之间沿X或Y轴方向进行剪切位移时的位移曲线;因为位移曲线形状类似余弦曲线,为了利用余弦曲线的对称性质,简化提取折射和散 射信息的数学表达,所以用余弦曲线拟合位移曲线S(θg),其表达式为:
其中或为光栅栅条平行于Y或X轴方向时,分析光栅相对分束光栅在X或Y轴方向的剪切角位移,xg或yg为分析光栅相对分束光栅产生的一维光束阵列在X或Y轴方向的剪切位移,D为分束光栅和分析光栅之间在光束传播方向上的距离,p为分析光栅的周期,也是位移曲线的周期,为无样品时位移曲线平均值,Smax和Smin分别为位移曲线的最大值和最小值,为无样品时位移曲线的可见度。图5(a)中位移曲线上a点对应亮场,d点对应暗场,b点对应左半亮场,c点对应右半亮场。图5(b)中位移曲线上a点对应亮场,d点对应暗场,b点对应下半亮场,c点对应上半亮场。亮场代表分束光栅产生的一维光束阵列几乎全部通过分析光栅,暗场代表分束光栅产生的一维光束阵列几乎全被分析光栅阻挡,半亮场代表分束光栅产生的一维光束阵列中,一半被分析光栅阻挡,一半通过分析光栅。
第二步:建立物函数数学模型。在建立物函数数学模型之前,先对样品中一点进行定义。在二维成像中,样品所在物面上一点(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)为样品的线性吸收系数。(2)式的物理意义为,吸收引起光线强度衰减,但不改变光线方向。(2)式还可以表示为分量表达式:
折射是一个能量守恒的过程,如图3所示,描述了样品对入射X射线折射作用的图像。根据图3,样品中一点(x,y)对通过该点光线的折射可以表达为:
上式左边表示入射光束,右边表示出射光束,其中表示光束角度矢量,为折射角矢量,其数学表达为:
其中δ(x,y,z)为样品折射率实部衰减率。(5)式的物理意义为,折射改变光线方向,但不改变光线强度。(5)式还可以写为分量表达式:
散射由面积元内部很多小颗粒的多重折射引起,也是一个能量守恒的过程,如图4所示,描述了样品对入射X射线散射作用的图像。散射和折射的不同之处在于,折射把样品物面上一个面积元作为一个整体来研究,即把样品物面上一个面积元作为一个微小棱镜,而散射则研究这个面积元内部的不均匀性质,即相当于研究微小棱镜内部的气泡、颗粒、微孔、微晶和杂质等不均匀。因此,对于每个面积元,只有一条折射光线和一个折射角,却有多条散射光线和多个散射角。换言之,散射是一个光束分散的过程。因为样品有一定厚度,在面积元内部沿着光束传播方向,各小颗粒分布是随机的,前后两个小颗粒产生的折射是相互独立的,小颗粒每次折射使入射光线偏离入射方向的角度是随机的,所以根据中心极限定理,散射角是以入射角(或者折射角)为中心服从二维正态统计分布,可以用方差来描述散射角分布范围。根据图4,一光线射入样品时,由于散射引起分散,出射光线分为两部分,散射光线和未散射光线,未散射光线仍然沿着入射方向传播,而散射光线偏离入射方向传播。随着光线在样品中穿行,散射事件的不断发生,散射光线不断产生,散射能量逐渐增强,而未散射光线逐渐消弱,未散射能量逐渐减弱,称为消光衰减。需要特别说明的是,每条光线都可能遇到面积元内部多个小颗粒的折射,需要把第一次被小颗粒折射产生第一次偏离的散射光线和该散射光线继续被后续小颗粒折射产生进一步偏离的散射光线分别考虑,这是因为散射能量是由一次折射产生的一次偏离决定的,而以后的多次折射产生的多次偏离仅仅使散射能量分布范围更大,增加散射角方差,而对增加或减少散射能量几乎不起作用。简言之,小颗粒一次折射决定散射能量和未散射能量的比例,小颗粒多次折射决定散射角方差。因此,消光衰减和吸收衰减一样遵循比尔定律。设入射光线能量为1,未散射光线继续沿着入射方向传播,其所携带能量,即消光衰减为exp(-Γ(x,y)),散射光线能量为1-exp(-Γ(x,y))。在散射中心对称的条件下,样品一点(x,y)对通过该点光线的散射可以表达为:
上式左边表示入射光束,右边表示出射光束,其中表示光束角度矢量。(8)式还可以写为分量表达式:
在(8)式或(9)式右边第一项中
其中γ(x,y,z)为样品的线性消光系数,第二项中σ2(x,y)为(x,y)点处样品整体厚度产生的散射角方差。因为样品整体厚度的散射角方差σ2(x,y),是光线传播路径上一系列厚度为Δzi薄片的微分散射角方差Δσ2(x,y,z)之和,所以样品整体厚度的散射角方差可以表示为:
其中ω(x,y,z)为线性散射系数。为了建立线性散射系数和线性消光系数之间的关系,把(11)式与(10)式进行比较,可得:
ω(x,y,z)=ε(x,y,z)γ(x,y,z), (12)
其中ε(x,y,z)为扩散因子。若样品是由散射性质相同的材料构成时,扩散因子ε(x,y,z)就是常量,则下式成立:
σ(x,y)=ε·Γ(x,y)。 (13)
此时就可以从一种信号获得另一种信号。换言之,若样品是由散射性质相同的材料构成时,则两个几何意义不同的散射信号可以归结为一个信号。
综合考虑上述三种作用,在样品散射中心对称的条件下,样品中一点(x,y)对通过该点光线的作用可以用物函数表达,
其分量表达式为:
根据(15)式,可知仅在X方向起作用的物函数为:
仅在Y方向起作用的物函数为:
根据(14)式或(15)式,可知出射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或Y轴方向作用的物函数Ox,y(x,y,θg)和位移曲线S(θg)的卷积推导而出:
其中,I0为无样品时分束光栅前的入射光光强,exp(-Μ(x,y))为吸收衰减像,θx,y(x,y)为折射角像,其中θg为分析光栅相对分束光栅沿X或Y轴方向的剪切角位移,V(x,y)为放入样品后位移曲线的可见度,又称为样品的可见度像,其表达式为:
为无样品时位移曲线的可见度,exp(-Γ(x,y))为样品的消光 衰减像,σ2(x,y)为样品的散射角方差像。
步骤(h2)中,数学表达采集投影数据的过程为:
设所述分析光栅相对所述分束光栅产生的一维光束阵列的剪切角位移θg为:
即在实验中把所述分束光栅产生的一维光束阵列和所述分析光栅之间的剪切位移固定在亮场位置,背景为亮场,放入样品,根据(18)式,样品相对于光束逐步或连续旋转,探测器逐幅采集到的亮场投影数据 可以表达为:
设所述分析光栅相对所述分束光栅产生的一维光束阵列的剪切角位移θg为:
即在实验中把所述分束光栅产生的一维光束阵列和所述分析光栅之间的剪切位移固定在暗场位置,背景为暗场,放入样品,根据(18)式,样品相对于光束逐步或连续旋转,探测器逐幅采集到的暗场投影数据 可以表达为:
设所述分析光栅相对所述分束光栅产生的一维光束阵列的剪切角位 移θg为:
即在实验中把所述分束光栅产生的一维光束阵列和所述分析光栅之间的剪切位移固定在右半亮场位置或上半亮场位置,背景为右半亮场或上半亮场,放入样品,根据(18)式,样品相对于光束逐步或连续旋转,探测器逐幅采集到的右半亮场投影数据可以表达为:
上半亮场投影数据可以表达为:
设所述分析光栅相对所述分束光栅产生的一维光束阵列的剪切角位移θg为:
即在实验中把所述分束光栅产生的一维光束阵列和所述分析光栅之间的剪切位移固定在左半亮场位置或下半亮场位置,背景为左半亮场或下半亮场,放入样品,(18)式,样品相对于光束逐步或连续旋转,探测器逐幅采集到的左半亮场投影数据可以表达为:
下半亮场投影数据可以表达为:
步骤(h3)中,提取样品半定量投影数据的过程为:
在忽略样品折射和散射的条件下,根据(20)式或(21)式,线性吸收系数投影数据的半定量表达式为:
或
在忽略样品吸收和散射的条件下,根据(22)式和(24)式,折射率实部衰减率平行于扇形光束扇面导数的投影数据的半定量表达式为:
或
根据(23)式和(25)式,折射率实部衰减率垂直于扇形光束扇面导数的投影数据的半定量表达式为:
或
在忽略样品吸收和折射条件下,
根据(20)式或(21)式,所述可见度像的半定量表达式为:
或
在弱散射条件下,
有:
把(34)式代入(19)式,得
把(32)式或(33)式代入(35)式,得线性散射系数投影数据的半定量表达式为:
或
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性消光系数投影数据的半定量表达式为:
或
在强散射条件下,
有:
和,
把(41)式代入(19)式,得:
把(32)式或(33)式代入(42)式,得线性消光系数投影数据的半定量表达式为:
或
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性散射系数投影数据的定量表达式为:
或
步骤(h4)中,提取样品定量投影数据的过程为:
根据(20)式和(21)式,或(22)式和(24)式,或(23)式和(25)式,获得线性吸收系数投影数据的定量表达式:
或
或从下列方程组获得:
将同一方向的亮场投影数据和暗场投影数据、或右半亮场投影数据/上半亮场投影数据和左半亮场投影数据/下半亮场投影数据,按照相应像素一一对准,并根据所述公式进行加法和对数运算;
根据(20)式、(21)式、(22)式和(24)式,获得折射率实部衰减 率平行于扇形光束扇面导数投影数据的定量表达式:
根据(20)式、(21)式、(23)式和(25)式,获得折射率实部衰减率垂直于扇形光束扇面导数投影数据的定量表达式
将同一方向的亮场投影数据、暗场投影数据、右半亮场投影数据/上半亮场投影数据和左半亮场投影数据/下半亮场投影数据,按照相应像素一一对准,并根据所述公式进行减法、除法和反正切运算;
在弱散射条件下,根据(34)式、(20)式、(21)式、(22)式和(24)式、或(23)式和(25)式,所述线性散射系数投影数据的定量表达式可从下列方程组获得:
或
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性散射系数投影数据的定量表达式可从下列方程组获得:
或
将同一方向的亮场投影数据、暗场投影数据、右半亮场投影数据/上半亮场投影数据和左半亮场投影数据/下半亮场投影数据,按照相应像素一一对准,并根据所述公式进行加法、减法、除法、乘方、开方和对数运算;
在强散射条件下,根据(41)式、(20)式、(21)式、(22)式和(24)式、或(23)式和(25)式,根据所述线性消光系数投影数据的定量表达式可从下列方程组获得:
或
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,根据所述线性散射系数投影数据的定量表达式:
或
将同一方向的亮场投影数据、暗场投影数据、右半亮场投影数据/上半亮场投影数据和左半亮场投影数据/下半亮场投影数据,按照相应像素一一对准,并根据所述公式进行加法、减法、除法、乘方、开方和对数运算。
根据关系式:
或,
同一方向的亮场投影数据、暗场投影数据、右半亮场投影数据/上半亮场投影数据和左半亮场投影数据/下半亮场投影数据中,其中任何一种都可以用其余三种来表示,因而在步骤(h4)中,可以把用四种背景投影数据表达的定量投影数据公式简化为用三种背景投影数据表达的定量投影数据公式。
当光栅栅条与样品转轴平行时,转角从0旋转到π的左半亮场投影数据和转角从π旋转到2π的右半亮场投影数据等价,或者转角从0旋转到π的右半亮场投影数据和转角从π旋转到2π的左半亮场投影数据等价,即以下两个关系式成立:
因而在步骤(h4)中,因而只需要采集从转角从0旋转到2π的右半亮场投影数据或左半亮场投影数据,就可获得转角从0旋转到π的右半亮场投影数据和左半亮场投影数据。换言之,用简单的转角从0旋转到2π的一种半亮场投影数据采集,替代复杂的转角从0旋转到π的两种半亮场投影数据采集。
步骤(h5)中,当所述X射线光束为平行束,建立重建算法公式的过程为:
设(x,y,z)为X射线光束坐标系,(x′,y′,z′)为样品坐标系,两套坐标系之间的关系为:
根据傅立叶中心切片定理,线性吸收系数的重建算法公式为
其中F-1代表一维逆傅立叶变换,为得耳塔函数;
根据傅立叶中心切片定理,折射率实部衰减率的重建算法公式为:
或,
其中i为虚数,其中F-1和分别代表一维逆傅立叶变换和二维逆傅立叶变换,(ρ,ν)为与空间坐标(x,y)对应的傅立叶频谱坐标;
根据傅立叶中心切片定理,折射率实部衰减率导数的重建算法公式为:
和,
其中 改变ψ角,可以重建不同ψ角方向的折射率实部衰减率导数;
根据傅立叶中心切片定理,在弱散射条件下,线性散射系数的重建算法公式为;
根据傅立叶中心切片定理,在强散射条件下,线性消光系数的重建算法公式为:
步骤(f6)中,当所述X射线光束为扇形光束,线性吸收系数、折射率实部衰减率、折射率实部衰减率导数、线性散射系数或线性消光系数的重建方法如下:
在扇形光束条件下,上述重建算法可以按照下列步骤转换为扇形光束重建算法。把所有扇形光束射线放在一起进行分组,把互相平行的射线分在一组,这样就把扇形光束重建问题化为平行束重建问题。设为扇形光束张角,样品旋转的投影光线集合中包含了平行束条件下样品旋转π的投影光线集合,通过重组,可以从转角区间为的扇形光束投影数据中,获得转角区间为[0,π]的平行束投影数据。因此,扇形光束重建问题可以转变为平行束重建问题。
以线性吸收系数重建为例说明,线性消光系数、线性散射系数和折射率实部衰减率及其导数的重建可如法炮制。如图6和图7所示,每一条扇形光束射线刚好对应一条平行束射线,这两条射线完全重合。每一条扇形光束射线的线性吸收系数的投影数据都刚好对应一条平行束射线的线性吸收系数的投影数据两者之间的坐标满足下列关系:
x=Lsinα (72)
其中,L为光源到旋转中心的距离,参考图8。当上述两个关系式成立时,平行束投影数据和扇形光束投影数据是相同的,有:
因此,下列方程组成立:
上述方程组中,等式左边是扇形光束投影数据,右边是平行光束投影数据。上述推理,对于探测器直接采集的投影数据也成立,因而下列方程组也成立:
在扇形光束投影数据重组为平行束投影数据的形式之后,就可以用平行束重建算法重建图像了。这种方法在投影数据重组时需要做插值运算,而插值运算会引入误差,因此这种方法有时不够精确。为了避免插值运算带来的误差,把对投影数据重组改为对重建算法的变量代换。从一个平行束重建算法的数学等式出发,其左边是待重建的图像,右边是积分表达式,积分号内有以平行束坐标变量表示的平行束投影数据接下来把平行束投影数据用等价的扇形光束投影数据来替换。在微积分中,变量替换时需要加入一个雅可比因子。雅可比因子是由一些偏导数算出的行列式,它是扇形光束坐标(α,y,β)的函数。在完成了平行束投影数据到扇形光束投影数据的替换,旧变量 到新变量的替换,并加入一个雅可比因子后,扇形光束重建算法就诞生了。
半定量重建算法公式如下:
在忽略样品折射和散射的条件下,根据(64)式和(26)式或(27)式,线性吸收系数的半定量的重建算法公式为:
其中F-1代表一维逆傅立叶变换,ρ为与空间坐标x对应的傅立叶频谱坐标,
或
在忽略样品吸收和散射的条件下,当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条垂直时,根据(65)式和(28)式或(29)式,利用折射率实部衰减率X方向导数投影数据半定量重建折射率实部衰减率的算法公式为:
其中i为虚数;根据(67)式和(28)式或(29)式,利用折射率实部衰减率X方向导数投影数据半定量重建折射率实部衰减率平行于扇形光束扇面导数的算法公式为:
改变ψ角,可以重建不同ψ角方向的折射率实部衰减率导数,
或
在忽略样品吸收和散射的条件下,当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行时,根据(66)式和(30)式或(31)式,利用折射率实部衰减率Y方向导数投影数据半定量重建折射率实部衰减率的算法公式为:
其中分别代表二维逆傅立叶变换,(ρ,ν)为与空间坐标(x,y)对应的傅立叶频谱坐标;根据(68)式和(30)式或(31)式,利用折射率实部衰减率Y方向导数投影数据半定量重建折射率实部衰减率垂直于扇形光束扇面导数的算法公式为:
其中
或
在忽略样品吸收和折射、且在弱散射条件下,根据(69)式和(36)式或(37)式,线性散射系数的半定量重建算法公式为:
其中
或
在忽略样品吸收和折射、且在强散射条件下,根据(70)式和(43)式或(44)式,线性消光系数的重建算法公式为:
其中
或
其中 α为扇形光束射线与光源焦点和样品转轴联线的夹角,β为z′轴与光源焦点和样品转轴联线的夹角,J(α,β)为雅可比因子。
定量重建算法公式如下:
线性吸收系数的重建算法公式为:
其中,
或
或从下列方程组获得:
当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条垂直时,利用折射率实部衰减率X方向导数投影数据定量重建折射率实部衰减率的算法公式为:
折射率实部衰减率平行于扇形光束扇面导数的重建算法公式为:
其中 改变ψ角,可以重建不同ψ角方向的折射率实部衰减率导数,表示平行于扇形光束扇面的折射角,可从下列方程组获得其定量表达式:
当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行时,利用折射率实部衰减率Y方向导数投影数据定量重建折射率实部衰减率的算法公式为:
折射率实部衰减率垂直于扇形光束扇面导数的重建算法公式为:
其中表示垂直于扇形光束扇面的折射角,可从下列方程组获得其定量表达式:
在弱散射条件下,线性散射系数的重建算法公式为:
其中为从以下方程组求得的解,
或
在强散射条件下,线性消光系数的重建算法公式为:
其中为从以下方程组求得的解,
或,
其中 α为扇形光束射线与光源焦点和样品转轴联线的夹角,β为z′轴与光源焦点和样品转轴联线的夹角,J(α,β)为雅可比因子。
当光栅栅条垂直于扇形光束扇面,分束光栅产生的一维光束阵列和分析光栅的剪切位移固定在半亮场位置时,样品在扇形光束中旋转2π的投 影光线集合、和在平行束中旋转2π的投影光线集合完全等价,两种光线集合之间的区别仅仅在于光线的排列顺序不同,通过重组,可以把转角区间为[0,2π]的扇形光束投影数据变换为转角区间为[0,2π]的平行束投影数据。根据(62)式和(63)式,参考图9,转角为β左半亮场投影数据和转角为β+π+2α的右半亮场投影数据等价,或者转角为β右半亮场投影数据和转角为β+π+2α的左半亮场投影数据等价,即以下两个关系式成立:
Claims (14)
1.一种用于人体医学检测的CT成像系统,其特征在于,包括:
光源装置,用于产生一维多缝X射线光源;
具有狭缝的扇形光束产生装置,用于在X射线光源产生的光束的照射下产生扇形光束,并照射分束光栅;
分束光栅,用于将所述扇形光束分束为一维周期性光束阵列;
样品台,用于承载人体,并沿着垂直于扇形光束扇面的方向移动人体;
分析光栅,用于产生成像所要求的背景光强,增强或抑制样品的折射信号或散射信号;
在光源光栅或所述分束光栅或所述分析光栅为吸收光栅时,其栅条厚度为至少使透过光强衰减到入射光强的10%所需的厚度;在所述分束光栅为相位光栅时,其栅条厚度需能够使透过光束获得π或π/2的相移;
探测器,用于探测不同的背景光强变化,采集所述样品在不同背景光强下的投影数据;
旋转装置,用于承载光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器,使光源装置、具有狭缝的扇形光束产生装置、分束光栅、分析光栅和探测器在保持相对位置不变的情况下围绕人体在与扇面平行的平面内旋转。
2.根据权利要求1所述的用于人体医学检测的CT成像系统,其特征在于,所述光源装置包括点光源或缝光源;或所述光源装置包括扩展光源和光源光栅,或所述光源装置为具有光源光栅互补结构的栅条靶;所述光源光栅用于将所述扩展光源分割成一维多缝光源,或所述栅条靶可以直接产生一维多缝光源。
3.根据权利要求1所述的用于人体医学检测的CT成像系统,其特征在于,所述样品台设置于分束光栅与光源装置之间并且紧邻分束光栅设置;或所述样品台设置于分束光栅与分析光栅之间并且紧邻分束光栅设置。
4.根据权利要求2所述的用于人体医学检测的CT成像系统,其特征在于,
所述光源光栅贴近所述光源放置;和/或,所述光源光栅的栅条宽度大于或等于缝宽,或所述栅条靶的栅条宽度小于或等于缝宽;
具有狭缝的扇形光束产生装置的狭缝宽度可调;和/或所述扇形光束产生装置为具有狭缝的板结构,所述板结构的厚度为至少使透过光强衰减到入射光强的1%所需的厚度;
所述光源光栅、所述分束光栅和所述分析光栅均为吸收光栅,或所述分束光栅为相位光栅,所述光源光栅和所述分析光栅为吸收光栅;和/或,
所述光源光栅或栅条靶的周期与所述分析光栅的周期形成针孔成像关系,针孔是分束光栅上的任意一条缝;和/或,
所述分束光栅和所述分析光栅之间的距离为0.1~5米;和/或,
所述分束光栅的周期为1~100微米;和/或,
所述分束光栅的栅条宽度和缝宽相等;和/或,
当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条垂直时,所述分束光栅为圆柱面光栅,圆柱的中心轴位于光源光栅或栅条靶中心;当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行时,所述分束光栅为平面光栅;
当所述分束光栅为圆柱面光栅时,所述分析光栅为与分束光栅平行的同心轴圆柱面光栅,产生一种背景光强;当所述分束光栅为平面光栅时,所述分析光栅为平行于分束光栅的平面光栅,并且分为上、中、下三部分,分别产生不同的背景光强;
所述分析光栅的周期等于所述光源中心光线对所述分束光栅周期的几何投影或几何投影的二分之一;和/或,
所述分析光栅的栅条宽度和缝宽相等;和/或,
所述探测器,当分析光栅产生一种背景光强时,为贴近所述分析光栅的一维线阵或二维面阵探测器,或所述探测器,当分析光栅产生三种不同的背景光强时,为三条独立的线阵探测器,分别贴近所述分析光栅上、中、下部分放置;和/或,
所述探测器包括多个探测单元构成的多个独立并行探测的一维线阵或二维面阵,
所述栅条靶为将靶光源与光源光栅集成为一体设置的结构。
5.一种用于人体医学检测的CT成像方法,其特征在于,包括:
调整光源装置,使所述光源装置产生一维多缝光源;
调整具有狭缝的扇形光束产生装置:使所述狭缝产生扇形光束照射分束光栅,并根据成像的不同需要,使狭缝限制的扇形光束或仅照射分析光栅三部分中的一部分,或照射分析光栅三部分中的两部分,或照射整个分析光栅;
调整分束光栅,使分束光栅平面垂直于所述光束中心传播方向,使分束光栅将所述光束分束为一维周期性光束阵列;
调整分析光栅,使所述分析光栅对准所述分束光栅产生的一维光束阵列;
探测器采集样品的投影数据:把分析光栅和所述分束光栅产生的一维光束阵列之间的剪切位移调整在背景光强满足成像要求的采集位置,放入样品,并使扇形光束围绕样品旋转,通过探测器采集预定旋转角度范围内所述样品在不同背景光强下的投影数据;
所述背景光强包括:亮场背景、暗场背景和/或半亮场背景;所述半亮场背景包括右半亮场背景和/或左半亮场背景,或上半亮场背景和/或下半亮场背景;
所述采集投影数据包括:所述样品在所述亮场背景下旋转的亮场投影数据、在所述暗场背景下旋转的暗场投影数据、和/或在所述半亮场背景下旋转的半亮场投影数据;所述半亮场投影数据包括右半亮场投影数据和/或左半亮场投影数据,或者包括上半亮场投影数据和/或下半亮场投影数据。
6.根据权利要求5所述的用于人体医学检测的CT成像方法,其特征在于,所述预定旋转角度范围为0度旋转到180+Δφ度的角度范围,Δφ为扇形光束角宽度,或所述预定旋转角度范围为0度旋转到360度的角度范围。
7.根据权利要求5所述的用于人体医学检测的CT成像方法,其特征在于,所述光源装置包括扩展光源和光源光栅,或所述光源装置为具有光源光栅互补结构的栅条靶,所述“调整光源装置,使所述光源装置产生一维多缝光源”具体为“调整扩展光源和光源光栅,光源光栅将所述扩展光源分割成一维多缝光源,或调整栅条靶产生的一维多缝光源,使每条缝光源都能产生光束照射分束光栅”。
8.根据权利要求5所述的用于人体医学检测的CT成像方法,其特征在于,所述调整分析光栅之后,还包括:在样品台上无样品的条件下,通过探测器阵列探测背景光强的变化,在垂直于扇形光束中心传播方向的平面内沿着垂直于光栅栅条的方向移动所述光源光栅或栅条靶或分束光栅或分析光栅,调整分析光栅和分束光栅产生的一维光束阵列之间的剪切位移,通过探测器的各探测单元采集随剪切位移变化的光强,形成光强随剪切位移变化的位移曲线,并用余弦函数曲线拟合位移曲线;
在探测器采集样品的投影数据之后还包括根据拟合的位移曲线,从所述采集的投影数据中提取所述样品的半定量或定量描述信息,然后进行三维重建的步骤。
9.根据权利要求8所述的用于人体医学检测的CT成像方法,其特征在于,从所述采集的投影数据中提取所述样品的半定量或定量描述信息,具体包括:
建立光栅剪切成像方程:用余弦函数曲线拟合测得的位移曲线,建立物函数数学模型、根据物函数和拟合位移曲线的卷积运算,建立光栅剪切成像方程;
求得所述采集投影数据的数学表达式:根据光栅剪切成像方程分别求得所述亮场投影数据、暗场投影数据和半亮场投影数据的数学表达式;
提取样品的半定量投影数据:分别把所述亮场投影数据、暗场投影数据和半亮场投影数据的数学表达式进行变形,求得线性吸收系数投影数据、折射率实部衰减率导数投影数据、线性散射系数投影数据或线性消光系数投影数据的半定量表达式;
或,
提取所述样品的定量投影数据:根据所述亮场投影数据、暗场投影数据和半亮场投影数据的数学表达式之间的定量关系,获得所述样品的线性吸收系数投影数据、折射率实部衰减率导数投影数据、线性散射系数投影数据或线性消光系数投影数据的定量表达式;
建立重建算法公式:根据傅立叶中心切片定理,建立线性吸收系数、折射率实部衰减率和折射率实部衰减率导数、线性散射系数或线性消光系数的重建算法公式;
三维重建:根据重建算法公式,重建样品的线性吸收系数、折射率实部衰减率和折射率实部衰减率导数、线性散射系数或线性消光系数的三维分布。
10.根据权利要求9所述的用于人体医学检测的CT成像方法,其特征在于,所述“提取样品的半定量投影数据”步骤中,
在忽略样品折射和散射的条件下,所述线性吸收系数投影数据的半定量表达式为:
或
在忽略样品吸收和散射的条件下,所述折射率实部衰减率平行于扇形光束扇面导数的投影数据的半定量表达式为:
或
所述折射率实部衰减率垂直于扇形光束扇面导数的投影数据的半定量表达式为:
或
在忽略样品吸收和折射、并在弱散射条件下,所述线性散射系数投影数据的半定量表达式为:
或
此时,若样品是由散射性质相同的材料构成,则所述线性消光系数投影数据的半定量表达式为:
或
在忽略样品吸收和折射、并在强散射条件下,所述线性消光系数投影数据的半定量表达式为:
或
此时,若样品是由散射性质相同的材料构成,则所述线性散射系数投影数据的半定量表达式为:
或
上述各式中, Smax和Smin分别为位移曲线的最大值和最小值,I0为无样品时分束光栅前的入射光强,为探测器探测到的样品在右半亮场背景光强下的投影数据,为探测器探测到的样品在左半亮场背景光强下的投影数据,为探测器探测到的样品在上半亮场背景光强下的投影数据,为探测器探测到的样品在下半亮场背景光强下的投影数据,为探测器探测到的样品在亮场背景光强下的投影数据,为探测器探测到的样品在暗场背景光强下的投影数据,D为分束光栅和分析光栅之间的距离,p为分析光栅和位移曲线的周期,ε为扩散因子。
11.根据权利要求9所述的用于人体医学检测的CT成像方法,其特征在于,所述“提取所述样品的定量投影数据”步骤中,
所述线性吸收系数投影数据的定量表达式为:
或
或从下列方程组获得:
所述折射率实部衰减率平行于扇形光束扇面导数投影数据的定量表达式从下列方程组获得:
所述折射率实部衰减率垂直于扇形光束扇面导数投影数据的定量表达式从下列方程组获得:
所述线性散射系数投影数据和所述线性消光系数投影数据的获取方法为:在弱散射条件下,所述线性散射系数投影数据的定量表达式从下列方程组获得:
或
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常量,所述线性消光系数投影数据的定量表达式可从下列方程组获得:
或
在强散射条件下,所述线性消光系数投影数据的定量表达式从下列方程组获得:
或
此时,若样品是由散射性质相同的材料构成,则扩散因子ε为常数,所述线性散射系数投影数据的定量表达式从下列方程组获得:
或
上述各式中, Smax和Smin分别为位移曲线的最大值和最小值,I0为无样品时分束光栅前的入射光强,为探测器探测到的样品在右半亮场背景光强下的投影数据,为探测器探测到的样品在左半亮场背景光强下的投影数据,为探测器探测到的样品在上半亮场背景光强下的投影数据,为探测器探测到的样品在下半亮场背景光强下的投影数据,为探测器探测到的样品在亮场背景光强下的投影数据,为探测器探测到的样品在暗场背景光强下的投影数据,D为分束光栅和分析光栅之间的距离,p为分析光栅和位移曲线的周期,ε为扩散因子。
12.根据权利要求10或11所述的用于人体医学检测的CT成像方法,其特征在于,所述的“建立重建算法公式”步骤中,当所述光束为扇形光束,预定旋转角度范围为0度旋转到360度的角度范围,半定量重建线性吸收系数、折射率实部衰减率、折射率实部衰减率导数、线性散射系数或线性消光系数的方法如下:
在忽略样品折射和散射的条件下,线性吸收系数的半定量重建算法公式为:
其中
或
在忽略样品吸收和散射的条件下,当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条垂直时,利用折射率实部衰减率X方向导数投影数据半定量重建折射率实部衰减率的算法公式为:
利用折射率实部衰减率X方向导数投影数据半定量重建折射率实部衰减率垂直样品转轴导数的算法公式为:
其中 改变ψ角,可以重建不同ψ角方向的折射率实部衰减率导数,
或
在忽略样品吸收和散射的条件下,当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行时,利用折射率实部衰减率Y方向导数投影数据半定量重建折射率实部衰减率的算法公式为:
用折射率实部衰减率Y方向导数投影数据半定量重建折射率实部衰减率平行样品转轴导数的算法公式为:
其中
或
在忽略样品吸收和折射、且在弱散射条件下,线性散射系数的半定量重建算法公式为:
其中
或
在忽略样品吸收和折射、且在强散射条件下,线性消光系数的半定量重建算法公式为:
其中
或
其中 α为扇形光束射线与光源焦点和样品转轴联线的夹角,β为z'轴与光源焦点和样品转轴联线的夹角,J(α,β)为雅可比因子。
13.根据权利要求10或11所述的用于人体医学检测的CT成像方法,其特征在于:当所述光束为扇形光束,预定旋转角度范围为0度旋转到360度的角度范围,当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条平行时,定量重建线性吸收系数、折射率实部衰减率、折射率实部衰减率导数、线性散射系数或线性消光系数的方法如下:
线性吸收系数的重建算法公式为:
其中
或
折射率实部衰减率的重建算法公式为:
折射率实部衰减率垂直于扇形光束扇面导数的重建算法公式为:
其中表示垂直于扇形光束扇面的折射角,从下列方程组获得:
在弱散射条件下,线性散射系数的重建算法公式为:
其中为从以下方程组求得的解,
在强散射条件下,线性消光系数的重建算法公式为:
其中为从以下方程求得的解:
其中
α为扇形光束射线与光源焦点和样品转轴联线的夹角,β为z'轴与光源焦点和样品转轴联线的夹角,J(α,β)为雅可比因子。
14.根据权利要求10或11所述的用于人体医学检测的CT成像方法,其特征在于,当所述光束为扇形光束,预定旋转角度范围为0度旋转到360度的角度范围,当具有狭缝的扇形光束产生装置上的狭缝与分束光栅栅条垂直时,定量重建线性吸收系数、折射率实部衰减率、折射率实部衰减率导数、线性散射系数或线性消光系数的方法如下:
线性吸收系数的重建算法公式为:
其中
或
或
折射率实部衰减率的重建算法公式为:
折射率实部衰减率导数的重建算法公式为:
其中 改变ψ角,重建不同ψ角方向的折射率实部衰减率导数;其中i为虚数,以下方程组获得:
在弱散射条件下,线性散射系数的重建算法公式为:
其中以下方程组获得:
在强散射条件下,线性消光系数的重建算法公式为:
其中以下方程组获得:
其中
α为扇形光束射线与光源焦点和样品转轴联线的夹角,β为z'轴与光源焦点和样品转轴联线的夹角,J(α,β)为雅可比因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310111218.XA CN103356223B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的 ct 成像系统及方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100966060 | 2012-04-01 | ||
CN201210096606 | 2012-04-01 | ||
CN201210096606.0 | 2012-04-01 | ||
CN201310111218.XA CN103356223B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的 ct 成像系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103356223A CN103356223A (zh) | 2013-10-23 |
CN103356223B true CN103356223B (zh) | 2015-07-08 |
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 | 可实现三维动态观测的光栅剪切成像装置和方法 |
CN201310111197.1A Expired - Fee Related CN103364839B (zh) | 2012-04-01 | 2013-04-01 | 基于光栅剪切成像的安检设备及方法 |
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 | 基于光栅剪切成像的医学检测设备和方法 |
CN201310110920.4A Expired - Fee Related CN103364418B (zh) | 2012-04-01 | 2013-04-01 | 光栅剪切二维成像系统及光栅剪切二维成像方法 |
CN201310111218.XA Expired - Fee Related CN103356223B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的 ct 成像系统及方法 |
CN201310111387.3A Expired - Fee Related CN103356208B (zh) | 2012-04-01 | 2013-04-01 | 用于人体医学检测的二维成像系统及方法 |
Family Applications Before (5)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310110905.XA Expired - Fee Related CN103365067B (zh) | 2012-04-01 | 2013-04-01 | 可实现三维动态观测的光栅剪切成像装置和方法 |
CN201310111197.1A Expired - Fee Related CN103364839B (zh) | 2012-04-01 | 2013-04-01 | 基于光栅剪切成像的安检设备及方法 |
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 | 基于光栅剪切成像的医学检测设备和方法 |
CN201310110920.4A Expired - Fee Related CN103364418B (zh) | 2012-04-01 | 2013-04-01 | 光栅剪切二维成像系统及光栅剪切二维成像方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
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 (24)
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射线相衬断层合成成像的探测器和成像系统 |
US10393890B2 (en) * | 2015-04-09 | 2019-08-27 | Shimadzu Corporation | X-ray imaging device |
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 | 中云智慧(北京)科技有限公司 | 核生化检测装置及安检装置 |
CN110337673B (zh) * | 2017-01-05 | 2023-10-20 | 通用电气公司 | 基于深度学习的在断层摄影重建中估计数据的方法和系统 |
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 | 李里 | 影像位移传感器及其测量方法 |
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 | 清华大学 | 光栅成像系统及其扫描方法 |
CN114460805A (zh) * | 2020-10-21 | 2022-05-10 | 中国科学院国家空间科学中心 | 一种基于高通滤波的遮挡物散射成像系统 |
CN113729593B (zh) * | 2021-09-28 | 2022-11-01 | 上海交通大学 | 基于多角度散射随机矩阵的3d内窥镜用血流成像方法 |
CN115684222B (zh) * | 2022-12-21 | 2023-04-11 | 济南汉江光电科技有限公司 | 一种快速低剂量的x射线多模态ct系统及成像方法 |
Family Cites Families (19)
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 |
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成像装置及实验方法 |
CN101011257B (zh) * | 2006-02-01 | 2011-07-06 | 西门子公司 | 产生投影或断层造影相位对比图像的焦点-检测器装置 |
EP1879020A1 (en) * | 2006-07-12 | 2008-01-16 | Paul Scherrer Institut | X-ray interferometer for phase contrast imaging |
ATE524056T1 (de) * | 2007-11-15 | 2011-09-15 | Suisse Electronique Microtech | Interferometervorrichtung und verfahren |
CN201191275Y (zh) * | 2007-11-23 | 2009-02-04 | 同方威视技术股份有限公司 | 一种x射线光栅相衬成像系统 |
CN101576515B (zh) * | 2007-11-23 | 2012-07-04 | 同方威视技术股份有限公司 | X射线光栅相衬成像系统及方法 |
US20100327175A1 (en) * | 2007-12-14 | 2010-12-30 | Yakov Nesterets | Phase-contrast imaging method and apparatus |
CN101467889B (zh) * | 2007-12-26 | 2010-08-25 | 中国科学院高能物理研究所 | 光栅剪切相位衬度ct成像数据采集和重建方法 |
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 | 富士フイルム株式会社 | 放射線撮影システム及び位置ずれ判定方法 |
JP2011224329A (ja) * | 2010-03-30 | 2011-11-10 | Fujifilm Corp | 放射線撮影システム及び方法 |
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 CN201310111197.1A patent/CN103364839B/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 CN201310110920.4A patent/CN103364418B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111218.XA patent/CN103356223B/zh not_active Expired - Fee Related
- 2013-04-01 CN CN201310111387.3A patent/CN103356208B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103365068A (zh) | 2013-10-23 |
CN103365067B (zh) | 2016-12-28 |
CN103365067A (zh) | 2013-10-23 |
CN103364418B (zh) | 2016-08-03 |
CN103365068B (zh) | 2016-03-09 |
CN103356223A (zh) | 2013-10-23 |
CN103356208A (zh) | 2013-10-23 |
CN103364839A (zh) | 2013-10-23 |
CN103356207B (zh) | 2016-08-03 |
CN103356208B (zh) | 2015-01-14 |
CN103356207A (zh) | 2013-10-23 |
CN103364839B (zh) | 2015-12-09 |
CN103364418A (zh) | 2013-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103356223B (zh) | 用于人体医学检测的 ct 成像系统及方法 | |
CN102221565B (zh) | X射线源光栅步进成像系统与成像方法 | |
CN101943668B (zh) | X射线暗场成像系统和方法 | |
CN105120755A (zh) | 用于光谱微分相衬锥形束ct和混合锥形束ct的方法和设备 | |
CN104622492A (zh) | 一种x射线光栅相位衬度成像装置和方法 | |
CN103364416B (zh) | X射线微分相衬显微成像系统及成像方法 | |
CN103858002A (zh) | 基于主成分分析的图像融合方法 | |
CN101726503A (zh) | X射线相衬层析成像 | |
Zhang et al. | Synchrotron radiation phase-contrast X-ray CT imaging of acupuncture points | |
CN103472074A (zh) | Ct成像系统和方法 | |
CN104574416A (zh) | 一种低剂量能谱ct图像去噪方法 | |
Brombal et al. | Monochromatic breast computed tomography with synchrotron radiation: phase-contrast and phase-retrieved image comparison and full-volume reconstruction | |
CN107144581A (zh) | 基于横向错位吸收光栅的x射线光栅差分相位衬度成像方法及装置 | |
CN102906783A (zh) | 边缘保护的噪声滤波 | |
CN105675631A (zh) | 一种快速扇束几何相位衬度ct成像装置和方法 | |
CN105852895A (zh) | 单次曝光的硬x射线光栅干涉仪的信息提取新方法 | |
CN101467889B (zh) | 光栅剪切相位衬度ct成像数据采集和重建方法 | |
CN104323790B (zh) | 同轴相衬成像方法及系统和相衬ct方法及系统 | |
Chen et al. | Image reconstruction for fan-beam differential phase contrast computed tomography | |
CN106618623B (zh) | 一次曝光的硬x射线光栅干涉仪的成像方法 | |
Acciavatti et al. | Oblique reconstructions in tomosynthesis. II. Super‐resolution | |
Pietersoone et al. | Combining Wave and Particle Effects in the Simulation of X-ray Phase Contrast—A Review | |
Kong et al. | Evaluation of an analytic reconstruction method as a platform for spectral cone-beam CT | |
Provinciali | X-ray phase imaging based on Hartmann Wavefront Sensor for applications on the study of neurodegenerative diseases | |
Jadick et al. | Cramér-Rao lower bound in the context of spectral x-ray imaging with propagation-based phase contrast |
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: 20150708 Termination date: 20180401 |