メーキャップ化粧品の光学的機能評価に向けた人工皮膚の開発
芝浦工業大学工学部機械工学科

山田純

In the current cosmetics industry, when developing make-up cosmetics and evaluating its appearance, a visual assessment
has been performed by applying the cosmetics to human’s skin. An evaluation using artificial skin using a resin material
has been attempted, but it is not used for final evaluation of cosmetics. Since the present artificial skin cannot completely
reproduce the characteristics of light propagation inside and surface of the living human skin, the appearance of the artificial
skin is different from the living skin. In this study, we proposed a method that can individually adjust three optical properties
and surface structures to create the artificial skin reproducing the appearance of the living human skin.
As for the absorption coefficient, it was shown that the mass concentration of dye and the absorption coefficient havelinearity, so that the absorption coefficient can be controlled. With regard to the scattering coefficient and phase function,
it was suggested that the scattering coefficient and phase function of the skin can be reproduced by controlling the particle
diameter and particle volume fraction based on Mie theory. From these results, it was shown that artificial skin reproducing
the optical properties of living human skin can be created by individually controlling the absorption coefficient by some dyes
and the scattering coefficient and scattering phase function by white particles.
Three machining methods such as cutting, sandblasting and shaking were tried on the A 5052 substrates to reproduce
the characteristics of surface reflection on the human skin. Subsequently to the developing the surface structure, the surface
structure data was measured using a laser microscope. These data indicated that it is possible to create a structure that can
reproduce the optical characteristics on the human skin surface by combined several machining methods appropriately.

1. 緒

た場合，実際の人の皮膚の持つ透明感などの質感を感じら

言れないことが課題とされていた．そこで，近年では皮膚現在の化粧品業界では，メーキャップ化粧品を開発し，内部での散乱光を二方向散乱面反射分布関数（BSSRDF：

その見栄えを評価する場合，実際に人の皮膚に化粧品を塗

Bidirectional Scattering Surface Reflectance Distribution

1 ‒ 3）．しFunction）を CG に適用したモデルが注目されるようになかし，この方法では被験者ごとの肌の違いや，同一被験者ってきた 6 , 7）．これらの研究結果から，皮膚─空気界面でであっても，肌の状態が変化するために，定量的な評価はの光反射性質だけでなく，皮膚内部での散乱及び吸収（光困難である．

物性）を考慮した反射光を再現することにより実際の人の布した際の仕上がりを視感によって評価している

これまでに，樹脂素材を使用した人工皮膚が製作されて皮膚に近い質感が再現されることが示された．

おり，それを用いた評価も試みられている．しかし，化粧人工皮膚作製に関する従来の研究は，化粧分野では，化品の最終評価には利用されていない．これは，現在ある人

粧品の塗りやすさの試験や評価 8）を目的としたものがある．

工皮膚では，皮膚内部の光伝播が十分に考慮されていない

また，医療分野では，
レーザー治療や画像診断技術の研究 9），からである．色味は似ていても，実際の皮膚の質感と異な

超音波診断や針・メスなどの医療器具の練習用 10）としてっていることが一つの理由である．

人工皮膚の開発が行われてきた．これらの人工皮膚では，人の皮膚の見え方については，Computer Graphics（CG）
分野を中心に盛んに研究が進められている

4 - 7）

人の皮膚の硬さや動摩擦係数など機械的性質の再現に関す．従来の研る研究が多く，光学特性に着目したものでも皮膚の見た目究として，皮膚の見え方を再現する手法の一つとして，二の再現が目的ではなく，使用するレーザーや x 線など限定

方 向 反 射 率 分 布 関 数（BRDF：Bidirectional Reflectance

的な波長での光学特性の再現とされている．一方，人の皮Distribution Function）を CG に適用したモデルが考えられ

膚の見た目を再現することを目的とした人工皮膚の作製にてきた

4 , 5）．しかし，この技術を用いて人の皮膚を再現し関する研究は少ない．Iva ら 11）は，人工皮膚で人の皮膚の分光反射率を再現することを目的に研究を行った．その結果，シリコーン樹脂に馬の血液と脂肪乳剤の主成分の一つ

Development of artificial skin for optical
performance evaluation of makeup
cosmetics
Jun Yamada
Department of Engineering, Shibaura
Institute of Technology

であるイントラリピッドを用いて，人の皮膚に近い分光反射率を持った独自の 9 層の人工皮膚を開発した．また，今後は，真皮層及び皮下組織も同様に色再現を行っていく必
要があることを結論として述べている．前述したように，
人の皮膚の見え方（質感）を再現した人工皮膚の作製には，− 75 −

コスメトロジー研究報告 Vol.27, 2019

分光反射率の調整だけでなく，皮膚─空気界面での光反射

具体的な手法として，透明な樹脂を母材とし，吸収係数

性質と皮膚内部での散乱及び吸収を再現することが不可欠

を染料で，散乱係数及び散乱位相関数を可視光において吸である．

収のない白色粒子で，独立に制御することにより光物性の

筆者ら

12）

は皮膚表面の構造が光反射性質に与える影響を再現を試みる．

調べ，皮膚の持つ皮溝と皮丘にある細かい凹凸の傾きの方向分布を決定することで，人の皮膚の空気─皮膚界面での光反射性質を再現できることを示した．また，筆者ら

2 . 2 . 空気 – 皮膚界面での光反射性質の制御方法

13 , 14）は，筆者らは，共焦点レーザー顕微鏡によって取得した皮膚

非侵襲で生きた人の皮膚の光物性を計測する手法及び計測

表面構造の形状データ（Fig. 1）を用いて，空気─皮膚界面装置の開発を行い，多数の日本人の皮膚の光物性の計測をでの光反射性質を予測できる解析モデルを構築してきた 12）．

行ってきた．本研究では，これらの計測データに基づき，

図 1 ⒜のような，皮膚表面形状データは，一定間隔で記

人工皮膚において人の皮膚の質感を再現するために，空気録された高さデータ群から構成される．この形状データの，─皮膚界面での光反射性質及び皮膚の光物性を再現する手

隣り合うデータを結んだ三角形の表面要素（図 1 ⒝）を考え，法の確立を目的とする．

その法線ベクトルの傾斜角度を算出する．そして，傾斜角2. 方度に関する確率密度関数としてまとめる．この確率密度関法

数は，M を表面要素の数，Mtot を一つの形状データに含ま

2 . 1. 光物性の制御方法

れる，全表面要素数，ΔΩe を表面要素の法線ベクトルが向人の皮膚では，外から皮膚に入射した光のうち，一部はく θe，φe 方向まわりの立体角として，（2）式で表される．

皮膚表面で反射されるが，残りは皮膚内部に浸透する．内……（2）部に浸透した光には，散乱されて再び皮膚の外に出ていく
もの，あるいは，吸収されて熱に変わるものがある．このこれによって算出された二方向反射率は，実験から得らことから，皮膚は散乱・吸収性を持つ連続媒体と考え，次れた皮膚表面の光反射性質と一致することが分かった．つの光輸送方程式により取り扱われることが多い

15）．

まり，表面構造の光反射性質は，その表面要素の傾斜角度…（1）2

ここで，I（W/m ･sr）は光の強度，Ω（−）は光の進行方向を表す単位ベクトル，s（m）はその方向の距離である．

に関する確率密度関数を算出することによって評価できる．
本研究では，人工皮膚作の表面形状データ（表面要素の傾斜角度に関する確率密度関数）が，実際の皮膚と同等となる人工皮膚の製作を目指す．

こ の 輸 送 方 程 式 に は， 吸 収 係 数 α（mm- 1）， 散 乱 係 数 σs

一般的に，人工皮膚の作製では，皮膚色に調整したシリ（mm- 1），散乱位相関数 p（-）といった光物性が現れる．こコーン樹脂やウレタン樹脂を形成容器に流し込み，硬化されら 3 つの光物性が皮膚内部での光伝播を特徴付ける．

せて作られる．よって，目的とする人工皮膚を実現するた

これまでの人工皮膚の製作では，上記の光物性の再現とめには，空気─皮膚界面での光反射性質を再現できる構造いう視点がなかったことに加えて，製作された人工皮膚が，を付与するための鋳型を作製する必要がある．本研究では，

実際の皮膚と同等であることを，光物性レベルで評価され鋳型となる金属表面に様々な加工方法を用いて細かい凹凸ることはなかった．本研究では，散乱係数，吸収係数の可を作製し，どのような加工方法が皮膚の表面構造の再現に視光域におけるスペクトルが実際の皮膚と同等となる人工適しているかを検証した．

皮膚の製作を目指す．

Fig. 1 （a）Surface morphology data of human skin and（b）
analytical model of surface structure.

− 76 −

メーキャップ化粧品の光学的機能評価に向けた人工皮膚の開発

3 . 結果と考察
3 . 1. 人工皮膚における光物性の再現
3 . 1 . 1 . 吸収係数の線形性
透明樹脂に添加する染料の質量濃度を変化させたサンプルを作製し，透過率 τ を分光光度計（株式会社日立
ハイテクサイエンス，U- 5100）を用いて測定し，以下のLambert-Beer の法則に基づき，染料の吸収係数を計測した．
……（3）
式中の sd（m）は光路長である．ここで，染料は母材の透明樹脂に溶解し，散乱を生じないものを用いる．そうする
ことで，減衰係数 β（mm- 1）のうちの散乱係数 σs が 0 となり，透過率から吸収係数が測定可能となる．本研究では，染料と透明樹脂との相性を考慮し，透明なウレタン樹脂（株式Fig. 2

Absorption coefficient of green dye samples.

Fig. 3

Reproduction of skin absorption coefficient.

会社ポリシス，PC- 00）と，日本化薬株式会社製；Yellow
GN，Red B， 有 本 工 業 株 式 会 社 製；Plast Blue，Plast
Green，Plast Orange の 5 種類の染料を使用した．
Fig. 2 に染料（Plast Green）の質量濃度と吸収係数のグラフを示す．Fig. 2 の結果から，質量濃度と吸収係数が線形性を持つことが示された．同様の計測を行った結果，全ての染料において，質量濃度と吸収係数の線形性が示された．
3 . 1 . 2 . 染料の足し合わせによる吸収係数の制御
3 . 1 . 1 . 項の結果から，染料の質量濃度を調整すること
で吸収係数を制御可能であることが確認された．そこで，前述した 5 種類の染料をウレタン樹脂に混ぜ，人の皮膚の吸収係数の再現する人工皮膚を作製し，そのときの光物性を，本研究室で開発した光物性計測装置 13）で計測した（Fig.
3）．この際，Fig. 2 の濃度に対する吸収係数の傾き αλ を全ての染料について算出し，非負の最小二乗法 16）を用いて，必要とされる染料の配合する質量濃度を求めた．その結果，人工皮膚の吸収係数は，予め計算した染料の足し合わせで得られることが確認できた．吸収係数の制御が本手法で可能である．
3 . 1 . 3 . Mie 理論による散乱係数の計算
媒質中の散乱係数は球形の粒子を用いた場合，Mie 理論
から導かれる粒子一つの散乱断面積と，粒子数密度を与えることで求められる．これまでの筆者らの研究 17）で，粒子の粒径分布を考慮した Mie 理論を適用することで，粒子の平均直径 DN（µm）と体積分率 fv（-）から，Fig. 4 に示す
ように，人の皮膚の散乱係数スペクトルに近い値を再現できると予測された．しかし，この結果を基に製作した人工
皮膚の散乱係数スペクトルは，大きさ，プロファイルとも
予測値とは異なった．この原因は，散乱位相関数の波長依存性にある．
3 . 1 . 4 . 散乱位相関数の影響を考慮した散乱係数の制御
皮膚の光物性を計測した手法（RSPM）13） では，波長− 77 −

Fig. 4

Calculated scattering coefficient by Mie theory.

コスメトロジー研究報告 Vol.27, 2019

λ= 633 nm の培養皮膚の散乱位相関数 7）を与えて，逆解析
を通じて散乱係数，吸収係数を推定している．しかし，実際の散乱位相関数には波長依存性があり，Mie 理論においても波長の関数として散乱位相関数が導かれる．先に述べた，人工皮膚で散乱係数を再現できなかった原因は，Mie理論での予測と計測における逆解析で，散乱位相関数が異Fig. 5なっていたためである．

Simulation flow.

そこで，Fig. 5 に示すモデルを用いて，Mie 理論から算出された散乱係数（σsa），吸収係数及び散乱位相関数を皮膚の物性値として割り当て，培養皮膚の散乱位相関数を用いて逆解析により推定される値（σse）を数値シミュレーションにより求めた（Fig. 6）．Fig. 6 の結果から，散乱位相関数が異なる場合，σsa と σse に大きな差が生じることが分かる．
また，このシミュレーションから，粒子径が 0 . 1 - 0 . 2 µmの Al 2 O 3 粒子を用いれば，σse が人の皮膚について計測された散乱係数スペクトルを再現できると考えられる．
本解析の検証として，DN= 0 . 6 µm の Al 2 O 3 粒子（住友化学株式会社，アドバンストアルミナ AA- 05）を使用したサンプルを作製し，そのときの光物性を RSPM で計測した（Fig. 7）
．Fig. 7 の結果より，シミュレーションから予想される結果と人工皮膚の計測結果の値と傾きはよく一致し，この手法から散乱係数・散乱位相関数の制御が可能である
ことが示された．

Fig. 6

Simulation of influence of scattering phase function.

Fig. 7Comparison of simulated result and artificial skin.

3 . 2 . 人工皮膚における空気- 皮膚界面での光反射性質の再現
3 . 2 . 1 . 鋳型の製作方法
鋳型となる形成容器の一面に機械加工を施し，その加工法に応じて，どのような光性質を持つ表面が作られるか検討する．Table 1 に，今回行った加工法を示す．鋳型は，アルミニウム合金で作製することを想定し，大きさ 50 ×
50 × 3 mm の A 5052 基板を対象としてその表面に加工実験を行う．
3 . 2 . 2 . 試作した鋳型の評価及び光性質の予測
Fig. 8 に，各加工法で処理された表面の拡大画像，Fig.
9 に，各表面構造のプロファイルの一例を示す．加工法 1（切削）では，約 15 µm 程度のピーク形状が見られる．加工法
2（サンドブラスト加工）では，不規則で細かい凹凸が形

Table 1 Machining to construct the surface structures on A 5052 substrates.
加工法 1：切削（Cutting）

カッターナイフによってランダムな方向に溝状の傷を付与する

加工法 2：
砥粒径 φ = 425‒600µm のアルミ系サンドブラスト用砂を，噴射ノズルから加工する A5052
サンドブラスト加工（Sandblasting） 基板表面までの距離を約 50mm に固定した状態で，噴射圧 0.3 MPa，8 分間吹き付ける．
加工法 3： 揺動（Shaking）

大きさ W73 × D73 × H80 の容器底部に両面テープで A5052 基板を貼付け，その中に粒
径 φ = 1.6mm の 球 状のステンレス投 射 材を入れて，手 動で上下に揺 動させることによってA5052 表面に凹凸構造を付与する．上下の揺動幅は約 110mm，振動数は 1 分間に約 200 回とした．

− 78 −

メーキャップ化粧品の光学的機能評価に向けた人工皮膚の開発

Fig. 8 Surface image of A5062 substrates after machining and humanskin surface.

成されているが，その深さは約 6 - 12 µm 程度と見積もられる．加工法 3（揺動）は，サンドブラスト加工と類似して，不規則で細かい凹凸が見られるが，その構造は，サンドブ
ラスト加工と比較してかなり小さく，ほとんど平滑であった．これは，手加工の方が，投射材の速度が遅かったためだと考えられる．
Fig. 10 に，各加工法によって形成された表面形状デー
タから算出した，傾斜角度に関する確率密度関数 psurf（θe ,φe）を示す．この図より，実際の皮膚の表面構造の psurf（θe ,φe）は，傾斜角 θe = 45 ˚ 付近でゆるやかに上昇する傾向が見られた．これは，皮溝の斜面に起因すると考えられる．同様に，切削加工の psurf（θe , φe）で θe = 75 ˚ 付近に小さなピ

Fig. 9 Surface profile of A5052 substrates after machining
and human skin.

ークが現れるのも，溝構造の斜面によるといえる．そして，前述の確率密度関数をもとに，各表面に，光を入射角
θin = 45 ˚ で入射させたときの光反射性質を予測した結果を，Fig. 11 に示す．全ての加工法において，皮膚表面の光反射性質と異なる結果となった．特に，揺動加工で作製した表面構造の光反射性質は，Fig. 11 に示されるように，正反射方向付近の二方向反射率が急激に上昇した．Fig. 10の結果より，これは，表面がほとんど傾斜のない構造であったことが原因だと考えられる．
これらの結果から，投射材形状や加工法により，表面構
造の特徴を制御できることが示された．よって，これらの適切な組み合わせが明らかになれば，実際の皮膚表面構造
を人工皮膚に転写せずとも，皮膚表面構造による光性質を再現する可能性があることが示された．

Fig. 10 Probability density function for machining surface
structures and human skin.

3 . 2 . 3 . 表面構造の違いが内部での光伝播に与える影響
次に，皮膚内部での光伝播も再現した人工皮膚作製に向け，Fig. 10 に示すような，表面傾斜角度の確率密度関数の違いが，皮膚内部での光伝播にどの程度影響を及ぼすかを，数値解析によって検討した．この光伝播解析は，Fig.
12 のような解析モデルに対して，Monte Carlo 法によって行う．これは，光を，エネルギーを持つ光束とし，解析モデルに入射させる試行を繰り返すことで近似的に解を求める解析方法である．皮膚表面での光挙動は，前述のように，表面形状データから算出される，傾斜角度の確率密度関数 psurf（θe , φe）よりに決定した．光の入射角度は θin = 45 ˚とし，光束が皮膚 - 空気界面に到達したとき，Fig. 5 の確

Fig. 11 Bidirectional reflectance for machining surface
structures and human skin obtained by numerical analysis.

− 79 −

コスメトロジー研究報告 Vol.27, 2019

率密度関数をもとに，これが出くわす表面要素，及び，その表面要素上で反射するか，あるいは透過して皮膚内部に進むかが決定される．解析モデルの大きさは，一般的な皮膚の厚みとされる 1 . 9 mm とし，その周囲は周期境界とした．そして，光束が，皮膚内部で吸収されるか，あるいは再び皮膚－空気界面に達して皮膚外に射出されるまで追跡する．
表面構造による光性質の違いが，皮膚内部での光伝播に与える影響を調べるために，実際の皮膚を含む各表面構造について，表面構造を透過した直後の二方向透過率を Fig.
13 に，皮膚内部で散乱されて，再び皮膚外に射出された光が到達した最大の深さと光量の関係をFig. 14 にまとめた．

Fig. 12 Analytical model to analyze the influence of thestructure of the skin surface on light propagation inside theskin.

Fig. 14 の結果より，光がより深くまで到達したのは，揺動加工によるものであった．この加工法によって作られた表面構造は，Fig. 9，10 の結果より，ほとんど傾斜がない．
そのため，光がこの構造表面を透過しても，その指向性が大きく変化せず，その結果として，深さ方向に進む光が減
少しなかったためである．逆に，皮膚表面を透過した光は，
その指向性が低下したため，皮膚の深くまで進む前に皮膚から射出されたと考えられる．これらの結果から，皮膚表面構造は，皮膚内部の光伝播にも影響を与えることが明らかとなった．

4. 総括

本研究では，生きた人の皮膚内部での光伝播及び空気 ‒Fig. 13 Light transmittance for each direction on thesurface structures obtained by the numerical analysis.

皮膚界面での光反射性質を再現した人工皮膚の作製を目的とし，吸収係数や散乱係数・散乱位相関数の 3 つの光物性を制御する方法及び表面構造を転写する鋳型の加工方法を検討した．
吸収係数については，人工皮膚の母材となる透明樹脂に溶解する染料を用いた場合，質量濃度と吸収係数が線形性を持ち，複数の染料を組み合わせることで，吸収係数を制御可能であることが示された．散乱係数については，散乱位相関数が Mie 理論と計測における逆解析とで異なっているため，その違いを考慮した新たな解析モデルを導入した．
その結果，この解析モデルから計算した値と，作製したサンプルの光物性の計測値を比較した結果，散乱係数スペク
トルの一致が確認された．これらのことから，吸収係数をFig. 14 Influence of the surface structure on the lightpropagation inside the skin.

染料で，散乱係数・散乱位相関数を白色粒子で，独立に制御することにより，人の皮膚の光物性を再現した人工皮膚が作製可能であることが示された．

本研究では，人工皮膚の光学的な性質を生きた人の皮膚

また，様々な加工方法を用いて，金属盤上に様々な表面と同等のものにするための手法を，内部と表面を個別に分構造を作製し，その構造が生じさせる光反射性質を分析しけて検証した．ここでは，作製した人工皮膚の質感の評価た結果，一種類の加工方法では人の皮膚の光反射性質を再までは行えておらず，その作製方法の検証のみであったが，現することができないが，複数の加工方法を組み合わせる

これらの内部と表面の製作方法を合わせて一つの人工皮膚ことで，人の皮膚の光反射性質を再現可能であることが示を作製することで，より生きた人の皮膚の質感を再現したされた．

物を作製可能になると考えられる．また，本研究では，光− 80 −

メーキャップ化粧品の光学的機能評価に向けた人工皮膚の開発

学的な性質を個別に制御する手法を提案しているため，本8） 菅沼薫 , 丹羽雅子：マッサージクリームの触感と物手法を利用することで，個人差が大きく，人によって異な理物性との関係 . 日本化粧品技術者会誌 , 24, 212-219,

る皮膚の光学的な性質に合わせて，化粧品の評価に使用す1991.
9） Jędrzejewska-Szczerska M, Wróbel MS, Galla S, Popov

るための様々な種類の人工皮膚の製作できる．

AP, Bykov AV, Tuchin VV, Cenian A. Investigation
謝辞of photothermolysis therapy of human skin diseases

本報告書の作成にあたり，共同研究者である河野貴裕氏

using optical phantoms. 1 8th Int. School on Quantum（芝浦工業大学）に多大な貢献をいただいた．ここに記してElectron.: Laser Phys. Applications, Int. Soc. Opt.
Photonics, 9447, 2015.

謝意を表す．

10）Aminzadeh R, Saviz M, Shishegar AA. Theoretical（引用文献）and experimental broadband tissue-equivalentphantoms at microwave and millimetre-wave

1） N i s h i m u r a H , T a k a s u k a Y , Y a m a m o t o M .

frequencies. Electron. Lett., 50 , 618 ‒ 620 , 2014 .

Optical properties of skin gloss and developmentof mizumizushii-looking makeup foundation. Int. J.

11）Iva S, Tanabe A, Maeda T, Funamizu H, Aizu Y.
Development of non-deterioration-type skin tissue

Cosmet. Sci., 29, 67, 2007.

phantom using silicone material. Opt. Rev., 2 1 ,2） Nishikata, N., Nishimura, H., Mohri, K., Nakamura, N.

353 ‒ 358 . 2014 .

A Natural-Looking Makeup. Cosmetics & Toiletries12）Nakamura K, Ogawa T, Takata S, Yamada J.

Magazine, 112, 39-56, 1997.
3） 上原孝一 , 南浩治 , 岩本啓 , 長田みゆき , 五十嵐祟訓 ,Numerical model for optical characteristics of the

中尾啓輔 , 大崎和友 , 植松隆史 , 野尻尚材 , 樫本明生 , 福human skin surface. Thermal Sci. Eng., 2 3 , 6 1 ‒ 6 9 ,田啓一：透過光制御による若顔印象ファンデーションの2015 .
13）Kono T, Yamada J. Development of measurement

開発 . 日本化粧品技術者会誌 , 44, 48-56, 2010.

instrument for optical properties of human skin in4） Weyrich T, Matusik W, Pfister H, Bickel B,

vivo. Jap. J. Thermophys. Prop., 31 , 72 ‒ 80 , 2017 .

Donner C, Tu C, McAndless J, Lee J, Ngan A, Jensen
HW, Gross M. Analysis of human faces using a14）Kono T, Toma S, Nozawa K, Yamada J. Scattering

measurement-based skin reflectance model. ACM

and Absorption Coefficients of Human Skin in JapaneseTrans. Graphics, 25, 1013-1024, 2006.

People. Proc. 1 st Pacific Rim Thermal Eng. Conf.,
Hawaii’ s Big Island, 15081 - 1 ‒ 5 , 2016 .

5） Marschner SR, Westin SH, Lafortune EPF, Torrance
KE, Greenberg DP. Image-based BRDF measurement15）Siegel R, Howel JR. Thermal Radiation Heat

including human skin. In Rendering Techniques’ 9 9,Transfer, 3 rd edn., Taylor and Francis, New York,

Springer, Vienna, 131-144, 1999.

686 ‒ 689 , 1992 .

6） Alotaibi S, Smith WAP. A Biophysical 3D Morphable16）Lawson CL, Hanson RJ. Solving Least Squares

Model of Face Appearance. Proc. 2 0 1 7 IEEE Conf.

Problems, Upper Saddle River, NJ: Prentice Hall, 161 ,1974 .

Comput. Vision and Pattern Recogn., 824-832, 2017.
7） Jensen HW, Marschner SR, Levoy M, Hanrahan P.

17）Sarasuk K, Kono T, Nakamura K, Yamada J. EffectA practical model for subsurface light transport. Proc.

of scattering particle size on extinction coefficient of

2 8th Annu. Conf. Comput. Gr. Interact. Techniques,

artificial human skin. Proc. 1 st Pacific Rim Thermal

511-518, 2001.

Eng. Conf., Hawaii’ s Big Island, 15111 - 1 ‒ 5 , 2016 .

− 81 −