一种基于自适应体素细分的光学断层成像方法
【技术领域】
[0001] 本发明属于光学技术领域，尤其设及一种基于自适应体素细分的光学断层成像方 法。
【背景技术】
[0002] 作为一种新兴的光学成像技术，光学Ξ维成像是通过融合生物体体表测量的多角 度光学信号、生物体的解剖结构和组织光学参数信息，基于精确的生物组织中的光传输模 型重建活体生物体体内祀向目标的位置和强度分布信息，其中，生物组织中光传输过程的 精确描述和祀向目标的准确快速重建是光学Ξ维成像方法实现的基础。在光学Ξ维成像技 术中，激发巧光断层成像、扩散光学层析断层成像和自发巧光断层成像是Ξ种主要的成像 模态。基于自发巧光断层成像模态，中科院自动化所在其专利申请文件"一种基于自适应有 限元的多光谱重建方法【申请号】2007 10302337.8,申请日：2007-12-26,授权号 CN200710302337.8，授权日2010.11.24)中提出了一种基于自适应有限元的多光谱重建方 法，该方法利用多种先验信息、多谱信息，结合自适应有限元的方法，提出后验可行光源区 域的策略，该方法能有效地提高自发巧光断层成像的重建质量和速度，在分子影像、重建算 法等领域有着重要的应用价值。但是该方法是基于四面体有限元网格，必须对组织进行四 面体剖分处理，用四面体数据构建前向矩阵，而构建四面体网格的过程比较复杂，重建流程 繁琐。而对组织采用六面体体素进行处理，则极大降低了数据处理的复杂性。六面体体素相 对于四面体有限元网格在存储空间和处理效率上更加具有优势，一般相同节点数量的情况 下，六面体体素数量是四面体体素数量的1/5,所W可W极大的降低数据存储空间。
[0003] 传统的对组织进行六面体体素处理过程，仅仅对组织进行单一类型的体素化，无 法精确的描述组织边界信息。
【发明内容】
[0004] 本发明的目的在于提供一种基于自适应体素细分的光学断层成像方法，旨在解决 传统对组织进行四面体有限元求解时，四面体剖分复杂，而对于组织进行六面体体素处理 过程，仅仅对组织进行单一类型的体素化处理，无法精确描述组织的边界信息的问题。
[0005] 本发明是运样实现的，一种基于自适应体素细分的光学断层成像方法，所述基于 自适应体素细分的光学断层成像方法包括:多模态的数据采集和自适应体素细分，对体内 光学信号重建;通过多模态分子影像系统获取投影数据和3DM抓分割器官获得的解剖结构 得到生物体表面的Ξ维巧光分布，利用获得的Ξ维分割图像进行体素的自适应细分;基于 自适应细分的离散体素结构和表面映射的巧光分布，建立系统方程，求解目标函数，实现对 生物体内祀向目标的准确、快速、高分辨率重建。
[0006] 进一步，所述基于自适应体素细分的光学断层成像方法包括W下步骤：
[0007] 步骤一，采集数据，利用分子影像系统，依次采集用于光学Ξ维成像的多角度巧光 数据、用于光学特性参数重建的多角度激光数据和用于获取动物体解剖结构的计算机断层 成像投影数据；
[0008] 步骤二，图像预处理和分割，采用3DMm)软件对预处理后的计算机断层成像投影数 据进行Ξ维重建，获得动物体Ξ维体数据;应用3DM邸软件中的人机交互式分割方法对获得 的动物体Ξ维体数据进行器官分割，获取动物体解剖结构；
[0009] 步骤Ξ，光学数据Ξ维表面映射，利用步骤二获取的动物体解剖结构和步骤姨获 取的巧光数据，应用非接触式光学断层成像方法利用生物体表面Ξ维能量重建技术获取动 物体表面的Ξ维巧光数据分布；
[0010] 步骤四，对组织进行自适应体素划分；
[0011] 步骤五，建立系统方程，经过步骤四自适应体素细分，在第k个离散体素上，利用多 级自适应有限元方法组装每个离散点上的子系统方程建立总的系统方程：
[0012] AkSk=巫k;
[0013] 其中，Ak是第k个离散体素上的系统矩阵，由离散化的体素结构和生物组织的光学 特性参数决定；
[0014] Sk是第k个离散体素上的祀向目标能量密度分布；
[001引 Ok是第k个离散体素边界节点上的光通量密度；
[0016] 步骤六，建立目标函数，根据第k个离散体素边界节点上的光通量密度计算值和测 量值之间的误差，结合祀向目标的稀疏分布约束，建立下列目标函数：
[0017]
[001引其中，Θ (Sk)是第k个离散体素上的目标函数；
[0019] Sf是第k个离散体素上的祀向目标能量密度的下限；
[0020] 5Γ是第k个离散体素上的祀向目标能量密度的上限；
[0021] ΦΓ是第k个离散体素边界节点上的光通量密度的测量值；
[0022] ||F| |i定义为求解矩阵F的h范数；
[0023] Ak是第k个离散体素上的正则化因子；
[0024] 步骤屯，求解目标函数。
[0025] 进一步，所述步骤一中多角度巧光数据的采集，控制成像体等间隔旋转，不大于 90%利用多模态分子影像系统中的CCD相机采集不少于四个角度的巧光数据;继续旋转，使 成像体转回到巧光数据初始采集的位置;多角度计算机断层成像投影数据的采集，控制成 像体等间隔旋转不大于1°，采集不少于360个角度的投影数据。
[00%]进一步，所述步骤四具体包括：
[0027] 首先，对组织进行初级体素离散化，遍历每个体素，并对每个体素所在组织进行编 号标记并标记出边界体素；
[0028] 将标记的边界体素利用八叉树的方法进行细分，判断细分后的次级体素是否是边 界体素，如果是边界体素，则对次级体素继续进行细分;如果不是边界体素或者已经达到最 小的分辨率，则停止细分。
[0029] 进一步，所述步骤屯具体包括：
[0030] 第一步，采用基于体素的方法求解建立的目标函数，根据第k个离散体素上形成的 系统矩阵的大小选择合适的优化方法，获得第k个离散体素上的祀向目标能量密度分布；
[0031] 第二步，利用祀向目标能量密度分布计算第k个离散体素边界节点上的光通量密 度；
[0032] 第Ξ步，判断第k个离散体素边界节点上的光通量密度的测量值和计算值之差，如 果小于给定阔值，则目标函数求解过程结束，获得光学Ξ维成像的祀向目标重建结果；否 贝1J，转向第一步。
[0033] 本发明的另一目的在于提供一种基因表达获取祀向目标的位置和强度信息的方 法，所述方法包含所述基于自适应体素细分的光学断层成像方法。
[0034] 本发明的另一目的在于提供一种肿瘤检测获取祀向目标的位置和强度信息的方 法，所述方法包含所述基于自适应体素细分的光学断层成像方法。
[0035] 本发明的另一目的在于提供一种药物研发与疗效评价获取祀向目标的位置和强 度信息的方法，所述方法包含所述基于自适应体素细分的光学断层成像方法。
[0036] 本发明提供的基于自适应体素细分的光学断层成像方法，采用基于自适应体素细 分的策略来对组织进行描述，并用自适应划分的六面体体素构建前向矩阵，不仅能更加准 确的描述组织的边界信息，同等数量的体素来表示组织边界，自适应体素的边界信息分辨 率是单一体素的3-5倍;而且最终能够获得更加精确的重建结果;提出了一种基于自适应体 素细分的光学断层成像方法，W实现对具有非规则解剖结构复杂生物体体内祀向目标的准 确、快速、高分辨率重建；同等数量的体素来对生物体内的祀向目标进行重建，自适应体素 的重建精度会比单一体素的提高30%-50%。自适应体素细分的过程，克服了均一体素无法 精确描述组织边界信息的缺点，采用自适应划分的多尺度体素来对组织进行描述，不仅实 现了对组织信息的精确定位，更加节省了数据存储空间，提高了数据处理效率;在同样节点 数目的情况下，自适应体素细分的方法的存储空间是四面体剖分的方法的Ξ分之一;如果 边界信息描述相同的情况下，自适应体素细分的方法比单一体素的存储空间降低了70% W 上。
[0037] 与现有技术相比本发明具有如下优点：
[0038] 第一，无需对组织进行四面体剖分，节省数据存储空间
[0039] 本发明是一种基于自适应体素细分的光学断层成像方法，克服了现有四面体有限 元重建技术中必须对组织进行四面体剖分才能重建W及数据存储方面的局限性，能够对具 有不规则解剖结构进行准确、快速成像。
[0040] 第二，自适应体素细分更准确的描述组织边界信息，提高重建准确性;如果是同等 数量的体素来对生物体内的祀向目标进行重建，自适应体素的重建精度会比单一体素的提 高 30%-50%。
[0041] 本发明由于对边界体素进行自适应细分，可W更加准确的描述前向传输矩阵，克 服了现有技术中体素无法精确描述边界信息的问题，能够有效地改善光学Ξ维成像的重建 精度。
【附图说明】
[0042] 图1是本发明实施例提供的基于自适应体素细分的光学断层成像方法流程图。
[0043] 图2是本发明实施例提供的基于自适应体素细分的光学断层成像方法详细步骤流 程图。
[0044] 图3是本发明实施例提供的体素细分及边界体素自适应细分示意图。
[0045] 图4是本发明实施例提供的边界体素自适应细分流程图。
【具体实施方式】
[0046] 为了使本发明的目的、技术方案及优点更加清楚明白，W下结合实施例，对本发明 进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用W解释本发明，并不用于 限定本发明。
[0047] 本发明通过多模态分子影像系统获取投影数据和3DM邸分割器官获得的解剖结构 得到生物体表面的Ξ维巧光分布，在此基础上利用获得的Ξ维分割图像进行体素的自适应 细分。基于自适应细分的离散体素结构和表面映射的巧光分布，建立系统方程，求解目标函 数，实现对生物体内祀向目标的准确、快速、高分辨率重建。
[0048] 下面结合附图对本发明的应用原理作详细的描述。
[0049] 如图1所示，本发明实施例的基于自适应体素细分的光学断层成像方法包括W下 步骤：
[0050] S101:利用分子影像系统，依次采集用于光学Ξ维成像的多角度巧光数据、用于光 学特性参数重建的多角度激光数据和用于获取动物体解剖结构的计算机断层成像投影数 据；
[0051] S102:采用3DMm)软件对预处理后的计算机断层成像投影数据进行Ξ维重建，获得 动物体Ξ维体数据;应用3DM抓软件中的人机交互式分割方法对获得的动物体Ξ维体数据 进行器官分割，获取动物体解剖结构；