(57)【特許請求の範囲】
【請求項１】
高次アンビソニックス（ＨＯＡ）フォーマットの３Ｄオーディオ信号のラウドネスレベルを調整する方法であって、前記３Ｄオーディオ信号は、０次の係数を有する少なくとも１つの係数チャネルと、０よりも高い次数の係数を有する１つ以上の係数チャネルとを有する、方法において、
前記３Ｄオーディオ信号の０次の係数を有する係数チャネルのラウドネスを測定するステップと、
前記０次の係数を有する係数チャネルの前記測定されたラウドネスに従って前記３Ｄオーディオ信号の全ての係数チャネルのためのマスタゲインをレベル調整するステップであって、レベル調整されたＨＯＡ信号が取得されるステップと
を有する方法。
【請求項２】
前記３Ｄオーディオ信号をレンダリングするステップであって、正規化されたレンダリング行列が使用されるステップと、
前記レンダリングされた信号のラウドネスをモニタするステップと
を更に有する請求項１に記載の方法。
【請求項３】
前記レンダリング行列は、Ｄｅがエネルギ保存レンダリング行列であるとして、
【数２９】
0006271586.tif 000032
に従って、又はＤが略エネルギ保存レンダリング行列であるとして、
【数３０】
0006271586.tif 000033
に従って、正規化係数Ｄｎにより正規化されるレンダリング行列を使用する、
請求項２に記載の方法。
【請求項４】
前記レンダリングは、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であり、且つ、Ｄが略エネルギ保存レンダリング行列であるとして、
【数３１】
0006271586.tif 000034
に従って、正規化係数
［外１］
0006271586.tif 000035
により正規化されるレンダリング行列を使用する、
請求項２に記載の方法。
【請求項５】
前記レベル調整するステップにより得られる信号をレンダリングし、該レンダリングされた信号を送信するステップを更に有し、
前記レンダリングは、Ｄが略エネルギ保存レンダリング行列であり、且つ、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であるとして、
【数３２】
0006271586.tif 000036
に従って、又は
【数３３】
0006271586.tif 000037
に従って、正規化係数により正規化されるエネルギ保存レンダリング行列を使用する、
請求項１乃至４のうちいずれか一項に記載の方法。
【請求項６】
前記レベル調整するステップにより得られる信号を送信し、受信し、レンダリングするステップを更に有し、
前記レンダリングは、送信及び受信の後に、Ｄが略エネルギ保存レンダリング行列であり、且つ、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であるとして、
【数３４】
0006271586.tif 000038
に従って、又は
【数３５】
0006271586.tif 000039
に従って、正規化係数により正規化されるエネルギ保存レンダリング行列を使用する、
請求項１乃至４のうちいずれか一項に記載の方法。
【請求項７】
前記係数チャネルのラウドネスを測定するステップ及び前記マスタゲインをレベル調整するステップは、送信側において、送信又は記憶のステップ前に行われる、
請求項１乃至６のうちいずれか一項に記載の方法。
【請求項８】
高次アンビソニックス（ＨＯＡ）フォーマットの３Ｄオーディオ信号のラウドネスレベル調整のための装置であって、前記３Ｄオーディオ信号は、０次の係数を有する少なくとも１つの係数チャネルと、０よりも高い次数の係数を有する１つ以上の係数チャネルとを有する、装置において、
前記０次の係数を有する係数チャネルのラウドネスを測定するラウドネス測定部と、
前記ラウドネス測定部から受け取った前記測定されたラウドネスに従って前記３Ｄオーディオ信号の全ての係数チャネルのゲインを調整し、レベル調整されたＨＯＡ信号を供給するマスタゲイン部と
を有する装置。
【請求項９】
正規化されたレンダリング行列が使用されて、前記レベル調整されたＨＯＡ信号をレンダリングし、該レンダリングされた信号のラウドネスをモニタし、該モニタされたラウドネスが前記ラウドネス測定部を制御するのに使用されるモニタ部
を更に有する請求項８に記載の装置。
【請求項１０】
前記ラウドネス測定部は、瞬時的なのラウドネス測定、短期のラウドネス測定、及び積分されたラウドネス測定を前記マスタゲイン部へ供給する、
請求項８又は９に記載の装置。
【請求項１１】
前記マスタゲイン部から得られた前記レベル調整されたＨＯＡ信号をレンダリングするレンダラであって、エネルギ保存レンダリング行列が使用され、該エネルギ保存レンダリング行列は、Ｄが略エネルギ保存レンダリング行列であり、且つ、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であるとして、
【数３６】
0006271586.tif 000040
に従って、又は
【数３７】
0006271586.tif 000041
に従って、正規化係数により正規化される、レンダラと、
前記レンダリングされた信号を送信する送信器と
を更に有する請求項８乃至１０のうちいずれか一項に記載の装置。
【請求項１２】
前記マスタゲイン部から得られた前記レベル調整されたＨＯＡ信号を送信する送信器と、
前記送信された信号を受信する受信器と、
前記受信された信号をレンダリングするレンダラであって、エネルギ保存レンダリング行列が使用され、該エネルギ保存レンダリング行列は、Ｄが略エネルギ保存レンダリング行列であり、且つ、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であるとして、
【数３８】
0006271586.tif 000042
に従って、又は
【数３９】
0006271586.tif 000043
に従って、正規化係数により正規化される、レンダラと
を更に有する請求項８乃至１０のうちいずれか一項に記載の装置。
【請求項１３】
前記マスタゲイン部から得られた前記レベル調整されたＨＯＡ信号を記憶する記憶部と、
前記記憶された信号を取り出す検索部と、
前記取り出された信号をレンダリングするレンダラであって、エネルギ保存レンダリング行列が使用され、該エネルギ保存レンダリング行列は、Ｄが略エネルギ保存レンダリング行列であり、且つ、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であるとして、
【数４０】
0006271586.tif 000044
に従って、又は
【数４１】
0006271586.tif 000045
に従って、正規化係数により正規化される、レンダラと
を更に有する請求項８乃至１０のうちいずれか一項に記載の装置。
【請求項１４】
レンダリング行列を使用することによって動作するラウドネス維持ＨＯＡレンダラを有し、前記レンダリング行列が正規化されることを特徴とするオーディオ再生装置。
【請求項１５】
前記レンダリング行列は、実質的にエネルギを保存し、
前記レンダリング行列は、Ｄが略エネルギ保存レンダリング行列であり、且つ、ｙが方向Ωｃの球面調和関数から成るＨＯＡテスト信号であるとして、
【数４２】
0006271586.tif 000046
に従って、又は
【数４３】
0006271586.tif 000047
に従って、あるいは、Ｄｅがエネルギ保存レンダリング行列であるとして、
【数４４】
0006271586.tif 000048
に従って、正規化部において正規化される、
請求項１４に記載のオーディオ再生装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、高次アンビソニックス（ＨＯＡ；Higher-Order  Ambisonics）コンテンツのためのラウドネスレベル調整に関する。特に、それは、ＨＯＡラウドネスレベルを測定する方法及びＨＯＡラウドネスレベルを測定する装置に関する。
【背景技術】
【０００２】
この項目は、以下で記載及び／又は請求される本発明の様々な態様に関連する可能性がある技術の様々な態様を読者に紹介することを目的とする。この議論は、本発明の様々な態様のより良い理解を促すために読者に背景情報を提供するのに役立つと思料される。然るに、それらの記述は、先行技術の承認としてではなく、この観点において読まれるべきであることが理解されるべきである。
【０００３】
例えばアンビソニックスなどの音場信号は、音場の表現を搬送する。アンビソニックス・フォーマットは、音場の球面調和関数展開に基づく。基本のアンビソニックス・フォーマット又はＢ−フォーマットは、次数０及び１の球面調和関数を使用し、一方、いわゆる高次アンビソニックス（ＨＯＡ）は、少なくとも２次の更なる球面調和関数も使用する。すなわち、ＨＯＡ信号は、例えば、次数０の信号（Ｗ−チャネル，Ｎ＝０）、次数１（Ｎ＝１）の１つ以上の信号、次数２（Ｎ＝２）の１つ以上の信号などのような、異なった次数Ｎの異なった部分信号を有する。復号化プロセスは、個々のラウドスピーカ信号を得るために必要とされる。オーディオシーンを合成するために、空間ラウドスピーカ配置を示すパン関数は、所与の音源の空間ローカリゼーションを得るために必要とされる。
【０００４】
復号器側で行われる１つのタスクは、再生レベルを定めることである。先行技術［１］に記載されており且つ図１に示されるように、各ラウドスピーカフィードの増幅器ゲインＧｌは、−１８ｄＢＦＳｒｍｓで入力されるデジタル全帯域ピンクノイズが７８±５ｄＢＡの音圧レベル（ＳＰＬ；Sound  Pressure  Level）を引き起こすように設定される。図１において、ピンクノイズ・テスト信号は、夫々のラウドスピーカについて個々に、増幅器１２におけるスピーカ増幅を調整することによって各ラウドスピーカ１３の音圧レベルをレベル調整するために使用される。デジタルのピンクノイズ・テスト信号は、Ｄ／Ａコンバータにおいてアナログ信号に変換される。ミキシング及び提示の現場でのＳＰＬレベル調整、並びにミキシング部屋でのコンテンツのラウドネスレベル調整は、プログラム又はアイテムを切り替えるときに、一定の知覚ラウドネスを可能にする。
【０００５】
［コンテンツ・ラウドネスレベル較正］
ミキシング設備及び提示現場の再生レベルがこのようにして設定される場合に、アイテム又はプログラムを切り替えることは、更なるレベル調整なしで可能であるべきである。チャネルに基づくコンテンツに関して、これは、コンテンツがミキシング場所で好ましいラウドネスレベルに調整される場合に簡単に達成される。好ましい聴取レベルの基準は、アイテム全体自体又はアンカー信号のラウドネスのいずれかであることができる。
【０００６】
基準としてアイテム全体自体を使用することは、コンテンツがファイルとして記憶される場合に、‘ショート形式のコンテンツ’にとって有用である。聴取による調整に加えて、ＥＢＵ  Ｒ１２８［２］に従うラウドネス・ユニット・フル・スケール（ＬＵＦＳ；Loudness  Units  Full  Scale）におけるラウドネスの調整は、コンテンツのラウドネス調整のために使用され得る。ＬＵＦＳの代替名は、ＩＴＵ−Ｒ  ＢＳ．１７７０［３］からの‘Loudness,  K-weighted,  relative  to  Full  Scale’（１ＬＵＦＳ＝１ＬＫＦＳ）である。あいにく、［２］における解決法は、最大で５チャネル・サラウンドまでのセットアップのためにしかコンテンツをサポートしない。２２チャネル・ファイルのラウドネス測定は、全２２チャネルが１の等しいチャネル重みによって因数分解される場合に、知覚ラウドネスと相関し得るが、依然として完全な聴取テストによる証拠又は裏付けがない。基準として例えばダイアログなどのアンカー信号を使用する場合に、レベルはこの信号に対して選択される。これは、例えばフィルム音響、ライブ録音及び放送などの‘ロング形式のコンテンツ’にとって有用である。好ましい聴取レベルを広げる更なる要件は、ここでは話し言葉の明りょう度である。
【０００７】
先と同じく、聴取による調整に加えて、コンテンツは、例えばＡＴＳＣ  Ａ／８５［４］において定義されるように、ラウドネス測定に関して正規化されてよい。コンテンツの最初の部分はアンカー部分として識別される。次いで、［３］において定義される測定が計算されるか、あるいは、それらの信号及び目標ラウドネスを達成するためのゲイン係数が決定される。ゲイン係数は、アイテム全体を調整するために使用される。あいにく、先と同じく、サポートされるチャネルの最大数は５に制限される。
【０００８】
ＩＴＵ−Ｒ  ＢＳ．１７７０［３］からの図２は、ＥＢＵ  Ｒ１２８［２］及びＡＴＳＣ  Ａ／８５［４］において使用されるラウドネス測定を示す。［２］は、コンテンツアイテム全体の測定されたラウドネスを−２３ｄＢＬＫＦＳへとゲイン調整することを提案する。［４］において、アンカー信号のラウドネスのみが測定され、コンテンツは、アンカー部分が−２４ｄＢＬＫＦＳの目標ラウドネスを達成するようにゲイン調整される。様々な入力信号Ｌ，Ｒ，Ｃ，Ｌｓ，ＲｓはＫ−フィルタ２１においてフィルタリングされ、各チャネルの電力は電力アベレージャ２２において平均化され、各チャネルは２３で重み付けされ、そして、重み付けされた信号は、測定されたラウドネス値を２５で得るよう２４で足し合わされる。
【０００９】
芸術的な考慮から、コンテンツはミキシングスタジオで調整されるべきである。これは、個人の聴取によってなされる。自動ラウドネス測定は、サポートとして、特定のラウドネスが超えられないことを示すために、使用され得る。
【００１０】
ＨＯＡ及びオーディオ・オブジェクト（ＡＯ；Audio  Object）に基づくコンテンツ、尚また、異なる数又は異なる位置のラウドスピーカにリミックスされる必要があるチャネルベースのコンテンツに関し、レンダリングは考慮に入れられるべきである。レンダラが満たすべき特性が存在する必要があり、そのようなレンダラは、ミキシングスタジオで且つ消費者の提示現場で使用されるべきである。
【発明の概要】
【発明が解決しようとする課題】
【００１１】
上記を鑑み、解決されるべき１つの課題は、復号されたオーディオデータの再生が、特に、異なるオーディオアイテム間を切り替えるときに、又は異なるラウドスピーカ・セットアップにレンダリングするときに、等しいラウドネスにより知覚されるべきことである。本発明は、少なくともこの課題を解決する。
【課題を解決するための手段】
【００１２】
原理上は、発明は、オーディオ信号のＷ−チャネル（すなわち、音場信号の０次の部分信号）のみのラウドネスを測定することと、Ｗ−チャネルの測定されたラウドネスに従って信号の全チャネルのためのマスタゲインを平均にすること（すなわち、マスタゲインのレベルを調整すること）とを有する。
【００１３】
発明の一実施形態に従って、ＨＯＡオーディオ信号のラウドネスレベルを調整する方法は、ＨＯＡオーディオ信号のＷ−チャネルのラウドネスを測定するステップと、所望のラウドネスレベルを得るために、Ｗ−チャネルの測定されたラウドネスに従ってＨＯＡ信号の全ての係数チャネルのためのＨＯＡマスタゲインを平均にするステップとを有する。
【００１４】
更に、発明の一実施形態において、ＨＯＡラウドネスレベル調整のための装置は、ＨＯＡ信号のＷ−チャネルのラウドネスを測定するＨＯＡラウドネス測定部と、ＨＯＡラウドネス測定部から受け取ったＷ−チャネルの測定されたラウドネスに従ってＨＯＡ信号の全ての係数チャネルのためのゲインを調整するＨＯＡマスタゲイン部とを有する。
【００１５】
有利に、Ｗ−チャネルのラウドネスを測定し、ＨＯＡマスタゲインを平均にすることは、送信側で、すなわち、ＨＯＡフォーマットのオーディオ信号を送信又は記憶する前に、行われ得る。発明の更なる有利な実施形態は、従属請求項、以下の記載及び図において開示される。
【図面の簡単な説明】
【００１６】
【図１】各ラウドスピーカの音圧レベルをレベル調整するためにピンクノイズ・テスト信号を使用する既知の解決法を示す。
【図２】既知のシステムにおいて使用されるＩＴＵ−Ｒ  ＢＳ．１７７０ラウドネス測定を示す。
【図３】発明の一実施形態に従う、特別な所与のラウドスピーカ・セットアップへのコンテンツのレンダリングのための構造を示す。
【図４】エネルギ保存レンダリング行列の実際的な例についての空間セットアップ及びエネルギ分布を示す。
【図５】次数Ｎ＝１〜６をレンダリングする２２．０チャネル・セットアップへレンダリングされる３つのテストアイテムの両耳ラウドネス値を示す。
【図６】発明に従うレンダリングにより５つの異なるスピーカ・セットアップへレンダリングされたテストアイテムの両耳ラウドネス値を示す。
【図７】レンダリング後の平均ラウドネスと比較されるＨＯＡ  Ｗ−チャネル信号のラウドネス値を示す。
【図８】例えば、ＨＯＡライブ放送のために有用な、本発明に従うシステムを示す。
【図９】チャネルベースの送信系のためのシステムを示す。
【発明を実施するための形態】
【００１７】
発明の例となる実施形態は、添付の図面を参照して記載される。
【００１８】
図３は、コンテンツがいずれかの特定のスピーカ・セットアップへレンダリングされるオーディオ再生装置又はその部分を示す。コンテンツのフォーマットは、高次アンビソニックス（ＨＯＡ）、オーディオ・オブジェクト（ＡＯ）又はチャネルベース（その場合に、それは異なるセットアップへレンダリングされる。）であることができる。ラウドスピーカフィードｗの再生は、特に、アイテムを切り替えるときに、又は異なるラウドスピーカ・セットアップへレンダリングするときに、等しいラウドネスにより知覚されるべきである。入力ベクトルｂは、ひと組の信号の単一の時間サンプルを、その単一の時間サンプルのレンダリングの前について記述する。ベクトルｗは、ひと組のスピーカフィードの単一の時間サンプルを、レンダラ３１におけるその単一の時間サンプルのレンダリング後について記述する。レンダラ３１の特性は、生成において且つ消費者側において同じである必要があり、あるいは、少なくとも、生成側において知られている必要がある。レンダリング後、ラウドスピーカフィードｗは、Ｄ／Ａコンバータ３２においてデジタルからアナログに変換され、増幅器３３において増幅され、そして、次いで、ラウドスピーカ３４に供給される。
【００１９】
以下で、所望のレンダリング特性は、３つ全てのコンテンツフォーマット（ＨＯＡ、ＡＯ、チャネルベース）から始まって、次いで、ＨＯＡの詳細に入りながら、導出される。ＨＯＡ信号は、例えば、次数０の信号（Ｗ−チャネル、Ｎ＝０）、次数１（Ｎ＝１）の１つ以上の信号、次数２（Ｎ＝２）の１つ以上の信号などのような、異なった次数Ｎの異なった部分信号を有することが知られる。
【００２０】
［信号エネルギ及び知覚ラウドネス］
理論解析を簡単にするために、我々は、次のように定義されるτ時間サンプルブロックの信号ｗ（図３）のエネルギに集中する（詳細は以下を参照）：
【数１】
0006271586.tif 000002
ここで、Ｗは、Ｌ個のラウドスピーカ・チャネル及びτ個のサンプルの行列であり、Ｗｌ，ｉは行列要素であり、ｌはスピーカのインデックスを表し、ｉはサンプルのインデックスを表す。エネルギＥＷは、［５］から知られるような２００Ｈｚを超える周波数についてのアンカー信号の知覚ラウドネスの公正な推定を与える。Ｋ−フィルタは２００Ｈｚよりも低い周波数を抑えるので［２］、ＥＷはラウドネス測定に略比例する。直接の比較のためのラウドネス測定及びエネルギ測定は、表１において以下で与えられる。明らかになるように、ＨＯＡ  Ｗ−チャネルのラウドネスレベル及び信号エネルギ（列“HOA(W)”）は、レンダリングされる信号（列“Rendered”）及び両耳信号（列“Binaural”）のラウドネスレベル及び信号エネルギに極めて近く、実質的に等しい。値は、見本となるように、２２．０セットアップ（ＮＨＫ）を参照しており、すなわち、２２個のラウドスピーカ及び４次のＨＯＡ信号（Ｌ＝２２，Ｎ＝４）による。
【００２１】
【表１】
0006271586.tif 000003
［レンダラの要件］
次の簡略化が利用される。すなわち、一般性を失わずに（すなわち、それは前提条件ではない。）、ＤＣ成分を含まないテスト信号ｘ（ｔ）が考えられる。Ａをこの信号のＲＭＳ値とし、ＥＡ＝Ａ２をそのエネルギとする。テスト信号は、方向Ω＝（θ，φ）Ｔの単一のソース信号であり、θは仰角であり、φは方位角である。エネルギ的考察は、最適な聴取エリア、いわゆるスイートエリアのために適用される。
【００２２】
テスト信号は、チャネルベースの信号、オーディオ・オブジェクト、又はＨＯＡ信号として表現され得る。以下の表はそれらの表現を提示し、ｂは抽象入力ベクトルであり、レンダリング後のチャネルベースの表現はベクトルｗによって表され、レンダリング後のエネルギはＥＷであり、レンダリングの等しいエネルギ（ラウドネス）のための条件はＥＡ＝ＥＷである。Ｄを、一般化されたレンダリング行列であるとする：
【数２】
0006271586.tif 000004
表２で、エネルギ保存要件ＥＡ＝ＥＷは、オーディオ・オブジェクトとして（上の行）、ＨＯＡ信号として（真ん中の行）、又はチャネルベースの表現として（下の行）符号化される指向性テスト信号について得られる。ベクトルｂは、レンダリングの前の符号化された入力である。ベクトルｗはレンダリング後の信号（Ｄ／Ａ前のラウドスピーカ信号）である。ＥＷはレンダリング後のエネルギである。ＥＡ＝Ａ２はテスト信号ｘ（ｔ）のエネルギであり、ｔｉは時間サンプルのインデックスである。更に以下で、ＥＷは、見本として、ＨＯＡの場合について導出される。
【００２３】
【表２】
0006271586.tif 000005
オーディオ・オブジェクト（ＡＯ）のレンダリングを考えると、オーディオ・オブジェクトのための考慮は、二乗されたレンダリング重みＣＷの和が１に等しいとのよく知られた要件に至る：
【数３】
0006271586.tif 000006
要件は、独立したＡＯの重ね合わせにも当てはまるが、相関ＡＯについては、チャネルベースの考慮が当てはまる。
【００２４】
チャネルベースのコンテンツのレンダリングのための要件は、以下によって与えられる：
【数４】
0006271586.tif 000007
ここで、ベクトルｃｂは、コンテンツ生成内で使用されるパニング重みを含み、行列Ｄは、一般化されたレンダリング行列（復号化行列とも呼ばれる。）を含む。
【００２５】
２つのタイプのコンテンツが、見本として考えられる：
１．  元のコンテンツは、独立したＡＯのパニングを用いてミキシングされており、二乗されたパニング重みの和は１に等しい。その場合に、レンダリング要件はＤＴＤ＝Ｉ（サイズＬｂ×Ｌｂの単位行列）になる。レンダリング行列のための結果として得られる要件は以下で論じられ、発明に従う解決法はこのようなコンテンツのために使用され得る。
２．  マイクロホン録音及びミキシングされたコンテンツに関し、パニング重みの性質は知られておらず、通常は、異なるラウドスピーカ・セットアップへミキシングする場合に、レンダリングされる信号のエネルギの予測は行われ得ない。よって、ダウンミキシング又は再レンダリング後のラウドネスの予測は一般に不可能であり、発明は適用可能でない。ここで、異なるアプローチのみが適用可能であると考えられる。すなわち、チャネル間の如何なる相関も取り除くための指向性及び無指向性部分の分離、並びに指向性コンテンツのためのオブジェクトに関連したレンダリング方法の使用である。その場合に、アンビエント・ミックスは、オリジナルの指向性対無指向性エネルギ比が保存されるように加えられる。従来の方法は、ミキシングスタジオで生成される新しい芸術的なミックスを作成するために使用されてきた。
【００２６】
エネルギを保存するＨＯＡレンダラに関し、ＨＯＡレンダリング要件は、次のとおりである：
【数５】
0006271586.tif 000008
ＨＯＡベースのコンテンツ内で、如何なる信号も平面波の重ね合わせとして表現される。平面波信号は、次の式によってＨＯＡへと符号化される：
【数６】
0006271586.tif 000009
この式で、ｙは、方向Ωの球面調和関数（ＳＨ；Spherical  Harmonics）から成る符号化ベクトルである。ベクトルｂ及びｙは、ＮがＨＯＡ次数を表すとして、（Ｎ＋１）２個の要素から成る。ＳＨの直交性に起因して、ベクトル積は、次のとおりである（Ｎ３Ｄにより正規化されたＳＨ）：
【数７】
0006271586.tif 000010
Ｉが（Ｎ＋１）２×（Ｎ＋１）２の単位行列であるとして、要件ＤＴＤ＝Ｉは［５］で論じられている。このような‘エネルギ保存’復号化行列を構築するためのプロシージャは、特異値分解（ＳＶＤ；Singular  Value  Decomposition）に基づく。［５］において、エネルギ保存は、次の式によって定義される：
【数８】
0006271586.tif 000011
エネルギ保存の要件を満足するレンダリング行列Ｄｎを導出するステップは、次のとおりである：
１．  エネルギ保存レンダリング行列Ｄｅを導出する。これは、スピーカの数ＬがＨＯＡ係数の数（Ｎ＋１）２よりも大きいか、又は少なくとも等しい限り、多くのスピーカ位置について可能である。
２．  エネルギ保存レンダリング行列をＨＯＡ次数と無関係にする。Ｄｅは、ＥＷ／Ｅｂ＝１であるように構築されるから、要件ＥＷ／Ｅａ＝１は、次の係数によるスケーリングをもたらす：
【数９】
0006271586.tif 000012
実際的なＨＯＡレンダラ構成に関し、次が当てはまる。実際上、スピーカの位置決めはしばしば最適でなく、Ｌ＜（Ｎ＋１）２である。その場合に、レンダリング行列設計は、理論に近づくことしかできない。偏差は、通常は、スピーカ密度が低い領域において最大であり、ＨＯＡ次数を増大させる場合に大きくなる。例は、図４において示される。良好なレンダリング行列設計は、スピーカ密度が低いエリアからのエネルギの減衰を特徴とする。すなわち、そのような方向から衝突する音響は、より小さい音声で知覚されるはずである。
【００２７】
図４は、エネルギ保存レンダリング行列ＤＨＯＡの実際の例を示す。上の並びには、異なるＨＯＡ次数についての様々なラウドスピーカ・セットアップが存在し、下の並びには、夫々の指向性エネルギ保存特性が存在する。球面をサンプリングする均等に割り付けられた方向を持った６７２４個の平面波は、ＨＯＡ符号化され、２２又は１２チャネル・セットアップへレンダリングされる。指向性エネルギ保存は、グレーレベルで示され、異なるセットアップ並びにそれらの最小及び最大値は次のとおりである：
図４ａ）は、ＨＯＡ次数Ｎ＝４を有する２２＿ＮＨＫ＿４５セットアップを示す：［−０．２ｄＢ；０．０ｄＢ］。
図４ｂ）は、ＨＯＡ次数Ｎ＝６を有する２２＿ＮＨＫ＿４５セットアップを示す：［−１．２ｄＢ；０．１ｄＢ］。
図４ｃ）は、ＨＯＡ次数Ｎ＝３を有する１２＿ＳｗｉｓｓＡｕｄｅｃ＿１１０セットアップを示す：［−１．１ｄＢ；０．２ｄＢ］。
【００２８】
明らかなように、信号エネルギは、Ｎ＝４を有する２２チャネル・セットアップについてむしろ均等に分布し、それにより、たった０．２ｄＢの偏差が現れる。Ｎ＝６の２２チャネル・セットアップ及びＮ＝３の１２チャネル・セットアップに関し、信号エネルギは依然としてむしろ均等に分布するが、既に１．３ｄＢの偏差を伴う（ラウドスピーカが利用可能でない方向に関してさえ）。信号エネルギの均等分布は、エネルギ保存レンダリングの利点である。
【００２９】
エネルギ保存は、ラウドスピーカが聴取位置の周りの完全な凸包に及ぶ場合にのみ可能であるから、レンダリング構成のステップは、次のようになる：
１．  略エネルギ保存レンダリング行列Ｄを構築する。２つの適切な方法は［５］、［６］において記載されている。エネルギ保存は、正面の領域において正確であるべきであり、背面下又は側面の領域においてそれほど正確でなくてよい。
２．  コーシー・シュワルツの不等式（以下参照）により、レンダリング行列は、次いで、次の式によって正規化される：
【数１０】
0006271586.tif 000013
これは、近エネルギ保存行列に関して式（６）に取って代わり、式（６）は、式（７）の特別な場合として見なされ得る。代替案は、モードベクトルｙ＝（Ｙ００（Ωｃ），Ｙ−１１（Ωｃ），・・，Ｙｎｍ（Ωｃ），・・，ＹＮＮ（Ωｃ））Ｔ及び方向Ωｃの球面調和関数Ｙｎｍ（Ωｃ）により、中心（Ωｃ＝（０，０）Ｔ）から衝突する単位エネルギのテスト信号を用いてレンダリング行列を正規化することである：
【数１１】
0006271586.tif 000014
モードベクトルｙは、方向Ωｃの球面調和関数から成るＨＯＡテスト信号であることができる。
【００３０】
リファレンスＨＯＡコンテンツは、エネルギ及びラウドネス測定による実際の検討において上記の検討を評価するために使用されてきた。ラウドネス測定は、ＨＯＡ０次の係数チャネル（Ｗ−チャネル）について、レンダリング後のマルチチャネル信号（５よりも多いチャネルを有するセットアップのための１のチャネル重みによる）について、並びにテクニカラーの聴取部屋の両耳室内インパルス応答（ＢＲＩＰ；Binaural  Room  Impulse  Response）によるマルチチャネル信号の両耳ダウンミックス及び中央のスピーカチャネルからの−２３ＬＫＦＳテスト信号による較正について、行われてきた。エネルギ測定（フルファイル）は、無音パッセージを伴うアイテムを除いて、ラウドネス測定と極めて密に整合する。それらのパッセージはＬＫＦＳ測定について無視されるので、値はここではより高い。これは、エネルギ測定が知覚ラウドネスの推定を与えるとの仮定を証明する。測定データは以下で与えられる。データの幾つかの評価は、以下で提示される。それらは、非公式の聴取テストによって確かめられる。ＬＫＦＳ測定精度は±２ｄＢである。
【００３１】
以下は、異なるＨＯＡ次数についてレンダリングするときのラウドネス維持に言及する。図５は、２２．０ラウドスピーカ・チャネル・セットアップのための異なるＨＯＡ次数についてレンダリングするときのラウドネス維持を示す。３つのテストアイテムＴ１，Ｔ２，Ｔ３の両耳ラウドネス測定は、発明に従う方法において、レンダリング次数Ｎ＝１，・・，６について２２．０チャネル・セットアップへレンダリングされる。テストアイテムの次数はＮ＝６であり、レンダラは、より低い次数についてレンダリングするときに、より高い次数の成分を無視する。第１のテストアイテムＴ１についての測定は、次数Ｎ＝１についてのｔ１１、次数Ｎ＝２についてのｔ１２等であり、第２及び第３のテストアイテムＴ２，Ｔ３についても同様である。明らかなように、異なるレンダリング次数を有する部分信号（すなわち、夫々のテスト信号の成分）は、それらのレンダリング次数と無関係に、実質的に同じラウドネスによりレンダリングされる。
【００３２】
以下は、異なるラウドスピーカ・セットアップへレンダリングするときのレンダリング維持に言及する。図６は、Ｎ＝４（すなわち、４次）を有する１４個のＨＯＡテストアイテムＬＴ１，・・，ＬＴ１４を５つの異なるラウドスピーカ・セットアップへレンダリングするときのラウドネス維持を示す。５つの列からなる各グループ、すなわち、各テストアイテムにおいて、列（左から右）は、２２．０（ＣｆＰ）、１２．０（Ａｕｒｏ３Ｄ）、９．０（Ａｕｒｏ３Ｄ）、５．０（ＩＴＵ）、２．０（±３０°）の選択されたラウドスピーカ・セットアップに対応する。明らかなように、５つの列からなる各グループにおける差は極めて小さい。すなわち、所与の信号は、ラウドスピーカ・セットアップ及びラウドスピーカの数と無関係に、実質的に同じラウドネスによりレンダリングされる。理論上のエネルギ保存は、スピーカの数が少なくなるにつれて劣化するので（特に、２Ｄセットアップ）、式（７）に従うレンダリング行列の正規化は、実際的に有用な改善されたレンダラをもたらす。
【００３３】
以下は、レンダリング後の１次ＨＯＡのＷ−チャネル及びアイテムのラウドネスに言及する。Ｎ３Ｄ正規化を用いると、０次の係数チャネルは、略等しいエネルギ及びラウドネスのレンダリングされた出力の調整されていないモノバージョンを含む。図７は、テクニカラーからの１３個のＮ＝４テストアイテムＬＭ１，・・，ＬＭ１３についてのこの相関、すなわち、レンダリング後の平均ラウドネスと比較されるＨＯＡ  Ｗ−チャネル信号のラウドネス測定を示す。モノチャネルＬＫＦＳ測定（夫々のテストアイテムの左手の列）は、図６の５つの異なるスピーカ・セットアップの両耳ＬＫＦＳ測定の平均（夫々のテストアイテムの右手の列）と比較される。驚くべき所見は、両方が実質的に等しいことである。よって、ＨＯＡ  Ｗ−チャネルは、レンダリング後の全体的なラウドネスの推定として使用され得る。この特徴は、マイクロホン・アレイを使用するときに、又はコンテンツ生成内でラウドネスをモニタするときに、ライブＨＯＡ放送録音のためにラウドネスレベルを設定するのに使用され得る。
【００３４】
発明に従って、知覚ラウドネスは、生成側でのミキシング設備及び消費者側での提示現場が同じラウドスピーカ・レベル較正、及び特別なエネルギ保存特性を備えたレンダラを使用する場合に、正規化され得る。エネルギ保存は、任意の方向からの指向性ソース信号のエネルギを保存するレンダラの特性を記述する。これは、全ての２Ｄオーディオ出力フォーマット、中でも少なくともオーディオ・オブジェクト、チャネルベース及びＨＯＡに当てはまる。
【００３５】
ＨＯＡコンテンツは、ミキシング設備において、例えば聴取によって、較正されたラウドネスレベルであり、エネルギ保存レンダリングはモニタリングのために使用される。適切なマスタゲインがＨＯＡコンテンツ（全ての係数チャネル）に適用され、次いで、ＨＯＡコンテンツは、分配のために記憶される。更に、エネルギ保存レンダリングは、提示現場でも使用される。
【００３６】
実際に、ＨＯＡレンダは、少なくともラウドスピーカ間のエリア（すなわち、ラウドスピーカ間の位置）においてエネルギ保存するよう構成され得る。良好なラウドネス維持結果をもたらすレンダリング行列の正規化は、次の式による正規化である：
【数１２】
0006271586.tif 000015
ＨＯＡ０次の係数チャネル（Ｗ−チャネル）のＬＫＦＳ測定は、その場合に、関連するチャネル表現のラウドネス測定又はそれらの関連する両耳ダウンミックスとおおよそ整合する。よって、それは、コンテンツ生成内で、特に、放送及び特にライブ放送内で、ラウドネスをモニタするために使用され得る。しかし、それはまた、記録されたコンテンツをモニタするために有利であり得る。
【００３７】
本発明は、一実施形態において、高次アンビソニックス（ＨＯＡ）のラウドネスレベルを調整する方法を提供し、それにより、音圧レベルを較正された３Ｄオーディオシステムでのレンダリング及び再生は、異なる現場において及び／又は異なるスピーカ・セットアップにより聴く場合に、等しいラウドネスにより知覚され得る。コンテンツのＨＯＡ符号化は特定される必要があり、レンダラはエネルギ保存の共通した特性を共有する必要がある。
【００３８】
一実施形態において、発明は、例えば、ＨＯＡマイクロホン録音により生成されるＨＯＡ信号、ＡＯ信号又はチャネルベース信号のライブ放送のために使用され得るＨＯＡコンテンツラウドネスをモニタするツールに関する。
【００３９】
以下で、ライブ放送のためのＨＯＡシステムが記載される。ＨＯＡマイクロホン・アレイは、スタジオでのポストプロセッシングなしで、実時間において３Ｄオーディオを捕捉することを可能にする。よって、そのようなコンテンツ捕捉システムは、例えば、スポーツ及び／又はコンサートイベントなどのようなライブ放送イベントに適している。以下は、３Ｄコンテンツのラウドネスを設定しモニタするシステムを記載する。
【００４０】
図８は、ＨＯＡシステム８０の一実施形態を表す。最初に、マイクロホン（又はマイクロホン・アレイ）８１は、記録すべきイベントまでの臨界距離内に位置付けられ、それにより、アナログ・セットアップ（例えば、適切な前置増幅）が保たれる。マイクロホン８１は、その出力を実時間ＨＯＡ符号器８２へ供給する。符号器８２は、符号化されたＨＯＡフォーマット信号を生成する。ＨＯＡフォーマット信号は、ＨＯＡラウドネスメータ８３、すなわち、ＨＯＡ０次チャネル（ＨＯＡ  Ｗ−チャネル）を評価し、ＬＫＦＳ（ＬＵＦＳ）ラウドネス測定を計算するユニットへ供給される。Ｗ−チャネルは、ＨＯＡ符号器８２、ＨＯＡラウドネスメータ８３又はＷ−チャネル抽出ユニット（図示せず。）内でＨＯＡ信号から取り出され得る。ＨＯＡラウドネスメータ８３は、瞬時的なラウドネス測定、短期のラウドネス測定、及び積分されたラウドネス測定ｑ８３［１］，［２］をＨＯＡマスタゲインユニット８４へ供給する。一実施形態において、ラウドネス測定はまた、手動により又は自動的に、ＨＯＡマスタゲインユニット８４へ供給されるより前に、任意の適応／選択ユニット８３ｘにおいて適応及び／又は選択されてよい。一実施形態において、ＨＯＡラウドネスメータ８３は、関連するエネルギ測定を実行する。このエネルギ測定は、ＨＯＡマスタゲインユニット８４においてゲインを設定するために且つラウドネスを調整するために使用される。
【００４１】
任意の第２のＨＯＡラウドネス計測ユニット８３ｍは、モニタリングの目的で、特に、より複雑なシステムのために、ＨＯＡマスタゲインユニット８４の後に取り付けられてよく、任意の第１のＨＯＡラウドネス維持レンダラ８５が、モニタリングのために使用され得る。良好な実際的な解決法は、初期マスタゲインを設定するためにスピーチ・アンカー・テスト記録を使用することである。
【００４２】
消費者装置における送信器による送信８６（又は記憶ユニットによる記憶）及び受信器による受信（又は検索ユニットによる記憶されたコンテンツの読み出し）の後、第２のＨＯＡラウドネス維持レンダラ８７は、オーディオ／コンテンツを所望のラウドスピーカ・セットアップ８８へレンダリングするために使用される。例えば、ステレオｓｔ、５．１サラウンド・サウンド５．１ｓｕｒ、他の３Ｄ音場ｓｆ又はヘッドホンｈｐなどの、異なるラウドスピーカ・セットアップ８８が使用されてよい。ＨＯＡフォーマットの１つの利点は、如何なる特定のラウドスピーカ・セットアップへの高品質でのそのレンダリングも最適化することが比較的簡単であることである。
【００４３】
原理上は、ＨＯＡシステム８０は、少なくとも２つの部分、すなわち、符号化部分及び復号化部分を含むものとして理解され得る。符号化部分は、符号器ユニット８２、ラウドネスメータ８３、マスタゲインユニット８４並びに、任意に、適応／選択ユニット８３ｘ、ラウドネス計測ユニット８３ｍ及び第１のラウドネス維持レンダラ８５のうちの１つ以上を含む。受信器における復号化部分は、少なくともラウドネス維持レンダリングユニット８７を含む。符号化部分及び復号化部分が記憶装置を通じて接続される場合は、その記憶装置は、それらのうちのいずれか一方の部分であってよく、あるいは、別個であってよい。
【００４４】
図９は、一実施形態において、チャネルベースの送信又は再生のためのシステム９０を表す。ＨＯＡ信号は、マイクロホン／マイクロホン・アレイ９１、ＨＯＡ符号器９２、ＨＯＡマスタゲインユニット９４、ＨＯＡラウドネス計測器９３及びＨＯＡラウドネス維持レンダラ９５を用いて、上述されたように捕捉されてラウドネス調整される。任意の第２のＨＯＡラウドネス計測ユニット９３ｍは、モニタリングのために、ＨＯＡマスタゲインユニット９４の後に取り付けられてよい。ＨＯＡラウドネス維持レンダラ９５は、この実施形態ではもはや任意でない。一実施形態において、ラウドネス測定はまた、手動により又は自動的に、ＨＯＡマスタゲインユニット９４へ供給されるより前に、任意の適応／選択ユニット９３ｘにおいて適応及び／又は選択されてよい。
【００４５】
ＨＯＡラウドネス維持レンダラ９５によって供給されるチャネルベースのオーディオは、次いで、送信ユニット又は記憶ユニット９６において送信又は記憶される。送信セットアップのスピーカ位置及び／又は元のコンテンツフォーマットを示すフラグは、メタデータとして信号に付加される。よって、送信セットアップのスピーカ位置、及び／又は元のコンテンツがＨＯＡであったことを伝えるフラグは、送信又は記憶される信号においてメタデータとして含まれる。送信／記憶及び消費者側でのコンテンツの受信の後、エネルギ保存チャネルベースレンダラ９９は、一実施形態において、消費者のラウドスピーカ・セットアップ９８へデータをレンダリングするために使用される。他の実施形態（図示せず。）では、エネルギ保存チャネルベースレンダラ９９は、送信又は記憶の前に、所定の又は個々のラウドスピーカ・セットアップ９８へデータをレンダリングするために使用される。何れの場合にも、チャネルベースレンダラ９９は、式（４）に従って作動する。異なるラウドスピーカ・セットアップが、上述されたように可能である。
【００４６】
代替的に、チャネルベースのコンテンツは、第２のＨＯＡ符号化ユニット（図示せず。例えば、チャネルベースレンダラ９９の代わり）において、送信されたメタデータを用いてＨＯＡへ再符号化され得る。このとき、第２のＨＯＡ符号化ユニットは、好適には、ラウドネス維持ＨＯＡレンダラである。
【００４７】
これより、ＨＯＡラウドネスメータ８３、９３が記載される。入力信号として、Ｎ３Ｄにより正規化された０次ＨＯＡ係数チャネルが使用される。処理は［１］，［２］と同じであり、処理はモノチャネルについてのみ記載される。様々な異なる実施形態において、ＨＯＡラウドネスメータ８３、９３は、瞬時的なラウドネス測定、短期のラウドネス測定及び積分されたラウドネス測定のうちの１つ以上を示すことができる。それらは、ラウドネス測定のために使用される時間フレームによって主として異なる。
【００４８】
ラウドネス維持レンダラが使用される条件下で、測定は、消費者のラウドスピーカ・セットアップでレンダリングした後に得られる（目標）ラウドネスに関する推測結果を導き出すことを可能にする。例えば、測定は、異なるコンテンツ又は異なるプログラムを切り替えるときにラウドネスのジャンプを伴わずに優れた聴取レベルを保つために使用され得る。音圧レベル（ＳＰＬ）較正システムによれば、等しいラウドネスが実現され得る。
【００４９】
ラウドネス維持レンダラは、次の要件を（少なくとも近似的に）満足することを必要とする：
・ＡＯベースのコンテンツについて：式（３）；
・チャネルベースのコンテンツについて：式（４）；及び
・ＨＯＡベースのコンテンツについて：式（５）。
【００５０】
一実施形態において、ＨＯＡレンダラは、式（６）、（７）、（８）に従うレンダリング行列正規化を使用する。次のエネルギ計算は、ＨＯＡレンダリング後の信号エネルギＥＷを導出する。テスト信号は、ＡのＲＭＳ値及び方向Ω＝（θ，φ）Ｔから衝突するエネルギＥＡ＝Ａ２を持ったＤＣ成分を含まない時間信号ｘ（ｔ）であり、θは仰角、φは方位角である。ｘがτ個のサンプルのベクトルを表すとする。ｘ＝（ｘ（ｔ１），・・・ｘ（ｔｉ），・・・ｘ（ｔτ））Ｔ。
【数１３】
0006271586.tif 000016
ＨＯＡ符号化は、次の式によって与えられる：
【数１４】
0006271586.tif 000017
なお、ベクトルｙ＝（Ｙ００＊（Ω），Ｙ１−１＊（Ω），・・）は、方向Ωの実数値の球面調和関数から作られる。次数のインデックスｏ及び時間サンプルのインデックスｉとして、行列要素Ｂｏ，ｉによれば：
【数１５】
0006271586.tif 000018
ＨＯＡ領域におけるエネルギは、行列Ｂの二乗されたフロベニウス・ノルムによって与えられる：
【数１６】
0006271586.tif 000019
そして、上記の符号化式によれば、次のようになる：
【数１７】
0006271586.tif 000020
ラウドスピーカ信号Ｗは、次の行列：
【数１８】
0006271586.tif 000021
によるレンダリング後に、次の式によって与えられる：
【数１９】
0006271586.tif 000022
なお：
【数２０】
0006271586.tif 000023
である。Ｗは、Ｂの列ベクトルｂｉから、次の式によって計算される：
【数２１】
0006271586.tif 000024
スピーカ信号のエネルギは、行列Ｗの二乗されたフロベニウス・ノルムによって与えられる：
【数２２】
0006271586.tif 000025
これは、列の和及び列ｂｉについてのユークリッド・ベクトル・ノルムを用いて書き換えられる：
【数２３】
0006271586.tif 000026
そして、
【数２４】
0006271586.tif 000027
によれば、次のとおりである：
【数２５】
0006271586.tif 000028
復号化行列の正規化は、以下において説明される。すなわち、如何にして式（７）及び、式（７）の特別な場合としての式（６）が導出されるかが説明される。コーシー・シュワルツの不等式は、フロベニウス・ノルムを用いて任意の行列について定義される：
【数２６】
0006271586.tif 000029
これから、レンダリング行列Ｄのための簡単な正規化は、次に従う：
【数２７】
0006271586.tif 000030
完璧なエネルギ保存は、ラウドスピーカの数がＨＯＡ係数の数よりも多いか又は等しい場合にのみ起こることができる（Ｌ≦（Ｎ＋１）２）。レンダリング行列Ｄが完璧にエネルギ保存している場合に、その特異値分解はＤ＝ＵＩＶになる。ここで、Ｕ，Ｖはユニタリ行列であり、Ｉはサイズ（Ｎ＋１）２×（Ｎ＋１）２の単位行列である。フロベニウス・ノルムは、ユニタリ変換の下で不変であり、この場合には、次のとおりである：
【数２８】
0006271586.tif 000031
これは式（６）を直接導く。
【００５１】
本発明の好適な実施形態に適用される本発明の基本的な新規の特徴が図示、記載及び指摘されてきたが、記載されている装置及び方法における、開示されている装置の形態及び詳細における、並びにそれらの動作における様々な省略及び置換は、本発明の主旨から逸脱することなしに当業者によって行われ得ることが理解されるであろう。同じ結果を達成するよう実質的に同じ方法において実質的に同じ機能を実行するそれらの要素の全ての組み合わせは、発明の適用範囲内にあることが明示的に意図される。１つの記載されている実施形態から他への要素の置換も十分に意図され且つ考えられている。本発明は単に一例として記載されているのであり、詳細の変更は発明の適用範囲から逸脱することなしに行われ得ることが理解されるであろう。明細書並びに（必要に応じて）特許請求の範囲及び図面において開示されている夫々の特徴は、独立して、又は何らかの適切な組み合わせにおいて、提供されてよい。特徴は、必要に応じて、ハードウェア、ソフトウェア、又は両者の組み合わせにおいて実装されてよい。接続は、適用可能である場合に、無線接続又は有線の、必ずしも直接又は専用でない接続として、実装されてよい。
【００５２】
［引用参考文献］
［１］ISO/IEC  JTC1/SC29/WG1  1/N13196,  ‘Calibration  of  22.2  multichannel  sound  reproduction’,  Audio  Subgroup,  October  2012,  Shanghai,  China
［２］EBU  Technical  Recommendation  R128,  “Loudness  Normalization  and  Permitted  Maximum  Level  of  Audio  Signals”,  Geneva,  2010,  [http://tech.ebu.ch/docs/r/r128.pdf]
［３］ITU-R  Recommendation  BS.1770-2,  “Algorithms  to  measure  audio  programme  loudness  and  true-peak  audio  level”,  Geneva,  2011
［４］ATSC  A/85,  “Techniques  for  Establishing  and  Maintaining  Audio  Loudness  for  Digital  Television”,  Advanced  Television  Systems  Committee,  Washington,  D.C.,  July  25,  2011
［５］F.  Zotter,  H.  Pomberger,  M.  Noisternig,  “Energy-Preserving  Ambisonic  Decoding”,  ACTA  ACUSTICA  UNITED  with  ACUSTICA,  Vol.98  (2012),  pages  37-47
［６］F.  Zotter,  M.  Frank,  “All-Round  Ambisonic  Panning  and  Decoding”,  J.  Audio  End.  Soc,  Vol.60,  No.10,  2012  October
【図１】
0006271586.tif 000049
【図２】
0006271586.tif 000050
【図３】
0006271586.tif 000051
【図４ａ）】
0006271586.tif 000052
【図４ｂ）】
0006271586.tif 000053
【図４ｃ）】
0006271586.tif 000054
【図５】
0006271586.tif 000055
【図６】
0006271586.tif 000056
【図７】
0006271586.tif 000057
【図８】
0006271586.tif 000058
【図９】
0006271586.tif 000059