JP6985056B2 - 医用画像処理装置及び医用画像処理プログラム - Google Patents

Description

本実施形態は、医用画像処理装置及び医用画像処理プログラムに関する。

心臓のＣＴスキャンにおける画質は、心臓の動きによる影響を受けることがある。心臓ＣＴスキャン取り込みの持続時間内の心臓の動きは、心臓ＣＴスキャンから導出された画像中のモーションアーチファクトの存在を招く可能性がある。

図１は、心臓のＣＴスキャンから導出された画像を示す。楕円２は、重度のモーションアーチファクトがある血管を示す。通常、著しいモーションアーチファクトのないスキャンから導出された画像は、血管がほぼ円形のように見えることが予想され得る。しかしながら、図１での楕円２中の血管は、かなり歪んでいる。

心臓の動きは、例えば毎分６０回（ｂｐｍ）を超える、高い心拍数で撮像されたスキャンの画質に、特に影響を与える可能性がある。高性能スキャナは、適度に高い心拍数では良好な画質を得ることが出来るが、もっと高い心拍数では画質の低下が依然として見られることがある。

高速回転でのスキャナは低速回転でのスキャナよりも、優れた結果を生み出す可能性があるが、モーションアーチファクトは高速回転が使用されたとしても、課題となり得る。

一定の患者については、例えば患者の心拍数を毎分６０回より下に低下させるため、患者の心拍数を減速させるために、β遮断薬などの薬物が投与されることがある。しかしながら、β遮断薬が禁忌である患者も中には居る。例えば病状のために、β遮断薬を摂ることが出来ない患者もいるかもしれない。

異なる時点の間に発生した動きを判断するために、時間において異なる点を表す画像データを比較することによって、心臓の動きを推定することが知られている。

特開２００９−２４０７２３号公報

しかしながら、心臓等の動きのある対象をＣＴスキャンする場合、動きによる画像への影響は以前として課題であり、さらなる改善の余地が求められる。

本実施形態は、上記事情に鑑み、適切な動き情報に基づいて動きの影響を低減させることができる医用画像処理装置及び医用画像撮像装置を提供するものである。

本実施形態に係る医用画像処理装置は、取得部と、生成部と、再構成部とを具備する。取得部は、Ｘ線管と被検体を介してＸ線管に対向する検出器とを、被検体の回りに相対的に回転させることでスキャン処理を実行するスキャン部を用いて、第一のスキャニング期間において収集された第一のデータセットと、第一のスキャニング期間と重なり合う第二のスキャニング期間に収集された第二のデータセットと、を取得する。生成部は、第一のスキャニング期間と第二のスキャニング期間とが重なる期間の少なくとも一部を含む期間において収集された第三のデータセットを第一のデータセットから除外して第一の部分データセットを生成し、第三のデータセットを第二のデータセットから除外して第二の部分データセットを生成し、第一の部分データセットと第二の部分データセットとを用いて、被検体の動きデータを生成する。再構成部は、動き情報と、第一のデータセットの少なくとも一部又は第二のデータセットの少なくとも一部と、を用いて、再構成処理を実行する。

モーションアーチファクトを示す心臓ＣＴ画像。 一実施形態による装置の概略図。 異なる画像に対応する回転角の範囲の概要図。 実施形態の概略的な方法を描いているフロー図。 部分再構成を収集する方法の概要を描いているフローチャート。

一実施形態による撮像データ処理装置１０が、図２に概略的に示されている。撮像データ処理装置１０は、計算装置１２、例えばパーソナルコンピュータ（ＰＣ：ｐｅｒｓｏｎａｌ ｃｏｍｐｕｔｅｒ）またはワークステーションを備え、計算装置１２は、ＣＴスキャナ１４、１つ以上のディスプレイ画面１６、そしてコンピュータキーボード、マウス、またはトラックボールなどの１つまたは複数の入力デバイス１８に接続される。

ＣＴスキャナ１４は、患者またはその他被検体の領域を描出した二次元または三次元ＣＴスキャンデータを得るように構成された任意のＣＴスキャナの可能性がある。本実施形態において、領域は心臓を具備する解剖学的領域である。他の実施形態では、領域は任意の適切な領域の可能性がある。その領域は、脳を含むことが出来る。またその領域は腹部を含むことが出来る。その領域は、任意の適切な血管（例えば冠動脈）または臓器（例えば肺または肝臓）を含むことが出来る。スキャンされる患者の領域は、測定ボリュームと呼ばれ得る。

ＣＴスキャナ１４は、Ｘ線源とガントリ上に備え付けられた受信機とを使用して、患者の領域をスキャンするよう構成されている。実施形態の中には、使用されるスキャンプロトコルは、連続的なボリューム取得であるものもある。ガントリは、本実施形態においては２７５ミリ秒で、一回転時間に患者の周りに３６０°フル回転を実行する。その他の状況では、患者の一心拍のうちにフル回転が完了することが出来る。

実施形態の中には、患者の周りをフル回転する度に、患者の領域のアキシャルスライスを描出したＣＴスキャンデータが提供される。その他の実施形態において、ＣＴスキャナ１４は、一回転における複数のスライスを取り込むよう構成されたマルチスライススキャナである。

代替実施形態で、ＣＴスキャナ１４は、例えばコーンビームＣＴスキャナ、ＭＲＩ（磁気共鳴イメージング：ｍａｇｎｅｔｉｃ ｒｅｓｏｎａｎｃｅ ｉｍａｇｉｎｇ）スキャナ、Ｘ線スキャナ、ＰＥＴ（ポジトロンエミッション断層撮像法：ｐｏｓｉｔｒｏｎ ｅｍｉｓｓｉｏｎ ｔｏｍｏｇｒａｐｈ）スキャナ、ＳＰＥＣＴ（単光子放出コンピュータ断層撮像法：ｓｉｎｇｌｅ ｐｈｏｔｏｎ ｅｍｉｓｓｉｏｎ ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ）スキャナ、または超音波スキャナ、またはハイブリッドスキャナ（例えば、ＣＴ−ＭＲまたはＣＴ−ＰＥＴスキャナ）など、任意のその他撮像モダリティにおけるスキャナによって、置き換えられたり補完されたりしてもよい。

以下の説明において、スキャンデータセット（例えば、ＣＴスキャンデータセット）という用語は、スキャナから受信される可能性のある生（再構成されていない）データを指すために使用される。ＣＴスキャンデータは、例えばＣＴスキャンの間中にスキャナによって収集される電圧データなど、ＣＴスキャンの間中にスキャナによって収集された測定を描出したものの可能性がある。ＣＴスキャンデータは、一枚またはそれ以上のアキシャルスライスを描出したデータを具備し得る。ある状況下で、ＣＴスキャンデータセットは、サイノグラムと呼ばれてもよい。

ＣＴスキャンデータセットは、データの複数のサブセットを具備する可能性があり、それぞれはスキャンの間中の異なる時間に対応し、従って異なるスキャン角に対応する。

本実施形態において、ＣＴスキャナ１４は、撮像データを収集するためのＣＴスキャンデータを再構成するように構成された、スキャナ再構成回路１５を具備する。再構成によって、スキャナ再構成回路１５は、生ＣＴデータを、空間の異なる点でのＸ線の減衰を表すボクセル輝度を具備する撮像データへと変換する。

以下の説明において、撮像データセットという用語は、再構成データ（撮像データとも呼ばれ得る）を指すために使用される。撮像データセットは、それぞれのボクセルが測定ボリュームにおける対応する位置を描出している状態で、例えば、ボクセルのアレイおよび関連する輝度を具備することが出来る。撮像データセットは、例えば表示用など、測定ボリュームの画像を生成するために、使用することが出来る。

ＣＴスキャナ１４は、例えばフィルタ逆投影など、任意の適切な方法を使用して、撮像データを収集するためにＣＴスキャンデータを再構成出来る。

本実施形態において、ＣＴスキャナ１４はＣＴスキャンデータのアキシャルスライスそれぞれに対する三つの撮像データセットを再構成するように構成されている。その他の実施形態において、ＣＴスキャナはＣＴスキャンデータのアキシャルスライスそれぞれに対し異なる数の撮像データセットを再構成するように構成されている。

アキシャルスライスそれぞれのために、ＣＴスキャンデータは、Ｘ線源のフル３６０°回転に対して収集される。所定のスライスに対する三つの撮像データセットのそれぞれは、所定のスライスに対するＣＴスキャンデータの各々の部分を使用して再構成される。それぞれの部分は、少なくとも半分の回転からのＣＴスキャンデータを含む、つまり、少なくとも１８０°の回転からのＣＴスキャンデータである。それぞれの部分は、時間においてオフセットされ、それによって回転角においてもオフセットされる。例えば、三つのデータセットは、Ｘ線源（６０°の回転）による１回転の６分の１に対応する時間の分だけの時間においてオフセットすることが出来る。

少なくとも１８０°回転からのデータは、完全な再構成を提供するために再構成されてもよく、こうして再構成された少なくとも１８０°回転からのデータは、アキシャルスライス全体の画像に対応する撮像データセットであってもよい。実施形態の中には、完全な再構成を提供するために、１８０°にスキャナのファンビームの幅を加えた角度範囲からのデータが再構成される。また実施形態の中には、１８０°を下回る参考例については、１８０°にスキャナのファンビームの幅を加えた角度範囲と置き換えられてもよい。

ある状況下で、１８０°の回転からのデータは、同時投影に対して使用することが出来る。リアルＣＴスキャナは、ファンビームスキャナの可能性があり、再構成において回転角度のより大きな範囲を使用出来る。しかし、いくつかの動き推定方法および／または動き補償方法は、ＣＴスキャナのファンビーム特質を無視する可能性があり、また同時投影を想定する可能性がある。例えば、心臓のスキャンは視野角が狭い（ｓｍａｌｌ）可能性があるので、ファンビーム角はかなり狭い可能性があり、平行なくらい十分近い可能性があり、その結果平行投影を使用して良い結果が達成され得る。

三つの撮像データセットのそれぞれは、回転角の異なる範囲からの測定を使用して収集されているフル再構成のそれぞれと共に、アキシャルスライスの完全な再構成を具備することが出来る。同じアキシャルスライスの三枚の個別の画像は、提供されてもよく、それぞれ画像が異なるスキャニング期間に対応している。

回転の間中に測定ボリュームにおいて動きが何も発生しなかった場合、三つのデータセットのそれぞれは実質的に同一である可能性が期待出来る。しかし、動きの存在は、異なる時間で（そして異なる角度で）収集された撮像データセットの間に食い違いを生じさせる可能性がある。

撮像データセットのそれぞれは、例えば図１に示されたようなモーションアーチファクトと同様のモーションアーチファクトなどを含むことが出来る。撮像データセットにおけるモーションアーチファクトは、ＣＴスキャンデータの関連のある部分が取得されたスキャニング期間の間中に発生している動きに起因する可能性がある。動きは、スキャナの１回転の半分よりも少ないタイムスケールで発生する可能性がある。例えば、心臓の少なくとも一部の形および／または位置は、第一の角度に対する計測の取得と、第一の角度から１８０°未満の角度分だけオフセットである第二の角度に対する計測の取得と、の間で変化し得る。

図３は、同じスライスの二枚の画像を描出した二つの撮像データセットＡ、Ｔが再構成されたガントリ回転角度の範囲の概略図である。第三の撮像データセットに対する角度の範囲、Ｂは同じスライスの第三の画像を描出しているが、図３には図示していない。第一の撮像データセットＡは、０°から１８０°までの角度の範囲に対する測定から再構成され、図３の領域３０と３２とを含む。第一の撮像データセットＡを再構成するために使用された角度の範囲は、垂直な網掛けを使用して表されている。第二の撮像データセットＴは、６０°から２４０°までの角度の範囲に対する測定から再構成され、図３の領域３２と３４とを含む。第一の撮像データセットＡを再構成するために使用された角度の範囲は、水平な網掛けを使用して表されている。６０°から１８０°までの角度の範囲に対する測定は、第一の撮像データセットＡの再構成において、また第二の撮像データセットＴの再構成において、両方とも使用される。図３の実施形態において、第三の撮像データセットＢは、１２０°から３００°までの角度の範囲に対する測定から再構成される（図示せず）。

本実施形態において、ＣＴスキャナ１４（例えば、Ａ、Ｔ、そしてＢ）によって再構成された撮像データセットは、メモリ２０に格納され、その後に計算装置１２に提供される。代替実施形態では、撮像データセットは、医用画像保管通信システム（ＰＡＣＳ）の一部を形成することが出来る遠隔データストア（図示せず）から供給される。メモリ２０または遠隔データストアは、任意の適切な形態の記憶装置を備え得る。

計算装置１２は、撮像データセットを自動的にまたは半自動的に処理するための処理リソースを提供し、中央処理装置（ＣＰＵ：ｃｅｎｔｒａｌ ｐｒｏｃｅｓｓｉｎｇ ｕｎｉｔ）２２を備える。本実施形態で、計算装置１２は、部分再構成回路２４と、レジストレーション回路２６と、画像生成回路２８と、を含む。

本実施形態において、部分再構成回路２４、レジストレーション回路２６、そして画像生成回路２８とは、それぞれ実施形態の方法を実行するために実行可能なコンピュータ可読指示を有するコンピュータプログラムを用いて、計算装置１２内で実施される。例えば、部分再構成回路２４、レジストレーション回路２６、そして画像生成回路２８とは、計算装置１２によって、例えばＣＰＵ２２（プロセッサ）によって実行可能である各々のコンピュータプログラムまたはアルゴリズムとしてそれぞれ実施することが出来る。しかしながら、他の実施形態で、様々なユニットは、１つまたは複数のＡＳＩＣ（特定用途向け集積回路）またはＦＰＧＡ（フィールドプログラマブルゲートアレイ）として実施することが出来る。

計算装置１２は、ハードドライブと、ＲＡＭとＲＯＭとデータバスとさまざまなデバイスドライバを含むオペレーティングシステムとグラフィックスカードを含むハードウェアデバイスとを含むＰＣの他の構成要素も含む。このような構成要素は、図が見やすいように、図２に示されていない。

図２のシステムは、図４を参照しながら以下で説明される、一実施形態の方法を実行するように構成されている。図４は、動き推定や動き補償再構成方法を具備する。三つの撮像データセットＡ、Ｔ、Ｂは上で述べられたが、簡単にするために図４の方法は、撮像データセットＲ、Ｔのペアが使用される実施形態に対して、以下で説明される。

実施形態の中に、単一のワープフィールドを作り出す、二つのボリュームＲ、Ｔから動きが推定されるものがある。ワープフィールドは、動き推定と考えられることが出来、再構成において使用出来る。実施形態には、三つの撮像データセットが使用されて、動き推定（例えば、図４を参考に以下で述べられる方法を使用して）が二度走らされるものがある。再構成ステージへと送られる、二つのワープフィールドが作り出される。その他の実施形態で、任意の適切な数の撮像データセットと任意の適切な数の動き推定とが使用出来る。

図３のフローチャートのステージ４０で、部分再構成回路２４は、メモリ２０から第一の撮像データセットＲを受信する。第一の撮像データセットは、第一のスキャニング期間の継続時間にわたってＣＴスキャナにより収集されたＣＴ測定を表すデータから、ＣＴスキャナ１４によって再構成された。第一のスキャニング期間において、ガントリは回転角度の第一の範囲（例えば、０°から１８０°）を回転する。その他の実施形態において、第一の撮像データセットＲは、ＣＴスキャナ１４からまたはデータストア、例えば遠隔データストアから受信出来る。

第一の撮像データセットＲは、本実施形態においては患者の心臓を具備する、患者の解剖学的領域を表すものである。本実施形態において、第一の撮像データセットＲは、患者の心臓を通過するアキシャルスライスを描出したものである。第一の撮像データセットＲは、アキシャルスライスに一致する、アキシャルスライスについてのフル再構成だと考えられ得る。

ステージ４２（ステージ４０と同時に起こってもよい）で、部分再構成回路２４は、ＣＴスキャナ１４から第二の撮像データセットを受信する。第二の撮像データセットＴは、第一の撮像データセットＲとしての患者の同じ領域の同じアキシャルスライスを描出している。第二の撮像データセットＴは、第二のスキャニング期間の継続時間を超えてＣＴスキャナにより収集されたＣＴ測定を表すデータから、ＣＴスキャナ１４によって再構成された。第二のスキャニング時間の間中、ガントリは回転角度の第二の範囲（例えば、６０°から２４０°）を回転する。

第一のスキャニング期間と第二のスキャニング期間とは、重なり合う。回転角度の第二の範囲は、回転角度の第一の範囲と重なり合う。

図３にもう一度戻り、本実施形態の第一の撮像データセットＲは、図３の領域３０と３２とに対応する、０°から１８０°までの角度範囲に対して取得された測定を使用して再構成される。第二の撮像データセットＴは、図３の領域３２と３４とに対応する、６０°から２４０°までの角度範囲に対して取得された測定を使用して再構成される。６０°から１８０°までの角度の範囲（領域３４）からのデータは、撮像データセットＲおよびＴと共通する。

図４のフローチャートのステージ４４で、部分再構成回路２４は、図３の領域３０（つまり、０°から６０°）によって表された角度の範囲に対する測定を描出したデータから再構成された、第一の部分再構成撮像データセットＲ’を収集するために、第一の撮像データセットを使用する。部分再構成回路２４は、図３の領域３４（つまり、１８０°から２４０°）によって表された角度の範囲に対する測定を描出したデータから再構成された、第二の部分再構成撮像データセットＴ’を収集するために、第二の撮像データセットを使用する。図３の領域３２（つまり、６０°から１８０°）によって表された角度の範囲に対する測定を描出したデータは、第一の部分再構成撮像データセットＲ’と第二の部分再構成撮像データセットＴ’との再構成から除外される。

図５は、完全な再構成（例えばＲ）から部分再構成（例えばＲ’）を収集する一方法を概要的に描いたフローチャートである。

本実施形態において、部分再構成回路２４は、図４のステージ４０で、ＣＴスキャナ１４から撮像データセットＲ（図５で入力画像６０として表されている）を収集する。本実施形態において、部分再構成回路２４は、ＣＴスキャナ再構成された撮像データセットＲからＣＴスキャンデータへのアクセスを有しない。

ステージ６２で、部分再構成回路２４は、第一のスキャニング期間の間中にスキャナによって収集された測定を表すデータを収集するために、撮像データセットＲを処理する。本実施形態において、部分再構成回路２４は、撮像データセットＲの離散フーリエ変換（ＤＦＴ）６４を作り出すために、二次元高速フーリエ変換（ＦＦＴ）を実行することで、撮像データセットＲを処理する。その他の実施形態において、任意の適切な処理方法が使用され得る。

中心スライス定理により、ＤＦＴ６４のラジアル角度は、ガントリ角度に一致し得る。所定のラジアル角度でＤＦＴ６４の原点を通り抜ける線に沿うデータは、対応するガントリ角度で取得された測定データを表すことができる。ＤＦＴ６４は、ある変換画像スライスであり、この変換画像スライスは、再構成画像の二次元ＤＦＴを演算すること（ステージ６２で実行されたように）と、またはサイノグラムのそれぞれ線の一次元ＤＦＴを演算すること、とのどちらかで、そして画像の中心を通り抜ける適切なラジアル角度での１ＤＤＦＴ結果を重なり合わせることで、生み出すことが出来る。従って、ＤＦＴ６４の中心を通る線は、サイノグラムの特定の線に一致する（そして、従って特定のガントリ角度に対応する）。

部分再構成回路２４は、撮像データセットＲ（本実施形態では、０°から１８０°）によって表された回転角度の範囲を具備する回転角度６６の範囲と、撮像データセットＴ（本実施形態では、６０°から２４０°）によって表された回転角度の範囲と、を具備する回転角度６６の範囲も、メモリ２０から（またはＣＴスキャナ１４から、または別のデータストアから）受信する。回転角度６６は、ＤＦＴ６４のどの部分をマスクで隠す（ｍａｓｋ ｏｕｔ）のかを決定するために、使用される。回転角度６６は、再構成ステージ５６において使用することも出来る。

ステージ６８で、部分再構成回路２４は、ＲとＴとに共通である角度を表すＲのＤＦＴ６４の部分をマスキングで隠すためのマスク７０を作る。ステージ６８は、ステージ６４に先立って、または後追いで、または同時に実行することが出来る。

マスク７０は、それぞれ画像スライスのフーリエ変換の部分をマスキングで隠すことにより、部分再構成を作るために使用出来る。ＤＦＴ６４を通る放射状のそれぞれ直線はサイノグラム（ＣＴスキャンデータ）の特定線に対応するので、ＤＦＴ６４の放射状の直線をマスキングで隠すことは、サイノグラムの直線により寄与する情報を削除することに対応する。マスキングすることは、動き情報を含むサイノグラムの部分から情報を残すこと、適切な角度に対する測定を表すサイノグラムの部分が重なり合うためにデータを消去すると考えられ得る。マスク７０は、所望の期間に対応するラジアル角度でＤＦＴ７４の原点を通り抜ける線を選択するように構成されている可能性がある。

ステージ７２で、部分再構成回路２４は、マスクで隠されたＤＦＴ７４を収集するために、マスク７０を撮像データセットＲのＤＦＴ６４に適用する。マスクで隠されたＤＦＴ７４は、Ｒに含まれるが、Ｔには含まれない角度に対する測定に対応するデータのみを具備する。

適切なマスクを画像Ｒのフーリエ変換に適用することで、部分再構成回路２４は、動き情報がある箇所に角度のみを残すことができ、ＲとＴとが重なり合いに対するガントリ回転角度の範囲から情報を除去すると考えることもできる。６０°から１８０°までの範囲におけるデータは、もはや含まれない。

ステージ７６で、部分再構成回路２４は、部分再構成撮像データセットＲ’（図５において部分的な画像７８と書かれている）を収集するために、マスクで隠されたＤＦＴ７４の二次元高速フーリエ逆変換（ＩＦＦＴ）を実行する。

一軸において低減された空間分解能という犠牲を払うことで、部分再構成撮像データセットＲ’は、当該Ｒ’が収集された撮像データセットＲよりも、より高い時間分解能を有することが出来る。低減された空間分解能を有する軸は、ＲとＴとの間のデータにおける重なり合いのために、ＲとＴとが位置合わせされるとしたら、有用な動き情報が何ら収集されない軸となることもある。

図５の方法は、Ｔには含まれるが、Ｒ（本実施形態においては、１８０°から２４０°までの角度である）には含まれない角度に対して収集された測定を表す、部分再構成撮像データセットＴ’を収集するために撮像データセットＴにも実行される。

なお、その他の実施形態で、図５を参考にした上述とは異なる方法を用いて、撮像データセットＲ、Ｔから部分再構成撮像データセットＲ’、Ｔ’とを収集することもできる。例えば、スキャンによって得られたＲＡＷデータ（サイノグラムデータ、投影データ）を入力データとして部分再構成撮像データセットＲ’、Ｔ’とを収集することも可能である。すなわち、部分再構成回路２４は、撮像データセットＲ、Ｔのそれぞれに対応するＲＡＷデータを入力し、ＲとＴとに共通である角度を表すＲのＲＡＷデータ（或いはＴのＲＡＷデータ）の部分をマスキングで隠すためのマスクを作る。部分再構成回路２４は、Ｒ、Ｔのそれぞれに対応するＲＡＷデータにマスクを適用し、マスクで隠されたＲ、Ｔのそれぞれに対応するＲＡＷデータに対して、例えば逆投影、高度な異方性カーネルを伴う撮像データセットの畳み込み等を実行し、図５のステージ７８に相当する部分的な再構成画像を取得する。以降、図４のステージ４６、４８以降の処理に従い、動き情報を用いた再構成処理を実現することができる。

また、再構成画像データとしての撮像データセットＲ、Ｔから順投影によりそれぞれのＲＡＷデータを生成し、得られたＲＡＷデータを用いて、部分再構成撮像データセットＲ’、Ｔ’とを収集することもできる。すなわち、部分再構成回路２４は、再構成画像データとしての撮像データセットＲ、Ｔを入力データとし、これらを順投影して撮像データセットＲ、Ｔのそれぞれに対応するＲＡＷデータを生成する。以降、ＲＡＷデータを入力データとした上記例と同様の処理によって、動き情報を用いた再構成処理を実現することができる。

撮像データセットＲ’、Ｔ’のそれぞれは、１８０°よりも小さい角度の範囲を描出（従って、撮像データセットは、フル再構成よりも寧ろ部分再構成である）。Ｒ’に対する角度（領域３０）は、Ｔ’に対する角度（領域３４）と正反対である。Ｒ’に対する角度は、Ｔ’に対する角度から１８０°分オフセットされる。

図４に戻って、ステージ４４の出力は、二つの部分再構成撮像データセットＲ’、Ｔ’である。ステージ４６で、部分再構成回路２４は、第一の部分再構成撮像データセットＲ’をレジストレーション回路２６に渡す。ステージ４８で、部分再構成回路２４は、第二の部分再構成撮像データセットＴ’をレジストレーション回路２６に渡す。

図４のステージ５０で、レジストレーション回路２６は、第一の部分再構成撮像データセットＲ’と第二の部分再構成撮像データセットＴ’とをレジストレーションする。本実施形態で、レジストレーション回路２６は、例えばＰｉｐｅｒ, Ｊ ｅｔ ａｌ, Ｏｂｊｅｃｔｉｖｅ ｅｖａｌｕａｔｉｏｎ ｏｆ ｔｈｅ ｃｏｒｒｅｃｔｉｏｎ ｂｙ ｎｏｎ−ｒｉｇｉｄ ｒｅｇｉｓｔｒａｔｉｏｎ ｏｆ ａｂｄｏｍｉｎａｌ ｏｒｇａｎ ｍｏｔｉｏｎ ｉｎ ｌｏｗ−ｄｏｓｅ ４Ｄ ｄｙｎａｍｉｃ ｃｏｎｔｒａｓｔ−ｅｎｈａｎｃｅｄ ＣＴ， Ｐｈｙｓｉｃｓ ｉｎ Ｍｅｄｉｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ ５７（６）， １７０１−１７１５ （２０１２）の方法を使用することで、非剛体レジストレーションを実行する。グローバル非剛体レジストレーション手順は、クラム−ヒル−ホークス計画を使用して計算されるワープフィールド（変形フィールド）、類似手段として相互情報量を使用して実行される。（Ｗｉｌｌｉａｍ Ｒ． Ｃｒｕｍ， Ｄｅｒｅｋ Ｌ． Ｇ． Ｈｉｌｌ， Ｄａｖｉｄ Ｊ． Ｈａｗｋｅｓ． Ｉｎｆｏｒｍａｔｉｏｎ Ｔｈｅｏｒｅｔｉｃ Ｓｉｍｉｌａｒｉｔｙ Ｍｅａｓｕｒｅｓ ｉｎ Ｎｏｎ−ｒｉｇｉｄ Ｒｅｇｉｓｔｒａｔｉｏｎ， Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＩＰＭＩ’２００３， ｐｐ．３７８−３８７）。本実施形態において、ワープフィールドは密度の高いワープフィールドであり、それぞれボクセルに対して定義された個々の置換ベクトルがある）。ワープフィールドは、時間の所定の期間を超えて、生体構造の移動を表している二次元または三次元ベクトルフィールドの可能性がある。任意のその他適切なレジストレーション手順が、代替的実施形態で使用され得る。例えば、第一の部分再構成撮像データセットＲ’と第二の部分再構成撮像データセットＴ’との差分画像を生成し、当該差分画像を基準としたレジストレーションを実行するようにしてもよい。

本実施形態において、ステージ５０の出力は、Ｒ’の中心時点とＴ’の中心時点との間の動きを表すワープフィールドである。Ｒ’の中心時点は、Ｒの中心時点と同じではなく、Ｔ’の中心時点とＴの中心時点と同じではないことに留意されたい。例えば、９０°の重なりの場合、ワープフィールドは、Ｒの中心時点の前の４５°の点とＴの中心時点の後の４５°の点との間の動きを表す（または逆も然りで、ＲまたはＴのどちらがまずスキャンされるかによる）。

ワープフィールドは、どのＲが収集されるかに対する第一のスキャニング期間とどのＴが収集されるかに対する第二のスキャニング期間との間の動きを表すと考えられ得る（第一のスキャニング期間と第二のスキャニング期間との間における時間で重なり合いに対して収集される動き情報は何もないということを考慮に入れる）。ワープフィールドは、動きの推定と考えられ得る。その他の実施形態において、動きの任意の適切な推定が、Ｒ’とＴ’とのレジストレーションから収集され得る。

ステージ５２で、レジストレーション回路２６は、ワープフィールドを画像生成回路２８へと渡す。

図４および図５の処理は、二つの画像データセットＲ、Ｔと単一のレジストレーションに関して述べられているが、その他の実施形態ではレジストレーションのその他の例があってもよい。例えば、実施形態には、レジストレーションが二例あり、再構成ステージに送られる二つのワープフィールドを生み出す。それぞれ例は、同じＲだが、異なる回転角度を伴う異なるＴボリュームを有する。

以前に知られていた方法には、単一の軸スライス（例えば、撮像データセットＲとＴで、それぞれが同じ軸スライスの完全な再構成であるなど））の二つ以上の異なる完全な再構成が、Ｒの測定時間とＴの測定時間との間に発生する動きの推定を収集するために、それぞれに対してレジストレーションされる可能性がある。

単一の軸スライスの異なる完全な再構成が、例えば、ＲとＴのスキャニング期間の重なり合いについて、スキャニング期間の重なり合いについて収集されたデータから再構成され得る。データの３６０°が所定スライスに対して収集された場合、その所定のスライスに対する完全な再構成のそれぞれが１８０°以上の角度範囲を表すデータから再構成され、その後それぞれ完全な再構成を収集するために使用されたデータは、重なり合うはずだ。

図３の領域３２によって表された角度の範囲Ｓは、撮像データセットＲとＴとで共通しているように、図３から見える可能性がある。領域３２（本実施形態において、６０°から１８０°）における角度の範囲を表す、撮像データセットＲとＴにおけるデータをレジストレーションすることで得られる有益な動き情報は、何にもないと考えられ得る。

それにより、本実施形態において、撮像データセットＲ（そして撮像データセットＴの部分ではない）の一部のみである角度の範囲を表す第一の部分再構成Ｒ’が収集され、撮像データセットＴ（そして撮像データセットＲの部分ではない）の一部のみである角度の範囲を表す第二の部分再構成Ｔ’が収集される。これらの角度の範囲におけるデータからの部分再構成は、現実の動き情報を提供するＲとＴとにおけるデータの全てを含むと考えられ得る。角度の範囲（領域３０と３４）は、矛盾なく１８０°同士であり、従って同じ空間的領域を表すと考えられ得る。

ある状況下では、部分再構成Ｒ’、Ｔ’は、乏しい空間分解能を有することが出来る（そして非常にディレクショナル）が、Ｒ’、Ｔ’が導出されたオリジナルボリュームＲ、Ｔとしてほぼ同じ動き情報を含むことが出来る。

ある状況下で、完全な再構成ＲとＴとをレジストレーションすることで収集出来るがそれよりも、更に正確な動きの推定が部分再構成のＲ’とＴ’をレジストレーションすることで、収集出来る。ボリュームＲとＴとはスキャンデータ（例えば、サイノグラムデータ）の重なり合う領域から再構成されるから、ＲとＴのボリュームデータの部分は高度に相関することもある。その他の実施形態において、ＲとＴのレジストレーションは動きを過剰または過少に推定する傾向の可能性がある。過剰または過少な動きの推定は、画像詳細の方位に依存することが出来る。過剰または過少な動きの推定は、予測するのが難しい可能性がある。さらに、例えば図１に示されたような完全な再構成において発生するモーションアーチファクトなど、モーションアーチファクトの存在下でレジストレーションを実行するのは難しい可能性もある。

ある状況下で、確度および／または整合性が増加した動き推定は、図４の方法を使用することで、収集することが出来る。動きは、ボリュームＲとＴの時間分解能よりも少ない期間を上回って推定され得る。

動きの推定は、測定ボリュームにおける任意の解剖学的構造を特定することなく収集され得る。動きの推定は、部分再構成のレジストレーションから単独で収集され得る。図４の方法は、血管が追跡される必要がないかもしれない。図４の方法は、完全に自動で実行されることが出来る。動き補償は、例えば血管近くの領域など、ボリュームのいくつかのパートに単に適用する代わりに、全体なボリュームに適用されることが出来る。

ある状況下で、ステージ４の方法は、いくつかの公知の方法を使用するよりも、動き推定のより速い方法を提供することが出来る。ステージ４の方法は、削減された時間における公知のいくつかの方法と同程度の結果を達成し得る。またある状況下では、確度はボリュームＲとＴの両方に共通する、制限する情報を除去することで改善されてもよい。

ステージ５４で、画像生成回路２８は、更なる撮像データセットを収集するために、測定ボリュームの測定を表すデータセットの再構成を実行する。再構成の処理において、画像生成回路２８は、測定を捉える間中に発生した動きに対して補償するために、推定された動きワープフィールド５２を使用する。推定動きに使用される再構成は、例えばＴａｎｇ ｅｔ ａｌ， Ａ ｃｏｍｂｉｎｅｄ ｌｏｃａｌ ａｎｄ ｇｌｏｂａｌ ｍｏｔｉｏｎ ｅｓｔｉｍａｔｉｏｎ ａｎｄ ｃｏｍｐｅｎｓａｔｉｏｎ ｍｅｔｈｏｄ ｆｏｒ ｃａｒｄｉａｃ ＣＴ， Ｐｒｏｃ． ＳＰＩＥ ９０３３， Ｍｅｄｉｃａｌ Ｉｍａｇｉｎｇ ２０１４： Ｐｈｙｓｉｃｓ ｏｆ Ｍｅｄｉｃａｌ Ｉｍａｇｉｎｇ， ９０３３０４ （１９ Ｍａｒｃｈ ２０１３）の方法を使用して、実行され得る。

本実施形態において、再構成されたデータのセットは、図５の処理が撮像データセットＴに適用された場合に、収集された撮像データセットＴのＤＦＴである。その他の実施形態において、任意の適切な方法が使用出来る。

本実施形態において、画像生成回路２８は、撮像データセットＴ（例えば、６０°から２４０°）に対するデータの角度範囲を１２の領域に分割し、それぞれの領域は１５°（約１２ミリ秒の時間に相当）の角度範囲を具備する。その他の実施形態において、任意の領域数が使用され得る。

画像生成回路２８は、例えば図５について上述の方法を使用することで、１２の領域のそれぞれに対するＤＦＴデータの個別の部分再構成を収集する。それぞれの部分再構成は、部分的な画像と呼ばれ得る。それぞれの部分的な画像は、一方向において高い空間分解能を有し、その他の方向においては非常に低い空間分解能を有することが出来る。ＴのＤＦＴのフル再構成画像が有するであろう時間分解能よりも、それぞれの部分的な画像の方がより良い時間分解能を有することができる。

個々の時間は、部分再構成のそれぞれに関連付けられる。例えば、部分再構成は、期間（例えば、１２ミリ秒辺りの期間）を超えて取得された測定を表すデータから再構成される可能性があり、部分再構成に関連付けられる時間は、データ取得の中間点の時間の可能性がある。

画像生成回路２８は、ワープフィールド５２を、部分再構成のそれぞれに関連けされた時間に対して補間する。画像生成回路２８は、それぞれが違ってもよい十二個の補間されたワープフィールドを収集する。

本実施形態において、補間は、補間を調整するために、追加のパラメータが追加された三次関数を使用することを具備する。本実施形態において、ワープフィールドは対象ボリュームの中間においてゼロに設定されている。その他の実施形態で、動きの推定を補間する任意の方法が使用することが出来る。

部分再構成のそれぞれに対し、画像生成回路２８は、ワープフィールド５２を当該部分再構成のそれぞれに関連した時間に補間することで得られた補間ワープフィールドに従って、部分再構成を変換する。

画像生成回路２８は、その後更なる撮像データセットを収集するために、変換された部分再構成を足し合わせる。

従って、本実施形態において、画像Ｔの完全な再構成は、Ｔの部分再構成を複数収集して、ワープフィールド５２を異なる部分再構成の異なる時間に補間して、補間されたワープフィールド５２を異なる部分再構成に変換するために適用して、それから変換された部分再構成を組み合わせる、ことで動きの推定に依存して調整される。動きの推定を再構成に適用することで、動きの影響が減る場合があるかもしれない。

実施形態には、ステージ５４で実行された再構成は、最終画像を収集するために、単一の軸スライスに対して収集された全てのＣＴデータの再構成であってもよい。ステージ５２の動き推定は、再構成ジオメトリを調整するために使用されてもよい。

ステージ５６で、画像生成回路２８は、ステージ５４で再構成された更なる撮像データセットを出力する。

更なる撮像データセットは、動きの影響がいくらか除去されてきたデータセットの可能性がある。撮像データセットＲ、Ｔのうちの少なくとも一方におけるモーションアーチファクトと比較した場合、モーションアーチファクトがいくらか減少することがある。例えば、より少ないモーションアーチファクトの影響は、対象フェーズデータセットＴにおいてよりも、新たなボリューメトリック撮像データセットにおいて現れる可能性があるし、またはモーションアーチファクトの深刻度は、対象フェーズデータセットＴにおけるモーションアーチファクトと比較した場合、減少することもある。

図４の方法は、ボリューメトリックＣＴデータセットにおけるそれぞれ軸スライスに対して、実行出来る。実施形態には、図４の方法は、マルチスライススキャナから収集されたマルチスライスボリュームに対して、実行することが出来る。

それぞれ軸スライスに対して、計算装置１２は、任意の適切な数の撮像データセットを受信出来る。計算装置１２は、重なり合うデータが撮像データセットから除外された、二つ以上の撮像データセットの部分再構成をレジストレーションすることで、一つ以上の動きの推定を収集することが出来る。それぞれスライスに対し、動きの推定または複数の推定は、動きが補償された、更なる撮像データセットを再構成する際に使用することが出来る。

実施形態には、図４の処理は三つの撮像データセットを使用して実行される。レジストレーション回路２６は、一つの最終再構成へと送る二つのワープフィールドを決定する。三つの撮像データセットのうちの一つは、双方のレジストレーションに対して参照（Ｒ）として、データの異なる部分がマスクで隠された状態で使用され得る。

実施形態の中には、撮像データセットを二つのみ（若しくは一つだけ）を使用して実行される動き補正は、三つ以上の撮像データセットを使用して実行される動き補正よりも、詳細な動き推定の提供が少なくなってしまう可能性がある。

実施形態の中には、動きの多重の推定が最終画像を再構成する際に使用され得る。例えば、部分再構成の異なる一組が、異なる時間に対しておよび／または異なる角度に対して、動き推定を収集するために使用され得る。ある状況下では、使用する推定をいくつにするかを決定する場合において、考慮することが出来る空間分解能と時間分解能との間のトレードオフが存在するだろう。

実施形態の中には、ＣＴスキャナ１４は、動きの推定を収集するために、計算装置１２によって受信され処理された、三つの撮像データセット（例えば、Ｔ、Ａ、そしてＢ）を再構成する。その他の実施形態で、任意の数の撮像データセットがボリューメトリックＣＴデータから再構成されている可能性があり、それぞれは異なる時点に対応する。例えば、対象フェーズデータセットと、一つ、三つ、四つ、または五つの参考フェーズデータセットは、ＣＴデータから再構成されている可能性がある。

図４を参考した上述の方法は、低減されたモーションアーチファクトで再構成データセットを生み出すために使用することが出来る。図４の方法は、冠動脈の動きの改善された推定を提供し得る。モーションアーチファクトを低減することにより、ＣＴスキャナの回転速度を上げることなく、画質を改善することが可能になるかもしれない。

ある状況下で、その他の動き補償方法によってもたらされる可能性がある濃い影のアーチファクトは、図５の方法を使用する動き補償において、もたらされない可能性がある。血管および／または石灰化の画像が改善され得る。濃い影のアーチファクトは、例えば肝静脈近くの濃い影が、改善され得る。

いくつかのＣＴシステムにとって、以前に収集が可能であった心臓画像よりも、より高い心拍数での高画質の心臓画像を収集することが、可能になるかもしれない。多重心臓信号についてスキャンする方法を使用しながら収集するよりも、より短いスキャン時間（従って、低い放射線量）で高画質の心臓の画像を収集することが可能かもしれない。より頻繁により高い画質を収集することが可能になるかもしれず、つまりこれは実行する必要のある可能性があるスキャンの数がより少なくなることを意味するかもしれない。より心拍数の高い患者を上手くスキャンすることが可能になるかもしれない。より少ない薬（例えば、β遮断薬）の使用で高い心拍数で、高い画質の画像を収集することが可能となるかもしれない。

動き推定に対する完全に再構成されたボリュームをレジストレーションする代わりに、部分再構成が、有益な動き情報を含むと考えられるデータのみから作り出される。部分再構成は、一軸における空間分解能を犠牲にして、完全に再構成されたボリュームよりもより良い時間分解能を有し得る。空間分解能が欠如した部分再構成に対する軸は、有益な動き情報が少ないか、全く含まない可能性がある。

図４と５を参考に上で述べられた実施形態において、計算装置１２は、撮像データセット（例えば、ＲとＴ）を受信し、部分的に再構成された撮像データセットを収集するために撮像データセット（例えば、Ｒ’とＴ’）を処理する。

その他の実施形態において、計算装置１２は、スキャン１４からスキャンデータ（例えば、ＣＴスキャンデータ）を受信する。計算装置１２は、部分再構成撮像データセットＲ’とＴ’の再構成の際に使用される予定のスキャンデータの部分を選択する。実施形態の中には、計算装置１２は、任意の完全な再構成（例えば、撮像データセットＲとＴとを再構成することなく）を第一に実行することなく部分再構成撮像データセットＲ’とＴ’とを再構成する。図３の例示的な角度に再び関して、一実施形態における部分再構成回路２４は、スキャンデータのセットを受信し、角度範囲３０（０°から６０°）において収集された測定を具備するスキャンデータの部分を選択する。角度範囲３０におけるスキャンデータは、第一のスキャニング期間（０°から１８０°のスキャンに撮られた時間）の間中に収集されたものだが、第二のスキャニング期間（６０°から２４０°のスキャンに撮られた時間）の間中に収集されたものではない。部分的再構成回路２４は、選択されたスキャンデータの部分から部分再構成撮像データセットＲ’を再構成する。

再構成回路２４は、角度範囲３４（１８０°から２４０°）において収集された測定を具備するスキャンデータの部分を選択する。角度範囲３０におけるスキャンデータは、第二のスキャニング期間（６０°から２４０°のスキャンに撮られた時間）の間中に収集されたものだが、第一のスキャニング期間（０°から１８０°のスキャンに撮られた時間）の間中に収集されたものではない。部分再構成回路２４は、スキャンデータの選択された部分から部分再構成撮像データセットＴ’を再構成する。レジストレーション回路２６は、部分再構成撮像データセットＲ’とＴ’からの動き推定を収集し、画像生成回路２８は、例えば図４を参考に上で述べられた方法を使用して、軸スライスに対するスキャンデータの完全な再構成である、撮像データセットを再構成する。

更なる実施形態において、スキャナ１４のスキャナ再構成回路１５は、部分再構成撮像データセットＲ’とＴ’とを収集するためにスキャンデータの部分を部分再構成回路２４に代わって、再構成する。実施形態の中には、スキャナ１４は、部分再構成撮像データセットを計算装置１２へと受け渡す。更なる実施形態において、図４の処理のうちのいくつかまたは全ては、スキャナ１４において、または任意の適切な装置において実行される。

本実施形態では、撮像データを提供するために、スキャナはＣＴデータを収集し、ＣＴデータを再構成する。その他の実施形態では、スキャナは、任意の適切なモダリティのデータを収集し（例えば、ＣＴ、コーンビームＣＴ、ＭＲ、ＰＥＴ、ＳＰＥＣＴ、Ｘ線、または超音波）、収集されたモダリティのデータから撮像データセットを再考するための適切な再構成方法を使用する。スキャナは、任意の適切な二次元または三次元撮像データセットを提供するために、スキャナからデータを再構成出来る。実施形態の中には、スキャナがハイブリッドスキャナ（例えば、ＣＴ−ＭＲまたはＣＴ−ＰＥＴスキャナ）で、図２の方法がハイブリッドスキャナからのデータのＣＴ部分に適用されるものがある。

実施形態には、スキャナがコーンビームＣＴスキャナのものがある。実施形態には、撮像データセットは、血管造影法撮像データセットのものがある。コーンビームＣＴは、その他のＣＴスキャン方法よりも、遅い回転速度を有してもよい。ある状況において、コーンビームＣＴスキャナは、回転の間に望ましくない動き（例えば、ぐらつくような動き）を経験する可能性がある。コーンビームＣＴ実施形態の中には、二つ以上のレジストレーションが動きの推定の際に、時間にわたって移動のより詳細なモデルを収集するために、実行される。

ボリューム測定は、例えば腹部領域など、ヒトまたは動物被検体の任意の適切な解剖学的領域であってよい。解剖学的領域は、例えば任意の臓器（例えば、心臓、脳、肺、または肝臓）や血管（例えば、冠動脈）など、任意の適切な解剖学的領域を具備することが出来る。

特定の実施形態は、医用撮像方法を提供してもよい。

患者の所定部分の二つの二次元または三次元ＣＴ画像データセットを受信すること、それぞれデータセットは、所定の期間を超えたサイノグラムデータから再構成され、二つのデータセットの期間が重なるようにして、それぞれ画像データセットに対し、その他データセットの期間と重ならない期間の範囲を決定すること、それぞれ画像データセットの部分再構成を計算すること、その他のデータセットの分け合わない情報のみを含むこと、二つの画像データセット間の動きの推定を収集するために、結果として生じる部分再構成の位置合わせ、を具備する。

部分再構成の計算は、それぞれデータセットのそれぞれ画像の二次元フーリエ変換（ＤＦＴ）を計算すること、ラジアル角度でそれぞれＤＦＴ画像の原点を通過する線のみを選択するマスクを重ならない期間に適用すること、結果として生じるマスクで隠された画像に二次元逆ＤＦＴを計算すること、を具備する。

データセットのうちの少なくとも一つは、モーションアーチファクトにより影響を受けていてもよい。患者の部分は、心臓であってよい。

特定の実施形態は、医用撮像方法を提供することが出来る。患者の所定の部分の少なくとも二つのＣＴ画像データセットを受信し、一つは対象フェーズとして呼ばれ、対象フェーズデータセットに対向する対象フェーズを除くそれぞれデータセットをレジストレーションし、患者の移動を記号化する動きフィールド（ｍｏｔｉｏｎ ｆｉｅｌｄ）のセットを収集し、動きフィールドとオリジナル画像データセットに基づいて、モーションアーチファクトによりより少なく影響を受けた更なる画像データセットを再構成すること、を具備する。

本実施形態では、特定のユニットについて説明してきた。いくつかの実施形態では、これらのユニットのうち１つまたは複数の機能は単一の処理リソースまたは他の構成要素によって提供可能であり、または、単一のユニットによって提供される機能は、組み合わされた２つ以上の処理リソースまたは他の構成要素によって提供可能である。単一のユニットへの言及は、そのユニットの機能を提供する複数の構成要素が互いに遠隔であるかどうかに関わりなく、そのような構成要素を包含し、複数のユニットへの言及は、それらのユニットの機能を提供する単一の構成要素を包含する。

他の実施形態として、第一のスキャニング期間の間中に測定ボリュームに相対的な医用スキャナの回転によって収集された前記測定ボリュームの少なくともいくつかの測定を描出した第一のデータセットを収集し、前記第一のスキャニング期間と重なり合う第二のスキャニング期間の間中に前記測定ボリュームに相対的な前記医用スキャナの回転によって収集された前記測定ボリュームの少なくともいくつかの測定を描出した第二のデータセットを収集し、前記第一のデータセットと前記第二のデータセットとに基づいた、前記第一のスキャニング期間と前記第二のスキャニング期間との間の動きの推定を収集するための手順を実行する、よう構成された処理回路を具備する医用撮像データ処理装置を提供し、前記収集は、前記第一のスキャニング期間と前記第二のスキャニング期間との間の重なりの間中に収集された前記測定を表す前記データの少なくともいくつかを、前記第一のデータセットおよび前記第二のデータセットから除外するようになっており、前記第一のスキャニング期間と前記第二のスキャニング期間のうちの少なくとも一方の間中に収集された測定を描出したデータセットは、前記測定ボリュームの少なくとも一部を描出した医用撮像データセットの再構成における使用に適している。

特定の実施形態は、第一のスキャニング期間の間中に測定ボリュームに相対的な医用スキャナの回転によって収集された前記測定ボリュームの少なくともいくつかの測定を描出した第一のデータセットを収集すること、前記第一のスキャニング期間と重なり合う第二のスキャニング期間の間中に前記測定ボリュームに相対的な前記医用スキャナの回転によって収集された前記測定ボリュームの少なくともいくつかの測定を描出した第二のデータセットを収集すること、前記第一のデータセットと前記第二のデータセットとに基づいた、前記第一のスキャニング期間と前記第二のスキャニング期間との間の動きの推定を収集するための手順を実行することと、を具備する医用撮像データ処理方法を提供し、前記収集は、前記第一のスキャニング期間と前記第二のスキャニング期間との間の重なりの間中に収集された前記測定を表す前記データの少なくともいくつかを、前記第一のデータセットおよび前記第二のデータセットから除外するようになっており、前記第一のスキャニング期間と前記第二のスキャニング期間のうちの少なくとも一方の間中に収集された測定を描出したデータセットは、前記測定ボリュームの少なくとも一部を描出した医用撮像データセットの再構成における使用に適している。

特定の実施形態について説明してきたが、これらの実施形態は、例として提示したにすぎず、本発明の範囲を限定することを意図したものではない。実際、本明細書で説明する新規な方法およびシステムはさまざまな他の形態で実施することが出来る。そのうえ、本明細書で説明する方法およびシステムの形態におけるさまざまな省略、置き換え、および変更は、本発明の趣旨から逸脱することなく行うことが出来る。添付の特許請求の範囲およびその等価物は、本発明の範囲に含まれるこのような形態または変形形態を包含することを意図するものである。

１０…撮像データ処理装置、１２…計算装置、１４…ＣＴスキャナ、１５…スキャナ再構成回路、１６…ディスプレイ画面、１８…入力デバイス、２０…メモリ、２２…ＣＰＵ、２４…部分再構成回路、２６…レジストレーション回路、２８…画像生成回路

Claims (13)

1. Ｘ線管と被検体を介して前記Ｘ線管に対向する検出器とを、前記被検体の回りに相対的に回転させることでスキャン処理を実行するスキャン部を用いて、第一のスキャニング期間において収集された第一のデータセットと、前記第一のスキャニング期間と重なり合う第二のスキャニング期間に収集された第二のデータセットと、を取得する取得部と、
   前記第一のスキャニング期間と前記第二のスキャニング期間とが重なる期間の少なくとも一部を含む期間において収集された第三のデータセットを前記第一のデータセットから除外して第一の部分データセットを生成し、前記第三のデータセットを前記第二のデータセットから除外して第二の部分データセットを生成し、前記第一の部分データセットと前記第二の部分データセットとを用いて、前記被検体の動きデータを生成する生成部と、
   前記動きデータと、前記第一のデータセットの少なくとも一部又は前記第二のデータセットの少なくとも一部と、を用いて、再構成処理を実行する再構成部と、
   を具備し、
   前記生成部は、
   前記第一のデータセットに基づく第一の再構成画像の第一のフーリエ変換データに、前記第三のデータセットに対応する角度を表すマスクを適用して、前記第一の部分データセットを生成し、
   前記第二のデータセットに基づく第二の再構成画像の第二のフーリエ変換データに、前記マスクを適用して、前記第二の部分データセットを生成する、
   医用画像処理装置。
2. 前記生成部は、前記動きデータとしてワープフィールドを生成する請求項１記載の医用画像処理装置。
3. 前記再構成部は、前記動きデータと、前記第一のデータセット又は前記第二のデータセットを用いて、前記再構成処理を実行する請求項１又は２記載の医用画像処理装置。
4. 前記再構成処理は、前記動きデータに依存する請求項１乃至３のうちいずれか一項記載の医用画像処理装置。
5. 前記再構成部は、
   前記第一のデータセットを前記第一のスキャニング期間に対応する回転角度を複数に分割して得られる複数の第一の分割データセットをそれぞれ部分再構成して複数の第一の部分再構成画像を再構成し、前記動きデータから前記複数の第一の部分再構成画像に関する複数の第一の分割動きデータを生成し、前記複数の第一の部分再構成画像と前記複数の第一の分割動きデータとに基づいて第一の最終画像を再構成する、又は、
   前記第二のデータセットを前記第二のスキャニング期間に対応する回転角度を複数に分割して得られる複数の第二の分割データセットをそれぞれ部分再構成して複数の第二の部分再構成画像を再構成し、前記動きデータから前記複数の第二の部分再構成画像に関する複数の第二の分割動きデータを生成し、前記複数の第二の部分再構成画像と前記複数の第二の分割動きデータとに基づいて第二の最終画像を再構成する、
   請求項１又は２記載の医用画像処理装置。
6. 前記取得部は、
   前記第一のスキャニング期間に対応するＲＡＷデータに基づく再構成画像を用いたフーリエ変換、投影の実行、高い異方性カーネルによる畳み込みのいずれかにより前記第一のデータセットを取得し、
   前記第二のスキャニング期間に対応するＲＡＷデータに基づく再構成画像を用いたフーリエ変換、投影の実行、高い異方性カーネルによる畳み込みのいずれかにより前記第二のデータセットを取得する、
   請求項１乃至５のうちいずれか一項記載の医用画像処理装置。
7. 前記取得部は、
   前記第一のスキャニング期間に対応するＲＡＷデータに基づく再構成画像を順投影することで前記第一のデータセットを取得し、
   前記第二のスキャニング期間に対応するＲＡＷデータに基づく再構成画像を順投影することで前記第二のデータセットを取得する、
   請求項１乃至５のうちいずれか一項記載の医用画像処理装置。
8. 前記第一のデータセットと前記第二のデータセットとは、時間的に連続するスキャン処理により取得される請求項１乃至７のうちいずれか一項記載の医用画像処理装置。
9. 前記第一のスキャニング期間及び前記第二のスキャニング期間は、前記回転に関して、少なくとも１８０°の回転の回転角度に対応する請求項１乃至８のうちいずれか一項記載の医用画像処理装置。
10. 前記第一のスキャニング期間及び前記第二のスキャニング期間は、前記回転に関して、一回転の期間に含まれる請求項１乃至９のうちいずれか一項記載の医用画像処理装置。
11. 前記第一のデータセットは、前記回転に関する第一の角度範囲に対応し、
    前記第二のデータセットは、前記回転に関する第二の角度範囲に対応し、
    前記第一の角度範囲及び前記第二の角度範囲は、同じサイズであること、一方に対して１８０°分オフセットされること、の少なくともいずれかを満たす請求項１乃至１０のうちいずれか一項記載の医用画像処理装置。
12. 前記第一のデータセット及び前記第二のデータセットは、投影データ、サイノグラムデータのうちの少なくとも一つを具備する請求項１乃至１１のうちいずれか一項記載の医用画像処理装置。
13. コンピュータに、
    Ｘ線管と被検体を介して前記Ｘ線管に対向する検出器とを、前記被検体の回りに相対的に回転させることでスキャン処理を実行するスキャン部を用いて、第一のスキャニング期間において収集された第一のデータセットと、前記第一のスキャニング期間と重なり合う第二のスキャニング期間に収集された第二のデータセットと、を取得させる取得機能と、
    前記第一のスキャニング期間と前記第二のスキャニング期間とが重なる期間の少なくとも一部を含む期間において収集された第三のデータセットを前記第一のデータセットから除外して第一の部分データセットを生成させ、前記第三のデータセットを前記第二のデータセットから除外して第二の部分データセットを生成させ、前記第一の部分データセットと前記第二の部分データセットとを用いて、前記被検体の動きデータを生成させる生成機能と、
    前記動き情データと、前記第一のデータセットの少なくとも一部又は前記第二のデータセットの少なくとも一部と、を用いて、再構成処理を実行させる再構成機能と、
    を実現させ、
    前記生成機能は、
    前記第一のデータセットに基づく再構成画像の第一のフーリエ変換データに、前記第三のデータセットに対応する角度を表すマスクを適用して、前記第一の部分データセットを生成し、
    前記第二のデータセットに基づく再構成画像の第二のフーリエ変換データに、前記マスクを適用して、前記第二の部分データセットを生成する、
    医用画像処理プログラム。

Images