删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

微分相位衬度计算机层析成像的感兴趣区域重建方法

本站小编 Free考研考试/2021-12-29

摘要:基于光栅干涉仪系统的X射线微分相位衬度计算机层析成像, 不仅可以重建物体的线性衰减系数, 还可以重建物体的相移系数和线性散射系数. 在实际应用时, 大面积光栅不易获得, 常常遇到样品大于光栅的情况. 当用小于样品的光栅对样品进行扫描时, 样品超出光栅成像视野的部分会导致微分相位投影信息被截断. 本文针对微分相位衬度计算机层析成像提出了一种相移系数的感兴趣区域重建方法. 该方法利用物体相移系数和线性衰减系数(即折射率实部减小量和折射率虚部)之间的近似线性关系; 通过重建相移系数的Lambda函数和线性衰减系数的Lambda逆函数的多项式组合, 近似重建物体感兴趣区域的相移系数. 数值模拟实验依据菲涅耳衍射积分理论, 进行计算机仿真X射线的传播过程和光栅成像过程. 实际实验利用上海同步辐射BL13W1站的Talbot光栅干涉仪系统, 分别对标准模体和生物样品进行光栅微分相位衬度计算机层析成像. 数值模拟和实际实验结果都验证了该方法的有效性.
关键词: 微分相位衬度/
计算机层析成像/
感兴趣区域重建/
相移系数

English Abstract


--> --> -->
基于吸收的传统计算机层析成像(computed tomography, CT)重建的是物体内部线性衰减系数的分布, 对吸收较弱的低原子序数样品难以获得高衬度CT图像. X射线穿过样品时, 不仅有振幅衰减, 还有相位改变, 即相移. 对于吸收较弱的低原子序数物体, 穿过物体后的X射线相移更明显[1,2], 因此探测相移的X射线相位衬度成像方法与CT技术相结合, 可以获得高衬度三维CT图像.
X射线相位衬度成像方法中被认为最有希望推广到临床应用的是基于光栅的微分相位衬度成像[3-5]. 利用光栅微分相位衬度成像可以获得样品的吸收图像、折射图像和散射图像. 当样品相对X射线旋转时, 可以获得各个角度的吸收投影图、折射投影图和散射投影图, 重建出被测样品的线性衰减系数、相移系数和线性散射系数的三维空间分布.
光栅微分相位衬度成像的视野是限制该方法广泛应用的主要原因之一. 由于光栅制作工艺的限制, 目前制作的光栅面积在几十毫米左右. 对于大于光栅成像视野的样品进行扫描时, 在超出视野范围部分, 微分相位投影信息被截断. 此时, 虽然难以对整个样品进行微分相位衬度计算机层析成像(differential phase contrast computed tomography, DPC-CT), 但是仍然有希望应用DPC-CT对样品的感兴趣区域进行成像. 这是DPC-CT的一个重要研究方向.
传统的吸收CT感兴趣区域重建可以分为精确重建方法和近似重建方法. 反投影滤波(back-projection filter, BPF) [6,7]类的感兴趣区域重建属于精确重建. 该方法通过对感兴趣区域的投影数据微分, 然后反投影到感兴趣区域, 最后沿着感兴趣区域的PI线进行有限区间的逆Hilbert变换, 重建感兴趣区域图像. 但是有限区间的逆Hilbert变换时需要满足沿Hilbert变换方向上至少包含一个端点, 所以BPF的方法不适用于不包含图像边界的内部感兴趣区域重建问题. Lambda感兴趣区域重建[8-10]属于近似重建. 该方法通过重建样品线性衰减系数的Lambda函数和Lambda逆函数的线性组合来近似重建出感兴趣区域图像.
光栅微分相位衬度成像获得的是被测样品的吸收、微分相位和散射投影. 要从微分相位投影重建感兴趣区域的相移系数, 一些吸收CT的感兴趣重建方法难以直接应用. Anastasio等[11]虽然证明了BPF方法可用于物体局部相移系数的重建, 但是BPF方法应用的限制条件依然存在. 迭代类算法也被用于DPC-CT的感兴趣区域重建, Cong等[12]在折射率n呈分段多项式分布的前提下, 提出了一种可用于相移系数重建的迭代算法; Pascal等[13]利用感兴趣区域上的先验信息, 提出了一种可以消除截断伪影的POCS算法; Yang等[14]提出了一种基于样条函数和希尔伯特变换的迭代算法, 迭代方法的计算比较复杂, 用途有限. Felsner等[15]在2018年提出了用吸收投影数据补全相位缺失投影的方法, 先根据线性衰减系数的重建图像进行材料分解, 再通过感兴趣区域的相位重建图像估计出每种材料的折射率实部的减小量, 然后结合每种材料对应的吸收投影图与截断的微分相位投影图外推出全部的微分相位投影图, 最后进行重建; 2020年Felsner等[16]对他们在2018年所提出的这个重建算法进行了改进, 直接根据感兴趣区域的微分相位投影图估计每种材料的折射率实部的减小量, 这种方法避免了由于截断的相位重建图像造成的误差.
本文先给出了平行束扫描模式下Lambda函数和Lambda逆函数的重建公式, 并利用相移系数和线性衰减系数(即折射率实部减小量和折射率虚部)之间的近似关系, 提出了一种针对微分相位衬度CT的感兴趣区域重建方法. 这种方法同时利用了微分相位投影和吸收投影两种信息, 通过重建相移系数的Lambda函数和线性衰减系数的Lambda逆函数的多项式组合来近似重建样品的相移系数. 其中, 相移系数的Lambda函数由微分相位投影进行一阶差分重建, 线性衰减系数的Lambda逆函数由吸收投影进行反投影重建.
2
2.1.数学准备
-->利用光栅微分相位衬度成像可以获得样品的吸收投影、微分相位投影和散射投影. 吸收投影的表达式为
${p_{{\rm{abs}}}} = \int {\mu {\rm{d}}l} = 2k\int {\beta {\rm{d}}l}, $
其中: $\mu $为样品的线性衰减系数, $\mu = 2 k\beta $; l表示射线路径; k为波数, $k = {{2{\text{π}}}}/{\lambda }$, $\lambda $为波长; $\beta $为样品的折射率虚部.
微分相位投影的表达式为
$\frac{{\partial {p_{{\rm{dpc}}}}}}{{\partial r}} = \frac{\partial }{{\partial r}}\int {( - \delta ){\rm{d}}l}, $
其中, ${p_{{\rm{dpc}}}}$表示相位投影, r表示射线的法向, $\delta $为样品的折射率实部减小量, 一般我们重建的是样品的相移系数$k\delta $. 散射投影就是散射角的二阶矩${\sigma ^2}$, 其表达式为
${\sigma ^2} = \int \alpha {\rm{d}}l, $
其中, $\alpha $为样品的线性散射系数.
本文所用的$\varLambda $算子是负拉普拉斯算子的平方根[9], 即
$\varLambda = \sqrt { - \varDelta } = \sqrt { - \left( {\frac{{{\partial ^2}}}{{\partial x_1^2}} + \frac{{{\partial ^2}}}{{\partial x_2^2}}} \right)}, $
其中, $({x_1}, {x_2})$表示二维直角坐标系下的坐标.
$\delta ({x_1}, {x_2})$表示点${{x}} = ({x_1}, {x_2})$处的折射率实部减小量. $\delta ({x_1}, {x_2})$的二维Fourier变换可以表示为
$\begin{split}&F\delta = \tilde \delta ({\omega _1},{\omega _2})= \\\;& \int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{ + \infty } {\delta ({x_1},{x_2}){{\rm{e}}^{ - {\rm{i2\pi }}({x_1},{x_2}) \cdot ({\omega _1},{\omega _2})}}{\rm{d}}{x_1}{\rm{d}}{x_2}} }, \end{split}$
其中, ${\omega _1} = \omega \cos \varphi $, ${\omega _2} = \omega \sin \varphi $, ${\omega _1}, \;{\omega _2}$是二维Fourier变换后的频域直角坐标, $\omega, \;\varphi $是极坐标.
$\tilde \delta ({\omega _1}, {\omega _2})$的二维逆Fourier变换可以表示为
$\begin{split}&{F^{ - 1}}\tilde \delta = \delta ({x_1},{x_2})= \\\;& \int_{ - \infty }^{ + \infty } {\int_{ - \infty }^{ + \infty } {\tilde \delta ({\omega _1},{\omega _2}){{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2}) \cdot ({\omega _1},{\omega _2})}}{\rm{d}}{\omega _1}{\rm{d}}{\omega _2}} } . \end{split}$
$\varLambda \delta $的二维Fourier变换为
$F(\varLambda \delta ({x_1},{x_2})) = |\omega | \cdot F\delta ({x_1},{x_2}). $
${\varLambda ^{{\rm{ - }}1}}\beta $的二维Fourier变换为
$F({\varLambda ^{ - 1}}\beta ({x_1},{x_2})) = \frac{1}{{|\omega |}} \cdot F\beta ({x_1},{x_2}). $
因为光栅相位衬度成像获得的是样品的吸收投影和散射投影, 所以可以直接利用吸收CT的重建算法[17-19]对样品的线性衰减系数和线性散射系数进行重建. 但光栅相位衬度成像获得的微分相位投影是相位投影的偏导数, 所以在对样品的相移系数进行重建时, 不能直接应用吸收CT的重建方法. 对于样品相移系数的重建, 一种方法是先将微分相位投影信息积分, 得到相位投影, 之后采用线性衰减系数的重建方法重建样品的相移系数; 另一种方法是直接利用微分相位投影信息重建, 将吸收重建的滤波反投影(filter back-projection, FBP)算法中所用的斜坡滤波器, 改为希尔伯特滤波器即可进行滤波反投影重建[20]. 对微分相位投影积分和希尔伯特滤波都要求投影数据没有截断, 因此这两种方法都无法直接应用于截断数据的感兴趣区域重建.
2
2.2.感兴趣区域重建方法
-->为了简明起见, 我们只讨论平行束扫描模式下的DPC-CT的重建方法, 其结论也可推广到扇束和锥束CT. 在平行束条件下, l表示射线路径, $l:{{x}} \cdot {{\phi}} = r$, 其中${{x}} = ({x_1}, {x_2}), \;{{\phi}} = (\cos \varphi, \sin \varphi )$, r为射线的法向, $\varphi $为射线旋转角度, $\varphi \in [0, {\text{π}}]$. 沿射线l的相位投影可以表示为
${p_{{\rm{dpc}}}}(r,\varphi ) = \int { - \delta ({x_1},{x_2}){\rm{d}}l} . $
${p_{{\rm{dpc}}}}(r, \varphi )$沿r方向的一维Fourier变换为
${F_r}{p_{{\rm{dpc}}}}(\omega,\varphi ) = \int_{ - \infty }^{ + \infty } {{p_{{\rm{dpc}}}}(r,\varphi )} {{\rm{e}}^{ - {\rm{i2\pi }}\omega \cdot r}}{\rm{d}}r. $
根据FBP重建公式
$\begin{split}\delta ({x_1},{x_2})=\;&\int_0^{\rm{\pi }} \int_{-\infty}^{+\infty} {F_r}{p_{{\rm{dpc}}}}(\omega,\varphi ) \cdot |\omega |\; \\ &\times {{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2})\!\cdot\! (\omega \cos \varphi,\omega \sin \varphi )}}{\rm{d}}\omega {\rm{d}}\varphi. \end{split}$
微分相位衬度成像$\varLambda \delta $的重建公式为
$\begin{split} \;&\varLambda \delta ({x_1},{x_2}) = {F^{ - 1}}(|\omega | \cdot F\delta ({x_1},{x_2})) \\ =\;& - \int_0^{\rm{\pi }} \int_{ - \infty }^{ + \infty } {F_r}{p_{{\rm{dpc}}}}(\omega,\varphi ) \cdot |\omega | \cdot |\omega |\\&\times {{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2}) \cdot (\omega \cos \varphi,\omega \sin \varphi )}}{\rm{d}}\omega {\rm{d}}\varphi \\ =\;& \frac{1}{{4{{\rm{\pi }}^2}}}\int_0^{\rm{\pi }} \int_{ - \infty }^{ + \infty } {F_r}\left(\frac{\partial }{{\partial {r^2}}}{p_{{\rm{dpc}}}}(r,\varphi )\right)\\&\times {{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2}) \cdot (\omega \cos \varphi,\omega \sin \varphi )}}{\rm{d}}\omega {\rm{d}}\varphi \\ = \;&\frac{1}{{4{{\rm{\pi }}^2}}}\int_0^{\rm{\pi }} \int_{ - \infty }^{ + \infty } {F_r}\left(\frac{\partial }{{\partial r}}\left(\frac{{\partial {p_{{\rm{dpc}}}}(r,\varphi )}}{{\partial r}}\right)\right)\\&\times{{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2}) \cdot (\omega \cos \varphi,\omega \sin \varphi )}}{\rm{d}}\omega {\rm{d}}\varphi . \end{split} $
$\beta ({x_1}, {x_2})$表示点x处的折射率虚部, 沿射线l的吸收投影可以表示为
${p_{{\rm{abs}}}}(r,\varphi ) = \int {\mu ({x_1},{x_2}){\rm{d}}l} = 2k\int {\beta ({x_1},{x_2}){\rm{d}}l}. $
微分相位衬度成像${\varLambda ^{ - 1}}\beta $的重建公式为
$\begin{split} \;& {\varLambda ^{ - 1}}\beta ({x_1},{x_2}) = \left({F^{ - 1}}\frac{1}{{|\omega |}} \cdot F\right)\beta ({x_1},{x_2}) \\ =\;&\frac{1}{{2k}}\int_0^{\rm{\pi }} \int_{ - \infty }^{ + \infty } {F_r}{p_{{\rm{abs}}}}(\omega,\varphi ) \cdot \frac{1}{{|\omega |}} \cdot |\omega | \\ \;& \times {{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2}) \cdot (\omega \cos \varphi,\omega \sin \varphi )}}{\rm{d}}\omega {\rm{d}}\varphi \\ = \;&\frac{1}{{2k}}\int_0^{\rm{\pi }} \int_{ - \infty }^{ + \infty } {F_r}({p_{{\rm{abs}}}}(\omega,\varphi ))\\ \;& \times {{\rm{e}}^{{\rm{i2\pi }}({x_1},{x_2}) \cdot (\omega \cos \varphi,\omega \sin \varphi )}}{\rm{d}}\omega {\rm{d}}\varphi. \end{split} $
根据折射率实部减小量和折射率虚部之间存在一种线性关系, $\delta = \varepsilon \beta $[21,22], 其中$\varepsilon $为一个与材质有关的常数. 因此,
${\varLambda ^{ - 1}}\delta ({x_1},{x_2})\;=\;\varepsilon {\varLambda ^{ - 1}}\beta ({x_1},{x_2}). $
相移系数$k\delta $则可以由$\varLambda \delta $${\varLambda ^{ - 1}}\beta $的多项式组合${P_n}$近似重建.
${P_n} = \sum\limits_{j = 1}^n {\sum\limits_{i = 0}^j {{a_{ji}}(\varLambda } \delta {)^{j - i}}{{({\varLambda ^{ - 1}}\beta )}^i}}, $
其中, n表示多项式组合的最高次数, ${a_{ji}}$为各个项的系数.
多项式系数${a_{ji}}$可以通过扫描已知材质模体进行求解, 假设模体包含M种材质, 每个材质的相移系数值为$k{\delta _m}$, 所在区域为${\varOmega _m}$, $m = 1, 2, \cdot \cdot \cdot, M$. 多项式的系数可以通过求解优化问题
${{a}} = \mathop {{\rm{arg min}}}\limits_{{a}} \sum\limits_{m = 1}^M {\sum\limits_{{\rm{(}}{x_1},{x_2}{\rm{)}} \in {\varOmega _{\rm{m}}}} {\left\| {k{\delta _m} - {P_n}{\rm{(}}{x_1},{x_2}{\rm{)}}} \right\|} }, $
得到, 其中$ {{a}} = ({a_{10}}, {a_{11}}, \cdots, {a_{ji}}) $.
2
3.1.模拟实验
-->本文依据菲涅耳衍射理论, 模拟了X射线光栅干涉仪相位步进法的成像过程. 首先采用单色平行束光源照射样品, 模拟了穿过物质后X射线的传播过程; 然后通过移动吸收光栅获得每一步所探测到的光强信息; 最后利用相位步进法提取得到样品的微分相位投影和吸收投影[23-26].
模拟实验中, 所用的X射线能量为30 keV, 探测器为$205$个探元的线阵探测器, 每个探元为0.08 mm. 成像示意图如图1所示, 样品的尺寸远大于光栅尺寸. 两个光栅的尺寸均为16.4 mm, 相位光栅的周期为0.016 mm, 吸收光栅的周期为0.008 mm, 两个光栅之间的距离为$7.743 \times {10^2}\;{\rm{mm}}$. 实验样品的直径为82 mm, 由大小不同的圆和椭圆组成, 模拟时离散的图像尺寸为$2048 \times 2048$, 每个像素大小为0.04 mm; 黑色部分的材质为PMMA, 折射率实部减小量$\delta $$2.968 \times {10^{ - 7}}$, 折射率虚部$\beta $$1.02 \times {10^{ - 10}}$, 其余白色部分的$\delta $$\beta $均为$0$; 红色虚线内为光栅成像能够覆盖的感兴趣区域(region of interest, ROI), 直径为16.4 mm. 图2(a)为信息提取所得到的吸收投影的正弦图, 图2(b)为微分相位投影的正弦图, 扫描角度为360°, 扫描角度数为720个.
图 1 模拟实验成像示意图
Figure1. Imaging schematic diagram of simulation experiment.

图 2 (a) 感兴趣区域吸收投影的正弦图; (b) 感兴趣区域微分相位投影的正弦图
Figure2. (a) Sinogram of absorption projection for the ROI; (b) sinogram of differential phase projection for the ROI.

对感兴趣区域正弦图采用本文方法进行重建, 重建结果见图3. 图3(a)为采用一次多项式近似重建的结果图, 图3(b)为采用二次多项式近似重建的结果图. 多项式各项的系数采用样品自身的相移系数值, 通过(16)式和(17)式获得, 一次多项式为${P_1} = {{1}}{{.29}}\varLambda \delta +{{ 0}}{{.28}}{\varLambda ^{ - 1}}\beta $, 二次多项式为${P_2} = $$- \; {{0}}{{.56}}\varLambda \delta - {{1}}{{.26}}{\varLambda ^{ - 1}}\beta +{{ 0}}{{.28\;(}}\varLambda \delta {)^2} +{{ 1}}{{.61}}{\varLambda ^{ - 1}}\beta \varLambda \delta\; +$$ {{ 2}}{{.05(}}{\varLambda ^{ - 1}}\beta {)^2} $. 图3(c)图3(a)图3(b)在绿色虚线位置处的剖线图, 蓝色曲线表示采用一次多项式近似重建的结果, 红色曲线表示采用二次多项式近似重建的结果, 黑色曲线表示理论值. 从图3(c)可以看出采用二次多项式近似重建的结果更接近理论值. 从表1的均方误差(mean-square error, MSE)和峰值信噪比(peak signal-to-noise ratio, PSNR)的值也可以看出二次多项式近似重建结果优于一次多项式近似重建结果.
图 3 相移系数重建结果 (a) 采用一次多项式近似的重建图像; (b)采用二次多项式近似的重建图像; (c) 图3(a)图3(b)在绿色虚线位置处的剖线图
Figure3. Reconstruction results of phase shift coefficient: (a) The reconstruction image using a first order polynomial approximation; (b) the reconstruction image using a second order polynomial approximation; (c) Fig.3 (a) and Fig.3 (b) in the green dotted line location profile chart.

方法MSEPSNR
一次多项式0.079610.9894
二次多项式0.027115.6647


表1一次多项式和二次多项式重建结果的MSE和PSNR
Table1.MSE and PSNR of reconstruction results of the first order polynomial and the second order polynomial.

2
3.2.实际实验
-->本文实验采用上海同步辐射BL13 W1站的Talbot光栅干涉仪系统, 数据采集中所用射线的光子能量为20 keV, 探测器为2048 × 2048个探元的面阵探测器, 探元大小为6.5 μm. 实验所用相位光栅的周期为2.396 μm, 成像面积直径为20 mm左右, 吸收光栅的周期为2.4 μm, 成像面积直径为20 mm左右, 两块光栅之间的距离为46.38 mm. 为了更好的验证本文方法重建的效果, 实验采用对完整的投影数据截断的方式, 获得感兴趣区域投影数据.
第一组实验采用的模体由LDPE, PMMA, PTFE三种成分不同的圆柱以及水组成, 将这三个圆柱和水放在一个外直径为10.2 mm的聚乙烯塑料试管内. 图4(a)是模体的照片. LDPE, PMMA, PTFE三种成分和水的折射率实部减小量$\delta $表2所示. 实验扫描角度为180°, 角度采样数540个. 采用8步相位步进法, 通过信息提取可获得样品的吸收投影数据和微分相位投影数据. 取一个断层的投影数据, 两点合并后的探测器单元个数为1024个, 截取的覆盖感兴趣区域的探测器单元个数为512个.
图 4 相移系数重建结果 (a)实验模体; (b)全局数据重建图像, 红色虚线内为感兴趣区域图像; (c)截断数据采用一次多项式近似的重建图像; (d)截断数据采用二次多项式近似的重建图像; (e) 图4(c)图4(d)在绿色虚线位置处的剖线图
Figure4. Reconstruction results of phase shift coefficient: (a) Experimental modle; (b) the reconstruction image of global data, the ROI image is in the red dotted line; (c) the reconstruction image of the truncated data using a first order polynomial approximation; (d) the reconstruction image of the truncated data using a second order polynomial approximation; (e) Fig.4 (c) and Fig.4 (d) in the green dotted line location profile chart.

材料水(H2O)PTFE
(C2F4)
PMMA
(C5O2H8)
LDPE
(C2H4)
$\delta $/10–75.269.656.305.46


表2水、PTFE、PMMA、LDPE的折射率实部减小量$\delta $
Table2.The decrement of the real part of the refractive index of water, PTFE, PMMA, and LDPE.

截断数据的感兴趣区域重建和全局数据重建的结果见图4. 全局数据重建采用的是将FBP重建算法中的滤波器改为Hilbert滤波器的方式[20], 重建结果见图4(b), 图中圆环为聚乙烯塑料试管, 内部不同灰度的小圆为不同成分的圆柱截面, 其余灰色部分为水, 红色虚线内为截断的感兴趣区域图像. 图4(c)图4(d)是采用本文方法对截断数据的重建结果, 图4(c)为采用一次多项式近似重建的结果图, 图4(d)为采用二次多项式近似重建的结果图. 多项式系数利用感兴趣区域内已知材质的相移系数值, 通过(16)式和(17)式获得, 一次多项式为${P_1} = 0.08\varLambda \delta +0.79{\varLambda ^{ - 1}}\beta $, 二次多项式为${P_2} = $$0.76\varLambda \delta \;+\; 0.35{\varLambda ^{ - 1}}\beta\; -\; 0.15\;{{\rm{(}}\varLambda \delta )^2} - 0.56{\varLambda ^{ - 1}}\beta \varLambda \delta\; +$$0.61{{\rm{(}}{\varLambda ^{ - 1}}\beta )^2} $. 图4(e)图4(c)图4(d)在绿色虚线位置处的剖线图, 其中蓝色曲线为一次多项式近似重建的结果, 绿色曲线为二次多项式近似重建的结果, 黑色曲线为理论值.
对比图4(b)图4(d)的重建结果, 可以看出本文方法可以很好地重建出感兴趣区域图像, 与全局数据重建的感兴趣区域图像相比, 对比度略低. 由于LDPE材质的$\delta $$5.46 \times {10^{ - 7}}$, 水的$\delta $$5.26 \times $$ {10^{ - 7}}$二者之间相差不大, 图4(b)中LDPE和水无法区分开, 在图4(c)图4(d)中LDPE和水稍有区别. 从图4(e)可以看出采用二次多项式近似重建结果比一次多项式近似重建结果更接近理论值.
第二组实验样品是一只仓鼠的前趾, 长度为5 mm左右, 实验参数与第一组实验相同. 选取光栅成像视野的1/4作为感兴趣区域, 探测器单元个数为256. 图5为全局吸收投影和微分相位投影的正弦图, 其中两条红色虚线间为截取的感兴趣区域正弦图, 图5(a)为吸收投影的正弦图, 图5(b)为微分相位投影的正弦图.
图 5 (a) 全局吸收投影的正弦图; (b)全局微分相位投影的正弦图. 两红色虚线间为截断的感兴趣区域正弦图
Figure5. (a) Sinogram of global absorption projection; (b) sinogram of the global differential phase projection. Between the two red dotted line for ROI of truncation sinogram.

截断数据的感兴趣区域重建和全局数据重建的结果见图6. 图6(a)为全局数据重建的结果图, 图6(b)为全局数据的感兴趣区域重建图像; 图6(c)为截断数据的感兴趣区域重建结果图, 重建过程中采用的是一次多项式${P_1} = \varLambda \delta + 0.58{\varLambda ^{ - 1}}\beta $, 由于生物软组织的主要元素为C, H, O, 其折射率实部减小量与水的比较接近, 因此多项式系数利用实验一中水的相移系数值, 通过(16)式和(17)式获得. 由于本文方法利用了微分相位投影和吸收投影两种信息, 对比图6(b)图6(c)可以看出本文方法重建图像比全局的微分相位重建图像细节更加丰富.
图 6 相移系数重建图像 (a) 全局数据重建图像, 红色虚线内为感兴趣区域图像; (b) 全局数据的感兴趣区域重建图像; (c) 截断数据采用一次多项式近似的重建图像
Figure6. Reconstruction image of phase shift coefficient: (a) The reconstruction image of global data, the ROI image is in the red dotted line; (b) the reconstruction image of the ROI from the global data; (c) the reconstruction image of the truncated data using a first order polynomial approximation.

本文通过重建相移系数的Lambda函数和线性衰减系数的Lambda逆函数的多项式组合, 来近似重建物体的相移系数分布. 由于Lambda函数和Lambda逆函数的重建具有局部性, 因此该方法可用于截断数据的感兴趣区域重建. 模拟实验和生物样品的重建结果都表明了该方法可以对超过光栅成像视野大小的样品进行感兴趣区域重建.
相关话题/图像 数据 实验 信息 计算机

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 大功率热平衡感应耦合等离子体数值模拟及实验研究
    摘要:感应耦合等离子体发生器是“临近空间高速目标等离子体电磁科学实验研究装置”的核心部件之一,常用于模拟高焓高速等离子体鞘套环境,为了研究大功率射频中压下感应耦合等离子体发生器的放电特性,采用数值模拟和实验相结合的方法研究其内部的传热与流动特性.本文基于局域热力学平衡条件,通过湍流场-电磁场-温度场 ...
    本站小编 Free考研考试 2021-12-29
  • 激光加速高能质子实验研究进展及新加速方案
    摘要:利用超强激光与等离子体相互作用来加速高能离子是激光等离子体物理及加速器物理领域的研究热点.经过了近20年的发展,激光离子加速已取得丰硕成果,催生了一批新的应用.本文概述了国内外激光离子加速所取得的标志性实验研究进展,围绕高能质子的产生这一关键问题进行了深入的探讨,介绍了近几年来发展的有潜力的新 ...
    本站小编 Free考研考试 2021-12-29
  • 飞秒超强激光驱动太赫兹辐射特性的实验研究
    摘要:强太赫兹源是太赫兹科学技术发展的关键,其中大能量强场太赫兹脉冲源在超快物态调控、新型电子加速器等领域具有重要的应用前景.超快超强激光与等离子体相互作用是近年来发展起来的一种新型的强场太赫兹辐射产生途径.本文报道了利用超强飞秒激光脉冲与金属薄膜靶作用产生太赫兹辐射的实验结果,研究了激光能量和离焦 ...
    本站小编 Free考研考试 2021-12-29
  • 态选择电荷交换实验测量以及对天体物理软X射线发射模型的检验
    摘要:在高温天体等离子体环境中,低能高电荷态离子与中性原子和分子之间的电荷交换是天体物理环境中软X射线发射的重要机制之一.电荷交换软X射线发射相关的天体物理建模需要大量的主量子数n和角量子数l分辨的态选择俘获截面数据,目前这类数据主要来自于经典或者半经典的原子碰撞理论模型.本文利用反应显微成像谱仪, ...
    本站小编 Free考研考试 2021-12-29
  • 超弹性球体入水过程空泡演化及球体变形实验
    摘要:超弹性材料是工程实际中常用的材料,具有很强的非线性力学性能.将超弹性材料应用于入水问题是一个新的跨学科研究方向.与传统的刚性球体入水现象不同,超弹性球体入水后极易发生变形.为了探究该大变形的入水流固耦合问题,本文采用高速摄像方法,对超弹性球体垂直入水问题开展实验研究.基于实验结果,对比分析了球 ...
    本站小编 Free考研考试 2021-12-29
  • 一种基于离散数据从局部到全局的网络重构算法
    摘要:网络的结构和功能彼此相互影响,网络上的功能往往体现为网络上的动力学过程,网络上的动力学过程通过网络中的行为表象数据进行体现.因此,根据网络上可观测的相关数据对网络结构进行重构将成为可能.本文拟解决如何根据网络上可观测的离散数据还原网络拓扑结构的问题,提出了在网络局部利用每一条离散数据对应节点的 ...
    本站小编 Free考研考试 2021-12-29
  • 原子系综中光学腔增强的Duan-Lukin-Cirac-Zoller写过程激发实验
    摘要:原子系综中的Duan-Lukin-Cirac-Zoller(DLCZ)过程是产生光与原子(量子界面)量子关联和纠缠的重要手段.当一束写光与原子发生作用时,将会产生斯托克斯(Stokes)光子的自发拉曼散射,并同时产生一个自旋波(spin-wave)存储在原子系综中,上述过程即为DLCZ量子记忆 ...
    本站小编 Free考研考试 2021-12-29
  • 基于YOLOv3框架的高分辨电镜图像原子峰位置检测
    摘要:高分辨电镜图像中原子峰位置的检测具有十分重要的现实意义,通过精确定量化原子峰位置可以分析物质在微观尺度上的结构形变、电极化矢量分布等重要信息.近年来深度学习技术在图像目标检测领域取得了巨大突破,这一技术可用在高分辨电镜图像处理上,因为原子位置的检测可以看作是一个目标检测问题.本文利用先进的机器 ...
    本站小编 Free考研考试 2021-12-29
  • 基于DNA折纸模板的铁原子阵列构建及其信息加密应用
    摘要:在后摩尔时代,突破原有技术极限,进行原子尺度的精准构筑,是当前的重大科学问题.DNA作为具有原子级精准度的生物大分子,能够进行程序性的分子识别,构筑原子数量与位置均严格确定的自组装结构,因此是进行原子制造的理想平台.本文提出基于DNA自组装折纸结构的精准定位能力,构筑铁原子阵列图案,并应用于对 ...
    本站小编 Free考研考试 2021-12-29
  • 周期调制结构平面薄膜电爆炸实验研究
    摘要:在金属层表面引入微结构以实现对Z箍缩等离子体形成和发展过程中不稳定性的调控具有重要研究价值.在“强光一号”装置上(峰值电流~1.4MA,上升时间~100ns),开展了针对具有一维周期性凹槽调制结构的金属薄膜的电爆炸实验研究.实验负载采用外推型平面结构,基底为30μm厚铝膜,刻蚀周期为2mm,刻 ...
    本站小编 Free考研考试 2021-12-29