計測自動制御学会東北支部 第 201 回研究集会 (2002.5.22)
資料番号 201-8

非侵襲小型連続血圧・心拍計測装置を用いた映像刺激の生体影響評価
Evaluation of Eﬀects of Visual Stimulation on the Human Using
Portale and Noninvasive Measuring Device for Blood Pressureand Heart Rate
○増田 達哉 ∗ ，吉澤誠 ∗∗ ，田中明 ∗ ，杉田 典大 ∗ ，阿部 健一 ∗ , 山家 智之 † , 仁田 新一 †○ Tatsuya Masuda∗ , Makoto Yoshizawa∗∗ , Akira Tanaka∗ , Norihiro Sugita∗ ,Ken-ichi Abe∗ , Tomoyuki Yambe† , Shin-ichi Nitta†*東北大学 大学院 工学研究科, **東北大学 情報シナジーセンター,† 東北大学 加齢医学研究所*Graduate School of Engineering, Tohoku University, **Information Synergy Center,Tohoku University, †Institute of Development, Aging and Cancer, Tohoku University
キーワード : 映像刺激 (visual stimulation)，生体影響 (biological eﬀect)，血圧 (blood pressure)，心拍数 (heartrate)，非侵襲 (noninvasiveness)

連絡先 : 〒 980-8579 仙台市青葉区荒巻字青葉 05 東北大学大学院 工学研究科 電気・通信工学専攻 吉澤研究室
増田 達哉，Tel.: (022)217-7074，Fax.: (022)263-9290，E-mail: masuda@abe.ecei.tohoku.ac.jp

1.

性を表す最大相互相関係数 ρmax に着目し ，これがはじめに映像刺激の生体影響評価法とし て，非侵襲的かつ近年の情報通信技術の飛躍的な発展により，コンピュータグ ラフィック (CG)によるデ ィジタル 映像等の多種多様の人工的な映像を視聴する機会が増加している．この人工的な映像刺激は，工学，医学，娯楽などに利用され ，今後のさらなる発展が期待される．しかし ，人工的な映像刺激は人体に何ら客観的・定量的評価が 可能なものとし て有望であることが 明らかとなった 2) ,3) ．し かし 従来の方法
では，ρmax を得るために大型で高価な計測装置を使用し なければ ならず，1 度に 1 人の被験者を対象とし た実験し か実行できなかった．このため従来の方法には次のような問題点がある．

かの精神的・肉体的影響を及ぼす可能性がある．実際に，ＴＶアニメーションの視聴により引き起こさ• 一人の 被験者に 対する一回の 実験が 長時間れ た発作や，人工仮想空間（ virtual environment:にわたるため被験者数の増加に限界がある．

VE ）におけ るいわゆる VE 酔いが 報告され ている• 各被験者の実験環境・実験条件を等し くする

1) ．

ことが 困難である．

これに関連し て，視覚刺激を主とし た外部刺激が 生体に 与え る影響に 関し ，著者らは ，Mayer 波
帯域における血圧変動と心拍数変動との間の関連• 一人の 被験者に 対する数人の 験者に よる実験のため ，通常の心 理状態が 保てず 映像刺激以外の要素に 対する反応が 混入する可能–1–

性がある．に ，小型の装置で計測できる心電図と光電脈波かこれ らに 対処する方法のひと つとし て ，1 回の実験をできるだけ多数の被験者に対し て同時に行うことが 考えられ る．これが 可能になれば ，被験者数の増加，実験環境・実験条件の均一化および心理的外乱の排除が 容易になると思われ る． そこで本研究では，ρmax を得るために従来使用し てきた大型で高価な計測装置の代わりに ，計測方法を変えて装置を小型化し 携帯可能にするとともに，できるだけ安価にすることによって，複数の被験ら簡単に得られ ，血圧と負の相関があるとされ ている脈波伝播時間（ pulse transmission time; PTT ）を利用する 4) ．PTTは，Fig.1のように ，心電図のR 波から 光セン サで簡単に 計測できる脈波の立ち上がりまでの時間である．PTTが 実際に血圧に反比例することが 確かめられれば ，心電図から心拍数も得られ るので ，光電脈波と心電計を計測する
ことにより，血圧と心拍数の相互相関係数を得ることができる．

者に対する同時実験ができるようにすることを目PTT的とし た．


)

まず，血圧計を使う代わりに ，安価な光電脈波(計から得られ る脈波と心電計から得られ る心電図を組合せて得られ る脈波伝播時間を血圧の代用に)できるかど うかを検討し た．被験者に筋肉負荷あ(

るいは 音負荷を行って血圧を変動させ ，このとき0の脈波伝播時間とト ノメト リ式血圧計による血圧との相関を求め，提案し た方法の妥当性を検討しFig. 1

た．次に ，光電脈波計と心電計から構成され る試作機を 8 台を完成させ ，8 人の被験者に対し 画像提示実験を行った．画像刺激によって変動を与えた心拍数・脈波伝播時間の周波数スペクトル解析を行い，複数同時実験に対する本装置の有効性を検討するとともに ，改良すべき点について明かにし た．

0.5

Time[sec]11.5

脈波伝播時間 (PTT)の定義

本研究で試作し た脈波・心電図計測装置を Fig.2に 示す．脈波とし ては ，赤外フォト・リフレ クタ(東芝製TLP907)を用いて指先の光電脈波を測定した．心電図は ，双極誘導による心電図検出回路を製作して測定した．ρmax の計算方法については 2.3
節で述べる．

2.

実験概要

2.1

実験装置



文献2) ,3) によると，0.1Hz 近傍の Mayer 波帯域における血圧と心拍数の間の規格化相互相関関数（ 相互相関係数 ）の最大値である ρmax は，映像刺激に対する人間の情動反応の変化に対応して変化する．

 


従来の方法で ρmax を非侵襲的に得るには心電計とト ノメト リ式血圧計のよ うな連続的に血圧を計測できる血圧計が 必要であった．これらの うち心電計には小型で安価なものもあるが ，連続血圧計はFig. 2大型で高価であるため，同時に実験ができる被験者の数は限られていた ．
そこで本研究では，血圧を直接計測する代わり–2–試作し た携帯型脈波・心電図計測装置の概観

2.2

実験内容

ノメト リ式血圧計のプ ローブ を装着し ている左手は極力動かさないように指示し た．

(1) 脈波伝播時間と血圧の相関性の検証実験
脈波伝播時間 PTTが実際に血圧の逆数と相関があるかど うかを検証するため，以下のような実験を行った．Fig.3に実験の概略を示す．

るためのものであり，コンピュータ・プ ログ ラム

MATLABで作られた次の式で定義される高音で不快な音を ，ヘッド ホンを通し て大きな音量で被験&$%また，音負荷とは ，被験者に精神的影響を与え者に聞かせることである．

sound[sin(2 ∗ π ∗ 500 ∗ sin(0.5 ∗ 2 ∗ π ∗ t)

 


"#$%


∗ sin(0.05 ∗ 2 ∗ π ∗ t)) ∗ sin(2 ∗ π ∗ t)]


 !

(1)

ここで，sound[]はオーディオ出力を行うためのMATLABの関数である．次のような画像提示実験

Fig. 3脈波伝播時間と血圧の相関性を検証する実験の概略を行い，心電計と光電脈波計で構成される装置（ 非侵襲小型連続血圧・心拍測定装置 ）を複数用いる

計測量は ，心電計から得られ る心電図波形，トノ メト リ式連続血圧計（ 日本コー リン JENTOW7700 ）から得られ る橈骨動脈圧波形，Fig.4に示しことにより，われわれが 提案し てきた ρmax 解析法が 可能かど うかについて検討し た．

(2) 画像提示実験
8人の被験者(20 歳〜22 歳の健康な男子学生)に対た光電脈波セン サからの光電脈波である．これらの波形は A/Dコンバータによってサンプ リング 周

し ，画像提示実験を 行った．画像はビデ オ・プ ロ期 1msで取得された．

ジェクタでスクリーンの背後から投影した．使用した画像は，Fig.6のような，人間に好感や快適な感 

情を抱かせると予想され る「ポジティブ 」と，不快感や嫌悪感を抱かせると予想され る「 ネガティブ 」の 2 種類のカテゴ リーから 成る静止画であ る．それぞれのカテゴ リーの画像は各20 枚である．Fig.7に実験プ ロト コルを示す．2 分間の「 画像提示中 」

 Fig. 4において ，各画像を 1 枚あたり 6 秒間ずつ連続し て提示した．
「 安静」では何の画像も提示し ていない．

光電脈波センサ

計測量は ，心電図波形と親指先で測った光電脈波5人の被験者(20 歳〜22 歳の健康な男子学生)に対である．この実験の場合，血圧は計測していない．

し ，Fig.5のようなプ ロトコルで，筋肉負荷と音負荷の 2 種類の負荷を与える試行を各 1 回ずつ行った．














10





Fig. 5

20(6 / )








20(6 / )筋肉負荷および 音負荷実験のプ ロト コル

Fig. 6
筋肉負荷とは ，強制的に血圧変動を作り出すためのものであり，右手で約4kgの錘を腕を曲げて支えることである．このとき，光電脈波セン サとト–3–

画像提示実験で使った静止画像の例



Fig. 7 



関関数である．

(2) 画像提示実験a) パワースペクトル

画像提示実験のプ ロト コル

まず，実験によって得られた光電脈波に基づく脈波伝播時間（ P T TP [ms] ）と心拍数（ HR [min-1] ）解析方法

2.3のデ ータを ，3 次のスプ ラ イン 補間を 用いて 時間(1) 脈波伝播時間と血圧の相関性の検証実験的に等間隔（ 60/128=0.469sec ）のデータに変換しa) 相関係数 r
解析では ，まず，各瞬時値デ ータを 心電図の R
波を 基 準とし て 拍数 単 位デ ータに 変 換し ，これを 3 次の スプ ラ イン 補間を 用いて 時間的に 等間隔（ 60/128=0.469sec ）のデータに変換した．次に，この時系列に対して周波数成分を Mayer波帯域に絞るため，帯域通過フィルタ (通過帯域 0.05Hz〜0.2Hz)を掛けた．さらに，P T TP あるいは P T TT のど ちら
かと ，ト ノ メト リ式血圧計からの橈骨動脈圧の拍内平均値 BP [mmHg]との間の相関係数 r を (2) 式のように計算し た．ただし ，比較のため，帯域通過フィルタの代わりに低域通過フィルタ(通過帯域DC〜0.2Hz)を掛け る場合も設けた．

K{x(k) − x̄}{y(k) − ȳ}
r =  k=1
KK22k=1 {x(k) − x̄}
k=1 {y(k) − ȳ}

た．次に，1ポ イントにつき前後 60 秒間，合計 120秒間のデ ータに ついて ，その 平均値を 差し 引き ，平均値からの変動成分を求めた．その変動成分について，高速フーリエ変換 FFTを用いてパワース
ペ クト ルを 求め ，0.04〜0.15Hz におけ るパワー ス
ペクトルの面積を全周波数におけるその面積で除し たものを LF 成分（ 全パワーに 占める LF 成分の割合 ），0.15〜0.45Hzにおけるパワースペクトルの面積を全周波数におけるその面積で除し たものをHF 成分（ 全パワーに 占める HF 成分の割合 ）として求めた．また，LF 成分と HF 成分の比 LF/HFを求めた．以上の解析を 1 秒毎に行い，それぞれの成分の時間変化を求めた．

b) 最大相互相関係数 ρmax
2.3(1)b)で述べた方法と同一の方法で P T TP から(2)HR までの ρmax を求めた．

ここで，x(k) と y(k) は対象とする２つの時系列のサンプル値であり，x̄ と ȳ はそれぞれの平均値，

K はデ ータ数であり，K = 1, 280（ 10 分間 ）であ

3.

結果(1) 脈波伝播時間と血圧の相関性の検証実験

る．

Fig.8は，筋肉負荷実験から得られたある被験者b) 最大相互相関係数 ρmax
相関係数 r と同様に ，対象とする 2 種類の信号にのフィルタ処理前の生の各時系列波形の例である．

対して 0.1Hz 前後の Mayer 波に対応する成分を抽出

同図から ，筋肉負荷を行うと心拍数と 血圧が 徐々するために，0.08Hzから 0.1Hzを通過域とする帯域に増加することがわかる．また ，ト ノ メト リ式血通過フィルタを 通し た ．このような処理を行った

圧波形から得られ る脈波伝播時間 P T TT と 光電脈後の 2つの信号をそれぞれ x(k),y(k) としたとき，1波から得られ る脈波伝播時間 P T TP は互いに よく試行10 分間のデータに対し ，1 秒毎にその前後60 秒似ており，それぞれが平均血圧のゆっくりとした推分，合計 120 秒間にわたり，(3) 式のように 計算さ移と大まかに逆相関し ていることがわかる．逆相れ る規格化相互相関関数（ 相互相関係数 ）ρxy (τ )関する理由は，平均血圧の上昇が血管を硬くし ，この最大値 ρmax を求めた．

れに伴って脈波の伝播時間が 短くなるからである

φxy (τ )ρmax = max ρxy (τ ) = τφxx (0) · φyy (0)(3)

ここで，φxy (τ ) は x(t) から y(t) への相互相関関数，φxx (0)，φyy (0)はそれぞれ x(t) , y(t) の自己相と推定され る．P T TT と P T TP は互いの形が似ているが ，それぞれの絶対値のレベルは異なる（ P T TTの平均値は P T TP の平均値より 200ms 程度短い ）．
この理由は ，2 者の計測地点に 差（ 血圧は 手首で–4–

0.05Hz 0.2Hz 0.05Hz

測る血圧，光電脈波は親指で測る ）があり，かつ，P T TT が 圧力波が 伝わる血圧脈波に 基づ くもので 

あるのに 対し ，P T TP は 血流の移動に 伴う容積脈Mayer 波であるという点にあると思われ る．

Fig.9は，Fig.8と同一のデ ータについて，Mayer

波帯域に周波数成分を限定するために帯域通過フィPTT Tルタ（ 0.05Hz〜0.2Hz ）を掛け，全試行区間(10分間)

について平均値 0 分散 1に規格化し た後の時系列波形である．筋肉負荷を 行うと ，Mayer 波帯域の各PTT P









            


波形の振幅が 小さくなっていることがわかる．これは ，筋肉負荷時には筋肉を収縮させるために血流量を増大させる必要があり，血圧との関連なしに ，すなわち開ループ 状態で心拍数が 上昇するた

Fig. 9Fig.8と 同一のデ ータに 対し ，Mayer 波帯域を通過域とする帯域通過フィルタ処理をし た後の各時系
列波形

め，血圧→心拍数→血圧→という閉ループ 状態で 

生じ る Mayer 波帯域のゆらぎ が 減少し たためであ !



 #  













  

PTT T



PTT PFig. 8

" #
$ %&'())*+  !



ると思われ る．

    

Fig. 10 ある被験者の筋肉負荷実験から得られたデータ(Fig.8と同一)を低域通過フィルタ（ DC〜0.2Hz ）処理をし たものに 関する，ト ノメト リ式血圧波形から 求ま 



る脈波伝播時間 P T TT と平均血圧（ BP ）との相関



筋肉負荷時のある被験者の各時系列波形の例

伝播時間 P T TP と平均血圧の相関を 示すものである．相関係数 r は 0.831 である．

（ フィルタ処理前 ）Fig.12は，同じデータにおいて，帯域通過フィル
フィル タ処理をし ない生の 波形の まま，P T TTタ（ 0.05Hz〜0.2Hz ）を 掛け，Mayer 波帯だけに 周あるいは P T TP と平均血圧との相関を求めても相波数成分を限定し た後の，光電脈波に基づく脈波関係数 rはいずれも低い．その理由は，脈波伝播時伝播時間 P T TP と平均血圧の相関を 示すものであ間を求める際の誤差の周波数成分が 高いところにる．相関係数 r は 0.870 である．

あることと，血圧の微小な速い変動は血管の硬さ

以上は ，平均血圧との相関は ，P T TT の場合よにあまり影響を及ぼさないことにあると思われる．

り P T TP の場合の方が強く，低域通過フィルタを掛Fig.10 は，ある被験者のデ ータ（ Fig.8 と同一 ）にけた場合よりも帯域通過フィルタを 掛けて Mayer

おいて，低域通過フィルタ（ DC〜 0.2Hz ）を掛けた波帯域に周波数成分を絞った方が 相関が 強いこと

後の，ト ノ メト リ式血圧波形に基づ く脈波伝播時を示し ている．

間 P T TT と 平均血圧の相関を 示すものである．相関係数 r は 0.680 である．

この傾向は他の被験者でも同様であった．また，
この傾向は ，筋肉負荷の場合に限らず音負荷実験

Fig.11 は ，同じ デ ータに おいて ，同一の 低域通の場合も同様であった．すなわち，5 人の被験者の過フィルタを 掛けた後の，光電脈波に基づ く脈波–5– 

  


!


!


!

!






















"#$ %&'())*+  

p<0.05!!p<0.05



Fig. 11 ある被験者の筋肉負荷実験から得られたデー  

タ(Fig.8と同一)を低域通過フィルタ（ DC〜0.2Hz ）処理

Fig. 13 脈波伝達時間と血圧の相関(5人の平均，Mayerをし たものに関する，光電脈波から求まる脈波伝播時波帯域に絞った場合)

間 P T TP と平均血圧（ BP ）との相関

LF成分と HF成分は，「 ポジティブ 」画像および「ネ       
  
  
!  " #$$%&ガ ティブ 」画像提示中に 減少し ，画像提示後の安静時に上昇する傾向にあることがわかる．
一方，HF成分に対するLF成分の比であるLF/HFも画像提示後の安静時に上昇する傾向にあった．このことは，安静時におけるLF 成分の上昇の度合いが HF 成分の 上昇の 度合いより大きいと いうことを意味している．試行開始直後の安静時に LF 成分

Fig. 12 ある被験者の筋肉負荷実験から得られたデータ (Fig.8 と 同一) を 帯域通過フィル タ（ 0.05Hz〜0.2Hz ）および LF/HFの値があまり大きくない理由は，試行開始直後のため被験者の心理的な緊張が 十分に緩んでいないからであると推測され る．

処理をし たものに 関する，光電脈波から 求まる脈波伝
播時間 P T TP と平均血圧（ BP ）との相関F

相関係数と回帰直線の傾きの平均値をまとめると，

L)aFig.13 のようになる．この図は ，筋肉負荷実験とPositivePositive

P T TP の回帰直線の傾きが 高い傾向があった．

伝播時間 P T TP の うち，デ ータの 欠損がなく有効であったものは 6 人分であった．

Fig.14 左列は，この有効な 6 人について，心拍数







Negative






Positive

NegativeFig. 14



Positive



Negative
H/ FL )c 
(2) 画像提示実験

 







Positive





画像提示実験で得られた 8 人の心拍数 HRと脈波Negative
F H ) bP T TT より高く，筋肉負荷実験より音負荷実験のNegative

た ，有意ではないが ，P T TP の回帰直線の傾きがPositive
F

数が P T TT より有意に 高いことを示し ている．ま





音負荷実験のど ちらについても，P T TP の相関係

Negative   







心拍数 HR および 脈波伝播時間 P T TP のa)LF成分，b)HF成分，c)LF/HF成分の経時的変化(6 人の平均)

変動 HRのパワースペクトルの，a)LF 成分，b)HF

Fig.14右列は，Fig.14左列と同じ 6人について，光成分，c)LF/HFの平均値と標準偏差の経時的変化
を 示すもので ある 5) 6) ．それぞれ ，個人毎にデ ー

電脈波に 基づ く脈波伝播時間 P T TP のパワースペ

タの初期値を１とし て規格化している．同図から，

クトルの，a)LF成分，b)HF成分，c)LF/HFの平均値と標準偏差の経時的変化を示すものである．そ–6–

れぞれ ，個人毎にデータの初期値を１とし て規格行った．その結果，データが 有効であった 6 人の被化し ている．

験者について，映像刺激により心拍数・脈波伝播

同図から ，LF 成分，HF 成分，LF/HFは ，心拍時間の LF 成分・HF 成分がともに 減少する傾向が数の場合とほぼ 同様の経時的変化を示すが ，LF 成確認できた．し かし ，試作機を使う場合，信号に分と HF 成分については「 ネガ ティブ 」画像提示後混入するアーティファクト や雑音が 多いため，現の上昇が 遅いことがわかる．

在のところ心拍数と脈波伝播時間の最大相関係数

ただし ，被験者 1 名を除いた他の 5 名の被験者にρmax は安定的に 求められていない．

ついては，相関係数ρ(τ )を最大とする遅れ時間τ が設定範囲外となり，ρmaxが 信頼できる値とはなら試作機の信号に混入するアーティファクトや雑音とし ては次のようなものが 考えられ る．

なかった．このため，ρmax の平均値は求めていな

1) 市販の心電計に比べて試作機の心電電極の雑い．τ が 設定範囲外となる主たる原因は ，試作機

音対策が 十分でなく，特に商用周波数雑音のの性能上，心拍数や P T TP に 雑音多くが 含まれ てレ ベルが 大きい．

いたこと ，および ，特に 画像提示中には Mayer 波2) 指先へのフォト リフレ クタの固定が 適切でな帯の信号成分（ 正弦波的ゆらぎ ）が 極端に少なくく，脈波の 感度が 低く，体動の 影響を 受けやなっていたことが 考えられ る．

すい．

3) 心電アンプ・脈波アンプ のダ イナミックレン

4.

ジが 狭く，信号が 飽和する場合がある．

考察(1) 脈波伝播時間と血圧の相関性の検証実験

し かし ，上記のいずれの場合も今後装置を改良す

以上の結果を元に脈波伝播時間P T TT が血圧の代ることにより改善され る可能性は存在する．

用となるかど うかを検討する．Fig.13 より，Mayer波帯域に周波数成分を絞った場合，ト ノ メト リ式
血圧計から得られ る圧脈波に基づく脈波伝播時間

5.

P T TT と 血圧と の 相関係数は 5 人の 被験者の 平均で 0.4 程度であった．これに対し て，指先における筋肉負荷と音負荷を用いた脈波伝播時間と血圧の相関性の検証実験の結果，光電脈波から得られ る容積脈波に基づ く脈波伝播

1) Mayer 波帯域に周波数成分を絞った場合，脈波時間 P T TP と血圧との相関係数は同様に 0.6 程度で伝播時間は血圧とほぼ 反比例する．

あった．また，血圧の代わりに P T TP を使って心拍2) 指先光電脈波に 基づ く脈波伝播時間の 方が ，数までの ρmax の経時的変動を求めたとき，P T TT

ト ノ メト リ式血圧に基づ く脈波伝播時間よりを 使うものより，血圧自体を使ったものに近いこ

血圧とよく相関する．ともわかった．したがって，視覚刺激の生体影響の評価指標とし て ，Mayer 波帯域の血圧から 心拍数おわりにことがわかった．

までの相互相関係数ρmax を計算するような方法を

脈波伝播時間と心拍数は，われわれが 考案し た

採用する場合，血圧の代わりに P T TP を用いるこ

光電脈波計と心電計から成る装置で取得可能であとはある程度可能であると思われ る．し かし ，両るため，上記1)と 2)は，Mayer 波帯域に周波数成分者の相関係数があまり高くないので ，ρmaxが 必ずを限定すれば ，脈波伝播時間で代用し た血圧変動

し も確実に計算できるとは限らないことに注意すから心拍数変動までの相互相関係数の最大値ρmaxる必要があることがわかった．

が 求められ る可能性があることを意味するもので(2) 画像提示実験ある．

光電脈波計と 心電計から 構成され る試作機を 8

また，この原理に基づ く試作機 8 台を完成させ ，台を完成させ，8 人の被験者に対し 映像提示実験を8人の被験者に対する画像提示実験，すなわち，複–7–数被験者に 対する同時実験を 行 うことが で きた．
その結果，データが 有効であった 6 人の被験者について ，画像提示刺激により心拍数・脈波伝播時間の LF 成分・HF 成分がともに 減少する傾向が 確認できた．し かし ，信号に 混入するアーティファク
トや雑音が 多いため，現在のところ ρmax は安定的に求められていない．本装置の実用化には ，雑音抑制など の改良が 必要であると考えられ る．
今後の課題とし て，非侵襲小型連続血圧・心拍測定分析装置については• 精度を高め，データ収集時のエラーを少なくする．

• データ収録部分を内蔵するように改良する．
• ρmax を安定的に求められ るよ うにする．
ことなど が 挙げ られ る．
また ，複数被験者による同時実験の今後の展開とし て，非侵襲小型連続血圧・心拍測定分析装置の標準化を進め，数十人から数百人規模の同時実験が 可能になるような装置への拡張性を検討する
ことが 考えられ る．

謝

辞

本研究は ，日本自転車振興会の 補助金を 受け，(財) 機械シ ステム振興協会，(社) 日本電子技術産業協会が 実施し たプ ロジェクト「 映像ディジ タル
コンテン ツ評価システムの開発に関するフィージビ リティスタデ ィ」の一部とし て行われ た．ここに謝意を表する．

参考文献
1) 中川千鶴，大須賀美恵子：ＶＥ酔い研究および 関連分野における研究の現状，ＶＲ学会論文誌，3(2)，31/39(1998)
2) 杉田典大，吉澤 誠，田中 明，阿部健一，山家智之，仁田新一：乗り物酔いし やすさと 3 次元映像
感受性の関連性の定量的評価，日本バーチャル リ
アリティ学会第６回大会論文集，121/124(2001)
3) Yoshizawa M., Sugita N., Tanaka A., Abe K.,
Yambe T., Nitta S.: Quantitative and physiological evaluation of three dimensional images, Proc.

–8–of the 7th International Conference on Virtual Systems and Multimedia, 864/871, Berkeley, U.S.A.
(2001)
4) Brian Gribbin, Andrew Steptoe, Peter Sleight:
Pulse wave velocity as a measure of blood pressure change,The Society for PsychophysiologicalResearch, Vol.13，No.1，86/90(1976)
5) Cerutti S., Baselli G., Bianchi A., Signorini M. G.:
Spectral techniques of analysis for blood pressureand heart rate signals, blood pressure andheart
rate variability, IOS Press, 39/55 (1992)
6) Pagani M., Federico L., Stefano G., Ornella R.,
Raﬀaellow F., Paolo P., Giulia S., Gabriella M.:
Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympathovagal interaction in man and conscious dog, Circ.
Res., 59(2), 178/193 (1986)