(57)【要約】本発明は、ノイズの多い構造において欠陥を潜在的に有する媒体を解析するための超音波医用撮像装置に関連する。Ｍ個の別個の連続励起により所与の深度に集束される超音波信号が放射され、画像形成モジュールが、媒体からの応答の受信後にその深度の画像を得ることを可能にする。本発明では、装置は、面積Ｎ＊Ｍの矩形応答行列を構成し、応答行列を特異値に分解し、更に、媒体における欠陥に対応する特異ゾーンを見つけるよう特異値に対応する特異ベクトルを使用するために、応答を利用するモジュールを有する。面積Ｎ＊Ｍの矩形応答行列の係数Ｋｎｍは、励起ｍ後の振動子ｎによって受信される媒体の応答を表す。本発明は、特に医用分野の超音波装置に適用される。
2006512129.tif 000002
【特許請求の範囲】
【請求項１】
ノイズの多い構造において欠陥を潜在的に有する有機媒体を解析する方法であって、
前記媒体は、Ｎ個の振動子からなるセットによって放射され、且つ、Ｍ個の別個の連続励起で所与の深度に集束される超音波信号によって励起され、それにより、前記媒体からの応答の受信後に、前記深度の画像を得る方法であって、
面積Ｎ＊Ｍの矩形応答行列を構成する段階と、
前記応答行列を特異値に分解する段階と、
前記媒体における欠陥に対応する特異ゾーンを見つけるよう前記特異値に対応する特異ベクトルを使用する段階と、
を有し、
前記面積Ｎ＊Ｍの矩形応答行列の係数Ｋｎｍは、励起ｍ後に振動子ｎによって受信される前記媒体の応答を表す、解析方法。
【請求項２】
応答行列Ｋｎｍは、複数の周波数に対して得られる請求項１記載の解析方法。
【請求項３】
Ｍ個の連続励起は、前記媒体の複数の深度に対して行われる請求項１又は２記載の解析方法。
【請求項４】
ノイズの多い構造において欠陥を潜在的に有する媒体を解析するための超音波医用撮像装置であって、
前記装置は、Ｍ個の別個の連続励起により所与の深度に集束される超音波信号を放射する振動子のセットと、
前記媒体からの応答の受信後に、前記深度の画像を得るための画像形成モジュールと、
面積Ｎ＊Ｍの矩形応答行列を構成し、前記応答行列を特異値に分解し、更に、前記媒体における欠陥に対応する特異ゾーンを見つけるよう前記特異値に対応する特異ベクトルを使用するために、前記応答を利用するモジュールと、
を有し、
前記面積Ｎ＊Ｍの矩形応答行列の係数Ｋｎｍは、励起ｍ後に振動子ｎによって受信される前記媒体の応答を表す、装置。
【請求項５】
応答行列Ｋｎｍは、複数の周波数に対して構成される請求項４記載の装置。
【請求項６】
Ｍ個の連続励起は、前記有機媒体の複数の深度に対して行われる請求項４又は５記載の装置。
【請求項７】
請求項４乃至６のうちいずれか一項記載の装置において使用されるプロセッサによって実行されるべきコンピュータプログラムであって、
請求項１乃至３のうちいずれか一項記載の有機媒体を解析する方法の段階を実行する命令のセットを有することを特徴とするコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ノイズの多い構造において欠陥を潜在的に有する有機媒体を解析する方法に係り、この媒体は、一組の振動子によって放射される超音波信号によって励起される。本発明は、特に、医用撮像及び超音波撮像機器内に導入可能な進歩した機能に係る。本発明は、特に、乳房の撮像及び乳房における微細石灰化の検出に係る。
【背景技術】
【０００２】
このような方法は、２００２年の８月のIEEE  Transactions  on  Ultrasonics,  Ferroelectrics,  and  Frequency  Control（第４９巻、第８号）においてE.  Kerbrat、C.  Prada、D.  Cassereau、及びM.  Finkによって公開された論文「Ultrasonic  Nondestructive  Testing  of  Scattering  Media  using  the  Decomposition  of  the  Time-Reversal  Operator」から公知である。公知の方法は、時間反転演算子の分解を用いて媒体を調べることを提案する。提案される方法では、振動子のアレイにおける第１のアレイが、短い励起によって励起され、媒体の応答の結果もたらされる信号は、そのアレイにおける全ての振動子によって受信される。この演算は、同じ励起を用いて各振動子に対して繰り返される。媒体のこれらの応答のフーリエ変換を生成することによって、四角の伝達行列Ｋが得られる。時間反転演算子は、Ｋ＊Ｋによって決定され、対角行列にされることが可能である。有意に非ゼロの固有値の数は、この方法によって検出される欠陥の数に等しい。次に、欠陥は、固有ベクトルの計算によって位置が突き止められる。
【０００３】
この文書の中で提案される方法は、方法に特有の多くの励起を必要とするという欠点を有する。これらの励起の特定の特性は、媒体内に存在する情報の一部のみしか考慮に入れることができないようにしてしまう。この方法は、従って、例えば、媒体の画像を得るために他の情報にアクセスを有することを可能にするために、媒体の他の高周波音波照射（insonification）と並列に、且つ、独立して使用されなければならない。更に、これらの特定の励起は、一般的な超音波撮像装置では行うことができず、従って、特定の装置の使用が必須となる。これは、結果が高速である、又は、更にはリアルタイムであることを必要とする適用のためにデータの容易且つ高速な獲得を必要とする医用撮像の分野においては問題である。更に、ここに開示される方法では、エネルギーは、媒体にほとんど伝達されず、これは、超音波の制限された伝搬をもたらす。この制限された伝搬は、正しい超音波画像を生成することを可能にしない。
【発明の開示】
【発明が解決しようとする課題】
【０００４】
本発明は、従来技術の方法の欠点を有さない、ノイズの多い構造において欠陥を潜在的に有する媒体を解析する方法であって、その媒体は超音波信号によって励起される方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の目的は、超音波信号は、Ｍ個の別個の連続励起によって所与の深度に集束され、それにより、媒体からの応答の受信後にその深度の画像を得るような導入段落に記載したような解析方法によって達成される。この方法は更に、面積Ｎ＊Ｍの矩形応答行列を構成する段階と、応答行列を特異値に分解する段階と、媒体における欠陥に対応する特異ゾーンを見つけるよう特異値に対応する特異ベクトルを使用する段階とを有し、面積Ｎ＊Ｍの矩形応答行列の係数Ｋｎｍは、励起ｍ後の振動子ｎによって受信される媒体の応答を表す。
【０００６】
本発明では、媒体は、例えば、一連の振動子又は振動子アレイの振動子の一連の幾何学的な一部分上に集中されて、超音波撮像において従来から使用される集束により励起される。
【０００７】
本発明は更に、媒体における欠陥を見つけるために振動子のそれぞれによって個々に受信した媒体のエコー応答を使用する。従って、本発明は、媒体の任意の特定の励起を必要とせず、また、従って、単一のデータ獲得を達成することを可能にする。更に、本発明は、反射特性の局所的な変化による特異ゾーンの発生源である反射体の検出及び位置の選定における認め得るほどの改善を得るために、小さな変更のみで超音波撮像装置に導入可能である。
【０００８】
第１の実施例では、応答の行列Ｋｎｍは、複数の周波数に対して得られる。実際は、様々な特異値は、全ての周波数に対して同じ強度で現れないことを以下に明記する。従って、幾つかの応答行列を構成することが有利であり得る。幾つかの応答行列、つまり、１つの応答行列を複数の周波数からの１つの周波数に対して構成することが有利である。
【０００９】
１つの有利な実施例では、Ｍ個の連続励起は、媒体の複数の深度に対して達成される。これは、一般的に、画像を得るために行われ、また、様々な深度における応答行列を構成するために有利に使用され、更に、媒体の拡大ゾーンでの特異ゾーンの検出を有するために行われる。
【００１０】
本発明の１つの好適な実施例では、特異ゾーンを見つけるために固有ベクトルを使用する段階は、媒体のバイナリ画像の形成をもたらす。即ち、値１は、コヒーレントな反射体の存在が検出されるゾーンに割当てられ、値０は、それ以外（例えば、ノイズの多いゾーン）に割当てられる。この実施例の目的は、本発明の方法によって与えられた情報の表示を可能にする、また、特に、例えば、媒体の高周波音波照射を適応させる、欠陥が検出されたゾーンを区別して処理する等のためにこの情報を利用することを可能にすることを目的とする。
【００１１】
本発明は、任意の超音波装置内に、ノイズの多いデータが結果としてもたらされる媒体におけるコヒーレントな反射体を検出する手段を導入することを可能にする。従って、本発明は更に、ノイズの多い構造において欠陥を潜在的に有する媒体を解析するための装置にも関連し、この装置は、Ｍ個の別個の連続励起により所与の深度に集束される超音波信号を放射する振動子のセットと、媒体からの応答の受信後に、深度の画像を得るための画像形成モジュールと、面積Ｎ＊Ｍの矩形応答行列を構成し、応答行列を特異値に分解し、更に、媒体における欠陥に対応する特異ゾーンを見つけるよう特異値に対応する特異ベクトルを使用するために、応答を利用するモジュールとを有し、面積Ｎ＊Ｍの矩形応答行列の係数Ｋｎｍは、励起ｍ後の振動子ｎによって受信される媒体の応答を表す。
【００１２】
従って、本発明は、医用撮像、特に、従来において得られる画像はノイズが多く、コヒーレントな反射変化は検出するのが困難であった超音波撮像の分野において有利に適用される。従って、本発明の装置は、従来通りの医用撮像ステーションである。
【００１３】
合成画像上では可視ではないコヒーレントな反射体を見つけることを可能にすることによって、本発明は更に、撮像結果の改良に寄与し、また、医用撮像の特定のケースにおけるより精密及び正確な診断を可能にし得る。この
本発明は更に、Ｍ個の別個の連続励起により所与の深度に集束される超音波信号を放射する振動子のセットと、媒体からの応答の受信後に、深度の画像を得るための画像形成モジュールとを有する撮像装置内に挿入可能な信号利用モジュールに関連し、この信号利用モジュールは、面積Ｎ＊Ｍの矩形応答行列を構成し、応答行列を特異値に分解し、更に、媒体における欠陥に対応する特異ゾーンを見つけるよう特異値に対応する特異ベクトルを使用することによって応答を利用する。面積Ｎ＊Ｍの矩形応答行列の係数Ｋｎｍは、励起ｍ後に振動子ｎによって受信される前記媒体の応答を表す。
【発明を実施するための最良の形態】
【００１４】
本発明は、図面に示す実施例を参照しながら更に説明するが、本発明は、これらの実施例に制限されるものではない。
【００１５】
以下に述べることは、参照符号に関する。同様のエンティティは、全図面において同一の文字で示す。幾つかの同様のエンティティが、単一の図面中に示される場合もある。その場合、類似のエンティティを区別するために、文字参照に、数字又は添え字が付けられる。この数字又は添え字は、便宜上省略する場合もある。このことは、説明及び請求項に適用する。
【００１６】
以下の説明は、当業者が、本発明を実施し且つ利用することができるよう提示する。好適な実施例の様々な代案が、当業者には明らかであり、また、本願に開示する本発明の一般的な原理は、他の実施例にも適用し得る。従って、本発明は、記載する実施例に制限されるものではなく、以下に説明する原理及び特徴に応じた最大限の範囲を有すると考えられる。
【００１７】
図１は、本発明の要件に応じて媒体ＭＩＤからくるエコー超音波信号の受信を説明する図である。媒体ＭＩＤは、集束された超音波ＦＯＣによって励起される。集束は、振動子アレイＡＲＲにおけるＰ個の振動子ＴＲ上で行われ、これらのＰ個の振動子の幾何学的中央に中心が合わされる。図１では、Ｐ＝４である。集束技術は更に、振動子アレイの任意の点に、超音波を集中することを可能にする。本発明では、媒体ＭＩＤにより返されたエコー信号は、次に、振動子ＴＲのアレイＡＲＲにおけるＮ個の振動子ＴＲのそれぞれにおいて記録される。医用撮像において従来から使用されている超音波獲得方法では、励起は、Ｐ個の振動子に、同じ集束であるが、図１に矢印ＳＣと示す走査方向において前の励起とオフセットにされて繰り返される。従って、従来の走査では、Ｍ個の獲得が行われる。獲得数Ｍは、様々であり得、そして、一般的に、振動子ＴＲのアレイＡＲＲにおける振動子ＴＲの数Ｎとは異なる。更に、従来の獲得は、一般的に、また、本発明の有利な実施例では、点Ｆ１、Ｆ２、及びＦ３と示す様々な集束深度に対して行われる。所与の深度における各獲得は、その深度における特定の情報を与える。しかし、これらの励起は、従来技術の提案される方法の実施には必須であるパルスに対する要素間応答（inter-element  response）にアクセスすることを可能にしない。
【００１８】
図２は、本発明を実施する超音波撮像装置の部分的な機能図である。この機能図は、より具体的には、振動子アレイにおける第１のＰ個の振動子に集束される広い周波数スペクトル波によって行われる第１の励起ｍ＝１を用いて本発明に従って行われた獲得を示す。例えば、スペクトルは、３乃至５ＭＨｚの周波数に中心が合わされ、全帯域幅の４０％の帯域幅を有する。エコー信号Ｓ［ｎ＝１，ｍ＝１］…Ｓ［ｎ＝Ｎ，ｍ＝１］は、励起ｍ＝１の後に、振動子ｎ（ｎ∈［１，Ｎ］）のそれぞれによって独立して受信され、そして、当業者には周知である手段に応じて超音波画像を生成するようビーム形成手段ＢＦに伝送される。この超音波画像は、図５を見て分かるように、一般的に、媒体ＭＩＤと、この媒体ＭＩＤの中にある物体ＯＢＪを観察することを可能にするが、一般的に、特に、スペックルノイズを示す、一様にノイズの多い画像における小規模の欠陥Ｘを区別することは可能にしない。この物体は、例えば、臓器であり、小さな欠陥は、病んでいる領域があることによる反射体であり得る。従って、乳房の場合、この小さな欠陥は、微細石灰化であり得る。本発明では、振動子アレイ全体における各振動子ｎによって独立して受信されたエコー信号Ｓ［１，ｍ＝１］…Ｓ［Ｎ，ｍ＝１］（Ｓｎ１とも示す）は、各振動子において受信された信号の時間部分を選択するセレクタＳＥＬにも伝送される。この時間部分は、一般的に、図１に示すような焦点Ｆ（１、２、又は３）の付近から来て受信される信号に対応する。この付近の大きさは、ユーザが有することを希望する、大きい深度における検出とこの検出の精度との間の妥協に依存する。次に、ｍ個の励起の走査が、医用超音波撮像において従来から使用される技術に従って、振動子アレイに沿って行われる。受信した信号Ｓ［１，ｍ］…Ｓ［Ｎ，ｍ］の時間部分は、Ｍ個の励起の媒体の走査の励起ｍのそれぞれに対し選択される。従って、ｍは、［１，Ｍ］に含まれる。これらの信号時間部分は、ｋｎｍ又はｋ［ｎ，ｍ］と示され、また、時間関数である。次に、これらの信号部分ｋ［ｎ，ｍ］は、信号の特定の利用のためのモジュールＰＥＭに伝送される。走査の各励起ｍに対し、モジュールＰＥＭは、特定の深度間隔に対応する媒体ＭＩＤからの振動子ｎによって受信された応答関数Ｓｎｍの部分ｋｎｍを、メモリ内に格納する。信号のフーリエ変換ｋｎｍ（ｔ）は、行列Ｋ＝（Ｋｎｍ（ω））１＜ｎ＜Ｎ；１＜ｍ＜Ｍを与え、これは、応答行列と呼ぶ。従って、媒体の集束された励起ｍの後に素子ｎによって受信された所与の励起周波数に対する媒体の応答をそれぞれ表す係数Ｋｎｍの行列が得られる。この行列は、矩形であり、また、一般的に、その要素（making）とは別個の要素によってスペクトルの各周波数に対し計算可能である。単一の周波数に対して１つの行列のみを獲得することが可能であるが、結果は、あまり正確ではないかもしれない。選択される周波数も、解像度及び媒体における減衰の制約に合うよう選択され得る。このような周波数は、実際には、本発明では、従来の媒体の超音波画像の獲得によって決定される。
【００１９】
モジュールＰＥＭは次に、応答関数Ｋの特異値への分解を計算する。事実上、この分解は、特に、特異系の解決（resolution）に対して使用され、また、実係数又は複素係数を有する面積ＮＭの矩形行列は、形式Ｋ＝ＵＤＶに分解され得、ここでは、Ｕは、面積ＮＮのユニタリ行列であり、Ｖは、面積ＭＭのユニタリ行列であり、Ｄは、面積ＮＭの対角行列である。面積ＮＭの行列Ｄの対角要素は、単に、行列Ｋ#Ｋの特異値の平方根であり、Ｋ#は、転置の共役である。
【００２０】
行列Ｋ#Ｋの特異ベクトルは、Ｕの列である。
【００２１】
第１の実施例では、モジュールＰＥＭは、複数の周波数に対して複数の応答行列を計算する。この場合、図３に示すような周波数ｆに応じた特異値ＶＰの振幅ＡＭＰを表すグラフが得られる。特異値ＶＰ１、ＶＰ２、及びＶＰ３は全て同じ周波数値に対して同じ相対強度で現れるわけではない。媒体における欠陥は、信号の反射における局所的な変化によってマーク付けされるので、反射体表現（reflector  term）で考えられ且つ決定されることが可能である。時間反転演算子の対角化の簡単な方法に対し、反射体の存在と非ゼロ特異値の存在との対応が、C.  Prada、M.  Finkによる「Eigenmodes  of  the  Time-Reversal  Operator:  A  Solution  to  Selective  Focusing  in  Multiple-Target  Media」（Wave  Motion  20/1994,  pp.  151-163）において調べられた。本発明は、この対応を得ることも可能であることが観察された。従って、本発明では、調べられた周波数領域において構成された矩形行列の非ゼロ特異値は、反射体の存在を明らかにする。従って、対応は、特異値＝反射体、また、最大特異値は、最大反射体に対応する。
【００２２】
第１の実施例では、複数の行列は、様々な周波数に対し構成され、行列Ｕの列である固有ベクトルの逆フーリエ変換が計算可能である。このことは、特異周波数ベクトルに対応する特異時間ベクトルを得ることを可能にする。これは、次には、圧力場を決定するよう単に、特異時間ベクトルを媒体内に伝搬することを可能にする。この圧力場の最高点は、媒体における欠陥に対応する。位相及び振幅における特異ベクトルの戻り伝搬は、例えば、受信集束技術を使用し、また、一般的に、媒体における音場をシミュレートするデジタル手段によって行われる。実際には、ソフトウェアが、この概念的な時間波動伝搬の機能を行う。例えば、図４では、この特異時間ベクトルの第１の伝搬技術は、高い圧力ゾーンＰＦＩの特定の方法でのカラーリングをもたらす。カラーレベルを有する画像も得ることが可能である。
【００２３】
ソフトウェアは更に、媒体の各深度に対する、従って、別個にされる伝搬行列を再構成するために提案されることも可能である。センサの平面から所与の深度の平面までの通過のための行列が得られ、従って、面積Ｙ上の反射体の位置を突き止めることを可能にする。有利な実施例を用いたこの様々な深度に対する再構成は、様々な深度において且つ別個にされるような方法で集束される波による媒体の高周波音波照射の結果であり得る。例えば、図１では、３つの深度Ｆ１、Ｆ２、Ｆ３の集束が達成され、そして、これらの３つの深度に対して伝搬行列が構成される。上に提案した特異ゾーンの位置を突き止める２つの技術が、１つのゾーンを分離し、また、実際には、図５に示すように物体ＯＢＪを含む媒体ＭＩＤの従来の超音波画像上に圧力場ＰＦＩを表示することを可能にする。従って、１つ以上の欠陥を見つけるために調べられたゾーンの超音波画像上にマーク付けが行われる。例えば、本発明によって、微細石灰化は、乳房内において検出可能である一方で、従来の超音波画像上のノイズ（「スペックル」）からは浮かび上がってこない。
【００２４】
本発明は、超音波獲得を用いる全てのタイプの医用撮像に導入可能である。本発明では、超音波撮像において使用される様々なタイプの集束による励起後の有機媒体によって与えられる応答を使用することが可能である。線密度（つまり、２つの連続する励起間の幾何学的間隔）も、本発明に関係なく適応され得、それにより、本発明によって与えられる結果を横方向により正確にし得る。従って、数ＭをＮより大きくすることが可能である。この場合、分解されるべき系は、退縮（degenerate）される。
【００２５】
本発明は、適応可能な高周波音波照射ビームを形成するよう用いられることが可能である。即ち、特異ベクトルは、扱いにくい領域に強い圧力場を送ることを可能にし、それにより、この領域に関するより正確な情報を有することを可能にする。
【００２６】
図６は、本発明の方法が内部で実施される装置を示す。本発明は、移動不可か、又は、モジュラー形式で実施可能であり、反射体検出モジュールが、従来の超音波装置に追加される。この検出モジュールは、超音波装置の振動子と関係なくエコー信号を受信し、また、例えば、上述したようなセレクタＳＥＬと利用モジュールＰＥＭを有する。図６では、本発明が不変である装置を示す。
【００２７】
この装置は、受信素子ＴＲを有するプローブＰＲＯＢを有し、このプローブは、従来の手段によって、データ処理装置ＬＡＢに接続される。公知の超音波撮像の技術に応じて戻りビーム及び画像を形成するモジュールＢＦに追加して、任意のデータ処理装置ＬＡＢは、上述したようなセレクタＳＥＬ及び動作モジュールＰＥＭを有する。装置ＬＡＢは、従来の表示機能によって、超音波装置によって従来通りに得られる画像に加えて、モジュールＰＥＭによって得られる情報から構成される画像を表示するディスプレイモジュールＤＩＳに接続される。組み合わせモジュールＣＭＢは、例えば、ビーム形成手段ＢＦによって得られたデータ、及び、動作モジュールＰＥＭによって得られたデータを組み合わせる。次に、組み合わせモジュールＣＭＢは、ディスプレイモジュールＤＩＳに接続される。従って、図５に示されるような、また、特異ゾーンの位置決めをする画像が、本発明によって得られる。特異ゾーンのグラフィカル表現の任意の手段（ゾーンを囲むバイナリ画像等）が、本発明の適用に区別なく使用することが可能である。ユーザインタフェースＵＩＦが、装置ＬＡＢを制御し且つそれをパラメータ化するために、装置ＬＡＢに有利に接続される。例えば、検出閾値、並びに、本発明の目的である反射体の位置の選定／検出の精度を決定可能な深度増加値がユーザによって変更可能である。
【００２８】
本発明は、ノイズの多い信号（「スペックル」）が得られる均一な媒体における任意の反射する欠陥の正確な場所を得ることを可能にし、それらの信号内では、公知の手段を用いてそのような欠陥を検出することは一般的に困難である。その適用の１つにおいて、本発明は、様々な方向において媒体を高周波音波照射し、より完全且つノイズの少ない画像が得られるよう結果を組み合わせることよりなる複合撮像に有利に関連する。
【００２９】
本発明の方法の段階において示された機能を満たす前述したモジュールは、従来の超音波装置における追加の適用として組み込まれるか、又は、本発明の機能を満たすよう従来の超音波装置に接続されることを意図する独立装置において使用されることが可能である。本発明の方法の段階において示された機能を、当業者がアクセス可能であるソフトウェア及び／又はハードウェア手段によって実施する方法は多く存在する。これは、図面が概略的である理由である。従って、図面は、様々なユニットによって行われる様々な機能を示すが、これは、幾つかの機能を実行可能にする単一のソフトウェア及び／又はハードウェア手段を排除するものではない。更に、これは、１つの機能を実行可能にするソフトウェア及び／又はハードウェア手段の組み合わせを排除するものではない。本発明は、提示する実施例に従って説明したが、当業者は、提示する実施例に対する変形が存在し、これらの変形は、本発明の精神及び範囲内であることを即座に認識するであろう。従って、多くの修正が、請求項によって決定される精神及び範囲から排除されることなく、当業者によって達成可能である。
【図面の簡単な説明】
【００３０】
【図１】本発明による媒体からくるエコー信号の受信と、矩形応答行列の形成を説明する図である。
【図２】本発明を実施する超音波撮像装置を示す部分機能図である。
【図３】本発明の第１の実施例に応じて得られる特異値のグラフである。
【図４】本発明の特異ゾーンの場所を示す図である。
【図５】超音波画像への本発明の適用を示す図である。
【図６】本発明の媒体を解析する装置を示す概略機能図である。
【図１】
2006512129.tif 000003
【図２】
2006512129.tif 000004
【図３】
2006512129.tif 000005
【図４】
2006512129.tif 000006
【図５】
2006512129.tif 000007
【図６】
2006512129.tif 000008
【国際調査報告】
2006512129.tif 000009
2006512129.tif 000010