(57)【特許請求の範囲】
【請求項１】
一連の画像において周期的運動を有している所望の構造の高品質画像を取得する方法であって、
（ａ）前記一連の画像を表示して、前記所望の構造を含んでいる選択された画像を基準画像Ｒとして設定するユーザ供給入力を受け取る工程と、
（ｂ）基準画像Ｒにおいて前記所望の構造を包含する関心領域ＲＯＩ１を設定する工程と、
（ｃ）前記一連の画像からの他の画像の各々とＲＯＩ１を相互相関させて、一連の相互相関画像を結果として得る工程と、
（ｄ）各々の相互相関画像の中央に位置する第２の関心領域ＲＯＩ２を画定する工程と、
（ｅ）各相互相関画像のＲＯＩ２内でＲＯＩ２のエッジから所定の距離よりも離隔した位置にピークを有している画像を「有効」画像として設定する工程と、
（ｆ）前記「有効」画像のみを処理して、より高品質の処理済み画像を結果として得る工程とを有している、前記高品質画像を取得する方法。
【請求項２】
前記の相互相関を行う工程は、
（ａ）前記基準画像Ｒを前記ＲＯＩ１内で逆フーリエ変換してｋ空間画像とする工程と、
（ｂ）前記ｋ空間画像をそのＸ座標に関してフリップさせる工程と、
（ｃ）前記ｋ空間画像をそのＹ座標に関してフリップさせて、カーネル・セットＫを結果として得る工程と、
（ｄ）前記ＲＯＩ１の外部で０を有し且つＲＯＩ１の内部で１を有するマスクＭを作成する工程と、
（ｅ）前記マスクＭを逆フーリエ変換してｋ空間画像とし、該マスクＭのｋ空間画像をそのＸ座標に関してフリップさせる工程と、
（ｆ）前記マスクＭのｋ空間画像をそのＹ座標に関してフリップさせて、加重係数セットＷを結果として得る工程と、
（ｇ）前記カーネル・セットＫ及び加重係数セットＷを、生データ画像Ｄにそれぞれ乗算して、積Ｋ及び積Ｗをそれぞれ結果として得る工程と、
（ｈ）前記積Ｋ及び積Ｗに対して２次元前方フーリエ変換を実行して、振幅及び位相情報（Ａ2，φ2）及び（Ａ1，φ1）をそれぞれ結果として得る工程と、
（ｉ）前記情報（Ａ2，φ2）及び（Ａ1，φ1）から濃淡度‖Ｉ2‖、‖Ｉ1‖をそれぞれ決定する工程と、
（ｊ）Ｉc＝‖Ｉ2‖／‖Ｉ1‖に従って相互相関画像Ｉcを決定する工程とを含んでいる請求項１に記載の方法。
【請求項３】
前記相互相関は、ｋ空間における乗算により実行されて、高速相互相関を結果として得る請求項１に記載の方法。
【請求項４】
（ａ）各相互相関画像の中心点からの有効画像の相互相関ピークの位置からオフセットを算出する工程と、
（ｂ）前記オフセットの大きさだけ前記有効画像を平行移動させて、前記オフセットを補正する工程と、
（ｃ）前記平行移動された画像を平均して、ぼやけが減少した処理済み画像を結果として得る工程とを更に含んでいる請求項１に記載の方法。
【請求項５】
（ａ）各相互相関画像の中心点からの有効画像の相互相関ピークの位置からオフセットを算出する工程と、
（ｂ）前記オフセットに従って画像をソートする工程と、
（ｃ）前記画像をソートされた順序で再生して、前記オフセットに従って変化する「平均化された映画」を提供する工程とを更に含んでいる請求項１に記載の方法。
【請求項６】
（ａ）前記画像は、取得中に呼吸している被検体の冠状血管アンジオグラフィ画像であり、
（ｂ）前記イメージングされる構造は、冠状血管であり、
（ｃ）前記オフセットは、呼吸運動によるものであり、
（ｄ）前記「平均化された映画」は、心臓の運動を除去したものである請求項５に記載の方法。
【請求項７】
周期的運動を有している所望の構造から取得される一連の画像から導出される高品質の処理済み画像を表示するシステムであって、
（ａ）フーリエｋ空間において複数の生画像を取得すると共に該画像を接続された構成要素へ供給するデータ・ソースと、
（ｂ）データ・ソース６に結合されており、ユーザが指定した画像を選択すると共に表示し、ユーザ選択の基準画像Ｒ上でＲＯＩ１を画定するように動作する対話式画像選択装置と、
（ｃ）前記対話式画像選択装置に結合されており、ＲＯＩ１を受け取ると共に、ＲＯＩ１の外部ですべて０を有し且つＲＯＩ１の内部ですべて１を有するマスクＭを作成するように動作するマスク装置と、
（ｄ）前記対話式画像選択装置、マスク装置及びデータ・ソースに結合されており、基準画像Ｒ、マスクＭ及び生データ画像Ｄから相関画像Ｉc を作成するｋ空間相関装置と、
（ｅ）前記ｋ空間相互相関装置に結合されており、ピークが存在しているならば、存在するピークについて、ピークの高さ、及び前記相互相関画像Ｉc の中心からの該ピークのオフセットを算出するように動作する最大値検出装置と、
（ｆ）前記最大値検出装置に結合されていると共に生データ・バッファから結合されており、前記オフセット、ピークの高さ及び生データ画像Ｄを受け取り、前記ピークの高さが所定の閾値よりも大きいか否か、また前記オフセットが、「有効」画像を画定する所定の閾値よりも小さいか否かを決定するように動作すると共に、「有効」画像に対して２次元フーリエ変換を実行するように動作する画像選択装置と、
（ｇ）複数の画像を記憶すると共に検索することが可能な画像メモリと、
（ｈ）前記画像選択装置に結合されており、前記有効画像の２次元フーリエ変換表現を受け取って、該画像を更に処理した後に、該処理済み画像を画像メモリに記憶させるように動作するルータと、
（ｉ）前記画像メモリに結合されており、画像メモリに記憶されている前記処理済み画像を表示する表示装置とを有している前記高品質の処理済み画像を表示するシステム。
【請求項８】
前記ルータは、当該ソート器に供給された有効画像をそのオフセットに従ってソートすることにより画像を処理するソート器を含んでいる請求項７に記載のシステム。
【請求項９】
前記ルータは、
（ａ）画像及びオフセットを受け取って、各々の画像をそのオフセットで平行移動させるシフト装置と、
（ｂ）画像データを記憶することが可能であり、処理中の画像平均を収容している連続処理用メモリと、
（ｃ）前記連続処理用メモリ及びシフト装置に結合されており、前記処理中の平均画像及び前記シフト装置からの前記平行移動後の画像を受け取ると共に平均し、結果を連続処理用メモリに記憶させる平均器とを含んでいる請求項７に記載のシステム。
【請求項１０】
（ａ）画像取得装置は、少なくとも２つのモードで動作し、その第１のモードは簡略化された画像の「迅速」取得であり、信号を受け取ったときに高分解能取得モードである第２のモードで動作し、
（ｂ）前記画像選択装置は、前記画像取得装置に接続されており、カレントの「有効」画像が検出されたときに画像取得装置に対して信号を送信して、画像取得装置を、元の生データＤを補足する追加情報を取得するように動作させる請求項７に記載のシステム。
【発明の詳細な説明】
【０００１】
【発明の背景】
１．発明の分野
本発明は画像処理に関し、より具体的には、サイクルによって画定される運動を有する構造のイメージングに関する。
２．関連する従来技術
運動している構造はイメージングが困難であり、特に、該構造が、空間内に設定されたイメージング平面内とその外との間を移動するような場合に困難である。このイメージングは、別の周期的運動が加わると特に困難になる。例えば、呼吸している被検体の構造のイメージングは、内部構造の周期的運動を生ずる。このイメージングは更に、構造が心臓に位置している又は心臓の付近に位置しているときには、心臓の拍動によって複雑化する。心臓の拍動は、圧縮運動と捩り運動との複合した運動であり、この運動が呼吸運動に重ね合わされる。
【０００２】
拍動する心臓の血管からの画像を取得する冠状血管アンジオグラフィ等では、運動している構造の画像を取得することがしばしば望ましい。画像を強調する目的のために、いくつかの画像を平均して、画像のノイズを減少させると共に視覚化の質を向上させている。
ＭＲ（磁気共鳴）冠状動脈アンジオグラフィは、典型的には、息を止めた状態での斜方単一スライス手法か、又は呼吸ゲート式３次元手法のいずれかを用いて実行されている。しかしながら、冠状血管疾患を持つ一部の患者にとっては、繰り返し息止めすることは可能でないことがあり、また現状のナヴィゲーション手法は、多岐にわたる患者の様々な呼吸パターンの範囲を網羅して作用するロバストな方法を一般的には提供していない。更に、これらのアプローチの両者とも、イメージングの開始後しばらくの間は、成功しているのか失敗しているのかが明確でないことが多い。本研究の目的は、冠状血管ＭＲＩのための適応平均化に基づくロバストで且つ息止めを要しない手法を開発することであった。
【０００３】
現在、従来可能であったよりも更に明瞭な運動する構造の画像を形成するイメージング方法が必要とされている。
【０００４】
【発明の概要】
被検体から、周期的運動を有している所望の構造の高品質の実時間画像が取得される。
いくつかの画像が取得されて、ユーザに対して表示される。ユーザは、所望の構造を含んでいる１つの画像を基準画像Ｒとして対話式で選択する。
【０００５】
所望の構造を包含している関心領域（ＲＯＩ１）が、基準画像Ｒ上で設定される。
次いで、ＲＯＩ１を、取得された相次ぐ画像の各々と相互相関させて、相互相関画像を得る。最も効率的な実施態様では、相互相関は、フーリエ領域すなわちｋ空間領域における乗算によって実行される。
【０００６】
相互相関画像の中央に、第２の関心領域を画定する。相互相関画像がこの関心領域内でピークを示しており、且つピークが関心領域のエッジから所定の距離よりも離隔しているならば、「有効（usable）」画像として設定される。
いくつかの有効画像を取得したら、これらの画像をシフトさせて平均してもよいし、又はそれらのオフセットによってソーティングして、同様なオフセットを有している画像のみを平均してもよい。
【０００７】
このようにして、オフセットについての補償が行われているので、ぼやけの少ない平均画像が形成される。また、平均化によって、より高品質の画像が得られる。
代替的な実施態様では、カレントの画像が「有効」画像であるときに、データ・ソースへ帰還が行われる。元の生データを補強するために追加情報が取得されて、より高分解能の画像が得られる。
【０００８】
【発明の目的】
本発明の目的は、周期的運動を有している構造についての高い信号対雑音比の画像を形成するシステムを提供することにある。
本発明のもう１つの目的は、イメージング中に被検体に対して息止めを要求することなく、被検体の構造についての高い信号対雑音比の画像を形成することにある。
【０００９】
本発明のもう１つの目的は、被検体が常態で呼吸している状態で、被検体の冠状血管構造についての高い信号対雑音比の画像を形成することにある。
【００１０】
【発明の詳しい説明】
運動している構造をイメージングする際には、様々な問題に遭遇する。これらのイメージングには多くの異なる研究分野があるが、最良の用途の１つは、冠状血管アンジオグラフィであろう。この用途について以下説明するが、本発明は、構造がイメージング平面内とその外との間を移動するような周期的運動を行う構造の任意のイメージングに利用することができる。
【００１１】
心臓の血管は、被検体の心臓の拍動によってイメージング平面内とその外との間を移動する可能性がある。また、被検体の呼吸によって、冠状血管が被検体の胸部の運動に伴ってイメージング平面内とその外との間を移動する可能性もある。
所望の冠状血管がイメージング平面内に位置しているような画像を識別して、これらの画像を用いて平均画像、又は一連の画像を形成すれば、品質は大幅に向上するであろう。更に、各々の画像における冠状血管のオフセットを取り出せば、このオフセットを用いて、その系列内の他の画像と平均する前に画像の位置を補正し、これにより、イメージング平面における血管の運動に起因するぼやけを回避することができよう。
【００１２】
ユーザは、実時間イメージング・ウィンドウを観察しながらシステムと対話して、所望の血管の周囲に関心領域（ＲＯＩ１）を画定する。このＲＯＩ１を用いてカーネルを作成し、新たなフレームが取得される都度、カーネルを新たなフレームの各々と相互相関させる。相互相関を高速化するために、フーリエ（ｋ空間）領域における乗算によって相互相関を実行する。
【００１３】
基準カーネルと緊密に整合するフレームは、比較的大きな中央ピークを持つ相互相関を生ずる。解剖学的構造が平面内で平行移動しているフレームは、ピークにオフセットを生じ、また血管が平面外へ移動したフレームは、中央領域には比較的小さな信号を有する。ピーク閾値が、フレームを排除する根拠として用いられる。回転又は変形がない場合には、平均化の前に画像を平行移動させることができ、相互相関画像の中心からのピークのオフセットは、元の画像に整列させるのに必要な平行移動量を表す。中心から所定の空間オフセットよりも大きなオフセットの位置にピークを有している画像もまた排除される。
【００１４】
図１〜図４では、被検体が常態で呼吸している状態で、呼吸サイクル又は心拍サイクルのいずれかによるゲート駆動を行わないで取得された心臓アンジオグラフィ画像を示している。各々の画像は、心拍サイクル中の特定の時期を表している。
血管が明瞭に視認できる画像の１つにおいて、血管の周囲にＲＯＩ１を設定した。図１では、ＲＯＩ１は、淡色の線である冠状動脈の大部分を包含している白のボックスである。この同じ動脈は図２及び図３においても視認することができるが、図４では、心臓は、この動脈がもはやイメージング平面内に存在していない位置へ動いてしまっている。
【００１５】
図１、図２、図３及び図４に対する相互相関を、それぞれ図５、図６、図７及び図８に示す。図５、図６及び図７は、基本画像（図１、図２及び図３）とカーネルとの間の緊密な対応性によって生じている顕著な相関ピークを示している（図５は、事実上、自己相関である。）。図３の冠状血管は、呼吸運動のため平面内で平行移動しており、図７における相関ピークも同様に上方にシフトしている。許容基準として閾値が設定されており、完全なデータ・セットから、許容されたフレームが平均された。図１とシフト後の図２及び図３との平均から図９が得られた。図９は、上記の各々の画像自体よりも明瞭である。
【００１６】
図１０及び図１１は合わせて、本発明の冠状血管アンジオグラフィ・システムの簡単化したブロック図を示す。ユーザ１がユーザ・インタフェイス１３と対話し、ユーザ・インタフェイス１３は画像取得装置３へ入力を供給する。画像取得装置３は、磁気共鳴（ＭＲ）スキャナであってもよいし、ｋ空間において生データを取得するような他のイメージング装置であってもよい。画像取得装置３は、冠状血管アンジオグラフィのスライス等の画像データを取得し、これらの画像データを生データ・バッファ５に記憶させる。この生データは、未処理のＭＲデータであり、フーリエ時間領域すなわち「ｋ空間」内にある。この生データは、複数のＭＲ画像を表すものである。画像取得装置及び生データ・バッファは一緒になって、データ・ソース６として動作する。
【００１７】
生の画像は、２次元高速フーリエ変換（２ＤＦＴ）装置７へ供給され、装置７は、これら生の画像を変換して画像空間領域に戻し、循環バッファ９を画像で充填する。
ユーザ１はユーザ・インタフェイス１３と対話して、表示装置１１に対話式で表示されている循環バッファ９内の画像から、基準画像Ｒとして１つの画像を選択する。
【００１８】
次いで、ユーザ１はユーザ・インタフェイス１３と対話して、イメージング対象の所望の構造、例えば冠状血管を包含している関心領域（ＲＯＩ１）を基準画像Ｒ上で選択する。
インタフェイス１３からのＲＯＩ１はマスク装置１５へ供給され、マスク装置１５は、ＲＯＩ１の外部のすべての区域が０であり且つＲＯＩ１の内部が１であるようなマスクＭを作成する。
【００１９】
マスクＭは、循環バッファ９からの選択された基準画像Ｒと共に乗算装置１７へ供給され、そこで乗算され、その結果としてマスク済み画像Ｉm を生じる。マスク済み画像Ｉm は、ＲＯＩ１の外部では本質的にゼロとなり、基準画像ＲのＲＯＩ１の内部の部分に相当するものとなる。
マスクＭは、逆フーリエ変換装置２１へ供給されて、マスクＭのｋ空間表現を形成する。
【００２０】
このマスクＭのｋ空間表現は、その座標値のすべてがフリップＸ装置２５においてｘ軸に関してフリップさせられ、この結果をフリップＹ装置２９に渡すと、フリップＹ装置２９がこれらの係数のすべてをＹ軸に関してフリップさせて、その結果として１組の加重係数Ｗを生じ、この１組の加重係数Ｗは、加重メモリ３３に記憶される。
【００２１】
マスク済み画像Ｉm についても、それぞれ逆ＦＦＴ装置１９、フリップＸ装置２３及びフリップＹ装置２７によって同様の処理が行われる。フリップＹ装置２７の出力は、カーネル・セットＫとなり、カーネル・セットＫはカーネル・メモリ３１に記憶される。
カーネル・セットＫ及び加重係数セットＷは両者とも、ｋ空間にある。カーネル・セットＫ及び加重係数セットＷを用いて、生データ・バッファ５からのｋ空間画像データＤに乗算する。即ち、生データ・バッファ５からのｋ空間画像スライスＤが、加重係数セットＷ及びカーネル・セットＫの各々と共に、それぞれ乗算器３５及び乗算器３７へ供給される。
【００２２】
乗算器３５及び３７の結果は、それぞれ前方２次元フーリエ変換（２ＤＦＴ）装置３９及び４１へ供給されて、それぞれ画像空間表現（Ａ1 ，φ1 ）、（Ａ2 ，φ2 ）を生成する。ここで、Ａは振幅情報であり、φは位相情報である。
濃淡度装置４３及び４５が、対応するそれぞれの前方２ＤＦＴ装置から位相及び振幅を受け取って、画像の２次元に広がる各々の点における画像の濃淡度（magnitude ）を表す画像濃淡度‖Ｉ1 ‖及び‖Ｉ2 ‖を作成する。各々の位置において、Ｉ2 をＩ1 で除算して、その位置における相関の強度であるＩc を得る。これにより、２次元相関マップが得られる。
【００２３】
図中、破線で囲まれた構成要素は、ｋ空間内での相関を実行するものであり、集合的にｋ空間相関装置３０と呼ばれる。
最大値検出装置４９が、２次元相互相関画像領域内の中央に位置している第２の関心領域（ＲＯＩ２）の全体にわたって、局所的な最大強度Ｉc であるピークが存在するか否かを決定する。ＲＯＩ２は、予め決定されていて記憶されていてもよいし、ユーザ・インタフェイス１３を介してユーザによって供給されてもよい。
【００２４】
相関マップ上のＲＯＩ２内での最大画像強度Ｉc の位置を画像の中心から減算して、オフセットを決定する。最大相関画像強度Ｉc の高さが、最大相関値Ｒとなる。
画像選択装置５１が、最大値検出装置４９から前述のオフセットにあるＲを受け取ると共に、生データ・バッファ５からｋ空間画像を受け取る。
【００２５】
画像選択装置５１は、相関マップ内にピークが存在するか否かを決定すると共に、このピークが相関画像のエッジから所定の距離の範囲内に位置しているか否かを決定するための所定のロジックを内部に有している。画像選択装置５１はまた、ピークは存在していないがＲＯＩのエッジに向かって定常的に増大する相関強度値Ｉc が存在しているか否かを決定することもできる。このことは、ＲＯＩ１が相関画像の冠状血管に重なり合っていないことを示す。この場合には、画像を廃棄して、他の画像を乗算器３５及び３７へ供給し、ここまでの処理を繰り返す。
【００２６】
しかしながら、画像選択装置５１が、ピークＩc が存在しており、且つピークＩc はＲＯＩのエッジから最小距離よりも離隔して位置していると決定したならば、画像選択装置５１は、この画像をよい画像、即ち「有効」画像であるものと決定し、前方２ＤＦＴ装置５２がこの画像を画像空間へ変換し、オフセット情報と共にこの画像をルータ装置５３へ供給する。
【００２７】
１つの実施態様では、ルータ５３は単純に、画像スライスを表す情報を各々保持することが可能な複数の位置を有しているセグメント型メモリ５５に画像を記憶させる。
第２の実施態様では、ルータ５３は、メモリ５５のビン（bin ）の数よりも多い画像及びオフセットを受け取る。すると、ルータ５３内に設けられているソート器５２が、セグメント型メモリ５５内のビンの各々に、所定のオフセット範囲によって、画像を記憶させる。あるビンに配置された各々の追加画像は、このビンに記憶されている以前の画像のすべての途中経過平均と共に平均化される。
【００２８】
これにより、「平均化された映画」が得られる。平均化された映画では、心臓の運動は本質的に除去されており、冠状血管が呼吸と共に上下に運動する様子が見られる。
更にもう１つの実施態様では、ルータ５３は、画像選択装置５１に接続されているシフト装置５７として実現され得る。シフト装置５７は、各画像及び各画像のオフセットを受け取り、各々の画像をそのオフセット分だけシフトさせて画像を正規化した後に、シフト後の画像を平均化器５９へ供給する。平均化器５９はまた、連続処理用メモリ６２に既に記憶されている内容、即ち処理中の画像平均を受け取る。次いで、平均化器は、これらの画像を平均し、平均後の画像を新たな処理中の平均として連続処理用メモリ６２に記憶させる。これらの実施態様の各々は、異なる用途及び結果を有する。
【００２９】
この適応平均化された画像は、実時間画像と比べて信号対雑音比（ＳＮＲ）について顕著な向上を示す。適応平均化は、冠状血管のイメージングに適している。
更にもう１つの代替的な実施態様では、画像取得装置３が、画像選択装置５１に結合されて、画像選択装置５１に応答するようにしてもよい。画像は、通常の方式で取得してもよいし、通常のイメージングの際の取得速度を高速化するために簡略化された方式で取得してもよい。但し、画像選択装置５１が、相関画像における高い相関ピークによって示されているようにカレント画像Ｄ内に所望の構造が存在しているものと決定したら、画像取得装置３に信号が送信されて、画像の追加情報を取得する。即ち、画像取得装置３は、高速取得の際にはｋ空間において低めの空間周波数をサンプリングするが、画像選択装置５１からの信号を受け取ると、高めの空間周波数をサンプリングすることになる。この追加情報は、画像選択装置５１へ直接渡されてもよいし、先ず生データ・バッファを通過させてもよく、既に取得された生画像データＤを補強することを目的とするものである。これにより、所望の構造が視認可能であるときにはより高分解能の画像を取得し、その他の場合には、より低分解能の「迅速（quick ）」画像を取得することが可能になる。
【００３０】
このシステムには、第１の表示装置１１と、第２の表示装置７１とが存在しているので、ユーザ１に対して少なくとも２つの異なる表現を同時に表示することができる。第１の表示装置１１は、被検体の拍動する心臓の実時間画像を連続表示してもよいし、接続を切断して、循環バッファ９に保存されている画像の任意のもの又はすべてを再生してもよい。
【００３１】
第２の表示装置は、第１の表示手段１１の実時間表示と比較させながら「平均化された映画」を再生してもよいし、連続処理用メモリ６１の平均画像を表示して、平均画像をユーザの入力に従って更新してもよい。また、更新の進行を連続的に表示することも可能である。
第２の表示装置７１はまた、メモリ５５内にある画像のすべて又は任意のものを再生してもよく、これらの画像を連続的に更新してもよいし、静止させていてもよい。
【００３２】
結    果
図１〜図４の冠状血管画像は、高性能勾配システム（ＧＥメディカル・システムズ製のシグナ・ホリゾン・エコースピード（Signa Horizon EchoSpeed） ）上でスパイラル・フルオロスコピー・パルス・シーケンスを用いて取得された。このスキャナは、ウルトラスパーク２（UltraSparc 2）から制御されており、画像の戻し転送、再構成及び表示を行って、拍動する心臓の実時間画像を形成した。６インタリーブ・スパイラル・パルス・シーケンスによって、１２８×１２８のマトリクス・サイズが形成され、１秒当たり５つの完全な新たな画像を形成し、１秒当たり９回の再構成を行った。適応冠状血管イメージングに先立って、１０秒〜２０秒の範囲内で最適な冠状動脈走査平面の位置を対話式で決定するために、スクリーン方式のツールを用いた。
【００３３】
適応平均化を行う実時間スパイラル・フルオロスコピーは、冠状血管のＭＲＩのためのロバストな手法としての可能性を示している。この手法は、息止めも、別個のナヴィゲーション・シーケンスも、ＥＣＧゲート駆動も必要としない。更に、この手法は、画像の開始時から、関心のある血管を連続して視覚化することを可能にする。冠状血管の走査平面の対話式での位置決定を、これらの走査平面の実時間イメージングの際の適応平均化とインタリーブさせることも容易に行うことができる。
【００３４】
新規の発明の現状で好ましいいくつかの実施の形態についてここに詳細に記載したが、当業者には今や、多くの改変及び変形が明らかになっていることであろう。従って、特許請求の範囲は、本発明の要旨の範囲内に含まれるこれらのようなすべての改変及び変形を網羅しているものとする。
【図面の簡単な説明】
【図１】被検体の心臓の断面についてのゲート駆動なし且つ息止めなしの実時間画像である。
【図２】被検体の心臓の断面についてのゲート駆動なし且つ息止めなしの実時間画像である。
【図３】被検体の心臓の断面についてのゲート駆動なし且つ息止めなしの実時間画像である。
【図４】被検体の心臓の断面についてのゲート駆動なし且つ息止めなしの実時間画像である。
【図５】関心領域１（ＲＯＩ１）の図１との相互相関である。
【図６】関心領域１の図２との相互相関である。
【図７】関心領域１の図３との相互相関である。
【図８】関心領域１の図４との相互相関である。
【図９】図１〜図３の画像の平均化から取得される被検体の断面画像である。
【図１０】本発明の実施態様の一部の簡略ブロック図である。
【図１１】本発明の実施態様の他の部分の簡略ブロック図である。
【図１】
0004316725.tif 000002
【図２】
0004316725.tif 000003
【図３】
0004316725.tif 000004
【図４】
0004316725.tif 000005
【図５】
0004316725.tif 000006
【図６】
0004316725.tif 000007
【図７】
0004316725.tif 000008
【図８】
0004316725.tif 000009
【図９】
0004316725.tif 000010
【図１０】
0004316725.tif 000011
【図１１】
0004316725.tif 000012