(57)【要約】  動きのある部位を対象として撮影を行い、診断目的に応じた断層像を再構成するX線CT装置において、範囲を絞り込んで動き情報を解析することにより、診断目的とする部位や位相での運動状態の解析精度を向上する為に、X線CT装置は、投影データとともに被検体の心電情報を取得し、投影データ及び心電情報に基づいて診断部位の動き情報を解析し、動き情報に基づいて再構成に適した運動状態(例えば、最も動き量の小さい状態)となる最適位相を決定し、最適位相にて断層像を再構成する。この際、着目する位相または部位において正確に最適位相を検出するため、動きの解析範囲を調整する。これにより、不都合な動き情報が解析対象に含まれないようにする。例えば、動き情報に対して所定の位相範囲や体軸方向範囲で適切な重みを重畳することにより解析範囲を絞り込む。
2013187461.tif 000002
【特許請求の範囲】
【請求項１】
動きのある診断部位を対象とし、被検体周囲の複数の方向からX線を照射して投影データを取得し、前記投影データに基づき被検体の断層像を再構成し、表示するX線CT装置であって、
前記診断部位の動き情報を算出する動き情報算出部と、
所定範囲の動き情報を解析対象とするための解析範囲重みを算出する解析範囲算出部と、
前記動き情報算出部により算出された動き情報に対して前記解析範囲重みを乗じることにより得た動き情報に基づいて、再構成に適した位相である最適位相を決定する最適位相決定部と、
前記最適位相における断層像を再構成する画像再構成部と、
を具備することを特徴とするX線CT装置。
【請求項２】
前記解析範囲重みは、任意の位相を基準とする位相方向の重み関数であることを特徴とする請求項1に記載のX線CT装置。
【請求項３】
被検体の心電情報を取得する心電情報取得部と、
心拍数毎の標準的な静止位相を格納する心拍数基準標準静止位相データベースと、を更に備え、
前記解析範囲算出部は、前記心電情報取得部により取得した心拍数に対応する前記静止位相を前記心拍数基準標準静止位相データベースから取得し、取得した静止位相を基準とする位相方向の重み関数を前記解析範囲重みとして算出することを特徴とする請求項2に記載のX線CT装置。
【請求項４】
撮影対象臓器について、標準的な静止位相を診断対象部位毎に格納する部位基準標準静止位相データベースを更に備え、
前記解析範囲算出部は、診断対象とする部位に対応する前記静止位相を前記部位基準標準静止位相データベースから取得し、取得した静止位相を基準とする位相方向の重み関数を前記解析範囲重みとして算出することを特徴とする請求項2に記載のX線CT装置。
【請求項５】
前記解析範囲重みは、前記投影データ計測時におけるX線強度に応じた形状とすることを特徴とする請求項2に記載のX線CT装置。
【請求項６】
前記最適位相は、最も動き量の小さい位相であることを特徴とする請求項1に記載のX線CT装置。
【請求項７】
前記解析範囲重みは、任意の撮影部位を基準とする体軸方向の重み関数であることを特徴とする請求項1に記載のX線CT装置。
【請求項８】
被検体の心電情報を取得する心電情報取得部を更に備え、
前記解析範囲重みは、前記心電情報取得部により取得した心電情報の変動が小さい撮影部位において重みの値が大きくなる形状とすることを特徴とする請求項7に記載のX線CT装置。
【請求項９】
前記解析範囲重みは、診断の重要性が高い撮影部位において重みの値が大きくなる形状とすることを特徴とする請求項7に記載のX線CT装置。
【請求項１０】
前記解析範囲重みは、静止位相算出の妨げになる構造物を含む撮影部位の重みの値が小さくなる形状とすることを特徴とする請求項7に記載のX線CT装置。
【請求項１１】
前記解析範囲重みは、心拍位相方向及び体軸方向を含む多次元の重み情報であることを特徴とする請求項1に記載のX線CT装置。
【請求項１２】
前記解析範囲重みを決定するためのパラメータを入力する入力部を更に備え、
前記解析範囲算出部は、前記入力部により入力されたパラメータに基づいて前記解析範囲重みを決定することを特徴とする請求項1に記載のX線CT装置。
【請求項１３】
前記動きのある診断部位は心臓であって、
前記入力部は、心臓の1心拍に相当する心電情報と、前記心電情報に重畳され、前記解析範囲重みに対応するオブジェクトとを表示する表示部を備えることを特徴とする請求項12に記載のX線CT装置。
【請求項１４】
前記入力部は、前記動きのある診断部位の評価用画像と、前記評価用画像の体軸方向位置に対応する前記解析範囲重みを設定するために用いられるオブジェクトとを表示する表示部を備えることを特徴とする請求項12に記載のX線CT装置。
【請求項１５】
前記入力部は、前記動きのある診断部位の体軸方向位置を示す軸と位相を示す軸とを有するエリアであって、各軸に対応する複数の画像が並べられたエリアと、前記体軸方向位置に対応する前記解析範囲重みを設定するために用いられるオブジェクトとを表示する表示部を備えることを特徴とする請求項12に記載のX線CT装置。
【請求項１６】
前記入力部は、前記動きのある診断部位の体軸方向位置を示す軸と位相を示す軸とを有するエリアであって、各軸に対応する複数の画像が並べられたエリアを表示する表示部を備え、前記エリア中の画像を選択することにより前記解析範囲重みを設定することを特徴とする請求項12に記載のX線CT装置。
【請求項１７】
動きのある診断部位を対象とし、被検体周囲の複数の方向からX線を照射して投影データを取得し、前記投影データに基づき被検体の断層像を再構成する画像再構成方法であって、
前記診断部位の動き情報を算出する動き情報算出ステップと、
所定範囲の動き情報を解析対象とするための解析範囲重みを算出する解析範囲算出ステップと、
前記動き情報算出ステップにより算出した動き情報に対して前記解析範囲重みを乗じることにより得た動き情報に基づいて、再構成に適した位相である最適位相を決定する最適位相決定ステップと、
前記最適位相における断層像を再構成する画像再構成ステップと、
を含むことを特徴とする画像再構成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動きのある部位の撮影に好適なX線CT装置及び画像再構成方法に関する。
【背景技術】
【０００２】
X線CT(Computed Tomography)装置によって身体の動きのある部位を撮影する場合、得られる断層像には動きに起因するアーチファクトが発生する。このアーチファクトを低減するためには、一般的には心電計や呼吸センサ等の生体センサを用いて生理的運動の計測を合わせて行い、得られた計測信号を用いて撮影を制御したり、画像を処理したりする。例えば、心臓を対象とする撮影では、心電計によって心電情報を計測しながら投影データを収集し、後から適切な心拍位相を選択して断層像を再構成する心電同期再構成法が行われる。心電同期再構成法の具体例としては、1心拍の投影データから1枚の断層像を再構成する心電同期ハーフ再構成法や、複数心拍の投影データを組み合わせて1枚の断層像を再構成する心電同期セグメント再構成法等が知られている。いずれの手法であっても、心臓の周期的な運動を解析し、最も動きの少ない心位相(以下、静止位相という)を求め、静止位相での投影データを用いて診断用の画像を再構成するため、動きの解析精度を向上する必要がある。【０００３】
静止位相の決定方法は、例えば、特許文献1に示す技術が開示されている。特許文献1には、心電情報とともに取得した投影データに基づいて、R波を基準とした相対位置で2％の間隔で各位相n(0％、2％、4％、・・・98％、100％)の断層像Tnを再構成し、再構成した各画像Tnについて、その画像Tnと2位相前の画像Tn-2とを差分して差分画像Xnを生成し、差分画像Xnの全画素または関心部位の絶対値総和や二乗和を心臓の動き量を示す指標値として求め、求めた動き量の時間変化の中から、例えば最も動き量の小さい位相を選択することが記載されている。また、動き量の時間変化を表示装置に表示して、操作者に手動で最適な心拍位相を指定させることが記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開2007-37782号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、心臓の動き量を示す曲線は撮影中の全時相及び全部位にわたって必ずしも理想的な曲線を得られるわけではなく、実際には動きの解析に不都合な部分を含むことがある。例えば、撮影中に被検体の体動が生じたり、心電計が不具合を起こしたり、不整脈が生じたりする場合には、本来の心臓の動きとは異なる動き量が含まれるため、このような不規則な動きの部分を除外して静止心位相を決定すべきである。また、例えばステントや石灰化といった高吸収体(大きなCT値となる物体)が撮影部位に含まれている場合に、上述のように差分画像の画素値の絶対値総和等を動き量とすると、動き量が本来の動きよりも大きく計算されてしまい、心運動解析の妨げになる。そのため、上述の特許文献1のように、心位相全体を対象として心臓の動き量を解析すると、最終的に算出される最適心拍位相の信頼性が低下する恐れがあった。【０００６】
本発明は、前述した問題点に鑑みてなされたものであり、その目的とすることは、動きのある部位を対象として撮影を行い、画像を再構成するX線CT装置において、範囲を絞り込んで動き情報を解析することにより、診断目的とする部位や位相での運動状態の解析精度を向上することが可能なX線CT装置及び画像再構成方法を提供することである。
【課題を解決するための手段】
【０００７】
前述した目的を達成するために第1の発明は、動きのある診断部位を対象とし、被検体周囲の複数の方向からX線を照射して投影データを取得し、前記投影データに基づき被検体の断層像を再構成し、表示するX線CT装置であって、前記診断部位の動き情報を算出する動き情報算出部と、所定範囲の動き情報を解析対象とするための解析範囲重みを算出する解析範囲算出部と、前記動き情報算出部により算出された動き情報に対して前記解析範囲重みを乗じることにより得た動き情報に基づいて、再構成に適した位相である最適位相を決定する最適位相決定部と、前記最適位相における断層像を再構成する画像再構成部と、を具備することを特徴とするX線CT装置である。【０００８】
また、第2の発明は、動きのある診断部位を対象とし、被検体周囲の複数の方向からX線を照射して投影データを取得し、前記投影データに基づき被検体の断層像を再構成する画像再構成方法であって、前記診断部位の動き情報を算出する動き情報算出ステップと、所定範囲の動き情報を解析対象とするための解析範囲重みを算出する解析範囲算出ステップと、前記動き情報算出ステップにより算出した動き情報に対して前記解析範囲重みを乗じることにより得た動き情報に基づいて、再構成に適した位相である最適位相を決定する最適位相決定ステップと、前記最適位相における断層像を再構成する画像再構成ステップと、を含むことを特徴とする画像再構成方法である。
【発明の効果】
【０００９】
本発明により、動きのある部位を対象として撮影を行い、断層像を再構成するX線CT装置において、範囲を絞り込んで動き情報を解析することにより、診断目的とする部位や位相での運動状態の解析精度を向上することが可能なX線CT装置及び画像再構成方法を提供できる。
【図面の簡単な説明】
【００１０】
【図１】X線CT装置1の全体構成図
【図２】本発明に係る画像再構成処理の流れを説明するフローチャート
【図３】心電同期セグメント再構成法について説明する図
【図４】動き情報の算出方法について説明する図。(a)はデータセット、(b)は体軸方向位置毎の動き量曲線(動き情報)、(c)は全体動き量曲線
【図５】位相方向解析範囲重みの例。(a)は心電波形、(b)は動き量曲線を示すグラフ、(c)、(d)は位相方向解析範囲重みの例
【図６】位相方向解析範囲重みを心拍数に応じて調整する例。(a)は低心拍の心電波形、(b)は低心拍の場合の位相方向解析範囲重み、(c)は高心拍の心電波形、(d)は高心拍の場合の位相方向解析範囲重み
【図７】心拍数基準標準静止位相データベースの一例
【図８】部位基準標準静止位相データベースの一例
【図９】位相方向解析範囲重みをX線出力分布に応じて調整する例。(a)は心電波形、(b)は心電情報に同期したX線出力分布の一例、(c)はX線出力に応じた位相方向解析範囲重みの一例
【図１０】体軸方向解析範囲重みの例。(a)は心電波形、(b)はデータセット、(c)は動き情報、(d)は体軸方向解析範囲重み、(e)は全体動き量曲線
【図１１】評価用画像に基づいて設定される体軸方向解析範囲重みの一例。(a)は評価用画像、(b)は動き情報、(c)は体軸方向解析範囲重み、(d)は全体動き量曲線
【図１２】評価用画像に基づいて体軸方向解析範囲重みを算出する場合の画像再構成処理の流れを説明するフローチャート
【図１３】心拍位相方向及び体軸方向を含む多次元の解析範囲重み(重みマップ)について説明する図。(a)は動き情報、(b)は重みマップ
【図１４】解析範囲設定画面の一例
【図１５】解析範囲設定画面の別の例
【図１６】操作者に解析用画像を提示した後に、解析範囲重みを算出する場合の画像再構成処理の流れを説明するフローチャート
【図１７】解析用画像データセットを表示した解析範囲設定画面の例
【図１８】解析用画像データセットを表示した解析範囲設定画面の別の例
【発明を実施するための形態】
【００１１】
以下、添付図面を参照しながら、本発明の好適な実施形態について詳細に説明する。
まず、図1を参照して、X線CT装置1の構成について説明する。
【００１２】
X線CT装置1は、被検体6に対するX線照射を行うとともに被検体6を透過したX線を検出するスキャナガントリ2と、被検体6を載置する寝台テーブル3と、スキャナガントリ2の回転動作及びX線照射を制御する測定制御装置202と、被検体6の心電情報を取得する心電計7と、X線CT装置1の各部を制御するための操作卓4とを備える。操作卓4は、システム制御装置401、画像処理装置403、表示装置405、及び入力装置406を備える。
【００１３】
スキャナガントリ2には、X線管201及びコリメータ203と、X線検出器205とが、回転盤207の開口部を介して対向配置されている。回転盤207のX線照射空間である開口部が設けられ、被検体6が載置される寝台テーブル3が搬入される。回転盤207は、測定制御装置202によって制御される回転盤駆動装置210から駆動伝達系を通じて伝達される駆動力によって被検体6の周囲を周回するよう駆動される。
【００１４】
X線管201はX線源であり、測定制御装置202に制御されて所定の強度のX線を連続的または断続的に照射する。測定制御装置202は、操作卓4のシステム制御装置401により決定されたX線管電圧及びX線管電流に従って、X線管201に印加または供給するX線管電圧及びX線管電流を制御する。コリメータ203は、X線管201から放射されたX線を、例えばコーンビーム(円錐形または角錐形ビーム)等のX線として被検体6に照射させるものであり、コリメータ203の開口幅は測定制御装置202により制御される。被検体6を透過したX線はX線検出器205に入射する。
【００１５】
X線検出器205は、例えばシンチレータとフォトダイオードの組み合わせによって構成されるX線検出素子群をチャネル方向(周回方向)に例えば1000個程度、列方向(体軸方向)に例えば1〜320個程度配列したものであり、被検体6を介してX線管201に対向するように配置される。X線検出器205はX線管201から放射されて被検体6を透過したX線を検出し、検出した透過X線データを図示しないデータ収集装置に出力する。データ収集装置は、X線検出器205の個々のX線検出素子により検出される透過X線データを収集し、ディジタルデータに変換し、投影データとして操作卓4の画像処理装置403に順次出力する。
【００１６】
測定制御装置202は、操作卓4のシステム制御装置401からの制御信号に従って、スキャナガントリ2内のX線源201、コリメータ203、回転板207の回転を制御する。
【００１７】
寝台テーブル3は、操作卓4のシステム制御装置401から送出される制御信号に従って、寝台テーブル3の高さを適切なものにし、また、体軸方向への前後動、体軸と垂直方向であって天板に平行な方向(左右方向)に移動する。これにより、被検体6がスキャナガントリ2の開口部(X線照射空間)に搬入及び搬出される。
【００１８】
操作卓4のシステム制御装置401は、CPU(Central Processing Unit)、ROM(Read Only Memory)、RAM(Random Access Memory)、ハードディスク等の記憶部等により構成されるコンピュータである。システム制御装置401は、スキャナガントリ2(測定制御装置202)、寝台テーブル3、及び心電計7を制御する。
【００１９】
システム制御装置401の記憶部には、画像処理装置403が生成する画像やX線CT装置1の機能を実現するためのプログラム等を記憶する。また、記憶部には後述する画像再構成処理において診断部位の動きの解析対象とする範囲を算出する際に使用される解析範囲重みの基となるデータや各種データベースが格納される。
【００２０】
システム制御装置401はスキャナガントリ2のデータ収集装置から入力された投影データと心電計7から入力された心電情報とを、撮影時刻に基づいて対応付ける。
【００２１】
画像処理装置403は、取得した投影データに対して、対数変換、感度補正等の前処理を行った後、断層像を再構成する。
【００２２】
断層像の再構成処理において、画像処理装置403は取得した投影データ及び心電情報に基づいて、まず、同一心拍位相で撮影位置の異なる断層像のデータセットを動き解析用の画像として再構成する。そして、異なる位相間で上述の画像を比較し、診断部位の動き情報を解析する。更に、解析した動き情報に基づいて再構成に適した運動状態(例えば、最も動き量の小さい状態；静止位相)となる最適位相を決定し、最適位相にて断層像を再構成する心電同期再構成処理を行う。
【００２３】
本発明において、画像処理装置403は、着目する位相または部位において正確に最適位相を検出するため、動きの解析範囲を調整する(解析範囲算出部)。これは、撮影中の体動や診断部位に含まれる構造物等の影響等による不都合な動き情報が解析対象に含まれないようにするための処理である。このため、画像処理装置403は、例えば、動き情報に対して所定の位相範囲や体軸方向範囲で適切な重み(以下、解析範囲重みという)を重畳する(動き情報算出部)。このように範囲を限定して求めた動き情報に基づいて最適位相を決定する(最適位相決定部)。
解析範囲重みの具体例については、後述する。
【００２４】
表示装置405は、液晶パネル、CRTモニタ等のディスプレイ装置と、ディスプレイ装置と連携して表示処理を実行するための論理回路で構成され、システム制御装置401に接続される。表示装置405は画像処理装置403から出力される再構成画像、並びにシステム制御装置401が取り扱う種々の情報を表示するものであり、操作者が表示内容を閲覧する。
【００２５】
入力装置406は、例えば、キーボード、マウス等のポインティングデバイス、テンキー、及び各種スイッチボタン等により構成され、操作者によって入力される各種の指示や情報をシステム制御装置401に出力する。操作者は、表示装置405及び入力装置406を使用して対話的にX線CT装置1を操作する。
【００２６】
心電計7は、被検体6に取り付けた電極を介して、心臓の心拍運動を反映した活動電位の時間変化を表す心電情報を計測し、例えば0.1秒間隔等の所定のサンプリングピッチでディジタル信号に変換する。心電計7により計測した心電情報は、システム制御装置401へ順次送出され、システム制御装置401により投影データに付加される。
【００２７】
次に、図2を参照して、X線CT装置1の動作について説明する。
本実施の形態のX線CT装置1は、図2のフローチャートに示す手順で画像再構成処理を実行する。すなわち、システム制御装置401のCPUは、記憶部から図2に示す画像再構成処理に関するプログラム及びデータを読み出し、このプログラム及びデータに基づいて処理を実行する。
【００２８】
まず、X線CT装置1は、操作者の操作に従って、動きのある部位の撮影を行う。以下の説明では、心臓を対象とする撮影を行うものとする。このとき、心電計7を用いて心電情報の取得も同時に行う(ステップS101)。システム制御装置401は、スキャナガントリ2から取得した投影データ及び心電計7から取得した心電情報を画像処理装置403へ送出する。
【００２９】
次に、画像処理装置403は、動きの解析範囲を算出する。これは、上述したように後段の処理(ステップS105の最適位相決定処理)において、着目する位相または部位において発生した不都合な動き情報の混入を防ぎ、正確に最適位相を検出するための処理である。具体的には、画像処理装置403は、ステップS104で求める動き情報に重畳する解析範囲重みを決定する(ステップS102)。解析範囲重みとは、例えば図5(c)、(d)等に示すように、位相毎に「0.0」〜「1.0」の値を持つ重み関数である。
解析範囲重みの詳細については後述する。
【００３０】
次に、X線CT装置1の画像処理装置403は、例えば、周知の心電同期セグメント再構成法等の手法を用いて、動き情報を解析するための解析用画像を再構成し、心拍位相の異なる複数のデータセットを作成する(ステップS103)。データセットとは、図4(a)に示すように、心拍位相(P1、P2、・・・Pn)毎に連続した体軸方向位置(z0、z1、・・・zx)の画像をまとめた断層像群のことである。図4では、心拍位相5％刻みで0％から100％(全位相範囲)のデータセットが作成された例を示しているが、本発明では全位相範囲のデータセットを作成する必要はなく、ステップS102の処理によって調整された解析範囲についてのみデータセットを作成すればよい。例えば、図5(c)、(d)に示すように、重みの値が「0」となる範囲については必ずしもデータセットを作成する必要はない。
【００３１】
ステップS103の解析用画像再構成処理において、1枚の断層像を作成するためには、少なくともビュー角度180°＋α(αは扇状X線のファン角度)分の投影データが必要である。例えば、心電同期セグメント再構成法では、取得した投影データを分割し、この分割投影データのうち複数の心拍から同じ心拍位相(例えば、R波の相対位置で80％位置)で得られたものを組み合わせて、180°＋α分の投影データを収集する。
【００３２】
図3は、心拍位相80％(R波の相対位置80％)の断層像を再構成するために必要な投影データの収集の仕方を説明する図である。
【００３３】
図3に示すように、画像処理装置403は、時相R0〜R1、R1〜R2、R2〜R3の3心拍中に撮影された投影データから同じ位相(例えば、R波の相対位置80％)で撮影角度の異なる分割投影データを180°＋αだけ収集する。図3の例では、時相R0〜R1間の投影データから角度0π〜1／3πの分割投影データ、時相R1〜R2間の投影データから角度1／3π〜2／3πの分割投影データ、時相R2〜R3間の投影データから角度2／3π〜πの分割投影データを収集している。このようにして収集された約180°分の投影データを用いて、解析用画像(断層像)を再構成する。この解析用画像を、各体軸方向位置(z0、z1、・・・zx)、及び各心拍位相(P1〜Pn、例えば、位相5％刻み)で再構成し、心拍位相毎に複数のデータセットを作成する。
【００３４】
解析用画像のデータセットが生成されると、次に画像処理装置403は、作成されたデータセットに基づいて動き情報を算出する(ステップS104)。
【００３５】
説明のために、まず、従来の全データ範囲を対象とする動き情報の算出方法を図4を参照して説明する。
【００３６】
画像処理装置403は、同一体軸位置で心拍位相が異なる複数の断層像を対象として、R波間隔の相対位置0％から100％の各位相P1〜Pnについて動き情報を算出する。
以下の説明では、体軸方向位置z1で算出された動き情報をt1、体軸方向位置z2で算出された動き情報をt2、・・・のように表記する(図4(b)参照)。
【００３７】
動き情報の算出方法は、様々な方法が知られているがいずれを用いてもよい。例えば、対象とする位相Piの画像と数％(例えば、3〜6％程度)異なる位相Pｊの画像とを差分し、その差分画像の各画素の絶対値総和や二乗和等を動き量として求め、動き量の位相方向の変化を示した動き量曲線を動き情報txとする。この処理を診断部位の全ての体軸方向位置(z1〜zx)で行うことにより、図4(b)の動き情報t1、t2、・・・txを算出する。そして、各体軸方向位置(z1〜zx)で得られた全ての動き情報t1〜txを体軸方向に加算し、診断対象臓器全体の動き量曲線T(以下、動き情報Tまたは全体動き情報Tという)とする。
【００３８】
そして、画像処理装置403は、全体動き情報Tに基づいて、診断用画像を再構成する所望の運動状態の位相(最適位相)を決定し(ステップS105)、決定した最適位相の断層像を再構成する(ステップS106)。最適位相としては、一般には、動きアーチファクトの影響を避けるために最も動き量の小さい位相(以下、静止位相という)を選択するが、必ずしもこれに限定されない。
【００３９】
また、ここで選択する位相は1つに限定せず、複数選択するようにしてもよい。例えば、心臓の収縮期と拡張期の各心拍位相範囲について、それぞれ最も動きの小さい位相(静止位相)を求め、各静止位相の診断用画像を再構成するようにしてもよい。
【００４０】
画像処理装置は、ステップS106で再構成した診断用画像を表示装置405に表示する(ステップS107)。診断用画像とは、診断目的に応じた再構成条件で作成された断層像や3次元画像等である。
【００４１】
以上が全データ範囲を対象として動き情報を求めた場合の画像再構成処理の全体の流れであるが、本発明では、ステップS104の動き情報算出処理において、ステップS102の解析範囲算出処理において求めた解析範囲重みを動き情報に重畳することにより、動きの解析対象とする範囲を調整する。これにより、信頼性の低い動き情報を除いて全体の動き情報を算出することが可能となる。
【００４２】
以下、解析範囲重みの例を説明する。
【実施例１】
【００４３】
動き情報の信頼性低下の要因として、撮影中に被検体の体動が生じた場合や、特定の心拍位相範囲にのみ影響が表れる構造物が撮影範囲に含まれる場合等が挙げられる。被検体に体動が生じると動き情報に非周期的な動き成分が混入する。また、特定の心拍位相範囲にのみ影響が表れる構造物が撮影範囲に含まれる場合には、その特定の心拍位相範囲における動き情報の値は標準的な曲線から大きく変動することがある。【００４４】
この問題を解決するためには、着目する心拍位相範囲において重みの値が大きくなり、それ以外の心拍位相範囲では重みの値が小さくなるような解析範囲重みを与える。以下、心拍位相方向にかける重み関数を「位相方向解析範囲重み」と呼ぶ。【００４５】
具体的には、ステップS102において画像処理装置403は、図5(c)、(d)に示すように、特定の位相を中心とし、位相方向に所定幅を有する形状の位相方向解析範囲重みを算出する。そして、ステップS104の動き情報算出処理において、動き情報に位相方向解析範囲重みを乗じ、全体動き情報Tを求める。【００４６】
図5(a)は心電波形、図5(b)は動き量曲線を示すグラフ、図5(c)、(d)は位相方向解析範囲重みの例を示す。【００４７】
例えば、着目する心拍位相範囲が70％である場合、図5(c)のように、心拍位相70％を中心として所定幅を持つ矩形の重み関数を位相方向解析範囲重みとして設定する。
図5(c)の例の重み関数は、心拍位相70％付近で「1.0」、その他の位相で「0」の値を持つ矩形形状のものであるが、重みの値はこれに限定されない。【００４８】
また、例えば、図5(d)に示すように、心拍位相70％をピークとする三角形状の重み関数を位相方向解析範囲重みとして設定するようにしてもよい。この場合、重み値の切り替わりが動き情報に与える影響を低減できる。【００４９】
矩形や三角の幅は、予め設定された標準的な位相範囲としてもよいし、操作者が適宜編集するようにしてもよい。【００５０】
なお、重み関数を動き情報に乗じるタイミングは以下の2通りのうち、いずれでもよい。【００５１】
(i)体軸方向位置(z1、z2、・・・zx)毎に求めた動き情報t1、t2、・・・txにステップS102で設定した重み関数をそれぞれ乗じ、それらを加算して体軸方向に加算して全体動き情報Tを算出する。【００５２】
(ii)体軸方向位置(z1、z2、・・・zx)毎に求められた動き情報t1、t2、・・・txを体軸方向に加算して求めた全体動き情報Tに、ステップS102で設定した重み関数を乗じる。【００５３】
(i)の場合は、体軸方向位置毎に解析範囲重みを乗じるため、精度よく動き情報を解析できるという効果がある。(ii)の場合は、全体動き情報Tに対して一度だけ解析範囲重みを乗じればよいため演算量を大幅に低減でき、処理を高速化できるという効果がある。【００５４】
(i)、(ii)のいずれの場合であっても、重みの値が「0」となる位相範囲についてデータセット(解析用画像)を作成する必要はないため、演算量を低減できる。【００５５】
動き情報に位相方向解析範囲重みをかけることにより、診断目的とする位相範囲に限定して最終的な動き情報Tを算出できる。このため、診断目的としない位相で生じた標準的でない動きの混入を防ぐことができ、正確に動き情報を算出することができる。そのため、動き情報の信頼性を向上でき、最適位相決定の精度が向上する。また、解析範囲を限定できるため、演算量を削減できる。また、図5の例の重み関数のように、特定の位相を中心とする矩形または三角形状の重みをかけることにより、特に、診断目的とする位相が予めわかっている場合等に、動きの解析範囲を絞り込むことができ、動き情報の信頼性を向上できる。
【実施例２】
【００５６】
心臓の運動状態は拍動周期の影響を受けることが知られている。そのため、静止心拍位相は被検体の心拍数によって異なる。【００５７】
そのため、図5に示したような心拍位相方向の解析範囲重みを設定する場合、撮影中の心拍数の変動の影響を避けるために、被検体の心拍数に応じて重み関数の形状を変更するようにしてもよい。【００５８】
図6は、心拍数に応じた位相方向解析範囲重みの調整の一例を示しており、(a)は低心拍の心電波形、(b)は低心拍の場合の位相方向解析範囲重み、(c)は高心拍の心電波形、(d)は高心拍の場合の位相方向解析範囲重みを示している。【００５９】
ステップS102の解析範囲算出処理において、画像処理装置403は、心電情報に応じて、図6に示すように、低心拍と高心拍とで重み関数の形状の中心位置を調整する。例えば、低心拍では心拍位相30％を中心とした所定幅を持つ矩形形状の重み関数とし、高心拍では心拍位相40％を中心とした所定幅を持つ矩形形状の重み関数とする。【００６０】
このように心拍数に応じて解析範囲重みを調整する場合、システム制御装置401の記憶部に心拍数と重み関数の形状の組み合わせを定義する心拍数基準標準静止位相データベース81を格納しておき、画像処理装置401は、取得した心電情報に基づき心拍数に対応する重み関数の形状を、心拍数基準標準静止位相データベース81を参照して決定するようにしてもよい。【００６１】
図7は、心拍数と重み関数の形状の組み合わせを定義した心拍数基準標準静止位相データベース81の一例である。【００６２】
図7の心拍数基準標準静止位相データベース81では、心拍数範囲毎に、収縮期と拡張期とでそれぞれ重み関数の矩形の先頭及び末尾の位相位置(％)を定義している。【００６３】
例えば、心拍数が40〜49［bpm］では、収縮期は20〜40［％］、拡張期は60〜80［％］の矩形重みとし、心拍数が50〜59［bpm］では、収縮期は25〜45［％］、拡張期は62〜82［％］の矩形重みとし、心拍数が60〜69［bpm］では、収縮期は30〜35［％］、拡張期は64〜84［％］の矩形重みとする。【００６４】
なお、この心拍数基準標準静止位相データベース81は、操作者が自在に編集できるようにすることが望ましい。【００６５】
また、静止心拍位相は臓器やその位置によっても異なることが一般的に知られている。
例えば、心臓の検査では、診断対象とする部位として、RCA(Right Coronary Artery；右冠動脈)、LAD(Left Descending Coronary Artery；左冠動脈前下行枝)、LCX(Left circumflex Coronary Artery；左冠動脈回旋枝)等があり、各部位でそれぞれ収縮期や拡張期の静止位相が異なっている。【００６６】
そこで、システム制御装置401の記憶部に診断対象となる部位と重み関数の形状の組み合わせを定義する部位基準標準静止位相データベース82を格納しておき、画像処理装置403は、診断目的とする部位に対応する重み関数の形状を、部位基準標準静止位相データベース82を参照して決定するようにしてもよい。【００６７】
図8は、診断部位と重み関数の形状の組み合わせを定義した部位基準標準静止位相データベース82の一例である。【００６８】
図8の部位基準標準静止位相データベース82では、診断部位毎に、収縮期と拡張期とでそれぞれ重み関数の矩形の先頭及び末尾の値(％)を定義している。【００６９】
例えば、診断部位「RCA：近位」では、収縮期は20〜40［％］、拡張期は60〜80［％］の矩形重みとし、診断部位「RCA：遠位」では、収縮期は25〜45［％］、拡張期は62〜82［％］の矩形重みとし、診断部位「LAD：近位」では、収縮期は30〜35［％］、拡張期は64〜84［％］の矩形重みとする。【００７０】
なお、図8の部位基準標準静止位相データベース82についても、操作者が編集できるようにすることが望ましい。
【実施例３】
【００７１】
心電情報を利用した心臓撮影において、撮影中の心電情報に同期し、X線出力を調整する手法がある。この手法は患者の被ばく線量の低減を目的とするものであり、例えば、診断対象となる心拍位相のみに限定してX線を曝射する手法や、診断対象となる心拍位相以外の領域のX線出力を抑制する手法等がある。後者の手法の場合、X線出力の小さい心拍位相で撮影された画像は画質が悪く、診断への有用性は低い。【００７２】
そこで、位相方向解析範囲重みをX線出力の分布に応じて調整するようにしても良い。【００７３】
図9にX線出力分布に応じた解析範囲重みの一例を示す。図9(a)は心電波形、(b)は心電情報に同期したX線出力分布の一例、(c)はX線出力に応じた位相方向解析範囲重みの一例である。【００７４】
図9(b)に示すように、X線出力が、診断目的とする拡張位相で100％、それ以外の位相で20％となるように抑制されている場合、位相方向解析範囲重みについても、図9(c)に示すように、X線出力に連動した重み形状となるようにする。【００７５】
これにより、診断への有用性の高い部位についての動き情報に限定して動きを解析することができるため、信頼性の高い動き情報を算出することができる。【００７６】
(実施例1)〜(実施例3)は、位相方向に乗じる解析範囲重みについて説明したが、解析範囲重みは、体軸方向に適用するようにしてもよい。
【実施例４】
【００７７】
心電同期再構成法では、心電計から取得した心電情報を基に画像再構成を行う。このため、被検体の健康状態によっては心電情報が不規則になり、心電計の不具合によって正しい心電情報が得られない場合がある。このような不具合が発生した時の動き情報は不正な情報を含むものとなるため、最終的に算出される全体動き情報に悪影響を及ぼす可能性がある。【００７８】
この問題に対応するため、ステップS102の解析範囲算出処理において、画像処理装置403は、不正な心電情報を取得した撮影位置における重みが小さくなるような重み関数を算出する。【００７９】
図10(a)に示すように、撮影中に不整脈が発生した場合、その心電波形を参照して作成される再構成画像は期待した心拍位相の画像とならない。そのため、X線CT装置1の画像処理装置403は、不整脈が生じた撮影時刻に対応する撮影位置(体軸方向位置)において、重みの値が小さくなるような形状の重み関数を設定し、この重み関数を各体軸方向位置の動き情報に対して乗じて全体の動き情報を算出する。【００８０】
図10(c)に示す重み関数は、不整脈が生じた撮影時刻に対応する体軸方向位置での重みの値を「0.2」程度、正常な心電波形を得た体軸方向位置での重みの値を「1.0」とする台形の形状とする。なお、動き関数の形状は、台形に限定されず、矩形や三角形状としてもよい。【００８１】
このように、心電情報の一部に明らかに信頼性の低い部分が含まれる場合は、信頼性の低い範囲について重みの値が小さくなるような解析範囲重みを設定する。これによって、不正または不規則な心電波形となった撮影部位の動き情報を除いて、全体動き情報を算出することができるため、信頼性の高い動き情報を算出することができる。
【実施例５】
【００８２】
図10の例では、心電情報を参照して体軸方向の解析範囲重みを算出するものとしたが、診断上、重要な部位における体軸位置の重みが大きくなるように、解析範囲重みを算出するようにしてもよい。【００８３】
例えば、心臓領域の冠状動脈の分岐部(図11(a)の91、92)は疾患が発生する確率が高いことが一般に知られており、診断における重要性が高い。また冠状動脈の疾患である粥状腫瘍(プラーク)(図11(a)の93)を含む体軸方向位置の画像は診断における重要性が高い。【００８４】
診断においてどの部位が重要であるかは、画像処理装置403が自動検出するものとしてもよいし、操作者が画像を評価して指定するものとしてもよい。【００８５】
以下、図12のフローチャートを参照して、診断上の重要性に応じて解析範囲重みを算出する処理を含む画像再構成処理の流れを説明する。【００８６】
図12の処理フローでは、図2と比較して、解析範囲を算出する処理(図12のステップS203)の前に、評価用画像90を再構成するステップ(S202)が含まれる点が異なる。【００８７】
すなわち、システム制御装置401は、図2のフローチャートのステップS101と同様に投影データ及び心電情報を取得し(ステップS201)、画像処理装置403へ送出する。そして、画像処理装置403は、評価用画像を再構成する(ステップS202)。
この評価用画像再構成処理では、画像処理装置403は、診断の重要性を判断するために最も適切な画像を評価用画像として再構成する。例えば、心臓を診断対象とする場合には、モーションアーチファクトの少ない画像を再構成すればよい。具体的には例えば、取得した投影データ及び心電情報に基づいて一般的な心電同期再構成処理を行い、拡張中期における心拍位相の断層像を作成する。【００８８】
ステップS202において、評価用画像90を作成すると、画像処理装置403は、作成した評価用画像90に基づいて診断上の重要性に応じた体軸方向の解析範囲重みを算出する(ステップS203)。【００８９】
すなわち、画像処理装置403は、作成した評価用画像90を解析し、冠状動脈の分岐部や、狭窄、粥状腫瘍等の位置を特定する。これらの部位の特定は、CAD(Computer Aided Diagnosis)と呼ばれる公知の画像診断支援方法等を適用すればよい。また、評価用画像90を表示装置405に表示し、操作者が目視により該当する部位を指定するようにしてもよい。【００９０】
そして画像処理装置403は、特定された部位における体軸方向位置の重みが大きくなるように解析範囲重みを算出する。【００９１】
例えば、図11(a)に示すような評価用画像90が再構成された場合、診断に重要な部位91、92、93を特定すると、画像処理装置403は特定した部位91、92、93を含む体軸方向位置で重みの値が大きくなるような形状の重み関数(図11(c))を算出する。【００９２】
その後、図2のステップS103と同様の手法で動き解析用の画像を再構成し(ステップS204)、各体軸方向位置(z1、z2、・・・z8)で動き情報t1、t2、・・・t8を算出し、これらの動き情報に対してステップS203で算出した体軸方向解析範囲重みを乗じ、加算することによって、全体動き情報を算出する(ステップS205)。【００９３】
そして、図2のステップS105〜ステップS107と同様に、算出した全体動き情報Tに基づいて所望の動き状態となる最適位相を決定し(ステップS206)、最適位相における診断用画像を再構成して(ステップS207)、表示装置405に表示する(ステップS208)。【００９４】
なお、図11及び図12の説明では、診断上の重要性の高い部位について重みを大きくする例について説明したが、同様の手法で解析に悪影響を及ぼす要因を含む体軸方向位置の重みを小さくすることで動き情報の信頼性を向上させるようにしてもよい。解析に悪影響を及ぼす要因としては、例えば、冠動脈内に留置されたステントや石灰化等の高吸収体(CT値が高い物体)の存在がある。隣接する位相の画像間の差分をとって動き情報を算出する場合、高吸収体を含む領域の動き量は他の領域と比較して大きく評価される。この場合、高吸収体の動き情報の影響が大きくなり、診断対象となる部位の動きの評価が正しく行えない場合がある。【００９５】
この場合も図12と同様の処理手順(ステップS201、S202、S203)で解析範囲重みを算出する。すなわち、画像処理装置403は、評価用画像90を解析し、動きの解析に悪影響を及ぼす要因(ステントや、石灰化等)の位置を、上述したように、例えばCADと呼ばれる公知の画像診断支援方法を適用したり、評価用画像を表示装置405に表示し、操作者が目視により該当する部位を入力したりして、特定する。そして、画像処理装置403は、特定された部位における体軸方向位置の重みが小さくなるように解析範囲重みを算出する。
【実施例６】
【００９６】
(実施例1)〜(実施例5)では、心拍位相方向の解析範囲重みと、体軸方向の解析範囲重みとを個別に設定する例について説明したが、双方を組み合わせた多次元的な解析範囲重みを作成することも可能である。【００９７】
図13(b)は、横軸を心拍位相、縦軸を体軸方向とした重みマップを示す。図13において、濃淡は重みの大きさに対応しており、白色ほど重みの値が大きいことを示している。【００９８】
このような重みマップを複数の重み情報群に乗じることで最終的な動き情報を算出する。【００９９】
図13のように多次元的な重みを作成することで、着目する身体部位と心拍位相の双方に最適化した動き情報を算出することが可能となる。【０１００】
図13は2次元の重みマップを示した例であるが、身体部位を体軸方向(z方向)のみならず奥行き方向(体軸と垂直な方向)を含む3次元情報に拡張すれば、3次元の重みマップを作成することも可能である。
【実施例７】
【０１０１】
次に、解析範囲算出(解析範囲重みの算出)に関するユーザインターフェースについて説明する。【０１０２】
図2のステップS102で実施する位相方向解析範囲重みの算出処理において、X線CT装置1の表示装置405及び入力装置406を利用して、操作者が直接、重みの範囲や形状のパラメータを入力するようにしてもよい。【０１０３】
この場合、画像処理装置403は、位相方向解析範囲重みの算出処理において、まず、表示装置405に図14に示すような解析範囲設定画面100を表示することが望ましい。【０１０４】
図14の解析範囲設定画面100には、1心拍に相当する心電情報101と、解析範囲を示す矩形オブジェクト102とが重ねて表示される。また、パラメータ入力のための中心位相入力欄104、解析範囲入力欄105等が設けられている。【０１０５】
矩形オブジェクト102は、その中心位置が重み形状の中心位置に該当し、また横幅が解析範囲に該当する。矩形オブジェクト102の中心位置や横幅は、マウス等の入力装置406からの入力指示に従って変化させることができるものとする。また、キーボード等の入力装置406(マウスポインタ103)を利用して、中心位相入力欄104や解析範囲入力欄105に数値が入力されると、入力された数値に応じて矩形オブジェクト102の中心位置や横幅を変更することも可能である。【０１０６】
画像処理装置403は、解析範囲設定画面100の矩形オブジェクト102によって設定された中心位相及び解析範囲に基づいて、例えば矩形の重み関数(図5(c)参照)を算出する。
【実施例８】
【０１０７】
また、例えば、図12のステップS203で実施する体軸方向解析範囲算出処理において、X線CT装置1の表示装置405及び入力装置406を利用して、操作者が直接、重みの範囲(解析範囲)や形状を指定するようにしてもよい。【０１０８】
この場合、画像処理装置403は、体軸方向解析範囲重みの算出処理において、まず、表示装置405に図15に示すような解析範囲設定画面110を表示することが望ましい。【０１０９】
図15の解析範囲設定画面110には、図12のステップS202で作成した評価用画像90(例えば、3次元画像)が表示され、その近傍に、体軸方向解析範囲重みの形状を設定するための線状オブジェクト112が並べて表示される。評価用画像90の体軸方向位置と、線状オブジェクト112の体軸方向位置は一致している。【０１１０】
線状オブジェクト112は、横軸を重みの値(「0.0」〜「1.0」)、縦軸を体軸方向として表示され、複数の線分が点オブジェクト113を介して接続されている。操作者が点オブジェクト113の位置をマウスポインタ114で移動することで、重み関数の形状を編集することができる。
【実施例９】
【０１１１】
図2、図12のフローチャートでは、予め動きの解析範囲を算出してから動き解析用の複数のデータセットを作成し、動き情報の算出を行った。しかし、処理手順はこれに限定されるものではなく、図16に示すフローチャートのように、まず、動き解析用の複数のデータセットを作成し、これらのデータセットの情報を操作者に提示(図17、図18参照)した後に、解析範囲重みを設定してもよい。【０１１２】
図16のフローチャートについて説明する。【０１１３】
まず、システム制御装置401は、撮影処理によって投影データ及び心電情報を取得し(ステップS301)、画像処理装置403へ送出する。そして、画像処理装置403は、動き解析用の画像を再構成し、複数のデータセットを作成し、表示する(ステップS302)。【０１１４】
図17は、複数のデータセットを表示した解析範囲設定画面120の一例である。【０１１５】
図17の解析範囲設定画面120には、解析用画像表示エリア121と、心電情報表示エリア122と、解析範囲重み設定用オブジェクト123とが設けられている。【０１１６】
解析用画像表示エリア121、心電情報表示エリア122、及び解析範囲重み設定用オブジェクト123の横軸は一致しており、いずれも撮影時刻(すなわち、体軸方向位置)を示す。【０１１７】
解析用画像表示エリア121には、ステップS302で作成した解析用画像の複数のデータセットが並べて表示される。図17の例では、縦方向に心拍位相5％刻みに、位相60％〜80％の解析用画像が並べられ、横方向には0.5秒刻み(10mm刻み)に0.0s〜3.5s(0.0mm〜70.0mm)の解析用画像が並べて表示される。【０１１８】
心電情報表示エリア122には、撮影中に得た心電情報が表示される。心電情報は、解析用画像表示エリア121の体軸方向に相当する縮尺で並べて表示される。【０１１９】
解析範囲重み設定用オブジェクト123は、複数の線分が点オブジェクト125を介して接続されたものである。操作者が点オブジェクト125の位置をマウスポインタ124で移動することで、重み関数の形状を編集することができる。解析範囲重み設定用オブジェクト123の横軸も、解析用画像表示エリア121の体軸方向に相当する縮尺で並べて表示される。【０１２０】
ステップS303の解析範囲算出処理において、操作者は、図17に示す解析範囲設定画面120の解析用画像や、心電情報を参照しながら、解析範囲重み設定用オブジェクト123の形状を編集する。【０１２１】
また、図17の解析範囲設定画面120に代えて、図18に示す解析範囲設定画面130を提示してもよい。【０１２２】
図18の解析範囲設定画面130には、図17と同様の解析用画像表示エリア131、及び心電情報表示エリア133が設けられている。解析用画像表示エリア131、及び心電情報表示エリア133の横軸は一致しており、いずれも撮影時刻(すなわち、体軸方向位置)を示す。【０１２３】
また、画面下部には、「選択範囲指定」モードまたは「除外範囲指定」モードのいずれかを選択するモード選択欄134、画像の選択状態を解除する際に操作される選択解除ボタン135、動きの解析を開始する際に操作される解析開始ボタン136等が設けられている。【０１２４】
図18の解析範囲設定画面130では、例えば「除外範囲指定」モードの場合、操作者が解析用画像や心電情報を目視しながら、マウスポインタ132により動き情報の解析から除外する心拍位相範囲と体軸方向範囲の断層像を選択する。「除外範囲指定」モードでは、選択した解析範囲の重みを「0.0」、それ以外を「1.0」として動き情報を算出する。逆に、「選択範囲指定」モードでは、選択した解析範囲の重みを「1.0」、それ以外を「0.0」として動き情報を算出するため、選択された断層像が動き情報の解析範囲に含まれる。【０１２５】
図18の例では、撮影時刻「2.0s〜2.5s」付近の心電情報に不整脈とみられる部分が含まれているため、その撮影時刻に相当する体軸方向範囲の解析用画像を「除外範囲指定」モードにて選択する(選択枠137)。すると撮影時刻「2.0s〜2.5s」で重み「0.0」、その他の撮影時刻で重み「1.0」の重み形状を設定することができる。【０１２６】
上述したように、図17または図18に示すような解析範囲設定画面120、130を利用して、体軸方向の解析範囲重みを設定すると、その後、画像処理装置403は、図2のステップS104〜ステップS107、図12のステップS205〜ステップS208と同様に、ステップS303で算出された解析範囲重みを動き情報に乗じた全体動き情報を算出し、所望の運動状態の位相を選択して診断用画像を再構成し、表示する(ステップS304〜ステップS307)。【０１２７】
図17や図18に示す解析範囲設定画面120、130を利用することにより、操作者は解析用画像や心電波形を目視しながら、解析用範囲重みを編集することができる。このため、より直観的に解析範囲重みを算出できる。【０１２８】
なお、図17や図18に示す解析範囲設定画面120、130では、複数の断層像を並べて選択対象としているが、3次元画像やMPR画像を並べて表示し、解析範囲の選択に利用してもよい。【０１２９】
以上説明したように、本発明に係るX線CT装置1は、動きのある診断部位について投影データを取得するとともに心電情報を取得し、取得した投影データ及び心電情報に基づいて動き解析用のデータセットを複数作成し、動き情報を算出する。そして、算出した動き情報に基づいて所望の運動状態における最適位相を決定し、最適位相における断層像を再構成する。この際、X線CT装置1は、所定範囲の動き情報を解析対象とするための解析範囲重みを算出し、動き情報に対して乗じることにより動きの解析範囲を調整する。【０１３０】
これにより、例えば、体動が生じた場合や、不整脈が生じた場合等、動きの解析に不都合な動き量が混入することを防ぎ、診断目的とする位相や部位における動き情報を正確に求めることが可能となる。これにより、断層像を再構成する最適位相を正確に決定することができる。【０１３１】
また、上述の解析範囲重みとして、任意の位相を基準とする位相方向に所定幅を持つ重み関数を設定することにより、例えば、撮影中に被検体の体動が生じた場合や、特定の心拍位相範囲にのみ影響が表れる構造物が撮影範囲に含まれる場合等に、余分の動き量が混入した位相を除いて動き情報を算出することができるため、動き情報の信頼性を向上できる。【０１３２】
また、心拍数によって位相方向の重みの範囲を動的に変更するようにすれば、撮影中の患者の心拍数に適した重みを設定できる。【０１３３】
また、静止心拍位相は臓器や撮影部位によって異なるため、部位によって位相方向の重みの範囲を変更するようにしてもよい。【０１３４】
また、解析範囲重みは、撮影時におけるX線強度に連動した形状とするようにしてもよい。この場合、X線強度を抑制して撮影された、診断上の有用性の低いとみなされる範囲を、解析対象から除いて動き情報を解析できる。【０１３５】
また、上述の解析範囲重みとして、任意の撮影部位を基準とする体軸方向の重み関数を設定すれば、例えば、心電計に不具合が生じた場合のように動き情報に悪影響を及ぼす撮影時刻(すなわち体軸方向位置)のデータを除いて動きを解析することができる。これにより、信頼性の高い動き情報を算出することができる。【０１３６】
特に、心電情報の変動が小さい撮影部位において重みの値が大きくなる形状の解析範囲重みを設定すれば、不整脈が生じた場合にその範囲を解析対象から除くことができる。【０１３７】
また、解析範囲重みは、診断の重要性が高い撮影部位において重みの値が大きくなる形状としてもよい。この場合は、診断の重要性の低い部位で生じた体動等の影響を受けにくくなるため、診断の重要性が高い撮影部位について動き情報を正確に算出することができる。【０１３８】
また、解析範囲重みは、静止位相算出の妨げになる構造物を含む撮影部位の重みの値が小さくなる形状としてもよい。この場合は、冠動脈内に留置されたステントや石灰化等の高吸収体(CT値が高い物体)が存在する場合に、動き量が実際より大きく算出されてしまい評価が正しく行えないという問題を解決できる。【０１３９】
また、解析範囲重みは、心拍位相方向及び体軸方向を含む多次元の重み情報としてもよい。この場合は、着目する身体部位と心拍位相の双方に最適化した動き情報を算出することが可能となる。【０１４０】
また、解析範囲重みを決定するためのパラメータを操作者が入力するようにしてもよい。この場合は、操作者の要望に応じて柔軟性の高い解析範囲の調整を行うことができる。【０１４１】
なお、上述の説明では、解析範囲重みは、心拍位相方向または体軸方向に乗じる例を説明したが、心拍位相方向解析範囲重み及び体軸方向解析範囲重みをともに算出し、組み合わせて使用することも可能である。【０１４２】
また、上述の各実施例で、投影データから解析用画像を再構成し、解析用画像を対象として動き情報を算出する例を説明したが、投影データを対象として動き情報を算出するようにしてもよい。【０１４３】
また、動き情報の解析及び心電同期再構成の説明において、心拍位相の表現を隣接するR波間の相対位置(％)を基準として説明したが、単独R波からの絶対時間を基準とする方法について適用してもよい。【０１４４】
また、対象を心臓として説明したが、本発明の適用対象は心臓に限定されるものではなく動きのある全ての部位の撮影に適用可能である。動きは周期的であっても非周期的であっても、本発明を適用できる。【０１４５】
その他、当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【０１４６】
1  X線CT装置、2  スキャナガントリ、3  寝台テーブル、201  X線源、202  測定制御装置、205  X線検出器、210  回転盤駆動装置、3  再構成演算器、4  操作卓、401  システム制御装置、403  画像処理装置、405  表示装置、406  入力装置、6  被検体、7  心電計
【図１】
2013187461.tif 000003
【図２】
2013187461.tif 000004
【図３】
2013187461.tif 000005
【図４】
2013187461.tif 000006
【図５】
2013187461.tif 000007
【図６】
2013187461.tif 000008
【図７】
2013187461.tif 000009
【図８】
2013187461.tif 000010
【図９】
2013187461.tif 000011
【図１０】
2013187461.tif 000012
【図１１】
2013187461.tif 000013
【図１２】
2013187461.tif 000014
【図１３】
2013187461.tif 000015
【図１４】
2013187461.tif 000016
【図１５】
2013187461.tif 000017
【図１６】
2013187461.tif 000018
【図１７】
2013187461.tif 000019
【図１８】
2013187461.tif 000020
【手続補正書】
【提出日】平成27年2月5日(2015.2.5)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
動きのある診断部位を対象とし、被検体周囲の複数の方向からX線を照射して投影データを取得し、前記投影データに基づき被検体の断層像を再構成し、表示するX線CT装置であって、
前記診断部位の動き情報を算出する動き情報算出部と、
所定範囲の動き情報を解析対象とするための解析範囲重みを算出する解析範囲算出部と、
前記動き情報算出部により算出された動き情報に対して前記解析範囲重みを乗じることにより得た動き情報に基づいて、再構成に適した位相である最適位相を決定する最適位相決定部と、
前記最適位相における断層像を再構成する画像再構成部と、
を具備することを特徴とするX線CT装置。
【請求項２】
前記解析範囲重みは、任意の位相を基準とする位相方向の重み関数であることを特徴とする請求項1に記載のX線CT装置。
【請求項３】
被検体の心電情報を取得する心電情報取得部と、
心拍数毎の標準的な静止位相を格納する心拍数基準標準静止位相データベースと、を更に備え、
前記解析範囲算出部は、前記心電情報取得部により取得した心拍数に対応する前記静止位相を前記心拍数基準標準静止位相データベースから取得し、取得した静止位相を基準とする位相方向の重み関数を前記解析範囲重みとして算出することを特徴とする請求項2に記載のX線CT装置。
【請求項４】
撮影対象臓器について、標準的な静止位相を診断対象部位毎に格納する部位基準標準静止位相データベースを更に備え、
前記解析範囲算出部は、診断対象とする部位に対応する前記静止位相を前記部位基準標準静止位相データベースから取得し、取得した静止位相を基準とする位相方向の重み関数を前記解析範囲重みとして算出することを特徴とする請求項2に記載のX線CT装置。
【請求項５】
前記解析範囲重みは、前記投影データの計測時におけるX線強度に応じた形状とすることを特徴とする請求項2に記載のX線CT装置。
【請求項６】
前記最適位相は、最も動き量の小さい位相であることを特徴とする請求項1に記載のX線CT装置。
【請求項７】
前記解析範囲重みは、任意の撮影部位を基準とする体軸方向の重み関数であることを特徴とする請求項1に記載のX線CT装置。
【請求項８】
被検体の心電情報を取得する心電情報取得部を更に備え、
前記解析範囲重みは、前記心電情報取得部により取得した心電情報の変動が小さい撮影部位において重みの値が大きくなる形状とすることを特徴とする請求項7に記載のX線CT装置。
【請求項９】
前記解析範囲重みは、診断の重要性が高い撮影部位において重みの値が大きくなる形状とすることを特徴とする請求項7に記載のX線CT装置。
【請求項１０】
前記解析範囲重みは、静止位相算出の妨げになる構造物を含む撮影部位の重みの値が小さくなる形状とすることを特徴とする請求項7に記載のX線CT装置。
【請求項１１】
前記解析範囲重みは、心拍位相方向及び体軸方向を含む多次元の重み情報であることを特徴とする請求項1に記載のX線CT装置。
【請求項１２】
前記解析範囲重みを決定するためのパラメータを入力する入力部を更に備え、
前記解析範囲算出部は、前記入力部により入力されたパラメータに基づいて前記解析範囲重みを決定することを特徴とする請求項1に記載のX線CT装置。
【請求項１３】
前記動きのある診断部位は心臓であって、
前記入力部は、心臓の1心拍に相当する心電情報と、前記心電情報に重畳され、前記解析範囲重みに対応するオブジェクトとを表示する表示部を備えることを特徴とする請求項12に記載のX線CT装置。
【請求項１４】
前記入力部は、前記動きのある診断部位の評価用画像と、前記評価用画像の体軸方向位置に対応する前記解析範囲重みを設定するために用いられるオブジェクトとを表示する表示部を備えることを特徴とする請求項12に記載のX線CT装置。
【請求項１５】
前記入力部は、前記動きのある診断部位の体軸方向位置を示す軸と位相を示す軸とを有するエリアであって、各軸に対応する複数の画像が並べられたエリアと、前記体軸方向位置に対応する前記解析範囲重みを設定するために用いられるオブジェクトとを表示する表示部を備えることを特徴とする請求項12に記載のX線CT装置。
【請求項１６】
前記入力部は、前記動きのある診断部位の体軸方向位置を示す軸と位相を示す軸とを有するエリアであって、各軸に対応する複数の画像が並べられたエリアを表示する表示部を備え、前記エリア中の画像を選択することにより前記解析範囲重みを設定することを特徴とする請求項12に記載のX線CT装置。
【請求項１７】
動きのある診断部位を対象とし、被検体周囲の複数の方向からX線を照射して投影データを取得し、前記投影データに基づき被検体の断層像を再構成する画像再構成方法であって、
前記診断部位の動き情報を算出する動き情報算出ステップと、
所定範囲の動き情報を解析対象とするための解析範囲重みを算出する解析範囲算出ステップと、
前記動き情報算出ステップにより算出した動き情報に対して前記解析範囲重みを乗じることにより得た動き情報に基づいて、再構成に適した位相である最適位相を決定する最適位相決定ステップと、
前記最適位相における断層像を再構成する画像再構成ステップと、
を含むことを特徴とする画像再構成方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
測定制御装置202は、操作卓4のシステム制御装置401からの制御信号に従って、スキャナガントリ2内のX線源201、コリメータ203、回転盤207の回転を制御する。
【国際調査報告】
2013187461.tif 000021
2013187461.tif 000022
2013187461.tif 000023
2013187461.tif 000024