JP6650245B2 - インパルス応答生成装置及びプログラム - Google Patents

Description

本発明は、インパルス応答生成装置及びプログラムに関する。

テレビ番組や音楽、映画などの音響制作において、音の広がり感や臨場感を与えるために、残響を付加する技術が用いられている。空間の残響をシミュレーションで模擬したＩＩＲ（Infinite impulse response、無限インパルス応答）フィルタ型のものや、実空間で測定したインパルス応答を用いて残響を生成するＦＩＲ（Finite impulse response、有限インパルス応答）フィルタ型のものなど、様々な残響付加装置が開発されている（例えば、特許文献１参照）。特に、ＦＩＲフィルタ型の残響付加装置では、インパルス応答を測定した空間の残響特性が忠実に再現できるという特長がある。

ところで、５．１チャネルサラウンド方式やスーパーハイビジョン用の２２．２マルチチャンネル方式といったマルチチャンネル音響方式は、平面上もしくは空間内に配置した複数のスピーカを用いて、様々な方向から到来する音を表現することができるという特長を持っている（例えば、非特許文献１参照）。一般的にこれらのコンテンツを制作する過程でも、残響付加装置を用いて豊かな臨場感と空間的な広がり感を付与することが多い。実空間において壁や床、天井など様々な方向から反射音が到来するように、マルチチャンネル音響方式における残響付加でも、ある任意の方向から到来する特徴的な反射音（方向別の残響）を、その方向に設置したスピーカから再生することで、インパルス応答を測定した空間の残響特性をできるだけ忠実に再現するのが理想である。特に直接音の後、１００ｍｓ（ミリ秒）くらいまでの間に到来する初期反射音は、音像の幅や方向感などの知覚に寄与することが分かっており、残響付加において初期反射音を方向別に再生することは、インパルス応答を測定した空間の印象を再現する上で重要である。そこで、指向性マイクロホンで測定したインパルス応答を用いて、方向別のインパルス応答を生成する技術が提案されてきた（例えば、特許文献２、３参照）。

一方、近接四点法（例えば、非特許文献２参照）やそれを拡張した方法（例えば、非特許文献３、４参照）など、複数のマイクロホンを用いて、音源や仮想音源（反射音）の位置を推定する様々な技術が提案されてきた。これらの手法では、各マイクロホンの相互相関やインテンシティ算出することにより、受音点における反射音の到来時刻、到来方向を予測することができる。

特開２０１４−４５２８２号公報 特開２０１３−２３８６４３号公報 特開２０１４−１１２７６７号公報

W.woszczyk，et al.，"Space Builder: An impulse response-based tool for immersive 22.2 channel ambiance design",日本，AES 40th international conference，2010年 山崎芳男，「近接4点法による音場の空間情報の可視化」，テレビジョン学会誌，1990年，Vol.44，No.3，p.253-258 西村昌浩,外2名,「BoSCマイクロホンによる仮想音源推定のための近接４点法の拡張アルゴリズムの提案」，2012年3月，日本音響学会講演論文集、p.789-790 小林万理恵他，外２名，「分散配置した近接４点法マイクロホンによる音源位置推定と音源分離」，日本音響学会講演論文集，2012年3月，p.685-686

上述した技術によって生成したインパルス応答には、厳密にはある任意の方向以外から到来する反射音も含まれている。指向性マイクロホンとは言っても、その指向性には幅があるため、近接する複数のマイクロホンで同一の反射音を拾ってしまうからである。このため、これらのインパルス応答を用いて残響付加を行う場合、ある反射音が複数のスピーカから再生され、リスニングポイントでの反射音の密度が高くなってしまい、測定音場の空間印象が厳密には再現されないという問題点がある。

本発明は、このような事情を考慮してなされたもので、残響付加において、測定音場の空間印象をより厳密に再現できる方向別のインパルス応答を生成するインパルス応答生成装置及びプログラムを提供する。

本発明の一態様は、同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出された各インパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源の情報を推定する仮想音源情報推定部と、前記仮想音源情報推定部により推定された前記仮想音源の情報と、前記インパルス応答との比較に基づいて、前記インパルス応答において、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音を判定する比較判定部と、前記インパルス応答から、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来すると前記比較判定部により判定された反射音以外の反射音を間引いて方向別のインパルス応答を生成する反射音除去部と、を備えることを特徴とするインパルス応答生成装置である。

本発明の一態様は、上述するインパルス応答生成装置であって、前記仮想音源の情報は、前記任意の位置における前記仮想音源の到来時刻及び到来方向を示し、前記比較判定部は、前記仮想音源の到来時刻と、前記仮想音源の到来方向に対応したマイクロホンにより収音した信号から算出された前記インパルス応答の波形構造が表す経時的な音の大きさの変化との比較に基づいて、前記インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音を判定する、ことを特徴とする。

本発明の一態様は、上述するインパルス応答生成装置であって、前記比較判定部は、複数の前記マイクロホンのうち前記仮想音源の到来方向に対応したマイクロホンを選択し、選択した前記マイクロホンにより収音した信号から算出された前記インパルス応答において、選択した前記マイクロホンと前記任意の位置との距離に応じた遅延を前記仮想音源の到来時刻から減算した時刻を含んだ所定の時刻範囲内で生じた局所的なピークを特定する処理を前記仮想音源のそれぞれについて行い、前記インパルス応答において前記処理により特定された局所的なピークを、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音と判定する、ことを特徴とする。

本発明の一態様は、コンピュータを、同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出された各インパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源の情報を推定する仮想音源情報推定手段と、前記仮想音源情報推定手段により推定された前記仮想音源の情報と、前記インパルス応答との比較に基づいて、前記インパルス応答において、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音を判定する比較判定手段と、前記インパルス応答から、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来すると前記比較判定手段により判定された反射音以外の反射音を間引いて方向別のインパルス応答を生成する反射音除去手段と、を具備するインパルス応答生成装置として機能させるためのプログラムである。

本発明によれば、残響付加において、測定音場の空間印象をより厳密に再現できる方向別のインパルス応答を生成することが可能となる。

単一指向性マイクロホンの指向特性を示す図である。 単一指向性マイクロホンにより測定されたインパルス応答を示す図である。 本発明の一実施形態によるインパルス応答生成装置の構成を示す機能ブロック図である。 同実施形態によるインパルス応答測定方法の例を示す図である。 同実施形態によるインパルス応答生成装置における方向別インパルス応答生成処理のフロー図である。 同実施形態によるインパルス応答測定方法の例を示す図である。 同実施形態によるインテンシティを用いた仮想音源情報の推定処理を説明するための図である。 同実施形態による相互相関係数を用いた仮想音源情報推定処理のフロー図である。 同実施形態による比較判定処理及び反射音除去処理の例を示す図である。

以下、図面を参照しながら本発明の実施形態を詳細に説明する。
まず、図１及び図２を用いて、複数の単一指向性マイクロホンによる収音を説明する。
図１は、単一指向性マイクロホンの指向特性を示す図である。マイクロホンＭ_ａ及びＭ_ｂは、単一指向性マイクロホンである。同図では、方向別の反射音を収音するため、マイクロホンＭ_ａをＤ１方向に向けて、マイクロホンＭ_ｂをＤ２方向に向けて配置している。マイクロホンＭ_ａの指向特性はＲ_ａであり、マイクロホンＭ_ｂの指向特性はＲ_ｂである。点線は、マイクロホンＭ_ａ、_Ｍｂが全指向性マイクロホンであった場合の収音範囲を示す。Ｄ１方向からは反射音Ａが、Ｄ２方向からは反射音Ｂが到来する。

図２は、図１に示す単一指向性マイクロホンにより測定されたインパルス応答を示す図である。インパルス応答は、ＴＳＰ（Time Stretched Pulse）信号などの音響測定用の信号を放射し、マイクロホンにより収音（測定）された信号を処理することによって算出されるものであり、音の大きさの経時的な変化を表す。本明細書では、あるマイクロホンにより収音された信号を処理することによって算出されたインパルス応答を、そのマイクロホンにより測定されたインパルス応答と記載する。具体的には、ある任意の位置のマイクロホンにより測定されたインパルス応答は、時刻ｔと、時刻ｔにおいてそのマイクロホンが収音した音の振幅とを対応付けた情報として得られる。図２（ａ）は、図１に示すマイクロホンＭ_ａにより測定されたインパルス応答ｈ_ａ（ｔ）を示し、図２（ｂ）は、図１に示すマイクロホンＭ_ｂにより測定されたインパルス応答ｈ_ｂ（ｔ）を示す。

本来、Ｄ１方向に向けて設置したマイクロホンＭ_ａのみでＤ１方向からの反射音Ａを収音し、Ｄ２方向に向けて設置したマイクロホンＭ_ｂのみでＤ２方向からの反射音Ｂを収音するのが理想である。しかし、図１に示すように、実際には指向性マイクロホンの指向特性には幅がある。そのため、反射音Ａ、ＢともマイクロホンＭ_ａ、Ｍ_ｂの両方で収音され、図２に示すインパルス応答ｈ_ａ（ｔ）、ｈ_ｂ（ｔ）のような波形となってしまう。つまり、インパルス応答ｈ_ａ（ｔ）では、反射音Ａより反射音Ｂの方がレベルは小さく、インパルス応答ｈ_ｂ（ｔ）では、反射音Ｂより反射音Ａの方がレベルは小さいものの、それぞれのマイクロホンに対応する方向（測定したい方向）以外から到来する反射音も含まれてしまう。

そのため、複数の指向性マイクロホンのそれぞれにより測定されたインパルス応答を使って残響付加を行うと、各マイクロホンに対応する方向（測定したい方向）以外から到来する反射音が複数のスピーカから再生されてしまい、測定空間の印象が変わってしまう可能性がある。そこで、本実施形態のインパルス応答生成装置は、マイクロホンＭ_ａにより収音された反射音Ｂや、マイクロホンＭ_ｂにより収音された反射音Ａのような、マイクロホンに対応する方向（測定したい方向）以外から到来する反射音を、測定されたインパルス応答から除去し、より正確な方向別のインパルス応答を生成する。

方向別のインパルス応答を生成するため、本実施形態では、まず、複数のマイクロホンを空間内の任意の位置に配置し、インパルス応答を測定する。インパルス応答生成装置は、各マイクロホンにより測定されたインパルス応答を用いて、マイクロホン群に囲まれた任意の位置における、空間内の仮想音源（反射音）の到来方向及び到来時刻を推定する。インパルス応答生成装置は、推定したそれらの情報と、各マイクロホンにより測定されたインパルス応答に含まれる反射音の到来時刻とを照らし合わせる。インパルス応答生成装置は、照らし合わせの結果に基づき、各マイクロホンに対応する方向（測定したい方向）以外から到来する反射音を、そのマイクロホンにより測定されたインパルス応答からの除去対象として特定する。インパルス応答生成装置は、各マイクロホンにより測定されたインパルス応答から、除去対象の任意の反射音を除去し、方向別のインパルス応答を生成する。なお、各マイクロホンにより収音されるべき到来方向は、そのマイクロホンの設置方向に加え、隣接するマイクロホンの設置方向等によっても異なり得る。

図３は、本発明の一実施形態によるインパルス応答生成装置１の構成を示す機能ブロック図であり、本実施形態と関係する機能ブロックのみを抽出して示してある。インパルス応答生成装置１は、例えば、コンピュータ装置により実現することができる。同図に示すように、インパルス応答生成装置１は、記憶部１１と、測定情報取得部１２と、仮想音源情報推定部１３と、比較判定部１４と、反射音除去部１５とを備えて構成される。

記憶部１１は、各マイクロホンにより測定されたインパルス応答など各種情報を記憶する。
測定情報取得部１２は、空間内の任意の位置に配置した複数のマイクロホンのそれぞれにより測定されたインパルス応答を取得し、記憶部１１に書き込む。以下では、ｎ本のマイクロホンのうちｉ番目のマイクロホンをＭ_ｉ（ｉは１以上ｎ以下の整数）と記載し、マイクロホンＭ_ｉにより測定されたインパルス応答をｈ_ｉ（ｔ）と記載する。また、インパルス応答ｈ_ｉ（ｔ）に含まれる局所的なピーク（包絡の極大点）における時刻を、時刻順にｔ_ｉ１，ｔ_ｉ２，…とする。

仮想音源情報推定部１３は、各マイクロホンにより測定されたインパルス応答からインテンシティや相互相関係数などを算出することにより、マイクロホン群（マイクロホンＭ_１〜Ｍ_ｎ）で囲まれた任意の点における、空間内の仮想音源の情報を推定する。推定の対象は、方向知覚への寄与が大きいと考えられるレベルの仮想音源である。方向知覚への寄与が大きいと考えられるレベルの仮想音源は、主に初期反射音部分、直接音から１００ｍｓ（ミリ秒）くらいまでの間に受音点に到来する仮想音源である。
比較判定部１４は、仮想音源情報推定部１３が推定した仮想音源の情報と元のインパルス応答の波形構造とを比較する。比較判定部１４は、残響付加の際に、ある反射音が複数のスピーカから重複して再生されないよう、比較結果に基づき、元のインパルス応答から除去すべき反射音（局所的なピーク）を判定する。
反射音除去部１５は、元のインパルス応答から、比較判定部１４が判定した反射音（局所的なピーク）を除去し、方向別のインパルス応答を生成する。この生成された方向別のインパルス応答を用いて、従来技術により残響付加が行われる。

図４は、インパルス応答測定方法の例を示す図である。インパルス応答の測定空間内に、インパルス応答測定用信号を再生するスピーカと複数のマイクロホンＭ_０〜Ｍ_ｎとを配置する。ここでは、マイクロホンＭ_０〜Ｍ_ｎは、単一指向性マイクロホンであるが、全指向性マイクロホンを使用してもよい。スピーカの位置や、マイクロホンの個数及び位置などは、測定空間の広さや残響の使用用途、再生方法に応じて設定する。例えば、２２．２マルチチャンネル方式の場合は、２４本もしくはＬＦＥ（Low Frequency Effect）チャンネルを除いた２２本のマイクロホンを使用することが考えられる。そして、適当な位置に設置したスピーカから、ＴＳＰ信号などのインパルス応答測定用の音響信号を放射し、マイクロホンＭ_０〜Ｍ_ｎからなるマイクロホン群で録音する。これにより各マイクロホンＭ_０、Ｍ_１、Ｍ_２、…、Ｍ_ｎのそれぞれにより測定されたインパルス応答ｈ_０（ｔ）、ｈ_１（ｔ）、ｈ_２（ｔ）、…、ｈ_ｎ（ｔ）が得られる。そして、マイクロホン群で囲まれた空間内の任意の点をＰとする。図４（ａ）では、マイクロホンＭ_０の設置位置が点Ｐであり、図４（ｂ）では、マイクロホンＭ_０〜Ｍ_ｎのいずれも設置されていない位置が点Ｐである。

図５は、インパルス応答生成装置１における方向別インパルス応答生成処理のフロー図である。
まず、インパルス応答生成装置１の測定情報取得部１２は、各マイクロホンＭ_０、Ｍ_１、Ｍ_２、…、Ｍ_ｎにより測定されたインパルス応答ｈ_０（ｔ）、ｈ_１（ｔ）、ｈ_２（ｔ）、…、ｈ_ｎ（ｔ）を取得する（ステップＳ１０）。例えば、測定情報取得部１２は、インパルス応答を各マイクロホンＭ_０〜Ｍ_ｎから受信してもよく、他のコンピュータ装置から受信してもよく、記録媒体から読み出してもよい。測定情報取得部１２は、取得したインパルス応答を記憶部１１に書き込む。

仮想音源情報推定部１３は、記憶部１１から読み出したインパルス応答ｈ_０（ｔ）、ｈ_１（ｔ）、ｈ_２（ｔ）、…、ｈ_ｎ（ｔ）を元に、マイクロホン群で囲まれた任意の点Ｐにおける、空間内の仮想音源の情報ｅ_１（ｒ_１,θ_１，φ_１，ｔ_１）、ｅ_２（ｒ_２,θ_２，φ_２，ｔ_２）、…、ｅ_ｍ（ｒ_ｍ,θ_ｍ，φ_ｍ，ｔ_ｍ）を予測する（ステップＳ２０）。ここで、ｒは振幅、θは方位角、φは仰角、ｔは到来時刻、ｍは予測された仮想音源の個数であり、ｒ、θ、φ及びｔの添え字は、ｍ個の仮想音源のうち何番目の仮想音源であるかを示す。振幅ｒは仮想音源のレベルを表し、方位角θ及び仰角φは点Ｐにおける仮想音源の到来方向を表す。方位角θ及び仰角φを表すために用いられる３次元座標空間は任意に設定し得る。ここでは、点Ｐの位置を原点０とし、ｘ軸をスピーカと点Ｐとを結ぶ直線、ｙ軸を地面と平行な平面上でｘ軸と垂直な方向、ｚ軸をｘｙ平面と垂直な方向とする３次元座標空間を用いる。また、ここでは、点ＰにマイクロホンＭ_０を設置した場合を考える。仮想音源情報推定部１３は、近接四点法やそれを拡張した方法など任意の既存の手法を用いてインテンシティや各チャンネル間の相互相関を算出することにより、仮想音源情報を予測する。算出手順の詳細については後述する。仮想音源情報推定部１３は、予測した仮想音源情報ｅ_１〜ｅ_ｍを記憶部１１に書き込む。

次に、比較判定部１４は、各インパルス応答ｈ_１（ｔ）〜ｈ_ｎ（ｔ）から除去すべき反射音を判定する比較判定処理を行う（ステップＳ３０）。まず、比較判定部１４は、記憶部１１からインパルス応答ｈ_０（ｔ）、ｈ_１（ｔ）、ｈ_２（ｔ）、…、ｈ_ｎ（ｔ）及び仮想音源情報ｅ_１（ｒ_１,θ_１，φ_１，ｔ_１）、ｅ_２（ｒ_２,θ_２，φ_２，ｔ_２）、…、ｅ_ｍ（ｒ_ｍ,θ_ｍ，φ_ｍ，ｔ_ｍ）を読み出す。比較判定部１４は、仮想音源情報推定部１３が推定した、点Ｐにおける仮想音源の情報ｅ_１〜ｅ_ｍから、各仮想音源の到来時刻（ｔ_１，ｔ_２，…，ｔ_ｍ）を抽出する。比較判定部１４は、各仮想音源の到来時刻（ｔ_１，ｔ_２，…，ｔ_ｍ）と、測定されたインパルス応答ｈ_１（ｔ）、ｈ_２（ｔ）、…、ｈ_ｎ（ｔ）のそれぞれに含まれる局所的なピークの到来時刻（ｔ_１１，ｔ_１２，…）、（ｔ_２１，ｔ_２２，…）、…、（ｔ_ｎ１，ｔ_ｎ２，…）とを比較する。比較判定部１４は、比較結果に基づいて、各インパルス応答から取り除くべき反射音を判定する。取り除くべき反射音とは、あるマイクロホンにより測定されたインパルス応答に含まれる反射音のうち、そのマイクロホンに対応する方向（マイクロホンを向けた方向）とは異なる方向から到来する反射音である。つまり、この取り除くべき反射音とは、その到来方向に向けて配置したマイクロホンだけではなく、その到来方向以外に向けて配置した他のマイクロホンにより測定されたインパルス応答において、重複して測定された反射音である。比較判定部１４は、取り除くべき反射音の判定にあたって、点Ｐから各マイクロホンＭ_１、Ｍ_２、…、Ｍ_ｎまでの距離に応じた遅延τ_１、τ_２、…、τ_ｎも考慮する。さらに、判定には、測定誤差等も加味する必要がある。そこで、比較判定部１４は、ある所定の窓幅ｗを設定する。比較判定部１４は、遅延を考慮した仮想音源の到来時刻と、インパルス応答においてその到来時刻を中心とした窓幅ｗの時刻範囲内で局所的なピークが生じた時刻とを比較して判定を行う。

反射音除去部１５は、各インパルス応答ｈ_１（ｔ）〜ｈ_ｎ（ｔ）から比較判定部１４が除去すべきと判定した反射音を除去する反射音除去処理を行い、方向別のインパルス応答を生成する（ステップＳ４０）。まず、反射音除去部１５は、記憶部１１から読み出したインパルス応答ｈ_１（ｔ）、ｈ_２（ｔ）、…、ｈ_ｎ（ｔ）のそれぞれから、比較判定部１４によりそのインパルス応答から除去すべきと判定された反射音を除去し、方向別のインパルス応答ｈ’_１（ｔ）、ｈ’_２（ｔ）、…、ｈ’_ｎ（ｔ）を生成する。

以下に、具体例を交えながら、図５のステップＳ２０〜ステップＳ４０の詳細な処理を説明する。
本来、三次元空間内の仮想音源情報を得るためには、同一平面内にない４点以上にマイクロホンを設置する必要があるが、ここでは原理を分かりやすく説明するため、同一平面内に複数のマイクロホンを置いて測定した場合を考える。

図６は、インパルス応答測定方法の例を示す図である。マイクロホンＭ_０を中心とし、それを取り囲むようにマイクロホンＭ_１〜Ｍ_４を、マイクロホンＭ_０から等距離かつ等間隔に設置する。隣り合うマイクロホンＭ_１〜Ｍ_４とマイクロホンＭ_０とがなす角度は９０度であり、マイクロホンＭ_０と各マイクロホンＭ_１〜Ｍ_４との距離はＬである。具体的には、マイクロホンＭ_１、Ｍ_２、Ｍ_３、Ｍ_４は、マイクロホンＭ_０を原点とした２次元空間のｘｙ平面上に設置される。マイクロホンＭ_０、Ｍ_１、Ｍ_２、Ｍ_３、Ｍ_４それぞれの設置方向は、ｘ軸に対してそれぞれ０度、０度、９０度、２７０度、１８０度である。

そして、スピーカから測定信号を再生し、各マイクロホンＭ_０〜Ｍ_４により収音された信号からインパルス応答ｈ_０（ｔ）〜ｈ_４（ｔ）を算出する。同図に示すように、各インパルス応答ｈ_ｉ（ｔ）には（ｉは１以上４以下の整数）、マイクロホンＭ_ｉに対応する方向（マイクロホンを向けた方向）から到来する反射音のレベル（実線）だけでなく、それ以外の方向から到来する反射音のレベル（破線）も含まれてしまう。指向性マイクロホンで測定すれば、全指向性マイクロホンで測定するよりも、マイクロホンに対応する方向（マイクロホンを向けた方向）以外から到来する反射音を抑制することができる。しかしながら、指向性マイクロホンとは言っても、その指向性には幅があるため、隣接する複数のマイクロホンで同じ反射音を拾ってしまう。そこで、インパルス応答生成装置１は、後述する処理によって、マイクロホン群の中心点Ｍ_０（マイクロホンＭ_０の位置）における仮想音源情報を推定し、その推定した仮想音源情報を元に、各マイクロホンＭ_１〜Ｍ_４のインパルス応答ｈ_１（ｔ）〜ｈ_４（ｔ）から重複する反射音を間引く。

次に、ステップＳ２０において仮想音源情報推定部１３が、測定されたインパルス応答から仮想音源情報を推定する処理を説明する。推定には様々な手法が用いられるが、ここでは例として、インテンシティを用いる手法と、相互相関係数を用いる手法とを示す。

まず、インテンシティを用いた仮想音源情報の推定処理を説明する。
図７は、インテンシティを用いた仮想音源情報の推定処理を説明するための図である。
最初に、仮想音源情報推定部１３は、インパルス応答ｈ_０（ｔ）及びインパルス応答ｈ_１（ｔ）から０−Ｘ方向（Ｍ_０−Ｍ_１の方向）の瞬時インテンシティＩ_ｘ（ｔ）を求める。さらに、仮想音源情報推定部１３は、インパルス応答ｈ_０（ｔ）及びインパルス応答ｈ_２（ｔ）とから０−Ｙ方向（Ｍ_０−Ｍ_２の方向）の瞬時インテンシティＩ_ｙ（ｔ）を求める。

瞬時インテンシティＩ_ｘ（ｔ）は、以下の式（１）により求められる。

式（１）において、ｐ_１（ｔ）、ｐ_２（ｔ）は、マイクロホンで測定された時刻ｔにおける音圧であり、それぞれ、インパルス応答ｈ_０（ｔ）、ｈ_１（ｔ）を用いる。また、ｄは、マイクロホン間の間隔であり、マイクロホンＭ_０とマイクロホンＭ_１の距離Ｌを用いる。ρ_０は、媒質の密度である。瞬時インテンシティＩ_ｘ（ｔ）の値の正負は、ｘ軸方向の到来方向を表す。
瞬時インテンシティＩ_ｙ（ｔ）も同様に求められる。ｐ_１（ｔ）、ｐ_２（ｔ）には、インパルス応答ｈ_０（ｔ）、ｈ_２（ｔ）を用い、ｄには、マイクロホンＭ_０とマイクロホンＭ_２の距離Ｌを用いる。瞬時インテンシティＩ_ｙ（ｔ）の値の正負は、ｙ軸方向の到来方向を表す。

なお、瞬時インテンシティＩ_ｘ（ｔ）は、ｘ軸上に設置した任意のマイクロホンにより測定されたインパルス応答の組合せ（例えば、インパルス応答ｈ_０（ｔ）及びｈ_４（ｔ））を用いて算出することができる。同様に、瞬時インテンシティＩ_ｙ（ｔ）は、ｙ軸上に設置した任意のマイクロホンにより測定されたインパルス応答の組合せ（例えば、インパルス応答ｈ_０（ｔ）及びｈ_３（ｔ））を用いて算出することができる。

続いて、仮想音源情報推定部１３は、ヒルベルト変換により瞬時インテンシティＩ_ｘ（ｔ）、瞬時インテンシティＩ_ｙ（ｔ）それぞれの包絡を取ったエンベロープインテンシティＩｈ_ｘ（ｔ）、Ｉｈ_ｙ（ｔ）を求める。仮想音源情報推定部１３は、エンベロープインテンシティＩｈ_ｘ（ｔ）、Ｉｈ_ｙ（ｔ）のレベルをベクトル合成し、合成インテンシティのレベルＬｅｖ＿Ｉｈ（ｔ）を求める。エンベロープインテンシティＩｈ_ｘ（ｔ）、Ｉｈ_ｙ（ｔ）のレベルをベクトル合成することにより、時刻ｔにおける仮想音源の到来方向も得られるが、同図では、レベルの経時的な変化のみを示している。仮想音源情報推定部１３は、合成インテンシティレベルＬｅｖ＿Ｉｈ（ｔ）の最初のピークレベルＬｅｖ＿Ｐｓ（直接音）から任意のレベルＬｅｖ＿Ｐｅまでの範囲に表れる局所的なピークＰ_１、Ｐ_２、…、Ｐ_ｍを検出する。仮想音源情報推定部１３は、検出したピークＰ_１、Ｐ_２、…、Ｐ_ｍにおける合成インテンシティから、室内（測定空間）における仮想音源の情報ｅ_１（ｒ_１,θ_１，φ_１，ｔ_１）、ｅ_２（ｒ_２,θ_２，φ_２，ｔ_２）、…、ｅ_ｍ（ｒ_ｍ,θ_ｍ，φ_ｍ，ｔ_ｍ）を得る。なお、到来時刻ｔ_１＜ｔ_２＜ｔ_３＜…＜ｔ_ｍとする。仮想音源の情報ｅ_１〜ｅ_ｍを到来時刻によって並べた情報をｅ（ｔ）とする。

上記は、マイクロホンＭ_０〜Ｍ_４が同一平面（ｘｙ平面）上に設置されている例であるため、仮想音源のＺ座標は０であり、仰角φ_１、φ_２、…、φ_ｍも０である。ｚ軸上にもマイクロホンを設置した場合、つまり三次元的にマイクロホンを設置した場合は、その設置したマイクロホンにより測定されたインパルス応答の組合せをさらに用いて上記と同様に合成インテンシティを算出することにより、各仮想音源の仰角を算出することができる。つまり、仮想音源情報推定部１３は、ｚ軸上に設置した２つのマイクロホンにより測定されたインパルス応答の組合せを用いて、0−Ｚ方向の瞬時インテンシティＩ_ｚ（ｔ）の包絡を取ったエンベロープインテンシティＩｈ_ｚ（ｔ）を算出する。仮想音源情報推定部１３は、エンベロープインテンシティＩｈ_ｘ（ｔ）、Ｉｈ_ｙ（ｔ）、Ｉｈ_ｚ（ｔ）のレベルをベクトル合成し、Ｌｅｖ＿ＰｓとＬｅｖ＿Ｐｅの範囲に表れる局所的なピークＰ_１、Ｐ_２、…、Ｐ_ｍを検出する。

次に、相互相関係数を用いた仮想音源情報の推定処理を説明する。
時刻ｔ_ａから時刻ｔ_ａ＋Δにおける２つの各インパルス応答ｈ_ｉ（ｔ）、ｈ_ｊ（ｔ）間の短時間相関係数Ｓ_ｉｊ（τ）は、次式（２）で与えられる。Δは分析窓長であり、τは時間差である。

異なる位置に設置された各マイクロホンに同一の反射音が到達する時刻にはずれがあるため、各マイクロホンにより測定されたインパルス応答には同一の反射音が異なる時刻に表れる。短時間相関係数Ｓ_ｉｊの値が大きいということは、同一の反射音が時間差τだけずれてｈ_ｉ（ｔ）とｈ_ｊ（ｔ）に出現することを表す。つまり、同一の反射音がτだけずれてマイクロホンＭ_ｉ、マイクロホンＭ_ｊに到達したことを表す。本実施形態では、ｘ軸上、ｙ軸上のそれぞれに設置した任意のマイクロホンにより測定されたインパルス応答の組合せを用いるため、ここでは、ｉ＝０，１，２、かつ、ｊ＝１，２（ｉ≠ｊ）とする。つまり、（ｉ，ｊ）の組合せには、（０，１）、（０，２）、（１，２）又は（２，１）の３つがある。以下では、ｉ＝０，１，２かつｊ＝１，２である場合を例に説明するが、ｘ軸上、ｙ軸上のそれぞれに設置した任意のマイクロホンにより測定されたインパルス応答の組合せ（例えば、ｉ＝０，３，４かつｊ＝３，４など）を用いることができる。

図８は、相互相関係数を用いた仮想音源情報推定処理のフロー図である。
まず、仮想音源情報推定部１３は、時刻ｔ_ａに初期値を設定する（ステップＳ１１０）。初期値は、例えば、インパルス応答ｈ_０（ｔ）において直接音が検出された時刻である。仮想音源情報推定部１３は、各（ｉ，ｊ）の組合せについて、相関係数Ｓ_ｉｊ（τ）に極大値を与える時間差τ＝τ_ｉｊを見つける（ステップＳ１２０）。仮想音源情報推定部１３は、各（ｉ，ｊ）の組合せごとに、相関係数Ｓ_ｉｊ（τ）がある一定の値（例えば、０．８）以上となる時間長（窓幅）Δ_ｉｊを求め、求めた中で最小の値の時間長Δ_{ｉｊ＿ｍｉｎ}を選択する（ステップＳ１３０）。時間長Δ_{ｉｊ＿ｍｉｎ}は、ある反射音が到達してから、次の反射音が到達するまでの最も短い時間に相当する。時間長Δ_ｉｊの最小値を選択することにより、窓幅の中に次の反射音が含まれないようにする。

次に、仮想音源情報推定部１３は、この時間長Δ_{ｉｊ＿ｍｉｎ}付近で補間処理を行う。ここでは、補間処理として、サンプリング周波数ｆｓを数十倍にあげる。本実施形態では、サンプリング周波数ｆｓを１６倍に上げる。仮想音源情報推定部１３は、補間処理により得られたインパルス応答ｈ_ｉ（ｔ）、ｈ_ｊ（ｔ）を用い、ステップＳ１２０と同様の処理を行い、相関係数Ｓ_ｉｊ（τ）に極大値を与え、かつ、より精度の高い時間差τ＝ｈτ_ｉｊを求める（ステップＳ１４０）。仮想音源情報推定部１３は、求めた時間差ｈτ_ｉｊを用いて仮想音源の座標及び角度を算出する（ステップＳ１５０）。

仮想音源とマイクロホンＭ_０の距離は、ｃを音速、仮想音源が放射された時刻をｔ_０とすると、（ｔ_ａ−ｔ_０）×ｃにより得られる。また、仮想音源とマイクロホンＭ_１の距離は、（ｔ_ａ＋ｈτ_０１−ｔ_０）×ｃにより、仮想音源とマイクロホンＭ_２の距離は、（ｔ_ａ＋ｈτ_０２−ｔ_０）×ｃにより得られる。マイクロホンＭ_０とｘ軸上のマイクロホンＭ_１との距離、及び、マイクロホンＭ_０とｙ軸上のマイクロホンＭ２との距離は既知（距離Ｌ）である。よって、これらのマイクロホン間の距離、及び、仮想音源と各マイクロホンＭ_０、Ｍ_１、Ｍ_２との距離とに基づいて、時刻ｔ_ａにおいて中心点Ｍ_０に到達した仮想音源の座標（Ｘ，Ｙ，Ｚ）が得られ、この座標から方位角及び仰角が得られる。

仮想音源情報推定部１３は、時刻ｔ_ａが分析終了の時刻に達していなければ（ステップＳ１６０：ＮＯ）、分析する窓のスタート地点である時刻ｔ_ａをΔ_{ｉｊ＿ｍｉｎ}だけずらして（ステップＳ１７０）、ステップＳ１２０からの処理を行う。仮想音源情報推定部１３は、時刻ｔ_ａが分析終了の時刻に達するまで（ステップＳ１６０：ＹＥＳ）、次々に仮想音源を検出していく。このような処理により、仮想音源情報推定部１３は、マイクロホン群の中心点Ｍ_０における仮想音源情報ｅ_１（ｒ_１,θ_１，φ_１，ｔ_１）、ｅ_２（ｒ_２,θ_２，φ_２，ｔ_２）、…、ｅ_ｍ（ｒ_ｍ,θ_ｍ，φ_ｍ，ｔ_ｍ）を推定することができる。

上記では、マイクロホンＭ_０〜Ｍ_４が同一平面（ｘｙ平面）上に設置されているため、仮想音源のＺ座標は０となり、従って、仰角φ_１、φ_２、…、φ_ｍも０となる。ｘｙ平面上にないｚ軸上のマイクロホンなど、同一平面上にないマイクロホンにより測定されたインパルス応答をさらに用いて上記と同様の処理を行うことで、各仮想音源のＺ座標が求められ、仰角も算出できる。

次に、図５のステップＳ３０における比較判定処理及び図５のステップＳ４０における反射音除去処理を説明する。
図９は、比較判定処理及び反射音除去処理の例を示す図である。ここでは、仮想音源情報推定の結果、４つの仮想音源情報ｅ_１（ｒ_１,θ_１，φ_１，ｔ_１）、ｅ_２（ｒ_２,θ_２，φ_２，ｔ_２）、ｅ_３（ｒ_３,θ_３，φ_３，ｔ_３）、ｅ_４（ｒ_４,θ_４，φ_４，ｔ_４）が推定された場合を例に説明する。比較判定部１４は、中心点Ｍ_０における仮想音源の情報ｅ_１〜ｅ_４から各仮想音源の到来時刻（ｔ_１，ｔ_２，ｔ_３，ｔ_４）を抽出する。比較判定部１４は、各仮想音源の到来時刻（ｔ_１，ｔ_２，ｔ_３，ｔ_４）と、測定されたインパルス応答ｈ_１（ｔ）、ｈ_２（ｔ）、ｈ_３（ｔ）、ｈ_４（ｔ）のそれぞれに含まれる局所的なピークの到来時刻（ｔ_１１，ｔ_１２，…），（ｔ_２１，ｔ_２２，…），（ｔ_３１，ｔ_３２，…），（ｔ_４１，ｔ_４２，…）とを比較し、取り除くべき反射音を判定する。取り除くべき反射音とは、あるインパルス応答において、そのインパルス応答が測定されたマイクロホンに対応する方向（マイクロホンを向けた方向）以外から到来する反射音である。

判定の流れを下記に示す。なお、ここで、ｍは予測した仮想音源の数（図９ではｍ＝４）であり、ｎはマイクロホンの数（図９ではｎ＝４）である。

（手順１）比較判定部１４は、仮想音源情報ｅ_ｋ（ｋは１以上ｍ以下の整数）ごとに、方位角θ_ｋに基づいて仮想音源の到来方向に対応するマイクロホンＭ_ｉ（ｉは１以上ｎ以下の整数）を特定し、特定したマイクロホンＭ_ｉにより測定されたインパルス応答ｈ_ｉ（ｔ）を選択する。具体的には、比較判定部１４は、設置方向が方位角θ_ｋであるマイクロホンにより測定されたインパルス応答を選択する。方位角θ_ｋが、いずれのマイクロホンＭ_１〜Ｍ_ｎの設置方向とも一致しない場合、比較判定部１４は、下記の（１）又は（２）のいずれかによりインパルス応答を選択する。

（１）比較判定部１４は、θ_ｋ方向の両隣に存在するマイクロホンのうち、設置方向がθ_ｋに近い方のマイクロホンのインパルス応答を選択する。この場合、図６に示すマイクロホンＭ_１、Ｍ_２、Ｍ_４、Ｍ_３で測定すべき反射音の方向は、例えば、方位角がｘ軸に対してそれぞれ０〜４５及び３１５〜３６０度、４５〜１３５度、１３５〜２２５度、２２５〜３１５度である。

（２）比較判定部１４は、θ_ｋ方向の両隣に存在するマイクロホンのそれぞれにより測定された２つのインパルス応答を選択する。この場合、マイクロホンＭ_１、Ｍ_２、Ｍ_４、Ｍ_３で測定すべき反射音の方向は、０〜９０及び２７０〜３６０度、０〜１８０度、９０〜２７０度、１８０〜３６０度である。

（手順２）比較判定部１４は、手順１により仮想音源情報ｅ_ｋについて選択したインパルス応答ｈ_ｉ（ｔ）において、ｔ_ｋ−τ_ｉ−ｗ／２≦ｔ≦ｔ_ｋ−τ_ｉ＋ｗ／２の範囲内で、局所的なピークが生ずる時刻ｔ＝ｔ_ｉｓを選択する。ここで、ｔ_ｋは、中心点Ｍ_０における仮想音源の到来時刻、τ_ｉは手順１において選択されたマイクロホンＭ_ｉから中心点Ｍ_０までの距離Ｌに応じた遅延、ｗは窓幅である。なお、窓幅ｗは予め設定されるが、隣接するマイクロホンとの距離が長い程大きな値とすることが望ましい。また、ｓは各インパルス応答における局所的なピークのインデックスであり、ｔ_ｉｓはインパルス応答ｈ_ｉ（ｔ）の波形構造においてｓ番目の局所的なピークが生ずる時刻である。

（手順３）比較判定部１４は、推定した仮想音源情報ごとに（手順１）及び（手順２）を繰り返し、各インパルス応答において選ばれなかった局所的なピークを、そのインパルス応答から除去すべき反射音であると判定する。反射音除去部１５は、各インパルス応答から、そのインパルス応答について比較判定部１４が除去すべきと判定した局所的なピークを除去する。

図９に示した例では、まず、比較判定部１４は、仮想音源情報ｅ_１（ｒ_１,θ_１，φ_１，ｔ_１）の到来方向を表すθ_１がマイクロホンＭ_１の設置方向であるため、マイクロホンＭ_１により測定されたインパルス応答ｈ_１（ｔ）を選択する。比較判定部１４は、インパルス応答ｈ_１（ｔ）において、ｔ_１−τ_１−ｗ／２≦ｔ≦ｔ_１−τ_１＋ｗ／２の範囲内で局所的なピークが生ずる時刻ｔ＝ｔ_１１を選択する。

次に、仮想音源情報ｅ_２（ｒ_２，θ_２，φ_２，ｔ_２）の到来方向を表すθ_２は、マイクロホンＭ_１の設置方向とマイクロホンＭ_２の設置方向の間の方位角であり、仮想音源はマイクロホンＭ_１とマイクロホンＭ_２の間から到来する。比較判定部１４は、マイクロホンＭ_１の設置方向よりもマイクロホンＭ_２の設置方向のほうがθ_２に近いため、マイクロホンＭ_２により測定されたインパルス応答ｈ_２（ｔ）を選択する（上記の手順１の（１）を適用）。比較判定部１４は、インパルス応答ｈ_２（ｔ）において、ｔ_２−τ_２−ｗ／２≦ｔ≦ｔ_２−τ_２＋ｗ／２の範囲内で局所的なピークが生ずる時刻ｔ＝ｔ_２２を選択する。

比較判定部１４は、仮想音源情報ｅ_３（ｒ_３，θ_３，φ_３，ｔ_３）及びｅ_４（ｒ_４，θ_４，φ_４，ｔ_４）についても同様に処理し、インパルス応答ｈ_３（ｔ）において時刻ｔ＝ｔ_３２、インパルス応答ｈ_４（ｔ）において時刻ｔ＝ｔ_４４を選択する。

反射音除去部１５は、インパルス応答ｈ_１（ｔ）の時刻ｔ_１１におけるピークＰ（ｔ_１１）、インパルス応答ｈ_２（ｔ）の時刻ｔ_２２におけるピークＰ（ｔ_２２）、インパルス応答ｈ_３（ｔ）の時刻ｔ_３２におけるピークＰ（ｔ_３２）、及び、インパルス応答ｈ_４（ｔ）の時刻ｔ_４４におけるピークＰ（ｔ_４４）以外の局所的なピークを除去する。これより、反射音除去部１５は、新たな方向別インパルス応答ｈ_１’（ｔ），ｈ_２’（ｔ），ｈ_３’（ｔ），ｈ_４’（ｔ）を生成する。

以上説明した実施形態によって、従来の技術よりも、測定音場に近似した方向別のインパルス応答を生成することが可能となる。
なお、平面上にない複数のマイクロホンを使用する場合、比較判定部１４は、方位角θ_ｋだけでなく、仰角φ_ｋをさらに加えて表される仮想音源の到来方向に基づいて、インパルス応答を選択してもよい。つまり、比較判定部１４は、方位角θ_ｋ及び仰角φ_ｋにより表される到来方向に主軸を向けたマイクロホン、あるいは、その到来方向の両隣に存在するマイクロホン又はそれらマイクロホンのうち主軸を向けた方向が到来方向に近いマイクロホンを特定し、特定したマイクロホンにより測定されたインパルス応答を選択する。
また、上記で示した例では中心点Ｍ_０と、それを取り囲むようにマイクロホン群を等間隔に設置しているが、必ずしもその必要はなく、任意のマイクロホン配置から、マイクロホン群で囲まれた空間内の任意の点における仮想音源情報を推定することができる。

上述した実施形態によれば、インパルス応答生成装置１は、空間内の異なる位置に設置された複数のマイクロホンが同一の音響信号を収音して得られた信号それぞれに基づき算出された各インパルス応答を用いて、複数のマイクロホンに囲まれた任意の位置における、空間内の仮想音源の情報を推定する。インパルス応答生成装置１は、仮想音源の情報とインパルス応答との比較に基づいて、任意のインパルス応答において、そのインパルス応答の算出の元となった信号が得られたマイクロホンに対応する方向から到来する反射音を判定する。インパルス応答生成装置１は、インパルス応答から、そのインパルス応答の算出の元となった信号が得られたマイクロホンに対応する方向以外から到来する余計な反射音を間引いて方向別のインパルス応答を生成する。これにより、インパルス応答生成装置１は、従来の手法よりも、測定音場に近い方向別のインパルス応答を生成することができる。

なお、上述のインパルス応答生成装置１は、内部にコンピュータシステムを有している。そして、インパルス応答生成装置１の動作の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータシステムが読み出して実行することによって、上記処理が行われる。ここでいうコンピュータシステムとは、ＣＰＵ及び各種メモリやＯＳ、周辺機器等のハードウェアを含むものである。

また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。

１ インパルス応答生成装置
１１ 記憶部
１２ 測定情報取得部
１３ 仮想音源情報推定部
１４ 比較判定部
１５ 反射音除去部

Claims (4)

1. 同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出されたインパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源それぞれの到来時刻及び到来方向を推定する仮想音源情報推定部と、
   前記仮想音源情報推定部により推定された前記仮想音源ごとに、複数の前記マイクロホンのうち前記任意の位置からの方向が前記仮想音源の前記到来方向に対応したマイクロホンを選択し、前記仮想音源の前記到来時刻と、当該仮想音源について選択された前記マイクロホンにより収音した前記信号から算出された前記インパルス応答における局所的なピークそれぞれが出現する時刻との差分に基づいて、当該インパルス応答における前記局所的なピークのうち、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音のピークを特定する特定処理を行い、前記特定処理により特定された前記インパルス応答における前記ピークを、当該インパルス応答の算出の元となった前記信号が得られた前記マイクロホンに対応する方向から到来する反射音と判定する比較判定部と、
   前記インパルス応答から、当該インパルス応答における前記ピークのうち前記比較判定部により反射音と判定された前記ピーク以外のピークを除去して方向別のインパルス応答を生成する反射音除去部と、
   を備えることを特徴とするインパルス応答生成装置。
2. 前記比較判定部は、前記仮想音源ごとの前記特定処理において、複数の前記マイクロホンのうち前記任意の位置からの方向が前記仮想音源の前記到来方向に対応したマイクロホンを選択し、選択した前記マイクロホンにより収音した信号から算出された前記インパルス応答における局所的なピークのうち、選択した前記マイクロホンと前記任意の位置との距離に応じた遅延を前記仮想音源の前記到来時刻から減算した時刻を含んだ所定の時刻範囲内で生じた前記ピークを、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音のピークと特定する、
   ことを特徴とする請求項１に記載のインパルス応答生成装置。
3. 前記仮想音源情報推定部は、前記任意の位置からの方向が異なる複数の前記マイクロホンそれぞれにより収音された前記信号に基づき算出された前記インパルス応答と、前記任意の位置に設置された前記マイクロホンにより前記音響信号を収音して得られた信号に基づき算出されたインパルス応答とを用いて算出した瞬時インテンシティに基づいて前記仮想音源の前記到来時刻及び前記到来方向を得る、又は、前記任意の位置からの方向が異なる複数の前記マイクロホンそれぞれにより収音された前記信号に基づき算出された前記インパルス応答と、前記任意の位置に設置された前記マイクロホンにより前記音響信号を収音して得られた信号に基づき算出されたインパルス応答との相関に基づいて前記仮想音源の前記到来時刻及び前記到来方向を得る、
   ことを特徴とする請求項１又は請求項２に記載のインパルス応答生成装置。
4. コンピュータを、
   同一の音響信号を空間内の異なる位置に設置された複数のマイクロホンにより収音して得られた信号それぞれに基づき算出されたインパルス応答を用いて、複数の前記マイクロホンに囲まれた任意の位置における、前記空間内の仮想音源それぞれの到来時刻及び到来方向を推定する仮想音源情報推定手段と、
   前記仮想音源情報推定手段により推定された前記仮想音源ごとに、複数の前記マイクロホンのうち前記任意の位置からの方向が前記仮想音源の前記到来方向に対応したマイクロホンを選択し、前記仮想音源の前記到来時刻と、当該仮想音源について選択された前記マイクロホンにより収音した前記信号から算出された前記インパルス応答における局所的なピークそれぞれが出現する時刻との差分に基づいて、当該インパルス応答における前記局所的なピークのうち、当該インパルス応答の算出の元となった信号が得られた前記マイクロホンに対応する方向から到来する反射音のピークを特定する特定処理を行い、前記特定処理により特定された前記インパルス応答における前記ピークを、当該インパルス応答の算出の元となった前記信号が得られた前記マイクロホンに対応する方向から到来する反射音と判定する比較判定手段と、
   前記インパルス応答から、前記比較判定手段により当該インパルス応答において反射音と判定された前記ピーク以外の前記ピークを除去して方向別のインパルス応答を生成する反射音除去手段と、
   を具備するインパルス応答生成装置として機能させるためのプログラム。

Images