一种基于体素的微扰荧光蒙特卡罗模拟方法技术

本发明专利技术涉及一种基于体素的微扰荧光蒙特卡罗模拟方法，将入射光源表征为设定数目光子的集合，确定光源的初始位置和方向及探测器的位置；投放激发光光子，追踪激发光光子在生物组织中传输，计算激发光光子转化为荧光光子比例，沿着激发光路径计算探测器上所接收到的荧光光子权重,保存光子的路径信息；当荧光团吸收系数微小变化时，利用保存的光子路径信息，来直接计算探测器上荧光光子权重，大大节约了模拟光子在生物组织中传输时间。本发明专利技术方法计算效率高，能够大大提高荧光断层成像重建的计算效率。

【技术实现步骤摘要】
一种基于体素的微扰荧光蒙特卡罗模拟方法
本专利技术属于数学仿真和生物医学工程领域，涉及一种基于体素的微扰荧光蒙特卡罗模拟方法。
技术介绍
大多数生物组织均为高散射的三维浑浊介质，建立一个高精度和高效率的计算方法对于荧光断层成像的定量精度有着重要的意义。蒙特卡罗是一种基于随机抽样过程的离散统计理论方法。相比较于其它方法，蒙特卡罗方法可模拟任意几何形状，边界条件和光学参数下的光子输运过程。由于其广泛的适用性，它作为一个模拟光子输运实际物理过程的最直接，最有效和最可信的方法。因而，它成为了评价其它特定应用方法的金标准。A.J.Welch最早描述了荧光在生物组织中激发和传播的规律，并提出了标准荧光蒙特卡罗法[1]。在半无限大的层状浑浊介质中，该方法的仿真结果已经被证明是精确的[2]。Liebert提出了运用于层状浑浊介质的微扰荧光蒙特卡罗方法[3]。在荧光断层成像重建中，荧光团吸收系数改变时，往往需要进行一次完整的标准荧光蒙特卡罗仿真来得到结果，计算量往往相当大；且通过CT和MRI等成像手段提供的生物组织结构信息均是基于体素的。A.J.Welch提出的方法会给计算机造成较大的运算负担，而Li本文档来自技高网...

【技术保护点】
一种基于体素的微扰荧光蒙特卡罗模拟方法，其特征在于包括以下步骤：(1)确定目标生物组织，对目标生物组织进行三维剖分，构建一个三维体素模型，建立一个三维数字矩阵，三维数字矩阵中每个元素与三维体素模型中体素一一对应，每个元素的数值标识一种生物组织，设置生物组织的光学特性参数：吸收系数、散射系数、荧光团吸收系数、折射率和各向异性因子；(2)将入射光源表征为设定数目光子的集合，确定光源的初始位置和方向及探测器的位置，投放激发光光子，追踪激发光光子在生物组织中传输，计算激发光光子转化为荧光光子比例，沿着激发光路径计算探测器上所接收到的荧光光子权重,保存光子的路径信息；(3)当生物组织内的荧光团吸收系数变...

【技术特征摘要】
2014.11.14 CN 20141064783751.一种基于体素的微扰荧光蒙特卡罗模拟方法，其特征在于包括以下步骤：(1)确定目标生物组织，对目标生物组织进行三维剖分，构建一个三维体素模型，建立一个三维数字矩阵，三维数字矩阵中每个元素与三维体素模型中体素一一对应，每个元素的数值标识一种生物组织，设置生物组织的光学特性参数：吸收系数、散射系数、荧光团吸收系数、折射率和各向异性因子；(2)将入射光源表征为设定数目光子的集合，确定光源的初始位置和方向及探测器的位置，投放激发光光子，追踪激发光光子在生物组织中传输，计算激发光光子转化为荧光光子比例，沿着激发光路径计算探测器上所接收到的荧光光子权重,保存激发光子的路径信息；所述激发光光子转化为荧光光子过程中，假定荧光是各向异性散射，荧光散射方向与激发光散射方向一致；假定所述激发光光子和荧光光子在生物组织中随机行走的路径相同；假定激发光在荧光区域连续激发出荧光；步骤(2)具体按以下步骤进行：(2.1)追踪激发光光子在生物组织中的传输，沿着激发光路径，利用公式计算激发光光子权重w(rs,r)；式中：rs为光源位置、r为荧光激发位置；为激发光的吸收系数、μaf为荧光团吸收系数；rj(j＝1,……pi)为激发光光子从rs到r的第j-1次和第j次散射事件间的路径长度；w0为激发光光子的初始权重、l(rj)为在rj微区所经历的路径长度；(2.2)如果激发光子在生物组织中的荧光区域被吸收并激发产生荧光，依据荧光团系数和量子效率来计算激发光转化为荧光的概率p(r)；p(r)＝η(1-exp(-μafl(r)))式中：η为量子效率；(2.3)沿着激发光路径，利用公式计算荧光光子权重w(r,rd)；

【专利技术属性】
技术研发人员：骆清铭，邓勇，罗召洋，江旭，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：湖北;42