(57)【要約】  本発明は、関心領域を含む立体測定領域を撮像するための超音波システムであって、前記立体測定領域にわたって可変周波数範囲内の超音波ビームをステアリング可能な複数のＣＭＵＴトランスデューサのアレイを有するプローブと、超音波ビームステアリングを制御し、前記立体測定領域の超音波画像データを提供する、前記アレイに結合されたビームフォーマと、前記周波数範囲内で、前記ＣＭＵＴトランスデューサの動作周波数を変化させる、前記ビームフォーマに結合されたトランスデューサ周波数コントローラであって、前記周波数コントローラは、前記立体測定領域内でステアリングされる前記超音波ビームについて、前記動作周波数を第１の周波数に設定する、トランスデューサ周波数コントローラと、前記超音波画像データに基づき超音波画像を生成する画像プロセッサとを備える超音波システムに関する。システムはさらに、前記超音波画像データに基づき関心領域を識別することを可能にする関心領域（ＲＯＩ）識別部であって、前記立体測定領域内の前記関心領域を示す識別データを生成する、ＲＯＩ識別部を備え、前記トランスデューサ周波数コントローラはさらに、前記識別データに基づき、前記関心領域内でステアリングされる前記超音波ビームについて、前記動作周波数を第２の周波数に変更し、前記第２の周波数は、前記第１の周波数よりも高い。
2018527096.tif 000002
【特許請求の範囲】
【請求項１】
関心領域を含む立体測定領域を撮像するための超音波システムであって、
前記立体測定領域にわたって可変周波数範囲内の超音波ビームをステアリング可能な複数のＣＭＵＴトランスデューサのアレイを有する、腔内撮像に適したプローブと、
超音波ビームステアリングを制御し、前記立体測定領域の超音波画像データを提供する、前記アレイに結合されたビームフォーマと、
前記可変周波数範囲内で、前記ＣＭＵＴトランスデューサの動作周波数を変化させる、前記ビームフォーマに結合されたトランスデューサ周波数コントローラであって、前記周波数コントローラは、前記立体測定領域内でステアリングされる超音波ビームについて、前記動作周波数を第１の周波数に設定する、トランスデューサ周波数コントローラと、
前記超音波画像データに基づき超音波画像を生成する画像プロセッサと、  
前記超音波画像データに基づき関心領域を識別することを可能にする関心領域（ＲＯＩ）識別部であって、前記立体測定領域内の前記関心領域を示す識別データを生成する、ＲＯＩ識別部と、
前記プローブ及び前記ＲＯＩ識別部に結合され、動作中、前記プローブを動かすように作用する駆動機構とを備え、
前記トランスデューサ周波数コントローラはさらに、前記識別データに基づき、前記関心領域内でステアリングされる超音波ビームについて、前記動作周波数を第２の周波数に変更し、前記第２の周波数は、前記第１の周波数よりも高く、前記駆動機構は、前記識別データに基づき前記プローブを動かし、前記プローブと前記ＲＯＩとの間の距離変更を可能にする、超音波システム。
【請求項２】
前記ビームフォーマは、前記立体測定領域内で比較的低い空間分解能を有し、前記関心領域内で比較的高い空間分解能を有する前記超音波画像データを提供する、請求項１に記載の超音波システム。
【請求項３】
前記画像プロセッサは、前記低い空間分解能のデータに基づき前記立体測定領域のワイドビューを生成し、前記高い空間分解能のデータに基づき前記関心領域の詳細ビューを生成する、請求項２に記載の超音波システム。
【請求項４】
前記立体測定領域の前記ワイドビュー、及び前記関心領域の前記詳細ビューの両方を表示する、前記画像プロセッサに結合された画像ディスプレイをさらに備える、請求項３に記載の超音波システム。
【請求項５】
前記ディスプレイは、前記ワイドビュー及び前記詳細ビューの両方を空間的レジストレーションで表示する、請求項４に記載の超音波システム。
【請求項６】
前記立体測定領域内の前記ＲＯＩの手動選択に応答する、前記ＲＯＩ識別部に結合されたユーザインターフェイスをさらに備える、請求項１に記載の超音波システム。
【請求項７】
前記ユーザインタフェイスは、前記ディスプレイに関連づけられたタッチスクリーンである、請求項６に記載の超音波システム。
【請求項８】
前記関心領域の前記詳細ビューはリアルタイム更新される、請求項３に記載の超音波システム。
【請求項９】
前記ディスプレイはさらに、前記立体測定領域内の前記プローブの現在位置を表示する、請求項４に記載の超音波システム。
【請求項１０】
前記トランスデューサ周波数コントローラは、前記動作周波数を、前記ＲＯＩ内でステアリングされる超音波ビームについて前記第２の周波数に、前記関心領域外でステアリングされる超音波ビームについて前記第１の周波数に同時に設定し、前記第２の周波数は前記第１の周波数よりも高い、請求項１に記載の超音波システム。
【請求項１１】
前記アレイは１次元アレイである、請求項１に記載の超音波システム。
【請求項１２】
前記アレイは２次元アレイである、請求項１に記載の超音波システム。
【請求項１３】
関心領域を含む立体測定領域の可変周波数超音波撮像方法であって、
前記立体測定領域にわたって、可変周波数範囲内の超音波ビームをステアリングするステップと、
前記超音波ビームのステアリングを制御し、前記立体測定領域の超音波画像データを提供するステップであって、前記制御は、前記ステアリングされるビームの周波数を第１の周波数に設定することを含む、ステップと、
超音波画像を生成するために前記超音波画像データを処理するステップと、  
前記立体測定領域内の前記超音波画像データにおける前記関心領域を示す識別データを生成するステップと、
プローブから前記関心領域までの距離を計算するステップと、
前記識別データに基づき、前記関心領域に対して前記プローブを動かすことによって前記距離を縮めるステップと、
前記関心領域内でステアリングされる超音波ビームの周波数を、前記第１の周波数よりも高い第２の周波数に変更するステップとを含む、方法。
【請求項１４】
前記立体測定領域内で比較的低い空間分解能を有し、前記関心領域内で比較的高い空間分解能を有する前記超音波画像データを提供するステップをさらに含む、請求項１３に記載の方法。
【請求項１５】
前記低い空間分解能のデータに基づき前記立体測定領域のワイドビューを生成し、前記高い空間分解能のデータに基づき前記関心領域の詳細ビューを生成するステップを含む、請求項１４に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、関心領域を含む立体測定領域を撮像するための超音波システムであって、前記立体測定領域にわたって可変周波数範囲内の超音波ビームをステアリング可能な複数のＣＭＵＴトランスデューサのアレイを有する、腔内撮像に適したプローブと、超音波ビームステアリングを制御し、前記立体測定領域の超音波画像データを提供する、前記アレイに結合されたビームフォーマと、前記周波数範囲内で、前記ＣＭＵＴトランスデューサの動作周波数を変化させる、前記ビームフォーマに結合されたトランスデューサ周波数コントローラであって、前記周波数コントローラは、前記立体測定領域内でステアリングされる前記超音波ビームについて、前記動作周波数を第１の周波数に設定する、トランスデューサ周波数コントローラと、前記超音波画像データに基づき超音波画像を生成する画像プロセッサとを備える超音波システムに関する。
【０００２】
本発明は、さらに、かかる超音波システムを使用する立体測定領域の可変周波数超音波撮像方法に関する。
【背景技術】
【０００３】
ＣＭＵＴトランスデューサプローブを有する超音波撮像システムは、ＷＯ２０１５０２８３１４Ａ１から公知である。このプローブは、以下のモードのいずれかで動作するように構成された複数のＣＭＵＴセルを有するアレイを備える。１つ目は、ＤＣバイアス電圧によって、ＣＭＵＴの動作中、セルのＣＭＵＴ膜をセルフロア上で自由に振動させる従来のモードであり、２つ目は、ＤＣバイアス電圧によって、ＣＭＵＴセルの動作中、セルのＣＭＵＴ膜をセルフロアまでつぶすつぶれモードである。ＤＣバイアス電圧の増加によって、つぶれモードでの動作中、ＣＭＵＴセルの周波数応答の中心周波数は増加し、ＤＣバイアス電圧の減少によって、つぶれモードでの動作中、ＣＭＵＴセルの周波数応答の中心周波数は減少する。ＤＣバイアス電圧は、身体の立体領域が撮像される周波数に応じて、異なる臨床用途に合わせて選択され得る。
【０００４】
ＣＭＵＴ技術の視点をさらに活用する新しい撮像技術に対するニーズが存在する。
【発明の概要】
【０００５】
本発明の目的は、超音波イメージングにおける改善された能力を可能にする超音波システムを提供することである。
【０００６】
この目的は、本発明によれば、超音波画像データに基づき関心領域を識別することを可能にする関心領域（ＲＯＩ）識別部と、プローブ及びＲＯＩ識別部に結合された駆動機構とを提供することによって達成される。
【０００７】
前記識別部は、前記立体測定領域内の前記関心領域を示す識別データを生成し、前記トランスデューサ周波数コントローラはさらに、前記識別データに基づき、前記関心領域内でステアリングされる前記超音波ビームについて、前記動作周波数を第２の周波数に設定し、前記第２の周波数は、前記第１の周波数よりも高く、前記駆動機構は、前記識別データに基づき前記プローブを動かし、前記プローブと前記ＲＯＩとの間の距離を縮めることを可能にする。
【０００８】
本発明は、識別された関心領域内の超音波ビームの周波数を高めることを可能にする新たな撮像技術を提供にあたり、ＣＭＵＴトランスデューサの可変周波数能力を利用する。ＲＯＩ識別部によって超音波データにおいてＲＯＩが識別されると、トランスデューサ周波数コントローラは、ＲＯＩが位置する立体特定領域の一部においてビーム周波数を増加させる。さらに、システムは、領域内の識別されたＲＯＩの位置に応じて立体測定領域に対するプローブの位置を自動的に調整することにより、立体測定領域に対して腔内プローブを移動させることによって、超音波撮像中にユーザにさらなる柔軟性を与える。識別されたＲＯＩとプローブとの間の距離が、選択された高い周波数における音響ビームの貫通深度よりも大きい場合、システムは、高いビーム周波数を有するＲＯＩの詳細ビューが生成され得るために、プローブをＲＯＩに近づけるよう駆動機構を構成し得る。
【０００９】
一実施形態では、前記ビームフォーマは、前記立体測定領域内で比較的低い空間分解能を有し、前記関心領域内で比較的高い空間分解能を有する前記超音波画像データを提供する。
【００１０】
この実施形態では、関心領域にわたって送信されるビーム周波数を高めることによって、ビームフォーマが、より高い周波数のＲＯＩ由来エコー信号を受信することを可能にし、よって、識別されたＲＯＩのより高い分解能の超音波データを提供する。従来技術のシステムと比較して、本発明の超音波システムは、超音波走査中に、立体測定領域に関するより詳細な超音波情報を受信することができる。
【００１１】
他の実施形態では、前記画像プロセッサは、前記低い空間分解能のデータに基づき前記立体測定領域のワイドビューを生成し、前記高い空間分解能のデータに基づき前記関心領域の詳細ビューを生成する。
【００１２】
音波減衰は、周波数の上昇とともに増加する。したがって、高い貫通深度を有するが低い空間分解能を有する立体測定領域のワイドビューと、ワイド視野内の詳細視野とを生成することが有益であり、ＲＯＩは、より高い空間分解能で撮像され得る。本発明の利点は、単一の超音波走査の間に同じＣＭＵＴトランスデューサアレイを使用して、両方の視野を生成可能なことである。
【００１３】
他の実施形態では、超音波システムはさらに、前記立体測定領域の前記ワイドビュー、及び前記関心領域の前記詳細ビューの両方を表示する、前記画像プロセッサに結合された画像ディスプレイを備える。
【００１４】
両視野はユーザに対して、別々の超音波画像として互いに隣接して表示されてもよいし、又は１つの超音波画像として空間的レジストレーションを使用して表示されてもよい。
【００１５】
他の実施形態では、超音波システムは、前記立体測定領域内の前記ＲＯＩの手動選択に応答する、前記ＲＯＩ識別部に結合されたユーザインターフェイスをさらに備える。
【００１６】
これにより、ユーザは、ＲＯＩ識別部によって識別されるべきＲＯＩを手動で選択する機会を得る。オプションで、ユーザインターフェイスは周波数コントロールに結合されてもよく、この場合、ユーザはさらに、立体測定領域内及び関心領域内でそれぞれステアリングされるビームの比較的低い周波数及び高い周波数を選択し得る。
【００１７】
他の実施形態では、アレイは２次元アレイ又は１次元アレイである。
【００１８】
アレイのデザインに応じて、超音波システムは、立体測定領域の３次元超音波画像又は２次元超音波画像（２Ｄスライス）を提供し得る。
【図面の簡単な説明】
【００１９】
【図１ａ】図１ａは、本発明の原理に係る立体測定領域の可変周波数撮像のための超音波システムを示す。
【図２】図２は、ＤＣバイアス電圧によって制御され、ｒ．ｆ．駆動信号によって駆動されるＣＭＵＴセルを示す。
【図３ａ−３ｄ】図３ａ〜図３ｄは、本発明の実施形態に適用されるつぶれモードのＣＭＵＴ動作の原理を示す。
【図４】図４は、つぶれモードで動作するＣＭＵＴセルの音響性能の等高線図を示す。
【図５】図５は、関心領域外の立体測定領域内でステアリングされる超音波ビームの比較的低い周波数と、関心領域内でステアリングされる超音波ビームの比較的高い周波数とを用いた立体測定領域のスキャンを示す。
【図６ａ−６ｂ】図６ａ〜ｂは、立体測定領域の超音波画像の表示を、関心領域の詳細ビューを含む立体測定領域のワイドビューとともに示す。
【図７ａ−７ｂ】図７ａ〜ｂは、立体測定領域に対して移動されるように適合された腔内プローブを使用する、可変ビーム周波数を有する立体測定領域のスキャンを示す。
【図８】図８は、本発明の原理に係る腔内プローブを使用して得られた超音波画像の表示を示す。
【図９】図９は、本発明に係る関心領域を含む立体測定領域を撮像するための超音波システムを示す。
【図１０】図１０は、本発明の可変周波数画像取得の基本原理のワークフローを示す。
【図１１】図１１は、本発明に係る可変周波数画像取得のためのワークフローを示す。
【発明を実施するための形態】
【００２０】
図１は、本発明の原理に係る立体測定領域の可変周波数撮像のための超音波システム１００を概略的かつ例示的に示す。プローブ１０は、例えばＣＭＵＴ（ｃａｐａｃｉｔｉｖｅ  ｍｉｃｒｏｍａｃｈｉｎｅｄ  ｕｌｔｒａｓｏｕｎｄ  ｔｒａｎｓｄｕｃｅｒｓ）などの可変周波数超音波トランスデューサのアレイ１４を含み得る。このアレイ１４は、２次元アレイ又は１次元アレイのいずれであってもよい。アレイのＣＭＵＴは、立体測定視野１３１（図５）（立体測定領域を含む）にわたって、可変周波数範囲内の超音波ビームを送信し、送信されたビームに応じたエコーを受信する。アレイ１４のトランスデューサは、アレイトランスデューサ１４のＣＭＵＴによって送信される超音波ビームのステアリングを制御するビームフォーマ６４に結合される。ビームフォーマは、さらに、トランスデューサによって受信されたエコーをビーム成形する。ビームは、トランスデューサアレイ１４から直進する（アレイに対して直交する）よう方向づけられてもよいし、又は、より大きな視野のために異なる角度に方向づけられてもよい。オプションで、超音波システムは、それぞれがトランスデューサのグループをビームフォーマ６４と結合する、複数のマイクロビームフォーマ（図示せず）を有し得る。上記複数のマイクロビームフォーム（サブアレイビームフォーマ）は、トランスデューサのグループからの信号を部分的にビーム成形し、それにより、プローブとメイン取得システムとを結合する信号チャネルの量を減少させる。マイクロビームフォーマは、好ましくは、集積回路の形態で作成され、プローブ１０のハウジング内であって、アレイトランスデューサの近くに配置される。プローブ１０は、さらに、プローブ１０の位置を示す信号をトランスデューサ位置検出器５４に供給する位置センサ５２を含み得る。センサ５２は、磁気、電磁気、無線周波数、赤外線、又は他のタイプのセンサであり得る。
【００２１】
マイクロビームフォーマによって生成された部分的にビーム成形された信号は、ビームフォーマ６４に送られ、各トランスデューサグループからの部分的にビーム成形された信号が結合され、完全にビーム成形された信号が生成される。超音波システム１００は、さらに、ＣＭＵＴアレイ１４及びビームフォーマ６４（又はオプションで複数のマイクロビームフォーマ）に結合されたトランスデューサ周波数コントローラ６２を備える。周波数コントローラ６２は、以下でより詳細に説明するように、アレイ１４内の各ＣＭＵＴトランスデューサの共振周波数を調整することによって、送受信される超音波ビームの周波数を制御する。完全にビーム成形された信号（すなわち、ビームに沿うコヒーレントエコー信号）は超音波画像データを表し、超音波画像データは、フィルタリング、振幅検出、ドップラー信号検出、及び信号プロセッサ６６による他の処理によって処理される。その後、超音波データは、画像プロセッサ６８によってプローブの座標系（例えば、ｒ、θ、φ）における超音波画像信号に処理される。超音波画像信号は、さらに、グラフィックプロセッサ７４によって所望の超音波画像形式（例えば、ｘ、ｙ、ｚデカルト座標）に変換され、ディスプレイ１８上に表示されてもよい。
【００２２】
画像プロセッサ６８には関心領域識別部７２が結合され、関心領域識別部７２は、超音波画像データの解析に基づいて、立体測定視野１３１内の関心領域８２’（ＲＯＩ）を示す識別データを生成するように構成される。画像プロセッサ６８及びＲＯＩ識別部７２はともに、１つの画像分析ユニット６８’の一部とすることができる。超音波撮像システム１００は、ユーザインターフェイス３８によって制御され得る。特に、ユーザインターフェイス３８は、ＲＯＩ識別部７２に接続されたり、又は画像分析ユニット６８’に直接接続されてもよく、これは、ディスプレイ上に表示された超音波画像に基づきＲＯＩ８２’を手動選択することを可能にする。さらに、ユーザは、ユーザインターフェイス３８を介して、ユーザがＲＯＩを撮像するために使用することを望む、アレイの可変周波数範囲内の所望の周波数を選択する。立体測定視野１３１内のＲＯＩ８２’の位置及びサイズ並びに所望のＲＯＩ撮像周波数などの上記ユーザ入力は、識別データの形で、画像分析ユニット６８’によってトランスデューサ周波数コントローラ６２に伝達される。本実施形態では、ＲＯＩ識別部と画像プロセッサ６８との間でユーザによって特定されたパラメータが交換され、画像プロセッサは、ＲＯＩ識別部によって提供される生成された識別データに基づき、ＲＯＩ８２’の座標、及び立体測定視野１３１内の特定されたＲＯＩを囲む立体測定領域１３２の座標を計算する。トランスデューサ周波数コントローラ６２は、ＲＯＩ識別部７２によって生成されて、画像プロセッサ６８によって処理される識別データに応答する。トランスデューサ周波数コントローラ６２は、ビームフォーマと共に、立体測定視野１３１内の特定されたＲＯＩを囲う立体測定領域１３２内で操作されるビームの周波数を調整する。ビームフォーマ及びトランスデューサ周波数コントローラは、周波数変更能力とビーム成形能力とを組み合わせた１つの可変周波数ビームフォーマユニット６４’としてデザインされてもよい。上記代替の実施形態では、マイクロビームフォーマが、トランスデューサ周波数コントローラ６２とともに、可変周波数ビームフォーマユニット６４’に組み合わせられ、プローブのハウジング内に配置されてもよい。
【００２３】
本発明によれば、超音波システムの撮像周波数の変化は、つぶれモードで動作するように構成されたＣＭＵＴトランスデューサを使用して提供される。ＣＭＵＴ技術は、バイアス電圧を変化させることによって撮像周波数を調整することを可能にする。この周波数範囲は広範囲におよび、この範囲に加えて、各周波数にも帯域幅があり、帯域幅の大部分は１００％に近い。この大きな周波数可変性は、広範囲の貫通度及び分解能にわたる撮像を可能にする。
【００２４】
本発明のＣＭＵＴトランスデューサアレイ１４は、複数のＣＭＵＴセル（トランスデューサ）を含む。各ＣＭＵＴセル１０３は、典型的には、シリコン基板１１２の上に隙間又はキャビティ１１８を介して吊り下げられた柔軟な膜又はダイアフラム１１４を備える。ダイアフラム１１４上には上部電極１２０が配置され、ダイアフラムと共に移動する。この例では、基板１１２の上面であるセルの床には、底部電極が配置される。電極１２０のデザインには他の形態も考えられ、例えば、電極１２０は、膜１１４に埋め込まれてもよく、又は追加層として膜１１４上に堆積されてもよい。この例では、非限定的な例として、底部電極１２２は円形に構成され、基板層１１２内に埋め込まれている。他の適切な構成では、例えば、底部電極１２２は他の電極形状及び他の場所を有し、例えば、基板層１１２上に配置され、底部電極１１２は隙間１１８に直接露出されるか、又は、上部電極１２０と底部電極１２２との間のショートを防止するために電気絶縁層若しくはフィルムによって隙間１１８から隔離される。さらに、膜層１１４は、基板層１１２の上面に対して固定され、膜層１１４と基板層１１２との間に球状又は円筒状のキャビティ１１８を形成するように構成及び寸法デザインされる。明瞭さのために、図２では、非限定的な例として、底部電極１２２は接地されている。例えば、接地された上部電極１２０、又は上部電極１２０及び底部電極１２２の両方が浮遊であるなどの他の構成も当然に同様に実現可能である。
【００２５】
セル１００及びそのキャビティ１１８は、別の幾何学的形状を示してもよい。例えば、キャビティ１１８は、長方形若しくは正方形の断面、六角形の断面、楕円形の断面、又は不規則な断面を示すことができる。本明細書では、ＣＭＵＴセル１０３の直径への言及は、セルの最大横方向寸法として理解されるべきである。
【００２６】
底部電極１２２は、キャビティ側の面が追加層（図示せず）によって絶縁されていてもよい。好ましい電気的絶縁層は、基板電極１２２の上かつ膜電極１２０の下に形成されたＯＮＯ（ｏｘｉｄｅ−ｎｉｔｒｉｄｅ−ｏｘｉｄｅ）誘電体層であるが、この層には任意の電気絶縁材料が考えられる。ＯＮＯ誘電体層は、有利なことに、デバイスの不安定性、出力音圧のドリフト及び減少をもたらす電極上の電荷蓄積を軽減する。
【００２７】
ＣＭＵＴ上のＯＮＯ誘電体層の作成例は、２００８年９月１６日に出願された「Ｃａｐａｃｉｔｉｖｅ  ｍｉｃｒｏｍａｃｈｉｎｅｄ  ｕｌｔｒａｓｏｕｎｄ  ｔｒａｎｓｄｕｃｅｒ」という名称の、Ｋｌｏｏｔｗｉｊｋらによる欧州特許出願ＥＰ２，３２６，４３２Ａ２に詳細に説明されている。ＯＮＯ誘電体層の使用は、吊り下げられた膜で動作するＣＭＵＴよりも電荷保持の影響を受けやすい、予めつぶれたＣＭＵＴで望ましい。開示される構成要素は、例えば、Ａｌ、Ｔｉ、窒化物（例えば、窒化シリコン）、酸化物（種々のグレード）、テトラエチルオキシシラン（ＴＥＯＳ）、ポリシリコンなどのＣＭＯＳ適合性材料から形成され得る。ＣＭＯＳ作成では、例えば、酸化物層及び窒化物層を化学気相成長法により形成し、メタライゼーション（電極）層をスパッタリングによって形成してもよい。
【００２８】
適切なＣＭＯＳプロセスは、ＬＰＣＶＤ及びＰＥＣＶＤであり、後者は、４００℃未満の比較的低い作業温度を有する。開示されるキャビティ１１８を作成するための例示的な技術は、膜層１１４の上面を追加する前に、膜層１１４の初期部分にキャビティを画定することを含む。他の製造の詳細事項は、米国特許第６，３２８，６９７号（Ｆｒａｓｅｒ）において見つけることができる。
【００２９】
図２に示すように、円柱状キャビティ１１８の直径は、円形に構成された電極プレート１２２の直径よりも大きい。電極１２０は、円形に構成された電極プレート１２２と同じ外径を有してもよいが、そのような合致は必要ではない。したがって、膜電極１２０は、下方の電極プレート１２２と位置合わせされるように、膜層１１４の上面に対して固定され得る。ＣＭＵＴセル１００の電極は、装置の容量性プレートを提供し、隙間１１８は、キャパシタのプレート間の誘電体である。ダイアフラムが振動すると、プレート間の誘電体ギャップの寸法が変化し、これにより、受信された音響エコーに対するＣＭＵＴセル１００の応答として感知される変化する静電容量が提供される。
【００３０】
電極間の間隔は、電圧源４５を用いて静電圧、例えばＤＣバイアス電圧を電極に印加することによって制御される。電圧源４５は、トランスデューサ周波数コントローラ６２内に実装され、その周波数制御能力を提供する。アレイ１４の複数のトランスデューサはそれぞれ、別個の電圧源を有してもよく、又はトランスデューサ周波数コントローラ６２内に実装されたいくつかの電圧源を共有してもよい。電圧源４５はまた、オプションで、ＣＭＵＴセル１０３の駆動電圧のＤＣ及びＡＣ成分又は刺激成分をそれぞれ提供するための別個のステージ１０２、１０４を備えることができる。第１のステージ１０２は、静（ＤＣ）電圧成分を生成するように構成され、第２のステージ１０４は、所定の交流周波数を有する交流可変電圧成分又は刺激を生成するように構成されてもよく、後者の信号は通常、全体駆動電圧とそのうちの上記静的成分との間の差である。印加駆動電圧の静的又はバイアス成分は、好ましくは、ＣＭＵＴセル１０３をつぶれ状態にするための閾値電圧以上である。これは、第１のステージ１０２が、全体電圧の特に低ノイズの静的成分を生成するために、比較的大きなコンデンサ、例えば平滑化コンデンサを含み得るという利点を有し、静的成分は通常、全体電圧を支配し、全体電圧信号のノイズ特性は、この静的成分のノイズ特性によって支配される。電圧源４５の他の適切な実施形態が明らかであろう。例えば、ある実施形態では、電圧源４５は３つの別個のステージを含み、すなわち、ＣＭＵＴ駆動電圧の静的ＤＣ成分を生成するための第１のステージ、駆動電圧の可変ＤＣ成分を生成するための第２のステージ、及び信号の周波数変更又は刺激成分を生成するための第３のステージ、例えばパルス回路などを含む。要するに、電圧源４５は、任意の適切な態様で実施され得る。
【００３１】
それ自体は知られているように、一定の閾値を超える静電圧を印加することによって、ＣＭＵＴセル１０３は、膜１１４が基板１１２上につぶれるつぶれ状態にされる。この閾値は、ＣＭＵＴセル１０３の厳密なデザインに依存し、ＤＣバイアス電圧の印加中に膜１１４がファンデルワール力によってセルの床に付着（接触）する該バイアス電圧として定義される。膜１１４と基板１１２との間の接触の量（面積）は、印加されるバイアス電圧に依存する。膜１１４と基板１１２との間の接触面積を増加させると、膜１１４の共振周波数が上昇する。これについては、図３ａ〜ｄを参照してより詳細に説明する。
【００３２】
つぶれモードのＣＭＵＴセル１０３の周波数応答は、つぶれ後のＣＭＵＴ電極に印加されるＤＣバイアス電圧を調整することによって変化させられ得る。結果として、より高いＤＣバイアス電圧が電極に印加されるにつれ、ＣＭＵＴセルの共振周波数は増加する。この現象の背後にある原理を、図３ａ及び図３ｂに示す。図３ａ及び図３ｃの断面図は、これを、膜１１４の外側支持部分と、膜が各図面においてキャビティ１１８の床に接触し始める点との間の距離Ｄ１及びＤ２によって一次元的に説明する。距離Ｄ１は、比較的低いバイアス電圧が印加された場合の図３ａの比較的長い距離であり、一方、図３ｃの距離Ｄ２は、より高いバイアス電圧が印加されたことによる、はるかに短い距離である。これらの距離は、両端で保持され、つま弾かれた長い糸と短い糸に類似する。長いゆるんだ糸は、はじかれると、より短く、張っている糸よりもはるかに低い周波数で振動する。同様に、図３ａのＣＭＵＴセルの共振周波数は、より高い引き下げバイアス電圧を受ける図３ｃのＣＭＵＴセルの共振周波数よりも低くなる。
【００３３】
実際には、ＣＭＵＴ膜の有効動作領域の関数であるので、この現象は、図３ｂ及び図３ｄの二次元図からも理解することができる。図３ａに示すように、膜１１４がＣＭＵＴセルの床にちょうど接触するとき、セル膜１１４の非接触（自由振動）部分の有効振動面積Ａ１は、図３ｂに示すように大きい。中心１７の小さい穴は、膜の中央接触領域を表す。大面積の膜は、比較的低い周波数で振動する。この領域１７は、ＣＭＵＴセルのフロアに潰れた膜１１４の領域である。しかし、図３ｃに示すように、より高いバイアス電圧によって膜がより深いつぶれ状態に引きさげられると、図３ｄに示すように、中央接触領域１７’がより大きくなり、自由振動領域Ａ２はより小さくなる。このより小さな領域Ａ２は、大きい領域Ａ１よりも高い周波数で振動する。したがって、ＤＣバイアス電圧が減少すると、つぶれたＣＭＵＴセルの周波数応答が減少し、ＤＣバイアス電圧が増加すると、つぶれたＣＭＵＴセルの周波数応答が増加する。
【００３４】
図４は、送信中の、ＡＣ変調、又は定周波数の周波数変調の形態の刺激を含む印加されるＤＣバイアス電圧の関数として、つぶれモードにおける典型的なＣＭＵＴセル１０３の音圧出力の等高線図を示す。対応するパルス長は、印加周波数の半分である。この等高線図から分かるように、ＣＭＵＴセル１０３が固定又は静的電圧、例えば固定値のＤＣバイアス電圧で動作するとき、わずかな周波数範囲でしか、最適な音響性能が得られない。しかし、等高線図に破線で示すように、バイアス電圧信号、及びバイアス電圧信号上の周波数変調を相関して変化させると、ＣＭＵＴセル１０３の最適な音響性能は、はるかに大きな周波数範囲にわたって達成され、これにより、ＣＭＵＴセル１０３を含む超音波プローブの送信モードで生成される超音波パルス（又はパルストレイン）の有効帯域幅を増加させることができる。したがって、周波数を、この例のように７〜１７ＭＨｚ、３〜１０ＭＨｚ、又はさらには２〜１５ＭＨｚのような広い周波数範囲で変化させることができる。
【００３５】
これは、つぶれ状態のＣＭＵＴセル１０３の共振周波数は、印加される（ＤＣ）バイアス電圧に依存することが説明されている図３ａ及び図３ｄを参照することによって理解することができる。適切な設定周波数を有する刺激を印加することによって特定の設定周波数の超音波パルスを発生させる際に印加されるバイアス電圧を調整することにより、パルス周波数毎にＣＭＵＴセル１０３の（近）最適音響性能を示す、異なる周波数のパルスを生成することができる。したがって、これは、撮像スペクトルの広い帯域幅にわたって、（近）最適な撮像分解能を保証する。
【００３６】
音波減衰は周波数の増加に伴って増加し、一方、超音波画像の解像度は周波数の増加と共に減少する。例えば、組織内の２サイクルパルスに関する典型的な奥行き（深度）及び軸方向分解能を下表に示す。
【表１】
2018527096.tif 000003
【００３７】
最適要件及び貫通（ｐｅｎｅｔｒａｔｉｏｎ）要件を合理的に満たすために、大部分の診断用途の周波数範囲は２〜１５ＭＨｚである。範囲の低い部分は、高い深度（例えば、関心領域が身体のより深い位置にある）又は高い減衰（例えば、経頭蓋研究において）の場合に有用である。周波数範囲の高い部分は、小さい貫通が必要な場合（例えば、乳房、甲状腺、若しくは表面血管の撮像又は小児撮像の場合）に有用である。ほとんどの大柄な患者の場合、３〜５ＭＨｚが十分な周波数であり、一方、細い患者及び小児では、５及び７．５ＭＨｚがよく使用され得る。１５ＭＨｚを超えるより高い周波数範囲は、ＩＶＵＳ、ＩＣＥ、ＦＬ−ＩＣＥなどの腔内（血管内）プローブを使用して高解像度画像を提供することができる。これらのプローブは、体腔、血管内でＲＯＩの近くに配置することができる。
【００３８】
本発明は、単一のＣＭＵＴトランスデューサアレイ１４を使用する、広周波数範囲における可変周波数超音波撮像のユニークな組み合わせを提供する。
【００３９】
図５は、プローブの位置がＲＯＩ８２’に対して固定される、本発明の基本原理を示す。プローブ１０は、立体測定視野１３１の超音波画像を取得するために使用される。トランスデューサ周波数コントローラ６２は、関心領域識別部７２に応答し、立体測定視野１３１内で操作される超音波ビームの比較的低い周波数、及び、特定されたＲＯＩ８２’の周囲の立体測定領域１３２内で操作される超音波ビームの比較的高い周波数を設定する。ＣＭＵＴエコーによって受信されたエコーは、ビームフォーマによって処理され、ＲＯＩ外の立体測定領域内で比較的低い空間分解能を有し、関心領域内で比較的高い空間分解能を有する立体測定領域の超音波画像データが提供される。これらの超音波データは画像プロセッサ６８で処理され、図６に示されるように、低い空間分解能のデータに基づく立体測定領域のワイドビュー８０、及び高空間分解能データに基づく関心領域８２の詳細ビュー１３２’が生成される。特定されたＲＯＩ８２の周囲の立体測定領域１３２の詳細ビュー１３２’は、プローブとＲＯＩとの間に位置する領域の画像１３３も含み得る。
【００４０】
図６は、ユーザに表示される２Ｄ超音波画像の表示９９を示し、ワイドビュー８０及び詳細ビュー１３２’は、互いに空間的に位置合わせされている（ｓｐａｔｉａｌ  ｒｅｇｉｓｔｒａｔｉｏｎ）。選択されたＲＯＩ８２’の表現８２は、詳細ビュー１３２’における増加した撮像周波数で表示される。比較的高い周波数を有する超音波ビームの貫通深度は、比較的低い周波数を有する超音波ビームの貫通深度に比べて低いので、比較的高い周波数範囲の上限周波数は、ＲＯＩが位置する深さ（プローブまでの距離）によって制限され、計算の際に画像プロセッサ６８によって考慮される。システム１００は、まず、比較的低いビーム周波数を有する立体測定視野の超音波データを取得し、立体測定領域の周辺コンテキストを提供し、さらに、特定後、ＲＯＩ８２に「ズームイン」する。ＲＯＩ８２の詳細ビュー１３２’は、図６ｃに示すように、先に取得され、コンテキストのために表示されたワイドビュー８０の隣にリアルタイムで更新され得る。
【００４１】
あるいは、ＲＯＩ８２の詳細ビュー１３２’及びワイドビュー８０は、互いに隣接して表示されてもよい。心臓撮像時の心臓学的用途において、超音波画像の表示及び取得は、ＥＣＧ同期（ＥＣＧ  ｇａｔｉｎｇ）によって心周期と同期させられてもよい。
【００４２】
ＣＭＵＴアレイ１４が線形アレイである場合、トランスデューサ周波数コントローラ６２は、ＲＯＩが高周波数で撮像される一方、他の要素は低周波数に維持されるよう、異なる周波数で個々のトランスデューサセル１０３を制御（駆動）することができる。線形アレイによって得られた模式的画像を図６ｂに示す。
【００４３】
リアルタイム低周波数ワイドビュー８０画像と同時に、埋め込みリアルタイム高周波数詳細ビュー１３２’画像が生成される。これは、周囲コンテキストも比較的低い深度でリアルタイムで撮像され（たとえより低い分解能であっても）、例えばＲＯＩの周辺で生じるツールの方向付け及びナビゲーションを可能にするという利点を有する。また、図６ａ及び図６ｃに示すように、ＣＭＵＴアレイ１４がフェーズドアレイである場合、同様の画像を得ることも可能である。フェーズドアレイの場合、画像を構成する各ラインについて、低周波数ラインを含むワイドビュー８０画像内に高周波数詳細ビュー１３２’画像が埋め込まれるように、全てのトランスデューサに適切な周波数が選択されるように、ビーム成形が実行される。ＲＯＩ８２の詳細ビュー１３２’及びワイドビュー８０の両方がリアルタイムで更新される場合、フェーズドアレイを含むシステムは、まず、低周波数で立体測定視野１３１ボリュームの全てのラインを取得し、続いて、より高い周波数で特定されたＲＯＩ８２を取り囲む立体測定領域１３２の全てのラインを取得することができる。取得されたビューは、さらに、インターリーブ又は補間されて１つの超音波画像にされてもよい。これは、図６ｃに示されている。他の取得ワークフローでは、ワイドビュー８０は詳細ビュー１３２’を過ぎて更新され、この場合、ユーザに表示される結果画像は図６ａに示されるようなものである。前者は、例えば、介入デバイスを追跡するための、全ボリュームのリアルタイムビューという利点を有する。後者は、より少ないラインが取得され、より高いフレームレートが達成され得るという利点を有する。
【００４４】
図７は、本発明の一実施形態を示し、プローブの位置を、立体測定視野１３１’内で変更することが可能である。プローブは、例えば、プローブがＲＯＩに対して容易に前進及び後進可能であるように、前向き又はエンドファイア構成で配置され得る。これは、例えばＥＰ１７４２５８０Ｂに記載されているように、ＩＶＵＳ（ｉｎｔｒａｖａｓｃｕｌａｒ  ｕｌｔｒａｓｏｕｎｄ）、ＩＣＥ（ｉｎｔｒａｃａｒｄｉａｃ  ｅｃｈｏｃａｒｄｉｏｇｒａｐｈｙ）、ＦＬ−ＩＣＥ（ｆｏｒｗａｒｄ  ｌｏｏｋｉｎｇ  ｉｎｔｒａｃａｒｄｉａｃ  ｅｃｈｏｃａｒｄｉｏｇｒａｐｈｙ）などの腔内プローブを提供することによって実現されてもよい。
【００４５】
腔内プローブは、立体測定領域を走査するためにスイープされる遠位先端におけるトランスデューサアレイを含み得る。ボリュームスイープは、１Ｄアレイの機械的移動、又は２Ｄアレイを用いたビームの電子的操作することによって提供され得る。トランスデューサアレイは、プローブの遠位先端に配置された流体チャンバ内に収容され、流体は、プローブと撮像される立体測定領域との間の適切な音響的結合を提供する。図９には、本発明の超音波システム１００’が概略的に示されている。システム１００’は、さらに、プローブ及びＲＯＩ識別部７２に（オプションで分析ユニット６８’に）結合された駆動機構２１を備え、駆動機構は識別（特定）データに基づき、撮像中にプローブ１０を動かすように作用する。駆動機構２１はまた、プローブの空間的位置を追跡する位置センサ５２から信号を受信し、よって、立体測定視野１３１’内でのプローブの移動を提供する。本実施形態は、ＲＯＩ８２’の撮像用の高周波数の上限に、より高い柔軟性を与える。ＲＯＩが識別された後、画像プロセッサ６８は、識別データに基づき、ＲＯＩ８２の座標、及び立体測定視野１３１内の識別されたＲＯＩを囲む立体測定領域１３２の座標を計算する。トランスデューサアレイ１４（又は実質的にはプローブ１０）とＲＯＩとの間の距離が、選択された高周波数のビームの貫通深度を超える場合、駆動機構２１は、ＲＯＩの「ズームイン」画像を取得するために、立体測定視野１３１’内のＲＯＩに近づくよう伝達される（図７ｂ）。
【００４６】
本発明は、（ＣＭＯＳ製造の進歩によって実現される）小型化ＣＭＵＴトランスデューサ、及び（つぶれ動作モードによって実現される）かかるトランスデューサの広い動作帯域を、駆動装置へのフィードバックループと組み合わせ、立体測定領域における自動ズームイン／アウトを実行可能な新世代の超音波システムをユーザに提供する。つぶれモードで動作するＣＭＵＴアレイの広い周波数帯域と、該アレイを含むプローブを物理的に並進移動させる手段との組み合わせは、より詳細な、よって改良された診断イメージングによって、高度な超音波イメージングにおける新たなユーザエクスペリエンスを可能にする。
【００４７】
図８は、ユーザに表示される２Ｄ超音波画像のディスプレイ９９を示す。詳細ビュー８２及びワイドビュー８０は、互いに隣接して、又は空間的に位置合わせされて示され得る。図８には後者の場合が示されており、線形アレイ及びフェーズドアレイによって得られた画像が、互いに隣接して配置されている。図６ａ〜ｂと比較すると、詳細ビュー８２は、プローブの位置がＲＯＩの位置に対して固定されている実施形態と比較して、より大きな貫通深度を有するものとして、ユーザの目に映る。ワイドビュー画像８０が、異なる時点で取得されたより高分解能の詳細ビュー画像８２でリアルタイム更新されるよう、詳細ビュー画像がプローブの進行（移動）中に連続的に取得され得る。また、立体測定領域内の移動しているプローブの現在位置が、ワイドビュー画像８０内に表示され得る。
【００４８】
画像プロセッサ６８は、ＲＯＩ識別パラメータ及びユーザ特定パラメータに基づいて、軸方向ノイズ、横方向スペックル、軸方向強度などの画質パラメータについて得られた超音波データを分析し得る。これらの品質パラメータは、ＲＯＩ画質の自動最適化のためのフィードバックループの一部となるよう、プローブを自動的に動かすための駆動機構への入力としても使用され得る。このような自動化は、プローブの細かい移動のために使用され、一方、全体的な運動はユーザインターフェイスを介して制御され得る。ユーザインターフェイス３８を介して、ユーザは、駆動機構の動作に関する追加制御を提供し得る。ユーザインターフェイスは、ユーザが、ＲＯＩ及び／又はプローブの動きを表示画像において手動で定めることを可能にする、ディスプレイに関連付けられたタッチスクリーンであってもよい。ＲＯＩに触れたり、「ピンチイン」又は「ピンチアウト」の動作をしたりすることで、プローブを特定の方向に物理的に移動させたり、又は所与のプローブ位置について貫通深度が十分な場合には詳細画像が取得されるように構成されてもよい。
【００４９】
他の実施形態では、比較的高い周波数で得られたＲＯＩのリアルタイム詳細３Ｄビューが、ワイドビュー２Ｄ画像内に埋め込まれる。これは、ワイドビュー２Ｄ画像の取得には、より少ない処理能力トランスデューサ利用が消費され、３Ｄ画像（又はバイプレーンＲＯＩ）を可能な限り高いフレームレートで取得できるという利点を有する。１次元の小さいアパーチャを有するアレイ（例えば、ＩＣＥ）の場合、本実施形態は、より良好なアパーチャ次元（ＩＣＥ軸方向及び横方向）に基づくワイドビュー撮像、及び全ての次元（例えば、ＩＣＥ：高さ（ｅｌｅｖａｔｉｏｎ）を含む）における詳細ＲＯＩ撮像を提供し、これは高周波数でより好ましくなる。
【００５０】
ＲＯＩ識別部は、超音波増強造影特徴によってマーキングされてもよいカテーテル、針、又はツールなどの特定の物体からの超音波データを使用して自動的にＲＯＩを識別し得る。これらの物体は、それらの形状及びアスペクト（又はマーカー若しくは位置センサ）に基づき画像解析ユニット６８’によって認識され、そして、ＲＯＩの座標が自動的に生成され得る。
【００５１】
他の実施形態では、最初に、比較的高い周波数のビームで関心ボリュームの画像が取得され、この関心ボリュームが、ユーザによってＲＯＩとして特定され得る。さらに、ユーザはユーザインターフェイスを介して、貫通深度がより大きいワイドビュー画像（ＲＯＩを含む）を得るために、ＲＯＩに使用されたものと比較して撮像周波数を低下させ得る。前述してきた実施形態と同様に、これらの視野は互いに隣接して表示されてもよいし、又は空間的レジストレーションによって表示されてもよい。
【００５２】
ＣＭＵＴに印加されるバイアス電圧の変化の最適速度を提供するために、可変周波数ビームフォーマユニット６４’の（又はオプションで、トランスデューサ周波数コントローラ６２の）集積回路（ＩＣ）エレクトロニクスに対して、別々の要件が課されてもよい。上記の場合のほとんどについて、現在のＩＣエレクトロニクス技術で十分であり得る。あるいは、さらに高速のバイアス電圧変化が必要な場合には、ＷＯ／２０１５／０８６４１３に記載されているような３端子ＣＭＵＴを使用することができる。
【００５３】
図１０には、可変周波数画像取得の基本原理のワークフロー２００が示されている。ステップ２０１において、立体測定視野１３１が撮像され、この視野は、ワイドビュー８０を含む。ステップ２０２において、ＲＯＩ８２が識別部によって検出され、自動検出は、例えば、特徴的な解剖学的特徴３０に基づいて、又はユーザ入力に基づいて実行され得る。ステップ２０３において、ＲＯＩのアウトラインがユーザに表示され得る。このステージにおいて、ユーザは、ユーザインターフェイス３８を介してシステム１００手動でインタラクトし、ＲＯＩのサイズ及び／又は位置を調整してもよい。さらに、ステップ２０４において、ユーザは、ＲＯＩの詳細ビューの所望の解像度（又は周波数）を選択し得る。画像プロセッサ６８はさらに、選択された解像度をトランスデューサ動作周波数に変換する。あるいは、このステップにおいて、画像プロセッサ６８は、プローブ１０（すなわち、プローブ内のトランスデューサアレイ１４）からＲＯＩまでの固定距離に基づきＲＯＩ８２を撮像可能な上限周波数を計算し得る。この情報はディスプレイ上に表示され得る。ステップ２０５において、システム１００は、より高い解像度でＲＯＩの詳細ビューを取得する。ステップ２０６において、ワイド視野及び詳細視野がユーザに表示される。
【００５４】
図１１には、本発明に係る可変周波数画像取得のためのワークフロー３００が示されている。ステップ３０１において、立体測定視野１３１が取得される。ステップ３０２において、ＲＯＩ８２が識別部によって検出される。ステップ３０３において、ＲＯＩのアウトラインがユーザに表示され得る。このステージにおいて、ユーザは、ユーザインターフェイス３８を介してシステム１００手動でインタラクトし、ＲＯＩのサイズ及び／又は位置を調整してもよい。並行して、ステップ３０７において、画像プロセッサ６８は、プローブからＲＯＩの最も遠い端部までの距離を計算する。さらに、ステップ３０４において、ユーザは、ＲＯＩの詳細ビューの所望の解像度（又は周波数）を選択し得る。ステップ３０９において、この情報に基づいて、画像プロセッサ６８は、選択された解像度（周波数）に対応する貫通深度を計算する。ステップ３０８において、プローブとＲＯＩとの間の距離が、貫通深度と比較される。計算された貫通深度がＲＯＩまでの距離よりも大きい場合、ワークフローはステップ３０５に進み、システム１００は、選択された解像度でＲＯＩの詳細ビューを取得する。計算された貫通深度がＲＯＩまでの距離よりも小さい場合、ワークフローはステップ３１０に進み、駆動機構は、プローブをＲＯＩの位置に近づける。移動距離は、ＲＯＩの位置と選択された解像度とによって決定される。移動距離が、撮像されるボリューム（物体）の解剖学的構造によって制限され、プローブがそれ以上移動できない場合、システム１００’は、解剖学的制限を考慮に入れてＲＯＩを取得可能な最適解像度を計算し、ユーザにフィードバックし得る。さらに、システム１００は、ステップ３０５において、選択された解像度又は提案された最適解像度でＲＯＩの詳細ビューを取得する。ステップ３０６において、ワイド視野及び詳細視野がユーザに表示される。
【００５５】
当業者は、本発明の原理が２Ｄ及び３Ｄ超音波イメージングの両方において実施可能であることを理解するであろう。
【００５６】
単一のユニットまたはデバイスが、請求項に記載される複数のアイテムの機能を果たし得る。複数の手段が互いに異なる従属請求項に記載されているからといって、これらの手段の組み合わせが好適に使用することができないとは限らない。
【００５７】
コンピュータプログラムは、他のハードウェアと共に又は他のハードウェアの一部として供給される光学記憶媒体又はソリッドステート媒体等の適切な媒体上で記憶及び／又は分配されてもよいし、インターネット又は他の有線若しくは無線テレコミュニケーションシステムを介して等の他の形態で分配されてもよい。
【００５８】
特許請求の範囲内のいかなる参照符号も、その範囲を限定するものと解釈されるべきではない。
【図１】
2018527096.tif 000004
【図２】
2018527096.tif 000005
【図３ａ】
2018527096.tif 000006
【図３ｂ】
2018527096.tif 000007
【図３ｃ】
2018527096.tif 000008
【図３ｄ】
2018527096.tif 000009
【図４】
2018527096.tif 000010
【図５】
2018527096.tif 000011
【図６Ａ】
2018527096.tif 000012
【図６Ｂ】
2018527096.tif 000013
【図６Ｃ】
2018527096.tif 000014
【図７Ａ】
2018527096.tif 000015
【図７Ｂ】
2018527096.tif 000016
【図８】
2018527096.tif 000017
【図９】
2018527096.tif 000018
【図１０】
2018527096.tif 000019
【図１１】
2018527096.tif 000020
【国際調査報告】
2018527096.tif 000021
2018527096.tif 000022
2018527096.tif 000023