(57)【要約】  医療装置（４００）は、少なくとも２つの経路（６０６、１１１４）に沿った超音波の速さを記述する超音波データを収集する超音波送信器（４４４、６０２）・受信器（４４６、６０４）システム（６００）を有する。医療装置は更に、医用画像データを収集する医用撮像システム（４０２）と、プロセッサに前記医用画像データを収集（１００、２００）させる命令を格納したメモリ（４６４）とを有する。前記命令はプロセッサに更に前記超音波データを収集させる（１０２、２０２）。前記命令はプロセッサに更に、前記医用画像データを少なくとも２つの組織型にセグメント化させる（１０４、２０４）。前記命令はプロセッサに更に、被検体内の前記少なくとも２つの経路に対応する少なくとも２つの距離を決定させる（１０６、２０６）。前記命令はプロセッサに更に、前記少なくとも２つの組織型内での超音波の速さを計算させる（１０８、２０８）。
2014517704.tif 000002
【特許請求の範囲】
【請求項１】
−  超音波データを収集する超音波送信器・受信器システムであり、前記超音波データは、被検体内の少なくとも２つの経路に沿った超音波の速さを記述する、超音波送信器・受信器システム；
−  前記被検体内の前記少なくとも２つの経路を記述する医用画像データを収集する医用撮像システム；
−  当該医療装置を制御するプロセッサ；
−  機械実行可能命令を格納したメモリであり、前記命令の実行が前記プロセッサに：
−  前記医用撮像システムを用いて前記医用画像データを収集させ；
−  前記超音波送信器・受信器システムを用いて前記超音波データを収集させ；
−  前記医用画像データをセグメント化させて少なくとも２つの組織型を特定させ；
−  前記被検体内の前記少なくとも２つの経路の各々に関して前記被検体内の少なくとも２つの距離を決定させ、ただし、前記少なくとも２つの距離は、前記超音波送信器・受信器システムによって生成された超音波が前記少なくとも２つの組織型内で進行した距離に等しく；
−  前記超音波データと、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離とに従って、前記少なくとも２つの組織型の各々内での超音波の速さを計算させる；
メモリ；
を有し、
当該医療装置は更に、前記被検体内の標的ゾーンを超音波破砕する高密度焦点式超音波システムを有し、前記高密度焦点式超音波システムは、超音波エネルギーを生成する超音波トランスデューサを有し、前記超音波トランスデューサは複数の超音波トランスデューサ素子を有し、前記複数の超音波トランスデューサ素子の各々の位相が制御可能であり、前記命令の実行は前記プロセッサに更に、前記セグメント化された医用画像データと、前記少なくとも２つの組織型の各々内での超音波の速さとに従って、一組のトランスデューサ素子位相を計算させる、
医療装置。
【請求項２】
当該医療装置はタンクを有し、前記タンクは、超音波伝導液と前記被検体の一部とを受け入れるための開口を有し、前記一部は前記標的ゾーンを含み、前記超音波送信器・受信器システムは、主として前記タンクを通る第１の平面内で、超音波を送信して受信するように適応され、前記第１の平面は前記開口の下方にあり、前記複数の超音波トランスデューサ素子は前記第１の平面の下方の第２の平面の下方にある、請求項１に記載の医療装置。
【請求項３】
前記超音波トランスデューサは前記第２の平面の上方に焦点を合わされる、請求項２に記載の医療装置。
【請求項４】
前記超音波トランスデューサは前記第２の平面の下方に焦点を合わされる、請求項２に記載の医療装置。
【請求項５】
前記少なくとも２つの経路は前記被検体を貫通している、請求項２乃至４の何れか一項に記載の医療装置。
【請求項６】
前記命令の実行は前記プロセッサに更に：
−  前記超音波送信器・受信器システムを用いて基準超音波データを収集させ、ただし、前記基準超音波データは、前記超音波伝導液を通る少なくとも２つの基準経路に沿って収集され、前記基準超音波データは超音波振幅測定結果を有し、前記超音波データは超音波振幅測定結果を有し；且つ
−  前記少なくとも２つ組織型の各々に関して減衰値を計算させ、ただし、前記減衰値は、前記超音波データの前記超音波振幅測定結果及び前記基準超音波データの前記超音波振幅測定結果と、前記少なくとも２つ経路の各々に関する前記少なくとも２つの距離と、に従って計算される；
請求項５に記載の医療装置。
【請求項７】
当該医療装置は更に、前記被検体及び／又は前記超音波送信器・受信器システムを位置決めする位置決めシステムを有し、前記命令の実行は前記プロセッサに更に：
−  前記セグメント化された医用画像データ内で骨組織を特定させ；且つ
−  前記少なくとも２つの経路が前記骨組織で反射されて前記受信器システムによって受信されるように、前記被検体及び／又は前記超音波送信器・受信器システムを位置決めさせる；
請求項１乃至４の何れか一項に記載の医療装置。
【請求項８】
前記命令の実行は前記プロセッサに更に、前記少なくとも２つの組織型の各々内での超音波減衰を、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離に従って計算させる、請求項１乃至７の何れか一項に記載の医療装置。
【請求項９】
前記命令の実行は前記プロセッサに更に：
−  前記高密度焦点式超音波システムに前記標的ゾーンを超音波破砕させるための超音波制御信号を、前記一組のトランスデューサ素子位相に従って生成させ；且つ
−  前記超音波制御信号を前記高密度焦点式超音波システムへと送信させる；
請求項１乃至８の何れか一項に記載の医療装置。
【請求項１０】
前記医用撮像システムは、撮像ゾーンから磁気共鳴データを収集する磁気共鳴撮像システムであり、前記標的ゾーンは前記撮像ゾーン内にあり、前記命令の実行は前記プロセッサに更に：
−  前記超音波トランスデューサから前記標的ゾーンまで生成される経路に沿って磁気共鳴サーモメトリデータを収集させ；
−  前記セグメント化された医用画像データと前記磁気共鳴サーモメトリデータとに従って、一組の温度補償されたトランスデューサ素子位相を計算させ；
−  前記一組の温度補償されたトランスデューサ素子位相に従って、温度補正された超音波制御信号を生成させ；且つ
−  前記温度補正された超音波制御信号を前記高密度焦点式超音波システムへと送信させる；
請求項９に記載の医療装置。
【請求項１１】
前記医用撮像システムは、磁気共鳴撮像システム、コンピュータ断層撮影システム、及び診断超音波システムのうちの何れか１つである、請求項１乃至９の何れか一項に記載の医療装置。
【請求項１２】
医療装置の作動方法であって、前記医療装置は、超音波データを収集する超音波送信器・受信器システムを有し、前記超音波データは、被検体内の少なくとも２つの経路に沿った超音波の速さを記述し、前記医療装置は更に、前記被検体内の前記少なくとも２つの経路を記述する医用画像データを収集する医用撮像システムを有し、当該方法は：
−  前記医用撮像システムを用いて前記医用画像データを収集するステップ；
−  前記超音波送信器・受信器システムを用いて前記超音波データを収集するステップ；
−  前記医用画像データを少なくとも２つの組織型にセグメント化するステップ；
−  超音波により、前記被検体内の前記少なくとも２つの経路の各々に関して前記被検体内の少なくとも２つの距離を決定するステップであり、前記少なくとも２つの距離は、前記超音波送信器・受信器システムによって生成された超音波が前記少なくとも２つの組織型内で進行した距離に相当する、ステップ；及び
−  前記超音波データと、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離とに従って、前記少なくとも２つの組織型の各々内での超音波の速さを計算するステップ；
を有し、
前記医療装置は更に、前記被検体内の標的ゾーンを超音波破砕する高密度焦点式超音波システムを有し、前記高密度焦点式超音波システムは、超音波エネルギーを生成する超音波トランスデューサを有し、前記超音波トランスデューサは複数の超音波トランスデューサ素子を有し、前記複数の超音波トランスデューサ素子の各々の位相が制御可能であり、当該方法は更に、前記セグメント化された医用画像データと前記少なくとも２つの組織型の各々内での超音波の速さとに従って一組のトランスデューサ素子位相を計算するステップを有する、
医療装置の作動方法。
【請求項１３】
医療装置を制御するプロセッサにより実行される機械実行命令を有するコンピュータプログラムであって、前記医療装置は、超音波データを収集する超音波送信器・受信器システムを有し、前記超音波データは、被検体内の少なくとも２つの経路に沿った超音波の速さを記述し、前記医療装置は更に、前記被検体内の前記少なくとも２つの経路を記述する医用画像データを収集する医用撮像システムを有し、前記命令の実行は前記プロセッサに：
−  前記医用撮像システムを用いて前記医用画像データを収集させ；
−  前記超音波送信器・受信器システムを用いて前記超音波データを収集させ；
−  前記医用画像データを少なくとも２つの組織型にセグメント化させ；
−  超音波により、前記被検体内の前記少なくとも２つの経路の各々に関して前記被検体内の少なくとも２つの距離を決定させ、ただし、前記少なくとも２つの距離は、前記超音波送信器・受信器システムによって生成された超音波が前記少なくとも２つの組織型内で進行した距離に相当し；
−  前記超音波データと、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離とに従って、前記少なくとも２つの組織型の各々内での超音波の速さを計算させ；且つ
超音波エネルギーを生成する超音波トランスデューサを有する高密度焦点式超音波システムを制御させ、前記超音波トランスデューサは複数の超音波トランスデューサ素子を有し、前記複数の超音波トランスデューサ素子の各々の位相が制御可能であり、前記セグメント化された医用画像データと、前記少なくとも２つの組織型の各々内での超音波の速さとに従って、一組のトランスデューサ素子位相が計算される；
コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用撮像に関し、特に、超音波と医用撮像とを用いた組織内の超音波の速度の計算に関する。
【背景技術】
【０００２】
焦点式超音波トランスデューサからの超音波は、人体の内部の領域を選択的に治療するために使用されることができる。超音波は高エネルギーの機械的振動として伝送される。これらの振動は減衰されるときに組織加熱を誘起するとともに、キャビテーションをもたらし得る。組織加熱及びキャビテーションはどちらも、或る臨床設定において、組織を破壊するために使用されることができる。しかしながら、超音波を用いて組織を加熱することは、キャビテーションより制御が容易である。超音波治療は、組織を切除するため、また、癌細胞の領域を選択的に消し去るために使用され得る。この技術は、子宮筋腫の治療に適用されており、子宮摘出術の必要性を減らしてきている。より低いパワー又はパルスモードにおいて、超音波は、遺伝物質又は薬剤を或る領域に選択的に送達するために使用されることができる。
【０００３】
超音波治療を行うため、特定の治療ボリューム（容積部）に超音波を集束させるよう、焦点式超音波トランスデューサを使用することができる。トランスデューサは典型的に、超音波を伝達することができる例えば脱気水などの媒体内に搭載される。そして、アクチュエータを用いて、超音波トランスデューサの位置が調整され、それにより、処置される組織領域が調整される。
【０００４】
特許文献１は、組織を通じて超音波エネルギーを集束させるシステムと、組織領域を撮像するイメージャとを示している。超音波エネルギーを集束させるために、位相及び振幅の補正係数が用いられている。それぞれの補正係数が、それらの組織型（tissue  type）内での組織型の音速の変化を補償する。それらの補正係数は、画像を分析してそれらの画像から組織特性を特定することによって決定されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】米国特許出願公開第２００４／０１２２３２３号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は、独立請求項にて、医療装置、方法、及びコンピュータプログラムを提供する。実施形態が従属請求項にて与えられる。
【課題を解決するための手段】
【０００７】
高密度焦点式超音波（High  Intensity  Focused  Ultrasound；ＨＩＦＵ）による乳癌の治療において、超音波は、例えば１４５１±３６ｍ／ｓ及び１５３８ｍ／ｓ±２２ｍ／ｓの超音波の波速によってそれぞれ特徴付けられる皮下脂肪及び腺組織など、幾つかの媒体を横切って進行する。ビーム経路内で、これら２つの組織型の１ｃｍ厚の変化は、ビームの焦点をずらすのに十分な２１０°の位相収差を誘起する。ＨＩＦＵ治療はしばしば、癌を位置特定し、治療を計画し、そして温熱療法（ハイパーサーミア）手順をモニタするため、磁気共鳴撮像（ＭＲＩ）を組み合わされる。
【０００８】
ＭＲＩ又はその他の撮像モダリティによって提供される解剖学的画像はまた、組織の波速及び減衰の先験性に基づいて位相及び振幅の逸脱補償を処理するために使用され得る。しかしながら、乳房組織の波速は（リフォーカシング（再焦点合わせ）を計算するのに必須のパラメータ）は患者及び年齢に依存する。本発明の実施形態は、医用撮像システムの内部に超音波の放射器（エミッタ）及び受信器（レシーバ）を付加することによって治療前又は治療中に組織型の波速及び減衰を特徴付ける装置及び／又は方法を提供し得る。超音波放射器及び受信器は、数少ない方向に沿って超音波の伝搬時間及び減衰を測定するために使用され、解剖学的画像が、それらの方向に沿っての各組織型の特徴付けを提供する。故に、波速及び減衰の値は、双方のデータを組み合わせることによって得られる。そのような組織特徴付けは、温熱療法の前又は最中に迅速に得られるものであり、より良好な治療効率及び精度のためにビームリフォーカシングを向上させることになる。
【０００９】
ＭＲ誘導式ＨＩＦＵによる乳癌の治療は、最近、臨床実用に入る大きな可能性を示している。しかしながら、超音波ビームの遠視野内に位置する胸郭、心臓及び肺への望ましくない組織損傷のかなりのリスクによって、幅広い臨床的支持が阻まれている。本発明の実施形態は、横方向の主ビーム方向を有するトランスデューサ設計の使用によって、これらのリスクを回避する手法を提供し得る。さらに、ビームの近視野における低めのエネルギー密度が、皮下組織層の望ましくない温度上昇に関する安全限界を超えることなく一層大きい腫瘍を治療することを可能にする。
【００１０】
大口径トランスデューサ設計のフォーカルポイント（焦点）の品質は、一様でない伝搬媒体の影響を遙かに受けやすい。乳房組織は、音響的観点において、６％もの音響伝搬速度の変動を示す腺組織、脂肪組織及び繊維組織の複合構造であるので、このような設計の使用は一般的には低質なビームをもたらし得る。
【００１１】
本発明の実施形態は、ビーム品質を回復する補正機能を提供し得る。全ての患者に対して同等の組織音響特性を仮定する補正関数の計算の基礎として高分解能ＣＴ撮像を使用する経頭蓋ＨＩＦＵアブレーションシステムの設計も同様の問題に遭遇している。そうとはいえ、それらの方法は、骨組織中のビーム伝搬に起因する収差の特徴付け及び補正を行うように最適化されているので、乳房組織に直接的に移転可能なものではない。また、マルチモダリティ（ＣＴ及びＭＲＩ）の使用は、システムの利用可能性、作業継続時間及び画像レジストレーション（整合）に関する問題を誘起する。故に、ＨＩＦＵを用いた乳癌の最適な非侵襲治療に大口径トランスデューサ設計を使用することを可能にする鍵となるのは、臨床的に意義のある時間フレーム（時間枠）で乳房組織の音響特性を特徴付ける効率的で非侵襲的な方法である。
【００１２】
ＭＲ撮像は、優れた軟組織コントラストを示し、故に、不均質組織構造を高い空間分解能でマッピングすることができる。この情報と、複数の異なる乳房組織型の体内（インビボ）波速測定の所定のテーブル（表）との組合せを用いて、リフォーカシング方法を処理し得る。しかしながら、皮下脂肪や腺組織の波速が患者及び年齢に依存するため、この方法には限界がある。
【００１３】
乳房組織の波速を定量化するために、飛行時間式超音波断層撮影プラットフォームが提案されている。しかしながら、それによる測定は、超音波の現状技術に制約されたコントラスト及び空間分解能で長いスキャン時間（品目のスライス当たり９分）を要する。本発明の実施形態は、超音波組織特性を特徴付ける手段を、（超音波波速測定以外の）１つ又は複数の撮像モダリティと組み合わせることによって、より患者特定的な、ＨＩＦＵ治療の組織不均質性補償を可能にし得る。
【００１４】
超音波放射器及び受信器を、患者の解剖学的画像と組み合わせて用いることで、超音波速さ及び超音波減衰を特徴付けることができる。解剖学的画像は、例えば：
・磁気共鳴撮像（ＭＲＩ）
・Ｘ腺写真
−  蛍光透視
−  Ｘ線投影
・核医学
−  シンチグラフィ
−  ＳＰＥＣＴ
−  陽電子放出断層撮影（ＰＥＴ）
・光音響撮像
・胸部サーモグラフィ
・断層撮影（トモグラフィ）
−  リニアトモグラフィ
−  ポリトモグラフィ
−  ゾノグラフィ
−  オルトパントモグラフィ（ＯＰＴ又はＯＰＧ）
−  コンピュータトモグラフィ（ＣＴ）
−  コンピュータアキシャルトモグラフィ（ＣＡＴ）
・超音波検査
などのうちの１つの撮像モダリティを用いて、あるいは複数の撮像モダリティを組み合わせて取得され得る。
【００１５】
超音波放射器及び受信器は、複数の組織層にわたっての伝搬遅延及び減衰を測定することができる。また、別の撮像モダリティは、超音波の伝搬方向に沿って、以下の特性：
・複数の組織層の各々の厚さ
・組織密度分布
・温度分布
の定量化を提供することができる。
【００１６】
超音波放射器及び受信器による測定と別の撮像モダリティとを組合せることは、最適なリフォーカシングを処理するのに必要な、患者特定的な、各組織型の超音波速さ及び減衰の見積もりを向上させることができる。
【００１７】
本発明の実施形態により、組織の特徴付けは、患者スクリーニング段階中に処理され、あるいはプランニング段階中に処理され、あるいは有意な手順遅延なくＨＩＦＵ治療中にオンザフライで処理され得る。超音波エコーの収集は、大抵のその他の医用撮像モダリティの収集時間より遙かに短い数ミリ秒を必要とするのみである。また、各システムが同一のプラットフォームに統合され、且つ収集システムが相互間での干渉を回避するように設計されている場合、超音波エコーの収集は、その他の撮像モダリティの使用と同時に行われ得る。あるいは、データは、何らかの順序で別々に収集されて、必要な全てのデータの収集後に処理工程に向けて共に組み合わされてもよい。
【００１８】
例えば、治療をプランニングするために最初に収集されたＭＲＩデータが、選択された標的領域に対する収差補正を処理するために再利用され得るので、それに関連する収集期間は追加期間として考えなくてよい。
【００１９】
しかしながら、生体構造（組織厚さ）を測定するとともに、同じセッション中に同時あるいは順次に超音波エコーを収集することは、それらが自動的に同じ基準フレーム内にあるという利点が存在する。組織厚さの測定とエコーの収集とが２つの異なるセッションで別々に行われる場合には、飛行時間を調べるために使用される超音波エコーと解剖学的情報とを共にレジストレーションするために、双方の測定中の患者位置を完全に確信する必要がある。
【００２０】
乳癌の治療は、乳房組織の容易な特徴付けのために超音波放射器と受信器とを直線に沿って互いに向き合わせることができるので、本発明の実施形態の好適な一用途である。また、大口径トランスデューサを用いた乳癌の治療は、組織収差に一層敏感であり、この用途を更に本質的なものにする。
【００２１】
‘コンピュータ読み取り可能記憶媒体’は、ここでは、コンピューティング装置のプロセッサによって実行可能な命令を格納し得る如何なる有形記憶媒体をも包含するものである。コンピュータ読み取り可能記憶媒体は、コンピュータ読み取り可能持続的記憶媒体とも呼ばれ得る。コンピュータ読み取り可能記憶媒体はまた、有形コンピュータ読み取り可能媒体とも呼ばれ得る。一部の実施形態において、コンピュータ読み取り可能記憶媒体はまた、コンピューティング装置のプロセッサによってアクセスされることが可能なデータを格納し得る。コンピュータ読み取り可能記憶媒体の例は、以下に限られないが、フロッピー（登録商標）ディスク、磁気ハードディスクドライブ、ソリッドステートハードディスク、フラッシュメモリ、ＵＳＢサムドライブ、ランダムアクセスメモリ（ＲＡＭ）、読み出し専用メモリ（ＲＯＭ）、光ディスク、磁気光ディスク、プロセッサのレジスタファイルを含む。光ディスクの例は、例えばＣＤ−ＲＯＭ、ＣＤ−ＲＷ、ＣＤ−Ｒ、ＤＶＤ−ＲＯＭ、ＤＶＤ−ＲＷ、又はＤＶＤ−Ｒといった、コンパクトディスク（ＣＤ）及びデジタル多用途ディスク（ＤＶＤ）を含む。コンピュータ読み取り可能記憶媒体なる用語はまた、ネットワーク又は通信リンクを介してコンピュータ装置によってアクセスされることが可能な様々な種類の記録媒体をも意味する。例えば、データは、モデム上、インターネット上、又はローカルエリアネットワーク上で取り出され得る。
【００２２】
‘コンピュータメモリ’又は‘メモリ’はコンピュータ読み取り可能記憶媒体の一例である。コンピュータメモリは、プロセッサにとって直接的にアクセス可能なメモリである。コンピュータメモリの例は、以下に限られないが、ＲＡＭメモリ、レジスタ、及びレジスタファイルを含む。
【００２３】
‘コンピュータストレージ’又は‘ストレージ’はコンピュータ読み取り可能記憶媒体の一例である。コンピュータストレージは不揮発性コンピュータ読み取り可能記憶媒体である。コンピュータストレージの例は、以下に限られないが、ハードディスクドライブ、ＵＳＢサムメモリ、フロッピー（登録商標）ドライブ、スマートカード、ＤＶＤ、ＣＤ−ＲＯＭ、及びソリッドステートハードドライブを含む。一部の実施形態において、コンピュータストレージはコンピュータメモリであってもよく、その逆もまた然りである。
【００２４】
‘プロセッサ’は、ここでは、プログラム又は機械実行可能命令を実行することができる電子部品を包含するものである。“プロセッサ”を有するコンピューティング装置への言及は、２つ以上のプロセッサ又は処理コアを含む場合があるものとして解釈されるべきである。プロセッサは例えばマルチコアプロセッサとし得る。プロセッサはまた、単一のコンピュータシステム内の、あるいは複数のコンピュータシステム間で分散された、複数のプロセッサの集合を意味し得る。コンピューティング装置なる用語も、各々が１つ以上のプロセッサを有する複数のコンピューティング装置の集合若しくはネットワークを意味する場合があるとして解釈されるべきである。多くのプログラムは、同一のコンピューティング装置内とし得る複数のプロセッサ、又は複数のコンピューティング装置にまたがって分散され得る複数のプロセッサ、によって実行される命令を有する。
【００２５】
‘ユーザインタフェース’は、ここでは、ユーザ又はオペレータがコンピュータ又はコンピュータシステムとインタラクトすることを可能にするインタフェースである。‘ユーザインタフェース’はまた、‘ヒューマンインタフェース装置’とも呼ばれ得る。ユーザインタフェースは、オペレータに情報又はデータを提供し、且つ／或いはオペレータから情報又はデータを受信し得る。ユーザインタフェースは、オペレータからの入力がコンピュータによって受信されることを可能にし得るとともに、コンピュータからの出力をユーザに提供し得る。換言すれば、ユーザインタフェースは、オペレータがコンピュータを制御あるいは操作することを可能にし得るとともに、コンピュータがオペレータの制御又は操作の効果を指し示すことを可能にし得る。ディスプレイ又はグラフィカルユーザインタフェース上でのデータ又は情報の表示は、オペレータに情報を提供することの一例である。キーボード、マウス、トラックボール、タッチパッド、ポインティングスティック、グラフィックタブレット、ジョイスティック、ゲームパッド、ウェブカム、ヘッドセット、ギアスティック、ステアリングホイール、ペダル、配線付きグローブ、ダンスパッド、リモートコントローラ、及び加速度計を介したデータの受信は全て、オペレータからの情報又はデータの受信を可能にするユーザインタフェースコンポーネントの例である。
【００２６】
‘ハードウェアインタフェース’は、ここでは、コンピュータシステムのプロセッサが外部のコンピューティング装置及び／又は機器とインタラクトする、あるいはそれらを制御する、ことを可能にするインタフェースを包含するものである。ハードウェアインタフェースは、プロセッサが外部コンピューティング装置及び／又は機器に制御信号又は命令を送信することを可能にし得る。ハードウェアインタフェースはまた、プロセッサが外部コンピューティング装置及び／又は機器とデータを交換することを可能にし得る。ハードウェアインタフェースの例は、以下に限られないが、ユニバーサルシリアルバス、ＩＥＥＥ１３９４ポート、パラレルポート、ＩＥＥＥ１２８４ポート、ＲＳ−２３２ポート、ＩＥＥＥ−４８８ポート、ブルートゥース接続、無線ローカルエリアネットワーク接続、ＴＣＰ／ＩＰ接続、イーサネット（登録商標）接続、制御電圧インタフェース、ＭＩＤＩインタフェース、アナログ入力インタフェース、及びデジタル入力インタフェースを含む。
【００２７】
‘医用画像データ’は、医用撮像システムを用いて収集された１次元、２次元又は３次元のデータを包含するものである。‘医用撮像システム’は、ここでは、患者の身体構造についての情報を収集し、１次元、２次元又は３次元の医用画像データのセットを構築するように適応された装置を包含するものである。医用画像データは、医師による診断に有用な視覚化を構築するために使用され得る。この視覚化はコンピュータを用いて実行されることができる。
【００２８】
磁気共鳴（ＭＲ）データは、ここでは、磁気共鳴撮像スキャン中に磁気共鳴装置のアンテナによって記録される、原子スピンにより放射される無線周波数信号の測定結果として定義される。磁気共鳴撮像画像は、ここでは、磁気共鳴撮像データ内に含まれる解剖学的データの、再構成された１次元、２次元又は３次元視覚化として定義される。この視覚化は、プロセッサ又はコンピュータを用いて実行されることができる。磁気共鳴画像は医用画像データの一例である。
【００２９】
‘超音波データ’は、ここでは、対象物中を伝送される超音波ビームを記述するデータを包含するものである。超音波データは、以下に限られないが、送信器による送信時間、放射される超音波の周波数、受信器による受信時間、送信と受信との間の遅延、送信振幅、反射振幅、及び受信器によって測定された超音波の振幅、のうちの如何なる組み合わせをも有していてもよい。例えば、超音波が骨又はその他の硬組織構造で跳ね返されたり跳ね飛ばされたりする場合、反射の振幅及び遅延が特に関心あるものとなり得る。
【００３０】
‘超音波の速さ’は、ここでは、固体媒体又は液体媒体中を超音波が伝播する速度を記述する量を包含するものである。超音波の速さなる用語は、超音波の波速、音の速さ、音速度、音速、超音波速、及び超音波速度と同義である。
【００３１】
‘超音波ウィンドウ’は、ここでは、超音波の波動又はエネルギーを伝送することが可能なウィンドウ（窓）を包含するものである。典型的に、薄膜又はメンブレインが超音波ウィンドウとして使用される。超音波ウィンドウは例えば、ＢｏＰＥＴ（二軸延伸ポリエチレンテレフタレート）の薄いメンブレインで製造され得る。
【００３２】
一態様において、本発明は、超音波データを収集する超音波送信器・受信器システムを有する医療装置を提供する。超音波データは、被検体内の少なくとも２つの経路に沿った超音波の速さを記述する。超音波送信器は超音波を放射し、受信器は超音波を受ける。例えば、超音波のパルスが、超音波送信器によって放射され、そして、超音波受信器によって受信される。受信された超音波パルスの振幅及び／又は遅延が、超音波受信器によって決定され得る。
【００３３】
医療装置は更に、被検体内の前記少なくとも２つの経路を記述する医用画像データを収集する医用撮像システムを有する。医療装置は更に、当該医療装置を制御するプロセッサを有する。医療装置は更に、機械実行可能命令を格納したメモリを有する。
【００３４】
機械実行可能命令の実行は、プロセッサに、医用撮像システムを用いて医用画像データを収集させる。例えば、プロセッサが制御信号を医用撮像システムに送信し、医用画像データが収集されるように制御する。機械実行可能命令の実行はプロセッサに更に、超音波送信器・受信器システムを用いて超音波データを収集させる。医用画像データと超音波データとが収集される順序は関係ない。例えば、先ず医用画像データが収集されてもよいし、先ず超音波データが収集されてもよいし、医用画像データ及び超音波データが同時に収集されてもよい。
【００３５】
機械実行可能命令の実行はプロセッサに更に、医用画像データを少なくとも２つの組織型にセグメント化させる。医用画像データは、被検体の生体構造を記述し得る。少なくとも２つの異なる種類の組織が、医用画像データ内で識別あるいはセグメント化される。医用画像データは、標準的なセグメンテーションアルゴリズムを用いてセグメント化されてよい。例えば、医用画像データの解剖学的データ内のランドマークが特定され得る。他の一例において、変形可能モデルが医用画像データにフィッティングされる。機械実行可能命令の実行はプロセッサに更に、被検体内の前記少なくとも２つの経路の各々に関して被検体内の少なくとも２つの距離を決定させる。該少なくとも２つの距離は、超音波送信器・受信器システムによって生成された超音波が前記少なくとも２つの組織型内で進行したそれぞれの距離に相当する。機械実行可能命令の実行はプロセッサに更に、超音波データと、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離とに従って、前記少なくとも２つの組織型の各々内での超音波の速さを計算させる。
【００３６】
この実施形態は、被検体の複数の組織型内での音又は超音波の速度を正確に決定するために使用され得るので、特に有利である。例えば、被検体の一部が２種類以上の組織を含む場合、異なる種類の組織内での音の速度は異なることが予想され得る。異なる組織領域における超音波の速度を正確に決定することにより、より正確に超音波を集束させ得る。これは、ワイドアパーチャトランスデューサの場合のように個々のフェイズドアレイトランスデューサ素子のビーム経路が大きく異なるときに特に有益である。複数の経路に沿って超音波が放射され、そして、医用画像データを用いて、複数の異なる組織型の各々内で超音波がどれだけ長くどれだけの距離だけ進行するかが決定される。これら距離の知見を、超音波が各経路を横断するのにどれだけの時間を要するかの知見と組み合わせることで、超音波の速度を計算することが可能になる。例えば、超音波データは、超音波パルスが超音波送信器によって送信されたのと、それが受信器によって受信されたのとの間の遅延を含み得る。
【００３７】
この実施形態はまた、組織内の超音波速度は診断的重要性を有するので有利となり得る。
【００３８】
他の一実施形態において、超音波送信器・受信器システムは第１の座標系を有し、医用撮像システムは第２の座標系を有する。第１の座標系は第２の座標系にレジストレーションされる。医用撮像システムと超音波送信器・受信器システムとによるデータ収集を相互に関連付けることができるので、これは有利である。これは、音速の正確な決定に必要な１つの共通座標系へのデータ融合を可能にする。
【００３９】
他の一実施形態において、医療装置は更に、被検体内の標的ゾーンを超音波破砕（ソニケーション）する高密度焦点式超音波システムを有する。高密度焦点式超音波システムは、超音波エネルギーを生成する超音波トランスデューサを有する。超音波トランスデューサは複数の超音波トランスデューサ素子を有する。複数の超音波トランスデューサ素子の各々の位相は制御可能である。例えば、超音波トランスデューサ素子の各々が交流電源によって駆動され得る。複数の超音波トランスデューサ素子の各々に供給される電力の位相及び／又は振幅が制御可能にされ得る。
【００４０】
命令の実行はプロセッサに更に、セグメント化された医用画像データと、前記少なくとも２つの組織型の各々内での超音波の速さとに従って、一組のトランスデューサ素子位相を計算させる。例えば、或る特定の位置を超音波破砕することが望ましい場合、この位置が、セグメント化された医用画像データ内で特定され得る。高密度焦点式超音波トランスデューサの位置と２つの組織型に対する超音波の速さとの間の関係の知見により、一組のトランスデューサ素子の各々の適切な位相が、高密度焦点式超音波が適切に集束されるように選定されることを可能にする。これを行わない場合には、標的ゾーンでの位相は、複数の異なる組織型で様々な異なる超音波速度が存在することによる影響を受けてしまい得る。この実施形態は、集束される超音波が標的ゾーン内に正確にフォーカシングされることを可能にするので、特に有利である。
【００４１】
他の一実施形態において、医療装置はタンクを有する。タンクは、超音波伝導液と被検体の一部とを受け入れるための開口を有する。タンクは例えば、超音波伝導液を充填されることが可能な容器又はコンテナとし得る。超音波伝導液は、ここでは、超音波エネルギーを伝達するように適応された如何なる液体をも包含するものである。超音波伝導液は、例えば、ゲル、オイル又は水とし得る。タンクが超音波伝導液で部分的に充たされ、そして、被検体の一部がタンク内に置かれる。該一部は標的ゾーンを含む。超音波送信器・受信器システムは、主としてタンクを通る第１の平面内で、超音波を送信して受信するように適応される。主としてとは、超音波送信器・受信器システムの超音波送信器及び受信器が、超音波送信器によって生成される超音波の大部分が第１の平面中を進行するように配置されることを意味する。超音波の大部分が、第１の平面中を進行するように生成される。第１の平面は前記開口の下方にある。前記複数の超音波トランスデューサ素子は第２の平面の下方にあり、第２の平面は第１の平面の下方にある。
【００４２】
他の一実施形態において、超音波トランスデューサは第２の平面の上方に焦点を合わされる。この実施形態は、超音波トランスデューサが上向きに、場合によりタンクの頂部又は外に、焦点合わせされるので有利である。
【００４３】
他の一実施形態において、超音波トランスデューサは第２の平面の下方に焦点を合わされる。この実施形態においては、ソニケーション領域の上方で超音波データが収集される。
【００４４】
他の一実施形態において、前記少なくとも２つの経路は被検体を貫通している。
【００４５】
他の一実施形態において、命令の実行はプロセッサに更に、超音波送信器・受信器システムを用いて基準超音波データを収集させる。基準超音波データは、超音波伝導液を通る少なくとも２つの基準経路に沿って収集される。基準超音波データは超音波振幅測定結果を有する。上記超音波データは超音波振幅測定結果を有する。命令の実行はプロセッサに更に、少なくとも２つ組織型の各々に関して減衰値を計算させる。減衰値は、超音波データの超音波振幅測定結果及び基準超音波データの超音波振幅測定結果と、前記少なくとも２つ経路の各々に関する前記少なくとも２つの距離とに従って計算される。この情報を用いて、各組織型における減衰を計算することができる。
【００４６】
他の一実施形態において、医療装置は更に、被検体及び／又は超音波送信器・受信器システムを位置決めする位置決めシステムを有する。命令の実行はプロセッサに更に、セグメント化された医用画像データ内で骨組織を特定させる。命令の実行はプロセッサに更に、前記少なくとも２つの経路が骨組織で反射されて受信器システムによって受信されるように、被検体及び／又は超音波送信器・受信器システムを位置決めさせる。例えば、プロセッサは、超音波送信器・受信器システムの位置を作動させる機械的作動（アクチュエイティング）システムに接続され得る。他の例では、プロセッサは、被検体の位置を制御することが可能な患者位置決めシステムを制御することができてもよい。この実施形態において、超音波は被検体に入り、骨によって反射される。例えば、このような一実施形態において、送信器及び受信器は、被検体が接触している超音波ウィンドウに対して位置決めされ得る。この実施形態は、超音波が通ることができない被検体領域の付近で組織中での超音波速度及び／又は組織内での超音波減衰を決定することができるので、非常に有利となり得る。例えば、超音波ビームを胸郭に通すことは困難である。
【００４７】
他の一実施形態において、命令の実行はプロセッサに更に、前記少なくとも２つの組織型の各々内での超音波減衰を、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離に従って計算させる。
【００４８】
他の一実施形態において、命令の実行はプロセッサに更に、前記一組のトランスデューサ素子位相に従って超音波制御信号を生成させる。ここでは、超音波制御信号とは、高密度焦点式超音波システムを制御することが可能なコマンド及び／又は信号のことである。超音波制御信号は、高密度焦点式超音波システムに標的ゾーンを超音波破砕させる。命令の実行はプロセッサに更に、超音波制御信号を高密度焦点式超音波システムへと送信させる。
【００４９】
他の一実施形態において、医用撮像システムは、撮像ゾーンから磁気共鳴データを収集する磁気共鳴撮像システムである。標的ゾーンは撮像ゾーン内にある。命令の実行はプロセッサに更に、超音波トランスデューサから標的ゾーンまで生成される経路に沿って磁気共鳴サーモメトリ（温度測定）データを収集させる。命令の実行はプロセッサに更に、セグメント化された医用画像データと磁気共鳴サーモメトリデータとに従って、一組の温度補償されたトランスデューサ素子位相を計算させる。
【００５０】
命令の実行はプロセッサに更に、一組の温度補償されたトランスデューサ素子位相に従って、温度補正された超音波制御信号を生成させる。命令の実行はプロセッサに更に、温度補正された超音波制御信号を高密度焦点式超音波システムへと送信させる。この実施形態は、前記少なくとも２つの組織型中での超音波速度が温度補償され得るので、特に有利である。身体組織の温度は比較的一定のままであるが、高密度焦点式超音波システムを用いて標的ゾーンを加熱するとき、ホットスポットが発達し得る。さらに、標的ゾーン内のソニケーションポイントも熱くなって超音波の速度を変化させ得る。これは、僅かながらも焦点ぼけ（デフォーカシング）を生じさせ得る。超音波トランスデューサの近視野における加熱も焦点ぼけを生じさせ得る。
【００５１】
他の一実施形態において、医用撮像システムは磁気共鳴撮像システムである。この実施形態において、医用画像データは磁気共鳴画像データセットとし得る。磁気共鳴撮像システムは磁気共鳴データを収集し、それを磁気共鳴画像データへと再構成し得る。
【００５２】
他の一実施形態において、医用撮像システムはコンピュータ断層撮影システムである。
【００５３】
他の一実施形態において、医用撮像システムは診断超音波システムである。
【００５４】
他の一実施形態において、前記少なくとも２つの経路は、被検体を通る複数あるいは多数の経路である。統計手法を用いて、被検体内の前記少なくとも２つの組織型を通る超音波の振幅及び／又は速度を正確に決定し得る。
【００５５】
他の一態様において、本発明は、医療装置の作動方法を提供する。当該方法はまた、コンピュータにより実行される方法を提供し得る。医療装置は、超音波データを収集する超音波送信器・受信器システムを有する。超音波データは、被検体内の少なくとも２つの経路に沿った超音波の速さを記述する。医療装置は更に、被検体内の前記少なくとも２つの経路を記述する医用画像データを収集する医用撮像システムを有する。当該方法は、医用撮像システムを用いて医用画像データを収集するステップを有する。当該方法は更に、超音波送信器・受信器システムを用いて超音波データを収集するステップを有する。当該方法は更に、医用画像データを少なくとも２つの組織型にセグメント化するステップを有する。当該方法は更に、超音波により、被検体内の前記少なくとも２つの経路の各々に関して被検体内の少なくとも２つの距離を決定するステップを有する。該少なくとも２つの距離は、超音波送信器・受信器システムによって生成された超音波が前記少なくとも２つの組織型内で進行した距離に相当する。当該方法は更に、超音波データと、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離とに従って、前記少なくとも２つの組織型の各々内での超音波の速さを計算するステップを有する。
【００５６】
医療装置は更に、被検体内の標的ゾーンを超音波破砕する高密度焦点式超音波システムを有する。高密度焦点式超音波システムは、超音波エネルギーを生成する超音波トランスデューサを有する。超音波トランスデューサは複数の超音波トランスデューサ素子を有する。複数の超音波トランスデューサ素子の各々の位相及び／又は振幅が制御可能である。命令の実行はプロセッサに更に、セグメント化された医用画像データと前記少なくとも２つの組織型の各々内での超音波の速さとに従って一組のトランスデューサ素子の位相及び／又は振幅を計算させる。
【００５７】
他の一態様において、本発明は、医療装置を制御するプロセッサにより実行される機械実行命令を有するコンピュータプログラムを提供する。当該コンピュータプログラムは例えば、コンピュータ読み取り可能媒体に格納され得る。医療装置は、超音波データを収集する超音波送信器・受信器システムを有する。超音波データは、被検体内あるいは被検体を貫通する少なくとも２つの経路に沿った超音波の速さを記述する。医療装置は更に、被検体内の前記少なくとも２つの経路を記述する医用画像データを収集する医用撮像システムを有する。
【００５８】
命令の実行はプロセッサに、医用撮像システムを用いて医用画像データを収集させる。命令の実行はプロセッサに更に、超音波送信器・受信器システムを用いて超音波データを収集させる。命令の実行はプロセッサに更に、医用画像データを少なくとも２つの組織型にセグメント化させる。命令の実行はプロセッサに更に、超音波により、被検体内の前記少なくとも２つの経路の各々に関して被検体内の少なくとも２つの距離を決定させる。該少なくとも２つの距離は、超音波送信器・受信器システムによって生成された超音波が前記少なくとも２つの組織型内で進行した距離に相当する。命令の実行はプロセッサに更に、超音波データと、前記少なくとも２つの経路の各々に関する前記少なくとも２つの距離とに従って、前記少なくとも２つの組織型の各々内での超音波の速さを計算させる。
【図面の簡単な説明】
【００５９】
以下、以下の図を含む図面を参照して、単なる例として、本発明の好適実施形態を説明する。
【図１】本発明の一実施形態に係る方法を例示するフロー図である。
【図２】本発明の一実施形態に係る方法の更なる例を示す図である。
【図３】人間の乳房内の相異なる組織型に関する超音波の波速の差を例示するプロットである。
【図４】本発明の一実施形態に係る医療装置を例示する図である。
【図５】図４の実施形態の高密度焦点式超音波システムをより詳細に例示する図である。
【図６】本発明の一実施形態に係る超音波送信器・受信器システムの一例を示す図である。
【図７】本発明の更なる一実施形態に係る超音波送信器・受信器システムを示す図である。
【図８】本発明の一実施形態に係る超音波トランスデューサ素子の構成を示す等角図である。
【図９】図８の実施形態の超音波トランスデューサ素子の構成の上面図である。
【図１０】図８の実施形態の超音波トランスデューサ素子の構成の側面図である。
【図１１】本発明の一実施形態に係る高密度焦点式超音波システムの超音波トランスデューサに対する超音波送信器・受信器システムの相対位置を例示する図である。
【図１２】本発明の更なる一実施形態に係る高密度焦点式超音波システムの超音波トランスデューサに対する超音波送信器・受信器システムの相対位置を例示する図である。
【図１３】本発明の更なる一実施形態に係る高密度焦点式超音波システムの超音波トランスデューサに対する超音波送信器・受信器システムの相対位置を例示する図である。
【図１４】本発明の更なる一実施形態に係る高密度焦点式超音波システムの超音波トランスデューサに対する超音波送信器・受信器システムの相対位置を例示する図である。
【図１５】本発明の更なる一実施形態に係る高密度焦点式超音波システムの超音波トランスデューサに対する超音波送信器・受信器システムの相対位置を例示する図である。
【発明を実施するための形態】
【００６０】
これらの図において似通った参照符号を付された要素は、等価な要素であるか、同じ機能を実行するかの何れかである。先述した要素は、機能が等価である場合、以下の図においては必ずしも説明しない。
【００６１】
図１は、本発明の一実施形態に係る方法を例示するフロー図を示している。ステップ１００にて、医用撮像システムを用いて医用画像データが収集される。次にステップ１０２にて、超音波送信器及び受信器を用いて超音波データが収集される。ステップ１００及び１０２は、どのような順序で行われてもよく、互いに同時に行われてもよい。次にステップ１０４にて、医用画像データが少なくとも２つの組織型にセグメント化（分割）される。ステップ１０６にて、該少なくとも２つの組織型の各々内で進行された距離が決定される。これは、超音波送信器及び受信器によって生成された超音波が進行した少なくとも２つの経路の距離をトレースすることによって成し遂げられ得る。最後にステップ１０８にて、該少なくとも２つの組織型の各々内での超音波の速さが計算される。この計算は、超音波データと、少なくとも２つの経路にて該少なくとも２つの組織型の各々内で進行された距離とに従って行われる。超音波データは、送信器及び受信器による超音波の送信と受信との間の遅延についての情報を含み得る。
【００６２】
図２は、本発明の一実施形態に係る方法の更なる例を示すフロー図を示している。ステップ２００−２０８は図１のステップ１００−１０８に相当する。ステップ２００にて、医用撮像システムを用いて医用画像データが収集される。ステップ２０２にて、超音波送信器及び受信器を用いて超音波データが収集される。被検体を通る少なくとも２つの経路でデータが収集される。ステップ２０４にて、医用画像データが少なくとも２つの組織型にセグメント化される。被検体を通る経路の数は組織型の数以上にされるべきである。ステップ２０６にて、該少なくとも２つの組織型の各々内で進行された距離が決定される。ステップ２０８にて、該少なくとも２つの組織型の各々内での超音波の速さが計算される。ステップ２１０にて、一組（セット）のトランスデューサ素子の位相が計算される。ステップ２１２にて、このトランスデューサ素子位相セットを用いて超音波制御信号が生成される。ステップ２１４にて、この超音波制御信号が高密度焦点式超音波システムへと送信される。
【００６３】
図３は、人間の乳房内の相異なる組織型に関する超音波の波速の差を例示するプロットを示している。ｘ軸は温度スケールである。これは、摂氏温度単位で測定され、平常の中核温度に対しての温度である。ｙ軸３０２はメートル／秒単位で波速すなわち超音波速度を示している。線３０４は、胸脂肪内での超音波速度である。線３０６は腺組織内での超音波速度である。このプロットにて見て取れることには、胸脂肪での速度３０４と腺組織での速度３０６とは異なる。温度が上昇するにつれて、この差も大きくなっている。
【００６４】
超音波測定手段を解剖学的画像以外の画像種類と組み合わせることができること、及び複数の撮像モダリティを共に組み合わせることができることを説明するために、波速の温度依存性を考察することができる。
【００６５】
大抵の生体組織では、５０℃までの温度範囲で、超音波速度は温度上昇につれて約２．２ｍｓ−１℃−１で増大する。対照的に、間質組織内の脂質では、人間の胸脂肪の場合でおよそ−３．１ｍ／（ｓ℃）の傾きで、速度は温度上昇につれて低下する。
【００６６】
文献において、胸腺及び胸脂肪の波速は、如何なる加温もない場合に、それぞれ１５３８ｍ／ｓ及び１４５１ｍ／ｓであると報告されている。故に、図３は、文献に基づき、腺組織及び脂肪組織に対してそれぞれ２．２ｍｓ−１℃−１及び−３．１ｍｓ−１℃−１の一定傾きを仮定して、波速の温度上昇関数を与えている。脂肪組織と腺組織との間での波速の差は、＋１６．４℃の温度上昇に伴って８７ｍ／ｓから１７４ｍ／ｓまで２倍に増大する。
【００６７】
乳房組織内での音響収差を補償するように設計される位相補正は波速の関数であるので、温度の関数でもある。乳房組織内の波速の差は１６．４℃という典型的な加温の間に２倍にも増大するので、温度の調整のための温度分布の定量化が、加温なしでの当初の波速の定量化と同様に重要である。
【００６８】
本発明に係る方法の実施形態は、波速の定量化調整を行うために、乳房組織内の温度分布を測定することを有し得る。空間的な温度分布は、プロトン共鳴周波数シフトに基づく１つ以上の温度マッピング（特に、腺組織に対して）及び／又はＴ１、Ｔ２又はＴ２スター緩和変化に基づく温度マップ（特に、脂肪組織に対して）を用いて、やはりＭＲＩによって測定されることができる。
【００６９】
各組織型の波速の温度依存性の患者特定適応化のため、ソニケーション（超音波破砕）と超音波測定との間で素早く切り換えることによって、あるいは干渉が起こらない場合に双方を同時に行うことによって、超音波測定を治療手順と同時に実行することが可能であり得る。オンザフライで収集される温度マップ及び超音波測定結果を用いて、各組織型について超音波組織特性をリアルタイムで見積もることができる。例えば、各組織層の、同等の温度を有する一層小さい領域への細分化を用いて、以前の方程式を再利用することができる。結果として、発生する各温度変化に対して、小領域各々内の波速及び減衰が得られる。しかしながら、多数の小領域は、多数の超音波放射器及び受信器を必要とすることになる。このシステムを単純なままにするため、この方程式を解くための代替手法は、超音波測定結果に最も合致する組織型ごとの波速の適切な温度依存性の反復的な探索である。
【００７０】
図４は、本発明の一実施形態に係る医療装置４００を例示する図を示している。この実施形態において、医療装置４００は、磁気共鳴撮像システム４０２と高密度焦点式超音波システム４０４とを有している。
【００７１】
図５は、高密度焦点式超音波システム４０４の更に詳細な図を示している。この実施形態においては、磁気共鳴撮像システム４０２が使用されている。磁気共鳴撮像システム４０２は例示を意図したものである。他の実施形態においては、その他の医用撮像システムも使用され得る。例えば、磁気共鳴撮像システム４０２がコンピュータ断層撮影システム又は診断超音波システムで置き換えられてもよい。図４及び図５を一緒に説明する。
【００７２】
磁気共鳴撮像システム４０２は磁石４０６を有する。磁石４０６は、円筒形の超伝導型磁石であり、それを貫くボア４０８を有している。磁石４０６は、超伝導コイルを備えた液体ヘリウム冷却式のクライオスタットを有する。永久磁石又は常伝導磁石を用いることも可能である。異なる複数種類の磁石の使用も可能であり、例えば、分割円筒磁石と所謂オープンマグネットとの双方を使用することも可能である。分割円筒磁石は、磁石のアイソプレーンへのアクセスを可能にするようにクライオスタットが２つの区画に分割されていることを除いて、標準的な円筒磁石と同様である。そのような磁石は、例えば、荷電粒子ビーム治療とともに使用され得るものである。オープンマグネットは、被検体を受け入れるのに十分な大きさの空間を間に有するように互いに上下に配置される２つの磁石区画を有する。２つの区画領域のこの構成は、ヘルムホルツコイルの構成と同様である。オープンマグネットは、被検体の閉じ込め度合いが低いので好まれている。円筒磁石のクライオスタットの内部に、一群の超伝導コイルが存在する。円筒磁石のボア４０８内に、磁気共鳴撮像を行うのに十分なように磁場が強く且つ均一にされる撮像ゾーン４１８がある。
【００７３】
磁石のボア内にはまた、磁気共鳴データの収集に使用される傾斜磁場コイル４１０がある。傾斜磁場コイル４１０は、磁石４０６の撮像ゾーン４１８内の磁気スピンを空間エンコードする。傾斜磁場コイル４１０は傾斜磁場コイル電源４１２に接続されている。傾斜磁場コイル４１０は代表的な例を意図したものである。典型的に、傾斜磁場コイル４１０は、３つの直交する空間方向での空間エンコーディングのために、３つの別々のコイルセットを含む。傾斜磁場コイル電源４１２は、傾斜磁場コイルに電流を供給する。傾斜磁場コイル４１０に供給される電流は、時間の関数として制御され、傾斜（ランプ）あるいはパルス化され得る。
【００７４】
撮像ゾーン４１８に隣接して、撮像ゾーン４１８内の磁気スピンの向きを操作し且つやはり撮像ゾーン４１８内のスピンからの無線送信を受信する無線周波数コイル４１４がある。無線周波数アンテナは複数のコイル素子を含み得る。無線周波数コイルは、チャンネル、アンテナ、又は無線周波数アンテナとも呼ばれ得る。無線周波数コイル４１４は無線周波数トランシーバ４１６に接続されている。無線周波数コイル４１４及び無線周波数トランシーバ４１６は、別々の送信コイルと受信コイル、及び別々の送信器と受信器で置き換えられてもよい。理解されるように、無線周波数コイル４１４及び無線周波数トランシーバ４１６は代表的なものである。無線周波数アンテナは専用送信アンテナ及び専用受信アンテナをも表すことが意図される。同様に、トランシーバ４１６は、別々の送信器及び受信器をも表し得る。
【００７５】
磁石４０６のボア４０８内で、被検体支持台４２２上で被検体４２０が横になっている。被検体４２０の一部が撮像ゾーン４１８内にある。高密度焦点式超音波システム４０４はタンク４２４を有している。タンク４２４は超音波伝導流体４２６で充たされている。流体の頂面すなわち液面がマーカー４２８によって指し示されている。被検体４２０の部分４３０が、被検体支持台４２２の開口４２３を通り抜けてタンク４２４内にある。部分４３０は延在して、超音波伝導流体４２６によって部分的に取り囲まれている。この例において、２つの超音波トランスデューサ素子４３２が流体４２６内にあるものとして図示されている。被検体４２０の部分４３０は乳房である。脂肪組織４３４の領域及び腺組織４３６の領域である２つの組織領域が存在している。破線４３８は、超音波トランスデューサ素子４３２から、流体４２６を通り、脂肪組織４３４を通り、且つ腺組織４３６を通って、標的ボリューム４４０内に位置するソニケーションボリューム４４２までの、超音波経路４３８を描いている。
【００７６】
超音波トランスデューサ素子の上方に、超音波送信器４４４及び超音波受信器４４６がある。超音波送信器４４４から超音波受信器４４６まで超音波ビーム４４８が進行する。鎖線によって指し示された第１の平面４５０及び第２の平面４５２が存在している。超音波送信器４４４及び超音波受信器４４６は主として第１の平面４５０内に位置している。超音波は主として第１の平面を通って進行する。第２の平面４５２は、第１の平面４５０の下方に位置している。超音波トランスデューサ素子４３２は第２の平面４５２の下方に位置するが、ソニケーションボリューム４４２は第２の平面４５２の上方に位置している。この実施形態は、ソニケーションボリューム４４２の近傍での乳房組織内超音波速度の直接測定を可能にするので有利である。実際には、一部の実施形態において、ソニケーションボリューム４４２も主として第１の平面４５０内に位置する。平面４５０は当然ながら２次元空間である。超音波ビーム４４８及びソニケーションボリューム４４２は３次元ボリュームである。本明細書では、平面がボリュームの一部を横断する場合、そのボリュームは主としてその平面内にあるとする。
【００７７】
高密度焦点式超音波システム４０４、傾斜磁場コイル電源４１２、及びトランシーバ４１６は、コンピュータシステム４５４のハードウェアインタフェース４５６に接続されるものとして図示されている。ハードウェアインタフェース４５６はプロセッサ４５８に接続されている。プロセッサ４５８は更に、ユーザインタフェース４６０、コンピュータストレージ４６２、及びコンピュータメモリ４６４に接続されている。
【００７８】
コンピュータストレージ４６２は、磁気共鳴撮像システム４０２によって収集された磁気共鳴データ４６６を含むように図示されている。磁気共鳴データ４６６から再構成された磁気共鳴画像データ４６８が存在している。磁気共鳴画像データ４６８は医用画像データの一例である。コンピュータストレージ４６２内にはまた、磁気共鳴画像データ４６８の画像セグメンテーション４７０が格納されている。画像セグメンテーション４７０は、少なくとも２つの組織型４３４、４３６の位置を記述する情報を含む。コンピュータストレージ内に温度マップ４７２が存在している。温度マップは、磁気共鳴データ４６６から再構成される。磁気共鳴データは、磁気共鳴サーモメトリデータを有する。コンピュータストレージはまた、距離データ４７４を含むものとして図示されている。距離データは、超音波ビーム４４８が第１及び第２の組織型４３４、４３６を通って進行した距離を記述する。第１の組織型は脂肪組織４３４であり、第２の組織型は腺組織４３６である。コンピュータストレージは、超音波データ４７６を含むものとして図示されている。超音波データは、超音波受信器４４６によって収集されたデータである。
【００７９】
コンピュータストレージは、超音波データ４７６から抽出されたパルス遅延データ４７８及び振幅データ４８０を含むものとして図示されている。パルス遅延データ４７８は、送信器４４４から受信器４４６まで進行する超音波ビーム４４８における遅延を記述する。振幅データ４８０は、送信器４４４から受信器４４６までの振幅の変化を記述する。コンピュータストレージは更に、超音波制御信号４８４を含むものとして図示されている。超音波制御信号は、高密度焦点式超音波システム４０４を動作させるための命令を含む。コンピュータストレージは、パルス遅延データ４７８と距離データ４７４とを用いて計算された超音波速度４８６を含むものとして図示されている。距離データ４７４及び振幅データ４８０は、２つの組織型４３４、４３６内での超音波減衰を計算するために使用され得る。
【００８０】
コンピュータメモリ４６４は、制御モジュール４９０を含むものとして図示されている。制御モジュール４９０は、医療装置４００の動作及び機能を制御するためのコンピュータ実行可能コードを含む。コンピュータメモリ４６４は更に、磁気共鳴データ４６６を収集するために磁気共鳴撮像システム４０２によって実行される操作を含むパルスシーケンス４９２を含むものとして図示されている。パルスシーケンス４９２を用いて磁気共鳴撮像システム４０２を制御する磁気共鳴撮像システム制御モジュール４９４が存在している。コンピュータメモリ４６４は更に、画像モジュール４９６を含むものとして図示されている。画像モジュール４９６は、磁気共鳴データ４６６を磁気共鳴画像データ４６８へと再構成するため、且つ／或いは磁気共鳴画像データ４６８からの画像セグメンテーション４７０の実行及び生成を行うためのコンピュータ実行可能コードを含む。画像モジュール４９６はまた、磁気共鳴データ４６６から温度マップ４７２を生成するためのコンピュータ実行可能コードを含み得る。
【００８１】
コンピュータメモリ４６４は更に、高密度焦点式超音波制御モジュール４９８を含むものとして図示されている。高密度焦点式超音波制御モジュール４９８は、高密度焦点式超音波システム４０４の動作及び機能を制御するためのコンピュータ実行可能コードを含む。高密度焦点式超音波制御モジュール４９８は、超音波制御信号４８４を生成するように適応される。コンピュータメモリ４６４は、画像分析モジュール５００を更に含むものとして図示されている。画像分析モジュールは、画像セグメンテーション４７０と磁気共鳴画像データ４６８とを用いて距離データ４７４を計算することができる。
【００８２】
コンピュータメモリ４６４は更に、速度計算モジュール５０２及び減衰計算モジュール５０４を含むものとして図示されている。速度計算モジュールは、距離データ４７４とパルス遅延データ４７８とから超音波速度４８６を計算するように適応されたコンピュータ実行可能コードを含む。減衰計算モジュール５０４は、振幅データ４８０と距離データ４７４とを用いて超音波減衰を計算するためのコンピュータ実行可能コードを含む。
【００８３】
図６は、本発明の一実施形態に係る超音波送信器・受信器システム６００の一例を示している。この図には４つの送信器６０２が示されている。各送信器６０２の真向かいに受信器６０４がある。破線６０６は送信器６０２から受信器６０４までの経路を描いている。経路６０６は、超音波伝導流体６０８、第１の組織型６１０及び第２の組織型６１２を通って進行している。第１及び第２の組織型６１０、６１２内での速度を決定するため、医用撮像システムは、経路６０６が位置するのと同じ平面内で医用画像データを収集し得る。各経路に対し、超音波伝導流体６０８内、第１の組織型６１０内及び第２の組織型６１２内で進行された距離を計算することができる。送信器６０２から受信器６０４への伝送間の遅延を決定することにより、超音波伝導流体内及び組織型６１０、６１２の各々内での速度が計算され得る。
【００８４】
ＭＲ−ＨＩＦＵ治療に先立ち、患者は、診断を確立するための様々な種類の医療検査を受ける。このスクリーニング段階は、その他の治療種類に対してのＨＩＦＵ治療の成功を推測するために必須である。超音波組織特性を特徴付ける手段をスクリーニング用撮像モダリティ（超音波の波速測定以外）と組み合わせて使用することにより、治療オプションが向上され得る。
【００８５】
乳癌のＨＩＦＵ治療では、結果として生じる位相収差は限られた範囲の値においてのみ補正され得るので、組織の波速及び形状の知識が不可欠である。故に、上述の方法を用いて得られる組織モデルに基づく音響場のシミュレーションにより、ＨＩＦＵによって患者を治療する可能性への答えが提供され得る。
【００８６】
このようなスクリーニング過程において、上述の方法と同じ方法を用いて各組織型での超音波の波速及び減衰を決定するため、超音波測定手段が、１つ以上の撮像モダリティ（超音波の波速測定以外）をＨＩＦＵシステムとは独立に組み合わせるシステム内に統合され得る。
【００８７】
例えば、図６にて説明したリング構造が、乳癌を検出するために通常使用されるＭＲＩマンモグラフィトラックテーブルに挿入され得る。ＭＲＩマンモグラフィトラックテーブルは、図４及び５に示した患者支持台４４２と同等のものである。この場合、超音波放射器及び受信器は、ＭＲコイル受信器ループ間に音響波を伝搬させるように配置され得る。
【００８８】
患者スクリーニングステップにおいて各組織型での超音波の波速及び減衰が定量化され、且つ患者がＨＩＦＵ治療に適格であると、得られた組織の特徴付けを治療手順に再利用することができる。各組織での超音波の波速及び減衰は、スクリーニング日と治療日とでほぼ確実に同じであるので、ビームをリフォーカシングするために適用すべき位相及び強度の補正を処理することには、セグメント化された組織層のレジストレーション又は組織層の新たなセグメンテーションのみが必要である。当然ながら、必要に応じて、超音波の波速及び減衰の測定を、これらの測定結果の信頼性を確実にすることが望まれる場合に、治療日に再収集することができる。しかしながら、これらの測定が時間のかかるものであり且つスクリーニングステップにて既に収集されている場合、総治療時間を短縮するために、治療日における超音波の波速及び減衰の定量化は省略され得る。
【００８９】
磁気共鳴−高密度焦点式超音波（ＭＲ−ＨＩＦＵ）プラットフォームは、ＭＲＩによってセグメント化された組織の波速を定量化すること担う超音波放射器／受信器システムを含むことができる。図６に示した例において、このシステムは、４つの方向に沿って乳房組織を特徴付けるために、剛体リングに取り付けられた、４つの受信器の正面に配置された４つの放射器からなっている。
【００９０】
超音波パルスの放射及び受信は、各方向ｉに沿った超音波進行時間Ｔｉ及び圧力振幅Ｐｉを提供する。同様に、先行するシステム校正ステップにて、胸部とトランスデューサとの間の結合媒体に従って水のみ又はオイルのみを用いて、ＴＲｅｆｉ及びＰＲｅｆｉの定量化も行うことができる。１５ｃｍの距離にわたる超音波進行時間は約０．１ｍｓであるので、１つの直線に沿った測定時間は非常に高速である。故に、治療手順時間の有意な遅れを引き起こすことなく、平均化目的で複数回の測定を繰り返すことができる。このような平均化は、遅延及び圧力振幅の測定の精度を向上させる。
【００９１】
ＭＲ画像のセグメンテーションアルゴリズムにより、脂肪組織Ｆ、腺組織Ｇ、及び水Ｗの描写が提供される。放射器から受信器まで各媒体の表面を横切って走る直線の交点が、長さＬＷｉ、ＬＦｉ、及びＬＧｉの定量化を提供する。しかしながら、長さＬＷｉの定量化は必須ではない。
【００９２】
進行時間は、水媒体、脂肪媒体及び腺媒体の波速ＣＷ、ＣＦ、ＣＧ及び減衰ＡＷ、ＡＦ、ＡＧによって、線形方程式：
【００９３】
【数１】
2014517704.tif 000003
に従って互いに結び付けられる。
【００９４】
双方の方程式の減算によって長さＬＷｉが相殺され、それにより、乳房組織での波速の定量化がより一層正確になる。ＭＲＩによる長さＬＷｉの測定は、磁化率の場の鋭い変化に関係するトランスデューサの近傍での画像歪みのために不正確である。ＴＲｅｆｉとＴｉとの間の伝搬時間の差：
【００９５】
【数２】
2014517704.tif 000004
は、２つの同等の波形の比較に相当し、これは相互相関によって非常に正確に評価されることができる。
【００９６】
結果として、波速ＣＦ及びＣＧは、先の方程式の反転によって得られる。
【００９７】
【数３】
2014517704.tif 000005
同様に、圧力振幅は、水媒体、脂肪媒体及び腺媒体の減衰αＷ、αＦ、αＧ及び放射圧力振幅Ｐ０によって、線形方程式：
【００９８】
【数４】
2014517704.tif 000006
に従って互いに結び付けられる。
【００９９】
双方の方程式の減算によって長さＬＷｉ及び圧力振幅Ｐ０が相殺され、それにより、組織での減衰の定量化がより一層正確になる。
【０１００】
【数５】
2014517704.tif 000007
結果として、減衰αＦ及びαＧは、先の方程式の反転によって得られる。
【０１０１】
【数６】
2014517704.tif 000008
この方法は、組織型の数が測定方向の数より少ないか等しいかである限り、３つ以上の組織型を特徴付けるように一般化されることができる。しかしながら、組織特徴の良好な測定精度は、２つの組織型と４つの方向とを有するこの例のように、方向の数が組織型の数より多いときに得られる。
【０１０２】
図７は、本発明の一実施形態に係る超音波送信器・受信器システム７００を示している。この例においては、単一の送信器７０２と複数の受信器７０４、超音波伝導流体７０６の領域、第１の組織型７０８、第２の組織型７１０、及び骨組織７１２が存在している。超音波は、送信器７０２から、経路７１４に沿って、超音波伝導流体、第１の組織型７０８、第２の組織型７１０を通って進行する。そして、超音波は骨７１２によって反射される。破線７１６は、受信器７０４へと戻る超音波の経路を示している。経路７１４及び７１６の平面内で収集される医用画像データを用いて、超音波伝導流体内、第１の組織型７０８内及び第２の組織型７１０内で進行される距離を計算し得る。これを、超音波の送信７０２と異なる受信器７０４による受信との間の遅延とともに用いて、これらの組織型内での超音波の速度を計算し得る。
【０１０３】
骨での超音波の反射を用いてその他の人体領域内の皮下脂肪及び筋肉を特徴付ける代替方法を考察し得る。骨転移並びに肝臓癌及び腎臓癌の治療において頻繁に骨構造に遭遇する。
【０１０４】
図８、９及び１０は、超音波トランスデューサ素子８０２の構成の等角図８００、上面図９００及び側面図１０００を示している。小さい円の各々が１つの超音波トランスデューサ素子である。破線８０４は、ソニケーションボリューム８０８までの超音波の経路を示している。参考までに、複数の形状８０６は、３つの異なる図８００、９００及び１０００で同じ組織構造を示している。
【０１０５】
図１１−１５は、高密度焦点式超音波システムの超音波トランスデューサに対する超音波送信器・受信器システムの様々な構成を示している。
【０１０６】
図１１において、領域１１００は被検体を表している。２つの超音波トランスデューサ素子が示されている。矢印１１０４及びその周囲の錐体状領域は、被検体１１００内のソニケーションボリューム１１０６に集束される治療超音波１１０４の経路を示している。ソニケーションボリューム１１０６を囲んで加熱領域１１０８が存在している。この実施形態において、超音波トランスデューサ素子１１０２は上向きに集束されている。超音波トランスデューサ素子１１０２の上方に位置する平面内に、被検体を通り抜ける経路１１１４に沿って超音波パルスを伝える送信器１１１０及び受信器１１１２が存在している。この実施形態においては、被検体１１００と送信器１１１０との間に空間１１１６が存在し、被検体１１００と受信器１１１２との間にも空間１１１６が存在している。
【０１０７】
図１２は、この実施形態においては空間が存在しないことを除いて同様の実施形態を示している。送信器１１１０は被検体１１００と直に接触している。受信器１１１２も被検体１１００と直に接触している。図１１の実施形態においては、空間１１１６は例えばゲルパッド又はその他の同様の超音波伝導媒体で充填され得る。
【０１０８】
図１３に示す実施形態において、超音波トランスデューサ素子１１０２はなおも、図１１及び１２においてと同じ関係を有している。しかしながら、この実施形態においては、別個の送信器１１１０は存在していない。超音波トランスデューサ素子１１０２が、図１１及び１２に示した送信器１１１０として機能する。２つの受信器１１１２が、この場合には治療超音波１１０４の経路に沿って位置して存在している。
【０１０９】
図１４には、異なる構成の超音波トランスデューサ素子が示されている。この実施形態において、超音波トランスデューサ素子１１０２は、被検体１１００の側方にあるものとして示されている。治療超音波１１０４の経路はここでも、被検体１１００内のソニケーションボリューム１１０６に集束されている。加熱領域１１０８がソニケーションボリューム１１０６を取り囲んでいる。治療超音波１１０４は、領域１１０８を立ち去って受信器１１１２まで進行する。図１３に示した実施形態と同様に、送信器は存在していない。超音波トランスデューサ素子１１０２が送信器１１１０として機能する。
【０１１０】
図１５に示す実施形態は、図１４の実施形態と同様である。この実施形態においては、２つの超音波トランスデューサ素子１１０２が存在している。これらは被検体１１００の隣り合う側方にある。治療超音波１１０４の経路は、双方の超音波トランスデューサ素子１１０２によってソニケーションボリューム１１０６に集束される。この実施形態においても、送信器は超音波トランスデューサ素子１１０２によって置き換えられている。受信器１１１２は各超音波トランスデューサ素子１１０２上に位置している。被検体を通り抜ける経路１１１４は治療超音波１１０４の経路と相等しい。
【０１１１】
図面及び以上の記載にて本発明を詳細に図示して説明したが、これらの図示及び説明は、限定的なものではなく、例示的あるいは典型的なものと見なされるべきであり、本発明は開示した実施形態に限定されるものではない。
【０１１２】
図面、明細書及び特許請求の範囲の学習から、請求項に記載の発明を実施使用とする当業者によって、開示した実施形態へのその他の変形が理解・達成され得る。請求項において、用語“有する”はその他の要素又はステップを排除するものではなく、不定冠詞“ａ”若しくは“ａｎ”は複数であることを排除するものではない。単一のプロセッサ又はその他のユニットが、請求項に記載された複数のアイテムの機能を果たしてもよい。特定の複数の手段が相互に異なる従属請求項に記載されているという単なる事実は、それらの手段の組合せが有利に使用され得ない、ということを指し示すものではない。コンピュータプログラムは、他のハードウェアとともに供給されるか、あるいは他のハードウェアの部分として供給されるかする例えば光記憶媒体又は半導体媒体などの好適な媒体上で格納／配信され得るが、例えばインターネット又はその他の有線若しくは無線の遠隔通信システムを介してなど、その他の携帯で配信されてもよい。請求項中の参照符号は、範囲を限定するものとして解されるべきでない。
【符号の説明】
【０１１３】
４００  医療装置
４０２  磁気共鳴撮像システム
４０４  高密度焦点式超音波システム
４０６  磁石
４０８  ボア
４１０  傾斜磁場コイル
４１２  傾斜磁場コイル電源
４１４  無線周波数コイル
４１６  トランシーバ
４１８  撮像ゾーン
４２０  被検体
４２２  被検体支持台
４２３  開口
４２４  タンク
４２６  流体
４２８  液面
４３０  被検体の部分
４３２  超音波トランスデューサ素子
４３４  脂肪組織
４３６  腺組織
４３８  超音波の経路
４４０  標的ゾーン
４４２  ソニケーションボリューム
４４４  超音波送信器
４４６  超音波受信器
４４８  超音波ビーム
４５０  第１の平面
４５２  第２の平面
４５４  コンピュータシステム
４５６  ハードウェアインタフェース
４５８  プロセッサ
４６０  ユーザインタフェース
４６２  コンピュータストレージ
４６４  コンピュータメモリ
４６６  磁気共鳴データ
４６８  磁気共鳴画像データ
４７０  画像セグメンテーション
４７２  温度マップ
４７４  距離データ
４７６  超音波データ
４７８  パルス遅延データ
４８０  振幅データ
４８２  トランスデューサ素子位相
４８４  超音波制御信号
４８６  超音波速度
４９０  制御モジュール
４９２  パルスシーケンス
４９４  磁気共鳴撮像システム制御モジュール
４９６  画像モジュール
４９８  高密度焦点式超音波制御モジュール
５００  画像分析モジュール
５０２  速度計算モジュール
５０４  減衰計算モジュール
６００  超音波送信器・受信器システム
６０２  送信器
６０４  受信器
６０６  経路
６０８  超音波伝導流体
６１０  第１の組織型
６１２  第２の組織型
７００  超音波送信器・受信器システム
７０２  送信器
７０４  受信器
７０６  超音波伝導流体
７０８  第１の組織型
７１０  第２の組織型
７１２  骨組織
７１４  超音波の経路
７１６  超音波の経路
８００  等角図
８０２  超音波トランスデューサ素子
８０４  超音波の経路
８０６  組織構造
８０８  ソニケーションボリューム
９００  上面図
１０００  側面図
１１００  被検体
１１０２  超音波トランスデューサ素子
１１０４  治療超音波の経路
１１０６  ソニケーションボリューム
１１０８  加熱領域
１１１０  送信器
１１１２  受信器
１１１４  被検体貫通経路
【図１】
2014517704.tif 000009
【図２】
2014517704.tif 000010
【図３】
2014517704.tif 000011
【図４】
2014517704.tif 000012
【図５】
2014517704.tif 000013
【図６】
2014517704.tif 000014
【図７】
2014517704.tif 000015
【図８】
2014517704.tif 000016
【図９】
2014517704.tif 000017
【図１０】
2014517704.tif 000018
【図１１】
2014517704.tif 000019
【図１２】
2014517704.tif 000020
【図１３】
2014517704.tif 000021
【図１４】
2014517704.tif 000022
【図１５】
2014517704.tif 000023
【国際調査報告】
2014517704.tif 000024
2014517704.tif 000025
2014517704.tif 000026