JP2022063800A - 肌の評価方法 - Google Patents
肌の評価方法 Download PDFInfo
- Publication number
- JP2022063800A JP2022063800A JP2020172228A JP2020172228A JP2022063800A JP 2022063800 A JP2022063800 A JP 2022063800A JP 2020172228 A JP2020172228 A JP 2020172228A JP 2020172228 A JP2020172228 A JP 2020172228A JP 2022063800 A JP2022063800 A JP 2022063800A
- Authority
- JP
- Japan
- Prior art keywords
- skin
- feature amount
- image
- data
- oxygen saturation
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 43
- XUMBMVFBXHLACL-UHFFFAOYSA-N Melanin Chemical group O=C1C(=O)C(C2=CNC3=C(C(C(=O)C4=C32)=O)C)=C2C4=CNC2=C1C XUMBMVFBXHLACL-UHFFFAOYSA-N 0.000 claims description 154
- 238000005259 measurement Methods 0.000 claims description 107
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 93
- 229910052760 oxygen Inorganic materials 0.000 claims description 93
- 239000001301 oxygen Substances 0.000 claims description 93
- 210000004204 blood vessel Anatomy 0.000 claims description 62
- 238000000034 method Methods 0.000 claims description 60
- 230000002792 vascular Effects 0.000 claims description 58
- 230000032683 aging Effects 0.000 claims description 49
- 238000013135 deep learning Methods 0.000 claims description 28
- 238000013527 convolutional neural network Methods 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 8
- 210000003491 skin Anatomy 0.000 description 218
- 239000010410 layer Substances 0.000 description 25
- 238000012549 training Methods 0.000 description 24
- 238000004458 analytical method Methods 0.000 description 22
- 238000009826 distribution Methods 0.000 description 20
- 108010054147 Hemoglobins Proteins 0.000 description 10
- 102000001554 Hemoglobins Human genes 0.000 description 10
- 238000000691 measurement method Methods 0.000 description 9
- 230000007423 decrease Effects 0.000 description 8
- 210000001519 tissue Anatomy 0.000 description 7
- 210000002615 epidermis Anatomy 0.000 description 6
- 210000004088 microvessel Anatomy 0.000 description 6
- 210000004003 subcutaneous fat Anatomy 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000010219 correlation analysis Methods 0.000 description 4
- 238000013136 deep learning model Methods 0.000 description 3
- 210000004207 dermis Anatomy 0.000 description 3
- 230000001678 irradiating effect Effects 0.000 description 3
- 238000009828 non-uniform distribution Methods 0.000 description 3
- 230000009759 skin aging Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 206010017711 Gangrene Diseases 0.000 description 2
- 208000004210 Pressure Ulcer Diseases 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000002344 surface layer Substances 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000037303 wrinkles 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
- 238000004483 ATR-FTIR spectroscopy Methods 0.000 description 1
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 206010011985 Decubitus ulcer Diseases 0.000 description 1
- 206010014970 Ephelides Diseases 0.000 description 1
- 208000003351 Melanosis Diseases 0.000 description 1
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 108010064719 Oxyhemoglobins Proteins 0.000 description 1
- 238000010220 Pearson correlation analysis Methods 0.000 description 1
- 208000012641 Pigmentation disease Diseases 0.000 description 1
- 206010064127 Solar lentigo Diseases 0.000 description 1
- 206010052428 Wound Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 241000212749 Zesius chrysomallus Species 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 238000000335 coherent Raman spectroscopy Methods 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013502 data validation Methods 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 201000001441 melanoma Diseases 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000037311 normal skin Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 206010033675 panniculitis Diseases 0.000 description 1
- 210000004694 pigment cell Anatomy 0.000 description 1
- 230000019612 pigmentation Effects 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000007665 sagging Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 231100000075 skin burn Toxicity 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Landscapes
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
【課題】容易に生体から得られる情報を利用して、生体の皮膚内部の構造に関する推定値を出力することができる、新規な肌の評価方法を提供する。【解決手段】被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力する。【選択図】なし
Description
本発明は、肌の評価方法に関する。
近年、測定手法の発達により、非侵襲的に皮膚の内部構造を測定することが可能になっている。
例えば、非特許文献1においては、組織透明化手法によって、血管の構造が明らかにされている。
例えば、非特許文献1においては、組織透明化手法によって、血管の構造が明らかにされている。
また、非侵襲的に測定した皮膚の内部構造に基づいて、皮膚の状態を診断する方法も知られている。
例えば特許文献1においては、光音響波を用いて皮膚内部の構造を非侵襲的に測定することにより、皮膚の熱傷、壊疽、褥瘡の診断を行う方法が開示されている。
例えば特許文献1においては、光音響波を用いて皮膚内部の構造を非侵襲的に測定することにより、皮膚の熱傷、壊疽、褥瘡の診断を行う方法が開示されている。
Shibata T, Kajiya K, Sato K, Yoon J, Kang HY. 3D microvascular analysis reveals irregularly branching blood vessels in the hyperpigmented skin of solar lentigo. Pigment cell & melanoma research.2018;31(6):725 7.
ところで、コラーゲン線維の太さや多寡、皮膚内部に保持される水分量等、皮膚の内部の形態的特徴や生理学的特徴、すなわち皮膚の内部の構造が、肌の見た目(しみ、しわ、くすみ等)に影響することは、知られている。
しかし、組織透明化手法による皮膚内部の構造の解析は、皮膚を摘出して行うため、血液や酸素が失われてしまい、組織が生存している状態については、測定することができない。また、解析のためのn数を十分に確保しにくい。
すなわち、上記のような手法では、生体における皮膚の内部の構造を測定することは制限され、このような手法で得られた指標に基づいた肌の評価は、生体の状態を反映しているものであるとは必ずしも言えなかった。
しかし、組織透明化手法による皮膚内部の構造の解析は、皮膚を摘出して行うため、血液や酸素が失われてしまい、組織が生存している状態については、測定することができない。また、解析のためのn数を十分に確保しにくい。
すなわち、上記のような手法では、生体における皮膚の内部の構造を測定することは制限され、このような手法で得られた指標に基づいた肌の評価は、生体の状態を反映しているものであるとは必ずしも言えなかった。
また、特許文献1に記載されているような、非侵襲的に測定した対象者の皮膚の内部構造に基づいて、皮膚の状態(熱傷、壊疽、褥瘡等の有無)を診断する方法は知られているものの、この方法では皮膚の内部構造の測定において高価な装置が必要となるという課題があった。
上記状況に鑑み、本発明は、高価な装置を用いることなく容易に生体から得られる情報を利用して、生体の皮膚内部の構造に関する推定値を出力することができる、新規な肌の評価方法を提供することを課題とする。
上記課題を解決する本発明は、被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力する、肌の評価方法である。
また、本発明の好ましい形態では、前記内部構造を示すデータは、非侵襲的に得られたデータである。
また、本発明の好ましい形態では、前記モデルは、前記肌画像を入力データとし、前記特徴量を目的変数として与えて深層学習を行うことにより作成される畳み込みニューラルネットワークモデルである。
また、本発明の好ましい形態では、入力された前記肌画像に基づいて、前記モデルの出力側から入力側に向かって分析することにより、入力された前記肌画像の領域ごとに、前記モデルによる推定結果への寄与を推定することを含む。
また、本発明の好ましい形態では、前記肌画像の領域ごとの前記推定結果への寄与を、前記肌画像に重ねて表示することを含む。
また、本発明の好ましい形態では、前記特徴量は、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚のメラニン構造に関する特徴量の少なくとも何れかである。
また、本発明の好ましい形態では、前記血管構造に関する特徴量は、皮膚上部における微小血管構造の量に関する特徴量である。
また、本発明の好ましい形態では、前記血管構造に関する特徴量は、血管の平均厚さ、血管の分岐の数及び皮膚上部における血管の体積から選ばれる。
本発明の好ましい形態では、前記酸素飽和度に関する特徴量は、血管の酸素飽和度である。
また、本発明の好ましい形態では、前記メラニン構造に関する特徴量は、メラニンのチューブ状の紋理構造の体積である。
また、本発明の好ましい形態では、前記内部構造を示すデータは、光音響法により得られた測定データ又は該測定データに基づいて生成されたデータである。
また、本発明の好ましい形態では、皮膚の血管構造に関する特徴量に基づいて、被験者の皮膚の酸素飽和度を推定することを含む。
また、本発明の好ましい形態では、皮膚のメラニン構造に関する特徴量に基づいて、皮膚の酸素飽和度を推定することを含む。
また、本発明の好ましい形態では、前皮膚の酸素飽和度に基づいて、皮膚の血管構造に関する特徴量及び/又は皮膚のメラニン構造に関する特徴量を推定することを含む。
また、本発明の好ましい形態では、前記モデルによって出力された前記推定値に基づいて、前記被験者の肌の加齢状態の指標を推定することを含む。
また、本発明の好ましい形態では、前記加齢状態は、前記特徴量及び年齢の組を複数与えることによって作成された相関モデルに、前記推定値を入力して、推定年齢を算出することによって評価される。
また、本発明は、非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力するようにコンピュータを機能させる、肌の評価プログラムである。
本発明によれば、肌画像に基づいて、容易に皮膚の内部構造に関する特徴量の推定値を出力する肌の評価方法が提供される。特に、光音響法等の非侵襲的手法により得られた測定データを学習データとして用いる場合には、生きた皮膚の内部構造に基づく学習の結果を反映した推定が可能となる。
本発明において、「皮膚」とは、表皮層、真皮層及び皮下脂肪層を含む生体組織を指し、組織の構造に特に着目する場合に用いる。
また、本発明において、「肌」は、生体組織としては「皮膚」と同義であるが、組織の構造に特に着目しない場合に用いる。
また、本発明の「評価」は、医療診断を含まない。
また、本発明において、「肌」は、生体組織としては「皮膚」と同義であるが、組織の構造に特に着目しない場合に用いる。
また、本発明の「評価」は、医療診断を含まない。
<1>肌の評価方法
本発明の肌の評価方法は、被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力する。
具体的には、被験者の顔の肌を(採取したサンプルではなく)生体において撮影して得られた肌画像を入力データとし、同一被験者の皮膚の内部構造を示すデータから抽出した特徴量を目的変数として複数与えて深層学習を行うことにより、畳み込みニューラルネットワークモデルを作成する。
そして、前記モデルに、特徴量の推定対象である被験者の肌画像を入力して、当該被験者に関する特徴量の推定値を出力する。
本発明の肌の評価方法は、被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力する。
具体的には、被験者の顔の肌を(採取したサンプルではなく)生体において撮影して得られた肌画像を入力データとし、同一被験者の皮膚の内部構造を示すデータから抽出した特徴量を目的変数として複数与えて深層学習を行うことにより、畳み込みニューラルネットワークモデルを作成する。
そして、前記モデルに、特徴量の推定対象である被験者の肌画像を入力して、当該被験者に関する特徴量の推定値を出力する。
すなわち、本実施形態では、まずモデルの作成において、光音響法に代表される非侵襲的な方法で得られた皮膚の内部構造を示すデータに基づき、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚のメラニン構造に関する特徴量を含む、特徴量を抽出する。
そして、同一被験者の生体において肌画像を撮影し、肌画像及び特徴量の組を学習データとして与えることで、肌画像に基づいて特徴量の推定値を出力するモデルを作成する。皮膚の内部構造を示すデータ及び特徴量については後述する。
そして、同一被験者の生体において肌画像を撮影し、肌画像及び特徴量の組を学習データとして与えることで、肌画像に基づいて特徴量の推定値を出力するモデルを作成する。皮膚の内部構造を示すデータ及び特徴量については後述する。
ここで発明者らは、上記光音響法を用いた手法により、皮膚の測定領域における血管の平均厚さと当該測定領域における酸素飽和度、及び、皮膚の測定領域における酸素飽和度と当該測定領域におけるメラニンの紋理構造の体積が、それぞれ互いに相関関係を有することを明らかにしている(図9及び図13参照)。
したがって、この相関関係を利用し、相互に他の皮膚構造の特徴量を推定することで、学習データの一部を作成してもよい。例えば、皮膚の内部構造を示すデータから皮膚の血管構造に関する特徴量を抽出し、皮膚の血管構造に関する特徴量に基づいて皮膚の酸素飽和度を推定して、それぞれ学習データにおける特徴量として用いてもよい。
また同様に相関関係を利用して、モデルにより出力された特徴量の推定値から、更に他の皮膚構造の特徴量を推定してもよい。
したがって、この相関関係を利用し、相互に他の皮膚構造の特徴量を推定することで、学習データの一部を作成してもよい。例えば、皮膚の内部構造を示すデータから皮膚の血管構造に関する特徴量を抽出し、皮膚の血管構造に関する特徴量に基づいて皮膚の酸素飽和度を推定して、それぞれ学習データにおける特徴量として用いてもよい。
また同様に相関関係を利用して、モデルにより出力された特徴量の推定値から、更に他の皮膚構造の特徴量を推定してもよい。
また本発明の好ましい形態では、入力された肌画像に基づいて、モデルを出力側から入力側に向かって分析することにより、入力された肌画像の領域ごとに、モデルによる推定結果への寄与を推定して肌画像に重ねて表示する。
分析手法としては、Grad-CAM法、GuidedBackPropagation法、またこれらを組み合わせたGuided Grad-CAM法等が好ましく例示される。これらの手法は、深層学習においてモデルが注目した領域を可視化する手法として知られている。入力画像に基づき、最終畳み込み層に流入するクラス固有勾配情報を使用して、逆伝播によって分析するという方法である。Guided Grad-CAM法では、出力から入力へ後ろ向きにトレースすることにより、入力画像の画素ごとに深層学習に関心のある特徴を抽出する技術である。
これにより、皮膚の内部構造に関する特徴量が、肌画像においてどのような表現と関係するかを可視化することができる。
分析手法としては、Grad-CAM法、GuidedBackPropagation法、またこれらを組み合わせたGuided Grad-CAM法等が好ましく例示される。これらの手法は、深層学習においてモデルが注目した領域を可視化する手法として知られている。入力画像に基づき、最終畳み込み層に流入するクラス固有勾配情報を使用して、逆伝播によって分析するという方法である。Guided Grad-CAM法では、出力から入力へ後ろ向きにトレースすることにより、入力画像の画素ごとに深層学習に関心のある特徴を抽出する技術である。
これにより、皮膚の内部構造に関する特徴量が、肌画像においてどのような表現と関係するかを可視化することができる。
<2>皮膚の内部構造を示すデータ及び特徴量の取得
本実施形態で学習データの作成に用いる皮膚の内部構造を示すデータは、非侵襲的な測定方法により得られるが、本発明はこれに限られず、侵襲的に得られた測定データに基づいてモデルの学習が行われてもよい。前記データは、公知の装置及び公知の方法を適宜用いて取得することができる。公知の非侵襲的な測定方法としては、例えば、高周波電流法、ATR-FTIR法、近赤外分光法、磁気共鳴画像法(MRI)、In vivo共焦点ラマン分光法、超音波断層撮影、光コヒーレンス・トモグラフィー、In vivo共焦点レーザー走査型顕微鏡、多光子顕微鏡、SHG顕微鏡、CARS顕微鏡、光超音波顕微鏡等を選択することができる。
また、上記非侵襲的な測定方法で測定することができれば、皮膚の内部構造を示すデータの種類は特に制限されず、例えば被験者の皮膚のある測定領域における血管構造に関するデータ、血流に関するデータ、酸素飽和度に関するデータ、メラニンのミクロ空間分布に関するデータ、水分に関するデータ、皮膚の粘弾性に関するデータ、特定の物質に関するデータ等を例示することができる。
また取得した皮膚の内部構造を示すデータは、生データのまま使用してもよく、一般的な統計学的加工(平均や標準偏差の取得等)や画像化等の加工を施してもよい。
本発明の好ましい形態では、非侵襲的な測定方法として、光音響法を採用する。また本発明の好ましい形態では、皮膚の内部構造を示すデータは、光音響法により取得した画像データである。
本実施形態で学習データの作成に用いる皮膚の内部構造を示すデータは、非侵襲的な測定方法により得られるが、本発明はこれに限られず、侵襲的に得られた測定データに基づいてモデルの学習が行われてもよい。前記データは、公知の装置及び公知の方法を適宜用いて取得することができる。公知の非侵襲的な測定方法としては、例えば、高周波電流法、ATR-FTIR法、近赤外分光法、磁気共鳴画像法(MRI)、In vivo共焦点ラマン分光法、超音波断層撮影、光コヒーレンス・トモグラフィー、In vivo共焦点レーザー走査型顕微鏡、多光子顕微鏡、SHG顕微鏡、CARS顕微鏡、光超音波顕微鏡等を選択することができる。
また、上記非侵襲的な測定方法で測定することができれば、皮膚の内部構造を示すデータの種類は特に制限されず、例えば被験者の皮膚のある測定領域における血管構造に関するデータ、血流に関するデータ、酸素飽和度に関するデータ、メラニンのミクロ空間分布に関するデータ、水分に関するデータ、皮膚の粘弾性に関するデータ、特定の物質に関するデータ等を例示することができる。
また取得した皮膚の内部構造を示すデータは、生データのまま使用してもよく、一般的な統計学的加工(平均や標準偏差の取得等)や画像化等の加工を施してもよい。
本発明の好ましい形態では、非侵襲的な測定方法として、光音響法を採用する。また本発明の好ましい形態では、皮膚の内部構造を示すデータは、光音響法により取得した画像データである。
ここで、皮膚の内部構造を示すデータの取得のための測定領域としては、任意の箇所を選択することができるが、日常生活において光を浴びやすい箇所(例えば、顔や手の甲)の皮膚を選択することが好ましい。また、日常生活において光を浴びにくい箇所(例えば、上腕の内側)の皮膚を選択することも、好ましい。
前者を選択する場合には、光による影響を考慮した総合的な肌の状態を評価することができ、後者を選択する場合には、光による影響を極力排除した(生理的な)肌の状態を評価することができる。
また、測定領域としては、視認できるシミやその他色素沈着、創傷等の異常がない正常な領域を設定する。
また、測定領域のサイズとしては、好ましくは5mm×5mm~15mm×15mmで設定することができる。
前者を選択する場合には、光による影響を考慮した総合的な肌の状態を評価することができ、後者を選択する場合には、光による影響を極力排除した(生理的な)肌の状態を評価することができる。
また、測定領域としては、視認できるシミやその他色素沈着、創傷等の異常がない正常な領域を設定する。
また、測定領域のサイズとしては、好ましくは5mm×5mm~15mm×15mmで設定することができる。
本発明において、皮膚の内部構造を示すデータとして画像データを使用する場合、特徴量を抽出するための画像データは、三次元データであることが好ましい。このような三次元データとしては、光音響法等によって得られた測定画像データそのものであってもよいし、測定画像データを加工した加工画像データであってもよい。
皮膚の内部構造を示すデータとして画像データを使用する場合、当該皮膚の画像データとしては、表皮層、真皮層、皮下脂肪層のすべての層を含む画像データを用いてもよいし、これを皮膚の深さ方向に任意に分割した画像データを用いてもよい。
例えば、皮膚の深さ方向に任意に分割した画像データとしては、表皮層から皮下脂肪層までを、上部(表皮から真皮乳頭層までの領域)、中間部(真皮網状層から皮下脂肪層の境界領域)、及び下部(皮下脂肪層の上層)に分割した画像データを用いることができる。
例えば、皮膚の深さ方向に任意に分割した画像データとしては、表皮層から皮下脂肪層までを、上部(表皮から真皮乳頭層までの領域)、中間部(真皮網状層から皮下脂肪層の境界領域)、及び下部(皮下脂肪層の上層)に分割した画像データを用いることができる。
本発明のモデルの作成において用いられる、皮膚の内部構造を示すデータから得られる皮膚の特徴量は、例えば、血管構造に関する特徴量、血流に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量、水分に関する特徴量を挙げることができる。
中でも好ましい形態として、血管構造に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量を挙げることができる。
このうち、血管構造に関する特徴量を用いる場合には、皮膚の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、メラニン構造に関する特徴量を用いる場合には、表皮層の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、酸素飽和度に関する特徴量を用いる場合には、表皮層の測定データを含むデータ、又は皮下脂肪層の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
中でも好ましい形態として、血管構造に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量を挙げることができる。
このうち、血管構造に関する特徴量を用いる場合には、皮膚の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、メラニン構造に関する特徴量を用いる場合には、表皮層の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、酸素飽和度に関する特徴量を用いる場合には、表皮層の測定データを含むデータ、又は皮下脂肪層の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
本発明の評価方法では、血管構造に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量のうち、1種類を用いてもよく、2種類以上を用いてもよい。
血管構造に関する特徴量としては、微細な血管構造に関する特徴量を用いることが好ましく、特に、皮膚の上部における微細な血管構造に関する特徴量を用いることが好ましい。このような特徴量として、血管の分岐の数、血管の曲率、血管の合計量(体積)、微細血管の数、血管の密度、血管の平均厚さ等が好ましく挙げられ、中でも、血管の分岐の数、血管の合計量(体積)、微細血管の数、血管の平均厚さを用いることが好ましい。
上記血管構造に関する特徴量は、上記非侵襲的な測定方法により測定した測定領域における血管構造に関するデータから、抽出することができる。
血管構造に関する特徴量を抽出するための画像データは、例えば、575nmと650nmの二波長を含むパルス光を測定部位に照射し、575nmの光音響シグナルから、650nmの光音響シグナルを差し引いた差分値を用いて、皮膚内部の血管構造を示す三次元画像データを生成することができる。例えば、株式会社アドバンテストのHadatomoZ(Hadatomoは登録商標)を用いて、このような三次元画像データ(測定画像データ)の生成を行うことができる。そして、生成した画像を解析することで、血管構造に関する特徴量を抽出することができる。
このような測定画像データを用いて、血管構造に関する特徴量を定量化してもよい。
好ましい形態では、この測定画像データを加工することにより、血管構造に関する特徴量の抽出の精度を上げることができる。
血管構造に関する特徴量を抽出するための画像データは、例えば、575nmと650nmの二波長を含むパルス光を測定部位に照射し、575nmの光音響シグナルから、650nmの光音響シグナルを差し引いた差分値を用いて、皮膚内部の血管構造を示す三次元画像データを生成することができる。例えば、株式会社アドバンテストのHadatomoZ(Hadatomoは登録商標)を用いて、このような三次元画像データ(測定画像データ)の生成を行うことができる。そして、生成した画像を解析することで、血管構造に関する特徴量を抽出することができる。
このような測定画像データを用いて、血管構造に関する特徴量を定量化してもよい。
好ましい形態では、この測定画像データを加工することにより、血管構造に関する特徴量の抽出の精度を上げることができる。
例えば、深層学習モデルを用いて、測定画像データを水平方向画像データ(平面画像データ)に分解し、この分解した測定平面画像データから血管構造を示す要素を抽出した抽出平面画像データを生成し、この抽出平面画像データを結合することで、血管構造が抽出された抽出三次元画像データを生成する方法があげられる。
測定平面画像データから抽出平面画像データを生成するためには、深層学習モデルを使用することができる。
この場合、データキュレーターが測定平面画像データを目視評価することで血管構造を抽出した注釈画像データを作成する。
訓練に必要な数の測定平面画像データと対応する注釈画像データを訓練用画像データとして用いて、深層学習をさせる。深層学習モデルとしては、エンコーダー・デコーダーモデル構造を有するネットワークであるUnetを用いることができる。
測定平面画像データから抽出平面画像データを生成するためには、深層学習モデルを使用することができる。
この場合、データキュレーターが測定平面画像データを目視評価することで血管構造を抽出した注釈画像データを作成する。
訓練に必要な数の測定平面画像データと対応する注釈画像データを訓練用画像データとして用いて、深層学習をさせる。深層学習モデルとしては、エンコーダー・デコーダーモデル構造を有するネットワークであるUnetを用いることができる。
酸素飽和度に関する特徴量としては、測定領域(例えば、皮膚の任意の層)の酸素飽和度、高酸素飽和度の微細血管が存在する領域の量などを用いることができる。
酸素飽和度に関する特徴量を抽出するための画像データは、例えば、575nmと586nmの二波長を含むパルス光を測定部位に照射し、酸素ヘモグロビンと脱酸素ヘモグロビンとの吸収係数の差を利用して、全体に対する酸素ヘモグロビンの割合を色分けなどで可視化した測定画像を生成することで、取得することができる。例えば、株式会社アドバンテストのHadatomoZ(Hadatomoは登録商標)を用いて、このような三次元画像データ(測定画像データ)の生成を行うことができる。そして、生成した画像を解析することで、酸素飽和度に関する特徴量を抽出することができる。
このような測定画像データを用いて、酸素飽和度に関する特徴量を定量化してもよい。
また、ある位置(領域)の各波長における光音響信号の振幅比に基づいて、当該位置の酸素飽和度を算出してもよい(実施例参照)。
このような測定画像データを用いて、酸素飽和度に関する特徴量を定量化してもよい。
また、ある位置(領域)の各波長における光音響信号の振幅比に基づいて、当該位置の酸素飽和度を算出してもよい(実施例参照)。
メラニン構造に関する特徴量は、メラニン顆粒の凝集により形成される構造に関する特徴量である。ここで、本発明においてメラニン構造に関連する特徴量は、視認できるようなシミやそばかす等に関連するメラニンの凝集などの構造ではなく、通常は視認できない1mm未満四方の空間における構造を意味する。
光音響法では、通常個々のメラニン顆粒は測定することができない。
この構造は、10~数百μmの範囲の大きさの構造である。また、この構造は、好ましくは、チューブ状の紋理構造である(図10(B)参照)。
メラニン構造に関する特徴量としては、前記構造の存在量に関する特徴量があげられ、体積、密度、数があげられる。
光音響法では、通常個々のメラニン顆粒は測定することができない。
この構造は、10~数百μmの範囲の大きさの構造である。また、この構造は、好ましくは、チューブ状の紋理構造である(図10(B)参照)。
メラニン構造に関する特徴量としては、前記構造の存在量に関する特徴量があげられ、体積、密度、数があげられる。
メラニン構造に関する特徴量を抽出するための画像データは、例えば、650nmのパルス光を測定部位に照射し、得られた光音響データに基づいた測定画像を生成することで、取得することができる。例えば、株式会社アドバンテストのHadatomoZ(Hadatomoは登録商標)を用いて、このような三次元画像データ(測定画像データ)の生成を行うことができる。そして、生成した画像を解析することで、メラニン構造に関する特徴量を抽出することができる。
このような測定画像データを用いて、メラニン構造に関する特徴量を定量化してもよい。
このような測定画像データを用いて、メラニン構造に関する特徴量を定量化してもよい。
<3>評価方法(内部構造の推定)
また、発明者らは、上記光音響法を用いた手法により、皮膚の測定領域における血管の平均厚さと当該測定領域における酸素飽和度、及び、皮膚の測定領域における酸素飽和度と当該測定領域におけるメラニンの紋理構造の体積が、それぞれ互いに相関関係を有することを明らかにした(図9及び図13参照)。
したがって、本発明は、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚におけるメラニン構造に関する特徴量の少なくとも何れかを解析することにより、相互に他の皮膚構造に関する特徴量を推定することを含む肌の評価方法を提供する。
また、発明者らは、上記光音響法を用いた手法により、皮膚の測定領域における血管の平均厚さと当該測定領域における酸素飽和度、及び、皮膚の測定領域における酸素飽和度と当該測定領域におけるメラニンの紋理構造の体積が、それぞれ互いに相関関係を有することを明らかにした(図9及び図13参照)。
したがって、本発明は、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚におけるメラニン構造に関する特徴量の少なくとも何れかを解析することにより、相互に他の皮膚構造に関する特徴量を推定することを含む肌の評価方法を提供する。
具体的には、皮膚の測定領域における血管の平均厚さが高い(より厚い)ほど、当該測定領域における血管の酸素飽和度が低いという関係がある(図9参照)。この関係は、特に皮膚の上部及び下部において顕著である。
したがって、皮膚の測定領域における血管の平均厚さを測定又は推定することによって、当該測定領域における血管の酸素飽和度を推定することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って減少するから、皮膚の測定領域における血管の平均厚さの測定データに基づいて、当該測定領域における血管の酸素飽和度を推定することによって、肌(特に、肌の加齢状態)を評価することができる。加齢状態の評価については後述する。
したがって、皮膚の測定領域における血管の平均厚さを測定又は推定することによって、当該測定領域における血管の酸素飽和度を推定することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って減少するから、皮膚の測定領域における血管の平均厚さの測定データに基づいて、当該測定領域における血管の酸素飽和度を推定することによって、肌(特に、肌の加齢状態)を評価することができる。加齢状態の評価については後述する。
また上記測定領域における血管の平均厚さと当該測定領域における酸素飽和度との関係に基づけば、皮膚の測定領域における酸素飽和度を測定又は推定することによって、当該測定領域における血管の平均厚さを推定することができる。
そして上記の通り、皮膚の測定領域における血管の平均厚さは、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域における血管の平均厚さを推定することによって、肌(特に、肌の加齢状態)を評価することができる。
そして上記の通り、皮膚の測定領域における血管の平均厚さは、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域における血管の平均厚さを推定することによって、肌(特に、肌の加齢状態)を評価することができる。
また、皮膚の測定領域における酸素飽和度が高いほど、当該測定領域に存在するメラニンの紋理構造の体積が少ないという関係がある(図13参照)。
したがって、皮膚の測定領域における酸素飽和度を測定又は推定することによって、当該測定領域において観察されるメラニンの紋理構造の体積を推定することができる。
そして上記の通り、皮膚の測定領域において観察されるメラニンの紋理構造の体積は、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域において観察されるメラニンの紋理構造の体積を推定することによって、肌(特に、肌の加齢状態)を評価することができる。
したがって、皮膚の測定領域における酸素飽和度を測定又は推定することによって、当該測定領域において観察されるメラニンの紋理構造の体積を推定することができる。
そして上記の通り、皮膚の測定領域において観察されるメラニンの紋理構造の体積は、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域において観察されるメラニンの紋理構造の体積を推定することによって、肌(特に、肌の加齢状態)を評価することができる。
また上記測定領域における酸素飽和度と、当該測定領域に存在するメラニンの紋理構造の体積との関係に基づけば、皮膚の測定領域において観察されるメラニンの紋理構造の体積を測定することによって、当該測定領域における酸素飽和度を推定することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って現象するから、皮膚の測定領域において観察されるメラニンの紋理構造の体積の測定データに基づいて、当該測定領域における酸素飽和度を推定することによって、肌(特に、肌の加齢状態)を評価することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って現象するから、皮膚の測定領域において観察されるメラニンの紋理構造の体積の測定データに基づいて、当該測定領域における酸素飽和度を推定することによって、肌(特に、肌の加齢状態)を評価することができる。
このように、皮膚の内部構造に関する特徴量を相互に推定することにより、本発明において肌画像から特徴量を推定するモデルの学習データの作成において、ある特徴量から他の特徴量に関する学習データを作成することができる。また、肌画像からある特徴量をモデルにより推定した後、他の特徴量を更に推定することができる。各特徴量やその抽出方法の好ましい形態は、<2>で述べた通りである。
<4>加齢状態の評価
本発明の肌の評価方法は、モデルによって出力された、被験者の皮膚の内部構造に関する特徴量の推定値に基づいて、被験者の肌の加齢状態の評価を行う工程を更に含んでいてもよい。
また、本発明の好ましい形態では、前記加齢状態は、前記特徴量及び年齢の組を複数与えることによって作成された相関モデルに、前記推定値を入力して、推定年齢を算出することによって評価される。
本発明によれば、肌画像に基づく推定値から肌の加齢状態を評価することができるので、簡便かつ非侵襲的に被験者の肌の加齢状態を評価することができる。
本発明の肌の評価方法は、モデルによって出力された、被験者の皮膚の内部構造に関する特徴量の推定値に基づいて、被験者の肌の加齢状態の評価を行う工程を更に含んでいてもよい。
また、本発明の好ましい形態では、前記加齢状態は、前記特徴量及び年齢の組を複数与えることによって作成された相関モデルに、前記推定値を入力して、推定年齢を算出することによって評価される。
本発明によれば、肌画像に基づく推定値から肌の加齢状態を評価することができるので、簡便かつ非侵襲的に被験者の肌の加齢状態を評価することができる。
本発明において、「肌の加齢状態」とは、第三者から視認可能な肌の加齢形質(例えば、しみ、しわ、くすみ等)の程度(例えば、大きさ、数、深さ等)、将来的に表出する可能性がある加齢形質の発生の予測、被験者の加齢印象(見た目年齢)等を例示することができる。
実施例に示す通り、前述した各特徴量は、年齢と関係を有することを、本発明者らは明らかにしている。したがって当該関係を利用して、前記特徴量に基づいて肌の加齢状態を評価することができる。
本実施形態では、学習済みの畳み込みニューラルネットワークモデルを利用して肌画像に基づき推定された特徴量から、加齢状態を評価する。モデルの学習に用いる学習データの作成において、皮膚の内部構造を示すデータから特徴量を抽出する方法は特に制限されず、非侵襲的な測定方法により得られた、皮膚の内部構造を示すデータに一般的な統計学的処理(例えば、平均、標準偏差等の算出)を施すことで特徴量を抽出する形態や、画像化等の加工を施した後に、当該画像を解析して特徴量を抽出する形態としてもよい。
本発明では、肌の加齢状態の評価においては、あらかじめ特徴量と年齢の相関関係を示す相関モデルを作成しておくことが、簡便性の問題から好ましい。これにより、肌画像から推定された特徴量を前記相関モデルに入力することで、推定年齢を算出することができる。なお年齢に限らず、特徴量と相関関係を有する任意の加齢状態に関する任意の指標を用いることができる。
以下に相関関係の解析と相関モデルの作成について説明する。
以下に相関関係の解析と相関モデルの作成について説明する。
本発明における加齢状態の評価に用いる特徴量(評価指標)は、以下のように決定することができる。
まず、非侵襲的な測定方法を用いて、年代の異なる複数被験者の皮膚の測定領域を測定して、皮膚の内部構造のデータを得る。皮膚の内部構造のデータは、好ましくは画像データであるか、又は得られたデータを画像化したものである。皮膚の内部構造のデータを画像データとする場合、当該画像は、皮膚断面画像等の二次元画像であっても、三次元画像であってもよい。
まず、非侵襲的な測定方法を用いて、年代の異なる複数被験者の皮膚の測定領域を測定して、皮膚の内部構造のデータを得る。皮膚の内部構造のデータは、好ましくは画像データであるか、又は得られたデータを画像化したものである。皮膚の内部構造のデータを画像データとする場合、当該画像は、皮膚断面画像等の二次元画像であっても、三次元画像であってもよい。
非侵襲的な測定方法としては、上記した各種装置及び方法を使用することができるが、光音響法を採用することが好ましい。
光音響法を採用することで、従来の方法では測定できなかった生体深部の構造を、より正確に測定することができる。
例えば、皮膚の内部構造の画像データを得る場合、年代の異なる複数の被験者の皮膚の測定対象部位にパルス光を照射して、発生した超音波を測定し、測定された皮膚内部の超音波のデータに基づいて、各年代の被験者についてそれぞれ皮膚内部の構造の画像データを生成する。
光音響法を採用することで、従来の方法では測定できなかった生体深部の構造を、より正確に測定することができる。
例えば、皮膚の内部構造の画像データを得る場合、年代の異なる複数の被験者の皮膚の測定対象部位にパルス光を照射して、発生した超音波を測定し、測定された皮膚内部の超音波のデータに基づいて、各年代の被験者についてそれぞれ皮膚内部の構造の画像データを生成する。
次に、得られた測定データに含まれるパラメータを抽出し、実年齢又は見た目年齢と相関のあるパラメータを探索する。
測定データからのパラメータの抽出方法は特に制限されず、公知の解析手法を用いることができる。例えば、測定データとして画像データを使用する場合、年代ごとの画像を目視で比較し、変化の見られた皮膚のパラメータを、抽出する方法を採用することができる。また、汎用の画像処理ソフトを用いた解析によって、色分布、輝度分布、空間周波数強度等を含む特徴を抽出することによって、パラメータを抽出することもできる。
測定データからのパラメータの抽出方法は特に制限されず、公知の解析手法を用いることができる。例えば、測定データとして画像データを使用する場合、年代ごとの画像を目視で比較し、変化の見られた皮膚のパラメータを、抽出する方法を採用することができる。また、汎用の画像処理ソフトを用いた解析によって、色分布、輝度分布、空間周波数強度等を含む特徴を抽出することによって、パラメータを抽出することもできる。
ここで、加齢状態の評価に用いる特徴量(評価指標)の決定においても、非侵襲的な方法で測定される皮膚の内部構造を示すデータから抽出されるパラメータ(特徴量)に替えて、モデルによって肌画像から推定された特徴量を用いてもよい。
次に、抽出したパラメータと、見た目年齢又は実年齢等の加齢状態を示す指標との相関を解析する。ここでの解析においても、モデルによって肌画像から推定された特徴量を用いてもよいし、非侵襲的な方法で測定される皮膚の内部構造を示すデータから得た特徴量を用いてもよい。なお、見た目年齢として、年齢を年代やスコア等に適宜変換した値を用いてもよい。
解析手法としては、回帰分析など、相関関係の解析に一般的に用いられる方法を採用することができる。
解析手法としては、回帰分析など、相関関係の解析に一般的に用いられる方法を採用することができる。
また、測定データとして画像データを使用する場合、非侵襲的な測定方法により得られた皮膚の内部構造を示す画像データに、見た目年齢又は実年齢等、加齢の度合いを示すデータのラベルを付したものを訓練データとした機械学習により、相関関係を解析することもできる。
解析の結果、年齢と相関があるとされた特徴量を、肌の評価指標として決定し、前記特徴量と年齢との相関関係を説明する相関モデルを作成する。相関モデルとしては、重回帰式、判別式、ニューラルネットワーク等があげられる。
発明者らの解析によって、具体的には、実年齢が高いほど、被験者の測定対象部位における血管の平均厚さが厚い(より血管が厚くなっている)という関係が見出された(図6参照)。
また、実年齢が高いほど、血管の体積が小さいこと、実年齢が高いほど皮膚の上部においては血管の分岐の数が少ないという関係が見いだされた。
また、実年齢が高いほど、測定領域における酸素飽和度が低く、特に皮膚の上部と下部においてこの関係が強い傾向が見出された(図8参照)。
また、実年齢が高いほど、血管の体積が小さいこと、実年齢が高いほど皮膚の上部においては血管の分岐の数が少ないという関係が見いだされた。
また、実年齢が高いほど、測定領域における酸素飽和度が低く、特に皮膚の上部と下部においてこの関係が強い傾向が見出された(図8参照)。
また、実年齢が高いほど、単位皮膚領域あたりに観察されるメラニンの紋理構造が多い(体積が大きい)という関係が見出された(図12参照)。
したがって、これらの特徴量は、本発明における加齢状態の評価方法に用いる特徴量(評価指標)として特に好ましく用いることができる。
本発明においては、上記特徴量のほかにも、肌の加齢状態の評価の指標として知られている項目を併せて測定し、これを加味して肌の加齢状態を評価してもよい。
例えば、公知の方法により、肌の水分量や粘弾性、加齢形質の多少等を非侵襲的に測定し、これらの測定結果と、上記特徴量の、光音響法による測定結果の少なくとも何れかの結果とを組み合わせ、被験者の肌の加齢状態を評価する形態とすることができる。
例えば、公知の方法により、肌の水分量や粘弾性、加齢形質の多少等を非侵襲的に測定し、これらの測定結果と、上記特徴量の、光音響法による測定結果の少なくとも何れかの結果とを組み合わせ、被験者の肌の加齢状態を評価する形態とすることができる。
特徴量から、肌を評価する方法としては、例えば、前述した相関分析に基づいて、あらかじめ特徴量と加齢状態の推定値(見た目年齢やスコア)との相関関係を示す式やモデルをコンピュータに記憶しておき、当該式やモデルに、特徴量を入力することで、加齢状態の推定値を出力する方法が好ましく用いられる。本発明においては、ここで入力する特徴量として、肌画像に基づいて推定された特徴量を用いることで、特徴量の推定を介して肌画像から加齢状態を推定することができる。
被験者(20代~60代の日本人女性、80名)の、シミや赤斑等が肉眼で観察される部分を除いた正常な頬部位(9mm×9mm)について、以下の方法により、光音響測定画像データを得て、様々な特徴量を抽出した。
更に同一被験者の顔画像から、光音響法による測定領域を含む領域の画像(1200px×1200px)を抽出して、更に頬部分の拡大皮膚画像(640px×480px)を30倍のビデオ顕微鏡を用いて取得し、肌画像と特徴量とを学習データとして深層学習を行った。
また、抽出した特徴量と、年齢との相関を解析することにより、加齢状態に関係のある特徴量を解析した。
更に同一被験者の顔画像から、光音響法による測定領域を含む領域の画像(1200px×1200px)を抽出して、更に頬部分の拡大皮膚画像(640px×480px)を30倍のビデオ顕微鏡を用いて取得し、肌画像と特徴量とを学習データとして深層学習を行った。
また、抽出した特徴量と、年齢との相関を解析することにより、加齢状態に関係のある特徴量を解析した。
<1>光音響測定画像の取得
光音響測定には、株式会社アドバンテストのHadatomoZ光音響測定装置を用いた。この装置は、二波長光音響イメージングおよび超音波イメージングが可能である。測定は、30μmの走査ステップで行った。
光音響測定には、株式会社アドバンテストのHadatomoZ光音響測定装置を用いた。この装置は、二波長光音響イメージングおよび超音波イメージングが可能である。測定は、30μmの走査ステップで行った。
測定には、2種類の二波長レーザー、(i)575nmおよび586nm、ならびに(ii)575nmおよび650nmを使用した。
血管のイメージングのためには、575nmの光音響信号から650nmの光音響信号を差し引いた値を使用した。
酸素飽和度分布のイメージングのためには、575nmの光音響信号および586nmの光音響信号から得られるオキシヘモグロビンとデオキシヘモグロビンとの間の吸収係数の差を使用した。
メラニン分布のイメージングのためには、650nmの光音響信号を使用した。
血管のイメージングのためには、575nmの光音響信号から650nmの光音響信号を差し引いた値を使用した。
酸素飽和度分布のイメージングのためには、575nmの光音響信号および586nmの光音響信号から得られるオキシヘモグロビンとデオキシヘモグロビンとの間の吸収係数の差を使用した。
メラニン分布のイメージングのためには、650nmの光音響信号を使用した。
<2>血管構造及びメラニン構造の解析のための画像生成
血管構造やメラニン構造に関する特徴量を抽出するために、MATLAB(登録商標)2020a(MathWorks、USA)、Python3.7、およびTensorflowを使用して、光音響測定により得られた画像データに基づいて、解析画像データを生成した。
光音響測定で得た三次元画像データをz軸上に積層した画像から水平方向(x‐y平面)画像に再構成した。再構成した画像から血管構造およびメラニン構造を抽出するために、エンコーダー・デコーダーモデル構造を有するネットワークであるUnetを、深層学習セグメント化方法として使用した。
血管構造およびメラニン構造が観察された測定領域の画像を、z軸の任意の複数の位置におけるx-y平面画像を切り出し、データキュレーターによって血管構造またはメラニン凝集構造と視覚的に見なされる領域についての注釈画像データ(訓練データ)として抽出画像を作成した。これらの抽出画像は、深層学習における出力データ、およびそのクラスラベルとした。これらの抽出画像は、血管構造およびメラニン構造について、それぞれ90例および30例であった。次に、元のx-y平面画像を入力データとして使用して血管およびメラニン構造を抽出するために、5回反復および2×2プーリングを有するネットワークを構築した。目的関数は、以下の式(1)によって計算されるダイス関数であると仮定した。
ここで、yは元のx-y平面画像の画素ベクトルを表し、y~は、注釈画像データの画素ベクトルを表す。
血管構造やメラニン構造に関する特徴量を抽出するために、MATLAB(登録商標)2020a(MathWorks、USA)、Python3.7、およびTensorflowを使用して、光音響測定により得られた画像データに基づいて、解析画像データを生成した。
光音響測定で得た三次元画像データをz軸上に積層した画像から水平方向(x‐y平面)画像に再構成した。再構成した画像から血管構造およびメラニン構造を抽出するために、エンコーダー・デコーダーモデル構造を有するネットワークであるUnetを、深層学習セグメント化方法として使用した。
血管構造およびメラニン構造が観察された測定領域の画像を、z軸の任意の複数の位置におけるx-y平面画像を切り出し、データキュレーターによって血管構造またはメラニン凝集構造と視覚的に見なされる領域についての注釈画像データ(訓練データ)として抽出画像を作成した。これらの抽出画像は、深層学習における出力データ、およびそのクラスラベルとした。これらの抽出画像は、血管構造およびメラニン構造について、それぞれ90例および30例であった。次に、元のx-y平面画像を入力データとして使用して血管およびメラニン構造を抽出するために、5回反復および2×2プーリングを有するネットワークを構築した。目的関数は、以下の式(1)によって計算されるダイス関数であると仮定した。
前記構築したネットワークを用いて、被験者から取得した光音響データから、血管構造およびメラニン構造が抽出されたx-y平面画像データを抽出し、この抽出した平面画像データを深さ方向で再び統合して、三次元データを作成した。
作成した三次元データに基づいて、長さ、厚さ、分岐の数、体積等のパラメータを測定した。皮膚の深さ方向に3等分した上部、中間部および下部のそれぞれについて、種々のパラメータを測定した。
<3>酸素飽和度の計算
酸素飽和度の空間分布SO2(r)は式(2)に示すように、酸化ヘモグロビンHbO2と脱酸素ヘモグロビンHbとの比率に基づいて算出することができる。
ここで、rは位置情報、εHbは脱酸素ヘモグロビンのモル吸収係数、εΔHbは酸化ヘモグロビンと脱酸素ヘモグロビンのモル吸収係数の差分値、μa(r)は光音響信号である。λ1とλ2はレーザー波長、波長λ1は575nm、波長λ2は586nmである。波長λ1およびλ2に対応する脱酸素ヘモグロビンおよび酸化ヘモグロビンのモル吸収係数が既知の値であると仮定すると、酸素飽和度の空間分布SO2(r)は、各波長における光音響信号の振幅比に基づいて算出することができる。皮膚の深さ方向に上部、中間部及び下部とわけ、そのそれぞれについて、酸素飽和度に関する種々のパラメータを測定した。
酸素飽和度の空間分布SO2(r)は式(2)に示すように、酸化ヘモグロビンHbO2と脱酸素ヘモグロビンHbとの比率に基づいて算出することができる。
<4>肌画像に基づく特徴量推定モデルの作成
顔画像の取得にはVISIA(Canfield、USA)を用いた。そして取得した顔画像から、光音響法による測定領域を含む領域の画像(1200px×1200px)を抽出し、更にそこから測定領域である頬部分の拡大皮膚画像(640px×480px)を、30倍のビデオ顕微鏡を用いて取得した。
顔画像の取得にはVISIA(Canfield、USA)を用いた。そして取得した顔画像から、光音響法による測定領域を含む領域の画像(1200px×1200px)を抽出し、更にそこから測定領域である頬部分の拡大皮膚画像(640px×480px)を、30倍のビデオ顕微鏡を用いて取得した。
深層学習は、MATLAB(登録商標)2020a及び深層学習ツールボックスを用いて行った。まず20代~60代の日本人女性80名の被験者から、特徴量の測定値によって、上位20名の上位群及び下位20名の下位群を抽出した。そして上位群及び下位群の合計40名の被験者の顔画像から抽出された、特徴量の測定領域を含む領域の画像及び拡大皮膚画像を入力データとし、測定値の上位群であるか下位群であるかを判別するように、畳み込みニューラルネットワークモデルに学習させた。ここで特徴量としては、光音響法により得られる、血管構造の体積、酸素飽和度及びメラニンの管状パターン構造の体積を採用し、そのそれぞれについて、光音響法による測定領域を含む領域の画像及び拡大皮膚画像から、上位群、下位群を判別するニューラルネットワークモデルを作成した。
更に本発明者らは、伝達学習のための訓練済みモデルであるGoogLeNetを用いて、サンプルサイズが小さいことによる影響を軽減した。
更に、反復的なトレーニングによって、予測が精度(95%以上)にまで向上した。また学習に用いていない検証データを用いた検討においても75%程度の精度を示すことが検証された。メラニン構造と皮膚画像との関係について深層学習を行った結果を代表例として図1に示す。
更に本発明者らは、伝達学習のための訓練済みモデルであるGoogLeNetを用いて、サンプルサイズが小さいことによる影響を軽減した。
更に、反復的なトレーニングによって、予測が精度(95%以上)にまで向上した。また学習に用いていない検証データを用いた検討においても75%程度の精度を示すことが検証された。メラニン構造と皮膚画像との関係について深層学習を行った結果を代表例として図1に示す。
<5>肌画像と皮膚の内部構造に関する特徴量の関係解析
本発明者らは、深層学習においてどのような画像特徴が識別に影響するかを可視化するため、Grad-CAM法を用いて可視化を行った。
本発明者らは、深層学習においてどのような画像特徴が識別に影響するかを可視化するため、Grad-CAM法を用いて可視化を行った。
血管構造と顔面及び皮膚の特徴との関係を図2(A)に示す。血管構造の平均厚さがより薄い被験者では、Grad-CAM画像に示されるヒートマップは均一に拡大する傾向があったが(図2(A)の上段画像)、血管構造の平均厚さがより厚い被験者では、局所的なヒートマップとなった(図2(A)の下段画像)。
Guided Grad-CAM画像で示されるように、血管構造の平均厚さがより厚い被験者において、雲様特徴の不均一な分布が観察された一方、微細な血管構造を有する被験者ではこのような特徴は観察されなかった。
更に、血管構造の平均厚さがより厚い例では、Guided Grad-CAM画像で示された領域と一致するカラーシェーディングとムラの存在が顔画像(入力画像)において感じられることが見出された。本研究で認められた血管構造の厚さは、皮膚のたるみとして知覚される顔面上の色むらにつながる可能性が示唆される。血管構造の平均厚さがより厚い被験者のGuided Grad-CAM画像中の矢印は、入力皮膚画像中のカラーシェーディング及び不均一性の存在と一致する雲様特徴の不均一な分布を示す。
Guided Grad-CAM画像で示されるように、血管構造の平均厚さがより厚い被験者において、雲様特徴の不均一な分布が観察された一方、微細な血管構造を有する被験者ではこのような特徴は観察されなかった。
更に、血管構造の平均厚さがより厚い例では、Guided Grad-CAM画像で示された領域と一致するカラーシェーディングとムラの存在が顔画像(入力画像)において感じられることが見出された。本研究で認められた血管構造の厚さは、皮膚のたるみとして知覚される顔面上の色むらにつながる可能性が示唆される。血管構造の平均厚さがより厚い被験者のGuided Grad-CAM画像中の矢印は、入力皮膚画像中のカラーシェーディング及び不均一性の存在と一致する雲様特徴の不均一な分布を示す。
次に、酸素飽和度と外観の関係を調べた(図2(B))。酸素飽和度の高い被験者のGrad-CAM画像に示されたヒートマップは、広く均一に広がる傾向があった(図2(B)の上段画像)。これに対し、酸素飽和度の低い被験者のGrad-CAM画像では、ヒートマップが分離する傾向があった(図2(B)の下段画像)。
更にGuided Grad-CAM画像に示されるように、酸素飽和度が低い被験者では雲の形で特徴の分布が不均一に広がるのに対し、酸素飽和度の高い被験者では均一に広がる。このことから、本研究で得られた酸素飽和度から、顔の比較的巨視的な色むらが導かれることが示唆された。Guided Grad-CAM画像中の矢印は、酸素飽和度が低い被験者において雲状に広がった特徴の不均一な分布を示し、これは入力画像のカラーシェーディング及び不均一性の存在と一致する。
更にGuided Grad-CAM画像に示されるように、酸素飽和度が低い被験者では雲の形で特徴の分布が不均一に広がるのに対し、酸素飽和度の高い被験者では均一に広がる。このことから、本研究で得られた酸素飽和度から、顔の比較的巨視的な色むらが導かれることが示唆された。Guided Grad-CAM画像中の矢印は、酸素飽和度が低い被験者において雲状に広がった特徴の不均一な分布を示し、これは入力画像のカラーシェーディング及び不均一性の存在と一致する。
また、メラニンの管状パターン構造と顔の皮膚画像との間の関係を調べた。本発明者らは、大量のメラニン構造を有する被験者のGuided Grad-CAM画像に示されるように、特徴が雲状に広く分布していることを見出した(図3の下段画像)。また入力画像の特徴がある領域では色むらが観察され、ここで見つかったメラニンの管状パターン構造が、スキンカメラで撮影された皮膚画像の色むらの視点(図中の矢印)として表現されていることが示唆される。このように本発明者らは、メラニンの管状パターン構造の存在が、皮膚上の不均一な色点の生成に関与する可能性を実証した。
<6>統計解析
上記で得られた種々のパラメータ、被験者の年齢について、JMPソフトウェアバージョン14.0(米国SAS Institute)を用いて解析した。各データ間の相関は、ピアソンの相関分析によって評価した。
上記で得られた種々のパラメータ、被験者の年齢について、JMPソフトウェアバージョン14.0(米国SAS Institute)を用いて解析した。各データ間の相関は、ピアソンの相関分析によって評価した。
<7>血管構造の解析
図4(A)は、575nmの光音響信号から得られた画像と650nmの超音波信号から得られた画像をマージした画像である。このマージ画像から、光音響測定によると、皮下の上部から下部の深さで血管構造を識別することができることがわかった。
図4(A)は、575nmの光音響信号から得られた画像と650nmの超音波信号から得られた画像をマージした画像である。このマージ画像から、光音響測定によると、皮下の上部から下部の深さで血管構造を識別することができることがわかった。
図4(B)(C)は、575nmの光音響信号から650nmの光音響信号を引いた差分値に基づいて抽出した、血管構造を示す代表的な測定画像である。(B)が20代被験者の画像であり、(C)が60代被験者の画像である。
特に、上面図を見ると明らかなとおり、図4より、20代被験者に見られた微小な血管構造が存在する領域(図4(B)中、白い三角形で示した血管)が、60代被験者では見られない(図4(C)中白い矢印で示した部分)ことがわかる。
特に、上面図を見ると明らかなとおり、図4より、20代被験者に見られた微小な血管構造が存在する領域(図4(B)中、白い三角形で示した血管)が、60代被験者では見られない(図4(C)中白い矢印で示した部分)ことがわかる。
図5は、前述の深層学習を用いて構築したネットワークにより、光音響測定により得られた血管構造を含むx-y平面の測定画像データと、これをネットワークに入力して得られた、血管構造の予測画像データの代表例である。測定画像データからデータキュレーターにより視覚的に抽出された血管構造の注釈画像データを併せて示す。(A)は、深層学習において訓練データとして用いられた測定画像データから得られた予測画像であり、(B)は、深層学習において訓練データとして用いられていない測定画像データから得られた予測画像である。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
(C)は、予測画像データ(x-y平面)をz軸方向に結合することで生成した三次元データの代表例である。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、いくつかのパラメータと年齢との間に相関関係が確認された。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、いくつかのパラメータと年齢との間に相関関係が確認された。
図6に示される通り、皮膚の各層の何れにおいても、血管構造の平均厚さ(Average thickness)は、年齢とともに増加することが観察された。
また、年齢の増加とともに、上部における血管構造の体積が減少し、下部における血管構造の体積が増加することが観察された(データは図示しない)。
また、年齢の増加とともに分岐の数は、減少することが観察された(データは図示しない)。
すなわち、加齢に伴い、皮膚上部においては、微小血管が減少し、血管構造全体の体積としても減少するのに対し、皮膚下部においては、微小血管は減少するものの血管構造全体の体積は減少しないことがわかり、加齢に伴う血管構造の変化は、皮膚の層(上部、中間部及び下部)によって異なることが示唆された。
また、年齢の増加とともに、上部における血管構造の体積が減少し、下部における血管構造の体積が増加することが観察された(データは図示しない)。
また、年齢の増加とともに分岐の数は、減少することが観察された(データは図示しない)。
すなわち、加齢に伴い、皮膚上部においては、微小血管が減少し、血管構造全体の体積としても減少するのに対し、皮膚下部においては、微小血管は減少するものの血管構造全体の体積は減少しないことがわかり、加齢に伴う血管構造の変化は、皮膚の層(上部、中間部及び下部)によって異なることが示唆された。
以上より、微小血管の量を示す特徴量については、加齢状態の評価のための指標として用いることができることが明らかとなった。
そして、その特徴量としては、血管構造の平均厚さ、分岐の数を用いることができることが明らかとなった。
また、皮膚上部における血管構造の体積についても、加齢状態の評価のための指標として用いることができることが明らかとなった。
また、皮膚上部における血管構造の平均厚さ、分岐の数、血管構造の体積の任意の組み合わせを用いて、加齢状態の評価のための指標として用いることができることが明らかとなった。
そして、その特徴量としては、血管構造の平均厚さ、分岐の数を用いることができることが明らかとなった。
また、皮膚上部における血管構造の体積についても、加齢状態の評価のための指標として用いることができることが明らかとなった。
また、皮膚上部における血管構造の平均厚さ、分岐の数、血管構造の体積の任意の組み合わせを用いて、加齢状態の評価のための指標として用いることができることが明らかとなった。
<8>酸素飽和度の解析
図7に、20代被験者及び60代被験者の光音響測定画像データを示す。光音響イメージングにより、酸素飽和度の空間分布が可視化できることが明らかとなった。
図7に、20代被験者及び60代被験者の光音響測定画像データを示す。光音響イメージングにより、酸素飽和度の空間分布が可視化できることが明らかとなった。
図7中、(A)が20代被験者の血管における酸素飽和度の空間分布を示す測定画像データであり、(B)が60代被験者の血管における酸素飽和度の空間分布を示す測定画像データである。
図7において、(A)と(B)とを比較すると、60代被験者の測定領域における酸素飽和度は、20代被験者の測定領域における血管の酸素飽和度に比べて、全体的に低くなっていた。また、20代被験者の測定画像データにおいては、微細血管が高い酸素飽和度で存在する領域((A)中、三角形で示す)が観察されるが、60代被験者の測定画像データにおいては、当該領域が存在しなかった((B)中、矢印で示す)。
図7において、(A)と(B)とを比較すると、60代被験者の測定領域における酸素飽和度は、20代被験者の測定領域における血管の酸素飽和度に比べて、全体的に低くなっていた。また、20代被験者の測定画像データにおいては、微細血管が高い酸素飽和度で存在する領域((A)中、三角形で示す)が観察されるが、60代被験者の測定画像データにおいては、当該領域が存在しなかった((B)中、矢印で示す)。
統計解析に供するために、前述した方法により酸素飽和度を算出し、先に測定した血管構造に関する特徴量、および年齢との相関関係を解析した。
図8に示す通り、加齢に伴い、皮膚の各層における酸素飽和度(Oxygen satulation)は減少する傾向があることが明らかになった。
すなわち、酸素飽和度は、肌の加齢状態の評価指標として使用できることが明らかになった。
すなわち、酸素飽和度は、肌の加齢状態の評価指標として使用できることが明らかになった。
また、図9に示す通り、皮膚の上部と下部においては、血管の平均厚さと酸素飽和度の間に負の相関があることが明らかになった。
したがって、血管の平均厚さから酸素飽和度を推定したり、反対に酸素飽和度の測定によって、血管の平均厚さを推定したりすることができる。
またこの解析により、微小な血管構造が加齢に伴って消失することで、酸素飽和度の低下を誘発することが示唆された。また、この現象は、少なくとも皮膚の上部及び下部で起こることが分かった。
またこの解析により、微小な血管構造が加齢に伴って消失することで、酸素飽和度の低下を誘発することが示唆された。また、この現象は、少なくとも皮膚の上部及び下部で起こることが分かった。
以上より、測定領域における酸素飽和度の分布を示す特徴量については、加齢状態の評価のための指標として用いることができることが明らかとなった。
そして、その特徴量としては、酸素飽和度、高酸素飽和度の微細血管領域の体積を用いることができることが明らかとなった。
特に、上部、及び下部における酸素飽和度は、加齢状態と関係のある血管の厚さの平均とも関係を有していることから、これらの層における酸素飽和度と血管構造に関する特徴量を組み合わせて評価に用いることで、より精度の高い評価ができることが示唆された。
そして、その特徴量としては、酸素飽和度、高酸素飽和度の微細血管領域の体積を用いることができることが明らかとなった。
特に、上部、及び下部における酸素飽和度は、加齢状態と関係のある血管の厚さの平均とも関係を有していることから、これらの層における酸素飽和度と血管構造に関する特徴量を組み合わせて評価に用いることで、より精度の高い評価ができることが示唆された。
<9>メラニン構造の解析
図10に、波長650nmのパルス光を用いた光音響法により取得した測定画像データを示す。
図10に、波長650nmのパルス光を用いた光音響法により取得した測定画像データを示す。
図10中、(A)が20代被験者のメラニン構造の分布を示す測定画像データであり、(B)が60代被験者のメラニン構造の分布を示す測定画像データである。
(A)に示す通り、20代被験者の肌においては、メラニン構造がそれほど豊富ではなかった。
これに対し、(B)に示す通り、60代被験者の肌においては、トロイダル空洞を伴う、チューブ状の紋理構造(矢印で示す)が観察された。
(A)に示す通り、20代被験者の肌においては、メラニン構造がそれほど豊富ではなかった。
これに対し、(B)に示す通り、60代被験者の肌においては、トロイダル空洞を伴う、チューブ状の紋理構造(矢印で示す)が観察された。
図11は、前述の深層学習を用いて構築したネットワークにより、光音響測定により得られたメラニン構造を含むx-y平面の測定画像データと、これをネットワークに入力して得られた、メラニン構造の予測画像データの代表例である。測定画像データからデータキュレーターにより視覚的に抽出されたメラニン構造の注釈画像データを併せて示す。(A)は、深層学習において訓練データとして用いられた測定画像データから得られた予測画像データであり、(B)は、深層学習において訓練データとして用いられていない測定画像データから得られた予測画像データである。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
(C)は、予測画像データ(x-y平面)をz軸方向に結合することで生成した三次元データの代表例である。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、メラニンの管状パターン構造の体積と年齢との間に相関関係が確認された。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、メラニンの管状パターン構造の体積と年齢との間に相関関係が確認された。
図12に示すとおり、加齢に伴って、皮膚の表層部に存在するメラニンのチューブ状の紋理構造の体積(volume of melanin structure)が増大することが明らかになった。
従来、正常皮膚におけるメラニンの三次元分布は研究されていない。また、本研究で実証された10~数百μmのスケールのメラニン構造についても報告はない。
この新たに発見された特徴的なメラニン構造に関する特徴量は、加齢状態の評価に有用であることが明らかになった。
従来、正常皮膚におけるメラニンの三次元分布は研究されていない。また、本研究で実証された10~数百μmのスケールのメラニン構造についても報告はない。
この新たに発見された特徴的なメラニン構造に関する特徴量は、加齢状態の評価に有用であることが明らかになった。
以上より、表層におけるメラニンの体積、特にメラニンのチューブ状の紋理構造の体積を評価指標として、加齢状態を評価することができることが明らかとなった。
また、前述した酸素飽和度と、メラニンの管状パターン構造の体積との相関関係を解析した。
図13に示す通り、酸素飽和度が減少するのに伴って、メラニンのチューブ状の紋理構造の体積が増加することが明らかになった。
したがって、メラニンのチューブ状の紋理構造の体積に基づいて、酸素飽和度の空間分布を推定できる。
また反対に、酸素飽和度に基づいて、メラニンのチューブ状の紋理構造の体積を推定できる。
図13に示す通り、酸素飽和度が減少するのに伴って、メラニンのチューブ状の紋理構造の体積が増加することが明らかになった。
したがって、メラニンのチューブ状の紋理構造の体積に基づいて、酸素飽和度の空間分布を推定できる。
また反対に、酸素飽和度に基づいて、メラニンのチューブ状の紋理構造の体積を推定できる。
本発明は、美容に関する肌の評価方法に利用することができる。
Claims (17)
- 被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力する、肌の評価方法。
- 前記内部構造を示すデータは、非侵襲的に得られたデータである、請求項1に記載の肌の評価方法。
- 前記モデルは、前記肌画像を入力データとし、前記特徴量を目的変数として与えて深層学習を行うことにより作成される畳み込みニューラルネットワークモデルである、請求項1又は2に記載の肌の評価方法。
- 入力された前記肌画像に基づいて、前記モデルの出力側から入力側に向かって分析することにより、入力された前記肌画像の領域ごとに、前記モデルによる推定結果への寄与を推定することを含む、請求項3に記載の肌の評価方法。
- 前記肌画像の領域ごとの前記推定結果への寄与を、前記肌画像に重ねて表示することを含む、請求項4に記載の肌の評価方法。
- 前記特徴量は、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚のメラニン構造に関する特徴量の少なくとも何れかである、請求項1~5の何れかに記載の肌の評価方法。
- 前記血管構造に関する特徴量は、皮膚上部における微小血管構造の量に関する特徴量である、請求項6に記載の肌の評価方法。
- 前記血管構造に関する特徴量は、血管の平均厚さ、血管の分岐の数及び皮膚上部における血管の体積から選ばれる、請求項6又は7に記載の肌の評価方法。
- 前記酸素飽和度に関する特徴量は、血管の酸素飽和度である、請求項6~8の何れかに記載の肌の評価方法。
- 前記メラニン構造に関する特徴量は、メラニンのチューブ状の紋理構造の体積である、請求項6~9の何れかに記載の肌の評価方法。
- 前記内部構造を示すデータは、光音響法により得られた測定データ又は該測定データに基づいて生成されたデータである、請求項6~10の何れかに記載の肌の評価方法。
- 皮膚の血管構造に関する特徴量に基づいて、被験者の皮膚の酸素飽和度を推定することを含む、請求項1~11の何れかに記載の肌の評価方法。
- 皮膚のメラニン構造に関する特徴量に基づいて、皮膚の酸素飽和度を推定することを含む、請求項1~11の何れかに記載の肌の評価方法。
- 皮膚の酸素飽和度に基づいて、皮膚の血管構造に関する特徴量及び/又は皮膚のメラニン構造に関する特徴量を推定することを含む、請求項1~11の何れかに記載の肌の評価方法。
- 前記モデルによって出力された前記推定値に基づいて、前記被験者の肌の加齢状態の指標を推定することを含む、請求項1~14の何れかに記載の肌の評価方法。
- 前記加齢状態は、前記特徴量及び年齢の組を複数与えることによって作成された相関モデルに、前記推定値を入力して、推定年齢を算出することによって評価される、請求項15に記載の肌の評価方法。
- 被験者の皮膚の内部構造を示すデータから得られる特徴量と、同一被験者の生体表面を撮影して得られる肌画像と、の複数の組を学習データとして、前記特徴量及び肌画像の関係を学習したモデルに、前記肌画像を入力することで、前記特徴量の推定値を出力するようにコンピュータを機能させる、肌の評価プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020172228A JP2022063800A (ja) | 2020-10-12 | 2020-10-12 | 肌の評価方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020172228A JP2022063800A (ja) | 2020-10-12 | 2020-10-12 | 肌の評価方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022063800A true JP2022063800A (ja) | 2022-04-22 |
Family
ID=81213130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020172228A Pending JP2022063800A (ja) | 2020-10-12 | 2020-10-12 | 肌の評価方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2022063800A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024085120A1 (ja) * | 2022-10-18 | 2024-04-25 | 株式会社資生堂 | 血管推定装置、血管推定方法、およびプログラム |
-
2020
- 2020-10-12 JP JP2020172228A patent/JP2022063800A/ja active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024085120A1 (ja) * | 2022-10-18 | 2024-04-25 | 株式会社資生堂 | 血管推定装置、血管推定方法、およびプログラム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rong et al. | A blood pressure prediction method based on imaging photoplethysmography in combination with machine learning | |
KR102597974B1 (ko) | 당뇨병성 순환계 합병증을 평가하기 위한 시스템 및 방법 | |
WO2015166503A1 (en) | Multimodal transcranial brain optical imaging | |
JP2002200050A (ja) | 肌色測定装置および肌色診断装置ならびに顔画像処理装置 | |
KR20180122622A (ko) | 조직 혈관 건강을 평가하는 방법 및 장치 | |
US11553844B2 (en) | Methods, systems and computer program products for calculating MetaKG signals for regions having multiple sets of optical characteristics | |
CN108135477A (zh) | 处理光学相干断层扫描 | |
D'Alessandro et al. | 3-D volume reconstruction of skin lesions for melanin and blood volume estimation and lesion severity analysis | |
Vyas et al. | Non-invasive estimation of skin thickness from hyperspectral imaging and validation using echography | |
US20230172565A1 (en) | Systems, devices, and methods for developing a model for use when performing oximetry and/or pulse oximetry and systems, devices, and methods for using a fetal oximetry model to determine a fetal oximetry value | |
Kainerstorfer et al. | Quantitative principal component model for skin chromophore mapping using multi-spectral images and spatial priors | |
JP2022063800A (ja) | 肌の評価方法 | |
Tomanic et al. | Estimating quantitative physiological and morphological tissue parameters of murine tumor models using hyperspectral imaging and optical profilometry | |
Else et al. | Effects of skin tone on photoacoustic imaging and oximetry | |
JP2022063799A (ja) | 肌の評価方法 | |
Granados-Castro et al. | Estimation of deoxygenated and oxygenated hemoglobin by multispectral blind linear unmixing | |
D’Alessandro et al. | Depth-dependent hemoglobin analysis from multispectral transillumination images | |
CN115474930A (zh) | 一种基于高光谱图像重建的无创血红蛋白检测方法 | |
JP7393654B2 (ja) | 唇年齢又は唇印象による唇の評価方法 | |
Ding et al. | Hyperspectral imaging with burn contour extraction for burn wound depth assessment | |
El-Sharkawy | Advancements in non-invasive hyperspectral imaging: Mapping blood oxygen levels and vascular health for clinical and research applications | |
Visscher | Imaging skin: past, present and future perspectives | |
Cano et al. | Deep learning assisted classification of spectral photoacoustic imaging of carotid plaques | |
Visscher et al. | From Image to Information: Image Processing in Dermatology and Cutaneous Biology | |
Benavides‐Lara et al. | High‐frequency photoacoustic and ultrasound imaging for skin evaluation: Pilot study for the assessment of a chemical burn |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20231004 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240430 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240521 |