[发明专利]全光学生物发光断层成像方法

摘要

本发明公开了一种全光学生物发光断层成像方法，解决了现有技术中需要融合结构成像技术，如计算机断层成像或磁共振成像，才能实现三维重建的问题。该方法采用白光信号图像获取生物体表面三维轮廓，结合生物体表面能量分布的定量重建，用于体内生物发光光源的反演，以实现低成本的全光学生物发光断层成像过程。具体实现步骤包括：数据采集与预处理、生物体表面三维轮廓重建、生物体表面能量分布重建、表面能量的定量校准和体内光源反演及三维显示。本发明的全光学生物发光断层成像方法具有灵敏度高、成本低的特点，同时避免了复杂繁琐、耗时耗力的分割和配准环节，可用于生物发光断层成像领域。

主权项

1.全光学生物发光断层成像方法，其实现的具体步骤是：步骤一，数据采集与预处理，采集多角度生物发光信号和白光信号图像，并依次进行背景噪声去除、感兴趣区域提取和坏点补偿等数据预处理；步骤二，生物体表面三维轮廓重建，(2.1)白光信号图像的阈值分割，应用下述公式对步骤一中采集的多角度白光信号图像进行阈值分割，将其转换为二值图像：I‾(x,y,θi)=1∀I(x,y,θi)≥Ith(i)0∀I(x,y,θi)<Ith(i)]]>式中，x是每幅白光信号图像中水平方向上的像素值；      y是每幅白光信号图像中竖直方向上的像素值；      θ_i是第i幅白光信号图像的拍摄视角；      I是经步骤一获取的白光信号图像；是经阈值分割之后的白光信号图像；是第i幅白光信号图像阈值分割时采用的阈值，其值可以根据白光信号图像的信噪比确定；(2.2)二值图像的图像边界提取；(2.3)确定生物体表面三维轮廓的笛卡尔坐标系；(2.4)计算二值图像边界点的坐标值，在步骤(2.3)确定的笛卡尔坐标系的基础上，应用下述公式计算步骤(2.2)中获得的图像边界点的坐标：X=mlcosθY=mlsinθZ=nlSIGN(flag)]]>式中，(X，Y，Z)是图像边界点的笛卡尔坐标；      m是当前像素点与Z轴之间的有效像素数量；      l是单个像素代表的物理长度；      θ是从0度图像算起的总的旋转角度；n是当前像素点与XOY平面之间的有效像素数量；SIGN()是符号函数；flag是标志位，如果当前点在XOY平面的上部，其值取1，否则其值取0；(2.5)建立生物体表面三维轮廓的三维体数据结构，应用步骤(2.4)中的公式，遍历步骤(2.2)中获得的所有图像边界，得到以点的形式构成的生物体表面的三维体数据结构；(2.6)建立生物体表面三维轮廓的三角面片数据结构；步骤三，生物体表面能量分布重建，利用步骤二重建的生物体表面三维轮廓和步骤一获得的生物发光信号图像，应用下式计算生物体表面能量分布：P=1π∫S∫ΩE(rd)ξ(rd,r)1|r-rvd|cosθscosθddΩdS]]>式中，P是生物体表面能量分布；      r是生物体表面三维轮廓S上的任意一点；      r_d是CCD相机上的任意一点；      r_vd是点r_d在虚拟探测平面上的成像点，虚拟探测平面定义为CCD相机的焦平面；      E(r_d)是CCD相机点r_d处的强度值；      ξ(r_d，r)是表面点r与CCD相机上点r_d之间的可视因子，其值由表面点r处出射的光线的方向、CCD相机中的薄透镜位置和直径尺寸确定；      cosθ_s是表面点r的法向量与表面点r到r_vd点之间连线的夹角余弦；      cosθ_d是r_vd点的法向量与r_vd点到表面点r之间连线的夹角余弦；      dΩ是CCD相机上点r_d处的微元面积；      dS是表面点r处的微元面积；步骤四，表面能量的定量校准，(4.1)利用积分球均匀光源系统，校准平面光源的出射流量密度与CCD相机所采集图像的灰度值之间的关系；(4.2)建立不同距离、不同视场情况下，生物体表面出射流量密度与CCD相机所采集图像的灰度值之间的定量关系：J=(G+r1te+r2)×r3(R-d)r4R]]>式中，J是生物体表面点的出射流量密度，单位为纳瓦每平方毫米；      G是CCD相机所采集图像的像素灰度值；      t_e是CCD相机采集生物发光信号时的曝光时间，单位为秒；      R是生物体表面点到有效透镜边缘的距离，单位为毫米；      d是生物体表面点到有效透镜中心的距离，单位为毫米；      r₁是18～23范围内的常数；r₂是7～12范围内的常数；r₃是57～62范围内的常数；r₄是3～8范围内的常数；(4.3)计算生物体表面的绝对能量分布；步骤五，体内光源反演及三维显示。