(57)【特許請求の範囲】
【請求項１】
ＤＯＡ推定値（１０５）から曖昧性を解消するための装置（１００）であって、
バイアスされたＤＯＡ推定値とバイアスされていないＤＯＡ推定値との間の関係を表わすバイアス情報（１０１）を使用することにより、前記ＤＯＡ推定値（１０５）を分析して複数の曖昧な分析パラメータ（１１５）を得るＤＯＡ推定値分析部（１１０）と、
前記複数の曖昧な分析パラメータ（１１５）における曖昧性を解消して、曖昧でない解消済みパラメータ（１２５)を得る曖昧性解消部（１２０）と、を備える装置。
【請求項２】
請求項１に記載の装置（２００）において、
前記ＤＯＡ推定値分析部（２１０）は、バイアス関数（２０１）に基づいて前記ＤＯＡ推定値（１０５）を複数の曖昧なＤＯＡ推定値（２１５）に関連付け、
前記バイアス関数（２０１）は前記バイアス情報（１０１）によって示され、
前記複数の曖昧なＤＯＡ推定値（２１５）は前記複数の曖昧な分析パラメータ（１１５）に対応しており、
前記曖昧性解消部（１２０）は、一つ又は複数の活性音源の先験的な空間情報（２１１）を使用することによって前記複数の曖昧なＤＯＡ推定値（２１５）から前記曖昧でない解消済みパラメータ（１２５）を決定する、装置。
【請求項３】
請求項１又は２に記載の装置（３００）において、
前記曖昧性解消部（１２０）は、複数の曖昧なＤＯＡ推定値（２１５）から曖昧でない解消済みパラメータ（３２５）を選択するセレクタ（３２０）として構成され、
前記セレクタ（３２０）が、前記複数の曖昧なＤＯＡ推定値（２１５）のうちの各々の曖昧なＤＯＡ推定値に対応する定位推定値と、一つ又は複数の活性音源の先験的な定位推定値（３１１）との間の距離の大きさに基づいて、前記曖昧でない解消済みパラメータ（３２５）の選択を実行する、装置。
【請求項４】
請求項１又は２に記載の装置（４００）において、
前記曖昧性解消部（１２０）は、複数の曖昧なＤＯＡ推定値（２１５）を前記曖昧でない解消済みパラメータ（４２５）へとマッピングするマッパー（４２０）として構成され、
前記マッパー（４２０）が、前記複数の曖昧なＤＯＡ推定値（２１５）の各々の曖昧なＤＯＡ推定値に対応する定位推定値、又は複数の曖昧なＤＯＡ推定値（２１５）から選択された選択済みＤＯＡ推定値（５１５）に対応する選択済み定位推定値と、一つ又は複数の活性音源の先験的な定位推定値（３１１）との間の距離の大きさに基づいて、前記複数の曖昧なＤＯＡ推定値（２１５）のマッピングを実行する、装置。
【請求項５】
請求項４に記載の装置（４００）において、
前記マッパー（４２０）は、時間につれて変数（５４５）を生成する数字発生器（５４０）と、いくつかの前記変数（５４５）と一つ又は複数の活性音源の所定のＤＯＡ（５２５）に対応する所定の先験的な定位推定値とを結合して修正されたＤＯＡ関連値（５３５）を得る結合部（５３０）とを備えており、
前記曖昧でない解消済みパラメータ（１２５）が前記修正されたＤＯＡ関連値（５３５）に相当する、装置。
【請求項６】
請求項２〜５のいずれか一項に記載の装置（６００）において、
前記装置は、一つ又は複数の活性音源の先験的な空間情報（２１１）を提供する空間情報供給部（６１０）をさらに備えており、
前記空間情報供給部（６１０）は、前記ＤＯＡ推定値（１０５）とは異なる入力パラメータ（６１１）に基づいて前記先験的な空間情報（２１１）を提供する、装置。
【請求項７】
請求項１に記載の装置（７００）において、
前記ＤＯＡ推定値分析部（７１０）は、前記複数の曖昧な分析パラメータ（１１５）の各々の曖昧な分析パラメータについて複数の分枝（７１５）を生成し、
前記複数の分枝（７１５）は前記複数の曖昧な分析パラメータ（１１５）における曖昧性を反映しており、
前記ＤＯＡ推定値分析部（７１０）は、前記ＤＯＡ推定値（１０５）から前記複数の分枝（７１５）を直接取得するか、又は複数の曖昧なＤＯＡ推定値（８１３）をさらに処理して前記複数の分枝（７１５）を取得し、
前記ＤＯＡ推定値分析部（７１０）は、スペクトル重み付け操作のための利得パラメータを計算し、
前記曖昧性解消部（７２０）は、前記曖昧な分析パラメータ（１１５）を表わす複数の分枝（７１５）から曖昧でない解消済みパラメータ（７２５）を決定する、装置。
【請求項８】
請求項７に記載の装置（７００）において、
前記曖昧性解消部（７２０）は、前記複数の分枝（７１５）によって表わされた前記複数の曖昧な分析パラメータ（１１５）を結合して前記曖昧でない解消済みパラメータ（７２５）を表わす結合済みパラメータを得る、結合部として構成されている、装置。
【請求項９】
請求項７に記載の装置（７００）において、
前記曖昧性解消部（７２０）は、前記複数の分枝（７１５）によって表わされた複数の曖昧な分析パラメータ（１１５）から曖昧でない解消済みパラメータ（７２５）を選択するセレクタとして構成されている、装置。
【請求項１０】
請求項７〜９のいずれか一項に記載の装置（８００）において、
前記ＤＯＡ推定値分析部（８１０）は、前記バイアス情報（１０１）によって示されるバイアス関数（２０１）に基づいて前記ＤＯＡ推定値（１０５）を複数の曖昧なＤＯＡ推定値（８１３）に関連付け、前記複数の曖昧なＤＯＡ推定値（８１３）をさらに処理して前記曖昧な分析パラメータ（１１５）を表わす前記複数の分枝（８１５）を得る、装置。
【請求項１１】
請求項７〜９のいずれか一項に記載の装置（９００）において、
前記ＤＯＡ推定値分析部（９１０）及び前記曖昧性解消部（７２０）は、少なくとも２つの処理工程が連続的に実行された場合に得られるであろう前記曖昧でない解消済みパラメータ（７２５）を表わす処理結果が提供されるように、単一の処理工程で前記曖昧でない解消済みパラメータ（７２５）を得る、装置。
【請求項１２】
請求項７に記載の装置（９００）において、
前記ＤＯＡ推定値分析部（９１０）は、前記複数の分枝（９１５）の各々の分枝について、前記ＤＯＡ推定値（１０５）を対応する曖昧な分析パラメータ（９１５）に関連付けるパラメータ処理ユニット（９１４）を備えている、装置。
【請求項１３】
請求項１〜１２のいずれか一項に記載の装置において、
前記装置は、空間オーディオ信号から前記ＤＯＡ推定値（１０５）を推定するためのＤＯＡ推定器（１１）をさらに備えており、前記空間オーディオ信号は無指向性の信号（Ｗ）及び複数の異なる指向性ダイポール信号を含んでいる、装置。
【請求項１４】
ＤＯＡ推定値（１０５）から曖昧性を解消するための方法であって、
バイアスされたＤＯＡ推定値とバイアスされていないＤＯＡ推定値との間の関係を表わすバイアス情報（１０１）を使用することにより、前記ＤＯＡ推定値（１０５）を分析して複数の曖昧な分析パラメータ（１１５）を得るステップと、
前記複数の曖昧な分析パラメータ（１１５）における曖昧性を解消して、曖昧でない解消済みパラメータ（１２５)を得るステップと、を含む方法。
【請求項１５】
コンピュータ上で実行されたときに、請求項１４に記載の方法を実行するためのプログラムコードを有するコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、到来方向推定値から曖昧性を解消する装置及び方法に関する。本発明のさらなる実施形態は、曖昧な空間オーディオパラメータ推定の強化に関する。
【背景技術】
【０００２】
本発明は、マイクロホン・フロントエンドを用いた空間オーディオ符号化に関する。本発明は、パラメトリック技法、即ちパラメトリック・サイド情報を伴う一つ又は複数のオーディオ信号によって空間サウンドが表現される技法に焦点を当てるものである。サイド情報は、少なくとも二つのマイクロホンに基づいて得られる。マイクロホンの配置及び信号処理構成は、通常は所望の仕様が厳格に満足されるような対象となるある周波数範囲をカバーするように設計されている。この対象領域を超えると、必要な空間的選択性からの顕著な逸脱が生じる可能性がある。極端な逸脱は、いわゆる空間エイリアシングを含む。音場分析において空間エイリアシングの最も重大な影響は、サウンドの到来方向の推定値が曖昧になることである。結果として、一つの方向についてのみ期待又は所望されるある挙動が他の一つ又は複数の方向においても観測されることになる。
【０００３】
空間オーディオ符号化の目的は、パラメトリック・サイド情報を伴う少数のオーディオチャネルを用いて空間音響シーンを再構成することである。そのような再構成で生じる一般的な課題は、不正確なパラメータ推定値の場合に、空間パラメータを利用する再構成処理又は任意の他の（例えば方向性フィルタリングのような）処理ユニットにおいて望ましくない結果を引き起こすことである。このような懸念は、上述した問題の影響を受ける周波数範囲であって空間エイリアシングが発生する周波数範囲に特に当てはまる。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】V. Pulkki. Spatial sound reproduction with directional audio coding. J. Audio Eng. Soc., 55(6):503_516, June 2007
【非特許文献２】C. Faller. Microphone front-ends for spatial audio coders. In 125th AES Convention, San Francisco, USA, Oct. 2008
【非特許文献３】ISO/IEC 23003-1:2007. Information technology - MPEG Audio technologies - Part 1: MPEG Surround. International Standards Organization, Geneva, Switzerland, 2007
【非特許文献４】E. Benjamin and T. Chen. The native B-format microphone: Part 1. In 119th AES Convention, Paper 6621, New York, Oct. 2005
【非特許文献５】J. Merimaa. Applications of a 3-D microphone array. In 112th AES Convention, Paper 5501, Munich, May 2002
【非特許文献６】G. W. Elko. Superdirectional microphone arrays. In S. L. Gay and J. Benesty, editors, Acoustic Signal Processing for Telecommunication, chapter 10, pages 181_235. Kluwer Academic Publishers, 2000
【非特許文献７】M. Kallinger, F. Kuech, R. Schultz-Amling, G. Del Galdo, J. Ahonen, and V. Pulkki. Analysis and Adjustment of Planar Microphone Arrays for Application in Directional Audio Coding. In 124th AES Convention, Paper 7374, Amsterdam, the Netherlands, May 2008
【非特許文献８】M. Kallinger, H. Ochsenfeld, G. Del Galdo, F. Kuech, D. Mahne, R. Schultz-Amling, and O. Thiergart. A spatial filtering approach for directional audio coding. In 126th AES Convention, Paper 7653, Munich, Germany, May 2009
【非特許文献９】O. Thiergart, R. Schultz-Amling, G. Del Galdo, D. Mahne, and F. Kuech. Localization of sound sources in reverberant environments based on directional audio coding parameters. In 128th AES Convention, Paper 7853, New York City, NY, USA, Oct. 2009
【発明の概要】
【発明が解決しようとする課題】
【０００５】
したがって、本発明の目的は、空間パラメータ処理のための改善された概念を提供することにある。
【課題を解決するための手段】
【０００６】
この目的は、請求項１に記載の装置、請求項１４に記載の方法、又は請求項１５に記載のコンピュータプログラムによって達成される。
【０００７】
本発明の実施形態によれば、ＤＯＡ（到来方向）推定値から曖昧性を解消する装置は、ＤＯＡ推定値分析部と曖昧性解消部とを備える。ＤＯＡ推定値分析部は、バイアス情報を用いてＤＯＡ推定値を分析することで、複数の曖昧な分析パラメータを得るよう構成される。ここでバイアス情報とは、バイアスされたＤＯＡ推定値とバイアスされていないＤＯＡ推定値との間の関係を表わすものである。さらに曖昧性解消部は、複数の曖昧な分析パラメータの中の曖昧性を解消することで、曖昧でない解消済みパラメータを得るよう構成されている。
【０００８】
本発明の根底にある基本的な考え方は、バイアスされたＤＯＡ推定値とバイアスされていないＤＯＡ推定値との間の関係を表わすバイアス情報を使用してＤＯＡ推定値を分析することで、複数の曖昧な分析パラメータが得られ、さらにそれら複数の曖昧な分析パラメータの中の曖昧性が解消されて曖昧でない解消済みパラメータが得られる場合に、空間パラメータ処理の改善を達成できるというものである。この方法は、空間エイリアシングによって引き起される不正確なＤＯＡパラメータの推定値に起因する曖昧性を解消できる。
【０００９】
さらなる実施形態によれば、ＤＯＡ推定値分析部は、バイアス情報によって示されるバイアス関数に基づいて、ＤＯＡ推定値を複数の曖昧なバイアスされていないＤＯＡ推定値に関連付ける。ここで、複数の曖昧なバイアスされていないＤＯＡ推定値は、複数の曖昧な分析パラメータに対応している。さらに、曖昧性解消部は、一つ又は複数の活性音源（active sources）の先験的な定位情報（a priori localization information）を使用して、複数の曖昧なバイアスされていないＤＯＡ推定値から曖昧でない解消済みパラメータを決定する。これにより、音源の位置についての先験的な情報から特定の又は補正済みＤＯＡパラメータを得ることが実質的に可能となる。
【００１０】
さらなる実施形態によれば、ＤＯＡ推定値分析部は、複数の曖昧な分析パラメータの各パラメータについて複数の分枝を生成し、それら複数の分枝は、複数の曖昧な分析パラメータの曖昧性を反映している。ここでＤＯＡ推定値分析部は、ＤＯＡ推定値から複数の分枝を直接的に取得するか、又は複数の曖昧なバイアスされていないＤＯＡ推定値をさらに処理することで複数の分枝を取得する。さらにＤＯＡ推定値分析部は、スペクトル重み付け操作の利得パラメータを計算する。代替的にＤＯＡ推定値分析部は、例えば加算、減算又は対数をとるなど、スペクトル重み付け操作と同様の代替手段によってスペクトル修正を行っても良い。さらに曖昧性解消部は、曖昧な分析パラメータを表わす複数の分枝から曖昧でない解消済みパラメータを決定する。このようにして、方向性オーディオ符号化(directional audio coding)の分野におけるレンダリング又は方向性フィルタリング操作などのスペクトル重み付け操作のために、適切な利得パラメータを得ることができる。さらに、この方法では先験的な定位情報を取得する必要がない。
【００１１】
さらなる実施形態においては、ＤＯＡ推定値分析部は、複数の分枝を得るための単一の処理ステップが、仮に二つの処理ステップが連続的に実行された場合に得られるであろう処理結果をもたらすように構成されている。特に、ＤＯＡ推定値分析部は、複数の分枝の各々の分枝のためにパラメータ処理ユニットを備えても良く、そのパラメータ処理ユニットは、ＤＯＡ推定値を処理結果によって示される通りにある所定の曖昧な分析パラメータに関連付けても良い。このようにして、複数の曖昧な分析パラメータの曖昧性が解消される場合に、より高い処理効率又は演算量の削減を実現できる。
【００１２】
本発明のさらなる実施形態は、ＤＯＡ推定値から曖昧性を解消する方法を提供する。
【００１３】
以下に、本発明の実施形態を添付図面を参照して説明する。
【図面の簡単な説明】
【００１４】
【図１】ＤＯＡ推定値から曖昧性を解消する装置の一実施形態のブロック図を示す。
【図２】先験的な定位情報を使用してＤＯＡ推定値から曖昧性を解消する装置の他の実施形態のブロック図を示す。
【図３】曖昧でない解消済みパラメータを選択するためのセレクタを備えた、ＤＯＡ推定値から曖昧性を解消する装置の他の実施形態のブロック図を示す。
【図４】複数の曖昧なバイアスされていないＤＯＡ推定値を曖昧でない解消済みパラメータへとマッピングするためのマッパーを備えた、ＤＯＡ推定値から曖昧性を解消する装置の他の実施形態のブロック図を示す。
【図５】数字発生器を備えた図４に係るマッパーの実施形態のブロック図を示す。
【図６】先験的な定位情報を供給する定位情報供給部を備えた、ＤＯＡ推定値から曖昧性を解消する装置のさらなる実施形態のブロック図を示す。
【図７】利得パラメータ計算を備えた、ＤＯＡ推定値から曖昧性を解消する装置のさらなる実施形態のブロック図を示す。
【図８】曖昧なバイアスされていないＤＯＡ推定値の生成部とパラメータプロセッサとを備えた、ＤＯＡ推定値から曖昧性を解消する装置のさらなる実施形態のブロック図を示す。
【図９】単一ステップで処理結果を得るための複数のパラメータ処理ユニットを備えた、ＤＯＡ推定値から曖昧性を解消する装置のさらなる実施形態のブロック図を示す。
【図１０】パラメトリック空間オーディオ符号器の概略図を示す。
【図１１】パラメトリック空間オーディオ復号器の概略図を示す。
【図１２】第５の無指向性マイクロホンを中央に有する、無指向性マイクロホンの正方形アレイの概略図を示す。
【図１３】本発明の一実施形態に係る、４〜７ｋＨｚの間の周波数における方位角推定値の例示的なグラフを示す。
【図１４】本発明の一実施形態に係る、３個のスピーカを有する設定のためのＤＯＡの関数としてのパンニング利得の例示的なグラフを示す。
【図１５】本発明の一実施形態に係る、方向性フィルタリング利得関数の例示的な指向性パターンを示す。
【図１６】空間パラメータの推定値から曖昧性を解消する装置の一実施形態の全体的な構造を示す。
【図１７】本発明のある実施形態に係る、バイアスされたＤＯＡ推定値を補正する構成を示す。
【図１８】本発明のある実施形態に係る、空間エイリアシング周波数を超える周波数領域におけるバイアスされていないＤＯＡの関数としての、推定されたＤＯＡの例示的なグラフを示す。
【図１９】本発明のさらなる実施形態に係る、バイアスされていないＤＯＡの関数として推定されたＤＯＡの例示的なグラフである。
【図２０】本発明の実施形態に係る、それぞれの確率密度関数(probability density functions)を有する二つの活性音源の例示的なグラフである。
【図２１】本発明の実施形態に係る、ＤＯＡの関数としての方向性フィルタリングのための例示的な設計関数を示す。
【図２２】曖昧なＤＯＡの関数として作動する処理ユニットを調節して補正済み出力パラメータを生成する、本発明の一実施形態に係る構成を示す。
【発明を実施するための形態】
【００１５】
図１は、ＤＯＡ推定値１０５の曖昧性を解消する装置１００の実施形態のブロック図を示す。ＤＯＡ推定値１０５は、ＤＯＡ推定器を使用することによって空間オーディオ信号から取得されても良い。ここで、空間オーディオ信号は無指向性の信号と複数の異なる指向性のダイポール信号とを含んでいてもよい。空間オーディオ信号から得られるＤＯＡ推定値１０５は、ＤｉｒＡＣ（方向性オーディオ符号化）パラメータに対応しても良く、サウンドのＤＯＡに関する角度
0005449624.tif 000002
によって表わすことができる。特に、ＤＯＡ推定値１０５は時間及び周波数に依存しても良く、即ち
0005449624.tif 000003
であっても良く、ここでｋは時間ブロックを表わし、ｍが周波数帯域を表わす。さらに、ＤＯＡ推定値１０５、即ち
0005449624.tif 000004
は、一つ又は複数の活性音源に由来する不正確な又は曖昧なＤＯＡ推定値に実質的に対応しても良い。図１を参照すると、装置１００はＤＯＡ推定値分析部１１０と曖昧性解消部１２０とを備えている。図１から明らかなように、ＤＯＡ推定値分析部１１０は、バイアス情報１０１を用いてＤＯＡ推定値１０５を分析し、複数の曖昧な分析パラメータを取得する。ここで、バイアス情報１０１はバイアスされたＤＯＡ推定値
0005449624.tif 000005
とバイアスされていないＤＯＡ推定値φとの間の関係
0005449624.tif 000006
を表わしている。特に、ＤＯＡ推定値がバイアスされたＤＯＡ推定値又はバイアスされていないＤＯＡ推定値を表現しても良いことに留意すべきである。さらに、曖昧性解消部１２０は、複数の曖昧な分析パラメータ１１５の曖昧性を解消し、曖昧でない解消済みパラメータ１２５を得るように構成されている。
【００１６】
図２は、先験的な定位又は空間情報２１１を使用してＤＯＡ推定値１０５の曖昧性を解消する装置２００のさらなる実施形態のブロック図を示す。図２に示すように、装置２００はＤＯＡ推定値１０５、即ち
0005449624.tif 000007
を、バイアス関数２０１、即ち
0005449624.tif 000008
に基づいて、複数の曖昧なバイアスされていないＤＯＡ推定値２１５、即ち
0005449624.tif 000009
へと関連付ける、ＤＯＡ推定値分析部２１０を備えている。図２に示した装置２００のＤＯＡ推定値分析部２１０が図１に示した装置１００のＤＯＡ推定値分析部１１０に対応しても良い。ここで、図２のＤＯＡ推定値分析部２１０によって使用されるバイアス関数２０１が図１のバイアス情報１０１によって表わされても良い一方で、図２のＤＯＡ推定値分析部２１０の出力における複数の曖昧なバイアスされていないＤＯＡ推定値２１５、即ち
0005449624.tif 000010
が図１の複数の曖昧な分析パラメータ１１５に対応しても良い。図２の実施形態において、曖昧性解消部１２０は、一つ又は複数の活性音源の先験的な定位情報２１１を使用して、複数の曖昧なバイアスされていないＤＯＡ推定値２１５、即ち
0005449624.tif 000011
から、曖昧でない解消済みパラメータ１２５、即ち
0005449624.tif 000012
を決定する。ここで、先験的な定位情報は一つ又は複数の活性音源のＤＯＡ２１１、即ち
0005449624.tif 000013
を表わしてもよい。さらに、ＤＯＡ推定から得られるＤＯＡ推定値１０５、即ち
0005449624.tif 000014
が一つ又は複数の活性音源に由来しても良いことに留意すべきである。
【００１７】
ここで、一般に２種類のＤＯＡがあることに留意すべきである。いくつかの実施形態において、ＤＯＡの
0005449624.tif 000015
が個々の音源の位置に相当し得る一方で、ＤＯＡ推定値
0005449624.tif 000016
が「瞬時の」ＤＯＡ（即ち個々の時間／周波数タイルに関連するＤＯＡ）に相当し得る。
【００１８】
図３は、セレクタ３２０を備えた、ＤＯＡ推定値１０５の曖昧性を解消する装置３００のさらなる実施形態のブロック図を示す。ここで、図３の装置３００は、基本的には図２の装置２００と同じブロックを備えている。したがって、類似の構成及び／又は機能を有する同じブロックを同じ番号で示す。しかしながら、図３の実施形態においては、曖昧性解消部１２０は、複数の曖昧なバイアスされていないＤＯＡ推定値２１５から曖昧でない解消済みパラメータ３２５、即ち
0005449624.tif 000017
を選択するセレクタ３２０として構成されている。いくつかの実施形態においては、セレクタ３２０は、複数の曖昧なバイアスされていないＤＯＡ推定値２１５のうちの一つの曖昧なバイアスされていないＤＯＡ推定値に対応する定位推定値と、一つ又は複数の活性音源の先験的な定位推定値３１１と、の間の距離の大きさに基づいて、曖昧でない解消済みパラメータ３２５、即ち
0005449624.tif 000018
の選択を行っても良い。ここで、図３に示したセレクタ３２０の出力における曖昧でない解消済みパラメータ３２５は、図２に示した曖昧性解消部１２０の出力における曖昧でない解消済みパラメータ１２５に対応しても良い。図３に示した先験的な定位推定値３１１は、一つ又は複数の活性音源のＤＯＡ
0005449624.tif 000019
のいずれか一つに対応しても良い。したがって、いくつかの実施形態においては、曖昧でない解消済みパラメータ３２５、即ち
0005449624.tif 000020
は、複数の曖昧なバイアスされていないＤＯＡ推定値２１５から選択され、一つ又は複数の活性音源に対する距離の大きさが最小となるような選択されたＤＯＡ推定値
0005449624.tif 000021
に対応しても良い。
【００１９】
図４は、マッパー４２０を備えた、ＤＯＡ推定値１０５からの曖昧性を解消する装置４００のさらなる実施形態のブロック図を示す。図４の装置４００は、基本的に図２の装置２００と同じブロックを備えているため、類似の構成及び／又は機能を有する同じブロックは、この場合も先と同様に同じ番号で示されている。図４の実施形態においては、曖昧性解消部１２０が複数の曖昧なバイアスされていないＤＯＡ推定値２１５を曖昧でない解消済みパラメータ４２５、即ち
0005449624.tif 000022
へとマッピングするためのマッパー４２０として構成されている。特に、マッパー４２０は、複数の曖昧なバイアスされていないＤＯＡ推定値２１５のうちの一つの曖昧なバイアスされていないＤＯＡ推定値に対応する定位推定値、又は複数の曖昧なバイアスされていないＤＯＡ推定値２１５から選択された選択済みＤＯＡ推定値５１５、即ち
0005449624.tif 000023
に対応する選択済み定位推定値と、一つ又は複数の活性音源の先験的な定位推定値３１１と、の間の距離の大きさに基づいて、複数の曖昧なバイアスされていないＤＯＡ推定値２１５のマッピングを行なっても良い。図３の実施形態と同様に、先験的な定位推定値３１１は、一つ又は複数の活性音源のＤＯＡ
0005449624.tif 000024
のいずれか一つに対応していても良い。
【００２０】
図５は、数字発生器５４０を備えた、図４に係るマッパー４２０の実施形態のブロック図を示す。特に、図５にはマッパー４２０のさらなる詳細が示されている。図５から明らかなように、マッパー４２０は、
0005449624.tif 000025
と、
0005449624.tif 000026
と、結合部５３０とを備えていてもよい。
0005449624.tif 000027
は、複数の曖昧なバイアスされていないＤＯＡ推定値２１５から、一つ又は複数の活性音源に対する距離の大きさが最小となるバイアスされていないＤＯＡ推定値５１５、即ち
0005449624.tif 000028
を選択しても良い。
0005449624.tif 000029
は、選択されたバイアスされていないＤＯＡ推定値５１５、即ち
0005449624.tif 000030
に対する距離の大きさが最小である一つ又は複数の活性音源のＤＯＡ５２５、即ち
0005449624.tif 000031
を決定しても良い。ここで、曖昧でない解消済みパラメータ
0005449624.tif 000032
が一つ又は複数の活性音源の決定されたＤＯＡ５２５、即ち
0005449624.tif 000033
に相当しても良い。
【００２１】
図５に示すマッパー４２０の数字発生器５４０が時間につれて変数５４５を生成する一方で、図５に示すマッパー４２０の結合部５３０が、数字発生器５４０によって生成されるいくつかの変数５４５と、一つ又は複数の活性音源の決定されたＤＯＡ５２５、即ち
0005449624.tif 000034
に対応する決定された先験的な定位推定値とを結合し、修正されたＤＯＡ関連値５３５、即ち
0005449624.tif 000035
を取得しても良い。ここで、曖昧でない解消済みパラメータ１２５、即ち
0005449624.tif 000036
が修正されたＤＯＡ関連値５３５、即ち
0005449624.tif 000037
に相当しても良い。時間とともに変数５４５を生成するために使用される数字発生器５４０は、例えば線形フィードバックシフトレジスタなどの疑似乱数発生器として構成されても良い。
【００２２】
図６は、定位又は空間情報供給部６１０を備えた、ＤＯＡ推定値１０５から曖昧性を解消する装置６００のさらなる実施形態のブロック図を示す。ここで、図６の装置６００は、基本的に図２の装置２００と同じブロックを備えている。したがって、類似の構成及び／又は機能を有する同じブロックは同じ番号で示されている。図６に示すように、定位情報供給部６１０は一つ又は複数の活性音源の先験的な定位情報２１１をもたらしても良い。図６の実施形態に示すように、定位情報供給部６１０は、ＤＯＡ推定値１０５、即ち
0005449624.tif 000038
とは異なる入力パラメータ６１１に基づいて、先験的な定位情報２１１をもたらしても良い。これは、一つ又は複数の活性音源のＤＯＡ
0005449624.tif 000039
の形式を持つ先験的な定位情報２１１が、基本的にＤＯＡ推定値１０５、即ち
0005449624.tif 000040
とは異なるデータから導出されることを意味する。
【００２３】
図７は、利得パラメータ計算を伴う、ＤＯＡ推定値１０５の曖昧性を解消する装置７００のさらなる実施形態のブロック図を示す。図７から明らかなように、装置７００はＤＯＡ推定値分析部７１０を備え、この推定値分析部７１０が複数の曖昧な分析パラメータ１１５の各々の曖昧な分析パラメータのための複数の分枝７１５を生成し、それら複数の分枝７１５が、複数の曖昧な分析パラメータ１１５における曖昧性を反映するよう構成されてもよい。図７の実施形態において、ＤＯＡ推定値分析部７１０は、複数の分枝７１５をＤＯＡ推定値１０５から直接的に取得しても良く、又は複数の曖昧なバイアスされていないＤＯＡ推定値８１３をさらに処理することで複数の分枝７１５を取得しても良い。特にＤＯＡ推定値分析部７１０は、方向性オーディオ符号化におけるレンダリング又は方向性フィルタリングの操作などのスペクトル重み付け操作のための利得パラメータを計算する。図７に示すように、利得パラメータの計算を含むＤＯＡ推定値分析部７１０は、ＤＯＡ推定値１０５の曖昧性を考慮することによって複数の分枝７１５を生成するために、バイアス情報１０１を受け取ることができる。さらに図７に示すように、装置７００の曖昧性解消部７２０は、曖昧な分析パラメータ１１５を表わす複数の分枝７１５から曖昧でない解消済みパラメータ７２５、即ちｆresを決定しても良い。
【００２４】
さらなる実施形態によれば、装置７００の曖昧性解消部７２０は、複数の分枝７１５によって表わされる複数の曖昧な分析パラメータを結合し、曖昧でない解消済みパラメータ７２５、即ちｆresを表わす結合済みのパラメータを得る結合部として構成されても良い。
【００２５】
さらなる実施形態によれば、装置７００の曖昧性解消部７２０は、複数の分枝７１５によって表わされる複数の曖昧な分析パラメータ１１５から曖昧でない解消済みパラメータ７２５、即ちｆresを選択するセレクタとして構成することもできる。
【００２６】
特に、図７に示す装置７００に含まれる、利得パラメータの計算を実行するためのＤＯＡ推定値分析部７１０と結合部又はセレクタの形態を持つ曖昧性解消部７２０とは、図１に示す装置１００のＤＯＡ推定値分析部１１０と曖昧性解消部１２０とにそれぞれ対応することができる。
【００２７】
基本的に、図７の装置７００は、特に先験的な定位情報に基づく操作を必要とせずに、複数の曖昧な分析パラメータから曖昧でない解消済みパラメータを特定するために使用することができる。図７の実施形態において、結合部又はセレクタとして構成される曖昧性解消部７２０は、例えば利得パラメータの最小又は最大値を選択しても良く、又は平均値を計算することによって利得パラメータを結合させても良い。
【００２８】
図８は、曖昧なバイアスされていないＤＯＡ推定値の生成部８１２とパラメータプロセッサ８１４とを備えた、ＤＯＡ推定値１０５から曖昧性を解消する装置８００のさらなる実施形態のブロック図を示す。ここで、図８の装置８００は、基本的に図７の装置７００に対応し、類似の構成及び／又は機能を有する同じブロックは同じ番号によって示されている。図８に示すとおり、装置８００のＤＯＡ推定値分析部８１０は、曖昧なバイアスされていないＤＯＡ推定値の生成部８１２と、パラメータプロセッサ８１４とを備えることができる。ここで、曖昧なバイアスされていないＤＯＡ推定値の生成部８１２は、バイアス情報１０１によって示されるバイアス関数２０１に基づいて、ＤＯＡ推定値１０５を複数の曖昧なバイアスされていないＤＯＡ推定値８１３に関連付けても良い。さらに、パラメータプロセッサ８１４は、複数の曖昧なバイアスされていないＤＯＡ推定値８１３をさらに処理することで、曖昧な分析パラメータ１１５を表わす複数の分枝８１５を取得しても良い。いくつかの実施形態によれば、曖昧な分析パラメータ１１５を表わす複数の分枝８１５は、三方向パンニング利得パラメータなどのレンダリング作業に使用される利得パラメータ（ｇp）のセット、又はＤｉｒＡＣ符号化の仕組みにおける方向性フィルタリングに適用できる利得パラメータ（Ｄ）にそれぞれ対応することができる。いくつかの実施形態においては、曖昧性解消部７２０を使用することによって、異なる分枝８１５の対応するレンダリング作業のための利得パラメータｇp又は方向性フィルタリング操作のための利得パラメータＤが例えばそれらの平均値を得るように線形結合されても良く、又は特定の利得パラメータ（例えば最小値又は最大値）が選択されても良い。曖昧性解消部７２０の出力において、曖昧でない解消済みパラメータ７２５、即ちｆresが得られる。
【００２９】
図９は、複数のパラメータ処理ユニット９１３を有する、ＤＯＡ推定値１０５の曖昧性を解消する装置９００のさらなる実施形態のブロック図を示す。ここで、図９の装置９００は図７の装置７００に対応でき、類似の構成及び／又は機能を有する同じブロックは同じ番号によって示されている。具体的には、図９の実施形態においては、装置９００は単一の処理ステップで複数の分枝９１５を取得できるＤＯＡ推定値分析部９１０を備えており、その推定値分析部は、二つの処理ステップが連続して実行された場合に取得できるであろう曖昧な分析パラメータ１１５を表わす処理結果をもたらすようにしても良い。
【００３０】
ここで、ＤＯＡ推定値分析部９１０によって得られる処理結果が、基本的に図８に示したＤＯＡ推定値分析部８１０によって得られる処理結果又は複数の分枝８１５に対応することに留意すべきである。ここで、複数のパラメータ処理ユニット９１３によって得られる各処理済みパラメータ又は出力パラメータは、基本的にパラメータプロセッサ８１４によって出力される対応する利得パラメータ（ｇp；Ｄ）に等しい。しかしながら、図９の実施形態においては、これらの同等なパラメータを単一の処理ステップを用いることによって効率的に得ることができる。したがって、装置９００のＤＯＡ推定値分析部９１０を、中間の処理ステップを必要とせずに処理結果を得るよう強化されたやり方でＤＯＡ推定値１０５を処理するための強化型プロセッサと称することもできる。
【００３１】
さらに、図９の実施形態において、ＤＯＡ推定値分析部９１０及び曖昧性解消部７２０もまた、曖昧でない解消済みパラメータ７２５を単一の処理ステップで取得するよう構成され、少なくとも二つの処理ステップが連続して実行された場合に取得可能な曖昧でない解消済みパラメータ７２５を表わす処理結果をもたらしても良い。
【００３２】
いくつかの実施形態において、曖昧でない解消済みパラメータ７２５は、対応するＤＯＡ推定値１０５に依存して事前決定又は事前計算できる。そのようにして得られたパラメータ及び対応するＤＯＡ推定値は、ルックアップ・テーブルに挿入することができる。この事前決定された情報を含むルックアップ・テーブルは、実際の処理の際に生じる演算量を回避又は少なくとも削減するために、効率的に使用することができる。
【００３３】
特に、ＤＯＡ推定器を使用することによって空間オーディオ信号から推定されたＤＯＡ推定値１０５を、最初にルックアップ・テーブルにおける入力値とし、次いでルックアップ・テーブル内でこの入力値に関連付けられた対応する曖昧でない解消済みパラメータ７２５を出力値として出力することができる。
【００３４】
代替的に、図９を参照すると、装置９００のＤＯＡ推定値分析部９１０は複数の分枝９１５の各分枝についてパラメータ処理ユニットを備えていてもよい。特に、複数のパラメータ処理ユニット９１３の各パラメータ処理ユニット９１４は、ＤＯＡ推定値１０５を対応する曖昧な分析パラメータに関連付けるよう構成されても良い。さらに、装置９００の曖昧性解消部７２０は、次に曖昧な分析パラメータ９１５を表わす処理結果の曖昧性を解消し、曖昧でない解消済みパラメータ７２５を最終的に得るために使用されても良い。
【００３５】
本発明のいくつかの実施形態は、パラメトリック空間オーディオ符号化及び方向性オーディオ符号化（ＤｉｒＡＣ）のいくつかの処理ユニットに関する。ＤｉｒＡＣ符号化の仕組みは、マイクロホン入力に基づいてさらなる考慮のための例示的な土台として機能する、効率的な空間符号化技法である。不正確な空間パラメータの推定値に対処する本発明の手法について、以下ではさらに詳細に説明する。
【００３６】
本発明の内容はパラメトリック空間オーディオ符号化に関するものである。まずパラメトリック空間オーディオ表現のいくつかの重要な例を説明する。マイクロホン信号に基づく空間オーディオコーデックの例として、非特許文献１に記載のＤｉｒＡＣという手法を主として検討する。しかしながら、本発明の基本的な原理は、他の空間オーディオ符号化の仕組みにも適用可能であり、そこでは、あるスペクトル領域における空間パラメータ推定値における曖昧性を例えば異なる領域からのパラメータ推定値によって解消でき、あるいはビデオ画像などの他の追加の情報源を考慮することによって解消できる。ＤｉｒＡＣの具体的詳細を検討する前に、すべてのパラメトリック空間符号化の仕組みに共通の基本的な符号器／復号器構造について検討する。
【００３７】
本発明の一実施形態によるパラメトリック空間オーディオ符号器１０の概略構造が図１０に示されている。図１０はパラメトリック空間オーディオ符号器を示す。符号器１０は、非特許文献１及び非特許文献２に記載のように、例えばマイクロホン信号などの複数のオーディオチャネル２を入力信号として受け取る。符号器１０の出力はダウンミックス信号３及び対応するサイド情報４である。ダウンミックス信号３は、一つ又は複数のオーディオチャネルからなる。サイド情報４は、観測された音場を表わすパラメトリック・メタデータを含んでいる。符号器１０の出力、即ちダウンミックス信号３とサイド情報４との組み合わせが、以下では空間オーディオストリーム又は空間オーディオ表現と呼ばれる。
【００３８】
対応するパラメトリック空間オーディオ復号器２０の概略構造が図１１に示されている。図１１はパラメトリック空間オーディオ復号器を示す。復号器２０は空間オーディオストリームを入力として受け取る。ダウンミックス信号７及びサイド情報６に含まれるメタデータに基づき、復号器２０は、所望の再構成設定に対応するスピーカ信号又はヘッドホン信号を計算する。例示的なスピーカの配置は、例えば非特許文献３に記載されている。
【００３９】
方向性オーディオ符号化（ＤｉｒＡＣ）は、周波数サブバンドにおけるサウンドの瞬時の到来方向（ＤＯＡ）と拡散性とを用いた音場のパラメータ表現を使用するものであり、人間の聴覚に関係がある特徴だけを考慮するものである。ＤｉｒＡＣの手法は、もし音場のＤＯＡが正しく再構成されるならば、両耳間の時間差（ＩＴＤ）及び両耳間のレベル差（ＩＬＤ）が正しく知覚されるという仮定に基づいている。したがって、もし音場の拡散性が正しく再構成されるならば、両耳間のコヒーレンス（ＩＣ）が正しく知覚されるものと仮定される。このようにして、再生側が方向及び拡散性のパラメータとモノラルのマイクロホン信号とだけしか必要とせずに、任意のスピーカセットを用いて、与えられた聴取位置における人間の空間オーディオの知覚に関係する特徴を生成することができる。
【００４０】
以下の全体を通して、点状の音源が広い周波数範囲においてサウンドのあるＤＯＡを生じさせると想定することに留意すべきである。例を挙げると、一名の活発な話し手がスピーチの全スペクトルにわたってサウンドの所定のＤＯＡを生じさせると期待される。２名の話し手が同時に活発である場合には、各々の周波数帯域のＤＯＡは、その特定の周波数帯域において「最高のサウンドエネルギーを放射」する、最も支配的な音源に対応するＤＯＡによって決定されると期待される。
【００４１】
以下では、ＤｉｒＡＣにおけるパラメータ推定のためのＤｉｒＡＣ分析を説明する。ＤｉｒＡＣにおいては、所望のパラメータ、即ち各周波数帯域ｍ及び時間ブロックインデックスｋにおけるサウンドのＤＯＡ、即ちｅDOA(k,m)（又は対応する角度φ(k,m)によって表わされる）と拡散性Ψ(k,m)とを、非特許文献１に記載のようにＢフォーマットのマイクロホン信号に基づいて、音場のエネルギー分析によって推定することができる。Ｂフォーマットのマイクロホン信号は、無指向性信号Ｗ(k,m)と、直交座標系のｘ及びｙ方向に対応する二つのダイポール信号（Ｕx(k,m)及びＵy(k,m)）とからなる。Ｂフォーマットの信号を、例えば非特許文献４に記載のような特定のＢフォーマットマイクロホンを使用するか、又は非特許文献５に記載のような無指向性のマイクロホンのアレイを使用して、直接的に測定することができる。後者においては、音響信号の波長が無指向性のマイクロホンの間隔よりもはるかに大きいという仮定が根底に存在する。この仮定が当てはまらない高い周波数においては、ＤＯＡの推定に関して系統的なバイアス(systematic bias)が存在する。非常に低い周波数においては、マイクロホン自身のノイズの影響が大きくなる。このように、高い周波数における問題は、非常に狭い間隔のアレイを設計することによって単純に解決できるわけではない。ＤＯＡ推定についての詳細及び高い周波数における上述の問題は、以下の段落において説明する。
【００４２】
ＤｉｒＡＣなどの空間オーディオ符号化における一つの根本的な基礎は、ＤＯＡが正しく推定されることである。この役割を果たすために、時間／周波数ドメイン表現における活性強度ベクトルの方向が調査される。
【数１】
0005449624.tif 000041
【００４３】
上述の式で、上付き文字Hは共役転置を表わしている。ρ0は空気の平均密度を表わし、ｃは音速である。最後に、Ｉa(k, m)の対向方向がサウンドのＤＯＡを向いている。
【数２】
0005449624.tif 000042
ここで、φ(k, m)はＤＯＡの方位角を表わしている。上付き文字Tは転置を表わしている。
【００４４】
二つのダイポール信号Ｕx(k, m)及びＵy(k, m)は、８の字形の指向性を有するマイクロホンによってもたらすことができる。コストを節約する代案は、非特許文献６に記載のような無指向性マイクロホンの一次差動アレイを使用することであると考えられる。しかしながら、どちらの選択肢においても、通常は非常に低い周波数及び高い周波数において所望の８の字の指向性から逸脱してしまう。高品質の指向性マイクロホンを使用すれば、この逸脱を一次差動アレイと比べてはるかに小さくすることができる。指向性マイクロホンにおける逸脱が測定値によって表示される一方で、アレイに基づくＤＯＡ推定のバイアスは、分析的に定式化することが可能である。具体的なアレイ形状についての研究は、非特許文献７に表されている。一例として、図１３では、ｄ＝４ｃｍの間隔を有する５個のマイクロホンからなるアレイ（図１２参照）について、実際のＤＯＡの関数としての不正確なＤＯＡ推定値１０５を示す。この例において、バイアスを以下のように表現することができる。
【数３】
0005449624.tif 000043
【００４５】
低い周波数ｆにおいては、
0005449624.tif 000044
という近似を適用でき、ほとんどバイアスされていない推定値
0005449624.tif 000045
が得られる。例えば５ｋＨｚのような所定の限界に近づくと、式（３）を使用して表現できるバイアスされていない推定値からの明確な逸脱が観測される。図１３から分かるように、いわゆるエイリアシング周波数ｆaliasまではマッピング曲線２０１を使用することによって、バイアスを除去することができる。この例では、
0005449624.tif 000046
である。この限界を超える領域、即ちｆ＞ｆaliasにおいては、推定値
0005449624.tif 000047
をもはや単一のバイアスされていないＤＯＡ推定値φに割り当てることができない。この例では３つの候補１１５、即ち
0005449624.tif 000048
が存在し、そのうちのただ一つが正確なＤＯＡである。非特許文献７によれば、曖昧でないＤＯＡ推定の上限周波数は下記の式によって示される。
【数４】
0005449624.tif 000049
【００４６】
これをＤｉｒＡＣの空間エイリアシング周波数と呼ぶ。falias未満での系統的なバイアスの補償は、非特許文献７に記載されている。しかしながら、本発明は、faliasを超える周波数領域での曖昧性を解消するための手法を含む。曖昧性の解消が系統的なバイアスの補正をも実行することに留意すべきである。
【００４７】
具体的には、図１２は第５の無指向性マイクロホンを中央に有する無指向性マイクロホンの正方形アレイを示す。隣り合うマイクロホンの間隔はｄ／２である。さらに、図１３が４〜７ｋＨｚの間の周波数における方位角の推定値
0005449624.tif 000050
を示す。灰色の線９がバイアスされていない角度推定を表わしている。７ｋＨｚにおいて、空間エイリアシングに起因する曖昧性を明らかに見て取ることができる。即ち、一つの例示的な推定値
0005449624.tif 000051
を単一のバイアスされていない推定値に割り当てることが不可能である。３つの候補、即ち
0005449624.tif 000052
が存在し、その中からバイアスされていない推定値を選択しなければならない。
【００４８】
次に、ＤｉｒＡＣにおける方向性フィルタリング及びレンダリングのための処理原理としてのスペクトル重み付けを説明する。ＤｉｒＡＣ分析の後、即ちすべてのＤｉｒＡＣパラメータ（ダウンミックス信号Ｗ(k, m)を伴う、φ(k, m)により表わされるサウンドのＤＯＡ及び拡散性Ψ(k, m)）が計算された後で、これらのパラメータを任意のスピーカ配置についてスピーカ信号をレンダリングするめに使用することができる。任意ではあるが、ＤｉｒＡＣのパラメータドメインにおいて方向性フィルタリングを実行することができる。レンダリング及び方向性フィルタリングの両者は、スペクトル重み付けの原理に基づく。それらは、恐らくは曖昧なＤＯＡ推定値によって影響される。各々の処理ユニットの基本的な機能、とりわけスペクトル重み付け関数の計算について、以下の二つの段落で概略的に説明する。
【００４９】
再生側におけるＤｉｒＡＣ合成について説明すると、各々のスピーカ信号がモノラルのダウンミックス信号、ＤｉｒＡＣのＤＯＡ及び拡散性パラメータに基づいて計算される。スピーカ信号は、直接サウンド、即ちコヒーレントなサウンドと、拡散サウンドとにそれぞれ対応する信号成分からなる。拡散サウンド成分はＤＯＡに依存しない。したがって、ここでは取り扱わない。ｐ番目のスピーカの信号Ｓp(k,m)のコヒーレント部分は、下記の式に従って計算される。
【数５】
0005449624.tif 000053
【００５０】
パンニング利得ｇp(k,m)又は利得パラメータ（８１５）は、サウンドのＤＯＡφ(k,m)又は
0005449624.tif 000054
と、スピーカの数と、所望の聴取位置に対するｐ番目のスピーカの位置とに依存する。図１４は、３つのスピーカを−６０°（ｇ1(k,m)）、０°（ｇ2(k,m)）及び６０°（ｇ3(k,m)）に有する設定について、ＤＯＡの関数としてのいくつかの例示的な利得関数を示す。換言すると、図１４は、−６０°(ｇ1(k,m))、０°(ｇ2(k,m))及び６０°(ｇ3(k,m))における３個のスピーカを有する設定について、ＤＯＡφ(k,m)の関数としてのパンニング利得の例を示す。各々の利得関数がスペクトル利得へのＤＯＡの割り当てを決定する。間違ったＤＯＡの使用は空間像の歪みにつながりかねない。
【００５１】
方向性フィルタリングについて説明すれば、ＤｉｒＡＣの分析と合成との間の処理の連鎖において、例えば空間フィルタリングを実行するための追加の処理ステップを適用するという選択肢が存在する。目標は、特定の方向から到来するサウンドを保ちつつ、他の方向からのサウンドを減衰させることである。原理的には、この技法をＤｉｒＡＣパラメータドメインにおける一種の「ビーム形成」と考えることができる。非特許文献８において、方向性フィルタリングと呼ばれるそのような技法が提案されており、基本的には、方向性フィルタリングは時間／周波数表現のダウンミックス信号Ｗ(k, m)に利得関数を適用して出力信号スペクトルを得る。
【数６】
0005449624.tif 000055
【００５２】
上式内の下付き文字dfは「方向性フィルタリング」が加えられたことを示す。利得関数Ｄ(k,m)又は利得パラメータ（８１５）は、ＤＯＡ又は
0005449624.tif 000056
ならびに拡散性の関数として設計されている。設計の原理を図１５によって説明することができる。
【００５３】
具体的には、図１５は方向性フィルタリングの利得関数の指向性パターンを示す。瞬時ＤＯＡφ(k,m)（極座標における方位角）を対応する利得（極座標における半径）へと割り当てることによって、利得関数の設計を行うことができる。拡散性が増大するにつれて指向性が拡大する。所望の音源が、利得関数が１に等しい３０°に位置すると仮定する。
【００５４】
極座標は、各時間ブロックｋ及び各周波数ビンｍにおけるＤＯＡφ(k,m)の関数としての、時間及び周波数依存性の利得Ｄ(k,m)についての設計関数を表わしている。拡散性Ψ(k,m)が増大するにつれて設計関数が拡大する。非特許文献８に係る設計関数それ自身は、周波数に対して不変であることに留意すべきである。
【００５５】
次に音源定位の技法を説明する。方向性フィルタリング及び本発明の実施形態のために、所望の音源の数及び位置を知ることが不可欠である。コンパクトなマイクロホンの配置を有する二次元の観点においては、位置情報を方位角によって表わすことができる。非特許文献９は、ＤｉｒＡＣパラメータによって音源の数及び位置を特定するためのロバストな方法を提案している。この原理は、各々の瞬時ＤＯＡ、即ち時間及び周波数依存性のＤＯＡが、実際の音源の可能性のある方向と推定値の所定のばらつきに起因する不確実性とを加えたものを表わす、という仮定に基づいている。ＤＯＡは、ランダム変数として取り扱われる。平均化及びクラスタリングのそれぞれの操作によってばらつきを減らし、活性音源の位置についての情報を得ることができる。
【００５６】
この方法は、ダウンミックス信号のパワースペクトル密度（ＰＳＤ）によって重み付けされ、低い拡散性を有する時間／周波数タイルにおいてのみ計算される、ＤＯＡのヒストグラムに基づく。実際に、低い拡散性パラメータは、ＤＯＡの信頼性の値と解釈することができる。この技法は、時間／周波数表現においてスピーチ信号が疎であることを頼りにする。さらにThiergartの手法は、定位のために使用されるパラメータについての統計的な情報ももたらすことができる。例えば、ＤＯＡのばらつきは上述した図５の実施形態の考慮対象であってもよい。
【００５７】
本発明の実施形態においては、種々の音源定位技法がＤｉｒＡＣ及び方向性フィルタリングに関連して同様に機能できることに留意すべきである。画像又は映像入力データに基づく顔又は体の定位など、音響のキューとは異なるキューも考慮することができる様々な音源定位方法が存在している。
【００５８】
本発明の実施形態は、空間オーディオ符号化技法の空間パラメータの推定において生じる非一意性を解消するための効率的な方法を表わす。ＤｉｒＡＣのＤＯＡ推定器などのいくつかのパラメータ推定器において、曖昧性はバイアスに関係している。このような場合、曖昧性の解消は先の実施形態において述べたようなバイアスの補償を暗黙的に含む。
【００５９】
さらに、本発明の手法の概念的構造を図１６に示す。例えば方位角の推定値１０５、即ち
0005449624.tif 000057
によって示されるＤＯＡパラメータの推定値が曖昧性に関して分析される。これは、アレイの形状についての測定値又は情報に基づいて実行することができる。後者の例が式（３）による関係によって与えられる。結果として、所望のＤＯＡについての多数の候補
0005449624.tif 000058
が得られる。通常は、現在のＤＯＡ推定値が次の処理ステップのための入力パラメータとして機能する。関数ｆ（・）がそのようなユニットを一般的な様相で表わしている。未だ適切なＤＯＡ候補についての判断が行なわれていないため、並列なＮ個の処理ユニットを考えることができる。最後に、出力パラメータ
0005449624.tif 000059
を適切な手法で結合することができる。即ち、出力パラメータのうちの一つだけを選択しても良く、出力パラメータの適切な結合を計算してもよい。可能であれば、曖昧性の解消を改善するために、この「選択又は結合」のステップに先験的な情報をさらに取り入れることができる。この情報は、例えば音源の位置についての可能性のある情報であってもよい。
【００６０】
具体的には、図１６は空間パラメータ推定値
0005449624.tif 000060
の曖昧性を解消する本発明の全体的な構造を示す。このような解消は、最も適切な候補を選択するか、又は連続的なパラメータ
0005449624.tif 000061
を結合することによって実行することができる。利用可能であれば、先験的な情報がこの選択を支援すべく使用されても良い。
【００６１】
図１６の概念的構造のより具体的な実現を考えた場合、本発明の実施形態に係る二つの方法を導き出すことができる。一実施形態においては、
0005449624.tif 000062
と設定することができる。即ち、いかなる後続の処理にかかわらず、ＤＯＡパラメータを補正することを目標とする。したがって、次の段階において適切なＤＯＡを選択することができる。ここで、例えば平均値をとるなどの多数のＤＯＡを結合することは、所望の結果につながる可能性が低い。適切なＤＯＡの選択は、音源の位置についての先験的な情報に基づく。実際には、多数のＤＯＡから活性音源の位置までの距離が最短であるＤＯＡを選択しても良い。この手順のさらに具体的な説明は、図２〜図４の実施形態で提示した。
【００６２】
第２の方法は、ＤＯＡが通常は分析されて次の処理ユニットのための入力パラメータとして役立てられるという事実に基づく。既に述べたように、そのようなユニットは出力としてスペクトル利得パラメータを生成することが多い。曖昧なＤＯＡｓによって、多数の出力パラメータが得られる。これらのパラメータは、図１６に
0005449624.tif 000063
によって示されている。ＤＯＡそのものとは対照的に、アプリケーションに応じて例えば最大値、最小値、平均値をとる等の手法で、出力パラメータを結合できることが分ってきた。この方法の利点は、先験的な情報を入手する必要がない点にある。この方法は、図７〜図９において詳細に説明した別の実施形態を構成している。
【００６３】
図１７に、ＤＯＡ推定値の補正を目的とするすべての手順の機能の概要を示す。ＤＯＡ推定値１０５、即ち
0005449624.tif 000064
が対応する推定器１１から得られる。曖昧性の考慮の後で、空間エイリアシング周波数を超える領域で、実際の補正済みＤＯＡの複数の候補
0005449624.tif 000065
が得られる。「適切な」候補の選択は、活性音源の物理的なＤＯＡへの近さに基づく。ここで、最終的な補正済みＤＯＡ１２５、即ち
0005449624.tif 000066
を決定するために二つの手法が存在する。候補
0005449624.tif 000067
のうちで、活性音源−ＤＯＡのうちの一つへの距離が最小である候補を選択することができ（図３参照）、又は「最も近い候補」を活性音源のＤＯＡへと直接マッピングすることができる。この選択肢は、図４の実施形態において説明済みである。
【００６４】
具体的には、図１７は、補正済み推定値の曖昧性によって反転(inversion)が複雑となる場合に、バイアスされたＤＯＡ推定値を補正する構成を示す。複数の補正済み推定値からの適切な候補の選択は、活性音源の定位に基づいている。補正済みＤＯＡ
0005449624.tif 000068
に基づく事後処理は任意に実行される。
【００６５】
図３の実施形態は、複数のＤＯＡから一つの候補を選択する方法を説明する。ここでは、活性音源の現在のＤＯＡについて考える。これらのＤＯＡは広帯域のパラメータであり、典型的には瞬時ＤＯＡと比べて急激に変化することがない。faliasを超える時間／周波数タイルにおけるＤＯＡ推定値の非一意性を解消するために、個々の時間／周波数タイルにおける観測／推定された瞬時ＤＯＡはある音源によって引き起されると仮定し、その音源の位置は先験的な知識３１１又は定位技法によって既知であると仮定する。
【００６６】
図１８に空間エイリアシング周波数を超える周波数における曖昧なバイアスされたパラメータ推定値を補正するための例示的な手順を示す。図１８では、例えば
0005449624.tif 000069
である瞬時ＤＯＡが分析される。曖昧性を考慮し、実際の瞬時ＤＯＡについて３つの候補２１５、即ち
0005449624.tif 000070
が得られる。これと並行して、音源定位技法により、Ｑ1＝７５°に位置する一つの活性音源が分析されている（定位推定値３１１）。補正済みＤＯＡの空間における３つの候補から、
0005449624.tif 000071
が［°］におけるＱ1までのユークリッド距離が最も小さい。したがって、
0005449624.tif 000072
という瞬時ＤＯＡを
0005449624.tif 000073
へと補正することができる。
0005449624.tif 000074
を依然としてランダム変数と見なせることができることに留意すべきである。
【００６７】
具体的には、図１８は、空間エイリアシング周波数を超える周波数におけるバイアスされていないＤＯＡφの関数としての推定されたＤＯＡ
0005449624.tif 000075
を示す。瞬時推定値
0005449624.tif 000076
が３つの補正済みの瞬時ＤＯＡ
0005449624.tif 000077
に割り当てられる。音源定位技法により、活性音源Ｑ1が発見されている。これにより、３つの候補
0005449624.tif 000078
のうちＱ1までの距離が最も小さい一つが、補正済みの瞬時ＤＯＡ
0005449624.tif 000079
へと割り当てられる。即ち、
0005449624.tif 000080
となる。
【００６８】
さらなる実施形態によれば、二つ以上の活性音源をもって瞬時補正を実行する可能性もある。二つ以上の音源が同時に活性である場合が存在しうる。上述した例示的な音源定位技法は、複数の音源の数及び位置を特定できる。上述したようなバイアス及び非一様性の解消の手順を、二つ以上の活性音源を有する動作に拡張できる。即ち、まず図３の実施形態において説明したとおりの手順が実行される。一つの推定されたＤＯＡ
0005449624.tif 000081
について、補正済みＤＯＡの複数の候補が得られる。次に、これらの候補について、一つの実際の音源までではなく、二つ以上の実際の音源までの距離を分析することができる。図１９では、二つの活性音源Ｑ1及びＱ2を有する例の略図を示す。空間エイリアシングを超える周波数における推定されたＤＯＡが活性音源のうちの一つに対応すると仮定していることから、補正済みＤＯＡを、候補のうちで補正済みＤＯＡの空間におけるいずれかの活性音源に対して最も近くに位置する候補へと割り当てることができる。
【００６９】
具体的には、図１９がバイアスされていないＤＯＡφの関数としての推定されたＤＯＡ
0005449624.tif 000082
を示している。音源定位技法が二つ以上の活性音源を分析する場合には、可能性のある補正済みの瞬時ＤＯＡのうちの一つを、活性音源のいずれか（ここではＱ1）に最も近く位置するバイアスされていない瞬時ＤＯＡへと割り当てる。
【００７０】
音源のＤＯＡの推定値をランダム変数として考察することにする。その場合、この手順は、各音源のＤＯＡ推定値が互いに独立し、同様かつ対称的に分布している場合に最適である。各音源のＤＯＡの分布がこれらの要件を満たさない場合、より適切な距離の値を使用しても良い。この場合について、以下で説明する。
【００７１】
活性音源のＤＯＡの一般的な分布を説明する。非特許文献９が提案する音源定位技法は、ノイズ及びバイアスによる影響を無視することができるすべての周波数帯域における、瞬時ＤＯＡの確率密度関数（ｐｄｆ）を利用している。上述した例は、対応する分布が互いに独立し、同様かつ対称的である場合、曖昧なＤＯＡの適切な割り当てをリターンする。各々の音源のｐｄｆについて実現可能な近似として、正規分布が使用されても良い。図２０は、二つの活性音源Ｑ1及びＱ2（３１１）を有する例示的な場面を概説している。それらのＤＯＡは、同様には分布していない。曖昧性の考慮の後で、補正済みの瞬時ＤＯＡのための３つの候補
0005449624.tif 000083
が得られる。さらなる実施形態によれば、いずれかの活性音源までのマハラノビス距離(Mahalanobis distance)が最小であるような
0005449624.tif 000084
を適切な候補として選択することができる。この場合、たとえＱ2までのユークリッド距離の方が近くても、Ｑ1が統計的な意味で最も近い音源となる。
【００７２】
具体的には、図２０において、それぞれｐｄｆを有し、Ｑ1＝６０°及びＱ2＝１００°に位置する二つの活性音源を示す。補正済みの瞬時ＤＯＡのための３つの候補
0005449624.tif 000085
が存在する。任意の活性音源（この場合にはＱ1）までのマハラノビス距離が最小であるため、
0005449624.tif 000086
が適切な候補として選択される。
【００７３】
マハラノビス距離はより一般的な距離の値の一例であり、正規分布において最適である。ＤＯＡが異なる様子で分布している場合には、他の適切な距離の値を使用しても良い。
【００７４】
図４の実施形態は、瞬時ＤＯＡを活性音源のＤＯＡへとマッピングする例を説明する。上述したように、瞬時ＤＯＡを、所定の平均の周囲に所定の分散を有する確率プロセスからのランダム変数としてモデル化することができる。瞬時ＤＯＡのバイアスを考慮することによって、関連するプロセスの平均及び分散が変更されるが、その確率的性質は維持される。曖昧性を解消する場合には、適切な候補を選択するための代替的な手法は、瞬時ＤＯＡ推定値を最も近い活性音源（例えばＱ1）のＤＯＡで置き換えることである。これは、瞬時ＤＯＡの確率的性質がこの種の量子化ステップによって有意に変更されることを意味する。次に続くアプリケーションが空間パラメータのより現実的な確率的性質の回復を必要とする場合には、活性音源のＤＯＡにランダム変数を加えることを提案する。例を挙げると、ある程度のランダム性がＤＯＡパラメータの推定値へと導入された場合に、ＤｉｒＡＣ合成の後で知覚される空間オーディオ品質が改善されることが分かる（図５参照）。このランダム変数η(k,m)がゼロ平均及び単位分散(unit variance)を有すると定義しても、一般性を失うことはない。瞬時ＤＯＡの分散を回復するために、ステアリング分散（steering variance）σ2stを次式で定義する。
【数７】
0005449624.tif 000087
【００７５】
ステアリング分散を、瞬時ＤＯＡ及びＤｉｒＡＣの拡散性パラメータΨの関数として定義することができる。両方のパラメータは、それ自身時間及び周波数依存性である。時間及び周波数のインデックスは読み易さのために省略されている。最終的に、
0005449624.tif 000088
の平均及び分散が一致しなければならない。ここで式(７)の
0005449624.tif 000089
、Ｑ1、及びσst；ηは、図５の実施形態に示した
0005449624.tif 000090
即ち５３５と、
0005449624.tif 000091
即ち５２５と、変数５４５とにそれぞれ対応しても良い。
【００７６】
ＤＯＡ推定値の分散が、例えば非特許文献９に記載のように所定の音源定位アルゴリズムによって分析され、その後、ディザリングの量の調節に使用できることに留意すべきである。
【００７７】
図７〜図９の実施形態は、後続の処理ユニットの調節に関する。空間オーディオ・コーデックのすべてのパラメータを、上述したように、例えばＤｉｒＡＣ合成又は方向性フィルタリングなどの連続する処理段階を操作するために使用することができる。瞬時ＤＯＡのようなパラメータが何らかの方法で損なわれる場合、連続する処理段階に先立ってこれらのパラメータを補正することが、最終的な高品質の出力を達成するための一つのやり方である。空間パラメータを補正する方法は、図２〜図４の実施形態において説明済みである。代案として、連続する処理段階を予想されるバイアス及び曖昧性へと調節することができる。図１６の一般的な構成に基づいた関係する手順を振り返る。この場合、関数ｆ（・）が例えば方向性フィルタリングの操作を表わしている。
【００７８】
基本的に方向性フィルタリングは、スペクトル利得を信号スペクトルで乗算することによって働く（式（６）参照）。図２１は、ゼロ位相の伝達関数の設計を示す。瞬時ＤＯＡが０と１との間の利得へとマッピングされる。この手順を、各々の時間ブロックｋ及び周波数帯域ｍについて実行することができる。
【００７９】
具体的には、図２１がＤＯＡの関数としての方向性フィルタリングのための例示的な設計関数を示す。この関数は、各々の時間ブロックｋ及び周波数帯域ｍにおいて使用される。この設計関数は、Ｑ1＝２８０°において所望の活性音源を有するという仮定のもとで「操縦」される。
【００８０】
設計規則が補正済みの瞬時ＤＯＡ
0005449624.tif 000092
に依存する点に留意することが重要である。しかしながら、曖昧性の考慮なしで一次的な推定値
0005449624.tif 000093
に基づいて設計規則を機能させることができる。
【００８１】
図１６に示すとおり、ＤＯＡパラメータ推定値１０５、即ち
0005449624.tif 000094
が曖昧性に関して分析される。その結果として、所望のＤＯＡについての多数の候補
0005449624.tif 000095
が得られる。各々の候補
0005449624.tif 000096
について、対応する利得関数又は利得パラメータ８１５、即ち
0005449624.tif 000097
が計算される。例えば実際の方向性フィルタリング操作に使用される所望の利得関数
0005449624.tif 000098
を、異なる候補関数
0005449624.tif 000099
によって提供される値を適切に結合することで取得することができる。これが「選択又は結合」というブロック（セレクタ又は結合部として構成された曖昧性解消部１２０）によって表わされている。適切な結合規則として、線形結合（例えば平均）や種々の利得値の選択（例えば最大値又は最小値）に対応する非線形結合などが挙げられる。
【００８２】
以上の検討から、方向性フィルタリングのアプリケーション又は空間オーディオ合成の段階において、現在の音源の位置についての先験的な情報は不要である。この場合、フィルタリングの操作及び曖昧性の検討を効率的な手法で融合できることが分かる。この手法が図２２に示され、以下に要約される。
【００８３】
上述のように、関数ｆ（・）が正確なＤＯＡφ(k,m)の関数としての方向性フィルタリング利得関数の設計を表わす。
【数８】
0005449624.tif 000100
【００８４】
φ(k,m)を対応する推定値
0005449624.tif 000101
によって置換すべきであると分っているため、空間エイリアシング周波数を超える周波数においては、方向性フィルタリング利得関数の不正確な挙動が予想される。ＤＯＡ推定におけるこのような曖昧性を補償するために、強化された設計関数９１０、即ち
0005449624.tif 000102
を導入することができ、それらの設計関数はＤＯＡ推定値の曖昧性及び潜在的なバイアスを反映するよう選択されたものである。図１６及び図２２に関し、
【数９】
0005449624.tif 000103
となり、ｉ＝Ｉ，・・・，Ｎである。
【００８５】
いくつかの実施形態において、所望の強化された利得関数
0005449624.tif 000104
を、種々の候補利得値
0005449624.tif 000105
を適切に結合することによって得ることができる。
【００８６】
具体的には、図２２が曖昧なＤＯＡの関数として動作する処理ユニットを調節して補正済み出力パラメータ
0005449624.tif 000106
を生成する構成を示す。曖昧性は、所望のアプリケーションに関して解消されるのであり、実際のＤＯＡに関して解消されるのではない。図１７に示したシステムとは対照的に、音源定位ユニットが不要である。
【００８７】
以下に、ＤＯＡ推定器１１に続く処理段階における曖昧性の解消がどのように音源定位ユニットを省略可能にするのかを説明するための例を示す。図２１について検討する。ここでは、活性音源がＱ1に位置すると仮定し、そのＱ1が１という利得に割り当てられている。約０．８へと割り当てられた
0005449624.tif 000107
がその最も近くに位置するため、図３の実施形態において説明した通りに処理する場合には、選択された候補になると考えられる。ここで可能性のある出力利得を結合することができる。いくつかの実施形態において、もし方向性フィルタリングの主要な目標が所望の音源信号の品質を維持することである場合には、すべての可能性のある利得のうちで最大のものを常に選択しても良い。
【数１０】
0005449624.tif 000108
【００８８】
反対に、他の実施形態によれば、干渉をできるだけ抑えたい場合には次のように最小の利得を選択しても良い。
【数１１】
0005449624.tif 000109
【００８９】
実際のアプリケーションの筋書きに応じて、例えば平均をとるなど、Ｎ個の関数ｆenh,i（・）を結合するための種々の方法も選択肢となりうる。図２１の例を考えると、実際のＤＯＡ候補
0005449624.tif 000110
の結合が、方向性フィルタリングにおいてＤＯＡの曖昧性を解消するための適切な手法を表わしてはいないことを指摘しておかなければならない。
【００９０】
関数ｆenh,i（・）の結合をアレイの形状及び所望の指向性パターンに基づいて効率的に事前計算できることに留意すべきである（図１５参照）。したがって、アルゴリズムの実行時に追加計算の負担が生じることがない。なぜなら、この手法の出力は推定されたＤＯＡの関数から構成されており、対応する補正済みバージョンから構成されてはいないからである。
【００９１】
これまでの実施形態を要約すると、本発明が提示する技法／概念は、以下のステップを含むことができる。最初に、マイクロホンアレイからの信号が記録されても良い。次いで、到来方向（ＤＯＡ）がこれらのマイクロホン信号から推定されても良い。次いで、ＤＯＡの曖昧性が活性音源の位置についての情報に基づいて解消されても良い。この情報は事前に入手可能であり、及び／又は作動時間の中で推定されても良い。代替的に、曖昧性をＤＯＡに関して解消するのではなく、例えばスペクトル利得など後続の処理段階のパラメータに関して解消してもよい。最後に、例えば方向性フィルタリングなどの後続の処理段階を適用することができる。
【００９２】
幾つかの態様を装置の文脈の中で説明してきたが、これら態様は対応する方法をもまた表すものであり、ブロックや装置は方法ステップ又は方法ステップの特徴に相当していることは明らかである。同様に、方法ステップの文脈の中で表された態様は、また対応するブロック、項目、又は対応する装置の特徴をも表している。方法ステップの一部又は全てを、例えばマイクロプロセッサ、プログラム可能なコンピュータ、又は電子回路などのハードウェア装置によって実行しても良い。いくつかの実施形態においては、方法ステップの最も重要な一つ又は複数をそのような装置によって実行しても良い。
【００９３】
本発明のパラメータを、デジタル記憶媒体に記憶することができ、又はインターネットなどの無線伝送媒体若しくは有線伝送媒体を介して伝送することができる。
【００９４】
所定の構成要件にも依るが、本発明の実施形態はハードウエア又はソフトウエアにおいて構成可能である。この構成は、その中に格納される電子的に読み取り可能な制御信号を有し、本発明の各方法が実行されるようにプログラム可能なコンピュータシステムと協働する（又は協働可能な）、デジタル記憶媒体、例えばフレキシブルディスク，ＤＶＤ，ＣＤ，ＲＯＭ，ＰＲＯＭ，ＥＰＲＯＭ，ＥＥＰＲＯＭ，フラッシュメモリなどを使用して実行することができる。したがって、デジタル記憶媒体はコンピュータ読み取り可能であってよい。
【００９５】
本発明に従う幾つかの実施形態は、電子的に読み取り可能な制御信号を有するデータキャリアを含んでも良く、その制御信号は上述した方法の１つを実行するようプログラム可能なコンピュータシステムと協働可能である。
【００９６】
一般的に、本発明の実施例はプログラムコードを有するコンピュータプログラム製品として構成することができ、このプログラムコードは当該コンピュータプログラム製品がコンピュータ上で作動するときに、本発明の方法の一つを実行するよう作動する。そのプログラムコードは例えば機械読み取り可能なキャリアに記憶されていても良い。
【００９７】
本発明の他の実施形態は、上述した方法の１つを実行するための機械読み取り可能なキャリア上に記憶されたコンピュータプログラムを含む。
【００９８】
換言すれば、本発明方法のある実施形態は、そのコンピュータプログラムがコンピュータ上で作動するときに上述した方法の１つを実行するためのプログラムコードを有する、コンピュータプログラムである。
【００９９】
本発明の他の実施形態は、上述した方法の１つを実行するために記録されたプログラムコードを含むデータキャリア（又はデジタル記憶媒体、又はコンピュータ読み取り可能な媒体）である。そのデータキャリア、デジタル記憶媒体、又は記録された媒体は、典型的には有形及び／又は非遷移型 (non-transitory)である。
【０１００】
本発明の他の実施形態は、上述した方法の１つを実行するためのコンピュータプログラムを表現する、データストリーム又は信号のシーケンスである。そのデータストリーム又は信号のシーケンスは、例えばインターネットを介したデータ通信接続を介して伝送されるように構成されても良い。
【０１０１】
他の実施形態は、上述した方法の１つを実行するように構成又は適用された、例えばコンピュータ又はプログラム可能な論理デバイスのような処理手段を含む。
【０１０２】
他の実施形態は、上述した方法の１つを実行するためのコンピュータプログラムがインストールされたコンピュータを含む。
【０１０３】
本発明の他の実施形態は、上述した方法の１つを実行するためのコンピュータプログラムを受信機へと（例えば、電子的又は光学的に）送信するように構成された、装置又はシステムを含む。受信機は、例えばコンピュータ、モバイルデバイス、メモリデバイス、などであってよい。その装置又はシステムは、例えばコンピュータプログラムを受信機へと送信するためのファイルサーバを備えても良い。
【０１０４】
幾つかの実施形態においては、（例えば書換え可能ゲートアレイのような）プログラム可能な論理デバイスが、上述した方法の幾つか又は全ての機能を実行するために使用されても良い。幾つかの実施形態では、書換え可能ゲートアレイは、上述した方法の１つを実行するためにマイクロプロセッサと協働しても良い。一般的に、そのような方法は、好適には任意のハードウエア装置によって実行される。
【０１０５】
上述した実施の形態は、本発明の原理を単に例示的に示したにすぎない。本明細書に記載した構成及び詳細について、修正及び変更が可能であることは、当業者にとって明らかである。従って、本発明は、本明細書に実施形態の説明及び解説の目的で提示した具体的詳細によって限定されるものではなく、添付した特許請求の範囲によってのみ限定されるべきである。
【０１０６】
要約すると、本発明は、空間エイリアシングによって引き起される曖昧性を解消する方法を含む。このような曖昧性を空間パラメータに関して解消することができ、又は空間パラメータを入力として使用する後続の処理ユニットにおいて解消することができる。
【符号の説明】
【０１０７】
１００，２００，３００，４００，６００，７００，８００，９００  曖昧性を解消する装置
１０１  バイアス情報
１０５  ＤＯＡ推定値
１１０  ＤＯＡ推定値分析部
１１５  曖昧な分析パラメータ
１２０  曖昧性解消部
１２５  曖昧でない解消済みパラメータ
【図１】
0005449624.tif 000111
【図２】
0005449624.tif 000112
【図３】
0005449624.tif 000113
【図４】
0005449624.tif 000114
【図５】
0005449624.tif 000115
【図６】
0005449624.tif 000116
【図７】
0005449624.tif 000117
【図８】
0005449624.tif 000118
【図９】
0005449624.tif 000119
【図１０】
0005449624.tif 000120
【図１１】
0005449624.tif 000121
【図１２】
0005449624.tif 000122
【図１３】
0005449624.tif 000123
【図１４】
0005449624.tif 000124
【図１５】
0005449624.tif 000125
【図１６】
0005449624.tif 000126
【図１７】
0005449624.tif 000127
【図１８】
0005449624.tif 000128
【図１９】
0005449624.tif 000129
【図２０】
0005449624.tif 000130
【図２１】
0005449624.tif 000131
【図２２】
0005449624.tif 000132