论文部分内容阅读
生物发光断层成像(BLT)是通过荧光测量值对生物组织发光源进行重建,得到发光源的三维分布。BLT包括前向问题和逆问题两部分,在BLT前向问题中,需要解决的主要问题是光传输方式的精确建模以及快速求解。在几种光传输模型中,SPN模型实用性越来越高。但是在SPN模型中,由于N的上升,导致该模型的求解计算速度不断降低。而在逆问题过程中,由于使用大规模荧光数据有助于提高重建质量,但同时也会消耗大量的计算时间,因此针对大规模数据,需要进行快速重建方面的研究。本文分别对BLT前向的SPN模型中计算代价和后向问题中重建速度两方面进行了探索。主要工作包括:1)基于GPU的BLT前向问题的并行加速求解策略和实现。在有限元方法下,对BLT前向过程进行计算时间分析可得,刚度矩阵生成和线性方程组求解两部分在整个前向过程中耗时最高,且这两部分都非常适合并行化处理,因此将SPN模型作为光传输模型,结合GPU的并行处理能力提出一种基于GPU/CPU双平台的BLT前向问题加速策略。将前向过程分解成刚度矩阵生成、线性方程组求解、数据交换和条件判断四部分,其中将前两部分传输给GPU端进行加速运算,而后两部分交给CPU执行。数字鼠仿真实验验证了在本策略下,针对不同网格数量下的SP5模型的BLT前向整体加速比最高能达到20倍左右。2)针对大规模数据导致BLT重建耗时大的问题,提出了一种基于交替方向乘子法(ADMM)结合随机对偶坐标下降方法(SDCA)的BLT重建算法。通过SDCA方法随机选取某一坐标系方向将目标函数分解成多个子函数,使每个坐标下都有对应的子函数,然后再计算每个子函数的最优解,然后在迭代过程中按照最优解与真实值之间的接近程度依次选取对应坐标下的子函数,减少无关数据在迭代过程中的使用率,提升收敛速度,实现快速重建目的。数字鼠仿体实验和真实老鼠实验验证了在保证BLT重建图像精度的同时,重建速度提高了5倍左右。