JP7341803B2 - X線ct装置、医用画像処理装置及びx線ctシステム - Google Patents

X線ct装置、医用画像処理装置及びx線ctシステム Download PDF

Info

Publication number
JP7341803B2
JP7341803B2 JP2019162621A JP2019162621A JP7341803B2 JP 7341803 B2 JP7341803 B2 JP 7341803B2 JP 2019162621 A JP2019162621 A JP 2019162621A JP 2019162621 A JP2019162621 A JP 2019162621A JP 7341803 B2 JP7341803 B2 JP 7341803B2
Authority
JP
Japan
Prior art keywords
ray
region
image
analysis
scatter diagram
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
JP2019162621A
Other languages
English (en)
Other versions
JP2020039872A (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.)
Canon Medical Systems Corp
National Cancer Center Japan
Original Assignee
Canon Medical Systems Corp
National Cancer Center Japan
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 Canon Medical Systems Corp, National Cancer Center Japan filed Critical Canon Medical Systems Corp
Priority to US16/563,569 priority Critical patent/US11295487B2/en
Priority to CN201910864900.3A priority patent/CN110881993B/zh
Publication of JP2020039872A publication Critical patent/JP2020039872A/ja
Application granted granted Critical
Publication of JP7341803B2 publication Critical patent/JP7341803B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/482Diagnostic techniques involving multiple energy imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/465Displaying means of special interest adapted to display user selection data, e.g. graphical user interface, icons or menus
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/467Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Human Computer Interaction (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Pulmonology (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

本発明の実施形態は、X線CT装置、医用画像処理装置及びX線CTシステムに関する。
従来、X線コンピュータ断層撮影(CT:Computed Tomography)装置により、異なる複数種類の管電圧で撮影を行なって画像を取得する手法がある。例えば、2種類の異なる管電圧を用いるデュアルエナジー(Dual-Energy:DE)収集では、2種類の異なる管電圧から得られた2つの投影データを、予め設定した2つの基準物質それぞれの投影データ(線積分データ)に分離し、分離した2つのデータそれぞれから、基準物質の存在率に基づく画像(基準物質画像)を再構成する技術が知られている。かかる技術では、2つの基準物質画像を用いて重み付け計算処理を行うことにより、単色X線画像や密度画像、実効原子番号画像等、様々な画像を取得することができる。
特開2009-261942号公報 WO2004/089218号公報
本発明が解決しようとする課題は、組織分解能を向上させたCT画像を提供することである。
実施形態のX線CT装置は、スキャン部と、解析部とを備える。スキャン部は、X線を利用したスキャンを実行することで、被検体の一部を含む領域について、第1のX線エネルギーに対応する第1データセット、及び前記第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを収集する。解析部は、前記第1データセット及び前記第2データセットに基づいて、前記領域内の各位置における複数の基準物質それぞれの含有量を表す散布図を生成し、前記散布図の少なくとも一部について重みを設定し、前記含有量及び前記重みに基づいて解析画像を生成する。
図1は、第1の実施形態に係る医用画像処理システムの構成の一例を示す図である。 図2は、第1の実施形態に係るX線CT装置の構成の一例を示すブロック図である。 図3は、第1の実施形態に係る表示情報の一例を示す図である。 図4は、第1の実施形態に係る散布図上の関心領域の指定を説明するための図である。 図5は、第1の実施形態に係る解析機能による重み付けの一例を説明するための図である。 図6は、第1の実施形態に係る解析画像の一例を示す図である。 図7は、第1の実施形態に係るX線CT装置の処理の流れを説明するためのフローチャートである。 図8は、第2の実施形態に係る重み付けの一例を説明するための図である。 図9は、第2の実施形態に係る解析画像の一例を示す図である。 図10は、第4の実施形態に係る解析画像の一例を示す図である。 図11は、第5の実施形態に係る解析機能による処理を説明するための模式図である。 図12は、第5の実施形態に係る医用画像処理装置の構成の一例を示す図である。
以下、図面を参照して、X線CT装置、医用画像処理装置及びX線CTシステムの実施形態について詳細に説明する。なお、本願に係るX線CT装置、医用画像処理装置及びX線CTシステムは、以下に示す実施形態によって限定されるものではない。
(第1の実施形態)
まず、第1の実施形態について説明する。第1の実施形態では、本願が開示する技術をX線CT装置に適用した場合の例を説明する。なお、以下、X線CT装置を含む医用画像処理システムを例に挙げて説明する。
図1は、第1の実施形態に係る医用画像処理システムの構成の一例を示す図である。図1に示すように、第1の実施形態に係る医用画像処理システム100は、X線CT(Computed Tomography)装置1と、画像保管装置2と、医用画像処理装置3とを備え、ネットワーク200を介して相互に接続される。なお、図1に示す例は、あくまでも一例であり、ネットワーク200にその他種々の装置(例えば、端末装置等)が接続される場合であってもよい。
例えば、第1の実施形態に係るX線CT装置1は、図1に示すように、ネットワーク200を介して、画像保管装置2と、医用画像処理装置3に接続される。なお、医用画像処理システム100は、ネットワーク200を介して、MRI装置、超音波診断装置及びPET(Positron Emission Tomography)装置等の他の医用画像診断装置や、端末装置などにさらに接続されてもよい。また、ネットワーク200は、院内で閉じたローカルネットワークにより構成されてもよいし、インターネットを介したネットワークにより構成されてもよい。
画像保管装置2は、各種の医用画像診断装置によって収集された画像データを保管する。例えば、画像保管装置2は、サーバ装置等のコンピュータ機器によって実現される。本実施形態では、画像保管装置2は、ネットワーク200を介してX線CT装置1から投影データやCT画像データ(ボリュームデータ)を取得し、取得した投影データやCT画像データを装置内又は装置外に設けられた記憶回路に記憶させる。
医用画像処理装置3は、ネットワーク200を介して各種の医用画像診断装置から画像データを取得し、取得した画像データを処理する。例えば、医用画像処理装置3は、ワークステーション等のコンピュータ機器によって実現される。本実施形態では、医用画像処理装置3は、ネットワーク200を介してX線CT装置1又は画像保管装置2から投影データやCT画像データを取得し、取得した投影データやCT画像データに対して各種画像処理を行う。医用画像処理装置3は、画像処理を行う前又は行った後のCT画像データをディスプレイ等に表示する。
X線CT装置1は、被検体のCT画像データを収集する。具体的には、X線CT装置1は、被検体を略中心にX線管及びX線検出器を旋回移動させ、被検体を透過したX線を検出して投影データを収集する。ここで、X線CT装置1は、異なる複数種類の管電圧で撮影を行なって、複数種類の投影データセットを取得する。例えば、X線CT装置1は、2種類の異なる管電圧を用いるデュアルエナジー(Dual-Energy:DE)収集により、1つの撮影対象部位に対して2つの投影データセットを収集する。そして、X線CT装置1は、収集した投影データに基づいて、CT画像データを生成する。
図2は、第1の実施形態に係るX線CT装置1の構成の一例を示す図である。図2に示すように、第1の実施形態に係るX線CT装置1は、架台装置10と、寝台装置30と、コンソール装置40とを有する。
ここで、図2においては、非チルト状態での回転フレーム13の回転軸又は寝台装置30の天板33の長手方向をZ軸方向とする。また、Z軸方向に直交し、床面に対し水平である軸方向をX軸方向とする。また、Z軸方向に直交し、床面に対し垂直である軸方向をY軸方向とする。なお、図2は、説明のために架台装置10を複数方向から描画したものであり、X線CT装置1が架台装置10を1つ有する場合を示す。
架台装置10は、X線管11と、X線検出器12と、回転フレーム13と、X線高電圧装置14と、制御装置15と、ウェッジ16と、コリメータ17と、DAS(Data Acquisition System)18とを有する。
X線管11は、熱電子を発生する陰極(フィラメント)と、熱電子の衝突を受けてX線を発生する陽極(ターゲット)とを有する真空管である。X線管11は、X線高電圧装置14からの高電圧の印加により、陰極から陽極に向けて熱電子を照射することで、被検体Pに対し照射するX線を発生する。例えば、X線管11には、回転する陽極に熱電子を照射することでX線を発生させる回転陽極型のX線管がある。
X線検出器12は、X線を検出する検出素子を複数有する。X線検出器12における各検出素子は、X線管11から照射されて被検体Pを通過したX線を検出し、検出したX線量に対応した信号をDAS18へと出力する。X線検出器12は、例えば、X線管11の焦点を中心とした1つの円弧に沿ってチャンネル方向(チャネル方向)に複数の検出素子が配列された複数の検出素子列を有する。X線検出器12は、例えば、チャンネル方向に複数の検出素子が配列された検出素子列が列方向(スライス方向、row方向)に複数配列された構造を有する。
例えば、X線検出器12は、グリッドと、シンチレータアレイと、光センサアレイとを有する間接変換型の検出器である。シンチレータアレイは、複数のシンチレータを有する。シンチレータは入射X線量に応じた光子量の光を出力するシンチレータ結晶を有する。グリッドは、シンチレータアレイのX線入射側の面に配置され、散乱X線を吸収するX線遮蔽板を有する。なお、グリッドはコリメータ(1次元コリメータ又は2次元コリメータ)と呼ばれる場合もある。光センサアレイは、シンチレータからの光量に応じた電気信号に変換する機能を有し、例えば、フォトダイオード等の光センサを有する。なお、X線検出器12は、入射したX線を電気信号に変換する半導体素子を有する直接変換型の検出器であっても構わない。
回転フレーム13は、X線管11とX線検出器12とを対向支持し、制御装置15によってX線管11とX線検出器12とを回転させる円環状のフレームである。例えば、回転フレーム13は、アルミニウムを材料とした鋳物である。なお、回転フレーム13は、X線管11及びX線検出器12に加えて、X線高電圧装置14やウェッジ16、コリメータ17、DAS18等を更に支持することもできる。更に、回転フレーム13は、図2において図示しない種々の構成を更に支持することもできる。
X線高電圧装置14は、変圧器(トランス)及び整流器等の電気回路を有し、X線管11に印加する高電圧を発生する高電圧発生装置と、X線管11が発生するX線に応じた出力電圧の制御を行うX線制御装置とを有する。高電圧発生装置は、変圧器方式であってもよいし、インバータ方式であってもよい。なお、X線高電圧装置14は、回転フレーム13に設けられてもよいし、図示しない固定フレームに設けられても構わない。
制御装置15は、CPU(Central Processing Unit)等を有する処理回路と、モータ及びアクチュエータ等の駆動機構とを有する。制御装置15は、入力インターフェース43からの入力信号を受けて、架台装置10及び寝台装置30の動作制御を行う。例えば、制御装置15は、回転フレーム13の回転や架台装置10のチルト、寝台装置30及び天板33の動作等について制御を行う。一例を挙げると、制御装置15は、架台装置10をチルトさせる制御として、入力された傾斜角度(チルト角度)情報により、X軸方向に平行な軸を中心に回転フレーム13を回転させる。なお、制御装置15は架台装置10に設けられてもよいし、コンソール装置40に設けられてもよい。
ウェッジ16は、X線管11から照射されたX線量を調節するためのフィルタである。具体的には、ウェッジ16は、X線管11から被検体Pへ照射されるX線が、予め定められた分布になるように、X線管11から照射されたX線を透過して減衰するフィルタである。例えば、ウェッジ16は、ウェッジフィルタ(wedge filter)やボウタイフィルタ(bow-tie filter)であり、所定のターゲット角度や所定の厚みとなるようにアルミニウム等を加工したフィルタである。
コリメータ17は、ウェッジ16を透過したX線の照射範囲を絞り込むための鉛板等であり、複数の鉛板等の組み合わせによってスリットを形成する。なお、コリメータ17は、X線絞りと呼ばれる場合もある。また、図2においては、X線管11とコリメータ17との間にウェッジ16が配置される場合を示すが、X線管11とウェッジ16との間にコリメータ17が配置される場合であってもよい。この場合、ウェッジ16は、X線管11から照射され、コリメータ17により照射範囲が制限されたX線を透過して減衰させる。
DAS18は、X線検出器12が有する各検出素子によって検出されるX線の信号を収集する。例えば、DAS18は、各検出素子から出力される電気信号に対して増幅処理を行う増幅器と、電気信号をデジタル信号に変換するA/D変換器とを有し、検出データを生成する。DAS18は、例えば、プロセッサにより実現される。
DAS18が生成したデータは、回転フレーム13に設けられた発光ダイオード(Light Emitting Diode: LED)を有する送信機から、光通信によって、架台装置10の非回転部分(例えば、固定フレーム等。図1での図示は省略している)に設けられた、フォトダイオードを有する受信機に送信され、コンソール装置40へと転送される。ここで、非回転部分とは、例えば、回転フレーム13を回転可能に支持する固定フレーム等である。なお、回転フレーム13から架台装置10の非回転部分へのデータの送信方法は、光通信に限らず、非接触型の如何なるデータ伝送方式を採用してもよいし、接触型のデータ伝送方式を採用しても構わない。
寝台装置30は、撮影対象の被検体Pを載置、移動させる装置であり、基台31と、寝台駆動装置32と、天板33と、支持フレーム34とを有する。基台31は、支持フレーム34を鉛直方向に移動可能に支持する筐体である。寝台駆動装置32は、被検体Pが載置された天板33を、天板33の長軸方向に移動する駆動機構であり、モータ及びアクチュエータ等を含む。支持フレーム34の上面に設けられた天板33は、被検体Pが載置される板である。なお、寝台駆動装置32は、天板33に加え、支持フレーム34を天板33の長軸方向に移動してもよい。
コンソール装置40は、メモリ41と、ディスプレイ42と、入力インターフェース43と、処理回路44とを有する。なお、コンソール装置40は架台装置10とは別体として説明するが、架台装置10にコンソール装置40又はコンソール装置40の各構成要素の一部が含まれてもよい。
メモリ41は、例えば、RAM(Random Access Memory)、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。メモリ41は、例えば、投影データやCT画像データを記憶する。また、メモリ41は、関心領域の情報を記憶する。なお、関心領域については、後に詳述する。また、例えば、メモリ41は、X線CT装置1に含まれる回路がその機能を実現するためのプログラムを記憶する。なお、メモリ41は、X線CT装置1とネットワークを介して接続されたサーバ群(クラウド)により実現されることとしてもよい。なお、メモリ41は、記憶部の一例である。
ディスプレイ42は、各種の情報を表示する。例えば、ディスプレイ42は、処理回路44によって生成された各種の画像を表示したり、操作者から各種の操作を受け付けるためのGUI(Graphical User Interface)を表示したりする。例えば、ディスプレイ42は、液晶ディスプレイやCRT(Cathode Ray Tube)ディスプレイである。ディスプレイ42は、デスクトップ型でもよいし、コンソール装置40本体と無線通信可能なタブレット端末等で構成されることにしても構わない。
入力インターフェース43は、操作者から各種の入力操作を受け付けて、受け付けた入力操作を電気信号に変換して処理回路44に出力する。例えば、入力インターフェース43は、ディスプレイ42によって表示された情報に対して指定操作を行うための入力操作を操作者から受け付ける。なお、指定操作については後述する。また、例えば、入力インターフェース43は、CT画像データを再構成する際の再構成条件や、CT画像データから後処理画像を生成する際の画像処理条件等の入力操作を操作者から受け付ける。
例えば、入力インターフェース43は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、音声入力回路等により実現される。なお、入力インターフェース43は、架台装置10に設けられてもよい。また、入力インターフェース43は、コンソール装置40本体と無線通信可能なタブレット端末等で構成されることにしても構わない。また、入力インターフェース43は、マウスやキーボード等の物理的な操作部品を備えるものだけに限られない。例えば、コンソール装置40とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を処理回路44へ出力する電気信号の処理回路も入力インターフェース43の例に含まれる。なお、入力インターフェース43は、受け付け部の一例である。
処理回路44は、X線CT装置1全体の動作を制御する。例えば、処理回路44は、スキャン制御機能441、前処理機能442、画像再構成機能443、解析機能444及び表示制御機能445を実行する。なお、スキャン制御機能441は、スキャン部の一例である。また、解析機能444は、解析部の一例である。
図2に示すX線CT装置1においては、各処理機能がコンピュータによって実行可能なプログラムの形態でメモリ41へ記憶されている。処理回路44は、メモリ41からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路44は、読み出したプログラムに対応する機能を有することとなる。
なお、図2においては、スキャン制御機能441、前処理機能442、画像再構成機能443、解析機能444及び表示制御機能445の各処理機能が単一の処理回路44によって実現される場合を示したが、実施形態はこれに限られるものではない。例えば、処理回路44は、複数の独立したプロセッサを組み合わせて構成され、各プロセッサが各プログラムを実行することにより各処理機能を実現するものとしても構わない。また、処理回路44が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。
以上、本実施形態に係るX線CT装置1の全体構成について説明した。かかる構成のもと、本実施形態に係るX線CT装置1は、組織分解能を向上させたCT画像を提供する。具体的には、X線CT装置1は、処理回路44による各種機能によって、複数種類の異なる管電圧での撮影を行なって収集した投影データ或いは投影データから再構成したCT画像データを用いて、組織ごとのコントラストを強調させたCT画像を生成する。
例えば、造影剤を用いない単純CT検査において撮影されたCT画像では、脳や、腫瘍、筋肉などの軟部組織はCT値が近い値となるため、CT値に基づくコントラストの違いによってCT画像における組織を明瞭に区別させることが困難となる場合がある。そこで、本願に係るX線CT装置では、複数種類の異なる管電圧での撮影から得られる基準物質画像を解析することで、画素ごとの基準物質の混合割合を取得し、取得した混合割合に基づいて、組織ごとのコントラストを強調させたCT画像を生成する。
なお、本実施形態では、「複数種類の異なる管電圧での撮影」は、2種類の異なる管電圧での「Dual Energyによる撮影」の他、3種類以上の異なる管電圧での「Multi Energyによる撮影」も含む。ここで、「Dual Energyによる撮影」の場合、以下の4つの撮影方法のいずれが用いられる場合でもよい。
例えば、第1の方法としては、1つのX線管を用いて第1の管電圧で撮影した後に、第1の管電圧と異なる第2の管電圧で撮影する「Slow-kV switching方式(2回転方式)」がある。また、例えば、第2の方法としては、回転中(スキャン中)のビュー毎に高速にX線管の管電圧を切り替えて撮影する「Fast-kV switching方式(高速スイッチング方式)」がある。この場合、管電圧の切り替えに同期してデータ収集装置がデータ収集を行ない、異なる管電圧のデータを1つのスキャン中に収集する。また、例えば、第3の方法としては、1つのX線管ではなく2のX線管を搭載した上でそれらを用いて異なる管電圧で撮影する「Dual Source方式(2管球方式)」がある。また、例えば、第4の方法としては、多層構造のX線検出器を用いる「積層型検出器方式」がある。例えば2層構造(浅い層の検出器、深い層の検出器)のX線検出器を用いる場合、浅い層の検出器で低エネルギーのX線が検出され、浅い層の検出器を通過した深い層の検出器で高エネルギーのX線が検出される。
以下、各種機能の詳細について説明する。なお、以下では、「Dual Energyによる撮影」を「Fast-kV switching方式」で実行することによって2種類の異なる管電圧での撮影を行ない、組織ごとのコントラストを強調させたCT画像を生成する場合について説明する。
スキャン制御機能441は、入力インターフェース43を介して操作者から受け付けた入力操作に基づいて、スキャンを制御する。具体的には、スキャン制御機能441は、入力操作に基づいて、X線高電圧装置14に制御信号を送信することで、高電圧発生装置からの出力電圧を制御する。また、スキャン制御機能441は、DAS18に制御信号を送信することで、DAS18によるデータ収集を制御する。
ここで、スキャン制御機能441は、X線を利用したスキャンを実行することで、被検体の一部を含む領域について、第1のX線エネルギーに対応する第1データセット、及び前記第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを収集する。例えば、スキャン制御機能441は、「Fast-kV switching方式」による「Dual Energyによる撮影」を制御する。
かかる場合には、スキャン制御機能441は、高電圧と低電圧とを切り替える制御信号をX線高電圧装置14に送信することで、X線高電圧装置14からX線管11への高電圧と低電圧の印加を制御する。また、スキャン制御機能441は、DAS18に対して制御信号を送信することで、検出した検出データが、高電圧のX線照射によるものであるか、或いは、低電圧のX線照射によるものであるかを識別させるように制御する。
前処理機能442は、DAS18から送信されたX線検出データに対して、対数変換処理や、オフセット補正、感度補正、ビームハードニング補正等の補正処理を行なうことで、投影データを生成する。なお、前処理前のデータ(検出データ)および前処理後のデータを総称して投影データと称する場合もある。例えば、前処理機能442は、第1の管電圧(例えば、高電圧)の検出データから投影データセット(以下、高エネルギー投影データセットと記載する)を生成する。また、前処理機能442は、第2の管電圧(例えば、低電圧)の検出データから投影データセット(以下、低エネルギー投影データセットと記載する)を生成する。
また、前処理機能442は、2種類の投影データセットを用いて、撮影の対象部位に存在する、予め決定された2つ以上の基準物質(水、ヨード、カルシウム、ハイドロキシアパタイト、脂肪等)を分離する。そして、前処理機能442は、2つ以上の基準物質のそれぞれに対応する2種類以上の単色X線の投影データセットを生成する。例えば、前処理機能442は、高エネルギー投影データセット及び低エネルギー投影データセットから、水とヨードの単色X線の投影データセットをそれぞれ生成する。なお、前処理機能442によって生成された投影データセットは、メモリ41によって記憶される。
画像再構成機能443は、メモリ41によって記憶された投影データセットから各種画像を生成し、生成した画像をメモリ41に格納する。例えば、画像再構成機能443は、投影データを種々の再構成法(例えば、FBP(Filtered Back Projection)などの逆投影法や、逐次近似法など)によって再構成することでCT画像データを再構成し、再構成したCT画像データをメモリ41に格納する。また、画像再構成機能443は、種々の画像処理を行うことにより、CT画像データからMPR画像などのCT画像を生成して、生成したCT画像をメモリ41に格納する。
例えば、画像再構成機能443は、メモリ41によって記憶された基準物質の単色X線の投影データセットを読み出し、基準物質画像データ(基準物質強調画像データ)を再構成する。一例を挙げると、画像再構成機能443は、水成分が強調された投影データセットに基づいて水成分の基準物質画像データを再構成し、ヨード成分が強調された投影データセットに基づいてヨード成分の基準物質画像データを再構成する。また、画像再構成機能443は、水成分の基準物質画像データ及びヨード成分の基準物質画像データに対してそれぞれ画像処理を実行することで、水成分の基準物質画像とヨード成分の基準物質画像を生成する。また、画像再構成機能443は、2つの基準物質画像データを用いて重み付け計算処理を行うことにより、所定のエネルギーにおける単色X線画像や密度画像、実効原子番号画像等、様々な画像を生成することができる。
また、例えば、画像再構成機能443は、メモリ41によって記憶された高エネルギー投影データセットと低エネルギー投影データセットとを読み出し、各投影データセットからCT画像データをそれぞれ再構成する。そして、画像再構成機能443は、CT画像データから高エネルギーに対応する多色X線画像と、低エネルギーに対応する多色X線画像とを生成することもできる。
解析機能444は、メモリ41によって記憶されたCT画像データや、CT画像に対して種々の解析処理を実行する。例えば、解析機能444は、第1データセット(例えば、高エネルギー投影データセット)及び第2データセット(例えば、低エネルギー投影データセット)に基づいて、領域内の各位置における複数の基準物質それぞれの含有量を表す散布図を生成し、散布図の少なくとも一部について重みを設定し、含有量及び重みに基づいて解析画像を生成する。
かかる場合には、解析機能444は、まず、第1データセットと第2データセットに基づいて、画素ごとに、複数の基準物質間での存在比率を推定する。例えば、解析機能444は、画像再構成機能443によって生成された高エネルギーに対応する多色X線画像と低エネルギーに対応する多色X線画像との間の画素値に基づいて、画素ごとに、水及びヨードにおける一方に対する他方の存在比率を推定する。すなわち、解析機能444は、Dual Energyによって撮影された領域の各画素について、水に対するヨードの存在比率と、ヨードに対する水の存在比率とをそれぞれ推定する。
そして、解析機能444は、推定した存在比率を軸とする散布図を生成し、生成した散布図における位置ごとに重みを設定することで、当該位置に対応する画素の輝度値を、重みが反映された輝度値に変換した解析画像を生成する。すなわち、解析機能444は、基準物質の存在比率を軸とした散布図を利用して、組織分解能を向上させたCT画像を生成する。なお、解析画像の生成については、後に詳述する。
表示制御機能445は、画像再構成機能443によって生成されたCT画像や、解析機能444によって生成された散布図、及び、解析画像などをディスプレイ42に表示させる。図3は、第1の実施形態に係る表示情報の一例を示す図である。ここで、図3においては、ヨード成分の基準物質画像と、水成分の基準物質画像と、単色X線画像と、散布図とを示す。例えば、表示制御機能445は、図3に示すように、脳の同一断面におけるヨード成分の基準物質画像(図中の左上)と、水成分の基準物質画像(図中の右上)と、単色X線画像(図中の右下)と、散布図(図中の左下)とをディスプレイ42に表示させるように制御する。
ここで、解析機能444によって生成される散布図は、例えば、図3の左下に示すように、水に対するヨードの存在比率を縦軸に示し、ヨードに対する水の存在比率を横軸に示すグラフによって示される。かかる場合には、例えば、解析機能444は、画素ごとに推定した水に対するヨードの存在比率及びヨードに対する水の存在比率に基づいて、各画素を示す点をグラフ中にプロットした散布図を生成する。すなわち、解析機能444は、画素ごとに、水とヨードの混合割合をグラフ中の座標で示した散布図を生成する。
ここで、図3に示す散布図は、2次元で表されているが、実際には、3次元の情報量を有する。具体的には、図3に示す散布図は、画素ごとの水とヨードの混合割合を示すものである。したがって、同一の組織(同一の物質)の場合、水とヨードの混合割合が略同じになり、同一の組織に含まれる画素は、散布図上の略同一の座標にプロットされることとなる。すなわち、図3に示す散布図は、奥行方向に複数の点が重なっており、縦軸、横軸及び奥行方向の3次元の情報量を有する。
また、図3に示す散布図は、表示制御機能445によって表示されるヨード成分の基準物質画像、水成分の基準物質画像及び単色X線画像と位置関係が対応づいている。すなわち、各画像における画素の位置と、散布図において対応する画素の位置とが対応づいており、例えば、各画像中の領域R1に含まれる画素群の散布図における位置が、矢印a1にて指示される点群に対応する。例えば、医師などの操作者が、入力インターフェース43を介して、いずれかの画像に領域R1を設定すると、表示制御機能445は、他の画像上の対応する位置に領域R1を表示させるとともに、矢印a1にて指示される点群の色を変えて表示させることも可能である。
同様に、表示制御機能445は、図3に示すように、各画像中の領域R2~領域R5に含まれる各画素群に対応する散布図中の点群を、矢印a2~矢印a5にて指示するように色を変えて表示させる。ここで、表示制御機能445は、対応する領域と点群とを同じ色で示し、かつ、領域R1~領域R5(矢印a1~矢印a5)間で異なる色を用いることで、それぞれを識別可能に表示させることができる。
上述したように、X線CT装置1は、散布図を利用して、組織分解能を向上させたCT画像を表示させる。ここで、解析機能444によって生成される散布図は、基準物質の混合割合をグラフ中の座標で示したものである。したがって、同一の組織の画素を示す点群は、散布図上の略同一位置にプロットされると考えられる。しかしながら、実際の散布図は、ノイズの影響などにより、図3に示すように、同一の組織に含まれる画素であっても、位置が分散されることが多い。
そこで、本願では、散布図上の位置に応じた重み付けを設定し、設定した重み付けをCT画像の輝度値に反映させることで、単純CT検査で撮影したCT画像であっても組織ごとのコントラストを強調させ、組織分解能を向上させたCT画像を提供する。
このような処理を行うために、解析機能444は、まず、コントラストを強調させる対象組織の散布図におけるおおよその位置を特定する。例えば、解析機能444は、入力インターフェース43を介して対象組織の位置を示す関心領域を指定する指定操作を受け付けることで、散布図における対象組織のおおよその位置を特定する。図4は、第1の実施形態に係る散布図上の関心領域の指定を説明するための図である。ここで、図4においては、左側に散布図を示し、右側に単色X線画像を示す。
例えば、図4に示すように、表示制御機能445が、ディスプレイ42に散布図と単色X線画像を表示させると、医師などの操作者が、入力インターフェース43を介して、散布図上に関心領域R11~関心領域R51などの領域を設定する。ここで、操作者による関心領域の設定は、散布図上に直接設定する場合でもよく、或いは、単色X線画像の位置を指定することで、各関心領域を設定する場合でもよい。
例えば、操作者が関心領域を散布図上に直接設定すると、表示制御機能445は、設定された関心領域に含まれる点群に対応する画素群を、単色X線画像上で強調して表示させる。一例を挙げると、操作者が関心領域R11を散布図上に直接設定すると、表示制御機能445は、図4の右側の図に示すように、単色X線画像において、関心領域R11に含まれる点群に対応する画素群の領域である領域R111の色を変えて表示させる。これにより、操作者は、散布図上で指定した関心領域に対応する部位を画像上で確認しながら、関心領域を設定することができる。
そして、コントラストを強調させる対象組織が複数ある場合には、例えば、図4に示すように、操作者は、散布図上に関心領域R21~関心領域R51を設定する。この場合にも、上記と同様に、表示制御機能445は、関心領域R21~関心領域R51に対応する画像上の領域R211~領域R511の色をそれぞれ変えて表示させることで、散布図上で指定した関心領域に対応する部位を画像上で確認させることができる。
また、操作者が画像上に領域を設定すると、表示制御機能445は、設定された領域に含まれる画素群に対応する点群を、散布図上で強調して表示させる。一例を挙げると、操作者が領域R111を単色X線画像上に設定すると、表示制御機能445は、散布図において、領域R111に含まれる画素群に対応する点群の色を変えて表示させる。これにより、操作者は、対象組織の画素群に対応する点群の散布図上での位置を把握することができる。操作者は、散布図上で色が変わっている点群を含むように、入力インターフェース43を介して関心領域を設定することで、例えば、領域R111に対応する散布図上の関心領域R11を設定することができる。
同様に、操作者が単色X線画像上に領域R211~領域R511を設定すると、表示制御機能445は、設定された領域R211~領域R511に対応する点群を散布図上でそれぞれ色を変えて表示させる。操作者は、散布図上で色が変わっている点群をそれぞれ含むように、複数の関心領域を設定することで、領域R211~領域R511に対応する散布図上の関心領域R21~関心領域R51を設定することができる。
上述したように、操作者によって散布図上の関心領域が設定されると、解析機能444は、設定された関心領域内を、コントラストを強調させる対象組織の散布図におけるおおよその位置として特定する。ここで、解析機能444は、操作者によって設定された散布図上の関心領域の位置及び形状をメモリ41に保存することができる。操作者によって設定された散布図上での関心領域の位置は、コントラストを強調させる対象組織の基準物質の混合割合を示すものである。したがって、同一被検体だけでなく、異なる被検体間でも同じ対象組織における基準物質の混合割合はあまり変わらないため、この対象組織の散布図におけるおおよその位置が同様の位置となる可能性が高い。
そこで、解析機能444は、操作者によって設定された散布図上の関心領域の位置及び形状をメモリ41に保存しておき、以後に実施する関心領域の設定の処理において、メモリ41に保存した関心領域の位置を、コントラストを強調させる対象組織の散布図におけるおおよその位置として特定することで、操作者の関心領域の設定の手間を省くとともに、どのような操作者が解析を行ったとしても正確な解析画像を生成することを可能にする。例えば、解析機能444は、メモリ41に関心領域の位置及び形状の情報を保存した被検体と同一の被検体の解析画像を再度生成する場合に、メモリ41に保存した形状の関心領域を保存した位置に配置し、配置した関心領域内を対象組織の散布図におけるおおよその位置として特定する。
なお、解析機能444は、保存した関心領域の位置の情報を用いた場合に、関心領域の調整を受け付けることもできる。かかる場合には、コントラストを強調させる対象組織の散布図におけるおおよその位置を特定する処理が開始されると、表示制御機能445は、保存された形状の関心領域を保存された位置に対応する散布図上の位置に配置して表示させる。さらに、表示制御機能445は、配置した関心領域内に含まれる点群に対応する画像内の画素群の色を変えた単色X線画像を表示させる。
操作者は、表示された散布図上の関心領域と単色X線画像とを参照しながら、入力インターフェース43を操作して、散布図上の関心領域の形状や位置を調整する。解析機能444は、調整後の関心領域内を対象組織の散布図におけるおおよその位置として特定する。例えば、解析機能444は、メモリ41に関心領域の位置及び形状の情報を保存した被検体とは異なる被検体の解析画像を生成する場合に、上述した関心領域の調整を受け付ける。
なお、解析機能444は、上記した関心領域を、被検体ごと、或いは、操作者ごとに保存することもできる。かかる場合には、例えば、解析機能444は、被検体情報や、操作者情報を取得し、取得した情報と関心領域の情報とを対応付けてメモリ41に保存する。また、解析機能444は、設定された関心領域を対象組織ごとに保存することもできる。かかる場合には、解析機能444は、関心領域が設定された対象組織の情報を取得し、取得した対象組織の情報と設定された関心領域の情報とを対応付けてメモリ41に保存する。なお、対象組織の情報は、例えば、操作者によって入力される。
また、解析機能444は、メモリ41に保存した関心領域の情報を適宜更新することもできる。例えば、解析機能444は、保存した関心領域について、上記した関心領域の調整が行われるごとに、メモリ41に保存する関心領域を調整後の領域に更新してもよい。また、対象組織ごとに関心領域が設定されている場合、例えば、解析機能444は、対象組織の情報が取得される関心領域が新たに設定されるごとに、メモリ41に保存済みの関心領域の情報と新たに設定された関心領域の情報とを用いて、さらに新たな関心領域を設定する。そして、解析機能444は、メモリ41に保存された当該対象組織に対応する関心領域を新たに設定した関心領域に更新する。
上述したように、対象組織の散布図におけるおおよその位置として関心領域が設定されると、解析機能444は、設定された関心領域に含まれる点群に対して異なる重み付けを行い、重み付けをCT画像の輝度値に反映させることで、対象組織のコントラストを変化させる。例えば、解析機能444は、関心領域に含まれる点群に対して、重心に基づく重み付けを行うことで、対象組織内にコントラストを付ける。
図5は、第1の実施形態に係る解析機能444による重み付けの一例を説明するための図である。ここで、図5は、図4に示す散布図において設定された関心領域R51内を拡大した図である。すなわち、図5において示す点群は、関心領域R51に含まれる点群である。例えば、解析機能444は、図5に示すように、点群の重心の濃度が濃く、重心から離れるほど薄くなるように重み付けを行う。
上述したように、散布図は、基準物質の混合割合を示しており、同一の組織(同一の物質)の画素を示す点は散布図上の略同一の座標にプロットされるはずであるが、ノイズなどの影響によって、同一の組織(同一の物質)の画素でも点が分散して示されている。しかしながら、ノイズなどの影響によって点が分散されたとしても、同一の組織(同一の物質)の画素を示す点は、散布図上の略同一の座標付近に集中して分布していると考えられる。一例を挙げると、脳実質に含まれる所定の物質を対象組織として設定された関心領域では、所定の物質の画素を示す点が、散布図上の略同一の座標付近に集中して分布していると考えられる。
そこで、本願の解析機能444は、関心領域に含まれる点群の重心を抽出し、抽出した重心を中心とした重み付けを行う。例えば、解析機能444は、関心領域に含まれる点群を対象として、散布図の座標ごとの点の数を計数する。そして、解析機能444は、計数した座標ごとの点の数が最多となる座標を、関心領域に含まれる点群の重心として抽出する。さらに、解析機能444は、抽出した座標からの距離に応じた重み付けを設定する。例えば、解析機能444は、抽出した座標からの距離が長くなるほど重みが小さくなるように、重み付けを設定する。
そして、解析機能444は、設定した重み付けをCT画像の画素値に反映させることで、解析画像を生成する。例えば、解析機能444は、重心の座標の点に対応する画素が最も高輝度(低輝度)となり、重心からの距離が長くなるにしたがって低輝度(高輝度)となるように、重み付けを行うことで、関心領域に含まれる点群に対応する画素のうち、同一の組織(同一の物質)の画素をより強調して表示させることができる。
例えば、解析機能444は、脳実質に含まれる所定の物質(例えば、灰白質を構成する物質)を対象組織として設定された関心領域内の点群から重心を抽出し、抽出した重心の座標の点に対応する画素が最も高輝度となり、重心からの距離が長くなるにしたがって低輝度となるように、重み付けを行う。これにより、解析機能444は、例えば、脳実質における灰白質が強調され、白質とのコントラストが明瞭となるCT画像(解析画像)や、灰白質における正常組織と病変組織とのコントラストが明瞭となる解析画像などを生成することができる。すなわち、解析機能444は、組織において同一の組織(物質の組成が同一となる組織)のコントラストを選択的に強調させることで、当該組織と他の組織とを識別可能に示した解析画像を生成することができる。なお、重心の座標の点に対応する画素の輝度値および重心からの距離に応じた輝度の低下度合い(輝度の勾配)は、予め任意に設定することができる。
ここで、解析画像を生成する際のCT画像は、どのようなものが用いられてもよい。例えば、解析機能444は、所定のエネルギーにおける単色X線画像や、多色X線画像、基準物質画像などを用いて、解析画像を生成することができる。なお、解析画像に用いられるCT画像は、例えば、診断内容に基づいて決められてもよい。
また、解析機能444は、散布図上に設定された複数の関心領域について、それぞれ重心を抽出し、抽出した重心に基づく重み付けをそれぞれ行うこともできる。例えば、解析機能444は、図4に示す関心領域R11~関心領域R51について、それぞれ重心を抽出して、抽出した重心に基づく重み付けをそれぞれ行う。そして、解析機能444は、設定した重み付けを、CT画像において対応する画素の画素値にそれぞれ反映させることで、解析画像を生成する。
このとき、解析機能444は、複数の関心領域について、それぞれ異なる重み付けを行う。すなわち、解析機能444は、各関心領域において抽出した重心に対して、それぞれ異なる重み付けを行う。一例を挙げると、関心領域R11~関心領域R51についてそれぞれ重み付けを行う場合、解析機能444は、各関心領域の重心の座標に対応する画素の輝度値がそれぞれ異なるように、各関心領域において抽出した重心に対して、それぞれ異なる重み付けを行う。
図6は、第1の実施形態に係る解析画像の一例を示す図である。例えば、解析機能444は、各関心領域の重心の座標に対応する画素の輝度値がそれぞれ異なるように、各関心領域において抽出した重心に対して、それぞれ異なる重み付けを行い、CT画像上の対応する画素の画素値を重み付けに基づいて補正することで、図6に示すような組織分解能が向上された解析画像を生成する。すなわち、解析機能444は、上述した処理により、脂肪や、筋肉、腫瘍、水、骨などについて、コントラストが強調されたCT画像を生成することができる。これは、造影剤を用いない単純CT検査によって撮影したCT画像をもとに生成することができる。換言すると、本実施形態に係るX線CT装置1は、単純CT検査の撮影によって、組織分解能が高いMRI画像に類似した解析画像を生成することができる。
なお、散布図における関心領域が複数の場合、各関心領域の重心の座標に対応する画素の輝度値をそれぞれ異なるようにすることが困難となる場合がある。そこで、X線CT装置1では、各関心領域における重み付けを、操作者が自由に変更できるGUIを設けることができる。かかる場合には、例えば、表示制御機能445が、解析機能444によって生成された解析画像とともに、関心領域R11~関心領域R51の重み付けをそれぞれ独立して変更することができるGUIをディスプレイ42に表示させる。操作者は、入力インターフェース43を介して、ディスプレイ42上のGUIを操作することで、関心領域R11~関心領域R51における重み付けをそれぞれ変更させる。ここで、操作者による重み付けの変更は、重心の座標の画素に適用する輝度値の変更と、重心からの距離に応じて変化する輝度の勾配とを含む。
解析機能444は、入力インターフェース43を介した各関心領域における重み付けの変更操作に応じて輝度値を変更した解析画像を生成する。表示制御機能445は、解析機能444によって解析画像が生成されるごとに、ディスプレイ42に表示させる解析画像を随時更新する。例えば、操作者は、ディスプレイ42に表示された解析画像を観察しながら重み付けを変更するためのGUIを操作することで、所望の解析画像を表示させることができる。これにより、X線CT装置1は、脳実質と水などのように画像上で重なる領域であっても、それぞれを区別して解析画像を表示させることができる。
なお、解析画像の表示形態は、図6に示す例に限られず、その他、種々の形態で表示されてもよい。例えば、表示制御機能445は、図3や図4に示すように、散布図や、基準物質画像と並べて解析画像を表示させることができる。ここで、散布図を並べて表示させる場合、表示制御機能445は、散布図に設定された関心領域を変更する操作を受け付けるGUIをさらに表示させることもできる。すなわち、操作者は、解析画像を参照しながら、散布図における関心領域や、関心領域ごとの重み付けを任意に変更することができる。
また、図6に示すように、複数の関心領域に対してそれぞれ設定した重み付けを、単一の解析画像に対して反映させるだけではなく、例えば、関心領域ごとに解析画像をそれぞれ生成して、生成した複数の解析画像を並べて表示させる場合でもよい。一例を挙げると、解析機能444は、関心領域R11に含まれる点群に対応する画素の画素値に重み付けを反映させた解析画像と、関心領域R21に含まれる点群に対応する画素の画素値に重み付けを反映させた解析画像と、関心領域R31に含まれる点群に対応する画素の画素値に重み付けを反映させた解析画像と、関心領域R41に含まれる点群に対応する画素の画素値に重み付けを反映させた解析画像と、関心領域R51に含まれる点群に対応する画素の画素値に重み付けを反映させた解析画像とを生成する。表示制御機能445は、生成された5つの解析画像を並べて表示させる。なお、表示制御機能445は、これらの解析画像についても、関心領域の重み付けをそれぞれ変更することができるGUIをディスプレイ42に表示させることができる。これにより、X線CT装置1は、例えば、水や、脂肪、筋肉、骨、などの組織(物質)ごとに強調された解析画像を生成して表示させることができる。
次に、X線CT装置1による処理の手順の一例を、図7を用いて説明する。図7は、第1の実施形態に係るX線CT装置1の処理の流れを説明するためのフローチャートである。なお、図7においては、操作者によって関心領域(ROI:Region of Interest)が設定される場合について示す。
ステップS101は、スキャン制御機能441に対応するステップである。ステップS102は、画像再構成機能443に対応するステップである。ステップS103は、画像再構成機能443及び解析機能444に対応するステップである。ステップS104~ステップS106及びステップS108は、解析機能444に対応するステップである。ステップS107は、表示制御機能445に対応するステップである。
まず、処理回路44は、異なるエネルギーで複数のデータセットを取得する(ステップS101)。次に、処理回路44は、基準物質画像を生成して(ステップS102)、散布図及び表示画像を生成して表示する(ステップS103)。その後、処理回路44は、ROIを受け付けたか否かを判定する(ステップS104)。
ここで、ROIの設定を受け付けると(ステップS104肯定)、処理回路44は、ROIに該当する画素において、各座標における画素数をカウントして、重心を設定する(ステップS105)。なお、処理回路44は、ROIの設定を受け付けるまで、待機状態である(ステップS104否定)。
ステップS105において重心を設定すると、処理回路44は、重心からの距離に応じて重み付けを変化させ(ステップS106)、変化させた重み付けに応じた濃淡をつけた表示画像(解析画像)を表示させる(ステップS107)。その後、処理回路44は、終了か否かを判定して(ステップS108)、終了操作を受け付けると(ステップS108肯定)、処理を終了する。なお、終了操作を受け付けるまで、処理回路44は、待機状態となる(ステップS108否定)。
上述したように、第1の実施形態によれば、スキャン制御機能441は、X線を利用したスキャンを実行することで、被検体の一部を含む領域について、第1のX線エネルギーに対応する第1データセット、及び第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを収集する。解析機能444は、第1データセット及び第2データセットに基づいて、領域内の各位置における複数の基準物質それぞれの含有量を表す散布図を生成し、散布図の少なくとも一部について重みを設定し、含有量及び重みに基づいて解析画像を生成する。従って、第1の実施形態に係るX線CT装置1は、同一の組織(物質)に対して重み付けして画像化した解析画像を表示させることができ、同一の組織(物質)を強調させ、組織分解能を向上させたCT画像を生成することを可能にする。
また、上述したように、第1の実施形態によれば、解析機能444は、散布図に設定された複数の関心領域について、異なる重みを設定する。従って、第1の実施形態に係るX線CT装置1は、組織ごとにコントラストをつけることができ、より組織分解能を向上させたCT画像を生成することを可能にする。
また、上述したように、第1の実施形態によれば、入力インターフェース43は、関心領域を変更する変更操作を受け付ける。解析機能444は、入力インターフェース43によって受け付けられた変更操作後の関心領域について、重みを設定する。従って、第1の実施形態に係るX線CT装置1は、対象組織に含める画素を適宜調整することができ、より観察しやすいCT画像を生成することを可能にする。
また、上述したように、第1の実施形態によれば、メモリ41は、散布図に設定された関心領域の情報を記憶する。解析機能444は、メモリ41に記憶された関心領域の情報に基づいて、新たに生成された散布図に関心領域を設定する。従って、第1の実施形態に係るX線CT装置1は、関心領域を設定する手間を省くとともに、どのような操作者によって操作された場合でも観察しやすいCT画像を生成することを可能にする。
また、上述したように、第1の実施形態によれば、解析機能444は、散布図において示される含有量の分布において、散布図に設定された関心領域における重心を抽出し、抽出した重心からの距離に応じた重みを設定する。従って、第1の実施形態に係るX線CT装置1は、同一の組織(物質)を容易に強調させることを可能にする。
また、上述したように、第1の実施形態によれば、解析機能444は、重心からの距離が長くなるほど重みが小さくなるように、重みを設定する。従って、第1の実施形態に係るX線CT装置1は、同一の組織(物質)を、精度よく強調させることを可能にする。
また、上述したように、第1の実施形態によれば、入力インターフェース43は、重みを変更する変更操作を受け付ける。解析機能444は、入力インターフェース43によって受け付けられた変更操作後の重みと含有量とに基づいて、解析画像を生成する。従って、第1の実施形態に係るX線CT装置1は、組織(物質)ごとにコントラストを調整することができ、より観察しやすいCT画像を提供することを可能にする。
また、上述したように、第1の実施形態によれば、解析機能444は、複数の基準物質それぞれの含有量を、領域内の各位置について推定し、散布図の少なくとも一部に対応する、領域の少なくとも一部について推定した含有量に基づく画素値を、重みで補正することで、解析画像を生成する。従って、第1の実施形態に係るX線CT装置1は、組織(物質)ごとにコントラストを調整することを可能にする。
(第2の実施形態)
上述した第1の実施形態では、関心領域ごとに重み付けをつける場合について説明した。第2の実施形態では、関心領域の位置ごとに重み付けをつける場合について説明する。第2の実施形態に係るX線CT装置1は、図2に示したX線CT装置1と同様の構成を有し、解析機能444による処理の一部が相違する。以下、第1の実施形態において説明した構成と同様の構成を有する点については、図2と同一の符号を付し、説明を省略する。
第2の実施形態に係る解析機能444は、散布図に設定された関心領域内の複数の位置について、異なる重みを設定する。例えば、解析機能444は、含有量の分布において、重心からの距離が閾値よりも小さい範囲に対して、同一の重みを設定する。図8は、第2の実施形態に係る重み付けの一例を説明するための図である。ここで、図8においては、図5に示す散布図(関心領域R51の拡大図)において位置ごとに重み付けを付ける場合について示す。また、図8において点群の分布の横に示すグラフは、縦軸が点の数を示し、横軸が点の位置を示すヒストグラムを示す。
例えば、解析機能444は、散布図における点群の分布に応じて、散布図に範囲を設定し、設定した範囲ごとに重み付けを行う。一例を挙げると、解析機能444は、図8に示すように、散布図に対して、ヒストグラムにおいてピーク位置を含む範囲aと、範囲bと、範囲cとを設定し、設定した範囲ごとに異なる重みを設定する。
例えば、解析機能444は、範囲aに含まれる点群に対して同一の重み付けを行い、範囲bと範囲cに含まれる点群に対しては、範囲a側の点群から重みが小さくなる重み付けをそれぞれ行う。これにより、解析機能444によって生成される解析画像は、範囲aに含まれる点群に対応する画素を強調したものとなる。上述したように、散布図においては、ノイズの影響などにより、同一の組織(物質)であっても散布図上の位置がずれる場合がある。上述した解析機能444の重み付けにより、ノイズの影響などによって点の位置がずれた同一の組織(物質)に画素を同一の輝度値で表現することができ、ノイズの低減された観察画像を提供することができる。
また、解析機能444は、範囲aに含まれる点群と、範囲bに含まれる点群と、範囲cに含まれる点群からそれぞれ重心を抽出して、抽出した重心に基づく重み付けをそれぞれ行うこともできる。一例を挙げると、解析機能444は、範囲aに含まれる点群を位置ごとにカウントして、重心の位置を抽出する。そして、解析機能444は、範囲aに含まれる点群を対象として、抽出した重心の濃度が濃く、重心から離れるほど薄くなるように重み付けを行う。同様に、解析機能444は、範囲bに含まれる点群を位置ごとにカウントして、重心の位置を抽出する。そして、解析機能444は、範囲bに含まれる点群を対象として、抽出した重心の濃度が濃く、重心から離れるほど薄くなるように重み付けを行う。また、解析機能444は、範囲cに含まれる点群を位置ごとにカウントして、重心の位置を抽出する。そして、解析機能444は、範囲cに含まれる点群を対象として、抽出した重心の濃度が濃く、重心から離れるほど薄くなるように重み付けを行う。
このように、解析機能444は、関心領域内の範囲ごとに異なる重み付けを行うことで、組織(物質)に対してより細かいコントラストを付けることができる。例えば、脳実質における浮腫に含まれる所定の物質を対象組織とした関心領域において、上記した重み付けを行うことで、解析機能444は、図9に示すように、矢印a6で指示する完全な浮腫領域と、矢印a5で指示する脳実質が残っている浮腫領域とを、それぞれ明確に示しつつ、異なる輝度値で示した解析画像を生成することができる。すなわち、同じ物質を含む組織であっても、含まれる物質の組成が異なる組織の場合、関心領域内で異なる位置に点群が形成される。本実施形態に係る解析機能444は、このような組織をそれぞれ識別可能に強調させることができる。なお、図9は、第2の実施形態に係る解析画像の一例を示す図である。
上述したように、第2の実施形態によれば、解析機能444は、散布図に設定された関心領域内の複数の位置について、異なる重みを設定する。従って、第2の実施形態に係るX線CT装置1は、組織(物質)に対してより細かいコントラストを付けることを可能にする。
また、第2の実施形態によれば、解析機能444は、含有量の分布において、重心からの距離が閾値よりも小さい範囲に対して、同一の重みを設定する。従って、第2の実施形態に係るX線CT装置1は、ノイズの低減されたCT画像を提供することを可能にする。
(第3の実施形態)
第3の実施形態では、散布図(画像)におけるノイズを特定し、特定したノイズを低減する場合について説明する。第3の実施形態に係るX線CT装置1は、図2に示したX線CT装置1と同様の構成を有し、解析機能444による処理の一部が相違する。以下、第1の実施形態において説明した構成と同様の構成を有する点については、図2と同一の符号を付し、説明を省略する。
第3の実施形態に係る解析機能444は、散布図におけるノイズを特定し、特定したノイズを低減した解析画像を生成する。例えば、解析機能444は、散布図において出現頻度が閾値よりも小さい点を、ノイズとして特定する。一例を挙げると、解析機能444は、散布図における点群において、同一の座標に分布する点の数を計数し、計数した値が閾値よりも少ない座標の点をノイズとして特定する。そして、解析機能444は、ノイズとして特定した点に対応する画素の輝度値を、例えば、重心の位置の点に対応する画素の輝度値と同じ値に補正することで、解析画像におけるノイズを低減する。
なお、上述した例では、ノイズとして特定した点に対応する画素の輝度値を、重心の位置の点に対応する画素の輝度値と同じ値に補正することでノイズを低減する例について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、ノイズとして特定した点を解析の対象外とし、重み付けを付与する画素から除く場合でもよい。
また、上述した例では、散布図を用いてノイズを特定する例について説明した。しかしながら、実施形態はこれに限定されるものではなく、画像からノイズを特定する場合でもよい。かかる場合には、解析機能444は、例えば、画像においてノイズに対応する画素を特定する。そして、解析機能444は、ノイズとして特定した画素に対して、その他のノイズ低減処理を実施する場合でもよい。
上述したように、第3の実施形態によれば、解析機能444は、散布図におけるノイズを特定し、特定したノイズを低減した解析画像を生成する。従って、第3の実施形態に係るX線CT装置1は、ノイズを低減させた解析画像を生成することを可能にする。
また、第3の実施形態によれば、解析機能444は、散布図において出現頻度が閾値よりも小さい点を、ノイズとして特定する。従って、第3の実施形態に係るX線CT装置1は、散布図から容易にノイズを特定することを可能にする。
(第4の実施形態)
上述した実施形態では、重み付けに基づいて、解析画像における輝度値を変化させる場合について説明した。第4の実施形態では、重み付けに基づいて、カラー画像を生成する場合について説明する。第4の実施形態に係るX線CT装置1は、図2に示したX線CT装置1と同様の構成を有し、解析機能444による処理の一部が相違する。以下、第1の実施形態において説明した構成と同様の構成を有する点については、図2と同一の符号を付し、説明を省略する。
第4の実施形態に係る解析機能444は、重み付けに応じた色を対応する画素に付与することで、カラー化させた解析画像を生成する。例えば、解析機能444は、関心領域ごとに色を対応付け、重み付けに基づいてカラーの濃度を変化させた解析画像を生成する。図10は、第4の実施形態に係る解析画像の一例を示す図である。なお、図10においては、図6に示す解析画像をカラー化させた場合について示す。
例えば、解析機能444は、図10の領域R6(脳脊髄液の領域)の画素に対応する関心領域に対して、所定の色(例えば、青等)を対応付けるとともに、関心領域内の点群に基づいて設定した重み付けに応じた色の濃淡をつけたカラー画像を生成する。すなわち、解析機能444は、領域R6に含まれる画素を、重み付けに応じて濃淡を変化させた青で示したカラー画像を生成する。表示制御機能445は、解析機能444によって生成されたカラー画像をディスプレイ42に表示させる。
ここで、解析機能444によるカラー画像の生成は、関心領域の選択に応じて実施させることができる。すなわち、操作者が入力インターフェース43を介して関心領域を変更するごとに、解析機能444は、変更後の関心領域を対象としたカラー画像を生成する。そして、表示制御機能445が、生成されたカラー画像をディスプレイ42に表示させる。これにより、操作者は、所望の組織(物質)に対して色を付けたカラー画像を観察することができる。
(第5の実施形態)
これまで第1~第4の実施形態について説明したが、上述した第1~第4の実施形態以外にも、種々の異なる形態にて実施されてよいものである。
上述した実施形態では、Dual-Energyによる撮影によって収集した高エネルギー投影データセットと低エネルギー投影データセットとを用いる場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、3種類以上の異なる管電圧での「Multi Energyによる撮影」によって収集した3種類以上の異なるエネルギーに対応する投影データセットが用いられる場合でもよい。
また、上述した実施形態では、組織ごとのコントラストを強調させることで、組織分解能が高いMRI画像に類似した解析画像を生成する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、仮想的な造影画像を生成することもできる。かかる場合には、例えば、解析機能444は、腫瘍や血管の富む組織等の造影剤の影響を受けやすい部分に対応する関心領域に含まれる点群に対応する画素群が最も高い輝度値を示すように重み付けすることにより、単純CT検査において撮影したCT画像を用いて、仮想的な造影画像を生成することができる。
また、上述した実施形態では、脳の解析画像を生成する例について説明したが、実施形態はこれに限定されるものではなく、その他、種々の部位の解析画像を生成する場合でもよい。
また、上述した実施形態では、点群における重心を重み付けの中心とする例について説明した。しかしながら、実施形態はこれに限定されるものではなく、点群の分布におけるピークの位置を重み付けの中心にしてもよい。
また、上述した実施形態では、単純CT検査の撮影で収集したデータセットを用いる例について説明した。しかしながら、実施形態はこれに限定されるものではなく、造影CT検査の撮影で収集したデータセットを用いてもよい。
また、上述した実施形態では、各組織(物質)のコントラストを強調するように、重み付けを行う例について説明した。しかしながら、実施形態はこれに限定されるものではなく、所定の組織(物質)に対応する画素が表示されないように重み付けを行う場合でもよい。一例を挙げると、解析機能444は、骨に対応する関心領域に含まれる点群に対して、輝度値が「0」となるように重み付けを行うことで、骨が非表示となる解析画像を生成する。
また、上述した実施形態では、複数の基準物質それぞれの含有量を表す散布図を生成する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、散布図を生成せずに基準物質の含有量に基づく基準物質の混合割合を用いて同一の組織(同一の物質)を抽出する場合でもよい。かかる場合には、解析機能444は、第1データセット及び第2データセットに基づいて、領域内の各位置における複数の基準物質それぞれの含有量を推定し、領域内の少なくとも一部について重みを設定し、含有量及び重みに基づいて解析画像を生成する。
具体的には、解析機能444は、異なるエネルギーに対応する複数のデータセットに基づいて、画素ごとに、複数の基準物質間での存在比率を推定する。そして、解析機能444は、散布図を生成せずに、各画素における存在比率に基づいて、解析画像を生成する。すなわち、解析機能444は、数値のみで処理を実行して、解析画像を生成する。例えば、解析機能444は、Dual Energyによって撮影された領域の各画素について、水に対するヨードの存在比率と、ヨードに対する水の存在比率とをそれぞれ推定する。そして、解析機能444は、推定した存在比率に基づいて、画素をグループに分類し、グループ間や、グループ内で重みを設定することで、各画素の輝度値を、重みが反映された輝度値に変換した解析画像を生成する。
図11は、第5の実施形態に係る解析機能444による処理を説明するための模式図である。なお、図11の上段は、画像データにおける各画素を示す。すなわち。(1,1)、(1,2)、(1,3)等は、それぞれ各画素の画像データにおける座標を示す。なお、図11においては、各画素を座標によって識別しているが、単に数字を用いて識別する場合でもよい。
例えば、解析機能444は、図11の中段の図に示すように、画像再構成機能443によって生成された高エネルギーに対応する多色X線画像の画素値と低エネルギーに対応する多色X線画像の画素値との関係に基づいて、画素ごとに、水に対するヨードの存在比率「Iodine(Water)」と、ヨードに対する水の存在比率「Water(Iodine)」を推定する。一例を挙げると、解析機能444は、画素(1,1)について、「Iodine(Water):10、Water(Iodine):815」を推定する。
上述したように全ての画素について水に対するヨードの存在比率とヨードに対する水の存在比率を推定すると、解析機能444は、推定した存在比率に基づいて、画素をグループに分類する。例えば、解析機能444は、存在比率の類似度に基づいて、画素をグループに分類する。なお、類似度の判定に用いる手法は、既存の手法が任意に用いられる。
一例を挙げると、解析機能444は、全画素のうち、「Iodine(Water):10、Water(Iodine):815」、「Iodine(Water):13、Water(Iodine):822」、「Iodine(Water):11、Water(Iodine):811」、・・・「Iodine(Water):13、Water(Iodine):772」、「Iodine(Water):13、Water(Iodine):822」・・に対応する各画素をグループAに分類する。また、解析機能444は、「Iodine(Water):60、Water(Iodine):1022」の画素及びこれに類似した各画素をグループBに分類し、「Iodine(Water):255、Water(Iodine):1506」の画素及びこれに類似した各画素をグループCに分類する。このように、解析機能444は、基準物質の含有量に基づいて、各画素をグループに分類する。
ここで、各グループは、水や、脂肪、筋肉、骨、などの組織(物質)にそれぞれ対応する。すなわち、識別の対象となる組織(物質)が何なるかに応じて、グループが変化することとなる。換言すると、識別の対象となる組織(物質)ごとに、基準物質の存在比率の範囲が設定される(例えば、脂肪の場合、Iodine(Water)がa~bまで、Water(Iodine)がc~dまで等)。そして、解析画像の生成の際にどの組織(物質)を対象とするかによって、グループ分けが変化する。
このように、各画素をグループに分類すると、解析機能444は、グループ間や、グループ内で重みを設定することで、各画素の輝度値を、重みが反映された輝度値に変換した解析画像を生成する。例えば、解析機能444は、グループAと、グループBと、グループCに対して異なる重み付けを行うことで、脂肪や、筋肉、腫瘍、水、骨などの組織(物質)について、コントラストを強調させた解析画像を生成することができる。
また、解析機能444は、グループ内で、画素ごとに異なる重み付けを行うことで、ノイズを低減したり、同一組織内で細かくコントラストをつけたりした解析画像を生成することができる。例えば、解析機能444は、存在比率の出現頻度に応じて異なる重みを設定する。一例を挙げると、解析機能444は、グループAおいて、出現頻度が高い「Iodine(Water):13、Water(Iodine):822」の画素と、出現頻度が低いその他の画素に対して異なる重み付けを行う。これにより、グループAに対応する組織において、ノイズを低減したり、細かくコントラストをつけたりすることが可能となる。
また、上述した実施形態では、X線CT装置1が各種処理を実行する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、医用画像処理装置3において各種処理が実行される場合であってもよい。図12は、第5の実施形態に係る医用画像処理装置3の構成の一例を示す図である。
例えば、図12に示すように、医用画像処理装置3は、入力インターフェース310と、ディスプレイ320と、記憶回路330と、処理回路340とを有する。
入力インターフェース310は、種々の設定などを行うためのトラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。
入力インターフェース310は、処理回路340に接続されており、操作者から受け付けた入力操作を電気信号に変換して処理回路340に出力する。なお、本明細書において入力インターフェース310は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を処理回路340へ出力する処理回路も入力インターフェースの例に含まれる。
ディスプレイ320は、処理回路340に接続され、処理回路340から出力される各種情報及び各種画像データを表示する。例えば、ディスプレイ320は、液晶モニタやCRT(Cathode Ray Tube)モニタ、タッチパネル等によって実現される。例えば、ディスプレイ320は、操作者の指示を受け付けるためのGUI(Graphical User Interface)や、種々の表示画像、処理回路350による種々の処理結果(例えば、散布図や、解析画像など)を表示する。
記憶回路330は、処理回路340に接続され、各種データを記憶する。例えば、記憶回路330は、RAM(Random Access Memory)、フラッシュメモリ等の半導体メモリ素子や、ハードディスク、光ディスク等によって実現される。例えば、記憶回路330は、X線CT装置1又は画像保管装置2から受信したCT画像データや、散布図、解析画像、関心領域の情報などを記憶する。また、記憶回路330は、処理回路340によって実行される各処理機能に対応するプログラムを記憶する。
処理回路340は、入力インターフェース310を介して操作者から受け付けた入力操作に応じて、医用画像処理装置3が有する各構成要素を制御する。例えば、処理回路340は、プロセッサによって実現される。本実施形態では、処理回路340は、入力インターフェース310から出力されるCT画像データを記憶回路330に記憶させる。また、処理回路340は、記憶回路330からCT画像データを読み出し、ディスプレイ320に表示する。
そして、処理回路350は、図12に示すように、制御機能341と、画像再構成機能342と、解析機能343と、表示制御機能344とを実行する。制御機能341は、医用画像処理装置3の全体を制御する。なお、制御機能341は、取得部の一例である。画像再構成機能342は、上述した画像再構成機能443と同様の処理を実行する。解析機能343は、上述した解析機能444と同様の処理を実行する。表示制御機能344は、上述した表示制御機能445と同様の処理を実行する。
上述した実施形態では、X線CT装置1及び医用画像処理装置3による独立した処理について説明した。しかしながら、実施形態はこれに限定されるものではなく、X線CT装置1と医用画像処理装置3とが協働するX線CTシステムとして機能する場合であってもよい。かかる場合には、X線CT装置1によって収集されたCT画像データが、医用画像処理装置3に転送され、医用画像処理装置3において、種々の処理が実行される。
上記説明において用いた「プロセッサ」という文言は、例えば、CPU、GPU(Graphics Processing Unit)、あるいは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサは、メモリ41又はメモリ53に保存されたプログラムを読み出し実行することで機能を実現する。
なお、図2においては、単一のメモリ41が各処理機能に対応するプログラムを記憶するものとして説明した。また、図12においては、単一の記憶回路330が各処理機能に対応するプログラムを記憶するものとして説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、複数のメモリ41を分散して配置するとともに、処理回路44が個別のメモリ41から対応するプログラムを読み出す構成としても構わない。また、例えば、複数の記憶回路330を分散して配置するとともに、処理回路340が個別の記憶回路330から対応するプログラムを読み出す構成としても構わない。また、メモリ41又は記憶回路330にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。
また、処理回路44及び処理回路340は、ネットワークを介して接続された外部装置のプロセッサを利用して、機能を実現することとしてもよい。例えば、処理回路44は、メモリ41から各機能に対応するプログラムを読み出して実行するとともに、X線CT装置1とネットワークNWを介して接続された外部のワークステーションやクラウドを計算資源として利用することにより、図2に示す各機能を実現する。また、例えば、処理回路340は、記憶回路330から各機能に対応するプログラムを読み出して実行するとともに、医用画像処理装置3とネットワークNWを介して接続された外部のワークステーションやクラウドを計算資源として利用することにより、図12に示す各機能を実現する。
上述した実施形態に係る各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。即ち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行われる各処理機能は、その全部又は任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現されうる。
また、上述した実施形態で説明した制御方法は、予め用意された制御プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この制御プログラムは、インターネット等のネットワークを介して配布することができる。また、この制御プログラムは、ハードディスク、フレキシブルディスク(FD)、CD-ROM、MO、DVD等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
以上説明した少なくとも1つの実施形態によれば、組織分解能を向上させたCT画像を提供することができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
1 X線CT装置
3 医用画像処理装置
44、340 処理回路
441 スキャン制御機能
442 前処理機能
342、443 画像再構成機能
343、444 解析機能
344、445 表示制御機能

Claims (15)

  1. X線を利用したスキャンを実行することで、被検体の一部を含む領域について、第1のX線エネルギーに対応する第1データセット、及び前記第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを収集するスキャン部と、
    前記第1データセット及び前記第2データセットに基づいて、前記領域内の各位置における複数の基準物質それぞれの含有量を表す散布図を生成し、前記散布図の少なくとも一部に含まれる前記含有量に対して重みを設定し、前記含有量及び前記重みを対応する画素の画素値に反映させた解析画像を生成する解析部と、
    を備えたX線CT装置。
  2. 前記解析部は、前記散布図に設定された複数の関心領域について、異なる重みを設定する、請求項1に記載のX線CT装置。
  3. 前記解析部は、前記散布図に設定された関心領域内の複数の位置について、異なる重みを設定する、請求項1に記載のX線CT装置。
  4. 前記関心領域を変更する変更操作を受け付ける受付部をさらに備え、
    前記解析部は、前記受付部によって受け付けられた変更操作後の関心領域について、重みを設定する、請求項2又は3に記載のX線CT装置。
  5. 前記散布図に設定された関心領域の情報を記憶する記憶部をさらに備え、
    前記解析部は、前記記憶部に記憶された関心領域の情報に基づいて、新たに生成された散布図に関心領域を設定する、請求項2~4のいずれか1つに記載のX線CT装置。
  6. 前記解析部は、前記散布図において示される前記含有量の分布において、前記散布図に設定された関心領域における重心を抽出し、抽出した重心からの距離に応じた重みを設定する、請求項2~5のいずれか1つに記載のX線CT装置。
  7. 前記解析部は、前記重心からの距離が長くなるほど重みが小さくなるように、前記重みを設定する、請求項6に記載のX線CT装置。
  8. 前記解析部は、前記含有量の分布において、前記重心からの距離が閾値よりも小さい範囲に対して、同一の重みを設定する、請求項6に記載のX線CT装置。
  9. 前記重みを変更する変更操作を受け付ける受付部をさらに備え、
    前記解析部は、前記受付部によって受け付けられた変更操作後の重みと前記含有量とに基づいて、前記解析画像を生成する、請求項1~8のいずれか1つに記載のX線CT装置。
  10. 前記解析部は、
    前記複数の基準物質それぞれの含有量を、前記領域内の各位置について推定し、
    前記散布図の少なくとも一部に対応する、前記領域の少なくとも一部について推定した前記含有量に基づく画素値を、前記重みで補正することで、前記解析画像を生成する、請求項1~9のいずれか1つに記載のX線CT装置。
  11. 前記解析部は、前記散布図におけるノイズを特定し、特定したノイズを低減した解析画像を生成する、請求項1~10のいずれか1つに記載のX線CT装置。
  12. 前記解析部は、前記散布図において出現頻度が閾値よりも小さい前記含有量を、前記ノイズとして特定する、請求項11に記載のX線CT装置。
  13. X線を利用したスキャンを実行することで、被検体の一部を含む領域について収集された、第1のX線エネルギーに対応する第1データセット、及び前記第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを取得する取得部と、
    前記第1データセット及び前記第2データセットに基づいて、前記領域内の各位置における複数の基準物質それぞれの含有量を表す散布図を生成し、前記散布図の少なくとも一部に含まれる前記含有量に対して重みを設定し、前記含有量及び前記重みを対応する画素の画素値に反映させた解析画像を生成する解析部と、
    を備えた医用画像処理装置。
  14. X線を利用したスキャンを実行することで、被検体の一部を含む領域について、第1のX線エネルギーに対応する第1データセット、及び前記第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを収集するスキャン部、
    を有するX線CT装置と、
    前記第1データセット及び前記第2データセットに基づいて、前記領域内の各位置における複数の基準物質それぞれの含有量を表す散布図を生成し、前記散布図の少なくとも一部に含まれる前記含有量に対して重みを設定し、前記含有量及び前記重みを対応する画素の画素値に反映させた解析画像を生成する解析部、
    を有する医用画像処理装置と
    を備えたX線CTシステム。
  15. X線を利用したスキャンを実行することで、被検体の一部を含む領域について、第1のX線エネルギーに対応する第1データセット、及び前記第1のX線エネルギーと異なる第2のX線エネルギーに対応する第2データセットを収集するスキャン部と、
    前記第1データセット及び前記第2データセットに基づいて、前記領域内の各位置における複数の基準物質それぞれの含有量を推定し、前記領域内の少なくとも一部に含まれる前記含有量に対して重みを設定し、前記含有量及び前記重みを対応する画素の画素値に反映させた解析画像を生成する解析部、
    を備えたX線CT装置。
JP2019162621A 2018-09-07 2019-09-06 X線ct装置、医用画像処理装置及びx線ctシステム Active JP7341803B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US16/563,569 US11295487B2 (en) 2018-09-07 2019-09-06 X-ray CT apparatus, medical image processing apparatus, and X-ray CT system
CN201910864900.3A CN110881993B (zh) 2018-09-07 2019-09-09 X射线ct装置、医用图像处理装置以及x射线ct***

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018167581 2018-09-07
JP2018167581 2018-09-07

Publications (2)

Publication Number Publication Date
JP2020039872A JP2020039872A (ja) 2020-03-19
JP7341803B2 true JP7341803B2 (ja) 2023-09-11

Family

ID=67875372

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019162621A Active JP7341803B2 (ja) 2018-09-07 2019-09-06 X線ct装置、医用画像処理装置及びx線ctシステム

Country Status (2)

Country Link
EP (1) EP3620112B1 (ja)
JP (1) JP7341803B2 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022032456A1 (zh) * 2020-08-10 2022-02-17 西安大医集团股份有限公司 图像重建方法、装置及计算机存储介质
CN112884858B (zh) * 2021-02-05 2024-04-26 明峰医疗***股份有限公司 Pet探测器的位置信息的分区方法及计算机可读存储介质
JP2024032518A (ja) * 2022-08-29 2024-03-12 キヤノン株式会社 画像処理装置、撮影システム、画像処理方法、及びプログラム
CN116646061B (zh) * 2023-04-28 2024-01-26 西安交通大学 分布式ct成像和智能诊疗***及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011110245A (ja) 2009-11-27 2011-06-09 Ge Medical Systems Global Technology Co Llc 画像表示装置、x線ct装置およびプログラム
WO2016158234A1 (ja) 2015-03-30 2016-10-06 株式会社日立製作所 画像生成装置、画像生成方法及びx線ct装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011110245A (ja) 2009-11-27 2011-06-09 Ge Medical Systems Global Technology Co Llc 画像表示装置、x線ct装置およびプログラム
WO2016158234A1 (ja) 2015-03-30 2016-10-06 株式会社日立製作所 画像生成装置、画像生成方法及びx線ct装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PETRONGOLO Michael et al.,Noise Suppression for Dual-Energy CT Through Entropy Minimization,IEEE Transactions on Medical Imaging,2015年11月,Vol.34, No.11,pp.2286-2297,DOI: 10.1109/TMI.2015.2429000

Also Published As

Publication number Publication date
JP2020039872A (ja) 2020-03-19
EP3620112B1 (en) 2023-04-19
EP3620112A1 (en) 2020-03-11

Similar Documents

Publication Publication Date Title
JP7341803B2 (ja) X線ct装置、医用画像処理装置及びx線ctシステム
US10052077B2 (en) Tomography imaging apparatus and method
JP7348376B2 (ja) 医用処理装置、x線ctシステム及び処理プログラム
JP7242284B2 (ja) 医用画像処理装置、x線診断装置及び医用画像処理プログラム
JP7467253B2 (ja) X線ctシステム及び医用処理装置
US10497152B2 (en) X-ray CT apparatus and reconstruction processing apparatus
CN110881993B (zh) X射线ct装置、医用图像处理装置以及x射线ct***
US20210196221A1 (en) Medical image processing apparatus, recording medium, and system
JP7466401B2 (ja) 医用画像診断装置
JP7483361B2 (ja) 医用画像処理装置、医用画像診断装置および医用画像処理プログラム
JP7353882B2 (ja) X線ctシステム及び医用処理装置
JP6373558B2 (ja) X線ct装置
JP2020103571A (ja) 医用処理装置及びx線診断システム
JP7467222B2 (ja) 医用情報処理装置、医用情報処理方法及び医用情報処理プログラム
JP2020096693A (ja) X線ctシステム及び処理プログラム
JP7321846B2 (ja) X線ctシステム及び医用処理装置
US20220383501A1 (en) Medical information display apparatus, medical image processing apparatus, medical information display method, and computer program product
JP6923414B2 (ja) X線コンピュータ断層撮影装置、医用画像処理装置、及びプログラム
US20240180503A1 (en) Photon counting computed tomography apparatus and photon-counting ct-scanning condition setting method
JP7233911B2 (ja) X線ctシステム及び処理プログラム
JP7236239B2 (ja) X線ctシステム及びctスキャン実行プログラム
JP2022065390A (ja) 医用画像処理装置、医用画像診断装置及びプログラム
JP2023012246A (ja) 医用画像処理装置、x線ct装置、医用画像処理方法及びプログラム
JP2024066111A (ja) 医用画像処理装置及び方法
JP2023160048A (ja) 医用画像処理装置、医用画像処理システム及び医用画像処理方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20190906

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220711

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230317

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230404

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230529

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230830

R150 Certificate of patent or registration of utility model

Ref document number: 7341803

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150