2.2 2(a) 2(b) 1 Torrance-Sparrow 2.3 2(c) (a) 2 (b) (c) RGB 0.07mm0.2mm c 2009 Information Processing Society of Japan

Similar documents
(a) 1 (b) 3. Gilbert Pernicka[2] Treibitz Schechner[3] Narasimhan [4] Kim [5] Nayar [6] [7][8][9] 2. X X X [10] [11] L L t L s L = L t + L s

IPSJ-CVIM

Ω + θr θi H α N β = + φ φ r x φ i 2 diffuse component specular component (a) 1 (b) BRDF (R) (G) (B) BRDF f BRDF (x, θ i, φ i, θ r, φ r ). (1) 1(b) N φ

情報処理学会研究報告 方向の組み合わせを変えながら反射光の強度を計測する この方法は単純であるが入射方向と反射方向の組み合わせ を同時に 1 つしか計測できないため すべての方向を計 測するには膨大な計測回数が必要となる問題がある その ため 計測を効率的に行う方法がこれまでに提案されてい る [3

Vol. 44 No. SIG 9(CVIM 7) ) 2) 1) 1 2) 3 7) 1) 2) 3 3) 4) 5) (a) (d) (g) (b) (e) (h) No Convergence? End (f) (c) Yes * ** * ** 1

20 Method for Recognizing Expression Considering Fuzzy Based on Optical Flow

2. Hašan [2] Dong [3] 3D Papas [4] Brunton [5] 3D [6] UV 3. UV LUT: Lookup Table LUT LUT 計測系 出力系 半透明材質の計測 UV インクの計測 要求材質の設計 ルックアップテーブル逆引き 印刷 2 印刷時の半透明

IPSJ SIG Technical Report Vol.2015-CVIM-196 No /3/6 1,a) 1,b) 1,c) U,,,, The Camera Position Alignment on a Gimbal Head for Fixed Viewpoint Swi

0 21 カラー反射率 slope aspect 図 2.9: 復元結果例 2.4 画像生成技術としての計算フォトグラフィ 3 次元情報を復元することにより, 画像生成 ( レンダリング ) に応用することが可能である. 近年, コンピュータにより, カメラで直接得られない画像を生成する技術分野が生

untitled

IPSJ SIG Technical Report Vol.2018-CVIM-210 No /1/18 1,a) [29] [14] [2] [25], [30] [3] [21] [22] a) ikeya.nobuhir

Vol.2018-CG-172 No.17 Vol.2018-DCC-20 No.17 Vol.2018-CVIM-214 No /11/8 情報処理学会研究報告 IPSJ SIG Technical Report プロジェクタ カメラシステムの同期遅延と 走査速度の制御にもとづく光伝

2 1 ( ) 2 ( ) i

3 2 2 (1) (2) (3) (4) 4 4 AdaBoost 2. [11] Onishi&Yoda [8] Iwashita&Stoica [5] 4 [3] 3. 3 (1) (2) (3)

1 UD Fig. 1 Concept of UD tourist information system. 1 ()KDDI UD 7) ) UD c 2010 Information Processing S

(a) (b) 2 2 (Bosch, IR Illuminator 850 nm, UFLED30-8BD) ( 7[m] 6[m]) 3 (PointGrey Research Inc.Grasshopper2 M/C) Hz (a) (b

IPSJ SIG Technical Report Vol.2013-CVIM-187 No /5/30 1,a) 1,b) 2,c) 1,d) 1. [1][2] 1 [3] 1 Presently with Osaka Uniersity 2 Presently with Micr

Vol1-CVIM-172 No.7 21/5/ Shan 1) 2 2)3) Yuan 4) Ancuti 5) Agrawal 6) 2.4 Ben-Ezra 7)8) Raskar 9) Image domain Blur image l PSF b / = F(

1 Table 1: Identification by color of voxel Voxel Mode of expression Nothing Other 1 Orange 2 Blue 3 Yellow 4 SSL Humanoid SSL-Vision 3 3 [, 21] 8 325

情報処理学会研究報告 IPSJ SIG Technical Report Vol.2013-CVIM-188 No /9/3 BRDF i

光学

IPSJ SIG Technical Report An Evaluation Method for the Degree of Strain of an Action Scene Mao Kuroda, 1 Takeshi Takai 1 and Takashi Matsuyama 1

す 局所領域 ωk において 線形変換に用いる係数 (ak 画素の係数 (ak bk ) を算出し 入力画像の信号成分を bk ) は次式のコスト関数 E を最小化するように最適化 有さない画素に対して 式 (2) より画素値を算出する される これにより 低解像度な画像から補間によるアップサ E(

3 3 3 Knecht (2-3fps) AR [3] 2. 2 Debevec High Dynamic Range( HDR) [4] HDR Derek [5] 2. 3 [6] 3. [6] x E(x) E(x) = 2π π 2 V (x, θ i, ϕ i )L(θ

Gaze Head Eye (a) deg (b) 45 deg (c) 9 deg 1: - 1(b) - [5], [6] [7] Stahl [8], [9] Fang [1], [11] Itti [12] Itti [13] [7] Fang [1],

IPSJ-CVIM

(MIRU2008) HOG Histograms of Oriented Gradients (HOG)

14 2 5


258 5) GPS 1 GPS 6) GPS DP 7) 8) 10) GPS GPS ) GPS Global Positioning System

,,.,.,,.,.,.,.,,.,..,,,, i

<95DB8C9288E397C389C88A E696E6462>

1(a) (b),(c) - [5], [6] Itti [12] [13] gaze eyeball head 2: [time] [7] Stahl [8], [9] Fang [1], [11] 3 -

IPSJ SIG Technical Report Vol.2012-CVIM-182 No /5/ RGB [1], [2], [3], [4], [5] [6], [7], [8], [9] 1 (MSFA: Multi-Spectrum Filt

Silhouette on Image Object Silhouette on Images Object 1 Fig. 1 Visual cone Fig. 2 2 Volume intersection method Fig. 3 3 Background subtraction Fig. 4

1 Kinect for Windows M = [X Y Z] T M = [X Y Z ] T f (u,v) w 3.2 [11] [7] u = f X +u Z 0 δ u (X,Y,Z ) (5) v = f Y Z +v 0 δ v (X,Y,Z ) (6) w = Z +

IPSJ SIG Technical Report Vol.2012-HCI-149 No /7/20 1 1,2 1 (HMD: Head Mounted Display) HMD HMD,,,, An Information Presentation Method for Weara

2. ICA ICA () (Blind Source Separation BBS) 2) Fig. 1 Model of Optical Topography. ( ) ICA 2.2 ICA ICA 3) n 1 1 x 1 (t) 2 x 2 (t) n x(t) 1 x(t

No.3 14

xx/xx Vol. Jxx A No. xx 1 Fig. 1 PAL(Panoramic Annular Lens) PAL(Panoramic Annular Lens) PAL (2) PAL PAL 2 PAL 3 2 PAL 1 PAL 3 PAL PAL 2. 1 PAL

& Vol.5 No (Oct. 2015) TV 1,2,a) , Augmented TV TV AR Augmented Reality 3DCG TV Estimation of TV Screen Position and Ro

3D VR CAD 3D CAD CAD [1] CAD 3DCG [2] [3] CAD 3D NC CG [4] Ccurve XY C curve α C curve [5], [6], [7], [8], [9] 2 [10] 1 [11], [12] 2.2 [13] Tcu

75 Author s Address: Possibility of Spatial Frequency Analysis of the Three-dimensional Appearance and Texture of Facial Skin

IPSJ SIG Technical Report Vol.2012-CG-148 No /8/29 3DCG 1,a) On rigid body animation taking into account the 3D computer graphics came

,,,,,,,,,,,,,,,,,,, 976%, i

光学

dsample.dvi

Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels).

1 Fig. 1 Extraction of motion,.,,, 4,,, 3., 1, 2. 2.,. CHLAC,. 2.1,. (256 ).,., CHLAC. CHLAC, HLAC. 2.3 (HLAC ) r,.,. HLAC. N. 2 HLAC Fig. 2

PowerPoint プレゼンテーション

untitled

OJT Planned Happenstance

IPSJ SIG Technical Report Vol.2010-CVIM-170 No /1/ Visual Recognition of Wire Harnesses for Automated Wiring Masaki Yoneda, 1 Ta

Miyazaki-3DForum dvi

Web Web [4] Web Web [5] Web 2 Web 3 4 Web Web 2.1 Web Web Web Web Web 2.2 Web Web Web *1 Web * 2*3 Web 3. [6] [7] [8] 4. Web 4.1 Web Web *1 Ama

<332D985F95B62D8FAC93638BA795DB90E690B62E706466>

塗装深み感の要因解析

IPSJ SIG Technical Report 1, Instrument Separation in Reverberant Environments Using Crystal Microphone Arrays Nobutaka ITO, 1, 2 Yu KITANO, 1

IPSJ SIG Technical Report Vol.2011-EC-19 No /3/ ,.,., Peg-Scope Viewer,,.,,,,. Utilization of Watching Logs for Support of Multi-

1214_KiyotaCalib_matsusita_fixed2.pdf

(3.6 ) (4.6 ) 2. [3], [6], [12] [7] [2], [5], [11] [14] [9] [8] [10] (1) Voodoo 3 : 3 Voodoo[1] 3 ( 3D ) (2) : Voodoo 3D (3) : 3D (Welc

2). 3) 4) 1.2 NICTNICT DCRA Dihedral Corner Reflector micro-arraysdcra DCRA DCRA DCRA 3D DCRA PC USB PC PC ON / OFF Velleman K8055 K8055 K8055

Abstract This paper concerns with a method of dynamic image cognition. Our image cognition method has two distinguished features. One is that the imag

2. CABAC CABAC CABAC 1 1 CABAC Figure 1 Overview of CABAC 2 DCT 2 0/ /1 CABAC [3] 3. 2 値化部 コンテキスト計算部 2 値算術符号化部 CABAC CABAC

1: A/B/C/D Fig. 1 Modeling Based on Difference in Agitation Method artisoc[7] A D 2017 Information Processing

IPSJ SIG Technical Report Pitman-Yor 1 1 Pitman-Yor n-gram A proposal of the melody generation method using hierarchical pitman-yor language model Aki

(4) ω t(x) = 1 ω min Ω ( (I C (y))) min 0 < ω < C A C = 1 (5) ω (5) t transmission map tmap 1 4(a) t 4(a) t tmap RGB 2 (a) RGB (A), (B), (C)

Table 1. Reluctance equalization design. Fig. 2. Voltage vector of LSynRM. Fig. 4. Analytical model. Table 2. Specifications of analytical models. Fig

untitled

IPSJ SIG Technical Report Vol.2012-CVIM-183 No /9/2 1,a) 2,b) Stephen Lin 2,c) 1,d) 1,e) 1. Woodham [25] Silver [20] I l n I(x) = ρ(x)n T l (1)

1 Web [2] Web [3] [4] [5], [6] [7] [8] S.W. [9] 3. MeetingShelf Web MeetingShelf MeetingShelf (1) (2) (3) (4) (5) Web MeetingShelf

IPSJ SIG Technical Report Vol.2009-CVIM-169 No /11/ Stereo by the horizontal rotary movement of the upswing fisheye camera Sat

芸術科学会論文誌 Vol. 14, No. 6, pp (2015) 1 はじめに 1 枚の人物写真に対して半透明化処理を行う場合 人 物領域を切り抜いて空となった背景部分の補間を行うこ とが考えられるが 画像 被写体人物の遮蔽によって背景部分の情報は画像中 には存在しない 補間す


第62巻 第1号 平成24年4月/石こうを用いた木材ペレット

2003/3 Vol. J86 D II No Fig. 1 An exterior view of eye scanner. CCD [7] CCD PC USB PC PC USB RS-232C PC

1

浜松医科大学紀要

Microsoft Word - toyoshima-deim2011.doc



IPSJ SIG Technical Report Vol.2013-GN-87 No /3/ Research of a surround-sound field adjustmen system based on loudspeakers arrangement Ak

Q [4] 2. [3] [5] ϵ- Q Q CO CO [4] Q Q [1] i = X ln n i + C (1) n i i n n i i i n i = n X i i C exploration exploitation [4] Q Q Q ϵ 1 ϵ 3. [3] [5] [4]

THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS TECHNICAL REPORT OF IEICE. UWB UWB

Perrett et al.,,,, Fig.,, E I, 76

IPSJ SIG Technical Report Vol.2014-CG-155 No /6/28 1,a) 1,2,3 1 3,4 CG An Interpolation Method of Different Flow Fields using Polar Inter

FUJII, M. and KOSAKA, M. 2. J J [7] Fig. 1 J Fig. 2: Motivation and Skill improvement Model of J Orchestra Fig. 1: Motivating factors for a

IPSJ SIG Technical Report Vol.2014-EIP-63 No /2/21 1,a) Wi-Fi Probe Request MAC MAC Probe Request MAC A dynamic ads control based on tra

IPSJ SIG Technical Report GPS LAN GPS LAN GPS LAN Location Identification by sphere image and hybrid sensing Takayuki Katahira, 1 Yoshio Iwai 1

MDD PBL ET 9) 2) ET ET 2.2 2), 1 2 5) MDD PBL PBL MDD MDD MDD 10) MDD Executable UML 11) Executable UML MDD Executable UML

<30315F985F95B65F90B490852E696E6464>

先端社会研究所紀要 第13号☆/2.矢崎

日本看護管理学会誌15-2

_念3)医療2009_夏.indd

Modal Phrase MP because but 2 IP Inflection Phrase IP as long as if IP 3 VP Verb Phrase VP while before [ MP MP [ IP IP [ VP VP ]]] [ MP [ IP [ VP ]]]


MA3-1 30th Fuzzy System Symposium (Kochi, September 1-3, 2014) Analysis of Comfort Given to Human by Using Sound Generation System Based on Netowork o

SERPWatcher SERPWatcher SERP Watcher SERP Watcher,

A Study on Throw Simulation for Baseball Pitching Machine with Rollers and Its Optimization Shinobu SAKAI*5, Yuichiro KITAGAWA, Ryo KANAI and Juhachi

[2] OCR [3], [4] [5] [6] [4], [7] [8], [9] 1 [10] Fig. 1 Current arrangement and size of ruby. 2 Fig. 2 Typography combined with printing

Transcription:

1 2 1 1 1000 2 3 2 Statistical Analysis of Human Skin Based on Reflection and Scattering Characteristics Youko Baba, 1 Tomohiro Mashita, 2 Yasuhiro Mukaigawa 1 and Yasushi Yagi 1 The characteristic of human skin depends on the age and gender. In this paper, we statistically analyze skin characteristics of about 1000 people in a wide age range to find the differences of skin characteristics in personality. Since skin is translucent, the incident ray on the skin surface is divided into reflection and scattering. That is, some rays reflect on the surface while the other rays scatter into the media. In our method, we first acquire reflection component and scattering component from facial images. Translucency of skin is analyzed by calculating the proportion of both components. Then, reflection component is divided into diffuse reflection and specular reflection. We analyze how oily the skin is by calculating the specular strength. Facial images among people over all ages were taken with measurement system including three projectors and two cameras. The characteristic of the facial image is calculated and the skin characteristic which depends on the age and gender is statistically analyzed. 1. ( CG) CG Debevec 1) 2) 3)4) 5) 6) 1000 2. 2.1 1 1 Institute of Scientific and Industrial Research, Osaka University 2 Cybermedia Center, Osaka University 1 c 2009 Information Processing Society of Japan

2.2 2(a) 2(b) 1 Torrance-Sparrow 2.3 2(c) (a) 2 (b) (c) RGB 0.07mm0.2mm 3. 3.1 2 2 c 2009 Information Processing Society of Japan

1 40 12 3 3 3 6 3 6 (a) (b) 3 3.2 3 3(a) RGB 10 1 25 3.3 Nayar 7) 1 3(b) 3 3 (i max ) (i min ) i max i min L r L s L r = i max i min (1) L s = 2i min (2) (1)(2) RGB RGB 20 20 3.4 9) 4 10) 4 4 3.3 1 ( 1) 1 ( 2) 2 1 1 2 2 1 2 1 1 2 3 c 2009 Information Processing Society of Japan

情報処理学会研究報告 図4 法線と鏡面反射率 源方向ベクトル 顔表面の法線ベクトル カメラの位置ベクトルの関係から鏡面反射成分が 強く観測される位置がわかる 図 4 のように光源方向の単位ベクトルを L カメラ方向の 単位ベクトルを C 法線単位ベクトルを N とする L + C の向きと N の向きが一致する 点で鏡面反射成分が最も強く観測される 本研究では L + C と N の為す角αが 15 度以下 (a) 機器構成 (b) 外観 の点を鏡面反射成分が強く観測される点とした 被験者の右頬の領域内で鏡面反射が強く観 図5 測される点の画素値を比較して最大値を求め 鏡面反射率とする 4. 計測システム 計測システム ロジェクタを異なる位置に 3 台設置することで様々な角度から光を投影することができる 4.1 機 器 構 成 これらのプロジェクタから 1 台ずつ適切なパターン光を投影して カメラで顔を撮影する 人間の肌の個性を計測するためには 幾何情報と反射特性を求める必要がある 幾何情報 パターン光を投影した画像を使って人間の顔の幾何情報と反射特性を得ることができる とは 3 次元形状と表面の法線であり 反射特性とは反射光と散乱光の強さや散乱光の空間的 4.2 キャリブレーション な広がりのことである 本研究では幾何情報と肌の反射特性についての情報をパターン光を 解析を行うにあたって幾何学的キャリブレーションと光学的キャリブレーションを行った 投影した画像から得る 幾何情報を得るためには X 方向と Y 方向のグレイコードパターン 鏡面反射成分を分離するのに必要な 3 次元座標を求めるために 空間コード化投影法に を使う 出射光を反射光と散乱光に分離するためには高周波パターンを使う また散乱光の よってキャリブレーションを行った 空間的な広がりを調べるためにはドットパターンを使う 本研究では統計的な解析のために 一方 光学的キャリブレーションでは 各プロジェクタとカメラを使って撮影した画像の 幅広い年代で多くのデータが必要であることと 被験者が静止する必要があることから撮影 色と明るさを一致させる 本研究では光源として複数のプロジェクタを用いて撮影を行う は短時間で終えなければならない よって各プロジェクタから投影するパターンと撮影する が 撮影した画像はプロジェクタの色や明るさの影響を受ける またプロジェクタやカメラ カメラは表 1 に示す必要最小限の組み合わせとした には個体差があるので 同じパターンを投影しても色と明るさが一致しない よって色と明 本研究ではこれらのカメラとプロジェクタの組み合わせを満たす図 5 のような計測シス るさはマクベスカラーチェッカを使って各 RGB 空間をアフィン変換することでマクベスカ テムを構築した 図 5(a) は計測システムの機器構成で 図 5(b) は実際の計測システムであ ラーチェッカの色空間へ値を変換した る この計測システムは箱の 1 面に顔を入れる穴が開いており その正面と右側に合計 2 台 の高感度カメラを設置している また正面と左右に光源となるプロジェクタを 3 台設置して いる カメラは Lumenera の LW-160C プロジェクタは EPSON の EMP-X5 である プ 4 c 2009 Information Processing Society of Japan

(a) (b) 6 2 R G B -0.15-0.12-0.05-0.29-0.05-0.56 0.55 0.07 0.03 5. 5.1 2009 3 3 1000 6 5.2 RGB 7 RGB 1 0 X 7 RGB 2 2 (c) 7 5.3 8 8(b) 8(c) 5 c 2009 Information Processing Society of Japan

情報処理学会研究報告 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 (a) 左頬 (a) 白パターン (b) 反射光成分 (c) 散乱光成分 㪅 㪅 㪅 (e) 鏡面反射成分 図8 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅㪌 㪅 㪅㪌 㪅 㪅㪌 (f) 唇 㪅㪌 㪅 (d) 額 (c) 鼻 㪅㪌 (d) 拡散反射成分 㪅㪌 㪅㪌 (e) あご 㪅 㪅㪌 㪅㪌 㪅㪌 㪅 㪅㪌 㪅 㪅㪌 (b) 右頬 㪅 㪅㪌 㪅 㪅 㪅㪌 㪅㪌 㪅 㪅 㪅㪌 㪅 (g) 左まぶた (h) 右まぶた 図 9 化粧による散乱光の減衰 (f) 解析部分の抽出 入射光の分離 る よって図 8(f) に示すように円形のシールを張った眼鏡を被験者にかけてもらい テン と考えられる プレートマッチングによって目の位置を抽出し その位置を基準に頬の位置を求める 5.4 化粧による散乱光の比率の変化 肌のてかりは全投影パターンでは鼻や唇に観測される てかりは反射光成分の中の鏡面反 射によって起きるので 反射光成分画像では観測されるが散乱光成分画像では観測されない 化粧は肌の反射 散乱特性を変化させる大きな要因になる 化粧をしている場合としてい ことが確認できる ない場合で 内部に入り込む光の量の変化ついて調べた 化粧をしている場合としていない また肌表面の凹凸をみると反射光成分は全体的にざらついていて 散乱光成分は滑らかで 場合の画像を撮影し 反射光成分と散乱光成分に分離し RGB ごとにその比率を調べた あることがわかる これは反射光成分が表面で反射された光の成分のみであるので表面の凹 被験者は 22 歳の女性で両頬 額 鼻 あご 下唇 両まぶたの 20 20 画素の平均の比率 凸の影響を受けるからである それに対して散乱光成分は光が肌の内部に入り込み 入射点 である 散乱光の反射光に対する比率を図 9 に示す の周辺に強い散乱光が到達するため 表面の細かい凹凸が見えにくくなるので表面が滑らか この結果より 化粧をすることで散乱光成分の比率が減ることがわかる つまり化粧をす に見える ることで肌の内部に入り込む光の量が減ると言える 部位によって比率の減少具合に差が出 また図 8(b) の画像から図 8(d) 8(e) の拡散反射成分と鏡面反射成分を分離する これら ている これは化粧をする際に肌全面に均一にムラなく化粧はせず 部位によって厚く塗る の分離した画像を使って解析を行う 解析を行うにあたり被験者の頬の位置を知る必要があ 部分と薄く塗る部分があることが原因と考えられる また赤い波長は緑と青の波長に比べ 6 c 2009 Information Processing Society of Japan

5.5 5.5.1 10 10(a) 11 11 10 10 20 5.4 10(a) RGB RGB 5.5.2 12 16bit 12 13 13 10 6. 3 1 (a) R 10 (c) B (b) G 1000 7 c 2009 Information Processing Society of Japan

11 12 13 4) K. Nitta, O. Matoba, and T. Yoshimura, Characteristics of Weight Function in a Steady-state Diffusion Optical Tomography, IEEJ Trans. FM, Vol. 127, No. 7, pp.397-401, 2007. 5) T. Weyrich, W. Matusik, H. Pfister, B. Bickel, C. Donner, C. Tu, J. McAndless, J. Lee, A. Ngan, H. W. Jensen, and M. Gross, Analysis of Human Faces using a Measurement-Based Skin Reflectance Model, Proc. SIGGRAPH2006, pp.1013-1024, 2006. 6) A. Ghosh, T. Hawkins, P. Peers, S. Frederiksen, and P. Debevec, Practical Modeling and Acquisition of Layered Facial Reflectance, Proc. SIGGRAPH Asia 2008. 7) S. K. Nayar, G. Krishnan, M. D. Grossberg, and R. Raskar, Fast Separation of Direct and Global Components of a Scene using High Frequency Illumination, Proc. SIGGRAPH2006 pp.935-944, 2006. 8) H.W.Jensen, S.R.Marschner, M.Levoy, and P.Hanrahan, A Practical Model for Subsurface Light Transport, Proc.SIGGRAPH2001,pp.511-518,2001. 9) S.K. Nayar, X. Fang and T.E. Boult, Removal of specularities using color and polarization, Proc. CVPR 93, pp.583-590, 1993. 10) S.Barsky and M.Petrou, The 4-Source Photometric Stereo Technique for Three- Dimensional Surfaces in the Presence of Highlights and Shadows, IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 25, Issue. 10,pp.1239-1252,2003. 1) P. Debevec, T. Hawkins, C. Tchou, H. Duiker, W. Sarokin and M. Sagar, Acquiring the Reflectance Field of a Human Face, Proc. SIGGRAPH 99, pp.145-156, 2000. 2) N. Tsumura, N. Ojima, K. Sato, M. Shiraishi, H. Shimizu, H. Nabeshima, S. Akazaki, K. Hori, and Y. Miyake, Image-based skin color and texture analysis/synthesis by extracting hemoglobin and melanin information in the skin, Proc. SIGGRAPH2003, pp.770-779, 2003. 3) V. Ntziachristos, A. G. Yodh, M. Schnall, and B. Chance, Concurrent MRI and diffuse optical tomography of breast after indocyanine green enhancement, Proc. National Academy of Science, vol.97, no.6, pp.2767-2772, 2000. 8 c 2009 Information Processing Society of Japan