一种光学投影断层成像中的吸收和散射系数重建方法
【专利摘要】本发明公开了一种光学投影断层成像中的吸收和散射系数重建方法，利用级数展开理论分解传输光强，结合OPT成像模型，利用测量信号中的弹道传播信号和散射信号的特性，采用统计迭代算法从散射分量中重建生物组织的散射系数，再利用衰减系数、散射系数与吸收系数之间的关系计算获得生物组织的吸收系数，最终实现生物组织吸收系数和散射系数空间分布的同时重建，应用于弱散射介质的投影断层重建。本发明实现了生物组织吸收系数和散射系数的三维重建，既有效解决了OPT成像存在的散射问题，又提高了OPT成像质量，同时更丰富了OPT技术提供的信息量。
【专利说明】
一种光学投影断层成像中的吸收和散射系数重建方法
技术领域
[0001] 本发明属于医学影像处理技术领域，尤其涉及一种光学投影断层成像中的吸收和 散射系数重建方法。
【背景技术】
[0002] 光学投影断层成像(Optical projection tomography,OPT)是一种新型的高分辨 率三维分子影像成像技术，其成像原理和X射线计算机断层成像的原理类似。OPT可获得扫 描样品结构像，又能利用荧光染料或荧光蛋白进行分子特异性标记，实现分子特性成像，并 且设备成本低，使用方便。
[0003] 0PT成像由于生物样品未经透明化处理，因此存在光子散射问题。常规的0PT成像 算法只考虑生物组织对光子的吸收特性，也只能重建光子的吸收系数不能重建光子的散射 系数。当采用0PT技术进行成像时，散射的影响不可忽略，并且会和光子的吸收特性混合在 一起，从而导致传统的0PT成像空间分辨率的降低以及重建结果的不准确。
【发明内容】
[0004] 本发明的目的在于提供一种应用于弱散射介质的光学投影断层成像中的吸收和 散射系数重建方法，旨在解决常规的0PT成像算法只考虑生物组织对光子的吸收特性，无法 同时重建光子的吸收系数和光子散射系数的问题。
[0005] 本发明是这样实现的，一种应用于弱散射介质的光学投影断层成像中的吸收和散 射系数重建方法利用级数展开理论分解传输光，将传输光分解为弹道传输分量、一次散射 传输分量和多次散射传输分量，结合本发明中弱散射介质的应用背景，将多次散射传输分 量忽略不计;结合0PT成像模型，分别测量弹道传输分量和一次散射传输分量;采用重建算 法从弹道传输分量中重建样品的衰减系数，从一次散射分量中重建样品的散射系数;再利 用衰减系数、散射系数与吸收系数之间的关系计算获得样品的吸收系数，最终实现样品吸 收系数和散射系数空间分布的同时重建。
[0006] 进一步，所述光学投影断层成像中的吸收和散射系数重建方法包括以下步骤：
[0007] 步骤一，光子弱散射生物组织传输建模；
[0008] 步骤二，光源和CCD相机相对放置时采集投影数据；
[0009] 步骤三，重建生物组织衰减系数；
[0010] 步骤四，光源和CCD相机偏移一定角度放置时采集投影数据；
[0011] 步骤五，重建生物组织散射系数；
[0012] 步骤六，计算样品的吸收系数。
[0013] 进一步，所述光子弱散射生物组织传输建模包括：
[0014] 表示上入射的光子到位置？且方向为;;时所传播的距离，则：
[0015] L| r,v) = ^(!gi；; =cxp(-￡' ' uXsv'yis)-；
[0016] 其中zP. b为？处的辐射度，表示单位立体角内、单位时间内、通过垂直于单位方向 矢量;;的单位面积上的平均功率通量密度，量纲是W/(m2. Sr) ;K〇为引入的光子弹道传播算 子，M F卜//a 0+A0表示生物组织在；:处的总的衰减系数，凡闪表示生物组织的吸收 系数，表示散射系数；
[0017] 再定义： pvify) f t /?
[0018] AX(r,v) =￡ cxp(-j -/^J^ p(\\\r)L(J: -t\\\r)diYdt ?
[0019] 其中K为引入的光子散射传播算子，为归一化的散射相位函数，表示光子从 方向$散射到〔方向的概率，满足]'/)&1/)^2' = 1，(^/表示单位方向矢量1；，对应的立体 角微元;定义m〇 = K()gin，:《" = ，贝lj有nn+i = Kmn(n彡0)，从而在f处的总的福射率为：
[0020] UK) = fjn,,=乞(K Y K hK 1:, =、! - K、_h, ' n-0. n--.{)
[0021] 其中mn表示经n次散射到达；:的辐射度分量；当光子在扩散区域传播时，K的谱半径 P(K)值接近1，当光子在弱散射区域传播时，P(K)》1，在此种情况下，当n-m时，快 速收敛；
[0022] 然后，在输出边界r+上描述探测器接收到的数据总量gQUt，即g ， 从而：
[0023] =K{)gm +^0gm +K2(l~KylK0gm\^ = A()gm + Algm+A2gin =Agm；
[0024] 其中A是描述光子传输的矩阵，Ao、Ai和A2分别描述弹道传输、一次散射传输和多次 散射传输部分，定义8〇 = 4呢^，81=4培^分别表不测量数值中的弹道传输分量和一次散射分 量，则知：
[0025] n," exp(-L，?",(、〒)办"G「1 ;
[0026] g, -1(11 "exiX-^JLI,(.s-v^^exiX-^' '(sv)ds) /js(/v)jt /iy\\^)g..(KldLre；
[0027]入射光方向为r，经一次散射后其方向为$，则上式中关于的积分仅在一 特定角度上有值，取I '系数k的取值由相位函数确定，同 .4. 7T 时定义似丨⑴=exp(-//,(.vi7)A) 'C〇2(t) = cxp(-￡ 办）分别表示光 子发生散射后和散射前的衰减量，则有：
[0028] g! = \ m^(j)- m2{t)^ids{tv)dt., r eT+〇 v 0
[0029]进一步，所述获取投影数据包括：
[0030] 首先，将光源和探测器相对放置，则光源以平行光对样品进行照射并且远心透镜 耦合的CCD相机接收平行入射的光线；
[0031] 其次，照射源对固定在电控旋转台上的样品进行水平投影断层扫描，照射源采用 激光器，并使用远心透镜将光线扩束为平行光照射样品；
[0032] 最后，使用CCD相机作为探测器采集样品的投影数据，采集的方法样品每旋转一个 角度采集一个投影数据，将每次采集的投影数据记录在计算机中 ±
[0033] 进一步，所述计算重建衰减系数包括：
[0034]在OPT成像中采用空间均匀分布的平行光对样品gin进行照射，通过CCD相机采集无 样品遮挡的照射光测得入射光强度;对" 厂广左右两 边同除以gin并取负对数，则：
[0036]采集到360度的测量数据Go后，采用精确高效的滤波反投影重建算法实现逆Radon 变换即计算出衰减系数，即yt = FBP(Go)。
[0037] 进一步，所述偏移一定角度采集投影数据包括：
[0038] 将OPT系统的探测器件偏移一定角度0放置，照射源对固定在电控旋转台上的样品 进行水平投影断层扫描，照射源采用激光器，并使用远心透镜将光线扩束为平行光照射样 品；
[0039] 使用探测器采集样品的投影数据，样品每旋转一个角度采集一个投影数据，将每 次采集的投影数据记录在计算机中。
[0040] 进一步，所述重建散射系数由公式良=| ( ' ％ (/) n(/) _(W)t//,g^ 含了OPT成像中散射的影响，当从某一确定的角度采集数据81时，交和i?的散射夹角确定，系 数k为一个确定的常数;两边同除以kgin，则有：
[0042] 由上式可知G , (0为散射系数延f方向的加权Radon变换，所加权值《 i(t)和《 2 (t)均是与衰减系数有关的函数，将<^(0离散化并以矩阵的形式表示如下：
[0043] ffys = Gi；
[0044]其中W表示离散化后的权值矩阵，yjPGi分别表示散射系数矢量和不同角度测量得 到的归一化测量矢量，利用带罚函数的加权最小二乘准则建立如下的目标函数：
[0046]其中表达式的第一项是似然函数的近似表达形式，第二项R(ys)为正则项，通常根 据图像的先验信息构造而成，0为正则化因子，矩阵C为协方差矩阵。以m表示CCD探测器检 测到的散射光子数，对应的协方差矩阵表不为： % 0 00 0、 0 ?2 0 0 0
[0047] €=0 0.?，.- 0 0 -- ?
[0048] 利用最优化方法对0(ys)的目标函数求解，即可求出散射系数：
[0049] ys = argminC> (ys) 〇
[0050] 进一步，所述计算吸收系数利用步骤三和步骤五的计算结果，计算样品的吸收系 数，利用关系式yt=y a+ys计算样品的吸收系数ya。
[0051] 本发明在透射式OPT成像方面，提出了同时重建吸收系数和散射系数的有效方法。 针对OPT成像存在弱散射的问题，结合OPT成像数据采集的自身特点，构建相应的数学模型， 并通过额外测量一组与入射光保持一定倾斜角度的数据的方法，对弹道传输分量和一次散 射传输分量进行分离，进而实现样品吸收系数和散射系数的三维重建，从而既可以有效解 决OPT成像存在的散射问题;提高了-OPT成像质量，同时更丰富了OPT技术提供的信息量，使 得0PT技术可以从吸收系数和散射系数两个角度描述生物的组织结构特性。
【附图说明】
[0052] 图1是本发明实施例提供的光学投影断层成像中的吸收和散射系数重建方法流程 图。
[0053]图2是本发明实施例提供的实施例的流程图。
[0054] 图3是本发明实施例提供的光学投影断层成像系统示意图。
[0055] 图4是本发明实施例提供的光线照射及射出样品的示意图。
[0056] 图5是本发明实施例提供的光源和CCD相机相对放置时信号采集示意图。
[0057] 图6是本发明实施例提供的光源和CCD相机偏移一定角度放置时信号采集示意图。
[0058] 图7是本发明实施中衰减系数与散射系数及吸收系数之间的关系示意图。
【具体实施方式】
[0059] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明 进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于 限定本发明。
[0060] 下面结合附图对本发明的应用原理作详细的描述。
[0061] 如图1所示，本发明实施例的光学投影断层成像中的吸收和散射系数重建方法包 括以下步骤：
[0062] S101:光子弱散射生物组织传输建模；
[0063] S102:光源和C⑶相机相对放置时采集投影数据；
[0064] S103:重建生物组织衰减系数；
[0065] S104:光源和(XD相机偏移一定角度放置时采集投影数据；
[0066] S105:重建生物组织散射系数；
[0067] S106:计算样品的吸收系数。
[0068] 下面结合具体实施例对本发明的应用原理作进一步的描述。
[0069] 如图2所示，本发明实施例的光学投影断层成像中的吸收和散射系数重建方法包 括以下步骤：
[0070] 步骤1，0PT光传输建模：
[0071] 首先，描述弱散射条件下的光子在生物组织内的传播，定义rib表示r-上入射的 光子到位置且方向为；时所传播的距离，如果仅考虑衰减对入射光的影响，则：
[0072] /.(/%v) = A'(lg," -oxp(-￡' 'jn,{sv)ds)-g!i：；
[0073] 其中if S)为f处的辐射度(Radiance)，表示单位立体角内、单位时间内、通过垂直 于单位方向矢量；的单位面积上的平均功率通量密度，它的量纲是W/(m 2. Sr) ;K〇为引入的 光子弹道传播算子，f卜凡D表示生物组织在?处的总的衰减系数，~表示 生物组织的吸收系数，A 表示散射系数。
[0074] 再定义：
[0075] KL(r,v^ =|_ cxp(-| //j.wOt/v)-//J p(v,v)L(F>
[0076] 其中K为引入的光子散射传播算子，/4夂.s' j为归一化的散射相位函数，表;^光子从 方向r?散射到;;方向的概率，满足丨/;^，'|/} /^ = 1，(1〇/表示单位方向矢量$对应的立体角 Am- 微元。定义_ = K〇gin，w,, = ，则有nn+i = Kmn(n彡0)，从而在f处的总的福射率为：
[0077] L(rS-) = X mn = X (^)BK0gk = (I-KY^K^ M = 0 H-{)
[0078] 其中mn表示经n次散射到达的辐射度分量。当光子在扩散区域传播时，K的谱半径P (K)值接近1，当光子在弱散射区域传播时，P(K)>>1，在此种情况下，当n-①时乂 快速收敛。
[0079] 然后，在输出边界r+上描述探测器接收到的数据总量gQUt，即^， 从而：
[0080] gout = K,gin u +KK〇gin u +K20-K)1 Ki)giii = 4)gin + Algiii + Algin = Ag.t ；
[0081 ]其中A是描述光子传输的矩阵，Ao、Ai和A2分别描述弹道传输、一次散射传输和多次 散射传输部分，定义8〇 = 4呢^，81=4培^分别表不测量数值中的弹道传输分量和一次散射分 量，则可知： /?T(r,v) .. 一
[0082] g{) = g.n exp(-j( ", 〇〒)办)，厂 e r ; /?r(F,y) M f
[0083] g\ exp(-|(//.(.svk/v)exp(-|( //.(si^)cis) tus(tv)j iii\\!)gii!cKldi,reP；
[0084]考虑到OPT成像中gin以平行光形式入射，并在一个特定的角度采集数据。假设其入 射光方向为铲，经一次散射后其方向为则上式中关于的积分仅在一特定角度上 有值，取J 系数k的取值可以由相位函数p(f,v')确定，同时定义 (0,(t) = 〇xp(-jn ^i(sv)cls)，似2(f) = exp(-j,:w '//,(.v〇/'s')分别表示光子发生散 射后和散射前的衰减量，则有：
[0085] gj = kgiu f (〇l(t)-(〇2(t)- u^(lv)dt, .? el"1 J 0
[0086] 步骤2,获取投影数据：
[0087] 首先，将光源和探测器相对放置，则光源以平行光对样品进行照射并且远心透镜 耦合的CCD相机只能接收平行入射的光线。
[0088] 其次，照射源对固定在电控旋转台上的样品进行水平投影断层扫描，照射源采用 激光器，并使用远心透镜将光线扩束为平行光照射样品。
[0089]最后，使用探测器采集样品的投影数据。本发明采集的方法是:样品每旋转一个 角度采集一个投影数据，将每次采集的投影数据记录在计算机中。
[0090] 步骤3,计算重建衰减系数：
[0091] 在OPT成像中采用空间均匀分布的平行光对样品gin进行照射，通过CCD相机采集无 样品遮挡的照射光测得入射光强度。对g9 =_ exp.(_- J"。< hF e r +左右两 边同除以gin并取负对数，则：
[0093]采集到360度的测量数据Go后，采用精确高效的滤波反投影（Filtering Backprojection，FBP)重建算法实现逆Radon变换即可计算出衰减系数，即yt = FBP(G〇)。 [0094]步骤4,偏移一定角度采集投影数据：
[0095]首先，将OPT系统的探测器件偏移一定角度0放置，照射源对固定在电控旋转台上 的样品进行水平投影断层扫描，照射源采用激光器，并使用远心透镜将光线扩束为平行光 照射样品。
[0096]然后，使用探测器采集样品的投影数据。样品每旋转一个角度采集一个投影数据， 将每次采集的投影数据记录在计算机中。
[0097]步骤5，重建散射系数
[0098]由公式私_二f ( %⑴?欣，（〇 ?//、（/〇<：//，可知gi包含了OPT成像中散射的 影响，当从某一确定的角度采集数据gi时，交和铲的散射夹角确定，系数k为一个可以确定的 常数。对(7)左右两边同除以kgin，则有：
L0100J 由上式口」知h (v)为散射糸数延v万冋的加权Radori变秧，所加权值《 Kt)和《 2 (t)均是与衰减系数有关的函数，利用步骤三中计算的衰减系数可以确定，从而散射系数可 以类似于吸收系数的重建，通过逆Radon变换计算得到。采用FBP算法实现这一过程，但是考 虑到散射分量相对微弱，信噪比相对比较低，因此本发明中拟采用基于统计迭代的方式实 现。为了描述方便，将q ((0离散化并以矩阵的形式表示如下：
[0101] ffys = Gi；
[0102] 其中W表示离散化后的权值矩阵，yjPGi分别表示散射系数矢量和不同角度测量得 到的归一化测量矢量。CCD探测到的散射光子的统计分布可以采用泊松模型来描述。基于该 统计模型，构造相应的最大似然函数，即可进行图像的三维重建。最大似然图像重建是一个 高维的求逆问题，通常具有不稳定性，通常的解决办法是在似然函数中引入正则项，从而可 以利用带罚函数的加权最小二乘(Penalizedweighted least-squares，PWLS)准则建立如 下的目标函数：
[0104] 其中表达式的第一项是似然函数的近似表达形式，第二项R(ys)为正则项，通常根 据图像的先验信息构造而成，0为正则化因子，矩阵C为协方差矩阵。以m表示CCD探测器检 测到的散射光子数，对应的协方差矩阵可以表不为： 、0 0 0 0'、 0 , 0 0 0 .
[0105] C~ 0 0 ? ? ? 0 0 ? ? ? 0 ..... nM/
[0106] 利用最优化方法对0(ys)的目标函数求解，即可求出散射系数：
[0107] Us = argmin〇 (ys)；
[0108]步骤6，计算吸收系数：
[0109]利用步骤3和步骤5的计算结果，计算样品的吸收系数。利用关系式^ =~+^计算 样品的吸收系数ya。
[0110] 下面结合附图3、附图4、附图5和附图6对本发明的重建结果做详细的描述。
[0111] 附图3为本发明的OPT成像示意图，在该系统中采用半导体激光器作为光源，并通 过远背光照明器将光束扩散为平行光对成像物体进行照射，在接收端采用与远心透镜耦合 的CCD相机进行信号接收;为了避免光线在生物组织体表产生折射效应，通常将样品浸泡在 BABB光学匹配液中，样品匹配液池采用与匹配液折射系数相当的硼硅酸盐玻璃制成。远背 光照明器使得入射光以平行光的形式对样品进行照射，远心透镜的使用保证CCD相机仅接 收平行入射的光线，发挥了类似于SPECT中射线准直器的作用，这既有助于高质量的OPT成 像，同时也可以在后续的数学建模中忽略折射的影响。
[0112]附图4是OPT成像的示意图，其中gin表示平行的入射光，gciut表示出射光线，射入和 射出边界分别用表示。
[0113] 未经运动伪影校正的重建结果。其中，扫描样品是固定在一个毛细玻璃管中的线 虫，旋转角度为360°，采集的投影数据的个数是500个。
[0114] 附图5是表示光源和CCD相机相对放置时信号采集的示意图，则光源以平行光对样 品进行照射并且远心透镜耦合的CCD相机只能接收平行入射的光线，因此在这种情形下， CCD相机可以接收到弹道传播的光子(如图中路径1 2所示），无法探测到一次散射的光子(如 图中路径li所示），对于多次散射的光子，可能会探测到(如图中路径13所示），也可能无法探 测到(如图中路径14所示），将此时探测器探测到数据记为gL。
[0115] 附图6是光源和CCD相机偏移一定角度放置时信号采集示意图。当将OPT系统中的 探测器件偏移一定角度探测时，只有满足特定散射角度的一次散射分量和部分多次散射分 量可以被探测到，将探测器探测到数据记为。
[0116] 附图7是衰减系数与散射系数ys及吸收系数ya之间关系的示意图，反应为公式即 14t 一 Ua+Us 〇
[0117] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精 神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种光学投影断层成像中的吸收和散射系数重建方法，其特征在于，所述光学投影 断层成像中的吸收和散射系数重建方法利用级数展开理论分解传输光，将传输光分解为弹 道传输分量、一次散射传输分量和多次散射传输分量，结合本发明中弱散射介质的应用背 景;结合OPT成像模型，分别测量弹道传输分量和一次散射传输分量;采用重建算法从弹道 传输分量中重建样品的衰减系数，从一次散射分量中重建样品的散射系数;再利用衰减系 数、散射系数与吸收系数之间的关系计算获得样品的吸收系数，最终实现样品吸收系数和 散射系数空间分布的同时重建。2. 如权利要求1所述的光学投影断层成像中的吸收和散射系数重建方法，其特征在于， 所述光学投影断层成像中的吸收和散射系数重建方法包括以下步骤： 步骤一，光子弱散射生物组织传输建模； 步骤二，光源和CCD相机相对放置时采集投影数据； 步骤三，重建生物组织衰减系数； 步骤四，光源和CCD相机偏移一定角度放置时采集投影数据； 步骤五，重建生物组织散射系数； 步骤六，计算样品的吸收系数。3. 如权利要求2所述的光学投影断层成像中的吸收和散射系数重建方法，其特征在于， 所述光子弱散射生物组织传输建模包括： )表示r -上Λ射的I子刹份罟：日玄-为二0计所佑播的昍商则：其中M?,为？处的辐射度，表示单位立体角内、单位时间内、通过垂直于单位方向矢量 ?的单位面积上的平均功率通量密度，量纲是W/(m2. Sr) ;KQ为引入的光子弹道传播算子， 卜)+凡.Ir)表示生物组织在？处的总的衰减系数，/(,D表示生物组织的吸收系 数，//,f 表示散射系数； 再定夕其中K为引入的光子散射传播算子，pi.W' 1为归一化的散射相位函数，表ZK光子从方向r 其中mn表示经η次散射到达?的辐射度分量；当光子在扩散区域传播时，K的谱半径P(K) 值接近1，当光子在弱散射区域传播时，P(K) >> 1，在此种情况下，当n-m时，快速 收敛； 然后，在输出边界Γ+上描述探测器接收到的数据总量8咖，即尽_=1^?|]^，从而：其中A是描述光子传输的矩阵，Ao、Al和Α2分别描述弹道传输、一次散射传输和多次散射 传输部分，定义go = Aogin，gl = Algin分别表亦测量数值中的弹道传输分量和一次散射分量， frlll/cn *4. 如权利要求2所述的光学投影断层成像中的吸收和散射系数重建方法，其特征在于， 所述获取投影数据包括： 首先，将光源和探测器相对放置，则光源以平行光对样品进行照射并且远心透镜耦合 的CCD相机接收平行入射的光线； 其次，照射源对固定在电控旋转台上的样品进行水平投影断层扫描，照射源采用激光 器，并使用远心透镜将光线扩束为平行光照射样品； 最后，使用探测器采集样品的投影数据，采集的方法样品每旋转一个角度采集一个投 影数据，将每次采集的投影数据记录在计算机中。5. 如权利要求2所述的光学投影断层成像中的吸收和散射系数重建方法，其特征在于， 所述计算重建衰减系数包括： 在OPT成像中采用空间均匀分布的平行光对样品gin进行照射，通过CCD相机采集无样品 遮挡的照射光测得入射光强度；Ξ右两边同 除以gin并取负对数，采集到360度的测量数据Go后，采用精确高效的滤波反投影重建算法实现逆Radon变换 即计算出衰减系数，即yt = FBP(G〇)。6. 如权利要求2所述的光学投影断层成像中的吸收和散射系数重建方法，其特征在于， 所述偏移一定角度采集投影数据包括： 将OPT系统的探测器件偏移一定角度Θ放置，照射源对固定在电控旋转台上的样品进行 水平投影断层扫描，照射源采用激光器，并使用远心透镜将光线扩束为平行光照射样品； 使用探测器采集样品的投影数据，样品每旋转一个角度采集一个投影数据，将每次采 集的投影数据记录在计算机中。7. 如权利要求2所沭的光学枵影断层成像中的吸收和散射系教重建方法，其特征在于， 所述重建散射系数由&含了 OPT成像中散 射的影响，当从某一确定的角度采集数据gi时，芬和铲的散射夹角确定，系数k为一个确定的 常数;两边同除以_夂由上式知G, (v~.)为散射糸数延交万冋的加权Radon变秧，所加权值ω : (t)和ω 2(t)均是 与衰减系数有关的函数，将Gi (P)离散化并以矩阵的形式表示如下： Wys=Gi; 其中W表示离散化后的权值矩阵，WdPG1分别表示散射系数矢量和不同角度测量得到的 归一化测量矢量，利用带罚函数的加权最小二乘准则建立如下的目标函数：其中表达式的第一项是似然函数的近似表达形式，第二项R(ys)为正则项，通常根据图 像的先验信息构造而成，β为正则化因子，矩阵C为协方差矩阵；以m表示CCD探测器检测到 的散射光子数，对应的协方差矩阵表不为：利用最优化方法对φ (ys)的目标函数求解，即求出散射系数： ys = arg η?ηΦ (ys) 〇8. 如权利要求2所述的光学投影断层成像中的吸收和散射系数重建方法，其特征在于， 所述计算吸收系数利用步骤三和步骤五的计算结果，计算样品的吸收系数，利用关系式W =μ3+μ5计算样品的吸收系数ya。
【文档编号】A61B5/00GK105894562SQ201610204224
【公开日】2016年8月24日
【申请日】2016年4月1日
【发明人】朱守平, 郎晋伟, 鲍翠平, 谢晖, 曹旭, 张建华, 梁继民
【申请人】西安电子科技大学