Improved AFEM algorithm for bioluminescence tomography based on dual-mesh alternation strategy

Adaptive finite element method (AFEM) is broadly adopted to recover the internal source in biological tissues. In this letter, a novel dual-mesh alternation strategy (dual-mesh AFEM) is developed for biolumi- nescence tomography. By comprehensively considering the error estimation of the finite element method solution on each mesh, two different adaptive strategies based on the error indicator of the reconstructed source and the photon flux density are used alternately in the process. Combined with the constantly adjusted permissible region in the adaptive process, the new algorithm can achieve a more accurate source location compared with the AFEM in the previous experiments.     

全变分约束迭代滤波反投影CT重建

重建算法是计算机断层成像(CT)技术的核心。在解析法CT重建过程中,结合先验信息和引入优化约束条件较为困难。通过对滤波反投影(FBP)原理及其重建图像与理想CT图像差值关系的分析,构造了以FBP为基础的迭代循环,解决了解析重建过程中先验信息的利用和优化约束条件的引入问题。为抑制迭代FBP产生的图像伪影,将全变分(TV)模型引入重建过程,建立了TV约束迭代滤波反投影CT重建方法。在数值模拟中,针对完善投影数据、稀疏投影数据、含金属投影数据和有限角投影数据等不同情况,重建出了与原始模型高度一致的CT图像,研究表明TV约束迭代滤波反投影方法是一种精度高、适应性较强的CT重建方法。     

Development of a DPET System Using in Tumor Therapy with Heavy Ions

Positron Emission Tomography (PET) is a functional imaging modality used in brain research to map in lively neurotransmitter and receptor activity and to investigate glucose utilization or blood flow patterns both in healthy and disease states. We are developing a Double-sector PET (DPET) system, which is employed as a tumor imaging detector for the heavy-ion therapy at Institute of Modern Physics (IMP).     

宽场光学相干断层成像系统的实验研究和模拟

改进和完善宽场光学相干断层成像技术(WFOCT)的主要方法是从光源、系统光路等方面进行改进.研究WFOCT系统的成像状态和环境状态对系统的横向分辨率的影响,借助ZEMAX光学设计软件模拟该系统的干涉效果,在实验中采用八步移相法来解析样品的二维图像信息,将解析出来的样品图像赋值到图像变量中进行存储显示.实验与模拟结果表明:模拟仿真获得的干涉信号含有对比度和平行度很高的干涉条纹,证明了利用ZEMAX软件模拟仿真干涉系统的可行性.     

基于非凸L_(1-2)正则子的锥束X射线发光断层成像EI北大核心CSCD

锥束X射线发光断层成像(CB-XLCT)是一种新型分子影像模态,对疾病的早期检测、靶向治疗以及药物研制等具有重要意义。然而,通过传统的压缩感知理论反演生物体内纳米目标的三维分布时,高维系统矩阵的强相关性会直接影响成像质量。基于非凸稀疏L_(1-2)正则子,将CB-XLCT的成像问题转化为一种新的稀疏重建模型。采用一种凸差分算法来解决非凸泛函最小化问题,在每一步凸差分子迭代中采用一种带自适应惩罚项的交替方向乘子法进行高效求解。设计了单目标数字鼠仿体、双目标数字鼠仿体以及真实在体老鼠实验验证提出算法的有效性和稳健性,并与五种常见正则子(L_(1/2),L_1,L_2,TV和L_0)进行对比和分析。实验结果表明,L_(1-2)正则子的成像性能最优,提出方法可以有效解决CB-XLCT的快速成像问题。     

Effects of cartilaginous rings on airflow and particle transport through simplified and realistic models of human upper respiratory tracts

In the present study, computational fluid dynamics （CFD） is used to investigate inspiratory and expiratory airflow characteristics in the human upper respiratory tract for the purpose of identifying the probable locations of particle deposition and the wall injury. Computed tomography （CT） scan data was used to reconstruct a three dimensional respiratory tract from trachea to first generation bronchi. To compare, a simplified model of respiratory tract based on Weibel was also used in the study. The steady state results are obtained for an airflow rate of 45 L/min, corresponding to the heavy breathing condition. The velocity distribution, wall shear stress, static pressure and particle deposition are compared for inspiratory flows in simplified and realistic models and expiratory flows in realistic model only. The results show that the location of cartilaginous rings is susceptible to wall injury and local particle deposition.     

基于块稀疏贝叶斯的生物发光断层重建

生物发光断层成像(BLT)是一种非侵入、高灵敏度的光学分子影像技术,可以通过探测生物体表面的光信号重建出生物体内部光源的三维分布情况。由于光在组织中传播时,散射占据主导作用,导致BLT重建问题的病态性,给光源重建带来巨大的挑战。在BLT重建中,基于光源稀疏分布的特征,稀疏正则化方法相比于传统的L2范数正则化取得了显著进展。更进一步,由于生物发光光源的分布具有的空间聚集特征,利用该特征将有助于进一步提高BLT重建的准确性。相比于传统的针对求解域中所有未知量进行稀疏重建的算法,探索了利用块稀疏进行生物发光断层成像重建的可行性,首先通过对系统矩阵进行相关系数分析将求解域划分成一系列数据块,然后利用块稀疏贝叶斯算法对生物发光光源的分布进行三维重建。通过仿真实验与小鼠活体实验,并与传统稀疏重建算法L1-LS进行了比较,结果表明该方法可以有效缓解BLT重建问题的病态性,抑制噪声,并且可提高重建结果的准确性。     

基于加权总差分最小化的中子稀疏投影计算机断层重建方法

为提升高噪声稀疏角度投影条件下中子计算机断层扫描(CT)质量,提出同时迭代重建方法(SIRT)与加权总差分最小化(WTDM)相结合的迭代重建方法(SIRT-WTDM)。在有无噪声情况下比较代数重建算法、联合代数重建算法及同时迭代重建算法的重建图像,证明了SIRT迭代重建具有较高的图像重建精度与较强的抗噪声性能,因此将SIRT作为高噪声中子投影图像CT迭代重建算法的保真项。考虑到对图像梯度稀疏性与连续性的约束,中子CT迭代重建方法的正则化约束项采用WTDM方法。由Shepp-Logan模体与真实冷中子层析扫描数据验证可知,在极端稀疏角度投影条件下,SIRT-WTDM可获得较好的重建效果。     

一种优化区域知识先验的稀疏角锥束X射线发光断层成像方法

锥束X射线发光断层（CB-XLCT）成像是一种可在生物体外对早期肿瘤进行有效检测的新型医学成像技术。稀疏角CB-XLCT成像加速了CB-XLCT技术的实时成像转化进程。然而,相对于传统多角度成像,稀疏角CB-XLCT成像的逆问题病态性明显加剧,这对传统成像方法的有效扩展提出了挑战。基于稀疏非凸L_p（01范数和L₀范数的代表算法验证所提方法的有效性和稳健性。实验结果表明,所提方法不仅可有效求解稀疏角CB-XLCT成像逆问题,还具有良好的可扩展性。