JP5821436B2 - 皮膚色素濃度測定方法 - Google Patents
皮膚色素濃度測定方法 Download PDFInfo
- Publication number
- JP5821436B2 JP5821436B2 JP2011194043A JP2011194043A JP5821436B2 JP 5821436 B2 JP5821436 B2 JP 5821436B2 JP 2011194043 A JP2011194043 A JP 2011194043A JP 2011194043 A JP2011194043 A JP 2011194043A JP 5821436 B2 JP5821436 B2 JP 5821436B2
- Authority
- JP
- Japan
- Prior art keywords
- skin
- concentration
- pigment
- image
- spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000049 pigment Substances 0.000 title claims description 75
- 238000000691 measurement method Methods 0.000 title description 5
- XUMBMVFBXHLACL-UHFFFAOYSA-N Melanin Chemical compound O=C1C(=O)C(C2=CNC3=C(C(C(=O)C4=C32)=O)C)=C2C4=CNC2=C1C XUMBMVFBXHLACL-UHFFFAOYSA-N 0.000 claims description 82
- 108010054147 Hemoglobins Proteins 0.000 claims description 58
- 102000001554 Hemoglobins Human genes 0.000 claims description 58
- 238000000862 absorption spectrum Methods 0.000 claims description 37
- 238000001228 spectrum Methods 0.000 claims description 35
- 238000000034 method Methods 0.000 claims description 34
- 230000008033 biological extinction Effects 0.000 claims description 21
- 230000014509 gene expression Effects 0.000 claims description 21
- 230000003595 spectral effect Effects 0.000 claims description 18
- 238000000611 regression analysis Methods 0.000 claims description 16
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 14
- 238000004088 simulation Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 210000003491 skin Anatomy 0.000 description 80
- 239000000975 dye Substances 0.000 description 30
- 230000000052 comparative effect Effects 0.000 description 18
- 238000010521 absorption reaction Methods 0.000 description 13
- 238000005259 measurement Methods 0.000 description 11
- 210000004207 dermis Anatomy 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 7
- 108010064719 Oxyhemoglobins Proteins 0.000 description 5
- 210000002615 epidermis Anatomy 0.000 description 5
- 108010036302 hemoglobin AS Proteins 0.000 description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 3
- BPYKTIZUTYGOLE-IFADSCNNSA-N Bilirubin Chemical compound N1C(=O)C(C)=C(C=C)\C1=C\C1=C(C)C(CCC(O)=O)=C(CC2=C(C(C)=C(\C=C/3C(=C(C=C)C(=O)N\3)C)N2)CCC(O)=O)N1 BPYKTIZUTYGOLE-IFADSCNNSA-N 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 239000002537 cosmetic Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 229910052736 halogen Inorganic materials 0.000 description 2
- 150000002367 halogens Chemical class 0.000 description 2
- 230000031700 light absorption Effects 0.000 description 2
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 1
- 101100366935 Caenorhabditis elegans sto-2 gene Proteins 0.000 description 1
- UPYKUZBSLRQECL-UKMVMLAPSA-N Lycopene Natural products CC(=C/C=C/C=C(C)/C=C/C=C(C)/C=C/C1C(=C)CCCC1(C)C)C=CC=C(/C)C=CC2C(=C)CCCC2(C)C UPYKUZBSLRQECL-UKMVMLAPSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 230000003796 beauty Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 150000001746 carotenes Chemical class 0.000 description 1
- 235000005473 carotenes Nutrition 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000000245 forearm Anatomy 0.000 description 1
- 238000005534 hematocrit Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 210000002780 melanosome Anatomy 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 230000036284 oxygen consumption Effects 0.000 description 1
- 206010033675 panniculitis Diseases 0.000 description 1
- 230000000886 photobiology Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
- 210000004003 subcutaneous fat Anatomy 0.000 description 1
- 210000004304 subcutaneous tissue Anatomy 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- NCYCYZXNIZJOKI-UHFFFAOYSA-N vitamin A aldehyde Natural products O=CC=C(C)C=CC=C(C)C=CC1=C(C)CCCC1(C)C NCYCYZXNIZJOKI-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Description
本発明は、皮膚色素濃度の測定方法及びこの方法を用いた皮膚色素画像の形成方法に関する。
皮膚の見た目の色は、皮膚に含まれる色素の種類と濃度に依存する。そこで、従来より皮膚の色の観測に基づいて皮膚に含まれるメラニン、酸化ヘモグロビン、還元ヘモグロビンの濃度を測定する方法が提案されている。
例えば、ランベルト−ベールの法則から、皮膚の吸光度スペクトルは、皮膚に含まれるメラニン、酸化ヘモグロビン及び還元ヘモグロビンの吸光度スペクトルの線形結合であるとのモデルをたて、皮膚のある領域に光を入射させると共に、同領域からの反射光を受光することにより平均値として得た所謂一点計測の皮膚の吸光度スペクトルを、メラニン、酸化ヘモグロビン、還元ヘモグロビンの吸光度スペクトルを用いて重回帰分析することにより、皮膚中のメラニン濃度、酸化ヘモグロビン濃度、又は還元ヘモグロビン濃度を求める方法が知られている（特許文献１）。
しかしながら、このモデルは、皮膚中の色素の吸収のみによって吸光度スペクトルが定まることを前提としており、皮膚を構成する表皮又は真皮における散乱や、表皮、真皮、皮下脂肪の各界面での散乱が考慮されていないため、色素濃度の測定精度が低い。
これに対し、皮膚における散乱を考慮したモンテカルロシミュレーションモデルで補われた重回帰分析法が提案されている（非特許文献１）。この分析法では、モンテカルロシミュレーションモデルによる吸光度スペクトルを目的変数とし、メラニン、酸化ヘモグロビン、還元ヘモグロビンの吸光係数スペクトルを説明変数として重回帰分析して、回帰パラメータの切片ａ0、メラニンの回帰パラメータａm、及び総ヘモグロビンの回帰パラメータａthを得(但し、ａth＝ａoh＋ａdh、ａoh＝酸化ヘモグロビンの吸収係数スペクトルの回帰パラメータ、ａdh＝還元ヘモグロビン吸収スペクトルの回帰パラメータ）、これらａ0、ａm及びａthと、メラニン濃度Ｃm及び総ヘモグロビン濃度Ｃthとの関係式Ｃm＝Ｆm(ａ0,ａm,ａth)、Ｃth＝Ｆth(ａ0,ａm,ａth)を求めておく。
一方、任意の皮膚の吸光度スペクトルを、メラニン、酸化ヘモグロビン及び還元ヘモグロビンの吸光係数スペクトルεm(λ)、εoh(λ)及びεdh(λ)を説明変数として重回帰分析して回帰パラメータの切片ａ0、メラニン吸光係数の回帰パラメータａm、及び総ヘモグロビン吸光係数の回帰パラメータａthを得、このａ0、ａm、ａthと前述の関係式Ｃm＝Ｆm(ａ0,ａm,ａth)、Ｃth＝Ｆth(ａ0,ａm,ａth)とからメラニン濃度Ｃmと総ヘモグロビン濃度Ｃthを得る。
この方法において、酸素飽和度ＣStO2＝ａoh／(ａoh＋ａdh)である。
この方法において、酸素飽和度ＣStO2＝ａoh／(ａoh＋ａdh)である。
このモンテカルロシミュレーションモデルを用いた従来の色素濃度測定方法は、皮膚の吸光度スペクトルに散乱の影響を考慮しているので、色素濃度の測定精度が高い。しかしながら、ある程度の広がりをもった皮膚領域の色素濃度を調べるために、モンテカルロシミュレーションモデルを用いた従来の色素濃度測定方法を、皮膚画像を用いた色素濃度の測定方法に適用すると、肌の凹凸や照明ムラ等により画像面内に陰影が存在する場合に陰影の影響を受けて色素濃度の測定精度が低下する。
Journal of Biomedical Optics 9(4): 700-710 (2004)
本発明は上述の従来の課題を解決しようとするものであり、皮膚画像から色素濃度を測定するにあたり、画像面内に陰影が存在する場合にもその影響を受けることなく、正確に測定できる方法を提供すること、また、皮膚の色素濃度分布を正確に表す皮膚色素画像を提供することを目的とする。
本発明者は、皮膚画像に陰影がない場合、皮膚の撮像により得た吸光度スペクトルを、メラニンの吸光係数スペクトルεm(λ)、酸化ヘモグロビンの吸光係数スペクトルεoh(λ)及び還元ヘモグロビンの吸光係数スペクトルεdh(λ)で重回帰分析して得られる回帰パラメータの切片ａ0には、各色素の濃度情報や皮膚散乱情報が反映されるために、この回帰パラメータの切片ａ0を用いて回帰パラメータと色素濃度との関係式を求めると、色素濃度を正確に求めることができるのに対し、皮膚画像に陰影がある場合には、回帰パラメータの切片ａ0を使用することなく、回帰パラメータと色素濃度との関係式を求めることにより、回帰パラメータと色素濃度との関係から陰影情報の影響を排除できることを見出し、そして、任意の皮膚の吸光度スペクトルを色素の吸光係数スペクトルで重回帰分析して得られる回帰パラメータと、前述の関係式とから皮膚の色素濃度を求めることにより、皮膚画像に形成された陰影の影響を受けることなく、色素濃度を正確に求められることを見出した。
即ち、本発明は、皮膚の色素濃度から分光反射率を算出するシミュレーションモデルに基づき、複数の色素を含む皮膚モデルの吸光度スペクトルを、各色素の吸光係数スペクトルεi(λ)（色素の種類をＮとしてi=1,2,…,N、 以下同じ）を用いて、各色素濃度Ｃi（i=1,2,…,N）の組み合わせが異なる複数の場合について求め、得られた複数の皮膚モデルの吸光度スペクトルを、各色素の吸光係数スペクトルεi(λ)（i=1,2,…,N）を説明変数として重回帰分析することにより、以下の重回帰式(1)で得られる回帰パラメータａi（i=1,2,…,N）の組を複数求め、
また、本発明は、被験者の皮膚の内部反射光画像の各ピクセルについて、上述の方法により皮膚の各色素濃度Ｃi（i=1,2,…,N）の少なくとも一つの色素濃度を求め、その色素濃度の分布画像を形成する皮膚色素画像の形成方法を提供する。
さらに、本発明は、こうして形成した皮膚色素画像を提供する。
さらに、本発明は、こうして形成した皮膚色素画像を提供する。
本発明によれば、皮膚画像に基づき、その画像に形成された陰影の影響を受けることなく、皮膚のメラニン、酸化ヘモグロビン、還元ヘモグロビン等の各色素濃度分布を正確に求めることができる。したがって、皮膚を観察したときの見え方と皮膚中の色素濃度との関係がわかる。また、ヘモグロビンの酸素飽和度も容易に分析することが可能となり、皮膚における酸素の消費状態、血行状態などを推定することが可能となる。さらに、化粧料を適用した皮膚について、本発明の皮膚色素画像を経時的にとり、例えば、メラニン濃度、酸化ヘモグロビン、還元ヘモグロビンの各濃度変化を調べることにより、化粧料の有効性の評価を行うことが可能となる。また、本発明は、店頭での美容のための肌診断などでも有用となる。
以下、図面に基づき、本発明を詳細に説明する。なお、各図中、同一符号は、同一又は同等の構成要素を表している。
本発明の皮膚色素濃度の測定方法では、被験者の皮膚の吸光度スペクトルＡ(λ)から被験者の皮膚中の色素濃度を求めるために、まず、皮膚のシミュレーションモデルに基づき、皮膚に複数の色素が種々の濃度Ｃi（i=1,2,…,N）で含まれている場合の吸光度スペクトルＡ(λ)を求め、この吸光度スペクトルＡ(λ)を各色素の吸光係数スペクトルεi(λ)（i=1,2,…,N）を説明変数として重回帰分析する。例えば、皮膚中の色素として、メラニン、酸化ヘモグロビン及び還元ヘモグロビンを考える場合に、皮膚のシミュレーションモデルに基づく吸光度スペクトルＡ(λ)を、メラニンの吸光係数スペクトルεm(λ)、酸化ヘモグロビンの吸光係数スペクトルεoh(λ)及び還元ヘモグロビンの吸光係数スペクトルεdh(λ)を説明変数として重回帰分析して以下の重回帰式(1)aを得る。
そして、この式(1)aで得られる回帰パラメータａm、ａoh、ａdhとその切片ａoのうち、切片ａoを使用することなく、回帰パラメータａm、ａoh、ａdhからメラニン濃度Ｃm、酸化ヘモグロビン濃度Ｃoh及び還元ヘモグロビン濃度Ｃdhを求める次式(2)aの関係式Ｆm、Ｆoh、Ｆdhを求めておくことを特徴とする。
重回帰式(1)aの切片ａoには、メラニン、酸化ヘモグロビン、還元ヘモグロビンの各色素の濃度情報が含まれると考えられ、さらに、波長依存性がないことから、真皮内の光散乱や真皮と皮下組織との界面での皮膚散乱の情報が含まれる、と考えられている。このため、メラニン、ヘモグロビン等の色素の種々の濃度ｃi（i=1,2,…,N）の組み合わせについて吸光度スペクトルを求め、色素濃度と重回帰式(1)の回帰パラメータａi（i=1,2,…,N）あるいはその切片ａoとの関係式Ｆを求めておき、一方、任意の吸光度スペクトルを重回帰分析することにより回帰パラメータａiあるいはその切片ａoを得、得られた回帰パラメータａiあるいは切片ａo、と前述の関係式Ｆとから色素濃度ｃiを得る場合に、モンテカルロシミュレーションモデルを用い、回帰パラメータとしてメラニンの吸光係数スペクトルεm(λ)の回帰パラメータａmと総ヘモグロビンの吸光係数スペクトルεth(λ)の回帰パラメータａthを使用する従来の色素濃度測定方法（非特許文献１）では、関係式Ｆに、回帰パラメータの切片ａoが含まれていた。これに対し、本発明はこの関係式として、回帰パラメータの切片ａoを含まない前述の式(2)、より具体的には、色素として、メラニン、酸化ヘモグロビン、及び還元ヘモグロビンを考える場合に前述の次式(2)aの関係式Ｆm、Ｆoh、Ｆdhを使用する。
これにより、吸光度スペクトルを撮像画像から求める場合に、その画像に陰影があっても、それらによる画像の暗さが皮膚内部での吸収・散乱の影響と誤測定されることなく、正確に色素濃度を測定することが可能となる。これは、回帰パラメータの切片ａoは、波長依存性がないことから、画像に陰影ができた場合、陰影情報がこの切片ａoに集約されるためと考えられる。
また、モンテカルロシミュレーションモデルを用いた従来の色素濃度測定方法（非特許文献１）において色素濃度Ｃm、Ｃoh及びＣdhを得る場合に、単に切片ａ０を除くだけでは、色素濃度がＣm、Ｃoh及びＣdhの３種類であるのに対して、重回帰パラメータはａmとａthの２種類のみなので、３対２の対応となるために精度が低下してしまうが、ａthの代わりにａｏｈとａｄｈを用いると、３種類の重回帰パラメータを対応させることができるので、精度の低下を避けることができる。
本発明においては、上述のように、式(2)の関係式
このシミュレーションモデルとしては、例えば、モンテカルロシミュレーション（L.-H. Wang, S.L. Jacques, and L.-Q.Zheng,Computer Methods and Programs in Biomedicine 47, 131-146 (1995).）を用いることができる。モンテカルロシミュレーションは、皮膚に入射した光が、実効散乱係数と吸収係数に従った確率で散乱、吸収され、拡散反射されるとするモデルである。この色素として、例えば図２に示すように、メラニン、酸化ヘモグロビン、還元ヘモグロビンを考える場合、実効散乱係数スペクトル（μｓ'(λ)）と、メラニン、酸化ヘモグロビン、還元ヘモグロビンの各色素の吸光係数スペクトルεm(λ)、εoh(λ)、εdh(λ)と、各色素の濃度Ｃm、Ｃoh、Ｃdhから分光反射率Ｒ(λ)を求めることを可能とする。
より具体的には、この皮膚モデルにおいて、
表皮の減衰係数は、Ｃm×εm(λ) で表され、
真皮の減衰係数は、Ｃoh×εoh(λ)＋Ｃdh×εdh(λ) で表される。
表皮の減衰係数は、Ｃm×εm(λ) で表され、
真皮の減衰係数は、Ｃoh×εoh(λ)＋Ｃdh×εdh(λ) で表される。
実効散乱係数スペクトル（μｓ'（λ））や各色素の吸光係数スペクトルεm(λ)、εoh(λ)、εdh(λ)は文献値を用いることができることから、任意に各色素の濃度Ｃm、Ｃoh、Ｃdhを設定することにより、設定した色素濃度での分光反射率Ｒ(λ)を求めることができる。
なお、本発明において、皮膚の色素濃度から分光反射率を算出するシミュレーションモデルとしては、モンテカルロシミュレーションに限らず、例えば、クベルカ-ムンクの式 (Kubelka-Munk formula)、Adding-Doubling法 (Adding-Doubling method)、拡散理論 (Diffusion theory)、擬似皮膚を使った実測値に基づくモデルなどを使用してもよい。
現象記述の正確さと得られる複数のスペクトル間との整合性の高さから、モンテカルロシミュレーションに基づくモデルが好ましい。
現象記述の正確さと得られる複数のスペクトル間との整合性の高さから、モンテカルロシミュレーションに基づくモデルが好ましい。
こうして得られた分光反射率Ｒ(λ)は、吸光度スペクトルＡ＝log(1/Ｒ)により、吸光度スペクトルに変換し、さらにメラニン、酸化ヘモグロビン、還元ヘモグロビンの各色素の吸光係数スペクトルεm(λ)、εoh(λ)、εdh(λ)を説明変数として式(1)aのように重回帰分析し、回帰パラメータａm、ａoh及びａdhと切片ａ0を求める。
この色素の濃度Ｃm、Ｃoh及びＣdhと回帰パラメータａm、ａoh及びａdhを、種々の色素濃度Ｃm、Ｃoh及びＣdhの組み合わせについて取得し、色素濃度Ｃm、Ｃoh及びＣdhと回帰パラメータａm、ａoh及びａdhとの関係式Ｆm、Ｆoh、Ｆdhを求める。
関係式の具体的な取得方法としては、例えば、次式のように多項式近似をすることができる。
色素濃度Ｃm、Ｃoh及びＣdhと回帰パラメータａm、ａoh及びａdhとの関係式Ｆm、Ｆoh、Ｆdhを求める手法としては、多項式近似に代えて、ルックアップテーブルによる対応付け、直線近似、指数関数近似、対数関数近似等を用いても良い。
本発明において、上述の図２に示した関係式Ｆm、Ｆoh、Ｆdhの求め方は、皮膚内にさらに他の色素が存在する場合、例えば、カロテン、ビリルビン、ＡＧＥｓなどが存在する場合にも適用することができ、式(2)を求めることができる。
一方、本発明の方法においては、上述の関係式Ｆi(ａ1,ａ2,…,ａN)とは別に、被験者の皮膚の吸光度スペクトルを取得する。吸光度スペクトルの取得方法は、皮膚の狭い領域に所定波長の光を入射させると共にその領域から受光した反射光の平均値を用いる、所謂一点計測の分光反射率測定装置を使用してもよいが、皮膚の代謝や血行状態の観察のし易さの点から皮膚の内部反射光の撮像により求めることが好ましい。
例えば、図１に示す撮像系のように、被験者１に、第１の偏光フィルタ２ａを通して光源３から光を入射させ、第２の偏光フィルタ２ｂを備えた波長可変フィルタ４を通してデジタルカメラ５で皮膚画像を撮る。この場合、第１の偏光フィルタ２ａと第２の偏光フィルタ２ｂの偏光方向を直交させることにより、第２の偏光フィルタ２ａで被験者１の皮膚からの表面反射光を遮光し、皮膚の色情報を担う内部反射光のみがデジタルカメラ５で受光されるようにする。
光源３としては、白色光を発するハロゲンランプ、メタルハロゲンランプ等を使用する。光源３は、平行光を発するものでも拡散光を発するものでもよい。本発明の方法によれば、取得画像内の陰影に関係なく皮膚中の色素濃度を正確に測定できるため、光源としては、被写体に陰影を形成するものでも使用することができる。好ましくは、入射角０〜４５°の範囲で被験者１の皮膚に光を出射するものを使用する。
波長可変フィルタ４は、コントローラ６の制御により透過波長を変更する。これにより、例えば、波長５００ｎｍから２０ｎｍ間隔で撮像波長をずらしながら、内部反射光画像を繰り返し撮る。この場合、白色板の明るさを基準に、波長領域ごとの感度を揃えることが好ましい。
デジタルカメラ５は、全測定波長領域に渡って感度を有するカメラであれば、モノクロカメラで充分であるが、カラーカメラでも差し支えない。
こうしてデジタルカメラ５で撮影した画像は、パーソナルコンピュータ等の演算装置７に保存する。演算装置７は、各ピクセルの画素値から、各ピクセルごとに反射率を算出する。撮影にカラーカメラを用いた場合は、画素値としてＲＧＢ各画素値の重み付き平均を用いればよく、測定波長に対する感度を有する特定の画素のみを用いても良い。演算装置７において、画素値から分光反射率を算出させるためには、画像処理が可能なプログラミング環境を用いて、波長ごとに白色板の画素値を反射率に変換するための係数を求め、その係数を該波長画像全体に適用すればよい。このようなプログラミング環境としては、ＭＡＴＬＡＢ（Ｍａｔｈｗｏｒｋｓ社製）等を使用することができる。
演算装置７は、分光反射率を吸光度スペクトルに変換し、さらにそれを各色素の吸光係数スペクトルεi(λ)（i=1,2,…,N）を説明変数として重回帰分析することにより回帰パラメータａi（i=1,2,…,N）を算出する。
また、この演算装置７には、前述の色素濃度Ｃi＝Ｆi(ａ1,ａ2,・・・ａN)と回帰パラメータａi（i=1,2,…,N）との関係式
さらに、被験者の皮膚画像の各ピクセルについて算出した各色素濃度Ｃi（i=1,2,…,N）の一つ又は任意の組み合わせを全ピクセルで表示させ、皮膚色素画像を形成することが好ましい。こうして得られる皮膚色素画像は、陰影の影響がなく、色素の正確な濃度分布を示すものとなる。
以下、実施例に基づいて本発明を具体的に説明する。
実施例１、比較例１
図１に示した撮像系で被験者の前腕内側の皮膚の内部反射光画像を、波長５００ｎｍ〜６００ｎｍの間を２０ｎｍ間隔で撮影して分光反射率を求めた。この場合、カメラの左右に配した光源が測定部位をできるだけむらなく照らすようにし設置し、画像に陰影ができないようにした。また、波長可変フィルタとして、ＶａｒｉＳｐｅｃ（ＣＲｉ社製）を用いた。
図１に示した撮像系で被験者の前腕内側の皮膚の内部反射光画像を、波長５００ｎｍ〜６００ｎｍの間を２０ｎｍ間隔で撮影して分光反射率を求めた。この場合、カメラの左右に配した光源が測定部位をできるだけむらなく照らすようにし設置し、画像に陰影ができないようにした。また、波長可変フィルタとして、ＶａｒｉＳｐｅｃ（ＣＲｉ社製）を用いた。
得られた分光反射率を、メラニンの吸光係数スペクトルεm(λ)、酸化ヘモグロビンの吸光係数スペクトルεoh(λ)及び還元ヘモグロビンの吸光係数スペクトルεdh(λ)を説明変数として重回帰分析することにより回帰パラメータａm、ａoh、ａdh及び切片ａm0を算出した。この場合、メラニンの吸光係数スペクトルεm(λ)、酸化ヘモグロビンの吸光係数スペクトルεoh(λ)及び還元ヘモグロビンの吸光係数スペクトルεdh(λ)は、次の文献に記載されている値を用いた。こうして得た回帰パラメータを図３(上段)に示す。
ＩＥＥＥ Ｔｒａｎｓ．Ｂｉｏｍｅｄ．Ｅｎｇ．３６（１９８９） １１４６‐１１５４
http://omlc.ogi.edu/news/jan98/skinoptics.html
http://omlc.ogi.edu/spectra/
ＳＰＩＥ Ｐｒｏｃ． ３２５２（１９９８） ７０‐８２
Ｊ．Ｉｎｖｅｓｔ Ｄｅｒｍａｔｏｌ ７７，１（１９８１） １３‐１９
http://omlc.ogi.edu/news/jan98/skinoptics.html
http://omlc.ogi.edu/spectra/
ＳＰＩＥ Ｐｒｏｃ． ３２５２（１９９８） ７０‐８２
Ｊ．Ｉｎｖｅｓｔ Ｄｅｒｍａｔｏｌ ７７，１（１９８１） １３‐１９
一方、計算ソフト「ＭＣＭＬ」（Computer Methods and Programs in Biomedicine 47:131-146 (1995)）を使用し、メラニン、酸化ヘモグロビン及び還元ヘモグロビンを含む皮膚モデルの分光反射スペクトルをモンテカルロシミュレーションモデル（Journal of Biomedical Optics 9(4):700-710 (2004)）から求め、さらに吸光度スペクトルに換算した。
この場合、皮膚モデルは、表皮厚0.006ｃｍ、真皮厚0.494ｃｍとし、表皮にメラニンが含まれ、真皮にヘモグロビンが含まれるとした。
表皮及び真皮の屈折率は1.4で一定とした。
分光波長は、500、520、540、560、580、600ｎｍの6条件とした。
メラニン濃度は、1〜10％の間で1％ごとの10条件とし、ヘモグロビン濃度は、0.2〜1.0％の間で0.2％ごとの5条件とし、酸素飽和度は、0〜100％の間で20％ごとの6条件とした。ここで、メラニン濃度は、吸光度がJournal of Photochemistry and Photobiology Vol.53 No.6 769-775で与えられたメラノソームの吸光度と等しくなるときの濃度を100%としたときの比率である。また、酸化、還元ヘモグロビン濃度は、それぞれヘマトクリット45%の血液の酸素飽和度１００％、０％のときの濃度を100%としたときの比率である。
表皮及び真皮の屈折率は1.4で一定とした。
分光波長は、500、520、540、560、580、600ｎｍの6条件とした。
メラニン濃度は、1〜10％の間で1％ごとの10条件とし、ヘモグロビン濃度は、0.2〜1.0％の間で0.2％ごとの5条件とし、酸素飽和度は、0〜100％の間で20％ごとの6条件とした。ここで、メラニン濃度は、吸光度がJournal of Photochemistry and Photobiology Vol.53 No.6 769-775で与えられたメラノソームの吸光度と等しくなるときの濃度を100%としたときの比率である。また、酸化、還元ヘモグロビン濃度は、それぞれヘマトクリット45%の血液の酸素飽和度１００％、０％のときの濃度を100%としたときの比率である。
そこで、6×10×5×6＝1800条件の反射率を計算し、10×5×6＝300条件の分光反射率を取得し、10×5×6＝300条件の色素濃度と回帰パラメータの対を取得し、2次の多項式近似を行うことにより回帰パラメータａm、ａoh及びａdhと色素濃度Ｃm、Ｃoh及びＣdhとの関係式Ｃm＝Ｆ1 m(ａm,ａoh,ａdh)、Ｃoh＝Ｆ1 oh(ａm,ａoh,ａdh)、Ｃoh＝Ｆ1 oh(ａm,ａoh,ａdh)を取得した（実施例１）。
また、同様にして、回帰パラメータ、ａm、ａth(ａth＝ａoh＋ａdh)、ａ0を用いて
関係式Ｃm＝Ｆ2 m(ａm,ａth,ａ0)、Ｃth＝Ｆ2 th(ａm,ａth,ａ0)を取得した（比較例１）。
関係式Ｃm＝Ｆ2 m(ａm,ａth,ａ0)、Ｃth＝Ｆ2 th(ａm,ａth,ａ0)を取得した（比較例１）。
図３(上段)は、被験者の肌画像から得られた切片ａ0及び回帰パラメータａm、ａoh及びａdh（実施例１、比較例１共通）である。モンテカルロシミュレーションから、切片ａ0がおよそ0.13から0.43の値をとることが予想されている。一方で、実施例１又は比較例１における画像中央付近では切片ａ0＝約2.3程度であった。このａ0の値は、前出のモンテカルロシミュレーションによる切片ａ0の値と比べて大きく乖離しているので、切片ａ0の値には肌の凹凸や照明の影響により意図しない陰影が反映されていると考えられる。
実施例１、比較例１のそれぞれで各色素濃度を求める関係式Ｆ1 m、Ｆ1 oh、Ｆ1 dh、Ｆ2 m、Ｆ2 thと、図３(上段)に示した回帰パラメータとから、実施例１、比較例１のそれぞれについてメラニン濃度と総ヘモグロビン濃度を求めた。このメラニン濃度と総ヘモグロビン濃度を画像の全ピクセルにおいて求め、表示させた。これを図５に示す。図５から、実施例１及び比較例１のいずれの方法でもメラニン濃度と総ヘモグロビン濃度が測定されているが、比較例１の画像はピンボケになっている。これは、実施例１及び比較例１の測定部位には、意図的に陰影を形成していないものの、肌状態や撮影条件等で生じた細かな陰影のため、比較例１の方法ではメラニン濃度や総ヘモグロビン濃度が正確には測定されていないことによる。
実施例２、比較例２
実施例１において画像を撮る際に、図４に示すように、カメラの左右に配置した２つある光源の内一方の光源３と被写体１との間に三日月型の遮蔽物８をおいてカメラ５で撮影される画像に意図的に明らかな陰影（三日月型の影）をつけ、この陰影付の画像を用いて実施例１、比較例１を繰り返した。なお、被験者の測定部位は実施例１と同じである。
実施例１において画像を撮る際に、図４に示すように、カメラの左右に配置した２つある光源の内一方の光源３と被写体１との間に三日月型の遮蔽物８をおいてカメラ５で撮影される画像に意図的に明らかな陰影（三日月型の影）をつけ、この陰影付の画像を用いて実施例１、比較例１を繰り返した。なお、被験者の測定部位は実施例１と同じである。
図３(下段)は、陰影付の被験者の肌画像から得られた回帰パラメータ（実施例２、比較例２共通）である。同図から、陰影付画像では、与えた陰影の影響が回帰パラメータの切片ａ0に現れており、このことから切片ａ0が色素濃度を正しく反映していないことがわかる。
図３（下段）の回帰パラメータと、実施例１及び比較例１で求めたそれぞれの関係式Ｆ1 m、Ｆ1 oh、Ｆ1 dh、Ｆ2 m、Ｆ2 thから実施例２及び比較例２のメラニン濃度及び総ヘモグロビン濃度を画像の全ピクセルについて求め、図５に示した。図５から、比較例２によれば、メラニン画像にも総ヘモグロビン画像にも遮蔽物８の陰影が映っており、色素濃度が正しく測定されていないのに対し、実施例２では陰影の影響が完全に除かれていることがわかる。
１ 被験者
２ａ 第１の偏光フィルタ
２ｂ 第２の偏光フィルタ
３ 光源
４ 波長可変フィルタ
５ デジタルカメラ
６ コントローラ
７ 演算装置
８ 遮蔽物
２ａ 第１の偏光フィルタ
２ｂ 第２の偏光フィルタ
３ 光源
４ 波長可変フィルタ
５ デジタルカメラ
６ コントローラ
７ 演算装置
８ 遮蔽物
Claims (7)
- 皮膚の色素濃度から分光反射率を算出するシミュレーションモデルに基づき、複数の色素を含む皮膚モデルの吸光度スペクトルを、各色素の吸光係数スペクトルεi(λ)（色素の種類をＮとしてi=1,2,…,N、 以下同じ）を用いて、各色素濃度Ｃi（i=1,2,…,N）の組み合わせが異なる複数の場合について求め、得られた複数の皮膚モデルの吸光度スペクトルを、各色素の吸光係数スペクトルεi(λ)（i=1,2,…,N）を説明変数として重回帰分析することにより、以下の重回帰式(1)で得られる回帰パラメータａi（i=1,2,…,N）の組を複数求め、
- 皮膚の色素濃度と分光反射率とのシミュレーションモデルとして、皮膚のモンテカルロシミュレーションを使用する請求項１記載の皮膚色素濃度測定方法。
- 被験者の皮膚の吸光度スペクトルを、被験者の皮膚の内部反射光画像から求める請求項１又は２記載の皮膚色素濃度測定方法。
- メラニン濃度Ｃm、酸化ヘモグロビン濃度Ｃoh、及び還元ヘモグロビン濃度Ｃdhの少なくとも一つの色素濃度を求める請求項１〜３のいずれかに記載の皮膚色素濃度測定方法。
- 被験者の皮膚の内部反射光画像の各ピクセルについて、請求項１〜３のいずれかに記載の方法により皮膚の少なくとも一つの色素濃度を求め、その色素濃度の分布画像を形成する皮膚色素画像の形成方法。
- メラニン濃度Ｃm、酸化ヘモグロビン濃度Ｃoh、及び還元ヘモグロビン濃度Ｃdhの少なくとも一つの色素濃度の分布画像を形成する請求項５記載の皮膚色素画像の形成方法。
- 演算装置を備えた皮膚色素画像の表示装置であって、該演算装置から出力される皮膚画像の各ピクセルが請求項５又は６記載の方法により形成され、各ピクセルが表示する色素濃度Ｃiが前記式(２)の関係式
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011194043A JP5821436B2 (ja) | 2011-09-06 | 2011-09-06 | 皮膚色素濃度測定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011194043A JP5821436B2 (ja) | 2011-09-06 | 2011-09-06 | 皮膚色素濃度測定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2013052182A JP2013052182A (ja) | 2013-03-21 |
| JP5821436B2 true JP5821436B2 (ja) | 2015-11-24 |
Family
ID=48129709
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2011194043A Active JP5821436B2 (ja) | 2011-09-06 | 2011-09-06 | 皮膚色素濃度測定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5821436B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6183030B2 (ja) * | 2013-07-25 | 2017-08-23 | 花王株式会社 | 皮膚色素濃度測定方法 |
| JP6755831B2 (ja) * | 2016-08-09 | 2020-09-16 | 花王株式会社 | 肌状態の評価方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0829333A (ja) * | 1994-07-15 | 1996-02-02 | Kubota Corp | 分光分析装置 |
| JP3706437B2 (ja) * | 1996-07-13 | 2005-10-12 | 株式会社堀場製作所 | 多成分水溶液の分析方法 |
| JP3297420B2 (ja) * | 2000-06-05 | 2002-07-02 | 財団法人塩事業センター | 成分濃度自動測定システム |
| JP2003149136A (ja) * | 2001-11-13 | 2003-05-21 | Shimadzu Corp | 光画像計測方法 |
| JP4872536B2 (ja) * | 2006-08-28 | 2012-02-08 | パナソニック電工株式会社 | 生体成分濃度測定方法 |
-
2011
- 2011-09-06 JP JP2011194043A patent/JP5821436B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2013052182A (ja) | 2013-03-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Spigulis et al. | Smartphone snapshot mapping of skin chromophores under triple-wavelength laser illumination | |
| US11751769B2 (en) | Systems and methods for assessing diabetic circulatory complications | |
| JP6755831B2 (ja) | 肌状態の評価方法 | |
| Chen et al. | In vivo real-time imaging of cutaneous hemoglobin concentration, oxygen saturation, scattering properties, melanin content, and epidermal thickness with visible spatially modulated light | |
| JP5903969B2 (ja) | 皮膚色素濃度測定方法 | |
| US11166652B2 (en) | Method and apparatus for assessing tissue vascular health | |
| Chatterjee et al. | Monte Carlo investigation of the effect of blood volume and oxygen saturation on optical path in reflectance pulse oximetry | |
| Denstedt et al. | Wavelet based feature extraction and visualization in hyperspectral tissue characterization | |
| Else et al. | Effects of skin tone on photoacoustic imaging and oximetry | |
| JP5821436B2 (ja) | 皮膚色素濃度測定方法 | |
| Urban et al. | Multimodal hyperspectral fluorescence and spatial frequency domain imaging for tissue health diagnostics of the oral cavity | |
| JP6183030B2 (ja) | 皮膚色素濃度測定方法 | |
| Huang et al. | Second derivative multispectral algorithm for quantitative assessment of cutaneous tissue oxygenation | |
| JP6131597B2 (ja) | 皮膚の評価方法及び皮膚の評価システム | |
| Kim et al. | Modeling and measuring extravascular hemoglobin: aging contusions | |
| da Cruz et al. | Optical properties of human skin phototypes and their correlation with individual angle typology | |
| Huang et al. | The image analysis of skin tissue irradiated with difference wavelengths of LED sources | |
| Zherebtsov et al. | Hyperspectral system for Imaging of skin chromophores and blood oxygenation | |
| JP6960323B2 (ja) | 肌状態の評価方法 | |
| Vasefi et al. | Toward in vivo diagnosis of skin cancer using multimode imaging dermoscopy:(II) molecular mapping of highly pigmented lesions | |
| Välisuo et al. | The colour of blood in skin: a comparison of Allen's test and photonics simulations | |
| Jakovels et al. | RGB imaging system for monitoring of skin vascular malformation's laser therapy | |
| JP7081576B2 (ja) | 皮膚組織における色素濃度の推定方法 | |
| Arimoto | Measurement of 2-D SpO2 distribution in skin tissue by multispectral imaging with depth selectivity control | |
| Livecchi et al. | Hyperspectral imaging with deep learning for quantification of tissue hemoglobin, melanin, and scattering |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140605 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150330 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150414 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150615 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150804 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150807 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150908 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150921 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 5821436 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |