JP6124295B2 - Method for evaluating the texture of bare skin or makeup skin - Google Patents

Method for evaluating the texture of bare skin or makeup skin Download PDF

Info

Publication number
JP6124295B2
JP6124295B2 JP2013096813A JP2013096813A JP6124295B2 JP 6124295 B2 JP6124295 B2 JP 6124295B2 JP 2013096813 A JP2013096813 A JP 2013096813A JP 2013096813 A JP2013096813 A JP 2013096813A JP 6124295 B2 JP6124295 B2 JP 6124295B2
Authority
JP
Japan
Prior art keywords
band
index value
skin
texture
image
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
Application number
JP2013096813A
Other languages
Japanese (ja)
Other versions
JP2014217456A (en
Inventor
陽平 石黒
陽平 石黒
彰 木嶋
彰 木嶋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Menard Cosmetic Co Ltd
Original Assignee
Nippon Menard Cosmetic Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nippon Menard Cosmetic Co Ltd filed Critical Nippon Menard Cosmetic Co Ltd
Priority to JP2013096813A priority Critical patent/JP6124295B2/en
Publication of JP2014217456A publication Critical patent/JP2014217456A/en
Application granted granted Critical
Publication of JP6124295B2 publication Critical patent/JP6124295B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)

Description

本発明は、素肌または化粧肌の視覚的質感の評価方法に関する。   The present invention relates to a method for evaluating the visual texture of bare skin or makeup skin.

スキンケア化粧品やメイクアップ化粧品による肌の視覚的質感(以下、「質感」という。)に対する作用を評価するために、素肌または化粧肌を撮像した画像から、光学的な指標値を得ることが行われていた。例えば、素肌または化粧肌を撮像した画像から、分散や歪度を質感の指標値とすることが行われていた(特許文献1)。また、質感は面上における輝度や色の不均一性によって生じ、どのような空間周波数の不均一性が含まれているかによってどのような質感であるかが大きく変化する。そのため、素肌を撮像した画像から鏡面反射成分を抽出し、さらに鏡面反射成分から得る特定の空間周波数成分を基に見かけの粗さを求める手法などが行われていた(特許文献2)。   In order to evaluate the effect of skin care cosmetics and makeup cosmetics on the visual texture of the skin (hereinafter referred to as “texture”), an optical index value is obtained from an image of bare skin or makeup skin. It was. For example, dispersion or distortion is used as a texture index value from an image of bare skin or makeup skin (Patent Document 1). Further, the texture is caused by unevenness of brightness and color on the surface, and the kind of texture varies greatly depending on what kind of spatial frequency non-uniformity is included. For this reason, a method has been used in which a specular reflection component is extracted from an image of bare skin and an apparent roughness is obtained based on a specific spatial frequency component obtained from the specular reflection component (Patent Document 2).

特開2009−131336号公報JP 2009-131336 A 特開2004−166801号公報JP 2004-166801 A

質感と光学的な指標値の関係は、その全てが明らかにされているわけではなく、光学的な指標値を用いた質感の評価方法においても、目的となる質感に対してより適切に評価し得る評価方法の開発が望まれていた。また、しばしば肌における不均一性は「ムラ」として美しさを損ねる要素とされてきたものの、肌は元来均一面ではなく、均一さは逆に人工的な不自然さを生じさせる。そのため、肌が美しさのために有するべき低周波帯域の不均一性や、キメとも呼ばれる肌が美しさのために有するべき高周波帯域の不均一性、美しさを損ねる不均一性である「ムラ」はそれぞれ適切に分類する必要があった。さらに、従来の質感評価方法においては、その容易化、低コスト化、省資源化等が望まれていた。   The relationship between the texture and the optical index value is not completely clarified. Even in the texture evaluation method using the optical index value, the target texture is evaluated more appropriately. Development of an evaluation method that can be obtained has been desired. Moreover, although unevenness in the skin is often regarded as an element that impairs beauty as “unevenness”, the skin is not originally a uniform surface, and the uniformity causes artificial unnaturalness. Therefore, the unevenness of the low frequency band that the skin should have for beauty, the non-uniformity of the high frequency band that the skin should have for beauty, and the nonuniformity that impairs the beauty. Each needed to be classified appropriately. Further, in the conventional texture evaluation method, it has been desired to facilitate, reduce costs, save resources, and the like.

本発明は、かかる事情に鑑みてなされたものであり、その目的は、素肌または化粧肌の視覚的質感を光学的な指標値を用いて定量的に評価することである。さらに詳細には、素肌または化粧肌の不均一性を適切に分類すること、並びに分類された各不均一性の関係から生じる質感を定量的に評価することである。   The present invention has been made in view of such circumstances, and an object thereof is to quantitatively evaluate the visual texture of bare skin or makeup skin using optical index values. More specifically, it is to appropriately classify the non-uniformity of bare skin or makeup skin, and to quantitatively evaluate the texture resulting from the relationship of each classified non-uniformity.

かかる課題を解決するために、本発明者らは、複数の頬を撮像した画像から空間周波数成分を取得し、各空間周波数成分における変化量を算出した結果、特定帯域における変化量が相対的に低いことによって見た目の「透け感」が生じることを見出した。さらにこの知見から、肌の不均一性は、「質感に寄与しない陰影」、「透け感を生じさせる低周波帯域の不均一性」、「透け感を損なう不均一性(ムラ)」、「透け感を生じさせる高周波帯域の不均一性(キメ)」の4帯域に分類できることを見出し、本発明を完成するに至った。なお、以降では、見た目に「透け感」が生じることを、単に「透ける」とも呼ぶ。   In order to solve such a problem, the present inventors acquired spatial frequency components from images obtained by imaging a plurality of cheeks, and calculated the amount of change in each spatial frequency component. It was found that the appearance of “translucency” appears due to the low. Furthermore, from this knowledge, skin non-uniformity is “shading that does not contribute to the texture”, “non-uniformity in the low frequency band that causes translucency”, “nonuniformity that impairs translucency (unevenness)”, “translucency” It has been found that it can be classified into four bands of “non-uniformity (texture) of the high frequency band causing a feeling”, and the present invention has been completed. Hereinafter, the appearance of “translucency” is also simply referred to as “transparency”.

本発明は、上述の課題の少なくとも一部を解決するためになされたものであり、以下の形態として実現することが可能である。   SUMMARY An advantage of some aspects of the invention is to solve at least a part of the problems described above, and the invention can be implemented as the following forms.

(1)本発明の一形態によれば、素肌または化粧肌の質感を評価する方法が提供される。この方法は、(a)素肌または化粧肌の撮像画像を得る工程と、(b)前記撮像画像の各空間周波数帯域における成分を得る工程と、(c)前記各空間周波数帯域における成分を少なくとも4帯域以上の複数帯域に分類する工程と、(d)前記複数帯域毎の変化量の指標値を得る工程と、(e)前記複数帯域の変化量の関係性から生じる質感の指標値を得る工程と、(f)前記質感の指標値を基に、素肌または化粧肌の質感を評価する工程と、を備える。この形態の方法によれば、素肌または化粧肌の撮像画像を空間周波数成分毎に分類し、指標値を得るので、各帯域の指標値から生じる質感を客観的に評価することができる。換言すると、スキンケア化粧品またはメイクアップ化粧品を使用した際に起こる各帯域の指標値の変化から生じる質感の変化を、客観的に評価することができる。なお、「前記複数帯域の変化量の関係性から生じる質感」としては、例えば、第1〜4帯域の4帯域に分類した際の、第2〜4帯域の変化量の関係から、第3帯域の変化量が相対的に低いことによって生じる透け感があげられる。また、前記工程(c)の分類は、4帯域に分類することがより好ましいが、本発明の趣旨を損なわない範囲で、更に細分化することもできる。   (1) According to one aspect of the present invention, a method for evaluating the texture of bare skin or makeup skin is provided. The method includes (a) obtaining a captured image of bare skin or makeup skin, (b) obtaining a component in each spatial frequency band of the captured image, and (c) at least four components in each spatial frequency band. A step of classifying into a plurality of bands equal to or higher than a band; (d) a step of obtaining an index value of the change amount for each of the plurality of bands; And (f) evaluating the texture of bare skin or makeup skin based on the texture index value. According to the method of this aspect, the captured image of bare skin or makeup skin is classified for each spatial frequency component and an index value is obtained, so that the texture generated from the index value of each band can be objectively evaluated. In other words, it is possible to objectively evaluate the change in texture caused by the change in the index value of each band that occurs when using skin care cosmetics or makeup cosmetics. The “texture resulting from the relationship between the change amounts of the plurality of bands” is, for example, from the relationship between the change amounts of the second to fourth bands when classified into the four bands of the first to fourth bands. There is a sense of see-through caused by a relatively low amount of change. The classification of the step (c) is more preferably classified into four bands, but can be further subdivided within a range not impairing the gist of the present invention.

(2)上記形態の方法において、前記工程(c)の分類は、各空間周波数成分を抜いた画像を用いた心理物理実験により、抜くことで透けて見える帯域(以下、「当該帯域」という。)を抽出し、第1帯域:1cycle/画像幅、第2帯域:2cycle/画像幅〜当該帯域未満、第3帯域:当該帯域、第4帯域:当該帯域よりも高い帯域、の4帯域として分類してもよい。この形態の方法によれば、人の質感知覚機構に準じて分類するので、適切に分類することができる。   (2) In the method of the above aspect, the classification of the step (c) is a band that can be seen through by a psychophysical experiment using an image from which each spatial frequency component is removed (hereinafter referred to as “the band”). ) Are extracted and classified as four bands: first band: 1 cycle / image width, second band: 2 cycle / image width to less than the band, third band: the band, and fourth band: a band higher than the band. May be. According to this form of the method, classification is performed according to the human texture perception mechanism, and therefore, classification can be performed appropriately.

(3)上記形態の方法において、前記工程(c)の分類は、顔面頬部の幅を基準として、第1帯域:1cycle/頬幅、第2帯域:2〜8cycle/頬幅、第3帯域:16〜32cycle/頬幅、第4帯域:64cycle〜/頬幅の4帯域として分類してもよい。この形態の方法によれば、短時間で解析結果を得ることができる。   (3) In the method of the above aspect, the classification in the step (c) is based on the width of the facial cheek part as follows: first band: 1 cycle / cheek width, second band: 2-8 cycles / cheek width, third band : 16 to 32 cycles / cheek width, fourth band: 64 cycles ~ / cheek width may be classified as four bands. According to this method, the analysis result can be obtained in a short time.

(4)上記形態の方法において、前記工程(e)は、第1〜4帯域の4帯域に分類した際の、第2〜4帯域の変化量の指標値をそれぞれSDとして、式1によって得る値を指標値としてもよい。式1では、第3帯域の変化量の指標値が第2または第4帯域の変化量の指標値よりも多い時にマイナスの値を取り、相対的に少ないほど高い値を取る。この形態の方法によれば、第3帯域の変化量が相対的に低いことによって生じる透け感を短時間で評価することができる。
(式1)指標値=(SD―SD)×(SD―SD
但し、SD>SDかつSD>SDの場合は評価不能とする
(4) In the method of the above aspect, in the step (e), the index values of the change amounts of the second to fourth bands when classified into the four bands of the first to fourth bands are expressed as SD 2 to 4 , respectively. A value obtained by 1 may be used as an index value. In Expression 1, a negative value is taken when the index value of the change amount of the third band is larger than the index value of the change amount of the second or fourth band, and a higher value is taken as it is relatively smaller. According to the method of this embodiment, it is possible to evaluate a sense of sheer caused by a relatively low change amount of the third band in a short time.
(Expression 1) Index value = (SD 2 −SD 3 ) × (SD 4 −SD 3 )
However, if SD 3 > SD 2 and SD 3 > SD 4 , the evaluation is impossible.

(5)上記形態の方法において、前記工程(e)は、、第1〜4帯域の4帯域に分類した際の、第2帯域の変化量の指標値を、肌が有するべき不均一性の強さの指標値としてもよい。この形態の方法によれば、人の質感知覚機構に準じて分類された肌が有するべき不均一性の強さを短時間で評価することができる。   (5) In the method of the above aspect, in the step (e), the index value of the change amount of the second band when the four bands of the first to fourth bands are classified, It may be an index value of strength. According to this form of the method, it is possible to evaluate in a short time the strength of non-uniformity that the skin classified according to the human texture perception mechanism should have.

(6)上記形態の方法において、前記工程(e)は、、第1〜4帯域の4帯域に分類した際の、第3帯域の変化量の指標値を、ムラの強さの指標値としてもよい。この形態の方法によれば、人の質感知覚機構に準じて分類されたムラの強さを短時間で評価することができる。   (6) In the method of the above aspect, in the step (e), the index value of the change amount of the third band when classified into the four bands of the first to fourth bands is used as the index value of the intensity of unevenness. Also good. According to this form of the method, the intensity of unevenness classified according to the human texture perception mechanism can be evaluated in a short time.

(7)上記形態の方法において、前記工程(e)は、、第1〜4帯域の4帯域に分類した際の、第4帯域の変化量の指標値を、キメの強さの指標値としてもよい。この形態の方法によれば、人の質感知覚機構に準じて分類されたキメの強さを短時間で評価することができる。   (7) In the method of the above aspect, in the step (e), the index value of the change amount of the fourth band when the four bands of the first to fourth bands are classified is used as the index value of the strength of the texture. Also good. According to this form of the method, the strength of the texture classified according to the human texture perception mechanism can be evaluated in a short time.

(8)上記形態の方法において、前記工程(b)は、ウェーブレット変換を用いてもよい。この形態の方法によれば、容易に各空間周波数成分を得ることができるので、短時間で解析結果を得ることができる。   (8) In the method of the above aspect, the step (b) may use wavelet transform. According to the method of this embodiment, each spatial frequency component can be easily obtained, so that an analysis result can be obtained in a short time.

(9)上記形態の方法において、前記工程(d)の指標値は、画素値から統計的に得られる値としてもよい。この形態の方法によれば、統計的手法により得られる値を指標値とするので、客観的な解析結果を得ることができる。   (9) In the method of the above aspect, the index value in the step (d) may be a value statistically obtained from the pixel value. According to the method of this embodiment, since the value obtained by the statistical method is used as the index value, an objective analysis result can be obtained.

本発明は、種々の態様で実現することも可能である。例えば、上述した評価方法の少なくとも一部の工程を実現する装置や、コンピュータープログラムや、そのコンピュータープログラムを記録した一時的でない記録媒体などの形態で実現することができる。   The present invention can be realized in various modes. For example, the present invention can be realized in the form of an apparatus that implements at least a part of the above-described evaluation method, a computer program, or a non-temporary recording medium that records the computer program.

本実施形態における質感の評価方法の手順を示すフローチャートである。It is a flowchart which shows the procedure of the evaluation method of the texture in this embodiment. ステップS115において、空間周波数成分を4帯域に分類する方法の手順を示すフローチャートである。It is a flowchart which shows the procedure of the method which classify | categorizes a spatial frequency component into 4 bands in step S115. ステップS210の処理結果を模式的に示す説明図である。It is explanatory drawing which shows the process result of step S210 typically. ステップS125における質感の指標値の一例を示す説明図である。It is explanatory drawing which shows an example of the parameter | index value of the texture in step S125. 第1実施例の撮像画像を示す説明図である。It is explanatory drawing which shows the captured image of 1st Example.

A.実施形態
A1.素肌または化粧肌の質感評価方法
図1は、本実施形態における、素肌または化粧肌の質感評価方法の手順を示すフローチャートである。本実施形態において、化粧肌とは、化粧下地や、ファンデーションや、おしろいや、チークなど、皮膚に塗布し得る任意の化粧料を、いずれかひとつまたは組み合わせて塗布した肌を意味する。
A. Embodiment A1. Bare or Makeup Skin Texture Evaluation Method FIG. 1 is a flowchart showing a procedure of a bare skin or makeup skin texture evaluation method in the present embodiment. In the present embodiment, the makeup skin means skin applied with any one or a combination of any cosmetics that can be applied to the skin, such as a makeup base, foundation, funnier, and teak.

図1に示すように、まず素肌または化粧肌の撮像画像を得る(ステップS105)。撮像装置としては、デジタルスチルカメラやデジタルビデオカメラなど、デジタル画像データを得ることができる任意のデバイスを採用することができる。但し、頬の質感に寄与する成分(例えば、キメ)を解析するのに十分な解像度と、頬全体の輝度や色の空間的な変化を解析するのに十分な撮像範囲の画像を得ることが可能なデバイスが好ましい。また撮像画像は、頬幅に合わせてトリミングすることが好ましい。なお、以降では、画素値として輝度を用いて説明するが、輝度に限らず、色相や彩度などの画素値を採用することができる。   As shown in FIG. 1, first, a captured image of bare skin or makeup skin is obtained (step S105). As the imaging apparatus, any device that can obtain digital image data, such as a digital still camera or a digital video camera, can be employed. However, it is possible to obtain an image having a resolution sufficient for analyzing a component (for example, texture) that contributes to the texture of the cheek and an imaging range sufficient for analyzing a spatial change in luminance and color of the entire cheek. Possible devices are preferred. The captured image is preferably trimmed to fit the cheek width. In the following description, luminance is used as the pixel value. However, not only the luminance but also pixel values such as hue and saturation can be adopted.

ステップS105で得られた撮像画像を空間周波数成分に分解する(ステップS110)。具体的な処理手順としては、先ず、撮像画像から各画素の輝度値データを取得し、ウェーブレット変換等によって各空間周波数成分を得る。撮像画像を頬幅に合わせてトリミングした画像を基にすることにより、頬幅を基準として2cycle/頬幅(n=0,1,2...)の空間周波数成分を得ることができる。ウェーブレット変換を用いた場合には最も単純なアルゴリズムであるHaarの基底を採用することができ、プログラムを作成し、パーソナルコンピューター等の演算装置に実装することによって実行することができる。 The captured image obtained in step S105 is decomposed into spatial frequency components (step S110). As a specific processing procedure, first, luminance value data of each pixel is acquired from a captured image, and each spatial frequency component is obtained by wavelet transform or the like. By using an image obtained by trimming the captured image according to the cheek width, a spatial frequency component of 2 n cycles / cheek width (n = 0, 1, 2,...) Can be obtained with the cheek width as a reference. . When the wavelet transform is used, Haar's basis, which is the simplest algorithm, can be adopted, and it can be executed by creating a program and mounting it on a computing device such as a personal computer.

得られる空間周波数成分を第1〜4帯域の4帯域に分類する(ステップS115)。本実施形態では、分類方法として、顔面の頬幅を基準として、第一帯域:1cycle/頬幅、第2帯域:2〜8cycle/頬幅、第3帯域:16〜32cycle/頬幅、第4帯域:64cycle〜/頬幅として分類する方法を採用する。   The obtained spatial frequency components are classified into four bands of the first to fourth bands (step S115). In the present embodiment, as a classification method, with reference to the cheek width of the face, the first band: 1 cycle / cheek width, the second band: 2-8 cycles / cheek width, the third band: 16-32 cycles / cheek width, the fourth Band: adopts a method of classifying as 64 cycle ~ / cheek width.

ステップS115で得られた各帯域成分を統合し、各帯域の変化量の指標値を取得する(ステップS120)。前記統合は、パーソナルコンピューターなどの演算装置によって加算処理を行うことで実行することができる。また前記指標値には、分散や標準偏差など、変化量を表す指標値を採用することができる。   Each band component obtained in step S115 is integrated, and an index value of the change amount of each band is acquired (step S120). The integration can be executed by performing an addition process using an arithmetic device such as a personal computer. In addition, as the index value, an index value representing the amount of change such as variance or standard deviation can be adopted.

ステップS120で得られた各帯域成分の変化量の指標値から、それらの関係によって生じる質感の指標値を取得する(ステップS125)。この工程は、前記指標値を基に、前述の式1を透け感の指標値とする算出方法を採用することができる。   From the index value of the change amount of each band component obtained in step S120, the texture index value generated by the relationship is acquired (step S125). In this step, a calculation method can be employed in which the above-described Equation 1 is used as the index value of the sense of transparency based on the index value.

前記質感の指標値としては、上述した透け感の指標値の他、各帯域の変化量の指標値をそのまま質感の指標値として採用することもできる。具体的には、第2帯域の変化量の指標値を肌が有するべき不均一性の強さの指標値とする方法、第3帯域の変化量の指標値を肌のムラの指標値とする方法、第4帯域の変化量の指標値を肌のキメの指標値とする方法を採用することができる。なお、第1帯域は、陰影の情報を含み、質感への寄与は低い。   As the texture index value, in addition to the above-described index value of translucency, the index value of the change amount of each band can be directly adopted as the texture index value. Specifically, the index value of the amount of change in the second band is used as an index value of the intensity of non-uniformity that the skin should have, and the index value of the amount of change in the third band is used as an index value of skin unevenness. It is possible to adopt a method that uses an index value for the amount of change in the fourth band as an index value for skin texture. Note that the first band includes shadow information and has a low contribution to the texture.

指標値が算出されると、得られた指標値に基づき、撮像した素肌または化粧肌の質感を評価する(ステップS130)。具体的には、上述した質感の指標値のうち、前記透け感の指標値を採用した場合、指標値がより高いほど、透け感がより高いと評価する。   When the index value is calculated, the texture of the captured skin or makeup skin is evaluated based on the obtained index value (step S130). Specifically, among the above-described index values of the texture, when the index value of the sense of sheer is adopted, it is evaluated that the sense of sheer is higher as the index value is higher.

以上のように素肌または化粧肌の質感を評価することができるので、化粧品を選定する際に、質感制御効果を考慮することができる。例えば、肌の質感改善効果を有する基礎化粧品を選定する場合に、連用前後に撮像した画像を用いて、上述した処理を行うことにより連用前後の肌の質感を把握して、質感改善効果の高さに基づき、いずれの基礎化粧品を選定するかを決定することができる。また別の例としては、メイクアップ化粧品を選定する場合に、複数のメイクアップ化粧品を製作し、これら試料を対象として上述した処理を行うことにより各試料によってつくられる肌の質感を把握して、かかる質感の制御効果に基づき、いずれのメイクアップ化粧品を選定するかを決定することができる。なお、各帯域成分の変化量は、メイクアップ化粧品の形態、メイクアップ化粧品に配合する原材料、塗布具の形態や材質、複数のメイクアップ化粧品の組み合わせ方法など種々の要素によって制御することができる。したがって、質感の制御効果の評価結果に基づき、前述の要素を調節することもできる。   Since the texture of bare skin or makeup skin can be evaluated as described above, the texture control effect can be taken into account when selecting cosmetics. For example, when selecting a basic cosmetic product that has an effect of improving the texture of the skin, it is possible to grasp the texture of the skin before and after continuous use by performing the above-described processing using the images taken before and after the continuous use, resulting in a high texture improvement effect. Based on this, it is possible to determine which basic cosmetic product is selected. As another example, when selecting a makeup cosmetic, a plurality of makeup cosmetics are manufactured, and by performing the above-described processing for these samples, the texture of the skin made by each sample is grasped, Based on the effect of controlling the texture, it can be determined which makeup cosmetic is selected. Note that the amount of change of each band component can be controlled by various factors such as the makeup cosmetic form, the raw materials to be blended in the makeup cosmetic, the form and material of the applicator, and the combination method of a plurality of makeup cosmetics. Therefore, the above-mentioned elements can be adjusted based on the evaluation result of the texture control effect.

A2.帯域の分類方法
図2は、上述したステップS115において、空間周波数成分を4帯域に分類する処理の手順を示すフローチャートである。上述したように、本実施形態のステップS115において、各空間周波数成分の分類は、顔面の頬幅を基準として第1帯域:1cycle/頬幅、第2帯域:2〜8cycle/頬幅、第3帯域:16〜32cycle/頬幅、第4帯域:64cycle〜/頬幅として分類する。かかる分類は、以下のようにして決定することができる。
A2. Band Classification Method FIG. 2 is a flowchart showing a processing procedure for classifying spatial frequency components into four bands in step S115 described above. As described above, in step S115 of the present embodiment, the classification of each spatial frequency component is based on the cheek width of the face as a reference: first band: 1 cycle / cheek width, second band: 2-8 cycles / cheek width, third The band is classified as 16 to 32 cycles / cheek width, and the fourth band: 64 cycles to / cheek width. Such classification can be determined as follows.

先ず、空間周波数成分ごとの変化量が一定のノイズ画像を作成する(ステップS205)。ノイズ画像の作成は、例えば、グレースケールのホワイトノイズ画像からウェーブレット変換により各空間周波数成分に分解し、各空間周波数成分の変化量として、0.0〜1.0のスケールにおける標準偏差を0.1とする規格化を行い、それらの空間周波数成分を再度統合することで作成する。なお、かかる画像の大きさとしては、例えば、縦及び横がいずれも256ピクセルの画像を採用することができる。また、このノイズ画像は、位相を変えて複数種類作成することが望ましい。   First, a noise image having a constant change amount for each spatial frequency component is created (step S205). The noise image is created by, for example, decomposing a grayscale white noise image into each spatial frequency component by wavelet transform, and setting the standard deviation on a scale of 0.0 to 1.0 as 0. It is created by standardizing to 1 and integrating those spatial frequency components again. As the size of the image, for example, an image having 256 pixels in both vertical and horizontal directions can be adopted. Further, it is desirable to create a plurality of types of noise images with different phases.

ステップS205で得たノイズ画像を基に、特定帯域成分だけを抜いた画像群を作成する(ステップS210)。例えば、1cycle/画像幅の成分だけを抜いた画像や、32cycle/画像幅の成分だけを抜いた画像などを作成する。特定帯域成分だけを抜いた画像の作成方法としては、例えば、ノイズ画像をウェーブレット変換によって各空間周波数成分に分解し、特定帯域成分だけを抜いて統合することで作成することができる。   Based on the noise image obtained in step S205, an image group from which only a specific band component is removed is created (step S210). For example, an image from which only a component of 1 cycle / image width is removed or an image from which only a component of 32 cycles / image width is removed is created. As a method of creating an image from which only a specific band component is removed, for example, it can be created by decomposing a noise image into each spatial frequency component by wavelet transform and extracting and integrating only the specific band component.

図3は、ステップS210において作成される画像を示す説明図である。各画像は、空間周波数成分ごとの変化量が一定のノイズから、各画像の下に示された空間周波数成分を抜いて作成された画像である。   FIG. 3 is an explanatory diagram showing the image created in step S210. Each image is an image created by removing the spatial frequency component shown below each image from noise with a constant change amount for each spatial frequency component.

特定帯域成分を抜いた画像群が得られると、被験者により、当該帯域が心理物理実験により特定される(ステップS215)。図3に示すように、特定空間周波数成分を抜いた画像のうち幾つかのものは、単一の層ではなく手前と奥に異なる複数の層が存在し、手前の層が透けているように知覚される。   When the image group from which the specific band component is removed is obtained, the band is specified by the subject through a psychophysical experiment (step S215). As shown in FIG. 3, in some images from which a specific spatial frequency component is removed, a plurality of different layers exist in the foreground and the back rather than a single layer so that the foreground layer is transparent. Perceived.

図4に示す例では、4人の被験者により、暗室に10分間順応した後、72dpiの提示装置から60cmの距離でランダムな順序で画像が提示され、透け感について「感じない(スコア0.0点)」、「やや感じる(スコア0.5点)」、「感じる(スコア1.0点)」の3つの回答から選択する心理物理実験を実施した。なお、回答は、同じ画像に対してノイズの位相を変えて6回の回答を行った。なお、ここでいう「透ける」とは、より厳密には、異なる別々の層が透過していることで重なり合っているように知覚されることと定義され、被験者にもそのように説明された。被験者はいずれも裸眼又は矯正された状態であり、正常の視力を有していた。   In the example shown in FIG. 4, after four subjects have adapted to a dark room for 10 minutes, images are presented in a random order at a distance of 60 cm from a 72 dpi presentation device. A psychophysical experiment was performed by selecting from the three answers of “point)”, “slightly feel (score 0.5 points)”, and “feel (score 1.0 points)”. The answers were made six times with the noise phase changed for the same image. It should be noted that “transparent” here is defined more strictly as being perceived as overlapping because different layers are transmitting, and was explained to the subject as well. All subjects were either naked or corrected and had normal visual acuity.

図4に示すように、16cycle/画像幅および32cycle/画像幅の帯域成分を抜いた画像において、透け感が増幅される。これは、16cycle/画像幅および32cycle/画像幅の成分を抜くことで、透け感が起こることを意味する。したがって、図4の例によると、ステップS210における抜くことで透ける帯域として、16〜32cycle/画像幅が決定され得る。なお、図4の例による2〜8cycle/画像幅の帯域は逆に抜くことで透け感が低下する。したがって、2〜8cycle/画像幅の帯域は透け感を生じさせるために必要な帯域であり、16〜32cycle/画像幅の帯域とは分類されるべきであることがわかる。   As shown in FIG. 4, the sense of see-through is amplified in an image from which band components of 16 cycles / image width and 32 cycles / image width are removed. This means that a sense of see-through occurs by removing components of 16 cycles / image width and 32 cycles / image width. Therefore, according to the example of FIG. 4, 16 to 32 cycles / image width can be determined as a transparent band by removing in step S210. It should be noted that the sense of transparency is reduced by removing the band of 2 to 8 cycles / image width according to the example of FIG. Therefore, it can be seen that the band of 2 to 8 cycles / image width is a band necessary for producing a sense of transparency, and the band of 16 to 32 cycles / image width should be classified.

抜くことで透ける帯域が特定されると、各空間周波数成分の分類を決定する(ステップS220)。具体的には、第1帯域:1cycle/画像幅、第2帯域:2cycle/画像幅〜当該帯域未満、第3帯域:当該帯域、第4帯域:当該帯域よりも高い帯域として分類する。したがって、図4の例によると、第1帯域:1cycle/画像幅、第2帯域:2〜8cycle/画像幅、第3帯域:16〜32cycle/画像幅、第4帯域:64cycle〜/画像幅として分類され得る。なお、上記分類を基に、ステップS210で得られる画像群から前述の式1により得られる透け感の指標値を算出すると、16cycle/画像幅および32cycle/画像幅の帯域成分を抜いた画像において高い値を取り、その他の画像においては0に近似した値かマイナスの値を取る。   When the transparent band is specified by extracting, the classification of each spatial frequency component is determined (step S220). Specifically, the first band: 1 cycle / image width, the second band: 2 cycle / image width to less than the band, the third band: the band, and the fourth band: a band higher than the band. Therefore, according to the example of FIG. 4, the first band: 1 cycle / image width, the second band: 2-8 cycles / image width, the third band: 16-32 cycles / image width, and the fourth band: 64 cycles / image width. Can be classified. Based on the above classification, when the index value of the translucency obtained by the above-described equation 1 is calculated from the image group obtained in step S210, it is high in an image from which band components of 16 cycles / image width and 32 cycles / image width are removed. It takes a value and takes a value close to 0 or a negative value in other images.

以上説明した実施形態における質感の評価方法に依れば、心理物理実験によって空間周波数成分の分類が決定されるので、人の質感知覚機構に準じた適切な分類ができる。したがって、適切に分類された空間周波数成分を基に光学的指標値を得ることにより、「透け感を生じさせる低周波帯域の不均一性」、「透け感を損なう不均一性(ムラ)」、「透け感を生じさせる高周波帯域の不均一性(キメ)」を人の質感知覚機構に準じて客観的に分類し、評価することができる。   According to the texture evaluation method in the embodiment described above, the classification of spatial frequency components is determined by psychophysical experiment, and therefore, appropriate classification according to the human texture perception mechanism can be performed. Therefore, by obtaining an optical index value based on appropriately classified spatial frequency components, "non-uniformity in the low frequency band that causes a sense of sheer", "non-uniformity that impairs the sheer feeling (unevenness)", It is possible to objectively classify and evaluate “non-uniformity (texture) in the high-frequency band causing a sense of see-through” according to the human texture perception mechanism.

さらに、従来明らかにされていなかった、複数の空間周波数成分の関係から生じる質感を定量的に評価することができる。具体的には、本発明の実施形態における質感の評価方法に依れば、抜くことで透ける帯域が心理物理実験によって特定される。換言すれば、他の帯域に比べ相対的に変化量が無いことで透ける第3帯域と、透け感を生じさせるために変化量が必要な帯域である第2帯域が特定されるため、光学的指標値である第2〜4帯域の変化量の関係性から、第3帯域の変化量の相対的な低さによって生じる透け感を定量的に評価することができる。   Furthermore, it is possible to quantitatively evaluate the texture that has not been clarified heretofore and that arises from the relationship between a plurality of spatial frequency components. Specifically, according to the texture evaluation method in the embodiment of the present invention, a band that can be seen through is specified by a psychophysical experiment. In other words, the third band that is transparent because there is no change amount relative to the other bands and the second band that is a band that requires a change amount in order to generate a sense of see-through are specified. From the relationship between the change amounts of the second to fourth bands, which are index values, it is possible to quantitatively evaluate the sense of transparency caused by the relative low change amount of the third band.

また、頬を撮像した画像を基に、ウェーブレット変換を行って各空間周波数成分を取得し、それを基に肌の質感の光学的指標値を得るので、短期間で解析結果を得ることができる。   Also, wavelet transform is performed on the basis of the image of the cheek image to acquire each spatial frequency component, and the optical index value of the skin texture is obtained on the basis thereof, so that the analysis result can be obtained in a short period of time. .

B.実施例
実施例では、前述のステップS105において、皮膚角質層の透過性を高める美容液による1カ月の連用試験を行い、連用試験開始前、1カ月の連用試験後でそれぞれ頬をデジタルカメラによってストロボ光下で撮像し、頬幅に合わせてトリミングした。
B. Example In the example, in the above-described step S105, a one-month continuous test with a cosmetic liquid that enhances the permeability of the skin stratum corneum is performed, and each cheek is strobed by a digital camera before the continuous test is started and after the one-month continuous test. Images were taken under light and trimmed to fit the cheek width.

図5は、実施例の連用試験前後の画像を用いた説明図である。図5(a)は、連用試験開始前の画像を示す。図5(b)は、1カ月の連用試験後の画像を示す。なお、この画像は9名の連用試験被験者を連用前後で撮像した画像のうち、最も皮膚の透け感を高める効果があらわれた画像として複数の専門家によって目視により選出されたものである。   FIG. 5 is an explanatory diagram using images before and after the continuous test of the example. FIG. 5A shows an image before the continuous test is started. FIG. 5B shows an image after a one-month continuous test. In addition, this image was visually selected by a plurality of experts as an image in which the effect of enhancing the translucency of the skin was most apparent among images obtained by photographing nine consecutive test subjects before and after continuous use.

撮像した2枚の画像をそれぞれ、ウェーブレット変換によって空間周波数成分を取得し(ステップS110)、各空間周波数成分を第1帯域:1cycle/頬幅、第2帯域:2〜8cycle/頬幅、第3帯域:16〜32cycle/頬幅、第4帯域:64cycle〜/頬幅として4帯域に分類した(ステップS115)。   Spatial frequency components are acquired for each of the two captured images by wavelet transform (step S110), and each spatial frequency component is determined for the first band: 1 cycle / cheek width, the second band: 2-8 cycles / cheek width, and third. Bands were classified into four bands as 16 to 32 cycles / cheek width and fourth band: 64 cycles to / cheek width (step S115).

前述の2枚の画像において取得した各帯域の輝度(L*値)標準偏差をそれぞれ各帯域の変化量の指標値として算出し(ステップS120)、前述の式1により得られる値を透け感の指標値として算出した(ステップS125)。その結果、連用試験開始前の値は−0.6×10−3であり、1カ月の連用試験後は1.5×10−2であった。したがって、ステップS130では、1カ月の連用試験後の方が透け感が高いと評価された。この画像は最も透け感を高める効果が表れた画像として選出された画像であることからも、上述した評価方法により得られた結果は、官能評価に沿ったものであった。 The luminance (L * value) standard deviation of each band acquired in the two images described above is calculated as an index value of the change amount of each band (step S120), and the value obtained by the above-described equation 1 is used as a transparent feeling. The index value was calculated (step S125). As a result, the value before the start of the continuous test was −0.6 × 10 −3 and 1.5 × 10 −2 after the continuous test for one month. Therefore, in step S130, it was evaluated that the sense of sheer was higher after one month of continuous testing. Since this image is an image selected as an image that has the most effect of enhancing the sense of transparency, the results obtained by the above-described evaluation method are in accordance with sensory evaluation.

C.変形例
C1.変形例1
上述した実施形態及び実施例では、ステップS110において、ウェーブレット変換を行うことにより各空間周波数成分を取得していたが、ウェーブレット変換に代えて、フーリエ変換を用いて各空間周波数成分を得ることもできる。すなわち、一般には、撮像画像を空間周波数解析可能な任意の手法を採用することができる。
C. Modification C1. Modification 1
In the above-described embodiments and examples, each spatial frequency component is acquired by performing wavelet transform in step S110. However, instead of wavelet transform, each spatial frequency component can also be obtained using Fourier transform. . In other words, in general, any method capable of spatial frequency analysis of a captured image can be employed.

C2.変形例2
上述した実施形態及び実施例では、ステップS120において、指標値を得るために4帯域に分類された各空間周波数成分を統合するとしたが、これを省略し、各帯域の指標値を、各帯域に該当する空間周波数成分の指標値から平均値を算出することによって得てもよい。
C2. Modification 2
In the above-described embodiments and examples, in step S120, the spatial frequency components classified into the four bands are integrated in order to obtain the index value. However, this is omitted, and the index value of each band is assigned to each band. You may obtain by calculating an average value from the index value of the applicable spatial frequency component.

C3.変形例3
上述した実施形態及び実施例において、肌の撮像(ステップS105)を、質感の評価とは別に既に実行されている場合には、質感の評価手順におけるステップS105を省略することもできる。この場合、ステップS105に代えて、「肌画像を取得する」工程を行う。肌画像の取得とは、例えば、ステップS110〜125を実行する演算装置に対して、肌の画像データを入力する工程を意味する。
C3. Modification 3
In the embodiments and examples described above, when the imaging of the skin (step S105) has already been performed separately from the texture evaluation, step S105 in the texture evaluation procedure can be omitted. In this case, a process of “acquiring a skin image” is performed instead of step S105. The acquisition of the skin image means, for example, a process of inputting skin image data to the arithmetic device that executes Steps S110 to S125.

本発明は、上述の実施形態や実施例、変形例に限られるものではなく、その趣旨を逸脱しない範囲において種々の構成で実現することができる。例えば、発明の概要の欄に記載した各形態中の技術的特徴に対応する実施形態、変形例中の技術的特徴は、上述の課題の一部又は全部を解決するために、あるいは、上述の効果の一部又は全部を達成するために、適宜、差し替えや、組み合わせを行うことが可能である。また、その技術的特徴が本明細書中に必要なものとして説明されていなければ、適宜、削除することができる。   The present invention is not limited to the above-described embodiments, examples, and modifications, and can be realized with various configurations without departing from the spirit thereof. For example, the technical features in the embodiments and the modifications corresponding to the technical features in each embodiment described in the summary section of the invention are to solve some or all of the above-described problems, or In order to achieve part or all of the effects, replacement or combination can be performed as appropriate. Further, if the technical feature is not described as necessary in the present specification, it can be appropriately deleted.

Claims (8)

素肌または化粧肌の視覚的質感の評価方法であって、
(a)素肌または化粧肌の撮像画像を得る工程と、
(b)前記撮像画像の各空間周波数帯域における成分を得る工程と、
(c)前記各空間周波数帯域における成分を、各空間周波数成分を抜いた画像を用いた心理物理実験により、抜くことで透けて見える帯域(当該帯域)を抽出し、第1帯域:1cycle/画像幅、第2帯域:2cycle/画像幅〜当該帯域未満、第3帯域:当該帯域、第4帯域:当該帯域よりも高い帯域、の4帯域に分類する工程と、
(d)前記4帯域の変化量の指標値を得る工程と、
(e)前記4帯域の変化量の指標値から生じる視覚的質感の指標値を得る工程と、
(f)前記視覚的質感の指標値を基に、素肌または化粧肌の視覚的質感を評価する工程と、
を備える、方法。
A method for evaluating the visual texture of bare skin or makeup skin,
(A) obtaining a captured image of bare skin or makeup skin;
(B) obtaining a component in each spatial frequency band of the captured image;
(C) Extracting a band (the band) that can be seen through by extracting the component in each spatial frequency band by a psychophysical experiment using an image from which each spatial frequency component is removed, and the first band: 1 cycle / image Classifying into four bands of width, second band: 2 cycle / image width to less than the band, third band: the band, and fourth band: a band higher than the band ;
(D) obtaining an index value of the change amount of the four bands ;
(E) obtaining a visual texture index value generated from the index value of the change amount of the four bands ;
(F) evaluating the visual texture of bare skin or makeup skin based on the visual texture index value;
A method comprising:
請求項1記載の方法において、
前記工程(c)の分類は、顔面頬部の幅を基準として、第1帯域:1cycle/頬幅、第2帯域:2〜8cycle/頬幅、第3帯域:16〜32cycle/頬幅、第4帯域:64cycle〜/頬幅の4帯域として分類する、方法。
The method of claim 1 , wherein
The classification of the step (c) is based on the width of the facial cheek part, the first band: 1 cycle / cheek width, the second band: 2-8 cycles / cheek width, the third band: 16-32 cycles / cheek width, 4 bands: A method of classifying as 4 bands of 64 cycles ~ / cheek width.
請求項1又は2記載の方法において、
前記工程(e)は、第2〜4帯域の変化量の指標値をそれぞれSD2〜4として、式1によって得る値を透け感の指標値とする、方法。
(式1)指標値=(SD―SD)×(SD―SD
但し、SD>SDかつSD>SDの場合は評価不能とする
The method according to claim 1 or 2 , wherein
The step (e) is a method in which the index value of the amount of change in the second to fourth bands is SD 2 to 4 , respectively, and the value obtained by Equation 1 is used as the index value of the sense of transparency.
(Expression 1) Index value = (SD 2 −SD 3 ) × (SD 4 −SD 3 )
However, if SD 3 > SD 2 and SD 3 > SD 4 , the evaluation is impossible.
請求項1又は2記載の方法において、前記工程(e)は、第2帯域の変化量の指標値を、肌が有するべき不均一性の強さの指標値とする、方法。 The method according to claim 1 or 2 , wherein the step (e) uses the index value of the amount of change in the second band as an index value of the strength of non-uniformity that the skin should have. 請求項1又は2記載の方法において、前記工程(e)は、第3帯域の変化量の指標値を、ムラの強さの指標値とする、方法。 3. The method according to claim 1 , wherein the step (e) uses an index value of the amount of change in the third band as an index value of unevenness intensity. 請求項1又は2記載の方法において、前記工程(e)は、第4帯域の変化量の指標値を、キメの強さの指標値とする、方法。 3. The method according to claim 1 , wherein the step (e) uses an index value of the amount of change in the fourth band as an index value of texture strength. 請求項1〜6のいずれか記載の方法において、
前記工程(b)は、ウェーブレット変換である、方法。
In the method in any one of Claims 1-6 ,
The method wherein step (b) is wavelet transform.
請求項1〜7のいずれか記載の方法において、
前記工程(d)の変化量の指標値は、画素値から統計的に得られる値である、方法。
In the method in any one of Claims 1-7 ,
The method, wherein the index value of the change amount in the step (d) is a value statistically obtained from the pixel value.
JP2013096813A 2013-05-02 2013-05-02 Method for evaluating the texture of bare skin or makeup skin Active JP6124295B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013096813A JP6124295B2 (en) 2013-05-02 2013-05-02 Method for evaluating the texture of bare skin or makeup skin

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013096813A JP6124295B2 (en) 2013-05-02 2013-05-02 Method for evaluating the texture of bare skin or makeup skin

Publications (2)

Publication Number Publication Date
JP2014217456A JP2014217456A (en) 2014-11-20
JP6124295B2 true JP6124295B2 (en) 2017-05-10

Family

ID=51936528

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013096813A Active JP6124295B2 (en) 2013-05-02 2013-05-02 Method for evaluating the texture of bare skin or makeup skin

Country Status (1)

Country Link
JP (1) JP6124295B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106573750B (en) * 2014-09-17 2018-09-28 三菱电机株式会社 Elevator device
JP6942012B2 (en) * 2017-09-12 2021-09-29 花王株式会社 Skin evaluation method
EP3664035B1 (en) * 2018-12-03 2021-03-03 Chanel Parfums Beauté Method for simulating the realistic rendering of a makeup product
JP7308443B2 (en) * 2019-03-19 2023-07-14 株式会社リコー Evaluation device, authenticity evaluation method and program

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4133248B2 (en) * 2002-11-18 2008-08-13 株式会社コーセー Skin gloss evaluation method
JP5146667B2 (en) * 2008-06-13 2013-02-20 花王株式会社 Evaluation method of skin smoothness
WO2011162050A1 (en) * 2010-06-21 2011-12-29 ポーラ化成工業株式会社 Age estimation method and gender determination method

Also Published As

Publication number Publication date
JP2014217456A (en) 2014-11-20

Similar Documents

Publication Publication Date Title
JP7235895B2 (en) Apparatus and method for visualizing cosmetic skin characteristics
US11605243B2 (en) Apparatus and method for determining cosmetic skin attributes
CA2678551C (en) Method and apparatus for simulation of facial skin aging and de-aging
JP6473401B2 (en) Skin gloss evaluation apparatus, gloss evaluation method, and gloss evaluation program
JP2009509708A (en) Method and apparatus for evaluating characteristics of defective part of skin and method for evaluating anti-aging effect of cosmetic product
US20080304736A1 (en) Method of estimating a visual evaluation value of skin beauty
JP6124295B2 (en) Method for evaluating the texture of bare skin or makeup skin
JP6530703B2 (en) Skin condition evaluation method
JP2010198382A (en) Simulation apparatus, simulation method and computer program
JP2008093048A (en) Elasticity evaluation method
JP5399874B2 (en) Image processing apparatus and image processing method
JP6730051B2 (en) Skin condition evaluation method
JP2008198156A (en) Hair image display method and display device
JP2007219911A (en) Method, system, and program for image processing
Kumari et al. Image fusion techniques based on pyramid decomposition
JP6393063B2 (en) Aging pore identification method and cosmetic evaluation method
JP7040889B2 (en) Image analysis method for distinguishing skin condition
JP5516772B2 (en) Display object
JP6099353B2 (en) Method for evaluating the effect of cosmetics to suppress the conspicuousness of dark colored areas of the skin
JP2023087699A (en) Level estimation method for face showing actual feeling of happiness, level estimation device for face showing actual feeling of happiness, and level estimation program for face showing actual feeling of happiness
JP2022169152A (en) Skin evaluation method, and powder evaluation method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160301

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170126

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170207

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170303

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: 20170328

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170329

R150 Certificate of patent or registration of utility model

Ref document number: 6124295

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250