(57)【要約】本発明はマイクロバブルを含む目標物から反射される超音波エコー信号を検出および撮像するための方法に関し、前記マイクロバブルは共振周波数によって特徴付けられ、前記方法は、前記目標物に少なくとも第１および第２の引き続く励起信号を照射するステップを有し、前記第１の励起信号は時間とともに増加する周波数掃引であり、前記第２の励起信号は時間とともに減少する周波数掃引であり、前記目標物からの前記第１および第２の励起信号のエコー信号を検出するステップと、前記エコー信号を組み合わせるステップと、をさらに有する。前記第１および第２の励起信号の最大周波数は前記共振周波数より低い。本発明は撮像装置にも関する。
2008543475.tif 000002
【特許請求の範囲】
【請求項１】
マイクロバブルを含む目標物から反射される超音波エコー信号を検出および撮像するための方法であって、
前記マイクロバブルは共振周波数によって特徴付けられ、
前記方法は、
前記目標物に少なくとも第１および第２の引き続く励起信号を照射するステップを有し、
前記第１の励起信号は時間とともに増加する周波数掃引であり、
前記第２の励起信号は時間とともに減少する周波数掃引であり、
前記目標物からの前記第１および第２の励起信号のエコー信号を検出するステップと、
前記エコー信号を組み合わせるステップと、
をさらに有し、
前記第１および第２の励起信号の最大周波数は前記共振周波数より低い方法。
【請求項２】
前記第１および第２の励起信号の最大周波数は同一である請求項１に記載の方法。
【請求項３】
前記第１および第２の励起信号は同一の周波数掃引を有する請求項１または２に記載の方法。
【請求項４】
前記第２の励起信号は、前記第１の励起信号の時間を反転した複製である請求項１から３のいずれか１項に記載の方法。
【請求項５】
前記第１および第２の励起信号の最大周波数は、前記マイクロバブルの共振周波数の９０％より高い請求項１から４のいずれか１項に記載の方法。
【請求項６】
前記第１および第２の励起信号の最大周波数は、前記マイクロバブルの共振周波数の９８％より高い請求項５に記載の方法。
【請求項７】
前記周波数掃引は時間の線形関数である請求項１から６のいずれか１項に記載の方法。
【請求項８】
前記周波数掃引は時間の非線形関数である請求項１から７のいずれか１項に記載の方法。
【請求項９】
マイクロバブルを含む目標物の超音波撮像のための装置であって、
前記マイクロバブルは共振周波数によって特徴付けられ、
前記装置は、
前記目標物に少なくとも第１および第２の引き続く励起信号を照射する照射手段（１，２）を備え、
前記第１の励起信号は時間とともに増加する周波数掃引であり、
前記第２の励起信号は時間とともに減少する周波数掃引であり、
前記目標物からの前記第１および第２の励起信号のエコー信号を検出する検出手段（６）と、
前記エコー信号を組み合わせる組み合わせ手段（７）と、
をさらに備え、
前記第１および第２の励起信号の最大周波数は前記共振周波数より低い装置。
【請求項１０】
前記第１および第２の励起信号の最大周波数は同一である請求項９に記載の装置。
【請求項１１】
前記第１および第２の励起信号は同一の周波数掃引を有する請求項９または１０に記載の装置。
【請求項１２】
前記第２の励起信号は、前記第１の励起信号の時間を反転した複製である請求項９から１１のいずれか１項に記載の装置。
【請求項１３】
前記第１および第２の励起信号の最大周波数は、前記マイクロバブルの共振周波数の９０％より高い請求項９から１２のいずれか１項に記載の装置。
【請求項１４】
前記第１および第２の励起信号の最大周波数は、前記マイクロバブルの共振周波数の９８％より高い請求項１３に記載の装置。
【請求項１５】
前記周波数掃引は時間の線形関数である請求項９から１４のいずれか１項に記載の装置。
【請求項１６】
前記周波数掃引は時間の非線形関数である請求項９から１４のいずれか１項に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概して、超音波コントラスト撮像技術を使用した身体の超音波診断に関し、特に、超音波コントラスト因子（ultrasound contrast agent（ＵＣＡ））を使用した方法および装置に関する。
【背景技術】
【０００２】
過去１０年以上、画像品質を向上させるために小さい気泡の形態でコントラスト因子が導入された。気泡は、他の器官から反射されるエコーから所望の器官を明瞭に識別するため、所望の器官から後方散乱されるエコーを増加させるために関心領域に注入される。小さいマイクロバブルの形態で溶けて気体は安定した。
【０００３】
超音波コントラスト撮像の有用性を拡張するために、研究は、効果的なＵＣＡおよび新たなコントラスト固有の撮像技術を開発することに活発に集中してきた。第２高調波に基づく技術は、心房のような多くの構造体の中でＵＣＡの検出を高める。それらは、気体のマイクロバブルおよび組織の間の超音波照射への応答の相違を利用する。超音波と相互作用するときコントラストバブルは非線形のまたは調和する振る舞いを示すのに対し、軟部組織は線形の反射体であることが知られている。
【０００４】
コントラストバブルの注入なしでも、第２高調波（２Ｈ）周波数において組織を撮像する技術が開発された。その特定の周波数は、画像品質を著しく向上させることが分かる。このモードは、造影剤を使用しない組織ハーモニック撮像（native tissue harmonic imaging）と呼ばれ、各種の臨床の応用において効果を立証してきた。最近の研究は、２Ｈ周波数のみでコントラストバブルおよび組織を撮像することは、いくつかの不利益と関係することを示した。そして、より高いハーモニック周波数は、コントラスト気泡と周囲の散乱媒体との間で検出または識別する能力を増加させる方法として示唆された。そのようなアプローチは、見込みのある生体外および生体内の結果を示した。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、かなりの研究の試みは、未だ、心筋灌流の評価においてＵＣＡの使用が約束される必要がある。パルス反転および電力変調のような新たなコントラスト撮像法は、左室造影（left ventricular opacification（ＬＶＯ））のために、現在、効果的に使用されているが、“撮像困難な”患者の超音波心臓検査法において心筋灌流を検出するために限られた成功を示した。
【０００６】
従って、新たな撮像法は、気体のマイクロバブルの固有の音響特性に基づいて研究されなければならない。最終的な灌流技術は、泡のエコーを増加させ、強いコントラスト対組織比（contrast-to-tissue ratio）を可能とする一方で、強い（線形または非線形の）組織のエコーの抑止を確実にすることを可能とするべきである。灌流していない組織とコントラストの灌流している組織との間の識別は、一般に、コントラスト対組織比と呼ばれる。
【課題を解決するための手段】
【０００７】
従って、本発明の目的は、既知の方法より向上したコントラスト因子検出のための方法を提供することである。本発明のもう１つの目的は、コントラスト因子検出のための向上した装置を提供することである。
【０００８】
従って、本発明は、請求項１による、マイクロバブルを含む目標物から反射される超音波エコー信号を検出および撮像するための方法を提供する。本発明は、さらに、請求項９による、超音波撮像のための装置を提供する。
【０００９】
本発明は、気体のマイクロバブルの共振の振る舞いを利用する。それらの共振周波数に近い周波数で励起された振動するマイクロバブルからのエコーが、組織のような振動しない構造体から放射されるエコーから検出およびえり分け（除去）される。
【００１０】
本発明の他の特徴および効果は、添付図面に関連して検討されるとき以下の説明においてさらに明らかになる。
【発明を実施するための最良の形態】
【００１１】
目標物に存在するガス状または気体のマイクロバブルは、注入されまたは吸収したＵＣＡ、または、ガス状の塞栓形成（embolism）、塞栓（emboli）の場合であり得る。
【００１２】
本発明による方法は、そのような気体のマイクロバブルの共振の振る舞いに基づく。気体のマイクロバブルに超音波が衝突すると、それは振動し、従って、超音波エネルギーを放射する。超音波の位相に応じて、泡は収縮または膨張する。振動する泡は、調和振動子に匹敵し、従って、自然の共振周波数によって特徴付けられ得る。泡の共振周波数はMinnaertによって定義され、数式：
【００１３】
【数１】
2008543475.tif 000003
【００１４】
によって与えられる。ここで、
Ｒは静止半径、
γはポリトロープ成分、
Ｐ０は周囲の圧力、
ρは媒体の密度である。
【００１５】
この数式は、泡の共振周波数がその大きさに反比例することを表わす。図１は、気体のマイクロバブルの共振周波数をその大きさの関数として表わす。直径６μｍの泡について、その共振周波数はおよそ１ＭＨｚである。
【００１６】
マイクロバブルがその共振周波数で励起されると、その振動は最大であり、従って、散乱エコー信号は強められる。図２は、直径３μｍの泡についてのこの実例におけるマイクロバブルの散乱の強さを示す散乱断面積を表わす。曲線は、その共振周波数に等しい周波数で生じる最大値を表わす。強いことに加えて、散乱信号の長さは、泡の共振の振る舞いのために増加される。図３、図４、図５に例が与えられる。図３は、泡の共振周波数より下の周波数を有する音響信号で励起されるとき、上側に励起信号、下側に直径３μｍの泡からの散乱信号を表わす。泡の応答は、送信信号とほぼ同一の長さを表わす。図４は送信信号および共振する泡の応答を表わす。第１に、散乱信号の高い振幅、第２に、マイクロバブルの共振の振る舞いのために著しく増加された信号の長さを理解することができる。従って、励起信号が通過した後、マイクロバブルは自由に振動を続ける。その共振周波数より上で励起された泡について、図５に結果が表され、ずっと短い信号応答に加えて、より低い散乱振幅が分かる。
【００１７】
従って、その共振周波数において、または、その共振周波数の近くにおいて励起されたマイクロバブルについて、それは超音波エネルギーをたいへん効果的に散乱し、ずっと長いパルス応答を表わす。
【００１８】
組織のような振動しない構造体について、振動しない散乱体から共振の振る舞いが観察されないので、音響応答は励起信号と同一の長さである。
【００１９】
ＵＣＡとして注入されるマイクロバブルは広範囲の大きさを示し、典型的に１〜５μｍで、１０μｍまで、または１μｍより小さい。これは、経験的に、ＵＣＡの泡についてそのような共振周波数を定義することを困難にし得る。それにもかかわらず、実験結果は、全てのマイクロバブルが同一の大きさであるかのように、全体のマイクロバブルは支配的な共振周波数とともに振る舞うことを示す。本発明による方法の続く説明において、ＵＣＡの共振周波数は、この支配的な共振周波数として定義される。そのような共振周波数は、典型的に、１０ＭＨｚ程度への１ＭＨｚの範囲内である。血管の超音波周波数は１〜５ＭＨｚの範囲である。血管内の超音波周波数は１０〜４０ＭＨｚの範囲であり、これに関するマイクロバブルの大きさは１μｍより小さい。気体のマイクロバブルのこれらの共振特性を利用するために、本発明による方法では、続発する２つの信号、両方がいくつかの周波数成分を含む、が送信される。従って、マイクロバブルを含む目標物から反射される超音波エコー信号を検出および撮像するために、目標物は少なくとも第１および第２の引き続く励起信号が照射され、第１の励起信号は時間とともに増加する周波数掃引または増加掃引周波数信号であり、一方、第２の励起信号は時間とともに減少する周波数掃引または減少掃引周波数信号である。
【００２０】
各々の信号の周波数帯域は、調査される泡の共振周波数より下から、調査される泡の共振周波数に近い周波数までの周波数を含む周波数掃引を有する。従って、第１および第２の励起信号の最大周波数はマイクロバブルの共振周波数より低い。そして、第１および第２の励起信号からのエコー信号は、目標物から到来して検出され、組み合わされる。
【００２１】
送信信号において使用される周波数掃引に応じて、泡は異なって応答する。第１信号、すなわち、図６に表わされているような時間とともに増加する周波数の周波数掃引を使用するとき、マイクロバブルは、まず、その共振周波数より下の周波数によって、そして、送信パルスの終了において高周波数成分によって、応答させられる。送信信号の高周波数成分はその共振周波数により近いので、泡は共振し、信号の終了において自由に振動する。これは引き伸ばされたパルス応答によって図８に表わされている。周波数が共振周波数に近づくに連れて、マイクロバブルはより多くのエネルギーを蓄積し、送信信号の影響は累積される。
【００２２】
第２信号、すなわち、図７に表わされているような減少掃引周波数信号の場合において、まず高周波数成分が送信され、これらの高周波数成分は、マイクロバブルの共振周波数に対応する。従って、（共振から離れた）低周波数成分がすぐに続き、開始された共振応答を打ち消すので、これらの高周波数成分による泡の共振の振る舞いが開始されるが、たいへん短い時間についてのみである。そして、自由な振動はすぐに打ち消され、泡は（共振から離れた）続く低周波数成分に応答し、従って、より短い全体のエコー応答を与える。信号の周波数が時間とともに減少するに連れて、累積した影響は増加掃引周波数信号の１つに匹敵し減少される。
【００２３】
図１０は、上側に、第１の信号、および、時間を反転した第２の複製を含む送信シーケンスを表わす。下側において、送信シーケンスへの予測された泡の応答が与えられている。
【００２４】
選択される周波数掃引はマイクロバブルおよびそれらの共振周波数に固有であるので、組織のような振動しないかつ共振しない散乱体は、マイクロバブルとはたいへん異なって応答し、増加掃引および減少掃引の励起の両方にたいへん類似した応答を示す。実際、共振周波数から離れた周波数に従うとき、他の散乱体（組織）およびマイクロバブルの両方は、（図３におけるような）短い振動する応答を表わし、それらの応答は、周波数が増加しようと減少しようと、励起信号とほぼ同一の長さを表わす。マイクロバブルの共振周波数により近づくと、周波数が増加しようと減少しようと、組織は類似の応答を示すが、マイクロバブルは上述したのと異なる。
【００２５】
マイクロバブルと他の散乱体との間の振る舞いにおける相違は、正しい組み合わせによって目標物を撮像することを可能とする。第１および第２の送信信号からのエコー信号が検出された後、これらのエコー信号および目標物を撮像するために組み合わせが実行される。実行される組み合わせは、第１および第２の信号のそれぞれの形状を考慮しなければならない。
【００２６】
本発明による方法の好ましい実施形態において、第１（増加掃引）および第２（減少掃引）の励起信号の最大周波数は同一である。さらなる実施形態において、両方の送信信号の周波数成分は同一であり、すなわち、第１および第２の励起信号は同一の周波数掃引を有する。図６および図７の例において、第１および第２の信号は互いの時間を反転した複製である。両方の送信信号は、同一の全体の周波数帯域および同一の音圧を有する。図６は増加する周波数掃引を有するチャープ信号であり、図７は正確に同一の信号であるが、減少する周波数掃引を有する。
【００２７】
両方の送信信号への泡の応答を計算するために、気泡力学を説明する理論モデルが使用された。このモデルは有名なRayleigh-Plessetの方程式に基づく。シミュレーションにおいて、直径３μｍの自由な泡が考察され、適用された音圧は１．８ＭＨｚにおいて５０ｋＰａであった。
【００２８】
図８および図９に結果が表わされている。図８は、図６に表わされている信号が照射されたとき、気体のマイクロバブルから散乱されたエコーを表わす。図９は、図７に表わされている信号が照射されたとき、泡の応答を表わす。送信信号において使用される周波数掃引に応じて、泡は異なって応答する。図６に表わされている信号を使用するとき、泡は共振し、増加掃引信号の終了において自由に振動するので、周波数掃引は時間とともに増加し、上述の累積した影響は、大きい振動を有するより長い応答（図８における引き伸ばされたパルス応答）をもたらす。減少掃引周波数信号を用いて、より低い周波数成分が続き、開始された共振応答を打ち消すので、これらの高周波数成分による泡の共振の振る舞いがまず開始されるが、たいへん短い時間についてのみである。図９は、そのようなより短い全体のエコー応答を表わす。図１０は、上側において、第１信号、および、第１信号の時間を反転した複製である第２信号を含む送信シーケンスを表わす。下側において、送信シーケンスへの予測された泡の応答が与えられている。
【００２９】
図１０におけるような送信信号を使用した、増加または減少する掃引励起への応答の間の相違を検出する最も簡単な方法の１つは、両方の応答の減算である。エコーの減算の結果は図１１に表わされている。マイクロバブルの増加および減少掃引信号への異なる応答のために、この場合において強い信号が得られる。同じ減算は、組織のような他の散乱体についてごくわずかの信号を示す。
【００３０】
英国特許第２４０４０２４号明細書における例について説明されているような、または、任意の他の適切な信号処理技術のような、よく知られた相関法のような他の検出法を使用することが可能である。
【００３１】
マイクロバブルの共振の振る舞いを最大化するために、第１および第２の信号の最大周波数は、マイクロバブルの共振周波数の９０％より高くなければならない。本発明による方法の向上した実施形態において、最大周波数はマイクロバブルの共振周波数の９８％より高い。周波数掃引は、選択された、時間に固有の関数とすることができ、例えば、周波数掃引は、時間の線形または非線形関数すなわち、周波数掃引は、線形または非線形変調である。図６および図７の例において、増加または減少掃引は、周波数掃引の上限および下限の平均値である基本周波数に中央揃えされた、時間の線形関数として選択されている。図１２において、線形変調および（非線形変調の例として）二次変調が表わされている。
【００３２】
増加および減少周波数信号期間は、典型的に、数μ秒、より詳しくは、図６および図７の例におけるように、４〜６μ秒の範囲である。適切ならば、１０μ秒までのより長いチャープ信号を使用することが可能である。引き続く第１および第２の信号の間の期間は、増加掃引信号および減少掃引信号からのエコー信号の間の任意の干渉を防止するために十分でなければならない。送信信号の振幅は、上述した例におけるように、１０ｋＰａ程度とすることが可能である。
【００３３】
泡から受信されるエコー信号は、（共振の振る舞いのために）長過ぎることがあり得、撮像のためにそのまま使用できないことがあり得る。しかしながら、撮像の目的により適合する信号を生成するために圧縮技術を使用することができる。そのような技術は、Misaridis TX、Gammelmark K、Jorgensen CH、Lindberg N、Thomsen AH、Pedersen MH、Jensen JAによる“Potential of coded excitation in medical ultrasound imaging”、Ultrasonics.、２０００年３月、38(1-8):183-9、および、Borsboom JM、Chin CT、de Jong Nによる“Nonlinear coded excitation method for ultrasound contrast imaging”、Ultrasound in Medical and Biology、２００３年２月、29(2):277-84において知られている。
【００３４】
図１３において、本発明による超音波撮像のための装置が表わされている。パルス生成器または送信源１は、上述した増加および減少周波数掃引を生成し、調査を受ける器官系においてそれらを目標物１１に向けるように構成される。パルス生成器１は、マイクロバブル１２を含む目標物１１に向けられる超音波信号を生成するためのアレイトランスデューサ２に接続される。送信源１およびアレイトランスデューサ２は照射手段である。目標物１１から受信される超音波エコー信号はアレイトランスデューサ２によって検知され、受信手段として動作する、それに接続された受信器６において受信される。受信されたエコー信号は、後処理回路７によって分析され、組み合わされ、処理される。
【００３５】
図６、図７、図１０、図１３に表わされている異なる増加および減少掃引信号は、ガウス包絡線によって特徴付けられる。これは、アレイトランスデューサ２のような典型的なアレイトランスデューサの周波数応答のためである。他の包絡線の形状も同様に使用することができる。
【００３６】
この方法の利点は、まず、コントラストマイクロバブルへのより良い感度である。この方法は、泡のエコーを増幅し、それによってコントラスト対組織比を増加させる一方で、組織から放射されるエコーの打ち消しを提供する。さらに、この方法は、十分な線形モードで動作することができ、従って、十分なトランスデューサ周波数帯域を利用することができる。超音波コントラスト撮像のためのそのような方法は、十分なトランスデューサ周波数帯域にわたって動作することができる。また、この方法は、第２高調波または優調和（superharmonic）のような非線形領域で動作することが可能であるが（国際公開第０３／０５９１６９号パンフレットを参照）、専用トランスデューサ以外のトランスデューサ帯域を乱切する。
【００３７】
また、この方法は、（異なる極性の第２パルスが続くパルス信号を送信し、第１および第２パルスへのハーモニック応答を検出し、組み合わせることにおいて、米国特許第５，７０６，８１９号明細書において定義されている）パルス反転、および／または、電力変調のような現在のコントラスト撮像法と組み合わせることが可能である。
【００３８】
動く組織において本発明による方法を使用するとき問題が生じ得る。実際、マイクロバブルを含む組織は、第１および第２の照射信号の間に動き得る。目標物を撮像するためにエコー信号を組み合わせるとき、結果として偽信号が現れ得る。そのような偽信号は、上述したように、いくつかの増加および減少掃引周波数信号を送信することによって効果的に減少させることができる。結果のエコー信号の特定の組み合わせによって、組織の動きを訂正することができ、従って、撮像の解像度を向上させる。そのような方法は、Peter N. BURNS、David HOPE SIMPSON、Michalakis A. AVERKIOUによる“NONLINEAR IMAGING” Ultrasound in Medical and Biology、第２６巻、別冊１、pp. S19-S22、２０００年において知られている。
【００３９】
また、気体の塞栓はガス状のマイクロバブルから成るので、本発明による方法は塞栓の検出に適用可能である。
【００４０】
理論的な発見を検証するために、光学的および音響学的な実験が行われた。１５ＭＨｚで動作するBrandarisカメラを使用して光学的な測定が行われた。図１４は、３．２μｍの半径で特徴付けられる実験用のコントラストマイクロバブルを使用して得られた結果を表わす。マイクロバブルは、１．７ＭＨｚの中心周波数の増加掃引周波数チャープが照射された。そして、同じマイクロバブルが減少掃引周波数を含む時間が反転されたチャープを用いて応答させられた。送信音圧は２００ｋＰａであった。結果は、図１４において時間の半径時間曲線関数として表わされている。増加掃引周波数チャープ（上側）は、広帯域周波数応答を示すマイクロバブルの高度な減衰振動を表わし、一方、反転したチャープは、より長い半径時間曲線によって特徴付けられる振動する泡の応答を生じさせる。
【００４１】
図１５は、組織の疑似模型に第１チャープ、および、その反転した複製が照射されたときの音響学的な測定から得られた結果を表わす。前の実験の送信パラメータが繰り返された。図１５の上側は、第１チャープへの組織の応答を表わし、下側は、反転したチャープが照射されたときの組織のエコーを表わす。この結果は、振動する構造体（気体物質）のみが送信されたチャープの周波数掃引に応じて異なって応答することを示す、両方の励起への組織の疑似物質の類似の応答を表わす。
【図面の簡単な説明】
【００４２】
【図１】気体のマイクロバブルの共振周波数をその大きさの関数として表わす。
【図２】直径３μｍの泡についてのこの実例における散乱断面積、すなわち、マイクロバブルの散乱の強さを表わす。
【図３】泡の共振周波数より下の周波数を有する音響信号で励起されるときの励起信号（上側）とともに３μｍの泡からの散乱信号（下側）を表わす。
【図４】泡の共振周波数に近い周波数を有する音響信号で励起されるときの励起信号（上側）とともに３μｍの泡からの散乱信号（下側）を表わす。
【図５】泡の共振周波数より上の周波数を有する音響信号で励起されるときの励起信号（上側）とともに３μｍの泡からの散乱信号（下側）を表わす。
【図６】本発明による方法において使用される典型的な増加掃引周波数信号を表わす。
【図７】図６における信号と同一の周波数帯域を含む、本発明による方法において使用される典型的な減少掃引周波数信号を表わす。
【図８】図６の励起信号へのマイクロバブルの応答を表わす。
【図９】図７の励起信号へのマイクロバブルの応答を表わす。
【図１０】送信された第１および第２の信号（上側）とともにマイクロバブルの応答（下側）を表わす。
【図１１】本発明による方法の組み合わせステップにおけるエコー信号の典型的な組み合わせを表わす。
【図１２】本発明による方法において使用される増加掃引周波数信号の線形変調および二次変調を表わす。
【図１３】本発明による装置を表わす。
【図１４】増加掃引周波数チャープ（上側）および減少掃引周波数チャープ（下側）を照射されたコントラストマイクロバブルの半径時間曲線を表わす。
【図１５】増加掃引周波数チャープ（上側）および減少掃引周波数チャープ（下側）への組織の疑似模型の応答を表わす。
【符号の説明】
【００４３】
１  放射源
２  アレイトランスデューサ
６  受信部
７  後処理回路
１１  目標物
１２  マイクロバブル
【図１】
2008543475.tif 000004
【図２】
2008543475.tif 000005
【図３】
2008543475.tif 000006
【図４】
2008543475.tif 000007
【図５】
2008543475.tif 000008
【図６】
2008543475.tif 000009
【図７】
2008543475.tif 000010
【図８】
2008543475.tif 000011
【図９】
2008543475.tif 000012
【図１０】
2008543475.tif 000013
【図１１】
2008543475.tif 000014
【図１２】
2008543475.tif 000015
【図１３】
2008543475.tif 000016
【図１４】
2008543475.tif 000017
【図１５】
2008543475.tif 000018
【国際調査報告】
2008543475.tif 000019
2008543475.tif 000020
2008543475.tif 000021