"二つのマイクロホンによって計測された二つの音響信号に基づいて聴覚臨場感を評価する聴覚臨場感評価装置であって、計測された前記二つの音響信号を周波数帯域ごとに分割する音響信号分割部と、分割された前記二つの音響信号を用いて、前記二つの前記音響信号に関する相互相関関数を前記周波数帯域ごとに算出する相互相関関数算出部と、算出された前記相互相関関数を、前記二つの音響信号のラグと、前記周波数帯域と、をフレームの座標軸とする二次元動画像の輝度、色相及び彩度のいずれかに変換することによって動画像を生成する動画像生成部と、生成された前記動画像の複数のフレームを用いて前記動画像の移動ベクトルを算出することによって、前記音像の移動ベクトルを算出する音像移動ベクトル算出部と、算出された前記音像の移動ベクトルに基づいて、前記音像の移動ベクトルが大きいほど聴覚臨場感評価値が大きくなるように当該聴覚臨場感評価値を算出する聴覚臨場感評価値算出部と、を備えることを特徴とする聴覚臨場感評価装置。前記動画像生成部は、前記相互相関関数を前記輝度に変換することを特徴とする請求項1に記載の聴覚臨場感評価装置。分割された前記二つの音響信号を用いて、前記二つの音響信号の音圧レベルを前記周波数帯域ごとに算出する音圧レベル算出部をさらに備え、前記動画像生成部は、算出された前記音圧レベルを、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数が変換されたもの以外のいずれかに変換することによって前記動画像を生成することを特徴とする請求項1に記載の聴覚臨場感評価装置。前記動画像生成部は、前記相互相関関数を前記輝度に変換するとともに、前記音圧レベルを前記彩度に変換することを特徴とする請求項3に記載の聴覚臨場感評価装置。前記音圧レベル算出部は、前記音圧レベルの差である音圧レベル差を前記周波数帯域ごとに算出し、前記動画像生成部は、算出された前記音圧レベル差を、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数及び前記音圧が変換されたもの以外に変換することによって前記動画像を生成することを特徴とする請求項3に記載の聴覚臨場感評価装置。前記動画像生成部は、前記相互相関関数を前記輝度に変換し、前記音圧レベルを前記彩度に変換するとともに、前記音圧レベル差を前記色相に変換することを特徴とする請求項5に記載の聴覚臨場感評価装置。前記相互相関関数算出部は、区間長の異なる複数の区間長別相互相関関数を算出し、算出された複数の区間長別相互相関関数に基づいて前記相互相関関数を算出することを特徴とする請求項1から請求項6のいずれか一項に記載の聴覚臨場感評価装置。前記音像移動ベクトル算出部は、前記動画像の三以上のフレームを用いて前記動画像の移動ベクトルを算出することを特徴とする請求項1から請求項7のいずれか一項に記載の聴覚臨場感評価装置。二つのマイクロホンによって計測された二つの音響信号に基づいて聴覚臨場感を評価する聴覚臨場感評価プログラムであって、コンピュータを、計測された前記二つの音響信号を周波数帯域ごとに分割する音響信号分割部、分割された前記二つの音響信号を用いて、前記二つの前記音響信号に関する相互相関関数を前記周波数帯域ごとに算出する相互相関関数算出部、算出された前記相互相関関数を、前記二つの音響信号のラグと、前記周波数帯域と、をフレームの座標軸とする二次元動画像の輝度、色相及び彩度のいずれかに変換することによって動画像を生成する動画像生成部、生成された前記動画像の複数のフレームを用いて前記動画像の移動ベクトルを算出することによって、前記音像の移動ベクトルを算出する音像移動ベクトル算出部、及び、算出された前記音像の移動ベクトルに基づいて、前記音像の移動ベクトルが大きいほど聴覚臨場感評価値が大きくなるように当該聴覚臨場感評価値を算出する聴覚臨場感評価値算出部、として機能させることを特徴とする聴覚臨場感評価プログラム。収音再生技術の発展に伴い、22.2マルチチャネル音響システム、WaveFieldSynthesis等といった高い臨場感を実現可能な音響システムが多数開発されており、音の品質の劣化度合いではなく、品質の向上度合いを客観的に評価することが求められている。ここで、KOzawa、YChujo、YSuzuki、TSone、Contentswhichyieldhighauditorypresenceinsoundreproduction、KanseiEngineeringInternational、2002年、Vol3、No4、p2530、及び、KOzawa、YChujo、ContentPresencevsSystemPresenceinAudioReproductionSystems、ProcoftheSecondInternationalSymposiumonUniversalCommunication、2008年、p5055に記載されているように、音響機器によって再生された音響信号による音像が、聴取者に対して相対的に移動している程、臨場感があると聴取者が感じることが明らかになっている。本願発明者は、鋭意研究の結果、音の品質の向上度合いを客観的に評価するためには、聴覚における臨場感、すなわち聴覚臨場感を定量化して客観的に評価すればよいのではないかという考えに至った。本発明は、前記した事情に鑑みて創案されたものであり、聴覚臨場感を客観的に評価することが可能な聴覚臨場感評価装置及び聴覚臨場感評価プログラムを提供することを課題とする。前記課題を解決するため、本発明の聴覚臨場感評価装置は、二つのマイクによって計測された二つの音響信号に基づいて聴覚臨場感を評価する聴覚臨場感評価装置であって、音響信号分割部と、相互相関関数算出部と、動画像生成部と、音像移動ベクトル算出部と、聴覚臨場感評価値算出部と、を備えることを特徴とする。かかる構成によると、音響信号分割部によって、計測された前記二つの音響信号を例えばある一定時間ごとにある区間長だけ切り出し、周波数帯域ごとに分割し、相互相関関数算出部によって、分割された前記二つの音響信号を用いて、前記二つの前記音響信号に関する相互相関関数を前記周波数帯域ごとに算出する。そして、動画像生成部によって、算出された前記相互相関関数を、前記二つの音響信号のラグと、前記周波数帯域と、をフレームの座標軸とする二次元動画像の輝度、色相及び彩度のいずれかに変換することによって動画像を生成し、音像移動ベクトル算出部によって、生成された前記動画像の複数のフレームを用いて前記動画像の移動ベクトルを算出することによって、前記音像の移動ベクトルを算出し、聴覚臨場感評価値算出部によって、算出された前記音像の移動ベクトルに基づいて、前記音像の移動ベクトルが大きいほど聴覚臨場感評価値が大きくなるように当該聴覚臨場感評価値を算出することができる。前記動画像生成部は、前記相互相関関数を前記輝度に変換する構成であってもよい。かかる構成によると、相互相関関数と相関する輝度成分を有するモノクロ動画像を生成し、生成されたモノクロ動画像の移動ベクトルから聴覚臨場感評価値を算出することができる。聴覚臨場感評価装置は、音圧レベル算出部をさらに備える構成であってもよい。かかる構成によると、音圧レベル算出部によって、分割された前記二つの音響信号を用いて、前記二つの音響信号の音圧レベルを前記周波数帯域ごとに算出し、動画像生成部によって、算出された前記音圧レベルを、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数が変換されたもの以外のいずれかに変換することによって前記動画像を生成することができる。すなわち、聴覚臨場感評価装置は、より多くのパラメータを用いて聴覚臨場感評価値を算出するので、聴覚臨場感の算出精度が向上する。前記動画像生成部は、前記相互相関関数を前記輝度に変換するとともに、前記音圧レベルを前記彩度に変換する構成であってもよい。かかる構成によると、相互相関関数と相関する輝度度成分及び音圧レベルと相関する彩度成分を有するカラー動画像を生成し、生成されたカラー動画像の移動ベクトルから聴覚臨場感評価値を算出することができる。また、前記音圧レベル算出部は、前記音圧レベルの差である音圧レベル差を前記周波数帯域ごとに算出し、前記動画像生成部は、算出された前記音圧レベル差を、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数及び前記音圧が変換されたもの以外に変換することによって前記動画像を生成する構成であってもよい。すなわち、聴覚臨場感評価装置は、より多くのパラメータを用いて聴覚臨場感評価値を算出するので、聴覚臨場感の算出精度が向上する。さらに、前記動画像生成部は、前記相互相関関数を前記輝度に変換し、前記音圧レベルを前記彩度に変換するとともに、前記音圧レベル差を前記色相に変換する構成であってもよい。かかる構成によると、相互相関関数と相関する輝度成分、音圧レベルと相関する彩度成分、及び音圧レベル差と相関する色相成分を有するカラー動画像を生成し、生成されたカラー動画像の移動ベクトルから聴覚臨場感評価値を算出することができる。前記相互相関関数算出部は、区間長の異なる複数の区間長別相互相関関数を算出し、算出された複数の区間長別相互相関関数に基づいて前記相互相関関数を算出する構成であってもよい。単一の区間長の相互相関関数のみを算出した場合には、偶々相関性が高く算出される等相互相関関数の値と実際の音像位置とが正確に一致しない場合があるが、かかる構成によると、複数の区間長を用いることで、音像位置の誤検出を防ぐことができる。前記音像移動ベクトル算出部は、前記動画像の三以上のフレームを用いて前記動画像の移動ベクトルを算出する構成であってもよい。かかる構成によると、移動ベクトルの算出精度が向上する。また、本発明は、コンピュータを前記した聴覚臨場感評価装置として機能させる聴覚臨場感評価プログラムとしても具現化可能である。本発明によれば、聴覚臨場感を客観的に評価することができる。以下、本発明の実施形態について、本発明を再生音場の評価に適用した場合を例にとり、適宜図面を参照しながら説明する。同様の部分には同一符号を付し、重複する説明を省略する。本発明において、音像は、聴取者が音源の存在を知覚する位置を指し、ステレオスピーカ等によって具現化された仮想的な音源の位置であったり、実際の音源の位置であったりする。また、本発明において、聴覚臨場感は、聴取者が聴覚において感じる臨場感を指す。従来、臨場感は、AV機器の性能を語る際の重要なキーワードであるにも関わらず、曖昧な意味で用いられていたのに対し、本発明の聴覚臨場感推定装置は、聴覚臨場感評価値を定量的かつ客観的に推定することができる。再生側の機器としては、ステレオスピーカ、5.1チャネルシステム、ヘッドホン等があるが、聴取者は、最終的に両耳に到達した音響信号によって聴覚臨場感を知覚する。したがって、本発明の聴覚臨場感評価装置は、聴取者の両耳と同様に、2チャネルの音響信号に基づいて音像の移動量を推定することとする。図1に示すように、本発明の実施形態に係る聴覚臨場感評価システム1は、再生側の機器として、記憶媒体再生装置10と、スピーカ群20と、を備えている。記憶媒体再生装置10は、図示しない記憶媒体に記憶されたデータを読み取り、ステレオスピーカ、5.1チャネルシステム、ヘッドホン等からなるスピーカ群20を介して音響信号を再生する。本発明の実施形態に係る聴覚臨場感評価システム1は、評価側の機器として、一対のマイクロホン30L,30Rと、聴覚臨場感評価装置40と、通知装置50と、を備える。<マイクロホン>一対のマイクロホン30L,30Rは、左右に配置された2チャネル構造を呈し、スピーカ群20で再生された音響信号を計測し、計測された2チャネルの音響信号を聴覚臨場感評価装置40へ出力する。本明細書では、スピーカ群20で再生された音波だけでなく、当該音波のマイクロホン30L,30Rによる計測結果も音響信号と記載する。<聴覚臨場感評価装置>聴覚臨場感評価装置40は、CPUCentralProcessingUnit、RAMRandomAccessMemory、ROMReadOnlyMemory、入出力回路等から構成されており、機能部として、音響信号分割部41と、相互相関関数算出部42と、音圧レベル算出部43と、動画像生成部44と、音像移動ベクトル算出部45と、聴覚臨場感評価値算出部46と、を備える。≪音響信号分割部≫音響信号分割部41は、一対のマイクロホン30L,30Rから出力された2チャネルの音響信号をある一定時間ごとにある区間長だけ取得し、取得された2チャネルの音響信号をM個のバンドパスフィルタによって周波数帯域ごとにM個に分割し、分割された音響信号を相互相関関数算出部42及び音圧レベル算出部43へ出力する。本実施形態において、音響信号分割部41は、5個のオクターブバンドパスフィルタを備えており、マイクロホン30Lから出力された音響信号Lを、中心周波数が125[Hz]である周波数帯域f1の第一の音響信号L1、中心周波数が250[Hz]である周波数帯域f2の第二の音響信号L2、中心周波数が500[Hz]である周波数帯域f3の第三の音響信号L3、中心周波数が1000[Hz]である周波数帯域f4の第四の音響信号L4、中心周波数が2000[Hz]である周波数帯域f5の第五の音響信号L5に分割する。同様に、音響信号分割部41は、マイクロホン30Rから出力された音響信号を、中心周波数が125[Hz]である周波数帯域f1の第一の音響信号R1、中心周波数が250[Hz]である周波数帯域f2の第二の音響信号R2、中心周波数が500[Hz]である周波数帯域f3の第三の音響信号R3、中心周波数が1000[Hz]である周波数帯域f4の第四の音響信号R4、中心周波数が2000[Hz]である周波数帯域f5の第五の音響信号R5に分割する。なお、音響信号の分割数は、5個に限定されず、対象とする音像の種類に応じて適宜変更可能である。例えば、自動車騒音のように広帯域な雑音成分を含む音を対象とする場合には、可聴周波数帯域20[Hz]~20[kHz]を1/1オクターブバンド幅で区切ったM=10程度とすることができる。また、楽器音のように高調波構造を有する音を対象とする場合には、1/10オクターブバンド幅で区切ったM=40程度とすることができる。また、汎用的なモデルのバンドパスフィルタとして、50チャネル程度の聴覚フィルタ聴覚末梢系の特性を考慮したバンドパスフィルタであり、代表的なものにROEXフィルタ、Gamma-toneフィルタ等があるを採用することによって、聴覚特性を反映した分割を行うことができる。≪相互相関関数算出部≫相互相関関数算出部42は、音響信号分割部41から出力された各音響信号L1~L5,R1~R5を取得し、取得された各音響信号L1~L5,R1~R5に基づいて、周波数帯域ごとの相互相関関数を算出し、算出された相互相関関数を動画像生成部44へ出力する。本実施形態において、相互相関関数算出部42は、算出周期ごと時刻t1,t2,...に、音響信号L1,R1の相互相関CC:CrossCorrelation関数を算出する。相互相関関数算出部42は、図2aに示すように、時刻t1を先頭として区切った区間長時間窓幅ともいうT1の音響信号L1,R1を抽出して抽出された音響信号L1,R1の相関値を算出するとともに、音響信号R1のみをラグmずつずらして区切って抽出し直して相関値を算出することを繰り返すことによって、図3aの下段に示すように、ラグmを変数とする相互相関関数を算出する。相互相関関数は、-1から1までの値をとり、音響信号L1,R1が一致した場合に1となり、音響信号L1,R1の波形が異なるほど値が小さくなる。本発明において、ラグmは、音響信号L1,R1における相互相関関数を算出するための区間のずれ時間差である。前記した時刻t1,t2,...を先頭として区切った区間長T1、p番目の算出周期、q番目の周波数帯域における時刻nの関数である音響信号Lp,qn,Rp,qnに対する相互相関関数rp,qmは、下記式によって算出される。ここで、σLp,q,σRp,qは、それぞれLp,qn,Rp,qnの標準偏差である。また、Nは、相互相関関数を算出する期間の長さTに対応する離散時間である。なお、区間長は、対象とする音像の種類に応じて適宜変更可能である。例えば、定常的な騒音を対象とする場合には、100~200[ms]程度とし、時間的な変化が大きい楽器音を対象とする場合には、10~20[ms]程度とすることができる。かかる区間長は、相互相関関数算出部42に予め記憶されていてもよく、聴覚臨場感評価装置40に接続されたキーボード、マウス等からなる入力装置図示せずをユーザが操作することによって設定される構成であってもよい。また、算出周期時刻t1,t2,...の間隔と、区間長とは、一致しても一致しなくてもよい。また、相互相関関数算出部42は、聴覚臨場感評価値算出部46によって算出された音像の移動量移動ベクトルの大きさに基づいて、算出周期を変更する構成であってもよい。例えば、音像の移動量が小さい場合には、算出周期を長く設定することによって演算量を削減することができる。この場合において、二次元動画像の一般的なフレームレートである1/30[s]を算出周期として相互相関関数を算出する場合であっても、区間長を1/30[s]よりも長く設定し、音響信号R1の抽出範囲をオーバーラップさせながら相互相関関数を算出することができる。相互相関関数算出部42には、音像の移動量と、相互相関関数の算出周期と、の関係がデータベースとして予め記憶されており、相互相関関数算出部42は、聴覚臨場感評価値算出部46によって算出された音像の移動量を用いてデータベースを検索することによって、音像の移動量の対応する算出周期を読み出し、読み出された算出周期に基づいて相互相関関数を算出する。また、相互相関関数算出部42は、区間長の異なる複数の区間長別相互相関関数を算出し、算出された複数の区間長別相互相関関数に基づいて相互相関関数を算出する構成であってもよい。例えば、相互相関関数算出部42には、区間長T1,T2,T3が予め記憶されており、相互相関関数算出部42は、図2bに示すように、区間長T1における区間長別相互相関関数を算出するとともに、区間長T1よりも長い区間長T2,T3における区間長別相互相関関数を算出し、算出された3個の区間長別相互相関関数の相加平均を算出することによって、相互相関関数を算出することができる。相互相関関数算出部42は、他の周波数帯域f2~f5に関しても同様に相互相関関数を算出する。≪音圧レベル算出部≫音圧レベル算出部43は、音響信号分割部41から出力された各音響信号L1~L5,R1~R5を取得し、取得された各音響信号L1~L5,R1~R5に基づいて、周波数帯域ごとの音圧レベル、及び、左右の音響信号の音圧の差である音圧レベル差を算出し、算出された音圧レベル及び音圧レベル差を動画像生成部44へ出力するとともに、算出された音圧レベルを聴覚臨場感評価値算出部46へ出力する。≪動画像生成部≫動画像生成部44は、相互相関関数算出部42から出力された周波数帯域ごとの相互相関関数と、音圧レベル算出部43から出力された音圧レベル及び音圧レベル差に基づいて、動画像より詳細には、二次元動画像を生成し、生成された動画像を音像移動ベクトル算出部45へ出力する。図5aに示すように、動画像は、時刻t1におけるフレーム100、時刻t2におけるフレーム200、時刻t3におけるフレーム300、時刻t4におけるフレーム400、...から構成されており、動画像の一のフレームは、縦軸を周波数周波数帯域f1~f5、横軸をラグmとする平面である。本実施形態では、一のフレームは、縦5×横8の画素から構成されており、動画像生成部44は、相互相関関数を対応する画素の輝度、音圧レベルを対応する画素の彩度、音圧レベル差を対応する画素の色相に変換することによって、動画像を生成する。例えば、相互相関関数のみを用いて動画像を生成する場合には、動画像生成部44は、相互相関関数が-1の場合に最低輝度となり、相互相関関数が1の場合に最高輝度となるように各画素の画素データを生成し図3b参照、かかる画素データによって構成された一連の画像からなる動画像を生成する。ここで、動画像生成部44は、図3aの下段に示す周波数帯域f1の相互相関関数から、図3cの下段に示す輝度を有する画素列110を生成し、図3aの上段に示す周波数帯域f2の相互相関関数から、図3cの上段に示す輝度を有する画素列120を生成する。また、動画像生成部44は、図4aに示す周波数帯域ごとの音響信号L1又はL2~L5,R1又はR2~R5の音圧レベル210,220の音圧レベル差230に基づいて、図4bに示すラグmごとの重みを算出し、算出された重みによって重み付けされた輝度を有する画素列110’を生成する構成であってもよい。相関の高い音響信号が左右の耳にある音圧レベル差で到来した場合、左右のレベル差によって、正面を0度として±90度のいずれかの方向に音像が知覚される。動画像生成部44は、かかる音圧レベル差と音像の知覚方向との関係を予め記憶しておき、音圧レベル差230を用いてかかる関係を参照することによって、何度の方向に音像があるかを推測することができる。動画像生成部44は、音像があると推定される方向に重みとして1を設定し、そこから遠ざかる方向につれて、重みとして小さい値を設定する。一方、相互相関関数の値が最大になるラグmが音像のある方向であるから、ラグmと音像の方向とは一対一に対応する。したがって、動画像生成部44は、重み付けによって、音圧レベル差230に基づいて音像があると推定された付近の相互相関関数の値はそのままにし、音圧レベル差230に基づいて音像がないと推定された付近の相互相関関数の値は小さくなるようにする。こうすることで、正弦波のように相互相関関数が1となるラグが数か所でてくるような場合でも、動画像生成部44は、音圧レベル差230に基づいて音像の方向を判定することができる。また、動画像生成部44は、相互相関関数を輝度又は色相に変換して動画像を生成する場合には、相互相関関数がマイナスの部分を0とみなして動画像を生成する構成であってもよい。これは、音像が1個の場合には、相互相関関数が最大となるラグに対応する位置が音像の位置と一致するためである。≪音像移動ベクトル算出部≫音像移動ベクトル算出部45は、動画像生成部44から出力された動画像の複数のフレームを用いて動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出し、算出された音像の移動ベクトルを聴覚臨場感評価値算出部46へ出力する。動画像の移動ベクトルの算出手法としては、オプティカルフロー、パターンマッチング、特徴量追従アルゴリズム等といった、画像情報処理のための公知のアルゴリズムを適宜採用することができる。音像移動ベクトル算出部45は、3個以上本実施形態では、時刻t1~t4に対応する4個のフレームを蓄積し、蓄積されたフレームの輝度、色相及び彩度に基づいて、周波数帯域ごとの移動ベクトルを算出する。詳細には、音像移動ベクトル算出部45は、図5bに示すように、時刻t1における画素列111ここでは、説明のため、1列の画素数を12個とする。、時刻t2における画素列211、時刻t3における画素列311を用いて、周波数帯域f1における移動ベクトルを算出するとともに、他の周波数帯域f2~f5についても同様の処理を行う。ここで、時刻t1,t2の2個のフレームを用いて移動ベクトルを算出する場合には、音像移動ベクトル算出部45は、時刻t1における画素列111の左から4個目の画素が、時刻t2においても画素列211の左から4個目に留まっていると誤算出するおそれがある実線の矢印。そのため、音像移動ベクトル算出部45は、時刻t3における画素列311をさらに用いることによって、時刻t1における画素列111の左から1~3個目の画素が、時刻t2における画素列211では左から5~7個目、時刻t3における画素列では左から10~12個目へと移動していると推定するため、時刻t1における画素列111の左から4個目の画素が、時刻t2においては画素列211の左から8個目へ移動していると算出することができる点線の矢印。このように、3個以上のフレームを用いて移動ベクトルを算出することによって、移動ベクトルの誤算出を防ぐことができる。また、音像移動ベクトル算出部45は、周波数特性が変動する音像に対して、周波数帯域ごとの音圧レベルの変化に応じて、周波数帯域を跨ぐように音像の移動ベクトルを算出してもよい。ここで、音像移動ベクトル算出部45は、音圧レベルの周波数成分の時間変化に基づいて、周波数、時間及び音圧レベルを軸とする三次元空間における音圧レベルのピーク尾根部分が周波数の軸上でシフトした場合に、音像の周波数特性が変動していると判定し、周波数帯域を跨ぐように音像の移動ベクトルを算出することができる。≪聴覚臨場感評価値算出部≫聴覚臨場感評価値算出部46は、音像移動ベクトル算出部45から出力された音像の移動ベクトルに基づいて、聴覚臨場感評価値を算出し、通知装置50へ出力する。本実施形態において、聴覚臨場感評価値算出部46は、閾値未満の音像の移動量を切り捨てるとともに、閾値以上の複数の音像の移動量を相加平均することによって、聴覚臨場感評価値を算出することができる。かかる聴覚臨場感評価値は、音像の移動ベクトルが大きいほど、当該聴覚臨場感評価値が大きくなるような値である。聴覚臨場感評価値算出部46は、動画像の周波数帯域ごとの移動ベクトル、すなわち、音像の移動ベクトルに基づいて、周波数帯域ごとの音像移動量を推定し、推定された周波数帯域ごとの音像移動量の最大値を採用する、音像移動量の相加平均を算出する、又は、音像移動量の積算値を算出することによって、全体としての音像移動量を推定する。本実施形態では、聴覚臨場感評価値算出部46は、音像移動量の積算値を採用する。すなわち、聴覚臨場感評価値算出部46は、図6aに示すように、前フレームと今フレームにおける各周波数帯域f1~f5の移動ベクトルV1a,V1b,V2~V5を算出する。続いて、聴覚臨場感評価値算出部46は、図6bに示すように、各周波数帯域f1~f5の移動ベクトルを統合する。ここでは、周波数帯域f2の移動ベクトルV2と周波数帯域f5の移動ベクトルV5が類似しているため統合されるとともに、周波数帯域f1の左の移動ベクトルV1aと周波数帯域f1の左の移動ベクトルV1bと周波数帯域f4の移動ベクトルV4とが類似しているため統合される。ここで、聴覚臨場感評価値算出部46は、二つの移動ベクトルの向きが同じで、大きさがほぼ同じ例えば、移動ベクトルの大きさの比率が一定範囲内であれば、二つの移動ベクトルが類似していると判定する構成であってもよく、さらに、二つの移動ベクトルの音圧レベル及び音圧レベル差の少なくとも一方がほぼ同じ比率が一定範囲内である場合に、二つの移動ベクトルが類似していると判定する構成であってもよい。続いて、聴覚臨場感評価値算出部46は、図6cに示すように、統合後の移動ベクトルV1b=V1a=V4,V3,V5の大きさの和積分値を算出し、算出結果を音像移動量とする。なお、聴覚臨場感評価値算出部46は、図6dに示すように、統合後の移動ベクトルV1b,V3,V5を音圧によって補正し、補正後の移動ベクトルV1b’,V3’,V5’の大きさの和積分値を算出し、算出結果を音像移動量とする構成であってもよい。例えば、音圧レベルが大きいほど聴覚臨場感評価値に対する音像の移動ベクトルの寄与度が高くなるように重み付けを行うことによって、聴覚臨場感評価値の推定精度を向上させることができる。続いて、聴覚臨場感評価値算出部46は、今フレームと次フレームにおける音像移動量を算出し、かかる算出を所定フレーム数に達するまで繰り返し、所定フレーム数分の音像移動量の和積算値を算出し、算出結果を最終的な音像移動量、すなわち、聴覚臨場感評価値とする。また、聴覚臨場感評価値算出部46は、フレーム数の異なる複数のフレーム数別聴覚臨場感評価値を算出し、算出された複数のフレーム数別聴覚臨場感評価値に基づいて聴覚臨場感評価値を算出する構成であってもよい。例えば、聴覚臨場感評価値算出部46は、3個のフレームにおけるフレーム数別聴覚臨場感評価値と、4個のフレームにおけるフレーム数別聴覚臨場感評価値と、を算出し、算出された2個のフレーム数別聴覚臨場感評価値の相加平均を算出することによって、聴覚臨場感評価値を算出することができる。また、音像が1個の場合には、聴覚臨場感評価値算出部46は、輝度すなわち、相互相関関数が最大となるラグmに基づいて、ある時刻における音像の方向方位角を推定することもできる。これは、音像が1個の場合には、相互相関関数が最大となるラグmが、音像から各マイクロホン30L,30Rまでの距離の差と関連するためである。また、聴覚臨場感評価値算出部46は、複数の周波数帯域の音圧レベルが同期して変動している場合には、当該複数の周波数帯域に同じ音像が存在するとみなして音像の移動ベクトルを統合することもできる。かかる構成によると、聴取者が知覚するのに似た音像の移動ベクトルを用いて聴覚臨場感評価値を算出することができる。<通知装置>通知装置50は、ディスプレイ、スピーカ等からなり、聴覚臨場感評価値算出部46から出力された聴覚臨場感評価値を表示したり音声出力したりすることによって、利用者へ通知する。<動作例>続いて、本発明の実施形態に係る聴覚臨場感評価装置40の動作例について説明する。まず、音響信号分割部41が、一対のマイクロホン30L,30Rから出力された2チャネルの音響信号を取得し、取得された2チャネルの音響信号を5個のバンドパスフィルタによって周波数帯域が異なる5個の音響信号L1~L5,R1~R5に分割する。続いて、相互相関関数算出部42が、各音響信号L1~L5,R1~R5に基づいて、周波数帯域ごとの相互相関関数を算出する。続いて、音圧レベル算出部43が、各音響信号L1~L5,R1~R5に基づいて、周波数帯域ごとの音圧レベル、及び、左右の音響信号の音圧の差である音圧レベル差を算出し、算出された音圧レベル及び音圧レベル差を動画像生成部44へ出力するとともに、算出された音圧レベルを聴覚臨場感評価値算出部46へ出力する。続いて、動画像生成部44が、相互相関関数と、音圧レベル算出部43から出力された音圧レベル及び音圧レベル差に基づいて、相互相関関数を輝度、音圧レベルを彩度、音圧レベル差を色相とみなした各画素の画素データを生成し、かかる画素データによって構成された一連の画像からなる動画像を生成する。続いて、音像移動ベクトル算出部45が、動画像生成部44から出力された動画像の複数のフレームを用いて動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出する。続いて、聴覚臨場感評価値算出部46が、音像移動ベクトル算出部45から出力された音像の移動ベクトルに基づいて、聴覚臨場感評価値を算出し、通知装置50へ出力する。<音像の移動ベクトルの算出例>続いて、本発明の実施形態に係る聴覚臨場感評価装置40の音像移動ベクトル算出部45が、移動する音像の移動ベクトルを実際に算出した例について説明する。ここで、マイクロホン30L,30Rとして、ダミーヘッド高研社製のSAMRAIを使用するとともに、スピーカ群20として、単一のスピーカを使用し、このスピーカを、ダミーヘッドを中心とした半径0.5[m]の円上で、左-90[deg]から右90[deg]まで、又は右から左まで移動させたときの音響信号をダミーヘッドが検出した。また、音響信号分割部41のサンプリング周波数は44100[Hz]であり、中心周波数125,250,500,1000,2000,4000[Hz]とするオクターブバンドパスフィルタを用いて音響信号分割部41を構成することにより、音響信号を6個の周波数帯域に分割した。相互相関関数算出部42は、相互相関関数の算出周期を1/30[s]、区間長を1/5[s]として相互相関関数を算出した。ここで、ダミーヘッドにおけるマイクロホン30L,30R間の距離は0.18[m]であり、音速は340[m/s]であるため、音響信号がマイクロホン30L,30Rに到達する時間差が最大となるスピーカが左又は右にある場合の時間差は、約0.7[ms]となる。そこで、ラグmとしては、0を中心とした±0.7[ms]に相当する部分にのみ着目した。また、音像移動ベクトル算出部45は、相互相関関数の最大値rp-1,qmmax,rp,qmmaxを有する画素同士を結ぶ手法を用いて動画像の移動ベクトル、すなわち、音像の移動ベクトルを算出した。≪単一音像の場合≫単一のスピーカから音響信号としてホワイトノイズを出力した場合、すなわち、単一の音源を用いて一の音像を生成した場合における、計測開始時点において動画像生成部44が生成した画像フレームを図7aに示し、現在途中経過時点において動画生成部44が生成した画像を図7bに示す。図7bにおける矢印は、開始時点から現在までの、周波数帯域ごとの音像の移動ベクトルである。図7abを参照すると、全ての周波数帯域において音像が左から右へ移動している様子が分かる。≪複数音像の場合≫左から右へ移動するスピーカからローパスノイズ遮断周波数700[Hz]を出力した音響信号と、右から左へ移動するスピーカからハイパスノイズ遮断周波数1000[Hz]を出力した音響信号と、を合成したものを音響信号分割部41に入力して分析を行った。この場合では、二つの音源を用いて二つの音像を生成している。この場合における、測開始時点において動画像生成部44が生成した画像フレームを図8aに示し、現在途中経過時点において動画生成部44が生成した画像を図8bに示す。図8bにおける矢印は、開始時点から現在までの、周波数帯域ごとの音像の移動ベクトルである。図8abを参照すると、低周波数帯域においては音像が左から右へ移動し、高周波数帯域においては音像が右から左へ移動している様子が分かる。本発明の実施形態に係る聴覚臨場感評価装置40は、複数の音像に対して処理を行う場合であっても、音像の方向を用いずに音像の移動量に基づいて聴覚臨場感評価値を算出するので、簡易な処理で聴覚臨場感を客観的に評価することができる。以上、本発明の実施形態について実施形態を参照して説明したが、本発明は前記実施形態に限定されず、本発明の要旨を逸脱しない範囲で適宜設計変更可能である。例えば、音圧レベル算出部43を備えずに、動画像生成部44が相互相関関数を輝度と変換したモノクロの動画像を生成する構成であってもよい。また、動画像生成部44が相互相関関数を輝度、音圧レベルを彩度としたカラーの動画像を生成する構成であってもよい。この場合には、動画像生成部44は、各画素の色相を所定値例えば、色相値50に設定することができる。また、相互相関関数、音圧レベル及び音圧レベル差を、輝度、彩度及び色相のいずれかに変換して動画像を生成する構成であってもよい。また、聴覚臨場感評価装置40は、記憶媒体に記憶されていない音響信号を再生する装置、再生装置ではない実際の音源に対しても、聴覚臨場感評価値を算出することができる。また、音像移動ベクトル算出部45は、複数の周波数帯域の画素列を用いて一つの移動ベクトルを算出することによって、統合された移動ベクトルを算出する構成であってもよく、移動ベクトルごとの画像の特徴量をパターンとして予め記憶しており、パターンマッチングを用いて移動ベクトルを算出する構成であってもよい。また、本発明は、コンピュータを聴覚臨場感評価装置40として機能させる聴覚臨場感評価プログラムとしても具現化可能である。"