JPWO2020003343A1 - 波源方向推定装置、波源方向推定方法、およびプログラム - Google Patents
波源方向推定装置、波源方向推定方法、およびプログラム Download PDFInfo
- Publication number
- JPWO2020003343A1 JPWO2020003343A1 JP2020526725A JP2020526725A JPWO2020003343A1 JP WO2020003343 A1 JPWO2020003343 A1 JP WO2020003343A1 JP 2020526725 A JP2020526725 A JP 2020526725A JP 2020526725 A JP2020526725 A JP 2020526725A JP WO2020003343 A1 JPWO2020003343 A1 JP WO2020003343A1
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- wave source
- estimation
- direction information
- source direction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 336
- 238000001228 spectrum Methods 0.000 claims description 206
- 230000010354 integration Effects 0.000 claims description 85
- 238000006243 chemical reaction Methods 0.000 claims description 72
- 238000005314 correlation function Methods 0.000 claims description 29
- 230000010365 information processing Effects 0.000 claims description 14
- 230000006870 function Effects 0.000 description 66
- 238000010586 diagram Methods 0.000 description 22
- 238000004891 communication Methods 0.000 description 7
- 239000006185 dispersion Substances 0.000 description 7
- 230000007423 decrease Effects 0.000 description 3
- 230000005236 sound signal Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000143306 Idaea sardoniata Species 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000001404 mediated effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
- G01S3/808—Systems for determining direction or deviation from predetermined direction using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/801—Details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
- G01S3/808—Systems for determining direction or deviation from predetermined direction using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
- G01S3/8083—Systems for determining direction or deviation from predetermined direction using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems determining direction of source
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/06—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being correlation coefficients
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/18—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R1/00—Details of transducers, loudspeakers or microphones
- H04R1/20—Arrangements for obtaining desired frequency or directional characteristics
- H04R1/32—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only
- H04R1/40—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers
- H04R1/406—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2430/00—Signal processing covered by H04R, not provided for in its groups
- H04R2430/20—Processing of the output signals of the acoustic transducers of an array for obtaining a desired directivity characteristic
- H04R2430/21—Direction finding using differential microphone array [DMA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Otolaryngology (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computational Linguistics (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
Description
まず、本発明の第1の実施形態に係る波源方向推定装置について図面を参照しながら説明する。以下においては、本実施形態の波源方向推定装置が、空気や水の振動波である音波の発生源(波源とも呼ぶ)の方向を推定する例について説明する。そのため、本実施形態の波源方向推定装置は、マイクロフォン(以下、マイクとも呼ぶ)によって電気信号に変換された音波の波源の方向を推定対象とする。なお、本実施形態の波源方向推定装置の推定対象は、音波の波源の方向に限定されず、振動波や電磁波などの任意の波動の波源の方向を推定することに使用できる。
図1は、本実施形態の波源方向推定装置10の構成を示すブロック図である。波源方向推定装置10は、複数の入力端子11、信号選択部12、相対遅延時間計算部13、複数の推定方向情報・周波数寄与率生成部15、重み計算部16、および統合部17を備える。
複数の入力端子11−1〜p(入力部とも呼ぶ)のそれぞれは、図示しないマイクに接続される。複数の入力端子11−1〜pのそれぞれには、異なる位置に配置された複数のマイクによって集音された音波(音信号とも呼ぶ)から変換された電気信号が入力信号として入力される。以下において、時刻tにおいてm番目の入力端子11−mに入力した入力信号をxm(t)と記載する(t:実数、m:1以上p以下の整数)。
信号選択部12は、入力端子11−1〜pに供給されるp個の入力信号の中から2つの入力信号を選択する。信号選択部12は、選択した2つの入力信号を推定方向情報・周波数寄与率生成部15−1〜Rのいずれかに出力し、それらの入力信号の供給元であるマイクの位置情報(以下、マイク位置情報とも呼ぶ)を相対遅延時間計算部13に出力する。
相対遅延時間計算部13には、信号選択部12からマイク位置情報が入力される。相対遅延時間計算部13は、予め設定された音源探索対象方向(音源方向とも呼ぶ)と、マイク位置情報とを用いて、信号選択部12が選択した全てのマイクペアごとに相対遅延時間を算出する。相対遅延時間とは、2つのマイクのマイク間隔と、音源方向とに基づいて一意に定まる音波の到達時間差のことである。例えば、音源方向は、所定の角度刻みで設定される。すなわち、相対遅延時間計算部13は、設定された音源方向について相対遅延時間を計算する。相対遅延時間計算部13は、算出した音源方向と相対遅延時間とのセットを推定方向情報・周波数寄与率生成部15に出力する。
推定方向情報・周波数寄与率生成部15−1〜Rのそれぞれには、信号選択部12によって選択された全マイクペアのうち1組のマイクペアの入力信号が入力される。また、推定方向情報・周波数寄与率生成部15−1〜Rのそれぞれには、入力信号の供給元のマイクペアの相対遅延時間が相対遅延時間計算部13から入力される。推定方向情報・周波数寄与率生成部15−1〜Rのそれぞれは、1組のマイクペアの入力信号と相対遅延時間とを用いて、そのマイクペアの入力信号間の推定方向情報と周波数寄与率とを生成する。推定方向情報・周波数寄与率生成部15−1〜Rは、推定方向情報を統合部17に出力し、周波数寄与率を重み計算部16に出力する。
変換部151には、信号選択部12から1組のマイクペアの入力信号(入力信号Aおよび入力信号B)が入力される。変換部151は、信号選択部12から供給された2つの入力信号を周波数領域の変換信号に変換する。変換部151は、入力信号を複数の周波数成分に分解するための変換を行う。例えば、変換部151は、フーリエ変換を用いて、入力信号を複数の周波数成分に分解する。変換部151は、クロススペクトル計算部152に変換信号を出力する。
クロススペクトル計算部152には、変換部151から変換信号が入力される。クロススペクトル計算部152は、変換部151から供給される変換信号を用いてクロススペクトルを計算する。クロススペクトル計算部152は、算出したクロススペクトルを平均計算部153に出力する。
平均計算部153には、クロススペクトル計算部152からクロススペクトルが入力される。平均計算部153は、クロススペクトル計算部152から供給されたクロススペクトルの平均(平均クロススペクトルとも呼ぶ)を計算する。平均計算部153は、算出した平均クロススペクトルを周波数別クロススペクトル生成部155に出力する。
分散計算部154には、平均計算部153から平均クロススペクトルが入力される。分散計算部154は、平均計算部153から供給された平均クロススペクトルを用いて分散を計算する。分散計算部154は、算出した分散を周波数別クロススペクトル生成部155に出力する。
周波数別クロススペクトル生成部155には、平均計算部153から平均クロススペクトルが入力され、分散計算部154から分散が入力される。周波数別クロススペクトル生成部155は、平均クロススペクトルと分散とを用いて平均クロススペクトルの周波数に対応するクロススペクトル(周波数別クロススペクトルとも呼ぶ)を計算する。周波数別クロススペクトル生成部155は、算出した周波数別クロススペクトルを逆変換部156に出力する。
周波数別基本クロススペクトル計算部551には、平均計算部153から平均クロススペクトルが入力される。周波数別基本クロススペクトル計算部551は、平均計算部153から供給される平均クロススペクトルを用いて、平均クロススペクトルの各周波数に対応するクロススペクトル(周波数別基本クロススペクトルとも呼ぶ)を計算する。周波数別基本クロススペクトル計算部551は、算出した周波数別基本クロススペクトルを乗算部553に出力する。周波数別基本クロススペクトルは、周波数成分ごとに相関関数を計算するために算出される。周波数別基本クロススペクトル計算部551は、ある周波数に対応する相関関数(周波数別相互相関関数とも呼ぶ)を後段で求めるための周波数別基本クロススペクトルを計算する。
カーネル関数スペクトル生成部552には、分散計算部154から分散が入力される。カーネル関数スペクトル生成部552は、分散計算部154から供給される分散を用いてカーネル関数スペクトルを算出する。カーネル関数スペクトルとは、カーネル関数をフーリエ変換し、その絶対値を取ったものである。なお、カーネル関数スペクトルには、カーネル関数をフーリエ変換し、その二乗を取ったものとしてもよい。また、カーネル関数スペクトルには、カーネル関数をフーリエ変換し、その絶対値の二乗を取ったものを用いてもよい。カーネル関数スペクトル生成部552は、算出したカーネル関数スペクトルを乗算部553に出力する。
Lは平均計算部153が平均クロススペクトルを求める際に平均化したフレーム数を表す。平均クロススペクトルの誤差は平均化フレーム数Lに反比例するため、式26および式27を用いることで、平均クロススペクトルの誤差(信頼性)を考慮して広がり制御パラメータを求めることができる。
乗算部553には、周波数別基本クロススペクトル計算部551から周波数別基本クロススペクトルが入力され、カーネル関数スペクトル生成部552からカーネル関数スペクトルが入力される。乗算部553は、周波数別基本クロススペクトル計算部551から供給される周波数別基本クロススペクトルと、カーネル関数スペクトル生成部552から供給されるカーネル関数スペクトルとの積を計算し、周波数別クロススペクトルを算出する。乗算部553は、算出した周波数別クロススペクトルを逆変換部156に出力する。
逆変換部156には、周波数別クロススペクトル生成部155の乗算部553から周波数別クロススペクトルが入力される。例えば、変換部151がフーリエ変換を用いて変換を行う場合、逆変換部156は逆フーリエ変換を用いて逆変換を行う。逆変換部156は、周波数別クロススペクトル生成部155から供給された周波数別クロススペクトルの逆変換を求める。
周波数別推定方向情報計算部157には、逆変換部156から周波数別相互相関関数が入力され、相対遅延時間計算部13から相対遅延時間が入力される。周波数別推定方向情報計算部157は、逆変換部156から供給された周波数別相互相関関数と、相対遅延時間計算部13から供給された相対遅延時間とを用いて、方向と相関値との対応関係を周波数別推定方向情報として求める。周波数別推定方向情報計算部157は、求めた周波数別推定方向情報を統合部17に出力する。
周波数寄与率計算部158には、周波数別推定方向情報計算部157から周波数別推定方向情報が入力される。周波数寄与率計算部158は、周波数別推定方向情報計算部157から供給された周波数別推定方向情報に基づいて、音源方向ごとの推定方向情報に対する周波数別の寄与率(周波数寄与率とも呼ぶ)を計算する。周波数寄与率計算部158は、算出した周波数寄与率を重み計算部16に出力する。
周波数統合部159には、周波数別推定方向情報計算部157から周波数別推定方向情報が入力される。周波数統合部159は、周波数別推定方向情報計算部157から供給された周波数別推定方向情報を統合する。すなわち、周波数統合部159は、周波数(周波数ビン番号)ごとに別々に算出された周波数別推定方向情報を一つに統合することによって、統合された推定方向情報を生成する。例えば、周波数統合部159は、個別に求めた複数の周波数別推定方向情報を合成したり、重ね合わせたりすることにより、一つの推定方向情報に統合する。周波数統合部159は、統合された推定方向情報を統合部17に出力する。
重み計算部16には、推定方向情報・周波数寄与率生成部15−1〜Rのそれぞれの周波数寄与率計算部158から周波数寄与率が入力される。重み計算部16は、推定方向情報・周波数寄与率生成部15−1〜Rから供給される周波数寄与率に基づいて、マイクペア間の周波数寄与率(単に寄与率とも呼ぶ)の類似度に応じた音源方向ごとの重みを計算する。重み計算部16は、算出した音源方向ごとの重みを統合部17に出力する。
統合部17には、推定方向情報・周波数寄与率生成部15−1〜Rから推定方向情報が入力され、重み計算部16から重みが入力される。統合部17は、推定方向情報・周波数寄与率生成部15−1〜Rから供給されるマイクペアごとの推定方向情報を全てのマイクペアについて統合し、統合した推定方向情報に重み計算部16から供給される重みをつけた統合推定方向情報を生成する。統合部17は、マイクペアごとの推定方向情報の総和や総乗を取ることで統合し、統合した推定方向情報に重みを掛けることによって統合推定方向情報を生成する。統合部17は、生成した統合推定方向情報を出力する。例えば、統合部17は、上位システム(図示しない)に統合推定方向情報を出力する。
次に、本実施形態の波源方向推定装置10の動作の概略について図面を参照しながら説明する。図9は、波源方向推定装置10の動作の概略について説明するためのフローチャートである。なお、図9のフローチャートに沿った説明においては、波源方向推定装置10を動作の主体として説明する。
次に、本発明の第2の実施形態に係る波源方向推定装置について図面を参照しながら説明する。
複数の入力端子21−1〜p(入力部とも呼ぶ)のそれぞれは、図示しないマイクに接続される(pは2以上の整数)。複数の入力端子21−1〜pのそれぞれには、異なる位置に配置されたマイクによって集音された音波(音信号とも呼ぶ)から変換された電気信号が入力信号として入力される。なお、入力端子21−1〜pは、第1の実施形態の入力端子21−1〜pと同様であるため、詳細な説明は省略する。
信号選択部22は、入力端子21−1〜pに供給されるp個の入力信号のうちの2つの入力信号を選択する。信号選択部22は、選択した2つの入力信号を周波数別推定方向情報生成部25−1〜Rに出力し、それらの入力信号の供給元であるマイクの位置情報(以下、マイク位置情報とも呼ぶ)を相対遅延時間計算部23に出力する(Rは1以上の整数)。なお、信号選択部22は、第1の実施形態の信号選択部12と同様であるため、詳細な説明は省略する。
相対遅延時間計算部23には、信号選択部22からマイク位置情報が入力される。相対遅延時間計算部23は、予め設定された音源探索対象方向(音源方向とも呼ぶ)と、マイク位置情報とを用いて、信号選択部22が選択した全てのマイクペアごとに相対遅延時間を計算する。相対遅延時間計算部23は、設定された音源方向と相対遅延時間とをセットにして周波数別推定方向情報生成部25に出力する。なお、相対遅延時間計算部23は、第1の実施形態の相対遅延時間計算部13と同様であるため、詳細な説明は省略する。
周波数別推定方向情報生成部25−1〜Rのそれぞれには、信号選択部22によって選択された全マイクペアのうち1組のマイクペアの入力信号が入力される。また、周波数別推定方向情報生成部25−1〜Rのそれぞれには、入力信号の供給元のマイクペアの相対遅延時間が相対遅延時間計算部23から入力される。周波数別推定方向情報生成部25−1〜Rのそれぞれは、1組のマイクペアの入力信号と相対遅延時間とを用いて、そのマイクペアごとに周波数別推定方向情報を生成する。周波数別推定方向情報生成部25−1〜Rは、生成したマイクペアごとの周波数別推定方向情報を重み生成部26および統合部27に出力する。
重み生成部26には、周波数別推定方向情報生成部25−1〜Rのそれぞれからマイクペアごとの周波数別推定方向情報が入力される。重み生成部26は、周波数別推定方向情報生成部25−1〜Rのそれぞれから供給されるマイクペアごとの周波数別推定方向情報を用いて、推定方向情報がピークをもつ音源方向ごとの重みを計算する。重み生成部26は、算出した音源方向ごとの重みを統合部27に出力する。
周波数統合部261−1〜Rのそれぞれには、周波数別推定方向情報生成部25−1〜Rのそれぞれから周波数別推定方向情報が入力される。周波数統合部261−1〜Rのそれぞれは、周波数別推定方向情報生成部25−1〜Rのそれぞれから供給される周波数別推定方向情報を統合し、推定方向情報を計算する。なお、周波数統合部261の動作は、第1の実施形態の周波数統合部159と同様であるため、詳細な説明は省略する。周波数統合部261−1〜Rのそれぞれは、算出した推定方向情報を波源方向候補選出部262−1〜Rのそれぞれに出力する。
波源方向候補選出部262−1〜Rのそれぞれには、周波数統合部261−1〜Rのそれぞれから推定方向情報が入力される。波源方向候補選出部262−1〜Rのそれぞれは、周波数統合部261−1〜Rのそれぞれから供給される推定方向情報からピークを検出し、ピークの方向を波源方向候補として選出する。波源方向候補選出部262−1〜Rのそれぞれは、選出した波源方向候補を波源方向候補統合部263に出力する。
波源方向候補統合部263には、波源方向候補選出部262−1〜Rのそれぞれから波源方向候補が入力される。波源方向候補統合部263は、波源方向候補選出部262−1〜Rのそれぞれによって選出された波源方向候補を統合する。波源方向候補統合部263は、統合した波源方向候補を波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれに出力する。
波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれには、周波数別推定方向情報生成部25−1〜Rのそれぞれからマイクペアごとの周波数別推定方向情報が入力される。また、波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれには、統合された波源方向候補が波源方向候補統合部263から入力される。波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれは、マイクペアごとの周波数別推定方向情報と、統合された波源方向候補とを用いて、波源方向候補ごとの推定方向情報に対する周波数寄与率である波源方向候補別周波数寄与率を計算する。波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれは、算出した波源方向候補別周波数寄与率を重み計算部265に出力する。第1の実施形態の周波数寄与率計算部158は、全ての音源方向について周波数寄与率を求める。それに対し、本実施形態の波源方向候補別周波数寄与率計算部264は、選出された波源方向候補について周波数寄与率を求める。
重み計算部265には、波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれから波源方向候補別周波数寄与率が入力される。重み計算部265は、波源方向候補別周波数寄与率計算部264−1〜Rのそれぞれから供給される波源方向候補別周波数寄与率を用いて、マイクペア間の周波数寄与率の類似度に応じた波源方向候補ごとの重みを計算する。重み計算部265は、算出した波源方向候補ごとの重みを統合部27に出力する。
統合部27には、周波数別推定方向情報生成部25−1〜Rのそれぞれから周波数別推定方向情報が入力され、重み生成部26から波源方向候補ごとの重みが入力される。統合部27は、周波数別推定方向情報生成部25−1〜Rから供給される周波数別推定方向情報を波源方向候補ごとに全てのマイクペアについて統合し、統合した波源方向候補ごとの周波数別推定方向情報に重みを付け、波源方向候補ごとの統合推定方向情報を算出する。統合部27は、マイクペアごとの周波数別推定方向情報の総和や相乗を取ることで統合し、統合した推定方向情報に重みを掛けることによって波源方向候補ごとの統合推定方向情報を生成する。例えば、統合部27は、以下の式45を用いて、波源方向候補Φごとの統合推定方向情報H(Φ,n)を算出する。
次に、本実施形態の波源方向推定装置20の動作の概略について図面を参照しながら説明する。図13は、波源方向推定装置20の動作の概略について説明するためのフローチャートである。なお、図13のフローチャートに沿った説明においては、波源方向推定装置20を動作の主体として説明する。
次に、本発明の第3の実施形態に係る波源方向推定装置について図面を参照しながら説明する。本実施形態の波源方向推定装置は、第1および第2の実施形態の波源推定装置を上位概念化した構成である。なお、本実施形態においては、波源を音源に限定せず、任意の波動の波源の方向を推定対象とする。
複数の入力端子31−1〜p(入力部とも呼ぶ)のそれぞれは、図示しないマイクに接続される(pは2以上の整数)。複数の入力端子31−1〜pのそれぞれには、異なる位置に配置されたセンサによって集音された波動(波動信号とも呼ぶ)から変換された電気信号が入力信号として入力される。なお、入力端子31−1〜pは、第1および第2の実施形態の対応する構成と同様であるため、詳細な説明は省略する。
信号選択部32は、入力端子31−1〜pに供給されるp個の入力信号のうちの2つの入力信号(入力信号のペアとも呼ぶ)を複数組選択する。信号選択部32は、選択した複数組の入力信号のペアを統合推定方向情報算出部35に出力し、それらの入力信号の供給元であるセンサの位置情報(以下、センサ位置情報とも呼ぶ)を相対遅延時間計算部33に出力する。なお、信号選択部32は、第1および第2の実施形態の対応する構成と同様であるため、詳細な説明は省略する。
相対遅延時間計算部33には、信号選択部32からセンサ位置情報が入力される。相対遅延時間計算部33は、センサ位置情報と音源探索対象方向(音源方向とも呼ぶ)とを用いて、信号選択部32が選択した全ての入力信号のペア間の相対遅延時間を計算する。相対遅延時間計算部33は、設定された音源方向と相対遅延時間とをセットにして統合推定方向情報算出部35に出力する。なお、相対遅延時間計算部33は、第1および第2の実施形態の対応する構成と同様であるため、詳細な説明は省略する。
統合推定方向情報算出部35には、複数組の入力信号のペアが信号選択部32から入力される。また、統合推定方向情報算出部35には、信号選択部32によって選択された全ての入力信号のペアの相対遅延時間が相対遅延時間計算部33から入力される。
次に、本発明の第4の実施形態に係る波源方向推定装置について図面を参照しながら説明する。本実施形態の波源方向推定装置は、第3の実施形態の波源方向推定装置に波源方向算出部を追加した構成を有する。
波源方向算出部46には、統合推定方向情報算出部45から統合推定方向情報が入力される。波源方向算出部46は、統合推定方向情報を用いて、波源方向を計算する。波源方向算出部46は、計算した波源方向を出力する。
ここで、各実施形態に係る波源方向推定装置の処理を実行するハードウェア構成について、図16の情報処理装置90を一例として挙げて説明する。なお、図16の情報処理装置90は、各実施形態の波源方向推定装置の処理を実行するための構成例であって、本発明の範囲を限定するものではない。
(付記1)
複数のセンサによって検出される波動に基づく電気信号を入力信号として取得する複数の入力手段と、
複数の前記入力信号のうち二つを組み合わせたペアを複数組選択する信号選択手段と、
前記ペアをなす二つの前記入力信号の供給元の前記センサへの波源方向ごとの波動の到達時間差を相対遅延時間として算出する相対遅延時間計算手段と、
前記ペアをなす前記入力信号と前記ペアの相対遅延時間とを用いて周波数別推定方向情報を前記ペアごとに生成し、生成した前記ペアごとの前記周波数別推定方向情報を用いて波源方向ごとの推定方向情報に対する周波数別の寄与率を算出し、算出した前記寄与率に基づいて異なる前記ペア間の前記寄与率の類似度に応じた波源方向ごとの重みを計算し、全ての前記ペアの前記推定方向情報に前記重みをつけて統合することによって統合推定方向情報を生成する統合推定方向情報算出手段とを備える波源方向推定装置。
(付記2)
前記信号選択手段は、
前記センサの間隔に基づいて、複数の前記入力信号から二つの前記入力信号を組み合わせた前記ペアを選択する付記1に記載の波源方向推定装置。
(付記3)
前記相対遅延時間計算手段は、
一つの前記ペアをなす二つの前記入力信号の供給元である前記センサにおける波源探索方向を基準として、前記信号選択手段によって選択された全ての前記ペアの前記相対遅延時間を、基準である前記波源探索方向の関数として計算する付記1または2に記載の波源方向推定装置。
(付記4)
前記統合推定方向情報算出手段は、
前記信号選択手段によって選択された全ての前記ペアのうち1組の前記ペアをなす二つの前記入力信号を取得するとともに、取得した前記ペアの前記相対遅延時間を取得し、前記ペアをなす二つの前記入力信号間の周波数別相互相関関数と前記ペアの前記相対遅延時間とを用いて、前記ペアをなす二つの前記入力信号間の前記推定方向情報と周波数別の前記寄与率とを生成する複数の推定方向情報・周波数寄与率生成手段と、
複数の前記推定方向情報・周波数寄与率生成手段から周波数ごとの前記寄与率を取得し、異なる前記ペア間の前記寄与率の類似度に応じて波源方向ごとの前記重みを計算する重み計算手段と、
複数の前記推定方向情報・周波数寄与率生成手段から前記推定方向情報を取得するとともに、前記重み計算手段から前記重みを取得し、複数の前記ペアごとの前記推定方向情報に前記重みをつけて統合して前記統合推定方向情報を生成する統合手段とを有する付記1乃至3のいずれか一項に記載の波源方向推定装置。
(付記5)
前記推定方向情報・周波数寄与率生成手段は、
前記ペアをなす二つの前記入力信号を周波数領域の変換信号に変換する変換手段と、
前記変換手段によって変換された前記変換信号を用いてクロススペクトルを計算するクロススペクトル計算手段と、
前記クロススペクトル計算手段によって算出される前記クロススペクトルを用いて平均クロススペクトルを計算する平均計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルを用いて分散を計算する分散計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルと、前記分散計算手段によって算出される分散とを用いて周波数別クロススペクトルを計算する周波数別クロススペクトル生成手段と、
前記周波数別クロススペクトル生成手段によって算出される前記周波数別クロススペクトルを逆変換して周波数別相互相関関数を計算する逆変換手段と、
前記ペアの前記周波数別相互相関関数と前記相対遅延時間とに基づいて波源方向と相関値との対応関係である周波数別推定方向情報を計算する周波数別推定方向情報計算手段と、
前記周波数別推定方向情報を統合して前記推定方向情報を計算する周波数統合手段と、
前記周波数別推定方向情報に基づいて波源方向ごとの前記推定方向情報に対する周波数別の前記寄与率を計算する周波数寄与率計算手段とを有する付記4に記載の波源方向推定装置。
(付記6)
前記重み計算手段は、
複数の前記ペア間において前記寄与率の類似度を計算し、算出した前記類似度に基づいて波源方向ごとの重みを計算する付記4または5に記載の波源方向推定装置。
(付記7)
前記統合推定方向情報算出手段は、
前記信号選択手段によって選択された全ての前記ペアのうち1組の前記ペアをなす二つの前記入力信号を取得するとともに、取得した前記ペアの前記相対遅延時間を取得し、前記ペアをなす二つの前記入力信号を用いて算出される周波数別相互相関関数と前記ペアの前記相対遅延時間とを用いて、前記ペアごとの周波数別推定方向情報を計算する複数の周波数別推定方向情報生成手段と、
複数の前記周波数別推定方向情報生成手段から前記ペアごとの前記周波数別推定方向情報を取得し、取得した前記ペアごとの前記周波数別推定方向情報に基づいて波源方向候補を計算し、算出した前記波源方向候補ごとに前記重みを計算する重み生成手段と、
複数の前記周波数別推定方向情報生成手段から前記ペアごとの前記周波数別推定方向情報を取得するとともに、前記重み生成手段から前記重みを取得し、複数の前記ペアごとの前記周波数別推定方向情報に前記重みをつけて統合して前記統合推定方向情報を生成する統合手段とを有する付記1乃至3のいずれか一項に記載の波源方向推定装置。
(付記8)
前記周波数別推定方向情報生成手段は、
前記ペアをなす二つの前記入力信号を周波数領域の変換信号に変換する変換手段と、
前記変換手段によって変換された前記変換信号を用いてクロススペクトルを計算するクロススペクトル計算手段と、
前記クロススペクトル計算手段によって算出される前記クロススペクトルを用いて平均クロススペクトルを計算する平均計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルを用いて分散を計算する分散計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルと、前記分散計算手段によって算出される分散とを用いて周波数別クロススペクトルを計算する周波数別クロススペクトル生成手段と、
前記周波数別クロススペクトル生成手段によって算出される前記周波数別クロススペクトルを逆変換して周波数別相互相関関数を計算する逆変換手段と、
前記ペアの前記周波数別相互相関関数と前記相対遅延時間とに基づいて波源方向と相関値との対応関係である周波数別推定方向情報を計算する周波数別推定方向情報計算手段とを有する付記7に記載の波源方向推定装置。
(付記9)
前記重み生成手段は、
前記周波数別推定方向情報を統合して前記ペアごとの前記推定方向情報を計算する複数の周波数統合手段と、
前記ペアごとの前記推定方向情報からピークを検出し、ピークとなる方向を前記ペアごとの波源方向候補として選出する複数の波源方向候補選出手段と、
複数の前記波源方向候補選出手段から前記ペアごとの前記波源方向候補を取得し、取得した前記ペアごとの前記波源方向候補を統合する波源方向候補統合手段と、
統合された前記波源方向候補と前記周波数別推定方向情報とに基づいて、前記波源方向候補ごとの前記推定方向情報に対する周波数別の前記寄与率を計算する複数の波源方向候補別周波数寄与率計算手段と、
前記波源方向候補ごとの前記推定方向情報に対する周波数別の前記寄与率を複数の前記波源方向候補別周波数寄与率計算手段から取得し、取得した前記波源方向候補ごとの前記推定方向情報に対する周波数別の前記寄与率に基づいて前記波源方向候補ごとの前記重みを計算する重み計算手段とを有する付記8に記載の波源方向推定装置。
(付記10)
前記重み計算手段は、
複数の前記ペア間の前記波源方向候補ごとの周波数別の前記寄与率の類似度を計算し、算出した前記類似度に基づいて前記波源方向候補ごとの前記重みを求める付記9に記載の波源方向推定装置。
(付記11)
前記統合推定方向情報算出手段によって算出される前記統合推定方向情報に基づいて波源方向を算出する波源方向算出手段を備える付記1乃至10のいずれか一項に記載の波源方向推定装置。
(付記12)
前記波源方向算出手段は、
一定時間ごとに前記統合推定方向情報が最大になる時刻に対応する方向を波源方向として算出する付記11に記載の波源方向推定装置。
(付記13)
複数の前記入力手段のそれぞれに対応して配置される前記センサを備える付記1乃至12のいずれか一項に記載の波源方向推定装置。
(付記14)
情報処理装置が、
複数のセンサによって検出される波動に基づく電気信号を入力信号として取得し、
複数の前記入力信号のうち二つを組み合わせたペアを複数組選択し、
前記ペアをなす二つの前記入力信号の供給元の前記センサへの波源方向ごとの波動の到達時間差を相対遅延時間として算出し、
前記ペアをなす前記入力信号と前記ペアの相対遅延時間とを用いて周波数別推定方向情報を前記ペアごとに生成し、
生成した前記ペアごとの前記周波数別推定方向情報を用いて波源方向ごとの推定方向情報に対する周波数別の寄与率を算出し、
算出した前記寄与率に基づいて異なる前記ペア間の前記寄与率の類似度に応じた波源方向ごとの重みを計算し、
全ての前記ペアの前記推定方向情報に前記重みをつけて統合することによって統合推定方向情報を生成する波源方向推定方法。
(付記15)
複数のセンサによって検出される波動に基づく電気信号を入力信号として取得する処理と、
複数の前記入力信号のうち二つを組み合わせたペアを複数組選択する処理と、
前記ペアをなす二つの前記入力信号の供給元の前記センサへの波源方向ごとの波動の到達時間差を相対遅延時間として算出する処理と、
前記ペアをなす前記入力信号と前記ペアの相対遅延時間とを用いて周波数別推定方向情報を前記ペアごとに生成する処理と、
生成した前記ペアごとの前記周波数別推定方向情報を用いて波源方向ごとの推定方向情報に対する周波数別の寄与率を算出する処理と、
算出した前記寄与率に基づいて異なる前記ペア間の前記寄与率の類似度に応じた波源方向ごとの重みを計算する処理と、
全ての前記ペアの前記推定方向情報に前記重みをつけて統合することによって統合推定方向情報を生成する処理とをコンピュータに実行させるプログラムを記録させたプログラム記録媒体。
11 入力端子
12 信号選択部
13 相対遅延時間計算部
15 推定方向情報・周波数寄与率生成部
16 重み計算部
17 統合部
20 波源方向推定装置
21 入力端子
22 信号選択部
23 相対遅延時間計算部
25 周波数別推定方向情報生成部
26 重み生成部
27 統合部
30 波源方向推定装置
31 入力端子
32 信号選択部
33 相対遅延時間計算部
35 統合推定方向情報算出部
41 入力端子
42 信号選択部
43 相対遅延時間計算部
45 統合推定方向情報算出部
46 波源方向算出部
151 変換部
152 クロススペクトル計算部
153 平均計算部
154 分散計算部
155 周波数別クロススペクトル生成部
156 逆変換部
157 周波数別推定方向情報計算部
158 周波数寄与率計算部
159 周波数統合部
251 変換部
252 クロススペクトル計算部
253 平均計算部
254 分散計算部
255 周波数別クロススペクトル生成部
256 逆変換部
257 周波数別推定方向情報計算部
261 周波数統合部
262 波源方向候補選出部
263 波源方向候補統合部
264 波源方向候補別周波数寄与率計算部
265 重み計算部
551 周波数別基本クロススペクトル計算部
552 カーネル関数スペクトル生成部
553 乗算部
Claims (15)
- 複数のセンサによって検出される波動に基づく電気信号を入力信号として取得する複数の入力手段と、
複数の前記入力信号のうち二つを組み合わせたペアを複数組選択する信号選択手段と、
前記ペアをなす二つの前記入力信号の供給元の前記センサへの波源方向ごとの波動の到達時間差を相対遅延時間として算出する相対遅延時間計算手段と、
前記ペアをなす前記入力信号と前記ペアの相対遅延時間とを用いて周波数別推定方向情報を前記ペアごとに生成し、生成した前記ペアごとの前記周波数別推定方向情報を用いて波源方向ごとの推定方向情報に対する周波数別の寄与率を算出し、算出した前記寄与率に基づいて異なる前記ペア間の前記寄与率の類似度に応じた波源方向ごとの重みを計算し、全ての前記ペアの前記推定方向情報に前記重みをつけて統合することによって統合推定方向情報を生成する統合推定方向情報算出手段とを備える波源方向推定装置。 - 前記信号選択手段は、
前記センサの間隔に基づいて、複数の前記入力信号から二つの前記入力信号を組み合わせた前記ペアを選択する請求項1に記載の波源方向推定装置。 - 前記相対遅延時間計算手段は、
一つの前記ペアをなす二つの前記入力信号の供給元である前記センサにおける波源探索方向を基準として、前記信号選択手段によって選択された全ての前記ペアの前記相対遅延時間を、基準である前記波源探索方向の関数として計算する請求項1または2に記載の波源方向推定装置。 - 前記統合推定方向情報算出手段は、
前記信号選択手段によって選択された全ての前記ペアのうち1組の前記ペアをなす二つの前記入力信号を取得するとともに、取得した前記ペアの前記相対遅延時間を取得し、前記ペアをなす二つの前記入力信号間の周波数別相互相関関数と前記ペアの前記相対遅延時間とを用いて、前記ペアをなす二つの前記入力信号間の前記推定方向情報と周波数別の前記寄与率とを生成する複数の推定方向情報・周波数寄与率生成手段と、
複数の前記推定方向情報・周波数寄与率生成手段から周波数ごとの前記寄与率を取得し、異なる前記ペア間の前記寄与率の類似度に応じて波源方向ごとの前記重みを計算する重み計算手段と、
複数の前記推定方向情報・周波数寄与率生成手段から前記推定方向情報を取得するとともに、前記重み計算手段から前記重みを取得し、複数の前記ペアごとの前記推定方向情報に前記重みをつけて統合して前記統合推定方向情報を生成する統合手段とを有する請求項1乃至3のいずれか一項に記載の波源方向推定装置。 - 前記推定方向情報・周波数寄与率生成手段は、
前記ペアをなす二つの前記入力信号を周波数領域の変換信号に変換する変換手段と、
前記変換手段によって変換された前記変換信号を用いてクロススペクトルを計算するクロススペクトル計算手段と、
前記クロススペクトル計算手段によって算出される前記クロススペクトルを用いて平均クロススペクトルを計算する平均計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルを用いて分散を計算する分散計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルと、前記分散計算手段によって算出される分散とを用いて周波数別クロススペクトルを計算する周波数別クロススペクトル生成手段と、
前記周波数別クロススペクトル生成手段によって算出される前記周波数別クロススペクトルを逆変換して周波数別相互相関関数を計算する逆変換手段と、
前記ペアの前記周波数別相互相関関数と前記相対遅延時間とに基づいて波源方向と相関値との対応関係である周波数別推定方向情報を計算する周波数別推定方向情報計算手段と、
前記周波数別推定方向情報を統合して前記推定方向情報を計算する周波数統合手段と、
前記周波数別推定方向情報に基づいて波源方向ごとの前記推定方向情報に対する周波数別の前記寄与率を計算する周波数寄与率計算手段とを有する請求項4に記載の波源方向推定装置。 - 前記重み計算手段は、
複数の前記ペア間において前記寄与率の類似度を計算し、算出した前記類似度に基づいて波源方向ごとの重みを計算する請求項4または5に記載の波源方向推定装置。 - 前記統合推定方向情報算出手段は、
前記信号選択手段によって選択された全ての前記ペアのうち1組の前記ペアをなす二つの前記入力信号を取得するとともに、取得した前記ペアの前記相対遅延時間を取得し、前記ペアをなす二つの前記入力信号を用いて算出される周波数別相互相関関数と前記ペアの前記相対遅延時間とを用いて、前記ペアごとの周波数別推定方向情報を計算する複数の周波数別推定方向情報生成手段と、
複数の前記周波数別推定方向情報生成手段から前記ペアごとの前記周波数別推定方向情報を取得し、取得した前記ペアごとの前記周波数別推定方向情報に基づいて波源方向候補を計算し、算出した前記波源方向候補ごとに前記重みを計算する重み生成手段と、
複数の前記周波数別推定方向情報生成手段から前記ペアごとの前記周波数別推定方向情報を取得するとともに、前記重み生成手段から前記重みを取得し、複数の前記ペアごとの前記周波数別推定方向情報に前記重みをつけて統合して前記統合推定方向情報を生成する統合手段とを有する請求項1乃至3のいずれか一項に記載の波源方向推定装置。 - 前記周波数別推定方向情報生成手段は、
前記ペアをなす二つの前記入力信号を周波数領域の変換信号に変換する変換手段と、
前記変換手段によって変換された前記変換信号を用いてクロススペクトルを計算するクロススペクトル計算手段と、
前記クロススペクトル計算手段によって算出される前記クロススペクトルを用いて平均クロススペクトルを計算する平均計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルを用いて分散を計算する分散計算手段と、
前記平均計算手段によって算出される前記平均クロススペクトルと、前記分散計算手段によって算出される分散とを用いて周波数別クロススペクトルを計算する周波数別クロススペクトル生成手段と、
前記周波数別クロススペクトル生成手段によって算出される前記周波数別クロススペクトルを逆変換して周波数別相互相関関数を計算する逆変換手段と、
前記ペアの前記周波数別相互相関関数と前記相対遅延時間とに基づいて波源方向と相関値との対応関係である周波数別推定方向情報を計算する周波数別推定方向情報計算手段とを有する請求項7に記載の波源方向推定装置。 - 前記重み生成手段は、
前記周波数別推定方向情報を統合して前記ペアごとの前記推定方向情報を計算する複数の周波数統合手段と、
前記ペアごとの前記推定方向情報からピークを検出し、ピークとなる方向を前記ペアごとの波源方向候補として選出する複数の波源方向候補選出手段と、
複数の前記波源方向候補選出手段から前記ペアごとの前記波源方向候補を取得し、取得した前記ペアごとの前記波源方向候補を統合する波源方向候補統合手段と、
統合された前記波源方向候補と前記周波数別推定方向情報とに基づいて、前記波源方向候補ごとの前記推定方向情報に対する周波数別の前記寄与率を計算する複数の波源方向候補別周波数寄与率計算手段と、
前記波源方向候補ごとの前記推定方向情報に対する周波数別の前記寄与率を複数の前記波源方向候補別周波数寄与率計算手段から取得し、取得した前記波源方向候補ごとの前記推定方向情報に対する周波数別の前記寄与率に基づいて前記波源方向候補ごとの前記重みを計算する重み計算手段とを有する請求項8に記載の波源方向推定装置。 - 前記重み計算手段は、
複数の前記ペア間の前記波源方向候補ごとの周波数別の前記寄与率の類似度を計算し、算出した前記類似度に基づいて前記波源方向候補ごとの前記重みを求める請求項9に記載の波源方向推定装置。 - 前記統合推定方向情報算出手段によって算出される前記統合推定方向情報に基づいて波源方向を算出する波源方向算出手段を備える請求項1乃至10のいずれか一項に記載の波源方向推定装置。
- 前記波源方向算出手段は、
一定時間ごとに前記統合推定方向情報が最大になる時刻に対応する方向を波源方向として算出する請求項11に記載の波源方向推定装置。 - 複数の前記入力手段のそれぞれに対応して配置される前記センサを備える請求項1乃至12のいずれか一項に記載の波源方向推定装置。
- 情報処理装置が、
複数のセンサによって検出される波動に基づく電気信号を入力信号として取得し、
複数の前記入力信号のうち二つを組み合わせたペアを複数組選択し、
前記ペアをなす二つの前記入力信号の供給元の前記センサへの波源方向ごとの波動の到達時間差を相対遅延時間として算出し、
前記ペアをなす前記入力信号と前記ペアの相対遅延時間とを用いて周波数別推定方向情報を前記ペアごとに生成し、
生成した前記ペアごとの前記周波数別推定方向情報を用いて波源方向ごとの推定方向情報に対する周波数別の寄与率を算出し、
算出した前記寄与率に基づいて異なる前記ペア間の前記寄与率の類似度に応じた波源方向ごとの重みを計算し、
全ての前記ペアの前記推定方向情報に前記重みをつけて統合することによって統合推定方向情報を生成する波源方向推定方法。 - 複数のセンサによって検出される波動に基づく電気信号を入力信号として取得する処理と、
複数の前記入力信号のうち二つを組み合わせたペアを複数組選択する処理と、
前記ペアをなす二つの前記入力信号の供給元の前記センサへの波源方向ごとの波動の到達時間差を相対遅延時間として算出する処理と、
前記ペアをなす前記入力信号と前記ペアの相対遅延時間とを用いて周波数別推定方向情報を前記ペアごとに生成する処理と、
生成した前記ペアごとの前記周波数別推定方向情報を用いて波源方向ごとの推定方向情報に対する周波数別の寄与率を算出する処理と、
算出した前記寄与率に基づいて異なる前記ペア間の前記寄与率の類似度に応じた波源方向ごとの重みを計算する処理と、
全ての前記ペアの前記推定方向情報に前記重みをつけて統合することによって統合推定方向情報を生成する処理とをコンピュータに実行させるプログラムを記録させたプログラム記録媒体。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2018/023971 WO2020003343A1 (ja) | 2018-06-25 | 2018-06-25 | 波源方向推定装置、波源方向推定方法、およびプログラム記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2020003343A1 true JPWO2020003343A1 (ja) | 2021-06-24 |
JP6933303B2 JP6933303B2 (ja) | 2021-09-08 |
Family
ID=68986143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020526725A Active JP6933303B2 (ja) | 2018-06-25 | 2018-06-25 | 波源方向推定装置、波源方向推定方法、およびプログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US11408963B2 (ja) |
JP (1) | JP6933303B2 (ja) |
WO (1) | WO2020003343A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220342026A1 (en) * | 2019-09-02 | 2022-10-27 | Nec Corporation | Wave source direction estimation device, wave source direction estimation method, and program recording medium |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008089312A (ja) * | 2006-09-29 | 2008-04-17 | Kddi Corp | 信号到来方向推定装置及び方法、並びに信号分離装置及び方法、コンピュータプログラム |
JP2010175431A (ja) * | 2009-01-30 | 2010-08-12 | Nippon Telegr & Teleph Corp <Ntt> | 音源方向推定装置とその方法と、プログラム |
JP2010193323A (ja) * | 2009-02-19 | 2010-09-02 | Casio Hitachi Mobile Communications Co Ltd | 録音装置、再生装置、録音方法、再生方法、及びコンピュータプログラム |
JP2011071686A (ja) * | 2009-09-25 | 2011-04-07 | Nec Corp | 映像音響処理装置、映像音響処理方法及びプログラム |
JP2015161551A (ja) * | 2014-02-26 | 2015-09-07 | 株式会社東芝 | 音源方向推定装置、音源方向推定方法およびプログラム |
WO2018003158A1 (ja) * | 2016-06-29 | 2018-01-04 | 日本電気株式会社 | 相関関数生成装置、相関関数生成方法、相関関数生成プログラムおよび波源方向推定装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4234746B2 (ja) * | 2006-09-25 | 2009-03-04 | 株式会社東芝 | 音響信号処理装置、音響信号処理方法及び音響信号処理プログラム |
JP6413741B2 (ja) * | 2014-12-16 | 2018-10-31 | 日本電気株式会社 | 振動発生源推定装置、方法およびプログラム |
JP6874766B2 (ja) * | 2016-06-15 | 2021-05-19 | 日本電気株式会社 | 波源方向推定装置、波源方向推定システム、波源方向推定方法および波源方向推定プログラム |
-
2018
- 2018-06-25 WO PCT/JP2018/023971 patent/WO2020003343A1/ja active Application Filing
- 2018-06-25 US US17/252,405 patent/US11408963B2/en active Active
- 2018-06-25 JP JP2020526725A patent/JP6933303B2/ja active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008089312A (ja) * | 2006-09-29 | 2008-04-17 | Kddi Corp | 信号到来方向推定装置及び方法、並びに信号分離装置及び方法、コンピュータプログラム |
JP2010175431A (ja) * | 2009-01-30 | 2010-08-12 | Nippon Telegr & Teleph Corp <Ntt> | 音源方向推定装置とその方法と、プログラム |
JP2010193323A (ja) * | 2009-02-19 | 2010-09-02 | Casio Hitachi Mobile Communications Co Ltd | 録音装置、再生装置、録音方法、再生方法、及びコンピュータプログラム |
JP2011071686A (ja) * | 2009-09-25 | 2011-04-07 | Nec Corp | 映像音響処理装置、映像音響処理方法及びプログラム |
JP2015161551A (ja) * | 2014-02-26 | 2015-09-07 | 株式会社東芝 | 音源方向推定装置、音源方向推定方法およびプログラム |
WO2018003158A1 (ja) * | 2016-06-29 | 2018-01-04 | 日本電気株式会社 | 相関関数生成装置、相関関数生成方法、相関関数生成プログラムおよび波源方向推定装置 |
Non-Patent Citations (2)
Title |
---|
KATO, MASANORI ET AL.: ""TDOA Estimation Based on Phase-Voting Cross Correlation and Circular Standard Deviation"", 2017 25TH EUROPEAN SIGNAL PROCESSING CONFERENCE(EUSIPCO), JPN6018034734, 2017, pages 1230 - 1234, XP033236133, ISSN: 0004554491, DOI: 10.23919/EUSIPCO.2017.8081404 * |
加藤 正徳 ほか: ""位相投票による相互相関と円周標準偏差に基づくTDOA推定"", 電子情報通信学会技術研究報告 IEICE TECHNICAL REPORT, vol. 第116巻, 第477号, JPN6018034733, March 2017 (2017-03-01), JP, pages 159 - 164, ISSN: 0004554490 * |
Also Published As
Publication number | Publication date |
---|---|
US20210263126A1 (en) | 2021-08-26 |
US11408963B2 (en) | 2022-08-09 |
WO2020003343A1 (ja) | 2020-01-02 |
JP6933303B2 (ja) | 2021-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9182475B2 (en) | Sound source signal filtering apparatus based on calculated distance between microphone and sound source | |
JP6289936B2 (ja) | 音源方向推定装置、音源方向推定方法およびプログラム | |
EP3484184A1 (en) | Acoustic field formation device, method, and program | |
KR101349268B1 (ko) | 마이크로폰 어레이를 이용한 음원 거리 측정 장치 | |
WO2013093565A1 (en) | Spatial audio processing apparatus | |
KR20170053623A (ko) | 사운드 소스들을 향상시키기 위한 방법 및 장치 | |
Padois | Acoustic source localization based on the generalized cross-correlation and the generalized mean with few microphones | |
Padois et al. | On the use of modified phase transform weighting functions for acoustic imaging with the generalized cross correlation | |
JPWO2018116943A1 (ja) | ノイズ抑圧装置、ノイズ抑圧方法、及びこれらを用いた受信装置、受信方法 | |
CN103688187A (zh) | 使用相位谱的声音源定位 | |
JP6862799B2 (ja) | 信号処理装置、方位算出方法及び方位算出プログラム | |
JP6933303B2 (ja) | 波源方向推定装置、波源方向推定方法、およびプログラム | |
Padois et al. | On the use of geometric and harmonic means with the generalized cross-correlation in the time domain to improve noise source maps | |
JP5986966B2 (ja) | 音場収音再生装置、方法及びプログラム | |
Schober et al. | Sound source localization using stochastic computing | |
JP7056739B2 (ja) | 波源方向推定装置、波源方向推定方法、およびプログラム | |
JP7286896B2 (ja) | 音源分離システム、音源位置推定システム、音源分離方法および音源分離プログラム | |
JP2005062096A (ja) | 話者位置検出方法、装置、プログラム、および記録媒体 | |
JP7276469B2 (ja) | 波源方向推定装置、波源方向推定方法、およびプログラム | |
JP6438786B2 (ja) | 信号中に含まれる正弦成分抽出装置、正弦成分抽出方法及びプログラム | |
JP2018170640A (ja) | マルチトーン信号検出装置、マルチトーン信号検出方法及びプログラム | |
JP5713933B2 (ja) | 音源距離測定装置、音響直間比推定装置、雑音除去装置、それらの方法、及びプログラム | |
US11276388B2 (en) | Beamforming system based on delay distribution model using high frequency phase difference | |
JP6585547B2 (ja) | 撮像システム、撮像方法 | |
WO2017154723A1 (ja) | 収録データ処理方法および収録データ処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201118 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20201118 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210720 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210802 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6933303 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |