JP5322548B2 - Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム - Google Patents

Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム Download PDF

Info

Publication number
JP5322548B2
JP5322548B2 JP2008238588A JP2008238588A JP5322548B2 JP 5322548 B2 JP5322548 B2 JP 5322548B2 JP 2008238588 A JP2008238588 A JP 2008238588A JP 2008238588 A JP2008238588 A JP 2008238588A JP 5322548 B2 JP5322548 B2 JP 5322548B2
Authority
JP
Japan
Prior art keywords
perfusion
time
data
interest
parameter
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.)
Active
Application number
JP2008238588A
Other languages
English (en)
Other versions
JP2010068958A (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.)
Toshiba Corp
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems 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 Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Priority to JP2008238588A priority Critical patent/JP5322548B2/ja
Priority to US12/560,047 priority patent/US9561011B2/en
Priority to CN2009101747675A priority patent/CN101675883B/zh
Publication of JP2010068958A publication Critical patent/JP2010068958A/ja
Application granted granted Critical
Publication of JP5322548B2 publication Critical patent/JP5322548B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/507—Clinical applications involving determination of haemodynamic parameters, e.g. perfusion CT

Description

本発明は、被検体から取得した血流のperfusion解析の結果を示す経時的な情報を作成して表示させるX線CT装置、医用画像処理装置および医用画像処理プログラムに関する。
従来から、X線CT (computed tomography)装置、MRI (magnetic resonance imaging)装置、核医学診断装置等の画像診断装置を用いて被検体の様々な臓器において血流のperfusion解析が行われている（例えば特許文献１参照）。従来のperfusion解析は、造影剤やトレーサを被検体に注入して得られるデータを用いて行われる。このperfusion解析によって血流の動態を観察することが可能となる。
perfusion解析の結果の表示方法としては、カラーマップが非常に多く用いられる。カラーマップを用いた観察は、血流動態および虚血部位の診断において積極的に利用されている。特に、頭部の血管や血管から栄養されている心臓や肝臓等の組織におけるperfusion解析は、組織の機能を把握する上で非常に重要な解析である。
図１は、従来のperfusion解析の結果得られる造影剤の時間濃度曲線(TDC: time density curve)を示す図であり、図２は、図１に示すようにTDCとして得られたperfusion解析結果の表示方法を示す図である。
図１の各グラフにおいて、縦軸は、被検体に注入された造影剤の濃度に対応するCT値(ハウンスフィールド値: HU)を示し、横軸は時間を示す。また、図１中の各曲線は、各ROI (ROI1, ROI2, ROI3, …)におけるTDCを表している。
従来のperfusion解析では、造影剤またはトレーサを注入して得られた複数の関心領域(ROI: region of interest)のX線CTデータから図１に示すようなROIごとのTDCが作成される。そして、作成されたTDCを基に様々なアルゴリズムによってROIごとのパーフュージョンパラメータが計算される。
perfusion解析の結果として得られたパーフュージョンパラメータは、図２に示すようにデータ値に応じて色分けされてカラーマップとして表示される。
特開２００６−２４７３８８号公報
しかしながら従来のperfusion解析の結果の表示方法であるカラーマップは、全ての時相のパーフュージョンデータを包括した結果として非ダイナミックに表示するものである。すなわち、カラーマップは、被検体から得られたパーフュージョンデータを時間軸方向に全て重ね合わせた結果を表示するものである。
このため、血流がどの方向からどのように流れてきたか等のパーフュージョンデータの経時的な情報を観察することができず、直感的に血流動態を把握することが困難になるという問題がある。
本発明はかかる従来の事情に対処するためになされたものであり、被検体から取得した血流のperfusion解析の結果を経時的な情報として表示させることが可能なX線CT装置、医用画像処理装置および医用画像処理プログラムを提供することを目的とする。
本発明に係るＸ線ＣＴ装置は、上述の目的を達成するために、被検体に造影剤を注入した状態でＸ線を曝射することによって検出されたＸ線検出データから前記被検体のパーフュージョンデータを生成するパーフュージョンデータ取得手段と、前記パーフュージョンデータから得られる複数の関心領域ごとの前記造影剤の濃度の時間変化に基づいて、前記造影剤の注入から各関心領域における濃度ピークまでの時間を示すパラメータであるパーフュージョンパラメータを前記複数の関心領域ごとに算出するパーフュージョンパラメータ算出手段と、前記造影剤の濃度の時間変化が算出される時間を複数の時間範囲に分割し、各時間範囲に対応する画像であって、前記複数のパーフュージョンパラメータにそれぞれ基づく複数の濃度ピーク時刻のうち前記各時間範囲に含まれる少なくとも１つの濃度ピーク時刻を抽出してその濃度ピーク時刻に対応する関心領域のみが示される画像を生成し、前記画像によって構成される複数の画像を順次表示させるパーフュージョンパラメータ表示手段と、を備えるものである。
また、本発明に係る医用画像処理装置は、上述の目的を達成するために、画像診断装置から被検体のパーフュージョンデータを取得するパーフュージョンデータ取得手段と、前記パーフュージョンデータから得られる複数の関心領域ごとのデータの値の時間変化に基づいて、前記造影剤の注入から各関心領域における濃度ピークまでの時間を示すパラメータであるパーフュージョンパラメータを前記複数の関心領域ごとに算出するパーフュージョンパラメータ算出手段と、前記データの値の時間変化が算出される時間を複数の時間範囲に分割し、各時間範囲に対応する画像であって、前記複数のパーフュージョンパラメータにそれぞれ基づく複数の濃度ピーク時刻のうち前記各時間範囲に含まれる少なくとも１つの濃度ピーク時刻を抽出してその濃度ピーク時刻に対応する関心領域のみが示される画像を生成し、前記画像によって構成される複数の画像を順次表示させるパーフュージョンパラメータ表示手段と、を備えるものである。
また、本発明に係る医用画像処理プログラムは、上述の目的を達成するために、コンピュータを、画像診断装置から被検体のパーフュージョンデータを取得するパーフュージョンデータ取得手段、前記パーフュージョンデータから得られる複数の関心領域ごとのデータの値の時間変化に基づいて、前記造影剤の注入から各関心領域における濃度ピークまでの時間を示すパラメータであるパーフュージョンパラメータを前記複数の関心領域ごとに算出するパーフュージョンパラメータ算出手段、および前記データの値の時間変化が算出される時間を複数の時間範囲に分割し、各時間範囲に対応する画像であって、前記複数のパーフュージョンパラメータにそれぞれ基づく複数の濃度ピーク時刻のうち前記各時間範囲に含まれる少なくとも１つの濃度ピーク時刻を抽出してその濃度ピーク時刻に対応する関心領域のみが示される画像を生成し、前記画像によって構成される複数の画像を順次表示させるパーフュージョンパラメータ表示手段、として機能させるものである。
本発明に係るX線CT装置、医用画像処理装置および医用画像処理プログラムにおいては、被検体から取得した血流のperfusion解析の結果を経時的な情報として表示させることができる。
本発明に係るX線CT装置、医用画像処理装置および医用画像処理プログラムの実施の形態について添付図面を参照して説明する。
（構成および機能）
図３は本発明に係るX線CT装置の実施の形態を示す構成図である。
Ｘ線CT装置１は、ガントリ部２およびコンピュータ装置３とから構成される。ガントリ部２は、Ｘ線管４、高電圧発生装置５、Ｘ線検出器６、データ収集部（DAS: data acquisition system）７、造影剤注入装置８および心電計９を有する。Ｘ線管４とＸ線検出器６とは、高速で且つ連続的に回転する図示しない回転リングに被検体Ｐを挟んで互いに対向する位置に搭載される。
造影剤注入装置８は、コンピュータ装置３からの制御信号により制御され、所定の条件に従って被検体Ｐに造影剤を注入する機能を有する。心電計９は、被検体Ｐに接着させた電極と接続される。心電計９は、被検体Ｐから電極を介して心電(ECG: electro cardiogram)信号を検出し、検出したECG信号から被検体Ｐの心電図を生成してコンピュータ装置３に与える機能を有する。
そして、高電圧発生装置５は、コンピュータ装置３からの制御信号により心電計９において生成された心電図に同期してＸ線管４に所定の管電圧の管電流を供給し、Ｘ線管４から所望のエネルギのＸ線を被検体Ｐに曝射できるように構成されている。Ｘ線検出器６は、２次元状に配置した複数のＸ線検出素子を備えている。そして、３６０度方向から被検体Ｐを透過したＸ線は、Ｘ線検出器６の各Ｘ線検出素子によって検出される。
DAS７は、Ｘ線検出器６の各Ｘ線検出素子から出力されたＸ線検出データに対して増幅、積分処理、A/D (analog to digital)変換処理、対数変換処理等のDAS処理を施すことによって各Ｘ線検出素子にそれぞれ対応する投影データを生成する機能と、生成した投影データをコンピュータ装置３に出力する機能を備えている。
コンピュータ装置３には、表示装置１０および入力装置１１が接続される。また、コンピュータ装置３の記憶装置にプログラム読み込ませて演算装置で実行することにより、コンピュータ装置３には各種機能が備えられる。ただし、コンピュータ装置３に機能を設ける代わりに対応する機能を有する回路を設けてもよい。具体的には、コンピュータ装置３は、Ｘ線検出器６からDAS７を経て取得した投影データから被検体Ｐの血流のパーフュージョンデータを求める機能と、求めたパーフュージョンデータの経時的な情報を表示装置１０に表示させる機能を有する。
そのために、コンピュータ装置３は、スキャン制御装置１２、パーフュージョンデータ生成部１３および医用画像処理装置１４として機能する。医用画像処理装置１４は、時間−濃度曲線作成部１５、パーフュージョンパラメータ算出部１６、時相分割部１７、マスク作成部１８およびパーフュージョンデータ表示処理部１９を有する。
スキャン制御装置１２は、心電計９からの心電図に基づいて、Ｘ線管４およびＸ線検出器６を回転させるためのガントリ部２の架台駆動装置、高電圧発生装置５、造影剤注入装置８にスキャン条件に従って制御信号を与えて制御する機能を有する。
パーフュージョンデータ生成部１３は、DAS７からの出力データである投影データを用いて血流のパーフュージョンデータとして造影X線CT画像データを生成する機能を有する。パーフュージョンデータは、血流部分に対応するROIごとの各時相におけるCT値またはCT値から得られる値を示すダイナミック造影X線CT画像データとして生成される。画像データの１ピクセルまたは隣接する複数のピクセルを１つのROIとすることもできる。各ピクセルをそれぞれROIとして設定すれば、より細部に亘る詳細な血流の動態表示が可能となる。
医用画像処理装置１４は、パーフュージョンデータ生成部１３において生成されたROIごとのパーフュージョンデータから時間に関するperfusionパラメータを算出し、perfusionパラメータの値がある時間幅の範囲にあるROIのみが選択的に描出されるようにマスクされた画像データを時間幅の範囲を変えつつ表示装置１０に順次表示させる機能を有する。
そのために、時間−濃度曲線作成部１５は、パーフュージョンデータ生成部１３において得られたROIごとのCT値またはCT値から得られる値の複数時相分のダイナミック造影CT画像データを取得して、各ROIにおけるTDCを作成する機能を有する。
尚、TDCの作成に用いる画像データのノイズを低減させるために、量子ノイズ低減フィルタ等のフィルタを画像データに対して施すようにしてもよい。また、TDCの平滑化のためにスプライン補間やγ関数を用いてTDCの曲線近似を行ってもよい。さらに、データ処理量の低減化のために、perfusionパラメータを算出すべき部分のみを抽出する処理をTDCの算出に先立ってダイナミック造影CT画像データに対して行ってもよい。この場合、データの抽出方法としては、動的輪郭線法やリジョングローイング法等の手法が知られている。
パーフュージョンパラメータ算出部１６は、時間−濃度曲線作成部１５から取得したTDCから各ROIにおける時間に関するperfusionパラメータを算出する機能を有する。時間に関するperfusionパラメータの例としては、TTP (time to peak)等のように時間を値とし、経時的に算出することが可能なパラメータが挙げられる。
時相分割部１７は、時間−濃度曲線作成部１５から取得したTDCの算出対象となった時間を予め任意に決定されたある時間範囲（間隔）で分割することにより複数の時相に対応する時間範囲を作成する機能を有する。
マスク作成部１８は、時相分割部１７から複数の時相に対応する時間範囲情報を、パーフュージョンパラメータ算出部１６から各ROIにおけるTTP等のperfusionパラメータを、それぞれ取得する機能と、複数の時相に対応する時間範囲のうちTTP等のperfusionパラメータの値が存在する時間範囲以外の時間範囲をマスクするためのマスク情報を作成する機能を有する。換言すれば、マスク作成部１８は、時間範囲ごとの複数のROIを含む画像データを作成する場合に、ある時間範囲内においてperfusionパラメータの値が存在するROIのみを抽出するためのマスクを作成する機能を備えている。マスクは、ユーザが容易に視認できるようなものであればよく、カラーマスクであってもモノクロのマスクであってもよい。
パーフュージョンデータ表示処理部１９は、マスク作成部１８において作成された複数の時間範囲ごとのマスクを用いて、複数の時間範囲ごとのマスク画像データを生成する機能と、生成した時系列の複数のマスク画像データを表示装置１０に順次表示させる機能を有する。マスクは、パーフュージョンパラメータの値が存在する時間範囲およびROIのみが選択的に抽出されるように作成されているため。時系列のマスク画像データは、各時間範囲（各時相）において、パーフュージョンパラメータの値が存在するROIを経時的に識別表示させるダイナミック画像データとなる。このため、ユーザは容易に各時間範囲においてパーフュージョンパラメータの値が存在するROIを視認することが可能となる。
図４は、図３に示す医用画像処理装置１４における時系列のマスク画像データの生成方法を説明する図である。
図４の各グラフにおいて、縦軸は被検体Ｐに注入された造影剤の濃度に対応するCT値(HU)を示し、横軸は時間を、斜め方向の軸は異なる複数のROI方向を、それぞれ示す。また、図４中の各曲線は、N個のROI (ROI1, ROI2, ROI3, …, ROI N)におけるTDC (ROI i)(iはN以下の自然数)をそれぞれ表している。
図４に示すように、時間−濃度曲線作成部１５によってROIごとのTDC (TDC(ROI1), TDC(ROI2), TDC(ROI3), …, TDC(ROIN))が作成される。また、パーフュージョンパラメータ算出部１６により、perfusionパラメータとしてROIごとのTTP (TTP(ROI1), TTP(ROI2), TTP(ROI3), …, TTP(ROIN))が算出される。一方、時相分割部１７によって各TDCの算出対象となった時間がある時間範囲DIVIDED TIME RANGEでそれぞれ分割される。
次に、マスク作成部１８において、各時間範囲においてTTPが存在するROI以外のROIをマスクするためのマスク情報MASK1, MASK2, MASK3, …が時間範囲ごとに作成される。例えば、ROI1のTTPが存在する時間範囲に対応するマスクは、ROI1を抽出するためのMASK3となる。そして、パーフュージョンデータ表示処理部１９は、マスク作成部１８において作成された複数の時間範囲ごとのマスクMASK1, MASK2, MASK3, …を用いて、複数の時相TIME PHASE 1, TIME PHASE 2, TIME PHASE 3, …にそれぞれ対応する時系列の複数のマスク画像データIMAGE (TIME PHASE 1), IMAGE (TIME PHASE 2), IMAGE (TIME PHASE 3), …を生成する。
尚、図４の例では、時間範囲ごとに異なるカラーが割当てられているため、マスク画像データは時相ごとに異なるカラーで表示されている。このようにカラースケールを用いてTTP等のperfusionパラメータを表示させることができる。ただし、グレースケールによりperfusionパラメータを輝度表示させてもよい。
パーフュージョンデータ表示処理部１９は、このように生成したマスク画像データの表示処理を行う機能も備えている。例えば、時系列の複数のマスク画像データが経時的に残像を残しながら表示されるようにする表示処理、残像を残さずに複数のマスク画像データの切換表示を行うための表示処理、時間的に隣接する２つのマスク画像データ間において補間処理を行って中間的なマスク画像データを作成することにより時間的にスムーズなダイナミック表示を行う表示処理、マスク画像データと３次元X線CT画像データとを重畳表示または同期表示させる表示処理、時間的に隣接する３次元X線CT画像データ間におけるサブトラクション画像データとマスク画像データ間とを重畳表示または同期表示させる表示処理等の様々な表示処理をマスク画像データに対して施すことができる。特に、３次元X線CT画像データのサブトラクション画像データは血管を強調させた画像データとなるため、マスク画像データとを同期させて表示すれば、組織にどのように血液が流れているのかが把握できるようになるという効果を期待することができる。
また、１スライスのみについての時系列の２次元マスク画像データではなく、３次元的に設定された複数のROIについての時系列の３次元マスク画像データを生成して、断面変換処理等の画像処理によって生成された２次元画像を表示装置１０に表示させることが実用的である。断面変換処理等の画像処理もパーフュージョンデータ表示処理部１９において行うことができる。
（動作および作用）
次にX線CT装置１の動作および作用について説明する。
図５は、図１に示すＸ線CT装置１により被検体Ｐの複数のROIについてパーフュージョンデータを取得してTTPが存在するROIのみを識別表示させる時系列の複数のマスク画像データを経時的な情報として順次表示させるための流れを示すフローチャートであり、図中Ｓに数字を付した符号はフローチャートの各ステップを示す。
まずステップＳ１において、スキャンの実行により被検体Ｐのパーフュージョンデータが収集される。すなわち、スキャン制御装置１２からの制御信号によって造影剤注入装置８から被検体Ｐに造影剤が注入される。一方、心電計９において検出されたECG信号に基づくスキャン制御装置１２からの制御信号によってＸ線管４およびＸ線検出器６が被検体Ｐ周りに回転し、高電圧発生装置５からECG信号に同期してＸ線管４に管電流が供給される。このため、Ｘ線管４からＸ線が被検体Ｐに曝射され、被検体Ｐを透過したＸ線がＸ線検出器６のＸ線検出素子により検出される。検出されたＸ線検出データは、DAS７に出力され、Ｘ線検出データから投影データが生成される。そして、投影データは、コンピュータ装置３に出力される。
次に、コンピュータ装置３のパーフュージョンデータ生成部１３は、DAS７から出力された投影データからROIごとの各時相における血流のパーフュージョンデータとして造影X線CT画像データを生成する。
次に、ステップＳ２において、医用画像処理装置１４の時間−濃度曲線作成部１５はパーフュージョンデータに基づいて各ROIにおけるTDCを作成する。
次に、ステップＳ３において、パーフュージョンパラメータ算出部１６は、各ROIにおけるTDCから時間に関するperfusionパラメータとして各ROIにおけるTTPを算出する。
一方、ステップＳ４において、時相分割部１７は、TDCの算出対象となった時間をある時間範囲でそれぞれ分割する。これにより、異なる時相に対応する複数の時間範囲が作成される。
次に、ステップＳ５において、マスク作成部１８は、各時相に対応する時間範囲内においてTTPが存在するROIのみを抽出するためのマスクを作成する。
次に、ステップＳ６において、パーフュージョンデータ表示処理部１９は、マスク作成部１８において作成された時間範囲ごとのマスクを用いてTTPが存在しないROIがマスクされた各時相に対応する複数のマスク画像データを経時的なパーフュージョンデータ情報として生成し、複数のマスク画像データを順次表示装置１０に表示させる。
図６は、図３に示すX線CT装置１により経時的に表示される複数時相のマスク画像の例を示す図である。
図６の各画像は、被検体Ｐの頭部における同一断面の画像である。図６に示すように経時的なパーフュージョンデータ情報を表す複数のマスク画像データを時相順に表示させることができる。図６は、各ピクセルをそれぞれROIに設定して１秒ごとに画像を生成した例を示している。このように、時系列のマスク画像データをダイナミック表示させることにより、ピクセル間における血流の移動を表示させることができる。つまり経時的にマスク画像データを表示させることにより、マスク領域および非マスク領域の移動により血流の流れる方向を視覚的に観察することができる。
以上のようなX線CT装置１は、被検体の各ROIにおいて取得した血流のperfusion解析の結果を示すTDCからTTP等の時間をデータ値とするperfusionパラメータを算出し、perfusionパラメータの値がある時間幅の範囲にあるROIのみが選択的に抽出されるようにマスクされた画像を時間幅の範囲を変えつつ順次表示させることによって血流動態を経時的に観察できるようにしたものである。
（効果）
このためX線CT装置１によれば、perfusion解析の結果を経時的に観察することが可能となる。また、血流がどの方向から流れてきたか等の血流の向きを含む血流動態を直感的に把握することも容易となる。特に、頭部や心臓等の部位においては、側副血行路によって栄養されている領域の認識が容易となることが期待される。さらに、３次元X線CT画像データや３次元ダイナミックX線CT画像データ間のサブトラクション画像データと各時相に対応するTTP等のパーフュージョンデータとを重畳表示させることにより、梗塞時における責任血管や支配領域を、より直感的に把握することが可能となる。
（変形例）
上述した実施形態では、被検体Ｐのパーフュージョンデータからperfusion解析の結果を経時的に表す画像データを生成して表示する医用画像処理装置１４をX線CT装置１に内蔵した例を示したが、X線CT装置１と医用画像処理装置１４とをネットワークを介して接続することもできる。さらに、X線CT装置１以外のMRI (Magnetic Resonance Imaging)装置、PET (positron emission computed tomography)やSPECT (single photon emission computed tomography)等の核医学診断装置、超音波診断装置等の画像診断装置に、画像診断装置において取得されたパーフュージョンデータからperfusion解析の結果を経時的に表す画像データを生成して表示する機能を有する医用画像処理装置１４を内蔵またはネットワークを介して接続することもできる。
また、上述した実施形態では、造影剤を用いた例について説明したが、造影剤を用いずに非造影で収集されたパーフュージョンデータからperfusion解析の結果を経時的に表す画像データを生成するようにしても良い。この場合には、パーフュージョンデータを表すMRエコーデータ等の収集データの信号強度の時間変化曲線から時間に関するパーフュージョンパラメータを算出することができる。
従来のperfusion解析の結果得られる造影剤の時間濃度曲線を示す図。 図１に示すようにTDCとして得られたperfusion解析結果の表示方法を示す図。 本発明に係るX線CT装置の実施の形態を示す構成図。 図３に示す医用画像処理装置における時系列のマスク画像データの生成方法を説明する図。 図１に示すＸ線CT装置により被検体の複数のROIについてパーフュージョンデータを取得してTTPが存在するROIのみを識別表示させる時系列の複数のマスク画像データを経時的な情報として順次表示させるための流れを示すフローチャート。 図３に示すX線CT装置１により経時的に表示される複数時相のマスク画像の例を示す図。
符号の説明
１ Ｘ線CT装置
２ ガントリ部
３ コンピュータ装置
４ Ｘ線管
５ 高電圧発生装置
６ Ｘ線検出器
７ データ収集部
８ 造影剤注入装置
９ 心電計
１０ 表示装置
１１ 入力装置
１２ スキャン制御装置
１３ パーフュージョンデータ生成部
１４ 医用画像処理装置
１５ 時間−濃度曲線作成部
１６ パーフュージョンパラメータ算出部
１７ 時相分割部
１８ マスク作成部
１９ パーフュージョンデータ表示処理部
Ｐ 被検体

Claims (4)

  1. 被検体に造影剤を注入した状態でX線を曝射することによって検出されたX線検出データから前記被検体のパーフュージョンデータを生成するパーフュージョンデータ取得手段と、
    前記パーフュージョンデータから得られる複数の関心領域ごとの前記造影剤の濃度の時間変化に基づいて、前記造影剤の注入から各関心領域における濃度ピークまでの時間を示すパラメータであるパーフュージョンパラメータを前記複数の関心領域ごとに算出するパーフュージョンパラメータ算出手段と、
    前記造影剤の濃度の時間変化が算出される時間を複数の時間範囲に分割し、各時間範囲に対応する画像であって、前記複数のパーフュージョンパラメータにそれぞれ基づく複数の濃度ピーク時刻のうち前記各時間範囲に含まれる少なくとも１つの濃度ピーク時刻を抽出してその濃度ピーク時刻に対応する関心領域のみが示される画像を生成し、前記画像によって構成される複数の画像を順次表示させるパーフュージョンパラメータ表示手段と、
    を備えるＸ線ＣＴ装置。
  2. 前記パーフュージョンパラメータ算出手段は、着目する領域内の各ピクセルをそれぞれ１つの関心領域として前記パーフュージョンパラメータを算出するように構成される請求項１記載のＸ線ＣＴ装置。
  3. 画像診断装置から被検体のパーフュージョンデータを取得するパーフュージョンデータ取得手段と、
    前記パーフュージョンデータから得られる複数の関心領域ごとのデータの値の時間変化に基づいて、前記造影剤の注入から各関心領域における濃度ピークまでの時間を示すパラメータであるパーフュージョンパラメータを前記複数の関心領域ごとに算出するパーフュージョンパラメータ算出手段と、
    前記データの値の時間変化が算出される時間を複数の時間範囲に分割し、各時間範囲に対応する画像であって、前記複数のパーフュージョンパラメータにそれぞれ基づく複数の濃度ピーク時刻のうち前記各時間範囲に含まれる少なくとも１つの濃度ピーク時刻を抽出してその濃度ピーク時刻に対応する関心領域のみが示される画像を生成し、前記画像によって構成される複数の画像を順次表示させるパーフュージョンパラメータ表示手段と、
    を備える医用画像処理装置。
  4. コンピュータを、
    画像診断装置から被検体のパーフュージョンデータを取得するパーフュージョンデータ取得手段、
    前記パーフュージョンデータから得られる複数の関心領域ごとのデータの値の時間変化に基づいて、前記造影剤の注入から各関心領域における濃度ピークまでの時間を示すパラメータであるパーフュージョンパラメータを前記複数の関心領域ごとに算出するパーフュージョンパラメータ算出手段、および
    前記データの値の時間変化が算出される時間を複数の時間範囲に分割し、各時間範囲に対応する画像であって、前記複数のパーフュージョンパラメータにそれぞれ基づく複数の濃度ピーク時刻のうち前記各時間範囲に含まれる少なくとも１つの濃度ピーク時刻を抽出してその濃度ピーク時刻に対応する関心領域のみが示される画像を生成し、前記画像によって構成される複数の画像を順次表示させるパーフュージョンパラメータ表示手段、
    として機能させる医用画像処理プログラム。
JP2008238588A 2008-09-17 2008-09-17 Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム Active JP5322548B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2008238588A JP5322548B2 (ja) 2008-09-17 2008-09-17 Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム
US12/560,047 US9561011B2 (en) 2008-09-17 2009-09-15 X-ray CT apparatus, medical image processing apparatus and medical image processing method
CN2009101747675A CN101675883B (zh) 2008-09-17 2009-09-17 X射线ct装置、医用图像处理装置和医用图像处理程序方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008238588A JP5322548B2 (ja) 2008-09-17 2008-09-17 Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム

Publications (2)

Publication Number Publication Date
JP2010068958A JP2010068958A (ja) 2010-04-02
JP5322548B2 true JP5322548B2 (ja) 2013-10-23

Family

ID=42007256

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008238588A Active JP5322548B2 (ja) 2008-09-17 2008-09-17 Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム

Country Status (3)

Country Link
US (1) US9561011B2 (ja)
JP (1) JP5322548B2 (ja)
CN (1) CN101675883B (ja)

Families Citing this family (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8270564B2 (en) * 2008-12-01 2012-09-18 Teratech Corporation Digital integration with detector correction
WO2012029928A1 (ja) * 2010-09-01 2012-03-08 株式会社 東芝 医用画像処理装置
EP2442249B1 (en) 2010-09-28 2019-11-20 Technische Universiteit Eindhoven Perfusion scanning detects angiogenesis from similarity in evolution of local concentrations of contrast agent
JP5660607B2 (ja) * 2010-10-28 2015-01-28 学校法人藤田学園 画像処理装置及びプログラム
US9226673B2 (en) 2011-01-10 2016-01-05 East Carolina University Methods, systems and computer program products for non-invasive determination of blood flow distribution using speckle imaging techniques and hemodynamic modeling
WO2012096878A2 (en) 2011-01-10 2012-07-19 East Carolina University Methods, systems and computer program products for non-invasive determination of blood flow distribution using speckle imaging techniques and hemodynamic modeling
WO2012174495A2 (en) * 2011-06-17 2012-12-20 Carnegie Mellon University Physics based image processing and evaluation process of perfusion images from radiology imaging
US8891845B2 (en) * 2011-09-28 2014-11-18 Telesystems Co., Ltd. Image processing apparatus and image processing method
US9299142B2 (en) * 2011-10-24 2016-03-29 Koninklijke Philips N.V. Perfusion imaging
WO2013111813A1 (ja) * 2012-01-27 2013-08-01 株式会社 東芝 医用画像処理装置
RU2015101806A (ru) * 2012-06-22 2016-08-10 Конинклейке Филипс Н.В. Временная разметка анатомической цели в ангиограммах
WO2014034940A1 (ja) 2012-08-31 2014-03-06 株式会社 東芝 医用診断画像処理装置
DE102012215721A1 (de) * 2012-09-05 2014-03-06 Siemens Aktiengesellschaft Vorrichtung, Verfahren und System zum Betreiben von bildgebenden Verfahren und Systemen
JP2014100555A (ja) * 2012-10-23 2014-06-05 Toshiba Corp 画像処理装置、医用画像診断装置及び画像処理プログラム
US10905391B2 (en) * 2012-11-23 2021-02-02 Imagia Healthcare Inc. Method and system for displaying to a user a transition between a first rendered projection and a second rendered projection
CN103839249B (zh) * 2012-11-23 2017-02-08 上海联影医疗科技有限公司 Ct肝灌注的图像后处理方法和ct肝灌注方法
CN105103191B (zh) 2013-04-03 2019-05-10 皇家飞利浦有限公司 血管分割
WO2014185424A1 (ja) 2013-05-13 2014-11-20 株式会社 東芝 医用画像解析装置
JP6054245B2 (ja) * 2013-05-13 2016-12-27 東芝メディカルシステムズ株式会社 医用画像解析装置
JP6054244B2 (ja) * 2013-05-13 2016-12-27 東芝メディカルシステムズ株式会社 医用画像解析装置
JP6100603B2 (ja) * 2013-05-13 2017-03-22 東芝メディカルシステムズ株式会社 医用画像撮影解析装置
WO2014199995A1 (ja) * 2013-06-11 2014-12-18 株式会社 東芝 Ｘ線コンピュータ断層撮影装置、およびスキャン開始タイミング決定方法
JP6566671B2 (ja) * 2014-04-10 2019-08-28 キヤノンメディカルシステムズ株式会社 Ｘ線診断装置及び画像処理装置
US9691433B2 (en) * 2014-04-18 2017-06-27 Toshiba Medical Systems Corporation Medical image diagnosis apparatus and medical image proccessing apparatus
CN105022719B (zh) * 2014-04-23 2019-06-28 Ge医疗系统环球技术有限公司 医学造影系统及方法
JP6537797B2 (ja) * 2014-09-29 2019-07-03 キヤノンメディカルシステムズ株式会社 医用画像処理装置及びｘ線診断装置
CA2963861C (en) 2014-10-14 2023-10-03 East Carolina University Methods, systems and computer program products for visualizing anatomical structures and blood flow and perfusion physiology using imaging techniques
US11553844B2 (en) 2014-10-14 2023-01-17 East Carolina University Methods, systems and computer program products for calculating MetaKG signals for regions having multiple sets of optical characteristics
CN107257655B (zh) 2014-10-14 2020-06-16 东卡罗莱娜大学 用于利用从多谱段血液流动和灌注成像获取的信号确定血液动力学状态参数的方法、系统和计算机程序产品
US10413253B2 (en) * 2014-11-21 2019-09-17 Samsung Electronics Co., Ltd. Method and apparatus for processing medical image
DE102015202494B4 (de) * 2015-02-12 2023-10-19 Siemens Healthcare Gmbh Auswertung einer dynamischen Kontrastmittelverteilung
US10818073B2 (en) 2015-03-10 2020-10-27 Wisconsin Alumni Research Foundation System and method for time-resolved, three-dimensional angiography with flow information
US10058256B2 (en) 2015-03-20 2018-08-28 East Carolina University Multi-spectral laser imaging (MSLI) methods and systems for blood flow and perfusion imaging and quantification
US10390718B2 (en) 2015-03-20 2019-08-27 East Carolina University Multi-spectral physiologic visualization (MSPV) using laser imaging methods and systems for blood flow and perfusion imaging and quantification in an endoscopic design
DE102017201330B4 (de) * 2017-01-27 2018-10-04 Siemens Healthcare Gmbh Verfahren zur Bildbearbeitung eines vierdimensionalen Angiographiedatensatzes, Bildverarbeitungseinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
JP7084193B2 (ja) 2018-04-10 2022-06-14 ザイオソフト株式会社 医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム
CN111429456B (zh) * 2020-05-06 2023-09-05 上海联影医疗科技股份有限公司 一种图像数据处理方法、装置、电子设备及存储介质
CN112862916B (zh) * 2021-03-11 2021-09-10 首都医科大学附属北京天坛医院 Ct灌注功能图量化参数处理设备及方法
JP2023032314A (ja) * 2021-08-26 2023-03-09 キヤノンメディカルシステムズ株式会社 核医学診断装置、データ処理方法及びプログラム

Family Cites Families (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7069068B1 (en) * 1999-03-26 2006-06-27 Oestergaard Leif Method for determining haemodynamic indices by use of tomographic data
US6516214B1 (en) * 2000-01-24 2003-02-04 The General Hospital Corporation Detection of stroke events using diffuse optical tomography
US6650928B1 (en) * 2000-11-27 2003-11-18 Ge Medical Systems Global Technology Company, Llc Color parametric and composite maps for CT perfusion
EP1302163B1 (en) * 2001-10-16 2006-07-05 Kabushiki Kaisha Toshiba Apparatus for calculating an index of local blood flows
JP3863414B2 (ja) * 2001-11-22 2006-12-27 株式会社東芝 超音波診断装置
US6692438B2 (en) * 2001-12-18 2004-02-17 Koninklijke Philips Electronics Nv Ultrasonic imaging system and method for displaying tissue perfusion and other parameters varying with time
JP4305720B2 (ja) * 2002-11-14 2009-07-29 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置
US6891918B2 (en) * 2002-11-27 2005-05-10 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for acquiring perfusion data
US7171257B2 (en) * 2003-06-11 2007-01-30 Accuray Incorporated Apparatus and method for radiosurgery
DE602005014214D1 (ja) * 2004-03-04 2009-06-10 Philips Intellectual Property
US7912269B2 (en) * 2004-03-31 2011-03-22 Kabushiki Kaisha Toshiba Medical image processing apparatus and method of processing medical image
WO2005114229A2 (en) * 2004-05-14 2005-12-01 The General Hospital Corporation Perfusion weighted mri with local arterial input functions
JP4795694B2 (ja) * 2004-06-11 2011-10-19 株式会社東芝 Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
JP2006130129A (ja) * 2004-11-08 2006-05-25 Toshiba Corp Ｘ線画像診断装置
US7853309B2 (en) * 2005-03-10 2010-12-14 Toshiba Medical Systems Corporation X-ray CT apparatus and myocardial perfusion image generating system
CN100584277C (zh) * 2005-04-14 2010-01-27 伯拉考开发股份有限公司 基于动画灌注成像的灌注评价
US20070016016A1 (en) * 2005-05-31 2007-01-18 Gabriel Haras Interactive user assistant for imaging processes
JP4777007B2 (ja) * 2005-08-03 2011-09-21 東芝メディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置
GB0609610D0 (en) * 2006-05-15 2006-06-21 Stiftelsen Universitetsforskni MR perfusion
JP5536974B2 (ja) * 2006-11-08 2014-07-02 株式会社東芝 Ｘ線診断装置及び画像処理装置
DE102007014133B4 (de) * 2007-03-23 2015-10-29 Siemens Aktiengesellschaft Verfahren zur Visualisierung einer Sequenz tomographischer Volumendatensätze der medizinischen Bildgebung
US8285014B2 (en) * 2007-04-06 2012-10-09 Siemens Aktiengesellschaft Measuring blood volume with C-arm computed tomography
US20080262344A1 (en) * 2007-04-23 2008-10-23 Brummett David P Relative value summary perfusion map
US8320647B2 (en) * 2007-11-20 2012-11-27 Olea Medical Method and system for processing multiple series of biological images obtained from a patient
JP5234905B2 (ja) * 2007-11-20 2013-07-10 東芝メディカルシステムズ株式会社 Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
JP4714228B2 (ja) * 2008-01-21 2011-06-29 株式会社東芝 脳組織内毛細血管の血流動態に関するインデックス演算方法、装置及び記憶媒体

Also Published As

Publication number Publication date
JP2010068958A (ja) 2010-04-02
CN101675883A (zh) 2010-03-24
US20100067767A1 (en) 2010-03-18
US9561011B2 (en) 2017-02-07
CN101675883B (zh) 2012-09-19

Similar Documents

Publication Publication Date Title
JP5322548B2 (ja) Ｘ線ｃｔ装置、医用画像処理装置および医用画像処理プログラム
EP1967140B1 (en) X-Ray CT-Apparatus, Myocardial Perfusion Information Generating System, X-Ray Diagnostic Method and Myocardial Perfusion Information Generating Method
JP5558672B2 (ja) 画像処理装置及びｘ線コンピュータ断層撮影装置
EP1604612B1 (en) X-ray CT apparatus for myocardial perfusion image generation
JP5148315B2 (ja) 医用画像処理装置、及び医用画像処理プログラム
JP5053606B2 (ja) 医用画像診断装置及び医用画像処理装置
WO2014034940A1 (ja) 医用診断画像処理装置
JP2018083056A (ja) 医用画像処理装置、医用画像処理方法及びｘ線ｃｔ装置
JP5596099B2 (ja) 医用画像処理装置、及び医用画像処理プログラム
JP2009125127A (ja) Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
JP5388638B2 (ja) 画像処理装置および画像処理プログラム
JP5156863B2 (ja) Ｘ線ｃｔ装置
JP2003204961A (ja) Ｘ線ｃｔ装置
JP2012000135A (ja) マルチモダリティ動画像診断装置
JP2010154982A (ja) Ｘ線コンピュータ断層撮影装置および画像処理装置
JP2006271484A (ja) 画像形成装置及び画像形成方法
JP5537072B2 (ja) 医用画像表示装置およびｘ線ｃｔ装置
JP5215444B2 (ja) 医用画像診断装置及び医用画像処理装置
JP2007289569A (ja) 医用画像処理装置及び医用画像処理方法
JP2013040829A (ja) ボリュームデータ処理装置及び方法
JP7024954B2 (ja) 医用画像処理装置及びｘ線ｃｔ装置
JP2007068726A (ja) 心機能診断装置
CN111340934A (zh) 用于生成组合式组织-血管表示的方法和计算机系统
JP2010022667A (ja) 脳血流解析装置
JP6878074B2 (ja) 医用画像診断装置及び医用画像処理装置

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20100422

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110916

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20111206

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20121127

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130128

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130326

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130527

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: 20130625

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130716

R150 Certificate of patent or registration of utility model

Ref document number: 5322548

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350