Nothing Special   »   [go: up one dir, main page]

JPWO2013047439A1 - Ｘ線ｃｔ装置及び画像補正方法 - Google Patents

Ｘ線ｃｔ装置及び画像補正方法 Download PDF

Info

Publication number
JPWO2013047439A1
JPWO2013047439A1 JP2013536265A JP2013536265A JPWO2013047439A1 JP WO2013047439 A1 JPWO2013047439 A1 JP WO2013047439A1 JP 2013536265 A JP2013536265 A JP 2013536265A JP 2013536265 A JP2013536265 A JP 2013536265A JP WO2013047439 A1 JPWO2013047439 A1 JP WO2013047439A1
Authority
JP
Japan
Prior art keywords
time
image
scan
scans
ray
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2013536265A
Other languages
English (en)
Other versions
JP6104166B2 (ja
Inventor
角村　卓是
卓是 角村
國分　博人
博人 國分
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Publication of JPWO2013047439A1 publication Critical patent/JPWO2013047439A1/ja
Application granted granted Critical
Publication of JP6104166B2 publication Critical patent/JP6104166B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/90—Dynamic range modification of images or parts thereof
    • G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/481—Diagnostic techniques involving the use of contrast agents
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/486—Diagnostic techniques involving generating temporal series of image data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54—Control of apparatus or devices for radiation diagnosis
    • A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54—Control of apparatus or devices for radiation diagnosis
    • A61B6/545—Control of apparatus or devices for radiation diagnosis involving automatic set-up of acquisition parameters
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pulmonology (AREA)
  • Cardiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

Prospective Triggering法を用いた造影撮影において、スライス間の時間的不連続に起因する造影濃度ムラを補正することが可能なX線CT装置を提供する。
X線CT装置1では、心臓の各体軸方向の位置で心電情報のR波を検出してから所定の時間後にX線を照射してスキャンし、心臓の特定の位相の投影データを得る。画像処理装置4は、収集した投影データに基づいて断層像を再構成し、上述の各スキャンにおける投影データ収集時刻を心電計5から収集した心電情報から算出し、算出した投影データ収集時刻に基づいてスキャン間時間を求め、スキャン間時間に応じた範囲の体軸方向位置の画像について濃度差が滑らかに変化するように補正する。

Description

本発明は、X線CT装置に関し、特に心電計から収集した被検者の心電情報を用いて撮影及び再構成するX線CT装置等に関する。
X線CT装置は、X線源とそれに対向するように配置したX線検出器とを被検体を中心に回転させながらX線を照射し、被検体を透過したX線を全周に亘って検出し、得られた投影データに基づいて画像再構成して断層像を得るものである。この画像再構成では、被検体または被検体内の臓器は静止していることが必要条件となっており、投影データを収集している間に被検体や臓器が動いているとその対象は画像再構成の際に正しく結像されず、ぼけた画像(モーションアーチファクト)となる。
特に心臓は常に拍動している臓器であるため、心臓の静止した画像を得るためには投影データを収集する時間、すなわち時間分解能を向上させる必要がある。そこでX線CT装置を用いた心臓撮影では、被検者に心電計を装着し、被検者の心電情報を用いて心臓の動きが少ない位相の投影データを収集して画像再構成することで比較的静止した心臓の画像を得る心電同期再構成法が行われている。
X線CT装置を用いた心臓撮影における撮影方法には、大きく分けてRetrospective Gating法とProspective Triggering法の2つがある。
Retrospective Gating法は、先に述べたように被検者の心電情報を収集しながらヘリカルスキャンし、得られた投影データの中から心電情報を用いて動きが少ない位相の投影データのみを収集して画像再構成する方法である。この方法はX線を連続的に照射してヘリカルスキャンするため、任意の位相の画像を作成することができる。しかし一方で複数の心拍から同位相の投影データを収集しなければならないため、テーブルの移動速度は遅くなり、被ばく量は通常のヘリカルスキャンの数倍となる。
Prospective Triggering法は、被検者の心電情報をモニタリングし、例えばR波を検出してから予め設定した特定の位相になる時間後にのみX線を照射してアキシャルスキャンする方法である。この場合、X線を連続的に照射せずに特定の位相になったときのみ照射するので、Retrospective Gating法と比べて非常に少ない被ばく量で済む(特許文献1)。
両者それぞれにメリットとデメリットがあるが、被検者への負担を考慮すると、被ばくの少ないProspective Triggering法で撮影するのが望ましい。
特開平9-24045号公報
しかしながら、上述のProspective Triggering法は、同じ体軸方向位置で1回しかスキャンしないため、次の体軸方向位置のスキャンは前のスキャンとは異なる心拍から投影データを収集することになる。そのため、スライス間で時間的に連続しない部分ができる。
したがってProspective Triggering法で冠動脈造影撮影を行うと、時間的に連続しない部分で造影剤の濃度が大きく変わり、MPR(Multi Plannar Reconstruction)画像を作成すると違和感のある画像となる。また、画像処理により冠動脈を抽出しようとしても、この造影濃度の急激な変化により正しく抽出できない場合がある。
本発明は、以上の問題点に鑑みてなされたものであり、Prospective Triggering法を用いた造影撮影において、スライス間の時間的不連続に起因する造影濃度ムラを補正することが可能なX線CT装置等を提供することを目的とする。
前述した目的を達成するために第1の発明は、被検者の心電情報を心電計から収集するとともに、収集した心電情報に同期したスキャンを実行し、スキャンで得た投影データを再構成することにより被検者の断層像を得るX線CT装置において、所定のタイミングでX線を照射してスキャンすることにより、心臓の各体軸方向位置で心臓の特定の位相の投影データを複数の時相で断続的に取得するデータ収集部と、前記データ収集部により取得した投影データに基づいて各体軸方向位置での断層像を再構成する再構成部と、あるスキャンから次のスキャンまでのスキャン間時間を算出し、算出した前記スキャン間時間の長さに応じて決定される体軸方向範囲の断層像について各断層像間で生じた濃度差を補正する補正部と、前記補正部により補正された断層像を基に生成される画像を表示する表示部と、を備えることを特徴とするX線CT装置である。
また、第2の発明は、被検者の心電情報を心電計から収集するとともに、収集した心電情報に同期したスキャンを実行し、スキャンで得た投影データを再構成することにより被検者の断層像を得るX線CT装置における画像補正方法であって、所定のタイミングでX線を照射してスキャンすることにより、心臓の各体軸方向位置で心臓の特定の位相の投影データを複数の時相で断続的に取得するデータ収集ステップと、取得した投影データに基づいて各体軸方向位置での断層像を再構成する再構成ステップと、あるスキャンから次のスキャンまでのスキャン間時間を算出し、算出した前記スキャン間時間の長さに応じて決定される体軸方向範囲の断層像について各断層像間で生じた濃度差を補正する補正ステップと、補正された断層像を基に生成される画像を表示する表示ステップと、を含むことを特徴とする画像補正方法である。
本発明のX線CT装置及び画像補正方法により、Prospective Triggering法を用いた造影撮影において、スライス間の時間的不連続に起因する濃度ムラを補正することが可能となり、違和感のないMPR画像または対象部位の抽出処理等の画像処理を容易に行うことができる。
X線CT装置1の構成図 画像処理装置4の構成図 X線CT装置1の動作手順を示すフローチャート Prospective Triggering Scanにおける心電波形とディレータイム、X線照射タイミングの関係を説明する概念図 X線の形状とテーブル移動量との関係を説明する概念図 スキャン間境界位置と濃度ムラについて説明する概念図 スキャン間時間を求める方法を説明する概念図 スキャン間時間が短い場合の補正範囲の一例 スキャン間時間が長い場合の補正範囲の一例 補正範囲の決定方法について説明する図 補正に用いる濃度線の別の例 スキャン間時間が短い場合のオーバーラップデータ寄与率を示す図 スキャン間時間が長い場合のオーバーラップデータ寄与率を示す図 補正範囲や補正の程度を示す表示例 補正範囲と補正の程度を示す別の表示例
以下、添付図面を参照しながら、本発明の好適な実施形態について詳細に説明する。
[実施形態]
まず、図1を参照して、X線CT装置1の構成について説明する。
X線CT装置1は、ガントリ100と、X線の照射を制御するX線制御器201と、ガントリ100の動作を制御するガントリ制御器202と、テーブル109の動作を制御するテーブル制御器203と、X線検出器105に入射したX線の強度を電気信号に変換し、収集するデータ収集装置(DAS；Data Acquisition System)204と、DAS204で変換した電気信号を取得し、取得したデータの補正及び画像再構成を行う再構成演算器3と、被検者6の心電情報を収集する心電計5と、心電計5で計測した心電情報を取得するとともに、X線制御器201、ガントリ制御器202、テーブル制御器203、及びDAS204を制御してガントリ100によるスキャンを実行し、また再構成演算器3にて再構成された断層像を取得する画像処理装置4と、を備える。
ガントリ100には、X線源101とX線検出器105とが被検者6の載置されるテーブル109を介して対向配置されている。
X線源101はX線制御器201に制御されて所定の強度のX線を照射する。X線検出器105はX線源101から照射され、被検者6を透過したX線を検出し、検出した透過X線の強度に従った電気信号を発生する。スキャンの際、X線源101及びX線検出器105はガントリ制御器202から入力される制御信号に従って回転中心の周りを回転する。
X線源101及びガントリ100の動作は、X線制御器201とガントリ制御器202によって制御される。X線制御器201は、X線源101に電力信号およびX線発生タイミング信号を供給する。ガントリ制御器202は、ガントリ100の回転速度及び位置を制御する。テーブル109はテーブル制御器203によって制御される。テーブル制御器203は、テーブル109の移動速度および位置を制御する。
X線検出器105に入射した透過X線は、DAS204によってデジタル信号に変換され、デジタルデータとして画像処理装置4に送られる。
再構成演算器3は画像処理装置4からデジタルデータを取得すると、感度補正、対数変換、オフセット補正等のデータ補正処理を行って投影データを生成し、投影データを用いて画像再構成処理を実行する。画像再構成処理によって再構成された画像データ(断層像)は画像処理装置4に入力され、データ記録装置403(図2)に保存されるとともに、画像表示装置405に表示される。
図2は画像処理装置4の構成図である。
画像処理装置4は主として、上述の各構成要素の動作を制御する中央処理装置(CPU401)と、画像処理装置4の制御プログラムが格納された主メモリ402と、画像データを格納するデータ記録装置403と、被検者6の画像データを一時記憶する表示メモリ404と、この表示メモリ404に一時記憶された画像データに基づく表示を行う画像表示装置405と、画像表示装置405上のソフトスイッチを操作するためのマウス、タッチパネル等のポインティングデバイス407及びコントローラ408と、各種パラメータ設定用のキーやスイッチを備えたキーボード等の外部入力装置406と、画像処理装置4をローカルエリアネットワーク、電話回線、インターネット等のネットワークに接続するためのネットワークアダプタ409と、上記各構成要素を接続するデータバス410とから構成される。データ記録装置403は、磁気ディスク等の記憶装置や、取り出し可能な外部メディアに対してデータの書き込みや読み出しを行う装置でもよい。画像処理装置4はネットワークアダプタ409及びネットワーク411を介して外部の画像データベース412と接続し、画像データベース412との間で画像データを送受信するようにしてもよい。
心電計5は、被検者6に取り付けた電極を介して、心臓の心拍運動を反映した活動電位の時間変化を表す心電情報を計測し、例えば0．1秒間隔等の所定のサンプリングピッチでデジタル信号に変換する。心電計5で得た心電情報は、画像処理装置4へ順次送出される。
次に、図3を参照して、X線CT装置1の動作について説明する。
本実施の形態のX線CT装置1は、図3のフローチャートに示す手順で心電同期撮影を実行する。すなわち、画像処理装置4のCPU401は、主メモリ402から図3に示す心電同期撮影処理に関するプログラム及びデータを読み出し、このプログラム及びデータに基づいて処理を実行する。
本処理では、冠動脈の造影撮影をProspective Triggering Scan法にて実行する。画像処理装置4は、心臓の静止位相の投影データを収集し、収集した投影データを用いて画像再構成する(ステップS1)。
Prospective Triggering Scanにおいて、X線CT装置1は、撮影前に予め被検者6の心拍情報に対してどの位相をスキャンするかの設定を任意の基準位相、例えばR波からのディレータイムの設定として受け付ける。スキャンの際、画像処理装置4は、被検者6の心電情報をモニタリングし、図4に示すようにR波を検出してからディレータイム後にX線を照射するよう制御する。
X線CT装置1は、R波を検出してから上述のディレータイム後に、体軸方向位置Z1に対してX線を照射し、スキャンする。その後、テーブル109を次の体軸方向位置Z2まで移動させ、R波を検出したら再度ディレータイム後にX線を照射し、スキャンする。
このとき例えば64スライス以上等の多列のX線検出器105を有するX線CT装置1では、X線の形状はコーン状となるが、再構成画像の体軸方向の連続性を保つため、図5に示すようにデータ収集時刻の異なる各スキャン間のテーブル移動量をガントリ100の回転中心におけるX線照射幅より短くして、各スキャン間でX線照射領域にオーバーラップ区間を設けることが望ましい。また、画像を再構成するスライス数(画像作成領域)もX線検出器105のスライス数未満となるように設定しておくことが望ましい。
再構成演算器3は、以上のようにして得られた投影データを基に画像再構成する。以上の処理により、心電情報の特定位相(R波検出後、所定のディレータイム経過後)の断層像を、複数の時相(データ収集時刻)にわたって断続的に得る。
このようにして得た断層像を用いてMPR画像を作成すると、図6に示す画像70のように、体軸方向で不自然な濃度差が生じることがある。これは、異なる心拍、すなわち異なるデータ収集時刻における投影データを用いて1枚のMPR画像70を作成するためである。各心拍で時間差があるため造影剤の濃度にも変化が生じ、異なる時相(データ収集時刻)の各スキャン間の境界スライス71、72、73、74で、濃度差が大きくなり、濃度ムラとなる。
こうした濃度ムラの発生を防ぐため、本発明では、ステップS2〜ステップS4の手順で補正処理を行う。
まず、画像処理装置4は、各スキャンのデータ収集時刻及びスキャン間時間を算出する(ステップS2)。
データ収集時刻とは、各スキャンにおける投影データの収集時刻である。また、画像処理装置4はあるスキャン(k回目)と次のスキャン(k＋1回目)との間の時間(スキャン間時間)を求める。
Prospective Triggering Scanでは、R波の検出後、所定のディレータイム経過をトリガとして体軸方向の各位置でスキャン(データ収集)するため、体軸方向で時間的に不連続なスライスが生じる。
スキャン間時間を求めるために、画像処理装置4のCPUは、まず、最初のスキャンのデータ収集時刻(スキャン開始時刻)から、各スキャンのデータ収集時刻までの時間を求める。
図7に示すように、投影データの各ビュー(View)と心電情報を時間的に対応させて投影データを収集し、例えば、最初にX線が照射されたビュー範囲(図7の斜線部)の中心から各スキャンでX線が照射されたビュー範囲の中心までのビュー数を求める。次に、求めたビュー数を時間に換算することによって、最初のスキャンから各スキャンまでに要する時間が求められる。最初のスキャンから各スキャンまでに要する時間(IT[ms])は、スキャナ1回転あたりの回転時間(ST[ms／rot])とその回転時間におけるスキャナ1回転あたりのビュー取り込み数(ビューレート：VR[view／rot])と先で求めたスキャン間のビュー数(N[view])と、から、以下の式(1)を用いて求めることができる。
Figure 2013047439
 なお、図7の例では、最初のスキャンから各スキャンまでの時間を求めるための基準位置をX線が照射されたビュー範囲の中心から中心までとしているが、これに限定されず、基準位置はビュー範囲の開始位置や終了位置であっても構わない。
スキャン間時間は、以上のようにして求めた各スキャンの最初のスキャンからの時間ITを互いに引くことで求められる。なお、スキャン間時間は、時間ではなく、R波の個数として表すようにしてもよい。
次に、画像処理装置4は、画像補正を行う断層像の体軸方向範囲(以下、補正範囲という)を算出する(ステップS3)。
画像処理装置4は、ステップS2で求めたスキャン間時間に応じて適切な補正範囲を決定する。
以下、図8及び図9を参照して補正範囲の求め方の一例を示す。
ステップS2で算出したスキャン間時間が短い場合、図8(A)に示すように、造影剤の濃度変化が比較的小さくなることが予想される。この場合、図8(B)に示すように、スキャン間境界位置80を中心として補正範囲を狭くとることで、造影剤の量が体軸方向に沿って滑らかに変化しているように見えるように画像補正する。
一方、ステップS2で算出したスキャン間時間が長い場合は、図9(A)に示すように造影剤の濃度変化が大きくなることが予想される。この場合、図9(B)に示すように、スキャン間境界位置90を中心として補正範囲を広くとることで、造影剤の量が体軸方向に沿って滑らかに変化しているように見えるように画像補正する。
なお、図8、図9に示すスキャン間境界位置80、90は、図6に示す境界スライス71、72、73、74のいずれかに対応するものである。
補正範囲を求める手順の詳細を、図10を参照して説明する。
まず、画像処理装置4は、ステップS1で得た元の断層像群のそれぞれ対応する画素の画素値から、図10(A)に示すように、体軸方向位置での造影濃度の違いを示す濃度曲線C1を得る。なお、図10(A)では説明を分かり易くするために濃度曲線を階段状としているが、必ずしも階段状となるものではない。
次に、画像処理装置4は、図10(B)に示すように、スキャン間境界位置80の濃度差81を求めて、その中央値(濃度中点82)を算出する。
次に、画像処理装置4は、図10(C)に示すように、濃度中点82を通る、ある傾きの濃度線83を算出する。
ここで、画像処理装置4は、スキャン間時間が長い場合は濃度線83の傾きを小さくし、スキャン間時間が短い場合は濃度線83の傾きを大きくするようにする。
そして、画像処理装置4は、図10(D)に示すように、濃度線83と濃度曲線C1との2つの交点84、85を算出し、図10(Ｅ)に示すように、2つの交点84、85間を、補正範囲86と決定する。
これにより、スキャン間時間が長く造影濃度差81が大きい場合は広い補正範囲となり、スキャン間時間が短く、造影濃度差81が小さい場合は狭い補正範囲となる。
なお、濃度線83は、図8〜10に示すような直線に限定されず、図11に示すような曲線としてもよい。その他、濃度線83の形状はこの限りではない。
また、ステップS1のスキャン時に各スキャンでX線照射幅をオーバーラップさせて投影データを得ている場合は、X線がオーバーラップする区間では、同じ体軸位置でデータ収集時刻の異なる2スキャン分の投影データを取得できる。
そのため、上述のステップS3で決定した補正範囲の断層像を補正する際に、上述のスキャン間時間の長さに応じて決定される寄与率とデータ収集時刻の異なる各投影データとを用いて画像を生成し、これにより各断層像間の濃度値ムラを補正するようにしてもよい。
そして、スキャン間時間が短い場合は上述の寄与率を小さくし、スキャン間時間が長い場合は上述の寄与率を大きくすることが望ましい。
すなわち、オーバーラップ区間では、データ収集時刻ｔ1の投影データとデータ収集時刻ｔ2の投影データとが収集されている。スキャン間境界位置80ではデータ収集時刻ｔ1の投影データとデータ収集時刻ｔ2の投影データを互いに50パーセントずつ用いて画像を生成する。また、スキャン間境界位置80から離れるに従って、時間的に遠い投影データの寄与率(投影データを利用する割合)を徐々に小さくなるようにして、画像を生成する。
このとき、上述のように、スキャン間時間が短い場合、造影剤の濃度変化は比較的小さい。したがって補正をあまりしなくても滑らかな濃度変化が得られる。これよりスキャン間の時間が短い場合は、図12に示すようにスキャン間境界位置80から離れた位置では時間的に遠い投影データの寄与率を少なくする。一方、スキャン間の時間が長い場合は、境界スライスから離れた位置でも補正を行わないと濃度差を抑えることができない。したがって図13に示すように境界スライスから離れた位置においても時間的に遠い投影データの寄与率を多くする。
補正範囲が決定されると、画像処理装置4はステップS3で求めた補正範囲のスキャン間境界部分の画像を補正する(ステップS4)。
なお、補正処理は、各断層像の特に造影部分についてのみ実行するようにしてもよい。
これにより着目される部分について高速に補正処理を行うことが可能となる。
次に、画像処理装置4は、ステップS3で求めた補正範囲を、補正後の画像95と併せて画像表示装置405に表示する(ステップS5)。
例えば、図14に示すように、ステップS3で求めた補正範囲の開始位置及び終了位置を補正後の画像95上にラインL1、L2で表示する。
また、例えば、図14に示すように、補正後の画像95上または当該画像95の近くに、画像位置と対応するように補正に用いた重み係数(オーバーラップ区間内の寄与率)を示すグラフ96で表示する。
或いは、図15に示すように、補正後の画像95上または画像95の近くに、画像位置と対応するように補正に用いた重み係数(濃度線の傾斜角度またはオーバーラップ区間内の寄与率)に対応する濃度や色を示す、補正バー97を表示する。
なお、これらの表示方法に限定されず、補正範囲、濃度線の傾斜角度、オーバーラップ区間内の寄与率等、補正範囲と補正の程度を示す指標をマーク、グラフ、図表、カラー表示等を用いて表示すればよい。また、補正範囲を表示することなく、補正後の画像95のみを画像表示装置405に表示しても良いし、補正後の画像95と補正前の画像とを並べて画像表示装置405に表示しても良い。
以上説明したように、本発明のX線CT装置1では、心電情報のR波を検出してから所定の時間後にX線を照射してスキャンすることにより、心臓の各体軸方向位置で心臓の特定の位相の投影データを複数の時相で断続的に取得し、取得した投影データに基づいて各体軸方向位置での断層像を再構成する(Prospective Triggering Scan)。そして、心電情報に基づいてあるスキャンから次のスキャンまでのスキャン間時間を算出し、算出したスキャン間時間の長さに応じて決定される体軸方向範囲の断層像について各断層像間で生じた濃度差を補正する。また、補正された断層像を基に生成される画像を表示するとともに、補正された体軸方向範囲を表示する。
したがって、Prospective Triggering法を用いた造影撮影において、スライス間の時間的不連続に起因する濃度ムラを補正することが可能となり、違和感のないMPR画像または対象部位の抽出処理等の画像処理を容易に行うことができるようになる。
また、スキャン間時間が短い場合は補正範囲を狭くし、スキャン間時間が長い場合は補正範囲を広くすることで、濃度変化を一様にすることができる。
また、あるデータ収集時刻におけるスキャンのX線照射領域と次のデータ収集時刻におけるスキャンのX線照射領域とにオーバーラップ区間を設けてスキャンする場合において、スキャン間時間の長さに応じて決定される寄与率でデータ収集時刻の異なる各投影データを用いて画像を生成するようにしてもよい。
この場合、スキャン間時間が短い場合は時間的に離れた投影データの寄与率を小さくし、スキャン間時間が長い場合は時間的に離れた投影データの寄与率を大きくすれば、各スキャンでの投影データを使用して信頼性の高い画像を作成しつつ、濃度ムラのない画像を得ることが可能となる。
また、造影された部分のみを補正対象とすれば、特に着目する部位について高速に補正処理を行える。
なお、本発明は、上述の実施形態に限定されるものではない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
1 X線CT装置、100 ガントリ、101 X線源、105 X線検出器、109 テーブル、201 X線制御器、202 ガントリ制御器、203 テーブル制御器、204 DAS、3 再構成演算器、4 画像処理装置、5 心電計、6 被検者、C1 濃度曲線、80、90 スキャン間境界位置、81 濃度差、82 濃度中点、83 濃度線、84、85 交点、95 補正後の画像、96、97 補正の指標
X線CT装置は、X線源とそれに対向するように配置したX線検出器とを被検者を中心に回転させながらX線を照射し、被検者を透過したX線を全周に亘って検出し、得られた投影データに基づいて画像再構成して断層像を得るものである。この画像再構成では、被検者または被検者内の臓器は静止していることが必要条件となっており、投影データを収集している間に被検者や臓器が動いているとその対象は画像再構成の際に正しく結像されず、ぼけた画像(モーションアーチファクト)となる。
しかしながら、上述のProspective Triggering法は、同じ体軸方向位置で1回しかスキャンしないため、次の体軸方向位置のスキャンは前のスキャンとは異なる心拍から投影データを収集することになる。そのため、スライス間で時間的に連続しない部分ができる。
したがってProspective Triggering法で冠動脈造影撮影を行うと、時間的に連続しない部分で造影剤の濃度が大きく変わり、MPR(Multi Planar Reconstruction)画像を作成すると違和感のある画像となる。また、画像処理により冠動脈を抽出しようとしても、この造影濃度の急激な変化により正しく抽出できない場合がある。

Claims (12)

  1. 被検者の心電情報を心電計から収集するとともに、収集した心電情報に同期したスキャンを実行し、スキャンで得た投影データを再構成することにより被検者の断層像を得るX線CT装置において、
    所定のタイミングでX線を照射してスキャンすることにより、心臓の各体軸方向位置で心臓の特定の位相の投影データを複数の時相で断続的に取得するデータ収集部と、
    前記データ収集部により取得した投影データに基づいて各体軸方向位置での断層像を再構成する再構成部と、
    あるスキャンから次のスキャンまでのスキャン間時間を算出し、算出した前記スキャン間時間の長さに応じて決定される体軸方向範囲の断層像について各断層像間で生じた濃度差を補正する補正部と、
    前記補正部により補正された断層像を基に生成される画像を表示する表示部と、
    を備えることを特徴とするX線CT装置。
  2. 前記表示部は、前記補正部により補正された体軸方向範囲を生成された画像とともに表示することを特徴とする請求項1に記載のX線CT装置。
  3. 前記補正部は、補正対象とする断層像の体軸方向範囲を、前記スキャン間時間が短い場合は狭くし、前記スキャン間時間が長い場合は広くすることを特徴とする請求項1に記載のX線CT装置。
  4. 前記データ収集部が、あるデータ収集時刻におけるスキャンのX線照射領域と次のデータ収集時刻におけるスキャンのX線照射領域とにオーバーラップ区間を設けてスキャンする場合において、
    前記補正部は、前記スキャン間時間の長さに応じて決定される寄与率とデータ収集時刻の異なる各投影データとを用いて画像を生成することにより各断層像間で生じた濃度差を補正することを特徴とする請求項1に記載のX線CT装置。
  5. 前記補正部は、前記スキャン間時間が短い場合は前記寄与率を小さくし、前記スキャン間時間が長い場合は前記寄与率を大きくすることを特徴とする請求項4に記載のX線CT装置。
  6. 前記補正部は、前記断層像の造影された部分のみを補正対象とすることを特徴とする請求項1に記載のX線CT装置。
  7. 被検者の心電情報を心電計から収集するとともに、収集した心電情報に同期したスキャンを実行し、スキャンで得た投影データを再構成することにより被検者の断層像を得るX線CT装置における画像補正方法であって、
    所定のタイミングでX線を照射してスキャンすることにより、心臓の各体軸方向位置で心臓の特定の位相の投影データを複数の時相で断続的に取得するデータ収集ステップと、 取得した投影データに基づいて各体軸方向位置での断層像を再構成する再構成ステップと、
    あるスキャンから次のスキャンまでのスキャン間時間を算出し、算出した前記スキャン間時間の長さに応じて決定される体軸方向範囲の断層像について各断層像間で生じた濃度差を補正する補正ステップと、
    補正された断層像を基に生成される画像を表示する表示ステップと、
    を含むことを特徴とする画像補正方法。
  8. 前記表示ステップは、前記補正ステップにより補正された体軸方向範囲を生成された画像とともに表示することを特徴とする請求項7に記載の画像補正方法。
  9. 前記補正ステップは、補正対象とする断層像の体軸方向範囲を、前記スキャン間時間が短い場合は狭くし、前記スキャン間時間が長い場合は広くすることを特徴とする請求項7に記載の画像補正方法。
  10. 前記データ収集ステップが、あるデータ収集時刻におけるスキャンのX線照射領域と次のデータ収集時刻におけるスキャンのX線照射領域とにオーバーラップ区間を設けてスキャンする場合において、
    前記補正ステップは、前記スキャン間時間の長さに応じて決定される寄与率とデータ収集時刻の異なる各投影データとを用いて画像を生成することにより各断層像間で生じた濃度差を補正することを特徴とする請求項7に記載の画像補正方法。
  11. 前記補正ステップは、前記スキャン間時間が短い場合は前記寄与率を小さくし、前記スキャン間時間が長い場合は前記寄与率を大きくすることを特徴とする請求項10に記載の画像補正方法。
  12. 前記補正ステップは、前記断層像の造影された部分のみを補正対象とすることを特徴とする請求項7に記載の画像補正方法。
JP2013536265A 2011-09-27 2012-09-24 Ｘ線ｃｔ装置及び画像補正方法 Active JP6104166B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2011210120 2011-09-27
JP2011210120 2011-09-27
PCT/JP2012/074397 WO2013047439A1 (ja) 2011-09-27 2012-09-24 Ｘ線ｃｔ装置及び画像補正方法

Publications (2)

Publication Number Publication Date
JPWO2013047439A1 true JPWO2013047439A1 (ja) 2015-03-26
JP6104166B2 JP6104166B2 (ja) 2017-03-29

Family

ID=47995471

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013536265A Active JP6104166B2 (ja) 2011-09-27 2012-09-24 Ｘ線ｃｔ装置及び画像補正方法

Country Status (4)

Country Link
US (1) US9129389B2 (ja)
JP (1) JP6104166B2 (ja)
CN (1) CN103764034B (ja)
WO (1) WO2013047439A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101664432B1 (ko) * 2014-02-12 2016-10-10 삼성전자주식회사 단층 영상 장치 및 그에 따른 단층 영상 디스플레이 방법
US9619889B2 (en) * 2014-10-21 2017-04-11 General Electric Company Methods and systems for normalizing contrast across multiple acquisitions
JP7051307B2 (ja) * 2016-05-09 2022-04-11 キヤノンメディカルシステムズ株式会社 医用画像診断装置
DE202019003376U1 (de) 2019-03-21 2019-09-13 Ziehm Imaging Gmbh Röntgensystem zur iterativen Bestimmung einer optimalen Koordinatentransformation zwischen überlappenden Volumina, die aus Volumendatensätzen von diskret abgetasteten Objektbereichen rekonstruiert wurden
CN117115577B (zh) * 2023-10-23 2023-12-26 南京安科医疗科技有限公司 一种心脏ct投影域最优相位识别方法、设备及介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001190547A (ja) * 1999-11-28 2001-07-17 Siemens Ag 周期的運動を行う身体範囲の検査方法
JP2003210456A (ja) * 2002-01-21 2003-07-29 Toshiba Corp 時系列画像の処理装置
JP2004057506A (ja) * 2002-07-29 2004-02-26 Toshiba Corp ディジタル画像処理装置及びｘ線診断装置
JP2004305358A (ja) * 2003-04-04 2004-11-04 Hitachi Medical Corp Ｘ線ｃｔ装置
JP2006021022A (ja) * 2004-06-11 2006-01-26 Toshiba Corp Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
WO2006077869A1 (ja) * 2005-01-18 2006-07-27 Hitachi Medical Corporation Ｘ線ｃｔ装置
JP2006326078A (ja) * 2005-05-27 2006-12-07 Hitachi Medical Corp 血流動態解析装置、ｘ線ｃｔ装置、ｍｒｉ装置、及び血流動態解析プログラム
US20110075906A1 (en) * 2009-09-29 2011-03-31 Thomas Allmendinger Scanning of a cyclically moved examination object using a contrast agent within the scope of a preliminary examination by means of a ct scanner

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU5653599A (en) * 1998-09-18 2000-04-10 Mitsubishi Plastics Inc. Image processing device and method, and recording medium
DE60212917T2 (de) * 2001-10-16 2007-03-01 Kabushiki Kaisha Toshiba Vorrichtung zur Berechnung eines Index von örtlichen Blutflüssen
AR047692A1 (es) * 2003-07-10 2006-02-08 Epix Medical Inc Imagenes de blancos estacionarios
CN100462050C (zh) * 2004-06-11 2009-02-18 株式会社东芝 X射线ct装置和心肌灌注图像产生系统
US7352840B1 (en) * 2004-06-21 2008-04-01 Radiation Monitoring Devices, Inc. Micro CT scanners incorporating internal gain charge-coupled devices
DE102005052368B4 (de) * 2005-10-31 2015-07-30 Bayer Pharma Aktiengesellschaft Röntgensystem zur Erstellung diagnostischer Röntgendarstellungen unter Applikation von Kontrastmitteln
WO2007138979A1 (ja) * 2006-05-25 2007-12-06 Hitachi Medical Corporation Ｘ線ｃｔ装置
JP5019930B2 (ja) * 2007-04-05 2012-09-05 富士フイルム株式会社 放射線断層画像取得装置
JP4966120B2 (ja) * 2007-07-18 2012-07-04 株式会社東芝 Ｘ線アンギオ撮影装置
JP5631698B2 (ja) * 2009-12-07 2014-11-26 株式会社東芝 医用画像処理装置及び医用画像処理方法
US9119590B2 (en) * 2010-02-18 2015-09-01 Los Angeles Biomedical Research Institute At Harbor-Ucla Medical Center Method for thoracic vertebral bone density measurement by thoracic quantitative computed tomography
US8718747B2 (en) * 2010-04-16 2014-05-06 Oslo Universitetssykehus Hf Estimating and correcting for contrast agent extravasation in tissue perfusion imaging
US8768031B2 (en) * 2010-10-01 2014-07-01 Mistretta Medical, Llc Time resolved digital subtraction angiography perfusion measurement method, apparatus and system

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001190547A (ja) * 1999-11-28 2001-07-17 Siemens Ag 周期的運動を行う身体範囲の検査方法
JP2003210456A (ja) * 2002-01-21 2003-07-29 Toshiba Corp 時系列画像の処理装置
JP2004057506A (ja) * 2002-07-29 2004-02-26 Toshiba Corp ディジタル画像処理装置及びｘ線診断装置
JP2004305358A (ja) * 2003-04-04 2004-11-04 Hitachi Medical Corp Ｘ線ｃｔ装置
JP2006021022A (ja) * 2004-06-11 2006-01-26 Toshiba Corp Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
WO2006077869A1 (ja) * 2005-01-18 2006-07-27 Hitachi Medical Corporation Ｘ線ｃｔ装置
JP2006326078A (ja) * 2005-05-27 2006-12-07 Hitachi Medical Corp 血流動態解析装置、ｘ線ｃｔ装置、ｍｒｉ装置、及び血流動態解析プログラム
US20110075906A1 (en) * 2009-09-29 2011-03-31 Thomas Allmendinger Scanning of a cyclically moved examination object using a contrast agent within the scope of a preliminary examination by means of a ct scanner

Also Published As

Publication number Publication date
US20140212016A1 (en) 2014-07-31
CN103764034B (zh) 2016-03-02
US9129389B2 (en) 2015-09-08
JP6104166B2 (ja) 2017-03-29
CN103764034A (zh) 2014-04-30
WO2013047439A1 (ja) 2013-04-04

Similar Documents

Publication Publication Date Title
US7684537B2 (en) X-ray CT apparatus
JP6618900B2 (ja) Ｘ線ｃｔ装置及び画像再構成方法
JP5643218B2 (ja) Ｘ線ｃｔ装置及びｘ線ｃｔ装置による画像表示方法
US9542762B2 (en) X-ray CT apparatus and image reconstruction method
US6434215B1 (en) EKG-less cardiac image reconstruction
JP2004275440A (ja) Ｘ線コンピュータ断層撮影装置
JP2009028065A (ja) Ｘ線ｃｔ装置
US7831011B2 (en) Computed tomography method and system
JP6509131B2 (ja) Ｘ線ｃｔ装置、画像処理装置、及び画像再構成方法
JP6104166B2 (ja) Ｘ線ｃｔ装置及び画像補正方法
JP2004121840A (ja) 周期的に運動する器官のｃｔ画像形成方法およびこの方法を実施するためのｃｔ装置
US11160523B2 (en) Systems and methods for cardiac imaging
JP2003204961A (ja) Ｘ線ｃｔ装置
CN102028493B (zh) 在使用造影剂的条件下对周期性运动检查对象的改善扫描
CN104755030A (zh) 借助动态准直进行剂量减小的ct拍摄
US9737279B2 (en) X-ray CT apparatus
WO2004071301A1 (ja) X線ct装置
JP4007928B2 (ja) Ｘ線ｃｔ装置
JP7024954B2 (ja) 医用画像処理装置及びｘ線ｃｔ装置
US11559272B2 (en) Methods and systems for computed tomography
JP2001161680A (ja) 同期心臓走査の間にｘ線照射量を制御するための方法及び装置
WO2014188936A1 (ja) Ｘ線ｃｔ装置及び撮影方法
JP5931642B2 (ja) Ｘ線ｃｔ装置
JP4338962B2 (ja) Ｘ線ｃｔ装置
US7023958B2 (en) Radiation image-acquiring apparatus, and radiation image-acquiring method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150827

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20160330

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20160427

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160823

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160914

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170214

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170228

R150 Certificate of patent or registration of utility model

Ref document number: 6104166

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350