CN105786762A - 一种人皮肤光谱的建模方法以及高拟合度的多个皮肤参数的数学建模方法 - Google Patents

Abstract

一种人皮肤光谱的建模方法以及高拟合度的多个皮肤参数的数学建模方法，本发明通过建立两个层级的分析模型，即皮肤光谱分析模型和皮肤参数数学模型，从而在皮肤光谱与光敏感的皮肤生物学参数之间建立起关联，实现由一组皮肤参数模拟出虚拟光谱，然后通过虚拟光谱与实际光谱利用信息处理技术进行优化迭代，能够找到一组最优解，达到指定的拟合度标准，实现皮肤生物学量化分析的目的。

Description

一种人皮肤光谱的建模方法以及高拟合度的多个皮肤参数的数学建模方法

技术领域

本发明涉及计算生物学，涉及建立皮肤光谱的分析模型的方法以及利用该模型构建皮肤生物学参数的建模方法，特别是涉及人皮肤光谱的建模方法以及高拟合度的多个皮肤参数的数学建模方法。

背景技术

光与物质相互作用会引起物质内部原子及分子能级的电子跃迁，使物质对光的吸收、反射、散射等在波长及强度信息上发生变化，光谱仪可用于检测并处理这类变化。与其他分析方法相比，光谱检测具有非破坏、高灵敏、高精准的特点，因而在各类材料的检测和鉴定方面得到广泛的应用。

生物学与光学的研究表明，人类皮肤是由多种生物学成分组成，其中的一些成分对光谱作用敏感，各成分具有特定的可见光光学性质，并且，这些皮肤成分含量已由生物学方法进行了测定。因此，形成一种籍于皮肤的可见光谱的数据处理方法实现皮肤的生物学成分量化分析是可行的。但在建模方面，要求可见光光学作用分析完整、皮肤参数尽量齐全，才能体现皮肤参数的光学性质细节，达到可见光全波段所需的量化分析精度要求。

目前，在皮肤参数量化分析领域主要存在以下手段，1)基于图像的皮肤分析。由于是平面的图像分析，其所能测得参数的与光谱的立体数据相比准确性有较大差距。所期望的其他的皮肤成分参数和皮肤结构参数无法获得。2)基于生物阻抗的皮肤水份检测。此种检测成分参数单一且精度不够。3)基于超声波的皮肤超声影像诊断系统。此种方式主要是对皮肤的结构组织进行定性观察，无法进行皮肤参数量化分析。4)基于X射线的三维断层成像的皮肤CT影像分析系统。此种方式主要是对皮肤的结构组织进行定性观察，无法进行皮肤参数量化分析。5)皮肤光谱检测系统。属于研究项目。此种方式是利用光谱进行某波段的直观比较观察，没有形成信息处理模型。

同时，建立皮肤渲染模型能够直观反映皮肤结构和组分。该类模型主要用于皮肤渲染，强调视觉效果，精度不是关键指标，因而将其用于皮肤参数量化分析存在问题有：①皮肤参数少，仅有4个成分参数变量，不能完整反映皮肤实际，用于量化分析的结果精度不够。②仅反映可见光的400-600nm波长区间，600nm以上波段无法拟合。(Kubelka-Munk模型由于未考虑水平方向散射，其结果会高于真实结果；Multipole模型由于其假设的离散点光源位置问题，在某些过薄皮肤层情况下无法模拟，该情况经常发生于590-700nm波长区间；Monte-Carlo模型结果具有不确定性，难以优化，且所需计算时间很长)③仅涉及皮肤的吸收相关的皮肤成分参数，散射相关的皮肤成分参数。④不包含皮肤结构参数。

因此，为了实现精确的皮肤生物学参数的量化目的，急需设计一种基于皮肤光谱进行皮肤参数量化分析的方案，可见光波段分析完整、皮肤生物学特征参数尽可能充足地被涉及，从而达到对不同人、不同部位差异尽可能多的分析维度和精度，并形成大数据处理的基础。

发明内容

本发明的目的是提供一种人皮肤光谱的建模方法以及高拟合度的多个皮肤参数的数学建模方法。填补了利用光谱模型进行皮肤生物学量化分析的空白，分析出可以描述皮肤光谱特征的一组皮肤参数，其虚拟出的皮肤光谱与实际皮肤光谱的拟合度非常高，提高分析的精确性。

为实现上述发明目的，本发明提供的技术方案是：

一种人皮肤光谱的建模方法，所述方法包括以下步骤：

步骤一根据可见光照射皮肤时吸收、反射、散射、透射的特征，将皮肤从上至下抽象为皮肤表面层、皮肤表皮层、皮肤真皮层、皮下组织层四层；

步骤二根据皮肤表皮层的光吸收、散射系数以及皮肤表皮层厚度变量参数，建立在皮肤表皮层中光反射和光透射方程，计算在皮肤表皮层的反射率和透射率；

步骤三根据皮肤真皮层的光吸收、散射系数以及皮肤真皮层厚度变量参数，建立在皮肤真皮层中光反射方程，计算在皮肤真皮层的反射率；

步骤四计算光在皮肤表皮层和皮肤真皮层的两层结构中的反射率；

步骤五计算光通过四层皮肤的总体结构中的反射率。

进一步地，在步骤一和步骤二之间还包括根据辐射传导方程和在皮肤粗糙表面层的反射系数，建立在皮肤粗糙表面层中光反射和光透射的方程，计算在皮肤表面层的反射率的步骤。

进一步地，所述在皮肤表面层中光反射和光透射的方程为[R_surf,T_surf]＝LSI(r_surf)，其中，R_surf为表面的光反射率，T_surf为表面的光透射率，LSI代表了用于计算皮肤反射率以及透射率的一系列计算公式，r_surf为皮肤表面的反射系数。

进一步地，在步骤二中，在所述皮肤表皮层中光反射和光透射由光从上方空气中进入和光从下方真皮层进入两种途径所决定，所述在皮肤表皮层中光反射和光透射的方程为和其中，为从空气进入表皮层的光反射率，为从空气进入表皮层的光透射率，为从真皮层进入表皮层的光反射率，为光从真皮层进入表皮层的光透射率，σ_α ^epi为表皮层吸收系数，为表皮层散射系数，d_epi为表皮层厚度，L_air→L_epi代表光由空气进入表皮层，L_derm→L_epi代表光由真皮层进入表皮层。

进一步地，在步骤三中，只考虑光反射光由光从上方表皮层中进入的途径，所述在皮肤真皮层中光反射方程为其中，为从表皮层进入真皮层的光反射率，为真皮层吸收系数，为真皮层散射系数，d_derm为真皮层厚度，L_epi→L_derm代表光由上方表皮层进入真皮层。

进一步地，在步骤四中，所述光在皮肤表皮层和皮肤真皮层的两层结构中的反射率计算公式为：

$\begin{array}{c}{R}_{\mathrm{int}er}=R_{epi}^{+}+T_{epi}^{+}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}\\ *R_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}\\ *R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+...\end{array}$

其中，R_inter为光在皮肤表皮层和皮肤真皮层的两层结构中的反射率。

进一步地，在步骤五中，所述总体结构中的反射率公式为R_skin＝R_surf+T_surfR_inter。

一种根据人皮肤光谱的建模方法建立的高拟合度的多个皮肤参数的数学建模方法，，包括以下步骤：

步骤一建立表示在皮肤表面层中的所述反射率和皮肤粗糙度以及皮肤光亮度之间关系的方程；

步骤二建立表示在皮肤表皮层中的所述吸收系数与表皮层黑色素体积分数、表皮层优黑素浓度、表皮层褐黑素浓度、表皮层水分体积分数、表皮层油脂体积分数、表皮层中胡萝卜素浓度、表皮层皮层厚度之间关系的方程；

步骤三建立表示在皮肤真皮层中的所述吸收率与真皮层水分体积分数、血液体积分数、血红蛋白浓度、血液中氧化血红蛋白体积分数、血液中脱氧血红蛋白体积分数、血液中一氧化碳血红蛋白体积分数、血液中高铁血红蛋白体积分数、血液中硫化血红蛋白体积分数、真皮层胡萝卜素浓度、真皮层胆红素浓度、血液中血小板体积分数、血液中血红蛋白体积分数、真皮层弹性蛋白体积分数之间关系的方程；

步骤四建立表示在皮肤真皮层中所述散射系数与皮肤胶原蛋白体积分数、皮肤胶原蛋白纤维半径或直径、皮肤胶原蛋白纤维束半径或直径之间关系的方程；

步骤五建立表示皮肤表皮层中光反射率、光透射率与表皮厚度之间关系的方程；

步骤六建立表示皮肤真皮层中光反射率与表皮厚度之间关系的方程。

进一步地，还包括步骤七：将步骤一至步骤四中计算的皮肤参数与样本实际的皮肤参数比较，求拟合度，判断所述多个皮肤参数数学建模方法是否精确。

进一步地，在步骤一中，所述方程为

其中f_surf表示皮肤表面对单一方向的反射率，d_ω0表示反射方向的微元。

进一步地，在步骤二中，所述方程为

$\begin{array}{c}{\mathrm{\σ}}_a^{epi}=f_{me}(c_{eu}{\mathrm{\σ}}_a^{eu}+c_{ph}{\mathrm{\σ}}_a^{ph})+(1-f_{me})f_{h_{2}o}^{epi}{\mathrm{\σ}}_a^{h_{2}o}+(1-f_{me})(1-f_{h_{2}o}^{epi})f_{lipid}{\mathrm{\σ}}_a^{lipid}+(1-\\ f_{me})(1-f_{h_{2}o}^{epi})(1-f_{lipid})(c_{\mathrm{\β}c}^{epi}{\mathrm{\σ}}_a^{\mathrm{\β}c}+{\mathrm{\σ}}_a^{baseline})\end{array}$

其中，表示表皮层吸收系数，f_me表示表皮层黑色素体积分数，c_eu表示表皮层优黑素浓度，表示优黑素吸收系数，c_ph表示褐黑素浓度，表示褐黑素吸收系数，表示表皮水分体积分数，表示水分吸收系数，f_lipid表示表皮层油脂体积分数，表示油脂吸收系数，表示表皮层胡萝卜素浓度，表示胡萝卜素吸收系数，表示皮肤基线吸收系数。

进一步地，在步骤三中，所述方程为

$\begin{array}{c}{\mathrm{\σ}}_a^{derm}=f_{h_{2}o}^{derm}{\mathrm{\σ}}_a^{h_{2}o}\\ +(1-f_{h_{2}o}^{derm})f_{blood}(S_{oxy}{c}_{blood}{\mathrm{\σ}}_a^{oxy}+(1-S_{oxy})S_{deoxy}{c}_{blood}{\mathrm{\σ}}_a^{deoxy}\\ +(1-S_{oxy})(1-S_{deoxy})S_{co}{c}_{blood}{\mathrm{\σ}}_a^{co}\\ +(1-S_{oxy})(1-S_{deoxy})(1-S_{co})S_{met}{c}_{blood}{\mathrm{\σ}}_a^{met}\\ +(1-S_{oxy})(1-S_{deoxy})(1-S_{co})(1-S_{met})S_{sulf}{c}_{blood}{\mathrm{\σ}}_a^{sulf}+c_{\mathrm{\β}c}^{derm}{\mathrm{\σ}}_a^{\mathrm{\β}c}\\ +c_{br}{\mathrm{\σ}}_a^{br}+f_{plt}{\mathrm{\σ}}_a^{plt})+(1-f_{h_{2}o}^{derm})(1-{\mathrm{Hf}}_{blood})f_{ela}{\mathrm{\σ}}_a^{ela}\\ +(1-f_{h_{2}o}^{derm})(1-{\mathrm{Hf}}_{blood})(1-f_{ela}){\mathrm{\σ}}_a^{baseline}\end{array}$

其中，表示真皮层吸收系数，表示真皮层水分体积分数，表示水分吸收系数，f_blood表示血液体积分数，S_oxy表示血液中含氧血红蛋白体积分数即血氧浓度，c_blood表示血红蛋白浓度，表示含氧血红蛋白吸收系数，S_deoxy表示血液中脱氧血红蛋白体积分数，表示脱氧血红蛋白吸收系数，S_co表示血液中一氧化碳血红蛋白体积分数，表示一氧化碳血红蛋白吸收系数，S_met表示血液中高铁血红蛋白体积分数，表示高铁血红蛋白吸收系数，S_sulf表示血液中硫化血红蛋白体积分数，表示硫化血红蛋白吸收系数，表示真皮层胡萝卜素浓度，表示胡萝卜素吸收系数，c_br表示胆红素浓度，表示胆红素吸收系数，f_plt表示血液中血小板体积分数，表示血小板吸收系数，H表示血液中血红蛋白体积分数，f_ela表示弹性蛋白体积分数，表示弹性蛋白吸收系数，表示皮肤基线吸收系数。

进一步地，在步骤四中，所述方程为

其中，σ_s表示真皮层散射系数，为真皮层瑞利散射系数，为真皮层米氏散射系数。

进一步地，在步骤五中，所述方程为和

其中，为从空气进入表皮层的光反射率，为从空气进入表皮层的光透射率，为从真皮层进入表皮层的光反射率，为光从真皮层进入表皮层的光透射率，σ_α ^epi为表皮层吸收系数，为表皮层散射系数，d_epi为表皮层厚度，L_air→L_epi代表光由空气进入表皮层，L_derm→L_epi代表光由真皮层进入表皮层。

进一步地，在步骤六中，所述方程为

$R_{derm}^{+}=LSI({\mathrm{\σ}}_a^{derm},{\mathrm{\σ}}_s^{derm},d_{derm},L_{epi}\→L_{derm})$

其中，为从表皮层进入真皮层的光反射率，为真皮层吸收系数，为真皮层散射系数，d_derm为真皮层厚度，L_epi→L_derm代表光由上方表皮层进入真皮层。

采用上述技术方案，本发明具有如下有益效果：

第一，本发明通过光在皮肤中的传导分析，构建了皮肤光谱模型和皮肤参数数学模型，建立了皮肤光谱与皮肤生物学参数之间的联系。利用本发明的数学模型计算出的一组皮肤生物学参数仿真的虚拟光谱与实际光谱高度拟合，模型准确可靠。

第二，本发明相关的皮肤的光谱吸收相关的参数19个，光谱散射相关的参数3个，皮肤表面漫反射相关的参数2个，皮肤厚度相关参数2个，共计26维特征向量，能够仿真皮肤光谱的细节，模拟真实的皮肤情况，精确度高，从而实现了皮肤生物学参数量化分析计算的目的，并可作为皮肤大数据的处理的基础。

第三，本发明所建立的皮肤表皮层、真皮层吸收系数与皮肤光吸收相关的生物学参数之间的数学模型中，引入了如下形式的“余量体积分数”进行归一化处理，使得在优化分析过程中能够得到合理的最优解。

${\mathrm{\&sigma;}}_a^{epi/derm}=\underset{i\&Element;Absorbant}{\&Sigma;}\&lsqb;\underset{j<i}{\mathrm{\&Pi;}}(1-f_j)\&rsqb;f_i{\mathrm{\&sigma;}}_a^i+\underset{i\&Element;Absorbant}{\mathrm{\&Pi;}}(1-f_i){\mathrm{\&sigma;}}_a^{baseline}$

该方法有3个作用：1、建立了皮肤某一吸收成分的体积分数对其他体积分数影响的数学关系，保证了皮肤各个吸收成分体积分数的物理约束；2、排除了搜索方向的矫正所造成的影响，保证了量化分析过程不会受到外部干扰；3、确保了皮肤光吸收相关的生物学参数量化分析结果的正确性。

第四，本发明的方法在可见光400-700nm的全波段内精确拟合，可以达到皮肤量化分析所需要的精度。

第五，本发明所基于的皮肤光谱信息可由无创的采集方式得到。

附图说明

图1为本发明拟合度的结构示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，下面结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的结构图及具体实施例仅用以解释本发明，并不用于限定本发明。

实施例1

基于辐射传导理论方程应用的需要，本发明构建了在可见光(400-700nm，作用于人的皮肤深度最大4mm)作用下皮肤的四层结构(皮肤表面、表皮层、真皮层、皮下组织)模型，抽象了光的四种作用形式(反射、透射、吸收、散射)在皮肤四层结构模型上的辐射传导路径分析。综合形成了皮肤光谱与反射、吸收、散射、透射之间的光学分析模型。该模型有两个特点：1、在皮肤的结构模型中，将皮肤的表皮层、真皮层厚度作为结构变量，突出了这两个变量对光谱的影响；2、该模型涵盖了光与皮肤的四种作用形式，是完整的光学模型。

在现有技术中，辐射传导方程(RTE，RadiativeTransferEquation)，该理论叙述了电磁波在介质中传播时，电磁波会：在“吸收”(Absorption)中损失能量，在“激发”(Emission)中获得能量，在“散射”(Scattering)中重新分配能量。

电磁波在介质传播中，其在位置x沿方向传播的能量可由辐射率表示。一束沿方向传播的能量在经过位置在x时，在经过一小段介质(dS(x))后，辐射率的变化量可写为：

$\begin{array}{c}dL(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})=-{\mathrm{\&sigma;}}_a(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})L(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})dS\left(x\right)+q(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})dS\left(x\right)-{\mathrm{\&sigma;}}_s(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})L(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})dS\left(x\right)\\ +\underset{\mathrm{\&Omega;}}{\&Integral;}{\mathrm{\&sigma;}}_s(x,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}^{\&prime;})p(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}^{\&prime;})L(x,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}^{\&prime;}){\mathrm{d\&omega;}}^{\&prime;}dS\left(x\right)\end{array}$

辐射传导方程为：

$\begin{array}{c}(\stackrel{\&RightArrow;}{\&dtri;}\&CenterDot;\stackrel{\&RightArrow;}{\mathrm{\&omega;}})L(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})=\frac{dL(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})}{dS\left(x\right)}=-({\mathrm{\&sigma;}}_a+{\mathrm{\&sigma;}}_s)(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})L(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})dS\left(x\right)+q(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}})dS\left(x\right)+\\ {\&Integral;}_{\mathrm{\&Omega;}}{\mathrm{\&sigma;}}_s(x,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}^{\&prime;})p(x,\stackrel{\&RightArrow;}{\mathrm{\&omega;}},{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}^{\&prime;})L(x,{\stackrel{\&RightArrow;}{\mathrm{\&omega;}}}^{\&prime;}){\mathrm{d\&omega;}}^{\&prime;}dS\left(x\right)\end{array}$

根据理论方程，建立皮肤中的光辐射传导方程，皮肤模型假设：

1)皮肤被抽象为多层结构，每一层拥有独立、多种吸收、散射介质；

2)皮肤在垂直于厚度延伸方向是无穷的，在平行于厚度延伸的方向可以有穷或无穷；

3)皮肤为平面平行结构，即皮肤的光学性质只在平行于厚度延伸方向改变，位于同一深度的皮肤的光学特性完全一致；

4)每层的散射、吸收介质都是均匀分布的，即皮肤任意一层中，任何位置的光学特性都完全相同，拥有相同的含量和吸收、散射系数；

5)皮肤中各种介质的吸收、散射均相互独立；

6)皮肤中散射介质为小尺度的胶原蛋白纤维及大尺度的胶原蛋白纤维束，其中小尺度胶原蛋白纤维的散射可以近似为球状瑞利散射，造成真皮层同时拥有球形瑞利和圆柱形米氏散射。

建立皮肤结构中光传播模型的目的是为了能够计算皮肤总体反射率，皮肤总体反射率的计算由3部分组成：1、皮肤表面反射和透射；2、皮肤内部多层结构的总体反射；3、皮肤总体结构的反射。

本发明人皮肤光谱的建模方法包括以下步骤，

步骤一根据可见光照射皮肤时吸收、反射、散射、透射的特征，将皮肤从上之下抽象为皮肤粗面层、皮肤表皮层、皮肤真皮层、皮下组织层四层。

表面层：无实际厚度，无限薄，位于最外层，上与外界环境、下与表皮层相连；

表皮层：为皮肤实际的第一层，拥有有限厚度，上连粗糙表面层、下与真皮层连接；

真皮层：为皮肤实际的第二层，拥有有限厚度，上连表皮层、下与皮下组织连接；

皮下组织层：皮下组织拥有无穷厚度，并不描述任何组成，其作用为吸收所有进入从真皮层透射进入皮下组织的光。

步骤二根据辐射传导方程和在皮肤表面层的光吸收、散射系数，建立在皮肤表面层中光反射和光透射的方程，计算在皮肤表面层的反射率和吸收率。

皮肤表面的反射和透射由表面的光反射模型计算方程为：

[R_surf,T_surf]＝LSI(r_surf)，(1)

公式(1)其中，r_surf为皮肤表面的反射系数，由表面粗糙度σ_surf和表面光亮度ρ_surf组成。

由于该层并不吸收任何光的能量，故其反射率和透射率拥有如下关系：

T_surf＝1-R_surf。(2)

步骤三根据皮肤表皮层的光吸收、散射系数，建立在皮肤表皮层中光反射和光透射方程，计算在皮肤表皮层的反射率和吸收率。

表皮层能够对进入该层的光吸收和散射，该层厚度会对在传播中光被吸收、散射的总量造成影响，且表皮上层(空气)和表皮下层(真皮)的折射率，以及光从何层进入表皮层也会对该层反射、透射造成影响，故该层的反射、透射率的计算需要同时计算上照明(光从空气进入，+号表示)和下照明(光从真皮层进入，-号表示)两种情况：

$\&lsqb;R_{epi}^{+},T_{epi}^{+}\&rsqb;=LSI({\mathrm{\&sigma;}}_a^{epi},{\mathrm{\&sigma;}}_s^{epi},d_{epi},L_{air}\&RightArrow;L_{epi});---\left(3\right)$

公式(3)，其中为表皮层吸收系数，为表皮层散射系数，d_epi为表皮层厚度，L_air→L_epi代表光从空气中进入表皮层。

$\&lsqb;R_{epi}^{-},T_{epi}^{-}\&rsqb;=LSI({\mathrm{\&sigma;}}_a^{epi},{\mathrm{\&sigma;}}_s^{epi},d_{epi},L_{derm}\&RightArrow;L_{epi}).---\left(4\right)$

公式(4)，其中为表皮层吸收系数，为表皮层散射系数，d_epi为表皮层厚度，L_derm→L_epi代表光从真皮层中进入表皮层。

步骤四根据皮肤真皮层的光吸收、散射系数，建立在皮肤真皮层中光反射方程，计算计算在皮肤真皮层的反射率。

真皮层也会对进入的光吸收和散射，真皮厚度对光被吸收、散射的总量造成影响，但由于假设了进入皮下组织的光完全被吸收，不会回到真皮层，故不用考虑下照明(光从皮下组织进入)以及真皮透射，公式如下：

$R_{derm}^{+}=LSI({\mathrm{\&sigma;}}_a^{derm},{\mathrm{\&sigma;}}_s^{derm},d_{derm},L_{epi}\&RightArrow;L_{derm}).---\left(5\right)$

公式(5)，其中为真皮层吸收系数，为真皮层散射系数，d_derm为真皮层厚度，L_epi→L_derm代表光从表皮层中进入真皮层。

步骤五计算光在皮肤表皮层和皮肤真皮层的两层结构中的反射率。由于皮肤内部是多层结构，故光能在表皮层和真皮层之间反射，光可能在被多次内部反射之后才离开皮肤，多层结构的反射率计算如下：

$R_{\mathrm{int}er}=R_{epi}^{+}+T_{epi}^{+}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*$

$\begin{array}{c}{R}_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}*\\ R_{derm}^{+}*T_{epi}^{-}+...\end{array}---\left(6\right)$

公式(6)，其中R_inter为皮肤表皮层和真皮层的两层结构的反射率，为光从空气进入表皮层的反射率，为光从空气进入表皮层的透射率，为光从真皮层进入表皮层的反射率，为光从真皮层进入表皮层的透射率，为光从表皮层进入真皮层的反射率。

步骤六计算光通过四层皮肤的总体结构中的反射率。计算完皮肤表面反射、透射，以及皮肤内部总体反射后，可以计算皮肤总体反射。光在照射到皮肤后，首先一部分光被皮肤表面反射，未被反射的光全部进入皮肤内部，并与皮肤内部相互作用(吸收、散射)，最后未被皮肤内部吸收的光通过反射逃离皮肤内部：

R_skin＝R_surf+T_surfR_inter(7)

公式(7)，其中R_skin为皮肤四层结构的总反射率，R_surf为皮肤表面的反射率，T_surf为皮肤表面的透射率，R_inter为皮肤表皮层和真皮层两层结构的反射率。

实施例2

在上述皮肤光谱特征模型的基础上，找到与光敏感有关的皮肤生物学参数，形成了与皮肤结构模型相对应的皮肤参数集，建立了皮肤参数、皮肤结构、光作用形式的关联。进而通过对三种光作用形式与皮肤参数的建模，建立了一组数学表达式，搭建了皮肤参数与光学性质之间的关系，在本发明中尽可能多地涉及相关皮肤生物学参数，使模型数据量充足，提高算法精确度即高拟合度的多个皮肤参数的数学建模方法，其步骤如下：

步骤一建立表示在皮肤粗糙表面层中的所述反射率和皮肤粗糙度以及皮肤光亮度之间关系的方程。

皮肤表面被抽象成无穷多的对称的微观V形沟槽，每个V形沟槽表面都对光造成镜面反射，并可以通过几何光学进行计算。最终皮肤表面的反射则由对所有的微观V形沟槽的反射进行统计得到。

V形沟槽的分布及每个表面反射强度对表面总反射均造成影响，而其V形沟槽相互之间也会对入射光、反射光造成遮挡，故将表面单一方向反射f_surf建模为表面反射强度ρ_surf、V形沟槽分布D，几何衰减G以及菲涅尔表面漫反射F的乘积：

$f_{surf}={\mathrm{\&rho;}}_{surf}\frac{F({\mathrm{\&omega;}}_o,h)G({\mathrm{\&omega;}}_i,{\mathrm{\&omega;}}_o,h)D(h,n,{\mathrm{\&sigma;}}_{surf})}{4\mathrm{\&pi;}(n\&CenterDot;{\mathrm{\&omega;}}_i)}---\left(8\right)$

公式(8),其中ω_i,ω_o是入射光、观察方向的反方向，n是宏观表面法向量，h是ω_i,ω_o的半向量。

微观V形沟槽是通过参数化随机模拟的，由所有沟槽平均深度—间隔比来控制，通常使用Beckmann分布，最终的皮肤表面向各个方向的总反射则是对单一方向反射在所有观察角度上积分：

R_surf＝∫_2πf_surfdω_o(9)

公式(9)，其中R_surf为皮肤表面反射率，f_surf为单一观察角度上的反射率，dω_o为观察角度的微元。

步骤二建立表示在皮肤表皮层中的所述吸收系数与表皮层黑色素体积分数、表皮层优黑素浓度、表皮层褐黑素浓度、表皮层水分体积分数、表皮层油脂体积分数、表皮层中胡萝卜素浓度之间关系的方程。

由皮肤的假设可知，皮肤中各个吸收成分相互独立，因此吸收系数可以表示为每个成分吸收的线性组合，表皮层吸收系数与对应成分参数的数学关系可以如下：

$\begin{array}{c}{\mathrm{\&sigma;}}_a^{epi}=f_{me}(c_{eu}{\mathrm{\&sigma;}}_a^{eu}+c_{ph}{\mathrm{\&sigma;}}_a^{ph})+(1-f_{me})f_{h_{2}o}^{epi}{\mathrm{\&sigma;}}_a^{h_{2}o}+(1-f_{me})(1-f_{h_{2}o}^{epi})f_{lipid}{\mathrm{\&sigma;}}_a^{lipid}\\ +(1-f_{me})(1-f_{h_{2}o}^{epi})(1-f_{lipid})(c_{\mathrm{\&beta;}c}^{epi}{\mathrm{\&sigma;}}_a^{\mathrm{\&beta;}c}+{\mathrm{\&sigma;}}_a^{baseline})\end{array}---\left(10\right)$

公式(10)，其中表示表皮层吸收系数，f_me表示表皮层黑色素体积分数，c_eu表示表皮层优黑素浓度，表示优黑素吸收系数，c_ph表示褐黑素浓度，表示褐黑素吸收系数，表示表皮水分体积分数，表示水分吸收系数，f_lipid表示表皮层油脂体积分数，表示油脂吸收系数，表示表皮层胡萝卜素浓度，表示胡萝卜素吸收系数，表示皮肤基线吸收系数。

步骤三建立表示在皮肤真皮层中的所述吸收系数与真皮层水分体积分数、血液体积分数、血红蛋白浓度、血液中氧化血红蛋白体积分数、血液中脱氧血红蛋白体积分数、血液中一氧化碳血红蛋白体积分数、血液中高铁血红蛋白体积分数、血液中硫化血红蛋白体积分数、真皮层胡萝卜素浓度、真皮层胆红素浓度、血液中血小板体积分数、血液中血红蛋白体积分数、真皮层弹性蛋白体积分数之间关系的方程。

$\begin{array}{c}{\mathrm{\&sigma;}}_a^{derm}=f_{h_{2}o}^{derm}{\mathrm{\&sigma;}}_a^{h_{2}o}+(1-f_{h_{2}o}^{derm})f_{blood}(S_{oxy}{c}_{blood}{\mathrm{\&sigma;}}_a^{oxy}+(1-\\ S_{oxy})S_{deoxy}{c}_{blood}{\mathrm{\&sigma;}}_a^{deoxy}+(1-S_{oxy})(1-S_{deoxy})S_{co}{c}_{blood}{\mathrm{\&sigma;}}_a^{co}+(1-S_{oxy})(1-\\ S_{deoxy})(1-S_{co})S_{met}{c}_{blood}{\mathrm{\&sigma;}}_a^{met}+(1-S_{oxy})(1-S_{deoxy})(1-S_{co})(1-\\ S_{met})S_{sulf}{c}_{blood}{\mathrm{\&sigma;}}_a^{sulf}+c_{\mathrm{\&beta;}c}^{derm}{\mathrm{\&sigma;}}_a^{\mathrm{\&beta;}c}+c_{br}{\mathrm{\&sigma;}}_a^{br}+f_{plt}{\mathrm{\&sigma;}}_a^{plt})+(1-f_{h_{2}o}^{derm})(1-\\ {\mathrm{Hf}}_{blood})f_{ela}{\mathrm{\&sigma;}}_a^{ela}+(1-f_{h_{2}o}^{derm})(1-{\mathrm{Hf}}_{blood})(1-f_{ela}){\mathrm{\&sigma;}}_a^{baseline}\end{array}---\left(11\right)$

公式(11)，其中表示真皮层吸收系数，表示真皮层水分体积分数，表示水分吸收系数，f_blood表示血液体积分数，S_oxy表示血液中含氧血红蛋白体积分数(血氧浓度)，c_blood表示血红蛋白浓度，表示含氧血红蛋白吸收系数，S_deoxy表示血液中脱氧血红蛋白体积分数，表示脱氧血红蛋白吸收系数，S_co表示血液中一氧化碳血红蛋白体积分数，表示一氧化碳血红蛋白吸收系数，S_met表示血液中高铁血红蛋白体积分数，表示高铁血红蛋白吸收系数，S_sulf表示血液中硫化血红蛋白体积分数，表示硫化血红蛋白吸收系数，表示真皮层胡萝卜素浓度，表示胡萝卜素吸收系数，c_br表示胆红素浓度，表示胆红素吸收系数，f_plt表示血液中血小板体积分数，表示血小板吸收系数，H表示血液中血红蛋白体积分数，f_ela表示弹性蛋白体积分数，表示弹性蛋白吸收系数，表示皮肤基线吸收系数。

步骤四建立表示在皮肤真皮层中所述散射系数与皮肤胶原蛋白体积分数、皮肤胶原蛋白纤维半径、皮肤胶原蛋白纤维束直径之间关系的方程；步骤五建立表示皮肤表皮层中光反射率、光透射率与表皮厚度之间关系的方程。

散射是光在与皮肤相互作用中，光的传播方向被改变而产生的现象，方向的改变并不遵循Snell定律(折射定律)。皮肤中的主要散射成分为胶原蛋白，并且皮肤具有两种不同的散射：瑞利散射以及米氏散射。瑞利散射主要由小尺度的胶原蛋白纤维结构产生，而米氏散射主要有大尺度的胶原蛋白纤维束(由多束具有相同走向的纤维绑定而成)产生，真皮的散射系数由3个参数表示：胶原蛋白体积分数，胶原蛋白纤维半径，胶原蛋白纤维束半径。根据生物学的研究，本皮肤模型中表皮的散射系数使用真皮散射系数近似，即两层拥有相同的散射系数。

${\mathrm{\&sigma;}}_s={\mathrm{\&sigma;}}_s^{Rayleigh}+{\mathrm{\&sigma;}}_s^{Mie}---\left(12\right)$

公式(12)，其中σ_s为真皮层散射系数，为真皮层瑞利散射系数，为真皮层米氏散射系数。

散射理论主要描述了介质的散射系数σ_s可以通过单个胶原蛋白截面散射系数σ_ss和其数量密度N的乘积获得，而胶原蛋白截面散射系数σ_ss则可以通过散射效率σ_sca乘以其截面面积A获得。

$\left\{\begin{array}{c}{\mathrm{\&sigma;}}_{ss}={\mathrm{\&sigma;}}_{sca}A\\ {\mathrm{\&sigma;}}_s={\mathrm{\&sigma;}}_{ss}N\end{array}\right.---\left(13\right)$

根据物理瑞利散射理论,米氏散射方程，计算米氏散射效率，计算各向异性指数，得到米氏散射系数，最终得到公式(12)。

虽然表皮层中并不含有胶原蛋白，但根据生物—物理学家的研究结果，表皮层散射特性与真皮层散射特性无异，故表皮层散射也由真皮层散射近似。

步骤五建立表示皮肤表皮层中光反射率、光透射率与表皮厚度之间关系的方程。

步骤六建立表示皮肤真皮层中光反射率与表皮厚度之间关系的方程。

实施例3

优选地，在本发明高拟合度的多个皮肤参数的数学建模方法中，还包括步骤七：将步骤一至步骤四中计算的皮肤参数与样本实际的皮肤参数比较，计算拟合度，判断上述多个皮肤参数的数学建模方法是否模拟精确，可以用于皮肤参数的量化分析。具体步骤包括1.测量目标皮肤各个波段的反射率样本，并进行预处理以增加平滑度和减少测量噪音；2.猜测一组初始皮肤生物学参数；3.根据所述步骤，分别计算皮肤表面、表皮层、真皮层的反射率及透射率；4.根据所述步骤，计算皮肤表皮层、真皮层以及皮下组织共三层结构的反射率；5.计算皮肤四层结构的反射率；6.比较计算得到的皮肤反射率和测量的皮肤反射率之间的距离；7.计算所选的参数组中，各个参数变化时计算得到的皮肤反射率和测量的皮肤反射率之间距离的变化；8.根据距离的变化选择下一组皮肤生物学参数，并重复步骤3至步骤8，直到得到拥有最小距离的一组生物学参数。其拟合曲线如图1所示，由拟合曲线可知，本方法拟合程度非常高。

以上所述实施例仅表达了本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (15)

1.一种人皮肤光谱的建模方法，其特征在于：所述方法包括以下步骤：

步骤一根据可见光照射皮肤时吸收、反射、散射、透射的特征，将皮肤从上至下抽象为皮肤表面层、皮肤表皮层、皮肤真皮层、皮下组织层四层；

步骤二根据皮肤表皮层的光吸收、散射系数以及皮肤表皮层厚度变量参数，建立在皮肤表皮层中光反射和光透射方程，计算在皮肤表皮层的反射率和透射率；

步骤三根据皮肤真皮层的光吸收、散射系数以及皮肤真皮层厚度变量参数，建立在皮肤真皮层中光反射方程，计算在皮肤真皮层的反射率；

步骤四计算光在皮肤表皮层和皮肤真皮层的两层结构中的反射率；

步骤五计算光通过四层皮肤的总体结构中的反射率。

2.根据权利要求1所述的人皮肤光谱的建模方法，其特征在于，还包括根据辐射传导方程和在皮肤表面层的反射系数，建立在皮肤表面层中光反射和光透射的方程，计算在皮肤表面层的反射率的步骤。

3.根据权利要求2所述的人皮肤光谱的建模方法，其特征在于，所述在皮肤表面层中光反射和光透射的方程为[R_surf,T_surf]＝LSI(r_surf)，其中，R_surf为表面的光反射率，T_surf为表面的光透射率，LSI代表了用于计算皮肤反射率以及透射率的一系列计算公式，r_surf为皮肤表面的反射系数。

4.权利要求1所述的人皮肤光谱的建模方法，其特征在于，在步骤二中，在所述皮肤表皮层中光反射和光透射由光从上方空气中进入和光从下方真皮层进入两种途径所决定，所述在皮肤表皮层中光反射和光透射的方程为

和

其中，为从空气进入表皮层的光反射率，为从空气进入表皮层的光透射率，为从真皮层进入表皮层的光反射率，为光从真皮层进入表皮层的光透射率，σ_α ^epi为表皮层吸收系数，为表皮层散射系数，d_epi为表皮层厚度，L_air→L_epi代表光由空气进入表皮层，L_derm→L_epi代表光由真皮层进入表皮层。

5.根据权利要求1所述的人皮肤光谱的建模方法，其特征在于，在步骤三中，只考虑光反射光由光从上方表皮层中进入的途径，所述在皮肤真皮层中光反射方程为

其中，为从表皮层进入真皮层的光反射率，为真皮层吸收系数，为真皮层散射系数，d_derm为真皮层厚度，L_epi→L_derm代表光由上方表皮层进入真皮层。

6.根据权利要求1所述的人皮肤光谱的建模方法，其特征在于，在步骤四中，所述光在皮肤表皮层和皮肤真皮层的两层结构中的反射率计算公式为：

$\begin{array}{c}{R}_{\mathrm{int}er}=R_{epi}^{+}+T_{epi}^{+}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}\\ *R_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+T_{epi}^{+}*R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*R_{epi}^{-}\\ *R_{derm}^{+}*R_{epi}^{-}*R_{derm}^{+}*T_{epi}^{-}+...\end{array}$

其中，R_inter为光在皮肤表皮层和皮肤真皮层的两层结构中的反射率。

7.根据权利要求1所述的人皮肤光谱的建模方法，其特征在于，在步骤五中，所述总体结构中的反射率公式为

R_skin＝R_surf+T_surfR_inter。

8.一种根据权利要求1所述的人皮肤光谱的建模方法，建立的高拟合度的多个皮肤参数的数学建模方法，其特征在于，包括以下步骤：

步骤一建立表示在皮肤表面层中的所述反射率和皮肤粗糙度以及皮肤光亮度之间关系的方程；

步骤二建立表示在皮肤表皮层中的所述吸收系数与表皮层黑色素体积分数、表皮层优黑素浓度、表皮层褐黑素浓度、表皮层水分体积分数、表皮层油脂体积分数、表皮层中胡萝卜素浓度、表皮层皮层厚度之间关系的方程；

步骤三建立表示在皮肤真皮层中的所述吸收率与真皮层水分体积分数、血液体积分数、血红蛋白浓度、血液中氧化血红蛋白体积分数、血液中脱氧血红蛋白体积分数、血液中一氧化碳血红蛋白体积分数、血液中高铁血红蛋白体积分数、血液中硫化血红蛋白体积分数、真皮层胡萝卜素浓度、真皮层胆红素浓度、血液中血小板体积分数、血液中血红蛋白体积分数、真皮层弹性蛋白体积分数之间关系的方程；

步骤四建立表示在皮肤真皮层中所述散射系数与皮肤胶原蛋白体积分数、皮肤胶原蛋白纤维半径或直径、皮肤胶原蛋白纤维束半径或直径之间关系的方程；

步骤五建立表示皮肤表皮层中光反射率、光透射率与表皮厚度之间关系的方程；

步骤六建立表示皮肤真皮层中光反射率与表皮厚度之间关系的方程。

9.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，还包括步骤七：将步骤一至步骤六中通过所述皮肤参数所虚拟的皮肤光谱与实际的皮肤光谱拟合。

10.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，在步骤一中，所述方程为

其中f_surf表示皮肤表面对单一方向的反射率，d_ω0表示反射方向的微元。

11.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，在步骤二中，所述方程为

$\begin{array}{c}{\mathrm{\&sigma;}}_a^{epi}=f_{me}\left(c_{eu}{\mathrm{\&sigma;}}_a^{eu}+c_{ph}{\mathrm{\&sigma;}}_a^{ph}\right)+\left(1-f_{me}\right)f_{h_{2}o}^{epi}{\mathrm{\&sigma;}}_a^{h_{2}o}+\left(1-f_{me}\right)\left(1-f_{h_{2}o}^{epi}\right)f_{lipid}{\mathrm{\&sigma;}}_a^{lipid}+(1-\\ f_{me}\left)\right(1-f_{h_{2}o}^{epi}\left)\right(1-f_{lipid}\left)\right(c_{\mathrm{\&beta;}c}^{epi}{\mathrm{\&sigma;}}_a^{\mathrm{\&beta;}c}+{\mathrm{\&sigma;}}_a^{baseline})\end{array}$

其中，表示表皮层吸收系数，f_me表示表皮层黑色素体积分数，c_eu表示表皮层优黑素浓度，表示优黑素吸收系数，c_ph表示褐黑素浓度，表示褐黑素吸收系数，表示表皮水分体积分数，表示水分吸收系数，f_lipid表示表皮层油脂体积分数，表示油脂吸收系数，表示表皮层胡萝卜素浓度，表示胡萝卜素吸收系数，表示皮肤基线吸收系数。

12.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，在步骤三中，所述方程为

其中，表示真皮层吸收系数，表示真皮层水分体积分数，表示水分吸收系数，f_blood表示血液体积分数，S_oxy表示血液中含氧血红蛋白体积分数即血氧浓度，c_blood表示血红蛋白浓度，表示含氧血红蛋白吸收系数，S_deoxy表示血液中脱氧血红蛋白体积分数，表示脱氧血红蛋白吸收系数，S_co表示血液中一氧化碳血红蛋白体积分数，表示一氧化碳血红蛋白吸收系数，S_met表示血液中高铁血红蛋白体积分数，表示高铁血红蛋白吸收系数，S_sulf表示血液中硫化血红蛋白体积分数，表示硫化血红蛋白吸收系数，表示真皮层胡萝卜素浓度，表示胡萝卜素吸收系数，c_br表示胆红素浓度，表示胆红素吸收系数，f_plt表示血液中血小板体积分数，表示血小板吸收系数，H表示血液中血红蛋白体积分数，f_ela表示弹性蛋白体积分数，表示弹性蛋白吸收系数，表示皮肤基线吸收系数。

13.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，在步骤四中，所述方程为

其中，σ_s表示真皮层散射系数，为真皮层瑞利散射系数，为真皮层米氏散射系数。

14.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，在步骤五中，所述方程为

和

$[R_{\mathrm{epi}}^{-},T_{\mathrm{epi}}^{-}]=\mathrm{LSI}({\mathrm{\&sigma;}}_a^{\mathrm{epi}},{\mathrm{\&sigma;}}_s^{\mathrm{epi}},d_{\mathrm{epi}},L_{\mathrm{derm}}\&RightArrow;L_{\mathrm{epi}}),$

其中，为从空气进入表皮层的光反射率，为从空气进入表皮层的光透射率，为从真皮层进入表皮层的光反射率，为光从真皮层进入表皮层的光透射率，σ_α ^epi为表皮层吸收系数，为表皮层散射系数，d_epi为表皮层厚度，L_air→L_epi代表光由空气进入表皮层，L_derm→L_epi代表光由真皮层进入表皮层。

15.根据权利要求8所述的高拟合度的多个皮肤参数的数学建模方法，其特征在于，在步骤六中，所述方程为

$R_{derm}^{+}=LSI({\mathrm{\&sigma;}}_a^{derm},{\mathrm{\&sigma;}}_s^{derm},d_{derm},L_{epi}\&RightArrow;L_{derm})$

其中，为从表皮层进入真皮层的光反射率，为真皮层吸收系数，为真皮层散射系数，d_derm为真皮层厚度，L_epi→L_derm代表光由上方表皮层进入真皮层。