JP6810619B2 - 肌状態評価方法 - Google Patents

Description

本発明は、肌状態を評価する技術に関して、特に、肌の凹凸形状に関する波形情報を分析する技術に関する。

皮膚表面の撮像情報に基づいて皮膚表面の状態を評価する手法が幾つか存在している。
特許文献１には、皮膚表面の三次元形状データに基づいて皮膚表面に関する原画像を生成し、その原画像をガウシアンフィルタ処理及び／又はノイズ除去処理を用いて補正し、その補正画像に基づいて抽出される皮丘又は皮溝に基づいて皮膚表面の状態を評価する手法が提案されている。この手法では、共焦点顕微鏡などによりレプリカ又は皮膚表面の画像データが当該三次元形状データとして取得されている。
特許文献２には、肌の表面凹凸形状の情報と肌の光学特性の情報とに基づいて、肌のつやを評価する手法が提案されている。この手法では、肌の表面をマイクロスコープで撮影することで得られる画像から肌の表面凹凸形状が取得されている。
非特許文献１では、皮膚表面から採取した皮表レプリカを顕微鏡に接続したテレビカメラで撮影し、得られた画像データに対してフーリエ変換を適用して得られるパワースペクトル分布を方向別又は周波数別に集計して皮表形態の特徴を数量化する手法が報告されている。

特開２０１３−１８８３２５号公報 特開２０１６−１８９８７７号公報

林 照次（他２名），"加齢及び化粧品による皮表形態の変化−画像解析及びフーリエ変換を応用した新しい方法による３次元解析−"，Ｊ．Ｓｏｃ．Ｃｏｓｍｅｔ．Ｃｈｅｍ．Ｊａｐａｎ．Ｖｏｌ．２３，Ｎｏ．１，１９８９，ｐ．４３−ｐ．５４

上述の提案手法は、皮膚表面の画像データに対して周波数解析、フィルタ処理などを施すことで、皮膚表面の状態の定量化を試みている。
ところが、皮膚表面は、スケールの異なる複数種の凹凸構造から形成されており、上述の提案手法では、このような複数種の凹凸構造を的確に評価することは難しい。例えば、皮膚表面形状の空間均一性といった観点での定量評価は困難である。

本発明は、このような課題に鑑みてなされたものであり、肌の凹凸形状に関する空間波形情報を用いて的確に肌状態を定量評価する技術を提供する。
本明細書において「肌」は、「皮膚」又は「皮膚表面」と意味的に区別せずに用いられる。

本発明の各態様では、上述した課題を解決するために、それぞれ以下の構成を採用する。

第一の態様は、肌状態評価方法に関する。第一の態様に係る肌状態評価方法は、肌の凹凸形状に関する空間波形情報を取得する工程と、前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を導出する導出工程と、前記導出された特性情報に基づいて、前記肌の状態を評価する評価工程と、を含む。

第二の態様は、肌状態評価装置に関する。第二の態様に係る肌状態評価装置は、肌の凹凸形状に関する空間波形情報を取得する取得手段と、前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を生成する生成手段と、前記生成された特性情報又は該特性情報に基づく前記肌の状態の評価情報を出力する出力処理手段と、を備える。

なお、本発明の別態様としては、上記第一態様に係る肌状態評価方法をコンピュータに実行させるプログラム、このようなプログラムを記録したコンピュータが読み取り可能な記録媒体があり得る。この記録媒体は、非一時的な有形の媒体を含む。

上記各態様によれば、肌の凹凸形状に関する空間波形情報を用いて的確に肌状態を定量評価する技術を提供することができる。

本実施形態に係る肌状態評価方法を示すフローチャートである。 第一実施形態に係る肌状態評価方法を示すフローチャートである。 第一実施形態におけるローレンツプロット及び座標変換を概念的に示す図である。 直線波形のローレンツプロット及び座標変換を示す図である。 三角波形及び三角波の変形波形のローレンツプロット及び座標変換を示す図である。 二種の指標値を座標軸とするグラフの例を示す図である。 第二実施形態に係る肌状態評価方法を示すフローチャートである。 第二実施形態におけるベクトル間の内積演算を概念的に示す図である。 第三実施形態に係る肌状態評価方法を示すフローチャートである。 工程（Ｓ４５）における類似度の算出例を概念的に示す図である。 第四実施形態に係る肌状態評価方法を示すフローチャートである。 第四実施形態におけるベクトル間の内積演算を概念的に示す図である。 それぞれ異なる空間周波数帯を選択したウェーブレット変換適用後の変換波形データ（凹凸パワー波形データ）の例を示す図である。 肌状態評価装置のハードウェア構成例を概念的に示す図である。 肌状態評価装置の処理構成例を概念的に示す図である。 皮膚外用剤適用前後の肌に関してローレンツプロットで得られるグラフを示す図である。 図１６に示される１６８Ｈｚから１２３８Ｈｚの周波数帯に対応する変換波形データ（凹凸パワー波形データ）から第一実施形態の手法で得られた５種の指標値を示す図である。 皮膚外用剤適用前後の同一被験者の肌に関する凹凸パワー波形データから第二実施形態の手法で得られたベクトル間の内積演算結果を示す図である。 被験者の肌と比較対象肌との類似性と第三実施形態の手法で導出された指標値間の類似度との関係を示す図である。 被験者の肌と比較対象肌との類似性と第四実施形態の手法で導出されたベクトル間の類似度との関係を示す図である。

以下、本発明の実施の形態について説明する。なお、以下に挙げる各実施形態はそれぞれ例示であり、本発明は以下の各実施形態の構成に限定されない。

まず、本実施形態に係る肌状態評価方法の概要について図１を用いて説明する。
図１は、本実施形態に係る肌状態評価方法を示すフローチャートである。
本実施形態に係る肌状態評価方法は、少なくとも、図１に示されるように工程（Ｓ１１）、工程（Ｓ１３）、及び工程（Ｓ１５）を含む。

工程（Ｓ１１）は、肌の凹凸形状に関する空間波形情報を取得する工程である。
ここで対象となる肌は、状態を評価したい被験者の任意の部位の肌の任意の領域である。
「肌の凹凸形状に関する空間波形情報」とは、肌の凹凸形状に関する空間変化（空間変動）を示す波形情報である。
「波形情報」又は「波形データ」とは、空間軸上に並ぶ振幅値の列データである。
「肌の凹凸形状に関する空間波形情報」は、高さの空間変化を示す凹凸波形データ（空間軸上に並ぶ高さ情報の列データ）であってもよいし、凹凸波形データから周波数解析などの演算により得られる値（パワー値などの振幅値）の空間変化を示す波形データであってもよい。
以降、前者の凹凸波形データをベース波形データと表記する場合がある。例えば、個々の画素値が皮膚表面の高さを示す画像データは凹凸波形データに該当する。
後者の例としては、凹凸波形データ（ベース波形データ）に対してウェーブレット変換といった空間情報を残す周波数解析を適用することにより、特定周波数成分の混合強度（パワー）の空間変化を示す波形データや、ベース波形データに対してフーリエ変換を適用することで特定周波数成分に限定した波形データがあり得る。このように、ベース波形データに対して何らかの演算を施された波形データを変換波形データと表記する場合がある。

工程（Ｓ１１）では、後述の肌状態評価装置のようなコンピュータが、レーザ顕微鏡などの肌の凹凸形状を測定した測定装置から凹凸波形データを当該空間波形情報として取得することができる。コンピュータは、他のコンピュータ又は可搬型記録媒体から凹凸波形データを取得してもよい。
また、当該コンピュータが、その取得された凹凸波形データに対して所定の演算を施すことにより、変換波形データを当該空間波形情報として取得することもできる。コンピュータは、他のコンピュータ又は可搬型記録媒体から当該変換波形データを取得してもよい。
ところで、肌の凹凸形状の測定装置は、肌から直接、高さ情報を測定してもよいし、肌から既存の手法で得られるレプリカ（非特許文献１参照）を用いて間接的に肌の高さ情報を測定してもよい。肌の凹凸形状の測定装置の測定原理は制限されない。

工程（Ｓ１３）は、工程（Ｓ１１）で取得された空間波形情報を用いて、対象空間波形の空間変動特性又はその対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報、より具体的には、対象空間波形の周期乱れを示す特性情報を導出する工程である。
ここで「対象空間波形の空間変動特性を示す特性情報」とは、対象となる空間波形の空間軸上での変動傾向を示す定量情報を意味する。また、「周期乱れを示す特性情報」とは、空間波形の乱れ、すなわち周波数や位相の乱れ（揺らぎ）を振幅から分離して表現した情報であり、本実施形態ではこれとは別に振幅の乱れ（揺らぎ）も同時に特性情報としている。この特性情報は、例えば、対象となる空間波形と、その空間波形を空間軸上で所定距離分シフトしたシフト波形との関係から得ることができる。この空間変動特性を示す特性情報の具体的な導出手法については後述する。

また、「対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報」とは、対象空間波形と比較対象肌の空間波形との関係性を示す定量情報を意味する。この特性情報は、例えば、両波形の類似度又は乖離度などで示され得る。この相互関係特性を示す特性情報の具体的な導出手法については後述する。

ここで、比較対象肌は、対象肌と同一人物の肌であっても、他人の肌であってもよい。例えば、比較対象肌は、或る美容施術又は化粧料を施す前の同一人物の肌である。この例によれば、美容施術又は化粧料を施す前後での被験者の肌を相対的に評価することができる。他の例として、同一人物の或る年齢時の肌が比較対象肌とされてもよい。この例によれば、その年齢時から今までの被験者の肌の変化を相対評価することができる。他の例として、同一人物の他の部位の肌が比較対象肌とされてもよい。この例によれば、被験者自身の肌を部位間で相対的に評価することができる。更に他の例として、理想とされる肌或いは被験者が憧れる人の肌が比較対象肌とされてもよい。また、対象肌と比較対象肌とでは同一部位が対象とされてもよいし、異なる部位が対象とされてもよい。

比較対象肌の空間波形情報は、対象肌に関する当該空間波形情報を取得する工程（Ｓ１１）の実行と同期して（並列又はその前後にて）、取得されてもよいし、対象肌に関する空間波形情報とは別時期に別途取得されてもよい。更に言えば、比較対象肌に関する空間波形情報は、対象肌に関するそれとは異なる方法で予め取得され、保持されていてもよい。工程（Ｓ１３）において、コンピュータが、他のコンピュータ又は可搬型記録媒体などから比較対象肌に関する空間波形情報を取得してもよい。

工程（Ｓ１５）は、工程（Ｓ１３）で導出された特性情報に基づいて、対象肌の状態を評価する工程である。
工程（Ｓ１５）で評価される肌の状態とは、肌の凹凸形状に関わる状態である。本実施形態では、肌の凹凸形状を形成する複数種の凹凸構造の空間均一性も肌状態として評価することができるため、評価される肌の状態としては、整った状態か否か、整った状態若しくは不揃い状態の程度、肌年齢など様々な指標で評価することができる。

工程（Ｓ１５）における具体的な評価手法については後述する。但し、この評価手法は、後述する例に限定されない。
例えば、コンピュータが、工程（Ｓ１３）で導出された特性情報に基づいて、所定のアルゴリズムにより対象肌の状態を評価してもよい。この場合、コンピュータは、評価結果（「最高レベル」、「良好」、「改善要」など）を提示することもできる。また、コンピュータが工程（Ｓ１３）で導出された特性情報を提示し、この特性情報を参照した者（被験者や評価者など）が、その特性情報を考慮して対象肌を評価してもよい。

概要として説明したように、本実施形態では、対象肌の凹凸形状に関する空間波形情報から、対象空間波形の空間変動特性を示す特性情報又は対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報が導出され、その特性情報に基づいて対象肌の状態が評価される。
これにより、本実施形態によれば、肌の凹凸形状に関わる肌状態を定量的に評価することができる。
更に、本実施形態において肌状態の評価の根拠とされる特性情報は、対象肌の凹凸形状の空間変動の情報を含む定量情報である。従って、本実施形態によれば、肌の凹凸形状の空間均一性といった観点でも肌状態を的確に定量評価することができる。

以下、本実施形態に係る肌状態評価方法についてより具体的に説明する。
［第一実施形態］
第一実施形態では、上述の空間変動特性を示す特性情報に基づく肌状態評価の一具体例を挙げる。以下、第一実施形態に係る肌状態評価方法を本評価方法と略称する場合もある。
図２は、第一実施形態に係る肌状態評価方法を示すフローチャートである。
本評価方法は、図２に示されるように、工程（Ｓ２１）から工程（Ｓ２５）を含む。工程（Ｓ２１）は、図１に示される工程（Ｓ１１）に対応し、工程（Ｓ２２）から工程（Ｓ２４）は、工程（Ｓ１３）に対応し、工程（Ｓ２５）は、工程（Ｓ１５）に対応する。
本評価方法は、例えば、後述の肌状態評価装置により実行される。但し、本評価方法は、各工程において人手が介入することを許容する。

工程（Ｓ２１）は、上述の工程（Ｓ１１）と同様であり、被験者の肌の凹凸形状に関する空間波形情報を取得する工程である。工程（Ｓ２１）では、上述したとおり、凹凸波形データを当該空間波形情報として取得してもよいし、その凹凸波形データから得られる変換波形データを当該空間波形情報として取得してもよい。また、工程（Ｓ２１）における当該空間波形情報の取得手法についても上述したとおりである。

工程（Ｓ２２）は、工程（Ｓ２１）で取得された空間波形情報を、対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換する工程である。
即ち、ここでの直交座標系は、空間波形情報により示される対象空間波形における空間軸上で所定距離離れた振幅値ペアの一方を横軸（ｘ軸）、他方を縦軸（ｙ軸）にした直交平面座標系となる。このように空間波形情報を直交座標系の座標点群に変換することをローレンツプロットと表記する場合もある。
ローレンツプロットに用いる当該所定距離は、シフト距離であり、抽出したい空間周波数成分に対応する適切な値に設定されることが望ましい。例えば、細かな空間変化を抽出したい場合には、所定距離Δｘは、大まかな空間変化を抽出したい場合に比べて、小さい値に設定すればよい。

工程（Ｓ２３）は、工程（Ｓ２２）で得られた直交座標系の座標点群を、当該直交座標系の直交軸であるｘ軸及びｙ軸で示される「ｙ＝ｘ」軸及び「ｙ＝−ｘ」軸に投影した変換座標点群に変換する工程である。

図３は、第一実施形態におけるローレンツプロット及び座標変換を概念的に示す図である。
図３で示される凹凸パワー波形データは、ベース波形データである凹凸波形データに対してウェーブレット変換を適用して得られる変換波形データである。この変換波形データが空間波形情報として取得され、この変換波形データがローレンツプロットにより直交座標系の座標点群に変換される。ローレンツプロットでは、変換波形データの振幅値（パワー値）が横軸（ｘ軸）に、所定距離Δｘシフトされたシフト波形データの振幅値（パワー値）が縦軸（ｙ軸）に取られる。つまり、当該直交座標系の一座標点のｙ軸要素である振幅値は、変換波形データにおいて、同座標点のｘ軸要素である振幅値から所定距離Δｘ離間した位置に現れる値である。
座標変換では、当該直交座標系の座標点群が「ｙ＝ｘ」軸及び「ｙ＝−ｘ」軸に投影された座標に変換される。図３では、「ｙ＝ｘ」軸が変換後の座標系においてＸ'軸と表記され、「ｙ＝−ｘ」軸が変換後の座標系においてＹ'軸と表記されている。

工程（Ｓ２４）は、工程（Ｓ２３）で得られた変換座標点群の分布状態を示す一種以上の指標値を当該対象空間波形の空間変動特性を示す特性情報として導出する工程である。
例えば、工程（Ｓ２４）では、当該変換座標点群における、「ｙ＝ｘ」軸（図３のＸ'軸）における原点からの距離の平均値若しくは標準偏差、「ｙ＝−ｘ」軸（図３のＹ'軸）における原点からの距離の平均値若しくは標準偏差、又は、それら平均値若しくはそれら標準偏差を用いて算出される比率若しくは分布面積の一つ以上を当該一種以上の指標値として導出することができる。

図３では、Ｘ'軸における原点からの距離の平均値が黒塗り星マークで示され、その距離の標準偏差がσｘで示され、Ｙ'軸における原点からの距離の平均値が白抜き星マークで示され、その距離の標準偏差がσｙで示されている。
比率としては、Ｘ'軸における平均距離とＹ'軸における平均距離との比率、又は、Ｘ'軸における標準偏差σｘとＹ'軸における標準偏差σｙとの比率（σｘ／σｙ）が算出され得る。
分布面積は、標準偏差σｘ及び標準偏差σｙを用いて、楕円面積として（π×σｘ×σｙ）で算出され得る。
但し、算出される一種以上の指標値は、このような例に限定されず、変換座標点群の分布状態を示す他の指標値であってもよい。例えば、距離の分散値又は分散値の比率が指標値として用いられてもよいし、分布面積は、長方形面積として（２σｘ×２σｙ）で算出されてもよい。工程（Ｓ２４）で算出される一種以上の指標値は、工程（Ｓ２３）で得られた変換座標点群の分布状態を示す複数種の指標値のうち、肌状態と相関の高い指標値とされることが望ましい。

本発明者らは、このようにローレンツプロットを用いて得られる指標値が、肌の凹凸形状に関する空間波形情報の特徴を示していることを見出している。
図４は、直線波形のローレンツプロット及び座標変換を示す図である。図４に示されるように、直線波形は、ローレンツプロット及び座標変換により一座標点に変換される。直線波形の振幅を増加させると、座標変換後の座標点はＸ'軸上を移動する。このことから、Ｘ'軸における平均距離は空間波形情報の振幅値（高さ、パワーなど）の平均情報を示すといえる。

図５は、三角波形及び三角波の変形波形のローレンツプロット及び座標変換を示す図である。図５に示されるように、三角波形はローレンツプロット及び座標変換によりＸ'軸に平行な直線状の複数座標点に変換される。この三角波形の振幅を増加させると、座標変換後の複数座標点の直線状分布はＸ'軸方向に延びる。このことから、Ｘ'軸における標準偏差σｘは空間波形情報の振幅（高さ、パワーなど）の幅を示すといえる。
また、三角波形に周期の乱れを含めると、座標変換後の複数座標点の分布がＹ'軸方向に延びる。更に、三角波形を不均一化すると、座標変換後の複数座標点の分布が散らばる。このことから、Ｙ'軸における平均距離は空間波形情報の周期の乱れの平均情報を示し、Ｙ'軸における標準偏差σｙは空間波形情報の周期の乱れの大きさを示すといえる。
更に、分布面積は、標準偏差σｘ及びσｙを用いるため、空間波形情報の振幅の幅及び周期の乱れの大きさを示すといえ、標準偏差の比率（σｘ／σｙ）は、空間波形情報の振幅の幅と周期の乱れの大きさとの比率を示すといえる。
なお、図５で示されるローレンツプロットでは所定距離Δｘが１周期に設定されている。
これらにより、ローレンツプロットを用いて得られる指標値が、肌の凹凸形状に関する空間波形情報の特徴を示すことが理論的に実証される。更に、本発明者らにより、これら指標値により人の肌の状態を定量評価できることは実証されている（実施例の項参照）。

工程（Ｓ２５）は、工程（Ｓ２４）で導出された一種以上の指標値を提示する工程である。
指標値の提示は、ディスプレイのような表示装置への表示で実現されてもよいし、プリンタへの印刷で実現されてもよいし、音声出力で実現されてもよい。
本評価方法では、被験者自ら又は評価者が、工程（Ｓ２５）で提示された指標値を参照して、被験者の肌状態を評価することができる。このとき、工程（Ｓ２５）では、被験者の肌に関する指標値と比較するための基準となる指標値又は肌状態の空間均一性の良し悪しの境界を示す閾値を更に提示することもできる。また、このような比較のための指標値又は閾値は、予め印刷媒体などに出力されていてもよい。

また、本評価方法では、工程（Ｓ２４）で導出された一種以上の指標値をそのような基準となる指標値や閾値などと比較することにより、被験者の肌状態を自動的に評価し、その評価結果を工程（Ｓ２５）で提示することもできる。例えば、指標種ごとに閾値を予め保持しておき、各種の指標値が閾値を超えるか否かにより、最高、良好、改善要、劣悪などと肌状態を評価し、その評価結果を提示することができる。

また、本評価方法では、工程（Ｓ２４）で複数種の指標値が導出される場合に、工程（Ｓ２５）において、各指標種を座標軸とする座標系に、被験者の肌に関する複数種の指標値に対応する点がプロットされたグラフを提示することもできる。
図６は、二種の指標値を座標軸とするグラフの例を示す図である。図６に示されるグラフでは、Ｙ'軸の平均距離が横軸（ｘ軸）とされ、Ｙ'軸の標準偏差σｙが縦軸（ｙ軸）とされている。空間波形に関するＹ'軸の平均距離は周期の乱れの平均情報を示し、Ｙ'軸の標準偏差は周期の乱れの大きさを示す。そのため、図６のグラフによれば、肌の凹凸構造の空間均一性に関する肌状態を視覚的に評価することができる。
例えば、図６に示されるように、被験者の肌（被験肌）に関する指標値と共に、同年代の平均的な肌状態や理想的な肌状態のような基準となる肌状態に関する指標値（基準指標値）がプロットされていてもよい。これにより、グラフ上におけるこの基準指標値との位置関係に応じて、被験者の肌状態を容易に評価することができる。

このように第一実施形態では、被験者の肌の凹凸形状に関する空間波形情報がローレンツプロットにより直交座標系の座標点群に変換され、更に座標変換される。そして、得られた変換座標点群の分布状態を示す一種以上の指標値が導出され、この指標値が被験者の肌状態を示す定量情報として提示される。この指標値が人の肌状態を示すことは本発明者らにより実証されている（実施例参照）。
従って、第一実施形態によれば、肌の凹凸形状の空間変動に関わる肌状態、例えば、凹凸の空間均一性といった肌状態を的確に定量評価することができる。

［第一実施形態の変形例］
上述の第一実施形態では、被験者の肌に関して導出される一種以上の指標値に基づいて、肌状態が評価されたが、指標値を用いることなく肌状態を評価することはできる。
例えば、図２の工程（Ｓ２４）の代わりに、工程（Ｓ２３）で得られた変換座標点群がプロットされたグラフを当該空間変動特性を示す特性情報として導出し、工程（Ｓ２５）においてこのグラフを提示してもよい。このグラフは、例えば、図３の下部に示される座標変換後のグラフである。
また、ローレンツプロットされたグラフを当該空間変動特性を示す特性情報として導出し、このグラフを提示してもよい。即ち、図２の工程（Ｓ２３）を省き、工程（Ｓ２４）の代わりに、工程（Ｓ２２）で得られた直交座標系の座標点群がプロットされたグラフを当該特性情報として導出する。この場合、図３の中央に示されるローレンツプロットされたグラフが提示される。
当該直交座標系の座標点群又はこの座標点群に対する座標変換により得られる変換座標点群がプロットされたグラフによれば、グラフの形状から視覚的にその指標値と同様の情報を把握することができる。即ち、このグラフによっても、被験者の凹凸形状に関わる肌状態を的確に定量評価することができる。

［第二実施形態］
第二実施形態では、上述の空間変動特性を示す特性情報に基づく肌状態評価の一具体例として、第一実施形態と異なる例を挙げる。以下、第二実施形態に係る肌状態評価方法を本評価方法と略称する場合もある。

図７は、第二実施形態に係る肌状態評価方法を示すフローチャートである。
本評価方法は、図７に示されるように、工程（Ｓ３１）から工程（Ｓ３４）を含む。工程（Ｓ３１）は図１に示される工程（Ｓ１１）に対応し、工程（Ｓ３２）及び工程（Ｓ３３）は工程（Ｓ１３）に対応し、工程（Ｓ３４）は工程（Ｓ１５）に対応する。
本評価方法は、例えば、後述の肌状態評価装置により実行される。但し、本評価方法は、各工程において人手が介入することを許容する。

工程（Ｓ３１）は、上述の工程（Ｓ１１）及び工程（Ｓ２１）と同様であり、被験者の肌の凹凸形状に関する空間波形情報を取得する工程である。
工程（Ｓ３２）は、工程（Ｓ３１）で取得された空間波形情報を空間軸上で所定距離シフトしたシフト波形情報を生成する工程である。
ここでの所定距離もローレンツプロットに用いる当該所定距離と同様にシフト距離であり、抽出したい周波数成分（凹凸構造のスケール）に対応する適切な値に設定されることが望ましい。例えば、細かな空間変化（小さい凹凸構造）を抽出したい場合には、所定距離Δｘは、大まかな空間変化（大きい凹凸構造）を抽出したい場合に比べて、小さい値に設定すればよい。

工程（Ｓ３３）は、工程（Ｓ３１）で取得された空間波形情報を成分とする元ベクトルと工程（Ｓ３２）で生成されたシフト波形情報を成分とするシフトベクトルとの演算により、対象空間波形の空間変動特性を示す特性情報を導出する工程である。
この演算としては、内積演算、ベクトル間のユークリッド距離演算などがある。内積演算によれば、内積値及びベクトル間の角度ｃｏｓθが当該空間変動特性を示す特性情報として算出され得る。また、ベクトル間のユークリッド距離が当該空間変動特性を示す特性情報として算出されてもよい。なお、当該演算により得られる値は、ベクトル間の類似度と呼ぶこともできる。

図８は、第二実施形態におけるベクトル間の内積演算を概念的に示す図である。
図８で示される凹凸パワー波形データは、図３と同様である。このような変換波形データが空間波形情報として取得され、この変換波形データが空間シフトされ、シフト波形データが生成される。例えば、凹凸パワー波形データの中の１００００個のパワー値を成分とする元ベクトルＡとシフト波形データの中の１００００個のパワー値を成分とするシフトベクトルＢとの内積演算が行われる。但し、各ベクトルの要素数（成分数）は、このような数に限定されない。
図８の例では、内積演算により内積値とｃｏｓθとが当該空間変動特性を示す特性情報として導出されている。

工程（Ｓ３４）は、工程（Ｓ３３）での演算結果、即ち、当該空間変動特性を示す特性情報を提示する工程である。提示される特性情報は、内積値又はｃｏｓθのように、一つの値であってもよいし、複数種の値であってもよい。
工程（Ｓ３４）での演算結果の提示手法は、第一実施形態における指標値の提示手法と同様である。
本評価方法では、被験者自ら又は評価者が、工程（Ｓ３４）で提示された演算結果を参照して、被験者の肌状態を評価することができる。このとき、工程（Ｓ３４）では、被験者の肌に関する演算結果と比較するための基準となる値又は肌の空間均一性の良し悪しの境界を示す閾値を更に提示することもできる。また、このような比較のための基準値又は閾値は、予め印刷媒体などに出力されていてもよい。

また、本評価方法では、工程（Ｓ３３）で導出された演算結果をそのような基準値や閾値などと比較することにより、被験者の凹凸形状に関わる肌状態を自動的に評価し、その評価結果を工程（Ｓ３４）で提示することもできる。例えば、閾値を予め保持しておき、演算結果が閾値を超えるか否かにより、最高、良好、改善要、劣悪などの良し悪しの程度や、肌年齢などにより肌状態を評価し、その評価結果を提示することができる。

このように第二実施形態では、被験者の肌の凹凸形状に関する空間波形情報を成分とするベクトルとその空間波形情報を所定距離シフトしたシフト波形情報を成分とするベクトルとの演算により、当該空間変動特性を示す特性情報が導出される。演算結果がベクトル間の類似度が高いことを示す場合、対象空間波形の空間変動が小さいことを示し、その類似度が低い場合、対象空間波形の空間変動が大きいことを示す。そして、その演算結果が被験者の肌状態を示す定量情報として提示される。このように提示される演算結果が人の肌状態を示すことは本発明者らにより実証されている（実施例参照）。
従って、第二実施形態によれば、肌の凹凸形状の空間変動に関わる肌状態、例えば、肌の凹凸構造の空間均一性といった肌状態を的確に定量評価することができる。

［第三実施形態］
上述の第一及び第二実施形態では、被験者の肌の凹凸形状に関する空間波形情報を処理して得られた情報に基づいて、被験者の肌の状態が評価された。第三実施形態では、被験者の肌に関する対象空間波形と比較対象肌に関する空間波形との間の相互関係特性を示す特性情報に基づいて、被験者の肌の状態が相対的に評価される。以下、第三実施形態に係る肌状態評価方法を本評価方法と略称する場合もある。

図９は、第三実施形態に係る肌状態評価方法を示すフローチャートである。
本評価方法は、図９に示されるように、工程（Ｓ４１）から工程（Ｓ４６）を含む。工程（Ｓ４１）は図１に示される工程（Ｓ１１）に対応し、工程（Ｓ４２）から工程（Ｓ４５）は工程（Ｓ１３）に対応し、工程（Ｓ４６）は工程（Ｓ１５）に対応する。
本評価方法は、例えば、後述の肌状態評価装置により実行される。但し、本評価方法は、各工程において人手が介入することを許容する。

工程（Ｓ４１）から工程（Ｓ４４）は、図２に示される工程（Ｓ２１）から工程（Ｓ２４）と同様である。
工程（Ｓ４５）は、被験者の肌に関して工程（Ｓ４４）で導出された一種以上の指標値と比較対象肌に関する一種以上の指標値との類似度を算出する工程である。
工程（Ｓ４５）で算出される類似度は、被験者の肌の対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報の一例である。

一種の指標値が導出される場合、被験者の肌に関する指標値と比較対象肌に関する指標値との差が当該類似度として算出されてもよい。また、複数種の指標値が導出される場合、当該類似度は、複数種の指標値を成分とするベクトル間のユークリッド距離として算出されてもよいし、それらベクトル間の角度として算出されてもよい。後者の場合、被験者の肌に関する複数種の指標値を成分とする特性ベクトルと比較対象肌に関する複数種の指標値を成分とする特性ベクトルとの内積演算を行い求められるｃｏｓθの値を類似度として算出することができる。

図１０は、工程（Ｓ４５）における類似度の算出例を概念的に示す図である。図１０の例では、被験者の肌（被験肌）に関する複数種の指標値を成分とする特性ベクトルＡと比較対象肌に関する複数種の指標値を成分とする特性ベクトルＢとの内積演算を行い求められるｃｏｓθの値が類似度として算出されている。特性ベクトルＡ及びＢの各成分数は、指標値として採用する上述の指標値の数となる。ｃｏｓθの値は、例えば、特性ベクトルＡと特性ベクトルＢとの内積値（Ａ・Ｂ）を特性ベクトルＡの長さ（｜Ａ｜）と特性ベクトルＢの長さ（｜Ｂ｜）との積で除算することで算出可能である。

比較対象肌に関する指標値は、被験者の肌に関する空間波形情報に対する工程（Ｓ４１）から工程（Ｓ４４）の実行と同期して（並列又はその前後にて）、比較対象肌に関する空間波形情報に対する工程（Ｓ４１）から工程（Ｓ４４）の実行により算出されてもよいし、被験者の肌に関する指標値とは別時期に別途算出されてもよい。更に言えば、比較対象肌に関する指標値は、被験者の肌に関するそれとは異なる方法で予め取得され、保持されていてもよい。工程（Ｓ４５）において、コンピュータが、他のコンピュータ又は可搬型記録媒体などから比較対象肌に関する指標値を取得してもよいし、当該指標値をユーザに入力させてもよい。
比較対象肌に関する指標値は、被験者の肌に関して工程（Ｓ４１）で取得された空間波形情報と同じ測定条件で測定された空間波形情報に基づいて算出されることが望ましい。その測定条件としては、肌の凹凸形状を測定する際の測定部の肌上の移動速度などである。測定条件を合わせることで、評価精度を向上させることができる。

工程（Ｓ４６）は、工程（Ｓ４５）で算出された類似度を提示する工程である。
類似度の提示は、ディスプレイのような表示装置への表示で実現されてもよいし、プリンタへの印刷で実現されてもよいし、音声出力で実現されてもよい。
本実施形態では、被験者自ら又は評価者が、工程（Ｓ４６）で提示された類似度を参照して、被験者の肌状態を比較対象肌との間で相対的に評価する。例えば、提示された類似度が高い程、被験者の肌状態は比較対象肌の状態により近いと評価でき、提示された類似度が低い程、被験者の肌状態は比較対象肌の状態からより乖離していると評価できる。
工程（Ｓ４６）では、類似度と共に、肌状態が近いと評価するための閾値と肌状態が遠い（異なる）と評価するための閾値とを提示することもでき、このような閾値は、予め印刷媒体などに出力されていてもよい。
ところで、被験者の肌状態と比較対象肌の状態との類似度合又は乖離（相違）度合が、肌に関する空間波形情報から得られる複数種の指標値間の類似度により定量的に評価できることは、本発明者らにより実証されている（実施例の項参照）。

このように第三実施形態では、第一実施形態と同様にローレンツプロットを用いて、被験者の肌に関する空間波形情報から一種以上の指標値が導出され、被験者の肌に関する一種以上の指標値と比較対象肌に関する一種以上の指標値との類似度が算出される。算出された類似度は、被験者の肌の対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報であり、被験者の肌と比較対象肌との類似度合又は乖離度合を示す定量情報である。
従って、第三実施形態によれば、肌の凹凸形状の空間変動に関わる肌状態、例えば、肌の凹凸構造の空間均一性といった被験者の肌状態を比較対象肌との比較により定量的に評価することができる。

［第三実施形態の変形例］
上述の第三実施形態では、工程（Ｓ４６）において類似度が提示されたが、類似度の提示に代えて、その類似度に対応する評価結果（「似ています」、「あまり似ていません」など）を提示してもよい。この場合、コンピュータが、予め保持される閾値などとその類似度とを比較することで、評価結果を導出し、この評価結果を提示することができる。

また、工程（Ｓ４５）における類似度の算出及び工程（Ｓ４６）における類似度の提示を実行することなく、次のような比較グラフを提示してもよい。具体的には、複数種の指標値が導出される場合、工程（Ｓ４６）では、各指標種を座標軸とする座標系に被験者の肌に関する複数種の指標値及び比較対象肌に関する複数種の指標値に対応する点がプロットされた比較グラフを提示することができる。この場合、コンピュータがその比較グラフを生成し、表示装置やプリンタなどにその比較グラフを出力し、その比較グラフを参照した者（被験者や評価者など）が、そのグラフ内の被験者の肌を示す点と比較対象肌を示す点との位置関係から被験者の肌を評価してもよい。また、コンピュータが、比較グラフの提示と共に、比較グラフ内の点の位置関係に対応する評価結果を提示することもできる。
例えば、この比較グラフは、図６に示されるグラフで、基準指標値の代わりに比較対象肌に関する指標値に対応する点がプロットされていればよい。

更に、工程（Ｓ４５）及び工程（Ｓ４６）に加えて、工程（Ｓ４４）が実行されなくてもよい。この場合、工程（Ｓ４３）で得られた変換座標点群がプロットされたグラフを生成し、比較対象肌に対しても同様の手法でそのグラフを生成する。そして、その両グラフが対比可能な状態で提示すればよい。また、変換座標点群がプロットされたグラフではなく、ローレンツプロットされたグラフを被験者の肌及び比較対象肌に関してそれぞれ生成し、その両グラフが対比可能な状態で提示されればよい。
言い換えれば、被験者の肌の対象空間波形の空間変動特性を示す特性情報として導出されたグラフと、比較対象肌に関して被験者の肌と同様に空間変動特性を示す特性情報として導出されたグラフとが対比可能な状態で示される比較情報を、被験者の肌の対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報として導出することができる。
当該直交座標系の座標点群又はこの座標点群に対する座標変換により得られる変換座標点群がプロットされたグラフによれば、グラフの形状から視覚的にその指標値と同様の情報を把握することができる。即ち、このグラフによっても、被験者の肌状態を比較対象肌との間で高精度に相対評価することができる。

［第四実施形態］
第四実施形態では、被験者の肌に関する対象空間波形と比較対象肌に関する空間波形との間の相互関係特性を示す特性情報に基づく肌状態評価の一具体例として、第三実施形態と異なる例を挙げる。以下、第四実施形態に係る肌状態評価方法を本評価方法と略称する場合もある。

図１１は、第四実施形態に係る肌状態評価方法を示すフローチャートである。
本評価方法は、図１１に示されるように、工程（Ｓ５１）から工程（Ｓ５３）を含む。工程（Ｓ５１）は図１に示される工程（Ｓ１１）に対応し、工程（Ｓ５２）は工程（Ｓ１３）に対応し、工程（Ｓ５３）は工程（Ｓ１５）に対応する。
本評価方法は、例えば、後述の肌状態評価装置により実行される。但し、本評価方法は、各工程において人手が介入することを許容する。

工程（Ｓ５１）は、上述の工程（Ｓ１１）、工程（Ｓ２１）などと同様であり、被験者の肌の凹凸形状に関する空間波形情報を取得する工程である。
工程（Ｓ５２）は、工程（Ｓ５１）で取得された空間波形情報を成分とする被験肌ベクトルと比較対象肌の空間波形情報を成分とする比較肌ベクトルとの類似度を、被験者の肌の対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報として導出する工程である。
当該類似度は、被験肌ベクトルと比較肌ベクトルとの間のユークリッド距離として算出されてもよいし、それらベクトル間の角度として算出されてもよい。後者の場合、両ベクトルの内積演算を行い求められるｃｏｓθの値を類似度として算出することができる。

図１２は、第四実施形態におけるベクトル間の内積演算を概念的に示す図である。
図１２の例では、被験者の肌に関する空間波形情報を成分とする被験肌ベクトルＡと比較対象肌に関する空間波形情報を成分とする比較肌ベクトルＢとの内積演算を行い求められるｃｏｓθの値が類似度として算出されている。この場合、類似度が１に近付く程、類似度が高く、類似度が０に近付く程、類似度が低くなる。被験肌ベクトルＡ及び比較肌ベクトルＢの各成分数（要素数）は、空間波形情報に含まれる振幅値群の一部又は全部とされる。

比較対象肌に関する比較肌ベクトルは、予め取得され、保持されていてもよい。例えば、工程（Ｓ５２）の実行時又はそれより前に、コンピュータが、他のコンピュータ又は可搬型記録媒体などから比較肌ベクトルを取得してもよいし、その比較肌ベクトルをユーザに入力させてもよい。
また、比較対象肌に関する空間波形情報が、被験者の肌に関する空間波形情報と共に（又は前後して）取得されてもよい。
比較対象肌に関する比較肌ベクトルは、被験者の肌に関して工程（Ｓ５１）で取得された凹凸形状に関する空間波形情報と同じ測定条件で測定された空間波形情報から生成されることが望ましい。その測定条件としては、肌の凹凸形状を測定する際の測定部の肌上の移動速度などである。測定条件を合わせることで、評価精度を向上させることができる。

工程（Ｓ５３）は、工程（Ｓ５２）において当該相互関係特性を示す特性情報として導出された類似度を提示する工程である。工程（Ｓ５３）は、第三実施形態における工程（Ｓ４６）と同様である。なお、被験者の肌状態と比較対象肌の状態との類似度合又は乖離（相違）度合が、肌の凹凸形状に関する空間波形情報を成分とするベクトル間の類似度により定量的に評価できることは、本発明者らにより実証されている（実施例の項参照）。

このように第四実施形態では、被験者の肌の凹凸形状に関する空間波形情報を成分とするベクトルと、比較対象肌の凹凸形状に関する空間波形情報を成分とするベクトルとの類似度が算出される。この算出された類似度は、被験者の肌の対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報であり、被験者の肌と比較対象肌との類似度合又は乖離度合を示す定量情報である。
従って、第四実施形態によれば、第三実施形態と同様に、被験者の肌の凹凸形状に関わる肌状態を比較対象肌との比較により定量的に評価することができる。

［変形例］
上述の第一、第二、及び第三実施形態では、ローレンツプロット又は波形情報の空間シフトにおいて、所定距離が用いられた。
これら各実施形態では、それぞれ異なる複数の所定距離を用いてローレンツプロット又は波形情報の空間シフトをそれぞれ行い、所定距離ごとに上述の一種以上の指標値をそれぞれ算出するようにしてもよい。このように所定距離Δｘを変えることで、肌の凹凸形状に関する空間波形において異なる成分の特徴を抽出できる。

第一実施形態及び第二実施形態では、被験者の肌に関して所定距離Δｘごとに算出された一種以上の指標値を用いて、被験者の肌の状態が評価されればよい。例えば、所定距離Δｘごとに指標値と閾値とを比較して、その比較結果を総合的に評価することで、被験者の肌の凹凸形状に関わる肌状態の評価結果とすることができる。
また、第三実施形態では、被験者の肌に関する指標値及び比較対象肌に関する指標値を、所定距離Δｘごとに比較することにより、被験者の肌状態を比較対象肌との間で相対的に評価することができる。例えば、上述の比較グラフが所定距離Δｘごとにそれぞれ提示されてもよい。また、被験者の肌と比較対象肌との指標値間の類似度が所定距離Δｘごとに算出及び提示されてもよい。
このようにすれば、小さい変化や大きい変化など着目すべき空間周波数成分に特化して、被験者の肌状態を高精度に評価することができる。

このような着目すべき空間周波数成分に特化した評価を行うためには、処理対象とする空間波形情報を特定の空間周波数成分に制限した波形情報としてもよい。
図１３は、それぞれ異なる空間周波数帯を選択したウェーブレット変換適用後の変換波形データ（凹凸パワー波形データ）の例を示す図である。図１３に示されるように、それぞれ異なる４つの空間周波数帯が選択された各変換波形データからローレンツプロットにより変換された座標点群はそれぞれ異なる分布を示している。即ち、着目すべき空間周波数成分を選択した空間波形情報を用いることで、肌の凹凸形状に関する空間波形において異なる成分の特徴を抽出できることが分かる。なお、図１３におけるローレンツプロットでは、所定距離Δｘとして５６０μｍが用いられた。

また、上述の各実施形態では、肌の凹凸形状に関する空間波形情報が用いられたが、測定部により肌の凹凸の高さの検出位置が特定できない場合、その高さの検出時に基づく肌の凹凸形状に関する時間変化を示す時間波形情報が用いられてもよい。この場合の時間波形情報は、時間軸上に並ぶ凹凸形状に関わる振幅値（高さ、パワーなど）のデータ列となる。

［肌状態評価装置］
上述の各実施形態及び変形例に係る肌状態評価方法は、図１４及び図１５例示される肌状態評価装置（以降、単に評価装置と略称する場合もある）により実行可能である。
図１４は、肌状態評価装置１０のハードウェア構成例を概念的に示す図である。
評価装置１０は、いわゆるコンピュータ（情報処理装置）であり、例えば、バスで相互に接続される、ＣＰＵ（Central Processing Unit）１１、メモリ１２、入出力インタフェース（Ｉ／Ｆ）１３、通信ユニット１４等を有する。評価装置１０を形成する各ハードウェア要素の数はそれぞれ制限されず、これらハードウェア要素は情報処理回路と総称することもできる。また、評価装置１０は、図１４に図示されないハードウェア要素を含んでもよく、そのハードウェア構成は制限されない。

ＣＰＵ１１は、一般的なＣＰＵ以外に、特定用途向け集積回路（ＡＳＩＣ）、ＤＳＰ（Digital Signal Processor）、ＧＰＵ（Graphics Processing Unit）等で構成してもよい。
メモリ１２は、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、補助記憶装置（ハードディスク等）である。
入出力Ｉ／Ｆ１３は、出力装置１５、入力装置１６等のユーザインタフェース装置と接続可能である。出力装置１５は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）ディスプレイのような、ＣＰＵ１１等により処理された描画データに対応する画面を表示する装置、印刷装置などの少なくとも一つである。入力装置１６は、キーボード、マウス等のようなユーザ操作の入力を受け付ける装置である。出力装置１５及び入力装置１６は一体化され、タッチパネルとして実現されてもよい。
通信ユニット１４は、他のコンピュータとの通信網を介した通信や、他の機器との信号のやりとり等を行う。通信ユニット１４には、可搬型記録媒体等も接続され得る。また、通信ユニット１４には、肌の触感に関わる物理量を検出可能な上述の検出部が接続されてもよい。

図１５は、肌状態評価装置１０の処理構成例を概念的に示す図である。
評価装置１０は、取得部２１、生成部２２、及び出力処理部２３を有する。これら処理モジュールは、ソフトウェア要素であり、例えば、メモリ１２に格納される肌状態評価プログラム（コンピュータプログラム）１８がＣＰＵ１１により実行されることにより実現される。このコンピュータプログラム１８は、例えば、ＣＤ（Compact Disc）、メモリカード等のような可搬型記録媒体やネットワーク上の他のコンピュータから入出力Ｉ／Ｆ１３又は通信ユニット１４を介してインストールされ、メモリ１２に格納されてもよい。
それら各ソフトウェア要素の動作により、評価装置１０は、上述の各実施形態及び変形例に係る肌状態評価方法を実行することができる。即ち、評価装置１０（ＣＰＵ１１）は、図１、図２、図７、図９、又は図１１に示される各工程を順次実行することができる。各工程の処理内容はそれぞれ上述の通りであるが、以下に、各工程の処理内容が適宜補足される。

取得部２１は、対象肌の凹凸形状に関する空間波形情報を取得する。即ち、取得部２１は、上述の工程（Ｓ１１）、工程（Ｓ２１）、工程（Ｓ３１）、工程（Ｓ４１）、又は工程（Ｓ５１）を実行する。
取得部２１は、肌の凹凸形状に関する空間波形情報を外部のコンピュータ、機器（検出部など）、可搬型記録媒体などから入出力Ｉ／Ｆ１３又は通信ユニット１４を介して取得することができる。
取得部２１は、測定装置により測定された肌の凹凸波形データを外部から取得し、この凹凸波形データに対して周波数解析などの演算を適用し、変換波形データを当該空間波形情報として取得してもよい。更に、取得部２１は、そのような変換波形データを外部から取得することもできる。

生成部２２は、取得部２１で取得された空間波形情報を用いて、対象空間波形の空間変動特性又は対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を生成する。即ち、生成部２２は、上述の工程（Ｓ１３）を実行する。
具体的には、生成部２２は、図２に示される工程（Ｓ２２）から工程（Ｓ２４）を実行して、工程（Ｓ２３）で得られる変換座標点群の分布状態を示す一種以上の指標値を対象肌の対象空間波形の空間変動特性を示す特性情報として算出することができる。また、生成部２２は、工程（Ｓ２２）の実行により得られる直交座標系の座標点群又は工程（Ｓ２２）及び工程（Ｓ２３）の実行により得られる変換座標点群がプロットされたグラフを当該対象空間波形の空間変動特性を示す特性情報として生成することもできる。また、生成部２２は、図７に示される工程（Ｓ３２）及び工程（Ｓ３３）を実行することにより、元ベクトルとシフトベクトルとの演算結果を当該対象空間波形の空間変動特性を示す特性情報として算出することもできる。
更に、生成部２２は、図９に示される工程（Ｓ４２）から工程（Ｓ４５）を実行することにより、対象肌に関して算出された一種以上の指標値と比較対象肌に関する一種以上の指標値との類似度を、上記相互関係特性を示す特性情報として算出することもできる。加えて、生成部２２は、図１１に示される工程（Ｓ５２）を実行することにより、対象肌に関する空間波形情報を成分とするベクトルと比較対象肌の空間波形情報を成分とするベクトルとの類似度を上記相互関係特性を示す特性情報として算出することもできる。

生成部２２により利用される、比較対象肌の凹凸形状に関する空間波形情報から得られる一種以上の指標値又は比較対象肌の凹凸形状に関する空間波形情報は、メモリ１２に予め格納されていてもよいし、通信ユニット１４を介して外部から通信により取得されてもよい。

出力処理部２３は、生成部２２により生成された特性情報又はその特性情報に基づく対象肌の状態の評価情報を出力する。即ち、出力処理部２３は、上述の工程（Ｓ１５）、工程（Ｓ２５）、工程（Ｓ３４）、工程（Ｓ４６）、又は工程（Ｓ５３）を実行する。
この評価情報は、例えば、対象肌（被験者の肌）の評価結果を示す情報である。例えば、当該評価情報は、被験者の肌の凹凸形状に関わる肌状態が良いか悪いか、肌の凹凸構造の空間均一性の良し悪し、被験者の肌状態が比較対象肌に近いか遠いかなどを示す。

出力処理部２３は、特性情報又は評価情報を出力装置１５に出力することができる。また、出力処理部２３は、特性情報又は評価情報を通信ユニット１４を介して外部のコンピュータに出力してもよいし、可搬型記録媒体に書き込んでもよい。出力処理部２３による特性情報及び評価情報の出力形態は制限されない。

以下に実施例を挙げ、上述の内容を更に詳細に説明する。以下の実施例の記載は、上述の各実施形態及び変形例に何ら限定を加えるものではない。

第一実施形態で述べた、肌の凹凸形状に関する空間波形情報からローレンツプロットを用いて算出された指標値により人の肌状態の違いが定量的に評価可能であることが、次のように検証された。
図１６は、皮膚外用剤適用前後の肌に関してローレンツプロットで得られるグラフを示す図である。図１６において、左側の各グラフは、皮膚外用剤適用前の被験者の肌に対応し（ＧＩ−１、ＧＩＩ−１、ＧＩＩＩ−１、ＧＩＶ−１）、右側の各グラフは、同一被験者に対して所定の皮膚外用剤を適用した後の肌に対応する（ＧＩ−２、ＧＩＩ−２、ＧＩＩＩ−２、ＧＩＶ−２）。また、それぞれ異なる空間周波数帯を選択したウェーブレット変換適用後の変換波形データ（凹凸パワー波形データ）から変換された直交座標系の座標点群がプロットされたグラフが縦に並んでいる。このローレンツプロットでは、所定距離Δｘとして５６０μｍが用いられた。
図１６によれば、皮膚外用剤適用前後において座標点群の分布形状が異なっている。特に、１６８Ｈｚから１２３８Ｈｚの周波数帯に対応するグラフでは、皮膚外用剤適用後に分布範囲が明らかに小さくなっている。これは、肌の凹凸形状の高さ及び周期の乱れが低減していること、即ち、皮膚外用剤により肌の凹凸形状がなめらかになりいわゆるキメが整ったことを示しているといえる。言い換えれば、皮膚外用剤により肌の凹凸構造の空間均一性が向上したことが示されている。
一方、皮膚外用剤の適用により被験者の肌のキメ状態は改善されたことが確認されている。
即ち、図１６によれば、肌の凹凸形状に関する空間波形情報からローレンツプロットを用いて得られるグラフにより人の肌状態の違いを可視化できることが実証されている。

図１７は、図１６に示される１６８Ｈｚから１２３８Ｈｚの周波数帯に対応する直交座標系の座標点群に基づいて算出された５種の指標値を示す図である。具体的には、図１７には、皮膚外用剤適用前後の肌の凹凸形状に関する空間波形情報から第一実施形態の手法で得られた５種の指標値、Ｘ'軸及びＹ'軸の平均距離、Ｘ'軸及びＹ'軸の標準偏差、並びに楕円面積（分布面積）が示されている。各グラフでは、皮膚外用剤適用前の被験者の肌に関する値（左側）と皮膚外用剤適用後の同一被験者の肌に関する値（右側）とが示されている。
図１７によれば、いずれの指標値においても、皮膚外用剤適用後の肌に関する指標値が皮膚外用剤適用前のそれよりも小さくなっている。これは、図１６の場合と同様に、皮膚外用剤により肌の凹凸構造の空間均一性が向上したことを示している。
このように、被験者の肌の凹凸形状の空間波形情報から得られた複数種の指標値により肌状態の違いを定量化可能であることが実証されている。
更に言えば、図１６及び図１７によれば、本実施形態が皮膚外用剤効果による肌状態の違いも可視化可能であることが確認される。

図１８は、皮膚外用剤適用前後の同一被験者の肌に関する凹凸パワー波形データから第二実施形態の手法で得られたベクトル間の内積演算結果を示す図である。図１８では、皮膚外用剤適用前の被験者の肌に関する凹凸パワー波形データＧ１と、同一被験者に対して所定の皮膚外用剤を適用した後の肌に関する凹凸パワー波形データＧ２とが示されている。そして、これら変換波形データを成分とする元ベクトルとその変換波形データを所定距離Δｘ（５６０μｍ）シフトしたシフトベクトルとの内積演算により、内積値及びｃｏｓθが算出されている。
図１８によれば、内積値及びｃｏｓθのいずれにおいても、皮膚外用剤適用前後で差が生じている。ｃｏｓθの値については、皮膚外用剤適用後（Ｇ２）のほうが皮膚外用剤適用前（Ｇ１）よりも大きくなっている。これは、皮膚外用剤適用後の肌の凹凸パワー波形データのほうが、皮膚外用剤適用前の肌のそれよりも空間変動が小さいことを示しており、肌の凹凸構造の空間均一性が皮膚外用剤により向上したことを示しているといえる。
このように、肌の凹凸形状に関する空間波形情報から第二実施形態の手法で得られた演算結果においてもそれぞれ皮膚外用剤に伴う肌状態の違いを表していることが実証される。

図１９は、被験者の肌と比較対象肌との類似性と第三実施形態の手法で導出された指標値間の類似度との関係を示す図である。図１９には、被験者の肌に関してローレンツプロットされた２つのグラフＧ１０及びＧ１１と、比較対象肌に関してローレンツプロットされたグラフＧ１２とが示されている。グラフＧ１１は、皮膚外用剤を適用した後の被験者の肌を示している。
上述したように、皮膚外用剤の適用により被験者の肌のキメ状態は改善されたことが確認されている。また、グラフＧ１２で示される比較対象肌は、キメが細かく整った肌である。つまり、被験者の実際の肌状態（凹凸形状に関わる状態）は、グラフＧ１０で示される状態よりもグラフＧ１１で示される状態のほうが、比較対象肌に近い状態であった。
一方、グラフＧ１０、Ｇ１１、及びＧ１２で示される各肌の凹凸形状に関して、第一実施形態及び第三実施形態の手法により、Ｘ'軸及びＹ'軸の平均距離、Ｘ'軸及びＹ'軸の標準偏差、並びに楕円面積の５種の指標値がそれぞれ算出された。そして、各肌の指標値間の類似度として各ベクトルのなす角度ｃｏｓθが算出された。
結果、グラフＧ１１とグラフＧ１２との指標値間の類似度がグラフＧ１０とグラフＧ１２との間の類似度よりも高い値を示した。これは、実際の肌のキメ状態の類似性と一致しているため、第三実施形態の手法で算出された複数種の特徴量間の類似度により、凹凸形状に関わる実際の肌状態の近さが定量化できることを実証している。

図２０は、被験者の肌と比較対象肌との類似性と第四実施形態の手法で導出されたベクトル間の類似度との関係を示す図である。図２０には、被験者の肌の凹凸波形データからウェーブレット変換により得られた変換波形データを示す二つのグラフＧ２０及びＧ２１と、比較対象肌の凹凸波形データからウェーブレット変換により得られた変換波形データを示すグラフＧ２２とが示されている。グラフＧ２１は、所定の皮膚外用剤を適用した後の被験者の肌を示している。
各肌の実際の状態（凹凸形状に関わる状態）は、図１９で述べたものと同様である。
一方で、グラフＧ２０、Ｇ２１、及びＧ２２で示される各肌に関する変換波形データが第四実施形態の手法によりベクトル化され、それらベクトル間の内積演算により得られるｃｏｓθがベクトル間の類似度として算出された。
結果、グラフＧ２１とグラフＧ２２とのベクトル間の類似度がグラフＧ２０とグラフＧ２２との間の類似度よりも高い値を示した。これは、実際の肌のキメ状態の類似性と一致しているため、第四実施形態の手法で算出されたベクトル間の類似度により、凹凸形状に関わる実際の肌状態の近さが定量化できることを実証している。

上述の内容の一部又は全部は、次のようにも特定され得る。但し、上述の内容が以下の記載に制限されるものではない。

＜１＞ 肌の凹凸形状に関する空間波形情報を取得する工程と、
前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を導出する導出工程と、
前記導出された特性情報に基づいて、前記肌の状態を評価する評価工程と、
を含む肌状態評価方法。

＜２＞ 前記空間変動特性又は前記相互関係特性は前記対象空間波形の周期乱れを示す、
＜１＞に記載の肌状態評価方法。
＜３＞ 前記導出工程では、
前記取得された空間波形情報を、前記対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換し、
前記直交座標系の前記座標点群を、前記直交座標系の直交軸であるｘ軸及びｙ軸で示される「ｙ＝ｘ」軸及び「ｙ＝−ｘ」軸に投影した変換座標点群に更に変換し、
前記変換座標点群の分布状態を示す一種以上の指標値を前記対象空間波形の前記空間変動特性を示す前記特性情報として導出する、
＜１＞又は＜２＞記載の肌状態評価方法。
＜４＞ 前記導出工程では、前記変換座標点群における、「ｙ＝ｘ」軸における原点からの距離の平均値若しくは標準偏差、「ｙ＝−ｘ」軸における原点からの距離の平均値若しくは標準偏差、又は、該平均値若しくは該標準偏差を用いて算出される比率若しくは分布面積の一つ以上を前記一種以上の指標値として導出する、
＜３＞記載の肌状態評価方法。
＜５＞ 前記導出工程では、前記変換座標点群の分布状態を示す複数種の指標値を導出し、
前記評価工程では、
各指標種を座標軸とする座標系に、前記肌に関する前記複数種の指標値に対応する点がプロットされたグラフを提示する、
＜３＞又は＜４＞記載の肌状態評価方法。
＜６＞ 前記導出工程では、
前記取得された空間波形情報を、前記対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換し、
前記直交座標系の座標点群又は該座標点群に対する座標変換により得られる変換座標点群がプロットされたグラフを前記対象空間波形の前記空間変動特性を示す前記特性情報として導出し、
前記評価工程では、前記グラフを提示する、
＜１＞から＜５＞のいずれか一項に記載の肌状態評価方法。
＜７＞ 前記導出工程では、
前記取得された空間波形情報を成分とする元ベクトルと該空間波形情報を空間軸上で所定距離シフトしたシフト波形情報を成分とするシフトベクトルとの演算により、前記対象空間波形の前記空間変動特性を示す前記特性情報を導出する、
＜１＞から＜６＞のいずれか一項に記載の肌状態評価方法。
＜８＞ 前記導出工程では、前記肌に関して導出された前記一種以上の指標値と前記比較対象肌に関する一種以上の指標値との類似度を、前記相互関係特性を示す前記特性情報として導出する、
＜３＞から＜５＞のいずれか一項に記載の肌状態評価方法。
＜９＞ 前記導出工程では、前記肌の前記対象空間波形の前記空間変動特性を示す前記特性情報として導出された前記グラフと、前記比較対象肌に関して前記肌と同様に前記空間変動特性を示す前記特性情報として導出されたグラフとが対比可能な状態で示される比較情報を、前記相互関係特性を示す前記特性情報として導出する、
＜６＞記載の肌状態評価方法。
＜１０＞ 前記導出工程では、前記取得された空間波形情報を成分とするベクトルと前記比較対象肌の空間波形情報を成分とするベクトルとの類似度を前記相互関係特性を示す前記特性情報として導出する、
＜１＞又は＜２＞記載の肌状態評価方法。
＜１１＞ 肌の凹凸形状に関する空間波形情報を取得する取得手段と、
前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を生成する生成手段と、
前記生成された特性情報又は該特性情報に基づく前記肌の状態の評価情報を出力する出力処理手段と、
を備える肌状態評価装置。
＜１２＞ 前記空間変動特性又は前記相互関係特性は前記対象空間波形の周期乱れを示す、
＜１１＞に記載の肌状態評価装置。
＜１３＞ 前記生成手段は、
前記取得された空間波形情報を、前記対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換し、
前記直交座標系の前記座標点群を、前記直交座標系の直交軸であるｘ軸及びｙ軸で示される「ｙ＝ｘ」軸及び「ｙ＝−ｘ」軸に投影した変換座標点群に更に変換し、
前記変換座標点群の分布状態を示す一種以上の指標値を前記対象空間波形の前記空間変動特性を示す前記特性情報として算出する、
＜１１＞又は＜１２＞記載の肌状態評価装置。
＜１４＞ 前記生成手段は、前記変換座標点群における、「ｙ＝ｘ」軸における原点からの距離の平均値若しくは標準偏差、「ｙ＝−ｘ」軸における原点からの距離の平均値若しくは標準偏差、又は、該平均値若しくは該標準偏差を用いて算出される比率若しくは分布面積の一つ以上を前記一種以上の指標値として算出する、
＜１３＞記載の肌状態評価装置。
＜１５＞ 前記生成手段は、前記変換座標点群の分布状態を示す複数種の指標値を算出し、
前記出力処理手段は、各指標種を座標軸とする座標系に、前記肌に関する前記複数種の指標値に対応する点がプロットされたグラフを出力する、
＜１３＞又は＜１４＞記載の肌状態評価装置。
＜１６＞ 前記生成手段は、
前記取得された空間波形情報を、前記対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換し、
前記直交座標系の座標点群又は該座標点群に対する座標変換により得られる変換座標点群がプロットされたグラフを前記対象空間波形の前記空間変動特性を示す前記特性情報として生成し、
前記出力処理手段は、前記グラフを出力する、
＜１１＞から＜１５＞のいずれか一項に記載の肌状態評価装置。
＜１７＞ 前記生成手段は、
前記取得された空間波形情報を成分とする元ベクトルと該空間波形情報を空間軸上で所定距離シフトしたシフト波形情報を成分とするシフトベクトルとの演算により、前記対象空間波形の前記空間変動特性を示す前記特性情報を算出する、
＜１１＞から＜１６＞のいずれか一項に記載の肌状態評価装置。
＜１８＞ 前記生成手段は、前記肌に関して導出された前記一種以上の指標値と前記比較対象肌に関する一種以上の指標値との類似度を、前記相互関係特性を示す前記特性情報として算出する、
＜１３＞から＜１５＞のいずれか一項に記載の肌状態評価装置。
＜１９＞ 前記出力処理手段は、前記生成手段により生成された前記グラフと、前記比較対象肌に関して前記肌と同様に前記空間変動特性を示す前記特性情報として生成されたグラフとを対比可能な状態で出力する、
＜１６＞記載の肌状態評価装置。
＜２０＞ 前記生成手段は、前記取得された空間波形情報を成分とするベクトルと前記比較対象肌の空間波形情報を成分とするベクトルとの類似度を前記相互関係特性を示す前記特性情報として算出する、
＜１１＞又は＜１２＞記載の肌状態評価装置。

１０ 肌状態評価装置（評価装置）
１１ ＣＰＵ
１２ メモリ
１３ 入出力Ｉ／Ｆ
１４ 通信ユニット
１５ 出力装置
１６ 入力装置
１８ 肌状態評価プログラム
２１ 取得部
２２ 生成部
２３ 出力処理部

Claims (10)

1. 肌の凹凸形状に関する空間波形情報を取得する工程と、
   前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を導出する導出工程と、
   前記導出された特性情報に基づいて、前記肌の状態を評価する評価工程と、
   を含み、
   前記導出工程では、
   前記取得された空間波形情報を、前記対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換し、
   前記直交座標系の前記座標点群を、前記直交座標系の直交軸であるｘ軸及びｙ軸で示される「ｙ＝ｘ」軸及び「ｙ＝−ｘ」軸に投影した変換座標点群に更に変換し、
   前記変換座標点群の分布状態を示す一種以上の指標値を前記対象空間波形の前記空間変動特性を示す前記特性情報として導出する、
   肌状態評価方法。
2. 前記導出工程では、前記変換座標点群における、「ｙ＝ｘ」軸における原点からの距離の平均値若しくは標準偏差、「ｙ＝−ｘ」軸における原点からの距離の平均値若しくは標準偏差、又は、該平均値若しくは該標準偏差を用いて算出される比率若しくは分布面積の一つ以上を前記一種以上の指標値として導出する、
   請求項１に記載の肌状態評価方法。
3. 前記導出工程では、前記変換座標点群の分布状態を示す複数種の指標値を導出し、
   前記評価工程では、
   各指標種を座標軸とする座標系に、前記肌に関する前記複数種の指標値に対応する点がプロットされたグラフを提示する、
   請求項１又は２に記載の肌状態評価方法。
4. 前記導出工程では、前記肌に関して導出された前記一種以上の指標値と前記比較対象肌に関する一種以上の指標値との類似度を、前記相互関係特性を示す前記特性情報として導出する、
   請求項１から３のいずれか一項に記載の肌状態評価方法。
5. 肌の凹凸形状に関する空間波形情報を取得する工程と、
   前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を導出する導出工程と、
   前記導出された特性情報に基づいて、前記肌の状態を評価する評価工程と、
   を含み、
   前記導出工程では、
   前記取得された空間波形情報を、前記対象空間波形における空間軸上で所定距離離れた振幅値ペアで形成される直交座標系の座標点群に変換し、
   前記直交座標系の座標点群又は該座標点群に対する座標変換により得られる変換座標点群がプロットされたグラフを前記対象空間波形の前記空間変動特性を示す前記特性情報として導出し、
   前記評価工程では、前記グラフを提示する、
   肌状態評価方法。
6. 前記導出工程では、前記肌の前記対象空間波形の前記空間変動特性を示す前記特性情報として導出された前記グラフと、前記比較対象肌に関して前記肌と同様に前記空間変動特性を示す前記特性情報として導出されたグラフとが対比可能な状態で示される比較情報を、前記相互関係特性を示す前記特性情報として導出する、
   請求項５に記載の肌状態評価方法。
7. 肌の凹凸形状に関する空間波形情報を取得する工程と、
   前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を導出する導出工程と、
   前記導出された特性情報に基づいて、前記肌の状態を評価する評価工程と、
   を含み、
   前記導出工程では、
   前記取得された空間波形情報を成分とする元ベクトルと該空間波形情報を空間軸上で所定距離シフトしたシフト波形情報を成分とするシフトベクトルとの演算により、前記対象空間波形の前記空間変動特性を示す前記特性情報を導出する、
   肌状態評価方法。
8. 肌の凹凸形状に関する空間波形情報を取得する工程と、
   前記取得された空間波形情報を用いて、対象空間波形の空間変動特性又は該対象空間波形と比較対象肌の空間波形との間の相互関係特性を示す特性情報を導出する導出工程と、
   前記導出された特性情報に基づいて、前記肌の状態を評価する評価工程と、
   を含み、
   前記導出工程では、前記取得された空間波形情報を成分とするベクトルと前記比較対象肌の空間波形情報を成分とするベクトルとの類似度を前記相互関係特性を示す前記特性情報として導出する、
   肌状態評価方法。
9. プロセッサ及びメモリを少なくとも備える肌状態評価装置であって、請求項１から８のいずれか一項の肌状態評価方法を実行する肌状態評価装置。
10. プロセッサ及びメモリを少なくとも備えるコンピュータに請求項１から８のいずれか一項の肌状態評価方法を実行させるコンピュータプログラム。

Images