(57)【要約】  チャネル間相関算出部（１０２）は、ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出する。ＤＭＡステレオ符号化部（１０４）とＤＭステレオ符号化部（１０５）は、チャネル間相関が閾値より大きい場合に共通の符号化モードを用いて、左チャネル信号及び右チャネル信号をそれぞれ符号化し、チャネル間相関が閾値以下の場合に左チャネル信号及び右チャネル信号に対して個別に判定された符号化モードを用いて、左チャネル信号及び右チャネル信号をそれぞれ符号化する。
2018221138.tif 000002
【特許請求の範囲】
【請求項１】
ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出する算出回路と、
前記チャネル間相関が閾値より大きい場合に、共通の符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化し、
前記チャネル間相関が前記閾値以下の場合に、前記左チャネル信号及び前記右チャネル信号に対して個別に判定された符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化する符号化回路と、
を具備する符号化装置。
【請求項２】
前記符号化回路は、左チャネルと右チャネルについて主要チャネルと非主要チャネルとを特定し、前記主要チャネルの符号化モードを判定するための第１のパラメータと、前記非主要チャネルの符号化モードを判定するための第２のパラメータとに対して重み付け加算を行い、前記重み付け加算によって得られる重み付けパラメータに基づいて前記共通の符号化モードを選択する、
請求項１に記載の符号化装置。
【請求項３】
前記第１のパラメータに対する第１の重み係数は、前記第２のパラメータに対する第２の重み係数より大きく、
前記チャネル間相関が小さいほど、前記第１の重み係数は大きい、
請求項２に記載の符号化装置。
【請求項４】
前記第１のパラメータに対する第１の重み係数は、前記第２のパラメータに対する第２の重み係数より大きく、
前記左チャネル信号と前記右チャネル信号との間のエネルギ差が大きいほど、前記第１の重み係数は大きい、
請求項２に記載の符号化装置。
【請求項５】
前記符号化回路は、現フレームで選択された前記共通の符号化モードが、過去のフレームで選択された前記共通の符号化モード、現フレームの前記第１のパラメータに基づいて判定される符号化モードと異なり、かつ、現フレームの前記第２のパラメータに基づいて判定される符号化モードの何れかと同一である場合、現フレームの前記共通の符号化モードを再選択する、
請求項２に記載の符号化装置。
【請求項６】
前記符号化回路は、現フレームの前記重み付けパラメータと、過去フレームの前記重み付けパラメータとを用いてスムージング処理を行い、前記スムージング処理後の重み付けパラメータに基づいて前記共通の符号化モードを再選択する、
請求項５に記載の符号化装置。
【請求項７】
前記符号化回路は、更に、前記チャネル間相関が、前記閾値よりも大きい第２の閾値よりも大きい場合、前記左チャネル信号及び前記右チャネル信号に対して、Mid/Sideステレオ符号化を行う、
請求項１に記載の符号化装置。
【請求項８】
前記算出回路は、前記左チャネル信号及び前記右チャネル信号のうち、一部の帯域の周波数スペクトルパラメータを用いて、前記チャネル間相関を算出する、
請求項１に記載の符号化装置。
【請求項９】
ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出するステップと、
前記チャネル間相関が閾値より大きい場合に、共通の符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化し、前記チャネル間相関が前記閾値以下の場合に、前記左チャネル信号及び前記右チャネル信号に対して個別に判定された符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化するステップと、
符号化方法。
【請求項１０】
前記符号化するステップにおいて、左チャネルと右チャネルについて主要チャネルと非主要チャネルとを特定し、前記主要チャネルの符号化モードを判定するための第１のパラメータと、前記非主要チャネルの符号化モードを判定するための第２のパラメータとに対して重み付け加算を行い、前記重み付け加算によって得られる重み付けパラメータに基づいて前記共通の符号化モードを選択する、
請求項９に記載の符号化方法。
【請求項１１】
前記第１のパラメータに対する第１の重み係数は、前記第２のパラメータに対する第２の重み係数より大きく、
前記チャネル間相関が小さいほど、前記第１の重み係数は大きい、
請求項１０に記載の符号化方法。
【請求項１２】
前記第１のパラメータに対する第１の重み係数は、前記第２のパラメータに対する第２の重み係数より大きく、
前記左チャネル信号と前記右チャネル信号との間のエネルギ差が大きいほど、前記第１の重み係数は大きい、
請求項１０に記載の符号化方法。
【請求項１３】
前記符号化するステップにおいて、現フレームで選択された前記共通の符号化モードが、過去のフレームで選択された前記共通の符号化モード、現フレームの前記第１のパラメータに基づいて判定される符号化モードと異なり、かつ、現フレームの前記第２のパラメータに基づいて判定される符号化モードの何れかと同一である場合、現フレームの前記共通の符号化モードを再選択する、
請求項１０に記載の符号化方法。
【請求項１４】
前記符号化するステップにおいて、現フレームの前記重み付けパラメータと、過去フレームの前記重み付けパラメータとを用いてスムージング処理を行い、前記スムージング処理後の重み付けパラメータに基づいて前記共通の符号化モードを再選択する、
請求項１３に記載の符号化方法。
【請求項１５】
前記符号化するステップにおいて、更に、前記チャネル間相関が、前記閾値よりも大きい第２の閾値よりも大きい場合、前記左チャネル信号及び前記右チャネル信号に対して、Mid/Sideステレオ符号化を行う、
請求項９に記載の符号化方法。
【請求項１６】
前記算出するステップにおいて、前記左チャネル信号及び前記右チャネル信号のうち、一部の帯域の周波数スペクトルパラメータを用いて、前記チャネル間相関を算出する、
請求項９に記載の符号化方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、符号化装置及び符号化方法に関する。
【背景技術】
【０００２】
近年、3GPP（3rd Generation Partnership Project）において、EVS（Enhanced Voice Services）コーデックが標準化された（例えば、非特許文献１を参照）。EVSコーデックは、モノラル音声音響信号を符号化するために設計されている。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】3GPP TS 26.445 V14.0.0, "Codec for Enhanced Voice services (EVS); Detailed algorithmic description (Release 14)", 2017-03
【非特許文献２】J.D.Johnston, A.J.Ferreira, “SUM-DIFFERENCE STEREO TRANSFORM CODING,” proc. IEEE ICASSP1992, pp.II-560 - II-572, 1992
【非特許文献３】E.Schuijers, W.Oomen, B.Brinker, and J. Breebaart, “Advances in Parametric Coding for High-Quality Audio”, in Preprint 5852, 114th AES convention, Amsterdam, Mar.2003.
【発明の概要】
【０００４】
EVSコーデックはステレオ信号の入出力をサポートしていないが、EVSコーデックのモノラル符号化を用いて、ステレオ信号の左チャネル、右チャネルをそれぞれ処理すれば、ステレオレンダリングシステムでも利用可能である。しかしながら、EVSコーデックのように多くの符号化モードを切り替えて符号化するマルチモードモノラルコーデックを用いてステレオ信号を符号化した場合、ステレオ信号の左チャネルと右チャネルとで異なる符号化モードを用いて符号化され、ステレオ再生時の音声品質を劣化させる恐れがある。なお、ステレオ信号のＬチャネル信号とＲチャネル信号とに分けて別々にモノラル符号化することを、「デュアルモノ符号化」と呼ぶこともある。
【０００５】
本開示の一態様は、マルチモードコーデックを用いてステレオ信号を符号化する場合でも、ステレオ再生時の音声品質の劣化を抑えることができる符号化装置及び符号化方法の提供に資する。
【０００６】
本開示の一態様に係る符号化装置は、ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出する算出回路と、前記チャネル間相関が閾値より大きい場合に共通の符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化し、前記チャネル間相関が前記閾値以下の場合に前記左チャネル信号及び前記右チャネル信号に対して個別に判定された符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化する符号化回路と、を具備する構成を採る。
【０００７】
本開示の一態様に係る符号化方法は、ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出し、前記チャネル間相関が閾値より大きい場合に共通の符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化し、前記チャネル間相関が前記閾値以下の場合に前記左チャネル信号及び前記右チャネル信号に対して個別に判定された符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化する。
【０００８】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラム、または、記録媒体で実現されてもよく、システム、装置、方法、集積回路、コンピュータプログラムおよび記録媒体の任意な組み合わせで実現されてもよい。
【０００９】
本開示の一態様によれば、マルチモードコーデックを用いてステレオ信号を符号化する場合でも、ステレオ再生時の音声品質の劣化を抑えることができる。
【００１０】
本開示の一態様における更なる利点および効果は、明細書および図面から明らかにされる。かかる利点および／または効果は、いくつかの実施形態並びに明細書および図面に記載された特徴によってそれぞれ提供されるが、１つまたはそれ以上の同一の特徴を得るために必ずしも全てが提供される必要はない。
【図面の簡単な説明】
【００１１】
【図１】EVSコーデックの一例を示す図
【図２】信号の分析パラメータと符号化モードとの対応関係の一例を示す図
【図３】デュアルモノ符号化の構成例を示す図
【図４】実施の形態１に係る符号化装置の一部の構成例を示すブロック図
【図５】実施の形態１に係る符号化装置の構成例を示すブロック図
【図６】実施の形態１に係る信号分析部及びＤＭＡステレオ符号化部の構成例を示すブロック図
【図７】実施の形態１に係る符号化モード選択処理の流れを示すフロー図
【図８】実施の形態１の変形例に係る符号化モード選択処理の流れを示すフロー図
【図９】実施の形態１の変形例に係る重み係数の選択処理の流れを示すフロー図
【図１０】実施の形態１の変形例に係るチャネル間エネルギ差と重み係数との対応関係の一例を示す図
【図１１】実施の形態２に係る信号分析部及びＤＭＡステレオ符号化部の構成例を示すブロック図
【図１２】実施の形態２に係る符号化モードの判定訂正処理の流れを示すフロー図
【図１３】実施の形態３に係る符号化装置の構成例を示すブロック図
【図１４】実施の形態３に係るチャネル間相関値の範囲と符号化モードとの対応関係の一例を示す図
【図１５】実施の形態４に係る信号分析部及びチャネル間相関算出部の構成例を示すブロック図
【図１６】実施の形態４に係る信号分析部及びチャネル間相関算出部の動作例を示す図
【図１７】実施の形態４の変形例２に係る信号分析部及びチャネル間相関算出部の構成例を示すブロック図
【発明を実施するための形態】
【００１２】
以下、本開示の実施の形態について図面を参照して詳細に説明する。
【００１３】
まず、マルチモードモノラル符号化システムの一例として，3GPP EVS符号化システムについて概説する（例えば、非特許文献１を参照）。
【００１４】
EVSコーデックでは、非特許文献１に記載されているように、複数の符号化技術（符号化モード）が採用されている（例えば、図１を参照）。EVSコーデックに採用された複数の符号化技術は、基本的に、以下の二つの原理に基づく。一つは線形予測（Linear Prediction：LP）ベースのアプローチであり、もう一つは周波数領域アプローチである。線形予測ベースの符号化では、CELP（Code Excited Linear Prediction）符号化技術に基づいて各ビットレート専用に最適化された符号化モード（例えば、ACELP（Algebraic CELP）等）が用いられる。また、周波数領域アプローチでは、HQ MDCT（High Quality Modified Discrete Cosine Transform）技術又はTCX（Transformed Code Excitation）技術などが採用されている。
【００１５】
EVSコーデックでは、入力された音声・音響信号に応じて、例えば、ACELP、HQ MDCT及びTCXの中から最も適した符号化モードが選択される。各符号化モードは各種信号を効率的に符号化できるように設計、調整されている。EVSコーデックでの符号化モード選択は、例えば、ビットレート、オーディオ信号の帯域幅、音声/音楽分類、選択された符号化モード、又はその他のパラメータ（特徴量）に基づいて行われる。図２は、一例として、ビットレート（[kbps]）、帯域幅（SWB（super wideband）、FB（fullband））、入力信号の種類（speech/audio）を示すパラメータと、各パラメータに応じて選択される符号化モード（ACELP、GSC、TCX、HQ MDCT）との対応関係を示す。
【００１６】
上述したように、EVSコーデックはモノラルコーデックだが、モノラルコーデックを用いてステレオ信号の各チャネルをそれぞれ処理すれば、ステレオレンダリングシステムでも利用可能である。図３は、一例として、ステレオ信号の各チャネル（左チャネル、右チャネル）の各々に対してモノラルコーデックを用いて処理するデュアルモノ符号化（dual mono encoder）の構成例を示す。
【００１７】
図３に示すように、ステレオ信号の左チャネル信号（以下、「Ｌ信号」と呼ぶ）及び右チャネル信号（以下、「Ｒ信号」と呼ぶ）は、モノラルコーデックによって個別に符号化される。この場合、ステレオ信号の左チャネルと右チャネルとで異なる符号化モードが選択され、符号化されることがある。具体的には、Ｌ信号及びＲ信号の特徴は、チャネル間の信号類似度によって変わってくるため、両方のチャネル信号がEVSコーデックのようなマルチモードコーデックで別々に処理された場合、両方のチャネルで異なる符号化モードがそれぞれ選択される場合が発生する。両方のチャネルで異なる符号化モードが選択されると、復号信号の主観品質が劣化し、ステレオ再生時に異音及び／又は歪となって聞こえたり、ステレオ定位感が乱れたりする原因となる場合がある。
【００１８】
そこで、本開示の各実施の形態では、多くの符号化モードを切り替えて符号化処理を行うマルチモードコーデックによりステレオ信号の両方のチャネル信号が別々に処理される場合でも、ステレオ再生時の音声品質の劣化（異音及び／又は歪み、定位感の乱れの発生）を抑える方法について説明する。
【００１９】
（実施の形態１）
［通信システムの概要］
本実施の形態に係る通信システムは、符号化装置（encoder）１００及び復号装置（decoder）（図示せず）を備える。
【００２０】
図４は、本実施の形態に係る符号化装置１００の一部の構成を示すブロック図である。図４に示す符号化装置１００において、チャネル間相関算出部１０２は、ステレオ信号を構成する左チャネル信号（Ｌ信号）と右チャネル信号（Ｒ信号）を用いて、左チャネルと右チャネルとの間のチャネル間相関（相互相関係数：Correlation Coefficient）を算出する。符号化部（ＤＭＡステレオ符号化部１０４及びＤＭステレオ符号化部１０５）は、チャネル間相関が閾値より大きい場合に共通の符号化モードを用いて、左チャネル信号及び右チャネル信号をそれぞれ符号化し、チャネル間相関が閾値以下の場合に左チャネル信号及び右チャネル信号に対して個別に判定された符号化モードを用いて、左チャネル信号及び右チャネル信号をそれぞれ符号化する。
【００２１】
［符号化装置の構成］
図５は、本実施の形態に係る符号化装置１００の構成例を示すブロック図である。図５において、符号化装置１００は、信号分析部１０１と、チャネル間相関算出部１０２と、切替スイッチ１０３と、ＤＭＡ（Dual Mono with mode alignment）ステレオ符号化部１０４と、ＤＭ（Dual Mono）ステレオ符号化部１０５と、多重化部１０６と、を含む構成を採る。
【００２２】
図５において、信号分析部１０１、チャネル間相関算出部１０２及び切替スイッチ１０３には、ステレオ信号を構成するＬ信号（Left channel）、及び、Ｒ信号（Right channel）が入力される。
【００２３】
信号分析部１０１は、入力されるＬ信号及びＲ信号に対して信号分析を行い、左チャネル及び右チャネルについて符号化モードの判定に必要なパラメータ（例えば、ビットレート、帯域幅、種類などの特徴量）をそれぞれ得る。信号分析部１０１は、得られた分析パラメータ（parameters）を切替スイッチ１０３に出力する。例えば、信号分析部１０１では、信号分析の際、チャネル信号の周波数領域変換処理、及び、エネルギ算出処理等が行われる。
【００２４】
チャネル間相関算出部１０２は、入力されるＬ信号及びＲ信号を用いて、例えば、次式（１）に従って、左チャネルと右チャネルとの間のチャネル間相関（相互相関係数）αを算出する。
【数１】
2018221138.tif 000003
【００２５】
式（１）において、R11及びR22は、Ｌ信号及びＲ信号のエネルギ（auto-correlation）を示す（例えば、R11がＬ信号に対応し、R22がＲ信号に対応）。また、R12は、Ｌ信号とＲ信号との間のクロススペクトルを示す。また、Framelengthはフレーム内の周波数スペクトルパラメータ（スペクトル係数）の数を示し、l(k)は、Ｌ信号におけるｋ番目のスペクトル係数を示し、R(k)は、Ｒ信号におけるｋ番目のスペクトル係数を示す。
【００２６】
また、チャネル間相関算出部１０２は、算出した相互相関係数αに基づいて、ステレオ信号（Ｌ信号及びＲ信号）に対するステレオ符号化モードを判定する。
【００２７】
ここで、ステレオ符号化モードには、例えば、図３に示すように、Ｌ信号及びＲ信号に対して符号化モードを個別に選択して符号化するモード（以下、「デュアルモノ符号化モード」又は「ＤＭステレオ符号化モード」と呼ぶ）、及び、後述するように、Ｌ信号及びＲ信号に対して共通の符号化モードを選択して符号化するモード（以下、「共通デュアルモノ符号化モード」又は「ＤＭＡステレオ符号化モード」と呼ぶ）がある。
【００２８】
具体的には、チャネル間相関算出部１０２は、相互相関係数αが閾値以下の場合にＤＭステレオ符号化モードと判定し、相互相関係数αが閾値より大きい場合にＤＭＡステレオ符号化モードと判定する。一例として、チャネル間相関算出部１０２は、相互相関係数αが０の場合（つまり、Ｌ信号とＲ信号とに相関が無い場合）にＤＭステレオ符号化モードと判定し、相互相関係数αが０より大きい場合（α＞０）にＤＭＡステレオ符号化モードと判定してもよい。
【００２９】
チャネル間相関算出部１０２は、相互相関係数α、ステレオ符号化モードの判定結果であるステレオモード判定フラグ（stereo mode decision）を、切替スイッチ１０３に出力する。
【００３０】
切替スイッチ１０３は、チャネル間相関算出部１０２から入力されるステレオモード判定フラグがＤＭＡステレオ符号化モードである場合、入力されるＬ信号、Ｒ信号、信号分析部１０１から入力される分析パラメータ、及び、相関算出部１０１から入力される相互相関係数αをＤＭＡステレオ符号化部１０４に出力する。一方、切替スイッチ１０３は、ステレオモード判定フラグがＤＭステレオ符号化モードである場合、Ｌ信号、Ｒ信号及び分析パラメータをＤＭステレオ符号化部１０５に出力する。
【００３１】
ＤＭＡステレオ符号化部１０４は、相互相関係数α、及び、分析パラメータを用いて、Ｌ信号及びＲ信号に対する共通の符号化モードを判定（選択）する。そして、ＤＭＡステレオ符号化部１０４は、判定した共通の符号化モードを用いて、Ｌ信号及びＲ信号をそれぞれ符号化し、生成された符号化ビットストリームを多重化部１０６へ出力する。なお、ＤＭＡステレオ符号化部１０４における符号化モードの選択方法の詳細については後述する。
【００３２】
ＤＭステレオ符号化部１０５は、分析パラメータを用いて、Ｌ信号及びＲ信号に対して個別に符号化モードを判定（選択）する。そして、ＤＭステレオ符号化部１０５は、判定した符号化モードを用いて、Ｌ信号及びＲ信号をそれぞれ符号化し、生成された符号化ビットストリームを多重化部１０６へ出力する（例えば、図３を参照）。
【００３３】
多重化部１０６は、ＤＭＡステレオ符号化部１０４又はＤＭステレオ符号化部１０５から入力される符号化ビットストリームを多重する。多重化されたビットストリームは、復号装置（図示せず）へ送信される。
【００３４】
なお、図５に示す符号化装置１００は、切替スイッチ１０３と、ＤＭＡステレオ符号化部１０４と、ＤＭステレオ符号化部１０５と、を備える代わりに、これらの構成部と同等の処理を行う符号化部を備える構成（図示せず）でもよい。すなわち、当該符号化部は、チャネル間相関算出部１０２からのチャネル間相関（相互相関係数α）に応じて、ステレオ符号化モード（ＤＭＡステレオ符号化又はＤＭステレオ符号化）を決定し、決定したステレオ符号化モードを用いてステレオ信号を構成するＬ信号及びＲ信号をそれぞれ符号化すればよい。
【００３５】
［ＤＭＡステレオ符号化部１０４の動作］
次に、ＤＭＡステレオ符号化部１０４における符号化モードの選択方法の詳細について説明する。
【００３６】
図６は、図５に示す信号分離部１０１及びＤＭＡステレオ符号化部１０４の構成を示すブロック図である。図６において、ＤＭＡステレオ符号化部１０４は、適応ミキシング部１４１と、符号化モード選択部１４２と、Lch符号化部１４３と、Rch符号化部１４４と、ビットストリーム生成部１４５と、を含む構成を採る。
【００３７】
図６に示すように、適応ミキシング部１４１には、信号分析部１０１（Lch信号分析部）においてＬ信号に対して信号分析を行って得られるLch分析パラメータ（Left channel parameters）が切替スイッチ１０３（図示せず）を介して入力される。同様に、図６に示すように、適応ミキシング部１４１には、信号分析部１０１（Rch信号分析部）においてＲ信号に対して信号分析を行って得られるRch分析パラメータ（Right channel parameters）が切替スイッチ１０３（図示せず）を介して入力される。
【００３８】
適応ミキシング部１４１は、チャネル間相関算出部１０２（図５を参照）から入力される相互相関係数αに基づいて、信号分析部１０１から入力されるLch分析パラメータ及びRch分析パラメータに対してミキシング（混合）を行い、ミキシング後の分析パラメータ（Mixed channel parameters）を符号化モード選択部１４２に出力する。換言すると、ミキシング後の分析パラメータは、Ｌ信号及びＲ信号に対する符号化モードの判定のための共通のパラメータ（特徴量）を表す。
【００３９】
符号化モード選択部１４２は、適応ミキシング部１４１から入力されるミキシング後の分析パラメータを用いて、Ｌ信号及びＲ信号の双方に共通して適用する符号化モードを選択する。符号化モード選択部１４２における符号化モードの選択方法は、ミキシング後の分析パラメータに応じて、例えば、図２で説明したEVSコーデック（モノラル符号化）における選択方法と同じ方法でもよい。符号化モード選択部１４２は、選択した符号化モードを示す符号化モード情報（coding mode decision）をLch符号化部１４３及びRch符号化部１４４に出力する。
【００４０】
Lch符号化部１４３は、符号化モード選択部１４２から入力される符号化モード情報に示される符号化モードを用いてＬ信号を符号化し、生成される符号化ビットストリームを、ビットストリーム生成部１４５へ出力する。
【００４１】
Rch符号化部１４４は、符号化モード選択部１４２から入力される符号化モード情報に示される符号化モードを用いてＲ信号を符号化し、生成される符号化ビットストリームを、ビットストリーム生成部１４５へ出力する。
【００４２】
ビットストリーム生成部１４５は、Lch符号化部１４３から入力される符号化ビットストリーム、及び、Rch符号化部１４４から入力される符号化ビットストリームを用いてステレオ符号化ビットストリームを生成し、多重化部１０６（図５を参照）へ出力する。
【００４３】
図７は、本実施の形態に係るＤＭＡステレオ符号化モードにおける符号化モードの選択処理の主な流れを示すフロー図である。
【００４４】
信号分析部１０１（Lch信号分析部及びRch信号分析部）は、Ｌ信号（左チャネル）及びＲ信号（右チャネル）のエネルギを算出する（ＳＴ１０１）。次に、適応ミキシング部１４１は、ＳＴ１０１で算出された各チャネルのエネルギを用いて、チャネル間エネルギ差Δを算出する（ＳＴ１０２）。
【００４５】
そして、適応ミキシング部１４１は、Ｌ信号（左チャネル）及びＲ信号（右チャネル）について、主要チャネル（dominant channel）と非主要チャネル（non-dominant channel）とを特定する（ＳＴ１０３）。
【００４６】
例えば、適応ミキシング部１４１は、ＳＴ１０２で算出したチャネル間エネルギ差Δに基づいて、主要チャネル及び非主要チャネルを特定してもよい。例えば、チャネル間エネルギ差Δを次式（２）で表す。
【数２】
2018221138.tif 000004
【００４７】
ここで、式（２）においてR11を左チャネルのエネルギとし、R22を右チャネルのエネルギとする場合、適応ミキシング部１４１は、チャネル間エネルギ差Δの正負に応じて主要チャネル及び非主要チャネルを特定する。具体的には、適応ミキシング部１４１は、エネルギ差Δが正の場合（Δ＞０。つまり、R11>R22）には左チャネルが主要チャネルであり、右チャネルが非主要チャネルであると特定する。一方、適応ミキシング部１４１は、エネルギ差Δが負の場合（Δ＜０。つまり、R11<R22）には左チャネルが非主要チャネルであり、右チャネルが主要チャネルであると特定する。なお、主要チャネル及び非主要チャネルの特定方法は上記方法に限定されるものではない。
【００４８】
次に、適応ミキシング部１４１は、相互相関係数αに基づいて、ＳＴ１０３で特定した主要チャネルの分析パラメータ及び非主要チャネルの分析パラメータに対する重み係数（ウェイト）を決定する（ＳＴ１０４）。そして、適応ミキシング部１４１は、主要チャネルの分析パラメータ及び非主要チャネルの分析パラメータに対して、ＳＴ１０４で決定した重み係数を用いて重み付け加算することにより、分析パラメータのミキシング（適応ミキシング）を行う（ＳＴ１０５）。
【００４９】
例えば、適応ミキシング部１４１は、次式（３）に従って分析パラメータのミキシング（重み付け加算）を行い、分析パラメータ（重み付けパラメータ）Mpを求める。
【数３】
2018221138.tif 000005
【００５０】
式（３）において、Dpは主要チャネルの符号化モードを判定するための分析パラメータを示し、NDpは非主要チャネルの符号化モードを判定するための分析パラメータを示す。また、W1は主要チャネルの分析パラメータに対する重み係数を示し、W2は非主要チャネルの分析パラメータに対する重み係数を示し、次式（４）で表される。
【数４】
2018221138.tif 000006
【００５１】
ただし、正規化相互相関係数（以下、単に「相互相関係数」と呼ぶ）αは、０＜α＜１である。
【００５２】
すなわち、重み係数W1の最小値は０．６となり、重み係数W2の最大値は０．４となる。これより、左チャネルと右チャネルとの間の相互相関係数αに依らず、重み係数W1は、重み係数W2より大きくなり、重み係数W1＞重み係数W2の関係となる。
【００５３】
つまり、適応ミキシング部１４１は、非主要チャネルの分析パラメータと比較して、主要チャネルの分析パラメータの重み係数を大きくして、分析パラメータMpを求める。これにより、重み付け加算によって得られる分析パラメータMpは、主要チャネルの分析パラメータがより強調された値となる。
【００５４】
また、左チャネルと右チャネルとの間のチャネル間相関を示す相互相関係数αが小さいほど、主要チャネルの分析パラメータに対する重み係数W1は大きくなり、非主要チャネルの分析パラメータに対する重み係数W2は小さくなる。
【００５５】
すなわち、式（４）に示す例では、常に主要チャネル側に大きな重み付けがなされることを保証しつつ、チャネル間相関（相互相関係数α）が高くなると、両チャネルの重み付けが均等に近づく。つまり、チャネル間相関が高い場合、両チャネルで算出される分析パラメータが類似するので、主要チャネルを特に強調する必要がないので、両チャネルの重み付けが均等に近づくような重み付けがなされる。一方、チャネル間相関が低い場合、両チャネルで算出される分析パラメータの差も大きくなる可能性が高いため、主要チャネルから求められる分析パラメータをより優先（強調）する重み付けがなされる。
【００５６】
このように、適応ミキシング部１４１は、チャネル間相関（相互相関係数α）に応じて、主要チャネルと非主要チャネルとの間の重み付けを調整して、分析パラメータをミキシングする。
【００５７】
一例として、相互相関係数α＝０．７の場合について説明する。この場合、重み係数W1及び重み係数W2は、次式（５）のように求められる。
【数５】
2018221138.tif 000007
【００５８】
また、適応ミキシング部１４１は、分析パラメータがｎ次元である場合、次式（６）に示すようにミキシング後の分析パラメータMpを求めてもよい。
【数６】
2018221138.tif 000008
【００５９】
式（６）において、ParaDTCX-HQは主要チャネルの分析パラメータを示し、ParaNDTCX-HQは非主要チャネルの分析パラメータを示す。
【００６０】
最後に、符号化モード選択部１４２は、ＳＴ１０５で求められた分析パラメータMpを用いて、Ｌ信号及びＲ信号の双方に共通の符号化モードを選択する（ＳＴ１０６）。符号化モード選択部１４２における符号化モードの選択方法は、図２で説明したEVSコーデック（モノラル符号化）における選択方法と同じ方法でもよい。
【００６１】
このように、本実施の形態では、符号化装置１００は、ステレオ信号のチャネル間相関がある場合、各チャネル信号の符号化に用いる符号化モードを共通化する。こうすることで、ステレオ信号の両方のチャネルで異なる符号化モードが選択された場合に復号信号の主観品質が劣化してしまうような状況でも、符号化装置１００は、ステレオ信号の両方のチャネルに対して共通の符号化モードを用いて符号化することで、復号信号の主観品質が劣化することを防止することができる。よって、本実施の形態によれば、複数の符号化モードを切り替えて符号化処理を行うマルチモードモノラルコーデックを用いてステレオ信号を符号化する場合でも、ステレオ再生時の音声品質の劣化を抑えることができる。
【００６２】
また、符号化装置１００は、共通の符号化モードを選択する際に、主要チャネル及び非主要チャネルを特定し、相互相関係数αに応じて、主要チャネルの分析パラメータを強調して、双方の分析パラメータをミキシングする。すなわち、本実施の形態によれば、符号化装置１００は、双方のチャネルのチャネル間相関に応じて分析パラメータの強調度合いを調整することにより、共通の符号化モードを適切に選択することができる。
【００６３】
一方、符号化装置１００は、ステレオ信号のチャネル間相関が無い場合、各チャネル信号の符号化に用いる符号化モードを個別に選択する。これにより、ステレオ信号の各チャネルで最適な符号化モードがそれぞれ選択される。
【００６４】
以上より、本実施の形態によれば、符号化装置１００は、ステレオ信号の両方のチャネルのチャネル間相関に応じて、各チャネルに対して適切な符号化モードを選択することができるので、音声品質を改善することができる。
【００６５】
［実施の形態１の変形例１］
実施の形態１では、符号化装置１００が相互相関係数αに基づいて各チャネルの分析パラメータに対する重み係数を決定する場合について説明したが、重み係数の決定方法はこれに限定されるものではない。変形例１では、一例として、相互相関係数αの代わりに、チャネル間エネルギ差に基づいて重み係数を決定する方法について説明する。
【００６６】
図８は、本実施の形態に係るＤＭＡステレオ符号化部１０４の主な処理の流れを示すフロー図である。なお、図８において、図７と同様の処理については同一符号を付し、その説明を省略する。
【００６７】
具体的には、図８に示すＳＴ１０４ａにおいて、適応ミキシング部１４１（図６を参照）は、ＳＴ１０２で算出したチャネル間エネルギ差Δに基づいて、ＳＴ１０３で特定した主要チャネルの分析パラメータ及び非主要チャネルの分析パラメータに対する重み係数（ウェイト）を決定する。
【００６８】
具体的には、適応ミキシング部１４１は、チャネル間エネルギ差Δが大きいほど、主要チャネルの分析パラメータに対する重み係数W1を大きくし、非主要チャネルの分析パラメータに対する重み係数W2を小さくする。つまり、適応ミキシング部１４１は、チャネル間エネルギ差Δが大きいほど、主要チャネルを優先（強調）するような重み付けを行う。
【００６９】
図９は、適応ミキシング部１４１における重み係数を決定する処理（図８のＳＴ１０４ａ）の一例を示すフロー図である。また、図１０は、チャネル間エネルギ差Δと重み係数（W1、W2）との対応関係の一例を示す図である。
【００７０】
適応ミキシング部１４１は、チャネル間エネルギ差Δが小さいか否か（例えば、Δ≦閾値thrLであるか否か）を判断する（ＳＴ１４１）。チャネル間エネルギ差Δが小さい場合（ＳＴ１４１：Ｙｅｓ）、適応ミキシング部１４１は、チャネル間エネルギ差Δが小さい場合（Δ：Low level）に対応する重み係数（図１０では、（W1=0.6、W2=0.4）を選択する（ＳＴ１４２）。
【００７１】
また、適応ミキシング部１４１は、チャネル間エネルギ差Δが中間レベルであるか否か（例えば、閾値thrL＜Δ≦thrMであるか否か）を判断する（ＳＴ１４３）。チャネル間エネルギ差Δが中間レベルである場合（ＳＴ１４３：Ｙｅｓ）、適応ミキシング部１４１は、チャネル間エネルギ差Δが中間レベルである場合（Δ：Moderate level）に対応する重み係数（図１０では、（W1=0.7、W2=0.3）を選択する（ＳＴ１４４）。
【００７２】
また、適応ミキシング部１４１は、チャネル間エネルギ差Δが大きいか否か（例えば、Δ＞thrMであるか否か）を判断する（ＳＴ１４５）。チャネル間エネルギ差Δが大きい場合（ＳＴ１４５：Ｙｅｓ）、適応ミキシング部１４１は、チャネル間エネルギ差Δが大きい場合（Δ：High level）に対応する重み係数（図１０では、（W1=0.8、W2=0.2）を選択する（ＳＴ１４６）。
【００７３】
チャネル間エネルギ差Δが大きいほど、ステレオ信号における主要チャネルの影響は、非主要チャネルに対して大きくなる可能性が高い。このため、図１０に示す例では、式（４）と同様、常に主要チャネル側に大きな重み付けがなされることを保証しつつ、チャネル間エネルギ差Δが大きいほど、主要チャネルから求められる分析パラメータをより優先（強調）する重み付けがなされる。
【００７４】
このように、変形例１では、適応ミキシング部１４１は、チャネル間エネルギ差Δに応じて、主要チャネルと非主要チャネルとの間の分析パラメータに対する重み付けを調整して、分析パラメータをミキシングする。
【００７５】
このように、符号化装置１００は、ステレオ信号における主要チャネルと非主要チャネルとのエネルギ差に応じて、分析パラメータのミキシングにおける、主要チャネルの分析パラメータの強調度合いを変更する。これにより、符号化装置１００は、チャネル間エネルギ差が大きい場合には、主要チャネルをより強調した分析パラメータを用いて共通の符号化モードを選択することができる。また、符号化装置１００は、チャネル間エネルギ差が小さい場合には、非主要チャネルがより反映された分析パラメータを用いて共通の符号化モードを選択することができる。通常、信号分析は、エネルギで正規化してから行われることが多い。そのような場合には分析パラメータがエネルギの大小を反映しなくなる。このため、エネルギ差に応じて主要チャネルのパラメータを強調することは、分析パラメータの領域でミキシングする場合に意味のあることである。
【００７６】
［実施の形態１の変形例２］
実施の形態１の説明で用いた値（例えば、式（４）に示すW1の最小値：０．６、図１０に示す重み係数等）は、一例であり、他の数値でもよい。
【００７７】
また、式（４）では、相互相関係数αに基づいて重み係数を求める一例を示しているが、これに限定されず、例えば、チャネル間相関（相互相関係数α）及びチャネル間エネルギ差Δの双方に基づいて重み係数を決定してもよい。
【００７８】
具体的には、適応ミキシング部１４１は、次式（７）に従って重み係数を算出してもよい。
【数７】
2018221138.tif 000009
【００７９】
ここで、βは、チャネル間エネルギ差Δに基づいて設定される値である。例えば、図１０におけるチャネル間エネルギ差Δと重み係数W1との対応関係と同様にして、チャネル間エネルギ差Δが大きいほど、βの値が大きくなってもよい。これにより、チャネル間エネルギ差Δが大きいほど、主要チャネルの分析パラメータに対する重み係数W1（最小値β）が大きくなる。
【００８０】
よって、適応ミキシング部１４１は、チャネル間相関によるチャネル間の信号類似度、及び、チャネル間エネルギ差の双方に応じて、主要チャネル及び非主要チャネルの強調度合い（優先度）を調整して、分析パラメータをミキシングすることができる。
【００８１】
（実施の形態２）
符号化モードの判定結果（選択結果）がフレーム間で頻繁に切り替わると、復号信号の主観品質の劣化につながることがある。そこで、本実施の形態では、フレーム間での符号化モードの判定結果が頻繁に切り替わることを抑える方法について説明する。
【００８２】
［符号化装置の構成］
本実施の形態に係る符号化装置は、実施の形態１に係る符号化装置１００と基本構成が共通するので、図５を援用して説明する。ただし、本実施の形態では、符号化装置１００は、図５に示すＤＭＡステレオ符号化部１０４の代わりに、図１１に示すＤＭＡステレオ符号化部１５０を備える。
【００８３】
図１１は、本実施の形態に係るＤＭＡステレオ符号化部１５０の構成例を示すブロック図である。
【００８４】
なお、図１１において、実施の形態１（図６）と同様の構成には同様の符号を付し、その説明を省略する。具体的には、図１１に示すＤＭＡステレオ符号化部１５０は、実施の形態１の構成（図６）と比較して、判定訂正部１５１を新たに備える。
【００８５】
また、本実施の形態では、信号分析部１０１（Lch信号分析部）は、実施の形態１の動作に加え、Lch分析パラメータに基づいて判定される符号化モード（例えば、図２を参照）を示すLch符号化モード判定結果（Left channel coding mode decision）を判定訂正部１５１に出力する。同様に、信号分析部１０１（Rch信号分析部）は、実施の形態１の動作に加え、Rch分析パラメータに基づいて判定される符号化モード（例えば、図２を参照）を示すRch符号化モード判定結果（Right channel coding mode decision）を判定訂正部１５１に出力する。
【００８６】
ＤＭＡステレオ符号化部１５０において、判定訂正部１５１は、過去のフレームにおいて適用された符号化モード、及び、信号分析部１０１から入力されるLch符号化モード判定結果、Rch符号化モード判定結果に基づいて、符号化モード選択部１４２から入力される符号化モード判定結果を訂正するか否かを判断する。
【００８７】
なお、ここでは、判定訂正部１５１に入力される符号化モードを「decision 1」と呼び、判定訂正部１５１から出力される符号化モードを「decision 2」と呼ぶ。
【００８８】
判定訂正部１５１は、符号化モード判定結果の訂正が不要と判断した場合、符号化モード判定結果を訂正せずにLch符号化部１４３及びRch符号化部１４４にそれぞれ出力する。一方、符号化モード判定結果の訂正が必要と判断した場合、符号化モード判定結果を訂正し、訂正後の符号化モード判定結果をLch符号化部１４３及びRch符号化部１４４にそれぞれ出力する。
【００８９】
図１２は、判定訂正部１５１における符号化モードの判定訂正処理の流れの一例を示すフロー図である。
【００９０】
図１２において、判定訂正部１５１は、符号化モード選択部１４２における現フレームの符号化モード判定結果（decision 1）が過去フレーム（例えば、１つ前のフレーム）において適用された符号化モードと同一であるか否かを判断する（ＳＴ１５１）。
【００９１】
符号化モード判定結果（decision 1）が過去フレームの符号化モードと同一である場合（ＳＴ１５１：Ｙｅｓ）、判定訂正部１５１は、符号化モード判定結果（decision 1）に対する訂正処理を行わずに処理を終了する（ＳＴ１５２）。
【００９２】
一方、符号化モード判定結果（decision 1）が過去フレームの符号化モードと同一ではない場合（ＳＴ１５１：Ｎｏ）、判定訂正部１５１は、過去フレーム（例えば、１つ前のフレーム）で用いられた符号化モードが、現フレームのLch符号化モード判定結果又は現フレームのRch符号化モード判定結果と同一であるか否かを判断する（ＳＴ１５３）。
【００９３】
ＳＴ１５３において，過去フレームで用いられた符号化モードが、現フレームのLch符号化モード判定結果又は現フレームのRch符号化モード判定結果と同一でない場合（ＳＴ１５３：Ｎｏ）、判定訂正部１５１は、符号化モード判定結果（decision 1）に対する訂正処理を行わずに処理を終了する（ＳＴ１５２）。
【００９４】
一方、過去フレームの符号化モードが、現フレームのLch符号化モード判定結果又は現フレームのRch符号化モード判定結果と同一である場合（ＳＴ１５３：Ｙｅｓ）、判定訂正部１５１は、現フレームの符号化モード判定結果及び過去フレームの符号化モードを用いて符号化モード判定結果（decision 1）の訂正処理（スムージング処理）を行う（ＳＴ１５４）。
【００９５】
すなわち、判定訂正部１５１は、現フレームで選択された共通の符号化モード（decision１）が、過去のフレームで選択された共通の符号化モードと異なり、かつ、過去のフレームで選択された共通の符号化モードが、現フレームのLch符号化モード判定結果か現フレームのRch符号化モード判定結果のいずれかと同じ場合に、現フレームの共通の符号化モードを再選択（訂正）する。
【００９６】
例えば、判定訂正部１５１は、次式（８）に従って、decision 1の判定処理において用いた分析パラメータMpを修正する。
【数８】
2018221138.tif 000010
【００９７】
式（８）において、Mp[-1]は１つ前のフレーム（過去フレーム）における分析パラメータMpを示し、Wは平滑化係数を示し、例えば、W=0.8としてもよい。なお、平滑化係数Wの値は０．８に限定されるものではない。また、スムージング処理において対象とする過去フレームは、式（８）に示すように１つ前のフレームに限らず、過去の複数フレームを対象としてもよい。
【００９８】
スムージング処理後に、判定訂正部１５１は、修正後の分析パラメータMpを用いて、符号化モードの再選択（再判定）を行う（ＳＴ１５５）。なお、符号化モードの再選択時における符号化モードの選択方法は、符号化モード選択部１４２における選択方法と同様でもよい。
【００９９】
このように、分析パラメータMpは、１つ前のフレーム及び現フレームに渡って平滑化される。また、式（８）に示すように、平滑化係数Wが大きいほど、修正後の分析パラメータMpは、過去フレームの分析パラメータMp[-1]により影響を受ける。すなわち、平滑化係数Wが大きいほど、修正後の分析パラメータMpに基づく符号化モードの再選択において、過去フレームで用いられた符号化モードが選択されやすくなる。
【０１００】
これにより、本実施の形態では、符号化モードの判定結果（選択結果）がフレーム間で頻繁に切り替わることを防止し、復号信号の主観品質の劣化を抑えることができる。
【０１０１】
（実施の形態３）
［符号化装置の構成］
図１３は、本実施の形態に係る符号化装置２００の構成を示すブロック図である。
【０１０２】
なお、図１３において、実施の形態１（図５）と同様の構成には同様の符号を付し、その説明を省略する。具体的には、図１３に示す符号化装置２００は、実施の形態１の構成（図５）に対して、ＤＭ−Ｍ／Ｓ（Mid/Side）変換部２０２、及び、Ｍ／Ｓステレオ符号化部２０４を新たに備える。
【０１０３】
符号化装置２００において、チャネル間相関算出部２０１は、算出したチャネル間相関（相互相関係数α）に基づいて、ＤＭステレオ符号化及びＤＭＡステレオ符号化に加え、Ｍ／Ｓステレオ符号化の中から、１つのステレオ符号化モードを選択する。チャネル相関算出部２０１は、選択した結果を示すステレオモード判定フラグを、ＤＭ−Ｍ／Ｓ変換部２０２、切替スイッチ２０３及び多重化部１０６に出力する。
【０１０４】
例えば、図１４に示すように、チャネル間相関算出部２０１は、相互相関係数αが０の場合にＤＭステレオ符号化モードと判定し、相互相関係数αが０より大きく、０．６以下の場合にＤＭＡステレオ符号化モードと判定し、相互相関係数αが０．６より大きい場合にＭ／Ｓステレオ符号化モードと判定してもよい。
【０１０５】
すなわち、チャネル間相関が高い場合（α：High。ここでは、0.6＜αの範囲）にはＭ／Ｓステレオ符号化が選択され、チャネル間相関が低い場合（α＝0）にはＤＭステレオ符号化が選択され、チャネル間相関が上記範囲の何れにも該当しない場合（α：Weak。ここでは、0＜α≦0.6）にはＤＭＡステレオ符号化が選択される。
【０１０６】
なお、図１４に示す相互相関係数αの範囲は一例であり、これに限定されるものではない。
【０１０７】
ＤＭ−Ｍ／Ｓ変換部２０２は、チャネル間相関算出部２０１から入力されるステレオモード判定フラグがＭ／Ｓステレオ符号化である場合には、Ｌ／Ｒ信号を後述するようにＭ／Ｓ信号に変換し、信号分析部１０１及び切替スイッチ２０３に出力する。ＤＭ−Ｍ／Ｓ変換部２０２は、ステレオモード判定フラグがＤＭステレオ符号化モード又はＤＭＡステレオ符号化モードの場合には、Ｌ／Ｒ信号をそのまま信号分析部１０１及び切替スイッチ２０３に出力する。
【０１０８】
切替スイッチ２０３は、実施の形態１（切替スイッチ１０３）の動作に加え、チャネル間相関算出部２０１から入力されるステレオモード判定フラグがＭ／Ｓステレオ符号化モードである場合、入力されるＬ信号、Ｒ信号、及び分析パラメータをＭ／Ｓステレオ符号化部２０４に出力する。
【０１０９】
Ｍ／Ｓステレオ符号化部２０４は、切替スイッチ２０３から入力されるＬ／Ｒの和信号、Ｌ／Ｒの差信号、及びそれぞれに対する分析パラメータを用いて、Ｍ／Ｓステレオ符号化を行う。Ｍ／Ｓステレオ符号化を行う場合には、ＤＭ−Ｍ／Ｓ変換部２０２において、ステレオ信号のＬ信号及びＲ信号が、双方のチャネルの和（sum）であるMidチャネルと、双方のチャネルの差（difference）であるSideチャネルとに変換されている。なお、Ｍ／Ｓステレオ符号化の詳細については、例えば、非特許文献２に記載された方法を用いてもよい。
【０１１０】
チャネル間相関が高い場合には、Ｍ／Ｓステレオ符号化は、ステレオ符号化と比較して、より効率的な符号化である。具体的には、チャネル間相関が高い場合には、双方のチャネルの差であるSideチャネルがゼロに近い値となるので、符号化情報の情報量を削減することができる。一方、チャネル間相関が低い場合には、Ｍ／Ｓステレオ符号化と比較して、デュアルモノ符号化によって符号化情報の情報量を削減することができる。また、チャネル間相関が高い場合には、音源が一つの点音源（例：一人の人が話しているようなケース）である可能性が高い。このような場合は、モノラル化した信号（Midチャネル信号）及びSideチャネル信号を用いてＬ／Ｒに振り分けるようにしたほうが安定したステレオ定位感が得られる。
【０１１１】
また、Ｍ／Ｓステレオ符号化では、上述したように、双方のチャネルの和及び差を符号化情報として生成するため、復号側（図示せず）では、フレーム毎の符号化情報（和及び差）に基づいて復号信号を復号する。つまり、和信号であるMidチャネル信号と差信号であるSideチャネル信号との和がＲチャネル信号となり、和信号（Midチャネル信号）と差信号（Sideチャネル信号）との差がＬチャネル信号となる。つまり、Midチャネル信号とSideチャネル信号の符号化モードが異なっていても、双方の信号がＬチャネルとＲチャネルの双方に反映されるため、符号化モードを必ずしも統一する必要がない。すなわち、Ｍ／Ｓステレオ符号化を用いれば、チャネル間で符号化モードが異なることによる、復号信号の主観品質の劣化を抑えることができる。
【０１１２】
このように、符号化装置２００は、チャネル間相関（相互相関係数α）に応じて、デュアルモノ符号化（ＤＭＡステレオ符号化又はＤＭステレオ符号化）及びＭ／Ｓステレオ符号化を切り替える。こうすることで、符号化装置２００は、チャネル間相関に応じて、適切な符号化モードを選択して、ステレオ信号を符号化することができるので、復号信号の主観品質を改善することができ、さらに、符号化情報を削減することができる。
【０１１３】
（実施の形態４）
本実施の形態では、チャネル間相関（相互相関係数α）を効率的に求める方法について説明する。
【０１１４】
本実施の形態に係る符号化装置は、実施の形態１に係る符号化装置１００と基本構成が共通するので、図５を援用して説明する。ただし、本実施の形態では、符号化装置１００は、図５に示すチャネル間相関算出部１０２の代わりに、図１５に示すチャネル間相関算出部３０１を備える。
【０１１５】
実施の形態１で説明した式（１）に示す相互相関係数αは、次式（９）で表される。
【数９】
2018221138.tif 000011
【０１１６】
すなわち、式（９）に示すように、相互相関係数αは、クロススペクトル成分（分子項の「Cross-Spectrum」）と、左チャネル及び右チャネルのエネルギ成分（分母項の「Left Channel Energy」及び「Right Channel Energy」）とに分けることができる。
【０１１７】
本実施の形態では、相互相関係数αの演算の際に、左チャネル及び右チャネルの全ての周波数スペクトルパラメータ（スペクトル係数）を用いるのではなく、一部の帯域の周波数スペクトルパラメータを用いることにより、相互相関係数αの演算量を削減する。
【０１１８】
図１５は、本実施の形態に係る信号分析部１０１及びチャネル間相関算出部３０１の構成例を示すブロック図である。
【０１１９】
信号分析部１０１は、Lch周波数領域変換部１１１と、Lchスペクトルバンドエネルギ算出部１１２と、Rch周波数領域変換部１１３と、Rchスペクトルバンドエネルギ算出部１１４と、を含む構成を採る。
【０１２０】
また、チャネル間相関算出部３０１は、エネルギ閾値算出部３１１と、主要帯域特定部３１２と、Lch主要帯域エネルギ算出部３１３と、Lch主要帯域スペクトル取得部３１４と、Rch主要帯域エネルギ算出部３１５と、Rch主要帯域スペクトル取得部３１６と、クロススペクトル算出部３１７と、相関演算部３１８と、を含む構成を採る。
【０１２１】
信号分析部１０１において、Lch周波数領域変換部１１１は、入力されるＬ信号を周波数領域変換し、Lch周波数スペクトルパラメータをLchスペクトルバンドエネルギ算出部１１２及びLch主要帯域スペクトル取得部３１４に出力する。
【０１２２】
Lchスペクトルバンドエネルギ算出部１１２は、Lch周波数領域変換部１１１から入力されるLch周波数スペクトルパラメータを複数のスペクトルバンドにグループ化し、各スペクトルバンドのエネルギを算出する。Lchスペクトルバンドエネルギ算出部１１２は、算出したLchバンドエネルギをエネルギ閾値算出部３１１、主要帯域特定部３１２及びLch主要帯域エネルギ算出部３１３に出力する。
【０１２３】
Rch周波数領域変換部１１３は、入力されるＲ信号を周波数領域変換し、Rch周波数スペクトルパラメータをRchスペクトルバンドエネルギ算出部１１４及びRch主要帯域スペクトル取得部３１６に出力する。
【０１２４】
Rchスペクトルバンドエネルギ算出部１１４は、Rch周波数領域変換部１１３から入力されるRch周波数スペクトルパラメータを複数のスペクトルバンドにグループ化し、各スペクトルバンドのエネルギを算出する。Rchスペクトルバンドエネルギ算出部１１４は、算出したRchバンドエネルギをエネルギ閾値算出部３１１、主要帯域特定部３１２及びRch主要帯域エネルギ算出部３１５に出力する。
【０１２５】
なお、図１５に示す信号分析部１０１における周波数領域変換及びスペクトルバンドエネルギ算出は、本チャネル間相関算出部の適用先であるコーデックにおいて行われる処理であるものとする。この場合、図１５に示す信号分析部１０１の各構成部は、本実施の形態に係るチャネル間相関算出のために新たに備えられる構成ではない。つまり、信号分析部１０１の処理量は増加しない。
【０１２６】
次に、チャネル間相関算出部３０１において、エネルギ閾値算出部３１１は、Lchスペクトルバンドエネルギ算出部１１２から入力されるLchバンドエネルギ、及び、Rchスペクトルバンドエネルギ算出部１１４から入力されるRchバンドエネルギを用いて、Lchエネルギ閾値、及び、Rchエネルギ閾値をそれぞれ算出する。エネルギ閾値算出部３１１は、算出したLch/Rchエネルギ閾値を主要帯域特定部３１２に出力する。
【０１２７】
主要帯域特定部３１２は、Lchスペクトルバンドエネルギ算出部１１２から入力されるLchバンドエネルギのうち、エネルギ閾値算出部３１１から入力されるLchエネルギ閾値より大きいエネルギを有するスペクトルバンドを、Lch主要帯域として特定する。同様に、主要帯域特定部３１２は、Rchスペクトルバンドエネルギ算出部１１４から入力されるRchバンドエネルギのうち、エネルギ閾値算出部３１１から入力されるRchエネルギ閾値より大きいエネルギを有するスペクトルバンドを、Rch主要帯域として特定する。主要帯域特定部３１２は、特定したLch主要帯域とRch主要帯域の総和、すなわちLch主要帯域またはRch主要帯域のいずれかに該当する帯域を「主要帯域」として、Lch主要帯域エネルギ算出部３１３及びLch主要帯域スペクトル取得部３１４及びRch主要帯域エネルギ算出部３１５及びRch主要帯域スペクトル取得部３１６に出力する。
【０１２８】
Lch主要帯域エネルギ算出部３１３は、Lchスペクトルバンドエネルギ算出部１１２から入力されるLchバンドエネルギのうち、主要帯域特定部３１２から入力される主要帯域に対応するバンドエネルギの総和を算出し、Lch主要帯域エネルギとして相関演算部３１８に出力する。
【０１２９】
Lch主要帯域スペクトル取得部３１４は、Lch周波数領域変換部１１１から入力されるLch周波数スペクトルパラメータのうち、主要帯域特定部３１２から入力される主要帯域に対応するLch周波数スペクトルパラメータを取り出し、Lch主要帯域スペクトルとしてクロススペクトル算出部３１７に出力する。
【０１３０】
Rch主要帯域エネルギ算出部３１５は、Rchスペクトルバンドエネルギ算出部１１４から入力されるRchバンドエネルギのうち、主要帯域特定部３１２から入力される主要帯域に対応するバンドエネルギの総和を算出し、Rch主要帯域エネルギとして相関演算部３１８に出力する。
【０１３１】
Rch主要帯域スペクトル取得部３１６は、Rch周波数領域変換部１１３から入力されるRch周波数スペクトルパラメータのうち、主要帯域特定部３１２から入力される主要帯域に対応するRch周波数スペクトルパラメータを取り出し、Rch主要帯域スペクトルとしてクロススペクトル算出部３１７に出力する。
【０１３２】
クロススペクトル算出部３１７は、Lch主要帯域スペクトル取得部３１４から入力されるLch主要帯域スペクトル、及び、Rch主要帯域スペクトル取得部３１６から入力されるRch主要帯域スペクトルを用いて、クロススペクトル（式（９）の分子項）を算出する。クロススペクトル算出部３１７は、算出したクロススペクトルを相関演算部３１８に出力する。
【０１３３】
相関演算部３１８は、Lch主要帯域エネルギ算出部３１３から入力されるLch主要帯域エネルギ、及び、Rch主要帯域エネルギ算出部３１５から入力されるRch主要帯域エネルギを用いて、左チャネル及び右チャネルのエネルギ（式（９）の分母項）を算出する。そして、相関演算部３１８は、算出したエネルギ（式（９）の分母項）と、クロススペクトル算出部３１７から入力されるクロススペクトル（式（９）の分子項）とを用いて、チャネル間相関（式（９）の相互相関係数α）を算出する。
【０１３４】
図１６は、チャネル間相関の算出処理に関する、信号分析部１０１及びチャネル間相関算出部３０１におけるＬ信号に対する処理の一例を示す。
【０１３５】
図１６に示すように、Lchスペクトルバンドエネルギ算出部１１２は、Lch周波数スペクトルパラメータlを、Nbands個のバンドにグループ化し、バンドkb（kb＝0〜（Nbands-1））のLchバンドエネルギLbandend(kb)を算出する。
【０１３６】
エネルギ閾値算出部３１１は、LchバンドエネルギLbandend(kb)を用いてLchエネルギ閾値l-を算出する。例えば、エネルギ閾値算出部３１１は、LchバンドエネルギLbandend(kb)の平均値、又は、非特許文献１に記載されたように、LchバンドエネルギLbandend(kb)の平均値及び標準偏差を用いて定義してもよい。
【０１３７】
例えば、バンドエネルギの平均Avgeneと標準偏差σbandeneとを用いる場合、エネルギ閾値thrは次式（１０）で表される。
【数１０】
2018221138.tif 000012
【０１３８】
また、バンドエネルギの平均Avgeneは次式（１１）で表される。
【数１１】
2018221138.tif 000013
【０１３９】
次に、主要帯域特定部３１２は、バンドkb（kb＝0〜（Nbands-1））のうち、LchバンドエネルギLbandend(kb)がLchエネルギ閾値l-より大きいバンドを主要帯域として特定する。図１６では、一例として、バンドkb（kb＝0〜（Nbands-1））のうち、kb＝0,1,2,5,6,7が主要帯域lidxとして特定されている。
【０１４０】
次に、Lch主要帯域エネルギ算出部３１３は、主要帯域lidxのバンドエネルギの総和をLchエネルギ（Left channel energy）として算出する。なお、LchバンドエネルギLbandend(kb)は信号分析部１０１で既に算出されているので、主要帯域エネルギ算出部３１３は、図１６に示すように、全バンドkbのエネルギの総和をLchエネルギとして算出してもよい。
【０１４１】
Lch主要帯域スペクトル取得部３１４は、Lch周波数スペクトルパラメータlのうち、Lch主要帯域lidxに含まれるLch周波数スペクトルパラメータL(lidx)を取得する。
【０１４２】
以上、Lchに対する処理について説明したが、信号分析部１０１及びチャネル間相関算出部３０１におけるＲ信号に対する処理についても図１６と同様に行えばよい（図示せず）。これにより、Ｒ信号に対して、Rchエネルギ（Right channel energy）、及び、Rch主要帯域ridxに含まれるRch周波数スペクトルパラメータR(ridx)が得られる。
【０１４３】
そして、クロススペクトル算出部３１７は、図１６に示すように、Lch主要帯域のLch周波数スペクトルパラメータL(lidx)、及び、Rch主要帯域のRch周波数スペクトルパラメータR(ridx)を用いてクロススペクトル（Cross-Spectrum）を算出する。
【０１４４】
ここで、idxlenは、主要帯域のバンド数（例えば、図１６の例ではidxlen=6）を示し、kは主要帯域内のスペクトルバンドのインデックス（例えば、図１６の例では、kb＝0,1,2,5,6,7に対してk=1〜6）を示す。
【０１４５】
最後に、相関演算部３１８は、Lchエネルギ（Left channel energy）、Rchエネルギ（Right channel energy）及びクロススペクトル（Cross-Spectrum）を用いて、式（９）に従ってチャネル間相関（α）を算出する。
【０１４６】
このように、本実施の形態によれば、チャネル間相関算出部３０１は、チャネル間相関を算出する際に、一部のスペクトルバンドを用いてチャネル間相関を算出する。また、チャネル間相関算出部３０１は、一部のスペクトルバンドとして、バンドエネルギがエネルギ閾値より大きい主要帯域を用いる。これにより、例えば、式（１２）に示すように、クロススペクトルの演算の対象を主要帯域の周波数スペクトルパラメータに限定することができる。よって、本実施の形態によれば、チャネル間相関の精度を維持しつつ、演算量を削減することができる。
【０１４７】
［実施の形態４の変形例１］
本実施の形態では、主要帯域特定部３１２においてLch及びRchの双方のバンドエネルギを用いて主要帯域を特定する場合について説明したが、主要帯域の特定方法はこれに限定されない。例えば、主要帯域特定部３１２は、Lch及びRchの中から主要チャネルを選択し、選択された主要チャネルのバンドエネルギを用いて、Lch及びRchの双方の主要帯域を特定してもよい。
【０１４８】
［実施の形態４の変形例２］
実施の形態４では、チャネル間相関算出部３０１において、主要帯域特定部３１２で選択されるスペクトルバンド（主要帯域）に含まれる周波数スペクトルパラメータを用いてチャネル間相関を求める場合について説明した。これに対して、変形例では、主要帯域の中から、主要なスペクトル成分をさらに選択して、チャネル間相関を求める場合について説明する。
【０１４９】
図１７は、変形例２に係るチャネル間相関算出部４０１の構成例を示すブロック図である。なお、図１７において、図１５と同様の構成には同一の符号を付し、その説明を省略する。図１７では、エネルギ閾値算出部３１１及び主要帯域特定部３１２は、Lch及びRchに対してそれぞれ備えられる。
【０１５０】
図１７において、Lch主要帯域分析部４１１は、Lch周波数領域変換部１１１から入力されるLch周波数スペクトルパラメータのうち、主要帯域特定部３１２−１から入力されるLch主要帯域内の周波数スペクトルパラメータの振幅（エネルギ）を算出し、Lch振幅閾値算出部４１２に出力する。
【０１５１】
Lch振幅閾値算出部４１２は、Lch主要帯域分析部４１１から入力される、主要帯域として特定されたスペクトルバンド内のLch周波数スペクトルパラメータの振幅値を用いて、平均振幅を算出する。Lch振幅閾値算出部４１２は、算出した平均振幅値をLch振幅閾値としてLch/Rch主要帯域スペクトル取得部４１５に出力する。
【０１５２】
また、Rch主要帯域分析部４１３及びRch振幅閾値算出部４１４は、Rchに対して、Lch主要帯域分析部４１１及びLch振幅閾値算出部４１２と同様の処理を行う。
【０１５３】
Lch/Rch主要帯域スペクトル取得部４１５は、Lch周波数領域変換部１１１から入力されるLch周波数スペクトルパラメータのうち、主要帯域に含まれ、かつ、Lch振幅閾値算出部４１２から入力されるLch振幅閾値より大きい振幅（エネルギ）を有するLch周波数スペクトルパラメータを選択し、Rch周波数領域変換部１１３から入力されるRch周波数スペクトルパラメータのうち、主要帯域に含まれ、かつ、Rch振幅閾値算出部４１４から入力されるRch振幅閾値より大きい振幅（エネルギ）を有するRch周波数スペクトルパラメータを選択する。そして、Lch/Rch主要帯域スペクトル取得部４１５は、LchとRchの少なくとも一方の周波数スペクトルパラメータが選ばれている周波数成分を相関演算に用いる、LchとRchに共通する周波数成分として選択する。Lch/Rch主要帯域スペクトル取得部４１５は、選択した周波数成分のLch周波数スペクトルパラメータ及びRch周波数スペクトルパラメータを相関演算部４１７に出力する。
【０１５４】
相関演算部４１７は、Lch/Rch主要帯域スペクトル取得部４１５から入力されるLch周波数スペクトルパラメータ及びRch周波数スペクトルパラメータを用いて、クロススペクトル（式（９）の分子項）を算出する。ここで、クロススペクトルの演算に用いる周波数スペクトルパラメータがLch主要帯域及びRch主要帯域内の特にエネルギの大きい成分に制限されているため、Lch主要帯域及びRch主要帯域内の全ての周波数スペクトルパラメータを用いる場合と比較して、演算量が削減される。
【０１５５】
また、相関演算部４１７は、相関算出部３１８と同様、式（９）の分母項も算出し、式（９）に示す相互相関係数αを算出する。
【０１５６】
このように、主要帯域特定部３１２で特定された主張帯域に含まれるスペクトル成分の数を更に限定することで、クロススペクトルの演算量を更に削減することができる。
【０１５７】
以上、本実施の形態の変形例１、２について説明した。
【０１５８】
なお、本実施の形態で説明した主要帯域を特定する方法は、スペクトルパラメータを符号化する種々の符号化方式に適応することができる。例えば、非特許文献３に示すようなBCC（Binaural Cue Coding）の原理を利用したパラメトリックステレオ符号化に適応することで、低ビットレート化、低演算量化を図ることができる。パラメトリックステレオ符号化では、チャネル間レベル差（ICLD：Inter Channel Level Difference）、チャネル間時間差（ICTD：Inter Channel Time Difference）、チャネル間コヒーレンス（ICC：Inter Channel Coherence）等のパラメータをサイド情報としてスペクトルバンド毎に符号化する。このとき、本実施の形態で説明したようなスペクトルバンドの選択及びスペクトル成分の選択を用いて、選択されたスペクトルバンド又はスペクトル成分のみを用いてICLD、ICTD、ICC等を計算すれば、サイド情報の算出に必要な演算量を減らすことができる。
【０１５９】
以上、本開示の各実施の形態について説明した。
【０１６０】
なお、上記実施の形態において、チャネル間エネルギ差Δ（例えば、式（２））を算出する際、主要チャネルの判定結果が安定するように、チャネル間エネルギ差の算出に、チャネルエネルギの瞬時値（現在のフレームにおけるチャネルエネルギ）ではなく、チャネルエネルギの長期平均を用いてもよい。例えば、符号化装置は、次式（１２）に従って、チャネル間エネルギ差Δを求め、求めたチャネル間エネルギ差Δを用いて主要チャネルの判定又は重み係数の取得を行ってもよい。これにより、符号化装置は、主要チャネルの判定又は重み係数の取得を精度良く行うことができる。
【数１２】
2018221138.tif 000014
【０１６１】
式（１２）において、Nはチャネルエネルギの長期平均の対象となるフレーム数を示し、framenocurは現フレームインデックスを示す。すなわち、（framenocur-m）は現フレームからｍフレーム前のフレームを表す。
【０１６２】
また、上記各実施の形態を組み合わせて適用してもよい。例えば、実施の形態３の符号化装置２００（図１３）において、ＤＭＡステレオ符号化部１０４の代わりに、実施の形態２に係るＤＭＡステレオ符号化部１５０（図１１）を備えてもよい。また、実施の形態３の符号化装置２００（図１３）において、チャネル間相関算出部１０２の代わりに、実施の形態４に係るチャネル間相関算出部３０１（図１５）又は４０１（図１７）を備えてもよい。
【０１６３】
また、上記実施の形態では、符号化モードとして、ACELP、TCX、HQ MDCT、GSC等を一例として用いる場合について説明したが、これらに限定されるものではない。
【０１６４】
また、本開示はソフトウェア、ハードウェア、又は、ハードウェアと連携したソフトウェアで実現することが可能である。上記実施の形態の説明に用いた各機能ブロックは、部分的に又は全体的に、集積回路であるＬＳＩとして実現され、上記実施の形態で説明した各プロセスは、部分的に又は全体的に、一つのＬＳＩ又はＬＳＩの組み合わせによって制御されてもよい。ＬＳＩは個々のチップから構成されてもよいし、機能ブロックの一部または全てを含むように一つのチップから構成されてもよい。ＬＳＩはデータの入力と出力を備えてもよい。ＬＳＩは、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。集積回路化の手法はＬＳＩに限るものではなく、専用回路、汎用プロセッサ又は専用プロセッサで実現してもよい。また、ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。本開示は、デジタル処理又はアナログ処理として実現されてもよい。さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適用等が可能性としてありえる。
【０１６５】
本開示の符号化装置は、ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出する算出回路と、前記チャネル間相関が閾値より大きい場合に共通の符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化し、前記チャネル間相関が前記閾値以下の場合に前記左チャネル信号及び前記右チャネル信号に対して個別に判定された符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化する符号化回路と、を具備する。
【０１６６】
本開示の符号化装置において、前記符号化回路は、左チャネルと右チャネルについて主要チャネルと非主要チャネルとを特定し、前記主要チャネルの符号化モードを判定するための第１のパラメータと、前記非主要チャネルの符号化モードを判定するための第２のパラメータとに対して重み付け加算を行い、前記重み付け加算によって得られる重み付けパラメータに基づいて前記共通の符号化モードを選択する。
【０１６７】
本開示の符号化装置において、前記第１のパラメータに対する第１の重み係数は、前記第２のパラメータに対する第２の重み係数より大きく、前記チャネル間相関が小さいほど、前記第１の重み係数は大きい。
【０１６８】
本開示の符号化装置において、前記第１のパラメータに対する第１の重み係数は、前記第２のパラメータに対する第２の重み係数より大きく、前記左チャネル信号と前記右チャネル信号との間のエネルギ差が大きいほど、前記第１の重み係数は大きい。
【０１６９】
本開示の符号化装置において、前記符号化回路は、現フレームで選択された前記共通の符号化モードが、過去のフレームで選択された前記共通の符号化モード、現フレームの前記第１のパラメータに基づいて判定される符号化モードと異なり、かつ、現フレームの前記第２のパラメータに基づいて判定される符号化モードの何れかと同一である場合、現フレームの前記共通の符号化モードを再選択する。
【０１７０】
本開示の符号化装置において、前記符号化回路は、現フレームの前記重み付けパラメータと、過去フレームの前記重み付けパラメータとを用いてスムージング処理を行い、前記スムージング処理後の重み付けパラメータに基づいて前記共通の符号化モードを再選択する。
【０１７１】
本開示の符号化装置において、前記符号化回路は、更に、前記チャネル間相関が、前記閾値よりも大きい第２の閾値よりも大きい場合、前記左チャネル信号及び前記右チャネル信号に対して、Mid/Sideステレオ符号化を行う。
【０１７２】
本開示の符号化装置において、前記算出回路は、前記左チャネル信号及び前記右チャネル信号のうち、一部の帯域の周波数スペクトルパラメータを用いて、前記チャネル間相関を算出する。
【０１７３】
本開示の符号化方法は、ステレオ信号を構成する左チャネル信号と右チャネル信号を用いて、左チャネルと右チャネルとの間のチャネル間相関を算出し、前記チャネル間相関が閾値より大きい場合に共通の符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化し、前記チャネル間相関が前記閾値以下の場合に前記左チャネル信号及び前記右チャネル信号に対して個別に判定された符号化モードを用いて、前記左チャネル信号及び前記右チャネル信号をそれぞれ符号化する。
【産業上の利用可能性】
【０１７４】
本開示の一態様は、マルチモード符号化技術を用いた音声通信システムに有用である。
【符号の説明】
【０１７５】
１００，２００  符号化装置
１０１  信号分析部
１０２，２０１，３０１，４０１  チャネル間相関算出部
１０３，２０３  切替スイッチ
１０４，１５０  ＤＭＡステレオ符号化部
１０５  ＤＭステレオ符号化部
１０６  多重化部
１４１  適応ミキシング部
１４２  符号化モード選択部
１４３  Lch符号化部
１４４  Rch符号化部
１４５  ビットストリーム生成部
１５１  判定訂正部
２０２  ＤＭ−Ｍ／Ｓ変換部
２０４  Ｍ／Ｓステレオ符号化部
３１１  エネルギ閾値算出部
３１２  主要帯域特定部
３１３  Lch主要帯域エネルギ算出部
３１４  Lch主要帯域スペクトル取得部
３１５  Rch主要帯域エネルギ算出部
３１６  Rch主要帯域スペクトル取得部
３１７  クロススペクトル算出部
３１８，４１７  相関演算部
４１１  Lch主要帯域分析部
４１２  Lch振幅閾値算出部
４１３  Rch主要帯域分析部
４１４  Rch振幅閾値算出部
４１５  Lch/Rch主要帯域スペクトル取得部
【図１】
2018221138.tif 000015
【図２】
2018221138.tif 000016
【図３】
2018221138.tif 000017
【図４】
2018221138.tif 000018
【図５】
2018221138.tif 000019
【図６】
2018221138.tif 000020
【図７】
2018221138.tif 000021
【図８】
2018221138.tif 000022
【図９】
2018221138.tif 000023
【図１０】
2018221138.tif 000024
【図１１】
2018221138.tif 000025
【図１２】
2018221138.tif 000026
【図１３】
2018221138.tif 000027
【図１４】
2018221138.tif 000028
【図１５】
2018221138.tif 000029
【図１６】
2018221138.tif 000030
【図１７】
2018221138.tif 000031
【国際調査報告】
2018221138.tif 000032
2018221138.tif 000033
2018221138.tif 000034