(57)【特許請求の範囲】
【請求項１】
所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて指定された領域から伝播される音の音圧波形を再生する再生装置であって、
前記計測した音源データを記憶するデータ記憶手段と、
前記データ記憶手段に記憶されている音源の水平角のデータと仰角のデータと音圧信号の大きさのデータとを用いて水平角と仰角とを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップを作成するマップ作成手段と、
前記音源マップを複数の空間領域に分割するマップ分割手段と、
前記各空間領域の周波数毎の音圧信号の大きさからオクターブバンド毎の音圧信号の大きさを前記空間領域毎に算出するとともに、前記オクターブバンド毎の音圧信号の大きさを前記空間領域毎に合成して前記空間領域におけるバンドパワー値をそれぞれ算出するバンドパワー値算出手段と、
前記算出されたオクターブバンド毎のバンドパワー値の比であるバンドパワー比を求めるバンドパワー比算出手段と、
オクターブバンドパスフィルターを備え、前記データ記憶手段に記憶された音源データを取り出して、オクターブバンド毎の音圧波形であるバンド波形をそれぞれ抽出するバンド波形抽出手段と、
前記抽出されたバンド波形の大きさと前記算出されたバンドパワー比とを用いて前記算出されたバンドパワー値をオクターブバンド毎に補正するバンドパワー値補正手段と、
前記補正されたバンドパワー値を用いて、前記複数の空間領域の内の指定された空間領域である指定領域から伝播される音の音圧波形を再生する音圧波形再生手段とを備えたことを特徴とする指定領域からの伝播音の再生装置。
【請求項２】
前記音圧波形再生手段で再生された音圧波形を入力して、指定領域から伝播される音を出力する伝播音再生手段を設けたことを特徴とする請求項１に記載の指定領域からの伝播音の再生装置。
【請求項３】
前記データ記憶手段に前記観測点から撮影した音源方向の映像の画像データを記憶するとともに、
前記音源マップと同様の複数の空間領域に分割された表示画面のうちの少なくとも前記指定領域を含む空間領域に前記画像データを表示する表示手段を設けたことを特徴とする請求項１または請求項２に記載の指定領域からの伝播音の再生装置。
【請求項４】
前記音源データを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対を有するマイクロフォン群と、前記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンと、前記２組のマイクロフォン対を構成するマイクロフォン間のそれぞれの位相差を求め、この求められた２組のマイクロフォン対の位相差の比から音源方向の水平角を推定するとともに、前記２組のマイクロフォン対を構成するマイクロフォン間の位相差と、前記第５のマイクロフォンと前記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差とを用いて前記音源方向の仰角を推定する音源方向推定手段とを備えた音源推定装置を用いて測定された音源データとしたことを特徴とする請求項１〜請求項３のいずれかに記載の指定領域からの伝播音の再生装置。
【請求項５】
所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて指定された領域から伝播される音を再生する方法であって、
前記音源の音源方向のデータである水平角のデータ及び仰角のデータと音圧信号の大きさのデータとを用いて水平角と仰角とを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップを作成するステップと、
前記音源マップを複数の空間領域に分割するステップと、
前記各空間領域の周波数毎の音圧信号の大きさからオクターブバンド毎の音圧信号の大きさを前記空間領域毎に算出するとともに、前記オクターブバンド毎の音圧信号の大きさを前記空間領域毎に合成して前記空間領域におけるバンドパワー値をそれぞれ算出するステップと、
前記算出されたオクターブバンド毎のバンドパワー値の比であるバンドパワー比を求めるステップと、
前記算出されたバンドパワー比を用いて前記算出されたバンドパワー値をオクターブバンド毎に基準化するステップと、
前記基準化されたバンドパワー値を用いて、前記複数の空間領域の内の指定された空間領域である指定領域から伝播される音の音圧波形を再生するステップとを備えたことを特徴とする指定領域からの伝播音の再生方法。
【請求項６】
所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて指定された領域から伝播される音を再生する方法であって、
前記音源の音源方向のデータである水平角のデータ及び仰角のデータと音圧信号の大きさのデータとを用いて水平角と仰角とを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップを作成するステップと、
前記音源マップを複数の空間領域に分割するステップと、
前記各空間領域の周波数毎の音圧信号の大きさからオクターブバンド毎の音圧信号の大きさを前記空間領域毎に算出するとともに、前記オクターブバンド毎の音圧信号の大きさを前記空間領域毎に合成して前記空間領域におけるバンドパワー値をそれぞれ算出するステップと、
前記算出されたオクターブバンド毎のバンドパワー値の比であるバンドパワー比を求めるステップと、
前記各空間領域の周波数毎の音圧信号から、オクターブバンド毎の音圧波形であるバンド波形をそれぞれ抽出するステップと、
前記抽出されたバンド波形の大きさと前記算出されたバンドパワー比とを用いて前記算出されたバンドパワー値をオクターブバンド毎に補正するステップと、
前記補正されたバンドパワー値を用いて、前記複数の空間領域の内の指定された空間領域である指定領域から伝播される音の音圧波形を再生するステップとを備えたことを特徴とする指定領域からの伝播音の再生方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数のマイクロフォンで採取した音の情報を用いて推定した音源方向のデータと音圧信号のデータとを用いて、指定された領域から伝播される音を再生する方法とその装置に関するものである。
【背景技術】
【０００２】
従来、音の到来方向を推定する方法としては、多数のマイクロフォンを等間隔に配置したマイクロフォンアレーを構築し、基準となるマイクロフォンに対する各マイクロフォンの位相差から音波の到来方向である音源方向を推定する、いわゆる音響学的手法が考案されている（例えば、非特許文献１参照）。
一方、複数のマイクロフォンから互いに交わる直線状に配置された複数のマイクロフォン対を構成し、対となる２つのマイクロフォン間の位相差に相当する到達時間差と、他の対となる２つのマイクロフォン間の到達時間差との比から音源方向を推定する方法が提案されている（例えば、特許文献１〜３参照）。
【０００３】
具体的には、図１２に示すように、４個のマイクロフォンＭ１〜Ｍ４を、互いに直交する２直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対（Ｍ１，Ｍ３）及びマイクロフォン対（Ｍ２，Ｍ４）を構成するように配置し、一方のマイクロフォン対（Ｍ１，Ｍ３）を構成するマイクロフォンＭ１，Ｍ３に入力する音圧信号の到達時間差Ｄ13と、他方のマイクロフォン対（Ｍ２，Ｍ４）を構成するマイクロフォンＭ２，Ｍ４に入力する音圧信号の到達時間差Ｄ24との比から、計測点と音源の位置との水平角θを推定する。また、マイクロフォンＭ１〜Ｍ４の作る平面上にない位置に第５のマイクロフォンＭ５を配置して４組のマイクロフォン対（Ｍ５, Ｍ１），（Ｍ５, Ｍ２），（Ｍ５, Ｍ３），（Ｍ５, Ｍ４）を構成して、前記２組のマイクロフォン対を構成するマイクロフォン間の到達時間差Ｄ13，Ｄ24と前記４組のマイクロフォン対を構成するマイクロフォン間の到達時間差Ｄ5j（ｊ＝１〜４）から、計測点と音源の位置との成す仰角φを推定する。
計測点から測った音源方向は前記推定された水平角θと仰角φとにより表わせる。音源方向は周波数毎に求められる。
【０００４】
これにより、マイクロフォンアレーを用いて音源方向を推定する場合に比較して、少ないマイクロフォン数で音源方向を正確に推定することができる。
また、ＣＣＤカメラ等の映像採取手段を設けて推定された音源方向の映像を撮影し、この映像の画像データと音源方向のデータとを合成して、映像中に推定した音源方向と音圧レベルとを図形で表示した音源推定用画像をディスプレイ等の表示画面に表示すれば、表示画面から音源を特定することも可能である。
図１３は、表示画面３７Ｄに音源方向推定画像を表示した一例を示す図で、横軸は音源の水平角θ、縦軸は音源の仰角φである。網目模様の丸印及び斜線を施した丸印の中心がそれぞれ音源Ｇ１〜Ｇ３の中心位置で、枠で囲った音源Ｇ１，Ｇ２は直接音の音源、Ｇ３は反射音の音源である。これにより、音源を視覚的に把握することができる。また、丸の大きさは音圧信号の大きさを表す。なお、図１３において、周波数は丸印Ｃの模様で区別しているが、表示画面３７Ｄをカラー画面とし周波数を色彩で区別してもよい。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００２−１８１９１３号公報
【特許文献２】特開２００６−３２４８９５号公報
【特許文献３】特開２００８−２２４２５９号公報
【非特許文献】
【０００６】
【非特許文献１】大賀寿郎，山崎芳男，金田豊；音響システムとディジタル処理，コロナ社，１９９５
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、前記従来の方法では、周波数毎の音源を視覚的に把握したり、特定の音源からの音を周波数毎に再生することはできるが、指定された領域から伝播される複数の周波数を含む音を容易にかつ短時間で再生することが困難であった。
【０００８】
本発明は、従来の問題点に鑑みてなされたもので、所定の観測点で計測した音源データを用いて指定された領域からの伝播音を短時間でかつ容易に再生することのできる方法とその装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本願の請求項１に記載の発明は、所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて指定された領域から伝播される音の音圧波形を再生する再生装置であって、前記計測した音源データを記憶するデータ記憶手段と、前記データ記憶手段に記憶されている音源の水平角のデータと仰角のデータと音圧信号の大きさのデータとを用いて水平角と仰角とを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップを作成するマップ作成手段と、前記音源マップを複数の空間領域に分割するマップ分割手段と、前記各空間領域の周波数毎の音圧信号の大きさからオクターブバンド毎の音圧信号の大きさを前記空間領域毎に算出するとともに、前記オクターブバンド毎の音圧信号の大きさを前記空間領域毎に合成して前記空間領域におけるバンドパワー値をそれぞれ算出するバンドパワー値算出手段と、前記算出されたオクターブバンド毎のバンドパワー値の比であるバンドパワー比を求めるバンドパワー比算出手段と、オクターブバンドパスフィルターを備え、前記データ記憶手段に記憶された音源データを取り出して、オクターブバンド毎の音圧波形であるバンド波形をそれぞれ抽出するバンド波形抽出手段と、前記抽出されたバンド波形の大きさと前記算出されたバンドパワー比とを用いて前記算出されたバンドパワー値をオクターブバンド毎に補正するバンドパワー値補正手段と、前記補正されたバンドパワー値を用いて、前記複数の空間領域の内の指定された空間領域である指定領域から伝播される音の音圧波形を再生する音圧波形再生手段とを備えたことを特徴とするものである。
これにより、容易にかつ短時間で指定領域から伝播される音の音圧波形を再生することができる。また、指定領域を複数個とすれば、大きな音もしくは特徴的な周波数を有する複数の空間領域からの音の音圧波形を再生することができる。
【００１０】
請求項２に記載の発明は、請求項１に記載の指定領域からの伝播音の再生装置において、前記音圧波形再生手段で再生された音圧波形を入力して、指定領域から伝播される音を出力する伝播音再生手段を設けたものである。
これにより、伝播音の音圧波形のみならず伝播音そのものも再生できるので、指定領域から伝播される音の特徴を確実に把握することができる。
請求項３に記載の発明は、請求項１または請求項２に記載の指定領域からの伝播音の再生装置において、前記データ記憶手段に前記観測点から撮影した音源方向の映像の画像データを記憶するとともに、前記音源マップと同様の複数の空間領域に分割された表示画面のうちの少なくとも前記指定領域を含む空間領域に前記画像データを表示する表示手段を設けたものである。
これにより、指定領域の音源についても視覚的に把握することができるので、どの音源からどのような音が伝播されるかを把握することができる。
【００１１】
請求項４に記載の発明は、請求項１〜請求項３のいずれかに記載の指定領域からの伝播音の再生装置であって、前記音源データを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対を有するマイクロフォン群と、前記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンと、前記２組のマイクロフォン対を構成するマイクロフォン間のそれぞれの位相差を求め、この求められた２組のマイクロフォン対の位相差の比から音源方向の水平角を推定するとともに、前記２組のマイクロフォン対を構成するマイクロフォン間の位相差と、前記第５のマイクロフォンと前記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差とを用いて前記音源方向の仰角を推定する音源方向推定手段とを備えた音源推定装置を用いて測定された音源データとしたことを特徴とする。
このような構成の音源推定装置を用いることにより、マイクロフォンアレーを用いた場合に比較して、少ないマイクロフォン数で音源方向と音圧信号の大きさとを正確に推定することができる。
【００１２】
請求項５に記載の発明は、所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて指定された領域から伝播される音を再生する方法であって、前記音源の音源方向のデータである水平角のデータ及び仰角のデータと音圧信号の大きさのデータとを用いて水平角と仰角とを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップを作成するステップと、前記音源マップを複数の空間領域に分割するステップと、前記各空間領域の周波数毎の音圧信号の大きさからオクターブバンド毎の音圧信号の大きさを前記空間領域毎に算出するとともに、前記オクターブバンド毎の音圧信号の大きさを前記空間領域毎に合成して前記空間領域におけるバンドパワー値をそれぞれ算出するステップと、前記算出されたオクターブバンド毎のバンドパワー値の比であるバンドパワー比を求めるステップと、前記算出されたバンドパワー比を用いて前記算出されたバンドパワー値をオクターブバンド毎に基準化するステップと、前記基準化されたバンドパワー値を用いて、前記複数の空間領域の内の指定された空間領域である指定領域から伝播される音の音圧波形を再生するステップとを備えたことを特徴とする。
これにより、複雑な計算を行うことなく、短時間で指定領域から伝播される音の音圧波形を再生することができる。
【００１３】
請求項６に記載の発明は、所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて指定された領域から伝播される音を再生する方法であって、前記音源の音源方向のデータである水平角のデータ及び仰角のデータと音圧信号の大きさのデータとを用いて水平角と仰角とを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップを作成するステップと、前記音源マップを複数の空間領域に分割するステップと、前記各空間領域の周波数毎の音圧信号の大きさからオクターブバンド毎の音圧信号の大きさを前記空間領域毎に算出するとともに、前記オクターブバンド毎の音圧信号の大きさを前記空間領域毎に合成して前記空間領域におけるバンドパワー値をそれぞれ算出するステップと、前記算出されたオクターブバンド毎のバンドパワー値の比であるバンドパワー比を求めるステップと、前記各空間領域の周波数毎の音圧信号から、オクターブバンド毎の音圧波形であるバンド波形をそれぞれ抽出するステップと、前記抽出されたバンド波形の大きさと前記算出されたバンドパワー比とを用いて前記算出されたバンドパワー値をオクターブバンド毎に補正するステップと、前記補正されたバンドパワー値を用いて、前記複数の空間領域の内の指定された空間領域である指定領域から伝播される音の音圧波形を再生するステップとを備えたことを特徴とする。
このように、音源の音圧信号を取り出してオクターブバンド毎の音圧波形であるバンド波形をそれぞれ抽出し、この抽出されたバンド波形の大きさと算出されたバンドパワー比とを用いて算出されたバンドパワー値をオクターブバンド毎に補正するようにすれば、指定領域から伝播される音の大きさについても精度良く再生することができる。
【００１４】
なお、前記発明の概要は、本発明の必要な全ての特徴を列挙したものではなく、これらの特徴群のサブコンビネーションもまた、発明となり得る。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施の形態に係る伝播音の再生装置の構成を示す図である。
【図２】音源マップの一例を示す図である。
【図３】音源マップの分割方法の一例を示す図である。
【図４】バンドパワー値の算出方法を示す図である。
【図５】バンドパワー値の算出方法を示す図である。
【図６】バンドパワー値の補正方法を示す図である。
【図７】指定領域から伝播される音の音圧波形を再生する方法を示す図である。
【図８】表示用マップの一例を示す図である。
【図９】本実施の形態に係る伝播音の再生方法を示すフローチャートである。
【図１０】音源推定装置の構成を示す図である。
【図１１】本発明による伝播音の再生方法の他の例を示すフローチャートである。
【図１２】従来のマイクロフォン対を用いた音源探査方法におけるマイクロフォンの配列を示す図である。
【図１３】音源方向推定画像を表示した表示画面の一例を示す図である。
【発明を実施するための形態】
【００１６】
以下、実施の形態を通じて本発明を詳説するが、以下の実施の形態は特許請求の範囲に係る発明を限定するものでなく、また、実施の形態の中で説明される特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
【００１７】
図１は本発明の実施の形態に係る伝播音の再生装置１０の構成を示す図で、伝播音の再生装置１０は、データ記憶手段１１と、マップ作成手段１２と、マップ分割手段１３と、バンドパワー値算出手段１４と、バンドパワー比算出手段１５と、バンド波形抽出手段１６と、バンドパワー値補正手段１７と、音圧波形再生手段１８と、伝播音再生手段１９と、画像処理手段２０と音源表示手段２１とを備える。
データ記憶手段１１はパーソナルコンピュータのハードディスクなどから構成され、マップ作成手段１２〜音圧波形再生手段１８の各手段と画像処理手段２０はパーソナルコンピュータのソフトウェアなどにより構成される。また、伝播音再生手段１９はスピーカーなどにより構成され、音源表示手段２１はパーソナルコンピュータのディスプレイなどから構成される。
データ記憶手段１１は、後述する音源推定装置３０を用いて所定の観測点で計測した１つもしくは複数の音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データと前記観測点から撮影した音源方向の映像の画像データとを記憶する。音源方向のデータは観測点と推定される音源との成す水平角θのデータと仰角φのデータである、なお、音圧信号の大きさは周波数毎に求められる。
【００１８】
マップ作成手段１２は、横軸を水平角θとし縦軸を仰角φとしたマップを作成するとともに、このマップ上に、データ記憶手段１１に記憶されている各音源の音源方向のデータ（θ，φ）と、音源方向が（θ，φ）で周波数がｆである音圧信号の大きさａf（θ，φ）とを表示して、図２に示すような音源マップ１２Ｍを作成する。なお、実際には、コンピュータのメモリー上に、音源方向のデータ（θ，φ）と音圧信号の大きさａf（θ，φ）のデータとを水平角θと仰角φとを座標軸としたマップに割り付けるだけで表示はしないが、表示した場合には、図２に示すように、従来と同じく中心が音源方向のデータ（θ，φ）で直径が音圧信号の大きさａf（θ，φ）で、模様が周波数を示す丸印Ｃで表示する。なお、図２では、図を見やすくするため、表示した周波数の種類を５種類とした。
マップ分割手段１３は、図３に示すように、音源マップ１２Ｍを複数の空間領域に分割する。
バンドパワー値算出手段１４は、図４に示すように、周波数毎の音圧信号の大きさａm,n,f（θ，φ）を合成してオクターブバンド毎の音圧信号の大きさａm,n,F（θ，φ）を空間領域Ｇm,n毎に求めた後、図５に示すように、オクターブバンド毎の音圧信号の大きさａm,n,F（θ，φ）を空間領域Ｇm,n毎に合成して空間領域Ｇm,nにおける音圧信号の大きさであるバンドパワー値ｐm,n（Ｆk）をオクターブバンド毎に算出する。同図において、丸印Ｄm,n,kは各空間領域Ｇm,n毎のバンドパワー値ｐm,n（Ｆk）を表示したもので、丸印Ｄm,n,kの大きさはバンドパワー値ｐm,n（Ｆk）の大きさを示し、模様はオクターブバンドの中心周波数Ｆkを示す。
Ｆkはオクターブバンドの中心周波数を示す。本例では、オクターブバンドとして１オクターブバンドを使用した。中心周波数Ｆ1〜Ｆ9は、それぞれ、３１．５Ｈｚ，６３Ｈｚ，１２５Ｈｚ，２５０Ｈｚ，５００Ｈｚ，１ｋＨｚ，２ｋＨｚ，４ｋＨｚ，８ｋＨｚである。
なお、帯域幅は、Ｆk−Ｆk/√２〜Ｆk＋Ｆk/√２である（ｋ＝１〜９）。
【００１９】
バンドパワー比算出手段１５は、バンドパワー値算出手段１４で算出されたオクターブバンド毎のバンドパワー値ｐm,n（Ｆk）の比であるパワー比Ｒm,n（Ｆk）をオクターブバンド毎に求める。
バンド波形抽出手段１６は、オクターブバンドパスフィルターを備え、データ記憶手段１１に記憶された観測された音の音圧信号の大きさＡを取り出してオクターブバンドパスフィルターを通すことで、オクターブバンド毎の音圧信号を抽出し、オクターブバンド毎のバンドパスの音の大きさＡ（Ｆk）を求める（ｋ＝１〜９）。
バンドパワー値補正手段１７は、図６に示すように、各空間領域Ｇm,nのオクターブバンド毎のバンドパスの音の大きさＡ（Ｆk）をパワー比Ｒm,n（Ｆk）を用いて各空間領域Ｇm,nに割り付けることで、バンドパワー値ｐm,n（Ｆk）を実際に観測された音の音圧信号の大きさＡで補正した補正バンドパワー値Ｐm,n（Ｆk）に変換する。
音圧波形再生手段１８は、図７に示すように、バンドパワー値補正手段１７で補正されたバンドパワー値Ｐm,n（Ｆk）を用いて、複数の空間領域Ｇm,nの内の指定された空間領域である指定領域ＧM,Nから伝播される音の音圧波形を再生する。
【００２０】
伝播音再生手段１９は、例えば、スピーカーなどから構成され、音圧波形再生手段１８で再生された音圧波形を入力して、指定領域から伝播される音を出力する。
画像処理手段２０は、水平角θと仰角φとを座標軸とした表示用マップを作成するとともにこの表示用マップを空間領域に分割し、各空間領域にデータ記憶手段１１に記憶されている観測点から撮影した音源方向の映像の画像データと補正されたバンドパワー値Ｐm,n（Ｆk）とを割り付ける。表示用マップの水平角θの範囲と仰角φの範囲及び分割方法は音源マップ１２Ｍと同一であるので、音源マップ１２Ｍをそのまま利用してもよい。
音源表示手段２１は表示画面２１Ｍを有し、表示画面２１Ｍ上に、図８（ａ），（ｂ）に示すような画像処理手段２０で作成された表示用マップを表示する。（ａ）図は主音源が１箇所の場合の例で、（ｂ）図は主音源が２箇所の場合の例である。
【００２１】
次に、本発明による伝播音の再生方法について、図９のフローチャートを参照して説明する。
まず、データ記憶手段１１から所定の観測点で計測した音源の音源方向と音圧信号の大きさとを含む音源データを取り出し、（ステップＳ１０）、図２に示すような、横軸が水平角θで縦軸を仰角φであるマップ上に、音源方向（θ，φ）と音圧信号の大きさａf（θ，φ）とを表示した音源マップ１２Ｍを作成する（ステップＳ１１）。
本例では、所定の観測点で計測した音源の音源方向と音圧信号の大きさとを含む音源データを、図１０に示すような音源推定装置３０を用いて求めた。具体的には、図２に示すように、図示しない工場の入り口を観測点として、工場内の工作機械Ｈのどの個所から騒音が発生するかを計測した。
音源推定装置３０は、図示しない騒音源からの騒音の音圧レベルを測定するために観測点に配置された計測用マイクロフォンＭ１〜Ｍ５と、音源位置近傍の映像を採取するためのＣＣＤカメラ（以下、カメラという）３１と、ローパスフィルタを備え、計測用マイクロフォンＭ１〜Ｍ５で採取された音響情報から所定の周波数以下の成分を取り出し増幅する増幅器３２と、増幅された音響情報（アナログ信号）をデジタル信号に変換するＡ／Ｄ変換器３３と、カメラ３１の映像情報（アナログ信号）をデジタル信号に変換するビデオ入出力ユニット３４と、Ａ／Ｄ変換された計測用マイクロフォンＭ１〜Ｍ５の音圧信号を用いて音源方向とこの音源からの音の大きさとを推定する音源方向推定手段３５と、上記映像信号に上記推定された音源方向を示す画像を付加した画像を生成する画像合成手段３６と、合成された画像を表示する音源位置表示手段３７とを備えている。
【００２２】
また、４０は計測用マイクロフォンＭ１〜Ｍ５を所定の位置に配列するためのマイクロフォンフレーム、４１は三脚から成る支持部材４２とこの支持部材４２の上部に配設された回転台４３とから成る測定用基台である。回転台４３により、マイクロフォンフレーム４０を回転させて、計測用マイクロフォンＭ１〜Ｍ５を水平面内で回転させることにより、音源方向を全方位にわたって推定することができる。
計測用マイクロフォンＭ１〜Ｍ４を、図１２に示すように、互いに直交する２直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対（Ｍ１，
Ｍ３）及びマイクロフォン対（Ｍ２，Ｍ４）を構成するように配置するとともに、マイクロフォンＭ５を、マイクロフォンＭ１〜Ｍ４の作る平面上にない位置に配置する。これにより、上記各マイクロフォン対（Ｍｉ, Ｍｊ）の位相差（時間遅れＤij）から、当該観測点から見た音源方向を推定することができる。
音の入射方向である水平角θと仰角φとは以下の式（１）及び式（２）で表わせる。
【数１】
0005693201.tif 000002
ここで、時間遅れＤijは、マイクロフォンＭiに到達する音圧信号と、このマイクロフォンＭiに対して対となるマイクロフォンＭjに到達する音圧信号との時間差であり、この対となる２つのマイクロフォンＭi及びＭjに入力される信号のクロススペクトルＰij（ｆ）を求め、更に、対象とする上記周波数ｆの位相角情報Ψ（ｒａｄ）を用いて、以下の式（３）により算出される。
【数２】
0005693201.tif 000003
音源方向推定手段３５では、Ａ／Ｄ変換器３３でＡ／Ｄ変換された計測用マイクロフォンＭ１〜Ｍ５で採取した音圧信号のデータを用いて、音源方向である水平角θと仰角φとを推定するとともに音圧信号の大きさを計測する。
なお、音源方向と音圧信号の大きさとは周波数毎に計測する。また、マイクロフォンＭ５に入力される信号の大きさＡ5を、観測される音の音圧信号の大きさＡとする。
一方、画像合成手段３６では、上記ビデオ入出力ユニット３４に入力された音源方向の映像信号に、推定された音源方向を示す画像を付加した画像を生成し、これを音源位置表示手段３７に送って表示する。そして、上記音源方向が示された画像から、音源の位置を決定する。
なお、本例では、音源方向のデータと音圧信号の大きさのデータと、カメラ３１で撮影されビデオ入出力ユニット３４でＡ／Ｄ変換された画像データとをデータ記憶手段１１に記憶させるようにしているが、画像合成手段３６で合成された画像のデータをデータ記憶手段１１に記憶しておいてもよい。この場合には、音源マップ１２Ｍの作成を省略することができる。
【００２３】
次に、音源マップ１２Ｍを複数の空間領域に分割（ステップＳ１２）した後、空間領域Ｇm,n毎にオクターブバンド毎の音圧信号の大きさであるバンドパワー値ｐm,n（Ｆk）を算出する（ステップＳ１３）。
図３は音源マップ１２Ｍの分割例を示す図で、水平角θの範囲を及び仰角φの範囲をそれぞれ−８０°＜θ≦＋８０°、−５０°＜φ≦＋５０°、水平角θの分割幅を４０°、仰角φの分割幅を２０°とすれば、音源マップ１２Ｍは２０個の空間領域Ｇm,n（ｍ＝１〜４、ｎ＝１〜５）に分割される。
バンドパワー値の算出は、まず図４に示すように、周波数毎の音圧信号の大きさａm,n,f（θ，φ）を合成してオクターブバンド毎の音圧信号の大きさａm,n,F（θ，φ）を求め、次に、図５に示すように、空間領域Ｇm,n毎にオクターブバンド毎の音圧信号の大きさａm,n,F（θ，φ）を合成して空間領域Ｇm,nにおけるバンドパワー値ｐm,n（Ｆk）をそれぞれ算出する。すなわち、空間領域Ｇm,nに複数のオクターブバンド毎の音圧信号があった場合には、それら複数の音圧信号を１つの音圧信号にまとめることでバンドパワー値ｐm,n（Ｆk）を算出する。このとき、算出されたバンドパワー値ｐm,n（Ｆk）が予め設定された閾値Ｋ（Ｆk）未満である場合には、空間領域Ｇm,nには中心周波数がＦkであるオクターブバンドの音圧信号がないとみなす。これにより、計算時間を短くすることができる。
【００２４】
バンドパワー値ｐm,n（Ｆk）の算出後には、バンドパワー値ｐm,n（Ｆk）の比であるパワー比Ｒm,n（Ｆk）をオクターブバンド毎に求める（ステップＳ１４）。
例えば、図５に示すように、予め設定された閾値Ｋ（Ｆ5）以上の大きさの５００Ｈｚバンドのバンドパワー値ｐm,n（Ｆ5）を有する音源が存在する空間領域はＧ2,2，Ｇ3,2，Ｇ3,3，Ｇ4,3の４つで、その他の空間領域Ｇm,nではｐm,n（Ｆ5）＝０である場合、空間領域Ｇm,n（Ｆ5）におけるパワー比Ｒm,n（Ｆ5）は、下記の式から算出される。
Ｒm,n（Ｆk）＝ｐm,n（Ｆ5）/｛ｐ1,1（Ｆ5）＋ｐ2,2（Ｆ5）＋……＋ｐ20,20（Ｆ5）｝
例えば、空間領域Ｇ2,2，Ｇ3,2，Ｇ3,3，Ｇ4,3のバンドパワー値がそれぞれｐ2,2（Ｆ5）＝１０，ｐ3,2（Ｆ5）＝６，ｐ3,3（Ｆ5）＝３，ｐ4,3（Ｆ5）＝１、その他の空間領域Ｇm,nではｐm,n（Ｆ5）＝０とすれば、空間領域Ｇ2,2，Ｇ3,2，Ｇ3,3，Ｇ4,3のパワー比Ｒm,n（Ｆ5）は、それぞれ、Ｒ2,2（Ｆ5）＝０．５，Ｒ3,2（Ｆ5）＝０．３，Ｒ3,3（Ｆ5）＝０．１５，Ｒ4,3（Ｆ5）＝０．０５となり、その他の空間領域Ｇm,nではｐm,n（Ｆ5）＝０となる。
【００２５】
ステップＳ１５では、データ記憶手段１１に記憶された音源の音圧信号の大きさＡからオクターブバンド毎の音圧波形であるバンド波形を抽出し、前記バンド波形の大きさであるオクターブバンド毎のバンドパスの音の大きさＡ（Ｆk）を求める。ここで、Ｂ（Ｆk，ｔ）をオクターブバンド毎の時刻歴バンドパス波形データとすると、音源の音圧信号ｆ（ｔ）は、ｆ（ｔ）＝Ａ（Ｆ1）・Ｂ（Ｆ1，ｔ）＋Ａ（Ｆ2）・Ｂ（Ｆ2，ｔ）＋……＋Ａ（Ｆ8）・Ｂ（Ｆ8，ｔ）＋Ａ（Ｆ9）・Ｂ（Ｆ9，ｔ）と近似することができる。
なお、このステップは、次のステップＳ１６の前であればどこで行ってもよい。
ステップＳ１６では、ステップＳ１４で算出されたバンドパワー値ｐm,n（Ｆk）を観測点で計測された音の大きさを用いて実際に観測されたバンドパワー値に補正する。この補正されたパワー値を以下補正バンドパワー値Ｐm,n（Ｆk）という。すなわち、オクターブバンド毎のバンドパスの音の大きさＡ（Ｆk）は、各空間領域Ｇm,nのオクターブバンド毎の補正バンドパワー値Ｐm,n（Ｆk）の二乗和の平均に等しいことから、図６に示すように、オクターブバンド毎のバンドパスの音の大きさＡ（Ｆk）をパワー比Ｒm,n（Ｆk）を用いて各空間領域Ｇm,nのバンドパワー値ｐm,n（Ｆk）に割り付ければ、バンドパワー値ｐm,n（Ｆk）を実際の音圧波形のバンドパワー値である補正バンドパワー値Ｐm,n（Ｆk）に補正することができる。
これら補正バンドパワー値Ｐm,n（Ｆk）を用いることで、各空間領域Ｇm,nから伝播される音の音圧波形を音の大きさを含めて再生することができる。
本例では、図７に示すように、複数の空間領域Ｇm,nのうちから音を再生する領域を指定し（ステップＳ１７）、この指定された空間領域である指定領域ＧM,Nから伝播される音の音圧波形のみを再生する（ステップＳ１８）。
具体的には、指定領域ＧM,Nのバンドパワー値ＰM,N（Ｆk）の平方根をａM,N（Ｆk）とし、オクターブバンド毎の時刻歴バンドパス波形データをＢ（Ｆk，ｔ）とすると、指定領域ＧM,Nから伝播される音の音圧波形ｆM,N（ｔ）は以下の式で表わせる。
ｆM,N（ｔ）＝ａM,N（Ｆ1）・Ｂ（Ｆ1，ｔ）＋ａM,N（Ｆ2）・Ｂ（Ｆ2，ｔ）＋……
……＋ａM,N（Ｆ8）・Ｂ（Ｆ8，ｔ）＋ａM,N（Ｆ9）・Ｂ（Ｆ9，ｔ）
【００２６】
ステップＳ１９では、水平角θと仰角φとを座標軸とし複数の空間領域Ｇm,nに分割された表示用マップ上に、観測点から撮影した音源方向の映像の画像データを割り付けた表示画面２１Ｍを表示するとともに、スピーカーを用いてステップＳ１８で求められた指定領域ＧM,Nから伝播される音の音圧波形ｆM,N（ｔ）を再生する。このとき、図８（ａ），（ｂ）に示すように、表示画面２１Ｍ上において指定領域ＧM,Nを枠で囲むなどして、スピーカーから出力される音と指定領域ＧM,Nに移っている建物あるいは機器などの発する音の音源と伝播音との対応が取れるようにすれば、音源とその音源の発する特徴的な音とを確実に把握することができる。
次に、音の再生を終了するか否か決定する（ステップＳ２０）。
別の指定領域ＧM’,N’から伝播される音を再生する場合にはステップＳ１７に戻り、指定領域ＧM’,N’ のバンドパワー値ＰM’,N’（Ｆk）を用いて指定領域ＧM’,N’から伝播される音の音圧波形ｆM’,N’（ｔ）を再生する。
【００２７】
このように、本実施の形態では、音源推定装置３０を用いて求めた所定の観測点で計測した音源の音源方向のデータと音圧信号の大きさのデータとを含む音源データを用いて、観測点から見た水平角θと仰角φとを座標軸としたマップ上に音源方向と音圧信号の大きさとを周波数毎に表示した音源マップ１２Ｍを作成した後、この音源マップ１２Ｍを複数の空間領域Ｇm,nに分割し、各空間領域Ｇm,nの周波数毎の音圧信号の大きさからオクターブバンドＦk毎の音圧信号の大きさであるバンドパワー値ｐm,n（Ｆk）をそれぞれ算出するとともに、バンドパワー値の比であるパワー比Ｒm,n（Ｆk）を求め、このパワー比Ｒm,n（Ｆk）とオクターブバンドパスフィルターに音圧信号を通して得られたオクターブバンド毎のバンドパスの音の大きさＡ（Ｆk）とを用いてバンドパワー値ｐm,n（Ｆk）を補正し、補正されたバンドパワー値Ｐm,n（Ｆk）を用いて、指定された空間領域である指定領域ＧM,Nから伝播される音の音圧波形を再生するようにしたので、容易にかつ短時間で指定領域ＧM,Nから伝播される音の音圧波形を再生することができる。
また、スピーカーなどの伝播音再生手段１９を設けて指定領域ＧM,Nから伝播される音を出力するとともに、指定領域ＧM,Nを含む空間領域に音源方向の画像データをディスプレイなどの音源表示手段２１に表示したので、指定領域ＧM,Nから伝播される音の特徴を確実に把握することができる。
【００２８】
なお、前記実施の形態では、音源を工場内の工作機械Ｈが発生する騒音としたが、観測対象はこれに限るものではなく、図１３に示した工場のある区域や橋梁などのような屋外の騒音発生源であってもよい。
また、前記例では、音源マップ１２Ｍを２０個の空間領域Ｇm,nに分割したが、分割数はこれに限るものではなく、観測対象により任意に設定すればよい。
また、音源データの取得については、マイクロフォンアレーを用いた他の音源推定装置を用いてもよいが、本例のように、音源推定装置３０を用いれば、少ないマイクロフォン数で音源方向と音圧信号の大きさとを正確に推定することができる。
また、上記例では、周波数毎の音圧信号の大きさａm,n,f（θ，φ）を合成してオクターブバンド毎の音圧信号の大きさａm,n,F（θ，φ）を求めてから指定領域ＧM,Nから伝播される音の音圧波形を再生するようにしたが、オクターブバンドを更に分割した１/２オクターブバンド毎に音圧信号の大きさａm,n,F（θ，φ）を求めたり、１/３オクターブバンド毎に音圧信号の大きさａm,n,F（θ，φ）を求めたりすれば、再生精度を更に向上させることができる。
【００２９】
なお、各空間領域Ｇm,nから伝播される音の音圧波形の音の大きさについては再現する必要がない場合、すなわち、音質のみが必要な場合には、バンド波形抽出手段１６とバンドパワー値補正手段１７とを省略し、バンドパワー比算出手段１５で算出したパワー比を用いて音声波形を再生すればよい。この場合のフローチャートを図１１に示す。ステップＳ１０〜ステップＳ１４までのステップは図９と同じである。
ステップＳ１４の直後のステップＳ２１では、パワー比Ｒm,n（Ｆk）を用いて算出されたバンドパワー値ｐm,n（Ｆk）をオクターブバンド毎に基準化する。ステップＳ１４で述べたように、パワー比Ｒm,n（Ｆk）の総計は１となるので、基準化した音の大きさをＡ0とすると、基準化されたバンドパワー値ｐ’m,n（Ｆk）はｐ’m,n（Ｆk）＝Ａ0・Ｒm,n（Ｆk）となる。
ステップＳ２２では、ステップＳ１７と同様に、複数の空間領域Ｇm,nのうちから音を再生する領域を指定する。
ステップＳ２３では、この指定された空間領域である指定領域ＧM,Nから伝播される音の音圧波形のみを再生する。指定領域ＧM,Nの基準化されたバンドパワー値ｐ’M,N（Ｆk）の平方根をａ’M,N（Ｆk）とすると、指定領域ＧM,Nから伝播される音の音圧波形ｆ’M,N（ｔ）は、ｆ’M,N（ｔ）＝ａ’M,N（Ｆ1）・Ｂ（Ｆ1，ｔ）＋ａ’M,N（Ｆ2）・Ｂ（Ｆ2，ｔ）＋……ａ’M,N（Ｆ8）・Ｂ（Ｆ8，ｔ）＋ａ’M,N（Ｆ9）・Ｂ（Ｆ9，ｔ）となる。ここで、Ｂ（Ｆk，ｔ）はオクターブバンド毎の時刻歴バンドパス波形データである。
【００３０】
ステップＳ２４では、水平角θと仰角φとを座標軸とし複数の空間領域Ｇm,nに分割された表示用マップ上に、観測点から撮影した音源方向の映像の画像データを割り付けた表示画面２１Ｍを表示するとともに、スピーカーを用いてステップＳ２３で求められた指定領域ＧM,Nから伝播される音の音圧波形ｆ’M,N（ｔ）を再生する。
その後、音の再生を終了するか否か決定し（ステップＳ２５）、別の指定領域ＧM’,N’から伝播される音を再生する場合にはステップＳ２２に戻り、パワー比Ｒm,n（Ｆk）を用いて算出された基準化されたバンドパワー値ｐ’m,n（Ｆk）を用いて指定領域ＧM’,N’から伝播される音の音圧波形ｆ’M’,N’（ｔ）を再生する。
【産業上の利用可能性】
【００３１】
以上説明したように、本発明によれば、指定領域から伝播される音の音圧波形を容易にかつ短時間で再生することができるので、観測された音源の特徴を効率よく把握することができる。
【符号の説明】
【００３２】
１０  伝播音の再生装置、１１  データ記憶手段、１２  マップ作成手段、
１２Ｍ  音源マップ、１３  マップ分割手段、１４  バンドパワー値算出手段、
１５  バンドパワー比算出手段、１６  バンド波形抽出手段、
１７  バンドパワー値補正手段、１８  音圧波形再生手段、１９  伝播音再生手段、
２０  画像処理手段、２１  音源表示手段、２１Ｍ  表示画面、
３０  音源推定装置、Ｍ１〜Ｍ５  計測用マイクロフォン、３１  ＣＣＤカメラ、
３２  増幅器、３３  Ａ／Ｄ変換器、３４  ビデオ入出力ユニット、
３５  音源方向推定手段、３６  画像合成手段、３７  音源位置表示手段、
４０  マイクロフォンフレーム、４１  測定用基台、４２  支持部材、４３  回転台。
【図１】
0005693201.tif 000004
【図２】
0005693201.tif 000005
【図３】
0005693201.tif 000006
【図４】
0005693201.tif 000007
【図５】
0005693201.tif 000008
【図６】
0005693201.tif 000009
【図７】
0005693201.tif 000010
【図８】
0005693201.tif 000011
【図９】
0005693201.tif 000012
【図１０】
0005693201.tif 000013
【図１１】
0005693201.tif 000014
【図１２】
0005693201.tif 000015
【図１３】
0005693201.tif 000016