JP6649749B2 - Image processing method for texture on skin surface and method for evaluating texture on skin surface using the same - Google Patents
Image processing method for texture on skin surface and method for evaluating texture on skin surface using the same Download PDFInfo
- Publication number
- JP6649749B2 JP6649749B2 JP2015221697A JP2015221697A JP6649749B2 JP 6649749 B2 JP6649749 B2 JP 6649749B2 JP 2015221697 A JP2015221697 A JP 2015221697A JP 2015221697 A JP2015221697 A JP 2015221697A JP 6649749 B2 JP6649749 B2 JP 6649749B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- skin surface
- texture
- image processing
- evaluating
- 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
Description
本発明は、皮膚表面のキメ(肌理)を評価するための画像処理方法と、その画像を用いた肌のキメの評価方法に関する。 The present invention relates to an image processing method for evaluating texture (texture) on the skin surface, and a method for evaluating skin texture using the image.
皮膚表面からの視覚的情報はシミ、ソバカス等の色の情報の他に形状による情報も与える。形状にはシワやタルミもあるが、キメ(肌理)も皮膚の恒常性、美観を表す重要な要素である。
キメの評価は種々の因子が関わっているので、一般的には熟練者が判断していたが、単位面積あたりの交点数を計測する方法(特許文献1)、画像に細線化処理等を行う方法(特許文献2)、短直線マッチング法での画像処理(特許文献)等、人を介さずに判断する方法が知られている。
特に、最近では、化粧品販売の現場で消費者の皮膚等を迅速に分析し、最適な化粧品を販売できるように求められている。
Visual information from the skin surface gives information on the shape in addition to color information such as spots and freckles. There are wrinkles and tarmi in the shape, but texture (texture) is also an important factor that expresses the homeostasis and beauty of the skin.
Since various factors are involved in the evaluation of texture, a skilled person generally judges it. However, a method of measuring the number of intersection points per unit area (Patent Document 1), thinning processing on an image, and the like are performed. A method of making a determination without human intervention, such as a method (Patent Document 2) and image processing by a short straight line matching method (Patent Document), is known.
In particular, recently, there has been a demand for the ability to quickly analyze the skin and the like of a consumer at a cosmetics sales site and to sell the most suitable cosmetics.
本発明の目的は、熟練した判断者を必要とせず、正確に、どこでも、迅速に皮膚表面のキメを評価できるようにするための画像処理方法とこれを用いた肌のキメの評価方法にある。 An object of the present invention is to provide an image processing method and an evaluation method for skin texture using the same, which do not require a skilled judge, accurately, anywhere, and can quickly evaluate the texture of the skin surface. .
本発明者らが鋭意検討した結果、画像処理方法としては、
皮膚表面のRGBで表現される画像ファイルからG(グリーン)を取り出し、取り出したG(グリーン)の画像に、ローリングボールアルゴリズム、ヒストグラム平坦化を実施した後二値化する方法が、本課題を解決すつ手段として最適であることを見出した。
に二値化した後、メディアンフィルタを実施するとさらに良好な結果が得られる場合がある。
上記の画像処理を行った画像に非線形判別分析を行った結果をもってスコア化することによって皮膚表面のキメの評価方法を見出した。
As a result of extensive studies by the present inventors, as an image processing method,
This method solves this problem by extracting G (green) from an image file expressed in RGB of the skin surface, performing a rolling ball algorithm, histogram flattening on the extracted G (green) image, and then binarizing the image. It has been found that it is most suitable as a means of fixing.
After the binarization, a better result may be obtained by performing a median filter.
A scoring method was performed based on the result of performing the non-linear discriminant analysis on the image subjected to the image processing described above, thereby finding a method for evaluating the texture of the skin surface.
詳しく説明すると、皮膚表面のRGBで表現される画像ファイルはデジタル式のマイクロスコープ等を用いることにより得ことができ、デジタル式マイクロスコープとしては、例えば、ショットモリテックス株式会社社製、チャームビュー、アイスコープR USB2.0等、スカラ株式会社社製、USBマイクロスコープ M3、ワイヤレススコープ AirMicro A1(U)、ワイヤレススコープ AirMicro A1等、株式会社フォルテシモ社等、USBマイクロスコープM3等が挙げられる。 More specifically, an image file expressed in RGB of the skin surface can be obtained by using a digital microscope or the like. USB Microscope M3, such as Scope R USB2.0, manufactured by Scalar Corporation, USB Microscope M3, Wireless Scope AirMicro A1 (U), Wireless Scope AirMicro A1, etc., Fortissimo Co., Ltd., etc.
上記のマイクロスコープ等で皮膚表面のRGBで表現される画像ファイルを得たら、この画像ファイルからG(グリーン)を取り出す。
取り出したG(グリーン)の画像に、ローリングボールアルゴリズムと、ヒストグラム平坦化を行う。
ローリングボールアルゴリズムは背景除去の1つの方法であり、ローリングボールの半径は画像の大きさ等によって最適な大きさを選択する。
ヒストグラム平坦化は、イコライゼーション処理ともいい、ヒストグラムの分布を均一になるように変換する処理である。
これらの処理後に二値化し、二値化した後、必要に応じて、ノイズ除去を実施する。
二値化の方法にはP-タイル法、モード法、判別分析法、最小誤差法があり、その中から選択して用いればよいが、判別分析法(大津の方法)がよく用いられる。
ノイズ除去にはいろいろ方法があり移動平均法、メディアンフィルター等がよく用いられるが、本発明にはメディアンフィルターがもっとも好ましかった。
メディアンフィルターとは、n×nの局所領域における濃度値を小さい順に並べ、 真ん中にくる濃度値を領域中央の画素の出力濃度とする処理である。
When an image file expressed in RGB of the skin surface is obtained with the above-mentioned microscope or the like, G (green) is extracted from this image file.
A rolling ball algorithm and histogram flattening are performed on the extracted G (green) image.
The rolling ball algorithm is one method of background removal, and the radius of the rolling ball selects an optimum size according to the size of the image and the like.
Histogram flattening is also called equalization processing, and is processing for converting the distribution of histograms so as to be uniform.
After these processes, binarization is performed, and after binarization, noise removal is performed as necessary.
The binarization method includes a P-tile method, a modal method, a discriminant analysis method, and a minimum error method, which may be selected and used, but the discriminant analysis method (Otsu's method) is often used.
There are various methods for removing noise, and a moving average method, a median filter, and the like are often used. In the present invention, a median filter is most preferred.
The median filter is a process in which density values in an n × n local area are arranged in ascending order, and a density value in the middle is set as an output density of a pixel in the center of the area.
上記のような方法で得た画像より、皮丘の面積とその変動係数、皮溝の幅や方向性等の各種のデータを得ることができる。
以上の処理は、公に利用可能なソフトウェア(ImageJ API/ライブラリー;http://rsb.info.nih.gov/ij/、NIH、MD)を使用して実施することができ、解析マクロ等を組むことによって、自動的に実施することができる。
From the image obtained by the above-described method, it is possible to obtain various data such as the area of the skin ridge and its variation coefficient, the width and directionality of the skin sulcus.
The above processing can be carried out using publicly available software (ImageJ API / library; http://rsb.info.nih.gov/ij/, NIH, MD), including analysis macros Can be implemented automatically.
さらに、本発明者は上記の画像より得たデータを用いて非線形判別分析により、熟練者しか行えなかった皮膚表面のキメの状態の分類を人を介することなく、自動的に行うことができた。
非線形判別分析とは、データ群を判別対象となる2つ以上のグループ(判別群)に判別するためのアルゴリズムであり、本発明では、独立変数として、「皮丘の平均面積」および「皮丘の面積の変動係数」を用い、この「皮丘の平均面積」と「皮丘の面積の変動係数」とにより被験者のデータの散布図を描いたときに、被験者のデータを2以上の判別群に分割するのに最適な関数(非線形判別関数)をそれぞれ求める。そして、未知データが入力された際に、その未知データがどのグループ群に属するかを非線形判別関数を基準にして決定することで、皮膚表面のキメの状態を判別できる。
判別群の群分けは目的によって変わるが、通常3〜20、よく利用されるのは、4〜10程度に群分けされることが多い。
非線形判別分析は、SAS(SAS社)、Stata(ライトストーン社)、SPSS(IBM社)、S_PLUS(数理システム社)等が利用できるし、また、フリーソフトである「R」でも分析できる。
Furthermore, the present inventor was able to automatically classify the texture of the skin surface, which could only be performed by a skilled person, without human intervention, by non-linear discriminant analysis using the data obtained from the above images. .
The non-linear discriminant analysis is an algorithm for discriminating a data group into two or more groups (discrimination groups) to be discriminated. In the present invention, as an independent variable, “average area of a crust” and “a crust” When a scatter diagram of the subject's data is drawn using the “average area of the crust” and the “coefficient of variation of the crust area” using the “coefficient of variation of the area of An optimal function (non-linear discriminant function) to divide into is calculated. Then, when unknown data is input, it is possible to determine the texture state of the skin surface by determining which group group the unknown data belongs to with reference to the nonlinear discriminant function.
The classification of the discrimination group varies depending on the purpose, but is usually 3 to 20, and is often used in about 4 to 10 in many cases.
For the non-linear discriminant analysis, SAS (SAS), Stata (Lightstone), SPSS (IBM), S_PLUS (Mathematical Systems) and the like can be used, and the analysis can be performed using free software "R".
次にに実施例を挙げて説明する。
1−1.女性の頬部をデジタルビデオスコープで皮膚表面をJPEG形式の画像ファイルを得た。(78名で実施した)
ImageJを用いて、上記RGBで表現される画像ファイルに以下の操作を行った。
1−2.RGBで表現される画像ファイルをRGBに分解し、G(グリーン)の画像を取り出した。
1−3.取り出したG(グリーン)の画像に、ローリングボールアルゴリズムを実施した。
1−4.ヒストグラム平坦化処理をした。
1−5.大津の方法によって二値化した。
1−6.メディアンフィルタ(3×3の画素の平均の値)を実施した。
以上の方法で、皮膚表面のキメを評価できる画像を得た。
Next, an example will be described.
1-1. A JPEG image file was obtained on the skin surface of the female cheek with a digital videoscope. (Conducted by 78 people)
Using ImageJ, the following operation was performed on the image file expressed in RGB.
1-2. An image file expressed in RGB was decomposed into RGB, and a G (green) image was extracted.
1-3. A rolling ball algorithm was performed on the extracted G (green) image.
1-4. The histogram was flattened.
1-5. It was binarized by Otsu's method.
1-6. A median filter (average value of 3 × 3 pixels) was performed.
By the above method, an image capable of evaluating the texture of the skin surface was obtained.
2−1.この画像より皮丘の面積の平均値と変動係数を得た。
この数値と、後述するスコアを教師データとして、ソフト「R」で非線形判別分析を行った。
デジタルビデオスコープで得たImageJで操作する前の画像を熟練者に、肌理がほとんどない状態(スコア1)から理想的な肌理の状態(スコア4)までの4つにクラス分けを行ってもらった。
2-1. From this image, the average value and the coefficient of variation of the area of the crust were obtained.
Non-linear discriminant analysis was performed by software “R” using this numerical value and a score described later as teacher data.
An image obtained by a digital videoscope before being operated with ImageJ was classified by an expert into four classes from a condition with almost no texture (score 1) to an ideal condition with a texture (score 4). .
非線形判別分析の結果、4つの判別群に分割するのに最適な非線形判別関数を得た。
確認のため上記で求めた非線形判別関数を用いて、上記の78のデータを4つにクラス分を再度行った結果、図3のように、熟練者による目視スコアとの一致率は69%で、両者の間で2段階以上異なる結果はなく、熟練者による目視に代わる手段として充分に利用できることがわかった。
As a result of the nonlinear discriminant analysis, an optimal nonlinear discriminant function for dividing into four discriminant groups was obtained.
For confirmation, the above-mentioned 78 data were re-classified into four classes using the nonlinear discriminant function obtained above, and as a result, as shown in FIG. 3, the coincidence rate with the visual score by an expert was 69%. There were no results that differed by more than two stages between the two, indicating that the two could be used sufficiently as a means to replace visual observation by a skilled person.
Claims (4)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014234069 | 2014-11-19 | ||
JP2014234069 | 2014-11-19 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016104124A JP2016104124A (en) | 2016-06-09 |
JP6649749B2 true JP6649749B2 (en) | 2020-02-19 |
Family
ID=56101999
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015221697A Active JP6649749B2 (en) | 2014-11-19 | 2015-11-12 | Image processing method for texture on skin surface and method for evaluating texture on skin surface using the same |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6649749B2 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109859257B (en) * | 2019-02-25 | 2021-12-17 | 北京工商大学 | Skin image texture evaluation method and system based on texture directionality |
CN110010228B (en) * | 2019-03-26 | 2022-12-23 | 广州艾颜佳美容美发设备有限公司 | Face skin perspective algorithm based on image analysis |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006061170A (en) * | 2004-07-30 | 2006-03-09 | Pola Chem Ind Inc | Discrimination method for skin |
JP5149527B2 (en) * | 2007-03-29 | 2013-02-20 | 常盤薬品工業株式会社 | How to display skin pigmentation |
JP2009018048A (en) * | 2007-07-12 | 2009-01-29 | Fujifilm Corp | Medical image display, method and program |
JP5651412B2 (en) * | 2010-09-08 | 2015-01-14 | 株式会社東芝 | Image processing apparatus, program, and image processing method |
JP6058902B2 (en) * | 2012-03-13 | 2017-01-11 | 株式会社 資生堂 | Skin condition analysis method, skin condition analysis apparatus, skin condition analysis system, program for executing the skin condition analysis method, and recording medium recording the program |
-
2015
- 2015-11-12 JP JP2015221697A patent/JP6649749B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2016104124A (en) | 2016-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6483826B2 (en) | Image analysis apparatus, image analysis method, image analysis system, image analysis program, and recording medium | |
TWI715877B (en) | Method and system for marking bone marrow cells | |
US20170140246A1 (en) | Automatic glandular and tubule detection in histological grading of breast cancer | |
US20170330315A1 (en) | Information processing apparatus, method for processing information, discriminator generating apparatus, method for generating discriminator, and program | |
WO2020062088A1 (en) | Image identification method and device, storage medium, and processor | |
CN109544583B (en) | Method, device and equipment for extracting interested area of leather image | |
Laga et al. | Image-based plant stornata phenotyping | |
CN113935666B (en) | Building decoration wall tile abnormity evaluation method based on image processing | |
JP2016507229A5 (en) | ||
Riana et al. | Extraction and classification texture of inflammatory cells and nuclei in normal pap smear images | |
CN103268492B (en) | A kind of corn grain type identification method | |
JP6649749B2 (en) | Image processing method for texture on skin surface and method for evaluating texture on skin surface using the same | |
JP2009075751A5 (en) | ||
JP6819445B2 (en) | Information processing equipment, control methods, and programs | |
JP2006061170A (en) | Discrimination method for skin | |
CN115362445A (en) | Data sorting device, data sorting method, and data sorting program | |
Nguyen et al. | A new method for splitting clumped cells in red blood images | |
JP5286574B2 (en) | Object detection recognition apparatus, object detection recognition method, and program | |
Ramaswamy et al. | A depth based approach to glaucoma detection using retinal fundus images | |
EP2776974B1 (en) | Methods and apparatus for image analysis using threshold compactness features | |
Hsu et al. | Vision-based detection of steel billet surface defects via fusion of multiple image features | |
CN115705748A (en) | Facial feature recognition system | |
CN108133210B (en) | Image format identification method and device | |
Jahani et al. | A novel method for centromere and length detection in microscopic images of human chromosomes | |
CN117474915B (en) | Abnormality detection method, electronic equipment and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20181010 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190617 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190813 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20191004 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191024 |
|
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: 20200114 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200117 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6649749 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |