JP7334048B2 - X線ctシステム及び方法 - Google Patents

X線ctシステム及び方法 Download PDF

Info

Publication number
JP7334048B2
JP7334048B2 JP2019053607A JP2019053607A JP7334048B2 JP 7334048 B2 JP7334048 B2 JP 7334048B2 JP 2019053607 A JP2019053607 A JP 2019053607A JP 2019053607 A JP2019053607 A JP 2019053607A JP 7334048 B2 JP7334048 B2 JP 7334048B2
Authority
JP
Japan
Prior art keywords
projection data
low
projection
data
sinogram
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
JP2019053607A
Other languages
English (en)
Other versions
JP2020099662A (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
Original Assignee
Canon 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 Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to EP19216170.1A priority Critical patent/EP3671646A1/en
Priority to CN201911346233.6A priority patent/CN111345835B/zh
Publication of JP2020099662A publication Critical patent/JP2020099662A/ja
Application granted granted Critical
Publication of JP7334048B2 publication Critical patent/JP7334048B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • 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/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4007Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units
    • A61B6/4014Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units arranged in multiple source-detector units
    • 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]
    • A61B6/035Mechanical aspects of 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/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/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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/10Image enhancement or restoration using non-spatial domain filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/60Image enhancement or restoration using machine learning, e.g. neural networks
    • 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/42Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4241Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20064Wavelet transform [DWT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/408Dual energy
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/441AI-based methods, deep learning or artificial neural networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Optics & Photonics (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

本発明の実施形態は、X線CTシステム及び方法に関する。
コンピュータ断層撮影(CT)システム及び方法は、特に医用撮像及び医用診断の為に、幅広く使用されている。CTシステムは、一般的に被検体の身体に関する一つ又は複数の断面的なスライスの画像を作成する。X線管などの放射源は、一側面から身体にX線を照射する。身体を通過してきた放射線の減衰は、検出器から受け取った電気信号を処理することで測定される。係る信号は、その後逆ラドン変換(又はその同等のもの)を実行することにより、身体の画像を再構成するために使用される。
エネルギー積分型検出器と光子計数検出器とは、デュアルエネルギーCT(dual-energy:DE)(「スペクトラルCT」と呼ばれることも)を生成するための個別のCT構成で使用することができる。DECTは、物質弁別を実行するために有利に使用することができ、それにより骨を身体における軟部組織と見分けることができ、より多くの臨床情報をもたらす。CTにおけるスペクトラル撮像に対して、様々な構成を使用することができる。一般的に、スペクトラルCT構成を、次の二つのタイプへと分類することができる。(i)異なるエネルギースペクトルを生成するX線源で(例えば、高速kV(kilo-voltage)切替、そしてデュアル線源構成)と、(ii)X線の個別のエネルギー(例えば、デュアルレイヤ検出器)に基づいて検出されたX線を見分ける検出器と、である。ここの、「kVp(kVのピーク値)」及び「kV」という言葉は、X線源からのX線放射のスペクトルの制御用の機構として、X線源(例えばX線管)へと適用される、管電圧の切り替えのコンテキストにおいて、同じ意味で使用される。
より具体的に、臨床的な意義の、DECT構成が三つ存在する―デュアルレイヤ検出器システム、デュアル線源システム、そして高速kV切替システムである。例えば、デュアルレイヤ検出器は、エネルギー分解検出器のうちの一タイプである。デュアルレイヤ検出器は、検出されたX線を異なるエネルギー帯へと分離する機能を実行するフィルタのように、個別のエネルギー積分型検出器の前面に配置された様々なX線エネルギーフィルタを使用する。
別のDECT構成が、デュアルX線源システムである。デュアルX線源システムにおいて、二つのX線源は向かい合って個別の検出器に配置され、各線源検出器ペアはそれ独自のCTシステムを形成している。係る各線源検出器ペアは、異なるX線スペクトルを用いて他のCTシステムと同時に操作する。例えば、二つのCTシステムは、ガントリを伴う同じ回転環状リングへと固定することができるが、互いに関し正しい角度で配置することもできる。従って、係る配置を以て、二つの異なるX線スペクトルで同時に、二つのCTスキャンを実行することができる。
第三のDECT構成において、単一の積分源は、患者の周りを回転するCTスキャナのビュー角として、高エネルギーX線スペクトルと低エネルギーX線スペクトルとを素早く切り替えるための、高速kV切替を使用するX線源(例えば、X線管)を用いて使用することができる。
これらの三つの各DECT構造(つまり、デュアルレイヤ、デュアル線源、そしてkV切替)は、それ自体に欠点がある。デュアルレイヤ検出器システムにおいて、シンチレータと光電子増倍管とを組み合わせることで、低エネルギーノイズや、帯域外のエネルギー抑制を達成するために不十分な最適化に悩まされる。その上、二つの読出しスペクトル間の著しい重複により、エネルギー分離が低下する。デュアル線源システムは、第二のシステムの追加されたコストや、またクロス散乱効果にも苦しむ。高速kV切替システムは、高エネルギー投影データ及び低エネルギー投影データと並行して収集するために使用される、極超短波発生器及びパラレルデータ収集システム(data acquisition systems:DAS)に対する必要条件が原因で、追加でのコストも相当負う。
従って、上記で特定された欠陥を解消する、より良いスペクトラルCTアプローチが望まれるのである。
米国特許第7995702号明細書
本発明が解決しようとする課題は、各DECT構造で特定された欠陥を解消するスペクトラルCTを実現することである。
実施形態に係るX線CTシステムは、取得部及び処理部を備えている。
前記取得部は、離散的な第1のビュー群に対応し、第1のエネルギーのX線の照射を伴って収集された第1の投影データセットと、前記第1のビュー群とは異なる離散的な第2のビュー群に対応し、前記第1のエネルギーとは異なる第2のエネルギーのX線の照射を伴って収集された第2の投影データセットと、を取得する。
前記処理部は、前記第1の投影データセット及び前記第2の投影データセットに基づいて、第3のビュー群に対応し、前記第1のエネルギーのX線に対応する第3の投影データセット、及び第4のビュー群に対応し、前記第2のエネルギーのX線に対応する第4の投影データセットを生成する学習済みモデルに対して、前記第1の投影データセット及び前記第2の投影データセットを入力することで、前記第3の投影データセット及び前記第4の投影データセットを生成する。
一実施形態に係る、CTスキャナにおけるX線源及び検出器の配置の概要的なダイアグラムを示している。 一実施形態に係る、高速kV切替のためのX線源の投影角度/ビューの例を示している。 一実施形態に係る、高速kV切替のための時間の関数として、kVのピーク値(kVp)の例を示している。 一実施形態に係る、疎なkV切替のためのX線源の投影角度/ビューの例を示している。 一実施形態に係る、疎なkV切替のための時間の関数として、kVpの例を示している。 一実施形態に係る、低及び高kV電圧設定で放射された、X線のエネルギースペクトルに対する確率密度関数のプロットを示している。 一実施形態に係る、ファントムに対する第一の物質成分(例えば水)の放射線濃度断面を示している。 一実施形態に係る、ファントムの図5Aにおける物質成分のサイノグラムを示している。 一実施形態に係る、ファントムに対する第二の物質成分(例えば骨)の放射線濃度断面を示している。 一実施形態に係る、ファントムの図6Aにおける物質成分のサイノグラムを示している。 一実施形態に係る、X線管に関し80kV設定に対するX線エネルギースペクトルに対応するサイノグラムを示している。 一実施形態に係る、X線管に関し135kV設定に対するX線エネルギースペクトルに対応するサイノグラムを示している。 一実施形態に係る、図7Aに図示のサイノグラムに対応する疎なビューサイノグラムを示している。 一実施形態に係る、図7Bに図示のサイノグラムに対応する疎なビューサイノグラムを示している。 一実施形態に係る、逐次近似応用被曝低減(adaptive iterative dose reduction:AIDR)3次元(three-dimensional:3D)再構成法を使用して再構成された画像を示しており、画像は、80kVpの低kV値に対する疎なkV切替投影データから再構成されている。 一実施形態に係る、AIDR-3D法を使用して再構成された画像を示しており、係る画像は、135kVpの高kV値に対する疎なkV切替投影データから再構成されている。 一実施形態に係る、図5Aと5Bとにおける画像の合計を示す。 一実施形態に係る、kV切替投影データ(例えば、疎なビューkV切替投影データ)のサイノグラム補填を実行するために、DL人工ニューラルネットワーク(artificial neural network:ANN)を学習且つ使用するためのフローダイアグラムの例を示している。 一実施形態に係る、3次元サイノグラム補填/補正を実行するための、二つの2次元ニューラルネットワークを使用する、2.5次元フィルタリングアプローチのフローダイアグラムを示している。 一実施形態に係る、3次元サイノグラム補填を実行するための、2.5次元アプローチの別のフローダイアグラムを示している。 一実施形態に係る、ANNベースのサイノグラム補填への周波数分割アプローチのフローダイアグラムを示している。 一実施形態に係る、サイノグラム復元を実行するための、DL ANNを学習且つ使用するためのフローダイアグラムの例を示している。 一実施形態に係る、画像強調を実行するための、DL ANNを学習且つ使用するためのフローダイアグラムの例を示している。 一実施形態に係る、高kV及び低kVサイノグラムに関し個別に単一のチャネルサイノグラム補填を実行するための、二つのDL ANNを学習且つ使用するためのフローダイアグラムの例を示している。 一実施形態に係る、損失エラー関数を最適化するための、DL-ANNネットワークの係数を逐次的に調整することによる、DL-ANNネットワークの学習のフローダイアグラムを示している。 一実施形態に係る、順伝播ANNの例を示している。 一実施形態に係る、「畳み込みニューラルネットワーク(convolutional neural network:CNN)」と呼ばれる、ANNのタイプの例を示している。 一実施形態に係る、疎なkV切替のためのCTスキャナの第一の例を示している。 一実施形態に係る、疎なkV切替のためのCTスキャナの第二の例を示している。 一実施形態に係る、再構成回路の構成の一例を示している。
以下、図面を参照しながら実施形態について述べる。実施形態は、kV(kilo-voltage)切替を用いて生成されたデュアルエネルギー(dual-energy:DE)投影データから、コンピュータ断層撮影(computed tomography:CT)画像を再構成するための、例えば、深層学習(deep learning:DL)ネットワークの使用に関する。また、実施形態は、例えば、エラーを減らし、且つ欠けている投影角度(ビュー)で投影値を埋めるためにサイノグラム補填(sinogram completion)で役立つDLネットワークの使用に関する。具体的には、実施形態に係るX線CTシステムは、取得部及び処理部を備えている。
ここで、取得部は、離散的な第1のビュー群に対応し、第1のエネルギーのX線の照射を伴って収集された第1の投影データセットと、当該第1のビュー群とは異なる離散的な第2のビュー群に対応し、当該第1のエネルギーとは異なる第2のエネルギーのX線の照射を伴って収集された第2の投影データセットと、を取得する。
処理部は、第1の投影データセット及び第2の投影データセットに基づいて、第3のビュー群に対応し、第1のエネルギーのX線に対応する第3の投影データセット、及び第4のビュー群に対応し、第2のエネルギーのX線に対応する第4の投影データセットを生成する学習済みモデルに対して、第1の投影データセット及び第2の投影データセットを入力することで、第3の投影データセット及び第4の投影データセットを生成する。
このような取得部及び処理部の詳細については後述する。
なお、処理部は、複数の検出器素子のX線放射の強度を表す投影データであって、第1の投影データセット及び第2の投影データセットを構成する投影データを取得部から取得する、としてもよい。投影データは、投影ビューの関数として低kV投影データと交互配置された高kV投影データを表しており、投影データにおいて高kV投影データの投影ビューは低kV投影データの投影ビューから分離されている。低kV投影データはX線源へと適用された第一の電圧を使用して収集されている。高kV投影データはX線源へと、第一の電圧に比べてより大きい第二の電圧を使用して収集されている。
また、処理部は、高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第一のニューラルネットワークである学習済みモデルを取得する、としてもよい。第一のニューラルネットワークは、入力データと対象データとを含む学習データセットを使用して学習されている。入力データはkV切替サイノグラムに基づいており、入力データにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されている。対象データは完全な低kVサイノグラムと完全な高kVサイノグラムとを含む完全なサイノグラムに基づいており、完全なサイノグラムのそれぞれが全ての投影ビューに対する投影データを含んでいる。
処理部は、高kV投影データと低kV投影データとに対するサイノグラム補填を実行する、としてもよい。処理部は、第一のニューラルネットワークに対して、低kV投影データに基づいた低kV入力と高kV投影データに基づいた高kV入力とを適用し、第一のニューラルネットワークは、サイノグラム補填からのエラーが低減された結果を出力している、よう構成されている、としてもよい。
補足すると、kV切替を用いて収集された投影データにおいて、高kV投影データは、低kV投影データと交互に配置されている。それは、一連の投影角度を通してCTスキャナが進む際に、高kV及び低kV電圧設定のうちのどちらか一方に対し所定の時間間隔で、X線源(例えば、X線管)にわたり適用される電圧が切り替えられるからである。このようにして、低kV電圧設定が適用される間の投影角度に対し、低kV投影データが収集されるが、高kVサイノグラムにおいて、高kV投影データに対する投影データは欠如している。従って、欠如している高kV投影データ(即ち、低kV電圧設定が適用された間の投影角度に対する高kV投影データ)を埋めるために、補間を使用することができる。同様に、欠如している低kV投影データを埋めるためにも、補間を使用することができる。この処理は、「サイノグラム補填(sinogram completion)」と呼ばれる。しかしサイノグラム補填は、完璧ではなく、またアーチファクトや画質の低下に繋がることがある。従って、サイノグラム補填を実行するための及び/又はサイノグラム補填における不完全を補正するための、改善された方法が望まれる。
ここに説明される方法と装置とは、サイノグラム補填に対して関連したアプローチにおいて、上記で述べたアーチファクトと画像劣化の提示とを解消する。これらの欠陥は、例えばkV切替投影データのより優れたサイノグラム補填を達成するために、高kV投影データの付加的な情報を使用することができる、二つのチャネル深層学習(deep learning:DL)人工ニューラルネットワーク(artificial neural network:ANN)を適用することにより、解消される。(i)サイノグラム補填ステップで又は直後で、(ii)基底/物質弁別ステップで又は直後で、(iii)画像再構成ステップで又は直後で、を含むCT画像再構成処理に沿った一つ又は複数のステップで、不完全なサイノグラム補填に関して補正するために、二つのチャネルDL-ANNを適用することができる。「投影データ」及び「サイノグラムデータ」の用語は、同じものを意味し、互いに置き換え可能に使用される。
ここで、(i)の場合、処理部は、例えば、低kV入力と高kV入力とを第一のニューラルネットワークへと適用することと一体的に、サイノグラム補填を実行するよう更に構成されてもよい。ここで、第一のニューラルネットワークへの高kV入力は、サイノグラム補填よりも前の高kV投影データである。第一のニューラルネットワークへの低kV入力は、サイノグラム補填よりも前の低kV投影データである。第一のニューラルネットワークは、サイノグラム補填を実行するために、高kV投影データと低kV投影データとの相補的な情報を使用するよう学習され、全ての投影ビューに対する投影データを含む完全な高kV投影データと、全ての投影ビューに対する投影データを含む完全な低kV投影データとを第3の投影データセット及び第4の投影データセットとして生成している。
あるいは、処理部は、例えば、第一のニューラルネットワークに対して、低kV入力と高kV入力との適用よりも前にサイノグラム補填を実行するよう更に構成され、全ての投影ビューに対する投影データを含む完全な高kV投影データと全ての投影ビューに対する投影データを含む完全な低kV投影データとを生成するサイノグラム補填を実行する。また、処理部は、高kV入力としての完全な高kV投影データと低kV入力としての完全な低kV投影データとを第一のニューラルネットワークに適用して、第一のニューラルネットワークから出力される完全な高kV投影データ及び完全な低kV投影データ内でサイノグラム補填エラーが低減された結果を出力する、としてもよい。
(ii)の場合、処理部は、例えば、高kV投影データと低kV投影データとを第一の基底投影データと第二の基底投影データとに弁別し、第一の基底投影データから第一の基底画像を再構成し、第二の基底投影データから第二の基底画像を再構成するよう更に構成された、としてもよい。
ここで、処理部は、第一の基底投影データ用の一方の入力チャネルと、第二の基底投影データ用の他方の入力チャネルとの二つのチャネルを有する第二のニューラルネットワークを取得し、第二のニューラルネットワークは、kV切替サイノグラムの基底弁別に基づいている入力データと、完全なサイノグラムの基底弁別に基づいている対象データと、を含む学習データセットを使用して学習されている、としてもよい。
あるいは、処理部は、高kV投影データと低kV投影データとを、物質弁別を使用して、第一の物質成分投影データと第二の物質成分投影データとに弁別するよう更に構成された、としてもよい。
なお、DL-ANNは、二つのチャネルに限らず、三つのチャネルを有してもよい。例えば、投影データは、投影ビューの関数として交互配置された高kV投影データと低kV投影データと中kV投影データを表す場合がある。この場合、中kV投影データは、例えば、第一の電圧と第二の電圧との間の切り替え途中の電圧を使用して収集されている。また、処理部は、投影データを取得し、高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルと、中kVデータ用の第三の入力チャネルとの三つの入力チャネルを有する第一のニューラルネットワーク(DL-ANN)を取得するよう構成されている場合がある。このような場合、三つのチャネルを有するDL-ANNが用いられる。
加えて、サイノグラムデータは、3次元で表示することができる。例えば、これらの3次元は、X線検出器アレイの2次元に対応する二つの空間次元を含むことができ、またサイノグラムデータの3次元は、CTスキャンの投影角度とすることができる。これらの3次元は、真の(true)3次元畳み込みニューラルネットワーク(convolutional neural network:CNN)を使用して明らかにすることができ、2次元(2D)畳み込みは、3次元畳み込みに比べてより早く計算することができる。その結果、2D CNNフィルタリングを第一の面に対して平行なスライスへとまず適用し、その後係る第一の面に対し直交である第二の面に対して平行なスライスへと2D CNNフィルタリングを適用することにより、全ての3Dに関して明らかにするために、2.5次元アプローチを使用することができる。つまり、一つの3D CNNを使用するよりも、二つの2D CNNを使用することにより3D効果を達成することができ、計算と時間との数を減らすことになる。
更に、DL-ANNネットワークの学習及び活用は、投影データを高及び低周波数成分へと二分することで、加速することもできる。係る高及び低周波数成分は、その後二つの個別のCNN―高周波数CNN及び低周波数CNN―へと適用される。
その上、特定の実施形態において、高kV投影データと低kV投影データとは、二つの個別のDL-ANNネットワークへと個別に適用することができる。例えば、DL-ANNは、チャネルのうちの一方が投影データであり、第二のチャネルはマスクである、二つのチャネルネットワークとすることができる。係るマスクは、サイノグラムにおけるどのピクセルが測定されたデータであり、またどれがサイノグラム補填処理により埋められるのかを示す。加えて又は代わりに、マスクに依存する重み付けを適用することにより、DL-ANNを学習するために、損失/コスト関数においてマスクを使用することができる。この場合、例えば、処理部は、[a]第一のニューラルネットワークを取得し、当該第一のニューラルネットワークは、kV切替低kVサイノグラムが収集された投影ビューに対してkV切替高kVサイノグラムが収集された投影ビューを表すkV切替マスクに対する第二のチャネルを有し、[b]第二のニューラルネットワークを取得し、第二のニューラルネットワークは、kV切替マスク用の第二のチャネルを有し、[c]第一のニューラルネットワークに対して、取得された投影データの高kV入力とkV切替マスクとを適用し、[d]第二のニューラルネットワークに対して、取得された投影データの低kV入力とkV切替マスクとを適用し、[e]第一のニューラルネットワークは、kV切替サイノグラムに対応する個別のkV切替マスクを含む入力データを含む第一の学習データセットを使用して学習されており、[f]第二のニューラルネットワークは、kV切替サイノグラムに対応する個別のkV切替マスクを含む入力データを含む第二の学習データセットを使用して学習されている、よう更に構成されていてもよい。
上記特性のおかげで、ここに説明される方法は、その他のデュアルエネルギーCT(dual-energy:DECT)システム及びアプローチに勝る利点を幾つか有する。例えば、背景技術で述べられたアプローチは、例えば、より大きなハードウェアコスト及び/又は再構成された画像の画質に対する低下を含む、様々な欠陥に悩まされる。例えば、ここに説明される方法は、高速kV切替よりも、疎なkV切替を使用して実現可能な高画質画像を作ることにより、ハードウェアコストを減らすことができる。即ち、ここに説明される方法は、高速kV切替に対して改善された画質を可能とするのみでなく、疎なkV切替に対しても可能とする。疎なkV切替を実現可能なものとすることで、疎なkV切替に対するハードウェアがよりコスト的に有効な場合があるので、ここに説明される方法は、ハードウェアコストを著しく下げることができる。
高速kV切替とデュアル/検出器-線源システムとの両方において、より高いコストが欠陥である。ここに説明される方法は、X線管に適用されるkVがまれに切り替えられる(つまり、kV切替が疎であり、結果的に低及び高kV値両方に対して疎なビュー投影データとなる)、疎なkV切替と相性が良い。従って、投影角度(ビュー)の各変化に対する低及び高kV値間を切り替えるよりも、ここで使用されるkV切替は疎であり、CTスキャナが他方のkV設定へと切り替わる前にいくつかの投影角度を通して回転しながら、所定のkV設定が維持されるようにして、kV切替の頻度が少なく実行されることを意味している。即ち、高kV設定へと切り替えられた後に、低kV設定へと再度切り替えられる前に、スキャナが多くの投影角度を通って回転し、且つそこで投影画像を収集しながら、X線源は高kV電圧を維持する。この低kV設定への再切替は、次のいくつかの投影角度を通し、その後維持され同様に続く。このようにして、単一のデータ収集システム(data acquisition system:DAS)を使用することができ、X線管にわたり電圧を、より遅い速度で且つkV設定切替間をより長い期間で切り替えるには、より簡素且つより安価なハードウェアで十分となる。
更に、ここに説明される方法は、デュアルレイヤ検出器の欠陥を解消する。それは、X線検出器でX線エネルギーをフィルタリングするよりも、X線源にわたって適用された電圧を調整/切替することにより、二つのエネルギースペクトルが達成されるからである。
高速kV切替及び疎なkV切替デュアルエネルギーCT(DECT)システムの両方を用いて、ここに説明される方法を使用することができる。しかし、疎なkV切替システムに関するコスト節約により、考察は疎なkV切替システムに主眼を置いている。疎なkV切替システムは、それぞれ低及び高エネルギースペクトルに対する、疎なビュー投影データを生成する。kV切替は、高速kV切替に比べ低頻度で実行されるので、高周波数発生器(高速kV切替に対して使用される極超短波発生器とは対照的)を使用して、ここに説明される方法を実行することができる。その上、単一の順次DAS(高速kV切替に対して使用されるパラレルDASとは対照的)を使用して、疎なkV切替を実行することができる。
DECTへの疎なkV切替アプローチに対する重要な問題のうちの一つに、疎なビュー投影データが再構成された画像の画質及び物質弁別に関して問題を提起するということがある。特に、疎なkV切替投影データに対して、ストリーク状アーチファクト、ビームハードニング、そして画質を低下させるようなその他の影響を受けにくい、効果的な再構成アルゴリズムを発達させるのが難しいということが判明している。一方で、フィルタ補正逆投影(filtered back-projection:FBP)法等の解析的再構成法が効果的なこともあるが、疎な投影データに適用された場合には、ストリーク状アーチファクトの影響に悩まされる再構成された画像を生成する。他方で、逐次再構成法が疎な投影データに適用された場合に、ストリーク状アーチファクトは減る可能性があるものの、各kV設定でのX線のより低い線量が原因で、空間分解能が低下し、且つノイズテクスチャが低下するという結果になる。kV設定は、X線管線源の陽極と陰極との間に適用されるピークキロ電圧のことである。「キロ電圧」は略してkV、そして「ピークキロ電圧」は略してkVpと、ここでは同義として使用される。
疎な切替投影データについてのもう一つの問題は、検出器ピクセルに対してX線により辿られた軌道が二つのkV設定間で異なるものの、サイノグラムドメイン物質弁別は、二つのkV設定に対して重複するX線軌道を必要とする点である。画像ドメイン物質弁別は、他方で、X線ビーム(例えば、ボウタイフィルタを通過した異なるパスが原因)のエネルギースペクトルにおけるビームハードニング補正及び空間的変動を含む、それ自体の難しさを一連で有する。このようにして、厳密なサイノグラム補填/補正は、サイノグラムドメインにおける厳密な物質弁別を可能とする重要な役割を果たす。そうでなければ、サイノグラム補填における不完全さが生き残る可能性があり、物質弁別処理及び画像再構成処理を通して悪化することもある。
次に図面を参照しながら、参照番号が数枚の図にわたって同一又は対応する部分を指し示す。図1は、エネルギー積分型検出器を有するコンピュータ断層撮影(CT)スキャナにおけるX線源及び検出器の構成を示している。図1が提供するのは、寝台116に横たわる画像化される被検体OBJ、X線源112、コリメータ/フィルタ114、そして個別の検出器素子のアレイ(つまりピクセル)を含む、画素化されるX線検出器103中の相対的な位置の非限定例である。
一実施形態において、X線源112とコリメータ/フィルタ114とは、ガントリに回転可能に接続されている回転コンポーネント110に固定して接続されている。例えば、係る回転コンポーネント110は、被検体OBJが寝台116上の空間に固定されたまま、ガントリ内で回転するよう構成された環状リングとすることができる。又は代替的に、ヘリカルスキャンにおいて、寝台は、X線源112及びX線検出器103がガントリのボアの周りを回転しながら、ガントリのボアに沿って平行移動することもできる。CTスキャナのガントリは、回転コンポーネント110のアイソセンタで中心とすることができるボアにおいて、開口アパーチャ115も含む。係る開口アパーチャ115は、被検体OBJがX線源からX線の投影面に位置されることを可能とする。特定の実施形態において、X線検出器103は、当該ガントリに回転可能に接続された別の回転コンポーネント130に固定して接続されている。回転/回転構成において、回転コンポーネント110と回転コンポーネント130とは、投影角度(つまりビュー)の進行で被検体OBJの投影データを取得するために、X線検出器103がX線源112と正反対に位置した状態を維持しながら共に回転することができる。一軸に沿って配置された投影角度での投影データと、その他の軸に沿って配置された投影データの空間的次元とを配置することで、サイノグラムは作り出される。
スペクトラルCTにおいて、多重エネルギー成分を有する投影データが、被検体OBJの投影測定を表すために使用される。これらの投影測定は、スペクトラルでないCTと同様の従来的なCT画像再構成法を可能にする、一連の角度で行われる。しかし、スペクトラルでないCTとは違い、スペクトラルCTは、投影測定を物質成分への弁別を可能とする、付加的な情報(つまり、スペクトラル減衰情報)を生成する。kエッジ法を除けば、コンプトン散乱及び光電子減衰がそれぞれ原因により、X線減衰のユニークなスペクトラル特徴に基づいて、物質数は大抵が二つである。つまり、二つの物質成分に対するX線減衰間のスペクトルの差は、係る物質成分が示す光電子減衰に対するコンプトン散乱の異なる割合から生じる(例えば、ヨードの様な高Z物質が原因のX線減衰は、水等の低いZ物質に比べて、光電子減衰に対するコンプトン散乱の異なる割合から成り立っている。)
サイノグラムドメイン(つまり画像再構成処理前)又は画像ドメイン(つまり画像再構成処理後)のどちらかに、スペクトラルドメインから物質ドメインへと投影データのマッピングを実行することができる。しかし、サイノグラム(投影)ドメインにおいて実行される場合に、係る投影データは、デュアルエネルギー成分のそれぞれに対する同一の(又はほぼ同一の)X線軌道を含まねばならない。kVp切替において、X線軌道間の差は、高kV設定が使用された投影ビューでは低kV投影データが欠如し、またその逆も然りという事実から生じる。
生体物質におけるX線の減衰は、二つの物理的プロセス(つまり光電子吸収とコンプトン散乱)によって占められている。このようにして、エネルギーの関数としての減衰係数は、以下の弁別によって近似することができる。
ここでμPE(E,x,y)は光電子減衰であり、μc(E,x,y)はコンプトン減衰である。代わりにこの減衰係数は、
となるように、高Z物質(つまり物質1)と低Z物質(つまり物質2)との弁別へと再整理することができる。ここでc1(x,y)及びc2(x,y)は、第一と第二の物質成分にそれぞれ対応する。物質弁別とは、測定された/再構成された減衰スペクトルと最も近似するc1(x,y)及びc2(x,y)から解く処理である。
図2A及び2Bは、高速kV切替の実施形態のダイアグラムを示している。X線管についてのkV設定は、各投影角度に対して変えられる。図2Aにおいて、高kV(低kV)設定を使用する投影画像の収集に対するX線源の位置が、白(黒)マルによって示されている。図2Bにおいて、X線管に適用された電圧は、白(黒)マルによって示された高kV(低kV)設定で収集された投影画像に対する時間及び電圧と共に、時間の関数として示されている。
図3A及び3Bは、疎なkV切替の実施形態のダイアグラムを示している。X線管についてのkV設定は、異なる投影角度で一連のN投影画像の後でのみ変更される。ここでNは、1よりも大きい数である。図3Aにおいて、高kV(低kV)設定を使用する投影画像の収集に対するX線源の位置は、白(黒)マルによって示されている。図3Bにおいて、X線管に適用された電圧は、白(黒)マルによって示された高kV(低kV)設定で収集された投影画像に対する時間及び電圧と共に、時間の関数として示されている。
図3A及び3Bの非限定例において、Nは3であるが、Nは2以上の大きい任意の整数とすることができる。更に、所定のkV設定で収集された一連の投影画像数は、CTスキャンを通して定数である必要はなく、また異なる間隔が異なるkV設定(例えば、所定のCTスキャンにおいて、低いkVp設定の時に比べて、高いkVp設定の時に、より多くの投影画像を収集することができるし、またその逆も然り)に対して適用することができる。
疎なkV切替設定の方法は、80kVの低kV設定及び135kVの高kV設定の非限定例を使用して、ここで描かれている。X線管に対して、X線スペクトルは、電子が陰極に衝突することにより突然止められ、電子の運動エネルギーを、制動放射線を経てX線へと変換する前に、電子を加速するために陽極及び陰極の間に適用された電圧(kV)で主に制御される。この処理により、異なるX線スペクトルは、X線管にわたり適用される電圧を変えることにより、生み出すことができる。図4は、80kVの低kV設定及び135kVの高kV設定を用いて生み出されたX線エネルギーの確率分布をそれぞれ示している。
ここに説明される方法は、二つの異なるX線スペクトルを生成するために、二つの電圧がX線源へと適用される、非限定的な例を使用して描かれている。しかし、ここに説明される方法は、三つ以上の異なるX線スペクトルを生成するために、X線源へと三つ以上の電圧設定が適用される使用にも適用可能、且つ使用も含む。例えば、kV切替は、高kV設定、中kV設定、そして低kV設定を使用して実行することができる。当該例において、例えば、高kV及び低kV設定を使用して描かれる場合に対して、下記に説明される二つのチャネルニューラルネットワークの代わりに、三つのチャネルネットワークを使用することができる。三つ以上のX線スペクトルを使用することにより、物質弁別に関する付加的な情報を提供するという利点がある。従って、三つ以上の電圧設定がX線源へと適用されるkV切替での実施形態は、ここに説明される方法の趣旨からは乖離することはない。
図5A及び6Aは、非限定例として使用される、ファントムの二つの物質成分に関する横断面ビューを示す。当該非限定例において、二つの物質成分は、水と骨とである。図5B及び図6Bは、図5A及び6Aに示された横断面に対応する水と骨とのサイノグラムを示す。これらの図5B及び6Bは、完全なサイノグラムである(つまり、データが何も欠如していないということ)。空間次元は横軸に沿って、且つ投影角度は縦軸に沿って、それぞれ表されている。物質成分サイノグラムは、80kV及び135kVのkV設定に対応する、図7A及び7Bに示されたサイノグラムを弁別することにより、取得することができる。
図7A及び7Bに示された完全なサイノグラムとは対照的に、図7C及び7Dは、各投影角度に対して、投影データが80kVのkV設定又は135kVのkV設定かのどちらかを使用して収集された、交互配置された疎なビュー投影データを示す。80kV及び135kVサイノグラムのそれぞれに対して欠如する投影角度で減衰値を埋めるために、サイノグラム補填法(つまり、補間又は修復)を使用することができる。
疎なビュー投影データは、高又は低kV設定かのどちらかに対応するとして投影角度の間隔を切り替えるよう設計し、その後高kV設定が適用される間隔に対する低kVサイノグラム値を削除(そして逆もまた同様)することで、完全なサイノグラムから近似することができる。高kV設定から低kV設定へ(逆もまた然り)の遷移領域において、切替処理の過渡効果(transient effects)を近似するために、高kV及び低kVサイノグラムの重畳を使用することができる。
上記で述べた通り、疎なkV切替を使用して収集された投影データを用いての画像再構成の問題は、ストリーク状アーチファクトが原因で、画質に低下が見られる傾向にあることである。例えば、図8A及び8Bは、80kVと135kVとに対する疎なkV切替投影データを用いてFBP再構成を使用して生成された、再構成された画像の例をそれぞれ示している。ストリーク状アーチファクトが鮮明に出ている。その上、厳密な検査により、ストリーク状アーチファクトが相補的(complementary)な傾向にあることを明らかとする。つまり、80kV画像における明るいストリークは、135kV画像における暗いストリークに対応し、逆もまた然りということである。これを描くために、図8Cは、80kV画像と135kV画像との合計を示す。かなりの程度にまで、80kV画像における明るいストリーク状アーチファクトが、135kV画像における暗いストリーク状アーチファクトに反作用する。これらのストリーク状アーチファクトは、図7Cと7Dとに示される通り、80kV及び135kV投影データに対する投影角度が相補的/交互配置であるという事実から生じるものとして、理解することができる。例えば、CTスキャナが高kV設定へと設定された投影角度は、CTスキャナが低kV設定へと設定された投影角度を補うものである。このようにして、個別の高及び低kV設定での二つの画像からの情報は、相補的となる。
ここに説明される方法は、不完全なサイノグラム補填から結果として生じる、未処理のエラーに関してサイノグラム補填及び/又は補正を実行するために、個別の低及び高kV投影データから組み合わされた情報を、如何に使用するかについて学習するために、DL-ANNを使用する。また、学習済みのDL-ANNを学習済みモデルとして使用する。
ここで、学習済みモデルと、学習の処理について順に述べる。
学習済みモデルは、CTガントリ等の取得部に通信可能な処理部に用いられる。処理部は、学習済みモデルに対して、取得部により取得された第1及び第2の投影データセットを入力することで、第3及び第4の投影データセットを生成する。詳しくは、取得部は、離散的な第1のビュー群に対応し、第1のエネルギーのX線の照射を伴って収集された第1の投影データセットと、当該第1のビュー群とは異なる離散的な第2のビュー群に対応し、当該第1のエネルギーとは異なる第2のエネルギーのX線の照射を伴って収集された第2の投影データセットと、を取得する。
処理部は、第1の投影データセット及び第2の投影データセットに基づいて、第3のビュー群に対応し、第1のエネルギーのX線に対応する第3の投影データセット、及び第4のビュー群に対応し、第2のエネルギーのX線に対応する第4の投影データセットを生成する学習済みモデルに対して、第1の投影データセット及び第2の投影データセットを入力することで、第3の投影データセット及び第4の投影データセットを生成する。
なお、第3のビュー群及び第4のビュー群は、それぞれ第1のビュー群及び第2のビュー群を含む、としてもよい。
第3のビュー群は、第1のビュー群より多く、第4のビュー群は、第2のビュー群より多い、としてもよい。
第1の投影データセット、第2の投影データセット、第3の投影データセット、及び第4の投影データセットは、それぞれサイノグラムである、としてもよい。
第3の投影データセット及び第4の投影データセットは、それぞれ画像再構成に必要な全てのビューに対応するサイノグラムである、としてもよい。
ここで、入力される第1及び第2の投影データセットと、生成される第3及び第4の投影データセットとしては、さまざまなバリエーションがある。さまざまなバリエーションのうち、以下の5通りの場合(1)~(5)を例に挙げて説明する。但し、さまざまなバリエーションは、以下の5通りの場合(1)~(5)に限定されない。また、以下の表を用いた説明中、記号「L」は、低kV(第1のエネルギーの例)に対応する。記号「H」は、高kV(第2のエネルギーの例)に対応する。なお、低kVが第2のエネルギーに対応し、高kVが第1のエネルギーに対応するように変形してもよい。各々の表中、1000個のビュー番号と、ビュー番号50個毎の範囲は、一例であり、これに限定されない。また、各々の表中、「第1、第2の投影データセット」の行における記号「L」に交差する列のビュー番号の複数の範囲(例、51-100, 151-200, ・・・ ,951-1000)は、第1のビュー群に対応する。「第1、第2の投影データセット」の行における記号「H」に交差する列のビュー番号の複数の範囲(例、1-50, 101-150, ・・・ ,901-950)は、第2のビュー群に対応する。
第1の場合(1)は、学習済みモデルに対して、交互配置された第1の投影データセット(L)及び第2の投影データセット(H)を入力することで、全ビューの第3の投影データセット(L)及び全ビューの第4の投影データセット(H)を生成する場合である。換言すると、離散的なビュー群の入力に対し、全ビューの生成物が得られる場合である。
第1の場合(1)には、第3の投影データセットは、第1のビュー群と第2のビュー群(全ビュー)に対応し、第1のエネルギー(L)に対応する投影データセットである。第4の投影データセットは、第1のビュー群と第2のビュー群(全ビュー)に対応し、第2のエネルギー(H)に対応する投影データセットである。
第2の場合(2)は、学習済みモデルに対して、交互配置された第1の投影データセット(L)及び第2の投影データセット(H)を入力することで、離散的な第3の投影データセット(L)及び離散的な第4の投影データセット(H)を生成する場合である。換言すると、離散的なビュー群の入力に対し、H,Lを反転した生成物が得られる場合である。
第2の場合(2)には、第3の投影データセットは、第2のビュー群(記号「H」の列)に対応し、第1のエネルギー(L)に対応する投影データセットである。第4の投影データセットは、第1のビュー群(記号「L」の列)に対応し、第2のエネルギー(H)に対応する投影データセットである。
第3の場合(3)は、学習済みモデルに対して、ブランクを介して交互配置された第1の投影データセット(L)及び第2の投影データセット(H)を入力することで、離散的な第3の投影データセット(L)及び離散的な第4の投影データセット(H)を生成する場合である。換言すると、離散的なビュー群の入力に対し、当該入力のビュー群にブランクのビュー群を含めた生成物が得られる場合である。
第3の場合(3)には、第3の投影データセットは、第のビュー群(記号「」の列)と、第1及び第2のビュー群とは異なるビュー群(例、中間のビュー群)とに対応し、第1のエネルギー(L)に対応する投影データセットである。第4の投影データセットは、第のビュー群(記号「」の列)と、第1及び第2のビュー群とは異なるビュー群(例、中間のビュー群)とに対応し、第2のエネルギー(H)に対応する投影データセットである。
第4の場合(4)は、学習済みモデルに対して、ブランクを介して交互配置された第1の投影データセット(L)及び第2の投影データセット(H)を入力することで、離散的な第3の投影データセット(L)及び離散的な第4の投影データセット(H)を生成する場合である。換言すると、離散的なビュー群の入力に対し、当該入力のビュー群とは異なるビュー群(ブランクのビュー群)の生成物が得られる場合である。
第4の場合(4)には、第3の投影データセットは、第1及び第2のビュー群とは異なるビュー群(例、中間のビュー群)とに対応し、第1のエネルギー(L)に対応する投影データセットである。第4の投影データセットは、第1及び第2のビュー群とは異なるビュー群(例、中間のビュー群)とに対応し、第2のエネルギー(H)に対応する投影データセットである。
第5の場合(5)は、学習済みモデルに対して、ブランクを介して交互配置された第1の投影データセット(L)及び第2の投影データセット(H)を入力することで、全ビューの第3の投影データセット(L)及び全ビューの第4の投影データセット(H)を生成する場合である。換言すると、離散的なビュー群の入力に対し、全ビュー(入力のビュー群より多い)の生成物が得られる場合である。
第5の場合(5)には、第3の投影データセットは、全ビュー(第1のビュー群と第2のビュー群より多い)に対応し、第1のエネルギー(L)に対応する投影データセットである。第4の投影データセットは、全ビュー(第1のビュー群と第2のビュー群より多い)に対応し、第2のエネルギー(H)に対応する投影データセットである。
以上が、さまざまなバリエーションのうちの5通りの場合の説明である。前述した通り、入力される第1及び第2の投影データセットと、生成される第3及び第4の投影データセットとのバリエーションは、前述した5通りの場合(1)~(5)に限定されない。
また、「離散的な第1のビュー群」及び「第1のビュー群とは異なる離散的な第2のビュー群」における「離散的」の定義は、管電圧の変化に対応していればよい。例えば、「離散的」の用語は、以下の場合(a)~(d)を含んでいる。ここで、管電圧が高kVの場合を「H」で表し、低kVの場合を「L」で表す。
(a)管電圧が、H→L→H→L→H→…で変化する場合(交互に変化する場合)。
(b)管電圧が、H→H→H→L→L→L→H→H→H→…で変化する場合(一部連続する場合)。
(c)管電圧が、H→M→L→M→…で変化する場合(高kVと低kVとの間の中kV(記号「M」)のビュー群を挟む場合)。前述したブランクのビュー群を含む場合も同様である。なお、中kVは、高kVと低kVとの間の切り替え途中の電圧である。但し、中kVは、これに限らず、高kVと低kVとの間の所定電圧としてもよい。いずれにしても、中kVを用いる場合、例えば、三つのチャネルを有するDL-ANNが用いられる。
(d)管電圧が、H→H→L→H→H→…で変化する場合(高kVと低kVの頻度(ビュー番号の範囲)が異なる場合)。
以上が、「離散的」の用語が含む場合(a)~(d)についての説明である。但し、「離散的」の用語は、これらの場合(a)~(d)に限定されない。
図9Aは、DL-ANN361の学習の処理310を含み、且つ最終的に高画質物質成分画像257を生成するのに疎なビュー投影データ251を補正するために、学習されたアーチファクト補正ネットワーク361を使用する処理202を含む、疎なkV切替CT画像再構成のための方法200のフローダイアグラムを示す。
処理310において、損失関数は、DL-ANNネットワーク361(例えば、DL-ANNネットワーク361のパラメータは、重み係数接続ネットワーク層、層におけるノードの活性化関数/ポテンシャルを含む場合がある)のパラメータを逐次的に調整/最適化するために使用される。学習されたネットワーク361を生成するために、停止基準が満たされるまで(例えば、所定のしきい値に対するネットワーク係数/パラメータの収束)、ネットワークパラメータの最適化は続く。損失関数は、完全なサイノグラム353(個別のkV設定のそれぞれで収集されたフルスキャン)と、疎なビュー収集357(例えば、疎なkV切替画像)がDL-ANNネットワーク361の現在のバージョンへと適用された場合に、ネットワークによる結果出力と、を比較する。疎なビューサイノグラム357は、投影角度の間隔が何ら減衰値を有さない、図7C及び7Dに示された様なサイノグラムとすることができる。代わりに、疎なビューサイノグラム357は、次のようなサイノグラムとすることもできる。それは、疎なビューサイノグラム357がDL-ANNネットワーク361へと適用されていることに先立って、不完全なサイノグラム補填法(例えば、補間又は修復)が既に適用されているものである。第一の場合に、DL-ANNネットワーク361は、欠如するサイノグラムデータを埋めるために学習され、第二の場合に、係るDL-ANNネットワーク361は、不完全なサイノグラム補填法の欠点に関して補正するために学習されている。
DL-ANNネットワーク361は、対象データ(即ち、完全なサイノグラム353)が二つのkV設定(即ち、図7A及び7Bに示されるような、高kV設定及び低kV設定)に個別に対応する完全なサイノグラムのペアを含む、二つのチャネルネットワークとすることができる。入力データ(即ち、疎なビューサイノグラム357)は、二つのkV設定へと個別に対応するサイノグラムのペアも含むが、対象データと比較される結果を生成するためにDL-ANNネットワーク361へと適用される入力データは、疎なビューサイノグラム(図7C及び7Dに示される通り)であるか、又は疎なビューサイノグラム357について、初期ではあるが不完全なサイノグラム補填法の実行結果であるか、のどちらかである。一般的に、処理310を実行するために使用される学習データセットは、個別の対象サイノグラムとペアになった多数の入力サイノグラムのセットを含むものである。完全なサイノグラムを「対象」又は「ラベル」と、そして疎なビューサイノグラムを「入力」と、それぞれ呼ぶことができる。学習データは、対応する対象及び入力の多数のセットを含むことができる。特定の実施形態において、オフライン学習は、所定のCTスキャン(例えば、CTスキャナが最初にインストールされ且つキャリブレートされた場合)より前に実行することができ、学習されたDL-ANNネットワーク361は、CTスキャンが収集され且つ画像再構成がステップ220で実行されるまで、メモリに格納することができる。一般的に、疎なビューサイノグラム357は、どのようにkV切替投影データ251が取得され、そしてその後方法200のステップ210において前処理されるのかと同じ方法で、又は整合するようにして準備される。このようなやり方で、ステップ210から結果として生じる前処理されたサイノグラムで生じる、特定のアーチファクト及び欠陥を補正するために、DL-ANNネットワーク361が学習される。
図9Aにおいて、高速か又は疎なkV切替かのどちらかを使用して生成されたCT投影データ251は、高画質な基底/成分画像257を生成するために、ステップ210、220、230、240を使用して処理される。基底/成分画像257は、物質成分画像(例えば、物質成分は、骨及び水の場合がある)であることが多い。しかし、物質成分基底の代わりに、その他の基底も使用することができる。例えば、基底弁別は、デュアルエネルギー画像を、コンプトン散乱減衰基底と光電効果減衰基底とに弁別することができる。基底弁別は、単なる物質弁別に比べてより一般的であるが、物質弁別の非限定例は、方法200の様々な変形例を描くために、ここで多く使用される。
特定の実施形態において、CT投影データ251は、ステップ210で前処理(例えば、信号プリコンディショニング、キャリブレーション補正、ベースライン補正、ビームハードニング補正等)されたCTスキャンから収集された、投影データとすることができる。特定の実施形態において、ステップ210で実行される前処理は、初期サイノグラム補填(例えば、補間)を含む。投影データ251は、様々なキャリブレーション及びジオメトリック因子(factor)を使用して補正される、サイノグラムとすることができる。例えば、前処理は、検出器オフセットに対するベースラインサブトラクション、検出器中の増幅器ゲイン及び量子効率における変形例に対する補正、X線フラックスの関数として非線形検出器応答に対する補正等の補正を含むことができる。その上、これらの補正は、キャリブレーションデータ、経験的に導出されたパラメータ、アプリオリ公知のパラメータに基づく場合がある。
処理202のステップ220において、ステップ210から生じるサイノグラムは、処理315から学習されたDL-ANNネットワーク361へと適用される。特定の実施形態において、処理315からのDL-ANNネットワーク361は、ステップ210で実行される初期サイノグラム補填における欠点に対して補正する。
上述の通り、高及び低kV疎なビューサイノグラム間の相補的な情報を利用するために、DL-ANNネットワーク361を二つのチャネルネットワークとすることができる。高及び低エネルギーデータは、互いに相補的であるから、両チャネルからの相互的、相補的な情報を解析する二つのチャネルネットワークを使用することで、各エネルギー成分(kV設定)を別々に捉える一つのチャネルネットワークを上回る、大幅な利益をもたらす。
従って、処理202のステップ220において、ネットワーク361は、低及び高エネルギー投影データ間の相補的な情報を利用する二つのチャネルネットワークとすることができる。特定の実施形態において、処理315で実行される学習は、ステップ210と同じ前処理技術を使用して生成された、入力画像を使用する。
処理202のステップ230において、基底/物質弁別は、ステップ220から生じるサイノグラムに関して実行される。一般的に、基底/物質成分サイノグラムを生成するために、任意の基底/物質弁別法を使用することができる。
処理202のステップ240において、逆投影法、フィルタ補正逆投影法、フーリエ変換ベースの画像投影法、逐次画像再構成法、逆行列画像再構成法、統計学的画像再構成法、又は当業者には公知のその他の再構成法を使用して、基底/物質成分サイノグラムに関して、画像再構成を実行することができる。例えば、再構成法は、ヘリカル再構成技術、コーンビーム再構成技術、フェルドカンプアルゴリズム、FBP再構成法、及び画像とサイノグラムドメインとのうちの一方又は両方におけるノイズ低減を伴う逐次近似応用被曝低減(adaptive iterative dose reduction:AIDR)三次元(3D)再構成法を使用することができる。再構成は、デノイジング、高減衰領域における光子不足(starvation)を最小化するための補正、量子ノイズを軽減するための補正、エッジ保存/強調補正、そしてフィルタリング(線形フィルタリングと非線形フィルタリング、そしてデノイジング)を含むことができる。
図9B及び9Cは、ステップ210からの前処理されたサイノグラム221を、学習されたDL-ANNネットワーク361へと適用するための、ステップ220に対するフローダイアグラムを示す。これに伴い、処理部は、例えば、[a1]高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第一のニューラルネットワークを取得し、[a2]高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第二のニューラルネットワークを取得し、[b]低kV投影データと高kV投影データとを第一の面に沿った個別の第一のスライスへと区切り、[c]部分的に補正された低kV投影データと部分的に補正された高kV投影データとを生成するために、個別の第一のスライスを第一のニューラルネットワークへと適用し、[d]部分的に補正された低kV投影データと部分的に補正された高kV投影データとを第一の面と直交である第二の面に沿った個別の第二のスライスへと区切り、[e]補正された低kV投影データと補正された高kV投影データとを生成するために、個別の第二のスライスを第二のニューラルネットワークへと適用する、よう更に構成された、としてもよい。
図9B及び9Cにおいて、DL-ANNネットワーク361は、サイノグラムを処理するために順番に使用される、二つのネットワーク(即ち、第一のニューラルネットワークと第二のニューラルネットワークとに)へと二分されている。サイノグラムは、三次元―検出器アレイの次元に対応する二つの空間次元(x、y)と、投影データが収集される投影角度θに対する三次元と―を有する。
3Dサイノグラムブロックを使用することは、三つの次元(ビュー、セグメント、チャネル)全てにおける特徴の連続性(feature continuity)により、欠如するkV測定を満たすので、好ましい。しかし、3Dサイノグラムブロックを学習し且つ使用することは、特徴空間(feature dimensions)が増えたことにより、kV切替サイノグラムがステップ220で適用された場合に、学習ステップ間の緩慢な収束及び時間のかかる計算という結果になる。このようにして、ステップ220は、二つのステップへと分けることが可能で、各ステップは、サイノグラムの三次元の異なるペアに対応する、異なる2Dネットワークを使用する。このニューラルネットワークの、「二.五次元(2.5D)アプローチ」と呼ぶことのできる2ステップアプリケーションは、サイノグラムの三つ全ての次元における情報を引き続き生かしながらも、フル3Dアプローチに関するコンピュータ的な負担を減らすことができる。
2.5Dアプローチは、デュアルkVサイノグラムをサイノグラムの2Dスライス(例えば、ビュー-チャネルスライス)へと分割することで始まり、第一のニューラルネットワークは、これらの第一の2Dスライスを使用する最適なサイノグラム補填を達成し、且つ部分的に補正/フィルタされた途中段階(intermediary)のサイノグラムを生成するために、学習される。この部分的に補正/フィルタされたサイノグラムは、更なるサイノグラム補填を補正/実行するために学習された第二のネットワークへとその後適用される、2Dスライス(例えば、チャネル-セグメントスライス)の第二のセットを生成するために、別の面に沿ってスライスされる。
ステップ222において、高及び低kVサイノグラムは、第一の面に対し平行な2D面へとスライスされる。例えば、図9Cにおいて3Dサイノグラムは、ビュー-チャネル面へとスライスされる/区切られる。ここで、空間次元xは「チャネル次元」と、そして空間次元yは「セグメント次元」と、それぞれ呼ばれる。
ステップ223において、2D畳み込みニューラルネットワーク(convolutional neural network:CNN)である第一のニューラルネットワークは、ステップ222からの各スライスへと適用され、係るスライスは高及び低kV設定に対する部分的にフィルタされたサイノグラム224を生成するために、再び組み合わせられる。第一のニューラルネットワークは、二つのチャネルにより提供された相補的な情報を利用するために、二つのチャネルネットワーク(即ち、低kVサイノグラムに対応する一つのチャネル及び高kVサイノグラムに対応する他方のチャネル)とすることもできる。
ステップ225において、部分的にフィルタされたサイノグラム224は、第二の面に対し平行な2D面へとスライスされる。例えば、図9Cにおける部分的にフィルタされたサイノグラム224は、チャネル-セグメント面に対してスライスされる/区切られる。
ステップ226において、2D CNNである第二のニューラルネットワークは、ステップ225からの各スライスへと適用され、係るスライスは高及び低kV設定に対する完全にフィルタされたサイノグラム227を生成するために、再び組み合わせられる。第二のニューラルネットワークは、二つのチャネルにより提供された相補的な情報を利用するために、二つのチャネルネットワーク(即ち、低kVサイノグラムに対応する一つのチャネル及び高kVサイノグラムに対応する他方のチャネル)とすることもできる。
2D畳み込みは、3D畳み込みに比べてずっと早く実行することができるので、二つの2D CNNは、一つの3D CNNに比べてより早く実行することができる。その上、学習の間、二つの2D CNNは、3D CNNに比べてずっと早く、別々に収束することがある。学習のため、完全なサイノグラムは、第一の及び第二の2D CNNのそれぞれの学習用の対象データとして、使用することができる。
ビュー-チャネル面にある第一のスライスと、チャネル-セグメント面にある第二のスライスとの選択は、説明の目的のために、非限定例として提供される。一般的に、第一のスライスと第二のスライスとに対し、任意の二つの面を選ぶことができる。
図9Dは、入力サイノグラムが、低周波数成分と高周波数成分とに二分される、ステップ220のフローダイアグラムを示している。これに伴い、処理部は、例えば、[a1]高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第一のニューラルネットワークを取得し、[a2]高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第二のニューラルネットワークを取得し、[b]低kV投影データを、高空間周波数低kV投影データと低空間周波数低kV投影データとを含む低kV周波数成分へと変換し、[c]高kV投影データを、高空間周波数高kV投影データと低空間周波数高kV投影データとを含む高kV周波数成分へと変換し、[d]第一のニューラルネットワークに対して、高kV周波数成分の高周波数成分と低kV周波数成分の高周波数成分とを適用し、[e]第二のニューラルネットワークに対して、高kV周波数成分の低周波数成分と低kV周波数成分の低周波数成分とを適用し、[f]補正された低kV投影データと補正された高kV投影データとを生成するために、個別の高周波数成分を、対応する低周波数成分と組み合わせるよう更に構成された、としてもよい。
ニューラルネットワークの学習において、低周波数成分は、高周波数成分に比べてより素早く収束する傾向にある。その上、低周波数成分は、より低い解像度へとダウンサンプルすることができ、計算を速め、低周波数成分に対する畳み込み層についての認知分野(perception fields)は、高周波数ニューラルネットワークに対する認知分野とは異なる(例えば、より大きい)場合がある。従って、サイノグラムを高及び低周波数成分に分け、且つサイノグラムの低周波数成分と共に使用するために低周波数成分ニューラルネットワークを学習することには、利益がある。これに伴い、処理部は、例えば、[g]ダウンサンプルされた投影データを生成するために、低空間周波数低kV投影データと低空間周波数高kV投影データとをダウンサンプリングすることと、[h]ダウンサンプルされた結果を生成するために、ダウンサンプルされた投影データを第二のニューラルネットワークへと適用することと、[i]補正された低空間周波数低kV投影データと補正された低空間周波数高kV投影データとを生成するために、ダウンサンプルされた結果をアップサンプルすることと、により低周波数成分を第二のニューラルネットワークへと適用するよう更に構成された、としてもよい。
従って、特定の実施形態において、ネットワークを低周波数成分と高周波数成分とに二分することにより、ネットワークの実行/収束及びネットワークの効率が改善される場合もある。その後、異なる周波数成分に対する二つの個別のネットワークを別々に学習することで、深層学習が行われる。例えば、この二分周波数アプローチは、第一にサンプルされたデュアルkVサイノグラム(例えば、ステップ220でのサイノグラム補填)を埋め且つ作り出すために、各高及び低kVサイノグラムを別々に補間することにより、実行することができる。第二に、低周波数成分を取得するために、ガウスのフィルタが各kVサイノグラムへと適用され、高周波数成分は、完全にサンプルされたデュアルkVサイノグラムとその低周波数成分との間の差としてその後取得される。第三に、二つの個別の深層学習ネットワークが学習される。(i)二つの低周波数成分サイノグラムに対する低周波数ニューラルネットワーク、(ii)二つの高周波数成分サイノグラムに対する高周波数ニューラルネットワーク、である。最後に、ステップ220におけるデュアルkVサイノグラム補填処理で、疎なデュアルkVサイノグラムが低及び高周波数成分へと分割され、その後個別の低及び高ニューラルネットワークへと適用される。個別のニューラルネットワークからの低及び高周波数の結果は、ステップ220から結果として生じるフィルタされたサイノグラム227を生成するために、その後再び組み合わされる。
ステップ228において、前処理されたサイノグラム221のそれぞれが、高及び低周波数成分へと弁別される。上述の通り、特定の実施形態において、この周波数二分は、低周波数成分サイノグラムを取得するために前処理されたサイノグラム221をローパスフィルタリングし、その後高周波数成分サイノグラムを取得するために対応する前処理されたサイノグラム221から低周波数成分サイノグラムを引くことで、実行することができる。代わりに、特定の実施形態において、周波数二分は、前処理されたサイノグラム221のそれぞれを高及び低周波数成分へと分けるために、ウェーブレット変換を使用して実行することもできる。加えて、更にその他の実施形態において、周波数二分は、開示の趣旨から乖離することなく、高及び低周波数成分を分けるために、任意のその他の変換又は周波数成分分離法を使用して実行する場合もある。言い換えると、処理部は、例えば、[s1]低空間周波数投影データを生成するために、投影データに対してローパスフィルタを適用し、且つ高空間周波数投影データを生成するために、低空間周波数投影データを投影データから引くことと、[s2]投影データを低空間周波数投影データと高空間周波数投影データとに分けるために、ウェーブレット変換を投影データへと適用することと、[s3]投影データを周波数ドメインへと変換し、且つ周波数ドメインにおける投影データを低空間周波数投影データと高空間周波数投影データとに分けることと、のうちの一つにより投影データを周波数成分へと変換するよう更に構成された、としてもよい。
ステップ229において、高及び低kVサイノグラム221の高周波数成分は、高周波ニューラルネットワークへと適用され、且つ高及び低kVサイノグラム221の低周波数成分は、低周波数ニューラルネットワークに併せて適用される。その後、個別のニューラルネットワークから結果として生じる高及び低周波数サイノグラムは、フィルタされた高及び低kVサイノグラム227を生成するために、組み合わされる(例えば、足すことで)。
図10は、DL-ANNネットワーク362が、ステップ230における基底弁別の後、且つステップ240における再構成前に実行されるステップ232で使用される、修正された方法200に対するフローダイアグラムを示している。基底弁別処理の非線形の本質が原因で、サイノグラム補填後の未処理エラーは、弁別の間に増幅することができ、再構成された基底画像における相当のエラーへとつながりかねない。これらのエラーを取り除くために、ここに説明された方法は、物質弁別されたドメイン(例えば、物質成分ドメイン)において、深層学習ベースのサイノグラム復元を実行することができ、それにより画像再構成前に基底弁別されたサイノグラムにおけるエラーを最小化する。例えば、このサイノグラム復元は、ステップ230から結果として生じる基底成分サイノグラムを、処理315を使用して学習されたDL-ANNネットワーク362へと適用することで、ステップ232で実行することができる。
ステップ220において、基底弁別の前のフルデュアルkVサイノグラムを作り出すために、低及び高kVサイノグラムについて、サイノグラム補填が実行される。ステップ220は、図9Aを参照して説明された通り、DL-ANNネットワーク361を使用して実行することができるが、又は補間アプローチ、修復アプローチ、或いはその他の方法を使用してサイノグラム補填を実行することもできる。
ステップ230において、物質/基底弁別は、上述の通り実行される。DL-ANNネットワーク362を学習するために使用される入力サイノグラム354は、ステップ230から結果として生じる低及び高kVサイノグラムと、同じ又は同様の処理を使用して、準備することができる。対象サイノグラムも、対象データがkV切替サイノグラムからというより完全なサイノグラムから準備されるので、サイノグラム補填無しで準備されることを除き、同様の方法で準備されている。
ステップ232において、ステップ230から結果として生じる基底成分サイノグラムは、学習されたDL-ANNネットワーク362へと適用される。
図11は、DL-ANNネットワーク363が、ステップ240での再構成後にステップ242で使用される、修正された方法200に対するフローダイアグラムを示している。再構成された基底画像の画質を更に改善するため、ここに説明される方法は、画像ドメインにおける深層学習ベースのポスト処理法を使用する場合もある。ステップ210-240の画像再構成処理は、上述の通り実行される。(例えば、ステップ220でDL-ANNネットワーク361を使用して任意に、且つステップ232でDL-ANNネットワーク362を使用して任意に)。
入力画像355は、同じ処理ステップを使用してkV切替データから準備され、また対象画像359は、サイノグラム補填に対してユニークである及び/又は不完全なサイノグラムデータから導出されるエラーの補正を除き、同じ処理ステップを使用して完全なサイノグラムデータから準備される。処理312において、DL-ANNネットワーク363は、入力画像と対象画像とを使用して学習される。これに伴い、処理部は、例えば、第一の基底画像用の一方の入力チャネルと、第二の基底画像用の他方の入力チャネルとの二つのチャネルを有する第二のニューラルネットワークを取得し、第二のニューラルネットワークは、入力データと対象データとを含む学習データセットを使用して学習されており、入力データは、kV切替サイノグラムの基底弁別から再構成された画像に基づいており、対象データは、完全なサイノグラムの基底弁別から再構成された画像に基づいているよう更に構成された、としてもよい。
ステップ242において、ステップ240からの基底成分画像257は、強調された基底成分画像258を生成するために、学習されたDL-ANNネットワーク363へと適用される。
図12は、二つの個別のDL-ANNネットワーク361Aと361Bとが、サイノグラム補填に対して実行/補正するためにステップ220A及び220Bで個別に使用される、修正された方法200に関するフローダイアグラムを示している。図12におけるフローダイアグラムは、低及び高kVサイノグラムが個別のチャネルへと適用される、二つのチャネルネットワークを使用して一緒に処理されているというよりも、別々のDL-ANNネットワーク361Aと361Bとを使用して、高kVサイノグラムと低kVサイノグラムとが別々に処理されることを除き、図9Aと同じである。つまり、処理310Aにおいて、DL-ANNネットワーク361Aは、入力データを生成するために高kV kV切替サイノグラムを使用し、且つ対象データとして高kVでの完全なサイノグラムを使用して、学習される。その上、処理310Bにおいて、DL-ANNネットワーク361Bは、入力データを生成するために低kV kV切替サイノグラムを使用し、且つ対象データとして低kVでの完全なサイノグラムを使用して、学習される。
DL-ANNネットワーク361Aと361Bとは、二つのkV切替サイノグラムのうちの一方又はその他のみを使用するので、DL-ANNネットワーク361Aと361Bとのそれぞれは、kV切替サイノグラムのその他によって提供される相補的な情報へとアクセスしない。それにも関わらず、DL-ANNネットワーク361A(361B)は、第二のチャネルとしてのサイノグラムと同じ次元であるマスクを適用することにより、二つのチャネルネットワークを使用することができる。前処理された高kV(低kV)サイノグラムは、第一のチャネルとして適用される。例えば、マスクは、オリジナルの高kVサイノグラムにおいて所定のピクセルが不在だったことを示すのに「1」の値を、そして存在したことを示すのに「0」を、それぞれ利用する場合がある。この付加的な情報は、より良いパフォーマンスを達成するために、DL-ANNネットワーク361A(361B)を助けることができる。マスクは、付加的な情報を提供するための第三のチャネルとして、DL-ANNネットワーク361、362、363のうちのその他の一つ以上に適用することができ、それによりパフォーマンスを改善する。この場合、処理部は、高kV投影データが収集された投影ビューに対して低kV投影データが収集された投影ビューを表すkV切替マスクを取得するよう更に構成され、第一のニューラルネットワークはkV切替マスクを受信するための第三の入力チャネルを有する、としてもよい。
加えて又は代わりに、ニューラルネットワーク(複数のニューラルネットワーク)のオフライン学習の間に使用される、損失関数における重みを生成するために、マスクを適用することができる。例えば、kV切替サイノグラムに初めは不在だったものの、サイノグラム補填処理の間に続いて埋められたピクセル値は、常に存在し且つサイノグラム補填処理によって変更されなかったサイノグラムの値に比べて、より多く重みづけることができる。この重み付けは、サイノグラムの埋められた領域におけるパフォーマンスを最適化する、DL-ANNネットワークに対する変更を優遇する(favoring)ことで、収束を改善することができる。
上述での、1及び0の値の二値マスクに加えて、距離マスクも使用する場合がある。距離マスクは、欠如する測定値と最も近い欠如していない測定値との間の距離を表すものかもしれない。より大きなマスク値は、最も近い測定された値からより遠く離れた欠如する測定値を示す。例えば、kV切替サイノグラムの測定された値は、マスクにおけるピクセルが「0」の値を割り当てることがある。上述の通り、如何に二値マスクが使用されるのか、と同様の方法で距離マスクを使用することができる。これに伴い、処理部は、アーチファクト低減ニューラルネットワークを学習させるよう更に構成されていてもよい。kV切替マスクは、kV切替サイノグラムにおける欠損値に対して、kV切替サイノグラムにおいて最も近い収集された値に対する近似を更に表す、ようにしてもよい。
要約すると、デュアルエネルギーCT(DECT)は、放射線科医により、高度な撮像ツールとして発展を遂げ、益々受け入れられつつある。従来的な単一のエネルギーCTスキャナと比べて、DECTスキャナは、複数の違った物質が単一のエネルギーCTの下では同じHU値を表すにも関わらず、それらを区別する能力を有する。DECTの実施形態は、主に次の三つのタイプへとカテゴライズができる―デュアル管kV切替、デュアルX線源、そしてデュアルレイヤ検出器、である。kV切替収集法に関して、最も安定した一つのkVで、各ビューで利用が可能である。サイノグラム空間における物質弁別を実行するために、二つのkVからの測定が収集されることを考えると、物質/基底弁別に先立ちサイノグラム補填が実行される。
モデルベースの逐次画像再構成等、不完全な測定の課題を解消するための関連するアプローチは、時間がかかるか、又は限られたアプリケーション(例えば、高速kV切替等高価なハードウェアが使用された場合)に対してのみ適切であるか、のどちらかである。深層学習(DL)は、画像修復や補填タスク等のアプリケーションに対して期待でき、関連する修復法に比べてより優れた結果を達成するであろう。関連する方法と比較して、DLは明確なモデリングを必要とせず、代わりに入力データにおける不備を補正するのに最適なアプローチを学習するための学習データに依存する。DLは従って、ピクセル間の内在する相関関係を、局所的にもグローバル的にも、捉える能力があり、対象(例えば、グラウンドトゥルース)とペアになった入力(例えば、破損した所見(corrupted observations))を含む学習データセットを使用してネットワークを学習することにより、高度なネットワークを構築することができる。計算的な時間やハードウェアコストを減らし、且つデュアルkV切替画質を更に改善するため、ここに説明される方法は、デュアルkV切替収集を伴うデュアルエネルギーCTにおける欠如する測定を補填するために、DLアプローチを使用する。
サイノグラム補填、物質/基底弁別、そして画像再構成の処理における様々なステップで、DL-ANNネットワークを使用することができる。例えば、DL-ANNネットワークは、サイノグラム補填から結果として生じるエラーや不完全さを減らすために、サイノグラム補填の間に、又は直後に使用することができる。その上、DL-ANNネットワークは、kV切替間に収集される不完全なサイノグラムから結果として生じる、未処理の不完全さが原因のエラーとアーチファクトとを減らす/補正するために、基底弁別及び/又は画像再構成の後の期間に(during after)、使用することもできる。これらのDL-ANNネットワークは、より優れたサイノグラム補填とCT画像から再構成されたCT画像に対する改善された画質とを達成するために、高及び低kVサイノグラムの相補的な情報を活用することができる。
上で記された特徴により、幾つかの利点を達成することができる。第一に、DL-ANNアプローチは、高速kV切替実行に対して最適なだけでなく、遅い(つまり、疎な)kV切替が採用された場合に、即ち、一つ以上の連続的なビューが個別のkV設定に対して欠如している場合にも最適である。つまり、上記特徴は、高速kV切替の実施形態と比較した場合に、より費用を抑えたハードウェアを用いて実行することができる、疎なkV切替の実施形態を可能とする。
第二に、DL-ANNアプローチにより実現される改善は、欠如しているデータをリカバーし、又は画像を再構成するための逐次再構成法を使用することとは対照的に、再構成が逆投影アプローチを使用して実行することが可能なデータを十分改善することができる。逆投影アプローチを使用することで、計算的な時間及びハードウェアコストを抜本的に減らすことが可能となる。
第三に、DL-ANNアプローチは、単一のエネルギーCTに比べて、意図的にDECTに向けて作られている。結果として、ここに説明される方法は、単一のエネルギーCTに向けたアプローチにおいて、見当たらなかった三つの利点を有する。第一に、特定の実施形態において、ここに説明される方法は、3D処理の利益を維持しながら、計算的な時間を減らす、二つのステップ、2.5Dサイノグラム処理戦略を採る。第二に、ここに説明される方法は、オフライン学習の効率及び収束を改善するために、周波数の二分を使用することもある。第三に、ここに説明される方法は、パフォーマンスを改善するために、サイノグラムマスクを用いて使用することもある。
次に、DL-ANNネットワークの学習の、より詳細な説明が提供される(例えば、処理315)。この説明は、完全なサイノグラム353である対象データと、疎なビューサイノグラム357である入力データと、の非限定例を使用して描かれているが、一般的には任意の対象データ及び入力データを使用することができる。以下の説明は、例えば、[a]サイノグラム補填ニューラルネットワークのネットワーク係数を初期化することと、[b]kV切替サイノグラムに基づいた高kV入力と低kV入力とのペアを含む入力データを取得することと、kV切替サイノグラムにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されており、kV切替低kVサイノグラムは、投影ビューの関数としてkV切替高kVサイノグラムと交互配置される、[c]完全な低kVサイノグラムと完全な高kVサイノグラムとを含む完全なサイノグラムに基づいた高kV対象と低kV対象とのペアを含む対象データを取得することと、完全なサイノグラムのそれぞれは全ての投影ビューに対する投影データを含む、[d]対象データと、入力データをサイノグラム補填ニューラルネットワークへ適用することからの出力との間の一致を表す損失関数を最適化するために、ネットワーク係数を逐次的にアップデートすることと、により、二つのチャネルを有するサイノグラム補填ニューラルネットワークを学習させるよう構成された処理部を具備するX線CTシステムに関する。
この処理部は、[e]kV切替高kVサイノグラムが収集された投影ビューに対してkV切替低kVサイノグラムが収集された投影データを表すkV切替マスクを取得することと、[f]損失関数を最適化するために、ネットワーク係数を逐次的にアップデートすることと、によりアーチファクト低減ニューラルネットワークを学習させるよう更に構成された、としてもよい。
また、処理部は、[g]アーチファクト低減ニューラルネットワークを学習させるよう更に構成され、kV切替マスクは、kV切替サイノグラムにおける欠損値に対して、kV切替サイノグラムにおいて最も近い収集された値に対する近似を更に表す、としてもよい。以下、このような処理部の処理315について具体的に説明する。
図13は、学習処理315の一実施形態のフローダイアグラムを示している。処理315において、低クオリティデータ357と高クオリティデータ353とは、DL-ANNネットワークを学習させるための学習データとして使用され、処理315のステップ319からの出力であるDL-ANNネットワークという結果になる。一般的に、データ357は、任意の不備呈示するサイノグラム又は画像の場合があり、係る「不備」とは、画像処理(例えば、ノイズ又はアーチファクト)を通して影響を受ける可能性のある、任意の望まれない特徴という場合がある。同様に、データ353は、「対象データ」や「不備低減データ」、「不備最小化データ」、又は「最適化されたデータ」と呼ばれることがあり、係る「不備」とは、画像357に比べてより少ない。オフラインDL学習処理315は、入力サイノグラム357から対象サイノグラム353に似ている画像を生み出すようDL-ANNネットワーク361を学習させるために、対応する対象サイノグラム353とペアにされている大量の入力サイノグラム357を使用して、DL-ANNネットワーク361を学習させる。
処理315は、学習データのセットが取得され、ネットワーク361は、エラー(例えば、損失関数により生み出された値)を減らすために逐次的にアップデートされる。言い換えると、対象データと、入力データをサイノグラム補填ニューラルネットワークへ適用することからの出力との間の一致を表す損失関数を最適化するために、ネットワーク係数が逐次的にアップデートされる。DL-ANNネットワークは、学習データにより暗示されたマッピングを推察し、コスト関数は、対象サイノグラム353と、DL-ANNネットワーク361の現在の形(incarnation)を入力サイノグラム357へと適用することにより生み出された結果と、の間のミスマッチに関連したエラー値を生み出す。例えば、特定の実施形態において、コスト関数は、平均二乗エラーを最小化するために、平均二乗エラーを使用する場合がある。多重層パーセプトロン(multilayer perceptrons:MLP)ニューラルネットワークの場合に、(確率論的)勾配降下法を使用して平均二乗エラーベースのコスト関数を最小化することで、ネットワークを学習させるために、逆伝播アルゴリズムを使用することができる。
処理315のステップ316において、DL-ANNネットワーク361の係数に対して、初期推測が生成される。例えば、初期推測は、画像化されている領域のアプリオリ知識、又は一つ又は複数の例示的なデノイジング法、エッジ保存検出法、及び/又はブロッブ検出法に基づく場合がある。加えて、初期推測は、LeCun初期化、Xavier初期化、Kaiming初期化のうちの一つに基づく場合もある。ステップ316の処理は、サイノグラム補填ニューラルネットワークのネットワーク係数を初期化することに対応する。
処理315のステップ316から319までは、DL-ANNネットワーク361を学習するための初期化法の非限定例を提供する。
ネットワーク361の現在のバージョンの適用後に、対象サイノグラム353(つまり、グラウンドトゥルース)と入力サイノグラム357との間の差(つまり、距離測定)の測定を表すために、エラーが計算される(例えば、損失関数又はコスト関数を使用して)。エラーは、上述のコスト関数を含む、画像データ間の任意の公知のコスト関数又は距離測定を使用して、計算することができる。その上、特定の実施形態において、エラー/損失関数は、一つ又は複数のヒンジ形損失関数(hinge loss)又は交差エントロピー損失関数を使用して計算することができる。
加えて、学習データに示された特定の例に対してネットワークを過学習させることを避けるために、損失関数を正則化アプローチと組み合わせることができる。正則化は、機械学習問題における過学習の回避に役立つ場合がある。あまりに長く学習された場合に、モデルが十分な表現力(representational power)を有すると仮定すると、ネットワークはそのデータセットに特定のノイズを学習するようになるが、これを「過学習」と呼ぶ。過学習の場合に、DL-ANNは汎化に乏しくなり、ノイズがデータセット間で異なるため、分散が大きいものとなる場合がある。最小合計エラーは、バイアスと分散との合計が最小の場合に発生する。従って、学習データにおけるノイズに特定の解というよりも、学習されたネットワークが一般的な解を表す可能性(likelihood)を最大化するのに最も簡単且つ可能な方法でデータを説明する、局所最小値に到達することが望ましい。この目標は、例えば、早期停止、重み正則化、ラソウ(lasso)正則化、リッジ正則化、又は弾性ネット正則化により、達成することができる。
特定の実施形態において、ネットワーク361は、逆伝播を使用して学習される。逆伝播は、ニューラルネットワークの学習用に使用することができ、また勾配降下最適化法と併せて使用される。前向きメッセージパッシング(forward pass)の間、アルゴリズムは、現在のパラメータθに基づいたネットワークの予測を計算する。これらの予測はその後、対応するグラウンドトゥルースラベル(即ち、高クオリティ画像353)と比較される予測により、損失関数へと入力される。後ろ向きメッセージパッシング(backward pass)の間、モデルは現在のパラメータに関する損失関数の勾配を計算する。その後、パラメータが、最小化された損失の方向における予め決めたサイズのサイズのステップ(step of size of a predefined size)を取ることによりアップデートされる(例えば、ネステロフのモメンタム法や様々な適応法等の加速された方法、ステップサイズは、損失関数を最適化するようより素早く収束するために選択することができる)。
逆投影が実行される最適化法は、勾配降下、バッチ勾配降下、確率論的勾配降下、ミニバッチ確率論的勾配降下のうちの一つ又は複数を使用することができる。更に、最適化法は、例えば、Adagradサブ勾配法、Adadelta又はAdagrad法のRMSPropパラメータアップデートバリエーション、そしてAdam適応最適化技術等のネステロフモーメンタム技術又は適応法を含む、深層学習における確率論的勾配降下のより速い収束速度という結果になる最適化アプローチにおける、一つ又は複数のモーメンタムアップデート技術を使用して、加速することができる。最適化法は、ヤコビ行列(Jacobian matrix)をアップデートステップへと組み込むことで、第二のオーダ法を適用することもできる。
前向きメッセージパッシングと後ろ向きメッセージパッシングとは、ネットワークの個別の層を通して、インクリメントに実行することができる。前向きメッセージパッシングにおいて、実行は、入力を第一の層を通してフィードすることにより開始し、その結果後続の層に対する出力活性化を作り出す。当該処理は、最後の層での損失関数が到達されるまで、繰り返される。後ろ向きメッセージパッシングの間、最後の層は、それ自体の学習可能なパラメータ(もしあれば)に関してと、以前の層に対する上流微分としての機能を果たす、それ自体の入力に関しても計算する。当該処理は、入力層が到達されるまで繰り返される。
図13に戻って、処理315のステップ317は、ネットワークにおける変化の関数としてエラーでの変化を計算することができる(例えば、エラー勾配)と判断し、当該エラーにおける変化は、DL-ANN361の重み/係数に対する後続の変化に対する方向及びステップサイズとを選択するために、使用することができる。この様な方法でエラーの勾配を計算することは、勾配降下最適化法の特定の実施形態と一貫性がある。特定のその他の実施形態において、当該ステップは、省略することができる、及び/又は、当業者にとって理解される通り、別の最適化アルゴリズム(例えば、シミュレートされた焼き鈍し法又は遺伝的アルゴリズムのような、非勾配降下最適化アルゴリズム)に従って、別のステップと置き換えることができる。
処理の315のステップ317において、係数の新たなセットが、DL-ANN361に対して決定される。例えば、重み/係数は、勾配降下最適化法又は過剰緩和加速法のように、ステップ317において計算した変化を使用してアップデートすることができる。
処理315のステップ318において、DL-ANN361のアップデートされた重み/係数を使用して、新たなエラー値が計算される。
ステップ319において、ネットワークの学習が完璧かどうかを判定するために、予め定めた停止基準が使用される。例えば、予め定めた停止基準は、新たなエラー及び/又は実行された総逐次数が所定値を上回るかどうかを評価することがある。例えば、新たなエラーが所定のしきい値を下回って低下したか、又は最大逐次数に到達したかのどちらかの場合に、停止基準を満たすことができる。停止基準が満たされなかった場合に、処理315において実行される学習処理は、新たな重みと係数とを使用してステップ317に戻ったり繰り返したりすることで、逐次ループの開始へと戻り続けるようになる(逐次ループは、ステップ317、318、319を含む)。停止基準が満たされた場合に、処理315において実行される学習処理は、完了である。
図14A及び14Bは、DL-ANNネットワーク361における層の間の相互連結の様々な例を示している。DL-ANNネットワーク361は、完全に連結された層、畳み込み層、そしてプーリング層を含むことができるが、これらは全て下記に説明される。DL-ANNネットワーク361の特定の推奨実施形態において、高レベルの推論(reasoning)を行う完全に連結された層は、損失関数に近づいて構造の更に下に位置される一方で、畳み込み層は入力層の近くに配置される。プーリング層は、畳み込みの後に挿入することができ、且つフィルタの空間的広がりを低下させる低減(reduction)を証明することができ、そしてつまり学習可能なパラメータ量とすることができる。活性化関数も、非線形性を取り込むために様々な層へと組み込まれ、複雑な予測関係を学習するためのネットワークを可能とする。活性化関数は、活性化関数を飽和状態にする(例えば、シグモイド又は双曲正接活性化関数)、又は矯正された活性化関数(例えば、上記で説明された第一の例及び第二の例に適用された、矯正された線形ユニット(Rectified Linear Unit:ReLU)とすることができる。DL-ANNネットワーク361の層は、上記で説明された第一及び第二の例でも例示された通り、バッチ規格化も組み込むことができる。
図14Aは、Nの入力、Kの隠れ層、そして三つの出力を有する、一般的な人工ニューラルネットワーク(artificial neural network:ANN)の例を示している。各層は、ノード(又はニューロンとも呼ばれる)でできており、各ノードは、出力を生成するために、入力の重み付けられた合計を実行し、当該重み付けられた合計の結果をしきい値と比較する。ANNは、可変しきい値、結合重み、又はノード数及び/又はノードの結合性等構造の詳細によって取得されたクラスのメンバーに対する関数のクラスを構成する。ANNにおけるノードは、ニューロン(又はニューロンノードとも)呼ぶことがあり、係るニューロンは、ANNシステムの異なる層の間に相互結合を有することができる。最も簡素なANNは三つの層を有し、オートエンコーダ(autoencoder)と呼ばれる。DL-ANNネットワーク361は、ニューロンの三つ以上の層を有し、また入力ニューロンと同数の出力ニューロン(「~」はxの上にある)を有すことがあり、ここでNは再構成された画像におけるピクセル数である。シナプス(つまり、ニューロン間の結合)は、計算におけるデータをマニピュレートする「重さ」と呼ばれる値(「係数」又は「重み付け係数」と同義で呼ばれる)を格納する。ANNの出力は、パラメータの次の三つのタイプに依存する。(i)ニューロンの異なる層間の相互結合パターン、(ii)相互結合の重さをアップデートする用の学習処理、(iii)ニューロンの重み付けられた入力を係るニューロンの出力活性化へと変換する活性化関数、である。
数学的に、ニューロンのネットワーク関数m(x)は、その他の関数の合成として更に定義することができる、その他の関数の合成n(x)として定義される。これは、図14Aに示されている通り、変数間の従属状態を描いている矢印を用いて、ネットワーク構造として便利に表すことができる。例えば、ANNは、非線形重み付けられた合計を使用することが可能であり、ここでm(x)=K(Σ(x))であり、またここでKは(「活性化関数」として通常呼ばれる)、双曲正接等のある所定の関数である。
図14A(そして同様に図14B)において、ニューロン(つまりノード)は、しきい値関数の周りの円で描かれている。図14Aに示された非限定例について、入力は線形関数の周りの円として描かれ、矢印はニューロン間の方向付けられた結合を指す。特定の実施形態において、DL-ANNネットワーク361は、順伝播型ネットワーク(feedforward network)である。
図14Bは、DL-ANNネットワーク361が畳み込みニューラルネットワーク(CNN)である非限定例を示している。CNNは、画像処理に対し有益な特性を有するANNのタイプであり、従って画像デノイジングのアプリケーションに対し特に関連性がある。ニューロン間の結合性パターンが画像処理における畳み込みを表すことのできる順伝播型ANNを、CNNは使用する。例えば、CNNは、「受容野(receptive fields)」と呼ばれる、入力画像の部分(portion)を処理する小さなニューロン収集の多層を使用することで、画像処理最適化に対して使用することができる。これらの収集の係る出力は、オリジナル画像のより良い表示を取得するために、これらの収集の出力が重なるようにしてその後タイル張り(tiled)にすることができる。この処理パターンは、畳み込み層とプーリング層との交代を有する多層にわたり、繰り返すことができる。
畳み込み層の後に続いて、CNNは、畳み込み層におけるニューロン集団の出力を組み合わせる、局所及び/又はグローバルプーリング層を含むことが可能である。更に、特定の実施形態において、CNNは、各層の終わりに又は各層の後に、適用された点別の非線形性を用いて、畳み込み層及び全結合層の様々な組み合わせも含むこともできる。
CNNは、画像処理に対する利点がいくつかある。フリーパラメータ数を減らし且つ生成を改善するために、入力の小さな領域について畳み込み演算が取り入れられる。CNNの特定の実施形態の重要な強みの一つは、畳み込み層における共有された重みの使用であり、つまり層における各ピクセルに対する係数として使用されるフィルタ(重みバンク)が同じフィルタであるということである。係る重要な強みは、メモリ使用量(memory footprint)を減らし、且つパフォーマンスを向上させる。その他の画像処理法と比較して、CNNは有利に比較的小さな前処理を使用する。これは、従来的なアルゴリズムでは手動で設計されたフィルタを学習することについて、信頼がおけることを意味する。特徴をデザインする際の予備知識及び人的努力に対する依存の欠如が、CNNに対する主たる強みである。
図15は、第三世代ジオメトリで配置されたエネルギー積分型検出器を有する、コンピュータ断層撮影(CT)スキャナの第一の実施形態を示している。ダイアグラムは、X線源112、コリメータ/フィルタ114、X線検出器103、そしてPCD1からPCDNまでの光子計数検出器の間の相対的位置を描いている。
図15に示されたX線源112及びX線検出器103の構成に加えて、投影データを取得するために、X線検出器及びX線源のその他のタイプと組み合わせとを使用することができる。
図15は、X線投影データを収集、格納、処理そして分布用の回路及びハードウェアも示している。係る回路及びハードウェアは、プロセッサ170、メモリ178、そしてデータ収集システム176を含む。
X線源112及びX線検出器103がガントリ140に収容され、且つ回転のコンポーネント110の円軌道の周りを回転する。X線検出器103における検出器素子は、送られたX線放射線を検出し、且つ検出器103が回転するときに検出された信号を出力する。一実施形態において、X線検出器103は、検出器面上の所定のチャネル及びセグメント方向に、密に配置されたエネルギー積分型検出器を有する。
一実施形態において、X線源112は、所定の高レベルエネルギーと所定の低レベルエネルギーとでのX線照射を放射するためのkV切替機能を実行するよう任意で構成された、単一のX線源である。例えば、X線源112は、適用される電圧に依存するエネルギースペクトルを有するX線を照射するよう構成されている。
X線検出器103は、シンチレータ素子と相互作用するX線放射線から結果として生じる、シンチレーションイベントからのシンチレーション光子の結果を検出するために、光電子増倍管又はアバランシェフォトダイオードを用いるシンチレーション素子等のエネルギー積分型検出器を使用することができる。係るシンチレータ素子は、結晶構造、有機液体、プラスチック、またはその他の公知のシンチレータの場合がある。例えば、エネルギー積分型X線検出器は、第一の電圧がX線源112へと適用される場合に、高kV投影データを出力し、且つ第一の電圧に比べて小さい第二の電圧がX線源へと適用された場合に、低kV投影データを出力するよう構成された複数の検出器素子を有する。
CTスキャナは、データチャネルも含む。当該データチャネルは、光子計数検出器及びX線検出器103からの投影測定結果を、データ収集システム176、プロセッサ170、メモリ178に送る。データ収集システム176は、検出器からの投影データの収集、デジタル化、ルーティング(routing)を制御する。データ収集システム176は、環状回転フレーム110の回転を制御する、放射線制御回路も含む。一実施形態において、データ収集システム176は、寝台116の移動やX線源112の操作、そしてX線検出器103の操作も制御するようになる。データ収集システム176は、一連の投影ビューを通して、X線源112とエネルギー積分型X線検出器とを回転させ、かつ第一の電圧と第二の電圧とを往復して、X線源112へと適用される電圧を変化させる、よう構成された制御部の一例である。この制御部は、例えば、X線源112へと適用された電圧が一連の投影ビューの二つ以上の最も近い投影ビューの間に変化しなかった後にのみ、X線源112へと適用される電圧を変化させることで、疎なビューkV切替を実行するよう構成された、としてもよい。データ収集システム176は、集中型システム又は代替的に分散型システムとすることができる。一実施形態において、データ収集システム176は、プロセッサ170に統合されている。この場合、プロセッサ170が制御部の一例となる。プロセッサ170は、投影データからの画像の再構成、投影データの前再構成処理、そして画像データの後再構成処理を含む機能を実行する。プロセッサ170は、ここに説明される機能及び方法も実行する。ガントリ140及びデータ収集システム176は、取得部の一例である。プロセッサ170は、処理部の一例である。処理部は、CTスキャナのプロセッサ176に限らず、取得部にネットワークを介して接続されたサーバ装置として実装してもよい。言い換えると、X線CTシステムは、単体の装置が取得部及び処理部を備えた構成で実現してもよく、取得部を備えた第1装置と、処理部を備えた第2装置とがネットワークを介して接続された構成で実現してもよい。このことは、全ての実施形態に共通している。
投影データの前再構成処理は、検出器キャリブレーション、検出器非線形性、極性効果、ノイズバランシング、そして物質弁別に対する補正を含むことができる。加えて、当該前再構成処理は、ステップ210における様々な処理を含むこともできる。
後再構成処理は、画像のフィルタリングやスムージング、ボリュームレンダリング処理、画像差分処理を必要に応じて含むことができる。加えて、後再構成処理は、処理202、310、320、330を含む、様々な方法200からのステップを含むことができる。
画像再構成処理は、フィルタ補正逆投影法、逐次画像再構成方法、又は確率論的画像再構成法を使用して、実行することができる。更に、画像再構成処理は、ステップ220を含むことがある。
プロセッサ170とデータ収集システム176との両方は、例えば投影データ、再構成された画像、キャリブレーションデータやパラメータ、そしてコンピュータプログラムを格納するために、メモリ178を利用できる。
プロセッサ170は、CPU及びネットワークコントローラを含むことができる。係るCPUは、特定用途向け集積回路(ASIC)、フィールドプログラマブルゲートアレイ(FPGA)、又はその他の複合プログラマブル論理デバイス(CPLD)など、個々の論理ゲートとして実行することができる。FPGA又はCPLDの実施形態は、VHDL、ベリログ、又は任意のその他のハードウェア記述言語でコード化することができ、そして係るコードは、FPGA又はCPLDにおいて直接電子メモリ内に格納することができるし、又は別箇の電子メモリとして格納することもできる。更に、メモリは、ROM、EPROM、EEPROM(登録商標)、又はFLASHメモリ等、不揮発性であってもよい。係るメモリは、静的又は動的RAM等揮発性でもよく、電子メモリの他にもFPGA又はCPLDとメモリとの相互作用を管理するためのマイクロコントローラやマイクロプロセッサ等プロセッサが提供されていてもよい。
代替的に、CPUは、ここに説明された機能を実行するコンピュータ読み取り可能命令のセットを含んでいるコンピュータプログラムを実行することができ、係るコンピュータプログラムは、任意の上述の非一時的電子メモリ及び/又はハードディスクドライブ、CD、DVD、FLASHドライブ、又はその他の任意の公知の格納媒体に格納されている。更に、コンピュータ読み取り可能命令は、ユーティリティアプリケーション、バックグラウンドデーモン、又はオペレーティングシステムのコンポーネント、又はそれらの組み合わせとして提供されてもよく、米国Intel社からのXenonプロセッサ又は米国AMD社からのOpteronプロセッサ等のプロセッサと、そしてMicrosoft VISTA、UNIX(登録商標)、Solaris、LINUX(登録商標)、Apple、MAC-OS等のオペレーティングシステムや、当業者にとって公知のその他オペレーティングシステムがプロセッサと一体となって実行する。更に、CPUは、命令を実行するために並行して協同的に動作する、マルチプルプロセッサとして実行することができる。
一実施形態において、ディスプレイ上に再構成された画像を映し出すことがある。当該ディスプレイは、LCDディスプレイ、CRTディスプレイ、プラズマディスプレイ、OLED、LED、又は当業者にとって既知の、任意のその他のディスプレイとすることができる。ネットワークコントローラは、例えば米国Intel社からのインテルイーサネット(登録商標)PROネットワークインターフェースカードとすることがあり、CTスキャナの様々な部分間とインターフェースを取ることがある。加えて、ネットワークコントローラは、外部ネットワークとインターフェースを取ることもある。理解されている通り、係る外部ネットワークは、インターネット等公衆ネットワークや、LAN又はWANネットワークなど私的ネットワークや、これらの任意の組み合わせでも良く、PSTN又はISDNサブネットワークを含んでも良い。外的ネットワークは、イーサネットワークのように有線で接続されていても良いし、又はEDGEや3G、4G等の無線セルラーシステムを含むセルラーネットワークのような無線でも良い。また無線ネットワークは、WiFi、Bluetooth(登録商標)、又は任意のその他の公知の通信の無線方式であっても良い。
メモリ178は、ハードディスクドライブ、CD-ROMドライブ、DVDドライブ、FLASHドライブ、RAM、ROM、又は当業者にとって既知の任意のその他の格納メディアであってもよい。
図16は、CT装置又はCTスキャナ100に含まれる放射線ガントリの第二の実施形態を描いている。図16に図示されるように、放射線ガントリ1000は側面から見て描かれており、X線管1001、環状フレーム1002、そして多列又は二次元アレイ型X線検出器1003を更に含む。X線管1001及びX線検出器1003は、環状フレーム1002上に被検体OBJを横切って正反対に取り付けられ、環状フレーム1002は回転軸RAの回りに回転可能に支持される。被検体OBJが図示された頁の奥の方向又は手前の方向の軸RAに沿って移動されながら、回転ユニット1007は環状フレーム1002を0.4秒/回転もの高速で回転させる。
本発明に係るX線コンピュータ断層撮影(CT)装置の第一の実施形態は、付随する図面を参照しながら以下に説明される。X線CT装置は、様々なタイプの装置を含むことに留意されたい。具体的には、X線管とX線検出器とが検査される予定の被検体の周辺を一緒に回る回転/回転型装置と、そして多数の検出器素子がリング状又は水平状に配置されており、X線管のみが検査される予定の被検体の周辺を回る固定/回転型装置とがある。本発明は、いずれのタイプにも適用可能である。今回は、回転/回転型タイプが例示される。
マルチスライスX線CT装置は高電圧発生器1009を更に含み、X線管1001がX線を生成するように、スリップリング1008を通して、高電圧発生器1009はX線管1001に印加される管電圧を生成する。X線は、被検体OBJに向かって放射され、被検体OBJの断面領域が円で表される。例えば、第一のスキャンにわたる平均的なX線エネルギーを有するX線管1001は、第二のスキャンにわたる平均的なX線エネルギーよりも小さい。このようにして、異なるX線エネルギーに対応して二回以上のスキャンを取得することができる。X線検出器1003は、被検体OBJを通り抜けてきた、照射されたX線を検出するために、被検体OBJを挟んでX線管1001から反対側の位置にある。X線検出器1003は、個々の検出器素子又は検出器ユニットを更に含む。
CT装置は、X線検出器1003から検出された信号を処理するための、その他のデバイスを更に含む。データ収集回路又はデータ収集システム(DAS)1004は、各チャネルに対するX線検出器1003からの出力信号を電圧信号に変換し、係る信号を増幅し、更にその信号をデジタル信号へと変換する。X線検出器1003及びDAS1004は、1回転当たりの所定全投影数(TPPR)を処理するよう構成されている。
上記で説明されたデータは、非接触データ送信装置1005を通して、放射線ガントリ1000外部のコンソール内に収容された、前処理回路1006へと送られる。前処理回路1006は、生データに対し感度補正など、特定の補正を実行する。ストレージ1012は、再構成処理直前のステージで、「投影データ」とも呼ばれる結果データを格納する。ストレージ1012は、再構成回路1014、入力インターフェース1015、ディスプレイ1016と共に、データ/制御バス1011を通して、処理回路1010に接続されている。処理回路1010は、CTシステムを駆動させるのに十分なレベルに達するまで電流を制限する電流調整器1013を制御する。「ストレージ」は、「メモリ」、「記憶部」又は「格納部」等と読み替えてもよい。「再構成回路」は、「再構成デバイス」又は「再構成処理部」等と読み替えてもよい。「ディスプレイ」は、「表示部」等と読み替えてもよい。放射線ガントリ1000は、取得部の一例である。再構成回路1014は、処理部の一例である。処理部は、X線CT装置の再構成回路1014に限らず、取得部にネットワークを介して接続されたサーバ装置として実装してもよい。
検出器は、どんな世代のCTスキャナシステムであっても、患者に対して回転及び/又は固定される。一実施形態において、上記のCTシステムは、第三世代ジオメトリシステムと第四世代ジオメトリシステムとが組み合わせられた例とする場合がある。第三世代ジオメトリシステムにおいて、X線管1001とX線検出器1003とは、環状フレーム1002上に正反対に取り付けられ、環状フレーム1002が回転軸RAの周りを回転する時に、被検体OBJの周りを回転する。第四世代ジオメトリシステムにおいて、検出器は患者の周辺に固定して取り付けられており、X線管は患者の周辺を回転する。代替的な実施形態において、放射線ガントリ1000は、Cアーム及びスタンドによって支持されている、環状フレーム1002上に配置された多数の検出器を有する。
ストレージ1012は、X線検出器1003でX線照射量を示す測定値を格納することができる。更に、ストレージ1012は、ここで説明された方法(例えば、方法200及びその変形例)を実行するための専用プログラムを格納することができる。
再構成回路1014は、ここに説明された方法(例えば、方法200のステップ220及びその変形例)の様々なステップを実行することができる。例えば、再構成回路1014は、前述した逐次近似応用被曝低減(AIDR)三次元(3D)再構成法を使用してもよい。この再構成法を使用する場合、再構成回路1014は、ノイズ低減、アーチファクト低減を目的として投影データ上と画像データ上との二段階でノイズモデルを用いたノイズ低減処理を実行し、S/N比及び質感の高い再構成画像を生成する。この場合、再構成回路1014は、図17に示すように、第1のノイズ取得部1014a、第2のノイズ取得部1014b、ノイズ推定部1014c、ノイズ低減処理部1014d、再構成画像生成部1014e、最適化処理部1014f、更新処理部1014g、ブレンド処理部1014hを有している。なお、これらの「・・・部」は、「・・・回路」等と読み替えてもよい。
第1のノイズ取得部1014aは、対数変換前の投影カウントデータに対して、システムモデル(装置ジオメトリ等の装置固有の特性を考慮したスキャナーモデル等)を利用してノイズ量を推定する。
第2のノイズ取得部1014bは、対数変換前の投影カウントデータに対して統計学的ノイズモデル(各X線量帯におけるフォトンノイズと電気的ノイズとを考慮した統計学的なノイズモデル)を利用してノイズ量を推定する。
ノイズ推定部1014cは、第1のノイズ取得部1014a、第2のノイズ取得部1014bのそれぞれにおいて推定された各ノイズ量に基づいて、対数変換前の投影カウントデータについての推定ノイズの総量を取得する。
ノイズ低減処理部1014dは、ノイズ推定部1014cにおいて推定されたノイズ量に従って、対数変換前の投影カウントデータに対しノイズ低減処理を実行する。また、ノイズ低減処理部1014dは、ノイズ低減された対数変換前の投影カウントデータに対数変換等を実行し出力する。
再構成画像生成部1014eは、ノイズ低減処理部1014dより受け取った(ノイズ低減処理及び対数変換後の)投影データを用いた再構成処理を実行し、再構成画像を生成する。この再構成画像生成部1014eが実行する再構成処理は、逐次近似再構成処理(フルIR処理)、フィルタードバックプロジェクション処理(FBP処理)、IR処理とFBP処理の双方を用いた処理(ハイブリッド処理)のいずれであってもよい。
最適化処理部1014fは、解剖学的モデルを用いて、再構成画像生成部1014eより受け取った再構成画像又は更新処理部1014gからフィードバックされる更新処理後の再構成画像に対し、三次元情報を用いて鮮鋭度の高い構造を維持すると同時に、画像データ上に存在するノイズを選択的に除去する最適化処理を実行する。また、最適化処理部1014fは、再構成画像生成部1014eより受け取った再構成画像に対し、所定の平滑化処理を実行する。
更新処理部1014gは、最適化処理部1014fから出力される、最適化処理後の画像データと平滑化処理後の画像データとを用いて、再構成画像の更新処理を実行する。更新後の再構成画像は再び最適化処理部1014fへと出力される。更新処理部1014gと最適化処理部1014fとの間で、所定の条件を満たすまで、上述した最適化処理、平滑化処理、更新処理が繰り返し実行される。
ブレンド処理部1014hは、更新処理部1014gから出力された再構成画像と、再構成画像生成部1014eから出力された(初期)再構成画像とを所定の比率でブレントし、S/N比を向上させつつ画像データ上でノイズ粒状性を維持し、より自然な質感を有する再構成画像を生成する。
なお、更新処理部1014gと最適化処理部1014fとの間における繰り返し処理の回数、ブレンド処理部1014hにおけるブレンド処理のブレンド比率は、撮像部位や臨床用途等を基準として決定することができる。
更に、再構成回路1014は、必要に応じてボリュームレンダリング処理や画像差分処理等、前再構成処理画像処理を実行することができる。
前処理回路1006によって実行された投影データの前再構成処理は、検出器キャリブレーション、検出器非直線性、極性効果のための補正を例えば含むことができる。更に、前再構成処理は、ステップ210を含むこともできる。
再構成回路1014によって実行される後再構成処理は、画像のフィルタリングやスムージング、ボリュームレンダリング処理、そして画像差分処理を、必要に応じて含むことができる。画像再構成処理は、方法200の様々なステップ(例えば、ステップ230、240)を実行でき、またDL-ANNネットワークのオフライン学習(例えば、処理310、320、330)も実行することができる。再構成回路1014は、例えば投影データ、再構成された画像、キャリブレーションデータやパラメータ、そしてコンピュータプログラムを格納するのに、ストレージを使うことがある。
再構成回路1014は、特定用途向け集積回路(ASIC)、フィールドプログラマブルゲートアレイ(FPGA)、又は複合プログラマブル論理デバイス(CPLD)など、個々の論理ゲートとして実行可能な、CPU(中央演算処理装置)を含むことができる。FPGA又はCPLDの実行は、VHDL、ベリログ、又は任意のその他のハードウェア記述言語でコード化することができ、そして係るコードはFPGA又はCPLDにおいて直接電子メモリ内に格納することができるし、又は別箇の電子メモリとして格納することもできる。更に、ストレージ1012は、ROM、EPROM、EEPROM(登録商標)、又はFLASHメモリ等、不揮発性メモリであってもよい。ストレージ1012は、静的又は動的RAM等揮発性でもよく、電子メモリだけでなくFPGA又はCPLDとメモリとの間の相互作用を管理するマイクロコントローラやマイクロプロセッサ等の処理部が設けられることもある。
代替的に、再構成回路1014におけるCPUは、ここで説明された機能を実行するコンピュータ読み取り可能命令のセットを含むコンピュータプログラムを実行することができ、係るコンピュータプログラムは、任意の上述の非一時的電子メモリ及び/又はハードディスクドライブ、CD、DVD、FLASHドライブ、又はその他の任意の既知の格納媒体に格納されている。更に、コンピュータ読み取り可能命令は、ユーティリティアプリケーション、バックグラウンドデーモン、又はオペレーティングシステムのコンポーネント、又はそれらの組み合わせで提供されてもよく、米国Intel社からのXenonプロセッサ又は米国AMD社からのOpteronプロセッサなどのプロセッサと、そしてMicrosoft VISTA、UNIX、Solaris、LINUX、Apple、MAC-OS等のオペレーティングシステムや、当業者にとっては既知のその他のオペレーティングシステムと一体となって実行する。更に、CPUは、命令を実行するために並行して協同的に動く、マルチプルプロセッサとして実行されてもよい。
一実施形態において、ディスプレイ1016上に再構成された画像を映し出すことができる。係るディスプレイ1016は、LCDディスプレイ、CRTディスプレイ、プラズマディスプレイ、OLED、LED、又は当業者にとって既知のその他のディスプレイであってもよい。
ストレージ1012は、ハードディスクドライブ、CD-ROMドライブ、DVDドライブ、FLASHドライブ、RAM、ROM、又は当業者にとって既知のその他の格納メディアであってもよい。
特定の実施形態が述べられてきた一方で、これらの実施形態は、一例として提示されたに過ぎず、本開示の範囲を限定することを意図したものではない。実際には、ここで説明された新たな方法、装置やシステムは、その他の様々な形態で具体化することが可能であり、更には、本開示の趣旨から乖離することなく、本開示に説明された方法、装置やシステムの形式で省略、置き換え、変更が可能である。
100 CTスキャナ
103 X線検出器
110 環状回転フレーム
112 X線源
114 コリメータ/フィルタ
140 ガントリ
170 プロセッサ
176 データ収集システム
178 メモリ
1000 放射線ガントリ
1001 X線管
1002 環状フレーム
1003 X線検出器
1004 データ収集回路(DAS)
1005 非接触データ送信機
1006 前処理回路
1007 回転ユニット
1008 スリップリング
1009 高電圧発生器
1010 処理回路
1011 データ/制御バス
1012 ストレージ
1013 電流調整器
1014 再構成回路
1014a 第1のノイズ取得部
1014b 第2のノイズ取得部
1014c ノイズ推定部
1014d ノイズ低減処理部
1014e 再構成画像生成部
1014f 最適化処理部
1014g 更新処理部
1014h ブレンド処理部
1015 入力インターフェース
1016 ディスプレイ

Claims (25)

  1. 離散的な第1のビュー群に対応し、第1のエネルギーのX線の照射を伴って収集された第1の投影データセットと、前記第1のビュー群とは異なる離散的な第2のビュー群に対応し、前記第1のエネルギーとは異なる第2のエネルギーのX線の照射を伴って収集された第2の投影データセットと、を取得する取得部と、
    前記第1の投影データセット及び前記第2の投影データセットに基づいて、第3のビュー群に対応し、前記第1のエネルギーのX線に対応する第3の投影データセット、及び第4のビュー群に対応し、前記第2のエネルギーのX線に対応する第4の投影データセットを生成する学習済みモデルに対して、前記第1の投影データセット及び前記第2の投影データセットを入力することで、前記第3の投影データセット及び前記第4の投影データセットを生成する処理部と、
    を備え、
    前記処理部は、
    複数の検出器素子のX線放射の強度を表す投影データであって、前記第1の投影データセット及び前記第2の投影データセットを構成する前記投影データを前記取得部から取得し、前記投影データは、投影ビューの関数として低kV投影データと交互配置された高kV投影データを表しており、前記投影データにおいて前記高kV投影データの投影ビューは前記低kV投影データの投影ビューから分離されており、前記低kV投影データはX線源へと適用された第一の電圧を使用して収集されており、前記高kV投影データは前記X線源へと、前記第一の電圧に比べてより大きい第二の電圧を使用して収集されており、
    高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第一のニューラルネットワークである前記学習済みモデルを取得し、前記第一のニューラルネットワークは、入力データと対象データとを含む学習データセットを使用して学習されており、前記入力データはkV切替サイノグラムに基づいており、前記入力データにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されており、前記対象データは完全な低kVサイノグラムと完全な高kVサイノグラムとを含む完全なサイノグラムに基づいており、前記完全なサイノグラムのそれぞれが全ての投影ビューに対する投影データを含んでいる、
    前記高kV投影データと前記低kV投影データとに対するサイノグラム補填を実行し、
    前記第一のニューラルネットワークに対して、前記低kV投影データに基づいた低kV入力と前記高kV投影データに基づいた高kV入力とを適用し、前記第一のニューラルネットワークは、前記サイノグラム補填からのエラーが低減された結果を出力している、
    よう構成されているX線CTシステム。
  2. 前記第3のビュー群及び前記第4のビュー群は、それぞれ前記第1のビュー群及び前記第2のビュー群を含む、請求項1に記載のX線CTシステム。
  3. 前記第3のビュー群は、前記第1のビュー群より多く、前記第4のビュー群は、前記第2のビュー群より多い、請求項1に記載のX線CTシステム。
  4. 前記第1の投影データセット、前記第2の投影データセット、前記第3の投影データセット、及び前記第4の投影データセットは、それぞれサイノグラムである、請求項1に記載のX線CTシステム。
  5. 前記第3の投影データセット及び前記第4の投影データセットは、それぞれ画像再構成に必要な全てのビューに対応するサイノグラムである、請求項1に記載のX線CTシステム。
  6. 前記処理部は、前記低kV入力と前記高kV入力とを前記第一のニューラルネットワークへと適用することと一体的に、前記サイノグラム補填を実行するよう更に構成され、
    前記第一のニューラルネットワークへの前記高kV入力は、前記サイノグラム補填よりも前の前記高kV投影データであり、
    前記第一のニューラルネットワークへの前記低kV入力は、前記サイノグラム補填よりも前の前記低kV投影データであり、
    前記第一のニューラルネットワークは、サイノグラム補填を実行するために、前記高kV投影データと前記低kV投影データとの相補的な情報を使用するよう学習され、全ての投影ビューに対する投影データを含む完全な高kV投影データと、全ての投影ビューに対する投影データを含む完全な低kV投影データとを前記第3の投影データセット及び前記第4の投影データセットとして生成している、請求項記載のX線CTシステム。
  7. 前記処理部は、前記第一のニューラルネットワークに対して、前記低kV入力と前記高kV入力との適用よりも前に前記サイノグラム補填を実行するよう更に構成され、
    全ての投影ビューに対する投影データを含む完全な高kV投影データと全ての投影ビューに対する投影データを含む完全な低kV投影データとを生成する前記サイノグラム補填を実行し、
    前記高kV入力としての前記完全な高kV投影データと前記低kV入力としての前記完全な低kV投影データとを前記第一のニューラルネットワークに適用して、前記第一のニューラルネットワークから出力される前記完全な高kV投影データ及び前記完全な低kV投影データ内でサイノグラム補填エラーが低減された結果を出力する、請求項記載のX線CTシステム。
  8. 前記処理部は、
    前記高kV投影データと前記低kV投影データとを第一の基底投影データと第二の基底投影データとに弁別し、
    前記第一の基底投影データから第一の基底画像を再構成し、
    前記第二の基底投影データから第二の基底画像を再構成するよう更に構成された、請求項記載のX線CTシステム。
  9. 前記処理部は、
    前記高kV投影データと前記低kV投影データとを、物質弁別を使用して、第一の物質成分投影データと第二の物質成分投影データとに弁別するよう更に構成された、請求項記載のX線CTシステム。
  10. 前記処理部は、
    前記第一の基底投影データ用の一方の入力チャネルと、前記第二の基底投影データ用の他方の入力チャネルとの二つのチャネルを有する第二のニューラルネットワークを取得し、
    前記第二のニューラルネットワークは、前記kV切替サイノグラムの基底弁別に基づいている入力データと、前記完全なサイノグラムの基底弁別に基づいている対象データと、を含む学習データセットを使用して学習されている、請求項記載のX線CTシステム。
  11. 前記処理部は、
    前記第一の基底画像用の一方の入力チャネルと、前記第二の基底画像用の他方の入力チャネルとの二つのチャネルを有する第二のニューラルネットワークを取得し、
    前記第二のニューラルネットワークは、入力データと対象データとを含む学習データセットを使用して学習されており、
    前記入力データは、前記kV切替サイノグラムの基底弁別から再構成された画像に基づいており、
    前記対象データは、前記完全なサイノグラムの基底弁別から再構成された画像に基づいているよう更に構成された、請求項記載のX線CTシステム。
  12. 前記処理部は、
    高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第二のニューラルネットワークを取得し、
    前記低kV投影データと前記高kV投影データとを第一の面に沿った個別の第一のスライスへと区切り、
    部分的に補正された低kV投影データと部分的に補正された高kV投影データとを生成するために、前記個別の第一のスライスを前記第一のニューラルネットワークへと適用し、
    前記部分的に補正された低kV投影データと前記部分的に補正された高kV投影データとを前記第一の面と直交である第二の面に沿った個別の第二のスライスへと区切り、
    補正された低kV投影データと補正された高kV投影データとを生成するために、前記個別の第二のスライスを前記第二のニューラルネットワークへと適用する、よう更に構成された、請求項記載のX線CTシステム。
  13. 前記処理部は、
    高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第二のニューラルネットワークを取得し、
    前記低kV投影データを、高空間周波数低kV投影データと低空間周波数低kV投影データとを含む低kV周波数成分へと変換し、
    前記高kV投影データを、高空間周波数高kV投影データと低空間周波数高kV投影データとを含む高kV周波数成分へと変換し、
    前記第一のニューラルネットワークに対して、前記高kV周波数成分の高周波数成分と前記低kV周波数成分の高周波数成分とを適用し、
    前記第二のニューラルネットワークに対して、前記高kV周波数成分の低周波数成分と前記低kV周波数成分の低周波数成分とを適用し、
    補正された低kV投影データと補正された高kV投影データとを生成するために、個別の前記高周波数成分を、対応する低周波数成分と組み合わせるよう更に構成された、請求項記載のX線CTシステム。
  14. 前記処理部は、
    低空間周波数投影データを生成するために、前記投影データに対してローパスフィルタを適用し、且つ高空間周波数投影データを生成するために、前記低空間周波数投影データを前記投影データから引くことと、
    前記投影データを低空間周波数投影データと高空間周波数投影データとに分けるために、ウェーブレット変換を前記投影データへと適用することと、
    前記投影データを周波数ドメインへと変換し、且つ前記周波数ドメインにおける前記投影データを低空間周波数投影データと高空間周波数投影データとに分けることと、
    のうちの一つにより前記投影データを周波数成分へと変換するよう更に構成された、請求項13記載のX線CTシステム。
  15. 前記処理部は、
    ダウンサンプルされた投影データを生成するために、前記低空間周波数低kV投影データと前記低空間周波数高kV投影データとをダウンサンプリングすることと、
    ダウンサンプルされた結果を生成するために、前記ダウンサンプルされた投影データを前記第二のニューラルネットワークへと適用することと、
    補正された低空間周波数低kV投影データと補正された低空間周波数高kV投影データとを生成するために、前記ダウンサンプルされた結果をアップサンプルすることと、
    により前記低周波数成分を前記第二のニューラルネットワークへと適用するよう更に構成された、請求項13記載のX線CTシステム。
  16. 前記処理部は、
    高kV投影データが収集された投影ビューに対して低kV投影データが収集された投影ビューを表すkV切替マスクを取得するよう更に構成され、
    前記第一のニューラルネットワークは前記kV切替マスクを受信するための第三の入力チャネルを有する、請求項記載のX線CTシステム。
  17. X線源へと適用される電圧に依存するエネルギースペクトルを有するX線を照射するよう構成された前記X線源と、
    第一の電圧が前記X線源へと適用される場合に、前記高kV投影データを出力し、且つ前記第一の電圧に比べて小さい第二の電圧が前記X線源へと適用された場合に、前記低kV投影データを出力するよう構成された複数の検出器素子を有するエネルギー積分型X線検出器と、
    一連の投影ビューを通して、前記X線源と前記エネルギー積分型X線検出器とを回転させ、かつ
    前記第一の電圧と前記第二の電圧とを往復して、前記X線源へと適用される電圧を変化させる、
    よう構成された制御部と、
    を備え、
    前記制御部は、
    前記X線源へと適用された前記電圧が前記一連の投影ビューの二つ以上の最も近い投影ビューの間に変化しなかった後にのみ、前記X線源へと適用される前記電圧を変化させることで、疎なビューkV切替を実行するよう構成された、請求項記載のX線CTシステム。
  18. 複数の検出器素子でX線照射の強度を表す投影データを取得し、前記投影データは、投影ビューの関数として低kV投影データと交互配置された高kV投影データを表しており、前記投影データにおいて前記高kV投影データの投影ビューは前記低kV投影データの投影ビューから分離されており、前記低kV投影データはX線源へと適用された第一の電圧を使用して収集されており、前記高kV投影データは前記X線源へと、前記第一の電圧に比べてより大きい第二の電圧を使用して収集されており、
    高kVデータ用の第一のチャネルを有する第一のニューラルネットワークを取得し、前記第一のニューラルネットワークは、入力データと対象データとを含む第一の学習データセットを使用して学習されており、前記入力データはkV切替サイノグラムに基づいており、前記入力データにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されており、前記対象データは、完全なサイノグラムのそれぞれが全ての投影ビューに対する投影データを含む、完全な高kVサイノグラムを含む前記完全なサイノグラムに基づいている、
    低kVデータ用の第一のチャネルを有する第二のニューラルネットワークを取得し、前記第二のニューラルネットワークは、入力データと対象データとを含む第二の学習データセットを使用して学習されており、前記入力データはkV切替サイノグラムに基づいており、前記入力データにおいて前記kV切替低kVサイノグラムの前記投影ビューが前記kV切替高kVサイノグラムの前記投影ビューから分離されており、前記対象データは、全ての投影ビューに対する投影データを含んでいる完全な低kVサイノグラムを含む、完全なサイノグラムに基づいている、
    前記高kV投影データに対するサイノグラム補填と、前記低kV投影データに対するサイノグラム補填とを実行し、
    前記第一のニューラルネットワークに対して、前記高kV投影データに基づいた高kV入力を適用し、前記第一のニューラルネットワークは、前記高kV投影データの前記サイノグラム補填からのエラーが低減された結果を出力している、
    前記第二のニューラルネットワークに対して、前記低kV投影データに基づいた低kV入力を適用し、前記第二のニューラルネットワークは、前記低kV投影データの前記サイノグラム補填からのエラーが低減された結果を出力している、
    よう構成された処理部を具備するX線CTシステム。
  19. 前記処理部は、
    前記第一のニューラルネットワークを取得し、前記第一のニューラルネットワークは、前記kV切替低kVサイノグラムが収集された投影ビューに対して前記kV切替高kVサイノグラムが収集された投影ビューを表すkV切替マスクに対する第二のチャネルを有し、
    前記第二のニューラルネットワークを取得し、前記第二のニューラルネットワークは、前記kV切替マスク用の第二のチャネルを有し、
    前記第一のニューラルネットワークに対して、前記取得された投影データの高kV入力と前記kV切替マスクとを適用し、
    前記第二のニューラルネットワークに対して、前記取得された投影データの低kV入力と前記kV切替マスクとを適用し、
    前記第一のニューラルネットワークは、前記kV切替サイノグラムに対応する個別のkV切替マスクを含む前記入力データを含む前記第一の学習データセットを使用して学習されており、
    前記第二のニューラルネットワークは、前記kV切替サイノグラムに対応する個別のkV切替マスクを含む前記入力データを含む前記第二の学習データセットを使用して学習されている、
    よう更に構成された、請求項18記載のX線CTシステム。
  20. 前記処理部は、
    前記投影データを取得し、前記投影データは、前記投影ビューの前記関数として交互配置された前記高kV投影データと前記低kV投影データと中kV投影データを表す、前記中kV投影データは、前記第一の電圧と前記第二の電圧との間の切り替え途中の電圧を使用して収集されており、
    前記高kVデータ用の前記一方の入力チャネルと、前記低kVデータ用の前記他方の入力チャネルと、中kVデータ用の第三の入力チャネルとの三つの入力チャネルを有する前記第一のニューラルネットワークを取得する、よう更に構成された、請求項記載のX線CTシステム。
  21. サイノグラム補填ニューラルネットワークのネットワーク係数を初期化することと、kV切替サイノグラムに基づいた高kV入力と低kV入力とのペアを含む入力データを取得することと、前記kV切替サイノグラムにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されており、前記kV切替低kVサイノグラムは、投影ビューの関数として前記kV切替高kVサイノグラムと交互配置される、
    完全な低kVサイノグラムと完全な高kVサイノグラムとを含む完全なサイノグラムに基づいた高kV対象と低kV対象とのペアを含む対象データを取得することと、前記完全なサイノグラムのそれぞれは全ての投影ビューに対する投影データを含む、
    前記対象データと、前記入力データを前記サイノグラム補填ニューラルネットワークへ適用することからの出力との間の一致を表す損失関数を最適化するために、前記ネットワーク係数を逐次的にアップデートすることと、
    により、二つのチャネルを有する前記サイノグラム補填ニューラルネットワークを学習させるよう構成された処理部を具備するX線CTシステム。
  22. 前記処理部は、
    前記kV切替高kVサイノグラムが収集された投影ビューに対して前記kV切替低kVサイノグラムが収集された投影データを表すkV切替マスクを取得することと、
    前記損失関数を最適化するために、前記ネットワーク係数を逐次的にアップデートすることと、
    によりアーチファクト低減ニューラルネットワークを学習させるよう更に構成された、請求項21記載のX線CTシステム。
  23. 前記処理部は、前記アーチファクト低減ニューラルネットワークを学習させるよう更に構成され、前記kV切替マスクは、前記kV切替サイノグラムにおける欠損値に対して、前記kV切替サイノグラムにおいて最も近い収集された値に対する近似を更に表す請求項22記載のX線CTシステム。
  24. kV切替を用いて収集された投影データのサイノグラム補填の方法であって、
    複数の検出器素子のX線放射の強度を表す投影データを取得することと、前記投影データは、投影ビューの関数として低kV投影データと交互配置された高kV投影データを表しており、前記投影データにおいて前記高kV投影データの投影ビューは前記低kV投影データの投影ビューから分離されており、前記低kV投影データはX線源へと適用された第一の電圧を使用して収集されており、前記高kV投影データは前記X線源へと、前記第一の電圧に比べてより大きい第二の電圧を使用して収集されており、
    高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第一のニューラルネットワークを取得することと、前記第一のニューラルネットワークは、入力データと対象データとを含む学習データセットを使用して学習されており、前記入力データはkV切替サイノグラムに基づいており、前記入力データにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されており、前記対象データは完全な低kVサイノグラムと完全な高kVサイノグラムとを含む完全なサイノグラムに基づいており、前記完全なサイノグラムのそれぞれが全ての投影ビューに対する投影データを含んでいる、
    前記高kV投影データと前記低kV投影データとに対するサイノグラム補填を実行することと、
    前記第一のニューラルネットワークに対して、前記低kV投影データに基づいた低kV入力と前記高kV投影データに基づいた高kV入力とを適用することと、前記第一のニューラルネットワークは、前記サイノグラム補填からのエラーが低減された結果を出力する、
    を具備する方法。
  25. 取得部と、
    処理部と、
    を備え、
    前記取得部は、
    複数の検出器素子のX線放射の強度を表す投影データを取得し、前記投影データは、投影ビューの関数として低kV投影データと交互配置された高kV投影データを表しており、前記投影データにおいて前記高kV投影データの投影ビューは前記低kV投影データの投影ビューから分離されており、前記低kV投影データはX線源へと適用された第一の電圧を使用して収集されており、前記高kV投影データは前記X線源へと、前記第一の電圧に比べてより大きい第二の電圧を使用して収集されており、
    前記処理部は、
    高kVデータ用の一方の入力チャネルと、低kVデータ用の他方の入力チャネルとの二つの入力チャネルを有する第一のニューラルネットワークを取得し、前記第一のニューラルネットワークは、入力データと対象データとを含む学習データセットを使用して学習されており、前記入力データはkV切替サイノグラムに基づいており、前記入力データにおいてkV切替低kVサイノグラムの投影ビューがkV切替高kVサイノグラムの投影ビューから分離されており、前記対象データは完全な低kVサイノグラムと完全な高kVサイノグラムとを含む完全なサイノグラムに基づいており、前記完全なサイノグラムのそれぞれが全ての投影ビューに対する投影データを含んでいる、
    前記高kV投影データと前記低kV投影データとに対するサイノグラム補填を実行し、
    前記第一のニューラルネットワークに対して、前記低kV投影データに基づいた低kV入力と前記高kV投影データに基づいた高kV入力とを適用し、前記第一のニューラルネットワークは、前記サイノグラム補填からのエラーが低減された結果を出力している、よう構成されている、X線CTシステム。
JP2019053607A 2018-12-20 2019-03-20 X線ctシステム及び方法 Active JP7334048B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP19216170.1A EP3671646A1 (en) 2018-12-20 2019-12-13 X-ray computed tomography (ct) system and method
CN201911346233.6A CN111345835B (zh) 2018-12-20 2019-12-20 X射线ct系统及方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US16/227,251 US11039806B2 (en) 2018-12-20 2018-12-20 Apparatus and method that uses deep learning to correct computed tomography (CT) with sinogram completion of projection data
US16/227,251 2018-12-20

Publications (2)

Publication Number Publication Date
JP2020099662A JP2020099662A (ja) 2020-07-02
JP7334048B2 true JP7334048B2 (ja) 2023-08-28

Family

ID=71099395

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019053607A Active JP7334048B2 (ja) 2018-12-20 2019-03-20 X線ctシステム及び方法

Country Status (3)

Country Link
US (2) US11039806B2 (ja)
JP (1) JP7334048B2 (ja)
CN (1) CN111345835B (ja)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019134905A1 (en) * 2018-01-02 2019-07-11 Koninklijke Philips N.V. Image reconstruction employing tailored edge preserving regularization
US11039806B2 (en) * 2018-12-20 2021-06-22 Canon Medical Systems Corporation Apparatus and method that uses deep learning to correct computed tomography (CT) with sinogram completion of projection data
US11176428B2 (en) * 2019-04-01 2021-11-16 Canon Medical Systems Corporation Apparatus and method for sinogram restoration in computed tomography (CT) using adaptive filtering with deep learning (DL)
CN117234356A (zh) * 2019-05-08 2023-12-15 义隆电子股份有限公司 用于触控装置的控制系统及方法
KR20210017185A (ko) * 2019-08-07 2021-02-17 한국전자통신연구원 심층 신경망을 기반으로 영상의 압축 포아송 잡음을 제거하는 방법 및 장치
US11455755B2 (en) * 2019-09-05 2022-09-27 Siemens Medical Solutions Usa, Inc. Methods and apparatus for neural network based image reconstruction
US11350895B2 (en) * 2019-11-29 2022-06-07 Wisconsin Alumni Research Foundation System and method for spectral computed tomography using single polychromatic x-ray spectrum acquisition
CN112051285B (zh) * 2020-08-18 2021-08-20 大连理工大学 X射线实时成像与ct断层扫描一体化智能无损检测系统
US20230326100A1 (en) * 2020-08-19 2023-10-12 Prismatic Sensors Ab Methods and systems related to x-ray imaging
EP3964131A1 (en) * 2020-09-03 2022-03-09 Koninklijke Philips N.V. Spectral x-ray material decomposition method
CN112651434B (zh) * 2020-12-21 2022-09-02 明峰医疗系统股份有限公司 一种基于深度学习的能谱ct基物质投影估计方法、系统及能谱ct
EP4057229A1 (de) * 2021-03-12 2022-09-14 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Verarbeitungsvorrichtung zum erhalten einer lückenfüller-sinogramminformation, computertomograph, verfahren und computerprogramm
US11823354B2 (en) * 2021-04-08 2023-11-21 GE Precision Healthcare LLC System and method for utilizing a deep learning network to correct for a bad pixel in a computed tomography detector
CN117355865A (zh) * 2021-04-13 2024-01-05 棱镜传感器公司 确定用于计算机断层扫描中的深度学习图像重建的置信度指示
KR102553068B1 (ko) * 2021-04-20 2023-07-06 연세대학교 원주산학협력단 다중 레벨 웨이브렛 인공지능 기반 선량 저 감화를 위한 컴퓨터 단층영상 재구성 방법
DE102021206417B3 (de) 2021-06-22 2022-10-06 Siemens Healthcare Gmbh Computerimplementierte Verfahren und Systeme zum Bereitstellen eines Korrekturalgorithmus für ein Röntgenbild und zum Korrigieren eines Röntgenbilds, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
US20230038970A1 (en) * 2021-07-22 2023-02-09 Leo Cancer Care, Inc. Scanner and method of image reconstruction
US12067651B2 (en) 2021-08-31 2024-08-20 Canon Medical Systems Corporation Projection based deep learning with frequency splitting for computed tomography
WO2023044114A1 (en) * 2021-09-17 2023-03-23 GE Precision Healthcare LLC Systems and methods for energy bin downsampling
KR20230052231A (ko) * 2021-10-12 2023-04-19 주식회사 유케어트론 측정스펙트럼 분석장치 및 분석방법
US11844639B2 (en) * 2021-10-29 2023-12-19 GE Precision Healthcare LLC System and method for utilizing dual energy imaging in a computed tomography imaging system
WO2023137461A1 (en) * 2022-01-13 2023-07-20 The Regents Of The University Of California Neural computed tomography reconstruction
US20240265591A1 (en) * 2023-02-03 2024-08-08 GE Precision Healthcare LLC Systems and methods for interpolation of dual-energy ct data
CN116612206B (zh) * 2023-07-19 2023-09-29 中国海洋大学 一种利用卷积神经网络减少ct扫描时间的方法及系统
CN116720554B (zh) * 2023-08-11 2023-11-14 南京师范大学 一种基于fpga技术的多段线性拟合的神经元电路实现方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004188187A (ja) 2002-11-27 2004-07-08 Ge Medical Systems Global Technology Co Llc アーティファクト低減を容易にする方法及び装置
JP2014503332A (ja) 2011-02-01 2014-02-13 コーニンクレッカ フィリップス エヌ ヴェ デュアルエネルギーct画像再構成のための方法及びシステム
US20160210749A1 (en) 2015-01-21 2016-07-21 Siemens Aktiengesellschaft Method and system for cross-domain synthesis of medical images using contextual deep network
WO2016175755A1 (en) 2015-04-28 2016-11-03 Siemens Healthcare Gmbh METHOD AND SYSTEM FOR SYNTHESIZING VIRTUAL HIGH DOSE OR HIGH kV COMPUTED TOMOGRAPHY IMAGES FROM LOW DOSE OR LOW kV COMPUTED TOMOGRAPHY IMAGES
WO2017223560A1 (en) 2016-06-24 2017-12-28 Rensselaer Polytechnic Institute Tomographic image reconstruction via machine learning

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7995702B2 (en) 2009-08-25 2011-08-09 General Electric Company System and method of data interpolation in fast kVp switching dual energy CT
WO2017214048A1 (en) 2016-06-06 2017-12-14 Perkinelmer Health Sciences, Inc. Systems and methods for automated sinogram completion, combination, and completion by combination
JP7202302B2 (ja) * 2017-01-05 2023-01-11 ゼネラル・エレクトリック・カンパニイ 断層撮影再構成に使用するためのデータのディープラーニングに基づく推定
WO2019063760A1 (en) * 2017-09-28 2019-04-04 Koninklijke Philips N.V. DISPERSION CORRECTION BASED ON DEEP LEARNING
US11062489B2 (en) * 2018-02-13 2021-07-13 Wisconsin Alumni Research Foundation System and method for multi-architecture computed tomography pipeline
KR102094598B1 (ko) * 2018-05-29 2020-03-27 한국과학기술원 뉴럴 네트워크를 이용한 희소 뷰 전산단층 촬영 영상 처리 방법 및 그 장치
WO2020014701A1 (en) * 2018-07-13 2020-01-16 Loyola University Chicago Phantoms and methods of calibrating dual energy imaging systems therewith
US11039806B2 (en) * 2018-12-20 2021-06-22 Canon Medical Systems Corporation Apparatus and method that uses deep learning to correct computed tomography (CT) with sinogram completion of projection data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004188187A (ja) 2002-11-27 2004-07-08 Ge Medical Systems Global Technology Co Llc アーティファクト低減を容易にする方法及び装置
JP2014503332A (ja) 2011-02-01 2014-02-13 コーニンクレッカ フィリップス エヌ ヴェ デュアルエネルギーct画像再構成のための方法及びシステム
US20160210749A1 (en) 2015-01-21 2016-07-21 Siemens Aktiengesellschaft Method and system for cross-domain synthesis of medical images using contextual deep network
WO2016175755A1 (en) 2015-04-28 2016-11-03 Siemens Healthcare Gmbh METHOD AND SYSTEM FOR SYNTHESIZING VIRTUAL HIGH DOSE OR HIGH kV COMPUTED TOMOGRAPHY IMAGES FROM LOW DOSE OR LOW kV COMPUTED TOMOGRAPHY IMAGES
WO2017223560A1 (en) 2016-06-24 2017-12-28 Rensselaer Polytechnic Institute Tomographic image reconstruction via machine learning

Also Published As

Publication number Publication date
CN111345835A (zh) 2020-06-30
US11039806B2 (en) 2021-06-22
US20210290193A1 (en) 2021-09-23
JP2020099662A (ja) 2020-07-02
US11864939B2 (en) 2024-01-09
CN111345835B (zh) 2024-08-09
US20200196972A1 (en) 2020-06-25

Similar Documents

Publication Publication Date Title
JP7334048B2 (ja) X線ctシステム及び方法
JP7309500B2 (ja) X線ctシステム及び方法
US11176428B2 (en) Apparatus and method for sinogram restoration in computed tomography (CT) using adaptive filtering with deep learning (DL)
JP7455622B2 (ja) 医用画像処理装置及び学習用画像の取得方法
US10925568B2 (en) Apparatus and method using physical model based deep learning (DL) to improve image quality in images that are reconstructed using computed tomography (CT)
JP2020168352A (ja) 医用装置及びプログラム
JP2020116377A (ja) 医用処理装置、医用処理方法及び記憶媒体
US11403791B2 (en) Apparatus and method using deep learning (DL) to improve analytical tomographic image reconstruction
JP7451256B2 (ja) X線システム、画像処理装置及びプログラム
JP7475979B2 (ja) X線システム及び撮像プログラム
EP3671646A1 (en) X-ray computed tomography (ct) system and method
CN115777114A (zh) 针对ct图像去噪的3d-cnn处理
JP7378404B2 (ja) 非スペクトルイメージングシステムを用いたスペクトルイメージング
WO2022268618A1 (en) Multi-energy x-ray imaging with anatomical intelligence
EP3764325B1 (en) System and program for avoiding focal spot blurring
Ström et al. Photon-counting CT reconstruction with a learned forward operator
EP3671647B1 (en) X-ray computed tomography (ct) system and method

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190718

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220117

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20221213

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20230106

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230131

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230403

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230816

R150 Certificate of patent or registration of utility model

Ref document number: 7334048

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150