(57)【要約】【課題】角質水分変化量を直接的に反映させた肌のシミュレーション画像を形成する。【解決手段】肌の元画像に対して角質水分量が変化した肌のシミュレーション画像の形成方法が、元画像の内部反射光画像と表面反射光画像を取得し、元画像の内部反射光画像の明るさを取得し、その明るさを、角質水分変化量（ΔＣ）の関数Ｆ（ΔＣ）から算出される調整量ΔＬで調整すると共に、元画像の内部反射光画像のＰＳＦ(Point-Spread Function)を取得し、そのＰＳＦを、角質水分変化量（ΔＣ）の関数Ｇ（ΔＣ）から算出される調整量ΔＰＳＦで調整し、一方、元画像の表面反射光画像の強度を取得し、その強度を角質水分変化量（ΔＣ）の関数Ｈ（ΔＣ）から算出される調整量ΔＳで調整し、明るさとＰＳＦを調整した内部反射光画像と、強度を調整した表面反射光画像とを再合成することからなる。【選択図】図１
2007105457.tif 000002
【特許請求の範囲】
【請求項１】
肌の元画像に対して角質水分量が変化した肌のシミュレーション画像の形成方法であって、
元画像の内部反射光画像と表面反射光画像を取得し、
元画像の内部反射光画像の明るさを取得し、その明るさを、角質水分変化量（ΔＣ）の関数Ｆ（ΔＣ）から算出される調整量（ΔＬ）で調整すると共に、元画像の内部反射光画像のＰＳＦ(Point-Spread Function)を取得し、そのＰＳＦを、角質水分変化量（ΔＣ）の関数Ｇ（ΔＣ）から算出される調整量（ΔＰ）で調整し、
元画像の表面反射光画像の強度を取得し、その強度を角質水分変化量（ΔＣ）の関数Ｈ（ΔＣ）から算出される調整量（ΔＳ）で調整し、
明るさとＰＳＦを調整した内部反射光画像と、強度を調整した表面反射光画像とを再合成してシミュレーション画像を形成する方法。
【請求項２】
元画像の表面反射光画像の強度の標準偏差を角質水分変化量（ΔＣ）の関数Ｈｄ（ΔＣ）から算出される調整量（ΔＳｄ）で調整することにより表面反射光画像の強度を調整する請求項１記載の方法。
【請求項３】
肌の内部反射光画像を色素成分画像と陰影成分画像に分離し、
内部反射光画像の明るさの調整とＰＳＦの調整を、それぞれ色素成分画像と陰影成分画像の一方又は双方について行う請求項１又は２記載の方法。
【請求項４】
複数の被験者の肌に所定の取り扱いを施し、その取り扱いの前後で角質水分量を測定することにより、その取り扱いと、それによる角質水分変化量との関係を求めておき、その関係に基づいて、当該被験者について肌に所定の取り扱いを施した場合の角質水分変化量を求め、その角質水分変化量から、被験者の肌に所定の取り扱いを施した場合の肌のシミュレーション画像を形成する請求項１〜３のいずれかに記載の方法。
【請求項５】
肌の元画像に対して角質水分量が変化した肌のシミュレーション画像を形成するシステムであって、
(1)肌の元画像、並びに、元画像の内部反射光画像及び表面反射光画像を形成することのできる画像形成手段、  
(2)元画像の内部反射光画像の明るさ（Ｌ0）を取得する手段、
(3)内部反射光画像のＰＳＦ（Ｐ0）を取得する手段、
(4)元画像の表面反射光画像の強度（Ｓ0）を取得する手段、並びに、
(5)角質水分変化量（ΔＣ）の関数Ｆ（ΔＣ）から内部反射光画像の明るさの調整量（ΔＬ）を算出し、この調整量（ΔＬ）で内部反射光画像の明るさ（Ｌ0）を調整すると共に、角質水分変化量（ΔＣ）の関数Ｇ（ΔＣ）から内部反射光画像のＰＳＦの調整量（ΔＰ）を算出し、この調整量（ΔＰ）で内部反射光画像のＰＳＦ（Ｐ0）を調整し、角質水分変化量（ΔＣ）の関数Ｈ（ΔＣ）から表面反射光画像の強度の調整量（ΔＳ）を算出し、この調整量（ΔＳ）で表面反射光画像の強度（Ｓ0）を調整し、明るさとＰＳＦを調整した内部反射光画像と強度を調整した表面反射光画像とを再合成してシミュレーション画像を形成する手段を有する演算手段、
を備えた肌のシミュレーション画像の形成システム。
【請求項６】
(4)元画像の表面反射光画像の強度（Ｓ0）を取得する手段と共に、該強度（Ｓ0）の標準偏差を取得する手段を備え、
(5)演算手段が、強度を調整した表面反射光画像として、表面反射光画像の強度の標準偏差を調整した画像を形成する請求項５記載の形成システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、角質水分量の変化量を肌のシミュレーション画像に反映させる方法とそのためのシステムに関する。
【背景技術】
【０００２】
従来より、被験者の画像を撮ってこれを元画像とし、画像処理により明度、色相、ＰＳＦ(Point-Spread Function)等の光学物性を変化させて元画像を修正することがなされている。
【０００３】
同様の画像処理により肌のシミュレーション画像を形成する方法も提案されており、例えば、肌に基礎化粧を施した場合の明度、色相等の光学物性値の変化に関するデータベースを構築しておき、そのデータベースに基づいて特定の基礎化粧を施した場合の光学物性値を得、その光学物性値を用いてシミュレーション画像を形成する方法（特許文献１）や、基礎化粧を施した場合のメラニン成分量やヘモグロビン成分量の変化に関するデータベースを構築しておき、そのデータベースに基づいて特定の基礎化粧を施したときのメラニン成分量やヘモグロビン成分量の変化量を得、これらの変化量に基づいて元画像の内部反射光の色素成分を制御し、シミュレーション画像を形成する方法（特許文献２）等がある。
【０００４】
一方、肌の角質水分量は、肌の透明感やはり等に影響することが知られており（特許文献３、特許文献４）、基礎化粧には保湿剤が広く使用されている。
【０００５】
【特許文献１】特開２００２−２０３２３８号公報
【特許文献２】特開２００２−２０００５０号公報
【特許文献３】特開２００４−２１５９９１号公報
【特許文献４】特開２００２−１０２１７７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、角質水分量が肌の見え方にどのように影響するかについては、大まかな傾向がわかっているにすぎず、数値として得られる角質水分変化量を直接的にシミュレーション画像に反映させることはこれまでになされていない。
【０００７】
そこで、本発明は、角質水分変化量を直接的に反映させた肌のシミュレーション画像を形成できるようにすることを目的とする。
【課題を解決するための手段】
【０００８】
本発明者らは、肌の角質水分変化量は、肌画像において、内部反射光画像の明るさとＰＳＦ、及び表面反射光画像の強度に影響し、これら３つの光学物性を角質水分変化量に基づいて調整すると、角質水分量が特定量変化した場合の画像を良好にシミュレーションできること、また更に、表面反射光画像の強度については、その標準偏差を角質水分変化量に基づいて調整すると、より良好なシミュレーション画像を形成できることを見出し、本発明を完成させた。
【０００９】
即ち、本発明は、肌の元画像に対して角質水分量が変化した肌のシミュレーション画像の形成方法であって、
元画像の内部反射光画像と表面反射光画像を取得し、
元画像の内部反射光画像の明るさを取得し、その明るさを、角質水分変化量（ΔＣ）の関数Ｆ（ΔＣ）から算出される調整量（ΔＬ）で調整すると共に、元画像の内部反射光画像のＰＳＦ(Point-Spread Function)を取得し、そのＰＳＦを、角質水分変化量（ΔＣ）の関数Ｇ（ΔＣ）から算出される調整量（ΔＰＳＦ）で調整し、
元画像の表面反射光画像の強度を取得し、その強度を角質水分変化量（ΔＣ）の関数Ｈ（ΔＣ）から算出される調整量（ΔＳ）で調整し、
明るさとＰＳＦを調整した内部反射光画像と、強度を調整した表面反射光画像とを再合成してシミュレーション画像を形成する方法を提供する。また更に、元画像の表面反射光画像の強度の標準偏差を角質水分変化量（ΔＣ）の関数Ｈｄ（ΔＣ）から算出される調整量（ΔＳｄ）で調整することにより表面反射光画像の強度を調整する態様を提供する。
【００１０】
また、本発明は、肌の元画像に対して角質水分量が変化した肌のシミュレーション画像を形成するシステムであって、
(1)肌の元画像、並びに、元画像の内部反射光画像及び表面反射光画像を形成することのできる画像形成手段、
(2)元画像の内部反射光画像の明るさ（Ｌ0）を取得する手段、
(3)内部反射光画像のＰＳＦ（Ｐ0）を取得する手段、
(4)元画像の表面反射光画像の強度（Ｓ0）を取得する手段、並びに、
(5)角質水分変化量（ΔＣ）の関数Ｆ（ΔＣ）から内部反射光画像の明るさの調整量（ΔＬ）を算出し、この調整量（ΔＬ）で内部反射光画像の明るさ（Ｌ0）を調整すると共に、角質水分変化量（ΔＣ）の関数Ｇ（ΔＣ）から内部反射光画像のＰＳＦの調整量（ΔＰ）を算出し、この調整量（ΔＰ）で内部反射光画像のＰＳＦ（Ｐ0）を調整し、角質水分変化量（ΔＣ）の関数Ｈ（ΔＣ）から表面反射光画像の強度の調整量（ΔＳ）を算出し、この調整量（ΔＳ）で表面反射光画像の強度（Ｓ0）を調整し、明るさとＰＳＦを調整した内部反射光画像と強度を調整した表面反射光画像とを再合成してシミュレーション画像を形成する手段を有する演算手段、
を備えた肌のシミュレーション画像の形成システムを提供する。また更に、(4)元画像の表面反射光画像の強度（Ｓ0）を取得する手段と共に、該強度（Ｓ0）の標準偏差を取得する手段を備え、
(5)演算手段が、強度を調整した表面反射光画像として、表面反射光画像の強度の標準偏差を調整した画像を形成するシステムを提供する。
【発明の効果】
【００１１】
本発明の肌のシミュレーション画像の形成方法によれば、元画像の内部反射光画像の明るさ及びＰＳＦ、並びに表面反射光画像の強度を、それぞれ角質水分変化量（ΔＣ）の関数から得られる調整量で調整するので、角質水分変化量（ΔＣ）を直接的に反映させた肌のシミュレーション画像を形成することができる。
【００１２】
ここで、内部反射光画像の明るさ及びＰＳＦ、並びに表面反射光画像の強度のいずれか一つでも調整をすることなくシミュレーション画像を形成すると、明るさ、マット感、透明性等が、角質水分量を実際にΔＣだけ変化させた肌の画像と異なって見える場合があるが、本発明によれば、内部反射光画像の明るさ及びＰＳＦ並びに表面反射光画像の強度の全てを角質水分変化量（ΔＣ）に基づいて調整するので、角質水分量を変化させた画像のシミュレーションを良好に行うことができる。
【００１３】
また、本発明の肌のシミュレーション画像の形成システムによれば、本発明の肌のシミュレーション画像の形成方法を良好に実施することができる。
【発明を実施するための最良の形態】
【００１４】
以下、図面を参照しつつ本発明を詳細に説明する。
【００１５】
図１は、本発明の一実施例の肌のシミュレーション画像の形成システム１００における処理の流れの説明図である。
【００１６】
このシステム１００は、画像形成手段１、演算手段２０及びディスプレイ３０を備え、さらに角質水分量計測装置４０、データベース５０及びプリンタ（図示せず）を備えている。
【００１７】
画像形成手段１は、白色通常光による肌画像（元画像）、集光照射による肌の内部散乱光画像、並びに肌の内部反射光画像及び表面反射光画像を形成するものである。かかる画像形成手段１は、任意の白色光源、集光光又は平行光を作ることのできるレンズ系、偏光板及びデジタルカメラから構成することができ、例えば、図２Ａに示したように、ハロゲン、メタルハライド等の白色光源２、白色光源２から導光した光を一端から出射し、点光源３として機能する光ファイバー（直径０．０１〜１ｍｍ）４、点光源３から出射された光を集光する集光レンズ５、集光レンズ５の後段に設けられたコリメータレンズ７、集光レンズ５とコリメータレンズ７との間に設けられた絞り６、コリメータレンズ７により形成された平行光を試料Ｓとする肌に集光する集光レンズ８、受光角（θ）が例えば４５°の向きに設置されたデジタルカメラ１０、点光源３とデジタルカメラの前にそれぞれ設けられた偏光板１１、１２を備えている。ここで、集光レンズ８、偏光板１１、１２は、着脱自在となっている。また、コリメータレンズ７と集光レンズ８との間にも、必要に応じて、光量調節及びレンズの収差低減のための絞りが設けられる（図示せず）。
【００１８】
この画像形成手段１において、図２Ａに示すように、集光レンズ８が光路に設置されている場合には、試料Ｓには集光光が照射されるが、図２Ｂに示すように、集光レンズ８が光路から外された場合には、試料Ｓには、コリメータレンズ７から平行光が照射されることとなる。そこで、この画像形成手段１において、集光照射による肌の内部散乱光画像は、図２Ａに示すように偏光板１１、１２を光路から外し、光路内に集光レンズ８を設置し、試料Ｓに集光光を照射し、その内部散乱光をデジタルカメラ１０で撮る。あるいは、点光源３の前の偏光板１１とデジタルカメラ１０の前の偏光板１２を、双方の偏光方向が直交するように装着して表面反射光成分が除去されるようにし、光路内に集光レンズ８を設置し、試料Ｓに集光光を照射し、その内部散乱光をデジタルカメラ１０で撮る。
【００１９】
一方、内部反射光画像は、図２Ｂに示すように、集光レンズ８を光路から外し、点光源３の前の偏光板１１とデジタルカメラ１０の前の偏光板１２を、双方の偏光方向が直交するように装着し、表面反射光成分を除去することにより形成し、表面反射光画像は、集光レンズ８を光路から外し、点光源３の前の偏光板１１とデジタルカメラ１０の前の偏光板１２を、双方の偏光方向が同じになるようにして装着することにより得た画像と、双方の偏光方向が直交するようにして得た画像との差分から形成する。
【００２０】
なお、内部反射光画像、表面反射光画像は、光源として図２Ｂに示すような平行光を使用することが望ましいが、レンズ等を使わず、通常の蛍光灯のような白色光源を使用してもよい。この場合にも、光源の前とカメラの前に偏光板を装着することにより内部反射光画像、表面反射光画像を得ることが出来る。
【００２１】
また、受光角（θ）が０°から９０°の範囲で所定の角度となるように、デジタルカメラ１０の向きを変えても良く、試料Ｓを回転させても良く、試料Ｓを照射する照射光の角度を変化させても良い。
【００２２】
点光源３として機能する光ファイバー４については、これを、光ファイバーを束ねたライトガイドに代えてもよい。ただし、光源の大きさを補正する操作が必要になることがある。
【００２３】
演算手段２０は、シミュレーション画像の形成に必要な演算機能及び画像処理機能を備えたパーソナルコンピュータからなっている。即ち、この演算手段２０は、元画像の内部反射光画像の明るさ（Ｌ0）の取得手段として、画像形成手段１で形成した内部反射光画像のＲＧＢ信号から、公知の方法によりＬ*ａ*ｂ*表色系のＬ*値やＸＹＺ表色系のＹ値等を算出する機能を有し、元画像の表面反射光画像の強度（Ｓ0）の取得手段として、画像形成手段１で形成した表面反射光画像のＲＧＢ信号の信号強度を算出する機能を有する。
【００２４】
また、演算手段２０は、元画像の内部反射光画像のＰＳＦの取得手段として、画像形成手段１で形成した集光照射による内部散乱光画像において内部散乱光の後方反射量を計測すること（ビデオリフレクトメトリ法）により減衰係数（μeff）を算出し、さらにＰＳＦを算出する機能を有する。より具体的には、例えば、集光照明による内部散乱光の動径方向の輝度分布を図３のように求め、必要に応じて、γ補正、暗電流補正等を行い、得られた輝度分布に対して、市販の科学的グラフ解析ソフトを使用し、次式（７）
【００２５】
【数１】
2007105457.tif 000003
（式中、ｒ：動径距離（独立変数）
ｙ：輝度（従属変数）
ｋ：パラメータ        ）
で非線形曲線にフィッティングし、減衰係数（μeff）を算出し、この減衰係数（μeff）からＰＳＦを取得する。
【００２６】
この他、元画像の内部反射光画像のＰＳＦを取得する手段としては、所謂ブラインドデコンボリューションによりＰＳＦを算出する機能を搭載してもよいが、ＰＳＦと角質水分変化量との対応の点から前述のように減衰係数（μeff）を算出し、さらにＰＳＦを算出する手法が好ましい。
【００２７】
演算手段２０は、さらにコンボリューション演算機能及びデコンボリューション演算機能を有する。ここで、コンボリューション（Convolution）は、畳み込み積分とも言われ、入力関数に対して重みづけを行い、出力関数を求める操作である。一方、デコンボリューション（Deconvolution）はコンボリューションの逆操作、即ち出力関数と重み関数とから入力関数を求める操作であり、劣化画像の復元に応用されている。
【００２８】
また、演算手段２０は、角質水分変化量（ΔＣ）に応じて、元画像の内部反射光画像の明るさ（Ｌ0）の調整量を算出する関数Ｆ（ΔＣ）及びＰＳＦ（Ｐ0）の調整量を算出する関数Ｇ（ΔＣ）、並びに表面反射光画像の強度（Ｓ0）の調整量を算出する関数Ｈ（ΔＣ）を有している。
【００２９】
これらの関数Ｆ（ΔＣ）、Ｇ（ΔＣ）、Ｈ（ΔＣ）は、肌に水分供給をしたり、肌を乾燥状態においたりすることにより肌の角質水分量を種々変化させた場合の内部反射光画像の明るさ（Ｌ）及びＰＳＦ（Ｐ）、並びに表面反射光画像の強度（Ｓ）と角質水分変化量（ΔＣ）との関係を多数取得し、そのデータ解析により得られるものであり、好ましくは、被験者の年齢、肌のタイプ（脂性か乾燥肌か）等に応じて関数を複数種持つことが好ましい。また、表面反射光画像の強度（Ｓ）については、その標準偏差と角質水分変化量（ΔＣ）との関数Ｈｄ（ΔＣ）を求め、それに基づいて表面反射光画像の強度の標準偏差を調整することが好ましい。
【００３０】
これらの関数Ｆ（ΔＣ）、Ｇ（ΔＣ）、Ｈ（ΔＣ）は、より具体的には、例えば、次の式（１）、式（３）、式（５）で表される。なお、このように内部反射光画像の明るさの調整量（ΔＬ）及びＰＳＦの調整量（ΔＰ）、並びに表面反射光画像の強度の調整量（ΔＳ）が、それぞれ角質水分変化量（ΔＣ）の関数で表されることは、本発明者が見出したものである。
【００３１】
【数２】
2007105457.tif 000004
（式中、ΔＣ：角質水分変化量
ΔＬ：内部反射光画像の明るさの調整量
ａ，ｂ：それぞれ係数  ）
【００３２】
【数３】
2007105457.tif 000005
（式中、ΔＣ：角質水分変化量
ΔＰ：内部反射光画像のＰＳＦの調整量
ｃ，ｄ：それぞれ係数  
μeff1：シミュレーション画像のμeff
μeff0：元画像のμeff
関数Ｇ’は、μeff1、μeff0のそれぞれのＰＳＦをＰ1(ｘ,ｙ)、Ｐ0(ｘ,ｙ)とし、
そのフーリエ変換をＦＰ1(ｕ,ｖ)、ＦＰ0(ｕ,ｖ)としたときの
ΔＦＰ(ｕ,ｖ)＝ＦＰ1(ｕ,ｖ)／ＦＰ0(ｕ,ｖ)
の逆フーリエ変換ΔＰ(ｘ,ｙ)を求める関数  ）
【００３３】
【数４】
2007105457.tif 000006
（式中、ΔＣ：角質水分変化量
ΔＳ：表面反射光画像の強度の調整量
ｅ，ｆ：それぞれ係数  ）
【００３４】
ここで、角質水分変化量（ΔＣ）は、角質水分量を変化させる前と後の角質水分量を所定の方法で測定した場合の測定値の差である。角質水分量の測定方法としては、コンダクタンス法、赤外法等をあげることができ、市販の測定器を使用することができる。なお、角質水分量の測定値は、測定方法によって異なる場合が多いので、本発明のシステムに備える角質水分量計測装置４０は、これらの関数の取得に使用した測定器と同種のものとすることが好ましい。
【００３５】
演算手段２０は、こうして算出される内部反射光画像の明るさの調整量（ΔＬ）及びＰＳＦの調整量（ΔＰ）、並びに表面反射光画像の強度の調整量（ΔＳ）に基づき、元画像の肌に対して角質水分量をΔＣだけ変化させた場合のシミュレーション画像の内部反射光画像と表面反射光画像を形成し、これを再合成してシミュレーション画像を完成する演算機能を有する。
【００３６】
より具体的には、画像処理により、シミュレーション画像の内部反射光画像の明るさ（Ｌ1）を次式（２）、ＰＳＦ（Ｐ1）を次式（４）、表面反射光画像の強度（Ｓ1）を次式（６ａ）又は（６ｂ）で得られる値に調整し、調整後の内部反射光画像と表面反射光画像を再合成する。
【００３７】
【数５】
2007105457.tif 000007
（式中、Ｌ1：シミュレーション画像の内部反射光画像の明るさ
Ｌ0：元画像の内部反射光画像の明るさ
ΔＬ：内部反射光画像の明るさの調整量  ）
【００３８】
【数６】
2007105457.tif 000008
（式中、
Ｉ1(ｕ,ｖ)：シミュレーション画像の内部反射光画像ｉ1(ｘ，ｙ)のフーリエ変換画像
Ｉ0(ｕ,ｖ)：元画像の内部反射光画像ｉ0(ｘ，ｙ)のフーリエ変換画像
ＦＰ0(ｕ,ｖ)：元画像の内部反射光画像のＰＳＦのフーリエ変換画像
ＦＰ1(ｕ,ｖ)：シミュレーション画像の内部反射光画像のＰＳＦのフーリエ変換画像
ΔＦＰ(ｕ,ｖ)：ＰＳＦの調整量ΔＰのフーリエ変換画像  ）
【００３９】
なお、式（４）による内部反射光画像のＰＳＦの調整は、元画像の内部反射光画像のフーリエ変換画像Ｉ0(ｕ,ｖ)を、まず元画像のＰＳＦのフーリエ変換画像ＦＰ0(ｕ,ｖ)でデコンボリューションしてＰＳＦの影響を取り除いた画像を形成し、それを算出されたＰＳＦのフーリエ変換画像ＦＰ1(ｕ,ｖ)でコンボリューションすることを意味しており、最終的にシミュレーション画像の内部反射光画像のフーリエ変換画像Ｉ1(ｕ,ｖ)を逆フーリエ変換して内部反射光画像ｉ1(ｘ，ｙ)を得、さらにその明るさを式（２）で調整してシミュレーション画像の内部反射光画像を得る。
【００４０】
【数７】
2007105457.tif 000009
（式中、Ｓ1：シミュレーション画像の表面反射光画像の強度
Ｓ0：元画像の表面反射光画像の強度
ΔＳ：表面反射光画像の強度の調整量
ｄ1：係数
ｄ2：係数
ｄ3：係数
ｇ：係数                              ）
【００４１】
これらの式（６ａ）と（６ｂ）は、シミュレーション画像における表面反射光画像の強度Ｓ1はΔＳの一次式でも、二次式でも表せることを意味している。いずれを使用するかは、求める精度によって適宜定める。
【００４２】
また、例えば、てかりの強い頬の肌の表面反射光画像とマットな頬の肌の表面反射光画像が、それぞれ図８の強度分布（ａ）、（ｂ）を示すように、表面反射光画像は、表面反射光の強度と度数のヒストグラムで表した場合の平均強度と標準偏差で特徴づけること（言い換えれば、明るさの平均値と標準偏差で特徴づけること）ができ、この標準偏差の角質水分変化量に基づく調整によっても、次式（８）のように、表面反射光画像の強度を調整することができる。
【００４３】
【数８】
2007105457.tif 000010
（式中、ΔＳｄ：表面反射光画像の強度の調整量
ｅｄ：係数
ｆｄ：係数）
【００４４】
標準偏差を調整した表面反射光画像は、市販の画像解析ソフト（例えば、Photoshop（Adobe社製））等により、肌の表面反射光画像のヒストグラムの標準偏差が所定量（Ｓｄ０＋ΔＳｄ；Ｓｄ０は表面反射光画像の強度の標準偏差の初期量を表す）になるようにコントラストを変化させる手法によっても得ることができる。
【００４５】
表面反射光画像の強度の調整は、好ましくは、平均強度と標準偏差の双方により行う。
【００４６】
一方、本発明のシステム１００においてデータベース５０は、多数の被験者の肌について、肌に水分を供給するか、あるいは肌を乾燥させる所定の取り扱いを施した場合に、その取り扱いの前後で肌の角質水分量、内部反射光画像の明るさ、ＰＳＦ、並びに表面反射光画像の強度を測定し、取り扱いとそれによる角質水分変化量（ΔＣ）および各画像の調整量ΔＬ、ΔＰ、ΔＳを求めるための係数ａ〜ｇ、ｄ1〜ｄ3との関係を、被験者の年齢、肌のタイプ（脂性肌か乾燥肌か）等ごとに蓄積し、データベース化したものである。このデータベース５０によれば、所定の角質水分変化量（ΔＣ）を指定することにより、その角質水分変化量（ΔＣ）をもたらす肌の取り扱い方法を得ることができ、また、所定の取り扱い方法を指定することにより、その取り扱い方法によりもたらされる角質水分変化量（ΔＣ）を得ることができる。
【００４７】
そして、このシステム１００では、元画像の肌に対して施す取り扱いおよび被験者の年齢、肌のタイプ（脂性肌か乾燥肌か）等を演算手段２０に入力すると、それに対応する角質水分変化量（ΔＣ）および各画像の調整量ΔＬ、ΔＰ、ΔＳを求めるための係数ａ〜ｇ、ｄ1〜ｄ3がデータベース５０から出力され、これに基づいて内部反射光画像の明るさの調整量（ΔＬ）及びＰＳＦの調整量（ΔＰ）、並びに表面反射光画像の強度の調整量（ΔＳ）が算出されるようになっている。
【００４８】
ディスプレイ３０には、画像形成手段１で形成された元画像、演算手段２０で形成されたシミュレーション画像、あるいはその形成の過程で使用される内部反射光画像、表面反射光画像等が、適宜切り替え、あるいは同時に表示される。
【００４９】
本発明の肌のシミュレーション画像の形成方法は、例えば、上述のシステム１００を用いて、図１のフローで次のように行われる。
【００５０】
まず、シミュレーション画像を形成する被験者の肌の通常光による画像（元画像）を、画像形成手段１により、白色光源下で偏光板を使用することなく撮る。また、被験者の肌の角質水分量を角質水分量計測装置４０で測定し、被験者の肌が乾燥状態にあるのか、湿潤状態にあるのか等を知り、例えば、乾燥状態にある場合には、所定の保湿化粧料の塗布を所定期間行った場合に増加する角質水分変化量（ΔＣ）および各画像の調整量ΔＬ、ΔＰ、ΔＳを求めるための係数ａ〜ｇ、ｄ1〜ｄ3をデータベース５０から得、この角質水分変化量（ΔＣ）に対応するシミュレーション画像を次にように形成する。即ち、まず、画像形成手段１において、偏光板１１、１２を用い、白色平行照明下で被験者の肌の画像を撮ることにより、その肌の内部反射光画像と表面反射光画像を得る。
【００５１】
また、元画像の内部反射光画像のＰＳＦを得るために、画像形成手段１で偏光板１１、１２を使用せずに肌に白色集光光を照射してその内部散乱光画像を撮り、演算手段２０で、内部散乱光の後方反射量を計測して図３に示すような動径方向の輝度分布を得、前述の式（７）で非線形曲線にフィッティングし、減衰係数（μeff）を算出する。また、肌に白色集光光を照射してその内部散乱光画像よりＰＳＦ（Ｐ0）を得る。
【００５２】
なお、内部反射光画像のＰＳＦはこの方法に限らず、例えば、所謂ブラインドデコンボリューションによりＰＳＦを算出してもよく、次の（ｉ）〜（iii）方法で取得してもよい。
（ｉ）測定対象物の等価散乱係数μs’（あるいは散乱係数μsと非等方散乱係数g）及び吸収係数μaを取得し、モンテカルロシミュレーション（Lihong Wang, Steven L. Jacques; http://omlc.ogi.edu/pubs/pdf/man_mcml.pdf）により、ＰＳＦを求める。
（ii）測定対象物を構成する散乱粒子径、粒子および媒体の屈折率を取得し、Mie散乱計算（Scott Prahl；http://omlc.bme.ogi.edu/software/mie/）により、等価散乱係数μs’（あるいは散乱係数μsと非等方散乱係数g）及び吸収係数μaを求め、(1)と同様にＰＳＦを求める。
（iii）測定対象物の等価散乱係数μs'（あるいは散乱係数μsと非等方散乱係数g）及び吸収係数μaを取得し、解析解（S. L. Jacques, A. Gutsche, J. A. Schwartz, L. H. Wang and F. K. Tittel: ゛Video-reflectometry to extract optical properties of tissue in vivo," Medical Optical Tomography: Functional Imaging and Monitoring, IS11 of SPIE Institute Series (SPIE, Bellingham, 1993) pp. 211-226）により、ＰＳＦを求める。
【００５３】
さらに、元画像の内部反射光画像の明るさ（Ｌ0）と表面反射光画像の強度（Ｓ0）を、それぞれ内部反射光画像と表面反射光画像のＲＧＢ信号から算出しておく。この場合、内部反射光画像の明るさとしては、Ｌ*ａ*ｂ*表色系のＬ*値やＸＹＺ表色系のＹ値等を使用する。
【００５４】
次に、元画像の内部反射光画像の明るさ（Ｌ0）とＰＳＦ（Ｐ0）をそれぞれ角質水分変化量（ΔＣ）で調整し、シミュレーション画像の内部反射光画像を得る。
【００５５】
ここで、明るさの調整は、前述の式（１）により明るさの調整量（ΔＬ）を求め、式（２）の明るさ（Ｌ1）に調整することにより行う。
【００５６】
また、ＰＳＦの調整は、前述の式（３）によりＰＳＦの調整量（ΔＰ）を求め、式（４）により調整後のＰＳＦ（Ｐ1）を有するシミュレーション画像を算出する。なお、元画像に対する明るさの調整とＰＳＦの調整は、いずれを先にしてもよい。
【００５７】
また、元画像の内部反射光に対する明るさとＰＳＦの調整は、元画像をさらに独立成分分析により色素成分画像（即ち、メラニン成分画像とヘモグロビン成分画像）と陰影成分画像に分離し、色素成分画像と陰影成分画像のいずれについても、明るさとＰＳＦの調整を行ってもよい。これにより、メラニン成分画像とヘモグロビン成分画像のそれぞれを、より適した明るさとＰＳＦに調整することができる。
【００５８】
なお、独立成分分析とは、皮膚の層構造をメラニンを主な色素成分として含有する表皮層と、ヘモグロビンを主な色素成分として含有する真皮層と、その他の色素成分を含有する皮下組織との積層構造であるとモデル化し、各層から独立的に信号が発せられ、それらが混合したものが画像信号になっていると考え、画像信号から各層の信号を分離抽出する分析方法である。この解析処理と画像処理の詳細はVol. 16, No. 9/  September 1999/ J. Opt. Soc. Am. A 2169に記載されており、パーソナルコンピュータに、市販の画像解析ソフト（例えば、Photoshop（Adobe社製））を搭載することにより行うことができる。
【００５９】
一方、元画像の表面反射光画像については、まず、式（５）で強度の調整量（ΔＳ）を求め、式（６ａ）又は式（６ｂ）で元画像の強度（Ｓ0）を調整し、シミュレーション画像の表面反射光画像を得る。
【００６０】
こうして得られたシミュレーション後の内部反射光画像と表面反射光画像を再合成することによりシミュレーション画像を得、ディスプレイ３０に表示する。
【００６１】
この肌のシミュレーション画像の形成方法及びシステムは、美容アドバイス、化粧料の推奨、化粧料の開発等で有用となる。例えば、美容アドバイザーは、化粧料の使用あるいは不使用による角質水分量の変化で肌の見え方が変わることを被験者に説明し易くなり、被験者は、化粧料の使用あるいは不使用により角質水分変化量が変化した状態の自分の肌画像を実感することが可能となり、化粧料の研究者は化粧料がもたらす角質水分変化量と、それによる肌の見え方との関係を活かして製品開発をすることが可能となる。
【実施例】
【００６２】
以下、実施例に基づき、本発明を具体的に説明する。
【００６３】
実施例１
２５名の被験者（２５〜３９歳、女性、乾燥肌）の顔の肌にスキンケア製剤を４週間連続使用するスキンケア試験をし、このスキンケア試験の前後で、各被験者の肌画像を偏光板を使用することなく白色光源下でデジタルカメラで撮ると共に、偏光板を用いて肌の内部反射光画像と表面反射光画像を撮った。さらに、頬部の角質水分量を水分測定器（ＩＢＳ社製  ＳＫＩＮＣＯＮ−２００）で測定した。この水分測定により、角質水分量（コンダクタンス）の平均値が４８μＳから７５μＳに上昇したことがわかった。
【００６４】
また、スキンケア試験の前後の肌について、(a)分光反射率の測定、(b)内部反射光画像の画像信号に基づく明るさ（Ｌ*ａ*ｂ*表色系のＬ*）の計測、(c)表面反射光画像の画像信号に基づく強度の計測、(d)集光光の照射による内部散乱光の動径方向の輝度分布曲線の取得と、それに基づく減衰係数（μeff）とＰＳＦの算出を行った。これらの結果を図４〜図７に示す。
【００６５】
図４から、４週間のスキンケア試験により、反射率が全体として高くなり、白色成分が増加しているといえるから、このスキンケア試験により色素成分の増減はなく、色素以外の要因で肌が明るくなったことがわかる。
【００６６】
図５から、スキンケア試験により内部反射光画像の明るさ（Ｌ*）が明るくなり、図６から１／μeff（即ち、光の広がりやすさ）が大きくなり、図７から表面反射光画像の強度が小さくなっていることがわかる。
【００６７】
以上により、スキンケア試験による角質層の水分増加により角質層の散乱性が変化して１／μeffが増加し、内部反射光が増えて内部反射光画像が明るくなり、また、表面反射光が減少してその強度が低下したと考えられる。
【００６８】
次に、上記２５名と異なる特定の被験者について、スキンケア試験前の画像を元画像として、角質水分変化量が特定量（ΔＣ）の場合のシミュレーション画像を次のように形成した。
【００６９】
まず、図５の結果から、内部反射光画像の明るさの調整量（ΔＬ）と角質水分変化量（ΔＣ）との関係式として、次式１’を得、
【００７０】
【数９】
2007105457.tif 000011
【００７１】
図６からスキンケア試験前後の減衰係数（μeff）を求め、これと角質水分変化量（ΔＣ）との関係式として次式（３’）を得、次いで前述のように内部反射光画像のＰＳＦの調整量（ΔＰ）を求めた。
【００７２】
【数１０】
2007105457.tif 000012
【００７３】
図７の結果から、表面反射光画像の強度の調整量（ΔＳ）と角質水分変化量（ΔＣ）との関係式として、次式５’を得た。
【００７４】
【数１１】
2007105457.tif 000013
【００７５】
そして、スキンケア試験前の内部反射光画像の明るさ（Ｌ0）を次式（２）で調整し、
【００７６】
【数１２】
2007105457.tif 000014
【００７７】
前述の次式（４）により、スキンケア試験後のＰＳＦ（Ｐ1）を有するシミュレーション画像の内部反射光画像を形成した。
【００７８】
【数１３】
2007105457.tif 000015
【００７９】
【数１４】
2007105457.tif 000016
【００８０】
こうして形成したシミュレーション画像の内部反射光画像と表面反射光画像を合成し、シミュレーション画像を形成した。
【００８１】
比較例１
内部反射光画像の明るさの調整を行わない以外は実施例１と同様にして、シミュレーション画像を形成した。
【００８２】
比較例２
内部反射光画像のＰＳＦの調整を行わない以外は実施例１と同様にして、シミュレーション画像を形成した。
【００８３】
比較例３
表面反射光画像の強度の調整を行わない以外は実施例１と同様にして、シミュレーション画像を形成した。
【００８４】
対比
実施例１及び比較例１〜３で得られたシミュレーション画像と、実際にスキンケア試験後に撮った被験者の肌画像とを対比した。その結果、比較例１のシミュレーション画像は実際のスキンケア試験後の肌画像に対して暗く、比較例２のシミュレーション画像はしわや毛穴が目立つ画像となり、比較例３のシミュレーション画像はつやが目立つとなった。これに対して、実施例１のシミュレーション画像は、明るさ、透明感、マット感のいずれについても実際のスキンケア試験後の肌画像を良好に再現できていた。
【産業上の利用可能性】
【００８５】
本発明は、角質水分量変化量を直接的に反映させたシミュレーション画像を形成するものであり、美容アドバイス、化粧料の推奨、化粧料の開発等で有用となる。
【図面の簡単な説明】
【００８６】
【図１】肌のシミュレーション画像の形成システムにおける処理の流れの説明図である。
【図２Ａ】画像形成手段の構成図である。
【図２Ｂ】画像形成手段の構成図である。
【図３】集光照明による内部散乱光の動径方向の輝度分布図である。
【図４】スキンケア試験前後における分光反射率の測定図である。
【図５】スキンケア試験前後の肌の内部反射光画像の明るさを示す図である。
【図６】スキンケア試験前後の肌の内部反射光画像の光の広がりやすさ（１／μeff）を示す図である。
【図７】スキンケア試験前後の肌の表面反射光画像の強度を示す図である。
【図８】てかりのつよい頬の肌の表面反射光画像における表面反射光強度と度数のヒストグラム（同図ａ）及びマットな頬の肌の表面反射光画像における表面反射光強度と度数のヒストグラム（同図ｂ）である。
【符号の説明】
【００８７】
１    画像形成手段
２    白色光源
３    点光源
４    光ファイバー
５    集光レンズ
６    絞り
７    コリメータレンズ
８    集光レンズ
１０    デジタルカメラ
１１    偏光板
１２    偏光板
２０    演算手段
３０    ディスプレイ
４０    角質水分量計測装置
５０    データベース
１００    肌のシミュレーション画像の形成システム
【図１】
2007105457.tif 000017
【図２Ａ】
2007105457.tif 000018
【図２Ｂ】
2007105457.tif 000019
【図３】
2007105457.tif 000020
【図５】
2007105457.tif 000021
【図６】
2007105457.tif 000022
【図７】
2007105457.tif 000023
【図４】
2007105457.tif 000024
【図８】
2007105457.tif 000025