ホーム > 特許ランキング > ポーラ化成工業株式会社
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-10-15
(45)【発行日】2024-10-23
(54)【発明の名称】肌の評価方法
(51)【国際特許分類】
A61B 5/00 20060101AFI20241016BHJP
A61B 5/107 20060101ALI20241016BHJP
A61B 8/13 20060101ALI20241016BHJP
A61B 5/1455 20060101ALI20241016BHJP
【ＦＩ】
A61B5/00 M
A61B5/107 800
A61B8/13
A61B5/1455
【請求項の数】
12
(21)【出願番号】P 2020172227
(22)【出願日】2020-10-12
(65)【公開番号】P2022063799
(43)【公開日】2022-04-22
【審査請求日】2023-10-04
(73)【特許権者】
【識別番号】000113470
【氏名又は名称】ポーラ化成工業株式会社
(74)【代理人】
【識別番号】100137338
【弁理士】
【氏名又は名称】辻田 朋子
(74)【代理人】
【識別番号】100196313
【弁理士】
【氏名又は名称】村松 大輔
(72)【発明者】
【氏名】水越 興治
(72)【発明者】
【氏名】▲濱▼中 祥弘
【審査官】増渕 俊仁
(56)【参考文献】
【文献】特開２０１０－２３３５８４（ＪＰ，Ａ）
【文献】特開２０１４－１０４１３２（ＪＰ，Ａ）
【文献】特開２０１７－０４７１９４（ＪＰ，Ａ）
【文献】国際公開第２０１９／２４４５８６（ＷＯ，Ａ１）
【文献】特開２０１９－１８７６８５（ＪＰ，Ａ）
【文献】特開２０１４－２２６４６４（ＪＰ，Ａ）
【文献】特開２０１４－１８０４５９（ＪＰ，Ａ）
【文献】特開平０７－０３９５５１（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１２／０３０００４９（ＵＳ，Ａ１）
【文献】韓国公開特許第１０－２０２０－００６９６９５（ＫＲ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／００－５／０１
Ａ６１Ｂ ５／０６－５／２２
Ａ６１Ｂ ８／００－８／１５
(57)【特許請求の範囲】
【請求項1】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
前記データは、三次元の画像データを含む、肌の評価方法。
【請求項2】
被験者の皮膚の内部構造を示すデータを取得する工程、取得した前記データから前記特徴量を抽出する工程、
抽出した前記特徴量を、予め用意した、前記特徴量と年齢との相関モデルに照合して、年齢の推定値を算出する工程、
を含む、請求項１に記載の肌の評価方法。
【請求項3】
前記データが、光音響法により得られた測定データ又は該測定データに基づいて生成されたデータである、請求項１又は２に記載の肌の評価方法。
【請求項4】
前記特徴量が、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚のメラニン構造に関する特徴量の少なくとも何れかである、請求項１～３の何れかに記載の肌の評価方法。
【請求項5】
前記血管構造に関する特徴量は、皮膚上部における微小血管構造の量に関する特徴量である、請求項４に記載の肌の評価方法。
【請求項6】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
前記特徴量が、皮膚の血管構造に関する特徴量であり、
前記血管構造に関する特徴量は、血管の平均厚さ、血管の分岐の数及び皮膚上部における血管の体積から選ばれる、肌の評価方法。
【請求項7】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
前記特徴量は、血管の酸素飽和度である、肌の評価方法。
【請求項8】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
前記特徴量は、メラニンのチューブ状の紋理構造の体積である、肌の評価方法。
【請求項9】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
前記特徴量が、皮膚の血管構造に関する特徴量であり、
前記血管構造に関する特徴量は、皮膚上部における微小血管構造の量に関する特徴量であり、
前記データは、光音響法により得られた測定データ又は該測定データに基づいて生成されたデータである、肌の評価方法。
【請求項10】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
非侵襲的に得られた被験者の皮膚の内部構造を示す画像データから得られる皮膚の血管構造に関する特徴量に基づいて、被験者の皮膚の酸素飽和度を推定することを含む、肌の評価方法。
【請求項11】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
非侵襲的に得られた皮膚の内部構造を示す画像データから得られる皮膚のメラニン構造に関する特徴量に基づいて、皮膚の酸素飽和度を推定することを含む、肌の評価方法。
【請求項12】
非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法であって、
非侵襲的に得られた皮膚の内部構造を示すデータから得られる皮膚の酸素飽和度に基づいて、皮膚の血管構造及び／又は皮膚のメラニン構造を推定することを含む、肌の評価方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、肌の評価方法に関する。
【背景技術】
【0002】
近年、測定手法の発達により、非侵襲的に皮膚の内部構造を測定することが可能になっている。
例えば、非特許文献１においては、組織透明化手法によって、血管の構造が明らかにされている。
例えば、非特許文献１においては、組織透明化手法によって、血管の構造が明らかにされている。
【0003】
また、非侵襲的に測定した皮膚の内部構造に基づいて、皮膚の状態を診断する方法も知られている。
例えば特許文献１においては、光音響波を用いて皮膚内部の構造を非侵襲的に測定することにより、皮膚の熱傷、壊疽、褥瘡の診断を行う方法が開示されている。
例えば特許文献１においては、光音響波を用いて皮膚内部の構造を非侵襲的に測定することにより、皮膚の熱傷、壊疽、褥瘡の診断を行う方法が開示されている。
【先行技術文献】
【非特許文献】
【0004】
【文献】Ｓｈｉｂａｔａ Ｔ， Ｋａｊｉｙａ Ｋ， Ｓａｔｏ Ｋ， Ｙｏｏｎ Ｊ， Ｋａｎｇ ＨＹ． ３Ｄ ｍｉｃｒｏｖａｓｃｕｌａｒ ａｎａｌｙｓｉｓ ｒｅｖｅａｌｓ ｉｒｒｅｇｕｌａｒｌｙ ｂｒａｎｃｈｉｎｇ ｂｌｏｏｄ ｖｅｓｓｅｌｓ ｉｎ ｔｈｅ ｈｙｐｅｒｐｉｇｍｅｎｔｅｄ ｓｋｉｎ ｏｆ ｓｏｌａｒ ｌｅｎｔｉｇｏ． Ｐｉｇｍｅｎｔ ｃｅｌｌ ＆ ｍｅｌａｎｏｍａ ｒｅｓｅａｒｃｈ．２０１８；３１（６）：７２５ ７．
【特許文献】
【0005】
【文献】特許第５８６０５３４号公報
【発明の概要】
【発明が解決しようとする課題】
【0006】
ところで、コラーゲン繊維の太さや多寡、皮膚内部に保持される水分量等、皮膚の内部の形態的特徴や生理学的特徴、すなわち皮膚の内部の構造が、肌の見た目（しみ、しわ、くすみ等）に影響することは、知られている。
しかし、組織透明化手法による皮膚内部の構造の解析は、皮膚を摘出して行うため、血液や酸素が失われてしまい、組織が生存している状態については、測定することができない。また、解析のためのｎ数を十分に確保しにくい。
すなわち、上記のような手法では、生体における皮膚の内部の構造を測定することは制限され、このような手法で得られた指標に基づいた肌の評価は、生体の状態を反映しているものであるとは必ずしも言えなかった。
しかし、組織透明化手法による皮膚内部の構造の解析は、皮膚を摘出して行うため、血液や酸素が失われてしまい、組織が生存している状態については、測定することができない。また、解析のためのｎ数を十分に確保しにくい。
すなわち、上記のような手法では、生体における皮膚の内部の構造を測定することは制限され、このような手法で得られた指標に基づいた肌の評価は、生体の状態を反映しているものであるとは必ずしも言えなかった。
【0007】
また、特許文献１に記載されているような、非侵襲的に測定した対象者の皮膚の内部構造に基づいて、皮膚の状態（熱傷、壊疽、褥瘡等の有無）を診断する方法は知られているものの、非侵襲的に測定した対象者の皮膚の内部構造に基づいて、当該対象者の肌の加齢状態について評価する方法については、検討されていなかった。
【0008】
上記状況に鑑み、本発明は、非侵襲的に測定した生体の皮膚内部の構造に基づく、新規な肌の評価方法、特に肌の加齢状態の評価方法を提供することを課題とする。
【課題を解決するための手段】
【0009】
上記課題を解決する本発明は、非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う、肌の評価方法である。
【0010】
本発明の好ましい形態では、肌の評価方法は、被験者の皮膚の内部構造を示すデータを取得する工程、取得したデータから前記特徴量を抽出する工程、抽出した特徴量を、予め用意した、前記特徴量と年齢との相関モデルに照合して、年齢の推定値を算出する工程、を含む。
【0011】
本発明の好ましい形態では、前記データが、画像データを含む。
【0012】
本発明の好ましい形態では、前記特徴量が、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚のメラニン構造に関する特徴量の少なくとも何れかである。
【0013】
また、本発明の好ましい形態では、前記血管構造に関する特徴量は、皮膚上部における微小血管構造の量に関する特徴量である。
【0014】
また、本発明の好ましい形態では、前記血管構造に関する特徴量は、血管の平均厚さ、血管の分岐の数及び皮膚上部における血管の体積から選ばれる。
【0015】
また、本発明の好ましい形態では、前記酸素飽和度に関する特徴量は、血管の酸素飽和度である。
【0016】
また、本発明の好ましい形態では、前記メラニン構造に関する特徴量は、メラニンのチューブ状の紋理構造の体積である。
前記各特徴量を用いることで、肌の加齢状態の評価を精度よく行うことができる。
前記各特徴量を用いることで、肌の加齢状態の評価を精度よく行うことができる。
【0017】
本発明の好ましい形態では、前記画像データは、光音響法により得られた測定データに基づいて生成される。
光音響法を用いることで、皮下脂肪層のような下部組織の構造も測定することが可能となる。
光音響法を用いることで、皮下脂肪層のような下部組織の構造も測定することが可能となる。
【0018】
また、本発明の好ましい形態では、非侵襲的に得られた被験者の皮膚の内部構造を示す画像データから得られる皮膚の血管構造に関する特徴量に基づいて、被験者の皮膚の酸素飽和度を推定することを含む。
【0019】
また、本発明の好ましい形態では、非侵襲的に得られた皮膚の画像データから得られる皮膚のメラニン構造に関する特徴量に基づいて、皮膚の酸素飽和度を推定することを含む。
【0020】
また、本発明の好ましい形態では、非侵襲的に得られた皮膚のデータから得られる皮膚の酸素飽和度に基づいて、皮膚の血管構造及び／又は皮膚のメラニン構造を推定することを含む。
本発明によれば、ある特徴量から、他の構造や性質を推定することにより、加齢状態を多面的に詳細に評価することができる。
本発明によれば、ある特徴量から、他の構造や性質を推定することにより、加齢状態を多面的に詳細に評価することができる。
【発明の効果】
【0021】
本発明によれば、生きた皮膚の状態に基づく肌の評価方法、特に加齢状態の評価方法が提供される。特に、光音響法により得られた測定データを用いる場合には、皮下脂肪層といった皮膚の下部の構造に基づく評価が可能となる。
【図面の簡単な説明】
【0022】
【図1】（Ａ）５７５ｎｍの光音響信号から得られた画像と、６５０ｎｍの超音波信号から得られた画像をマージした画像である。（Ｂ）光音響法により得られた測定データから生成された、２０代被験者の皮膚の上部における血管構造を表す図である。（Ｃ）光音響法により得られた測定データから生成された、６０代被験者の皮膚の上部における血管構造を表す図である。
【図2】（Ａ）深層学習において訓練データとして用いられた測定画像データ、注釈画像データ及び、深層学習による血管構造の予測画像データである。（Ｂ）深層学習において訓練データとして用いられていない測定画像データ、注釈画像データ及び、深層学習による血管構造の予測画像データである。（Ｃ）予測画像データ（ｘ－ｙ平面）をｚ軸方向に結合することで生成した三次元データの代表例である。
【図3】被験者の年齢と、皮膚の各層における血管構造の平均厚さとの関係を示す図である。
【図4】（Ａ）２０代被験者の測定領域における、酸素飽和度の空間分布を示す図である。（Ｂ）６０代被験者の測定領域における、酸素飽和度の空間分布を示す図である。
【図5】被験者の年齢と、皮膚の各層における酸素飽和度との関係を示す図である。
【図6】皮膚の各層における血管構造の平均厚さと、対応する層における酸素飽和度との関係を示す図である。
【図7】（Ａ）２０代被験者の測定領域における、メラニン構造の分布を示す図である。（Ｂ）６０代被験者の測定領域における、メラニン構造の分布を示す図である。
【図8】（Ａ）深層学習において訓練データとして用いられた測定画像データ、注釈画像データ及び、深層学習によるメラニン構造の予測画像データである。（Ｂ）深層学習において訓練データとして用いられていない測定画像データ、注釈画像データ及び、深層学習によるメラニン構造の予測画像データである。（Ｃ）予測画像データ（ｘ－ｙ平面）をｚ軸方向に結合することで生成した三次元データの代表例である。
【図9】被験者の年齢と、皮膚の上部に存在するメラニンのチューブ状の紋理構造との関係を示す図である。
【図10】皮膚の各層における酸素飽和度と、対応する層におけるメラニンのチューブ状の紋理構造の体積との関係を示す図である。
【発明を実施するための形態】
【0023】
本発明において、「皮膚」とは、表皮層、真皮層及び皮下脂肪層を含む生体組織を指し、組織の構造に特に着目する場合に用いる。
また、本発明において、「肌」は、生体組織としては「皮膚」と同義であるが、組織の構造に特に着目しない場合に用いる。
また、本発明の「評価」は、医療診断を含まない。
また、本発明において、「肌」は、生体組織としては「皮膚」と同義であるが、組織の構造に特に着目しない場合に用いる。
また、本発明の「評価」は、医療診断を含まない。
【0024】
＜１＞肌の評価方法（加齢状態の評価）
本発明の肌の評価方法は、非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う工程を含む。
また、本発明の好ましい形態では、被験者の皮膚の内部構造を示すデータを取得する工程、取得したデータから特徴量を抽出する工程、抽出した特徴量を、予め用意した、前記特徴量と年齢との相関モデルに照合して、年齢の推定値を算出する工程、を含む。
本発明によれば、皮膚を侵襲的に測定したデータに基づいて肌の加齢状態を評価することができるので、簡便に被験者の肌の加齢状態を評価することができる。
本発明の肌の評価方法は、非侵襲的に得られた被験者の皮膚の内部構造を示すデータから得られる特徴量に基づいて、被験者の肌の加齢状態の評価を行う工程を含む。
また、本発明の好ましい形態では、被験者の皮膚の内部構造を示すデータを取得する工程、取得したデータから特徴量を抽出する工程、抽出した特徴量を、予め用意した、前記特徴量と年齢との相関モデルに照合して、年齢の推定値を算出する工程、を含む。
本発明によれば、皮膚を侵襲的に測定したデータに基づいて肌の加齢状態を評価することができるので、簡便に被験者の肌の加齢状態を評価することができる。
【0025】
本発明において、「肌の加齢状態」とは、見た目上の加齢の程度、加齢印象（見た目年齢）等を例示することができる。
【0026】
本発明において用いる皮膚の内部構造を示すデータは、非侵襲的な測定方法により得られる。前記データは、公知の装置及び公知の方法を適宜用いて取得することができる。公知の非侵襲的な測定方法としては、例えば、高周波電流法、ＡＴＲ－ＦＴＩＲ法、近赤外分光法、磁気共鳴画像法（ＭＲＩ）、Ｉｎ ｖｉｖｏ共焦点ラマン分光法、超音波断層撮影、光コヒーレンス・トモグラフィー、Ｉｎ ｖｉｖｏ共焦点レーザー走査型顕微鏡、多光子顕微鏡、ＳＨＧ顕微鏡、ＣＡＲＳ顕微鏡、光超音波顕微鏡等を選択することができる。
また、上記非侵襲的な測定方法で測定することができれば、皮膚の内部構造を示すデータの種類は特に制限されず、例えば被験者の皮膚のある測定領域における血管構造に関するデータ、血流に関するデータ、酸素飽和度に関するデータ、メラニンのミクロ空間分布に関するデータ、水分に関するデータ、皮膚の粘弾性に関するデータ、特定の物質に関するデータ等を例示することができる。
また取得した皮膚の内部構造を示すデータは、生データのまま使用してもよく、一般的な統計学的加工（平均や標準偏差の取得等）や画像化等の加工を施してもよい。
本発明の好ましい形態では、非侵襲的な測定方法として、光音響法を採用する。また本発明の好ましい形態では、皮膚の内部構造を示すデータは、光音響法により取得したデータ、例えば画像データである。
また、上記非侵襲的な測定方法で測定することができれば、皮膚の内部構造を示すデータの種類は特に制限されず、例えば被験者の皮膚のある測定領域における血管構造に関するデータ、血流に関するデータ、酸素飽和度に関するデータ、メラニンのミクロ空間分布に関するデータ、水分に関するデータ、皮膚の粘弾性に関するデータ、特定の物質に関するデータ等を例示することができる。
また取得した皮膚の内部構造を示すデータは、生データのまま使用してもよく、一般的な統計学的加工（平均や標準偏差の取得等）や画像化等の加工を施してもよい。
本発明の好ましい形態では、非侵襲的な測定方法として、光音響法を採用する。また本発明の好ましい形態では、皮膚の内部構造を示すデータは、光音響法により取得したデータ、例えば画像データである。
【0027】
ここで、皮膚の内部構造を示すデータの取得のための測定領域としては、任意の箇所を選択することができるが、日常生活において光を浴びやすい箇所（例えば、顔や手の甲）の皮膚を選択することが好ましい。また、日常生活において光を浴びにくい箇所（例えば、上腕の内側）の皮膚を選択することも、好ましい。
前者を選択する場合には、光老化の要素を含む総合的な加齢状態を評価することができ、後者を選択する場合には、光老化の要素を極力排除した（生理的な）加齢状態を評価することができる。
また、測定領域としては、視認できるシミやその他色素沈着、創傷等の異常がない正常な領域を設定する。
また、測定領域のサイズとしては、好ましくは５ｍｍ×５ｍｍ～１５ｍｍ×１５ｍｍで設定することができる。
前者を選択する場合には、光老化の要素を含む総合的な加齢状態を評価することができ、後者を選択する場合には、光老化の要素を極力排除した（生理的な）加齢状態を評価することができる。
また、測定領域としては、視認できるシミやその他色素沈着、創傷等の異常がない正常な領域を設定する。
また、測定領域のサイズとしては、好ましくは５ｍｍ×５ｍｍ～１５ｍｍ×１５ｍｍで設定することができる。
【0028】
本発明において、皮膚の内部構造を示すデータとして画像データを使用する場合、特徴量を抽出するための画像データは、三次元データであることが好ましい。このような三次元データとしては、光音響法等によって得られた測定画像データそのものであってもよいし、測定画像データを加工した加工画像データであってもよい。
評価において、予め用意した相関モデルを用いる場合には、当該相関モデルを作成する際に解析した画像データを用いて特徴量を抽出する。加工画像データの作成の形態については、各特徴量の抽出方法として後述する。
評価において、予め用意した相関モデルを用いる場合には、当該相関モデルを作成する際に解析した画像データを用いて特徴量を抽出する。加工画像データの作成の形態については、各特徴量の抽出方法として後述する。
【0029】
皮膚の内部構造を示すデータとして画像データを使用する場合、当該皮膚の画像データとしては、表皮層、真皮層、皮下脂肪層のすべての層を含む画像データを用いてもよいし、これを皮膚の深さ方向に任意に分割した画像データを用いてもよい。
例えば、皮膚の深さ方向に任意に分割した画像データとしては、表皮層から皮下脂肪層までを、上部（表皮から真皮乳頭層までの領域）、中間部（真皮網状層から皮下脂肪層の境界領域）、及び下部（皮下脂肪層の上層）に分割した画像データを用いることができる。
例えば、皮膚の深さ方向に任意に分割した画像データとしては、表皮層から皮下脂肪層までを、上部（表皮から真皮乳頭層までの領域）、中間部（真皮網状層から皮下脂肪層の境界領域）、及び下部（皮下脂肪層の上層）に分割した画像データを用いることができる。
【0030】
本発明の肌の評価に用いることができる、皮膚の内部構造を示すデータから得られる皮膚の特徴量は、例えば、血管構造に関する特徴量、血流に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量、水分に関する特徴量を挙げることができる。
中でも好ましい形態として、血管構造に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量を挙げることができる。
このうち、血管構造に関する特徴量を用いる場合には、皮膚の上部の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、メラニン構造に関する特徴量を用いる場合には、表皮層の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、酸素飽和度に関する特徴量を用いる場合には、表皮層の測定データを含むデータ、又は皮下脂肪層の測定データを含むデータを用いることが好ましい。また、当該データは、酸素飽和度の高低を色分けして示す画像データであってもよい。
中でも好ましい形態として、血管構造に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量を挙げることができる。
このうち、血管構造に関する特徴量を用いる場合には、皮膚の上部の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、メラニン構造に関する特徴量を用いる場合には、表皮層の測定データを含むデータを用いることが好ましい。また、当該データは、画像データであることが好ましい。
また、酸素飽和度に関する特徴量を用いる場合には、表皮層の測定データを含むデータ、又は皮下脂肪層の測定データを含むデータを用いることが好ましい。また、当該データは、酸素飽和度の高低を色分けして示す画像データであってもよい。
【0031】
本発明の評価方法では、血管構造に関する特徴量、酸素飽和度に関する特徴量、メラニン構造に関する特徴量のうち、１種類を用いてもよく、２種類以上を用いてもよい。
【0032】
血管構造に関する特徴量としては、微細な血管構造に関する特徴量を用いることが好ましく、特に、皮膚の上部における微細な血管構造に関する特徴量を用いることが好ましい。このような特徴量として、血管の分岐の数、血管の曲率、血管の合計量（体積）、微細血管の数、血管の密度、血管の平均厚さ等が好ましく挙げられ、中でも、血管の分岐の数、血管の合計量（体積）、微細血管の数、血管の平均厚さを用いることが好ましい。
【0033】
上記血管構造に関する特徴量は、上記非侵襲的な測定方法により測定した測定領域における血管構造に関するデータから、抽出することができる。
血管構造に関する特徴量を抽出するための画像データは、例えば、５７５ｎｍと６５０ｎｍの二波長を含むパルス光を測定部位に照射し、５７５ｎｍの光音響シグナルから、６５０ｎｍの光音響シグナルを差し引いた差分値を用いて、皮膚内部の血管構造を示す三次元画像データを生成することができる。例えば、株式会社アドバンテストのＨａｄａｔｏｍｏＺ（Ｈａｄａｔｏｍｏは登録商標）を用いて、このような三次元画像データ（測定画像データ）の生成を行うことができる。そして、生成した画像を解析することで、血管構造に関する特徴量を抽出することができる。
このような測定画像データを用いて、血管構造に関する特徴量を定量化してもよい。
好ましい形態では、この測定画像データを加工することにより、血管構造に関する特徴量の抽出の精度を上げることができる。
血管構造に関する特徴量を抽出するための画像データは、例えば、５７５ｎｍと６５０ｎｍの二波長を含むパルス光を測定部位に照射し、５７５ｎｍの光音響シグナルから、６５０ｎｍの光音響シグナルを差し引いた差分値を用いて、皮膚内部の血管構造を示す三次元画像データを生成することができる。例えば、株式会社アドバンテストのＨａｄａｔｏｍｏＺ（Ｈａｄａｔｏｍｏは登録商標）を用いて、このような三次元画像データ（測定画像データ）の生成を行うことができる。そして、生成した画像を解析することで、血管構造に関する特徴量を抽出することができる。
このような測定画像データを用いて、血管構造に関する特徴量を定量化してもよい。
好ましい形態では、この測定画像データを加工することにより、血管構造に関する特徴量の抽出の精度を上げることができる。
【0034】
例えば、深層学習モデルを用いて、測定画像データを水平方向画像データ（平面画像データ）に分解し、この分解した測定平面画像データから血管構造を示す要素を抽出した抽出平面画像データを生成し、この抽出平面画像データを結合することで、血管構造が抽出された抽出三次元画像データを生成する方法があげられる。
測定平面画像データから抽出平面画像データを生成するためには、深層学習モデルを使用することができる。
この場合、データキュレーターが測定平面画像データを目視評価することで血管構造を抽出した注釈画像データを作成する。
訓練に必要な数の測定平面画像データと対応する注釈画像データを訓練用画像データとして用いて、深層学習をさせる。深層学習モデルとしては、エンコーダー・デコーダーモデル構造を有するネットワークであるＵｎｅｔを用いることができる。
測定平面画像データから抽出平面画像データを生成するためには、深層学習モデルを使用することができる。
この場合、データキュレーターが測定平面画像データを目視評価することで血管構造を抽出した注釈画像データを作成する。
訓練に必要な数の測定平面画像データと対応する注釈画像データを訓練用画像データとして用いて、深層学習をさせる。深層学習モデルとしては、エンコーダー・デコーダーモデル構造を有するネットワークであるＵｎｅｔを用いることができる。
【0035】
酸素飽和度に関する特徴量としては、測定領域（例えば、皮膚の任意の層）の酸素飽和度、高酸素飽和度の微細血管が存在する領域の量などを用いることができる。
【0036】
酸素飽和度に関する特徴量を抽出するための画像データは、例えば、５７５ｎｍと５８６ｎｍの二波長を含むパルス光を測定部位に照射し、酸素ヘモグロビンと脱酸素ヘモグロビンとの吸収係数の差を利用して、全体に対する酸素ヘモグロビンの割合を色分けなどで可視化した測定画像を生成することで、取得することができる。例えば、株式会社アドバンテストのＨａｄａｔｏｍｏＺ（Ｈａｄａｔｏｍｏは登録商標）を用いて、このような三次元画像データ（測定画像データ）の生成を行うことができる。そして、生成した画像を解析することで、酸素飽和度に関する特徴量を抽出することができる。
このような測定画像データを用いて、酸素飽和度に関する特徴量を定量化してもよい。
また、ある位置（領域）の各波長における光音響信号の振幅比に基づいて、当該位置の酸素飽和度を算出してもよい（実施例参照）。
このような測定画像データを用いて、酸素飽和度に関する特徴量を定量化してもよい。
また、ある位置（領域）の各波長における光音響信号の振幅比に基づいて、当該位置の酸素飽和度を算出してもよい（実施例参照）。
【0037】
メラニン構造に関する特徴量は、メラニン顆粒の凝集により形成される構造に関する特徴量である。ここで、本発明においてメラニン構造に関連する特徴量は、視認できるようなシミやそばかす等に関連するメラニンの凝集などの構造ではなく、通常は視認できない１ｍｍ未満四方の空間における構造を意味する。
光音響法では、通常個々のメラニン顆粒は測定することができない。
この構造は、１０～数百μｍの範囲の大きさの構造である。また、この構造は、好ましくは、チューブ状の紋理構造である（図７（Ｂ）参照）。
メラニン構造に関する特徴量としては、前記構造の存在量に関する特徴量があげられ、体積、密度、数があげられる。
光音響法では、通常個々のメラニン顆粒は測定することができない。
この構造は、１０～数百μｍの範囲の大きさの構造である。また、この構造は、好ましくは、チューブ状の紋理構造である（図７（Ｂ）参照）。
メラニン構造に関する特徴量としては、前記構造の存在量に関する特徴量があげられ、体積、密度、数があげられる。
【0038】
メラニン構造に関する特徴量を抽出するための画像データは、例えば、６５０ｎｍのパルス光を測定部位に照射し、得られた光音響データに基づいた測定画像を生成することで、取得することができる。例えば、株式会社アドバンテストのＨａｄａｔｏｍｏＺ（Ｈａｄａｔｏｍｏは登録商標）を用いて、このような三次元画像データ（測定画像データ）の生成を行うことができる。そして、生成した画像を解析することで、メラニン構造に関する特徴量を抽出することができる。
このような測定画像データを用いて、メラニン構造に関する特徴量を定量化してもよい。
このような測定画像データを用いて、メラニン構造に関する特徴量を定量化してもよい。
【0039】
実施例に示す通り、前述した各特徴量は、年齢と関係を有することを、本発明者らは明らかにしている。したがって当該関係に基づいて、非侵襲的に取得した皮膚内部の画像データから得られる前記特徴量に基づいて、肌の加齢状態を評価することができる。
【0040】
取得した皮膚の内部構造を示すデータから特徴量を抽出する方法は特に制限されず、非侵襲的な測定方法により得られた、皮膚の内部構造を示すデータに一般的な統計学的処理（例えば、平均、標準偏差等の算出）を施すことで特徴量を抽出する形態や、画像化等の加工を施した後に、当該画像を解析して特徴量を抽出する形態としてもよい。
【0041】
本発明では、肌の加齢状態の評価においては、あらかじめ特徴量と年齢の相関関係を示す相関モデルを作成しておくことが、簡便性の問題から好ましい。
以下に相関関係の解析と相関モデルの作成について説明する。
以下に相関関係の解析と相関モデルの作成について説明する。
【0042】
本発明の評価方法に用いる特徴量（評価指標）は、以下のように決定することができる。
まず、非侵襲的な測定方法を用いて、年代の異なる複数被験者の皮膚の測定領域を測定して、皮膚の内部構造のデータを得る。皮膚の内部構造のデータは、好ましくは画像データであるか、又は得られたデータを画像化したものである。皮膚の内部構造のデータを画像データとする場合、当該画像は、皮膚断面画像等の二次元画像であっても、三次元画像であってもよい。
まず、非侵襲的な測定方法を用いて、年代の異なる複数被験者の皮膚の測定領域を測定して、皮膚の内部構造のデータを得る。皮膚の内部構造のデータは、好ましくは画像データであるか、又は得られたデータを画像化したものである。皮膚の内部構造のデータを画像データとする場合、当該画像は、皮膚断面画像等の二次元画像であっても、三次元画像であってもよい。
【0043】
非侵襲的な測定方法としては、上記した各種装置及び方法を使用することができるが、光音響法を採用することが好ましい。
光音響法を採用することで、従来の方法では測定できなかった生体深部の構造を、より正確に測定することができる。
例えば、皮膚の内部構造の画像データを得る場合、年代の異なる複数の被験者の皮膚の測定対象部位にパルス光を照射して、発生した超音波を測定し、測定された皮膚内部の超音波のデータに基づいて、各年代の被験者についてそれぞれ皮膚内部の構造の画像データを生成する。
光音響法を採用することで、従来の方法では測定できなかった生体深部の構造を、より正確に測定することができる。
例えば、皮膚の内部構造の画像データを得る場合、年代の異なる複数の被験者の皮膚の測定対象部位にパルス光を照射して、発生した超音波を測定し、測定された皮膚内部の超音波のデータに基づいて、各年代の被験者についてそれぞれ皮膚内部の構造の画像データを生成する。
【0044】
次に、得られた測定データに含まれるパラメータを抽出し、実年齢又は見た目年齢と相関のあるパラメータを探索する。
測定データからのパラメータの抽出方法は特に制限されず、公知の解析手法を用いることができる。例えば、測定データとして画像データを使用する場合、年代ごとの画像を目視で比較し、変化の見られた皮膚のパラメータを、抽出する方法を採用することができる。また、汎用の画像処理ソフトを用いた解析によって、色分布、輝度分布、空間周波数強度等を含む特徴を抽出することによって、パラメータを抽出することもできる。
測定データからのパラメータの抽出方法は特に制限されず、公知の解析手法を用いることができる。例えば、測定データとして画像データを使用する場合、年代ごとの画像を目視で比較し、変化の見られた皮膚のパラメータを、抽出する方法を採用することができる。また、汎用の画像処理ソフトを用いた解析によって、色分布、輝度分布、空間周波数強度等を含む特徴を抽出することによって、パラメータを抽出することもできる。
【0045】
次に、抽出したパラメータと、見た目年齢又は実年齢との相関を解析する。なお、見た目年齢として、年齢を年代やスコア等に適宜変換した値を用いてもよい。
解析手法としては、回帰分析など、相関関係の解析に一般的に用いられる方法を採用することができる。
解析手法としては、回帰分析など、相関関係の解析に一般的に用いられる方法を採用することができる。
【0046】
また、測定データとして画像データを使用する場合、非侵襲的な測定方法により得られた皮膚の内部構造を示す画像データに、見た目年齢又は実年齢等、加齢の度合いを示すデータのラベルを付したものを訓練データとした機械学習により、相関関係を解析することもできる。
【0047】
解析の結果、年齢と相関があるとされたパラメータを、肌の評価指標として用いることができる特徴量として決定し、前記特徴量と年齢との相関関係を説明する相関モデルを作成する。相関モデルとしては、重回帰式、判別式、ニューラルネットワーク等があげられる。
【0048】
発明者らの解析によって、具体的には、実年齢が高いほど、被験者の測定対象部位における血管の平均厚さが厚い（より血管が厚くなっている）という関係が見出された（図３参照）。
また、実年齢が高いほど、血管の体積が小さいこと、実年齢が高いほど皮膚の上部においては血管の分岐の数が少ないという関係が見いだされた。
また、実年齢が高いほど、測定領域における酸素飽和度が低く、特に皮膚の上部と下部においてこの関係が強い傾向が見出された（図５参照）。
また、実年齢が高いほど、血管の体積が小さいこと、実年齢が高いほど皮膚の上部においては血管の分岐の数が少ないという関係が見いだされた。
また、実年齢が高いほど、測定領域における酸素飽和度が低く、特に皮膚の上部と下部においてこの関係が強い傾向が見出された（図５参照）。
【0049】
また、実年齢が高いほど、単位皮膚領域あたりに観察されるメラニンの紋理構造が多い（体積が大きい）という関係が見出された（図９参照）。
【0050】
したがって、これらの特徴量は、本発明の肌の評価方法に用いる特徴量（評価指標）として特に好ましく用いることができる。
【0051】
本発明においては、上記特徴量のほかにも、肌の加齢状態の評価の指標として知られている項目を併せて測定し、これを加味して肌の加齢状態を評価してもよい。
例えば、公知の方法により、肌の水分量や粘弾性、加齢形質の多少等を非侵襲的に測定し、これらの測定結果と、上記特徴量の、光音響法による測定結果の少なくとも何れかの結果とを組み合わせ、被験者の肌の加齢状態を評価する形態とすることができる。
例えば、公知の方法により、肌の水分量や粘弾性、加齢形質の多少等を非侵襲的に測定し、これらの測定結果と、上記特徴量の、光音響法による測定結果の少なくとも何れかの結果とを組み合わせ、被験者の肌の加齢状態を評価する形態とすることができる。
【0052】
特徴量から、肌を評価する方法としては、例えば、前述した相関分析に基づいて、あらかじめ特徴量と加齢状態の推定値（見た目年齢やスコア）との相関関係を示す式やモデルをコンピュータに記憶しておき、当該式やモデルに、光音響法により得られた測定データ（又はこれから得られた画像データ）を入力することで、加齢状態の推定値を出力する方法が好ましく用いられる。
【0053】
＜２＞評価方法（内部構造の推定）
また、発明者らは、上記光音響法を用いた手法により、皮膚の測定領域における血管の平均厚さと当該測定領域における酸素飽和度、及び、皮膚の測定領域における酸素飽和度と当該測定領域におけるメラニンの紋理構造の体積が、それぞれ互いに相関関係を有することを明らかにした（図６及び図１０参照）。
したがって、本発明は、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚におけるメラニン構造に関する特徴量の少なくとも何れかを解析することにより、相互に他の皮膚構造を推定することを含む肌の評価方法を提供する。
また、発明者らは、上記光音響法を用いた手法により、皮膚の測定領域における血管の平均厚さと当該測定領域における酸素飽和度、及び、皮膚の測定領域における酸素飽和度と当該測定領域におけるメラニンの紋理構造の体積が、それぞれ互いに相関関係を有することを明らかにした（図６及び図１０参照）。
したがって、本発明は、皮膚の血管構造に関する特徴量、皮膚の酸素飽和度に関する特徴量、及び皮膚におけるメラニン構造に関する特徴量の少なくとも何れかを解析することにより、相互に他の皮膚構造を推定することを含む肌の評価方法を提供する。
【0054】
具体的には、皮膚の測定領域における血管の平均厚さが高い（より厚い）ほど、当該測定領域における血管の酸素飽和度が低いという関係がある（図６参照）。この関係は、特に皮膚の上部及び下部において顕著である。
したがって、皮膚の測定領域における血管の平均厚さを測定することによって、当該測定領域における血管の酸素飽和度を推定することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って減少するから、皮膚の測定領域における血管の平均厚さの測定データに基づいて、当該測定領域における酸素飽和度を推定することによって、肌（特に、肌の加齢状態）を評価することができる。
したがって、皮膚の測定領域における血管の平均厚さを測定することによって、当該測定領域における血管の酸素飽和度を推定することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って減少するから、皮膚の測定領域における血管の平均厚さの測定データに基づいて、当該測定領域における酸素飽和度を推定することによって、肌（特に、肌の加齢状態）を評価することができる。
【0055】
また上記測定領域における血管の平均厚さと当該測定領域における酸素飽和度との関係に基づけば、皮膚の測定領域における酸素飽和度を測定することによって、当該測定領域における血管の平均厚さを推定することができる。
そして上記の通り、皮膚の測定領域における血管の平均厚さは、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域における血管の平均厚さを推定することによって、肌（特に、肌の加齢状態）を評価することができる。
そして上記の通り、皮膚の測定領域における血管の平均厚さは、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域における血管の平均厚さを推定することによって、肌（特に、肌の加齢状態）を評価することができる。
【0056】
また、皮膚の測定領域における酸素飽和度が高いほど、当該測定領域に存在するメラニンの紋理構造の体積が少ないという関係がある（図１０参照）。
したがって、皮膚の測定領域における酸素飽和度を測定することによって、当該測定領域において観察されるメラニンの紋理構造の体積を推定することができる。
そして上記の通り、皮膚の測定領域において観察されるメラニンの紋理構造の体積は、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域において観察されるメラニンの紋理構造の体積を推定することによって、肌（特に、肌の加齢状態）を評価することができる。
したがって、皮膚の測定領域における酸素飽和度を測定することによって、当該測定領域において観察されるメラニンの紋理構造の体積を推定することができる。
そして上記の通り、皮膚の測定領域において観察されるメラニンの紋理構造の体積は、加齢に伴って増加するから、皮膚の測定領域における酸素飽和度の測定データに基づいて、当該測定領域において観察されるメラニンの紋理構造の体積を推定することによって、肌（特に、肌の加齢状態）を評価することができる。
【0057】
また上記測定領域における酸素飽和度と、当該測定領域に存在するメラニンの紋理構造の体積との関係に基づけば、皮膚の測定領域において観察されるメラニンの紋理構造の体積を測定することによって、当該測定領域における酸素飽和度を推定することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って現象するから、皮膚の測定領域において観察されるメラニンの紋理構造の体積の測定データに基づいて、当該測定領域における酸素飽和度を推定することによって、肌（特に、肌の加齢状態）を評価することができる。
そして上記の通り、皮膚の測定領域における酸素飽和度は、加齢に伴って現象するから、皮膚の測定領域において観察されるメラニンの紋理構造の体積の測定データに基づいて、当該測定領域における酸素飽和度を推定することによって、肌（特に、肌の加齢状態）を評価することができる。
【0058】
各特徴量やその抽出方法の好ましい形態は、＜１＞で述べた通りである。
【実施例】
【0059】
被験者（２０代～６０代の日本人女性、８０名）の、シミや赤斑等が肉眼で観察される部分を除いた正常な頬部位（９ｍｍ×９ｍｍ）について、以下の方法により、光音響測定画像データを得て、様々な特徴量を抽出した。
そして、抽出した特徴量と、年齢との相関を解析することにより、加齢状態に関係のある特徴量を解析した。
そして、抽出した特徴量と、年齢との相関を解析することにより、加齢状態に関係のある特徴量を解析した。
【0060】
＜１＞光音響測定画像の取得
光音響測定には、株式会社アドバンテストのＨａｄａｔｏｍｏＺ光音響測定装置を用いた。この装置は、二波長光音響イメージングおよび超音波イメージングが可能である。測定は、３０μｍの走査ステップで行った。
光音響測定には、株式会社アドバンテストのＨａｄａｔｏｍｏＺ光音響測定装置を用いた。この装置は、二波長光音響イメージングおよび超音波イメージングが可能である。測定は、３０μｍの走査ステップで行った。
【0061】
測定には、２種類の二波長レーザー、（ｉ）５７５ｎｍおよび５８６ｎｍ、ならびに（ｉｉ）５７５ｎｍおよび６５０ｎｍを使用した。
血管のイメージングのためには、５７５ｎｍの光音響信号から６５０ｎｍの光音響信号を差し引いた値を使用した。
酸素飽和度分布のイメージングのためには、５７５ｎｍの光音響信号および５８６ｎｍの光音響信号から得られるオキシヘモグロビンとデオキシヘモグロビンとの間の吸収係数の差を使用した。
メラニン分布のイメージングのためには、６５０ｎｍの光音響信号を使用した。
血管のイメージングのためには、５７５ｎｍの光音響信号から６５０ｎｍの光音響信号を差し引いた値を使用した。
酸素飽和度分布のイメージングのためには、５７５ｎｍの光音響信号および５８６ｎｍの光音響信号から得られるオキシヘモグロビンとデオキシヘモグロビンとの間の吸収係数の差を使用した。
メラニン分布のイメージングのためには、６５０ｎｍの光音響信号を使用した。
【0062】
＜２＞血管構造及びメラニン構造の解析のための画像生成
血管構造やメラニン構造に関する特徴量を抽出するために、ＭＡＴＬＡＢ（登録商標） ２０２０ａ（ＭａｔｈＷｏｒｋｓ、ＵＳＡ）、Ｐｙｔｈｏｎ３．７、およびＴｅｎｓｏｒｆｌｏｗを使用して、光音響測定により得られた画像データに基づいて、解析画像データを生成した。
光音響測定で得た三次元画像データをｚ軸上に積層した画像から水平方向（ｘ-ｙ平面）画像に再構成した。再構成した画像から血管構造およびメラニン構造を抽出するために、エンコーダー・デコーダーモデル構造を有するネットワークであるＵｎｅｔを、深層学習セグメント化方法として使用した。
血管構造およびメラニン構造が観察された測定領域の画像を、ｚ軸の任意の複数の位置におけるｘ-ｙ平面画像を切り出し、データキュレーターによって血管構造またはメラニン凝集構造と視覚的に見なされる領域についての注釈画像データ（訓練データ）として抽出画像を作成した。これらの抽出画像は、深層学習における出力データ、およびそのクラスラベルとした。これらの抽出画像は、血管構造およびメラニン構造について、それぞれ９０例および３０例であった。次に、元のｘ-ｙ平面画像を入力データとして使用して血管およびメラニン構造を抽出するために、５回反復および２×２プーリングを有するネットワークを構築した。目的関数は、以下の式（１）によって計算されるダイス関数であると仮定した。
ここで、ｙは元のｘ-ｙ平面画像の画素ベクトルを表し、y~は、注釈画像データの画素ベクトルを表す。
血管構造やメラニン構造に関する特徴量を抽出するために、ＭＡＴＬＡＢ（登録商標） ２０２０ａ（ＭａｔｈＷｏｒｋｓ、ＵＳＡ）、Ｐｙｔｈｏｎ３．７、およびＴｅｎｓｏｒｆｌｏｗを使用して、光音響測定により得られた画像データに基づいて、解析画像データを生成した。
光音響測定で得た三次元画像データをｚ軸上に積層した画像から水平方向（ｘ-ｙ平面）画像に再構成した。再構成した画像から血管構造およびメラニン構造を抽出するために、エンコーダー・デコーダーモデル構造を有するネットワークであるＵｎｅｔを、深層学習セグメント化方法として使用した。
血管構造およびメラニン構造が観察された測定領域の画像を、ｚ軸の任意の複数の位置におけるｘ-ｙ平面画像を切り出し、データキュレーターによって血管構造またはメラニン凝集構造と視覚的に見なされる領域についての注釈画像データ（訓練データ）として抽出画像を作成した。これらの抽出画像は、深層学習における出力データ、およびそのクラスラベルとした。これらの抽出画像は、血管構造およびメラニン構造について、それぞれ９０例および３０例であった。次に、元のｘ-ｙ平面画像を入力データとして使用して血管およびメラニン構造を抽出するために、５回反復および２×２プーリングを有するネットワークを構築した。目的関数は、以下の式（１）によって計算されるダイス関数であると仮定した。
【数1】
【0063】
前記構築したネットワークを用いて、被験者から取得した光音響データから、血管構造およびメラニン構造が抽出されたｘ-ｙ平面画像データを抽出し、この抽出した平面画像データを深さ方向で再び統合して、三次元データを作成した。
【0064】
作成した三次元データに基づいて、長さ、厚さ、分岐の数、体積等のパラメータを測定した。皮膚の深さ方向に上部、中間部及び下部とわけ、そのそれぞれについて、種々のパラメータを測定した。
【0065】
＜３＞酸素飽和度の計算
酸素飽和度の空間分布ＳＯ２（ｒ）は式（２）に示すように、酸化ヘモグロビンＨｂＯ２と脱酸素ヘモグロビンＨｂとの比率に基づいて算出することができる。
ここで、ｒは位置情報、εＨｂは脱酸素ヘモグロビンのモル吸収係数、εΔＨｂは酸化ヘモグロビンと脱酸素ヘモグロビンのモル吸収係数の差分値、μａ（ｒ）は光音響信号である。λ１とλ２はレーザー波長、波長λ１は５７５ｎｍ、波長λ２は５８６ｎｍである。波長λ１およびλ２に対応する脱酸素ヘモグロビンおよび酸化ヘモグロビンのモル吸収係数が既知の値であると仮定すると、酸素飽和度の空間分布ＳＯ２（ｒ）は、各波長における光音響信号の振幅比に基づいて算出することができる。皮膚の深さ方向に上部、中間部及び下部とわけ、そのそれぞれについて、酸素飽和度に関する種々のパラメータを測定した。
酸素飽和度の空間分布ＳＯ２（ｒ）は式（２）に示すように、酸化ヘモグロビンＨｂＯ２と脱酸素ヘモグロビンＨｂとの比率に基づいて算出することができる。
【数2】
【0066】
＜４＞統計解析
上記で得られた種々のパラメータ、被験者の年齢について、ＪＭＰソフトウェアバージョン１４．０（米国ＳＡＳ Ｉｎｓｔｉｔｕｔｅ）を用いて解析した。各データ間の相関は、ピアソンの相関分析によって評価した。
上記で得られた種々のパラメータ、被験者の年齢について、ＪＭＰソフトウェアバージョン１４．０（米国ＳＡＳ Ｉｎｓｔｉｔｕｔｅ）を用いて解析した。各データ間の相関は、ピアソンの相関分析によって評価した。
【0067】
＜５＞血管構造の解析
図１（Ａ）は、５７５ｎｍの光音響信号から得られた画像と６５０ｎｍの超音波信号から得られた画像をマージした画像である。このマージ画像から、光音響測定によると、皮膚の上部から下部の深さで血管構造を識別することができることがわかった。
図１（Ａ）は、５７５ｎｍの光音響信号から得られた画像と６５０ｎｍの超音波信号から得られた画像をマージした画像である。このマージ画像から、光音響測定によると、皮膚の上部から下部の深さで血管構造を識別することができることがわかった。
【0068】
図１（Ｂ）（Ｃ）は、５７５ｎｍの光音響信号から６５０ｎｍの光音響信号を引いた差分値に基づいて抽出した、血管構造を示す代表的な測定画像である。（Ｂ）が２０代被験者の画像であり、（Ｃ）が６０代被験者の画像である。
特に、上面図を見ると明らかなとおり、図１より、２０代被験者に見られた微小な血管構造が存在する領域（図２（Ｂ）中、白い三角形で示した血管）が、６０代被験者では見られない（図２（Ｃ）中白い矢印で示した部分）ことがわかる。
特に、上面図を見ると明らかなとおり、図１より、２０代被験者に見られた微小な血管構造が存在する領域（図２（Ｂ）中、白い三角形で示した血管）が、６０代被験者では見られない（図２（Ｃ）中白い矢印で示した部分）ことがわかる。
【0069】
図２は、前述の深層学習を用いて構築したネットワークにより、光音響測定により得られた血管構造を含むｘ－ｙ平面の測定画像データと、これをネットワークに入力して得られた、血管構造の予測画像データの代表例である。測定画像データからデータキュレーターにより視覚的に抽出された血管構造の注釈画像データを併せて示す。（Ａ）は、深層学習において訓練データとして用いられた測定画像データから得られた予測画像であり、（Ｂ）は、深層学習において訓練データとして用いられていない測定画像データから得られた予測画像である。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
【0070】
（Ｃ）は、予測画像データ（ｘ－ｙ平面）をｚ軸方向に結合することで生成した三次元データの代表例である。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、いくつかのパラメータと年齢との間に相関関係が確認された。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、いくつかのパラメータと年齢との間に相関関係が確認された。
【0071】
図３に示される通り、皮膚の各層の何れにおいても、血管構造の平均厚さ（Ａｖｅｒａｇｅ ｔｈｉｃｋｎｅｓｓ）は、年齢とともに増加することが観察された。
また、年齢の増加とともに、上部における血管構造の体積が減少し、下部における血管構造の体積が増加することが観察された（データは図示しない）。
また、年齢の増加とともに分岐の数は、減少することが観察された（データは図示しない）。
すなわち、加齢に伴い、皮膚上部においては、微小血管が減少し、血管構造全体の体積としても減少するのに対し、皮膚下部においては、微小血管は減少するものの血管構造全体の体積は減少しないことがわかり、加齢に伴う血管構造の変化は、皮膚の層（上部、中間部及び下部）によって異なることが示唆された。
また、年齢の増加とともに、上部における血管構造の体積が減少し、下部における血管構造の体積が増加することが観察された（データは図示しない）。
また、年齢の増加とともに分岐の数は、減少することが観察された（データは図示しない）。
すなわち、加齢に伴い、皮膚上部においては、微小血管が減少し、血管構造全体の体積としても減少するのに対し、皮膚下部においては、微小血管は減少するものの血管構造全体の体積は減少しないことがわかり、加齢に伴う血管構造の変化は、皮膚の層（上部、中間部及び下部）によって異なることが示唆された。
【0072】
以上より、微小血管の量を示す特徴量については、加齢状態の評価のための指標として用いることができることが明らかとなった。
そして、その特徴量としては、血管構造の平均厚み、分岐の数を用いることができることが明らかとなった。
また、皮膚上部における血管構造の体積についても、加齢状態の評価のための指標として用いることができることが明らかとなった。
また、皮膚上部における血管構造の平均厚み、分岐の数、血管構造の体積の任意の組み合わせを用いて、加齢状態の評価のための指標として用いることができることが明らかとなった。
そして、その特徴量としては、血管構造の平均厚み、分岐の数を用いることができることが明らかとなった。
また、皮膚上部における血管構造の体積についても、加齢状態の評価のための指標として用いることができることが明らかとなった。
また、皮膚上部における血管構造の平均厚み、分岐の数、血管構造の体積の任意の組み合わせを用いて、加齢状態の評価のための指標として用いることができることが明らかとなった。
【0073】
＜６＞酸素飽和度の解析
図４に、２０代被験者及び６０代被験者の光音響測定画像データを示す。光音響イメージングにより、酸素飽和度の空間分布が可視化できることが明らかとなった。
図４に、２０代被験者及び６０代被験者の光音響測定画像データを示す。光音響イメージングにより、酸素飽和度の空間分布が可視化できることが明らかとなった。
【0074】
図４中、（Ａ）が２０代被験者の血管における酸素飽和度の空間分布を示す測定画像データであり、（Ｂ）が６０代被験者の血管における酸素飽和度の空間分布を示す測定画像データである。
図３において、（Ａ）と（Ｂ）とを比較すると、６０代被験者の測定領域における酸素飽和度は、２０代被験者の測定領域における血管の酸素飽和度に比べて、全体的に低くなっていた。また、２０代被験者の測定画像データにおいては、微細血管が高い酸素飽和度で存在する領域（（Ａ）中、三角形で示す）が観察されるが、６０代被験者の測定画像データにおいては、当該領域が存在しなかった（（Ｂ）中、矢印で示す）。
図３において、（Ａ）と（Ｂ）とを比較すると、６０代被験者の測定領域における酸素飽和度は、２０代被験者の測定領域における血管の酸素飽和度に比べて、全体的に低くなっていた。また、２０代被験者の測定画像データにおいては、微細血管が高い酸素飽和度で存在する領域（（Ａ）中、三角形で示す）が観察されるが、６０代被験者の測定画像データにおいては、当該領域が存在しなかった（（Ｂ）中、矢印で示す）。
【0075】
統計解析に供するために、前述した方法により酸素飽和度を算出し、先に測定した血管構造に関する特徴量、および年齢との相関関係を解析した。
【0076】
図５に示す通り、加齢に伴い、皮膚の各層における酸素飽和度（Ｏｘｙｇｅｎ ｓａｔｕｒａｔｉｏｎ）は減少する傾向があることが明らかになった。
すなわち、酸素飽和度は、肌の加齢状態の評価指標として使用できることが明らかになった。
すなわち、酸素飽和度は、肌の加齢状態の評価指標として使用できることが明らかになった。
【0077】
また、図６に示す通り、皮膚の上部と下部においては、血管の平均厚さと酸素飽和度の間に負の相関があることが明らかになった。
【0078】
したがって、血管の平均厚さから酸素飽和度を推定したり、反対に酸素飽和度の測定によって、血管の平均厚さを推定したりすることができる。
またこの解析により、微小な血管構造が加齢に伴って消失することで、酸素飽和度の低下を誘発することが示唆された。また、この現象は、少なくとも皮膚の上部及び下部で起こることが分かった。
またこの解析により、微小な血管構造が加齢に伴って消失することで、酸素飽和度の低下を誘発することが示唆された。また、この現象は、少なくとも皮膚の上部及び下部で起こることが分かった。
【0079】
以上より、測定領域における酸素飽和度の分布を示す特徴量については、加齢状態の評価のための指標として用いることができることが明らかとなった。
そして、その特徴量としては、酸素飽和度、高酸素飽和度の微細血管領域の体積を用いることができることが明らかとなった。
特に、皮膚の上部、及び下部における酸素飽和度は、加齢状態と関係のある血管の厚さの平均とも関係を有していることから、これらの層における酸素飽和度と血管構造に関する特徴量を組み合わせて評価に用いることで、より精度の高い評価ができることが示唆された。
そして、その特徴量としては、酸素飽和度、高酸素飽和度の微細血管領域の体積を用いることができることが明らかとなった。
特に、皮膚の上部、及び下部における酸素飽和度は、加齢状態と関係のある血管の厚さの平均とも関係を有していることから、これらの層における酸素飽和度と血管構造に関する特徴量を組み合わせて評価に用いることで、より精度の高い評価ができることが示唆された。
【0080】
＜７＞メラニン構造の解析
図７に、波長６５０ｎｍのパルス光を用いた光音響法により取得した測定画像データを示す。
図７に、波長６５０ｎｍのパルス光を用いた光音響法により取得した測定画像データを示す。
【0081】
図７中、（Ａ）が２０代被験者のメラニン構造の分布を示す測定画像データであり、（Ｂ）が６０代被験者のメラニン構造の分布を示す測定画像データである。
（Ａ）に示す通り、２０代被験者の肌においては、メラニン構造がそれほど豊富ではなかった。
これに対し、（Ｂ）に示す通り、６０代被験者の肌においては、トロイダル空洞を伴う、チューブ状の紋理構造（矢印で示す）が観察された。
（Ａ）に示す通り、２０代被験者の肌においては、メラニン構造がそれほど豊富ではなかった。
これに対し、（Ｂ）に示す通り、６０代被験者の肌においては、トロイダル空洞を伴う、チューブ状の紋理構造（矢印で示す）が観察された。
【0082】
図８は、前述の深層学習を用いて構築したネットワークにより、光音響測定により得られたメラニン構造を含むｘ－ｙ平面の測定画像データと、これをネットワークに入力して得られた、メラニン構造の予測画像データの代表例である。測定画像データからデータキュレーターにより視覚的に抽出されたメラニン構造の注釈画像データを併せて示す。（Ａ）は、深層学習において訓練データとして用いられた測定画像データから得られた予測画像データであり、（Ｂ）は、深層学習において訓練データとして用いられていない測定画像データから得られた予測画像データである。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
これより、学習データとして用いられていない測定画像データからも、解析者が作成した正解データと高い類似性をもって予測画像データを生成することができることが確認された。
【0083】
（Ｃ）は、予測画像データ（ｘ－ｙ平面）をｚ軸方向に結合することで生成した三次元データの代表例である。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、メラニンの管状パターン構造の体積と年齢との間に相関関係が確認された。
この三次元データから抽出された各種パラメータ、被験者の年齢の相関解析の結果、メラニンの管状パターン構造の体積と年齢との間に相関関係が確認された。
【0084】
図９に示すとおり、加齢に伴って、皮膚の表層部に存在するメラニンのチューブ状の紋理構造の体積（ｖｏｌｕｍｅ ｏｆ ｍｅｌａｎｉｎ ｓｔｒｕｃｔｕｒｅ）が増大することが明らかになった。
従来、正常皮膚におけるメラニンの三次元分布は研究されていない。また、本研究で実証された１０～数百μｍのスケールのメラニン構造についても報告はない。
この新たに発見された特徴的なメラニン構造に関する特徴量は、加齢状態の評価に有用であることが明らかになった。
従来、正常皮膚におけるメラニンの三次元分布は研究されていない。また、本研究で実証された１０～数百μｍのスケールのメラニン構造についても報告はない。
この新たに発見された特徴的なメラニン構造に関する特徴量は、加齢状態の評価に有用であることが明らかになった。
【0085】
以上より、表層におけるメラニンの体積、特にメラニンのチューブ状の紋理構造の体積を評価指標として、加齢状態を評価することができることが明らかとなった。
【0086】
また、前述した酸素飽和度と、メラニンの管状パターン構造の体積との相関関係を解析した。
図１０に示す通り、酸素飽和度が減少するのに伴って、メラニンのチューブ状の紋理構造の体積が増加することが明らかになった。
したがって、メラニンのチューブ状の紋理構造の体積に基づいて、酸素飽和度の空間分布を推定できる。
また反対に、酸素飽和度に基づいて、メラニンのチューブ状の紋理構造の体積を推定できる。
図１０に示す通り、酸素飽和度が減少するのに伴って、メラニンのチューブ状の紋理構造の体積が増加することが明らかになった。
したがって、メラニンのチューブ状の紋理構造の体積に基づいて、酸素飽和度の空間分布を推定できる。
また反対に、酸素飽和度に基づいて、メラニンのチューブ状の紋理構造の体積を推定できる。
【産業上の利用可能性】
【0087】
本発明は、美容に関する肌の評価方法に利用することができる。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
