(57)【要約】高調波ノイズ重み付けの大きさを最適化することができるように高調波ノイズ重み付け（ＨＮＷ）係数（εｐ）の値を選択する必要性に対処するために、高調波ノイズ重み付けをデジタル音声符号器において行なう方法及び装置を本明細書に開示する。動作状態の間、受信音声を解析して（５０３）ピッチ周期を求める。次に、ＨＮＷ係数をピッチ周期に基づいて選択し（５０５）、知覚ノイズ重み付けフィルタ（Ｃ（ｚ））を高調波ノイズ重み付け（ＨＮＷ）係数（εｐ）に基づいて決定する。
2007513364.tif 000002
【特許請求の範囲】
【請求項１】
高調波ノイズ重み付けをデジタル音声符号器において行なう方法であって、
音声入力ｓ（ｎ）を受信するステップと、
ピッチ周期（Ｄ）を前記音声入力から求めるステップと、
高調波ノイズ重み付け係数εｐを前記ピッチ周期に基づいて決定するステップと、
知覚ノイズ重み付け関数ＷＨ（ｚ）を前記高調波ノイズ重み付け係数に基づいて決定するステップと、
を備える方法。
【請求項２】
請求項１に記載の方法において、εｐはＤの減少関数である、方法。
【請求項３】
請求項２に記載の方法において、次式が成り立ち、
【数１】
2007513364.tif 000012
上式において、
εｍａｘは前記高調波ノイズ重み付け係数の最大許容値であり、
εｍｉｎは前記高調波ノイズ重み付け係数の最小許容値であり、
Ｄｍａｘは最大ピッチ周期であり、該周期を超えると、前記高調波ノイズ重み付け係数がεｍｉｎに設定され、
Δは前記高調波ノイズ重み付け係数の傾きである、方法。
【請求項４】
高調波ノイズ重み付けをデジタル音声符号器において行なう方法であって、
音声入力ｓ（ｎ）を受信するステップと、
閉ループ・ピッチ遅延（τ）を前記音声入力から求めるステップと、
高調波ノイズ重み付け係数εｐを前記閉ループ・ピッチ遅延に基づいて決定するステップと、
知覚ノイズ重み付け関数ＷＨ（ｚ）を前記高調波ノイズ重み付け係数に基づいて決定するステップと、
を備える方法。
【請求項５】
請求項４に記載の方法において、εｐはτの減少関数である、方法。
【請求項６】
請求項５に記載の方法において、
次式が成り立ち、
【数２】
2007513364.tif 000013
上式において、
εｍａｘは前記高調波ノイズ重み付け係数の最大許容値であり、
εｍｉｎは前記高調波ノイズ重み付け係数の最小許容値であり、
τｍａｘは最大閉ループ・ピッチ遅延であり、該ピッチ遅延を超えると、前記高調波ノイズ重み付け係数がεｍｉｎに設定され、
Δは前記高調波ノイズ重み付け係数の傾きである、方法。
【請求項７】
装置であって、
音声（ｓ（ｎ））を入力として有し、かつピッチ周期（Ｄ）を前記音声に基づいて出力するピッチ解析回路と、
Ｄを入力として有し、かつ高調波ノイズ重み付け係数（εｐ）をＤに基づいて出力する高調波ノイズ係数発生器と、
εｐを入力として有し、かつｅ（ｎ）がｓ（ｎ）とｓ（ｎ）の推定値との間の差に基づいて表わされる場合に、εｐを使用して重み付け誤差信号ｅ（ｎ）を生成する知覚誤差重み付けフィルタと、
を備える装置。
【請求項８】
装置であって、
閉ループ・ピッチ遅延（τ）を入力として有し、かつ高調波ノイズ重み付け係数（εｐ）をτに基づいて出力する高調波ノイズ係数発生器と、
εｐを入力として有し、かつｅ（ｎ）がｓ（ｎ）とｓ（ｎ）の推定値との間の差に基づいて表わされる場合に、εｐを使用して重み付け誤差信号ｅ（ｎ）を生成する知覚誤差重み付けフィルタと、
を備える装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は概して信号圧縮システムに関し、より詳細には、符号励起線形予測（ＣＥＬＰ：Code Excited Linear Prediction）型音声符号化システムに関する。
【背景技術】
【０００２】
デジタル音声信号及びオーディオ信号の圧縮は公知である。圧縮では通常、信号を通信チャネルを通して効率的に送信するか、あるいは圧縮信号を固体メモリ装置またはコンピュータ・ハード・ディスクのようなデジタル・メディア装置に保存する必要がある。多くの圧縮（または「符号化」）技術があるが、デジタル音声符号化に非常に一般的に用いられてきた一つの方法が符号励起線形予測（ＣＥＬＰ）であり、この方法は、一連の「合成的解析（analysis-by-synthesis ）」符号化アルゴリズムの一つである。合成的解析は普通、符号化プロセスを意味し、このプロセスでは、デジタル・モデルのパラメータを使用して一連の候補信号を合成し、これらの候補信号を入力信号と比較して歪みを解析する。次に、最小歪み、または誤差成分を生じさせる一連のパラメータを送信するか、あるいは保存する。一連のパラメータを最終的に使用して元の入力信号の推定値を再構成する。ＣＥＬＰは特定の合成的解析法であり、この解析法では、一つ以上の励起コードブックを使用し、これらのコードブックは基本的に、コードブックからコードブック・インデックスに応じて取り出される複数組みのコードベクトルを含む。これらのコードベクトルは、「試行錯誤（trial and error ）」プロセスにおける音声合成装置への入力として使用され、このプロセスでは、誤差基準を候補コードベクトルの各々に関して評価し、そして最小誤差を生じさせる候補を選択する。
【０００３】
例えば、図１は先行技術によるＣＥＬＰ符号器１００のブロック図である。ＣＥＬＰ符号器１００では、音声サンプルｎ（ｓ（ｎ））を含む入力信号を線形予測符号化（ＬＰＣ：Linear Predictive Coding）解析ブロック１０１に入力し、このブロックにおいて、線形予測符号化を使用して短時間スペクトルの包絡線を推定する。結果として得られるスペクトル・パラメータ（またはＬＰパラメータ）は伝達関数Ａ（ｚ）によって表わされる。スペクトル・パラメータをＬＰＣ量子化ブロック１０２に入力し、ＬＰＣ量子化ブロックはスペクトル・パラメータを量子化して、マルチプレクサ１０８における使用に適する量子化スペクトル・パラメータＡｑを生成する。次に、量子化スペクトル・パラメータＡｑがマルチプレクサ１０８に転送され、そしてマルチプレクサは符号化ビット・ストリームを量子化スペクトル・パラメータ、及び自乗誤差最小化／パラメータ量子化ブロック１０７により決定される一連のパラメータ、τ、β、κ、及びγに基づいて生成する。この技術分野の当業者であれば分かることであるが、τ、β、κ、及びγは閉ループ・ピッチ遅延（closed loop pitch delay ）、適応コードブック・ゲイン、固定コードブック・ベクトル・インデックス、及び固定コードブック・ゲインとしてそれぞれ定義される。
【０００４】
量子化スペクトル・パラメータまたはＬＰパラメータはまた、該当する伝達関数１／Ａｑ（ｚ）を有するＬＰＣ合成フィルタ１０５に自律的に転送される。ＬＰＣ合成フィルタ１０５はまた、合成励起信号ｕ（ｎ）を第１コンバイナ１１０から受信し、入力信号の推定値ｓ＾（ｎ）を量子化スペクトル・パラメータＡｑ及び合成励起信号ｕ（ｎ）に基づいて生成する。合成励起信号ｕ（ｎ）は次のようにして生成される。適応コードブック・コードベクトルｃτは、適応コードブック（ＡＣＢ：Adaptive Codebook ）１０３からインデックス・パラメータτに基づいて選択される。次に、適応コードブック・コードベクトルｃτに対して、ゲイン・パラメータβに基づいて重み付けを行ない、そして重み付け済み適応コードブック・コードベクトルｃτが第１コンバイナ１１０に転送される。固定コードブック・コードベクトルｃｋは、固定コードブック（ＦＣＢ：Fixed Codebook）１０４からインデックス・パラメータκに基づいて選択される。次に、固定コードブック・コードベクトルｃｋに対して、ゲイン・パラメータγに基づいて重み付けを行ない、そしてこの重み付け済み固定コードブック・コードベクトルも第１コンバイナ１１０に転送される。次に、第１コンバイナ１１０は合成励起信号ｕ（ｎ）を、重み付け済み適応コードブック・コードベクトルｃτを重み付け済み固定コードブック・コードベクトルｃｋと合成することにより生成する。（読者が理解し易いように、変数もこれらの変数のｚ変換形で表示している。変数のｚ変換は該当する大文字で表わされる。例えばｅ（ｎ）のｚ変換はＥ（ｚ）と表わされる）。
【０００５】
ＬＰＣ合成フィルタ１０５は入力信号推定値ｓ＾（ｎ）を第２コンバイナ１１２に転送する。第２コンバイナ１１２はまた、入力信号ｓ（ｎ）を受信し、そして入力信号の推定値ｓ＾（ｎ）を入力信号ｓ（ｎ）から減算する。入力信号ｓ（ｎ）と入力信号推定値ｓ＾（ｎ）との間の差を知覚誤差重み付けフィルタ（perceptual error weighting filter)１０６に入力し、知覚誤差重み付けフィルタは知覚重み付け誤差信号ｅ（ｎ）をｓ＾（ｎ）とｓ（ｎ）との間の差、及び重み付け関数ｗ（ｎ）に基づいて、次式が成り立つように生成する。
【０００６】
【数１】
2007513364.tif 000003
次に、知覚重み付け誤差信号ｅ（ｎ）が自乗誤差最小化／パラメータ量子化ブロック１０７に転送される。自乗誤差最小化／パラメータ量子化ブロック１０７は誤差信号ｅ（ｎ）を使用して、入力信号ｓ（ｎ）の最良の推定値ｓ＾（ｎ）を生成する最適な一連のτ、β、κ、及びγを求める。
【０００７】
図２は、符号器１００からの送信を受信する先行技術による復号器２００のブロック図である。この技術分野の当業者であれば理解できることであるが、符号器１００が生成する符号化ビット・ストリームを復号器２００のデマルチプレクサが使用して、最適な一連のパラメータ、すなわちτ、β、κ、及びγを、符号器１００が実行する合成プロセスと同じプロセスで復号化する。従って、符号器１００が生成する符号化ビット・ストリームを復号器２００が誤差を生じることなく受信する場合、復号器２００が出力する音声ｓ＾（ｎ）は、符号器１００が生成する入力音声推定値ｓ＾（ｎ）の正確な複製として再構成することができる。
【０００８】
図１に戻って同図を参照すると、重み付けフィルタＷ（ｚ）はヒトの耳が周波数マスク現象を示すことを利用して、信号及びノイズの周波数が近接するとした場合に、同時に生じるノイズが相対的に強い信号によってマスクされるようにしていることが示される。ＩＥＥＥ  Ｔｒａｎｓ．  Ｏｎ  Ｖｅｈｉｃｕｌａｒ  Ｔｅｃｈｎｏｌｏｇｙ，  ｐｐ．  ８０８−８１６，  Ａｕｇ．  １９９４に掲載されたＳａｌａｍｉ  Ｒ．，  Ｌａｆｌａｍｍｅ  Ｃ．，  Ａｄｏｕｌ  Ｊ−Ｐ，  Ｍａｓｓａｌｏｕｘ  Ｄ．，による「パーソナル通信システム用トール・クオリティ８Ｋｂ／ｓ音声符号器（A toll quality 8 Kb/s speech coder for personal communication system）」と題された論文に記載されているように、Ｗ（ｚ）はＬＰＣ係数ａｉから生成され、次式により与えられる。
【０００９】
【数２】
2007513364.tif 000004
上式において、次式
【００１０】
【数３】
2007513364.tif 000005
により表される関係があり、ｐはＬＰＣの次数である。重み付けフィルタはＬＰＣスペクトルから生成されるので、重み付けフィルタは「スペクトル重み付け」とも呼ばれる。
【００１１】
上述の手順では、信号周期性も基本周波数及び基本周波数の倍数の周波数でのスペクトル・ピークに影響することを考慮に入れていない。種々の技術が、これらの基本周波数の高調波のノイズをマスクする方法を利用することを提案している。例えば、「デジタル音声符号器及び高調波ノイズ重み付けを利用する方法（Digital speech coder and method utilizing harmonics noise weighting ）」と題されたＧｅｒｓｏｎ及びＪａｓｉｕｋの米国特許第５，５２８，７２３号、及びＰｒｏｃ．  ＩＥＥＥ  ＩＣＡＳＳＰ，  ｐｐ．  ２０５−２０８，  １９９３に掲載されたＧｅｒｓｏｎ  Ｉ．  Ａ．，  Ｊａｓｉｕｋ  Ｍ．Ａ．，による「ＣＥＬＰ型音声符号器の性能を向上させる技術（Techniques for improving the performance of CELP type speech coders ）」と題された論文においては、重み付けフィルタにおいて高調波ノイズにマスクをかける方法が提案されている。上記文献が示しているように、高調波ノイズ重み付けを、スペクトル重み付けフィルタに高調波ノイズ重み付けフィルタＣ（ｚ）によって変形を加えることにより取り入れ、この高調波ノイズ重み付けは次式で表わされる。
【００１２】
【数４】
2007513364.tif 000006
上式において、Ｄはピッチ周期またはピッチ・ラグ、またはピッチ遅延に対応し、ｂｉはフィルタ係数であり、そして０≦εｐ＜１は高調波ノイズ重み付け係数である。高調波ノイズ重み付けを取り入れた重み付けフィルタは次式により表わされる。
【００１３】
【数５】
2007513364.tif 000007
【発明の開示】
【発明が解決しようとする課題】
【００１４】
高調波ノイズ重み付けの大きさは通常、積εｐｂｉに依存する。ｂｉは遅延に依存するので、高調波ノイズ重み付けの大きさは遅延の関数である。上記の先行技術文献は、高調波ノイズ重み付け係数（εｐ）の異なる値を異なる所定時間に使用することができる、すなわち、εｐを時間変化パラメータ（例えば、サブフレーム毎に変化することができる）とすることができることを提案しているが、先行技術はεｐを選択する方法を提示していない。従って、高調波ノイズ重み付けを、最適かつ動的にεｐの適切な値を決定するデジタル音声符号器において行なって、高調波ノイズ重み付けの大きさの最適化を可能にする方法及び装置が必要になる。上記の先行技術文献は、高調波ノイズ重み付け係数（εｐ）の異なる値を異なる時間で使用することができる（例えば、εｐはサブフレーム毎に変化することができる）ことを提案しているが、先行技術はεｐを変える方法を提示していないか、あるいはこのような方法がいつ、またはどのようにして利点をもたらすかについて提案していない。従って、高調波ノイズ重み付けを、最適かつ動的にεｐの適切な値を決定するデジタル音声符号器において行って、知覚重み付け全体を向上させることができる方法及び装置が必要になる。
【課題を解決するための手段】
【００１５】
高調波ノイズ重み付け（ＨＮＷ：harmonic noise weighting）の大きさを最適化することができるように高調波ノイズ重み付け係数（εｐ）の値を選択する必要性に対処するために、高調波ノイズ重み付けをデジタル音声符号器において行なう方法及び装置を本明細書に開示する。動作状態の間、受信音声を解析してピッチ周期を求める。次に、ＨＮＷ係数をピッチ周期に基づいて選択し、知覚ノイズ重み付けフィルタ（Ｃ（ｚ））を高調波ノイズ重み付け（ＨＮＷ）係数（εｐ）に基づいて決定する。ピッチ周期（Ｄ）が長い場合、基本周波数の高調波のピークは非常に近接するので、隣接高調波成分の谷間は隣接ピーク領域をマスクする位置に配置することができる。従って、大きな高調波ノイズ重み付け係数を大きな値のＤに関して付与する必要がない。
【００１６】
ＨＮＷ係数はピッチ周期の関数であるので、より高精度のノイズ重み付けを行なうことができ、従って音声歪みは聴取者にとってはさほど気にならない。
本発明は高調波ノイズ重み付けをデジタル音声符号器において行なう方法を含む。本方法は、音声入力ｓ（ｎ）を受信するステップと、ピッチ周期（Ｄ）を音声入力から求めるステップと、高調波ノイズ重み付け係数εｐをピッチ周期に基づいて決定するステップと、を含む。次に、知覚ノイズ重み付け関数ＷＨ（ｚ）を高調波ノイズ重み付け係数に基づいて決定する。
【００１７】
本発明は更に、高調波ノイズ重み付けをデジタル音声符号器において行なう方法を含む。本方法は、音声入力ｓ（ｎ）を受信するステップと、閉ループ・ピッチ遅延（τ）を音声入力から求めるステップと、高調波ノイズ重み付け係数εｐを閉ループ・ピッチ遅延に基づいて決定するステップと、を含む。次に、知覚ノイズ重み付け関数ＷＨ（ｚ）を高調波ノイズ重み付け係数に基づいて決定する。
【００１８】
本発明は更にピッチ解析回路を備える装置を含み。ピッチ解析回路は音声ｓ（ｎ）を入力として有し、ピッチ周期（Ｄ）を音声に基づいて出力する。装置は高調波ノイズ重み付け係数発生器を備え、発生器はＤを入力として有し、高調波ノイズ重み付け係数（εｐ）をＤに基づいて決定する。装置は更に知覚ノイズ重み付けフィルタを備え、フィルタはεｐを入力として有し、εｐを使用して重み付け誤差信号ｅ（ｎ）を生成する。ｅ（ｎ）はｓ（ｎ）とｓ（ｎ）の推定値との間の差に基づいて表わされる。
【００１９】
最後に、本発明は高調波ノイズ重み付け係数発生器を備える装置を含み、高調波ノイズ重み付け係数発生器は閉ループ・ピッチ遅延（τ）を入力として有し、高調波ノイズ重み付け係数（εｐ）をτに基づいて出力する。装置は知覚誤差重み付けフィルタを備え、知覚誤差重み付けフィルタはεｐを入力として有し、εｐを使用して重み付け誤差信号ｅ（ｎ）を生成する。ｅ（ｎ）はｓ（ｎ）とｓ（ｎ）の推定値との間の差に基づいて表わされる。
【発明を実施するための最良の形態】
【００２０】
次に、同様の参照番号が同様の構成要素を指す添付の図を参照すると、図３は、本発明の好適な実施形態によるＣＥＬＰ符号器３００のブロック図を示している。図示されているように、ＣＥＬＰ復号器３００は先行技術に示されるＣＥＬＰ復号器に類似するが、ピッチ解析回路３１１及びＨＮＷ係数発生器３０９が追加されている点が異なる。更に、知覚誤差重み付けフィルタ３０６は、ＨＮＷ係数をＨＮＷ係数発生器３０９から受信するように構成されている。符号器３００の動作は次のようにして行なわれる。
【００２１】
入力音声ｓ（ｎ）をピッチ解析回路３１１に入力し、この回路では、ｓ（ｎ）を解析してピッチ周期（Ｄ）を求める。この技術分野の当業者であれば理解できることであるが、ピッチ周期（更に、ピッチ・ラグ、遅延、またはピッチ遅延と呼ばれる）は通常、タイム・ラグであり、このタイム・ラグで、過去の入力音声が現在の入力音声と最大の相関を有する。
【００２２】
一旦、ピッチ周期（Ｄ）が求まると、ＤをＨＮＷ係数発生器３０９に入力し、この発生器では、特定の音声に関するＨＮＷ係数（εｐ）を決定する。上述したように、高調波ノイズ重み付け係数はピッチ周期Ｄの関数として動的に変化する。高調波ノイズ重み付けフィルタは次式により表わされる。
【００２３】
【数６】
2007513364.tif 000008
上述したように、大きな値のＤに関して小さな高調波ノイズ重み付け（Ｃ（ｚ））を行なうことが望ましい。εｐをＤの減少関数として選択する（式（７）参照）ことにより、大きな値のピッチ遅延に関して小さな大きさの高調波ノイズ重み付けを確実に行なうことができる。εｐ（Ｄ）に関する多くの関数があるが、本発明の好適な実施形態では、εｐ（Ｄ）は式（７）により与えられ、そして図４にグラフとして示される。
【００２４】
【数７】
2007513364.tif 000009
上式において、
εｍａｘは高調波ノイズ重み付け係数の最大許容値であり、
εｍｉｎは高調波ノイズ重み付け係数の最小許容値であり、
Ｄｍａｘは最大ピッチ周期であり、この周期を超えると、高調波ノイズ重み付け係数がεｍｉｎに設定され、
Δは高調波ノイズ重み付け係数の傾きである。
【００２５】
一旦、εｐ（Ｄ）が発生器３０９によって決定されると、εｐ（Ｄ）をフィルタ３０６に供給して重み付けフィルタＷＨ（ｚ）を生成する。上述のように、ＷＨ（ｚ）はＷ（ｚ）及びＣ（ｚ）の積である。誤差ｓ（ｎ）−ｓ＾（ｎ）を重み付けフィルタ３０６に供給して重み付け誤差信号ｅ（ｎ）を生成する。先行技術による符号器におけるように、誤差重み付けフィルタ３０６は誤差信号ｅ（ｎ）を入力信号と推定入力信号との間の差に基づいて生成する、すなわち次式に従って生成する。
【００２６】
【数８】
2007513364.tif 000010
重み付けフィルタＷＨ（ｚ）は、ヒトの耳が周波数マスク現象を示すことを利用して、同時に生じるノイズが、信号及びノイズの周波数が近接すると仮定した場合に、強い方の信号によってマスクされるようにする。ｅ（ｎ）の値に基づいて、自乗誤差最小化／パラメータ量子化回路３０７はτ、κ、γ、βの値を生成し、これらの値はチャネルを伝送されるか、あるいはデジタル・メディア装置に保存される。
【００２７】
上述したように、ＨＮＷ係数はピッチ周期の関数なので、精度の高いノイズ重み付けを行なうことができ、従って音声歪みは聴取者にはさほど気にならない。
図５は符号器３００の動作を示すフロー・チャートである。論理フローはステップ５０１から始まり、このステップでは、音声入力（ｓ（ｎ））をピッチ解析回路３１１が受信する。ステップ５０３では、ピッチ解析回路３１１はピッチ周期（Ｄ）を求め、ＤをＨＮＷ係数発生器３０９に出力する。ＨＮＷ係数発生器３０９はＤを使用して高調波ノイズ重み付け係数（εｐ）をＤに基づいて決定し、そしてεｐを知覚誤差重み付けフィルタ３０６に出力する（ステップ５０５）。最後に、ステップ５０７では、フィルタ３０６はεｐを使用して知覚ノイズ重み付け関数ＷＨ（ｚ）を生成する。
【００２８】
本発明について特定の実施形態を参照しながら特定の形で示し、記載してきたが、この技術分野の当業者であれば、形式及び詳細における種々の変更を実施形態に、本発明の技術思想及び技術範囲から逸脱しない範囲において加え得ることが理解できるであろう。例えば、ＷＨ（ｚ）をεｐに基づいて生成するために特定の公式を与えたが、ＷＨ（ｚ）をεｐに基づいて生成する他の手段を利用することも考えられる。例えば、式（６）のＣ（ｚ）の定義を行なうための加算項は、εｐを乗じる前に更に変更することができる。更に、別の実施形態では、εｐをτに基づいて表わすことができ、この場合、τ（図６参照）を式（７）のＤの代わりに用いる。上述したように、τは閉ループ・ピッチ遅延として定義され、この場合、εｐはτの減少関数である。従って、式（７）は次式のように表わされる。
【００２９】
【数９】
2007513364.tif 000011
上式において、
εｍａｘは高調波ノイズ重み付け係数の最大許容値であり、
εｍｉｎは高調波ノイズ重み付け係数の最小許容値であり、
τｍａｘは最大閉ループ・ピッチ遅延であり、この遅延を超えると、高調波ノイズ重み付け係数がεｍｉｎに設定され、
Δは高調波ノイズ重み付け係数の傾きである。
【図面の簡単な説明】
【００３０】
【図１】先行技術による符号励起線形予測（ＣＥＬＰ）符号器のブロック図。
【図２】先行技術によるＣＥＬＰ復号器のブロック図。
【図３】本発明の好適な実施形態によるＣＥＬＰ復号器のブロック図。
【図４】εｐ対ピッチ・ラグ（Ｄ）の関係を表わすグラフ。
【図５】本発明の高調波ノイズ重み付け方法を取り入れるためにＣＥＬＰ符号器が実行するステップを示すフロー・チャート。
【図６】本発明の別の実施形態によるＣＥＬＰ符号器のブロック図。
【図１】
2007513364.tif 000014
【図２】
2007513364.tif 000015
【図３】
2007513364.tif 000016
【図４】
2007513364.tif 000017
【図５】
2007513364.tif 000018
【図６】
2007513364.tif 000019
【国際調査報告】
2007513364.tif 000020