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 PDF

Info

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
Application number
JP2015221697A
Other languages
Japanese (ja)
Other versions
JP2016104124A (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.)
Mikimoto Pharmaceutical Co Ltd
Original Assignee
Mikimoto Pharmaceutical 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 Mikimoto Pharmaceutical Co Ltd filed Critical Mikimoto Pharmaceutical Co Ltd
Publication of JP2016104124A publication Critical patent/JP2016104124A/en
Application granted granted Critical
Publication of JP6649749B2 publication Critical patent/JP6649749B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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.

特開2001−170028号公報JP 2001-170028 A 特開2006−61170号公報JP 2006-61170 A 特開2008−61892号公報JP 2008-61892 A

本発明の目的は、熟練した判断者を必要とせず、正確に、どこでも、迅速に皮膚表面のキメを評価できるようにするための画像処理方法とこれを用いた肌のキメの評価方法にある。   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.

実施例の画像の一例 画像A=デジタルビデオスコープで得た画像 画像B=画像Aより取り出したG(グリーン)の画像 画像C=画像Bにローリングボールアルゴリズム(半径50ピクセル)を実施した画像 画像D=画像Cにヒストグラム平坦化処理を実施した画像 画像E=画像Dを大津の方法で二値化した画像 画像F=画像Eにメディアンフィルタを実施した画像Image A = Image obtained with a digital video scope Image B = Image of G (green) extracted from Image A Image C = Image obtained by performing rolling ball algorithm (radius 50 pixels) on Image B Image D = Image C subjected to histogram flattening processing Image E = Image D obtained by binarizing image D using the Otsu method Image F = Image obtained by applying a median filter to image E 熟練者による目視スコアの代表例と実施例によって処理した画像Representative examples of visual scores by skilled workers and images processed by the examples 熟練者による目視スコアと実施例による非線形判別分析結果との比較Comparison between the visual score of the expert and the result of the nonlinear discriminant analysis by the examples

Claims (4)

皮膚表面のRGBで表現される画像ファイルからG(グリーン)を取り出し、取り出したG(グリーン)の画像に、ローリングボールアルゴリズム、ヒストグラム平坦化を実施した後二値化する皮膚表面のキメの画像処理方法。 G (green) is extracted from the image file of the skin surface expressed in RGB, and the extracted G (green) image is subjected to a rolling ball algorithm, histogram flattening, and then binarized image processing of the skin surface texture Method. さらに二値化した後、メディアンフィルタを実施する請求項1の皮膚表面のキメの画像処理方法。   2. The image processing method according to claim 1, wherein a median filter is performed after binarization. 請求項1乃至請求項2のいずれかの方法で処理した画像に非線形判別分析を行った結果をもってスコア化することを特徴とする皮膚表面のキメの評価方法。   A method for evaluating texture on the skin surface, characterized in that an image processed by the method according to any one of claims 1 and 2 is scored based on a result of performing a nonlinear discriminant analysis. 前記画像ファイルは、デジタル式マイクロスコープを用いて得られた画像ファイルである請求項1乃至請求項2のいずれかに記載の皮膚表面のキメの画像処理方法。3. The image processing method according to claim 1, wherein the image file is an image file obtained by using a digital microscope.
JP2015221697A 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 Active JP6649749B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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