JP6942012B2 - Skin evaluation method - Google Patents
Skin evaluation method Download PDFInfo
- Publication number
- JP6942012B2 JP6942012B2 JP2017175149A JP2017175149A JP6942012B2 JP 6942012 B2 JP6942012 B2 JP 6942012B2 JP 2017175149 A JP2017175149 A JP 2017175149A JP 2017175149 A JP2017175149 A JP 2017175149A JP 6942012 B2 JP6942012 B2 JP 6942012B2
- Authority
- JP
- Japan
- Prior art keywords
- skin
- spatial frequency
- frequency band
- evaluation
- information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Description
本発明は、肌(皮膚)の表面形状に基づく肌評価技術に関する。 The present invention relates to a skin evaluation technique based on the surface shape of skin (skin).
肌表面の画像又は三次元形状データを用いて皮膚表面を定量評価する各種手法が提案されている。
特許文献1には、皮膚表面又はそのレプリカに光を照射して得られる皮膚表面画像の輝度を、多段階の階層(階調)を有するデジタル信号値に変換し、複数の測定点において、その輝度デジタル信号値を読み取り、各測定点間における輝度デジタル信号値差を統計処理し、輝度差の平均値、標準偏差、メジアン及びモードを皮膚特性値として算出する皮膚表面形状測定方法が開示されている。
特許文献2には、皮膚表面におけるしわやキメの視覚的粗さ及び皮膚化粧料が持つしわ・キメの隠ぺい効果を数値化可能とする皮膚表面解析方法が開示されている。本解析方法では、サンプル皮膚表面の微細明暗分布が強調された二次元サンプル画像を撮像し、この二次元サンプル画像を加工処理した後、一次元明暗プロファイルを抽出し、この一次元明暗プロファイルに基づいて明暗ピーク間距離の平均値、標準偏差及び変動係数を算出する。
特許文献3には、第三者の目に映る、肌の見た目の美しさを、客観的な数値に置換する技術が開示されている。本技術では、肌のレプリカ標本の凹凸状態を計測した3次元形状データから得られる表面粗さパラメータであるSm(凹凸の平均間隔)及び/又はPmr(負荷長さ率)及び/又Ra(算術平均粗さ)と、肌の写真を主観的判断して得た肌の美しさのスコアとの相関関係から得られる回帰式(予測式)に、鑑別すべき肌のレプリカ標本のSm及び/又はPmr及び/又Raを適用して、肌の美しさのスコアを得る。
Various methods for quantitatively evaluating the skin surface using images of the skin surface or three-dimensional shape data have been proposed.
In Patent Document 1, the brightness of a skin surface image obtained by irradiating a skin surface or a replica thereof with light is converted into a digital signal value having a multi-step hierarchy (gradation), and the brightness is converted into a digital signal value having a multi-step hierarchy (gradation). A skin surface shape measuring method that reads a luminance digital signal value, statistically processes the luminance digital signal value difference between each measurement point, and calculates the average value, standard deviation, median, and mode of the luminance difference as skin characteristic values is disclosed. There is.
Patent Document 3 discloses a technique for replacing the beauty of the appearance of the skin seen by a third party with an objective numerical value. In this technique, Sm (average interval of unevenness) and / or Pmr (load length ratio) and / or Ra (arithmetic), which are surface roughness parameters obtained from three-dimensional shape data obtained by measuring the unevenness state of a skin replica sample. Sm and / or Sm and / or of the skin replica sample to be differentiated by the regression equation (prediction equation) obtained from the correlation between the average roughness) and the skin beauty score obtained by subjectively judging the skin photograph. Apply Pmr and / or Ra to get a skin beauty score.
しかしながら、上述の各手法には、肌の表面形状の解析について改善の余地がある。
肌の表面形状は、寸法の異なる複数の凹凸要素からなる多重構造を持つことが知られている。例えば、肌表面には、シワ、毛穴、皮丘と皮溝とから形成されるキメ、皮丘上に形成された凹凸形状(以降、二次レリーフと表記する)などが存在し、それら凹凸要素はそれぞれ寸法が異なる。上述の各手法では、このような肌表面の多重構造を考慮した処理が行われていないため、算出された肌表面に関する指標値には改善の余地が残る。
本発明は、肌表面の形状に基づいて肌を評価する技術に関する。
However, each of the above methods has room for improvement in the analysis of the surface shape of the skin.
It is known that the surface shape of the skin has a multiple structure composed of a plurality of uneven elements having different dimensions. For example, on the surface of the skin, there are wrinkles, pores, texture formed from skin hills and skin grooves, uneven shape formed on the skin hill (hereinafter referred to as secondary relief), and the like. Have different dimensions. In each of the above-mentioned methods, since the treatment considering the multiple structure of the skin surface is not performed, there is room for improvement in the calculated index value related to the skin surface.
The present invention relates to a technique for evaluating skin based on the shape of the skin surface.
本発明の態様では、上述した課題を解決するために、以下の構成を採用する。
本発明の態様は、肌評価方法に関する。肌評価方法は、肌又は肌のレプリカから計測された肌表面の凹凸情報を取得する工程と、前記取得された凹凸情報から所定の空間周波数帯域の凹凸情報を抽出する分解工程と、前記抽出された所定の空間周波数帯域の凹凸情報から凹凸形状に関する統計値を算出する算出工程と、前記算出された統計値を用いて前記肌を評価する評価工程と、を含む。
なお、本発明の別態様は、例えば、上述の肌評価方法を実行する肌評価装置、コンピュータ(情報処理装置)に関するものであり、上述の肌評価方法をコンピュータに実行させるプログラムに関するものであり、このようなプログラムを記録したコンピュータが読み取り可能な記憶媒体に関するものである。この記録媒体は、非一時的な有形の媒体を含む。
In the aspect of the present invention, the following configuration is adopted in order to solve the above-mentioned problems.
Aspects of the present invention relate to skin evaluation methods. The skin evaluation method includes a step of acquiring unevenness information on the skin surface measured from the skin or a replica of the skin, a decomposition step of extracting unevenness information in a predetermined spatial frequency band from the acquired unevenness information, and the extraction. It includes a calculation step of calculating a statistical value regarding the uneven shape from the unevenness information of a predetermined spatial frequency band, and an evaluation step of evaluating the skin using the calculated statistical value.
Another aspect of the present invention relates to, for example, a skin evaluation device and a computer (information processing device) that execute the above-mentioned skin evaluation method, and relates to a program that causes a computer to execute the above-mentioned skin evaluation method. It relates to a computer-readable storage medium on which such a program is recorded. This recording medium includes a non-temporary tangible medium.
上記態様によれば、肌表面の形状に基づいて高精度に肌を評価する技術を提供することができる。 According to the above aspect, it is possible to provide a technique for evaluating the skin with high accuracy based on the shape of the skin surface.
以下、本発明の好ましい実施形態の例(以降、単に本実施形態と略称する)について説明する。以下に挙げる各実施形態はそれぞれ例示であり、本発明は以下の各実施形態の構成に限定されない。 Hereinafter, an example of a preferred embodiment of the present invention (hereinafter, simply abbreviated as the present embodiment) will be described. Each of the following embodiments is an example, and the present invention is not limited to the configuration of each of the following embodiments.
[第一実施形態]
第一実施形態に係る肌評価方法について図1を用いて説明する。
図1は、第一実施形態に係る肌評価方法を示すフローチャートである。
第一実施形態に係る肌評価方法は、図1に示されるように、工程(S11)、工程(S13)、工程(S15)、及び工程(S17)を含む。
[First Embodiment]
The skin evaluation method according to the first embodiment will be described with reference to FIG.
FIG. 1 is a flowchart showing a skin evaluation method according to the first embodiment.
As shown in FIG. 1, the skin evaluation method according to the first embodiment includes a step (S11), a step (S13), a step (S15), and a step (S17).
工程(S11)は、対象肌又は対象肌のレプリカから計測される等の肌表面の凹凸情報を取得する工程である。
工程(S11)で取得される凹凸情報は、肌表面の高さ情報を少なくとも含んでいればよく、高さ情報を含む三次元情報、肌の断面形状を示す波形情報など、その具体的なデータ形式等は制限されない。当該凹凸情報は、肌表面の位置情報と高さ情報との二次元情報である凹凸波形データであってもよいし、高さを輝度で示す画像データであってもよい。例えば、共焦点レーザ顕微鏡等を用いて対象肌又は対象肌のレプリカを計測することで、肌表面の三次元形状を示すデータを得ることもでき、そのデータ中において所望のラインを指定することで、そのライン上の凹凸形状プロファイル(凹凸波形データ)を得ることもできる。凹凸情報を取得することが可能な機器、すなわち工程(S11)で用いることができる機器は、共焦点レーザ顕微鏡に限らず、光学系の三次元測定器、接触式又は非接触式の表面粗さ測定器等、様々な測定器が存在するが、空間周波数分解能が1/20(cycle/μm)以上、高さ方向の分解能(検出能)が0.05μm以下の機器を用いることが好ましい。
The step (S11) is a step of acquiring unevenness information on the skin surface such as measured from the target skin or a replica of the target skin.
The unevenness information acquired in the step (S11) may include at least the height information of the skin surface, and specific data such as three-dimensional information including the height information and waveform information indicating the cross-sectional shape of the skin. The format is not restricted. The unevenness information may be unevenness waveform data which is two-dimensional information of the position information and height information of the skin surface, or may be image data indicating the height by brightness. For example, by measuring the target skin or a replica of the target skin using a confocal laser microscope or the like, it is possible to obtain data showing the three-dimensional shape of the skin surface, and by specifying a desired line in the data. , It is also possible to obtain a concave-convex shape profile (concave-convex waveform data) on the line. The device capable of acquiring unevenness information, that is, the device that can be used in the step (S11) is not limited to the confocal laser microscope, but is a three-dimensional measuring device of an optical system, a contact type or a non-contact type surface roughness. There are various measuring instruments such as measuring instruments, but it is preferable to use an instrument having a spatial frequency resolution of 1/20 (cycle / μm) or more and a resolution (detectability) in the height direction of 0.05 μm or less.
図2は、肌表面の凹凸情報の例を示す図である。
図2の例では、被験者の頬から肌表面形状が転写されたレプリカを採取し、そのレプリカからレーザ顕微鏡により計測された画像(三次元情報)が上部に示されている。更に、その画像上で所望のラインが指定され、そのライン上における位置情報及び高さ情報からなる凹凸波形データが下部に示されている。
図2の例では、肌表面の凹凸情報が肌のレプリカから計測されたが、肌から直接計測されてもよい。
FIG. 2 is a diagram showing an example of unevenness information on the skin surface.
In the example of FIG. 2, a replica in which the skin surface shape is transferred from the cheek of the subject is collected, and an image (three-dimensional information) measured by a laser microscope from the replica is shown at the upper part. Further, a desired line is specified on the image, and uneven waveform data including position information and height information on the line is shown at the bottom.
In the example of FIG. 2, the unevenness information on the skin surface was measured from the skin replica, but it may be measured directly from the skin.
肌のレプリカは、評価対象となる部位の肌表面から採取される。当該レプリカは、肌表面の凹凸形状が転写されていればよく、その生成手法は限定されない。例えば、シリコン印象材などの剤を肌表面に押し当てることで、レプリカを生成することができる。
また、図2の例では、肌表面の凹凸情報として、高さ情報を含む三次元情報(画像)から取得された凹凸波形データが用いられるが、その三次元情報(画像)がそのまま肌表面の凹凸情報として利用されてもよい。
Skin replicas are collected from the skin surface of the site to be evaluated. The replica only needs to have the uneven shape of the skin surface transferred, and the generation method thereof is not limited. For example, a replica can be produced by pressing an agent such as a silicone impression material against the skin surface.
Further, in the example of FIG. 2, as the unevenness information of the skin surface, the unevenness waveform data acquired from the three-dimensional information (image) including the height information is used, but the three-dimensional information (image) is used as it is on the skin surface. It may be used as unevenness information.
工程(S13)は、工程(S11)で取得された凹凸情報から所定の空間周波数帯域の凹凸情報を抽出する分解工程である。
ここで「所定の空間周波数帯域」とは、工程(S13)の実行前に既に決められている空間周波数帯域を意味する。
「空間周波数帯域」とは、肌表面の凹凸形状に適した単位長(例えば、1μm)当たりの高さ変動の回数の範囲を意味する。
所定の空間周波数帯域として、肌の評価指標と相関が認められる空間周波数帯域が設定されることが望ましい。また、肌の表面形状の多重構造を形成する或る凹凸要素の寸法(スケール)に対応する空間周波数帯域が設定されてもよい。肌は部位によって凹凸形状のスケールが異なるため、測定したい部位の形状特徴に合わせて、適宜、空間周波数帯域を選定する。例えば、顔の肌表面の凹凸形状におけるキメや皮丘といった凹凸要素に対応する空間周波数帯域として、1/800(cycle/μm)から1/60(cycle/μm)の少なくとも一部の帯域が設定されてもよいし、二次レリーフの凹凸要素に対応する空間周波数帯域として、1/100(cycle/μm)から1/20(cycle/μm)の少なくとも一部の帯域が設定されてもよいし、それら凹凸要素よりも更に大きい寸法の凹凸要素に対応する空間周波数帯域として、1/400(cycle/μm)以下が設定されてもよい。
The step (S13) is a decomposition step of extracting unevenness information in a predetermined spatial frequency band from the unevenness information acquired in the step (S11).
Here, the "predetermined spatial frequency band" means a spatial frequency band that has already been determined before the execution of the step (S13).
The “spatial frequency band” means a range of the number of height fluctuations per unit length (for example, 1 μm) suitable for the uneven shape of the skin surface.
As a predetermined spatial frequency band, it is desirable to set a spatial frequency band in which a correlation with the skin evaluation index is recognized. Further, a spatial frequency band corresponding to the size (scale) of a certain uneven element forming a multiple structure of the surface shape of the skin may be set. Since the scale of the uneven shape of the skin differs depending on the part, the spatial frequency band is appropriately selected according to the shape characteristics of the part to be measured. For example, at least a part of the band from 1/800 (cycle / μm) to 1/60 (cycle / μm) is set as the spatial frequency band corresponding to the uneven elements such as texture and skin hills in the uneven shape of the skin surface of the face. However, at least a part of the band from 1/100 (cycle / μm) to 1/20 (cycle / μm) may be set as the spatial frequency band corresponding to the uneven element of the secondary relief. , 1/400 (cycle / μm) or less may be set as the spatial frequency band corresponding to the concavo-convex element having a size larger than those concavo-convex elements.
工程(S13)において所定の空間周波数帯域の凹凸情報を抽出する手法には、様々な手法を利用することができる。例えば、工程(S11)で取得された凹凸情報をフーリエ変換等を用いて空間周波数領域に変換し、空間周波数領域においてフィルタを適用した後、逆フーリエ変換等を用いて元の空間領域に戻すことで、所定の空間周波数帯域の凹凸情報を抽出することができる。また、ウェーブレット変換や2CRフィルタ、位相補償形ディジタルフィルタ(ガウシアンフィルタ)など他の周波数解析手法を用いて所定の空間周波数帯域の凹凸情報を抽出することもできる。 Various methods can be used as a method for extracting unevenness information in a predetermined spatial frequency band in the step (S13). For example, the unevenness information acquired in the step (S11) is converted into a spatial frequency domain by using a Fourier transform or the like, a filter is applied in the spatial frequency domain, and then the original spatial domain is restored by using an inverse Fourier transform or the like. Therefore, it is possible to extract unevenness information in a predetermined spatial frequency band. Further, it is also possible to extract unevenness information in a predetermined spatial frequency band by using another frequency analysis method such as a wavelet transform, a 2CR filter, or a phase compensation type digital filter (Gaussian filter).
工程(S15)は、工程(S13)で抽出された所定の空間周波数帯域の凹凸情報から凹凸形状に関する統計値を算出する算出工程である。
工程(S15)で算出される統計値は、肌表面の高さ方向の偏差に基づくパラメータであり、平均、標準偏差、分散、歪度、若しくは尖度、又は、それらのいずれか複数の組合せである。
平均として例えば算術平均粗さ(Ra)が算出される。算術平均粗さ(Ra)は、粗さ曲線の基準長さにおける平均線からの偏差の絶対値の平均値で示される。
標準偏差或いは分散として例えば二乗平均平方根粗さ(Rq)が算出される。二乗平均平方根粗さ(Rq)は、粗さ曲線の基準長さにおける平均線からの偏差の二乗を平均して平方根を取った値で示される。
歪度としては例えば二乗平均平方根粗さ(Rq)の三乗によって無次元化した基準長さにおける凹凸波形の三乗平均であるスキューネス(Rsk)が算出される。このスキューネス(Rsk)によれば、肌表面の凹凸の対称性が示される。
尖度としては例えば二乗平均平方根粗さ(Rq)の四乗によって無次元化した基準長さにおける凹凸波形の四乗平均であるクルトシス(Rku)が算出される。このクルトシス(Rku)によれば、肌表面の凹凸の尖り度が示される。
The step (S15) is a calculation step of calculating a statistical value regarding the shape of the unevenness from the unevenness information of the predetermined spatial frequency band extracted in the step (S13).
The statistical value calculated in the step (S15) is a parameter based on the deviation in the height direction of the skin surface, and is a combination of mean, standard deviation, variance, skewness, kurtosis, or any one or more of them. be.
For example, the arithmetic mean roughness (Ra) is calculated as an average. The arithmetic mean roughness (Ra) is indicated by the average value of the absolute values of the deviations from the average line at the reference length of the roughness curve.
For example, the root mean square roughness (Rq) is calculated as the standard deviation or variance. The root mean square roughness (Rq) is indicated by the value obtained by averaging the square of the deviation from the average line at the reference length of the roughness curve and taking the square root.
As the skewness, for example, skewness (Rsk), which is the cube average of the uneven waveform at the reference length dimensionless by the cube of the root mean square roughness (Rq), is calculated. According to this skewness (Rsk), the symmetry of the unevenness of the skin surface is shown.
As the kurtosis, for example, Kurtosis (Rku), which is the fourth average of the uneven waveforms at the reference length dimensionless by the square of the root mean square roughness (Rq), is calculated. According to this Kurtosis (Rku), the sharpness of the unevenness on the skin surface is shown.
工程(S17)は、工程(S15)で算出された統計値を用いて対象肌を評価する評価工程である。
本発明者らは、様々な手法で取得可能な肌表面の凹凸情報の中の、或る特定の空間周波数帯域の凹凸情報に着目して、その凹凸情報における或る特定種の統計値と肌の或る特定の評価指標との間に高い相関があることを見出した。これは、肌の評価指標が肌の表面形状の多重構造を形成する凹凸要素群の中の或る特定の凹凸要素の状態に強く影響を受けているということもできる。
例えば、後述の実施例で記載するように、肌の触感と、空間周波数帯域(1/500(cycle/μm)以上1/80(cycle/μm)未満)の凹凸形状の歪度との相関が高いことが見出されている。また、空間周波数帯域(1/500(cycle/μm)未満)の算術平均粗さ(Ra)、二乗平均平方根粗さ(Rq)と肌の見た目との相関が高いことが見出されている。
The step (S17) is an evaluation step of evaluating the target skin using the statistical values calculated in the step (S15).
The present inventors pay attention to the unevenness information of a specific spatial frequency band in the unevenness information of the skin surface that can be obtained by various methods, and the statistical value of a specific kind and the skin in the unevenness information. We found that there was a high correlation with a certain evaluation index of. It can also be said that the evaluation index of the skin is strongly influenced by the state of a specific uneven element in the uneven element group forming the multiple structure of the surface shape of the skin.
For example, as described in Examples described later, there is a correlation between the tactile sensation of the skin and the skewness of the uneven shape in the spatial frequency band (1/500 (cycle / μm) or more and less than 1/80 (cycle / μm)). It has been found to be expensive. Further, it has been found that there is a high correlation between the arithmetic mean roughness (Ra) and the root mean square roughness (Rq) of the spatial frequency band (less than 1/500 (cycle / μm)) and the appearance of the skin.
これにより、工程(S17)では、上述の所定の空間周波数帯域の凹凸形状の統計値を用いることで、その統計値と高い相関を示す評価指標(例えば、触感、見た目など)で、対象肌を評価することができる。例えば、空間周波数帯域(1/500(cycle/μm)以上1/80(cycle/μm)未満)の凹凸情報から算出された歪度を用いて、対象肌の触感を評価することができる。また、空間周波数帯域(1/500(cycle/μm)未満)の凹凸情報から算出された算術平均粗さ(Ra)、二乗平均平方根粗さ(Rq)を用いて、対象肌の見た目の印象を評価することができる。もちろん、見た目及び触感のように複数の評価指標で肌が評価されてもよい。また、肌の評価指標は、これら触感及び見た目のみに限定されない。 As a result, in the step (S17), by using the above-mentioned statistical value of the uneven shape of the predetermined spatial frequency band, the target skin is set with an evaluation index (for example, tactile sensation, appearance, etc.) showing a high correlation with the statistical value. Can be evaluated. For example, the tactile sensation of the target skin can be evaluated using the skewness calculated from the unevenness information of the spatial frequency band (1/500 (cycle / μm) or more and less than 1/80 (cycle / μm)). In addition, the arithmetic average roughness (Ra) and the root mean square roughness (Rq) calculated from the unevenness information of the spatial frequency band (less than 1/500 (cycle / μm)) are used to give an impression of the appearance of the target skin. Can be evaluated. Of course, the skin may be evaluated by a plurality of evaluation indexes such as appearance and touch. Moreover, the evaluation index of the skin is not limited to these tactile sensations and appearances.
工程(S17)の実行主体は、人であってもよいし、コンピュータであってもよい。
例えば、工程(S11)、工程(S13)及び工程(S15)がコンピュータにより実行され、工程(S15)で算出された統計値がそのコンピュータから表示装置や印刷装置等を介して出力される。この場合、人がその出力された統計値を参照し、その統計値と比較表とを見比べて、対象肌の評価を行うことができる。比較表は、例えば、肌の或る評価指標の評価結果の分類(例えば、見た目が良い、標準及び悪いといった3分類、きれい及びきれいではないといった2分類など)ごとに統計値の範囲が設定された表であってもよい。この比較表は、年齢ごと、性別ごと、模範となり得る有名人ごとにそれぞれ生成することもでき、予め印刷されていてもよいし、コンピュータから出力されてもよい。
他の例として、工程(S17)がコンピュータにより実行される場合には、コンピュータが工程(S15)で算出された指標値から評価スコアを算出し、その評価スコアを出力してもよい。評価スコアは、良し悪しの2値であってもよいし、3値以上の数値であってもよいし、評価値に対応する文字列(「良い」、「悪い」、「標準」、「きれい」など)であってもよい。コンピュータによる評価スコアの算出手法については後述する。
The execution subject of the process (S17) may be a person or a computer.
For example, the steps (S11), steps (S13), and steps (S15) are executed by a computer, and the statistical values calculated in the steps (S15) are output from the computer via a display device, a printing device, or the like. In this case, a person can refer to the output statistical value, compare the statistical value with the comparison table, and evaluate the target skin. In the comparison table, for example, a range of statistical values is set for each classification of evaluation results of a certain evaluation index of skin (for example, 3 classifications such as good appearance, standard and bad appearance, 2 classifications such as clean and not clean). It may be a table. This comparison table can be generated for each age, gender, and potential celebrity, and may be pre-printed or output from a computer.
As another example, when the step (S17) is executed by the computer, the computer may calculate the evaluation score from the index value calculated in the step (S15) and output the evaluation score. The evaluation score may be a binary value of good or bad, a numerical value of 3 values or more, and a character string corresponding to the evaluation value ("good", "bad", "standard", "clean". , Etc.). The calculation method of the evaluation score by the computer will be described later.
このように、第一実施形態では、肌又は肌のレプリカから計測された肌表面の凹凸情報が取得され、その凹凸情報から抽出された所定の空間周波数帯域の凹凸情報の統計値を用いて、肌が評価される。つまり、第一実施形態によれば、肌の表面形状の多重構造を形成する凹凸要素群の中の或る特定の凹凸要素の特徴を限定的に参照することができる。
従って、多重構造を有する肌表面形状に関して、或る特定の凹凸要素に着目することにより、肌を高精度に評価することができる。
As described above, in the first embodiment, the unevenness information of the skin surface measured from the skin or the skin replica is acquired, and the statistical value of the unevenness information of the predetermined spatial frequency band extracted from the unevenness information is used. The skin is evaluated. That is, according to the first embodiment, the characteristics of a specific uneven element in the uneven element group forming the multiple structure of the surface shape of the skin can be referred to in a limited manner.
Therefore, the skin can be evaluated with high accuracy by paying attention to a specific uneven element with respect to the skin surface shape having a multiple structure.
[第二実施形態]
上述の第一実施形態では、特定の空間周波数帯域の凹凸情報が用いられたが、第二実施形態では、当該肌表面の凹凸情報から多重分解された複数の空間周波数帯域の凹凸情報が用いられる。以下、第二実施形態に係る肌評価方法について図3を用いて説明する。
図3は、第二実施形態に係る肌評価方法を示すフローチャートである。
第二実施形態に係る肌評価方法は、図3に示されるように、工程(S31)、工程(S33)、工程(S35)、及び工程(S37)を含む。
工程(S31)は、上述の工程(S11)と同様である。
[Second Embodiment]
In the first embodiment described above, the unevenness information of a specific spatial frequency band is used, but in the second embodiment, the unevenness information of a plurality of spatial frequency bands multiple-decomposed from the unevenness information of the skin surface is used. .. Hereinafter, the skin evaluation method according to the second embodiment will be described with reference to FIG.
FIG. 3 is a flowchart showing a skin evaluation method according to the second embodiment.
As shown in FIG. 3, the skin evaluation method according to the second embodiment includes a step (S31), a step (S33), a step (S35), and a step (S37).
The step (S31) is the same as the above-mentioned step (S11).
工程(S33)は、工程(S31)で取得された凹凸情報を複数の空間周波数帯域へと分解する工程である。即ち、工程(S33)で分解する所定の空間周波数帯域は少なくとも2帯域とされる。
ここで、肌又は肌のレプリカから計測された凹凸情報を複数の空間周波数帯域の凹凸情報に分解することを「多重分解」と表記する場合もある。
当該複数の空間周波数帯域は、工程(S33)の実行前に予め決められている。当該複数の空間周波数帯域には、第一実施形態で述べた所定の空間周波数帯域が含まれていてもよい。また、当該複数の空間周波数帯域には、他の空間周波数帯域と一部が重複する空間周波数帯域が含まれていてもよい。
当該複数の空間周波数帯域は、肌の表面形状の多重構造を形成するいずれか複数の凹凸要素の代表寸法(代表スケール)に基づいて定めることもできる。また、空間周波数帯域の区分けを変えながら肌の評価指標との相関を調査し、その評価指標との相関の高さに基づいて、当該複数の空間周波数帯域が定められてもよい。
The step (S33) is a step of decomposing the unevenness information acquired in the step (S31) into a plurality of spatial frequency bands. That is, the predetermined spatial frequency band decomposed in the step (S33) is at least two bands.
Here, decomposing the unevenness information measured from the skin or the skin replica into the unevenness information of a plurality of spatial frequency bands may be referred to as "multiple decomposition".
The plurality of spatial frequency bands are predetermined before the execution of the step (S33). The plurality of spatial frequency bands may include a predetermined spatial frequency band described in the first embodiment. Further, the plurality of spatial frequency bands may include a spatial frequency band that partially overlaps with other spatial frequency bands.
The plurality of spatial frequency bands can also be determined based on the representative dimensions (representative scale) of any of the plurality of uneven elements forming the multiple structure of the surface shape of the skin. Further, the correlation with the skin evaluation index may be investigated while changing the division of the spatial frequency band, and the plurality of spatial frequency bands may be determined based on the height of the correlation with the evaluation index.
図4は、肌表面の多重構造を示す図である。図5は、第二実施形態における肌表面の凹凸情報の多重分解を示す図である。
図5の例では、工程(S31)で取得された肌表面の凹凸情報が図5の最上部に凹凸波形データとして示されている。工程(S33)でその凹凸波形データが多重分解されることにより、4つの空間周波数帯域の凹凸波形データ(凹凸情報)が抽出されている。具体的には、1/28(cycle/μm)以上の空間周波数帯域、1/80(cycle/μm)以上1/28(cycle/μm)未満の空間周波数帯域、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域、1/500(cycle/μm)未満の空間周波数帯域の4つの空間周波数帯域に多重分解されている。なお、図5の各グラフは、空間周波数帯域に対応する波長範囲で示されている。
FIG. 4 is a diagram showing a multi-layered structure on the skin surface. FIG. 5 is a diagram showing multiple decompositions of unevenness information on the skin surface in the second embodiment.
In the example of FIG. 5, the unevenness information of the skin surface acquired in the step (S31) is shown as unevenness waveform data at the uppermost part of FIG. By multiplexing the uneven waveform data in the step (S33), the concave-convex waveform data (concavo-convex information) of the four spatial frequency bands is extracted. Specifically, a spatial frequency band of 1/28 (cycle / μm) or more, a spatial frequency band of 1/80 (cycle / μm) or more and less than 1/28 (cycle / μm), 1/500 (cycle / μm) It is multiplexed into four spatial frequency bands: a spatial frequency band of less than 1/80 (cycle / μm) and a spatial frequency band of less than 1/500 (cycle / μm). Each graph of FIG. 5 is shown in a wavelength range corresponding to the spatial frequency band.
ここで、肌表面の多重構造を形成する凹凸要素には、図4に示されるように、皮丘及び皮溝からなるキメ、キメよりも大寸法のキメの集合体や毛穴、及びキメよりも小寸法の二次レリーフが少なくとも含まれる。これら3つの凹凸要素を考慮した場合、図5に例示される4つの凹凸波形のうち、1/28(cycle/μm)未満1/80(cycle/μm)以上の空間周波数帯域の凹凸波形を二次レリーフの代表寸法に対応付け、1/80(cycle/μm)未満1/500(cycle/μm)以上の空間周波数帯域の凹凸波形をキメの代表寸法に対応付け、1/500(cycle/μm)未満の空間周波数帯域の凹凸波形はキメの集合体や毛穴の代表寸法に対応付けることができる。
これにより、図5に示される多重分解の例では、肌表面の多重構造を形成する各凹凸要素に対応する凹凸情報が抽出されていると考えることができる。
Here, as shown in FIG. 4, the uneven elements forming the multiple structure on the skin surface include textures composed of skin hills and skin grooves, aggregates and pores of textures having a size larger than the textures, and textures. At least a small secondary relief is included. When these three concave-convex elements are taken into consideration, of the four concave-convex waveforms exemplified in FIG. 5, two concave-convex waveforms in the spatial frequency band of less than 1/28 (cycle / μm) and 1/80 (cycle / μm) or more are obtained. Corresponding to the representative dimensions of the next relief, the uneven waveform in the spatial frequency band of less than 1/80 (cycle / μm) and 1/500 (cycle / μm) or more is associated with the representative dimensions of the texture, 1/500 (cycle / μm). Concavo-convex waveforms in the spatial frequency band less than) can be associated with representative dimensions of texture aggregates and pores.
As a result, in the example of multiple decomposition shown in FIG. 5, it can be considered that the unevenness information corresponding to each unevenness element forming the multiple structure on the skin surface is extracted.
工程(S35)は、工程(S33)で多重分解されたいずれか複数の空間周波数帯域の凹凸情報から複数種の統計値を算出する工程である。
工程(S35)で算出される複数種の統計値は、第一実施形態と同様に、肌表面の高さ方向の偏差に基づくパラメータである。ここでは、平均、標準偏差、分散、歪度、又は尖度のいずれか複数の組合せが複数種の統計値として算出される。
また、算出される当該複数種の統計値は、工程(S33)で多重分解された全ての空間周波数帯域の凹凸情報から算出されてもよいし、一部の複数の空間周波数帯域の凹凸情報から算出されてもよい。つまり、工程(S35)では、平均、標準偏差、分散、歪度、又は尖度を多重分解された各帯域のうち一部の帯域において求めてもよい。
また、各統計値は、統計値の種別ごとに異なる空間周波数帯域の凹凸情報からそれぞれ算出されてもよい。例えば、所定の第一空間周波数帯域の凹凸情報から歪度が算出され、第一空間周波数帯域とは少なくとも一部の帯域が異なる所定の第二空間周波数帯域の凹凸情報から平均、分散(又は標準偏差)、尖度が算出されてもよい。但し、算出される当該複数種の統計値の一部は、共通の空間周波数帯域の凹凸情報から算出されてもよい。
The step (S35) is a step of calculating a plurality of types of statistical values from the unevenness information of any one or more of the plurality of spatial frequency bands decomposed in the step (S33).
The plurality of types of statistical values calculated in the step (S35) are parameters based on the deviation in the height direction of the skin surface, as in the first embodiment. Here, a plurality of combinations of mean, standard deviation, variance, skewness, or kurtosis are calculated as a plurality of kinds of statistical values.
Further, the calculated statistical values of the plurality of types may be calculated from the unevenness information of all the spatial frequency bands multiple-decomposed in the step (S33), or from the unevenness information of some of the plurality of spatial frequency bands. It may be calculated. That is, in the step (S35), the mean, standard deviation, variance, skewness, or kurtosis may be obtained in a part of the multiple-decomposed bands.
Further, each statistical value may be calculated from the unevenness information of the spatial frequency band different for each type of statistical value. For example, the skewness is calculated from the unevenness information of the predetermined first spatial frequency band, and the average, variance (or standard) is calculated from the unevenness information of the predetermined second spatial frequency band whose band is different from that of the first spatial frequency band. Deviation) and skewness may be calculated. However, a part of the calculated statistical values of the plurality of types may be calculated from the unevenness information of the common spatial frequency band.
本発明者らは、実施例で例示するように、肌の複数の評価指標の各々について、どの空間周波数帯域のどの種の統計値と相関が高いかを調査し、肌の評価指標と高い相関を示す所定の複数の空間周波数帯域の複数種の統計値の存在を明らかにした。
これにより、第二実施形態では、このように見出された肌評価指標と各空間周波数帯域の凹凸情報の統計値との相関性に基づいて、どの空間周波数帯域の凹凸情報からどの種の統計値を算出すべきかが定められることが望ましい。
具体的には、実施例によれば、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域の凹凸波形データから歪度が算出され、1/500(cycle/μm)未満の空間周波数帯域の凹凸波形データから平均、分散(又は標準偏差)が算出される。この例で示されるように、工程(S35)で複数種の統計値の算出に用いられるいずれか複数の空間周波数帯域は、上限周波数と下限周波数との比率が相互に異なることが望ましい。
As illustrated in Examples, the present inventors investigated which kind of statistical value in which spatial frequency band is highly correlated with each of a plurality of skin evaluation indexes, and highly correlated with the skin evaluation index. We have clarified the existence of multiple types of statistical values in a plurality of predetermined spatial frequency bands.
As a result, in the second embodiment, based on the correlation between the skin evaluation index found in this way and the statistical value of the unevenness information of each spatial frequency band, what kind of statistics are obtained from the unevenness information of which spatial frequency band. It is desirable to determine whether the value should be calculated.
Specifically, according to the embodiment, the skewness is calculated from the uneven waveform data of the spatial frequency band of 1/500 (cycle / μm) or more and less than 1/80 (cycle / μm), and the skewness is calculated to be 1/500 (cycle / μm). The mean and variance (or standard deviation) are calculated from the uneven waveform data in the spatial frequency band of less than μm). As shown in this example, it is desirable that the ratios of the upper limit frequency and the lower limit frequency of any one or more of the spatial frequency bands used in the calculation of the plurality of kinds of statistical values in the step (S35) are different from each other.
工程(S37)は、工程(S35)で算出された複数種の統計値を用いて対象肌を評価する工程である。工程(S37)における評価手法や実行主体は、複数種の統計値を用いることを除き、第一実施形態における工程(S17)と同様である。例えば、上述の比較表が統計値の種別ごとに設けられていればよい。
また、工程(S37)では、複数種の統計値を用いて一つの評価指標で対象肌を評価してもよいし、複数の評価指標で対象肌を評価してもよい。前者の場合、例えば、母集団の肌サンプルデータを用いた重回帰分析により予め取得された重回帰式であって、複数種の統計値(例えば、歪度及び尖度)を説明変数とし肌の評価スコア(例えば、滑らかさ度合)を目的変数とする重回帰式に、工程(S35)で算出された複数種の統計値を代入することで、肌の評価スコアを算出することができる。
後者の場合、例えば、統計値の種別ごとに肌の評価指標の種別が対応付けられていることで、複数種の指標値を用いて複数の指標で対象肌を評価することができる。
The step (S37) is a step of evaluating the target skin using a plurality of types of statistical values calculated in the step (S35). The evaluation method and the execution subject in the step (S37) are the same as those in the step (S17) in the first embodiment except that a plurality of types of statistical values are used. For example, the above-mentioned comparison table may be provided for each type of statistical value.
Further, in the step (S37), the target skin may be evaluated by one evaluation index using a plurality of types of statistical values, or the target skin may be evaluated by a plurality of evaluation indexes. In the former case, for example, it is a multiple regression equation obtained in advance by multiple regression analysis using the skin sample data of the population, and multiple kinds of statistical values (for example, skewness and kurtosis) are used as explanatory variables for the skin. The skin evaluation score can be calculated by substituting the plurality of types of statistical values calculated in the step (S35) into the multiple regression equation with the evaluation score (for example, the degree of smoothness) as the objective variable.
In the latter case, for example, by associating the types of skin evaluation indexes with each type of statistical value, it is possible to evaluate the target skin with a plurality of indexes using a plurality of types of index values.
このように、第二実施形態では、対象肌又はそのレプリカから計測された肌表面の凹凸情報が複数の空間周波数帯域に多重分割され、いずれか複数の空間周波数帯域の凹凸情報から算出される複数種の統計値を用いて対象肌が評価される。つまり、第二実施形態によれば、多重構造を有する肌の表面形状を複数の凹凸要素の特徴に分解してそれらを個別に参照することができる。そのように個別に参照した凹凸要素の各特徴に基づいて肌評価することで、肌の評価精度を向上させることができる。 As described above, in the second embodiment, the unevenness information of the skin surface measured from the target skin or a replica thereof is divided into a plurality of spatial frequency bands, and the plurality of irregularities calculated from the unevenness information of any one of a plurality of spatial frequency bands. Target skin is evaluated using species statistics. That is, according to the second embodiment, the surface shape of the skin having a multiple structure can be decomposed into the features of a plurality of uneven elements and can be referred to individually. By evaluating the skin based on each feature of the uneven element referred to individually in this way, the evaluation accuracy of the skin can be improved.
第一実施形態及び第二実施形態に係る肌評価方法は、次のような肌評価装置により実行され得る。以下、肌評価装置について図6及び図7を用いて説明する。
図6は、肌評価装置10のハードウェア構成例を概念的に示す図であり、図7は、肌評価装置10の処理構成例を概念的に示す図である。
The skin evaluation method according to the first embodiment and the second embodiment can be executed by the following skin evaluation device. Hereinafter, the skin evaluation device will be described with reference to FIGS. 6 and 7.
FIG. 6 is a diagram conceptually showing a hardware configuration example of the
肌評価装置10は、いわゆるコンピュータであり、例えば、バスで相互に接続される、CPU(Central Processing Unit)11、メモリ12、入出力インタフェース(I/F)13、通信ユニット14等を有する。肌評価装置10を形成する各ハードウェア要素の数はそれぞれ制限されず、これらハードウェア要素は情報処理回路と総称することもできる。また、肌評価装置10は、図6に図示されないハードウェア要素を含んでもよく、そのハードウェア構成は制限されない。
The
CPU11は、一般的なCPU以外に、特定用途向け集積回路(ASIC)、DSP(Digital Signal Processor)、GPU(Graphics Processing Unit)等で構成してもよい。
メモリ12は、RAM(Random Access Memory)、ROM(Read Only Memory)、補助記憶装置(ハードディスク等)である。
入出力I/F13は、出力装置15、入力装置16等のユーザインタフェース装置と接続可能である。出力装置15は、LCD(Liquid Crystal Display)やCRT(Cathode Ray Tube)ディスプレイのような、CPU11等により処理された描画データに対応する画面を表示する装置、印刷装置などの少なくとも一つである。入力装置16は、キーボード、マウス等のようなユーザ操作の入力を受け付ける装置である。出力装置15及び入力装置16は一体化され、タッチパネルとして実現されてもよい。
通信ユニット14は、他のコンピュータとの通信網を介した通信や、他の機器との信号のやりとり等を行う。通信ユニット14には、可搬型記録媒体等も接続され得る。また、通信ユニット14には、対象肌又はそのレプリカから凹凸情報を計測する計測器(レーザ顕微鏡など)が接続されてもよい。
また、肌評価装置10は、その計測器を内蔵する機器であってもよい。
In addition to a general CPU, the
The
The input / output I / F 13 can be connected to a user interface device such as an
The
Further, the
肌評価装置10は、取得部21、分解部22、算出部23、出力処理部24等を有する。これらは、ソフトウェア要素であり、例えば、メモリ12に格納される評価プログラム19がCPU11にロードされ実行されることにより実現される。この評価プログラム19は、例えば、CD(Compact Disc)、メモリカード等のような可搬型記録媒体やネットワーク上の他のコンピュータから入出力I/F13又は通信ユニット14を介してインストールされ、メモリ12に格納されてもよい。
The
取得部21は、工程(S11)又は工程(S31)を実行する。
取得部21は、対象肌又はそのレプリカから計測された肌表面の凹凸情報を計測器、外部のコンピュータ、可搬型記録媒体などから入出力I/F13又は通信ユニット14を介して取得することができる。また、取得部21が計測器として実現されている場合、取得部21が計測しその凹凸情報を生成してもよい。
The
The
分解部22は、工程(S13)又は工程(S33)を実行する。分解部22による所定の空間周波数帯域の凹凸情報の抽出手法及び多重分解手法については上述したとおりである。
算出部23は、工程(S15)又は工程(S35)を実行する。算出部23による統計値の算出手法などについても上述したとおりである。
The disassembling
The
出力処理部24は、算出部23により算出された一種以上の統計値を用いて肌の評価情報を出力する。例えば、出力処理部24は、工程(S17)又は工程(S37)を実行する。出力処理部24により出力される評価情報の具体的内容及びその生成手法は制限されない。
例えば、出力処理部24は、母集団の肌サンプルデータを用いた回帰分析により取得された単回帰式又は重回帰式を保持していてもよい。単回帰式は、いずれか一種の統計値を説明変数とし肌の評価スコアを目的変数とする回帰式であればよく、重回帰式は、複数種の統計値を説明変数とし肌の評価スコアを目的変数とする回帰式であればよい。例えば、出力処理部24は、歪度及び尖度を説明変数とし肌の滑らかさ度合を目的変数とする重回帰式を保持することができる。出力処理部24は、このような回帰式に算出部23により算出された統計値を代入することで、肌の評価スコアを算出し、その評価スコアを肌の評価情報として出力することができる。もちろん、出力される評価情報は、数値ではなく或いは数値に加えて、その数値に対応する評価文字列(「良い」、「悪い」など)や図柄などであってもよい。
また、出力処理部24は、算出部23により算出された統計値と上述の比較表とを肌の評価情報として出力することもできる。この場合には、実際に肌の評価を行うのは、出力された評価情報を参照した人である。
The
For example, the
Further, the
[変形例]
上述の実施形態では、対象肌に関する凹凸形状の一以上の統計値を用いて肌が評価されたが、統計値の変化に基づいて肌の状態変化(状態改善、状態悪化など)を評価することもできる。
実施例で例示されるように、本発明者らは、統計値の変化と肌の或る評価指標における状態変化と間に相関があることを見出している。
そのため、皮膚外用剤の継続使用に伴い、対象肌又はそのレプリカから計測された凹凸情報をそれぞれ取得し、上述の実施形態と同様に一以上の統計値を算出し、その統計値の差に基づいて、肌状態の変化状況を評価することができる。例えば、予めサンプル肌に基づいて肌状態の変化状況(改善、悪化、変化無し)ごとに統計値差の範囲を求めておき、被験者の対象肌に関して得られた統計値差がどの範囲に属するかによって、改善した、悪化した、又は変化無しというように対象肌の変化状況を決定することができる。肌状態の変化状況は、改善、悪化、変化無しといった3値のみでなく、改善及び悪化の2値、或いは、4値以上の変化スコアとして算出されてもよい。
[Modification example]
In the above-described embodiment, the skin is evaluated using one or more statistical values of the uneven shape of the target skin, but the skin condition change (condition improvement, condition deterioration, etc.) is evaluated based on the change in the statistical value. You can also.
As illustrated in the examples, the present inventors have found that there is a correlation between changes in statistical values and changes in the state of skin at a certain evaluation index.
Therefore, with the continuous use of the external preparation for skin, the unevenness information measured from the target skin or its replica is acquired, one or more statistical values are calculated in the same manner as in the above-described embodiment, and the difference between the statistical values is used. Therefore, it is possible to evaluate the change status of the skin condition. For example, the range of the statistical value difference is obtained in advance for each change state (improvement, deterioration, no change) of the skin condition based on the sample skin, and which range the statistical value difference obtained for the target skin of the subject belongs to. It is possible to determine the state of change of the target skin such as improved, deteriorated, or no change. The change state of the skin condition may be calculated not only as a three-value such as improvement, deterioration, and no change, but also as a two-value of improvement and deterioration, or a change score of four or more values.
よって、上述の実施形態に係る肌評価方法を用いて、化粧料、医薬品、医薬部外品などの皮膚外用剤又は医療機器や美容機器等の効果を評価する方法を実現することもできる。言い換えれば、皮膚外用剤を適用する前後の肌又は肌のレプリカから計測された肌表面の凹凸情報をそれぞれ取得する工程と、取得された凹凸情報から所定の空間周波数帯域の凹凸情報を抽出する分解工程と、抽出された所定の空間周波数帯域の凹凸情報から凹凸形状に関する統計値を算出する算出工程と、皮膚外用剤を適用する前後に関して算出された当該統計値の差を算出する工程と、算出された差を用いて皮膚外用剤を適用する前後の肌を評価する評価工程と、を含む皮膚外用剤の効果を評価する方法を実現することができる。 Therefore, it is also possible to realize a method for evaluating the effects of external preparations for skin such as cosmetics, pharmaceuticals, quasi-drugs, medical devices, beauty devices, etc., by using the skin evaluation method according to the above-described embodiment. In other words, the process of acquiring the unevenness information of the skin surface measured from the skin or the skin replica before and after applying the external preparation for skin, and the decomposition of extracting the unevenness information of a predetermined spatial frequency band from the acquired unevenness information. A step, a calculation step of calculating a statistical value regarding the uneven shape from the extracted unevenness information of a predetermined spatial frequency band, and a step of calculating the difference between the statistical values calculated before and after applying the external preparation for skin, and calculation. It is possible to realize an evaluation step of evaluating the skin before and after applying the external preparation for skin using the difference obtained, and a method for evaluating the effect of the external preparation for skin including.
また、上述の説明で用いた複数のフローチャートでは、複数の工程(処理)が順番に記載されているが、各工程の実行順序は、その記載の順番に制限されない。図示される工程の順番を内容的に支障のない範囲で変更することができる。 Further, in the plurality of flowcharts used in the above description, a plurality of steps (processes) are described in order, but the execution order of each step is not limited to the order of description. The order of the illustrated steps can be changed within a range that does not hinder the contents.
上述の内容の一部又は全部は、次のようにも特定され得る。但し、上述の内容が以下の記載に制限されるものではない。 Part or all of the above content may also be specified as follows. However, the above contents are not limited to the following descriptions.
<1>肌又は肌のレプリカから計測された肌表面の凹凸情報を取得する工程と、
前記取得された凹凸情報から所定の空間周波数帯域の凹凸情報を抽出する分解工程と、
前記抽出された所定の空間周波数帯域の凹凸情報から凹凸形状に関する統計値を算出する算出工程と、
前記算出された統計値を用いて前記肌を評価する評価工程と、
を含む肌評価方法。
<1> The process of acquiring the unevenness information of the skin surface measured from the skin or the replica of the skin, and
A decomposition step of extracting unevenness information in a predetermined spatial frequency band from the acquired unevenness information, and
A calculation step of calculating statistical values related to the uneven shape from the extracted uneven information of a predetermined spatial frequency band, and
An evaluation process for evaluating the skin using the calculated statistical values, and
Skin evaluation method including.
<2>前記算出工程で算出される前記統計値は、少なくとも歪度又は尖度を含み、
前記評価工程では、前記算出工程で算出された少なくとも歪度又は尖度を用いて前記肌を評価する、
<1>に記載の肌評価方法。
<3>前記分解工程では、前記所定の空間周波数帯域を少なくとも2帯域とし、
前記算出工程では、平均、標準偏差、分散、歪度、又は尖度を各帯域ごとに求める、
<1>に記載の肌評価方法。
<4>前記分解工程では、前記所定の空間周波数帯域を含む異なる複数の空間周波数帯域に、前記取得された凹凸情報を多重分解し、
前記算出工程では、前記多重分解されたいずれか複数の空間周波数帯域の凹凸情報から複数種の統計値を算出し、
前記評価工程では、前記算出された複数種の統計値を用いて前記肌を評価する、
<1>から<3>のいずれか一つに記載の肌評価方法。
<5>前記分解工程では、異なる所定の第一空間周波数帯域及び所定の第二空間周波数帯域を含む前記複数の空間周波数帯域に、前記取得された凹凸情報を多重分解し、
前記算出工程では、前記所定の第一空間周波数帯域の凹凸情報から歪度を算出し、前記所定の第二空間周波数帯域の凹凸情報から尖度を算出し、
前記評価工程では、前記算出された歪度及び尖度を用いて、前記肌を評価する、
<4>に記載の肌評価方法。
<6>前記評価工程では、母集団の肌サンプルデータを用いた重回帰分析により予め取得された重回帰式であって、前記歪度及び前記尖度を説明変数とし肌の滑らかさ度合を目的変数とする重回帰式を用いて、前記肌の滑らかさ度合を決定する、
<5>に記載の肌評価方法。
<7>前記算出工程で前記複数種の統計値の算出に用いられる前記いずれか複数の空間周波数帯域は、上限周波数と下限周波数との比率が相互に異なる、
<4>から<6>のいずれか一つに記載の肌評価方法。
<8>肌又は肌のレプリカから計測された肌表面の凹凸情報を取得する取得手段と、
前記取得された凹凸情報から所定の空間周波数帯域の凹凸情報を抽出する分解手段と、
前記抽出された所定の空間周波数帯域の凹凸情報から凹凸形状に関する統計値を算出する算出手段と、
前記算出された統計値を用いて前記肌の評価情報を出力する出力処理手段と、
を備える肌評価装置。
<9>前記算出手段は、前記統計値として凹凸形状の少なくとも歪度又は尖度を算出し、
前記出力処理手段は、前記算出手段により算出された少なくとも歪度又は尖度を用いて前記肌の評価情報を生成する、
<8>に記載の肌評価装置。
<10>前記分解手段は、前記所定の空間周波数帯域を含む異なる複数の空間周波数帯域に、前記取得された凹凸情報を多重分解し、
前記算出手段は、前記多重分解されたいずれか複数の空間周波数帯域の凹凸情報から複数種の統計値を算出し、
前記出力処理手段は、前記算出された複数種の統計値を用いて前記肌の評価情報を生成する、
<8>又は<9>に記載の肌評価装置。
<11>前記分解手段は、異なる所定の第一空間周波数帯域及び所定の第二空間周波数帯域を含む前記複数の空間周波数帯域に、前記取得された凹凸情報を多重分解し、
前記算出手段は、前記所定の第一空間周波数帯域の凹凸情報から歪度を算出し、前記所定の第二空間周波数帯域の凹凸情報から尖度を算出し、
前記出力処理手段は、前記算出された歪度及び尖度を用いて、前記肌の評価情報を生成する、
<10>に記載の肌評価装置。
<12>前記出力処理手段は、母集団の肌サンプルデータを用いた重回帰分析により予め取得された重回帰式であって、前記歪度及び前記尖度を説明変数とし肌の滑らかさ度合を目的変数とする重回帰式を用いて、前記肌の滑らかさ度合を決定する、
<11>に記載の肌評価装置。
<13>前記算出手段で前記複数種の統計値の算出に用いられる前記いずれか複数の空間周波数帯域は、上限周波数と下限周波数との比率が相互に異なる、
<10>から<12>のいずれか一つに記載の肌評価装置。
<2> The statistical value calculated in the calculation step includes at least skewness or kurtosis.
In the evaluation step, the skin is evaluated using at least the skewness or kurtosis calculated in the calculation step.
The skin evaluation method according to <1>.
<3> In the decomposition step, the predetermined spatial frequency band is set to at least two bands.
In the calculation step, the mean, standard deviation, variance, skewness, or kurtosis is obtained for each band.
The skin evaluation method according to <1>.
<4> In the decomposition step, the acquired unevenness information is multiple-decomposed into a plurality of different spatial frequency bands including the predetermined spatial frequency band.
In the calculation step, a plurality of types of statistical values are calculated from the unevenness information of any of the plurality of spatial frequency bands obtained by the plural decomposition.
In the evaluation step, the skin is evaluated using the calculated plurality of types of statistical values.
The skin evaluation method according to any one of <1> to <3>.
<5> In the decomposition step, the acquired unevenness information is multiple-decomposed into the plurality of spatial frequency bands including different predetermined first spatial frequency bands and predetermined second spatial frequency bands.
In the calculation step, the skewness is calculated from the unevenness information of the predetermined first spatial frequency band, and the kurtosis is calculated from the unevenness information of the predetermined second spatial frequency band.
In the evaluation step, the skin is evaluated using the calculated skewness and kurtosis.
The skin evaluation method according to <4>.
<6> In the evaluation step, the multiple regression equation obtained in advance by the multiple regression analysis using the skin sample data of the population is used, and the degree of smoothness of the skin is aimed at by using the skewness and the kurtosis as explanatory variables. The degree of skin smoothness is determined by using a multiple regression equation as a variable.
The skin evaluation method according to <5>.
<7> In any one of the plurality of spatial frequency bands used for calculating the plurality of types of statistical values in the calculation step, the ratio of the upper limit frequency and the lower limit frequency is different from each other.
The skin evaluation method according to any one of <4> to <6>.
<8> An acquisition means for acquiring information on the unevenness of the skin surface measured from the skin or a replica of the skin, and
Decomposition means for extracting unevenness information in a predetermined spatial frequency band from the acquired unevenness information, and
A calculation means for calculating statistical values related to the shape of unevenness from the extracted information on unevenness in a predetermined spatial frequency band, and
An output processing means that outputs the evaluation information of the skin using the calculated statistical value, and
A skin evaluation device equipped with.
<9> The calculation means calculates at least the skewness or kurtosis of the uneven shape as the statistical value.
The output processing means generates evaluation information of the skin using at least the skewness or kurtosis calculated by the calculation means.
The skin evaluation device according to <8>.
<10> The decomposition means multiplexizes the acquired unevenness information into a plurality of different spatial frequency bands including the predetermined spatial frequency band.
The calculation means calculates a plurality of types of statistical values from the unevenness information of any of the plurality of spatial frequency bands obtained by the plural decomposition.
The output processing means generates the evaluation information of the skin using the calculated plurality of kinds of statistical values.
The skin evaluation device according to <8> or <9>.
<11> The decomposition means multiplexizes the acquired unevenness information into the plurality of spatial frequency bands including different predetermined first spatial frequency bands and predetermined second spatial frequency bands.
The calculation means calculates the skewness from the unevenness information of the predetermined first spatial frequency band, calculates the kurtosis from the unevenness information of the predetermined second spatial frequency band, and calculates the kurtosis.
The output processing means uses the calculated skewness and kurtosis to generate the skin evaluation information.
The skin evaluation device according to <10>.
<12> The output processing means is a multiple regression equation previously acquired by multiple regression analysis using population skin sample data, and the degree of smoothness of the skin is determined by using the skewness and the kurtosis as explanatory variables. The degree of skin smoothness is determined by using the multiple regression equation as the objective variable.
The skin evaluation device according to <11>.
<13> In any of the plurality of spatial frequency bands used in the calculation means for calculating the plurality of types of statistical values, the ratio of the upper limit frequency and the lower limit frequency is different from each other.
The skin evaluation device according to any one of <10> to <12>.
以下に実施例を挙げ、上述の内容を更に詳細に説明する。但し、以下の実施例の記載は、上述の内容に何ら限定を加えるものではない。 Examples will be given below, and the above contents will be described in more detail. However, the description of the following examples does not impose any limitation on the above contents.
本実施例では、上述の実施形態の効果を検証すべく、複数人の被験者の肌をサンプル肌として、各サンプル肌のレプリカからそれぞれ計測された凹凸情報に基づいて、各サンプル肌の見た目(外観)及び触感と肌の凹凸形状に関する統計値との関係について調査された。 In this embodiment, in order to verify the effect of the above-described embodiment, the skins of a plurality of subjects are used as sample skins, and the appearance (appearance) of each sample skin is based on the unevenness information measured from each sample skin replica. ) And the relationship between the tactile sensation and the statistical values regarding the uneven shape of the skin were investigated.
各サンプル肌のレプリカは、親水性ビニルシリコン印象剤を各被験者の鼻横から頬にかけての部位に押し付けることで生成された。このため、各サンプル肌のレプリカには、各被験者の対象部位の皮膚表面の凹凸形状が反転された状態で転写された。
このように採取された各サンプル肌のレプリカの肌転写面の三次元形状をレーザ顕微鏡で10倍の倍率でそれぞれ計測し、各サンプル肌に関し高さ情報(凹凸情報)を色で表す肌画像がそれぞれ生成された。更に、各サンプル肌の肌画像に対して約7mmのラインを指定して、そのライン上の高さと位置との関係を示す凹凸波形データが取得された。
Replicas of each sample skin were generated by pressing a hydrophilic vinyl silicone impression agent against the area from the side of the nose to the cheeks of each subject. Therefore, it was transferred to the replica of each sample skin in a state where the uneven shape of the skin surface of the target site of each subject was inverted.
The three-dimensional shape of the skin transfer surface of the replica of each sample skin collected in this way is measured with a laser microscope at a magnification of 10 times, and a skin image showing height information (unevenness information) for each sample skin in color is obtained. Each was generated. Further, a line of about 7 mm was designated for the skin image of each sample skin, and uneven waveform data showing the relationship between the height and the position on the line was acquired.
図8は、3つのサンプル肌の肌画像と各肌画像から比較例として取得された凹凸形状に関する統計値との関係を示す図である。
比較例では、上述の実施形態のような多重分解を行うことなく、全体の凹凸波形データから凹凸形状に関する統計値が算出された。具体的には、サンプル肌A、B及びCの各肌画像に対して20本のラインをそれぞれ指定することで、サンプル肌ごとに20個の凹凸波形データを取得し、各凹凸波形データの統計値の平均値が図8に示されている。統計値としては、図8に示されるとおり、算術平均粗さ(Ra)、二乗平均平方根粗さ(Rq)、歪度(スキューネス)(Rsk)、尖度(クルトシス)(Rku)が算出された。
比較例として算出された統計値と各サンプル肌の肌画像とを比較すると、いずれの統計値においても肌の凹凸形状を的確には指標化できていないことが分かる。特に、肌画像で示されているように、サンプル肌Aとサンプル肌Bとは明らかに表面形状が異なっているにも関わらず、いずれの統計値もそのような表面形状の違いを表しているとは言えない。
FIG. 8 is a diagram showing the relationship between the skin images of the three sample skins and the statistical values regarding the uneven shape obtained from each skin image as a comparative example.
In the comparative example, the statistical value regarding the uneven shape was calculated from the entire uneven waveform data without performing the multiple decomposition as in the above-described embodiment. Specifically, by designating 20 lines for each skin image of sample skins A, B, and C, 20 uneven waveform data are acquired for each sample skin, and statistics of each uneven waveform data are obtained. The average value of the values is shown in FIG. As statistical values, as shown in FIG. 8, arithmetic mean roughness (Ra), root mean square roughness (Rq), skewness (skewness) (Rsk), and kurtosis (Rku) were calculated. ..
Comparing the statistical values calculated as comparative examples with the skin images of each sample skin, it can be seen that the uneven shape of the skin cannot be accurately indexed in any of the statistical values. In particular, as shown in the skin image, although the surface shapes of the sample skin A and the sample skin B are clearly different, all the statistical values show such a difference in the surface shape. It can not be said.
本実施例では、上述の複数人の被験者に共通の皮膚外用剤を利用させ、皮膚外用剤の使用前の肌と使用後の肌とをそれぞれサンプル肌とし、各サンプル肌のレプリカから上述したとおりレーザ顕微鏡により凹凸波形データがそれぞれ取得された。そして、各サンプル肌の凹凸波形データを4つの空間周波数帯域に多重分解することで、4つの各空間周波数帯域の凹凸波形データを抽出し、空間周波数帯域ごとに4種の指標値(Ra、Rq、Rsk、Rku)が算出された。その上で、皮膚外用剤の使用前のサンプル肌と使用後のサンプル肌とに関してそれぞれ算出された指標値の差が被験者ごとに算出された。結果として、各被験者について、4つの空間周波数帯域の各々における4種の指標値差がそれぞれ取得された。
4つの空間周波数帯域は、1/28(cycle/μm)以上の空間周波数帯域、1/80(cycle/μm)以上1/28(cycle/μm)未満の空間周波数帯域、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域、1/500(cycle/μm)未満の空間周波数帯域である。
In this embodiment, the above-mentioned plurality of subjects are made to use a common skin external preparation, and the skin before and after the use of the skin external preparation is used as sample skin, respectively, and as described above from the replica of each sample skin. Concavo-convex waveform data was acquired by a laser microscope. Then, by multiplexing the uneven waveform data of each sample skin into four spatial frequency bands, the uneven waveform data of each of the four spatial frequency bands is extracted, and four kinds of index values (Ra, Rq) are extracted for each spatial frequency band. , Rsk, Rku) were calculated. Then, the difference between the index values calculated for the sample skin before the use of the external preparation for skin and the sample skin after the use was calculated for each subject. As a result, for each subject, four kinds of index value differences in each of the four spatial frequency bands were acquired.
The four spatial frequency bands are a spatial frequency band of 1/28 (cycle / μm) or more, a spatial frequency band of 1/80 (cycle / μm) or more and less than 1/28 (cycle / μm), and 1/500 (cycle / μm). It is a spatial frequency band of μm) or more and less than 1/80 (cycle / μm), and a spatial frequency band of less than 1/500 (cycle / μm).
一方で、本実施例では、当該複数の被験者における皮膚外用剤の使用前の肌と使用後の肌とについて官能評価により見た目(外観)及び触感がスコア化された。
見た目の評価は、専門の評価者がサンプル肌を目視することで、小じわの目立ち具合を10段階(小じわスコアと表記される)で官能評価する方法で行われた。ここでの「小じわ」とは、保湿を加えると消えるようなシワを意味する。
触感の評価は、専門の評価者がサンプル肌を触診し、滑らかさ度合を5段階(触感スコアと表記される)で官能評価する方法で行われた。
更に、同一被験者にける皮膚外用剤の使用前の肌の評価スコアと使用後の肌の評価スコアとの間の差が算出された。具体的には、使用前の肌の小じわスコアと使用後の肌の小じわスコアとの差が算出され、使用前の肌の触感スコアと使用後の肌の触感スコアとの差が算出された。
On the other hand, in this example, the appearance (appearance) and tactile sensation of the skin before and after the use of the external preparation for skin in the plurality of subjects were scored by sensory evaluation.
The appearance was evaluated by a method in which a professional evaluator visually evaluated the sample skin to perform a sensory evaluation of the conspicuousness of fine wrinkles on a scale of 10 (referred to as a fine wrinkle score). The term "fine wrinkles" here means wrinkles that disappear when moisturized.
The tactile sensation was evaluated by a method in which a professional evaluator palpates the sample skin and sensory evaluates the degree of smoothness on a scale of 5 (referred to as tactile sensation score).
Furthermore, the difference between the evaluation score of the skin before the use of the external preparation for skin and the evaluation score of the skin after the use of the external preparation for the same subject was calculated. Specifically, the difference between the fine wrinkle score of the skin before use and the fine wrinkle score of the skin after use was calculated, and the difference between the tactile sensation score of the skin before use and the tactile sensation score of the skin after use was calculated.
図9(a)は、4つの空間周波数帯域の凹凸成分及び全ての凹凸成分の各々に関する算術平均粗さ(Ra)と小じわとの関係を示すグラフであり、図9(b)は、4つの空間周波数帯域の凹凸成分及び全ての凹凸成分の各々に関する二乗平均平方根粗さ(Rq)と小じわとの関係を示すグラフである。図9(a)及び図9(b)には、皮膚外用剤の使用後の肌に関する情報が示されている。
図9において、横軸は4つの空間周波数帯域の凹凸成分(対応する波長範囲を表記)と全ての凹凸成分(「全体」と表記)を示し、縦軸は各統計値を示す。「良い」の棒グラフは小じわが標準より目立たない被験者の群を示し、「標準」の棒グラフは小じわの目立ちが標準であった被験者の群を示し、「悪い」の棒グラフは小じわが標準より目立つ被験者の群を示す。各棒グラフの長さが被験者群ごとの統計値の平均を示し、各棒グラフに重畳されている縦線は被験者群ごとの統計値の標準誤差を示す。
FIG. 9A is a graph showing the relationship between the arithmetic mean roughness (Ra) and fine wrinkles for each of the unevenness components of the four spatial frequency bands and all the unevenness components, and FIG. 9B is a graph showing four. It is a graph which shows the relationship between the mean square root roughness (Rq) and fine wrinkles about each of the unevenness component of a spatial frequency band and all unevenness components. 9 (a) and 9 (b) show information about the skin after use of the external preparation for skin.
In FIG. 9, the horizontal axis shows the unevenness component of the four spatial frequency bands (indicated by the corresponding wavelength range) and all the unevenness components (indicated as “whole”), and the vertical axis shows each statistical value. The "good" bar graph shows the group of subjects whose fine wrinkles are less noticeable than the standard, the "standard" bar graph shows the group of subjects whose fine wrinkles are more noticeable than the standard, and the "bad" bar graph shows the subjects whose fine wrinkles are more noticeable than the standard. Shows a group of. The length of each bar graph indicates the average of the statistical values for each subject group, and the vertical line superimposed on each bar graph indicates the standard error of the statistical values for each subject group.
図9(a)及び図9(b)によれば、1/500(cycle/μm)未満の空間周波数帯域の凹凸成分に関する、算術平均粗さ(Ra)、二乗平均平方根粗さ(Rq)が小じわの目立ちやすさと相関があることが確認される。
これにより、1/500(cycle/μm)未満の空間周波数帯域の凹凸成分に関する、算術平均粗さ(Ra)及び二乗平均平方根粗さ(Rq)の統計値差により、被験者の肌の見た目の変化を推定できることが実証されている。
According to FIGS. 9 (a) and 9 (b), the arithmetic mean roughness (Ra) and the root mean square roughness (Rq) with respect to the uneven component of the spatial frequency band of less than 1/500 (cycle / μm) are It is confirmed that there is a correlation with the conspicuity of fine wrinkles.
As a result, the appearance of the subject's skin changes due to the difference in statistical values of the arithmetic mean roughness (Ra) and the root mean square roughness (Rq) for the uneven component in the spatial frequency band of less than 1/500 (cycle / μm). Has been proven to be able to be estimated.
図10は4つの空間周波数帯域の凹凸成分及び全ての凹凸成分の各々に関する歪度(Rsk)と触感(滑らかさ)との関係を示すグラフである。図10には、皮膚外用剤の使用後の肌に関する情報が示されている。
図10において、横軸は4つの空間周波数帯域の凹凸成分(対応する波長範囲を表記)と全ての凹凸成分(「全体」と表記)を示し、縦軸は各統計値を示す。「良い」の棒グラフは触感(滑らかさ)が標準より優れる被験者の群を示し、「標準」の棒グラフは触感(滑らかさ)が標準であった被験者の群を示し、「悪い」の棒グラフは触感(滑らかさ)が標準より悪い被験者の群を示す。各棒グラフの長さが被験者群ごとの統計値の平均を示し、各棒グラフに重畳されている縦線は被験者群ごとの統計値の標準誤差を示す。
FIG. 10 is a graph showing the relationship between the skewness (Rsk) and the tactile sensation (smoothness) for each of the uneven components of the four spatial frequency bands and all the uneven components. FIG. 10 shows information about the skin after the use of the external preparation for skin.
In FIG. 10, the horizontal axis indicates the unevenness component of the four spatial frequency bands (indicated by the corresponding wavelength range) and all the unevenness components (indicated as “whole”), and the vertical axis indicates each statistical value. The "good" bar graph shows the group of subjects with better than standard tactile sensation (smoothness), the "standard" bar graph shows the group of subjects with standard tactile sensation (smoothness), and the "bad" bar graph shows the group of subjects with standard tactile sensation (smoothness). Shows a group of subjects whose (smoothness) is worse than normal. The length of each bar graph indicates the average of the statistical values for each subject group, and the vertical line superimposed on each bar graph indicates the standard error of the statistical values for each subject group.
図10によれば、1/80(cycle/μm)未満1/500(cycle/μm)以上の空間周波数帯域の凹凸成分に関する歪度が触感スコアと相関があることが確認される。
これにより、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域の凹凸成分に関する歪度により、被験者の肌の触感(滑らかさ)を推定できることが実証されている。
更に言えば、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域の凹凸成分に関する歪度が小さい程、触感(滑らかさ度合)は良くなるということも新たに見出されている。但し、歪度が大きいからといって必ずしも触感が良いと判断できない可能性もある。よって、肌の凹凸形状に他の統計値(例えば、平均、分散(及び標準偏差)、尖度など)を更に考慮したほうがより高精度に触感を評価できる場合がある。
According to FIG. 10, it is confirmed that the skewness of the unevenness component in the spatial frequency band of less than 1/80 (cycle / μm) and 1/500 (cycle / μm) or more correlates with the tactile score.
This demonstrates that the tactile sensation (smoothness) of the subject's skin can be estimated from the skewness of the unevenness component in the spatial frequency band of 1/500 (cycle / μm) or more and less than 1/80 (cycle / μm). ..
Furthermore, the smaller the skewness of the uneven component in the spatial frequency band of 1/500 (cycle / μm) or more and less than 1/80 (cycle / μm), the better the tactile sensation (smoothness). Has been found. However, even if the skewness is large, it may not always be judged that the tactile sensation is good. Therefore, it may be possible to evaluate the tactile sensation with higher accuracy by further considering other statistical values (for example, mean, variance (and standard deviation), kurtosis, etc.) for the uneven shape of the skin.
図11(a)は4つの空間周波数帯域の凹凸成分及び全ての凹凸成分の各々に関する歪度(Rsk)の変化と触感の変化との関係を示すグラフであり、図11(b)は4つの空間周波数帯域の凹凸成分及び全ての凹凸成分の各々に関する尖度(Rku)の変化と触感の変化との関係を示すグラフである。
図11において、横軸は4つの空間周波数帯域の凹凸成分(対応する波長範囲を表記)と全ての凹凸成分(「全体」と表記)を示し、縦軸は各統計値の差(使用後の統計値から使用前の統計値を減算した値)を示す。「+」の棒グラフは触感(滑らかさ)が改善した被験者群を示し、「±」の棒グラフは触感の変化がなかった被験者群を示し、「−」の棒グラフは触感が悪化した被験者群を示す。各棒グラフの長さが被験者群ごとの統計値差の平均を示し、各棒グラフに重畳されている縦線は被験者群ごとの統計値差の標準誤差を示す。
FIG. 11 (a) is a graph showing the relationship between the change in skewness (Rsk) and the change in tactile sensation for each of the uneven components of the four spatial frequency bands and all the uneven components, and FIG. 11 (b) shows four. It is a graph which shows the relationship between the change of the kurtosis (Rku) and the change of the tactile sensation about each of the uneven component of the spatial frequency band and all the uneven components.
In FIG. 11, the horizontal axis shows the unevenness component of the four spatial frequency bands (indicated by the corresponding wavelength range) and all the unevenness components (indicated as “whole”), and the vertical axis shows the difference between the statistical values (after use). The value obtained by subtracting the statistical value before use from the statistical value) is shown. The "+" bar graph indicates the subject group with improved tactile sensation (smoothness), the "±" bar graph indicates the subject group with no change in tactile sensation, and the "-" bar graph indicates the subject group with deteriorated tactile sensation. .. The length of each bar graph indicates the average of the statistical value differences for each subject group, and the vertical lines superimposed on each bar graph indicate the standard error of the statistical value differences for each subject group.
図11(a)及び図11(b)によれば、1/80(cycle/μm)未満1/500(cycle/μm)以上の空間周波数帯域の凹凸成分に関する歪度の変化、及び、1/500(cycle/μm)未満の空間周波数帯域の凹凸成分に関する尖度の変化が触感スコアの変化と相関があることが確認される。
歪度(Rsk)の統計値差については、触感が改善した被験者群が負の値を示しており、触感が悪化した被験者群が正の値を示している。これは、触感が改善した被験者群については、皮膚外用剤の使用後に歪度が小さくなっていることを示し、触感が悪化した被験者群については、逆に歪度が大きくなっていることを示している。即ち、肌表面の凹凸形状の歪度が小さい程、滑らかさは向上するという論理に即した結果が得られている。
一方で、尖度(Rku)の統計値差については、歪度の逆の状態となっている。即ち、触感が改善した被験者群については、皮膚外用剤の使用後に尖度が大きくなっていることを示し、触感が悪化した被験者群については、逆に尖度が小さくなっていることを示している。
これにより、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域の凹凸成分に関する歪度の変化、及び、1/500(cycle/μm)未満の空間周波数帯域の凹凸成分に関する尖度の変化により、被験者の肌の触感の変化を高精度に推定できることが実証されている。
更に言えば、1/500(cycle/μm)以上1/80(cycle/μm)未満の空間周波数帯域の凹凸成分に関する歪度が小さい程、触感(滑らかさ度合)は良くなり、1/500(cycle/μm)未満の空間周波数帯域の凹凸成分に関する尖度が大きい程、触感(滑らかさ度合)は良くなるということも新たに見出されている。但し、尖度が大きいからといって必ずしも触感が良いと判断できない可能性もある。よって、肌の凹凸形状に関する尖度で触感を評価する場合には、他の統計値(例えば、Rqや分散値、歪度など)を更に考慮したほうがより高精度に触感を評価できる場合がある。
更に、当該尖度の変化[Rku(500μm−)]及び当該歪度の変化[Rsk(80−500μm)]と触感スコアの変化との関係を重回帰分析することにより、触感スコアの変化を当該尖度の変化と当該歪度の変化とで推定できることが確認された。その重回帰分析により次のような重回帰式が算出された(決定係数:0.146)。
触感スコアの変化=4.01×[Rsk(80−500μm)]−0.73×[Rku(500μm−)−0.25
According to FIGS. 11 (a) and 11 (b), the change in skewness with respect to the unevenness component in the spatial frequency band of less than 1/80 (cycle / μm) and 1/500 (cycle / μm) or more, and 1 /. It is confirmed that the change in kurtosis with respect to the uneven component in the spatial frequency band of less than 500 (cycle / μm) correlates with the change in the tactile score.
Regarding the difference in the statistical values of the skewness (Rsk), the group of subjects with improved tactile sensation showed a negative value, and the group of subjects with deteriorated tactile sensation showed a positive value. This indicates that the skewness of the subjects with improved tactile sensation decreased after the use of the external preparation for skin, and conversely, the skewness of the subjects with deteriorated tactile sensation increased. ing. That is, the result is based on the logic that the smaller the skewness of the uneven shape on the skin surface, the better the smoothness.
On the other hand, the difference in the statistical value of kurtosis (Rku) is the opposite of the skewness. That is, for the group of subjects with improved tactile sensation, it was shown that the kurtosis was increased after the use of the external preparation for skin, and for the group of subjects with deteriorated tactile sensation, conversely, it was shown that the kurtosis was decreased. There is.
As a result, the change in skewness related to the unevenness component of the spatial frequency band of 1/500 (cycle / μm) or more and less than 1/80 (cycle / μm), and the spatial frequency band of less than 1/500 (cycle / μm) It has been demonstrated that the change in the tactile sensation of the subject's skin can be estimated with high accuracy by the change in the kurtosis related to the unevenness component.
Furthermore, the smaller the skewness of the uneven component in the spatial frequency band of 1/500 (cycle / μm) or more and less than 1/80 (cycle / μm), the better the tactile sensation (smoothness), and the better the tactile sensation (smoothness). It has also been newly found that the greater the kurtosis of the uneven component in the spatial frequency band of less than cycle / μm), the better the tactile sensation (degree of smoothness). However, even if the kurtosis is large, it may not always be judged that the tactile sensation is good. Therefore, when evaluating the tactile sensation based on the kurtosis of the uneven shape of the skin, it may be possible to evaluate the tactile sensation with higher accuracy by further considering other statistical values (for example, Rq, dispersion value, skewness, etc.). ..
Furthermore, by performing multiple regression analysis of the relationship between the change in kurtosis [Rku (500 μm-)] and the change in skewness [Rsk (80-500 μm)] and the change in tactile score, the change in tactile score is measured. It was confirmed that it can be estimated from the change in kurtosis and the change in skewness. The following multiple regression equation was calculated by the multiple regression analysis (coefficient of determination: 0.146).
Change in tactile score = 4.01 x [Rsk (80-500 μm)] -0.73 x [Rku (500 μm-) -0.25
本実施例では、1/28(cycle/μm)未満1/80(cycle/μm)以上の空間周波数帯域の凹凸波形を二次レリーフの代表寸法に対応付け、1/80(cycle/μm)未満1/500(cycle/μm)以上の空間周波数帯域の凹凸波形をキメの代表寸法に対応付け、1/500(cycle/μm)未満の空間周波数帯域の凹凸波形はキメの集合体や毛穴の代表寸法に対応付けている。これらの空間周波数帯域を分けて、平均、標準偏差、分散、歪度、若しくは尖度で表現することで、各寸法の肌の形状を評価することができることが検証された。 In this embodiment, the uneven waveform in the spatial frequency band of less than 1/28 (cycle / μm) and 1/80 (cycle / μm) or more is associated with the representative dimension of the secondary relief, and is less than 1/80 (cycle / μm). Concavo-convex waveforms in the spatial frequency band of 1/500 (cycle / μm) or more are associated with typical texture dimensions, and uneven waveforms in the spatial frequency band of less than 1/500 (cycle / μm) are representative of texture aggregates and pores. Corresponds to dimensions. It was verified that the shape of the skin of each dimension can be evaluated by dividing these spatial frequency bands and expressing them by mean, standard deviation, variance, skewness, or kurtosis.
10 肌評価装置
11 CPU
12 メモリ
13 入出力I/F
14 通信ユニット
15 出力装置
16 入力装置
19 評価プログラム
21 取得部
22 分解部
23 算出部
24 出力処理部
10
12 Memory 13 I / O I / F
14
Claims (6)
肌又は肌のレプリカから計測された肌表面の凹凸情報を取得する工程と、
前記取得された凹凸情報を多重分解して、相互に異なる複数の空間周波数帯域の凹凸情報を抽出する分解工程と、
前記多重分解されたいずれか複数の空間周波数帯域の凹凸情報から、凹凸形状に関する統計値であって空間周波数帯域ごとに異なる種別の複数種の統計値を算出する算出工程と、
前記肌の評価のために前記算出された複数種の統計値を出力する出力工程と、
を含む肌評価方法。 A skin evaluation method performed by a computer equipped with at least a CPU (Central Processing Unit) and memory.
The process of acquiring the unevenness information of the skin surface measured from the skin or skin replica,
A decomposition step of multiplying the acquired unevenness information to extract unevenness information of a plurality of spatial frequency bands different from each other, and a decomposition step.
A calculation step of calculating a plurality of types of statistical values related to the uneven shape, which are different types for each spatial frequency band, from the unevenness information of any one of the plurality of spatial frequency bands that have been multiple-decomposed.
An output process that outputs the calculated plurality of types of statistical values for the evaluation of the skin, and
Skin evaluation method including.
請求項1に記載の肌評価方法。 Wherein the plurality of types of statistical values output by said output step includes a pre-Symbol skewness and kurtosis calculated for each spatial frequency band different from each other in the calculating step at least,
The skin evaluation method according to claim 1.
を更に含む請求項2に記載の肌評価方法。 A multiple regression equation obtained in advance by multiple regression analysis using the skin sample data of the mother group, with a multiple regression equation for the purpose variable smoothness degree of the skin and explanatory variables the skewness and the kurtosis The process of determining the degree of smoothness of the skin,
The skin evaluation method according to claim 2, further comprising.
前記取得された凹凸情報を多重分解して、相互に異なる複数の空間周波数帯域の凹凸情報を抽出する分解手段と、
前記多重分解されたいずれか複数の空間周波数帯域の凹凸情報から、凹凸形状に関する統計値であって空間周波数帯域ごとに異なる種別の複数種の統計値を算出する算出手段と、
前記算出された複数種の統計値を用いて前記肌の評価情報を出力する出力処理手段と、
を備える肌評価装置。 An acquisition method for acquiring information on the unevenness of the skin surface measured from the skin or a replica of the skin,
Decomposition means that multiplexes the acquired unevenness information and extracts unevenness information of a plurality of spatial frequency bands that are different from each other.
A calculation means for calculating a plurality of types of statistical values related to the uneven shape, which are different types for each spatial frequency band, from the unevenness information of any of the plurality of spatial frequency bands that have been multiple-decomposed.
An output processing means for outputting the skin evaluation information using the calculated plurality of types of statistical values, and
A skin evaluation device equipped with.
前記出力処理手段は、前記算出手段により相互に異なる空間周波数帯域ごとに算出された歪度及び尖度を少なくとも用いて前記肌の評価情報を生成する、
請求項4に記載の肌評価装置。 The calculation means calculates at least the skewness and kurtosis of the uneven shape as the plurality of types of statistical values.
The output processing means generates the evaluation information of the skin using at least the skewness and the kurtosis calculated for each spatial frequency band different from each other by the calculation means.
The skin evaluation device according to claim 4.
請求項5に記載の肌評価装置。 The output processing means is a multiple regression equation previously acquired by multiple regression analysis using population skin sample data, with the skewness and kurtosis as explanatory variables and the smoothness of the skin as the objective variable. The degree of smoothness of the skin is determined by using the multiple regression equation.
The skin evaluation device according to claim 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017175149A JP6942012B2 (en) | 2017-09-12 | 2017-09-12 | Skin evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017175149A JP6942012B2 (en) | 2017-09-12 | 2017-09-12 | Skin evaluation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019050855A JP2019050855A (en) | 2019-04-04 |
| JP6942012B2 true JP6942012B2 (en) | 2021-09-29 |
Family
ID=66013187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017175149A Active JP6942012B2 (en) | 2017-09-12 | 2017-09-12 | Skin evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6942012B2 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022118890A1 (en) * | 2020-12-02 | 2022-06-09 | 株式会社資生堂 | Method for generating data on tactile perception of human skin, device for generating data on tactile perception of human skin, method for evaluating tactile perception of human skin, device for evaluating tactile perception of human skin, device for presenting tactile perception of human skin, and method for presenting tactile perception of human skin |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4133248B2 (en) * | 2002-11-18 | 2008-08-13 | 株式会社コーセー | Skin gloss evaluation method |
| JP4808975B2 (en) * | 2005-02-15 | 2011-11-02 | 株式会社 資生堂 | Simulation apparatus and simulation method for evaluating and simulating skin unevenness |
| JP5146667B2 (en) * | 2008-06-13 | 2013-02-20 | 花王株式会社 | Evaluation method of skin smoothness |
| JP6124295B2 (en) * | 2013-05-02 | 2017-05-10 | 日本メナード化粧品株式会社 | Method for evaluating the texture of bare skin or makeup skin |
-
2017
- 2017-09-12 JP JP2017175149A patent/JP6942012B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2019050855A (en) | 2019-04-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ahn et al. | Correlation between a Cutometer® and quantitative evaluation using Moire topography in age‐related skin elasticity | |
| CN102870136B (en) | Age estimation method | |
| JP5146667B2 (en) | Evaluation method of skin smoothness | |
| CN105263399B (en) | Transparent sensation evaluating device, transparent appraisal method | |
| US10743769B2 (en) | Skin evaluation apparatus, skin evaluation method, and skin evaluation program | |
| JP5753055B2 (en) | Skin image analysis apparatus and skin image analysis method | |
| JP7477678B2 (en) | Skin distortion measurement method | |
| JP2017051474A (en) | Skin gloss evaluation device, gloss evaluation method and gloss evaluation program | |
| JP2009082338A (en) | Skin discrimination method using entropy | |
| CN118314284B (en) | Image processing method and system for medical three-dimensional high-resolution display | |
| JP6124295B2 (en) | Method for evaluating the texture of bare skin or makeup skin | |
| JP6008698B2 (en) | Face image analyzing apparatus and face image analyzing method | |
| JP6942012B2 (en) | Skin evaluation method | |
| Bousseau et al. | Gloss perception in painterly and cartoon rendering | |
| JP2014033944A (en) | Method and apparatus of evaluating skin condition | |
| JP6850194B2 (en) | Skin condition evaluation method and skin condition evaluation device | |
| JP7731228B2 (en) | Skin image analysis method | |
| Wilson et al. | The prevalence of lip vermilion morphological traits in a 15‐year‐old population | |
| CN106030659B (en) | Ageing analysis method and ageing analysis device | |
| Kravchenko et al. | High‐fidelity iridal light transport simulations at interactive rates | |
| JP6810619B2 (en) | Skin condition evaluation method | |
| Clarke et al. | Measuring visual discomfort—a novel two-step method for reducing criterion effects when measuring subjective responses | |
| Wei et al. | Objective evaluation of optical illusion skirt based on image texture features | |
| JP2017158921A (en) | Skin condition evaluation method | |
| JP6854702B2 (en) | Rough skin evaluation method and rough skin evaluation device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200605 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20210319 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210323 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20210517 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210701 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210713 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210813 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210831 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210907 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 6942012 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |