(57)【特許請求の範囲】
【請求項１】
二つのマイクロホンによって計測された二つの音響信号に基づいて聴覚臨場感を評価する聴覚臨場感評価装置であって、
計測された前記二つの音響信号を周波数帯域ごとに分割する音響信号分割部と、
分割された前記二つの音響信号を用いて、前記二つの前記音響信号に関する相互相関関数を前記周波数帯域ごとに算出する相互相関関数算出部と、
算出された前記相互相関関数を、前記二つの音響信号のラグと、前記周波数帯域と、をフレームの座標軸とする二次元動画像の輝度、色相及び彩度のいずれかに変換することによって動画像を生成する動画像生成部と、
生成された前記動画像の複数のフレームを用いて前記動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出する音像移動ベクトル算出部と、
算出された前記音像の移動ベクトルに基づいて、前記音像の移動ベクトルが大きいほど聴覚臨場感評価値が大きくなるように当該聴覚臨場感評価値を算出する聴覚臨場感評価値算出部と、
を備えることを特徴とする聴覚臨場感評価装置。
【請求項２】
前記動画像生成部は、前記相互相関関数を前記輝度に変換する
ことを特徴とする請求項１に記載の聴覚臨場感評価装置。
【請求項３】
分割された前記二つの音響信号を用いて、前記二つの音響信号の音圧レベルを前記周波数帯域ごとに算出する音圧レベル算出部をさらに備え、
前記動画像生成部は、算出された前記音圧レベルを、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数が変換されたもの以外のいずれかに変換することによって前記動画像を生成する
ことを特徴とする請求項１に記載の聴覚臨場感評価装置。
【請求項４】
前記動画像生成部は、前記相互相関関数を前記輝度に変換するとともに、前記音圧レベルを前記彩度に変換する
ことを特徴とする請求項３に記載の聴覚臨場感評価装置。
【請求項５】
前記音圧レベル算出部は、前記音圧レベルの差である音圧レベル差を前記周波数帯域ごとに算出し、
前記動画像生成部は、算出された前記音圧レベル差を、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数及び前記音圧が変換されたもの以外に変換することによって前記動画像を生成する
ことを特徴とする請求項３に記載の聴覚臨場感評価装置。
【請求項６】
前記動画像生成部は、前記相互相関関数を前記輝度に変換し、前記音圧レベルを前記彩度に変換するとともに、前記音圧レベル差を前記色相に変換する
ことを特徴とする請求項５に記載の聴覚臨場感評価装置。
【請求項７】
前記相互相関関数算出部は、区間長の異なる複数の区間長別相互相関関数を算出し、算出された複数の区間長別相互相関関数に基づいて前記相互相関関数を算出する
ことを特徴とする請求項１から請求項６のいずれか一項に記載の聴覚臨場感評価装置。
【請求項８】
前記音像移動ベクトル算出部は、前記動画像の三以上のフレームを用いて前記動画像の移動ベクトルを算出する
ことを特徴とする請求項１から請求項７のいずれか一項に記載の聴覚臨場感評価装置。
【請求項９】
二つのマイクロホンによって計測された二つの音響信号に基づいて聴覚臨場感を評価する聴覚臨場感評価プログラムであって、
コンピュータを、
計測された前記二つの音響信号を周波数帯域ごとに分割する音響信号分割部、
分割された前記二つの音響信号を用いて、前記二つの前記音響信号に関する相互相関関数を前記周波数帯域ごとに算出する相互相関関数算出部、
算出された前記相互相関関数を、前記二つの音響信号のラグと、前記周波数帯域と、をフレームの座標軸とする二次元動画像の輝度、色相及び彩度のいずれかに変換することによって動画像を生成する動画像生成部、
生成された前記動画像の複数のフレームを用いて前記動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出する音像移動ベクトル算出部、及び、
算出された前記音像の移動ベクトルに基づいて、前記音像の移動ベクトルが大きいほど聴覚臨場感評価値が大きくなるように当該聴覚臨場感評価値を算出する聴覚臨場感評価値算出部、
として機能させることを特徴とする聴覚臨場感評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音響信号に基づいて聴覚における臨場感を評価する聴覚臨場感評価装置及び聴覚臨場感評価プログラムに関する。
【背景技術】
【０００２】
複数話者が存在する場合の音声認識の前処理、ロボットの自律運動制御等に有用な技術として、音源方向の推定に関する研究が行われてきた（非特許文献１参照）。
かかる音源方向の推定技術は、多チャンネルマイクロホンによって構成されるマイクロホンアレー信号処理を行うものであり、複数のマイクロホンを配置するため、大規模な装置が必要であった。
【０００３】
一方、人間は左右の耳で聞こえる２チャネルの音響信号を用いて十分な精度の音源方向の推定を実現しており、これを模して、２チャネルの音響信号を用いて音源方向の推定を行う研究も行われている（非特許文献２参照）。
かかる２チャネルの音響信号を用いた音源方向の推定技術は、装置の小型化を図ることができるが、情報量が少ないため、推定の精度を確保するためには複雑な信号処理が必要であった（非特許文献３参照）。
【０００４】
また、従来の２チャネルの音響信号を用いた音源方向の推定技術は、一つの音源の方向しか推定することができず、また、推定対象が静止音源に限られていた（特許文献１参照）。
すなわち、従来の音源方向の推定技術は、高精度であることを目指すが故に装置として大規模になりがちで実時間処理が困難である上、一つの音源の方向しか推定することができないという問題を有していた。
【０００５】
一方、符号化された音に対する客観的な評価法として、ＰＥＡＱ（Perceptual Evaluation of Audio Quality）が開発されており、かかる評価法は、ＩＴＵ（International Telecommunications Union：国際電気通信連合）によって規格化されている（ITU-R BS. 1387）。
かかる評価法は、ヒトの聴覚末梢系を模倣しているが、聴覚中枢系の部分に関しては、簡易なニューラルネットワークで代用している。
【０００６】
また、ヒトが感じる音声品質の客観的な評価法として、特許文献２に記載された評価法が存在する。かかる評価法は、基準となる音に関する劣化度合いを客観的に評価する手法である。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開平５−８７９０３号公報
【特許文献２】特開２００４−１７２７５３号公報
【非特許文献】
【０００８】
【非特許文献１】Arthur N. Popper、Richard R. Fay、Eds., 「Sound Source Localization」、Springer Handbook of Auditory Research、2005年、Springer、New York
【非特許文献２】永田仁史、「２チャネル信号に基づく到来音方向推定の計算モデル」、日本比較生理生化学会、2010年、第27巻、第1号、p.10-18
【非特許文献３】Nicoleta Roman、DeLiang Wang、「Binaural Tracking of Multiple Moving Sources」、IEEE Transactions on Audio、Speech, and Language Processing、2008年、Vol. 16、No. 4、p. 728-739
【発明の概要】
【発明が解決しようとする課題】
【０００９】
収音再生技術の発展に伴い、２２．２マルチチャネル音響システム、Wave Field Synthesis等といった高い臨場感を実現可能な音響システムが多数開発されており、音の品質の劣化度合いではなく、品質の向上度合いを客観的に評価することが求められている。
【００１０】
ここで、K. Ozawa、Y. Chujo、Y. Suzuki、T. Sone、「Contents which yield high auditory-presence in sound reproduction」、Kansei Engineering International、2002年、Vol. 3、No. 4、p. 25-30、及び、K. Ozawa、Y. Chujo、「Content Presence vs. System Presence in Audio Reproduction Systems」、Proc. of the Second International Symposium on Universal Communication (ISUC2008)、2008年、p. 50-55に記載されているように、音響機器によって再生された音響信号による音像が、聴取者に対して相対的に移動している程、臨場感があると聴取者が感じることが明らかになっている。
【００１１】
本願発明者は、鋭意研究の結果、音の品質の向上度合いを客観的に評価するためには、聴覚における臨場感、すなわち聴覚臨場感を定量化して客観的に評価すればよいのではないかという考えに至った。
【００１２】
本発明は、前記した事情に鑑みて創案されたものであり、聴覚臨場感を客観的に評価することが可能な聴覚臨場感評価装置及び聴覚臨場感評価プログラムを提供することを課題とする。
【課題を解決するための手段】
【００１３】
前記課題を解決するため、本発明の聴覚臨場感評価装置は、二つのマイクによって計測された二つの音響信号に基づいて聴覚臨場感を評価する聴覚臨場感評価装置であって、音響信号分割部と、相互相関関数算出部と、動画像生成部と、音像移動ベクトル算出部と、聴覚臨場感評価値算出部と、を備えることを特徴とする。
【００１４】
かかる構成によると、音響信号分割部によって、計測された前記二つの音響信号を例えばある一定時間ごとにある区間長だけ切り出し、周波数帯域ごとに分割し、相互相関関数算出部によって、分割された前記二つの音響信号を用いて、前記二つの前記音響信号に関する相互相関関数を前記周波数帯域ごとに算出する。そして、動画像生成部によって、算出された前記相互相関関数を、前記二つの音響信号のラグと、前記周波数帯域と、をフレームの座標軸とする二次元動画像の輝度、色相及び彩度のいずれかに変換することによって動画像を生成し、音像移動ベクトル算出部によって、生成された前記動画像の複数のフレームを用いて前記動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出し、聴覚臨場感評価値算出部によって、算出された前記音像の移動ベクトルに基づいて、前記音像の移動ベクトルが大きいほど聴覚臨場感評価値が大きくなるように当該聴覚臨場感評価値を算出することができる。
【００１５】
前記動画像生成部は、前記相互相関関数を前記輝度に変換する構成であってもよい。かかる構成によると、相互相関関数と相関する輝度成分を有するモノクロ動画像を生成し、生成されたモノクロ動画像の移動ベクトルから聴覚臨場感評価値を算出することができる。
【００１６】
聴覚臨場感評価装置は、音圧レベル算出部をさらに備える構成であってもよい。
【００１７】
かかる構成によると、音圧レベル算出部によって、分割された前記二つの音響信号を用いて、前記二つの音響信号の音圧レベルを前記周波数帯域ごとに算出し、動画像生成部によって、算出された前記音圧レベルを、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数が変換されたもの以外のいずれかに変換することによって前記動画像を生成することができる。すなわち、聴覚臨場感評価装置は、より多くのパラメータを用いて聴覚臨場感評価値を算出するので、聴覚臨場感の算出精度が向上する。
【００１８】
前記動画像生成部は、前記相互相関関数を前記輝度に変換するとともに、前記音圧レベルを前記彩度に変換する構成であってもよい。かかる構成によると、相互相関関数と相関する輝度度成分及び音圧レベルと相関する彩度成分を有するカラー動画像を生成し、生成されたカラー動画像の移動ベクトルから聴覚臨場感評価値を算出することができる。
【００１９】
また、前記音圧レベル算出部は、前記音圧レベルの差である音圧レベル差を前記周波数帯域ごとに算出し、前記動画像生成部は、算出された前記音圧レベル差を、前記輝度、前記色相及び前記彩度のうち、前記相互相関関数及び前記音圧が変換されたもの以外に変換することによって前記動画像を生成する構成であってもよい。すなわち、聴覚臨場感評価装置は、より多くのパラメータを用いて聴覚臨場感評価値を算出するので、聴覚臨場感の算出精度が向上する。
【００２０】
さらに、前記動画像生成部は、前記相互相関関数を前記輝度に変換し、前記音圧レベルを前記彩度に変換するとともに、前記音圧レベル差を前記色相に変換する構成であってもよい。かかる構成によると、相互相関関数と相関する輝度成分、音圧レベルと相関する彩度成分、及び音圧レベル差と相関する色相成分を有するカラー動画像を生成し、生成されたカラー動画像の移動ベクトルから聴覚臨場感評価値を算出することができる。
【００２１】
前記相互相関関数算出部は、区間長の異なる複数の区間長別相互相関関数を算出し、算出された複数の区間長別相互相関関数に基づいて前記相互相関関数を算出する構成であってもよい。単一の区間長の相互相関関数のみを算出した場合には、偶々相関性が高く算出される等相互相関関数の値と実際の音像位置とが正確に一致しない場合があるが、かかる構成によると、複数の区間長を用いることで、音像位置の誤検出を防ぐことができる。
【００２２】
前記音像移動ベクトル算出部は、前記動画像の三以上のフレームを用いて前記動画像の移動ベクトルを算出する構成であってもよい。かかる構成によると、移動ベクトルの算出精度が向上する。
【００２３】
また、本発明は、コンピュータを前記した聴覚臨場感評価装置として機能させる聴覚臨場感評価プログラムとしても具現化可能である。
【発明の効果】
【００２４】
本発明によれば、聴覚臨場感を客観的に評価することができる。
【図面の簡単な説明】
【００２５】
【図１】本発明の実施形態に係る聴覚臨場感評価システムを示すブロック図である。
【図２】一対のマイクロホンによって計測された音響信号を示すグラフである。
【図３】（ａ）は相互相関関数を示すグラフであり、（ｂ）は、相互相関関数の値と輝度との関係を示す図であり、（ｃ）は、相互相関関数から生成された画素列を示す図である。
【図４】（ａ）は音圧レベル差を示すグラフであり、（ｂ）は音圧レベル差による重みを示すグラフであり、（ｃ）は画素列に対する重み付けを示す図である。
【図５】（ａ）は動画像を示す図であり、（ｂ）は移動ベクトルを示す図である。
【図６】（ａ）〜（ｄ）は、音像移動量の算出手法を説明するための図である。
【図７】（ａ）（ｂ）は、単一音像の移動ベクトルの算出例を説明するための図である。
【図８】（ａ）（ｂ）は、複数音像の移動ベクトルの算出例を説明するための図である。
【発明を実施するための形態】
【００２６】
以下、本発明の実施形態について、本発明を再生音場の評価に適用した場合を例にとり、適宜図面を参照しながら説明する。同様の部分には同一符号を付し、重複する説明を省略する。本発明において、「音像」は、聴取者が音源の存在を知覚する位置を指し、ステレオスピーカ等によって具現化された仮想的な音源の位置であったり、実際の音源の位置であったりする。また、本発明において、「聴覚臨場感」は、聴取者が聴覚において感じる臨場感を指す。従来、臨場感は、ＡＶ機器の性能を語る際の重要なキーワードであるにも関わらず、曖昧な意味で用いられていたのに対し、本発明の聴覚臨場感推定装置は、聴覚臨場感評価値を定量的かつ客観的に推定することができる。
【００２７】
再生側の機器としては、ステレオスピーカ、５．１チャネルシステム、ヘッドホン等があるが、聴取者は、最終的に両耳に到達した音響信号によって聴覚臨場感を知覚する。したがって、本発明の聴覚臨場感評価装置は、聴取者の両耳と同様に、２チャネルの音響信号に基づいて音像の移動量を推定することとする。
【００２８】
図１に示すように、本発明の実施形態に係る聴覚臨場感評価システム１は、再生側の機器として、記憶媒体再生装置１０と、スピーカ群２０と、を備えている。記憶媒体再生装置１０は、図示しない記憶媒体に記憶されたデータを読み取り、ステレオスピーカ、５．１チャネルシステム、ヘッドホン等からなるスピーカ群２０を介して音響信号を再生する。
【００２９】
本発明の実施形態に係る聴覚臨場感評価システム１は、評価側の機器として、一対のマイクロホン３０Ｌ，３０Ｒと、聴覚臨場感評価装置４０と、通知装置５０と、を備える。
【００３０】
＜マイクロホン＞
一対のマイクロホン３０Ｌ，３０Ｒは、左右に配置された２チャネル構造を呈し、スピーカ群２０で再生された音響信号を計測し、計測された２チャネルの音響信号を聴覚臨場感評価装置４０へ出力する。本明細書では、スピーカ群２０で再生された音波だけでなく、当該音波のマイクロホン３０Ｌ，３０Ｒによる計測結果も音響信号と記載する。
【００３１】
＜聴覚臨場感評価装置＞
聴覚臨場感評価装置４０は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read-Only Memory）、入出力回路等から構成されており、機能部として、音響信号分割部４１と、相互相関関数算出部４２と、音圧レベル算出部４３と、動画像生成部４４と、音像移動ベクトル算出部４５と、聴覚臨場感評価値算出部４６と、を備える。
【００３２】
≪音響信号分割部≫
音響信号分割部４１は、一対のマイクロホン３０Ｌ，３０Ｒから出力された２チャネルの音響信号をある一定時間ごとにある区間長だけ取得し、取得された２チャネルの音響信号をＭ個のバンドパスフィルタによって周波数帯域ごとにＭ個に分割し、分割された音響信号を相互相関関数算出部４２及び音圧レベル算出部４３へ出力する。
【００３３】
本実施形態において、音響信号分割部４１は、５個のオクターブバンドパスフィルタを備えており、マイクロホン３０Ｌから出力された音響信号Ｌを、中心周波数が１２５［Ｈｚ］である周波数帯域ｆ１の第一の音響信号Ｌ１、中心周波数が２５０［Ｈｚ］である周波数帯域ｆ２の第二の音響信号Ｌ２、中心周波数が５００［Ｈｚ］である周波数帯域ｆ３の第三の音響信号Ｌ３、中心周波数が１０００［Ｈｚ］である周波数帯域ｆ４の第四の音響信号Ｌ４、中心周波数が２０００［Ｈｚ］である周波数帯域ｆ５の第五の音響信号Ｌ５に分割する。
【００３４】
同様に、音響信号分割部４１は、マイクロホン３０Ｒから出力された音響信号を、中心周波数が１２５［Ｈｚ］である周波数帯域ｆ１の第一の音響信号Ｒ１、中心周波数が２５０［Ｈｚ］である周波数帯域ｆ２の第二の音響信号Ｒ２、中心周波数が５００［Ｈｚ］である周波数帯域ｆ３の第三の音響信号Ｒ３、中心周波数が１０００［Ｈｚ］である周波数帯域ｆ４の第四の音響信号Ｒ４、中心周波数が２０００［Ｈｚ］である周波数帯域ｆ５の第五の音響信号Ｒ５に分割する。
【００３５】
なお、音響信号の分割数は、５個に限定されず、対象とする音像の種類に応じて適宜変更可能である。例えば、自動車騒音のように広帯域な雑音成分を含む音を対象とする場合には、可聴周波数帯域（２０［Ｈｚ］〜２０［ｋＨｚ］）を１／１オクターブバンド幅で区切ったＭ＝１０程度とすることができる。また、楽器音のように高調波構造を有する音を対象とする場合には、１／１０オクターブバンド幅で区切ったＭ＝４０程度とすることができる。また、汎用的なモデルのバンドパスフィルタとして、５０チャネル程度の聴覚フィルタ（聴覚末梢系の特性を考慮したバンドパスフィルタであり、代表的なものにＲＯＥＸフィルタ、Ｇａｍｍａ−ｔｏｎｅフィルタ等がある）を採用することによって、聴覚特性を反映した分割を行うことができる。
【００３６】
≪相互相関関数算出部≫
相互相関関数算出部４２は、音響信号分割部４１から出力された各音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５を取得し、取得された各音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５に基づいて、周波数帯域ごとの相互相関関数を算出し、算出された相互相関関数を動画像生成部４４へ出力する。
【００３７】
本実施形態において、相互相関関数算出部４２は、算出周期ごと（時刻ｔ１，ｔ２，…）に、音響信号Ｌ１，Ｒ１の相互相関（ＣＣ：Cross Correlation）関数を算出する。
相互相関関数算出部４２は、図２（ａ）に示すように、時刻ｔ１を先頭として区切った区間長（時間窓幅ともいう）Ｔ１の音響信号Ｌ１，Ｒ１を抽出して抽出された音響信号Ｌ１，Ｒ１の相関値を算出するとともに、音響信号Ｒ１のみをラグｍずつずらして区切って抽出し直して相関値を算出することを繰り返すことによって、図３（ａ）の下段に示すように、ラグｍを変数とする相互相関関数を算出する。相互相関関数は、−１から＋１までの値をとり、音響信号Ｌ１，Ｒ１が一致した場合に＋１となり、音響信号Ｌ１，Ｒ１の波形が異なるほど値が小さくなる。本発明において、ラグｍは、音響信号Ｌ１，Ｒ１における相互相関関数を算出するための区間のずれ（時間差）である。
【００３８】
前記した時刻ｔ１，ｔ２，…を先頭として区切った区間長Ｔ１、ｐ番目の算出周期、ｑ番目の周波数帯域における時刻ｎの関数である音響信号Ｌｐ，ｑ（ｎ），Ｒｐ，ｑ（ｎ）に対する相互相関関数ｒｐ，ｑ（ｍ）は、下記式によって算出される。
【００３９】
【数１】
0005658588.tif 000002
【００４０】
ここで、σＬｐ，ｑ，σＲｐ，ｑは、それぞれＬｐ，ｑ（ｎ），Ｒｐ，ｑ（ｎ）の標準偏差である。また、Ｎは、相互相関関数を算出する期間の長さＴに対応する離散時間である。
【００４１】
なお、区間長は、対象とする音像の種類に応じて適宜変更可能である。例えば、定常的な騒音を対象とする場合には、１００〜２００［ｍｓ］程度とし、時間的な変化が大きい楽器音を対象とする場合には、１０〜２０［ｍｓ］程度とすることができる。かかる区間長は、相互相関関数算出部４２に予め記憶されていてもよく、聴覚臨場感評価装置４０に接続されたキーボード、マウス等からなる入力装置（図示せず）をユーザが操作することによって設定される構成であってもよい。
【００４２】
また、算出周期（時刻ｔ１，ｔ２，…の間隔）と、区間長とは、一致しても一致しなくてもよい。
【００４３】
また、相互相関関数算出部４２は、聴覚臨場感評価値算出部４６によって算出された音像の移動量（移動ベクトルの大きさ）に基づいて、算出周期を変更する構成であってもよい。例えば、音像の移動量が小さい場合には、算出周期を長く設定することによって演算量を削減することができる。この場合において、二次元動画像の一般的なフレームレートである１／３０［ｓ］を算出周期として相互相関関数を算出する場合であっても、区間長を１／３０［ｓ］よりも長く設定し、音響信号Ｒ１の抽出範囲をオーバーラップさせながら相互相関関数を算出することができる。相互相関関数算出部４２には、音像の移動量と、相互相関関数の算出周期と、の関係がデータベースとして予め記憶されており、相互相関関数算出部４２は、聴覚臨場感評価値算出部４６によって算出された音像の移動量を用いてデータベースを検索することによって、音像の移動量の対応する算出周期を読み出し、読み出された算出周期に基づいて相互相関関数を算出する。
【００４４】
また、相互相関関数算出部４２は、区間長の異なる複数の区間長別相互相関関数を算出し、算出された複数の区間長別相互相関関数に基づいて相互相関関数を算出する構成であってもよい。例えば、相互相関関数算出部４２には、区間長Ｔ１，Ｔ２，Ｔ３が予め記憶されており、相互相関関数算出部４２は、図２（ｂ）に示すように、区間長Ｔ１における区間長別相互相関関数を算出するとともに、区間長Ｔ１よりも長い区間長Ｔ２，Ｔ３における区間長別相互相関関数を算出し、算出された３個の区間長別相互相関関数の相加平均を算出することによって、相互相関関数を算出することができる。
【００４５】
相互相関関数算出部４２は、他の周波数帯域ｆ２〜ｆ５に関しても同様に相互相関関数を算出する。
【００４６】
≪音圧レベル算出部≫
音圧レベル算出部４３は、音響信号分割部４１から出力された各音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５を取得し、取得された各音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５に基づいて、周波数帯域ごとの音圧レベル、及び、左右の音響信号の音圧の差である音圧レベル差を算出し、算出された音圧レベル及び音圧レベル差を動画像生成部４４へ出力するとともに、算出された音圧レベルを聴覚臨場感評価値算出部４６へ出力する。
【００４７】
≪動画像生成部≫
動画像生成部４４は、相互相関関数算出部４２から出力された周波数帯域ごとの相互相関関数と、音圧レベル算出部４３から出力された音圧レベル及び音圧レベル差に基づいて、動画像（より詳細には、二次元動画像）を生成し、生成された動画像を音像移動ベクトル算出部４５へ出力する。
【００４８】
図５（ａ）に示すように、動画像は、時刻ｔ１におけるフレーム１００、時刻ｔ２におけるフレーム２００、時刻ｔ３におけるフレーム３００、時刻ｔ４におけるフレーム４００、…から構成されており、動画像の一のフレームは、縦軸を周波数（周波数帯域ｆ１〜ｆ５）、横軸をラグｍとする平面である。本実施形態では、一のフレームは、縦５×横８の画素から構成されており、動画像生成部４４は、相互相関関数を対応する画素の輝度、音圧レベルを対応する画素の彩度、音圧レベル差を対応する画素の色相に変換することによって、動画像を生成する。
【００４９】
例えば、相互相関関数のみを用いて動画像を生成する場合には、動画像生成部４４は、相互相関関数が−１の場合に最低輝度となり、相互相関関数が＋１の場合に最高輝度となるように各画素の画素データを生成し（図３（ｂ）参照）、かかる画素データによって構成された一連の画像からなる動画像を生成する。ここで、動画像生成部４４は、図３（ａ）の下段に示す周波数帯域ｆ１の相互相関関数から、図３（ｃ）の下段に示す輝度を有する画素列１１０を生成し、図３（ａ）の上段に示す周波数帯域ｆ２の相互相関関数から、図３（ｃ）の上段に示す輝度を有する画素列１２０を生成する。
【００５０】
また、動画像生成部４４は、図４（ａ）に示す周波数帯域ごとの音響信号Ｌ１（又はＬ２〜Ｌ５），Ｒ１（又はＲ２〜Ｒ５）の音圧レベル２１０，２２０の音圧レベル差２３０に基づいて、図４（ｂ）に示すラグｍごとの重みを算出し、算出された重みによって重み付けされた輝度を有する画素列１１０’を生成する構成であってもよい。
【００５１】
相関の高い音響信号が左右の耳にある音圧レベル差で到来した場合、左右のレベル差によって、正面を０度として±９０度のいずれかの方向に音像が知覚される。動画像生成部４４は、かかる音圧レベル差と音像の知覚方向との関係を予め記憶しておき、音圧レベル差２３０を用いてかかる関係を参照することによって、何度の方向に音像があるかを推測することができる。動画像生成部４４は、音像があると推定される方向に重みとして１を設定し、そこから遠ざかる方向につれて、重みとして小さい値を設定する。
【００５２】
一方、相互相関関数の値が最大になるラグｍが音像のある方向であるから、ラグｍと音像の方向とは一対一に対応する。したがって、動画像生成部４４は、重み付けによって、音圧レベル差２３０に基づいて音像があると推定された付近の相互相関関数の値はそのままにし、音圧レベル差２３０に基づいて音像がないと推定された付近の相互相関関数の値は小さくなるようにする。こうすることで、正弦波のように相互相関関数が１となるラグが数か所でてくるような場合でも、動画像生成部４４は、音圧レベル差２３０に基づいて音像の方向を判定することができる。
【００５３】
また、動画像生成部４４は、相互相関関数を輝度又は色相に変換して動画像を生成する場合には、相互相関関数がマイナスの部分を０とみなして動画像を生成する構成であってもよい。これは、音像が１個の場合には、相互相関関数が最大となるラグに対応する位置が音像の位置と一致するためである。
【００５４】
≪音像移動ベクトル算出部≫
音像移動ベクトル算出部４５は、動画像生成部４４から出力された動画像の複数のフレームを用いて動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出し、算出された音像の移動ベクトルを聴覚臨場感評価値算出部４６へ出力する。動画像の移動ベクトルの算出手法としては、オプティカルフロー、パターンマッチング、特徴量追従アルゴリズム等といった、画像情報処理のための公知のアルゴリズムを適宜採用することができる。
【００５５】
音像移動ベクトル算出部４５は、３個以上（本実施形態では、時刻ｔ１〜ｔ４に対応する４個）のフレームを蓄積し、蓄積されたフレームの輝度、色相及び彩度に基づいて、周波数帯域ごとの移動ベクトルを算出する。詳細には、音像移動ベクトル算出部４５は、図５（ｂ）に示すように、時刻ｔ１における画素列１１１（ここでは、説明のため、１列の画素数を１２個とする。）、時刻ｔ２における画素列２１１、時刻ｔ３における画素列３１１を用いて、周波数帯域ｆ１における移動ベクトルを算出するとともに、他の周波数帯域ｆ２〜ｆ５についても同様の処理を行う。ここで、時刻ｔ１，ｔ２の２個のフレームを用いて移動ベクトルを算出する場合には、音像移動ベクトル算出部４５は、時刻ｔ１における画素列１１１の左から４個目の画素が、時刻ｔ２においても画素列２１１の左から４個目に留まっていると誤算出するおそれがある（実線の矢印）。そのため、音像移動ベクトル算出部４５は、時刻ｔ３における画素列３１１をさらに用いることによって、時刻ｔ１における画素列１１１の左から１〜３個目の画素が、時刻ｔ２における画素列２１１では左から５〜７個目、時刻ｔ３における画素列では左から１０〜１２個目へと移動していると推定するため、時刻ｔ１における画素列１１１の左から４個目の画素が、時刻ｔ２においては画素列２１１の左から８個目へ移動していると算出することができる（点線の矢印）。このように、３個以上のフレームを用いて移動ベクトルを算出することによって、移動ベクトルの誤算出を防ぐことができる。
【００５６】
また、音像移動ベクトル算出部４５は、周波数特性が変動する音像に対して、周波数帯域ごとの音圧レベルの変化に応じて、周波数帯域を跨ぐように音像の移動ベクトルを算出してもよい。ここで、音像移動ベクトル算出部４５は、音圧レベルの周波数成分の時間変化に基づいて、周波数、時間及び音圧レベルを軸とする三次元空間における音圧レベルのピーク（尾根部分）が周波数の軸上でシフトした場合に、音像の周波数特性が変動していると判定し、周波数帯域を跨ぐように音像の移動ベクトルを算出することができる。
【００５７】
≪聴覚臨場感評価値算出部≫
聴覚臨場感評価値算出部４６は、音像移動ベクトル算出部４５から出力された音像の移動ベクトルに基づいて、聴覚臨場感評価値を算出し、通知装置５０へ出力する。本実施形態において、聴覚臨場感評価値算出部４６は、閾値未満の音像の移動量を切り捨てるとともに、閾値以上の複数の音像の移動量を相加平均することによって、聴覚臨場感評価値を算出することができる。かかる聴覚臨場感評価値は、音像の移動ベクトルが大きいほど、当該聴覚臨場感評価値が大きくなるような値である。
【００５８】
聴覚臨場感評価値算出部４６は、動画像の周波数帯域ごとの移動ベクトル、すなわち、音像の移動ベクトルに基づいて、周波数帯域ごとの音像移動量を推定し、推定された周波数帯域ごとの音像移動量の最大値を採用する、音像移動量の相加平均を算出する、又は、音像移動量の積算値を算出することによって、全体としての音像移動量を推定する。
【００５９】
本実施形態では、聴覚臨場感評価値算出部４６は、音像移動量の積算値を採用する。すなわち、聴覚臨場感評価値算出部４６は、図６（ａ）に示すように、前フレームと今フレームにおける各周波数帯域ｆ１〜ｆ５の移動ベクトルＶ１ａ，Ｖ１ｂ，Ｖ２〜Ｖ５を算出する。
【００６０】
続いて、聴覚臨場感評価値算出部４６は、図６（ｂ）に示すように、各周波数帯域ｆ１〜ｆ５の移動ベクトルを統合する。ここでは、周波数帯域ｆ２の移動ベクトルＶ２と周波数帯域ｆ５の移動ベクトルＶ５が類似しているため統合されるとともに、周波数帯域ｆ１の左の移動ベクトルＶ１ａと周波数帯域ｆ１の左の移動ベクトルＶ１ｂと周波数帯域ｆ４の移動ベクトルＶ４とが類似しているため統合される。
【００６１】
ここで、聴覚臨場感評価値算出部４６は、二つの移動ベクトルの向きが同じで、大きさがほぼ同じ（例えば、移動ベクトルの大きさの比率が一定範囲内）であれば、二つの移動ベクトルが類似していると判定する構成であってもよく、さらに、二つの移動ベクトルの音圧レベル及び音圧レベル差の少なくとも一方がほぼ同じ（比率が一定範囲内）である場合に、二つの移動ベクトルが類似していると判定する構成であってもよい。
【００６２】
続いて、聴覚臨場感評価値算出部４６は、図６（ｃ）に示すように、統合後の移動ベクトルＶ１ｂ（＝Ｖ１ａ＝Ｖ４），Ｖ３，Ｖ５の大きさの和（積分値）を算出し、算出結果を音像移動量とする。
【００６３】
なお、聴覚臨場感評価値算出部４６は、図６（ｄ）に示すように、統合後の移動ベクトルＶ１ｂ，Ｖ３，Ｖ５を音圧によって補正し、補正後の移動ベクトルＶ１ｂ’，Ｖ３’，Ｖ５’の大きさの和（積分値）を算出し、算出結果を音像移動量とする構成であってもよい。例えば、音圧レベルが大きいほど聴覚臨場感評価値に対する音像の移動ベクトルの寄与度が高くなるように重み付けを行うことによって、聴覚臨場感評価値の推定精度を向上させることができる。
【００６４】
続いて、聴覚臨場感評価値算出部４６は、今フレームと次フレームにおける音像移動量を算出し、かかる算出を所定フレーム数に達するまで繰り返し、所定フレーム数分の音像移動量の和（積算値）を算出し、算出結果を最終的な音像移動量、すなわち、聴覚臨場感評価値とする。
【００６５】
また、聴覚臨場感評価値算出部４６は、フレーム数の異なる複数のフレーム数別聴覚臨場感評価値を算出し、算出された複数のフレーム数別聴覚臨場感評価値に基づいて聴覚臨場感評価値を算出する構成であってもよい。
例えば、聴覚臨場感評価値算出部４６は、３個のフレームにおけるフレーム数別聴覚臨場感評価値と、４個のフレームにおけるフレーム数別聴覚臨場感評価値と、を算出し、算出された２個のフレーム数別聴覚臨場感評価値の相加平均を算出することによって、聴覚臨場感評価値を算出することができる。
【００６６】
また、音像が１個の場合には、聴覚臨場感評価値算出部４６は、輝度（すなわち、相互相関関数）が最大となるラグｍに基づいて、ある時刻における音像の方向（方位角）を推定することもできる。これは、音像が１個の場合には、相互相関関数が最大となるラグｍが、音像から各マイクロホン３０Ｌ，３０Ｒまでの距離の差と関連するためである。
【００６７】
また、聴覚臨場感評価値算出部４６は、複数の周波数帯域の音圧レベルが同期して変動している場合には、当該複数の周波数帯域に同じ音像が存在するとみなして音像の移動ベクトルを統合することもできる。かかる構成によると、聴取者が知覚するのに似た音像の移動ベクトルを用いて聴覚臨場感評価値を算出することができる。
【００６８】
＜通知装置＞
通知装置５０は、ディスプレイ、スピーカ等からなり、聴覚臨場感評価値算出部４６から出力された聴覚臨場感評価値を表示したり音声出力したりすることによって、利用者へ通知する。
【００６９】
＜動作例＞
続いて、本発明の実施形態に係る聴覚臨場感評価装置４０の動作例について説明する。まず、音響信号分割部４１が、一対のマイクロホン３０Ｌ，３０Ｒから出力された２チャネルの音響信号を取得し、取得された２チャネルの音響信号を５個のバンドパスフィルタによって周波数帯域が異なる５個の音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５に分割する。
【００７０】
続いて、相互相関関数算出部４２が、各音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５に基づいて、周波数帯域ごとの相互相関関数を算出する。
【００７１】
続いて、音圧レベル算出部４３が、各音響信号Ｌ１〜Ｌ５，Ｒ１〜Ｒ５に基づいて、周波数帯域ごとの音圧レベル、及び、左右の音響信号の音圧の差である音圧レベル差を算出し、算出された音圧レベル及び音圧レベル差を動画像生成部４４へ出力するとともに、算出された音圧レベルを聴覚臨場感評価値算出部４６へ出力する。
【００７２】
続いて、動画像生成部４４が、相互相関関数と、音圧レベル算出部４３から出力された音圧レベル及び音圧レベル差に基づいて、相互相関関数を輝度、音圧レベルを彩度、音圧レベル差を色相とみなした各画素の画素データを生成し、かかる画素データによって構成された一連の画像からなる動画像を生成する。
【００７３】
続いて、音像移動ベクトル算出部４５が、動画像生成部４４から出力された動画像の複数のフレームを用いて動画像の移動ベクトルを算出することによって、音像の移動ベクトルを算出する。
【００７４】
続いて、聴覚臨場感評価値算出部４６が、音像移動ベクトル算出部４５から出力された音像の移動ベクトルに基づいて、聴覚臨場感評価値を算出し、通知装置５０へ出力する。
【００７５】
＜音像の移動ベクトルの算出例＞
続いて、本発明の実施形態に係る聴覚臨場感評価装置４０の音像移動ベクトル算出部４５が、移動する音像の移動ベクトルを実際に算出した例について説明する。
【００７６】
ここで、マイクロホン３０Ｌ，３０Ｒとして、ダミーヘッド（高研社製のSAMRAI）を使用するとともに、スピーカ群２０として、単一のスピーカを使用し、このスピーカを、ダミーヘッドを中心とした半径０．５［ｍ］の円上で、左（−９０［ｄｅｇ］）から右（９０［ｄｅｇ］）まで、又は右から左まで移動させたときの音響信号をダミーヘッドが検出した。また、音響信号分割部４１のサンプリング周波数は４４１００［Ｈｚ］であり、中心周波数１２５，２５０，５００，１０００，２０００，４０００［Ｈｚ］とするオクターブバンドパスフィルタを用いて音響信号分割部４１を構成することにより、音響信号を６個の周波数帯域に分割した。
【００７７】
相互相関関数算出部４２は、相互相関関数の算出周期を１／３０［ｓ］、区間長を１／５［ｓ］として相互相関関数を算出した。ここで、ダミーヘッドにおけるマイクロホン３０Ｌ，３０Ｒ間の距離は０．１８［ｍ］であり、音速は３４０［ｍ／ｓ］であるため、音響信号がマイクロホン３０Ｌ，３０Ｒに到達する時間差が最大となる（スピーカが左又は右にある）場合の時間差は、約０．７［ｍｓ］となる。そこで、ラグｍとしては、０を中心とした±０．７［ｍｓ］に相当する部分にのみ着目した。
【００７８】
また、音像移動ベクトル算出部４５は、相互相関関数の最大値ｒｐ−１，ｑ（ｍ）ｍａｘ，ｒｐ，ｑ（ｍ）ｍａｘを有する画素同士を結ぶ手法を用いて動画像の移動ベクトル、すなわち、音像の移動ベクトルを算出した。
【００７９】
≪単一音像の場合≫
単一のスピーカから音響信号としてホワイトノイズを出力した場合、すなわち、単一の音源を用いて一の音像を生成した場合における、計測開始時点において動画像生成部４４が生成した画像（フレーム）を図７（ａ）に示し、現在（途中経過時点）において動画生成部４４が生成した画像を図７（ｂ）に示す。図７（ｂ）における矢印は、開始時点から現在までの、周波数帯域ごとの音像の移動ベクトルである。図７（ａ）（ｂ）を参照すると、全ての周波数帯域において音像が左から右へ移動している様子が分かる。
【００８０】
≪複数音像の場合≫
左から右へ移動するスピーカからローパスノイズ（遮断周波数７００［Ｈｚ］）を出力した音響信号と、右から左へ移動するスピーカからハイパスノイズ（遮断周波数１０００［Ｈｚ］）を出力した音響信号と、を合成したものを音響信号分割部４１に入力して分析を行った。この場合では、二つの音源を用いて二つの音像を生成している。この場合における、測開始時点において動画像生成部４４が生成した画像（フレーム）を図８（ａ）に示し、現在（途中経過時点）において動画生成部４４が生成した画像を図８（ｂ）に示す。図８（ｂ）における矢印は、開始時点から現在までの、周波数帯域ごとの音像の移動ベクトルである。図８（ａ）（ｂ）を参照すると、低周波数帯域においては音像が左から右へ移動し、高周波数帯域においては音像が右から左へ移動している様子が分かる。
【００８１】
本発明の実施形態に係る聴覚臨場感評価装置４０は、複数の音像に対して処理を行う場合であっても、音像の方向を用いずに音像の移動量に基づいて聴覚臨場感評価値を算出するので、簡易な処理で聴覚臨場感を客観的に評価することができる。
【００８２】
以上、本発明の実施形態について実施形態を参照して説明したが、本発明は前記実施形態に限定されず、本発明の要旨を逸脱しない範囲で適宜設計変更可能である。例えば、音圧レベル算出部４３を備えずに、動画像生成部４４が相互相関関数を輝度と変換したモノクロの動画像を生成する構成であってもよい。また、動画像生成部４４が相互相関関数を輝度、音圧レベルを彩度としたカラーの動画像を生成する構成であってもよい。この場合には、動画像生成部４４は、各画素の色相を所定値（例えば、色相値５０％）に設定することができる。また、相互相関関数、音圧レベル及び音圧レベル差を、輝度、彩度及び色相のいずれかに変換して動画像を生成する構成であってもよい。
【００８３】
また、聴覚臨場感評価装置４０は、記憶媒体に記憶されていない音響信号を再生する装置、再生装置ではない実際の音源に対しても、聴覚臨場感評価値を算出することができる。また、音像移動ベクトル算出部４５は、複数の周波数帯域の画素列を用いて一つの移動ベクトルを算出することによって、統合された移動ベクトルを算出する構成であってもよく、移動ベクトルごとの画像の特徴量をパターンとして予め記憶しており、パターンマッチングを用いて移動ベクトルを算出する構成であってもよい。また、本発明は、コンピュータを聴覚臨場感評価装置４０として機能させる聴覚臨場感評価プログラムとしても具現化可能である。
【符号の説明】
【００８４】
１      聴覚臨場感評価システム
３０Ｌ，３０Ｒ  マイクロホン
４０    聴覚臨場感評価装置
４１    音響信号分割部
４２    相互相関関数算出部
４３    音圧レベル算出部
４４    動画像生成部
４５    音像移動ベクトル算出部
４６    聴覚臨場感評価値算出部
【図１】
0005658588.tif 000003
【図２】
0005658588.tif 000004
【図６】
0005658588.tif 000005
【図３】
0005658588.tif 000006
【図４】
0005658588.tif 000007
【図５】
0005658588.tif 000008
【図７】
0005658588.tif 000009
【図８】
0005658588.tif 000010