(57)【要約】  本発明は、媒体信号に埋め込まれた少なくともプレゼンスデータシンボルの識別を可能とする方法、装置、及びコンピュータプログラムに関する。送信された媒体信号を取得する（ステップ２２）。その媒体信号は一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられた、ホスト信号を歪めたバージョンを含む。一組のディザ値のうちの一部のディザ値区間を設ける。各区間は１つのディザ値に対応する（ステップ２６）。各区間について１つのヒストグラムを決定する（ステップ３０）。ヒストグラムは送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、対応する区間中の１つのディザ値を有する。別のヒストグラムを結合し、量子化ステップサイズ推定するために結合したヒストグラムに基づいてリスケーリングファクタを決定する（ステップ３２）。
2007516680.tif 000002
【特許請求の範囲】
【請求項１】
媒体に埋め込まれたデータシンボルがあることを少なくとも識別可能とする方法であって、
送信された媒体信号を取得するステップであって、前記媒体信号は、一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられたホスト信号を歪めたバージョンを含むところのステップと、
前記一組のディザ値内に、それぞれが１つまたは少数のディザ値に対応する複数のディザ値区間を設けるステップと、
各ディザ値区間について１つのヒストグラムを決定するステップであって、ヒストグラムは前記送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、前記対応するディザ値区間中の１つのディザ値を有するところのステップと、
前記区間に対応する別々のヒストグラムを結合して単一のヒストグラムにするステップと、
前記量子化ステップサイズを推定するために、前記結合されたヒストグラムに基づいてリスケーリングファクタを決定するステップと、を有することを特徴とする方法。
【請求項２】
請求項１に記載の方法であって、
前記送信された媒体信号のオフセットを推定するステップと、前記信号から前記推定されたオフセットを取り除くステップと、をさらに有することを特徴とする方法。
【請求項３】
請求項１に記載の方法であって、
データが埋め込まれた前記元の媒体信号を少なくとも近似的に回復するために、前記リスケーリングファクタを用いて前記信号をリスケーリングするステップをさらに有することを特徴とする方法。
【請求項４】
請求項３に記載の方法であって、
前記埋め込まれたデータを検出して取り出すために、前記リスケーリングされた信号を処理するステップをさらに有することを特徴とする方法。
【請求項５】
媒体信号に埋め込まれたデータシンボルがあることを少なくとも識別可能とする装置であって、
一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられたホスト信号を歪めたバージョンを含む送信された媒体信号を取得するように構成された信号取得部と、
前記一組のディザ値内に、それぞれが１つまたは少数のディザ値に対応する複数のディザ値区間を設け、
各ディザ値区間について１つのヒストグラムを決定し、ここでヒストグラムは前記送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、前記対応するディザ値区間中の１つのディザ値を有し、
前記区間に対応する別々のヒストグラムを単一のヒストグラムに結合し、
前記量子化ステップサイズを推定するために、前記結合されたヒストグラムに基づいてリスケーリングファクタを決定するように構成された信号歪み判定部と、を有することを特徴とする装置。
【請求項６】
請求項５に記載の方法であって、
前記信号歪み判定部は、前記送信された媒体信号のオフセットを推定するようにさらに構成され、  前記装置は、前記信号から前記推定されたオフセットを取り除くように構成されたユニットをさらに有することを特徴とする装置。
【請求項７】
請求項５に記載の装置であって、
データが埋め込まれた前記元の媒体信号を少なくとも近似的に回復するため、前記送信された媒体信号で前記リスケーリングファクタにかけるように構成された乗算部をさらに有することを特徴とする装置。
【請求項８】
請求項７に記載の装置であって、
前記埋め込まれたデータを検出または取り出すために、前記リスケーリングされた信号を処理するように構成された透かし検出部をさらに有することを特徴とする装置。
【請求項９】
媒体信号に埋め込まれたデータシンボルがあることを少なくとも識別可能とするコンピュータプログラムであって、
前記プログラムをコンピュータに読み込んだ時、前記コンピュータに、
一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられたホスト信号を歪めたバージョンを含む、送信され取得された媒体信号に対して、
前記一組のディザ値内に、それぞれが１つまたは少数のディザ値に対応する複数のディザ値区間を設け、
各ディザ値区間について１つのヒストグラムを決定し、ここでヒストグラムは前記送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、前記対応するディザ値区間中の１つのディザ値を有し、
前記区間に対応する別々のヒストグラムを単一のヒストグラムに結合し、
前記量子化ステップサイズを推定するために、前記結合されたヒストグラムに基づいて少なくともリスケーリングファクタを決定させることを特徴とするコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、透かしを入れた媒体信号、及び媒体信号に埋め込まれた透かしの検出の分野に関する。本発明は、特に、ホスト信号を量子化することにより媒体信号に埋め込まれたデータシンボルの少なくとも識別を可能とする方法、装置、及びコンピュータプログラム製品に関する。
【背景技術】
【０００２】
コンテント所有者の権利を侵害や詐欺から守るため、媒体信号に透かしをいれることは周知である。
このような透かし方法として、いわゆるＱＩＭ（Quantization  Index  Modulation）透かし法とＤＣ−ＱＩＭ（Distortion-Compensated  Quantization  Index  Modulation）透かし法がある。透かしはサンプリンググリッド上のサンプル値を選択することにより、いわゆるホスト信号に埋め込まれる。透かし埋め込み後にも元の信号と統計的に類似した信号となるように、透かし信号にディザを加えることが一般的である。
【０００３】
しかし、このように透かしを埋め込んだホスト信号をレシーバに送信すると、使用されるスケーリングファクタがリスケーリングされ、スケーリングファクタは量子化ステップサイズと比例する。そうすると、量子化サイズが正しくないので、透かしの検出が困難または不可能となる。このようなスケーリングファクタのリスケーリングは、例えば使用した送信チャネルの特性等、いろいろな理由で起こる。１つの理由は、チャネルがアナログであり、そのゲインファクタが分からないことである。他の理由は、悪意者による信号への攻撃である。
【０００４】
国際特許出願公報第０３／０５２６８９号には、リスケーリングされている可能性がある受信ホスト信号からスケーリングファクタを回復する方法が記載されている。この文献によると、受信信号サンプルは３つの異なるディザ値区間に分割される。次に、各区間における信号サンプルの部分確率密度関数を用いてヒストグラムを作り、この３つのヒストグラムを分析して元のスケーリングファクタ及びＤＣオフセットを求める。このように、分析する前にディザ値を変えたグループで信号サンプルを供給する。つまり、信号サンプルの情報の一部はスケーリングファクタの決定には使用されないことがある。この文献には、スケーリングファクタの回復でパイロット信号を使用することも記載されている。
【０００５】
しかし、このスケーリングファクタの決定のために、異なる信号サンプルの情報をより多く使用することができる、受信媒体信号のスケーリングファクタを決定する代替的な方法があれば有利である。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
よって、本発明の目的の１つは、先行技術において使用するより多い受信信号サンプル中の情報を使用する、受信媒体信号のスケーリングファクタを決定する代替的方法を提供することである。
【課題を解決するための手段】
【０００７】
本発明の第１の態様によると、この目的は、媒体信号に埋め込まれたデータシンボルがあることを少なくとも識別することを可能とする方法により達成される。該方法は、一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられたホスト信号を歪めたバージョンを含む、送信された媒体信号を取得するステップと、前記一組のディザ値内に、それぞれが１つまたは少数のディザ値に対応する複数のディザ値区間を設けるステップと、各ディザ値区間について１つのヒストグラムを決定するステップであって、ここでヒストグラムは前記送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、前記対応するディザ値区間中の１つのディザ値を有するステップと、前記区間に対応する別々のヒストグラムを単一のヒストグラムに結合するステップと、前記量子化ステップサイズを推定するために、前記結合されたヒストグラムに基づいて少なくともリスケーリングファクタを決定させるステップと、を有する。
【０００８】
本発明の第２の態様によると、この目的は、媒体信号に埋め込まれたデータシンボルがあることを少なくとも識別することを可能とする装置によっても達成される。該装置は、
一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられたホスト信号を歪めたバージョンを含む送信された媒体信号を取得するように構成された信号取得部と、前記一組のディザ値内に、それぞれが１つまたは少数のディザ値に対応する複数のディザ値区間を設け、各ディザ値区間について１つのヒストグラムを決定し、ここでヒストグラムは前記送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、前記対応するディザ値区間中の１つのディザ値を有し、前記区間に対応する別々のヒストグラムを単一のヒストグラムに結合し、前記量子化ステップサイズを推定するために、前記結合されたヒストグラムに基づいて少なくともリスケーリングファクタを決定するように構成された信号歪み判定部と、を有する。
【０００９】
本発明の第３の態様によると、この目的は、媒体信号に埋め込まれたデータシンボルがあることを少なくとも識別することを可能とするコンピュータプログラム製品によっても達成される。該製品はコンピュータ読み取り可能媒体を有し、該媒体上に、前記プログラムをコンピュータに読み込んだ時、前記コンピュータに、一定の量子化ステップサイズを用いて量子化することによりデータシンボルが埋め込まれ、そのディザに一組のディザ値が加えられたホスト信号を歪めたバージョンを含む、送信され取得された媒体信号に対して、前記一組のディザ値内に、それぞれが１つまたは少数のディザ値に対応する複数のディザ値区間を設け、各ディザ値区間について１つのヒストグラムを決定し、ここでヒストグラムは前記送信された媒体信号の一組の信号サンプルのすべてのサンプル値に対して決定され、前記対応するディザ値区間中の１つのディザ値を有し、前記区間に対応する別々のヒストグラムを単一のヒストグラムに結合し、前記量子化ステップサイズを推定するために、前記結合されたヒストグラムに基づいて少なくともリスケーリングファクタを決定させる、コンピュータプログラムコード手段、を有する。
【００１０】
請求項２と６は、送信された媒体信号のオフセットの推定にも関する。
【００１１】
請求項３と７は、スケーリングファクタを用いた、送信された媒体信号のリスケーリングに関する。
【００１２】
請求項４と８は、埋め込まれたデータを検出し、これを取り出すための、リスケーリングされた信号の処理に関する。
【００１３】
本発明は、量子化ステップサイズを推定するために、異なる信号サンプル値のヒストグラム中のすべての情報を用いるという利点を有する。ヒストグラムの数は多いので、ヒストグラム間の周期性を見つけるのはさらに容易であり、これにより量子化ステップサイズのロバスト（robust）な決定が可能となる。さらに、本発明は、量子化ステップサイズを推定するためにパイロット信号を受け取る必要はない。
【００１４】
本発明の基本的アイデアは、埋め込まれた情報の検出器が、それが有する、サンプル位置に対応するディザ値に関する知識を最適に利用することである。特に、その検出器は、どのサンプル位置が同じディザ値を有するかに関する知識を使用する。ディザ値に対応するすべてのサンプル位置に対して、（粗い）ヒストグラムを構成する。次に、これらの粗いヒストグラムを結合して単一の振る舞いのよいヒストグラムにして、それからスケーリングファクタを求める。
【００１５】
本発明の上記その他の態様を、以下に説明する実施形態を参照して明らかにし、説明する。
【００１６】
ここで、添付した図面を参照して、本発明をより詳しく説明する。
【００１７】
［実施形態の詳細な説明］
本発明は、受信したホスト信号の透かしを検出してそれを取り出すために、受信した透かしを埋め込んだ信号の元のスケーリングファクタを推定することに関する。本発明をよりよく理解するため、ＱＩＭ（Quantization  Index  Modulation）とＤＣ−ＱＩＭ（Distortion  Compensated  Quantization  Index  Modulation）の基本概念についてここで説明する。
【００１８】
ＱＩＭによると、透かしはサンプリンググリッド上の量子化値を選択することによりホスト信号に埋め込まれる。このアプローチの場合、透かしは次式（１）で表せる：
【００１９】
【数１】
2007516680.tif 000003
ここで、ｑは透かしを埋め込んだ信号、ｓは元の信号、Ｄは量子化ステップサイズ、ｄは区間０≦ｄ＜１にあるディザ値、ｂは透かし情報、Ｎは透かし信号中の情報シンボル数である。実際のアプリケーションでは多くの場合Ｎは２である。これはいわゆる２値アルファベットを示している。
【００２０】
ディザは広い範囲の値を取る。本発明の一実施形態においては、０と１の間の１２８通りの値を取る。ディザ値は固定された疑似乱数変数であり、情報シンボルとは独立である。ディザは後でスケーリングファクタを決定するために使用されるが、これについては後で説明する。透かしを埋め込んだ信号にディザを設ける理由は、透かしを埋め込んだ信号ｑは元の信号ｓとは事実上区別でき、その区別は透かしを埋め込んだ信号と元の信号のヒストグラムを調べて判断できる。
【００２１】
以上、ＱＩＭ透かし法について説明した。この変形がＤＣ−ＱＩＭ法であり、基本的な量子化アプローチが歪み補償部で改良されている。
【００２２】
この方法による透かしを表す式は以下の式（２）の通りである：
【００２３】
【数２】
2007516680.tif 000004
ここで、ａは０と１の間で適当に選択したファクタである。
【００２４】
これらの透かし法に関するより詳細な情報は、国際特許出願公報第０３／０５２６８９号及び論文「Quantization  Index  Modulation:  A  Class  of  Provably  Good  Methods  for  Digital  Watermarking  and  Information  Embedding」（B.  Chen、G.  W.  Wornell著、IEEE  Transactions  on  Information  Theory、Vol.  47,  No.  4,  May  2001）に記載されている。これらの文献はここに参照援用する。
【００２５】
これら２つの方法のいずれかにより透かしを埋め込まれた信号は、チャネルを通じて透かし検出器に送信されると、受信信号ｒ中に埋め込まれたメッセージｂを検出することは非常に容易である。検出は次式（３）により行う：
【００２６】
【数３】
2007516680.tif 000005
すなわち、Ｄ／Ｎで値ｒを割り、リスケーリングされたディザＮｄを引き、その結果を丸め、その丸めた値にモジュロＮ演算をする。
【００２７】
ここで、受信信号ｒは、ノイズｎを含む透かしを埋め込んだ信号ｑに対応すると仮定した。
【００２８】
実際には、ノイズがのっているだけでなく、元の信号がリスケーリングされていることがある。このリスケーリングは使用した送信チャネルの特定に起因し、その送信チャネルがアナログチャネルである場合に起こる。また、リスケーリングは、透かしを埋め込んだ信号に対する悪意者による攻撃の結果であってもよい。本発明は、受信信号のリスケーリングファクタの識別に関し、元のスケーリングを推定し、その後回復した信号の透かしを識別して取り出すために信号をスケーリングすることは簡単である。
【００２９】
変化後の信号と元の信号の関係は次式（４）でモデル化できる：
【００３０】
【数４】
2007516680.tif 000006
ここで、Ｒは乗法的劣化（multiplicative  degradation）であり、ＯはＤＣ値のオフセットを表す。次式（５）により、
【００３１】
【数５】
2007516680.tif 000007
新しい信号ｒ′を決定し、それから埋め込まれた透かしを回復する。
【００３２】
本発明を説明するため図１と図２を参照する。図１はデータ信号の識別を可能とする装置を示し、図２はデータ信号の識別を可能とする方法を示すフローチャートである。装置１０は、受信部１２を有し、チャネル（図示せず）を通して透かしが埋め込まれた媒体信号を受信する（ステップ２２）。チャネルはアナログチャネルでもよく、その場合、受信部はラジオ回路を含む。受信信号ｒは信号歪み判定部１４に送られる（ステップ２４）。信号歪み判定部１４は本発明のコアとなる要素であり、この信号歪み判定部１４において、幾つかのステップが実行される。次にこのステップを説明する。
【００３３】
受信信号ｒは一範囲に値を有する幾つかのサンプルからできており、透かし検出器側でディザ区間が既知であると仮定すると、信号歪み判定部１４は、ディザが取り得る多数のディザ区間（例えば１２８の異なる区間）をまず設定する（ステップ２６）。そして、受信信号の信号サンプル値をディザ区間に従ってグループ分けする。すなわち、１つのディザ区間内のディザを有する信号サンプル値をすべてまとめて１つのグループとする。ここで、各ディザ値区間は１つのディザ値のみに基づく。このように、各ディザ区間にいくつかの異なるサンプルグループが入る（ステップ２８）。各ディザ区間の信号サンプルすべてについてヒストグラムを作り（ステップ３０）、各ディザ区間に対して１つのヒストグラムが作られ、または計算されるようにする。これらのヒストグラムは、単一のヒストグラムに結合される（ステップ３１）。そして、リスケーリングファクタ１／Ｒ及びオフセットＯを、この結合したヒストグラムに基づき推定する（ステップ３２）。推定のしかたは少し後でより詳細に説明する。ここでは、（一般的には粗な）別々のヒストグラムを結合して単一の（振る舞いのよい）ヒストグラムにして、必要なパラメータをそのヒストグラムから求める。次に、決定されたオフセットＯは、信号歪み判定部１４から加算部１６に送られる。加算部１６は信号ｒも受け取る。次に、加算部１６は、信号ｒにオフセット−Ｏを加算、すなわち、受信信号からオフセットＯを減算して（ステップ３４）、このように調節した受信信号を乗算部１８に送る。リスケーリングファクタ１／Ｒは、信号歪み判定部１４により推定された後、乗算部１８に送られる。次に、透かしを埋め込んだ信号の元のレベルまたは元の量子化ステップサイズを回復するために、信号をスケーリングするために、リスケーリングファクタ１／Ｒを調節された受信信号と乗算する（ステップ３６）。次に、このようにリスケーリングされた信号ｒ′を透かし検出器２０に送り、透かし検出を実行する（ステップ３８）。
【００３４】
ここで、図３を参照して、スケーリングファクタをどのように決定するかを説明する。図３は、２５６信号レベル、１２８ディザ値区間、Ｄ＝１６、Ｎ＝２である画像について、本発明の原理に従って求めたいくつかのヒストグラムのマトリックスを示す。ヒストグラムは、このようにそれぞれの列を占める。図から分かるように、結合されたヒストグラムは、左から右に伸びた、少し上を向いた幾つかのラインとなる。この角度は、ディザにより信号のレベルが上がることによるものである。このように、各ヒストグラムによりこれらのラインが得られる。ラインが上がる角度を調べると、受信信号のスケーリングファクタを決定し、そのスケーリングファクタにより信号をリスケーリングすることができる。図中、縦方向にハニング窓を適用した。これがマトリックスの上部と下部が黒くなっている理由である。インデックス１から１２８まで、位相がＮＤ＝２Ｄだけシフトしている。この図から、乗法的劣化（multiplicative  degradation）が約０．９５であることが分かる。さらに、各ディザ区間には明りょうな周期性がある。
【００３５】
数学的には、図１０のヒストグラムは、次式（６）で表すことができる：
【００３６】
【数６】
2007516680.tif 000008
ここで、ｈ（ｋ，ｎ）はディザ区間ｋのサンプルｎの頻度であり、Ｋはディザ区間の数、Ｏはオフセット、Ｒはスケーリングファクタである。マトリックスの各列は他の列をシフトしたものであり、信号の周期性は周波数領域で容易に回復できるので、高速フーリエ変換（ＦＦＴ）を縦軸ｎに沿って適用するが、その前に、ラップアラウンド（wrap  around）効果を減らすために窓関数ｗ（ｎ）を適用する。
【００３７】
窓関数を適用すると、次式（７）となる：
【００３８】
【数７】
2007516680.tif 000009
ここで、ｗ（ｎ）は十分速く崩壊し、近似的にｈ′（ｋ＋ｍ，ｎ）＝ｈ′（ｋ，ｎ−ｍ＊Ｄ／ｋ）、すなわち、窓関数適用後のシフト分散となると仮定した。次に、次式（８）により、縦方向にフーリエ変換する：
【００３９】
【数８】
2007516680.tif 000010
ここで、一組の信号値のサイズはＭである。それゆえ、図３の各列について、フーリエ変換の位相は異なるが、振幅は異ならないことが分かる。（すべてのヒストグラムをそれぞれ継ぎ足した）結合されたヒストグラムｈ（ｎ）のＨ′（ｖ）は、変数ｋにわたって次式の｜Ｈ′（ｋ，ｖ）｜を平均することにより推定できる：
【００４０】
【数９】
2007516680.tif 000011
これは、図４に概略を示した。図４は｜Ｈ′（ｖ）｜をｖの関数として示している。この図から、｜Ｈ′（ｖ）｜が強いピークを有し、そのピークは元のヒストグラムの周期的振る舞いに対応し、ピーク間の距離の周期はＲＤであることが分かる。図中の最初の強いピークは、基本周期ＲＤに対応し、周波数領域では基本周波数Ｍ＊Ｎ／（Ｒ＊Ｄ）である。図中、この周波数は３３と推定され、これからファクタＲ＝０．９６と決定することができる。
【００４１】
この基本周波数は、オフセットＯを決定するためにも使用できる。この基本周波数の（ｋの関数としての）位相的振る舞いは、次式（１０）で与えられる：
【００４２】
【数１０】
2007516680.tif 000012
図５には、この位相関係をプロットし、最良の線形近似（直線）も示した。最小のディザインデックスから最大のディザインデックスまでで、２つの完全なサイクルがカバーされていることが明らかに分かる。オフセット値は０．１７＊Ｄと推定される。
【００４３】
本装置中のユニットは、ソフトウェアでもハードウェアでも提供してもよい。ソフトウェアの場合、信号歪み判定部は、加算部及び乗算部とともに、１つ以上のプロセッサと、本発明による方法を実行するプログラムコードを含む対応するプログラムメモリとで提供してもよい。プログラムコードは、コンピュータプログラム製品で提供してもよい。図６にＣＤ−ＲＯＭディスク４０であるコンピュータプログラム製品の一例を示した。これは一例に過ぎず、その他いろいろなタイプのコンピュータプログラム製品も可能である。さらに、プログラムコードは、例えばインターネットを介して、サーバからダウンロードしてもよい。
【００４４】
本発明により幾つかの利点が得られる。量子化ステップサイズを推定するために、異なる信号サンプル値についてのヒストグラム中のすべての乗法を使用する。ヒストグラムの数は多いので、ヒストグラム間の周期性を見つけるのはさらに容易であり、これにより量子化ステップサイズ及びオフセットのロバスト（robust）な決定が可能となる。さらに、本発明は、量子化ステップサイズを推定するためにパイロット信号を受け取る必要はない。さらに、本方法は、スケーリングファクタがサンプル値の範囲によって一定でない場合にも容易に拡張することができる。
【００４５】
本発明をいろいろと変形することができる。媒体信号を本装置が受信する必要はない。それよりも、本装置がメモリから信号をフェッチしてもよい。さらに、本装置は、透かし検出器を含んでいる必要はない。この検出は他の装置で実行することができる。使用するディザ値の数は、変化してもよいし、必ずしも１２８に限定されてはいない。それゆえ、本発明は添付した特許請求の範囲によって限定されるだけである。
【図面の簡単な説明】
【００４６】
【図１】本発明による装置を示すブロック図である。
【図２】本発明による方法を示すフローチャートである。
【図３】本発明による方法で計算される、受信画像信号を示すヒストグラムである。
【図４】図３のヒストグラムのフーリエ変換を示すグラフである。
【図５】図４の変換されたヒストグラムの位相関係を示すグラフである。
【図６】本発明による方法を実行するコンピュータプログラムコードを含むコンピュータプログラム製品を示す図である。
【図１】
2007516680.tif 000013
【図２】
2007516680.tif 000014
【図３】
2007516680.tif 000015
【図４】
2007516680.tif 000016
【図５】
2007516680.tif 000017
【図６】
2007516680.tif 000018
【国際調査報告】
2007516680.tif 000019
2007516680.tif 000020