(57)【要約】本発明の態様は、第３イメージング・モダリティ（例えば、超音波システム）の使用を介して得られる２つのイメージング・システム（例えば、ＣＴシステムおよびＭＲＩシステム）の画像のレジストレーションを検証および／または取得するシステムおよび方法を供する。【選択図】図１Ａ〜Ｃ
2018534038.tif 000002
【特許請求の範囲】
【請求項１】
第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像のレジストレーションを検証するための方法であって、
（ａ）第１イメージング・システムを用いて、解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得ること、
（ｂ）第３イメージング・システムを用いて、第３イメージング・システムと解剖学的対象との間の距離を測定すること、
（ｃ）第２イメージング・システムを用いて、解剖学的対象の第２画像を得ること、
（ｄ）前記レジストレーションを用いて、第１画像および第２画像を位置合せすること、ならびに
（ｅ）（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）前記測定距離に基づいて、前記レジストレーションの誤差を算出することを含んで成る、方法。
【請求項２】
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムが、磁気共鳴イメージング（ＭＲＩ）システム、コンピュータ・トポグラフィ（ＣＴ）システムおよび超音波トランスデューサ・システムをそれぞれ含んで成る、請求項１に記載の方法。
【請求項３】
第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を得ることをさらに含んで成る、請求項１に記載の方法。
【請求項４】
第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置が、飛行時間法に基づいて決定される、請求項３に記載の方法。
【請求項５】
第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出することをさらに含んで成る、請求項３に記載の方法。
【請求項６】
第３イメージング・システムの位置を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換することをさらに含んで成る、請求項５に記載の方法。
【請求項７】
解剖学的対象の第２画像を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換することをさらに含んで成る、請求項６に記載の方法。
【請求項８】
第３イメージング・システムの変換された位置および第１イメージング・システムの座標系における変換された第２画像に基づいて、第３イメージング・システムと解剖学的対象との間の距離を算出することをさらに含んで成る、請求項７に記載の方法。
【請求項９】
ステップ（ｂ）で得られる測定距離を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換することをさらに含んで成る、請求項８に記載の方法。
【請求項１０】
前記レジストレーションの誤差が、前記算出距離から変換された測定距離のずれに基づいて決定される、請求項９に記載の方法。
【請求項１１】
前記レジストレーションの誤差と予め決められた閾値とを比較すること、および該比較に基づいて、前記レジストレーションの妥当性を決定することをさらに含んで成る、請求項１０に記載の方法。
【請求項１２】
第３イメージング・システムから送信され、また第３イメージング・システムによって受信される信号に基づいて、第３イメージング・システムと解剖学的対象との間の距離が測定される、請求項１に記載の方法。
【請求項１３】
第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像のレジストレーションを検証するためのシステムであって、
解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得るための第１イメージング・システム、
解剖学的対象の第２画像を得るための第２イメージング・システム、ならびに
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムと連通するコントローラを含んで成り、
コントローラが、
第３イメージング・システムと解剖学的対象との間の距離を測定するようになっており、
第１画像および第２画像を位置合せするようになっており、ならびに
（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）前記測定距離に基づいて、前記レジストレーションにおける誤差を算出するようになっている、システム。
【請求項１４】
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムが、磁気共鳴イメージング（ＭＲＩ）システム、コンピュータ・トポグラフィ（ＣＴ）システムおよび超音波トランスデューサ・システムをそれぞれ含んで成る、請求項１３に記載のシステム。
【請求項１５】
コントローラが、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を決定するようにさらに構成されている、請求項１３に記載のシステム。
【請求項１６】
コントローラが、飛行時間法に基づいて、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を決定するようにさらに構成されている、請求項１５に記載のシステム。
【請求項１７】
コントローラが、第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システム座標系を関連付ける変換を算出するようにさらに構成されている、請求項１５に記載のシステム。
【請求項１８】
コントローラが、第３イメージング・システムの位置を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようにさらに構成されている、請求項１７に記載のシステム。
【請求項１９】
コントローラが、解剖学的対象の第２画像を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようにさらに構成されている、請求項１８に記載のシステム。
【請求項２０】
コントローラが、第１イメージング・システムの座標系における第３イメージング・システムの変換された位置および変換された第２画像に基づいて、第３イメージング・システムと解剖学的対象との間の距離を算出するようにさらに構成されている、請求項１９に記載のシステム。
【請求項２１】
コントローラが、ステップ（ｂ）で得られる測定距離を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようにさらに構成されている、請求項２０に記載のシステム。
【請求項２２】
コントローラが、前記算出距離から変換された測定距離のずれに基づいて、前記レジストレーションの誤差を決定するようにさらに構成されている、請求項２１に記載のシステム。
【請求項２３】
コントローラが、前記レジストレーションの誤差と予め決められた閾値とを比較し、該比較に基づいて、該レジストレーションの妥当性を決定するようにさらに構成されている、請求項２２に記載のシステム。
【請求項２４】
コントローラが、第３イメージング・システムから送信され、また第３イメージング・システムによって受信される信号に基づいて、第３イメージング・システムと解剖学的対象との間の距離を測定するようにさらに構成されている、請求項１３に記載のシステム。
【請求項２５】
第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像を位置合せするための方法であって、
（ａ）第１イメージング・システムを用いて、解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得ること、
（ｂ）第２イメージング・システムを用いて、解剖学的対象の第２画像を得ること、
（ｃ）第３イメージング・システムを用いて、解剖学的対象の第３画像を得ること、
（ｄ）第２画像および第３画像を位置合せすること、
（ｅ）（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）位置合せされた第２画像および第３画像に基づいて、第１イメージング・システムの座標系、第２イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出することを含んで成る、方法。
【請求項２６】
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムが、磁気共鳴イメージング（ＭＲＩ）システム、コンピュータ・トポグラフィ（ＣＴ）システムおよび超音波トランスデューサ・システムをそれぞれ含んで成る、請求項２５に記載の方法。
【請求項２７】
第２画像および第３画像のレジストレーションが、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を第３イメージング・システムの座標系における座標に変換することを含んで成る、請求項２５に記載の方法。
【請求項２８】
第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を第３イメージング・システムの座標系における内側解剖学的対象に関連する座標にフィッティングすることをさらに含んで成る、請求項２７に記載の方法。
【請求項２９】
第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を得ることをさらに含んで成る、請求項２５に記載の方法。
【請求項３０】
第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置は、飛行時間法に基づいて決定される、請求項２９に記載の方法。
【請求項３１】
第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出することをさらに含んで成る、請求項２９に記載の方法。
【請求項３２】
第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換が、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換に基づいて算出される、請求項２５に記載の方法。
【請求項３３】
第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像を位置合せするためのシステムであって、
解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得るための第１イメージング・システム、
解剖学的対象の第２画像を得るための第２イメージング・システム、
解剖学的対象の第３画像を得るための第３イメージング・システム、ならびに
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムと連通するコントローラを含んで成り、
コントローラが、
第２画像および第３画像を位置合せするようになっており、
（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）位置合せされた第２画像および第３画像に基づいて、第１イメージング・システムの座標系、第２イメージング・システムのシステムおよび第３イメージング・システムの座標系を関連付ける変換を算出するようになっている、システム。
【請求項３４】
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムが、磁気共鳴イメージング（ＭＲＩ）システム、コンピュータ・トポグラフィ（ＣＴ）システムおよび超音波トランスデューサ・システムをそれぞれ含んで成る、請求項３３に記載のシステム。
【請求項３５】
コントローラが、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を第３イメージング・システムの座標系における座標に変換するようにさらに構成されている、請求項３３に記載のシステム。
【請求項３６】
コントローラが、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を、第３イメージング・システムの座標系における内側解剖学的対象に関連する座標にフィッティングするようにさらに構成されている、請求項３５に記載のシステム。
【請求項３７】
コントローラが、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を決定するようにさらに構成されている、請求項３３に記載のシステム。
【請求項３８】
コントローラが、飛行時間法に基づいて、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を決定するようにさらに構成されている、請求項３７に記載のシステム。
【請求項３９】
コントローラが、第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出するようになっている、請求項３７に記載のシステム。
【請求項４０】
コントローラが、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出するようにさらに構成されている、請求項３３に記載のシステム。
【請求項４１】
処置シーケンスの間にて動く解剖学的特徴を検出するための方法であって、
（ａ）処置シーケンスの前に、（ｉ）解剖学的特徴および第１イメージング・システムの少なくとも一部分の第１画像を得ること、ならびに（ｉｉ）解剖学的特徴と第１イメージング・システムの少なくとも一部分との間の最短距離を算出するように第１画像を処理すること、ならびに
（ｂ）処置シーケンスの間にて、（ｉ）解剖学的特徴と第１イメージング・システムの少なくとも一部分との間の最短距離を測定すること、（ｉｉ）測定された最短距離とステップ（ａ）で得られる算出された最短距離とを、それらからのずれを決定するように比較すること、ならびに（ｉｉｉ）該ずれに基づいて、解剖学的特徴の動きを決定することを含んで成る、方法。
【請求項４２】
第１画像が、第２イメージング・システムを用いて得られ、前記方法が、処置シーケンスの前に、第３イメージング・システムを用いて解剖学的特徴の第２画像を得ることをさらに含んで成る、請求項４１に記載の方法。
【請求項４３】
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムが、超音波トランスデューサ・システム、磁気共鳴イメージング（ＭＲＩ）システムおよびコンピュータ・トポグラフィ（ＣＴ）システムをそれぞれ含んで成る、請求項４２に記載の方法。
【請求項４４】
第１イメージング・システムの座標系における第１イメージング・システムの少なくとも一部分の位置を得ることをさらに含んで成る、請求項４３に記載の方法。
【請求項４５】
第１イメージング・システムの少なくとも一部分の位置が、飛行時間法に基づいて決定される、請求項４４に記載の方法。
【請求項４６】
第１画像および第１イメージング・システムの座標系における第１イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第２イメージング・システムの座標系を関連付ける変換を算出することをさらに含んで成る、請求項４４に記載の方法。
【請求項４７】
算出された変換に基づいて、第１イメージング・システムの位置を、第１イメージング・システムの座標系から第２イメージング・システムの座標系に変換することをさらに含んで成る、請求項４６に記載の方法。
【請求項４８】
第１画像および第２画像を位置合せすることをさらに含んで成る、請求項４７に記載の方法。
【請求項４９】
解剖学的特徴の第２画像を、第３イメージング・システムの座標系から第２イメージング・システムの座標系に変換することをさらに含んで成る、請求項４８に記載の方法。
【請求項５０】
前記最短距離が、第１イメージング・システムの変換された位置および第２イメージング・システムの座標系における変換された第２画像に基づいて算出される、請求項４９に記載の方法。
【請求項５１】
解剖学的特徴と第１イメージング・システムの少なくとも一部分との間の前記最短距離が、第１イメージング・システムから送信され、また第１イメージング・システムによって受信される信号に基づいて測定される、請求項５０に記載の方法。
【請求項５２】
前記測定された最短距離を、第１イメージング・システムの座標系から第２イメージング・システムの座標系に変換することをさらに含んで成る、請求項５１に記載の方法。
【請求項５３】
前記ずれと予め決められた閾値とを比較すること、および該比較に基づいて、解剖学的特徴の動きを決定することをさらに含んで成る、請求項５２に記載の方法。
【請求項５４】
処置シーケンスの間にて動く解剖学的特徴を検出するためのシステムであって、
処置シーケンスの前に、解剖学的特徴および第２イメージング・システムの少なくとも一部分の第１画像を得るための第１イメージング・システム、ならびに
第１および第２システムと連通するコントローラを含んで成り、
コントローラが、
（ａ）処置シーケンスの前に、解剖学的特徴と第２イメージング・システムの少なくとも一部分との間の最短距離を算出するように第１画像を処理するようになっており、ならびに
（ｂ）処置シーケンスの間にて、（ｉ）解剖学的特徴と第２イメージング・システムの少なくとも一部分との間の最短距離を測定するようになっており、（ｉｉ）測定された最短距離とステップ（ａ）で得られる算出された最短距離とを、それらからのずれを決定するように比較するようになっており、ならびに（ｉｉｉ）該ずれに基づいて解剖学的特徴の動きを決定するようになっている、システム。
【請求項５５】
処置シーケンスの前に解剖学的特徴の第２画像を得るための第３イメージング・システムをさらに含んで成る、請求項５４に記載のシステム。
【請求項５６】
第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムが、磁気共鳴イメージング（ＭＲＩ）システム、超音波トランスデューサ・システムおよびコンピュータ・トポグラフィ（ＣＴ）システムをそれぞれ含んで成る、請求項５５に記載のシステム。
【請求項５７】
コントローラが、第２イメージング・システムの座標系における第２イメージング・システムの少なくとも一部分の位置を得るようにさらに構成されている、請求項５６に記載のシステム。
【請求項５８】
コントローラが、飛行時間法に基づいて、第２イメージング・システムの少なくとも一部分の位置を決定するようにさらに構成されている、請求項５７に記載のシステム。
【請求項５９】
コントローラが、第１画像および第２イメージング・システムの座標系における第２イメージング・システムの少なくとも一部分の位置に基づいて、第２イメージング・システムの座標系および第１イメージング・システムの座標系を関連付ける変換を算出するようにさらに構成されている、請求項５７に記載のシステム。
【請求項６０】
コントローラが、算出された変換に基づいて、第２イメージング・システムの位置を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようにさらに構成されている、請求項５９に記載のシステム。
【請求項６１】
コントローラが、第１画像および第２画像を位置合せするようにさらに構成されている、請求項６０に記載のシステム。
【請求項６２】
コントローラが、解剖学的特徴の第２画像を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようにさらに構成されている、請求項６１に記載のシステム。
【請求項６３】
コントローラが、第２イメージング・システムの変換された位置および第１イメージング・システムの座標系における変換された第２画像に基づいて、前記最短距離を算出するようにさらに構成されている、請求項６２に記載のシステム。
【請求項６４】
コントローラが、第２イメージング・システムから送信され、また第２イメージング・システムによって受信される信号に基づいて、前記最短距離を測定するようにさらに構成されている、請求項６３に記載のシステム。
【請求項６５】
コントローラが、測定された最短距離を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようにさらに構成されている、請求項６４に記載のシステム。
【請求項６６】
コントローラが、前記ずれと予め決められた閾値とを比較し、該比較に基づいて、解剖学的特徴の動きを決定するようにさらに構成されている、請求項６５に記載のシステム。
【請求項６７】
複数の処置期間を有する処置シーケンスの間にて動く解剖学的特徴を検出するための方法であって、
処置シーケンスの間にて、解剖学的特徴とイメージング・システムの少なくとも一部分との間の最短距離を測定すること、
現在の処置期間における測定された最短距離と前の処置期間における測定された最短距離とを、それらからのずれを決定するように比較すること、ならびに
前記ずれに基づいて、解剖学的特徴の動きを決定することを含んで成る、方法。
【請求項６８】
解剖学的特徴と前記イメージング・システムの少なくとも一部分との間の前記最短距離が、該イメージング・システムから送信され、また該イメージング・システムによって受信される信号に基づいて測定される、請求項６７に記載の方法。
【請求項６９】
前記ずれと予め決められた閾値とを比較すること、および該比較に基づいて、解剖学的特徴の動きを決定することをさらに含んで成る、請求項６８に記載の方法。
【請求項７０】
前記イメージング・システムが、超音波トランスデューサ・システムを含んで成る、請求項６７に記載の方法。
【請求項７１】
複数の処置期間を有する処置シーケンスの間にて動く解剖学的特徴を検出するためのシステムであって、
イメージング・システムであって、処置シーケンスの間にて、解剖学的特徴と該イメージング・システムの少なくとも一部分との間の最短距離を測定するためのイメージング・システム、ならびに
前記イメージング・システムと連通するコントローラを含んで成り、
コントローラが、
現在の処置期間における測定された最短距離と前の処置期間における測定された最短距離とを、それらからのずれを決定するように比較するようになっており、および
前記ずれに基づいて、解剖学的特徴の動きを決定するようになっている、システム。
【請求項７２】
コントローラが、前記イメージング・システムから送信され、また該イメージング・システムによって受信される信号に基づいて、解剖学的特徴と該イメージング・システムの少なくとも一部分との間の最短距離を測定するようにさらに構成されている、請求項７１に記載のシステム。
【請求項７３】
コントローラが、前記ずれと予め決められた閾値とを比較するようになっており、該比較に基づいて、解剖学的特徴の動きを決定するようにさらに構成されている、請求項７２に記載のシステム。
【請求項７４】
前記イメージング・システムが、超音波トランスデューサ・システムを含んで成る、請求項７１に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はイメージングに関し、より詳細には、種々のイメージング・モダリティを用いて得られる画像を位置合せするための、および／または画像レジストレーションを検証するためのシステムおよび方法に関する。
【背景技術】
【０００２】
内臓の医療イメージングは、医療従事者が処置の決定をするために用いることができる重要な解剖学的情報および診断情報を供する。医療画像は、コンピュータ・トポグラフィ（またはコンピュータ断層撮影、computed topography、ＣＴ）、磁気共鳴イメージング（ＭＲＩ）、または超音波イメージングなどの様々な非侵襲的イメージング手順を用いて得ることができる。ＣＴシステムは、対象の解剖学的部位を通してＸ線を送り、Ｘ線の減衰係数に基づいて、対象の断面画像（「スライス（slice）」）を再構成することができる。結果として、ＣＴシステムは、骨構造に関する詳細を見ること、肺および胸の疾患を診断すること、および癌を検出することに非常に適している。ＣＴイメージングの利点は、例えば、完全なスキャンの短いスキャン時間（典型的には５分未満）、低コスト（ＭＲＩ装置の約半額）、体内の骨組織を正確に描出する能力、高速イメージング速度（各スキャン時間は３０秒未満）に起因する小さいモーション・アーチファクトなどを含む。しかしながら、ＣＴシステムは患者に照射し、その結果生じるリスクを負う。このため、ＣＴスキャンは妊婦や子供には推奨しない。
【０００３】
高周波音波を身体に通すことを含む超音波イメージングは、広く使用されている別の技術である。超音波は、軟組織を十分に貫通し、その短波長に起因して、数ミリメートルの寸法のスポットに焦点をあてることができる。典型的な超音波検査では、トランスデューサ・プローブが皮膚上または身体開口部の内側に直接配置される。薄い層のゲルを皮膚に適用して、トランスデューサ・プローブと皮膚との間の直接的な接触を供し、それによって体内への超音波エネルギーの効率的な伝達を可能にすることができる。内側解剖学的構造の超音波画像は、それらの構造によって反射された波、特に、例えば、反射信号の振幅および位相、および／または超音波が身体を通過するのにかかる時間から構築することができる。超音波画像はリアルタイムで捕捉されるので、体内の内臓の動きおよび血管を流れる血液も表すことができる。さらに、超音波イメージングは、高い時間分解能、音響散乱体（例えば、石灰化および気泡）に対する高い感度、優れた視覚化、低コスト、携帯性、および電離放射線被ばくが無く、したがって、一般的に、妊婦および子供に非常に安全であるとされている。
【０００４】
ＭＲＩは、対象組織を視覚化するために使用される、さらに別のイメージング・モダリティである。ＭＲＩは、以下でさらに記載するように、処置中の超音波集束をガイドするために超音波と共に使用することができる。つまり、ＭＲＩは、患者を均一な静磁場に配置し、組織内の水素核のスピンを整列させることを含む。次に、適当な周波数の高周波（ＲＦ）電磁パルス（「共振周波数」）を印加することによって、スピンを反転させて一時的にその整列を壊し、応答信号を生じさせることができる。異なる組織は、異なる応答信号を生成し、ＭＲ画像におけるこれらの組織の間のコントラストをもたらす。共振周波数および応答信号の周波数は磁場強度に依存するため、磁界強度を位置に依存させるために勾配磁場を均一磁場に重ねることによって、応答信号の起源および周波数を制御することができる。時間的に変化する勾配磁場を用いることにより、組織のＭＲＩ「スキャン」を得ることができる。
【０００５】
多くのＭＲＩプロトコルは、２または３の互いに垂直な方向に時間依存性勾配を利用する。勾配磁場およびＲＦパルスの相対強度およびタイミングは、パルス・シーケンスで特定される。時間依存磁場勾配は、例えば、脳腫瘍を視覚化し、患者の頭蓋骨に対するその位置を決定するために、ＭＲＩ応答信号の組織依存性と組み合わせて利用することができる。ＭＲＩは、（患者を動かさずに）多断面をイメージングする能力、高い信号対ノイズ比、軟組織の形態および機能におけるわずかな変化に対する高感度、ならびに放射線被爆が無いことを含む利点を有する。しかしながら、ＭＲＩは、長いスキャン時間（典型的には３０分〜数時間）、骨格構造の低解像度画像および他の高周波（ＲＦ）装置の動作による干渉に起因して患者の動きに敏感である。
【０００６】
各イメージング技術は、それ自身の長所と短所を有し、異なる種類の情報を供することができるので、異なるイメージング技術を組み合わせることが実施において有利である。例えば、ＣＴスキャンおよびＭＲＩを組み合わせることにより、軟組織間のわずかな違いだけでなく、骨構造に関して良好に詳細を供することができる。ＭＲＩおよび超音波イメージングの組合せは、術中の脳神経外科手術および胸部生検ガイダンスにおけるより良い診断のための追加の診断情報を供することが示されている。さらに、周囲の健康な組織に重大な損傷を引き起こすことなく、病的な（例えば、癌性の）組織を加熱および除去するために、超音波エネルギーを処置的に使用することができるため、ＭＲＩおよび超音波の組合せは、治療的医療処置の間、イメージング能力を供する。超音波集束システムは、一般的に、１以上の超音波ビームを生成するために、音響トランスデューサ（acoustic transducer）表面またはトランスデューサ表面のアレイを含む。トランスデューサ・アレイでは、個々の表面または「素子（element）」は、典型的に、個々に制御可能である（すなわち、それらの振動位相および／または振幅を互いに独立して設定することができ、ビームを所望の方向に操作し、所望の距離に焦点をあてることを可能にする）。超音波システムは、超音波ベースの処置を監視するのに役立つ、トランスデューサ・アレイに統合された、または別個の検出器の形態で供される受信素子（receiving elements）も含むことが多い。例えば、受信素子は、皮膚の火傷を避けるために除去される必要がある皮膚における気泡に起因して、トランスデューサと対象組織との間の界面（または境界面、interfaces）によって反射された超音波を検出し得る。受信素子はまた、過熱組織におけるキャビテーション（または流体管壁、cavitation）（すなわち、組織の液体において形成される泡の崩壊による空洞の形成）を検出するために使用され得る。
【０００７】
集束超音波トランスデューサ・システム（focused ultrasound transducer system）は、ＭＲ画像における対象組織（腫瘍など）に対するトランスデューサの位置および向きを決定するためのＭＲトラッキング・コイルまたは他のマーカーを含んでもよい。所要のトランスデューサ素子の位相および振幅の算出に基づいて、トランスデューサ・アレイは、超音波を対象に集束させるように駆動される。超音波集束自体は、ＭＲＩまたは音響共振力イメージング（ＡＲＦＩ）を用いて視覚化し得、そのような視覚化を用いて焦点位置を調整し得る。これらの方法は、一般的に超音波の磁気共鳴ガイド・フォーカシング（ＭＲｇＦＵＳ）と称される。
【０００８】
２以上のイメージング・システムを首尾よく統合し、および／または異なるシステムによって供される情報を組み合わせるためには、異なるイメージング座標系で得られる画像データを位置合せする必要がある。レジストレーションに対する従来のアプローチは、典型的に、複雑な算出手順を含み、処置目的のために十分な精度を供しない場合がある。このため、このような画像レジストレーションを検証せずに使用すると、エネルギーの不正確な適用、失敗または長時間の処置につながる場合がある。したがって、処置用途をサポートするために、十分に高速で信頼できる様式で異なるモダリティを用いて得られる画像のレジストレーションを確立し、検証する必要がある。
【発明の概要】
【０００９】
本発明の態様は、異なるモダリティ（または処理方式、modalities）を用いて、また異なる座標系で得られる画像レジストレーション（または画像位置合わせ、image registration）の精度を評価および定量化するためのシステムおよび方法を供する。種々の態様では、２つの画像座標系（例えば、ＣＴシステムおよびＭＲＩシステム）に関連する画像レジストレーションの精度は、第３イメージング（または画像化、imaging）・モダリティ（例えば、超音波システム）の使用を介して評価される。例えば、超音波座標系およびＭＲＩ座標系間の空間的関係は、超音波システムの少なくとも一部分（例えば、トランスデューサ素子のいくつか）を含むＭＲ画像を、ＭＲ画像内に含まれる超音波システムのその部分の既知の空間配置と共に得ることによって決定することができる。次いで、対象とする解剖学的部位（または身体構造における部位もしくは領域、anatomic region）の先に得られたＣＴ画像の座標は、調査（study）中の画像レジストレーションを用いてＭＲＩシステムにおける座標に変換されてもよい。結果として、対象とする解剖学的部位および超音波システム間の距離は、ＭＲＩ座標系において算出されてもよい。この距離は、超音波システムを用いて音響的に測定されてもよく、超音波座標系の測定された距離は、その後、ＭＲＩ座標系に変換されてもよい。ＭＲＩ座標系における算出された距離と測定された距離とを比較することによって、それらの間の不一致に基づいて誤差ベクトルを割り当てることができる。この誤差ベクトルは、調査中の画像レジストレーションの精度を示す。すなわち、誤差ベクトルが小さいほど、算出された距離および測定された距離のより小さいずれ（または偏差、deviation）、つまりより高いレジストレーションの精度を示す。当然のことながら、イメージング・モダリティのこの例示的な使用は、例示的な目的のみのためであり、あらゆる３つのモダリティを、用途に適切なもの（記載された役割のいずれかで機能するあらゆるモダリティ）として用いてもよい。
【００１０】
ある態様では、本発明はまた、第３モダリティ（例えば、超音波システム）を用いて、２つのイメージング・モダリティ（例えば、ＣＴシステムおよびＭＲＩシステム）の異なる座標系で得られる画像を位置合せするためのアプローチを（それらのレジストレーションを検証すること（verifying）とは対照的に）供する。例えば、超音波座標系は、超音波システムの少なくとも一部分のＭＲ画像および上記のような超音波システムの既知の空間配置の使用を介して、最初にＭＲＩ座標系に位置合せされてもよい。超音波座標系およびＣＴ座標系間のレジストレーションは、対象とする解剖学的対象の超音波画像およびＣＴ画像を介して確立してもよい。続いて、ＭＲＩ座標系およびＣＴ座標系は、超音波座標系およびＭＲＩ座標系を関連付ける画像レジストレーション、ならびに超音波座標系およびＣＴ座標系を関連付ける画像レジストレーションを用いて位置合せされてもよい。
【００１１】
本発明は、ＭＲＩ座標系およびＣＴ座標系間の画像レジストレーションの検証および／または取得のための超音波システムの使用を参照して記載されたが、そのような詳細は、本発明の範囲における制限としてみなされるべきではないことに留意されたい。例えば、上記のように、ＭＲＩ座標系は、超音波座標系およびＣＴ座標系を位置合せするために使用されてもよく、ＣＴシステムは、超音波システムおよびＭＲＩシステムを位置合せするために使用されてもよい。さらに、他のイメージング座標系間の画像レジストレーションは、本明細書に記載のアプローチを実行することによって評価および／または取得されてもよい。
【００１２】
したがって、一態様では、本発明は、第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像レジストレーションを検証する方法に関する。種々の態様では、この方法は、（ａ）第１イメージング・システムを用いて、解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得ること、（ｂ）第３イメージング・システムを用いて、第３イメージング・システムと解剖学的対象との間の距離を測定すること、（ｃ）第２イメージング・システムを用いて、解剖学的対象の第２画像を得ること、（ｄ）そのレジストレーションを用いて第１画像および第２画像を位置合せすること、ならびに（ｅ）（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）測定距離に基づいて、レジストレーションの誤差を算出することを含む。一態様では、第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムは、ＭＲＩシステム、ＣＴシステムおよび超音波トランスデューサ・システムをそれぞれ含む。
【００１３】
この方法は、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を得ることをさらに含んでもよく、その位置は、飛行時間法（time-of-flight method）に基づいて決定されてもよい。さらに、本方法は、第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出することを含んでもよい。一態様では、第３イメージング・システムの位置は、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換される。解剖学的対象の第２画像は、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換される。さらに、第３イメージング・システムおよび解剖学的対象間の距離は、第３イメージング・システムの変換された位置および第１イメージング・システムの座標系における変換された第２画像に基づいて算出される。
【００１４】
第３イメージング・システムおよび解剖学的対象間の距離は、第３イメージング・システムから送信され、第３イメージング・システムによって受信される信号に基づいて測定されてもよい。種々の態様では、方法は、第３イメージング・システムの座標系から第１イメージング・システムの座標系への測定された距離を変換することをさらに含む。次いで、算出された距離から変換された測定距離のずれに基づいて、レジストレーションの誤差が決定される。一態様では、レジストレーションの誤差は予め決められた閾値と比較され、レジストレーションの妥当性（validity）がその比較に基づいて決定される（または判定される、determined）。
【００１５】
別の態様では、本発明は、第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像のレジストレーションを検証するためのシステムに関する。ある態様では、システムは、解剖学的対象の第１画像および第３イメージング・システムの少なくとも一部分を得るための第１イメージング・システム、解剖学的対象の第２画像を得るための第２イメージング・システムならびに第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムと連通するコントローラを含む。一態様では、コントローラは、第３イメージング・システムおよび解剖学的対象間の距離を測定するようになっており、第１画像および第２画像を位置合せするようになっており、ならびに（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）測定距離に基づいて、レジストレーションにおける誤差を算出するようになっている。第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムは、ＭＲＩシステム、ＣＴシステムおよび超音波トランスデューサ・システムをそれぞれ含んでもよい。
【００１６】
コントローラは、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を決定するようにさらに構成されてもよく、その位置は、飛行時間法に基づいて決定されてもよい。さらに、コントローラは、第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システム座標系を関連付ける変換を算出するようになっていてもよい。一態様では、コントローラは、第３イメージング・システムの位置を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換する。コントローラは、解剖学的対象の第２画像を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換する。さらに、コントローラは、第１イメージング・システムの座標系における第３イメージング・システムの変換された位置および変換された第２画像に基づいて、第３イメージング・システムおよび解剖学的対象間の距離を算出するようになっている。
【００１７】
さらに、コントローラは、第３イメージング・システムから送信され、第３イメージング・システムによって受信される信号に基づいて、第３イメージング・システムおよび解剖学的対象間の距離を測定するようになっていてもよい。コントローラは、第３イメージング・システムの座標系から第１イメージング・システムの座標系に測定距離を変換してもよい。さらに、コントローラは、算出された距離からの変換された測定距離のずれに基づいてレジストレーションの誤差を決定し、レジストレーションの誤差を予め決められた閾値と比較し、その比較に基づいてレジストレーションの妥当性を決定する。
【００１８】
本発明の別の態様は、第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像を位置合せするための方法に関する。種々の態様では、この方法は、（ａ）第１イメージング・システムを用いて、解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得ること、（ｂ）第２イメージング・システムを用いて、解剖学的対象の第２画像を得ること、（ｃ）第３イメージング・システムを用いて、解剖学的対象の第３画像を得ること、（ｄ）第２画像および第３画像を位置合せすること、（ｅ）（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）位置合せされた第２画像および第３画像に基づいて、第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムの座標系を関連付ける変換を算出することを含む。一態様では、第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムは、ＭＲＩシステム、ＣＴシステム、および超音波トランスデューサ・システムをそれぞれ含む。
【００１９】
第２画像および第３画像のレジストレーションは、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を第３イメージング・システムの座標系における座標に変換することを含んでもよい。例えば、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標は、第３イメージング・システムの座標系における内側解剖学的対象に関連する座標にフィッティングされてもよい。
【００２０】
さらに、この方法は、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を得ることを含んでもよい。一態様では、その位置は、飛行時間法に基づいて決定される。この方法は、第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出することを含んでもよい。ある態様では、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換は、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換に基づいて算出される。
【００２１】
さらに別の態様では、本発明は、第１イメージング・システムおよび第２イメージング・システムを用いて得られる内側解剖学的対象の画像を位置合せするためのシステムに関する。ある態様では、システムは、解剖学的対象および第３イメージング・システムの少なくとも一部分の第１画像を得るための第１イメージング・システム、解剖学的対象の第２画像を得るための第２イメージング・システム、解剖学的対象の第３画像を得るための第３イメージング・システムならびに第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムと連通するコントローラを含む。ある態様では、コントローラは、第２画像および第３画像を位置合せするようになっており、（ｉ）第１画像における第３イメージング・システムの少なくとも一部分および解剖学的対象の位置、ならびに（ｉｉ）位置合せされた第２画像および第３画像に基づいて、第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムの座標系を関連付ける変換を算出するようになっている。第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムは、例えば、ＭＲＩシステム、ＣＴシステム、および超音波トランスデューサ・システムをそれぞれ含んでもよい。
【００２２】
コントローラは、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を、第３イメージング・システムの座標系における座標に変換するようにさらになっていてもよい。例えば、コントローラは、第２イメージング・システムの座標系における内側解剖学的対象に関連する座標を、第３イメージング・システムの座標系における内側解剖学的対象に関連する座標にフィッティングさせてもよい。
【００２３】
さらに、コントローラは、第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置を決定し得る。一態様では、その位置は、飛行時間法に基づく。さらに、コントローラは、第１画像および第３イメージング・システムの座標系における第３イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出するようになっていてもよい。コントローラは、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換に基づいて、第１イメージング・システムの座標系および第３イメージング・システムの座標系を関連付ける変換を算出するようになっていてもよい。
【００２４】
本発明のさらに別の態様は、処置シーケンス（treatment sequence）の間にて動く解剖学的特徴（anatomic feature）を検出する方法に関する。種々の態様では、この方法は、（ａ）処置シーケンスの前に、（ｉ）解剖学的特徴および第１イメージング・システムの少なくとも一部分の第１画像を得ること、ならびに（ｉｉ）解剖学的特徴と第１イメージング・システムの少なくとも一部分との間の最短距離を算出するように第１画像を処理すること、ならびに（ｂ）処置シーケンスの間にて、（ｉ）解剖学的特徴と第１イメージング・システムの少なくとも一部分との間の最短距離を測定すること、（ｉｉ）測定された最短距離とステップ（ａ）で得られる算出された最短距離とを、それらからのずれを決定するように比較すること、ならびに（ｉｉｉ）そのずれに基づいて、解剖学的特徴の動きを決定することを含む。一態様では、第１画像が、第２イメージング・システムを用いて得られ、方法が、処置シーケンスの前に、第３イメージング・システムを用いて解剖学的特徴の第２画像を得ることをさらに含む。第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムは、超音波トランスデューサ・システム、ＭＲＩシステム、およびＣＴシステムをそれぞれ含む。
【００２５】
この方法は、第１イメージング・システムの座標系における第１イメージング・システムの少なくとも一部分の位置を得ることを含んでもよく、その位置は、飛行時間法に基づいて決定されてもよい。さらに、この方法は、第１画像および第１イメージング・システムの座標系における第１イメージング・システムの少なくとも一部分の位置に基づいて、第１イメージング・システムの座標系および第２イメージング・システムの座標系を関連付ける変換を算出することを含む。一態様では、この方法は、算出された変換に基づいて、第１イメージング・システムの位置を、第１イメージング・システムの座標系から第２イメージング・システムの座標系に変換することをさらに含む。
【００２６】
ある態様では、この方法は、第１画像および第２画像を位置合せすることを含む。さらに、この方法は、解剖学的特徴の第２画像を、第３イメージング・システムの座標系から第２イメージング・システムの座標系に変換することを含む。次いで、その最短距離は、第１イメージング・システムの変換された位置および第２イメージング・システムの座標系における変換された第２画像に基づいて算出されてもよい。
【００２７】
種々の態様では、解剖学的特徴と第１イメージング・システムの少なくとも一部分との間の最短距離は、第１イメージング・システムから送信され、第１イメージング・システムによって受信される信号に基づいて測定される。一態様では、測定された最短距離は、第１イメージング・システムの座標系から第２イメージング・システムの座標系に変換される。さらに、この方法は、ずれと予め決められた閾値とを比較し、その比較に基づいて解剖学的特徴の動きを決定することをさらに含む。
【００２８】
別の態様では、本発明は、処置シーケンスの間にて動く解剖学的特徴を検出するためのシステムに関する。種々の態様では、システムは、処置シーケンスの前に、解剖学的特徴および第２イメージング・システムの少なくとも一部分の第１画像を得るための第１イメージング・システム、ならびに第１および第２システムと連通するコントローラを含む。一態様では、コントローラは、（ａ）処置シーケンスの前に、解剖学的特徴と第２イメージング・システムの少なくとも一部分との間の最短距離を算出するように第１画像を処理するようになっており、ならびに（ｂ）処置シーケンスの間にて、（ｉ）解剖学的特徴と第２イメージング・システムの少なくとも一部分との間の最短距離を測定するようになっており、（ｉｉ）測定された最短距離とステップ（ａ）で得られる算出された最短距離とを、それらからのずれを決定するように比較するようになっており、ならびに（ｉｉｉ）該ずれに基づいて解剖学的特徴の動きを決定するようになっている。一態様では、第１イメージング・システム、第２イメージング・システムおよび第３イメージング・システムは、ＭＲＩシステム、超音波トランスデューサ・システムおよびＣＴシステムをそれぞれ含む。
【００２９】
コントローラは、第２イメージング・システムの座標系における第２イメージング・システムの少なくとも一部分の位置を得るようにさらになっていてもよい。その位置は、飛行時間法に基づいて決定されてもよい。さらに、コントローラは、第１画像および第２イメージング・システムの座標系における第２イメージング・システムの少なくとも一部分の位置に基づいて、第２イメージング・システムの座標系および第１イメージング・システムの座標系を関連付ける変換を算出し得る。一態様では、コントローラは、算出された変換に基づいて、第２イメージング・システムの位置を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換するようになっている。
【００３０】
ある態様では、コントローラは、第１画像および第２画像を位置合せするようにさらに構成されている。さらに、それは、解剖学的特徴の第２画像を、第３イメージング・システムの座標系から第１イメージング・システムの座標系に変換してもよい。次いで、コントローラは、第２イメージング・システムの変換された位置および第１イメージング・システムの座標系における変換された第２画像に基づいて最短距離を算出し得る。
【００３１】
種々の態様では、コントローラは、第２イメージング・システムから送信され、第２イメージング・システムによって受信される信号に基づいて最短距離を測定するようにさらに構成されている。一態様では、コントローラは、測定された最短距離を、第２イメージング・システムの座標系から第１イメージング・システムの座標系に変換する。さらに、コントローラは、そのずれと予め決められた閾値とを比較し、その比較に基づいて解剖学的特徴の動きを決定する。
【００３２】
さらに別の態様では、本発明は、複数の処置期間を有する処置シーケンスの間にて動く解剖学的特徴を検出する方法に関する。種々の態様では、この方法は、処置シーケンスの間にて、解剖学的特徴とイメージング・システム（例えば、超音波トランスデューサ・システム）の少なくとも一部分との間の最短距離を測定すること、現在の処置期間における測定された最短距離と前の処置期間における測定された最短距離とを、それらからのずれを決定するように比較すること、ならびにそのずれに基づいて、解剖学的特徴の動きを決定することを含む。解剖学的特徴とイメージング・システムの少なくとも一部分との間の最短距離は、イメージング・システムから送信され、またイメージング・システムによって受信される信号に基づいて測定されてもよい。さらに、この方法は、ずれと予め決められた閾値とを比較し、その比較に基づいて解剖学的特徴の動きを決定することを含んでもよい。
【００３３】
別の態様では、本発明は、複数の処置期間を有する処置シーケンスの間にて動く解剖学的特徴を検出するためのシステムに関する。種々の態様では、システムは、イメージング・システム（例えば、超音波トランスデューサ・システム）であって、処置シーケンスの間にて、解剖学的特徴とそのイメージング・システムの少なくとも一部分との間の最短距離を測定するためのイメージング・システム、ならびにイメージング・システムと連通するコントローラを含む。一態様では、コントローラは、現在の処置期間における測定された最短距離と前の処置期間における測定された最短距離とを、それらからのずれを決定するように比較するようになっており、およびそのずれに基づいて、解剖学的特徴の動きを決定するようになっている。コントローラは、イメージング・システムから送信され、イメージング・システムによって受信される信号に基づいて、解剖学的特徴とイメージング・システムの少なくとも一部分との間の最短距離を測定するようにさらになっていてもよい。さらにコントローラは、ずれと予め決められた閾値とを比較し、その比較に基づいて解剖学的特徴の動きを決定してもよい。
【００３４】
種々の態様において、本明細書に記載の方法は、ヒトまたは動物の体の処置方法ではない場合がある。例えば、この方法は、処置設定の範囲外（例えば、方法のキャリブレーション（または較正、calibration）時または方法がブランクに適用される場合）で実施されてもよい。本出願を通しての本発明の方法の検討は、非処置的方法に関連してもよい。
【００３５】
本明細書で使用される場合、用語「およそ（approximately）」、「おおよそ（roughly）」および「実質的に（substantially）」は、約１０％を意味し、ある態様では約５％を意味する。本明細書を通して使用される「１の例（one example、またはan example）」、「一態様（one embodiment、またはan embodiment）」に対する言及は、例に関連して記載される特定の特徴、構造または特性が、現在の技術の少なくとも１の例において含まれることを意味する。したがって、本明細書全体の様々な箇所における「一例において（in one example、またはin an example）」または「一態様（one embodiment、またはan embodiment）」という語句は必ずしもすべて同じ例を指しているわけではない。さらに、特定の特徴、構造、ルーチン、ステップまたは特性は、技術の１以上の例において、いかなる適切な方法で組み合わせてもよい。本明細書で供される見出しは、便宜上のものであり、請求する技術の範囲または意味を、限定または解釈することを意図するものではない。
【図面の簡単な説明】
【００３６】
図面において、同様の参照符号は、一般的に異なる図を通して同じ部分を指す。また、図面は必ずしも原寸に比例するものではなく、その代わりに本発明の原理を例示することに重点が置かれている。以下の記載では、本発明の種々の態様を以下の図面を参照して記載する。
【００３７】
【図１Ａ】図１Ａは、本発明の種々の態様による例示的なＭＲＩシステムを模式的に示す。
【図１Ｂ】図１Ｂは、本発明の種々の態様による例示的なＣＴシステムを模式的に示す。
【図１Ｃ】図１Ｃは、本発明の種々の態様による例示的な超音波システムを模式的に示す図である。
【図２Ａ】図２Ａは、本発明の種々の態様による超音波座標系をＭＲＩ座標系に相関させる（または関係付ける、correlating）アプローチを示す図である。
【図２Ｂ】図２Ｂは、本発明の種々の態様によるＣＴシステムをＭＲＩ座標系に相関させるアプローチを示す図である。
【図２Ｃ】図２Ｃは、本発明の種々の態様による各トランスデューサ素子と対象の反射面における最も近い点との間の距離を算出するアプローチを示す図である。
【図３Ａ】図３Ａは、本発明の種々の態様による各トランスデューサ素子と対象の反射面における最も近い点との間の距離を音響的に測定するためのアプローチを示す図である。
【図３Ｂ】図３Ｂは、本発明の種々の態様による超音波システムから送信され、超音波システムによって受信される信号の波形を示す図である。
【図３Ｃ】図３Ｃは、本発明の種々の態様による各トランスデューサ素子と対象の反射面における最も近い点との間の距離を決定するための信号処理を示す。
【図４】図４は、本発明の種々の態様による従来のアプローチを用いて得られる画像レジストレーションの精度を評価するアプローチを示す図である。
【図５Ａ】図５Ａは、本発明の種々の態様による処置中の患者の動きを検出するためのアプローチを示す図である。
【図５Ｂ】図５Ｂは、本発明の種々の態様による処置中の患者の動きを検出するためのアプローチを示す図である。
【図６Ａ】図６Ａは、本発明の種々の態様による同じトランスデューサ・アレイおよび異なる時間における目標位置を有する距離測定の誤差関数の例示的な値を示す図である。
【図６Ｂ】図６Ｂは、本発明の種々の態様による同じトランスデューサ・アレイおよび異なる時間における目標位置を有する距離測定の誤差関数の例示的な値を示す図である。
【図７Ａ】図７Ａは、本発明の種々の態様による異なるトランスデューサ・アレイおよび目標位置を有する距離測定の誤差関数の例示的な値を示す図である。
【図７Ｂ】図７Ｂは、本発明の種々の態様による異なるトランスデューサ・アレイおよび目標位置を有する距離測定の誤差関数の例示的な値を示す図である。
【図７Ｃ】図７Ｃは、本発明の種々の態様による異なるトランスデューサ・アレイおよび目標位置を有する距離測定の誤差関数の例示的な値を示す図である。
【図７Ｄ】図７Ｄは、本発明の種々の態様による異なるトランスデューサ・アレイおよび目標位置を有する距離測定の誤差関数の例示的な値を示す図である。
【図８Ａ】図８Ａは、本発明の種々の態様による２つのイメージング座標系のレジストレーションに対するアプローチを示す。
【図８Ｂ】図８Ｂは、本発明の種々の態様による２つのイメージング座標系のレジストレーションに対するアプローチを示す。
【図８Ｃ】図８Ｃは、本発明の種々の態様による２つのイメージング座標系のレジストレーションに対するアプローチを示す。
【発明を実施するための形態】
【００３８】
図１Ａは、例示的なＭＲＩ装置１０２を示す。装置１０２は、シリンダー形状の電磁石１０４を含むことができ、電磁石１０４は、電磁石１０４のボア１０６内に必要な静磁場を生成する。医療処置中、患者は、可動支持テーブル（movable support table）におけるボア１０６の内側に位置付けられる。患者内の対象部位（region of interest）１１０（例えば、患者の頭部）は、電磁石１０４が実質的に均一場を生成するイメージング部位１１２内に位置付けられてもよい。１組のシリンダー形状の勾配磁場コイル（magnetic field gradient coils）１１３を、ボア１０６内に、また患者を取り囲むように供されてもよい。勾配磁場コイル１１３は、予め決められた時間で、また互いに直交する３方向において、予め決められた大きさの勾配磁場を生成する。勾配磁場により、異なる空間位置が異なる歳差運動周波数に関連させることができ、それによりＭＲ画像にその空間解像度を与えることができる。イメージング部位（imaging region）１１２を取り囲むＲＦ送信コイル（RF transmitter coil）１１４は、患者の組織に磁気共鳴（ＭＲ）応答信号を放射するために、イメージング部位１１２にＲＦパルスを放射する。未処理のＭＲ応答信号は、ＲＦコイル１１４によって感知され、次いで、ＭＲコントローラ１１６に渡され、次いで、ユーザーに表示され得るＭＲ画像を算出する。あるいは、別個のＭＲトランスミッタ・コイルおよびレシーバ・コイルを用いてもよい。ＭＲＩ装置１０２を用いて得られる画像は、異なる組織間の視覚的コントラストおよび従来のＸ線技術で視覚化することができない患者の解剖学的構造の詳細な内側視野を放射線医師および医師に供し得る。
【００３９】
ＭＲＩコントローラ１１６は、パルス・シーケンス、すなわち磁場勾配の相対的なタイミングおよび強度ならびにＲＦ励起パルスならびに応答検出期間を制御し得る。ＭＲ応答信号は、従来の画像処理システムを用いて増幅され、調整され、ならびに生データへとデジタル化され、さらに、当業者に既知の方法によって画像データのアレイに変換される。画像データに基づいて、処置部位（例えば、腫瘍）が同定される。画像処理システムは、ＭＲＩコントローラ１１６の部分であってもよいし、ＭＲＩコントローラ１１６と連通する別個のデバイス（例えば、画像処理ソフトウェアを含む汎用コンピュータ）であってもよい。ある態様では、１つ以上の超音波システム１２０または１以上のセンサ１２２は、以下でさらに記載するように、ＭＲＩ装置１０２のボア１０６内で移動する。
【００４０】
図１Ｂは、本発明の態様による例示的なＣＴシステム１３０を示す。ＣＴシステム１３０は、支持構造（support structure）１３４、Ｘ線源１３６およびＸ線源１３６の同じ側または反対側にあるＸ線検出器１３８を有するスキャナ１３２を含む。スキャナ１３２は、スキャンされるべき患者のためのテーブル１４０を受けるようになっている。テーブル１４０は、イメージング・シーケンスの間にてスキャンされるイメージング容積またはイメージング平面に患者を適切に位置付けるように、スキャナ１３２における開口部を通って動かすことができる。ＣＴシステム１３０は、放射線源コントローラ１４３、データ取得コントローラ１４４およびテーブル・コントローラ１４５をさらに含む。放射線源コントローラ１４３は、Ｘ線源１３６から患者に向かってＸ線放射線の放出のタイミングを調整する。データ取得コントローラ１４４は、患者を介して透過された、または患者から反射されたＸ線強度データを測定するように、Ｘ線検出器１３８と通信する。得られる強度データは、ＣＴ画像再構成のためにデータ取得コントローラ１４４またはシステム・コントローラ１４６によって処理される。テーブル・コントローラ１４５は、適用されるイメージング・プロトコルに応じて、イメージング・シーケンス中またはイメージング・シーケンス間のテーブル位置を調整するように機能する。コントローラ１４３、１４４、１４５、１４６は、別個のユニットであってもよいし、単一のユニットとして統合されてもよい。
【００４１】
さらに、システム・コントローラ１４６は、ユーザー・インターフェース１４７およびメモリー・デバイス１４８に結合されてもよい。ユーザー・インターフェース１４７は、システム・コントローラ１４６と統合されてもよく、一般的には、イメージング・シーケンスを開始し、そのようなシーケンスを制御し、およびイメージング・シーケンスの間にて得られるデータを操作するためのユーザー・ワークステーション（user workstation）を含むであろう。メモリー・デバイス１４８は、イメージング・システムに対して局所的であってもよく、またはシステムから部分的にまたは完全に離れていてもよく、ＣＴ画像再構成のための生データ、部分処理データまたは完全処理データを受けるように構成されてもよい。
【００４２】
典型的なＣＴスキャンの間、Ｘ線スキャナ１３２は患者の周りを予め決められた速度で回転し、Ｘ線源１３６は体内を通る狭いＸ線ビームを放出する。患者の身体に対して様々な角度で患者を透過または患者から反射したＸ線強度は、検出器１３８によって測定される。患者の２次元断層画像（すなわち、「スライス」）は、Ｘ線の検出されたビーム角度および強度を含む。異なる角度で得られる複数の「スライス」を処理して、ＣＴ画像を再構成することができる。ある態様では、患者のＣＴ画像は、計画段階で（すなわち、熱処理の前に）取得され、処置計画を作成する時間を可能にする。
【００４３】
図１Ｃは、本発明のいくらかの態様による例示的な超音波システム１５０を示す。図示するように、超音波システムは、ハウジング１５４の表面にアレイ１５３で配置された複数の超音波トランスデューサ素子１５２を含む。このアレイは、トランスデューサ素子１５２の単一の行またはマトリクスを含んで成ることができる。トランスデューサ素子１５２は、調整なしに配置されてもよい（すなわち、規則的な間隔が置かれる必要はなく、または規則的なパターンで配置される必要もない）。アレイは、図示するように湾曲した（例えば、球状または放物線状の）形状を有してもよく、または１以上の平面または他の形状の部分を含んでもよい。その寸法は、用途に応じて、ミリメートルおよび数十センチメートルの間で変化し得る。トランスデューサ素子１５２は、圧電セラミック素子であってもよい。圧電複合材、または一般的に、電気エネルギーを音響エネルギーに変換することができるいかなる材料を用いてもよい。素子１５２間の機械的結合を減衰させる（または弱める、damp）ために、これらは、シリコーンゴムまたは任意の他の適切な減衰材（damping material）を用いてハウジング１５４に取り付けることができる。
【００４４】
トランスデューサ素子１５２は、別個に制御可能であり、すなわち、それぞれ、他のトランスデューサの振幅および／または位相と依存しない振幅および／または位相で超音波を放射することができる。トランスデューサ・コントローラ１５６が、トランスデューサ素子１５２を駆動するために供される。ｎトランスデューサ素子の場合、コントローラ１５６は、増幅および位相遅延回路を含んで成るｎ制御回路を含んでもよく、それぞれの制御回路がトランスデューサ素子の１つを駆動する。コントローラ１５６は、典型的には０．１ＭＨｚ〜４ＭＨｚの範囲のＲＦ入力信号をｎ制御回路のためのｎチャネルに分割してもよい。集束された超音波ビームを集中して生成するように、アレイの個々のトランスデューサ素子１５２を同じ周波数、一方で異なる位相および異なる振幅で駆動するようになっていてもよい。トランスデューサ・コントローラ１５６は、所望の焦点位置１５８に対する必要な位相および振幅を算出するように、ソフトウェア、ハードウェア、ファームウェア、ハードワイヤリング、またはそれらの任意の組合せで実行され得る、算出機能を供することが望ましい。一般的に、コントローラ１５６は、例えば、周波数ジェネレータ、増幅器および位相遅延回路を含むビーム形成器、ならびに算出を実行し、個々のトランスデューサ素子１５２の位相および振幅をビーム形成器に通信するコンピュータ（例えば、汎用コンピュータ）のような種々の別個の装置を含んでもよい。そのようなシステムは、容易に入手可能であるか、過度の実験をすることなく実施することができる。
【００４５】
超音波イメージングを実行するために、コントローラ１５６は、トランスデューサ素子１５２を、イメージングされる部位に音響信号を送信し、患者の体内の様々な構造および器官からの反射信号を受信するように駆動する。各トランスデューサ素子１５２に印加されるパルスを適切に遅延させることによって、集束された超音波ビームを所望のスキャン・ラインに沿って送信することができる。患者の体内の所与のポイントから反射された音響信号は、異なる時間にトランスデューサ素子１５２によって受信される。次いで、トランスデューサ素子は、受信した音響信号をビーム形成器に供給される電気信号に変換する。各トランスデューサ素子１５２からの遅延信号は、ビーム形成器によって集約され（summed）、所与のスキャン・ラインに沿った反射されたエネルギーレベルの表示であるスキャナ信号を供する。このプロセスは、患者の体の予め決められた部位の画像を生成するための信号を供するように、複数のスキャン・ラインに対して繰り返される。典型的には、スキャン・パターンはセクター・スキャンであり、スキャン・ラインは超音波トランスデューサの中心で始まり、異なる角度に向けられる。線形、曲線またはあらゆる他のスキャン・パターンを利用することもできる。
【００４６】
超音波システム１５０は、ＭＲＩ装置１０２のボア１０６内に設けられてもよく、またはＭＲＩ装置１０２の近傍に位置付けられてもよい。組み合わされたＭＲＩ超音波システムは、処置および／または安全の目的のための超音波の適用を監視することができることが知られている。超音波システム１５０およびＭＲＩ装置１０２の相対的な位置を決定するために、超音波システム１５０は、それに関連するＭＲトラッカー（MR trackers）１６０をさらに含んでもよく、ＭＲトラッカーは、システム１５０に対して固定された位置および向きに配置される。トラッカー１６０は、例えば超音波システムのハウジングに組み込まれるか、または超音波システムのハウジングに取り付けられる。ＭＲトラッカー１６０および超音波システム１５０の相対的な位置および向きが分かっている場合、結果として生じる画像において、ＭＲトラッカー１６０を含むＭＲスキャンは、ＭＲＩ座標（ＭＲＩ装置１０２の座標系）における超音波システム１５０の位置を暗示的に明らかにする。超音波座標系をＭＲＩ座標系に関連付けることを補助するために、ある態様では、超音波システム１５０の少なくとも一部分（例えば、いくらかのトランスデューサ素子）を含むＭＲ画像が取得される。このＭＲ画像は、以下にさらに記載するように、ＭＲＩ座標系におけるトランスデューサ素子の位置（全てが得られるＭＲ画像に含まれる必要はない）を決定するように、超音波トランスデューサ素子の既知の空間配置と組み合わせて用いてもよい。
【００４７】
患者に穿刺する間、超音波は、典型的には、組織（例えば、密度、構造および結果として超音波伝播特性が異なる、骨、筋肉または脂肪）の複数の層を走査する（encounter）。組織における不均質性および異方性に起因して、超音波波面は歪んでいることが多い。さらに、異なるトランスデューサ素子からの信号は、トランスデューサ素子と、イメージングまたは処置される部位との間の材料の異なる厚さおよび輪郭、場合によっては空気充填または液体充填ポケットを走査する場合があり、異なる位相シフトおよび減衰をもたらす。したがって、多層組織構造の構造、密度、および／または厚さの知識は、トランスデューサ素子に課される適切な位相シフトおよび増幅ファクターによる影響を補整するために、また集束特性の劣化を回避するために、重要である。ＭＲＩは、一般的に、軟組織（例えば、脳）の高感度画像を供するが、ＣＴスキャンは、骨構造（例えば、頭蓋骨）に関するより詳細な画像を生成する。したがって、ＭＲＩ装置１０２およびＣＴシステム１３０から得られる画像情報を組み合わせることは有益である。そのようにするために、画像のレジストレーション、つまりＭＲＩ座標系およびＣＴ座標系のレジストレーションが当然に必要となる。
【００４８】
画像レジストレーションへの多数の従来のアプローチが利用可能である。種々の態様において、任意の所望のアプローチを用いて得られる画像レジストレーションの精度は、２つのイメージング・モダリティ（例えば、ＭＲＩ装置１０２と組み合わせた超音波システム１５０）を用いて検証および評価される。図２Ａは、超音波座標系１５０をＭＲＩ座標系１０２に関連付けるアプローチ２００を示す。第１ステップ２０２において、空間的パラメータ（例えば、超音波座標系におけるトランスデューサ素子１５２を特徴付ける配向および位置）は、適切なアプローチを用いて得られる。例えば、各トランスデューサ素子１５２は、焦点１５８または他の位置に位置付けられたセンサ（図示せず）に向かってパルスを放射してもよい。各トランスデューサ素子１５２とセンサとの間の距離は、飛行時間法に基づいて決定されてもよい。この決定された距離を用いて、超音波座標系１５０における各トランスデューサ素子の実際の位置および／または向きを決定することができる。トランスデューサ素子の空間的配置は、メモリに記憶してもよい。第２ステップ２０４において、対象とする解剖学的対象（例えば、脳および頭蓋骨を含む患者の頭部）および超音波システムの少なくとも一部分（例えば、トランスデューサ・アレイ１５３の少なくともいくらかの素子１５２）のＭＲ画像が取得される。第３ステップ２０６において、解剖学的対象および／またはＭＲ画像において取得されるトランスデューサ素子１５２の部分の空間的パラメータに関連するＭＲＩ座標が決定される。第４ステップ２０８において、（第１ステップ２０２で得られる）超音波座標系２０２および（第２ステップ２０４で得られる）ＭＲＩ座標系２０４におけるトランスデューサ素子１５２の空間的パラメータに基づいて、２つの座標系間の画像変換マトリックス（image transformation matrix）が生成される。第５ステップ２１０において、トランスデューサ素子の全体のグループ（または少なくともそのＭＲ画像で得られる部分よりも大きな部分）に関連する空間的パラメータを、超音波座標系からＭＲＩ座標系に変換するために、変換マトリックスが適用される。これらの操作の結果として、トランスデューサ素子１５２に関連する空間的パラメータはＭＲＩ座標系で表される。
【００４９】
図２Ｂは、解剖学的対象のＣＴイメージングのデータセットをＣＴ座標系１３０からＭＲＩ座標系１０２に変換するアプローチ２２０を示す。第１ステップ２２２において、ＣＴイメージングのデータセットは、処置計画段階の間、または処置段階の前に取得される。ある態様では、得られるＣＴデータセットは、解剖学的対象に関するより詳細を抽出するように処理される（ステップ２２４）。例えば、ＣＴ画像は、頭蓋骨のような骨構造を含んでもよく、頭蓋骨の表面を決定するようにＣＴデータセットの得られる密度値に閾値法（thresholding technique）および／または他の既知のイメージング処理を適用することによって、頭蓋骨をセグメント化してもよい。その後、従来の提案されたアプローチを利用して、解剖学的対象のＣＴ画像とＭＲ画像との間の画像レジストレーションを得ることができる（ステップ２２６）。得られるＭＲＩ−ＣＴ画像レジストレーションに基づいて、解剖学的対象の抽出された詳細を含むＣＴデータセットを、ＣＴ座標系からＭＲＩ座標系に変換してもよい（ステップ２２８）。次いで、ＭＲＩ装置１０２およびＣＴスキャナ１３０を用いて得られる情報をＭＲＩ座標系で組み合わせて、解剖学的対象に関する改善されたイメージング品質または付加的な情報を供することができる（ステップ２３０）。このアプローチは、関心のある対象が、軟組織および骨組織の両方を含む場合に特に有用である。ＭＲＩは軟組織（例えば、脳）の高感度画像を供し、ＣＴスキャンは骨構造（例えば、頭蓋骨）についても詳細を供する。
【００５０】
トランスデューサ素子の空間的パラメータおよびＣＴスキャンを用いて得られる情報の両方がＭＲＩ座標系の座標に変換されると、ＭＲＩ−ＣＴ画像レジストレーションの精度を評価し得る。関心のある解剖学的対象が患者の頭部である、図２Ｃを参照すると、各トランスデューサ素子１５２と、頭蓋骨表面におけるトランスデューサ素子に最も近いポイント２３４との間の距離２３２は、ＭＲＩ座標系１０２において算出することができる。算出された距離２３２は、レジストレーション頭蓋骨距離ベクトル「ＲＳＤＶ」として示される。
【００５１】
図３Ａは、飛行時間法によるアプローチ３００を使用した上記距離の音響測定（すなわち、トランスデューサ素子１５２から送信され、またはトランスデューサ素子１５２によって受信される信号のタイミングおよび超音波システムが動作する媒体における（例えば、水または空気における）既知の音速に基づく音響測定）を示す。例えば、一態様では、トランスデューサ素子のすべて（または５０％よりも多く）が送信モードに設定され、頭蓋骨に向かって短い音響パルス（〜５０μｓ）を放出する（ステップ３０２）。パルス放出の実質的に直後（例えば、数マイクロ秒後）に、トランスデューサ素子の少なくともいくらかは、頭蓋骨から反射された信号を受信するように受信モードに切り替えられる（ステップ３０４）。一態様では、すべての受信トランスデューサ素子によって受信された信号が同時に記録される（ステップ３０６）。別の態様では、トランスデューサ素子は、複数のサブセットにグループ分けされる。上記のように、各サブセットが順次活性化されて（activated）、信号を放出および受信し、ならびに所望のトランスデューサ素子のセットからの信号が含まれるまで、受信信号を順次記録する。記録された信号は、各トランスデューサ素子１５２と、頭蓋骨におけるトランスデューサ素子１５２に最も近いポイントとの間の距離ベクトル（または音響距離ベクトル「ＡＤＶ」）を算出するために処理することができる（ステップ３０８）。
【００５２】
図３Ｂは、典型的な超音波信号の記録された波形を示す。波形３１２は、波形３１４を有する反射された信号の位相測定のリファレンスとして使用され得る送信された信号を表す。記録された信号を処理してＡＤＶを判定する。例えば、図３Ｃも参照して、信号は、超音波信号からノイズ信号を除去するために、１００ＫＨｚの範囲で中央の超音波周波数（例えば、６５０ＫＨｚ）でフィルタリングしてもよい（ステップ３１６）。次に、フィルタリングされた信号にエンベロープ検出アプローチ（envelope-detection approach）（例えば、ヒルベルト変換）を適用して、反射信号３１４の振幅エンベロープ３１８を決定してもよい（例えば、図２Ｃに示す７０μｓと２００μｓとの間の時間枠（time window）において）（ステップ３２０）。次いで、振幅エンベロープ３１８の有効時間枠３２２を決定してもよい（ステップ３２４）。例えば、振幅エンベロープ３１８の開始時間および停止時間は、信号振幅がそれぞれ、予め決められた振幅閾値３２６を上回るか、または下回る場合に決定されてもよい。次いで、反射信号３１４の平均波形時間を算出することができる（例えば、決定された有効時間枠にわたって信号を積分する）（ステップ３２８）。さらに、有効時間枠３２０における反射信号３１４に関連する位相は、送信信号３１２のリファレンス位相に基づいて得られてもよい（ステップ３３０）。反射信号３１４に関連する決定された平均波形時間および位相に基づいて、反射信号３１４と送信信号３１２との間の経過時間を算出することができる（ステップ３３２）。次いで、各トランスデューサ素子１３２とトランスデューサ素子１３２に最も近い反射体との間のＡＤＶは、媒体における経過時間および音速に基づいて決定することができる（ステップ３３４）。測定されたＡＤＶは、図２に示すように得られる超音波−ＭＲＩ変換マトリックスを用いて、超音波座標系からＭＲＩ座標系に変換されてもよい。この例では、頭蓋骨の外面が超音波の最大フラクションを反射するので（他の中間層の頭蓋骨表面からの反射は無視できる）、最も近い反射体は頭蓋骨外面における最も近いポイントである。一般的に、このアプローチを用いて測定されるＡＤＶは、高い精度（中心周波数で波長の１／１０のオーダー（すなわち、０．２５ｍｍ）の精度）を有し、したがって、従来得られる画像レジストレーションを評価するために、および／または患者の動きを検出するために適している。
【００５３】
図４を参照して、算出された距離（ＲＳＤＶ）は音響的に測定された距離（ＡＤＶ）と比較されて、従来のアプローチを用いて得られる画像レジストレーションの精度を評価することができる。一態様では、ＡＤＶ測定のために活性化される個々のトランスデューサ素子ごとに、ＲＳＤＶおよびＡＤＶの比較が行われる。比較は、２つのベクトル（例えば、以下にさらに記載する誤差関数）を比較することができる適切なアプローチを用いて、ＲＳＤＶとＡＤＶとの間の相関スコアおよび／または不一致として定量化されてもよい。ＡＤＶ測定のために活性化されたすべての（または、少なくともいくらかの）トランスデューサ素子の定量化された比較値は、次いで、共に加算され、すべての（または、少なくともいくらかの）トランスデューサ素子のＲＳＤＶとＡＤＶとの間の「全体的な（entire）」相関／不一致を表す。ＭＲＩ座標系におけるＡＤＶの音響測定およびＲＳＤＶの算出の両方が高精度で決定できるため、ずれおよび／または相関スコアは、ＭＲＩ−ＣＴ画像レジストレーションの精度を示し、より小さいずれまたはより高い相関スコアは、ＭＲＩ−ＣＴ画像レジストレーションのより高い精度を示す。
【００５４】
距離ベクトルＲＳＤＶおよびＡＤＶは、適切なアプローチを用いて比較してもよい。例えば、２つの距離ベクトル間の単純な誤差関数を以下のように算出してもよい。
【数１】
2018534038.tif 000003
式中において、Ｘｒｅｇｉ、ＹｒｅｇｉおよびＺｒｅｇｉは、各トランスデューサ素子ｉに対するＲＳＤＶの値を表し、Ｘｒｅｆｉ、ＹｒｅｆｉおよびＺｒｅｆｉは、各トランスデューサ素子ｉに対するＡＤＶの値を表す。ｘｓｈｉｆｔ、ｙｓｈｉｆｔおよびｚｓｈｉｆｔは、誤差関数の最小値を得るために変化し得るグローバル・シフト値（global shift value）を表す。ＮｅｌはＡＤＶ測定のために活性化されるトランスデューサ素子の数を表す。誤差関数は、画像レジストレーションの精度を評価するために使用される。例えば、誤差関数の最小値が見出せない場合、または最小値が予め決められた閾値を超える場合、画像レジストレーションは不正確である（または、少なくとも医療目的のためには十分に正確ではない）。ある態様では、各トランスデューサ素子１５２は、レジストレーションの検証目的で、３ｍｍ未満の平均誤差を有することが可能であり、したがって、１０００個のトランスデューサ素子を有する例示的な超音波システムでは、画像レジストレーションの検証を決定するための閾値は１０４ｍｍ２に設定される。この値は、各トランスデューサ素子の許容誤差に基づいて調整することができる。トランスデューサ素子について距離が誤って算出される場合（例えば、測定ノイズなど）があるが、前述のアプローチは、多数のトランスデューサ素子の使用（例えば、１０００のオーダーでの使用）により画像レジストレーションを評価するため、誤って算出された素子の少数は、誤差関数の最小値に大きな影響を与え得ない。したがって、ＲＳＤＶおよびＡＤＶが不一致または破損している場合にのみ、高い誤差値が発生する。
【００５５】
種々の態様において、熱処理中のＡＤＶの時間発展ならびに／またはＡＤＶおよびＲＳＤＶの比較を用いて、患者の動きを検出することができる。本明細書で使用する「処置中（during treatment）」という用語は、処置セッションの全体的な時間を意味し、一般的に各超音波処理の前、間および後の時間を含む。図５Ａを参照すると、種々の態様では、超音波処理の前、間および／または後のたびに、新しいＡＤＶがリアルタイムで測定される。新しいＡＤＶは、トランスデューサ・アレイのサブセットおよび／またはより低い品質で測定してもよいため、ＡＤＶの測定期間は、中断することなく処置期間とインターリーブされる（または交互配置される、interleaved）のに十分に短くし得る（例えば、２００μｓ）。あるいは、トランスデューサ・アレイ１５３の１以上のサブセットは、ＡＤＶ測定専用であり、アレイ１５３の残りの部分は、処置目的のために超音波を集束する。さらに、ＡＤＶ測定のために頭蓋骨から反射された信号を受信するために、素子の別個の超音波トランスデューサ・アレイ（または、トランスデューサ・アレイ１５３の素子の１以上のサブセット）が供されてもよい。別個であれば、受信トランスデューサは、超音波トランスデューサ・アレイ１５３の近傍に設けられてもよく、またはそのハウジング１５４に統合されてもよい。さらに、トランスデューサ・アレイは、ＭＲＩ装置１０２のボア１０６内に設けられてもよく、またはそれの近傍に位置付けられてもよい。
【００５６】
次いで、新たに測定されたＡＤＶは、処置中のより早い時期（例えば、超音波処理が行われる前、複数の超音波処理前、または１の超音波処理前）に得られるＡＤＶと比較してもよい。新しいＡＤＶが以前に得られるＡＤＶから大幅に逸脱した場合、患者はその間に移動した可能性が高い。種々の態様では、ずれが閾値（例えば、１０００のトランスデューサ素子もしくは２．２ｍｍの平均誤差を有する各トランスデューサ素子に対して５０００ｍｍ２、またはパーセンテージで５％超もしくは１０％超）を超える場合、患者の動きが有意であるとみなされ、その動きを確認および／または補整するために是正措置を実施してもよい。
【００５７】
ある態様では、処置中のＡＤＶの時間発展を監視して、患者の動きを予測する。例えば、１０００のトランスデューサ素子を含む処置に対して、各超音波処理中にＡＤＶが５５０ｍｍ２まで増加する場合、患者は１０回目の超音波処理によって臨床的に許容される範囲を超えて動くことが予想される。したがって、患者および超音波システムの不位置合せに起因する健常な非対象組織への損傷を避けるために、９回目の超音波処理の終了時に処置を中断してもよい。
【００５８】
図５Ｂに示すように、別の態様では、動き検出は、ＡＤＶおよびＳＲＤＶの比較に基づいて判定される。例えば、画像レジストレーションの精度が検証された後、新しいＡＤＶが測定され、検証された画像レジストレーションを用いて算出されたＲＳＤＶと比較される。再び、ＡＤＶとＲＳＤＶとの間の差が閾値（例えば、１０００個のトランスデューサ素子に対して５０００ｍｍ２）を下回る場合、患者の動きは無視できるか、または臨床的に許容可能な限界内であるとみなすことができる。しかしながら、ＡＤＶおよびＲＳＤＶの差が閾値を超える場合、患者は大きく動いた可能性があり、是正措置を取ってもよい。
【００５９】
したがって、ＡＤＶとＲＳＤＶ、または以前に得られるＡＤＶとを比較することにより、対象の画像再構成の必要なしに疑わしい患者の動きを識別し、動きをリアルタイムで検出することができる。式（１）の誤差関数が最小値を取る場合、新たに測定されたＡＤＶは、算出されたＲＳＤＶ（または、以前に得られるＡＤＶ）とベストマッチする（すなわち、それらの最小ずれを有する）。したがって、最小値に関連する値および／または位置（すなわち、ｘｓｈｉｆｔ、ｙｓｈｉｆｔおよびｚｓｈｉｆｔの値）を監視することによって、患者の動きを検出し得る。動き検出の感度は、許容ずれ閾値を変更することによって調整し得る。図６Ａおよび図６Ｂは、異なる時間におけるトランスデューサ・アレイおよび頭蓋骨の同じ位置に基づく同じ画像レジストレーション（それによる同じＲＳＤＶ）を使用した、２の繰返しＡＤＶ測定に対する主軸における誤差関数の例示的な値を示す。図６Ａおよび図６Ｂは、患者の画像レジストレーションおよび／または動きの妥当性を決定するための誤差関数の最小値および／または位置を使用するアプローチが信頼でき、繰り返し可能であることを示している。患者の動きがＡＤＶプロファイルに基づいて（すなわち、新たに測定されたＡＤＶを以前に得られるＡＤＶと比較することによって）決定された場合で、患者の動きがない場合、誤差関数の最小値および位置のより高い類似性が予期される。
【００６０】
図７Ａ〜図７Ｄは、同じ画像レジストレーションを使用するが、異なる頭蓋骨位置を用いた４のＡＤＶ測定に対する主軸における誤差関数の例示的な値を示す。最小値および位置は、頭蓋骨の位置によって変化することが示されている。したがって、本発明は、画像レジストレーションの検証および／または患者の動きの検出のための、高感度で正確なアプローチを供する。
【００６１】
上記のように、比較によって有意な動きが生じている可能性があることが示された場合には、是正措置を実行してもよい。是正措置は、動きを確認するために新しいＭＲ画像を取得すること、超音波処置を一時停止すること、動きを補整するために超音波システムおよび／または患者の位置を調整することなどを含んでもよい。例えば、解剖学的対象の新たに得られるＭＲ画像は、動きを確認するために最後の有効なＭＲ画像と照合してもよい。これは、アレイ１５３におけるトランスデューサのうちのいくらかのみがＡＤＶの測定専用であり、それにより（トランスデューサ・アレイ全体を用いて測定されたＡＤＶと比較して）より低い感度を示す場合に特に有用である。ＭＲ画像比較によって動きが生じたことが明らかであれば、超音波システムおよび／または患者の位置が調整されるまで、超音波処置を中断してもよい。しかしながら、ＭＲ画像比較により、動きが生じなかったか、または臨床的に有意でないことが示された場合、処置は計画どおりに進行する。
【００６２】
ある態様では、組み合わされた超音波およびＭＲＩシステムは、処置目的のために十分な精度でＭＲＩおよびＣＴ座標系間のイメージングのレジストレーションを供する。図８Ａおよび図８Ｂに示す例示的な態様を参照すると、第１ステップ８０２において、超音波座標系は、最初に、超音波システムの少なくとも一部分のＭＲ画像を、上記のような図２Ａに関連する超音波座標系におけるトランスデューサ・アレイの空間配置を特定する情報と併せて用いてＭＲＩ座標系に位置合わせされる。この超音波−ＭＲＩのレジストレーション（例えば、変換マトリックス（transformation matrix）として表現することができる）は、超音波座標系からＭＲＩ座標系にトランスデューサ素子の空間配置を変換することを可能にする。第２ステップ８０４において、関心のある解剖学的対象の超音波画像が取得される。第３ステップ８０６において、関心のある解剖学的対象のＣＴ画像が取得される（例えば、処置計画段階中に、または処置開始前に先に得られる画像データを保存するコンピュータ・メモリからアップロードされる）。第４ステップ８０８では、得られる超音波画像がＣＴ画像と比較されて、超音波画像およびＣＴ画像が位置合せされる。超音波−ＣＴ画像のレジストレーションは、例えば、ＣＴ画像データを、相関させるアプローチまたはあらゆる他の適切なアプローチを用いて超音波画像データにフィッティングすることによって確立することができる。第５ステップ８１０において、ＣＴ座標系およびＭＲＩ座標系間の変換（例えば、変換マトリックス）は、まず、（第４ステップ８０８で得られる超音波−ＣＴ画像レジストレーションを用いて）ＣＴ座標系から超音波座標系へのＣＴイメージング・データの変換を算出し、続いて、（第１ステップ８０２で得られる超音波−ＭＲＩ画像レジストレーションを用いて）変換されたデータを超音波座標系からＭＲＩ座標系に変換する。次いで、得られるＭＲＩ−ＣＴ座標変換は、２つのイメージング・システムを用いて得られる画像を組み合わせることを可能にし、それにより、正確かつ効率的な処置計画のための、軟組織および骨構造に関する詳細を供する。
【００６３】
別の態様では、図８Ｃを参照して、第２ステップ８０４で得られる解剖学的対象の超音波画像は、第１ステップ８０２で得られる変換マトリックスを用いて、超音波座標系からＭＲＩ座標系に変換される（ステップ８１２）。次いで、ＭＲＩ座標系における解剖学的対象の変換された画像データは、ＣＴ座標系における解剖学的対象の先に得られたＣＴ画像データと比較され、ＭＲＩ−ＣＴ画像レジストレーションが算出される（ステップ８１４）。これは、例えば、解剖学的対象のＣＴ画像データを、ＭＲＩ座標系における解剖学的対象の変換された画像データにフィッティングさせることによって達成することができる。したがって、ＭＲＩ−ＣＴ画像レジストレーションは、本明細書で記載されるように適宜算出され得る。
【００６４】
本発明は、ＭＲＩ座標系およびＣＴ座標系を関連付ける座標変換を評価するため、および／もしくは取得するため、ならびに／または患者の動きを検出するための超音波システムの使用を参照して記載したが、それは本発明の範囲を限定することを意図するものではない。例えば、ＭＲＩ座標系は、超音波座標系およびＣＴ座標系を位置合せするために使用されてもよく、同様に、ＣＴシステムは、超音波システムおよびＭＲＩシステムを位置合せするために使用されてもよい。さらに、上記のアプローチを用いてあらゆる２つのイメージング・システムを関連付けるイメージングのレジストレーションを検証および／もしくは取得のため、ならびに／または患者の動きをリアルタイムでの検出のために、上記のイメージング・モダリティのいずれかの代わりに他のイメージング・モダリティを用いてもよい。
【００６５】
さらに、本明細書に記載された種々の態様の特徴は、必ずしも相互に排他的ではなく、本発明の主旨および範囲から逸脱することなく、種々の組合せおよび順列とすることが（そのような組合せおよび順列が本明細書に示されていなくても）できることを理解されたい。実際、本明細書に記載されているものの変形、修正および他の態様が、本発明の主旨および範囲から逸脱することなく当業者に想到されるであろう。
【００６６】
一般的に、２以上のイメージング・システム間のイメージングのレジストレーションを評価するためおよび／もしくは算出するため、ならびに／または患者の動きを検出するための機能は、ＭＲＩのコントローラ、超音波システムのコントローラおよび／またはＣＴのコントローラと統合されているか、または別個の外部コントローラによって供されているかどうかにかかわらず、ハードウェア、ソフトウェアまたはその両方の組合せで実施される１以上のモジュールで構成されてもよい。機能が１以上のソフトウェア・プログラムとして供される態様では、プログラムは、ＦＯＲＴＲＡＮ、ＰＡＳＣＡＬ、ＪＡＶＡ（登録商標）、Ｃ、Ｃ＋＋、Ｃ＃、ＢＡＳＩＣ、様々なスクリプト言語および／またはＨＴＭＬなどのいくらかの高レベル言語のいずれかで書かれてもよい。さらに、ソフトウェアは、対象コンピュータに存するマイクロプロセッサに向けられたアセンブリ言語で実装することができる。例えば、ソフトウェアが、ＩＢＭ ＰＣまたはＰＣクローンにおいて動作するように構成されている場合、インテル８０ｘ８６アセンブリ言語で実装してもよい。ソフトウェアは、限定されないが、フロッピー（登録商標）ディスク、ジャンプドライブ、ハードディスク、光ディスク、磁気テープ、ＰＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、フィールド−プログラマブル・ゲートアレイまたはＣＤ−ＲＯＭを含む製品で具体的に表現されてもよい。ハードウェア回路を使用する態様は、例えば、１以上のＦＰＧＡ、ＣＰＬＤまたはＡＳＩＣプロセッサを用いて実装してもよい。
【図１Ａ】
2018534038.tif 000004
【図１Ｂ】
2018534038.tif 000005
【図１Ｃ】
2018534038.tif 000006
【図２Ａ】
2018534038.tif 000007
【図２Ｂ】
2018534038.tif 000008
【図２Ｃ】
2018534038.tif 000009
【図３Ａ】
2018534038.tif 000010
【図３Ｂ】
2018534038.tif 000011
【図３Ｃ】
2018534038.tif 000012
【図４】
2018534038.tif 000013
【図５Ａ】
2018534038.tif 000014
【図５Ｂ】
2018534038.tif 000015
【図６Ａ】
2018534038.tif 000016
【図６Ｂ】
2018534038.tif 000017
【図７Ａ】
2018534038.tif 000018
【図７Ｂ】
2018534038.tif 000019
【図７Ｃ】
2018534038.tif 000020
【図７Ｄ】
2018534038.tif 000021
【図８Ａ】
2018534038.tif 000022
【図８Ｂ】
2018534038.tif 000023
【図８Ｃ】
2018534038.tif 000024
【国際調査報告】
2018534038.tif 000025
2018534038.tif 000026
2018534038.tif 000027
2018534038.tif 000028
2018534038.tif 000029