基于双网格的有限投影的荧光分子断层成像重建方法
【技术领域】
[0001] 本发明属于分子影像领域，更进一步设及光学分子成像技术领域中的基于双网格 策略的有限投影角度的巧光分子断层成像重建方法，可用于在体巧光分子断层成像的逆问 题重建。
【背景技术】
[0002] 巧光分子断层成像下简称FMT)是近年发展起来的一种新型成像模态。它利用 外部光源激发巧光探针(巧光团、巧光染料等），使其发射光子，利用巧光采集装置收集巧光 信号，结合数学模型，可W获得重建目标内巧光探针的位置和浓度，W此实现对活体状态下 的生物过程进行细胞和分子水平的定性和定量研究，被广泛用于疾病早期诊断、疗效监测、 新药研发等领域。
[0003] 巧光分子断层成像属于逆向问题，具有严重病态性。由于光的强散射特性，使得光 子在其内部的传输不再沿直线传播，而是经过大量无规则的散射过程。另外，采集的巧光数 据是成像目标表面值，数量有限，运使得逆问题是一个未知数远远大于方程数的数学问题， 进一步增加了求解问题的不适定性。多点激发和多角度投影能够增加测量数据，一定程度 上缓解问题的不适定性。但由此也带来了新的问题，一方面会增加数据采集时间不利于巧 光分子断层成像的实时成像和实时重建，另一方面盲目增加测量数据会带来数据冗余。再 者，求解逆问题一般采用有限元方法，由此需要对成像目标进行网格离散求解。网格尺寸过 大，不能准确逼近成像目标;网格尺寸过小，重建问题节点数目增加，未知数增加，由此反而 增加问题的病态性。
【发明内容】
[0004] 为了解决巧光分子断层成像的实时重建问题，W及存在的病态性和解不唯一性， 本发明提出了一种基于双网格策略的有限投影角度的巧光分子断层成像重建方法。为了实 现实时重建，采用单点激发，有限角度测量。为了提高重建质量，采用双网格策略，即用网格 尺寸不同的网格一一粗网格和细网格离散成像目标，分别求解。其中粗网格上的重建结果 为细网格重建提供巧光目标的可行区域，建立表面的测量数据与可行区域内巧光目标分布 的线性关系，将线性关系转化为11范数极小化问题求解，从而获得重建目标内部的巧光目 标的Ξ维分布与浓度。
[000引1.本发明的技术方案是:一种基于双网格策略的有限投影角度的巧光分子断层 成像重建方法，其特征在于:包括W下步骤： (1) 获取有限个角度的测量数据； (2) 获得重建目标的解剖结构信息和光学特性参数； (3) 重建粗网格，基于粗网格的重建结果为细网格重建提供巧光目标的可行区域； (4) 在所述的可行区域下重建细网格，实现巧光目标的Ξ维分布。
[0006]所述的步骤（1)实现过程是:利用多模态光学分子影像系统，依次采集巧光数据， 用于光学特性参数重建的激光数据w及用于获得成像对象的解剖结构的计算机断层成像 数据，其中巧光数据为有限个角度的巧光数据。
[0007]所述步骤(2)中巧光数据的获取为有限个角度，一般为不大于四个角度的投影数 据。
[000引所述的步骤(2)包括： (la)利用3DM邸软件将计算机断层成像数据进行Ξ维重建并分割得到成像对象的解剖 结构； (化)利用激光数据和步骤(la)中的解剖结构，采用扩散光学层析成像方法重建生物体 各个组织的光学特性参数。
[0009] 所述的步骤(3)包括： (3a)利用Amira软件对成像目标进行离散获得网格尺寸较粗的网格； (3b)基于光传输模型和有限元理论，将步骤(2)获得的解剖结构信息和光学特性参数 作为先验信息，建立有限角度的巧光数据与重建对象内部巧光目标分布的线性关系； (3c)利用正则化理论，将上述线性关系转化为11范数极小化问题：
主是正则化参数，可通过不完全变量截断共辆梯度法求解，获得粗网格上巧光目标的 Ξ维分布； (3d)对步骤(2c)中的结果，将节点上分布值大于某个阔值的所有节点设定为巧光目标 可行的区域，为细网格重建提供巧光目标的可行区域。
[0010] 所述步骤(3b)中是建立表面巧光数据与成像目标内所有节点之间的线性关系，步 骤(3c)是全域重建。
[0011] 所述的步骤(4)包括： (4a)利用Amira软件对成像目标进行离散获得网格尺寸较细的网格； (4b)结合步骤(3d)给出的巧光目标可行区域，基于光传输模型和有限元理论，建立有 限角度的巧光数据与可行区域内巧光目标分布的线性关系； (4c)重复步骤(3c)，获得细网格上巧光目标的Ξ维分布； (4d)显示结果。
[0012] 所述步骤(4b)中是建立表面巧光数据与可行区域内所有节点之间的线性关系，步 骤(4c)是局部域重建。
[001引本发明具有W下优点： 第一，本发明采用的有限点激发，有限角度测量，既缩短了数据采集时间，也减小了重 建问题规模，有利于实现实时测量和实时重建； 第二，本发明利用双网格策略，采用粗细网格离散成像目标，粗网格重建结果为细网格 重建提供巧光目标的可行区域，降低了细网格上重建问题的病态性，提高了重建结果质量。
[0014]
【附图说明】
[0015] 图1为基于本发明的巧光断层成像重建方法的流程图； 图2为用于仿真实验的数字鼠模型； 图3为粗网格重建为细网格重建提供的巧光目标可行区域； 图4为用本发明的重建算法获得的重建结果，其中(a)图显示为z=16.4mm截面的巧光目 标巧光产额值分布图；（b)图显示为真实巧光目标与重建的代表巧光目标的四面体的Ξ维 展示。
[0016]
【具体实施方式】
[0017] 下面结合附图对本发明作进一步的详细描述，应指出的是，所描述的实施例仅旨 在便于对本发明的理解，而对其不起任何限定作用。
[0018] 参见图1，本发明巧光断层成像重建方法的流程如下： (1)获取有限角度的巧光数据： la)激发光源对固定在电控旋转台上的重建目标进行有限角度的透射式断层成像； 透射式断层成像，将激光器与光学检测仪器放置在成像目标的两侧，激光照射重建目 标激发巧光目标发出巧光，巧光穿透成像目标被激光器对面的光学检测仪器检测到。
[0019] 有限个角度透射式断层扫描，用电脑控制电控旋转台等间隔旋转一定角度，一般 大于90° (本例中选120°)，激光器发射点状激光照射成像目标，一般转一个角度激发一次， 采集一次巧光数据。
[0020] lb)使用光学检测仪器获得测量数据，获得光通量密度Φκ ; 在步骤la)中，将测量数据应用非接触式光学断层成像方法（授权公告号 CN101692971B)中描述的生物体表面Ξ维能量重建技术获取成像目标体表面的Ξ维巧光数 据分布。
[0021] (2)获得重建对象的解剖结构信息W及光学特性参数： 2a)重建对象的解剖结构信息 对计算机断层成像投影数据进行Ξ维重建，并用3DM抓软件预处理获得成像目标的Ξ 维体数据;采用3DM邸软件中的人机交互式半自动化分割方法对体数据进行组织分割，获得 成像目标的解剖结构； 2b)获取光学特性参数 利用解剖结构信息和应用基于生物组织特异性的光学Ξ维重建方法（申请号 201110148500.6,申请日2011.06.02)中描述的基于区域的扩散光学层析成像算法获得成 像目标内各个组织的光学特性参数。
[0022] (3)利用Amira软件对成像目标进行离散获得网格尺寸较粗的网格，网格尺寸粗细 由生成网格的