(57)【要約】本発明は、オーディオデータの処理に関する。本発明は、（ａ）三次元空間で伝わり基準点から第１距離（Ｐ）に配置された音源から生じる音を表す信号を符号化して、前記基準点に対応する原点の球面調和関数基数で表される成分による音の表現を得ること、（ｂ）音響再生の場合、再生点（ＨＰｉ）と聴き手が通常位置する聴覚点（Ｐ）との間の距離を定義する第２距離（Ｒ）に基づき、フィルタ処理による近接音場効果の補償を前記成分に適用すること、を特徴とする。
2006506918.tif 000002
【特許請求の範囲】
【請求項１】
音響データを処理する方法であって、
ａ）三次元空間で伝播し基準点（Ｏ）から第１距離（ρ）に位置する音源から生じる少なくとも１つの音を表す信号が、前記基準点（Ｏ）に対応する原点の球面調和関数の基数で表される成分（Ｂｍｎσ）による音の表現を得るように符号化され、
ｂ）近接音場効果の補償が、再生装置による音の再生の場合、再生点（ＨＰｉ）と聴覚点（Ｐ）との間の距離を実質的に定義する第２距離（Ｒ）に依存するフィルタ処理によって前記成分（Ｂｍｎσ）に適用される方法。
【請求項２】
請求項１に記載の方法であって、前記音源は、基準点（Ｏ）から遠方に設置され、
連続次数ｍの成分が、球面調和関数の前記基数での音の表現用に取得され、
フィルタ（１／Ｆｍ）が適用され、その係数は、各々次数ｍの成分に適用されるが、ベキｍの多項式の逆数形式で解析的に表され、その変数は、再生装置のレベルに近接音場効果を補償するように、音響周波数及び前記第２距離（Ｒ）に反比例する方法。
【請求項３】
請求項１に記載の方法であって、前記音源は、前記第１距離（ρ）に想定される仮想音源であり、
連続次数ｍの成分は、球面調和関数の前記基数で音を表現するために取得され、
広域フィルタ（Ｈｍ）が適用され、その係数は、各々次数ｍの成分に適用されるが、分数の形式で解析的に表され、
分子は、ベキｍの多項式であり、その変数は、仮想音源の近接音場効果をシミュレートするように、音響周波数及び前記第１距離（ρ）に反比例し、
分母は、ベキｍの多項式であり、その変数は、低音響周波数の仮想音源の近接音場効果を補償するように、音響周波数及び前記第２距離（Ｒ）に反比例する方法。
【請求項４】
上記請求項の１つに記載の方法であって、
ステップａ）及びｂ）において符号化及びフィルタ処理されたデータは、前記第２距離を表すパラメータ（Ｒ／ｃ）と共に再生装置に送信される方法。
【請求項５】
請求項１乃至３の１つに記載の方法であって、再生装置は、記憶媒体を読み取るための手段を含み、
ステップａ）及びｂ）において符号化及びフィルタ処理されたデータは、前記第２距離を表すパラメータ（Ｒ／ｃ）と共に、再生装置が読み取る記憶媒体に記憶される方法。
【請求項６】
請求項４及び５の１つに記載の方法であって、
前記聴覚点（Ｐ）から第３距離（Ｒ２）に配置された複数のスピーカを含む再生装置による音再生の前に、前記第２（Ｒ１）及び第３距離（Ｒ２）に依存する係数を有する適応フィルタ（Ｈｍ（Ｒ１／ｃ，Ｒ２／ｃ））が、符号化及びフィルタ処理されるデータに適用される方法。
【請求項７】
請求項６に記載の方法であって、
前記適応フィルタ（Ｈｍ（Ｒ１／ｃ，Ｒ２／ｃ））の係数は、各々次数ｍの成分に適用されるが、分数の形式で解析的に表され、
分子は、ベキｍの多項式であり、その変数は音響周波数及び前記第２距離（Ｒ）に反比例し、
分母は、ベキｍの多項式であり、その変数は音響周波数及び前記第３距離（Ｒ２）に反比例する方法。
【請求項８】
請求項２、３及び７の１つに記載の方法であって、ステップｂ）を実行するために、
偶数次数ｍの成分に関しては、次数２のカスケード状セル形式の音響デジタルフィルタと、
奇数次数ｍの成分に関しては、次数２のカスケード状セルと次数１の追加セルとの形式の音響デジタルフィルタと、が提供される方法。
【請求項９】
請求項８に記載の方法であって、次数ｍの成分の場合、音響デジタルフィルタの係数は、前記ベキｍの多項式のベキ根の数値により定義される方法。
【請求項１０】
請求項２、３、７、８及び９の１つに記載の方法であって、前記多項式は、ベッセル多項式である方法。
【請求項１１】
請求項１、２、及び４乃至１０の１つに記載の方法であって、
三次元空間で伝播する少なくとも１つの音を表す前記信号を得るように、実質的に前記基準点（Ｏ）に対応する中心を有する球の表面上に実質的に配置されるアレイ状の音響変換器を含むマイクが提供される方法。
【請求項１２】
請求項１１に記載の方法であって、
広域フィルタが、一方では、前記第２距離（Ｒ）の関数として近接音場効果を補償するように、他方では、変換器から生じる信号を等化して前記変換器の方向性の重み付けを補償するように、ステップｂ）において適用される方法。
【請求項１３】
請求項１１及び１２の１つに記載の方法であって、
球面調和関数の前記基数で音を表すために選択された成分の総数に依存する複数の変換器が提供される方法。
【請求項１４】
上記請求項の１つに記載の方法であって、
ステップａ）では、成分の総数は、再生時、音の再生が忠実であり成分の総数と共に増加する寸法を有する知覚点（Ｐ）周辺の空間領域を得るように、球面調和関数の基数から選択される方法。
【請求項１５】
請求項１４に記載の方法であって、
成分の前記総数に少なくとも等しい数のスピーカを含む再生装置が提供される方法。
【請求項１６】
請求項１乃至５及び８乃至１３の１つに記載の方法であって、
聴き手からの選択した距離に配置された少なくとも第１及び第２スピーカを含む再生装置が提供され、
聴き手からの所定の基準距離（Ｒ）に位置する音源の空間位置の認識キューが、この聴き手に対して取得され、
ステップｂ）の補償は、前記基準距離が実質的に第２距離として適用される方法。
【請求項１７】
請求項４及び５の１つと組み合わせて選択される請求項１乃至３及び８乃至１３の１つに記載の方法であって、
聴き手からの選択した距離に配置された少なくとも第１及び第２スピーカを含む再生装置が提供され、
聴き手からの所定の基準距離（Ｒ２）に位置する音源の空間位置の認識キューが、この聴き手に対して取得され、
再生装置による音再生の前に、第２距離（Ｒ）及び実質的に基準距離（Ｒ２）に依存する係数を有する適応フィルタ（Ｈｍ（Ｒ／ｃ，Ｒ２／ｃ））が、ステップａ）及びｂ）において符号化及びフィルタ処理されるデータに適用される方法。
【請求項１８】
請求項１６及び１７の１つに記載の方法であって、
再生装置は、聴き手のそれぞれの耳用に２つのヘッドホンを備えたヘッドセットを含み、
各ヘッドホンに対して別々に、ステップａ）及びｂ）の符号化及びフィルタ処理は、第１距離（ρ）として、再生される音源の位置（Ｍ）から各耳をそれぞれ離間する距離（ｒＲ，ｒＬ）で、各ヘッドホンに供給すべきそれぞれの信号に関して適用される方法。
【請求項１９】
上記請求項の１つに記載の方法であって、ステップａ）及びｂ）において行列系が形成され、前記系には、少なくとも、
球面調和関数の基数での前記成分を含む行列（Ｂ）と、
ステップｂ）のフィルタ処理係数に対応する係数を有する対角行列（対角（１／Ｆｍ）と、が含まれ、前記行列は乗算され補償成分
【数１】
2006506918.tif 000072
の結果行列が得られる方法。
【請求項２０】
請求項１９に記載の方法であって、
再生装置は、聴覚点（Ｐ）から同一の距離（Ｒ）に実質的に配置される複数のスピーカを含み、
ステップａ）及びｂ）で符号化及びフィルタ処理された前記データを復号化するために、また、前記スピーカへの供給に適する信号を形成するために、
行列系が、再生装置に特有な、前記結果行列
【数２】
2006506918.tif 000073
及び所定の復号行列（Ｄ）を含んで形成され、
行列（Ｓ）が、補償成分
【数３】
2006506918.tif 000074
と前記復号行列（Ｄ）の乗算によりスピーカ供給信号を表す係数を含んで得られる方法。
【請求項２１】
実質的に球の表面に配置されたアレイ状音響変換器が備えられたマイクを含む集音装置であって、前記装置には、更に、処理ユニットが含まれ、前記処理ユニットは、
変換器が発する各信号を受信し、
前記球（Ｏ）の中心に対応する原点の、球面調和関数の基数で表された成分（Ｂｍｎσ）による音の表現を得るように前記信号に符号化を適用し、
一方では、球の半径（ｒ）に対応する距離に依存し、他方では、基準距離（Ｒ）に依存するフィルタ処理を前記成分（Ｂｍｎσ）に適用するように、構成されることを特徴とする装置。
【請求項２２】
請求項２１に記載の装置であって、
前記フィルタ処理は、一方では、前記変換器の方向性の重み付けを補償するように、変換器から生じる信号を球の半径の関数として等化すること、また、他方では、音の再生の場合、再生点（ＨＰｉ）と聴覚点（Ｐ）との間の距離を実質的に定義する選択された基準距離（Ｒ）の関数として近接音場効果を補償すること、に存することを特徴とする装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、オーディオデータの処理に関する。
【背景技術】
【０００２】
特に、専用音響シミュレーション及び／又は再生等の三次元空間における音波の伝播に関係する手法は、音響及び心理的音響現象のシミュレーションに適用されるオーディオ信号処理法を実現する。このような処理方法は、音場の空間的符号化、その送信、及び一組のスピーカ又はステレオ式ヘッドセットのヘッドホンでのその空間化再生を行う。
【０００３】
空間化音の手法の中で特徴的な手法には、互いに補完的であるが一般的に同一系内で双方共実現される２つの処理カテゴリである。
【０００４】
一方では、第１処理カテゴリは、空間効果、即ち、より一般的には環境効果を合成するための方法に関する。１つ又は複数の音源（発信信号、位置、方位、方向性等）の記述から、また（空間幾何学的配置あるいは所望の音響知覚を含む）空間効果モデルに基づき、一組の基本的音響現象（直接波、反射波、又は回折波）あるいは巨視的音響現象（反響及び拡散音場）を計算し記述し、三次元空間において、選択した聴覚点に位置する聴き手のレベルで空間効果を伝え得る。次に、通常、反射（主受信波の再放射で活性化し空間位置属性を有する“二次”音源）に関係する一組の信号及び／又は遅延反響に関係する一組の信号（拡散音場に対する非相関信号）を計算する。
【０００５】
他方では、第２カテゴリの方法は、音源の位置的又は方向的表現に関する。これらの方法は、（一次及び二次音源を含む）上述した第１カテゴリの方法によって決定される信号に適用され、それらに関係する空間的記述（音源の位置）の関数として適用される。特に、この第２カテゴリによるこのような方法では、スピーカ又はヘッドホンに分配される信号を取得して、最終的に、聴き手周辺の所定位置それぞれに配置された音源の聴覚的印象を聴き手に与え得る。この第２カテゴリによる方法は、聴き手による音源位置認識が、三次元空間に分布するため、“三次元音像の生成法”と称される。一般的に、第２カテゴリによる方法には、基本的音響現象を空間的に符号化する第１ステップが含まれ、これによって三次元空間における音場の表現が生成される。第２ステップでは、後で用いるために、この表現が送信又は記憶される。復号化の第３ステップでは、復号信号が再生装置のスピーカ又はヘッドホンに伝えられる。
【０００６】
本発明は、どちらかと言えば、上記第２カテゴリに包含され、特に、音源の空間的符号化やこれら音源の三次元音響表現の仕様に関する。本発明は、また、１つ又は複数の三次元アレイのマイクによる集音時、自然音場の“音響”符号化について、“仮想的”音源の符号化（音源がシミュレートされるゲームや空間化会議等の用途）にも同様に該当する。
【０００７】
考え得る音響空間化手法の中では、“アンビソニック”法が好まれる。アンビソニック符号化は、詳細に後述するが、球面調和関数の基数で１つ又は複数の音波に関係する信号を表示する（特に、仰角及び方位角を含む球座標で音又は複数の音の方向を特徴付ける）ことにある。また、これらの信号を表しこの球面調和関数の基数で表される成分は、近接音場で放射される波に関して、この音場を発する音源と球面調和関数の基数の原点に対応する点との間の距離にも依存する。特に、この距離依存性は、後で分かるように、音響周波数の関数として表される。
【０００８】
このアンビソニック手法は、特に、仮想音源のシミュレーションの点では、極めて多くの機能を提供し、また、一般的には、次の利点を示す。即ち、
音響現象の実体を合理的に伝え、現実的で得心させる埋没的空間表現を提供する。
音響現象の表現は、適応性がある。即ち、様々な状況に適応し得る空間分解能を提供する。具体的には、この表現は、符号化信号送信時の処理能力制約条件及び／又は再生装置の制限条件の関数として送信及び利用し得る。
アンビソニック表現は、柔軟性があり、また、音場の回転のシミュレーションが可能であり、あるいは、再生時、多様な幾何学的配置のあらゆる再生装置にアンビソニック信号の復号化を適応し得る。
【０００９】
公知のアンビソニック法では、仮想音源の符号化は、本質的に方向性を有する。符号化機能は、球面座標の仰角及び方位角に依存する球面調和関数によって表される音波の入射角に依存する利得を計算することである。特に、復号化の際、スピーカは、再生時、遠方配置されると仮定する。これにより、再生波面の形状が歪む（即ち、湾曲する）。具体的には、上述したように、球面調和関数の基数での信号音の成分は、近接音場の場合、実際、音源の距離及び音響周波数にも依存する。より厳密には、これらの成分は、上記距離及び音響周波数に反比例する変数を有する多項式の形式で数学的に表し得る。従って、アンビソニック成分は、それらの理論的表現の点で、低周波で発散し、特に、それらが、有限の距離にある音源によって放射される近接音場音を表す場合、音響周波数が減少しゼロになると、無限大になる。この数学的現象は、アンビソニック表現の分野では、次数１の場合については、用語“低音域ブースト”によって既に知られている。特に、ゲルゾン（M.A.GERZON）による“聴覚局在化の一般メタ理論”、第９２回ＡＥＳ会議予稿３３０６，１９９２年、頁５２において公知である。この現象は、高ベキの多項式を含む高球面調和関数次数の場合、特に重要になる。
【００１０】
以下の文献、ゾンタッヒ（SONTACCHI）及びヘルドリッヒ（HOELDRICH）による“距離符号化を用いた３Ｄ音場の詳細調査”（２００１年１２月６日〜８日、アイルランド、リメリックにおけるデジタルオーディオ効果（DAFX-01）に関するCOST_G-6会議の議事録）は、アンビソニック表現に近い表現の範囲内で波面の湾曲を考慮するための手法を開示し、その原理は以下に存する。即ち、
ＷＦＳ（“波音場合成”の意）型の（シミュレートされる）仮想集音から生じる信号に（高次数の）アンビソニック符号化を適用すること。
１つの領域に渡って音場を再現することであって、領域境界に渡るその値、従って、ホイヘンスーフレネルの原理に基づき再現すること。
【００１１】
しかしながら、この文献に提示された手法は、高次数に対してアンビソニック表現を用いるという事実により期待が持てるが、いくつかの問題が生じる。即ち、
ホイヘンスーフレネルの原理の適用を可能にする全ての面の計算に必要なコンピュータ資源だけでなく、必要な計算時間が過剰である。
狭い間隔で配置された仮想マイク格子が選択されない場合、マイク間の距離のために“空間的エイリアシング”と称する処理アーティファクトが生じることによって、処理が更に面倒になる。
この手法は、集音時、実音源が存在すると、センサがアレイ状に配置される実例への移行が困難である。
再生時、三次元音響表現が、再生装置の固定半径に暗黙裡に束縛される。これは、ここでは、初期アレイ状マイクと同じ寸法のアレイ状スピーカ上でアンビソニック復号化を行わなければならないためであり、この文献は、他のサイズの再生装置に符号化又は復号化を適応する手段を提示してしない。
【００１２】
結局、この文献は、水平アレイ状センサを提示し、これによって、対象の音響現象が、ここでは、水平方向にのみ伝播すると仮定することによって、他の方向のあらゆる伝播を除外し、また従って、通常の音場の物理的実体を表さない。
【００１３】
更に一般的に、この手法は、任意の種類の音源、特に、近接音場源ではなく、むしろ遠方音源（平面波）の充分な処理が可能であり、このことは、数多くの用途において制限的及び人為的状況に対応する。
【発明の開示】
【発明が解決しようとする課題】
【００１４】
本発明の目的は、任意の種類の音場、特に近接音場における音源の効果を符号化、送信、及び再生によって処理するための方法を提供することである。
【００１５】
本発明の他の目的は、方向だけでなく距離に関しても、仮想音源の符号化を可能にする方法を提供し、任意の再生装置に適応可能な復号方式を定義することである。
【００１６】
本発明の他の目的は、特に、三次元アレイ状のマイクを用いて自然音場の集音を行うために、（低周波を含む）任意の音響周波数の音を処理する強力な方法を提供することである。
【課題を解決するための手段】
【００１７】
この目的のために、本発明は、音響データを処理する方法を提案する。本方法においては、
ａ）三次元空間で伝播し基準点から第１距離に位置する音源から生じる少なくとも１つの音を表す信号が、前記基準点に対応する原点の球面調和関数の基数で表される成分による音の表現を得るように符号化され、
ｂ）近接音場効果の補償が、再生装置による音の再生の場合、再生点と聴覚点との間の距離を実質的に定義する第２距離に依存するフィルタ処理によって前記成分に適用される。
【００１８】
第１実施形態では、前記音源は、基準点から遠方に設置され、
連続次数ｍの成分が、球面調和関数の前記基数での音の表現用に取得され、
フィルタが適用され、その係数は、各々次数ｍの成分に適用されるが、ベキｍの多項式の逆数形式で解析的に表され、その変数は、再生装置のレベルに近接音場効果を補償するように、音響周波数及び前記第２距離に反比例する。
【００１９】
第２実施形態では、前記音源は、前記第１距離に想定される仮想音源であり、
連続次数ｍの成分は、球面調和関数の前記基数で音を表現するために取得され、
広域フィルタが適用され、その係数は、各々次数ｍの成分に適用されるが、分数の形式で解析的に表され、
分子は、ベキｍの多項式であり、その変数は、仮想音源の近接音場効果をシミュレートするように、音響周波数及び前記第１距離に反比例し、
分母は、ベキｍの多項式であり、その変数は、低音響周波数の仮想音源の近接音場効果を補償するように、音響周波数及び前記第２距離に反比例する。
【００２０】
好適には、ステップａ）及びｂ）において符号化及びフィルタ処理されたデータを前記第２距離を表すパラメータと共に再生装置に送信する。
【００２１】
補足又は変形例として、再生装置は、記憶媒体を読み取るための手段を含み、再生装置が読み取る記憶媒体に、前記第２距離を表すパラメータと共に、ステップａ）及びｂ）において符号化及びフィルタ処理されたデータを記憶する。
【００２２】
利点として、前記聴覚点から第３距離に配置された複数のスピーカを含む再生装置による音再生の前に、前記第２及び第３距離に依存する係数を有する適応フィルタが、符号化及びフィルタ処理されるデータに適用される。
【００２３】
ある実施形態では、前記適応フィルタの係数は、各々次数ｍの成分に適用されるが、分数の形式で解析的に表され、
分子は、ベキｍの多項式であり、その変数は音響周波数及び前記第２距離に反比例し、
分母は、ベキｍの多項式であり、その変数は音響周波数及び前記第３距離に反比例する。
【００２４】
利点として、ステップｂ）を実行するために、
偶数次数ｍの成分に関しては、次数２のカスケード状セル形式の音響デジタルフィルタと、
奇数次数ｍの成分に関しては、次数２のカスケード状セルと次数１の追加セルとの形式の音響デジタルフィルタと、が提供される。
【００２５】
本実施形態では、次数ｍの成分の場合、音響デジタルフィルタの係数は、前記ベキｍの多項式のベキ根の数値により定義される。
【００２６】
ある実施形態では、前記多項式は、ベッセル多項式である。
【００２７】
信号音の集音時、三次元空間で伝播する少なくとも１つの音を表す前記信号を得るように、実質的に前記基準点に対応する中心を有する球の表面上に実質的に配置されるアレイ状の音響変換器を含むマイクが利点として提供される。
【００２８】
本実施形態では、広域フィルタが、一方では、前記第２距離の関数として近接音場効果を補償するように、他方では、変換器から生じる信号を等化して前記変換器の方向性の重み付けを補償するように、ステップｂ）において適用される。
【００２９】
好適には、球面調和関数の前記基数で音を表すために選択された成分の総数に依存する複数の変換器を提供する。
【００３０】
有利な特徴によれば、ステップａ）では、成分の総数は、再生時、音の再生が忠実であり成分の総数と共に増加する寸法を有する知覚点周辺の空間領域を得るように、球面調和関数の基数から選択される。
【００３１】
好適には、更に、成分の前記総数に少なくとも等しい数のスピーカを含む再生装置が提供される。
【００３２】
バイノーラル又はトランスオーラル合成による再生装置の枠内における変形例として、
聴き手からの選択した距離に配置された少なくとも第１及び第２スピーカを含む再生装置が提供され、
聴き手からの所定の基準距離に位置する音源の空間位置の予想認識キューが、この聴き手に対して取得され、いわゆる“トランスオーラル”又は“バイノーラル合成”手法を適用し、
ステップｂ）の補償は、前記基準距離が実質的に第２距離として、適用される。
【００３３】
２つのヘッドホンを有するように適合される再生装置の変形例では、
聴き手からの選択した距離に配置された少なくとも第１及び第２スピーカを含む再生装置が提供され、
聴き手からの所定の基準距離に位置する音源の空間位置の認識キューが、この聴き手に対して取得され、
再生装置による音再生の前に、第２距離及び実質的に基準距離に依存する係数を有する適応フィルタが、ステップａ）及びｂ）において符号化及びフィルタ処理されるデータに適用される。
【００３４】
特に、バイノーラル合成による再生装置の枠内では、
再生装置は、聴き手のそれぞれの耳用に２つのヘッドホンを備えたヘッドセットを含み、
好適には、各ヘッドホンに対して別々に、ステップａ）及びｂ）の符号化及びフィルタ処理は、第１距離として、再生空間で再生される音源の位置と各耳とをそれぞれ離間する距離で、各ヘッドホンに供給すべきそれぞれの信号に関して適用される。
【００３５】
好適には、ステップａ）及びｂ）において行列系が形成され、前記系には、少なくとも、
球面調和関数の基数での前記成分を含む行列と、
ステップｂ）のフィルタ処理係数に対応する係数を有する対角行列と、
が含まれ、前記行列は乗算され補償成分の結果行列が得られる。
【００３６】
好適には、再生時、
再生装置は、聴覚点から同一の距離に実質的に配置される複数のスピーカを含み、
ステップａ）及びｂ）で符号化及びフィルタ処理された前記データを復号化するために、また、前記スピーカへの供給に適する信号を形成するために、
行列系が、再生装置に特有な、前記補償成分の結果行列及び所定の復号行列を含んで形成され、
行列が、結果行列と前記復号行列の乗算によりスピーカ供給信号を表す係数を含んで得られる。
【００３７】
また本発明は、実質的に球の表面に配置されたアレイ状音響変換器が備えられたマイクを含む集音装置を目的とする。本発明によれば、本装置は、更に、処理ユニットが含まれ、この処理ユニットは、
変換器が発する各信号を受信し、
前記球の中心に対応する原点の、球面調和関数の基数で表された成分による音の表現を得るように前記信号に符号化を適用し、
また、一方では、球の半径に対応する距離に依存し、他方では、基準距離に依存するフィルタ処理を前記成分に適用するように、構成される。
【００３８】
好適には、処理装置によって行われるフィルタ処理は、一方では、前記変換器の方向性の重み付けを補償するように、変換器から生じる信号を球の半径の関数として等化すること、また、他方では、前記基準距離の関数として近接音場効果を補償することにある。
【発明を実施するための最良の形態】
【００３９】
本発明の他の利点及び特徴は、本明細書の以下の詳細な説明を解釈し添付図面を検討すると明らかになる。
【００４０】
最初に、音響空間化用広域システムを例示する図１を参照する。仮想場面をシミュレートするためのモジュール１ａは、信号の仮想音源として、例えば、モノラル音響等の音響オブジェクトを三次元空間の選択位置で定義し、また、音の方向を定義する。更に、仮想空間の幾何学的仕様を提供して、音の反響をシミュレートし得る。処理モジュール１１は、聴き手に関して、これら音源の内の１つ又は複数の管理（この聴き手に関する音源の仮想位置の定義）を適用する。それは、遅延及び／又は標準のフィルタ処理を適用することによって、反響等をシミュレートするための空間効果プロセッサを実現する。このように再現された信号は、音源の基本的寄与の空間的符号化用モジュール２ａに送信される。
【００４１】
これと並行して、自然集音は、実音源（モジュール１ｂ）に関して選択的に配置された１つ又は複数のマイクによって録音の枠組み内で行い得る。マイクによって集音された信号は、モジュール２ｂによって符号化される。集音及び符号化された信号は、中間表現フォーマット（モジュール３ｂ）に基づき変換された後、モジュール１ａによって生成されモジュール２ａによって符号化される（仮想音源から生じる）信号と、モジュール３によって混合し得る。その後、混合信号は、後で再生するために、送信（矢印ＴＲ）されるか又は媒体に記憶される。その後、スピーカを含む再生装置６で再生するために復号モジュール５に供給される。場合によっては、復号化ステップ５は、復号モジュール５の上流に処理モジュール４が提供されることで、音場を処理するステップ、例えば、回転が、先行することもある。
【００４２】
再生装置は、再生時、特に、三次元空間における音の方向を確実に認識するように、例えば、三次元（多重チャネル）構成において、球面上に配置された多数のスピーカの形態を取り得る。この目的のために、聴き手は、一般的に、アレイ状のスピーカによって形成され上述した聴覚点に対応する球の中心に自分自身を置く。変形例として、再生装置のスピーカは、平面（二次元パノラマ構成）に配置され、スピーカは、特に、円上に配置され、聴き手は、この円の中心に通常配置し得る。他の変形例では、再生装置は、“サラウンド”型（５．１）装置の形態を取り得る。最後に、有利な変形例では、再生装置は、再生される音をバイノーラル合成するための２つのヘッドホンを備えたヘッドセットの形態を取ることができ、これによって聴き手は、更に詳細に分かるように、三次元空間での音源の方向を認識し得る。また、三次元空間で認識するために２つのスピーカを備えたこのような再生装置は、聴き手からの選択した距離に２つのスピーカを配置したトランスオーラル再生装置の形態も取り得る。
【００４３】
次に図２を参照して、基本音源の三次元音響再生用の空間的符号化及び復号化について説明する。音源１乃至Ｎから生じる信号並びにその位置（実位置又は仮想位置）が、空間的符号化モジュール２に送信される。その位置は、入射角（聴き手から見た音源の方向）又はこの音源と聴き手との間の距離によって、同様に適切に定義し得る。このように符号化された複数の信号により、広域音場の多重チャンネル表現を得ることができる。符号化された信号は、図１で述べたように、音響再生装置６に送信され（矢印ＴＲ）三次元空間において音が再生される。
【００４４】
次に図３を参照して、音場の三次元空間の球面調和関数によるアンビソニック表現について以下説明する。如何なる音源もない原点Ｏを中心とする領域（半径Ｒの球）について考える。原点Ｏから球の点への各ベクトル
【数１】
2006506918.tif 000003
が、方位θｒ、高さδｒ、及び（原点Ｏからの距離に対応する）半径ｒで記述される球座標系を採用する。
【００４５】
この球内部の圧力音場
【数２】
2006506918.tif 000004
（Ｒが球の半径である場合、ｒ＜Ｒ）は、周波数領域において級数として表現し得るが、級数項は、角度関数
【数３】
2006506918.tif 000005
と、動径関数ｊｍ（ｋｒ）との重み付け積であり、従って、ｆが音響周波数、ｃが伝播媒体中の音速である場合、ｋ＝２πｆ／ｃである伝播項に依存する。
【００４６】
圧力音場は、
【数４】
2006506918.tif 000006
と表し得る。
【００４７】
重み付け係数の集合
【数５】
2006506918.tif 000007
は、暗黙裡に周波数に依存し、従って、対象領域の圧力音場を記述する。このため、これらの係数は、“球面調和関数成分”と呼ばれ、また、球面調和関数
【数６】
2006506918.tif 000008
の基数での音（又は圧力音場）に対する周波数表現を表す。
【００４８】
角度関数は、“球面調和関数”と呼ばれ、
【数７】
2006506918.tif 000009
によって定義される。ここで、Ｐｍｎ（ｓｉｎδ）は、度数ｍ、次数ｎのルジャンドル関数であり、δｐ，ｑは、クロネッカの記号である（ｐ＝ｑの時、１、そうでない場合、０に等しい）。
【００４９】
球面調和関数は、直交基数を形成し、高調波成分間において、一般的には、２つの関数ＦとＧとの間のスカラ積が、それぞれ
【数８】
2006506918.tif 000010
及び
【数９】
2006506918.tif 000011
によって定義される。
【００５０】
球面調和関数は、図４に示すように、次数ｍ並びに指数ｎ及びσの関数として拘束される実関数である。明部及び暗部は、球面調和関数の正及び負の値にそれぞれ対応する。次数ｍが高くなればなるほど、角度関数が（また、従って、関数間の相違が）大きくなる。動径関数ｊｍ（ｋｒ）は、球面ベッセル関数であり、その係数は、図５に次数ｍの幾つかの値に対して示す。
【００５１】
球面調和関数の基数によるアンビソニック表現の解釈は、以下のように与え得る。同様な次数ｍのアンビソニック成分は、最終的に、原点Ｏ（図３に表す球の中心）付近における圧力音場の次数ｍの“導関数”又は“モーメント”を表す。
【００５２】
特に、
【数１０】
2006506918.tif 000012
は、圧力のスカラ量を示し、他方、
【数１１】
2006506918.tif 000013
は、原点Ｏでの圧力勾配（又は、特定の速度）に関係する。これらの最初の４つの成分Ｗ，Ｘ，Ｙ，及びＺは、（次数０の成分Ｗの場合）全方位性のマイクを用いて、また（後続の他の３成分の場合）双方向性のマイクを用いて、自然集音中に得られる。極めて多くの音響変換器を用いることによって、適切な処理、特に等化によって、更に（１より大きくより高次数ｍの）アンビソニック成分を取得し得る。
【００５３】
（１より大きい）より高次数の追加成分を考慮し、従って、アンビソニック記述の角度分解能を大きくすることによって、原点Ｏを中心として、音波の波長について、より広い近隣領域に渡って圧力音場の近似値へのアクセスが得られる。従って、角度分解能（球面調和関数の次数）と、表現し得る半径方向の範囲（半径ｒ）との間には密接な関係が存在することを理解されたい。つまり、図３の原点Ｏから空間的に離れると、アンビソニック成分の数が多くなり（次数Ｍが高い）、また、これらアンビソニック成分の集合による音の表現が良くなる。しかしながら、音のアンビソニック表現は、原点Ｏから離れると、不充分になることも理解されたい。この影響は、特に（短波長の）高音響周波数に対して重要になる。従って、最大限の数のアンビソニック成分を得ることによって、音の再生が忠実でありその大きさが成分の総数と共に大きくなる知覚点周辺の空間領域を生成できるようにすることが関心の対象である。
【００５４】
以下、空間化音響符号化／送信／再生システムへの適用例について説明する。
【００５５】
実際には、アンビソニックシステムは、上述したように、球面調和関数の成分の部分集合を考慮する。後者が指数ｍ＜Ｍのアンビソニック成分を考慮する時、次数Ｍの系について語ることになる。スピーカを備えた再生装置によって再生を処理する時、これらのスピーカが水平面に配置される場合、指標ｍ＝ｎの調和関数だけが利用されることを理解されたい。他方、再生装置が、球の表面に渡り配置される（“多重チャネルの”）スピーカを含む場合、存在するスピーカと同じ数の調和関数を利用することが原則として可能である。
【００５６】
基準Ｓは、平面波によって搬送される圧力信号を示し、図３の球の中心に対応する点Ｏ（球座標の基数の原点）で集音される。波の入射角は、方位角θ及び仰角δによって記述される。この平面波に関連する音場成分の式は、関係式
【数１２】
2006506918.tif 000014
によって与えられる。
【００５７】
原点Ｏから距離ρにある近接音場源を符号化（シミュレート）する場合、近接音場が、第１近似値に合わせて球面波を放射することを考慮して、フィルタ
【数１３】
2006506918.tif 000015
が、波面形状を“湾曲する”ように適用される。音場の符号化成分は、
【数１４】
2006506918.tif 000016
となる。また、上記フィルタ
【数１５】
2006506918.tif 000017
に対する表現は、関係式
【数１６】
2006506918.tif 000018
によって与えられる。ここで、ω＝２πｆは、波の角周波数であり、ｆは、音響周波数である。
【００５８】
最終的に、これら後者の２つの関係式〔Ａ４〕及び〔Ａ５〕は、（シミュレートされる）仮想音源、及び近接音場における実音源の双方に対して、アンビソニック表現の音成分が、ここでは、ベッセル多項式であるベキｍの多項式の形態で数学的に（特に解析的に）表され、その変数（ｃ／２ｊωρ）が、音響周波数に反比例することを示す。
【００５９】
従って、以下のことを理解されたい。即ち、
平面波の場合、符号化は、実数の有限の利得によってのみ原信号と異なる信号を生成し、このことは、純粋に方向的符号化（関係式〔Ａ３〕）に対応する。また、
球面波（近接音場源）の場合、関係式〔Ａ５〕で表されるように、追加フィルタ
【数１７】
2006506918.tif 000019
が、アンビソニック成分に関する表現に、周波数に依存する複素振幅比を導入することによって距離キューを符号化する。
【００６０】
この追加フィルタは、“積分器”型であり、音響周波数がゼロに向かって減少するにつれて、増幅効果が大きくなり発散する（無限である）ことに留意されたい。図６は、各次数ｍに対して、低周波での利得の増加を示す（ここで、第１距離ρ＝１ｍである）。従って、任意の信号音にそれらを適用しようとする時、不安定で発散性のフィルタを取り扱っている。この発散は、高い次数値ｍの場合、更に決定的である。
【００６１】
特に、関係式〔Ａ３〕、〔Ａ４〕、及び〔Ａ５〕から、図６に示すように、近接音場で仮想音源をモデル化すると、高い次数ｍの場合特に決定的なように、低周波で発散性のアンビソニック成分を呈することが理解されるであろう。この発散は、低周波において、上述した“低音域ブースト”現象に対応する。また、それは、実音源に対する集音においても出現する。
【００６２】
特にこの理由によって、アンビソニック手法は、特に高次数ｍの場合、最先端技術では、音響処理において（理論的なもの以外の）具体的用途の経験がない。
【００６３】
再生時、アンビソニック表現で符号化される波面形状に適応するように、近接音場の補償が必要になることを特に理解されたい。図７において、再生装置は、上述した例では、聴覚点Ｐから、同一の距離Ｒに配置された複数のスピーカＨＰｉを含む。この図７では、
スピーカＨＰｉが位置する各点は、上述した再生点に対応し、
点Ｐは、上述した聴覚点であり、
これらの点は、上述した第２距離Ｒだけ離間されており、
一方、上述した図３では、
点Ｏは、上述した基準点に対応し、球面調和関数の基数の原点を形成し、
点Ｍは、基準点Ｏから、上述した第１距離ρに位置する（実又は仮想）音源の位置に対応する。
【００６４】
本発明によれば、近接音場の事前補償は、実際の符号化段階で導入され、この補償は、解析的形態のフィルタ
【数１８】
2006506918.tif 000020
を含み、上記アンビソニック成分
【数１９】
2006506918.tif 000021
に適用される。
【００６５】
本発明により提供される利点の１つによれば、図６で効果が現れる増幅
【数２０】
2006506918.tif 000022
は、符号化に引き続き適用されるフィルタ
【数２１】
2006506918.tif 000023
の減衰を通して補償される。特に、この補償フィルタ
【数２２】
2006506918.tif 000024
の係数は、音響周波数と共に増加し、特に、低周波の場合、ゼロに近づく。利点として、この事前補償は、符号化から直ぐに行われ、送信されるデータが、低周波に対して発散しないことを保証する。
【００６６】
補償フィルタに関わる距離Ｒの物理的意味を示すために、例示により、信号音の集音時の初期実平面波を考慮する。この遠方音源の近接音場効果をシミュレートするために、関係式〔Ａ４〕に示したように、関係式〔Ａ５〕の第１フィルタを適用する。この時、距離ρは、近接仮想音源Ｍと、図３の球基数の原点を表す点Ｏとの間の距離を表す。このように、近接音場シミュレーション用の第１フィルタは、上述した距離ρにある仮想音源の存在をシミュレートするために適用される。しかしながら、一方では、上述したように、このフィルタの係数の項は、低周波で発散し（図６）、他方では、上記距離ρは、再生装置のスピーカと聴覚点Ｐとの間の距離を必ずしも表さない（図７）。本発明によれば、符号化時、事前補償が、上述したように
【数２３】
2006506918.tif 000025
型のフィルタを含んで適用され、これによって、一方では、拘束信号の送信が可能になり、他方では、図７に示すように、スピーカＨＰｉを用いて音響再生のために、符号化から直ちに距離Ｒを選択し得る。特に、認識されたいことは、集音時、原点Ｏから距離ρに置かれた仮想音源をシミュレートした場合、再生時（図７）（スピーカＨＰｉから距離Ｒの）聴覚点Ｐに位置する聴き手は、聞き取る際、聴覚点Ｐから距離ρに配置され集音時シミュレートされる仮想音源に対応する音源Ｓの存在を認識する。
【００６７】
従って、符号化段階で（距離Ｒに配置された）スピーカの近接音場の事前補償は、距離ρに配置された仮想音源のシミュレートされる近接音場効果と組み合わせ得る。符号化時、一方では、近接音場のシミュレーションによる、他方では、近接音場の補償による総フィルタは、最終的に利用されるが、このフィルタの係数は、関係式
【数２４】
2006506918.tif 000026
によって解析的に表し得る。
【００６８】
関係式〔Ａ１１〕によって与えられる総フィルタは、安定しており、また、図８に示すように、本発明に基づく空間的アンビソニック符号化における“距離符号化”部分を構成する。これらフィルタの係数は、周波数に対する単調な伝達関数に対応し、高周波で値１に、低周波で、値（Ｒ／ρ）ｍに近づく。図９において、フィルタのエネルギースペクトル
【数２５】
2006506918.tif 000027
は、（距離ρ＝１ｍでここに配置された）仮想音源の音場効果による符号化成分の増幅を伝達し、（距離Ｒ＝１．５ｍに配置された）スピーカの音場が事前補償される。従って、デシベル単位の増幅は、ρ＜Ｒの時（図９の場合）正であり、ρ＞Ｒの時（ρ＝３ｍ及びＲ＝１．５ｍの図１０の場合）負である。空間化再生装置では、聴覚点とスピーカＨＰｉとの間の距離Ｒは、実際には、１又は数メートルのオーダである。
【００６９】
また図８では、通常の方向パラメータθ及びδとは別に、符号化に含まれる距離に関するキューが、送信されることを理解されるであろう。このように、球面調和関数
【数２６】
2006506918.tif 000028
に対応する角度関数は、方向的符号化に対して保持される。
【００７０】
しかしながら、本発明の骨子の範囲内において、アンビソニック成分に、それらの次数ｍの関数として適用される総フィルタ（近接音場補償、及び場合によっては、近接音場のシミュレーション）
【数２７】
2006506918.tif 000029
を更に提供して、図８に示すように、距離符号化を実現する。音響デジタル領域におけるこれらのフィルタの実施形態は、後で詳述する。
【００７１】
これらのフィルタは、まさしく距離符号化（ｒ）から直ちに、また、方向符号化（θ，δ）の前でさえ、適用し得ることを特に留意されたい。このように、上述したステップａ）及びｂ）は、共に同一の広域ステップに持ち込むことができたり、入れ替えたり（方向符号化及び補償フィルタ処理の後、距離符号化）できることを理解されたい。従って、本発明に基づく方法は、ステップａ）及びｂ）の連続した時間的具体化に限定されない。
【００７２】
図１１Ａは、総次数Ｍ＝１５の系及び３２個のスピーカ上での再生に対して、（距離パラメータは、図９と同じ状態の）水平面において球面波を補償して、近接音場の再現の（上方から見た）視覚化を表す。図１１Ｂには、再生空間では、図７の聴覚点Ｐに対応する集音空間の点から距離ρに位置する近接音場源からの初期音波の伝播を表す。図１１Ａにおいて留意されたいことは、（図式化した頭部で表される）聴き手は、図１１Ｂの聴覚点Ｐから距離ρに位置する同一の地理的位置に仮想音源を厳密に特定し得ることである。
【００７３】
このように、符号化された波面の形状は、復号化及び再生後の状態に適合することが実際に検証される。しかしながら、図１１Ａに示したような点Ｐの右側での干渉は、顕著であり、この干渉は、スピーカ（従って、考慮されるアンビソニック成分）の数が、スピーカによって範囲が定められる面全体に渡って含まれる波面の完璧な再現には不充分であるという事実によるものである。
【００７４】
以下、例として、本発明の骨子の範囲内で本方法を実現するための音響デジタルフィルタの取得について説明する。
【００７５】
上述したように、符号化から直ちに補償される近接音場効果をシミュレーションしようとする場合、形態
【数２８】
2006506918.tif 000030
のフィルタが、音のアンビソニック成分に適用される。
【００７６】
関係式〔Ａ５〕によって与えられる近接音場のシミュレーションの表現から、遠方音源（ρ＝∞）に対して、関係式〔Ａ１１〕は、単に
【数２９】
2006506918.tif 000031
になることが明白である。
【００７７】
従って、この後者の関係式〔Ａ１２〕から明白なことは、シミュレートされる音源が遠方音場（遠方音源）で放射する事例は、関係式〔Ａ１１〕で公式化されるように、フィルタに対する一般表現の特定の事例に過ぎないことである。
【００７８】
音響デジタル処理の分野では、連続時間アナログ領域におけるこのフィルタの解析的表現からデジタルフィルタを定義する有利な方法は、“双一次変換”から成る。
【００７９】
関係式〔Ａ５〕は、ラプラス変換の形態でまず表され、これは、
【数３０】
2006506918.tif 000032
に対応する。ここで、τ＝ρ／ｃ（ｃは、媒体中の音速であり、通常、空気中で３４０ｍ／ｓ）である。
【００８０】
双一次変換は、サンプリング周波数ｆｓに対して、関係式〔Ａ１１〕を、ｍが奇数の場合、
【数３１】
2006506918.tif 000033
の形態で、また、ｍが偶数の場合、
【数３２】
2006506918.tif 000034
の形態で、提示することにある。ここで、ｚは、上記関係式〔Ａ１３〕に関して、
【数３３】
2006506918.tif 000035
によって定義され、
【数３４】
2006506918.tif 000036
である。ここで、ｘ＝ａの場合、α＝４ｆｓＲ／ｃであり、ｘ＝ｂの場合、α＝４ｆｓρ／ｃである。
【００８１】
Ｘｍ，ｑは、ベッセル多項式のｑ個の連続根
【数３５】
2006506918.tif 000037
であり、ｍが奇数の時、それらの実部、（コンマで区切られた）それらの係数、及びそれらの（実数）値のそれぞれの形態で、様々な次数ｍについて、以下の表１に示す。
表１：ＭＡＴＬＡＢ（Ｃ）計算ソフトウェアを用いて計算されるベッセル多項式の値Ｒｅ〔Ｘｍ，ｑ〕，｜Ｘｍ，ｑ｜，（及びｍが奇数の時、Ｒｅ〔Ｘｍ，ｍ〕）
【表１】
2006506918.tif 000038
2006506918.tif 000039
【００８２】
このように、デジタルフィルタは、上記により与えられた関係式〔Ａ１４〕を用いて、カスケード状の次数２（ｍが偶数）のセル及び追加セル（ｍが奇数）を提供することによって、表１の値を用いて配置される。
【００８３】
従って、デジタルフィルタは、後述するように容易にパラメータ化し得る無限大のインパルス応答形態で具体化される。留意されたいことは、有限のインパルス応答形態の実例は、予見でき、解析的公式から伝達関数の複素スペクトルを計算し、そして、そこから逆フーリエ変換によって有限のインパルス応答を演繹することにある。その後、畳み込み演算がフィルタ処理に対して適用される。
【００８４】
このように、符号化時、この近接音場の事前補償を導入することによって、修正アンビソニック表現（図８）は、
【数３６】
2006506918.tif 000040
の形態で定義され、周波数領域において表される信号を、送信可能な表現として採用する。
【００８５】
上述したように、Ｒは、補償近接音場効果に関連する基準距離であり、ｃは、音速（通常、空気中で３４０ｍ／ｓ）である。この修正アンビソニック表現は、（図１の矢印ＴＲに近接して“囲まれた”データを送信することによって概略的に表される）同じ拡張特性を持ち、また、通常のアンビソニック表現と同じ音場回転変換（図１のモジュール４）に従う。
【００８６】
以下、アンビソニック受信信号の復号化の場合に実行される動作を示す。
【００８７】
最初に、復号化動作は、上述した基準距離Ｒとは異なる半径Ｒ２の任意の再生装置に適応し得ることを示す。この目的のために、前述したような
【数３７】
2006506918.tif 000041
型のフィルタが適用されるが、ρ及びＲの代わりに、距離パラメータはＲ及びＲ２である。特に、パラメータＲ／ｃのみが、符号化と復号化との間に記憶（及び／又は送信）される必要があることを留意されたい。
【００８８】
図１２において、そこに示すフィルタ処理モジュールは、例えば、再生装置の処理ユニットに備えられる。受信アンビソニック成分は、符号化時、第２の距離としての基準距離Ｒ１が事前補償されている。しかしながら、再生装置は、聴覚点Ｐから第３の距離Ｒ２に配置された複数のスピーカを含み、この第３の距離Ｒ２は、上記第２の距離Ｒ１とは異なる。そして、図１２のフィルタ処理モジュールは、
【数３８】
2006506918.tif 000042
の形態で、データ受信時、距離Ｒ２での再生の事前補償を距離Ｒ１に適合する。勿論、上述したように、再生装置は、パラメータＲ１／ｃも受信する。
【００８９】
更に本発明は、音場（実及び／又は仮想音源）の幾つかのアンビソニック表現を混合することを可能にし、その基準距離Ｒは、異なる（場合によっては、無限基準距離が遠方音源に対応する）ことを認識されたい。好適には、最短の基準距離でのこれら全ての音源の事前補償は、アンビソニック信号を混合する前に、フィルタ処理され、これによって、再生時、音響解放の正確な定義を取得し得る。
【００９０】
再生時、（光プロジェクタが光学系の選択方向を照明するように）空間的に選択された方向に対する音響強化効果での所謂“音響集束”処理の枠組み内で、（アンビソニック成分の重み付けによる）音響集束の行列処理を含み、近接音場事前補償を集束処理と組み合わせるようにして、距離符号化を有利に適用する。
【００９１】
以下、再生時、スピーカの近接音場を補償するアンビソニック復号化方法について説明する。
【００９２】
成分
【数３９】
2006506918.tif 000043
から、アンビソニック形式に基づき、及び、図７の再生点Ｐに対応する聴き手の“理想的”配置を提供する再生装置のスピーカを用いることによって符号化される音場を再現するために、各スピーカによって放射される波は、以下のように、再生装置の中心でアンビソニック音場の事前の“再符号化”処理によって定義される。
【００９３】
この“再符号化”に対して、簡単のために、音源が遠方音場で放射するとまず考える。
【００９４】
また図７では、指数ｉ及び入射角（θｉ及びδｉ）のスピーカによって放射される波は、信号Ｓｉで供給される。このスピーカは、その寄与
【数４０】
2006506918.tif 000044
を通して成分Ｂ’ｍｎの再現に加わる。
【００９５】
指数ｉのスピーカに関連する符号化係数のベクトルｃｉは、関係式
【数４１】
2006506918.tif 000045
によって表される。
【００９６】
Ｎ個のスピーカの集合から発する信号のベクトルＳは、式
【数４２】
2006506918.tif 000046
によって与えられる。
【００９７】
（最終的に“再符号化”行列に対応する）これらＮ個のスピーカに対する符号化行列は、関係式
【数４３】
2006506918.tif 000047
によって表される。ここで、各項ｃｉは、上記関係式〔Ｂ１〕に基づくベクトルを表す。
【００９８】
従って、アンビソニック音場Ｂ’の再現は、関係式
【数４４】
2006506918.tif 000048
によって定義される。
【００９９】
このように、関係式〔Ｂ４〕は、再生前、再符号化動作を定義する。最終的に、復号化は、このように、一般的関係式、
Ｂ’＝Ｂ・・・〔Ｂ６〕
を定義するように、
【数４５】
2006506918.tif 000049
の形態で、再生装置によって受信される原アンビソニック信号を再符号化信号
【数４６】
2006506918.tif 000050
と比較することにある。
【０１００】
これには、特に、関係式
Ｓ＝Ｄ．Ｂ・・・〔Ｂ７〕
を満たす復号化行列Ｄの係数を決定することが含まれる。
【０１０１】
好適には、スピーカの数は、復号化されるアンビソニック成分の数以上であり、また、復号化行列Ｄは、再符号化行列Ｃの関数として、
Ｄ＝ＣＴ（ＣＣＴ）−１・・・〔Ｂ８〕
の形態で表し得る。ここで、記号ＣＴは、行列Ｃの転置行列に対応する。
【０１０２】
留意されたいことは、各周波数領域に対して異なる基準を満たす復号化の定義が可能であり、これによって、聞き取り条件の関数として、特に、再生時、図３の球の中心Ｏでの位置決めの制約に関して、最適な再生を提供し得ることである。この目的のために、各アンビソニック成分での段階的周波数等化によって、簡単なフィルタ処理が有利に提供される。
【０１０３】
しかしながら、元々符号化された波の再現を得るために、スピーカに対する遠方音場の仮定の修正が必要である。つまり、上述した再符号化行列Ｃでそれらの近接音場効果を表す必要があり、また、この新しいシステムを反転して復号器を定義する必要がある。この目的のために、（図７の点Ｐから同一の距離Ｒに配置された）スピーカの同心性を仮定すると、全てのスピーカは、
【数４７】
2006506918.tif 000051
型の各アンビソニック成分では、同じ近接音場効果
【数４８】
2006506918.tif 000052
を有する。対角行列の形態で近接音場項を導入することによって、上述した関係式〔Ｂ４〕は、
【数４９】
2006506918.tif 000053
となる。
【０１０４】
上述した関係式〔Ｂ７〕は、
【数５０】
2006506918.tif 000054
となる。
【０１０５】
従って、行列化動作には、各成分
【数５１】
2006506918.tif 000055
の近接音場を補償するフィルタ処理動作が先行し、これは、上述したように、関係式〔Ａ１４〕を基準にして、デジタル形態で実現し得る。
【０１０６】
実際には、“再符号化”行列Ｃは、再生装置に特有であることを思い起こされたい。その係数は、初期的には、パラメータ化によって、また、再生装置が所定の励起に反応する音響特性によって決定し得る。同様に、復号化行列Ｄは、再生装置に特有である。その係数は、関係式〔Ｂ８〕によって決定し得る。前の表記法を続けると、
【数５２】
2006506918.tif 000056
が、事前補償アンビソニック成分の行列である場合、これら後者は、
【数５３】
2006506918.tif 000057
として、
【数５４】
2006506918.tif 000058
の行列形態で、再生装置に送信し得る。
【０１０７】
その後、再生装置は、
【数５５】
2006506918.tif 000059
で、スピーカＨＰｉに供給するための信号Ｓｉを形成するように、事前補償されたアンビソニック成分に復号化行列Ｄを適用することによって、行列形態
【数５６】
2006506918.tif 000060
（送信成分の列ベクトル）で受信されるデータを復号化する。
【０１０８】
また図１２では、復号化動作が、基準距離Ｒ１と異なる半径Ｒ２の再生装置に適応しなければならない場合、本来の上述した復号化に先立つ適合用のモジュールは、半径Ｒ２の再生装置にそれを適合させるように、各アンビソニック成分
【数５７】
2006506918.tif 000061
をフィルタ処理し得る。その後、本来の復号化動作が、上述したように、関係式〔Ｂ１１〕について行われる。
【０１０９】
バイノーラル合成への本発明の適用について、以下に述べる。
【０１１０】
バイノーラル合成装置の２つのヘッドホンを備えたヘッドセットを有する聴き手を表す図１３Ａを参照する。聴き手の２つの耳は、空間のそれぞれの点ＯＬ（左耳）及びＯＲ（右耳）に置かれる。聴き手の頭の中心は、点Ｏに置かれ、また、聴き手の頭の半径は、値ａである。音源は、聴き手の頭の中心から距離ｒ（また、それぞれ右耳から距離ｒＲ及び左耳からｒＬ）に位置する空間の点Ｍにおいて、聴覚的に知覚されなければならない。更に、点Ｍに配置された音源の方向は、ベクトル
【数５８】
2006506918.tif 000062
によって定義される。
【０１１１】
一般的に、バイノーラル合成は、以下のように定義される。
【０１１２】
聴き手は、各々特有な耳の形状を有する。この聴き手による空間音の知覚は、この聴き手に特有な耳の形状（特に、耳介の形状及び頭の寸法）の関数として出生時から学習によって行われる。空間音の知覚は、音が一方の耳に他方の耳より先に到達するという事実によってとりわけ顕在化し、これにより、バイノーラル合成を適用する再生装置の各ヘッドホンによって放射される信号間に遅延τが生じる。
【０１１３】
再生装置は、同一の聴き手に対して、自分の頭の中心から同じ距離Ｒにおいて、自分の頭周辺の音源を掃引することによって、初期的にパラメータ化される。従って、この距離Ｒは、上述した”再生点“と聴覚点（ここでは聴き手の頭の中心Ｏ）との間の距離であり得ることを理解されたい。
【０１１４】
以下において、指標Ｌは、左耳に接するヘッドホンによって再生される信号と関連し、指標Ｒは、右耳と接するヘッドホンによって再生される信号と関連する。図１３Ｂにおいて、遅延は、別個のヘッドホン用に信号を生成する各経路に対する初期信号Ｓに適用し得る。これらの遅延τＬ及びτＲは、最大遅延τｍａｘに依存する。前述したように、ａが聴き手の頭の半径に対応し、ｃが音速に対応する場合、最大遅延τｍａｘは、ここでは、比率ａ／ｃに対応する。特に、これらの遅延は、点Ｏ（頭の中心）から点Ｍ（図１３Ａでは、音が再生される音源の位置）までの距離と各耳からこの点Ｍまでの距離との差の関数として定義される。利点として、それぞれの利得ｇＬ及びｇＲは、更に、各経路に適用され、点Ｏから点Ｍまでの距離と各耳から点Ｍまでの距離との比に依存する。各経路２Ｌ及び２Ｒに適用されるそれぞれのモジュールは、本発明の骨子の範囲内での近接音場事前補償ＮＦＣ（“近接音場補償”の意）によるアンビソニック表現で、各経路の信号を符号化する。従って、本発明の骨子の範囲内の方法の実行によって、音源Ｍから生じる信号は、それらの方向（方位角θＬ及びθＲ並びに仰角δＬ及びδＲ）によってだけでなく、音源Ｍから各耳ｒＬ及びｒＲを離間する距離の関数としても、定義し得ることを理解されたい。このように符号化された信号は、各経路５Ｌ及び５Ｒ用のアンビソニック復号モジュールを含む再生装置に送信される。従って、アンビソニック符号化／復号化は、複製形態で、（ここでは“Ｂフォーマット”型の）バイノーラル合成による再生において、各経路（左ヘッドホン、右ヘッドホン）に対して、近接音場補償と共に適用される。近接音場補償は、各耳と再生される音源の位置Ｍとの間の距離ｒＬ及びｒＲを第１距離ρとして、各経路に対して行われる。
【０１１５】
以下、アンビソニック表現での集音の文脈内において、本発明の骨子の範囲内の補償の適用について説明する。
【０１１６】
図１４において、マイク１４１は、音響圧力を集音可能であり電気信号Ｓ１，・・・，ＳＮを再現可能な複数の変換器カプセルを含む。カプセルＣＡＰｉは、所定の半径ｒの球（ここでは、例えば、ピンポン玉等の剛体球）上に配置される。カプセルは、球の全面において一定の間隔で離間される。実際には、カプセルの数Ｎは、アンビソニック表現の所望の次数Ｍの関数として選択される。
【０１１７】
以下、剛体球上に配置されるカプセルを含むマイクの文脈内で、アンビソニックの文脈において符号化から直ちに近接音場効果を補償する方法について述べる。従って、分かることは、近接音場の事前補償が、上述したように、仮想音源シミュレーションに対してだけでなく、集音時にも適用され、また、より一般的に、アンビソニック表現を含む全種類の処理と近接音場事前補償を組み合わせることによって適用し得ることである。
【０１１８】
（受信音波の回折を生じ易い）剛体球が存在する場合、上述の関係式〔Ａ１〕は、
【数５９】
2006506918.tif 000063
となる。
【０１１９】
球ハンケル関数の導関数ｈ−ｍは、繰り返し則、即ち、
【数６０】
2006506918.tif 000064
に従う。
【０１２０】
関係式
【数６１】
2006506918.tif 000065
によって与えられる投射及び等化演算を実行することによって、球の表面における圧力音場から初期音場のアンビソニック成分
【数６２】
2006506918.tif 000066
を演繹する。
【０１２１】
この式では、ＥＱｍは、重み付けＷｍを補償する等化器フィルタであり、重み付けＷｍは、カプセルの方向性に関係し、また更に、剛体球による回折を含む。
【０１２２】
このフィルタＥＱｍの式は、以下の関係式
【数６３】
2006506918.tif 000067
によって与えられる。
【０１２３】
この等化フィルタの係数は、安定でなく、無限利得が、超低周波で得られる。更に、球面調和関数成分は、それら自体、音場が平面波、つまり、前述したように、遠方音源から生じる波の伝播に限定されない場合、有限の振幅ではないことに留意することが適当である。
【０１２４】
更に、固体球に埋め込むカプセルを提供するよりもむしろ、カージオイド型カプセルを提供する場合、遠方音場の方向性は、式
Ｇ（θ）＝α＋（１−α）ｃｏｓθ・・・［Ｃ５］
によって与えられる。
【０１２５】
“音響的に透過な”支持体上に搭載されるこれらのカプセルを考えることによって、補償される重み付け項は、
【数６４】
2006506918.tif 000068
となる。
【０１２６】
関係式〔Ｃ６〕によって与えられるこの重み付けの解析的な逆数に対応する等化フィルタの係数は、超低周波に対して発散性があることが再度明白になる。
【０１２７】
一般的に、任意の種類の方向センサの場合、センサの方向性に関する重み付けＷｍを補償するフィルタＥＱｍの利得は、低い音響周波数に対して無限であることが示される。図１４では、近接音場事前補償は、関係式
【数６５】
2006506918.tif 000069
によって与えられる等化フィルタＥＱｍに対する実際の式に有利に適用される。
【０１２８】
従って、信号Ｓ１乃至ＳＮは、マイク１４１から回収される。適宜、これらの信号の前置等化は、処理モジュール１４２によって適用される。モジュール１４３は、行列形態で、アンビソニック文脈でこれらの信号を表し得る。モジュール１４４は、マイク１４１の球半径ｒの関数として表されるアンビソニック成分に関係式〔Ｃ７〕のフィルタを適用する。近接音場補償は、第２距離として基準距離Ｒに対して行われる。従って、モジュール１４４によってフィルタ処理される符号化信号は、場合によって、基準距離を表すパラメータＲ／ｃと共に送信し得る。
【０１２９】
従って、それぞれ、近接音場仮想音源の生成、実音源から生じる信号音の集音、又は（スピーカの近接音場効果を補償するための）再生、に関連する様々な実施形態では、本発明の骨子の範囲内の近接音場補償が、アンビソニック表現を含む全種類の処理に適用し得ることが明白である。この近接音場補償は、音源の方向及びその距離が有利に考慮されなければならない多数の音の文脈にアンビソニック表現を適用し得る。更に、アンビソニック文脈内で全種類の音響現象（近接又は遠方音場）を表現する可能性が、アンビソニック成分の有限の実数値に対する制限のために、この事前補償によって保証される。
【０１３０】
勿論、本発明は、一例として上述した実施形態に限定されず、他の変形例に適用される。
【０１３１】
従って、近接音場事前補償は、符号化時、近接音源に対しても遠方音源と同様取り込み得ることを理解されたい。後者の場合（遠方音源、及び平面波の受信）では、上述の距離ρは、前述のフィルタＨｍに対する式を実質的に修正することなく、無限であると見なされる。従って、一般的に、遅延拡散音場（遅延反響）のモデル化に使用可能な非相関信号を提供する空間効果プロセッサを用いた処理は、近接音場事前補償と組み合わせ得る。これらの信号は、同様なエネルギーであると見なされ、また、全方位性の成分
【数６６】
2006506918.tif 000070
に対応する拡散音場（図４）の割当てに対応すると見なし得る。次に、（選択された次数Ｍの）様々な球面調和関数成分は、各アンビソニック成分に対する利得修正を適用することによって構成でき、また、スピーカの近接音場補償が、（図７に示すように、基準距離Ｒが聴覚点からスピーカを離間する状態で）適用される。
【０１３２】
勿論、本発明の骨子の範囲内の符号化の原理は、単極音源（実又は仮想）及び／又はスピーカ以外の放射モデルに一般化できる。具体的には、任意の放射形状（特に空間に広がった音源）は、基本的点音源の連続的分布の統合によって表現し得る。
【０１３３】
更に、再生の文脈では、任意の再生の文脈に近接音場補償を適応し得る。この目的のために、伝達関数（音響が再生される空間での実際の伝播を考慮した各スピーカに対する近接音場球面調和関数成分の再符号化）を計算するだけでなく、この再符号化を反転して復号化の再定義も行い得る。
【０１３４】
以上、アンビソニック成分を含む行列系を適用した復号化方法について述べた。変形例では、高速フーリエ変換（円又は球）による一般的処理を行い、計算時間、及び復号化処理に必要な（メモリに関する）計算処理資源を制限する。
【０１３５】
図９及び１０で述べたように、近接音場源の距離ρに対する基準距離Ｒの選択は、音響周波数の様々な値に対する利得の相違をもたらすことが分かる。事前補償による符号化の方法は、音響デジタル圧縮と組み合わせて、各周波数サブバンドに対する利得の量子化と調整が可能になることが分かる。
【０１３６】
利点として、本発明は、全種類の音響空間化システムに適用し、特に、“仮想現実”型の用途（三次元空間の仮想場面を介したナビゲーション、三次元音響空間化によるゲーム、インターネット網上で発する“チャット”型の会話）に、インターフェイスの音響装備に、音楽を記録、混合、及び再生するための音響編集ソフトウェアに、適用する。また、更に、ミュージカル又は映画撮影の集音や、あるいは、例えば、音響装備された“ウェブカメラ”用などのインターネットにおける音響雰囲気を送信するための三次元マイクを用いた集音にも適用する。
【図面の簡単な説明】
【０１３７】
【図１】空間化再生装置による符号化、送信、復号化、及び再生により、仮想音源のシミュレーションによって、信号音を集音し生成するためのシステムを示す概略図。
【図２】強度的に、また、信号が生じる音源の位置的に双方で定義される信号の符号化をより厳密に表す図。
【図３】球座標でのアンビソニック表現に含まれるパラメータを示す図。
【図４】様々な次数の球面調和関数
【数６７】
2006506918.tif 000071
の球座標の基準枠における三次元計量による表現を示す図。
【図５】次数ｍの連続値の球面ベッセル関数であって、圧力音場のアンビソニック表現に用いられる動径関数ｊｍ（ｋｒ）の係数における変化を示すグラフ。
【図６】特に低周波での様々な連続次数ｍに対する近接音場効果による増幅を表す図。
【図７】上記聴覚点（参照Ｐ）、上記第１距離（参照ρ）、及び上記第２距離（参照Ｒ）に複数のスピーカＨＰｉを含む再生装置を表す概略図。
【図８】本発明に基づく、方向符号化及び距離符号化によるアンビソニック符号化に含まれるパラメータを表す概略図。
【図９】仮想音源の第１距離ρ＝１ｍ、及び第２距離Ｒ＝１．５ｍに位置するスピーカの事前補償に対してシミュレートされる補償及び近接音場フィルタのエネルギースペクトルを表す図。
【図１０】仮想音源の第１距離ρ＝３ｍ、及び距離Ｒ＝１．５ｍに位置するスピーカの事前補償に対してシミュレートされる補償及び近接音場フィルタのエネルギースペクトルを表す図。
【図１１Ａ】本発明に基づく水平面における球面波の補償による近接音場の再現を表す図。
【図１１Ｂ】図１１Ａと比較して、音源Ｓから生じる初期波面を表す図。
【図１２】聴覚点から第３距離Ｒ２に配置された複数のスピーカを含む再生装置に合わせて、受信され事前補償されたアンビソニック成分を、第２距離として基準距離Ｒに対する符号化に適応するためのフィルタ処理モジュールを表す概略図。
【図１３Ａ】再生時、音源が近接音場で放射する状態で、バイノーラル合成を適用する再生装置を用いて、聴き手に対する音源Ｍの配置を表す概略図。
【図１３Ｂ】アンビソニック符号化／復号化が組み合わせられる図１３Ａのバイノーラル合成の枠組みにおいて近接音場効果による符号化及び復号化のステップを表す概略図。
【図１４】本発明に基づき、アンビソニック符号化、等化、及び近接音場補償によって、例示により、球上に配置された複数の圧力センサを含むマイクから生じる信号の処理を表す概略図。
【符号の説明】
【０１３８】
１０・・・仮想場面の記述、１１・・・音源の管理、２ａ・・・基本的寄与の空間的符号化、４・・・音場の操作（回転）、５・・・復号化、１ｂ・・・自然集音、２ｂ・・・音響符号化、３ｂ・・・中間表現フォーマット、６・・・バイノーラル／トランスオーラル。
【図１】
2006506918.tif 000075
【図２】
2006506918.tif 000076
【図３】
2006506918.tif 000077
【図４】
2006506918.tif 000078
【図５】
2006506918.tif 000079
【図６】
2006506918.tif 000080
【図７】
2006506918.tif 000081
【図８】
2006506918.tif 000082
【図９】
2006506918.tif 000083
【図１０】
2006506918.tif 000084
2006506918.tif 000085
【図１２】
2006506918.tif 000086
2006506918.tif 000087
【図１３Ｂ】
2006506918.tif 000088
【図１４】
2006506918.tif 000089
【国際調査報告】
2006506918.tif 000090
2006506918.tif 000091
2006506918.tif 000092
2006506918.tif 000093
2006506918.tif 000094
2006506918.tif 000095