本发明专利技术公开了一种基于稳态辐射率间质测量装置的多波长光学参数反构方法，包括：设定光源‑探测器的距离，以1°为间隔，测量水平面360°范围内辐射率的分布；选取波长，并从测量信号中，提取出该波长下的辐射率分布，利用本发明专利技术所选择的辐射传输方程在P9近似下的解析解作为光子传输正向模型拟合实验所测量的辐射率，最终重建该波长下待测样品的光学参数。本发明专利技术仅需在单个光源‑探测器距离下进行稳态辐射率间质测量，有望实现微创检测；本发明专利技术能对低约化反照率样本(如前列腺增生组织)的光学参数进行准确反构；本发明专利技术所提供方法应用前景广阔，除生物体组织方面的应用，可用于牛奶、液态仿体等光学参数的测量。
【技术实现步骤摘要】
本专利技术属于组织光学研究中的光学参数检测领域，具体涉及一种基于稳态辐射率间质测量装置的多波长组织光学参数反构方法。
技术介绍
前列腺癌是男性常见的恶性肿瘤之一，由于生活方式的改变及人口老龄化，其发病率及死亡率显著上升，已成为影响中国男性健康的重要疾病之一。近年来，光动力疗法(PhotodynamicTherapy，PDT)成为前列腺治疗的新手段[1]。在PDT治疗中，PDT的光剂量是制定治疗方案及在线优化的关键因素。光剂量是指通过测量组织吸收系数、散射系数的分布决定光通量密度的分布，由于从光学参数计算光通量密度分布的方法已较成熟，因此光剂量的关键就是获得组织的光学参数的分布。此外，由于前列腺是内部器官，因此本专利技术考虑的是间质光学参数测量。从光源方面着眼，间质光学参数的测量可分为使用连续光照射的稳态方法、使用光强度调制的频域法及使用脉冲光源的时间分辨方法，其中，稳态测量系统简单、快捷、廉价，有望应用与临床医学。从测量物理量着眼，间质光学参数测量方法可分为基于光子密度的测量方法及基于辐射率的测量方法。常用的基于光子密度的测量方法需要测量多个光源-探测器距离下的光子密度分布，将会加剧内部器官的损伤[1]。与此相对应的，基于辐射率的间质光学测量参数方法[2]仅需要在固定的光源-探测器距离下在测量多个角度下的辐射率分布，有望实现无创或微创测量。此外，与光子密度测量方式相比，辐射率含有的角度信息对光学参数的变化更加敏感，且具有较强的抗噪能力[3]。因此，本专利技术选择基于稳态辐射率的间质光学参数测量。光学参数反构的通常思路为，依据辐射率和生物组织光学参数之间的正向模型，将不同角度下测量的辐射率代入反构算法得到生物组织光学参数的估计值，因此准确的光子传输模型显得尤为重要。在辐射率测量基本框架下，常用光子传输正向模型包括，Monte-Carlo模拟模型、辐射传输方程(radiativetransferequation，RTE)模型、扩散近似模型、PN近似模型。其中，Monte-Carlo模型能够准确地模拟光在组织中的传输过程，但其由于其耗时过长，不适用于临床中快速检测的需求；RTE被认为是描述光子在生物组织中传输的准确模型，但其解析解很难直接得到；扩散近似模型最为简单，但在间质辐射率测量方式中，当光源-探测器距离较小，组织(如前列腺增生组织)的约化返照率较低时，不满足该模型的要求[4]；PN近似模型为RTE在N阶球谐函数展开下的解析解[5]，其解具体为无限媒质中的各项同性点光源下的辐射率。经计算知，PN近似解的误差随阶数N的增大而减小，但当N大于9时，误差几乎不随N的增大而减小，且误差很小。与此同时，PN近似解的计算时间随阶数N的增大而增加。综合考虑，本专利技术选择辐射传输方程的P9近似下的解析解作为光子传输正向模型。[参考文献][1]T.C.Zhu,J.C.Finlay,S.M.Hahn,\Determinationofthedistributionoflight,opticalproperties,drugconcentration,andtissueoxygenationin-vivoinhumanprostateduringmotexafinlutetium-mediatedphotodynamictherapy,\JournalofPhotochemistryandPhotobiologyB:Biology,79(3):231-241(2005)。[2]S.Grabtchak,L.G.Montgomery,W.M.Whelan,\Feasibilityofinterstitialnear-infraredradiancespectroscopyplatformforexvivocanineprostatestudies:opticalpropertiesextraction,hemoglobinandwaterconcentration,andgoldnanoparticlesdetection,\JournalofBiomedicalOptics,19(5):057003(2014)。[3]O.Barajas,A.M.Ballangrud,G.G.Milleretal,\MonteCarlomodelingofangularradianceintissuephantomsandhumanprostate:PDTlightdosimetry,\PhysicsinMedicineandBiology,42(9):1675–1687(1997)。[4]E.L.HULL,T.H.FOSTER,\Steady-statereflectancespectroscopyintheP3approximation,\JournaloftheOpticalSocietyofAmericaA,18(3):584-599(2001)。[5]A.Liemert,A.Kienle,\AnalyticalGreen’sfunctionoftheradiativetransferradiancefortheinfinitemedium,\PhysicalReview,83(3):036605(2011)。
技术实现思路
针对多波长下组织体光学参数间质测量的需求，本专利技术提出一种基于稳态辐射率间质测量装置的多波长组织光学参数反构方法。为了解决上述技术问题，本专利技术提出的一种基于稳态辐射率间质测量装置的多波长光学参数反构方法，其中，所采用的稳态辐射率间质测量装置的结构是，包括光源模块、探测模块、位移模块及控制模块；所述光源模块包括连接至光源光纤一端的卤钨灯；所述光源光纤的另一端设有各向同性散射头，用于将所述卤钨灯发出的光转为点光源；所述探测模块包括侧照光纤和光谱仪，所述侧照光纤的一端耦合至所述光谱仪，所述侧照光纤的另一端为与光纤轴线呈45°角的探测斜面，用于探测与光纤轴线垂直的入射的光信号；所述位移模块包括手动精密位移台、电动平移台和电动旋转台，所述电动平移台的上方固定有一金属立柱，所述手动精密位移台固定在所述金属立柱的上方，所述光源光纤通过一金属固定管与所述手动精密位移台固定；所述电动旋转台固定在一金属支架上，所述侧照光纤通过一金属固定管垂直的固定在所述电动旋转台上；所述手动精密位移台用于调整光源光纤的高度，使各向同性散射头与探测斜面保持在同一水平高度；所述电动平移台用于在水平方向移动光源光纤；所述电动旋转台用于旋转侧照光纤，确定探测斜面的朝向以改变光子接收方向；所述控制模块包括步进电机控制箱和上位计算机；所述上位计算机一方面发出指令至步进电机控制箱，从而控制电动平移台及电动旋转台的运动；另一方面，用于控制光谱仪探测光信号；该多波长光学参数反构方法包括以下步骤：步骤一、利用手动精密位移台、电动平移台，调整光源光纤的空间位置，从而确定各向同性散射头与探测斜面之间的距离，即光源-探测器之间的距离r；通过电动旋转台调整探测斜面的初始探测方向；设定电动旋转台的单位旋转角度为1°；步骤二、所述卤钨灯发出白光，经过光源光纤传输至各向同性散射头，光经过各向同性散射后，以点光源的形式传输至置于样品池中的待测样品中；光经过待测样品的吸收散射后，经过探测斜面传本文档来自技高网...
【技术保护点】
一种基于稳态辐射率间质测量装置的多波长光学参数反构方法，其中，所采用的稳态辐射率间质测量装置的结构是，包括光源模块、探测模块、位移模块及控制模块；所述光源模块包括连接至光源光纤(2)一端的卤钨灯(1)；所述光源光纤(2)的另一端设有各向同性散射头(3)，用于将所述卤钨灯(1)发出的光转为点光源；所述探测模块包括侧照光纤(5)和光谱仪(6)，所述侧照光纤(5)的一端耦合至所述光谱仪(6)，所述侧照光纤(5)的另一端为与光纤轴线呈45°角的探测斜面(4)，用于探测与光纤轴线垂直的入射的光信号；所述位移模块包括手动精密位移台(7)、电动平移台(8)和电动旋转台(9)，所述电动平移台(8)的上方固定有一金属立柱，所述手动精密位移台(7)固定在所述金属立柱的上方，所述光源光纤(2)通过一金属固定管与所述手动精密位移台(7)固定；所述电动旋转台(9)固定在一金属支架上，所述侧照光纤(5)通过一金属固定管垂直的固定在所述电动旋转台(9)上；所述手动精密位移台(7)用于调整光源光纤(2)的高度，使各向同性散射头(3)与探测斜面(4)保持在同一水平高度；所述电动平移台8用于在水平方向移动光源光纤(2)；所述电动旋转台(9)用于旋转侧照光纤(5)，探测斜面(4)的朝向以改变光子接收方向；所述控制模块包括步进电机控制箱(10)和上位计算机(11)；所述上位计算机|(11)一方面发出指令至步进电机控制箱(10)，从而控制电动平移台(8)及电动旋转台(9)的运动；另一方面，用于控制光谱仪(6)探测光信号；其特征在于，该多波长光学参数反构方法包括以下步骤：步骤一、利用手动精密位移台(7)、电动平移台(8)，调整光源光纤(20的空间位置，从而确定各向同性散射头(3)与探测斜面(4)之间的距离，即光源‑探测器之间的距离r；通过电动旋转台(9)调整探测斜面(4)的初始探测方向；设定电动旋转台(9)的单位旋转角度为1°；步骤二、所述卤钨灯(1)发出白光，经过光源光纤(2)传输至各向同性散射头(3)，光经过各向同性散射后，以点光源的形式传输至置于样品池(12)中的待测样品中；光经过待测样品的吸收散射后，经过探测斜面(4)传入侧照光纤(5)，被光谱仪(6)接收；光谱仪(6)测量到的信号传入上位计算机(11)，并存储；步骤三、以步骤一设定的单位旋转角度旋转电动旋转台9，重复步骤二，直至实现在水平面内360°范围内辐射率的测量；步骤四、波长λ下组织光学参数的反构，步骤如下：首先，选取待反构光学参数所在波长为λ，提取出波长λ下的辐射率M＝[m0,m1,...,mk,...,m359]，其中，mk为第k个角度下的辐射率测量值；然后，设定最优化目标函数为：f(μ~a,μ~s)=Σk=0359(mkm0-lk(μ~a,μ~s)l0(μ~a,μ~s))2---(1)]]>式(1)中；分别为样品池中待测样品的吸收系数、散射系数的估计值；是光源‑探测器距离为r、探测角度为k/180×π处P9近似下辐射传输方程的解，即P9近似下的辐射率，其表达式为：lk(μ~a,μ~s)=Σn=0N2n+14πφn(μ~a,μ~s,r)Pn(cos(k·π180))---(2)]]>式(2)中：n＝0,1,...,N，N＝9为近似阶数，k＝0,1,...,359为第k个探测角度；Pn为n阶勒让德多项式；φn(μa,μs,r)为第n阶光子密度计算值，其表达式为：φn(μa,μs,r)=14πΣi=1N+12υi(μa,μs)Gi(ωi(μa,μs))ln(υi(μa,μs))κn(υi(μa,μs)r)---(3)]]>其中，κn(x)的表达式如下：κ0(x)=exp(x)xκ1(x)=(1-1x)exp(x)xκn(x)=κn-2(x)-2n-1xκn-1(x)---(4)]]>定义递推式：Bn+1(ω(μa,μs))＝(2n+1)σnBn(ω(μa,μs))+ω(μa,μs)n2Bn‑1(ω(μa,μs))，其中，σn≡μa+μs(1‑gn)，g为组织体的各向异性因子，g取值为0.9；当，B1(ω(μa,μs))＝1，B2(ω(μa,μs))＝3σ1，时BN+1(ω(μa,μs))=R(ω(μa,μs))=Σl=0N-12alωl(μa,μs)---(5)]]>当，B0(ω(μa,μs))＝1，B1...
【技术特征摘要】
1.一种基于稳态辐射率间质测量装置的多波长光学参数反构方法，其中，所采用的稳态辐射率间质测量装置的结构是，包括光源模块、探测模块、位移模块及控制模块；所述光源模块包括连接至光源光纤(2)一端的卤钨灯(1)；所述光源光纤(2)的另一端设有各向同性散射头(3)，用于将所述卤钨灯(1)发出的光转为点光源；所述探测模块包括侧照光纤(5)和光谱仪(6)，所述侧照光纤(5)的一端耦合至所述光谱仪(6)，所述侧照光纤(5)的另一端为与光纤轴线呈45°角的探测斜面(4)，用于探测与光纤轴线垂直的入射的光信号；所述位移模块包括手动精密位移台(7)、电动平移台(8)和电动旋转台(9)，所述电动平移台(8)的上方固定有一金属立柱，所述手动精密位移台(7)固定在所述金属立柱的上方，所述光源光纤(2)通过一金属固定管与所述手动精密位移台(7)固定；所述电动旋转台(9)固定在一金属支架上，所述侧照光纤(5)通过一金属固定管垂直的固定在所述电动旋转台(9)上；所述手动精密位移台(7)用于调整光源光纤(2)的高度，使各向同性散射头(3)与探测斜面(4)保持在同一水平高度；所述电动平移台8用于在水平方向移动光源光纤(2)；所述电动旋转台(9)用于旋转侧照光纤(5)，探测斜面(4)的朝向以改变光子接收方向；所述控制模块包括步进电机控制箱(10)和上位计算机(11)；所述上位计算机|(11)一方面发出指令至步进电机控制箱(10)，从而控制电动平移台(8)及电动旋转台(9)的运动；另一方面，用于控制光谱仪(6)探测光信号；其特征在于，该多波长光学参数反构方法包括以下步骤：步骤一、利用手动精密位移台(7)、电动平移台(8)，调整光源光纤(20的空间位置，从而确定各向同性散射头(3)与探测斜面(4)之间的距离，即光源-探测器之间的距离r；通过电动旋转台(9)调整探测斜面(4)的初始探测方向；设定电动旋转台(9)的单位旋转角度为1°；步骤二、所述卤钨灯(1)发出白光，经过光源光纤(2)传输至各向同性散射头(3)，光经过各向同性散射后，以点光源的形式传输至置于样品池(12)中的待测样品中；光经过待测样品的吸收散射后，经过探测斜面(4)传入侧照光纤(5)，被光谱仪(6)接收；光谱仪(6)测量到的信号传入上位计算机(11)，并存储；步骤三、以步骤一设定的单位旋转角度旋转电动旋转台9，重复步骤二，直至实现在水平面内360°范围内辐射率的测量；步骤四、波长λ下组织光学参数的反构，步骤如下：首先，选取待反构光学参数所在波长为λ，提取出波长λ下的辐射率M＝[m0,m1,...,mk,...,m359]，其中，mk为第k个角度下的辐射率测量值；然后，设定最优化目标函数为：f(μ~a,μ~s)=Σk=0359(mkm0-lk(μ~a,μ~s)l0(μ~a,μ~s))2---(1)]]>式(1)中；分别为样品池中待测样品的吸收系数、散射系数的估计值；是光源-探测器距离为r、探测角度为k/180×π处P9近似下辐射传输方程的解，即P9近似下的辐射率，其表达式为：lk(μ...
【专利技术属性】
技术研发人员：赵会娟，刘玲玲，蒋景英，张丽敏，高峰，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。