[发明专利]基于生物组织特异性的光学三维成像方法

权利要求书

1.基于生物组织特异性的光学三维成像方法，包括：

(1)采集数据

利用多模态分子影像系统，依次采集用于光学三维成像的多角度荧光数据、用于光学特性参数重建的多角度激光数据和用于获取生物体解剖结构的计算机断层成像投影数据；

(2)预处理

2a)利用多模态分子影像系统中的预处理软件对采集的多角度激光数据和荧光数据进行去除背景噪声、提取感兴趣区域、补偿坏点预处理；

2b)利用多模态分子影像系统中的预处理软件对计算机断层成像投影数据进行补偿坏点坏线、亮场校正、几何校正预处理；

(3)获取解剖结构

利用3DMED软件对预处理后的计算机断层成像投影数据进行三维重建，获得生物体三维体数据；利用3DMED软件中的人机交互式分割方法对获得的生物体三维体数据进行器官分割，获取生物体解剖结构；

(4)获取表面光学数据

对步骤(3)获取的生物体解剖结构和步骤(2)获取的荧光和激光数据，应用非接触式光学断层成像方法中的生物体表面三维能量重建技术获取生物体表面的三维荧光和激光数据分布；

(5)重建光学特性参数

对步骤(3)获取的生物体解剖结构和步骤(4)获取的生物体表面三维激光数据分布，应用基于区域的扩散光学层析成像方法重建生物体各个组织的光学特性参数；

(6)建立光传输模型

6a)根据生物体在解剖结构和组织光学特性参数方面的差异，将生物组织划分为高散射特性组织、低散射特性组织、空腔组织三大类；

6b)根据扩散近似方程、简化球谐波近似方程、朗伯源特性方程的适用范围，采用相应的方程描述不同特性生物组织中的光传输过程；

6c)构造外推边界和折射率不匹配边界条件耦合不同特性生物组织的光传输方程，建立统一的、描述光在整个生物体中传输过程的混合数学模型；

(7)建立系统方程

7a)如果是在第一级离散网格上建立系统方程，利用Amira软件对生物体的高散射特性和低散射特性组织进行初始离散；否则，利用编写的程序对上一级离散网格进行手动调整；

7b)在第k级离散网格上，利用多级自适应有限元方法对步骤6c)建立的混合数学模型进行离散化，组装每个离散点上的子系统方程建立总的系统方程：

A_kS_k＝Φ_k

其中，A_k是第k级离散网格上的系统矩阵，依赖于生物体内三类特性生物组织的分布和生物组织的光学特性参数；

S_k是第k级离散网格上的靶向目标能量密度分布；

Φ_k是第k级离散网格边界节点上的光通量密度；

(8)建立目标函数

根据第k级离散网格边界节点上的光通量密度计算值和测量值之间的误差，结合靶向目标的稀疏分布约束，建立下列目标函数：

minΘ(Sk)Skinf≤Sk≤Sksup={12||AkSk-Φkm||1+λk||Sk||1}]]>

其中，Θ(S_k)是第k级离散网格上的目标函数；

是第k级离散网格上的靶向目标能量密度的下限；

是第k级离散网格上的靶向目标能量密度的上限；

是第k级离散网格边界节点上的光通量密度的测量值；

‖F‖₁定义为求解矩阵F的l₁范数；

λ_k是第k级离散网格上的正则化因子；

(9)求解目标函数

9a)采用基于任务导向的混合优化方法求解建立的目标函数，根据第k级离散网格上形成的系统矩阵的大小选择合适的优化方法，获得第k级离散网格上的靶向目标能量密度分布；

9b)利用靶向目标能量密度分布计算第k级离散网格边界节点上的光通量密度；

9c)判断第k级离散网格边界节点上的光通量密度的测量值和计算值之差，如果小于给定阈值，则目标函数求解过程结束，获得光学三维成像的靶向目标重建结果，转向步骤(10)；否则，转向步骤9d)；

9d)根据步骤9a)获得的靶向目标能量密度分布和步骤9b)获得的边界节点上的光通量密度计算值，调整第k+1级离散网格，转向步骤(7)；

(10)显示结果。