(57)【要約】【課題】システムが、表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求める。【解決手段】表示面は任意の形状および姿勢を有することができる。表示面上で座標が既知である位置が特定される。各位置は、表示面のスルーホールに設置される光ファイバによって光センサと光学的に結合される。既知の較正パターンが投影される間に、各較正パターンについて各位置において光の強度を直接検知する。強度は、位置とプロジェクタの出力画像の画素との間の対応関係を求めるために用いられ、投影画像を表示面に合わせてワーピングできるようにする。【選択図】図１
2007536766.tif 000002
【特許請求の範囲】
【請求項１】
任意形状の表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求める方法であって、
既知の較正パターンセットを前記表示面に投影すること、
各較正パターンについて前記表示面上の複数の位置のそれぞれにおいて光の強度を直接検知することであって、位置毎に１つの別個の光センサが関連付けられている、前記検知すること、および
前記位置における強度を相関させ、前記複数の位置と前記プロジェクタの出力画像の画素との間の対応関係を求めること
を含む、任意形状の表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求める方法。
【請求項２】
各位置は既知の座標を有する、
請求項１記載の方法。
【請求項３】
前記較正パターンはＧｒａｙコードの形態である、
請求項１記載の方法。
【請求項４】
前記対応関係は、前記プロジェクタのパラメータを求めるために用いられる、
請求項１記載の方法。
【請求項５】
前記パラメータは、前記プロジェクタの内部および外部パラメータならびに非線形歪みを含む、
請求項４記載の方法。
【請求項６】
前記対応関係に応じて前記プロジェクタに対する入力画像をワーピングすること、および
前記ワーピングした入力画像を前記表示面に投影し、歪みがないように見せることをさらに含む、
請求項１記載の方法。
【請求項７】
前記プロジェクタは平坦な前記表示面と無造作に位置合わせされる、
請求項１記載の方法。
【請求項８】
前記表示面は平坦である、
請求項１記載の方法。
【請求項９】
前記表示面は二次面である、
請求項１記載の方法。
【請求項１０】
鑑賞者および前記プロジェクタは前記表示面に対して同じ側に位置する、
請求項１記載の方法。
【請求項１１】
前記表示面は平坦であり、位置の数は４である、
請求項８記載の方法。
【請求項１２】
前記光センサはフォトトランジスタである、
請求項１記載の方法。
【請求項１３】
前記光センサは、光ファイバによって対応する前記位置と結合される、
請求項１２記載の方法。
【請求項１４】
前記強度は０または１に量子化される、
請求項１記載の方法。
【請求項１５】
前記対応関係に応じて前記プロジェクタに対する入力画像シーケンスをワーピングすること、および
前記ワーピングした入力画像シーケンスを前記表示面に投影し、映像として歪みがないように見せることをさらに含む、
請求項１記載の方法。
【請求項１６】
前記表示面および前記プロジェクタは、前記対応関係を求め、画像シーケンスをワーピングし、前記ワーピングした入力画像シーケンスを投影する間、互いに対して移動する、
請求項１５記載の方法。
【請求項１７】
前記表示面は実世界のオブジェクトの３Ｄモデルの外面である、
請求項１記載の方法。
【請求項１８】
前記表示面は前記３Ｄモデルを乗せる背景幕を含む、
請求項１記載の方法。
【請求項１９】
前記光は赤外線である、
請求項１記載の方法。
【請求項２０】
各較正画像はペアとして投影され、前記ペアの２番目の画像は前記較正画像の反転である、
請求項１記載の方法。
【請求項２１】
前記対応関係は、前記プロジェクタを再配置するために用いられる、
請求項１記載の方法。
【請求項２２】
前記対応関係は、前記表示面を変形させるために用いられる、
請求項１記載の方法。
【請求項２３】
任意形状の表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求めるシステムであって、
座標が既知である複数の位置を有する表示面と、
複数の既知の較正パターンと、
各較正パターンについて前記表示面上の前記複数の位置のそれぞれにおいて光の強度を直接検知する手段と、
前記位置における強度を相関させて、前記複数の位置と前記プロジェクタの出力画像の画素との間の対応関係を求める手段と
を備える、任意形状の表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求めるシステム。
【請求項２４】
各位置は、光ファイバによって別個の光センサと光学的に結合される、
請求項２３記載のシステム。
【請求項２５】
前記光ファイバは、前記表示面のスルーホールに配置される、
請求項２４記載のシステム。
【請求項２６】
任意形状の表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求める方法であって、
前記表示面に投影される複数の較正パターンの各々について前記表示面上の複数の位置のそれぞれにおいて光の強度を直接検知することであって、位置毎に１つの別個の光センサが関連付けられている、前記検知すること、および
前記位置における強度を相関させ、前記複数の位置と前記プロジェクタの出力画像の画素との間の対応関係を求めること
を含む、任意形状の表示面上の位置とプロジェクタの出力画像の画素との間の対応関係を求める方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、包括的にはプロジェクタの較正に関し、特に任意形状の表示面に対するプロジェクタの較正に関する。
【背景技術】
【０００２】
ポータブルなデジタルプロジェクタが今や普及している。これらのプロジェクタは、大型の画像および映像を表示することができる。通常、プロジェクタは卓上に位置付けられるか、投影ブースに配置されるか、あるいは天井に取り付けられる。
【０００３】
従来技術において、歪みのない画像を生成するには、プロジェクタの光軸は平坦な表示面に直角でなければならなかった。また、水平な画像を得るには、プロジェクタの横軸は水平でなければならなかった。上記の制約が満たされても、投影環境の物理的制約によっては、投影画像を投影面の所定のターゲット画像領域と完全に位置合わせすることは依然として難しいか、あるいは不可能でさえある。プロジェクタを無造作に設置した場合、画像の補正が必要となる。
【０００４】
平坦な表示面に対する完全な補正は、歪みを最小にするために３つの位置自由度、２つのスカラ自由度、および３つの回転自由度を考慮する必要がある。これらの補正は、表示面が任意の多様体である場合には不十分である可能性がある。以下において、多様体という用語は特に、３次元で任意の形状および姿勢を有する位相接続面（topological connected surface）を指す。姿勢とは向きおよび位置を意味する。
【０００５】
投影する画像を歪めて、投影画像を正確に位置合わせされ歪みがないように見せることが可能である。しかしこれは、プロジェクタを表示面に対して入念に較正することを要求する。この較正プロセスは、手動で行うと時間がかかる退屈なものとなる可能性があり、高品質な画像を維持するためには頻繁に行わなければならない。プロジェクタまたは表示面のいずれかまたは両方が投影中に動く動的な表示環境の場合、これは極めて難しい。
【０００６】
ほとんどの従来技術の自動較正技法は、補正できる自由度の数が著しく制限されており、通常は１つまたは２つのキーストン補正度のみである。それらの技法はまた、平坦な表示面に制限されている。位置、サイズ、回転、キーストン歪みならびに不規則な表面を自動補正することができる従来技術の技法は、部屋、表示面、および較正カメラの絶対的または相対的な幾何学的形状データの知識に頼っている。カメラを較正に用いる場合、表示面は、較正パターンをカメラへ反射するために反射性でなければならない。いくつかの技法は、プロジェクタにチルトセンサを設置する改良を必要とする。
【発明の開示】
【発明が解決しようとする課題】
【０００７】
このような技法の欠点として、幾何学的な較正データが入手できない場合あるいは場所で、またはプロジェクタに関係のない変更（表示面の位置変更あるいは形状変更または較正カメラの交換等）が行われる場合にプロジェクタを使用できないことがある。表示面が非反射性である場合、あるいは表示面が高反射性である場合（紛らわしい鏡面ハイライトを生じる）、カメラベースの較正システムは機能しない。また、カメラベースのシステムでは、カメラ画像の画素を投影画像中の対応画素に相関させることが難しい。
【０００８】
したがって、任意形状の表面に対してプロジェクタを全自動で較正する方法が必要とされている。
【課題を解決するための手段】
【０００９】
本発明は、任意形状の表示面に対してプロジェクタを較正する方法およびシステムを提供する。この較正は、プロジェクタの位置および回転、画像サイズ、およびキーストン歪み、ならびに非平坦な表面の幾何学的形状を補正する。
【００１０】
本発明は、おそらくは任意形状の表示面上の位置とプロジェクタ画像の画素との間の対応関係を求める方法およびシステムを提供する。例えば、本システムは、オブジェクトにおいて、プロジェクタの左側部分によって照明される部分と、プロジェクタの右側部分によって照明される部分とを分類するために用いることができる。
【００１１】
本発明によるシステムは、表示面またはその付近に取り付けられる別個の光センサを用いる。本方法は、表面に直接投影される光を測定する。これは、表面から反射された光を間接的に測定し、それによりさらなる混乱をもたらすカメラベースのシステムとは区別される。また、各センサは、投影画像の１画素に対応する。カメラベースのシステムでは、異なる光学特性、異なる幾何学的形状、異なる解像度、ならびに異なる内因性および外因性パラメータを少なくとも含むいくつかの理由から、カメラ画素とプロジェクタ画素との間の対応関係を求めることが難しい。
【００１２】
個々のセンサが、投影画像の強度を各位置において直接測定する。１つまたは複数の投影パターンを用いて、本システムは、プロジェクタ画像のどの画素がどの検知位置を照明しているかを推定する。
【００１３】
表示面の２Ｄまたは３Ｄ形状および幾何学的形状が既知であり、かつこの幾何学的形状内の光センサの位置が既知である場合、どのプロジェクタ画素がどのセンサを照明するかに関する情報を用いて、表示面に対してプロジェクタを較正することができる。
【００１４】
取得した較正パラメータを用いて、投影する入力画像を歪め、投影された出力画像を、表示面上で歪みがないように見せる。較正パラメータはまた、表示面に対するプロジェクタの姿勢を求める、内部および外部の幾何学的パラメータを求める、プロジェクタと表示面の間の距離を求める、幾何学的形状が既知である表示面に対するプロジェクタ光線の入射角を求める、表面領域をプロジェクタによって照明されるセグメントと照明されないセグメントとに分類する、プロジェクタの接線方向（radial）歪みを計算する、複数のプロジェクタからの表示面において重複する画像間の関係を求める、および表示面の変形を求めるといった他の目的にも用いることができる。
【発明を実施するための最良の形態】
【００１５】
システム構造
図１に示すように、プロジェクタ１００が平坦な表示面１０１と無造作に位置合わせされる。ここで、鑑賞者（viewer）１０およびプロジェクタ１００は表示面１０１に対して同じ側に位置する。したがって、投影画像は表示面によって鑑賞者へ反射される。無造作な位置合わせにより、プロジェクタの出力画像すなわち映像１０２は、表示面の所望の画像領域１０３と完全には一致しない可能性がある。したがって、入力画像１１０を歪ませて、出力画像として投影されたときに画像領域１０３に合うようにすることが必要である。
【００１６】
したがって、表示面１０１は、２Ｄまたは３Ｄの座標が既知である４つの位置１０４を含む。表示面１０１のサイズおよび位相に応じてさらなる位置を用いることもできることに留意すべきである。任意の投影角度および平坦な表示面について出力画像を矩形の画像領域１０３に合わせるために必要な位置の最小数は４である。
【００１７】
光センサは、既知の位置１０４における光エネルギーの強度を直接測定する。これは、投影画像が表示面によって反射された後でその画像を間接的に測定するカメラベースのシステムと対照的である。直接測定にはいくつかの利点がある。すなわち、カメラベースのプロジェクタの較正と異なり、本システムは、より複雑な幾何学的形状を有する反射光に基づく強度の測定値を処理する必要がない。
【００１８】
一実施形態において、センサは、表面またはその付近の位置１０４に取り付けられるフォトダイオードまたはフォトトランジスタである。別法として、図５に示すように、光センサ５０１は、光ファイバ５０２によって表面の位置１０４と結合される。表面は、スルーホール５０３を含み、光ファイバ５０２のための光路またはルートを提供する。スルーホールは、直径１ミリメートル以下とすることができる。極細の光ファイバを作成する方法はよく知られている。これは、検知位置のサイズをプロジェクタ画素のサイズ以下まで小さくするのに役立つ。本発明の目的で、各検知位置は実質的に、投影される出力画像の画素に対応する。本実施形態は、プロジェクタ１００によって拡大（augment）される小型の３Ｄモデルを取り付ける（instrument）ためにも役立つ。
【００１９】
位置１０４は、画像領域と位置の間の幾何学的関係が既知である限り、画像領域１０３から独立させることができる。これは、表面が平坦であるかまたはパラメータによって定義される、例えば表面が二次面（二次曲面）または他の高次面（高次曲面）（パラメータによって記述できない表面を含む）である場合は簡単である。
【００２０】
較正モジュール（プロセッサ）１０５は、センサ５０１の各々からセンサデータを取得する。好適な実施形態において、センサデータは、Ａ／Ｄ変換後に、０および１のビットに量子化される。０は光が検知されなかったことを示し、１は光が検知されたことを示す。これを可能にするために、光強度を閾値処理することができる。利点として、２値の強度読み取り値は、周囲の背景照明にあまり影響を受けない。ただし、強度は、グレースケールで測定することもできることを理解すべきである。本明細書中に記載する様々な構成要素間のリンクは有線または無線とすることができる。較正モジュールは、ＰＣまたはラップトップコンピュータの形態とすることができる。
【００２１】
図４に示すように、較正モジュール１０５はまた、較正パターン４０１および４０２のセットを生成してプロジェクタ１００へ送出することができる。これらのパターンは以下でより詳細に説明する。較正パターンは、表示面１０１および位置１０４上に投影される。各パターンについて各位置で測定された光強度に基づいて、較正モジュール１０５は、映像処理モジュール１０６に伝えられるワーピング関数（Ｗ）１１１の較正パラメータを求める。較正パラメータは、プロジェクタの内部および外部（ときには内因性および外因性としても知られる）パラメータおよび非線形歪みを反映する。
【００２２】
映像処理モジュール１０６は、映像源１０７によって生成された入力画像１１０を歪め、出力画像１０２を表示面１０１に投影するときにその歪みをなくし、画像領域１０３と位置合わせする。用途によっては、較正パラメータおよびワーピング関数を映像源１０７に直接伝えることが有用である場合がある。
【００２３】
較正モジュール１０５、映像処理モジュール１０６、および映像源１０７、ならびにプロジェクタ１００は、より少ない数の別個の構成要素、例えばプロジェクタサブアセンブリを有する１つのプロセッサモジュールに統合することができる。光センサおよび画像生成ハードウェアを除いて、本システムの機能の大部分はソフトウェアで実施することができる。しかし、ソフトウェアは全てハードウェア回路で実施することもできる。
【００２４】
図２は、複雑で非平坦な画像領域１０３、例えば自動車モデルの外面を示す。このモデルは、フルサイズ、または拡大縮小したものとすることができる。好適な実施形態において、このモデルは、プラスチックまたは紙で作られ、高範囲の色をレンダリングするために白で塗られたモデルカーである。このモデルは、「路面」と「景色」を形成する背景幕の正面に置くことができる。背景幕にもセンサを取り付けることができる。この目的は、外面を実際に塗り直すことなく、モデルを様々な配色で見せることである。背景幕は、車が景色の中の道路を走っているように見えるように照明されることができる。したがって、潜在的な顧客は、購入を決める前に、シミュレートされた環境でモデルを見ることができる。この場合、４つより多くの検知位置を用いる。任意の投影角度および非平坦な表示面について出力画像を表示領域に合わせるために必要な位置の最小数は６である。
【００２５】
本発明は、光学検知位置１０４を含む平坦および非平坦な表示面１０１に対してプロジェクタを較正することを可能にする。較正システムは、画像の位置合わせおよび歪みを補償して、投影される出力画像１０２を表示面の画像領域１０３に合わせることができる。
【００２６】
較正方法
図３は、本発明による較正方法を示す。較正パターン４０１および４０２のセットを順次投影する３００。これらのパターンは、一意の光エネルギーシーケンスを検知位置１０４へ送出する。センサはセンサデータ３１１を取得する３０１。センサデータを復号化して３０２、位置１０４の座標データ３１２を求める。座標データは、ワーピング関数３１３を計算する３０３ために用いられる。ワーピング関数は、入力画像をワーピングして、歪めた出力画像３１４を生成するために用いる。次にこの歪めた画像３１４を投影して３０５、表示面１０１の画像領域１０３と位置合わせすることができる。歪めた画像は、位置の座標から直接生成することもできることに留意すべきである。
【００２７】
較正パターン
図４に示すように、好ましい較正パターン４０１および４０２は、Grayに対して１９５３年３月に発行された米国特許第２，６３２，０５８号に記載される一連のバイナリ符号化マスクに基づく。これらはＧｒａｙコードとして知られる。Ｇｒａｙコードは、機械位置エンコーダに多く用いられる。利点として、Ｇｒａｙコードは、１ビットのみに影響する僅かな位置変化を検出することができる。従来のバイナリコードを用いた場合、最高ｎビットまで変化する可能性があるため、センサ要素間の僅かな位置合わせ誤差が大きく誤った読み取り値を生じる可能性がある。Ｇｒａｙコードにはこの問題がない。Ａ、Ｂ、Ｃ、Ｄ、Ｅとして示す最初の５つのレベル４００は、後の各パターンと前のパターンとの間の関係を、縦の間隔がより細く分割されたものとして示す。４００の５つのレベルは、右側の５つの画像ペア（Ａ、Ｂ、Ｃ、Ｄ、Ｅとして示す）とそれぞれ関連する。各画像ペアは、符号化方式を用いて像平面の横軸４０１および縦軸４０２を分割する方法を示す。この細分化プロセスは、各ビンのサイズがプロジェクタ画素の解像度よりも小さくなるまで続く。他のパターンを用いることもでき、例えばパターンはＧｒａｙの正弦波の形態とすることもできることに留意すべきである。
【００２８】
所定の順序で投影されると、較正パターンは、一意の光エネルギーパターンを各位置１０４へ送出する。パターンは、［ｌｏｇ２（ｎ）］個のパターンしか必要とせずに、位置１０４の画素間の位置付けを区別する。ここで、ｎは投影画像の画素数である。
【００２９】
未処理の強度値を、パターンセットについて、各位置における光の有無［０，１］に対応する２進数シーケンス３１１に変換する。このビットシーケンスを次に、各位置の座標に対応する出力画像の画素の横座標および縦座標に適宜復号化する。
【００３０】
較正パターンの数は、位置の数およびその座標とは無関係である。表示面は、特にそれが任意の複雑な多様体である場合、任意数の検知位置を含むことができる。検知位置は表面に固定されているため、計算は大幅に簡単になる。実際に、全較正を一瞬で行うことができる。
【００３１】
計算の単純さおよび速度により、動的な較正が可能になる。言い換えれば、表示面が形状または位置を変化させている間に、画像または映像を表示面に投影しながら動的に較正を行うことができる。実際に、表面の形状は、検知データ３１１に動的に適合させることができる。また、赤外線センサ、または高速の瞬間的な潜像を用いることによって較正パターンを不可視にすることもできることに留意すべきである。したがって、較正パターンは、基礎となる表示プログラムの表示を妨害しない。
【００３２】
別法として、較正パターンは画像ペアとすることができ、ステガノグラフィのように、１つの画像の直後にその相補的な否定（negation）または反転を投影して、パターンを実質的に不可視にする。これには、光強度測定を微分して（be differential）周囲の背景光の貢献（影響）を少なくするという利点もある。
【００３３】
ワーピング関数
同一平面上の４つの位置について、較正モジュール１０５は、次のワーピング関数を求める。
ｐｓ＝ｗ＊ｐｏ
ここで、ｗはワープ行列（warp matrix）であり、ｐｓは検知位置の座標であり、ｐｏは、表示面の各位置と位置合わせされる出力画像中の対応画素の座標である。
【００３４】
検知したｐｓの値および既知のｐｏの値を用いて、ワープ行列ｗの対応関係を解くことができる。これは、ホモグラフィとして知られる、或る任意の四辺形領域を別の任意の四辺形にワーピングするための従来の技法である。正式には、ホモグラフィは、１つの座標枠における３Ｄ平面の画像を第２の座標枠の画像にマッピングする３×３の８自由度の投影変換Ｈである。ホモグラフィの計算方法はよく知られている。Sukthankar等著「カメラホモグラフィツリーを用いた大型マルチプロジェクタディスプレイのスケーラブルな位置合わせ（Scalable Alignment of Large-Format Multi-Projector Displays Using Camera Homography Trees）」（Proceedings of Visualization, 2002）を参照のこと。
【００３５】
通常、画素ｐｏは出力画像の角にある。４つより多くの位置を用いる場合、平面の最良一致プロセスを用いることができる。４つより多くの位置を用いることによって、較正の質が向上する。基本的に、本発明は、強度の測定値を用いて、それらの位置を出力画像中の対応画素に相関させる。
【００３６】
ワープ行列ｗは映像処理モジュール１０６に伝えられる。ワーピング関数は、入力画像を歪めて、位置、スケール、回転、およびキーストン歪みを補正し、結果として得られる出力画像を、歪みがなく画像領域と位置合わせされるように見せる。
【００３７】
非平坦な表面
任意の多様体は、面法線がプロジェクタの光軸に対して鈍角である位置を含む可能性がある。これらの位置に対応するセンサは、プロジェクタからの直接照明を受けず、較正プロセス中に不可視となる可能性がある。したがって、検知位置は、プロジェクタによって直接照明されるように選択されなければならない。
【００３８】
図３に示すような、任意の多様体に画像を投影するための一般化された技法は、Raskar他によって、２００１年８月１５日付で出願された「投影画像を用いて実オブジェクトを動画化するシステムおよび方法（System and Method for Animating Real Objects with Projected images）」と題する米国特許出願第０９／９３０，３２２号（参照により本明細書中に援用される）に記載されている。この技法は、表示面の幾何学的形状が分かっていること、および最低６つの較正点を用いることを要求する。投影画像は、用途に応じて、表示面の特徴を強調する（enhance）、拡大する、または隠す（disguise）ために用いることができる。出力画像を歪ませる代わりに、較正データを用いて、プロジェクタを新たな位置に機械的に移動させて歪みを補正することもできる。また、表示面自体を移動または変形させて歪みを補正することも可能である。上記を様々に組み合わせて、例えば出力をワーピングするとともにプロジェクタを移動させる、または出力をワーピングするとともに表示面を移動させることも可能である。これらは全て、システムを較正した状態に保ちながら動的に行うことができる。
【００３９】
本方法の主目的は、入力画像を歪めるかまたはワーピングするために用いることができるプロジェクタパラメータを求めて、ワーピングした出力画像を、表示面上で歪みがないように見せることである。しかし、較正パラメータは、他の目的、例えば表示面に対するプロジェクタの姿勢（内部および外部の幾何学的パラメータを含む）を求めるためにも用いることができる。姿勢は、他の用途、例えば画像を強調した環境における照明の計算に、または合成オブジェクトをシーンに挿入するために用いることができる。
【００４０】
パラメータはまた、プロジェクタと表示面の間の距離を求める、幾何学的形状が既知である表面に対するプロジェクタ光線の入射角を求める（例えば、３Ｄレンダリングプログラムにおいて照明の計算を行うか、または入力強度を変化させて、表示面における画像強度を均一に見せる）、表面領域をプロジェクタによって照明されるセグメントと照明されないセグメントとに分類する、プロジェクタの接線方向歪みを求める、および表示面の変形を求めるために用いることができる。
【００４１】
本発明はまた、複数のプロジェクタを同時に較正するために用いることができる。ここで、複数のプロジェクタは、重複する画像を表示面に投影する。これは、表示面が非常に大きい（例えばプラネタリウム）、あるいは表示面を多くの側面から見る場合に有用である。
【００４２】
本発明を好適な実施形態の例として説明してきたが、本発明の精神および範囲内で様々な他の適応および修正を行うことができることが理解される。したがって、添付の特許請求の範囲の目的は、本発明の真の精神および範囲に入る変形および修正をすべて網羅することである。
【図面の簡単な説明】
【００４３】
【図１】本発明による光センサを備える平坦な表示面に対してプロジェクタを較正するシステムの概略図である。
【図２】非平坦な表示面に対してプロジェクタを較正するシステムの概略図である。
【図３】光センサを備える表示面に対してプロジェクタを較正する方法のフロー図である。
【図４】本発明により用いられるＧｒａｙコード較正パターンを示す図である。
【図５】別個の光センサを有する表示面の側面図である。
【図１】
2007536766.tif 000003
【図２】
2007536766.tif 000004
【図３】
2007536766.tif 000005
【図４】
2007536766.tif 000006
【図５】
2007536766.tif 000007
【国際調査報告】
2007536766.tif 000008
2007536766.tif 000009
2007536766.tif 000010