特許 6395811 高次アンビソニックス表現を圧縮および圧縮解除する方法および装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6395811

(24)【登録日】2018年9月7日

(45)【発行日】2018年9月26日

(54)【発明の名称】高次アンビソニックス表現を圧縮および圧縮解除する方法および装置

(51)【国際特許分類】

   G10L 19/008 20130101AFI20180913BHJP

   G10L 25/51 20130101ALI20180913BHJP

【ＦＩ】

   G10L19/008 100

   G10L25/51 400

【請求項の数】22

【全頁数】36

(21)【出願番号】特願2016-509473(P2016-509473)

(86)(22)【出願日】2014年4月24日

(65)【公表番号】特表2016-520864(P2016-520864A)

(43)【公表日】2016年7月14日

(86)【国際出願番号】EP2014058380

(87)【国際公開番号】WO2014177455

(87)【国際公開日】20141106

【審査請求日】2017年3月1日

(31)【優先権主張番号】13305558.2

(32)【優先日】2013年4月29日

(33)【優先権主張国】EP

(73)【特許権者】

【識別番号】510185767

【氏名又は名称】ドルビー・インターナショナル・アーベー

(74)【代理人】

【識別番号】100107766

【弁理士】

【氏名又は名称】伊東 忠重

(74)【代理人】

【識別番号】100070150

【弁理士】

【氏名又は名称】伊東 忠彦

(74)【代理人】

【識別番号】100091214

【弁理士】

【氏名又は名称】大貫 進介

(72)【発明者】

【氏名】クルーガー，アレクサンダー

(72)【発明者】

【氏名】コルドン，スフエン

【審査官】 山下 剛史

(56)【参考文献】

【文献】 特開２０１２−１３３３６６（ＪＰ，Ａ）

【文献】 国際公開第２０１２／０５９３８５（ＷＯ，Ａ１）

【文献】 国際公開第２０１１／１１７３９９（ＷＯ，Ａ１）

【文献】 欧州特許出願公開第２０９４０３２（ＥＰ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ１０Ｌ １９／００−１９／２６

Ｈ０４Ｓ １／００− ７／００

ＡＥＳ Ｅ−Ｌｉｂｒａｒｙ

(57)【特許請求の範囲】

【請求項1】

所定数（Ｉ）の知覚符号化処理を使用して音場のＨＯＡと称する高次アンビソニックス表現をＨＯＡ係数列の入力される時間フレーム

[この文献は図面を表示できません]

を用いて圧縮する方法であって、
前記方法は、フレーム単位で行われ、
−現在のフレーム

[この文献は図面を表示できません]

に対して、支配的な方向のセット

[この文献は図面を表示できません]

および対応する検出された方向性信号のインデックスのデータセット

[この文献は図面を表示できません]

を推定するステップ（１３）と、
−前記現在のフレームのＨＯＡ係数列を分解するステップ（１４、１５）であって、非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））であって、支配的な方向推定値の前記セット

[この文献は図面を表示できません]

に含まれる各々の方向と前記方向性信号のインデックスの各々の遅延したデータセット

[この文献は図面を表示できません]

とを用いた、前記非所定数（Ｍ）が前記所定数（Ｉ）よりも小さい、該非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））と、前記所定数（Ｉ）と前記非所定数（Ｍ）との差に対応する低減された数のＨＯＡ係数列によって表現される残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））と、対応する前記低減された数の残差のアンビエントＨＯＡ係数列のインデックスのデータセット

[この文献は図面を表示できません]

と、に分解する、該分解するステップ（１４、１５）と、
−前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））および前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））のＨＯＡ係数列を前記所定数（Ｉ）に対応する数のチャンネルに割り当てるステップ（１６）であって、前記割り当てのために、前記方向性信号のインデックスの前記遅延されたデータセット

[この文献は図面を表示できません]

および前記低減された数の残差のアンビエントＨＯＡ係数列のインデックスの前記データセット

[この文献は図面を表示できません]

が使用される、該割り当てるステップと、
−関連するフレーム（Ｙ（ｋ−２））の前記チャンネルを知覚符号化するステップ（１７）であって、符号化された圧縮されたフレーム

[この文献は図面を表示できません]

が得られる、該知覚符号化するステップ（１７）と、
を含む、前記方法。

【請求項2】

所定数（Ｉ）の知覚符号化処理を使用して音場のＨＯＡと称する高次アンビソニックス表現をＨＯＡ係数列の入力される時間フレーム

[この文献は図面を表示できません]

を用いて圧縮する装置であって、当該装置は、フレーム単位の処理を実行し、
−現在のフレーム

[この文献は図面を表示できません]

に対して、支配的な方向のセット

[この文献は図面を表示できません]

および対応する検出された方向性信号のインデックスのデータセット

[この文献は図面を表示できません]

を推定するように構成された手段（１３）と、
−前記現在のフレームのＨＯＡ係数列を分解するように構成された手段（１４、１５）であって、非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））であって、支配的な方向推定値の前記セット

[この文献は図面を表示できません]

に含まれる各々の方向と、前記方向性信号のインデックスの各々の遅延したデータセット

[この文献は図面を表示できません]

とを用いた、前記非所定数（Ｍ）が前記所定数（Ｉ）よりも小さい、該非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））と、前記所定数（Ｉ）と前記非所定数（Ｍ）との差に対応する低減された数のＨＯＡ係数列によって表現される残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））と、対応する前記低減された数の残差のアンビエントＨＯＡ係数列のインデックスのデータセット

[この文献は図面を表示できません]

と、に分解するように構成され、割り当てのために、前記方向性信号のインデックスの前記遅延されたデータセット

[この文献は図面を表示できません]

および前記低減された数の残差のアンビエントＨＯＡ係数列の前記データセット

[この文献は図面を表示できません]

が使用される、該手段（１４、１５）と、
−前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））および前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））のＨＯＡ係数列を前記所定数（Ｉ）に対応する数のチャンネルに割り当てるように構成された手段（１６）であって、該割り当てにより前記割り当てを記述する選択されたアンビエントＨＯＡ係数列のインデックスのパラメータ

[この文献は図面を表示できません]

を取得し、該パラメータが圧縮解除側での対応する再配分に使用可能である、前記手段（１６）と、
−関連するフレーム（Ｙ（ｋ−２））の前記チャンネルを知覚符号化するように構成された手段（１７）であって、符号化された圧縮されたフレーム

[この文献は図面を表示できません]

が得られる、該手段（１７）と、
を備える、前記装置。

【請求項3】

前記非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））は、知覚に関連する基準に従って、
−前記圧縮のための所定の所与の数のチャンネルを用いることにより、対応して圧縮解除されるＨＯＡ表現が最小の知覚可能な誤差をもたらし、前記基準で考慮される誤差は、
−−複数の異なる数の前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））と前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））に対して複数の異なる数のＨＯＡ係数列とを
使用することから生ずるモデル化誤差と、
−−前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））の知覚符号化に伴う量子化雑音と、
−−前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））の個々のＨＯＡ
係数列を符号化することに伴う量子化雑音と、であり、
−前記３つの誤差の結果として生ずる合計誤差は、当該合計誤差の知覚可能性に関して、複数のテスト方向および複数の臨界帯域に対して考慮され、
−前記最小の知覚可能な誤差を達成するために、平均の知覚可能な誤差または最大の知覚可能な誤差を最小化するように前記非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））が選択される、ように決定される、請求項１に記載の方法。

【請求項4】

前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））を表現する低減された数のＨＯＡ係数列の選択は、３つの場合を区別する基準に従って行われ、前記３つの場合は、
−前記現在のフレーム（ｋ）のＨＯＡ係数列の数がその前のフレーム（ｋ−１）のＨＯＡ係数列の数と同一であるときに、前記前のフレームと同一のＨＯＡ係数列が選択される場合と、
−前記現在のフレーム（ｋ）のＨＯＡ係数列の数がその前のフレーム（ｋ−１）のＨＯＡ係数列の数よりも小さいときに、方向性信号によって占められている前記現在のフレーム内のチャンネルに割り当てられた前記前のフレームに存在した、前記前のフレームからのＨＯＡ係数列が非アクティブ化される場合と、
−前記現在のフレーム（ｋ）のＨＯＡ係数列の数がその前のフレーム（ｋ−１）のＨＯＡ係数列の数よりも大きいときに、前記前のフレームで選択されていたＨＯＡ係数列が前記現在のフレームでも選択され、追加的なＨＯＡ係数列が当該追加的なＨＯＡ係数列の知覚的な重要性または最も高い平均パワーに従って選択可能である場合と、
である、請求項１または３に記載の方法。

【請求項5】

前記割り当て（１６）の実行で、
−アクティブな方向性信号は、前記知覚符号化（１７）のための連続的な信号を取得するために、そのアクティブな方向性信号がチャンネル・インデックスを保持するように、所与のチャンネルに割り当てられ、
−前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））のＨＯＡ係数列は、最小の数（Ο_ＲＥＤ）の当該ＨＯＡ係数列が常に対応する数（Ο_ＲＥＤ）の前のチャンネルが含まれるように、割り当てられ、
−前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））の追加的なＨＯＡ係数列を割り当てるために、該追加的なＨＯＡ係数列が前記前のフレーム（ｋ−１）でも選択されていたかが判定され、
−−前記追加的なＨＯＡ係数列が前記前のフレーム（ｋ−１）でも選択されていた場合には、前記ＨＯＡ係数列の知覚符号化（１７）されるチャンネルに対する割り当て（１６）が前記前のフレームに対する割り当てと同一となり、
−−前記追加的なＨＯＡ係数列が前記前のフレーム（ｋ−１）では選択されておらず、
ＨＯＡ係数列が新たに選択されている場合には、前記ＨＯＡ係数列は、まずインデックスに関して昇順に配列され、当該昇順で、方向性信号によってまだ占められていない知覚符号化（１７）されるチャンネルに割り当てられる、
請求項１、３、および４のいずれか１項に記載の方法。

【請求項6】

Ο_ＲＥＤは、前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ,ＲＥＤ}（ｋ−２））を表現するＨＯＡ係数列の数であり、前記割り当て（１６）を記述するパラメータは、前記残差のアンビエントＨＯＡ成分を表現するΟ_ＲＥＤ個のＨＯＡ係数列に追加して使用される追加の数のＨＯＡ係数列に対応する長さを有するビットアレイに配列され、前記ビットアレイ内で各ｏ番目のビットは、（Ο_ＲＥＤ＋ｏ）番目の追加のＨＯＡ係数列が前記残差のアンビエントＨＯＡ成分を表現するために使用されているかどうかを示す、請求項１および３〜５のいずれか１項に記載の方法。

【請求項7】

前記割り当て（１６）を記述するパラメータは、非アクティブな方向性信号の数に対応する長さを有する割り当てベクトルに配列され、該割り当てベクトルの要素は、前記残差のアンビエントＨＯＡ成分の追加的なＨＯＡ係数列のうちのいずれが非アクティブな方向性信号を有するチャンネルに割り当てられるかを示す、請求項１および３〜５のいずれか１項に記載の方法。

【請求項8】

前記現在のフレームのＨＯＡ係数列の分解（１４）は、前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））から元のＨＯＡ表現の部分を予測するために圧縮解除側で使用可能なパラメータ（ζ（ｋ−２））をさらに供給する、請求項１および３〜７のいずれか１項に記載の方法。

【請求項9】

前記割り当て（１６）は、割り当てベクトル（γ（ｋ））を供給し、該割り当てベクトルの要素は、前記残差のアンビエントＨＯＡ成分に対する追加的なＨＯＡ係数列のうちのいずれが非アクティブな方向性信号を有するチャンネルに割り当てられるかについての情報を表現する、請求項５〜８のいずれか１項に記載の方法。

【請求項10】

請求項１に記載の方法に従って圧縮された高次アンビソニックス表現を圧縮解除する方法であって、
−チャンネルの知覚復号されたフレーム

[この文献は図面を表示できません]

を得るために、現在の符号化圧縮されたフレーム

[この文献は図面を表示できません]

を復号するステップ（３１）と、
−方向性信号のインデックスの前記データセット

[この文献は図面を表示できません]

と選択されたアンビエントＨＯＡ係数列のインデックスの前記データセット

[この文献は図面を表示できません]

を使用して、方向性信号

[この文献は図面を表示できません]

の前記対応するフレームと前記残差のアンビエントＨＯＡ成分

[この文献は図面を表示できません]

の前記対応するフレームとを再形成するために、チャンネルの前記知覚復号されたフレーム

[この文献は図面を表示できません]

を再配分するステップ（３２）と、
−検出された方向性信号のインデックスの前記データセット

[この文献は図面を表示できません]

および支配的な方向推定値の前記セット

[この文献は図面を表示できません]

を使用して、方向性信号

[この文献は図面を表示できません]

の前記フレームと前記残差のアンビエントＨＯＡ成分

[この文献は図面を表示できません]

の前記フレームとからＨＯＡ表現の現在の圧縮解除されたフレーム

[この文献は図面を表示できません]

を再合成するステップ（３３）と、を含み、
均一に分布した方向に対する方向性信号が前記方向性信号

[この文献は図面を表示できません]

から予測され、その後に、前記現在の圧縮解除されたフレーム

[この文献は図面を表示できません]

は、方向性信号

[この文献は図面を表示できません]

の前記フレーム、前記予測された信号、および前記残差のアンビエントＨＯＡ成分

[この文献は図面を表示できません]

から再合成される、前記方法。

【請求項11】

請求項１に記載の方法に従って圧縮された高次アンビソニックス表現を圧縮解除する装置であって、
−チャンネルの知覚復号されたフレーム

[この文献は図面を表示できません]

を得るために、現在の符号化圧縮されたフレーム

[この文献は図面を表示できません]

を復号するように構成された手段（３１）と、
−検出された方向性信号のインデックスの前記データセット

[この文献は図面を表示できません]

と選択されたアンビエントＨＯＡ係数列のインデックスの前記データセット

[この文献は図面を表示できません]

を使用して、方向性信号

[この文献は図面を表示できません]

の前記対応するフレームと前記残差のアンビエントＨＯＡ成分

[この文献は図面を表示できません]

の前記対応するフレームとを再形成するために、チャンネルの前記知覚復号されたフレーム

[この文献は図面を表示できません]

を再配分するように構成された手段（３２）と、
−検出された方向性信号のインデックスの前記データセット

[この文献は図面を表示できません]

および支配的な方向推定値の前記セット

[この文献は図面を表示できません]

を使用して、方向性信号

[この文献は図面を表示できません]

の前記フレームと前記残差のアンビエントＨＯＡ成分

[この文献は図面を表示できません]

の前記フレームとからＨＯＡ表現の現在の圧縮解除されたフレーム

[この文献は図面を表示できません]

を再合成するように構成された手段（３３）と、を備え、
均一に分布した方向に対する方向性信号が前記方向性信号

[この文献は図面を表示できません]

から予測され、その後に、前記現在の圧縮解除されたフレーム

[この文献は図面を表示できません]

は、方向性信号

[この文献は図面を表示できません]

の前記フレーム、前記予測された信号、および前記残差のアンビエントＨＯＡ成分

[この文献は図面を表示できません]

から再合成される、前記装置。

【請求項12】

均一に分布した方向に対する方向性信号の前記予測は、前記方向性信号

[この文献は図面を表示できません]

から、前記予測のための前記供給されたパラメータ（ζ（ｋ−２））を使用して行われる、
請求項８に記載の方法。

【請求項13】

前記再配分（３２）において、検出された方向性信号のインデックスの前記データセット

[この文献は図面を表示できません]

および選択されたアンビエントＨＯＡ係数列のインデックスの前記データセット

[この文献は図面を表示できません]

の代わりに、受信した割り当てベクトル（γ（ｋ））が使用され、該割り当てベクトルの要素は、前記残差のアンビエントＨＯＡ成分の追加的なＨＯＡ係数列のうちのいずれが非アクティブな方向性信号を有するチャンネルに割り当てられるかについての情報を表現する、請求項１０または１２に記載の方法。

【請求項14】

前記非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））は、知覚に関連する基準に従って、
−前記圧縮のための所定の所与の数のチャンネルを用いることにより、対応して圧縮解除されるＨＯＡ表現が最小の知覚可能な誤差をもたらし、前記基準で考慮される誤差は、
−−複数の異なる数の前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））と前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））に対して複数の異なる数のＨＯＡ係数列とを
使用することから生ずるモデル化誤差と、
−−前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））の知覚符号化に伴う量子化雑音と、
−−前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））の個々のＨＯＡ
係数列を符号化することに伴う量子化雑音と、であり、
−前記３つの誤差の結果として生ずる合計誤差は、当該合計誤差の知覚可能性に関して、複数のテスト方向および複数の臨界帯域に対して考慮され、
−前記最小の知覚可能な誤差を達成するために、平均の知覚可能な誤差または最大の知覚可能な誤差を最小化するように前記非所定数（Ｍ）の方向性信号（Ｘ_ＤＩＲ（ｋ−２））が選択される、ように決定される、請求項２に記載の装置。

【請求項15】

前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））を表現する低減された数のＨＯＡ係数列の選択は、３つの場合を区別する基準に従って行われ、前記３つの場合は、
−前記現在のフレーム（ｋ）のＨＯＡ係数列の数がその前のフレーム（ｋ−１）のＨＯＡ係数列の数と同一であるときに、前記前のフレームと同一のＨＯＡ係数列が選択される場合と、
−前記現在のフレーム（ｋ）のＨＯＡ係数列の数がその前のフレーム（ｋ−１）のＨＯＡ係数列の数よりも小さいときに、方向性信号によって占められている前記現在のフレーム内のチャンネルに割り当てられた前記前のフレームに存在した、前記前のフレームからのＨＯＡ係数列が非アクティブ化される場合と、
−前記現在のフレーム（ｋ）のＨＯＡ係数列の数がその前のフレーム（ｋ−１）のＨＯＡ係数列の数よりも大きいときに、前記前のフレームで選択されていたＨＯＡ係数列が前記現在のフレームでも選択され、追加的なＨＯＡ係数列が当該追加的なＨＯＡ係数列の知覚的な重要性または最も高い平均パワーに従って選択可能である場合と、
である、請求項２または１４に記載の装置。

【請求項16】

前記割り当て（１６）の実行で、
−アクティブな方向性信号は、前記知覚符号化（１７）のための連続的な信号を取得するために、そのアクティブな方向性信号がチャンネル・インデックスを保持するように、所与のチャンネルに割り当てられ、
−前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））のＨＯＡ係数列は、最小の数（Ο_ＲＥＤ）の当該ＨＯＡ係数列が常に対応する数（Ο_ＲＥＤ）の前のチャンネルが含まれるように、割り当てられ、
−前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２））の追加的なＨＯＡ係数列を割り当てるために、該追加的なＨＯＡ係数列が前記前のフレーム（ｋ−１）でも選択されていたかが判定され、
−−前記追加的なＨＯＡ係数列が前記前のフレーム（ｋ−１）でも選択されていた場合には、前記ＨＯＡ係数列の知覚符号化（１７）されるチャンネルに対する割り当て（１６）が前記前のフレームに対する割り当てと同一となり、
−−前記追加的なＨＯＡ係数列が前記前のフレーム（ｋ−１）では選択されておらず、
ＨＯＡ係数列が新たに選択されている場合には、前記ＨＯＡ係数列は、まずインデックスに関して昇順に配列され、当該昇順で、方向性信号によってまだ占められていない知覚符号化（１７）されるチャンネルに割り当てられる、
請求項２、１４および１５のいずれか１項に記載の装置。

【請求項17】

Ο_ＲＥＤは、前記残差のアンビエントＨＯＡ成分（Ｃ_{ＡＭＢ,ＲＥＤ}（ｋ−２））を表現するＨＯＡ係数列の数であり、前記割り当て（１６）を記述するパラメータは、前記残差のアンビエントＨＯＡ成分を表現するΟ_ＲＥＤ個のＨＯＡ係数列に追加して使用される追加の数のＨＯＡ係数列に対応する長さを有するビットアレイに配列され、前記ビットアレイ内で各ｏ番目のビットは、（Ο_ＲＥＤ＋ｏ）番目の追加のＨＯＡ係数列が前記残差のアンビエントＨＯＡ成分を表現するために使用されているかどうかを示す、請求項２および１４〜１６のいずれか１項に記載の装置。

【請求項18】

前記割り当て（１６）を記述するパラメータは、非アクティブな方向性信号の数に対応する長さを有する割り当てベクトルに配列され、該割り当てベクトルの要素は、前記残差のアンビエントＨＯＡ成分の追加的なＨＯＡ係数列のうちのいずれが非アクティブな方向性信号を有するチャンネルに割り当てられるかを示す、請求項２および１４〜１６のいずれか１項に記載の装置。

【請求項19】

前記現在のフレームのＨＯＡ係数列の分解（１４）は、前記方向性信号（Ｘ_ＤＩＲ（ｋ−２））から元のＨＯＡ表現の部分を予測するために圧縮解除側で使用可能なパラメータ（ζ（ｋ−２））をさらに供給する、請求項２及び１４〜１８のいずれか１項に記載の装置。

【請求項20】

前記割り当て（１６）は、割り当てベクトル（γ（ｋ））を供給し、該割り当てベクトルの要素は、前記残差のアンビエントＨＯＡ成分に対する追加的なＨＯＡ係数列のうちのいずれが非アクティブな方向性信号を有するチャンネルに割り当てられるかについての情報を表現する、請求項１６〜１９のいずれか１項に記載の装置。

【請求項21】

均一に分布した方向に対する方向性信号の前記予測は、前記方向性信号

[この文献は図面を表示できません]

から、前記予測のための前記供給されたパラメータ（ζ（ｋ−２））を使用して行われる、
請求項１９に記載の装置。

【請求項22】

前記再配分（３２）において、検出された方向性信号のインデックスの前記データセット

[この文献は図面を表示できません]

および選択されたアンビエントＨＯＡ係数列のインデックスの前記データセット

[この文献は図面を表示できません]

の代わりに、受信した割り当てベクトル（γ（ｋ））が使用され、該割り当てベクトルの要素は、前記残差のアンビエントＨＯＡ成分の追加的なＨＯＡ係数列のうちのいずれが非アクティブな方向性信号を有するチャンネルに割り当てられるかについての情報を表現する、請求項１１または２１に記載の装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、方向性信号成分およびアンビエント信号成分を別々に処理することによって高次アンビソニックス表現を圧縮および圧縮解除する方法および装置に関する。

【背景技術】

【0002】

高次アンビソニックス（ＨＯＡ）は、波面合成法（ＷＦＳ）や２２．２のようなチャンネルに基づくアプローチといった他の技術が存在する一方で、三次元音声を表現する１つの可能性を提供している。チャンネルに基づく方法と対照的に、ＨＯＡ表現には、特定のラウドスピーカの設定とは独立しているという利点がある。しかしながら、この柔軟性を得るためには、特定のラウドスピーカの設定でＨＯＡ表現を再生するための復号処理が必要となる。通常、必要なラウドスピーカの数が大変多くなるＷＦＳのアプローチと比較して、ＨＯＡは極めて少ない数のラウドスピーカのみで構成される設定にすることができる。ＨＯＡのさらなる利点は、ヘッドフォンへのバイノーラル・レンダリングにも変更を必要とすることなく同じ表現を利用できる点にある。

【0003】

ＨＯＡは、切断球面調和関数（ＳＨ）展開による複素調和平面波振幅の空間密度の表現に基づいている。各展開係数は角周波数の関数であり、これを、時間領域関数によって同等に表現することができる。したがって、一般性を失うことなく、完全なＨＯＡ音場表現は、実際には、“Ο”個の時間領域関数から構成されるものと考えることができる。ここで、Οは、展開係数の数を表している。これらの時間領域関数は同等の意味を有するものとして以下のＨＯＡ係数列またはＨＯＡチャンネルを参照する。

【0004】

ＨＯＡ表現の空間解像度は、展開の最大次数Ｎの増加とともに向上する。残念ながら、展開係数の数“Ο”は、次数Ｎに対して二乗的に増加し、特にΟ＝（Ｎ＋１）²となる。例えば、次数Ｎ＝４を使用した一般的なＨＯＡ表現には、Ο＝２５の個数のＨＯＡ（展開）係数が必要となる。上記の点を考慮して、ＨＯＡ表現の伝送のための合計ビットレートは、所望の単一チャンネルのサンプリング・レートｆ_ｓおよびサンプル毎のビットの数Ｎ_ｂが与えられると、Ο・ｆ_ｓ・Ｎ_ｂによって求めることができる。したがって、サンプル毎にＮ_ｂ＝１６の個数のビットを使用してｆ_ｓ＝４８ｋＨｚのサンプリング・レートでの次数Ｎ＝４のＨＯＡ表現を伝送すると、結果として、ビットレートは、１９．２メガビット／秒となるが、これは、多くの実用的なアプリケーション、例えば、ストリーミングでは極めて高いビットレートである。

【0005】

ＨＯＡ音場表現の圧縮は、欧州特許出願第１２３０６５６９号および欧州特許出願第１２３０５５３７号において提案されている。例えば、Ｅ．Ｈｅｌｌｅｒｕｄ、Ｉ．Ｂｕｒｎｅｔｔ、Ａ．ＳｏｌｖａｎｇおよびＵ．Ｐ．Ｓｖｅｎｓｓｏｎの「ＡＡＣを用いた高次アンビソニックスの符号化」１２４回ＡＥＳコンベンション、アムステルダム、２００８年、において行われているような、ＨＯＡ係数列を個々に知覚符号化することの代わりに、特に音場分析を行い、所与のＨＯＡ表現を方向性成分および残差アンビエント成分に分解することによって、知覚符号化される信号の数を減少させる試みが行われている。一般的には、方向性成分は、一般的な平面波関数とみなすことができる少数の支配的な方向性信号によって表現されるものとされる。残差のアンビエントＨＯＡ成分の次数が低減される。その理由は、支配的な方向性信号を抽出した後には、より低次のＨＯＡ係数が最も関連する情報を保持していると考えられるからである。

【発明の概要】

【0006】

総括すると、そのような処理を行うことによって、知覚符号化されるＨＯＡ係数列の初期数（Ｎ＋１）^２は、Ｄ個の支配的な方向性信号の所定数と、切断次数Ｎ_ＲＥＤ＜Ｎを用いて残差のアンビエントＨＯＡ成分を表現する（Ｎ_ＲＥＤ＋１）^２個のＨＯＡ係数列の数とに低減される。それによって、符号化される信号の数が決まり、すなわち、Ｄ＋（Ｎ_ＲＥＤ＋１）^２となる。特に、この数は、時間フレームｋにおけるアクティブな支配的な方向性音源の実際に検出された数Ｄ_ＡＣＴ（ｋ）≦Ｄとは独立している。これは、時間フレームｋにおいて、アクティブな支配的な方向性音源の実際に検出された数Ｄ_ＡＣＴ（ｋ）が方向性信号の最大許容数Ｄよりも小さい場合、知覚符号化される支配的な方向性信号のいくつかまたは全てさえもが零となることを意味している。つまり、これはこの複数のチャンネルが音場の関連情報を捕捉するために全く使用されないことを意味する。

【0007】

この状況で、欧州特許出願第１２３０６５６９号および欧州特許出願第１２３０５５３７号における処理の別の想定される弱点は、各時間フレーム内の支配的な方向性信号の数を決定するための基準である。その理由は、音場の連続的な知覚符号化に関してアクティブな支配的な方向性信号の最適な数を決定する試みが行われていないからである。例えば、欧州特許出願第１２３０５５３７号においては、支配的な音源の数が単純なパワー基準を使用して、すなわち、最大の固有値に属する係数間の相関行列の部分空間の次元を求めることによって推定される。欧州特許出願第１２３０６５６９号においては、支配的な方向性音源のインクリメンタル検出が提案されている。ここで、各々の方向からの平面波関数のパワーが最初の方向性信号に対して十分に高い場合には、方向性音源が支配的であると考慮される。欧州特許出願第１２３０６５６９号および欧州特許出願第１２３０５５３７号の場合のようなパワーに基づく基準を使用すると、音場の知覚符号化に関して最適であるとは云えない方向性−アンビエント分解となることもある。

【0008】

本発明によって解決される課題は、現在のＨＯＡオーディオ信号コンテンツに対して、所定の低減された数のチャンネルに、方向性信号およびアンビエントＨＯＡ成分に対する係数をどのように割り当てるかを決定することによって、ＨＯＡ圧縮を改善することにある。この課題は、請求項１および３に開示されたそれぞれの方法によって解決される。これらの方法を利用する装置は、請求項２および４において開示されている。

【0009】

本発明は、２つの態様において、欧州特許出願第１２３０６５６９号で提案されている圧縮処理を改善する。第１に、知覚符号化される所与の数のチャンネルによってもたらされる帯域幅が良好に利用される。支配的な音源信号が検出されない時間フレームでは、支配的な方向性信号に対して当初より確保されているチャンネルは、アンビエント成分についての追加的な情報を捕捉するために、残差のアンビエントＨＯＡ成分の追加的なＨＯＡ係数列の形式で使用される。第２に、所与のＨＯＡ音場表現を知覚符号化するために所与の数のチャンネルを利用するという目的を念頭に置くと、ＨＯＡ表現から抽出される方向性信号の数を決定するための基準は、その目的に対して適応化される。方向性信号の数は、復号され再構築されたＨＯＡ表現によって知覚される誤差が最も小さくなるように決定される。その基準は、方向性信号を抽出することと残差のアンビエントＨＯＡ成分を記述するためにＨＯＡ係数列をより少なく使用することとから生ずるモデル化誤差と、方向性信号を抽出することなく、その代わりに残差のアンビエントＨＯＡ成分を記述するために追加的なＨＯＡ係数列を使用することから生ずるモデル化誤差とを比較する。その基準は、さらに、その双方の場合に対して、方向性信号および残差のアンビエントＨＯＡ成分のＨＯＡ係数列の知覚符号化によってもたらされる量子化雑音の空間パワー分布を考慮する。

【0010】

上述した処理を実施するために、ＨＯＡ圧縮を開始する前に、信号（チャンネル）の合計数Ｉが定められる。この合計数Ｉは、当初のΟ個のＨＯＡ係数列の数と比較して低減させられたものである。アンビエントＨＯＡ成分は、最小の数Ο_ＲＥＤ個のＨＯＡ係数列によって表現されるものと仮定される。場合によっては、その最小の数が零となることもある。残りのＤ＝Ｉ−Ο_ＲＥＤ個のチャンネルは、方向性信号抽出処理が判定する知覚的に意味のよりあるものに依存して、方向性信号またはアンビエントＨＯＡ成分の追加的な係数列のいずれかを含むものとされる。方向性信号またはアンビエントＨＯＡ成分係数列のいずれかの残りのＤ個のチャンネルに対する割り当ては、フレーム単位で変更可能であるものと仮定される。受信機側での音場の再構築のために、この割り当てについての情報は、追加の副情報として送信される。

【0011】

原理的には、本発明の圧縮方法は、所定数の知覚符号化処理を使用して、ＨＯＡと称する音場の高次アンビソニックス表現をＨＯＡ係数列の入力される時間フレームを用いて圧縮するのに適している。この方法は、フレーム単位で行われ、
−現在のフレームに対して、支配的な方向のセットおよび対応する検出された方向性信号のインデックスのデータセットを推定するステップと、
−上記現在のフレームのＨＯＡ係数列を分解するステップであって、非所定数の方向性信号であって、支配的な方向推定値の上記セットに含まれる各々の方向と上記方向性信号のインデックスの各々のデータセットとを用いた、上記非所定数が上記所定数よりも小さい、上記非所定数の方向性信号と、上記所定数と上記非所定数との差に対応する低減された数のＨＯＡ係数列によって表現される残差のアンビエントＨＯＡ成分と、対応する上記低減された数の残差のアンビエントＨＯＡ係数列のインデックスのデータセットと、に分解する、上記分解するステップと、
−上記方向性信号および上記残差のアンビエントＨＯＡ成分のＨＯＡ係数列を上記所定数に対応する数のチャンネルに割り当てるステップであって、上記割り当てのために、上記方向性信号のインデックスの上記データセットおよび上記低減された数の残差のアンビエントＨＯＡ係数列のインデックスの上記データセットが使用される、上記割り当てるステップと、
−関連するフレームの上記チャンネルを知覚符号化するステップであって、符号化された圧縮されたフレームが得られる、上記知覚符号化するステップと、を含む。

【0012】

原理的には、本発明の圧縮装置は、所定数の知覚符号化処理を使用して音場のＨＯＡと称する高次アンビソニックス表現をＨＯＡ係数列の入力される時間フレームを用いて圧縮するのに適している。
上記装置は、フレーム単位の処理を実行し、
−現在のフレームに対して、支配的な方向のセットおよび対応する検出された方向性信号のインデックスのデータセットを推定するように構成された手段と、
−上記現在のフレームのＨＯＡ係数列を分解するように構成された手段であって、非所定数の方向性信号であって、支配的な方向推定値の上記セットに含まれる各々の方向と、上記方向性信号のインデックスの各々のデータセットとを用いた、上記非所定数が上記所定数よりも小さい、上記非所定数の方向性信号と、上記所定数と上記非所定数との差に対応する低減された数のＨＯＡ係数列によって表現される残差のアンビエントＨＯＡ成分と、対応する上記低減された数の残差のアンビエントＨＯＡ係数列のインデックスの対応するデータセットと、に分解するように構成された、上記手段と、
−上記方向性信号および上記残差のアンビエントＨＯＡ成分のＨＯＡ係数列を上記所定数に対応する数のチャンネルに割り当てるように構成された手段であって、上記割り当てのために、上記方向性信号のインデックスの上記データセットおよび上記低減された数の残差のアンビエントＨＯＡ係数列のインデックスの上記データセットが使用される、上記手段と、
−関連するフレームの上記チャンネルを知覚符号化するように構成された手段であって、符号化された圧縮されたフレームが得られる、上記手段と、を含む。

【0013】

原理的には、本発明の圧縮解除方法は、上述の圧縮方法に従って圧縮された高次アンビソニックス表現を圧縮解除するのに適している。この圧縮解除方法は、
−チャンネルの知覚復号されたフレームを得るために、現在の符号化圧縮されたフレームを復号するステップと、
−検出された方向性信号のインデックスの上記データセットと上記選択されたアンビエントＨＯＡ係数列のインデックスの上記データセットを使用して、方向性信号の上記対応するフレームと残差のアンビエントＨＯＡ成分の上記対応するフレームとを再形成するために、チャンネルの上記知覚復号されたフレームを再配分するステップと、
−検出された方向性信号のインデックスの上記データセットおよび支配的な方向性推定値の上記セットを使用して、方向性信号の上記フレームと上記残差のアンビエントＨＯＡ成分の上記フレームとからＨＯＡ表現の現在の圧縮解除されたフレームを再合成するステップと、を含み、
均一に分布した方向に対する方向性信号が上記方向性信号から予測され、その後に、上記現在の圧縮解除されたフレームは、方向性信号の上記フレーム、上記予測された信号、および上記残差のアンビエントＨＯＡ成分から再合成される。

【0014】

原理的には、本発明の圧縮解除装置は、上述の圧縮方法に従って圧縮された高次アンビソニックス表現を圧縮解除するのに適している。この装置は、
−チャンネルの知覚復号されたフレームを得るために、現在の符号化圧縮されたフレームを復号するように構成された手段と、
−検出された方向性信号のインデックスの上記データセットと選択されたアンビエントＨＯＡ係数列のインデックスの上記データセットを使用して、方向性信号の上記対応するフレームと上記残差のアンビエントＨＯＡ成分の上記対応するフレームとを再形成するために、チャンネルの上記知覚復号されたフレームを再配分するように構成された手段と、
−検出された方向性信号のインデックスの上記データセットおよび支配的な方向性推定値の上記セットを使用して、方向性信号の上記フレームと上記残差のアンビエントＨＯＡ成分の上記フレームとから、上記ＨＯＡ表現の現在の圧縮解除されたフレームを再合成するように構成された手段と、を含み、
均一に分布した方向に対する方向性信号が上記方向性信号から予測され、その後に、上記現在の圧縮解除されたフレームは、方向性信号の上記フレーム、上記予測された信号、および上記残差のアンビエントＨＯＡ成分から再合成される。

【0015】

本発明の追加的な実施形態は、各々の従属請求項に開示されており、有利なものである。

【図面の簡単な説明】

【0016】

【図1】ＨＯＡ圧縮のブロック図である。

【図2】支配的な音源方向の推定のブロック図である。

【図3】ＨＯＡ圧縮解除のブロック図である。

【図4】球面座標システムを示す図である。

【図5】複数の異なるアンビソニックス次数Ｎおよび角度θ∈[0,π]に対する正規化された分散関数ν_N(Θ)を示す図である。

【発明を実施するための形態】

【0017】

本発明の例示的な実施形態は、添付図面を参照して説明される。
Ａ．改良されたＨＯＡ圧縮
本発明に係る圧縮処理は、欧州特許出願第１２３０６５６９号に基づいており、図１に示されている。ここで、信号処理ブロックは、欧州特許出願第１２３０６５６９号に対して変更が加えられ、または新たに導入されており、その信号処理ブロックは太字のボックスで示されており、本出願における「

[この文献は図面を表示できません]

」（方向推定値とされたもの）および「Ｃ」は、それぞれ、欧州特許出願第１２３０６５６９号の「Ａ」（方向推定値の行列）および「Ｄ」に対応する。

【0018】

ＨＯＡ圧縮のために、長さＬのＨＯＡ係数列の重複しない入力フレームＣ（ｋ）を用いたフレーム単位の処理が使用される。ここで、ｋは、フレームのインデックスを表す。フレームは、下記の式（１）に特定されたＨＯＡ係数列に関して定義される。

【数1】

[この文献は図面を表示できません]

ここで、Ｔ_ｓは、サンプリング期間を表す。

【0019】

図１のステップまたはステージ１１／１２は、任意に行われ、ＨＯＡ係数列の重複しないｋ番目のフレームおよび（ｋ−１）番目のフレームを下記の式に従って連結して長いフレーム

[この文献は図面を表示できません]

にすることを含む。

【数2】

[この文献は図面を表示できません]

この長いフレームは、隣接する長いフレームと５０％重複し、長いフレームは、支配的な音源方向の推定に連続的に使用される。

[この文献は図面を表示できません]

の表記と同様に、チルダ記号は、以下の説明において、各々の量が長い重複するフレームを指すことを示すために使用される。ステップ／ステージ１１／１２が存在しない場合には、チルダ記号は特別な意味を持たない。

【0020】

原理的には、支配的な音源の推定ステップまたはステージ１３は、欧州特許出願第１３３０５１５６号に提案されているように行われるが、重要な変更を有する。この変更は、検出される方向の数の決定、すなわち、何個の方向性信号がＨＯＡ表現から抽出されるとするかに関する。これは、アンビエントＨＯＡ成分の良好な近似計算のために、追加的なＨＯＡ係数列を使用することよりも方向性信号を抽出することの方が知覚的に関連性が高い場合にのみ、追加的なＨＯＡ係数列を使用する代わりに方向性信号を抽出しようとする考えから成し遂げられるものである。Ａ．２の項目でこの技術についての詳細な説明を行う。

【0021】

支配的な音源の推定により、検出された方向性信号のインデックスのデータセット

[この文献は図面を表示できません]

と、対応する方向推定値のセット

[この文献は図面を表示できません]

とが得られる。Ｄは、ＨＯＡ圧縮を開始する前に設定しなければならない方向性信号の最大数を示している。

【0022】

ステップまたはステージ１４において、ＨＯＡ係数列の現在の（長い）フレーム

[この文献は図面を表示できません]

が、セット

[この文献は図面を表示できません]

内に含まれる方向に属する複数の方向性信号Ｘ_ＤＩＲ（ｋ−２）と、残差のアンビエントＨＯＡ成分Ｃ_ＡＭＢ（ｋ−２）とに分解される（欧州特許出願第１３３０５１５６号に提案されているように）。滑らかな信号を得るために、重畳加算処理の結果として２つのフレーム分の遅延が導入される。Ｘ_ＤＩＲ（ｋ−２）は、合計Ｄ個のチャンネルを含むものの、このうち、アクティブな方向性信号に対応するチャンネルのみが零でないと仮定される。このチャンネルを特定するインデックスは、データセット

[この文献は図面を表示できません]

内において出力されるものと仮定される。さらに、ステップ／ステージ１４における分解によって、方向性信号から元のＨＯＡ表現の部分を予測するために圧縮解除側で使用されるいくつかのパラメータζ（ｋ−２）を供給する（より詳細には欧州特許出願第１３３０５１５６号参照）。ステップまたはステージ１５において、アンビエントＨＯＡ成分Ｃ_ＡＭＢ（ｋ−２）の係数の数はインテリジェントに低減され、Ο_ＲＥＤ＋Ｄ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個の非零のＨＯＡ係数列のみを含むようになる。ここで、

[この文献は図面を表示できません]

は、データセット

[この文献は図面を表示できません]

の組の数、すなわち、フレームｋ−２内のアクティブな方向性信号の数を示す。アンビエントＨＯＡ成分は、最小の数Ο_ＲＥＤ個のＨＯＡ係数列によって常に表現されると仮定されるため、この問題は、実際には、想定されるΟ−Ο_ＲＥＤ個のＨＯＡ係数列から残りのＤ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個のＨＯＡ係数列を選択することに集約される。滑らかな低減されたアンビエントＨＯＡ表現を取得するために、この選択は、前のフレームｋ−３で行った選択と比較して、変更が可能な限り少なくなるように行われる。

【0023】

特に、以下の３つの場合を区別すべきである。

【0024】

ａ）Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）＝Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−３）：この場合、フレームｋ−３の場合と同様に、同一のＨＯＡ係数列が選択されるものと想定される。

【0025】

ｂ）Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）＜Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−３）：この場合、現在のフレーム内のアンビエントＨＯＡ成分を表現するために、この前のフレームｋ−３よりも多いＨＯＡ係数列を使用することができる。ｋ−３において選択済のそのＨＯＡ係数列は、現在のフレーム内でも選択されるものと仮定される。異なる基準に従って追加的なＨＯＡ係数列を選択可能である。例えば、最高の平均パワーを有するＨＯＡ係数列をＣ_ＡＭＢ（ｋ−２）内で選択するか、あるいは、それぞれの知覚的な重要性に関してＨＯＡ係数列を選択する。

【0026】

ｃ）Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）＞Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−３）：この場合、現在のフレーム内のアンビエントＨＯＡ成分を表現するために、最後のフレームｋ−３に存在するＨＯＡ係数列よりも少ないＨＯＡ係数列を使用することができる。ここで解決すべき課題は、既に選択済のＨＯＡ係数列のうち、どれを非アクティブ化しなければならないかである。合理的な解決法は、フレームｋ−３で、信号を割り当てるステップまたはステージ１６でチャンネル

[この文献は図面を表示できません]

に割り当てられたＨＯＡ係数列を非アクティブ化することである。

【0027】

追加的なＨＯＡ係数列がアクティブ化または非アクティブ化されるときのフレーム境界での不連続を回避するために、各々の信号を平滑的にフェード・インまたはフェード・アウトさせるとよい。

【0028】

Ο_ＲＥＤ＋Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個の低減された個数の最終的なアンビエントＨＯＡ表現は、Ｃ_{ＡＭＢ，ＲＥＤ}（ｋ−２）によって示される。選択されたアンビエント係数列のインデックスは、データセット

[この文献は図面を表示できません]

内に出力される。

【0029】

ステップ／ステージ１６において、Ｘ_ＤＩＲ（ｋ−２）に含まれるアクティブな方向性信号およびＣ_{ＡＭＢ，ＲＥＤ}（ｋ−２）に含まれるＨＯＡ係数列は、個々の知覚符号化のためにＩ個のチャンネルのフレームＹ（ｋ−２）に割り当てられる。より詳細に信号の割り当てを記述すると、フレームＸ_ＤＩＲ（ｋ−２）、Ｙ（ｋ−２）およびＣ_{ＡＭＤ，ＲＥＤ}（ｋ−２）は、下記のように、個々の信号ｘ_DIR，d（ｋ−２）（ｄ∈｛１，… ，Ｄ｝）、ｙ_i（ｋ−２）（ｉ∈｛１，… ，I｝）およびｃ_{AMB, RED, ο}（ｋ−２）（ο∈１，… ，Ο）によって構成されるものと仮定される。

【数3】

[この文献は図面を表示できません]

【0030】

連続する知覚符号化のために連続した信号を取得するために、それぞれのチャンネルのインデックスを保持するようにアクティブな方向性信号が割り当てられる。これを下記の式のように表すことができる。

【数4】

[この文献は図面を表示できません]

【0031】

アンビエント成分のＨＯＡ係数列は、最小の数のΟ_ＲＥＤ個の係数列がＹ（ｋ−２）の最後のΟ_ＲＥＤ個の信号に常に含まれるように、すなわち、下記の式に従って割り当てられる。

【数5】

[この文献は図面を表示できません]

【0032】

追加的なＤ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個のアンビエント成分のＨＯＡ係数列については、これらが前のフレームでも選択されていたかどうかを区別すべきである。
ａ）追加的なＤ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個のアンビエント成分のＨＯＡ係数列が送信されるものとして前のフレーム内でも選択されていた場合、すなわち、各々のインデックスもまた、データセット

[この文献は図面を表示できません]

に含まれる場合には、これらの係数列のＹ（ｋ−２）における信号への割り当ては、前のフレームに対する割り当てと同じである。この処理は、滑らかな信号ｙ_ｉ（ｋ−２）を確保するものであり、ステップまたはステージ１７における連続的な知覚符号化にとって好ましいものである。
ｂ）そうではなく、いくつかの係数列が新たに選択されている場合、すなわち、これらのインデックスがデータセット

[この文献は図面を表示できません]

に含まれているが、データセット

[この文献は図面を表示できません]

に含まれていない場合には、これらはまず、インデックスに関して昇順に配列され、この順番で方向性信号によってまだ占められていないＹ（ｋ−２）のチャンネル

[この文献は図面を表示できません]

に割り当てられる。

【0033】

この特定の割り当ては、ＨＯＡ圧縮解除処理の間に信号の再配分および合成が、どのアンビエントＨＯＡ係数列がＹ（ｋ−２）個のどのチャンネルに含まれているかについての情報無しに行えるようになるという利点を提供する。代わりに、データセット

[この文献は図面を表示できません]

および

[この文献は図面を表示できません]

の情報のみで、ＨＯＡ圧縮解除の間に割り当てを再構築することができる。

【0034】

この割り当て処理によって、割り当てベクトル

[この文献は図面を表示できません]

ももたらされることが有利である。この要素γ_ο（ｋ）（ο＝１，… ，Ｄ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２））は追加的なＤ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個のアンビエント成分のＨＯＡ係数列の各々のインデックスを表す。換言すれば、割り当てベクトルγ（ｋ）の要素により、追加的なΟ−Ο_ＲＥＤ個のアンビエントＨＯＡ成分のＨＯＡ係数列のうちのいずれがＤ−Ｎ_{ＤＩＲ，ＡＣＴ}（ｋ−２）個の非アクティブな方向性信号のチャンネルに割り当てられるかについての情報が得られる。このベクトルは、ＨＯＡ圧縮解除のために行われる再配分処理の初期化（項目Ｂ参照）を可能にするために、追加的に、フレームレートによる送信よりも低い頻度ではあるが送信されることがある。知覚符号化ステップ／ステージ１７は、フレームＹ（ｋ−２）のＩ個のチャンネルを符号化し、符号化されたフレーム

[この文献は図面を表示できません]

を出力する。

【0035】

ステップ／ステージ１６でベクトルγ（ｋ）が送信されないフレームについては、圧縮解除側で、データ・パラメータ・セット

[この文献は図面を表示できません]

および

[この文献は図面を表示できません]

がベクトルγ（ｋ）の代わりに再配分を行うために使用される。

【0036】

Ａ．１ 支配的な音源方向の推定
図１の支配的な音源方向に対する推定ステップ／ステージ１３が図２により詳細に描かれている。これは、本質的に、欧州特許出願第１３３０５１５６号に記載された内容に従って行われるが、決定的な違いがある。その決定的な違いは、支配的な音源の数を決定する手法である。支配的な音源の数は、所与のＨＯＡ表現から抽出される方向性信号の数に対応する。この数は重要であり、その理由は、より多くの方向性信号を使用すること、あるいはその代わりに、より多くのＨＯＡ係数列を使用してアンビエントＨＯＡ成分をより良好にモデル化することのいずれかによって、所与のＨＯＡ表現がより良好に表現されているかを制御するためにこの数が使用されるからである。

【0037】

支配的な音源方向の推定は、入力されるＨＯＡ係数列の長いフレーム

[この文献は図面を表示できません]

を使用して、支配的な音源方向の予備サーチで、ステップまたはステージ２１において開始する。予備的な方向推定値

[この文献は図面を表示できません]

と共に、個々の音源によって形成されるものとされる、予備的な方向推定値に対応する方向性信号

[この文献は図面を表示できません]

およびＨＯＡ音場成分

[この文献は図面を表示できません]

を欧州特許出願第１３３０５１５６号に記載された内容に従って算出する。

【0038】

ステップまたはステージ２２において、予備的な方向推定値、方向性信号、およびＨＯＡ音場成分は、抽出される方向性信号の数

[この文献は図面を表示できません]

を決定するために入力されるＨＯＡ係数列のフレーム

[この文献は図面を表示できません]

と共に使用される。結果として、

[この文献は図面を表示できません]

の方向性推定値

[この文献は図面を表示できません]

、これと対応する方向性信号

[この文献は図面を表示できません]

、およびＨＯＡ音場成分

[この文献は図面を表示できません]

が破棄される。代わりに、

[この文献は図面を表示できません]

の方向推定値

[この文献は図面を表示できません]

のみが、次に、既に見つかっている音源に対して割り当てられる。

【0039】

ステップまたはステージ２３において、結果として得られる方向軌跡は、音源動きモデルに従ってスムージング（滑らかに）され、音源のいずれがアクティブであるとされるかが決定される（欧州特許出願第１３３０５１５６号参照）。この最後の処理により、アクティブな方向性音源のインデックスのセット

[この文献は図面を表示できません]

とこれに対応する方向推定値のセット

[この文献は図面を表示できません]

とが得られる。

【0040】

Ａ．２ 抽出される方向性信号の数の決定
ステップ／ステージ２２において方向性信号の数を決定するために、知覚的に最も関連する音場情報を捕捉するために利用される所与の合計数のＩ個のチャンネルが存在する状況が想定される。したがって、全体としてのＨＯＡ圧縮／圧縮解除品質にとって、より多くの方向性信号を使用すること、あるいは、アンビエントＨＯＡ成分のより良好なモデル化のためにより多くのＨＯＡ係数列を使用することのいずれかによって、現在のＨＯＡ表現がより良好に表現されるかというという課題を考慮して、抽出される方向性信号の数が決定される。抽出される方向性音源の数を決定するための基準をステップ／ステージ２２において導出するために、どの基準が人間の知覚に関連しているか、ＨＯＡ圧縮が、特に、以下の２つの処理によって行われることが考慮される。
−アンビエントＨＯＡ成分を表現するためのＨＯＡ係数列の低減（これは、関連するチャンネルの数の低減を意味する）
−方向性信号およびアンビエントＨＯＡ成分を表現するためのＨＯＡ係数列の知覚符号化

【0041】

抽出された方向性信号の数Ｍ（０≦Ｍ≦Ｄ）に依存して、１番目の処理により、下記の式に従って近似計算が行われる。

【数6】

[この文献は図面を表示できません]

【数7】

[この文献は図面を表示できません]

ここで、

【数8】

[この文献は図面を表示できません]

は、Ｍ個の個々に考慮される音源によって形成されるとするＨＯＡ音場成分

[この文献は図面を表示できません]

から構成される方向性成分のＨＯＡ表現を示し、

[この文献は図面を表示できません]

は、Ｉ−Ｍ個の非零ＨＯＡ係数列のみを用いたアンビエント成分のＨＯＡ表現を示している。

【0042】

２番目の処理からの近似計算を下記の式によって表現することができる。

【数9】

[この文献は図面を表示できません]

【数10】

[この文献は図面を表示できません]

ここで、

[この文献は図面を表示できません]

および

[この文献は図面を表示できません]

は、それぞれ、知覚復号処理の後に合成された方向性成分およびアンビエントＨＯＡ成分を示している。

【0043】

基準の形成
抽出される方向性信号の数

[この文献は図面を表示できません]

は、合計近似誤差（ここで

[この文献は図面を表示できません]

である）

【数11】

[この文献は図面を表示できません]

が人間の知覚の点で可能な限り顕著とならないように選択される。これを確実にするために、個々のバーク尺度臨界帯域に対する合計誤差の方向性パワー分布は、所定の数Ｑ個のテスト方向Ω_q （ｑ＝１，… ，Ｑ）で考慮される。このテスト方向は、単位球面上でほぼ均一に分布する。より具体的に述べると、ｂ番目の臨界帯域（ｂ＝１，… ，Ｂ）に対する方向性パワー分布は、下記のベクトルによって表現される。

【数12】

[この文献は図面を表示できません]

ベクトルの成分

[この文献は図面を表示できません]

は、方向Ω_q、ｂ番目のバーク尺度臨界帯域、およびｋ番目のフレームに関連する合計誤差

[この文献は図面を表示できません]

のパワーを示す。合計誤差

[この文献は図面を表示できません]

の方向性パワー分布

[この文献は図面を表示できません]

は、元のＨＯＡ表現

[この文献は図面を表示できません]

による下記の方向性知覚マスキングパワー分布と比較される。

【数13】

[この文献は図面を表示できません]

次に、各テスト方向Ω_qおよび臨界帯域ｂに対して、合計誤差の知覚レベル

[この文献は図面を表示できません]

が算出される。知覚レベルは、ここで、本質的に、合計誤差

[この文献は図面を表示できません]

の方向性パワーと方向性マスキングパワーとの比率として下記の式に従って定義される。

【数14】

[この文献は図面を表示できません]

【0044】

「１」を減算し、連続的な最大値を求める処理が行われ、誤差パワーがマスキング閾値未満である限り確実に知覚レベルが零になるようにする。最終的に、抽出される方向性信号の数

[この文献は図面を表示できません]

は、全ての臨界帯域に亘る誤差知覚レベルの最大値の全てのテスト方向に対する平均値が最小になるように、すなわち、下記の式に従って選択される。

【数15】

[この文献は図面を表示できません]

【0045】

なお、代替的には、式（１５）において誤差知覚レベルの最大値を平均化処理によって置き換えることができる。

【0046】

方向性知覚マスキングパワー分布の算出
元のＨＯＡ表現

[この文献は図面を表示できません]

による方向性知覚マスキングパワー分布

[この文献は図面を表示できません]

の算出のために、元のＨＯＡ表現

[この文献は図面を表示できません]

は、テスト方向Ω_q （ｑ＝１，… ，Ｑ）から到来する一般的な平面波

[この文献は図面を表示できません]

によって表現されるようにするために、空間領域に変換される。行列

[この文献は図面を表示できません]

内の一般的な平面波信号

[この文献は図面を表示できません]

を

【数16】

[この文献は図面を表示できません]

のように配列すると、空間領域への変換は、下記の処理によって表現される。

【数17】

[この文献は図面を表示できません]

ここで、Ξは、テスト方向Ω_q （ｑ＝１，… ，Ｑ）に対して以下の式によって定義されるモード行列を示す。

【数18】

[この文献は図面を表示できません]

ここで、Ｓ_q：＝

【数19】

[この文献は図面を表示できません]

元のＨＯＡ表現

[この文献は図面を表示できません]

による、方向性知覚マスキングパワー分布

[この文献は図面を表示できません]

の要素

[この文献は図面を表示できません]

は、個々の臨界帯域ｂに対する一般的な平面波関数

[この文献は図面を表示できません]

のマスキングパワーに対応する。

【0047】

方向性パワー分布の算出
以下の説明において、方向性パワー分布

[この文献は図面を表示できません]

を算出するための以下の２つの代替策が示される。

【0048】

ａ．１つの可能性は、項目Ａ．２の最初に記載されている２つの処理を実行することによって、所望のＨＯＡ表現

[この文献は図面を表示できません]

の近似値

[この文献は図面を表示できません]

を実際に算出することである。次に、合計近似誤差

[この文献は図面を表示できません]

が式（１１）に従って算出される。次に、合計近似誤差

[この文献は図面を表示できません]

が、テスト方向Ω_q （ｑ＝１，… ，Ｑ）から到来する一般的な平面波

[この文献は図面を表示できません]

によって表現されるために、空間領域に変換される。一般的な平面波信号を以下のように表される行列

[この文献は図面を表示できません]

内に配置すると、

【数20】

[この文献は図面を表示できません]

空間領域への変換は、下記の処理によって表現される。

【数21】

[この文献は図面を表示できません]

合計近似誤差

[この文献は図面を表示できません]

の方向性パワー分布

[この文献は図面を表示できません]

の要素

[この文献は図面を表示できません]

は、個々の臨界帯域ｂ内で一般的な平面波関数

[この文献は図面を表示できません]

のパワーを算出することによって取得される。

【0049】

ｂ．代替的な解決法は、

[この文献は図面を表示できません]

の代わりに近似値

[この文献は図面を表示できません]

のみを算出することである。この方法には、個々の信号の複雑な知覚符号化を直接行う必要がないという利点がある。この代わりに、個々のバーク尺度臨界帯域内の知覚量子化誤差のパワーを知ることで十分である。この目的のため、式（１１）に定義された合計近似誤差を、以下の３つの近似誤差の合計として記述することができる。

【数22】

[この文献は図面を表示できません]

【数23】

[この文献は図面を表示できません]

【数24】

[この文献は図面を表示できません]

この３つの近似誤差は、互いに独立しているものと仮定することができる。この独立性のため、合計誤差

[この文献は図面を表示できません]

の方向性パワー分布は、３つの個々の誤差

[この文献は図面を表示できません]

、

[この文献は図面を表示できません]

、および

[この文献は図面を表示できません]

の方向性パワー分布の合計として表現することができる。

【0050】

以下、個々のバーク尺度臨界帯域に対する３つの誤差の方向性パワー分布をどのように算出するかについて記載する。

【0051】

ａ．誤差

[この文献は図面を表示できません]

の方向性パワー分布を算出するために、まず、下記の式によって、空間領域への変換が行われる。

【数25】

[この文献は図面を表示できません]

ここで、近似誤差

[この文献は図面を表示できません]

は、したがって、テスト方向Ω_q （ｑ＝１，… ，Ｑ）から到来する一般的な平面波

[この文献は図面を表示できません]

によって表現され、これは、下記の式に従って、行列

[この文献は図面を表示できません]

内に配列される。

【数26】

[この文献は図面を表示できません]

結果として、近似誤差

[この文献は図面を表示できません]

の方向性パワー分布

[この文献は図面を表示できません]

の要素

[この文献は図面を表示できません]

は、個々の臨界帯域ｂ内で、一般的な平面波関数

[この文献は図面を表示できません]

のパワーを算出することによって取得される。

【0052】

ｂ．誤差

[この文献は図面を表示できません]

の方向性パワー分布

[この文献は図面を表示できません]

を算出するために、方向性信号

[この文献は図面を表示できません]

を知覚符号化することによって、この誤差が方向性ＨＯＡ成分

[この文献は図面を表示できません]

に導入されることに留意すべきである。さらに、方向性ＨＯＡ成分が式（８）によって与えられることを考慮すべきである。そして、簡略化のために、ＨＯＡ成分

[この文献は図面を表示できません]

が、空間領域内で、Ο個の一般的な平面波関数

[この文献は図面を表示できません]

によって、等価的に表現されるものと仮定する。これは、単なるスケーリングによって、すなわち、下記の式に従って方向性信号

[この文献は図面を表示できません]

から形成される。

【数27】

[この文献は図面を表示できません]

ここで、

[この文献は図面を表示できません]

は、スケーリング・パラメータを示している。各々の平面波方向

[この文献は図面を表示できません]

は、単位球面上で均一に分布し、

[この文献は図面を表示できません]

が方向推定値

[この文献は図面を表示できません]

と対応するように、回転されるものと仮定される。したがって、スケーリング・パラメータ

[この文献は図面を表示できません]

は「１」である。

【0053】

回転された方向

[この文献は図面を表示できません]

に対して

[この文献は図面を表示できません]

をモード行列として定義し、

【数28】

[この文献は図面を表示できません]

に従ってベクトル内の全てのスケーリング・パラメータ

[この文献は図面を表示できません]

を配列すると、ＨＯＡ成分

[この文献は図面を表示できません]

を下記の式のように記述することができる。

【数29】

[この文献は図面を表示できません]

【0054】

結果として、真の方向性ＨＯＡ成分

【数30】

[この文献は図面を表示できません]

と、

【数31】

[この文献は図面を表示できません]

【数32】

[この文献は図面を表示できません]

によって知覚復号された方向性信号

[この文献は図面を表示できません]

（ｄ＝１，… ，Ｍ）が合成されたものとの間の誤差

[この文献は図面を表示できません]

（式（２３）参照）は、下記の式で表される知覚符号化誤差

【数33】

[この文献は図面を表示できません]

の点で個々の方向性信号において下記の式によって表現することができる。

【数34】

[この文献は図面を表示できません]

【0055】

テスト方向Ω_q （ｑ＝１，… ，Ｑ）に対して、空間領域内の誤差

[この文献は図面を表示できません]

の表現は、下記の式によって与えられる。

【数35】

[この文献は図面を表示できません]

【0056】

ベクトルの要素β^（ｄ）（ｋ）を

[この文献は図面を表示できません]

と表し、個々の知覚符号化誤差

[この文献は図面を表示できません]

が互いに独立しているものと仮定することにより、式（３５）から、知覚符号化誤差

[この文献は図面を表示できません]

の方向性パワー分布

[この文献は図面を表示できません]

の要素

[この文献は図面を表示できません]

は、下記の式によって算出することができる。

【数36】

[この文献は図面を表示できません]

[この文献は図面を表示できません]

は、方向性信号

[この文献は図面を表示できません]

におけるｂ番目の臨界帯域内の知覚量子化誤差のパワーを表現するように想定されている。このパワーは、方向性信号

[この文献は図面を表示できません]

の知覚マスキングパワーに対応するものとすることができる。

【0057】

ｃ．アンビエントＨＯＡ成分のＨＯＡ係数列の知覚符号化の結果として得られる誤差

[この文献は図面を表示できません]

の方向性パワー分布

[この文献は図面を表示できません]

を算出するために、各ＨＯＡ係数列が独立して符号化されるものとする。したがって、各バーク尺度臨界帯域内の個々のＨＯＡ係数列内に導入される誤差は、相関性がないとすることができる。これは、誤差

[この文献は図面を表示できません]

の係数間相関行列は、各バーク尺度臨界帯域に対して対角である、すなわち、下記の式で表される。

【数37】

[この文献は図面を表示できません]

要素

[この文献は図面を表示できません]

は、

[この文献は図面を表示できません]

内のｏ番目の符号化されたＨＯＡ係数列におけるｂ番目の臨界帯域内の知覚量子化誤差のパワーを表現するものとする。これは、ｏ番目のＨＯＡ係数列

[この文献は図面を表示できません]

の知覚マスキングパワーに対応するものと仮定することができる。したがって、知覚符号化誤差

[この文献は図面を表示できません]

の方向性パワー分布は、下記の式によって算出される。

【数38】

[この文献は図面を表示できません]

【0058】

Ｂ．改良されたＨＯＡ圧縮解除
対応するＨＯＡ圧縮解除処理が図３に示されており、このＨＯＡ圧縮解除処理は、以下のステップまたはステージを含む。

【0059】

ステップまたはステージ３１において、

[この文献は図面を表示できません]

内の復号された信号を取得するために、

[この文献は図面を表示できません]

内に含まれるＩ個の信号の知覚復号処理が行われる。

【0060】

信号再配分ステージまたはステージ３２において、

[この文献は図面を表示できません]

内の知覚復号された信号は、方向性信号のフレーム

[この文献は図面を表示できません]

およびアンビエントＨＯＡ成分のフレーム

[この文献は図面を表示できません]

を再形成するために再配分される。インデックスのデータセット

[この文献は図面を表示できません]

および

[この文献は図面を表示できません]

を使用して、ＨＯＡ圧縮に対して行われる割り当て処理を再現することによって、どのように信号を再配分するかについての情報が取得される。これは、再帰的な処理であるため（項目Ａ参照）、例えば、送信に不具合が発生しているような場合に再配分処理を初期化できるようにするために、追加的に送信される割り当てベクトルγ（ｋ）を使用することができる。

【0061】

合成ステップまたはステージ３３において、（欧州特許出願第１２３０６５６９号の図２ｂおよび図４に関連して記載されている処理に従って、）方向性信号のフレーム

[この文献は図面を表示できません]

、対応する方向のセット

[この文献は図面を表示できません]

と共にアクティブな方向性信号のインデックスのセット

[この文献は図面を表示できません]

、方向性信号からのＨＯＡ表現の部分を予測するためのパラメータζ（ｋ−２）、および低減されたアンビエントＨＯＡ成分のＨＯＡ係数列のフレーム

[この文献は図面を表示できません]

を使用して、所望の合計ＨＯＡ表現の現在のフレーム

[この文献は図面を表示できません]

が再合成される。

[この文献は図面を表示できません]

は、欧州特許出願第１２３０６５６９号における

[この文献は図面を表示できません]

に対応し、

[この文献は図面を表示できません]

および

[この文献は図面を表示できません]

は、欧州特許出願第１２３０６５６９号における

[この文献は図面を表示できません]

に対応する。ここでアクティブな方向性信号のインデックスは、

[この文献は図面を表示できません]

の行列要素においてマーク付けされる。すなわち、均一に分布する方向に対する方向性信号は、予測のための受信済のパラメータ（ζ（ｋ−２）)を使用して方向性信号

[この文献は図面を表示できません]

から予測される。その後、現在の圧縮解除されたフレーム

[この文献は図面を表示できません]

が、方向性信号

[この文献は図面を表示できません]

のフレーム、予測された部分および低減されたアンビエントＨＯＡ成分

[この文献は図面を表示できません]

から再合成される。

【0062】

Ｃ． 高次アンビソニックスの基礎
高次アンビソニックス（ＨＯＡ）は注目されるコンパクトな領域内の音場の記述に基づいていており、音源が存在しないものと仮定される。その場合、注目領域内の時間ｔおよび位置ｘでの音圧ｐ（ｔ，ｘ）の空間時間的な挙動は、均質媒質の波動方程式によって物理的に完全に求められる。以下の内容は、図４に示された球面座標システムに基づいている。使用されている座標システムにおいて、ｘ軸は前方の位置を指し、ｙ軸は左側を指し、ｚ軸は上方を指す。空間内の位置ｘ＝（ｒ，θ，φ）^Ｔは、半径ｒ＞０（すなわち、座標原点への距離）、極軸ｚから測定される傾斜角θ∈［０，π］、さらに、ｘ軸からの、ｘ−ｙ平面内で反時計周りに測定される、方位角φ∈［０，２π］によって表される。さらに、（・）^Ｔは、転置を表す。

【0063】

Ｆ_ｔ（・）によって表される時間に対する音圧のフーリエ変換、すなわち、

【数39】

[この文献は図面を表示できません]

は下記の式に従った一連の球面調和関数に拡張される（Ｅ.Ｇ. Ｗｉｌｌｉａｍｓ著“Ｆｏｕｒｉｅｒ Ａｃｏｕｓｔｉｃｓ（フーリエ・アコースティックス））”、応用数理科学、第９３巻、アカデミックプレス社、１９９９年参照）。ここで、ωは角周波数を表し、ｉは虚数単位を表す。

【数40】

[この文献は図面を表示できません]

式（４０）において、ｃ_ｓは音速を示し、ｋは角波数を示し、この角波数ｋはｋ＝ｗ／ｃ_ｓによって角周波数ωに関連している。さらに、ｊ_ｎ（・）は、第１種球ベッセル関数を表しており、

[この文献は図面を表示できません]

は、Ｃ．１の項目で定義されている次数ｎおよび位数ｍの実数値の球面調和関数を示している。展開係数

[この文献は図面を表示できません]

は、角波数ｋのみに依存する。上述した内容において、音圧は、空間的に帯域制限されているものと暗黙的に仮定されている。したがって、球面調和関数の級数が次数インデックスｎに対して上限Ｎで打ち切られ、これは、ＨＯＡ表現の次数と呼ばれる。

【0064】

音場が相異なる角周波数ωの調和平面波の無限個の重ね合わせによって表現され、角の組（θ，φ）によって特定される全ての想定可能な方向から到来する場合には、各々の平面波複素振幅関数Ｃ（ω，θ，φ）は、下記の球面調和展開によって表すことができることが分かる（Ｂ. Ｒａｆａｅｌｙ著、“Ｐｌａｎｅ−ｗａｖｅ Ｄｅｃｏｍｐｏｓｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｕｎｄ Ｆｉｅｌｄ ｏｎ ａ Ｓｐｈｅｒｅ ｂｙ Ｓｐｈｅｒｉｃａｌ Ｃｏｎｖｏｌｕｔｉｏｎ（球面畳み込みによる球面上の音場の平面波分解）”、米国音響学会誌４(１１６)、２１４９−２１５７頁、２００４年参照）。

【数41】

[この文献は図面を表示できません]

ここで、展開係数

[この文献は図面を表示できません]

は、展開係数

[この文献は図面を表示できません]

と下記の式によって関連する。

【数42】

[この文献は図面を表示できません]

【0065】

個々の係数

[この文献は図面を表示できません]

が角周波数ωの関数であると仮定すると、逆フーリエ変換（

[この文献は図面を表示できません]

）によって示される）を適用することにより、下記の時間領域関数をもたらす。

【数43】

[この文献は図面を表示できません]

これは、各次数ｎおよび位数ｍに対して、下記の単一のベクトルｃ（ｔ）にまとめられる。

【数44】

[この文献は図面を表示できません]

ベクトルｃ（ｔ）内の時間領域関数

[この文献は図面を表示できません]

の位置インデックスは、ｎ（ｎ＋１）＋１＋ｍによって与えられる。ベクトルｃ（ｔ）内の要素の総計は、Ο＝（Ｎ＋１）^２によって与えられる。

【0066】

最終的なアンビソニックス形式は、サンプリング周波数ｆ_ｓを使用して、下記のｃ（ｔ）のサンプリングされたバージョンをもたらす。

【数45】

[この文献は図面を表示できません]

ここで、Ｔ_ｓ＝１／ｆ_ｓは、サンプリング期間を示す。ｃ（ｌＴ_ｓ）の要素は、アンビソニックス係数として参照される。時間領域信号

[この文献は図面を表示できません]

は実数値であり、したがって、アンビソニックス係数は実数値である。

【0067】

Ｃ．１ 実数値の球面調和関数の定義
実数値の球面調和関数

[この文献は図面を表示できません]

は、下記の式によって与えられる。

【数46】

[この文献は図面を表示できません]

ここで

【数47】

[この文献は図面を表示できません]

関連するルジャンドル関数Ｐ_ｎ，ｍ（ｘ）は、下記の式で定義される。

【数48】

[この文献は図面を表示できません]

ここで、ルジャンドル多項式Ｐ_ｎ（ｘ）を用い、上述した、Ｅ．Ｇ．Ｗｉｌｌｉａｍｓ著の文献の場合とは異なり、コンドン-ショートレーの位相項（−１）^ｍを用いない。

【0068】

Ｃ．２ 高次アンビソニックスの空間解像度
方向Ω_０＝（θ_０，φ_０）^Ｔから到来する一般的な平面波関数ｘ（ｔ）は、下記の式によってＨＯＡにおいて表現される。

【数49】

[この文献は図面を表示できません]

平面波振幅の対応する空間密度

[この文献は図面を表示できません]

は、下記の式によって与えられる。

【数50】

[この文献は図面を表示できません]

【数51】

[この文献は図面を表示できません]

【0069】

式（５１）から理解されるように、これは、一般的な平面波関数ｘ（ｔ）と空間分散関数ν_Ｎ（Θ）との積であり、空間分散関数ν_Ｎ（Θ）は、下記の式の特性を有するΩとΩ_０との間の角度Θのみに依存するように示されている。

【数52】

[この文献は図面を表示できません]

想定のとおり、無限次元の極限、つまり、Ｎ→∞である場合において、空間分散関数は
ディラックのデルタ関数δ（・）、すなわち、下記のように変化する。

【数53】

[この文献は図面を表示できません]

【0070】

しかしながら、有限次元Ｎの場合には、方向Ω_０からの一般的な平面波の寄与は、近隣の方向ににじみ、このにじみの度合いは次数の増加に伴い減少する。Ｎの複数の異なる値に対する正規化された関数ν_Ｎ（Θ）のプロットが図５に示されている。

【0071】

任意の方向Ωでの平面波振幅の空間密度の時間領域の挙動は、他の任意の方向での平面波振幅の空間密度の時間領域の挙動の倍数となることが指摘される。特に、時間ｔに対して、何らかの所定方向Ω_１およびΩ_２についての関数ｃ（ｔ，Ω_１）およびｃ（ｔ，Ω_２）は、高い相関性がある。

【0072】

Ｃ．３ 球面調和関数変換
平面波振幅の空間密度がΟ個の空間方向Ω_ｏ（１≦ο≦Ο）で離散化される場合、空間方向Ω_ｏは単位球面上でほぼ均一に分布するのだが、Ο個の方向性信号ｃ（ｔ，Ω_ｏ）が取得される。これらの信号をベクトルにまとめると、下記の式で表され、

【数54】

[この文献は図面を表示できません]

式（５０）を使用してこのベクトルを、下記のような単純な行列乗算によって式（４４）に定義される連続的なアンビソニックス表現ｃ（ｔ）から計算可能であることを検証できる。
ｃ_ＳＰＡＴ（ｔ）＝Ψ^Ｈｃ（ｔ） （５５）
ここで、（・）^Ｈは、複素共役転置を示し、Ψは、下記の式によって定義されるモード行列を表す。

【数55】

[この文献は図面を表示できません]

ここで、

【数56】

[この文献は図面を表示できません]

【0073】

方向Ω_ｏは単位球面上にほぼ均一に分布しているため、一般的には、モード行列は、可逆である。したがって、連続的なアンビソニックス表現は、方向性信号ｃ（ｔ，Ω_ｏ）から下記の式によって計算することができる。

【数57】

[この文献は図面を表示できません]

【0074】

双方の式は、アンビソニックス表現と空間領域との間の変換および逆変換を構成する。本願において、これらの変換は、球面調和関数変換および逆球面調和関数変換と呼ばれる。

【0075】

なお、方向Ω_ｏは単位球面上でほぼ均一に分布するため、近似計算

【数58】

[この文献は図面を表示できません]

が利用可能となり、式（５５）において、Ψ^Ｈの代わりにΨ^−１を使用することが正当化される。

【0076】

上述した関係の全てが離散時間領域にも有効であることは有利である。

【0077】

本発明の処理を単一のプロセッサまたは電子回路、または、並列に動作する複数のプロセッサまたは電子回路、および／または、本発明の処理の複数の異なる部分に対して動作する、複数のプロセッサまたは電子回路で実行することができる。

【図1】

[この文献は図面を表示できません]

【図2】

[この文献は図面を表示できません]

【図3】

[この文献は図面を表示できません]

【図4】

[この文献は図面を表示できません]

【図5】

[この文献は図面を表示できません]