JP7346429B2 - スペクトルコンピュータ断層撮影(ct)スキャナによって生成される画質が向上したバーチャル非造影画像 - Google Patents

スペクトルコンピュータ断層撮影(ct)スキャナによって生成される画質が向上したバーチャル非造影画像 Download PDF

Info

Publication number
JP7346429B2
JP7346429B2 JP2020541482A JP2020541482A JP7346429B2 JP 7346429 B2 JP7346429 B2 JP 7346429B2 JP 2020541482 A JP2020541482 A JP 2020541482A JP 2020541482 A JP2020541482 A JP 2020541482A JP 7346429 B2 JP7346429 B2 JP 7346429B2
Authority
JP
Japan
Prior art keywords
contrast
image
spectral
training
generate
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
JP2020541482A
Other languages
English (en)
Other versions
JP2021511874A (ja
JPWO2019149711A5 (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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2021511874A publication Critical patent/JP2021511874A/ja
Publication of JPWO2019149711A5 publication Critical patent/JPWO2019149711A5/ja
Application granted granted Critical
Publication of JP7346429B2 publication Critical patent/JP7346429B2/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
    • 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/60Image enhancement or restoration using machine learning, e.g. neural networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/92Dynamic range modification of images or parts thereof based on global image properties
    • 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/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/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Landscapes

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

Description

以下は、概してイメージングに関し、より具体的には、スペクトルコンピュータ断層撮影(CT)スキャナで収集されたデータで生成されたバーチャル非造影画像の画質向上に関する。
CTスキャナは一般に、検出器アレイの反対側で回転式ガントリ上に取り付けられる単一広帯域X線管を含む。X線管は、X線管と検出器アレイとの間にある検査領域の周りを回転し、検査領域を横断する多色放射線を放出する。検出器アレイは、検査領域を横断する放射線を検出し、投影データを生成する。投影データは再構成されてボリュメトリック画像データが生成され、これを使用して1つ以上の2次元画像を生成することができる。画像は、相対的な放射線密度に対応するグレースケール値で表されるピクセルを含む。これらの値は、スキャンされた被験者の減衰特性を反映し、一般に被験者内の解剖学的構造といった構造を示す。
造影CTスキャンは、造影剤が投与された被験者をスキャンすることを含む。造影スキャンを含む手順には、非造影スキャンも含まれている。例えば冠動脈疾患の場合、造影画像は、血管の狭窄及び/又は漏出を検出するのに役立つが、ヨウ素が存在するために血管の石灰化を検出するのにはあまり適していない。血管石灰化の検出のために別の非造影スキャンが行われる。別の例では、腎腫瘤の場合、患者は非造影スキャンと造影腎撮影スキャンとを受け、非造影腎撮影画像と造影腎撮影画像との減衰値が比較されて、腫瘤が亢進しているかどうかが判断される。残念ながら、1回ではなく2回のスキャンを行うと、患者の総X線被ばく量及びスキャン時間が増加し、患者スループットが低下する。
このようなCTスキャナでは、被験者及び/又は物体による放射線の吸収は被験者及び/又は物体を通過する光子のエネルギーに依存するため、検出された放射線はスペクトル情報も含む。このようなスペクトル情報は、被験者及び/又は物体の組織及び/又は材料の元素組成又は材料組成(例えば原子番号)を示す情報といった追加情報を提供することができる。残念ながら、投影データは、検出器アレイによって出力される信号はエネルギースペクトル全体で積分されるエネルギーフルエンスに比例するため、スペクトル特性を反映しない。このようなCTスキャナは、本明細書では非スペクトルCTスキャナとも呼ばれる。スペクトルイメージング用に構成されたCTスキャナは、スペクトル情報を利用する。このようなCTスキャナは、本明細書ではスペクトルCTスキャナと呼ばれる。
スペクトルCTスキャナは、被験者内の材料(例えば骨、軟組織、脂肪、注入された造影剤)のエネルギー依存減衰特性の違いを利用して、材料分解を通じて材料固有の画像を提供する。例示的なスペクトルCTスキャナは、デュアルエネルギー用に構成されており、投影データは、2つの異なる(例えばより高い及びより低い)エネルギーレベルで収集される。例示的なデュアルエネルギー構成には、1)1つのエネルギーレベルで放射X線を放出する1つのX線管と、低エネルギーX線と高エネルギーX線とをそれぞれ検出する2層のX線検出器、2)高速kVスイッチングを有する1つのX線管と、単層検出器、及び、3)互いに角度的に(例えば90度で)オフセットしている2つのX線管/単層検出器のペアが含まれる。
造影スキャンからの材料固有の画像は、「バーチャル(virtual)」非造影(VNC)画像を生成するように処理される。VNC画像は、非造影スキャンからの「真の」非造影(TNC)画像の近似を表す。したがって、スペクトルCTスキャナは、1回の造影スキャンから造影画像とVNC画像との両方を生成することができる。一例では、これは、両タイプのスキャンを必要とする検査のために、非造影スキャンと造影スキャンとの両方を行う必要を軽減する。この結果、患者の総X線被ばく量及びスキャン時間が低減し、患者スループットを増加させることができる。
しかし、VNC画像の画質は、TNC画像に比べて低下する傾向がある。例えばVNC画像では、信号が、グラウンドトゥルースよりも小さくなったり、TNC画像に比べて欠落したりすることがある。文献では、TNC画像には両方が含まれているのに、同じVNC画像において1つの腎臓結石が存在し、別の腎臓結石が欠落している例が示されている。更に、VNC画像は、TNC画像に比べてノイズ及びアーチファクトが多くなる傾向がある。ノイズが多いと感度が低下し、検出における偽陽性率が増加する可能性がある。アーチファクトの例には、ビーム硬化及びリムアーティファクトが含まれる。残念ながら、このような画質低下は診断に影響を与える可能性がある。
本明細書に説明する態様は、上記問題及びその他に対処する。
一態様では、スペクトルイメージングシステムが、放射X線を放出する放射線源と、放射X線を検出し、スペクトルデータを生成する検出器アレイとを含む。スペクトルイメージングシステムは更に、画質が向上したバーチャル非造影画像を生成するようにトレーニングされたニューラルネットワークを含むコンピュータ実行可能命令を含むバーチャル非造影画像向上モジュールを格納するメモリを含む。ニューラルネットワークは、トレーニング用スペクトルデータと、非造影スキャンから生成されるトレーニング用非造影画像とでトレーニングされる。スペクトルイメージングシステムは更に、トレーニングされたニューラルネットワークでスペクトルデータを処理して画質が向上したバーチャル非造影画像を生成するプロセッサを含む。
別の態様では、スペクトルイメージングシステムが、放射X線を放出する放射線源と、放射X線を検出し、スペクトルデータを生成する検出器アレイとを含む。スペクトルイメージングシステムは更に、ニューラルネットワークを含むコンピュータ実行可能命令を含むバーチャル非造影画像向上モジュールを格納するメモリを含む。スペクトルイメージングシステムは更に、トレーニング用スペクトルデータと、非造影スキャンから生成されるトレーニング用非造影画像とで画質が向上したバーチャル非造影画像を生成するようにニューラルネットワークをトレーニングするプロセッサを含む。
別の態様では、コンピュータ可読記憶媒体が、コンピューティングシステムのプロセッサによって実行されると、当該プロセッサに、放射線源で放射X線を放出させ、検出器アレイで、放出された放射X線を検出させ、スペクトルデータを生成させ、画質が向上したバーチャル非造影画像を生成するようにトレーニングされたニューラルネットワークをトレーニングさせ、画質が向上したバーチャル非造影画像を生成するように、トレーニングされたニューラルネットワークでスペクトルデータを処理させるコンピュータ可読命令でエンコードされる。
当業者であれば、添付の説明を読んで理解すると、本願の更に他の態様を認識するであろう。
本発明は、様々なコンポーネント及びコンポーネントの構成、また、様々なステップ及びステップの構成の形を取ってよい。図面は、好適な実施形態を例示するに過ぎず、本発明を限定するものと解釈されるべきではない。
図1は、スキャナによって生成されたVNC画像の画質を向上させるようにトレーニングされたニューラルネットワークを含むVNC画像向上モジュールを有する例示的なスペクトルコンピュータ断層撮影スキャナを概略的に示す。 図2は、トレーニング用VNC画像のセット及びトレーニング用TNC画像のセットを用いるニューラルネットワークのトレーニングの一例を概略的に示す。 図3は、図2のトレーニングされたニューラルネットワークを、スペクトルコンピュータ断層撮影スキャナによって生成されたVNC画像と共に使用して、画質が向上したVNC画像を生成する一例を概略的に示す。 図4は、トレーニング用デュアルエネルギー画像のセット及びトレーニング用TNC画像のセットを用いるニューラルネットワークのトレーニングの一例を概略的に示す。 図5は、図4のトレーニングされたニューラルネットワークを、スペクトルコンピュータ断層撮影スキャナによって生成された造影デュアルエネルギー画像と共に使用して、画質が向上したVNC画像を生成する一例を概略的に示す。 図6は、トレーニング用デュアルエネルギー投影データのセット及びTNC画像のセットを用いるニューラルネットワークのトレーニングの一例を概略的に示す。 図7は、図6のトレーニングされたニューラルネットワークを、スペクトルコンピュータ断層撮影スキャナによって生成された造影デュアルエネルギー投影データと共に使用して、画質が向上したVNC画像を生成する一例を概略的に示す。 図8は、トレーニング用VNC画像及びトレーニング用TNC画像を用いてニューラルネットワークをトレーニングし、トレーニングされたニューラルネットワークを使用して画質が向上したVNC画像を生成する例示的な方法を示す。 図9は、トレーニング用デュアルエネルギー画像及びトレーニング用TNC画像を用いてニューラルネットワークをトレーニングし、トレーニングされたニューラルネットワークを使用して画質が向上したVNC画像を生成する例示的な方法を示す。 図10は、トレーニング用デュアルエネルギー投影データ及びTNC画像を用いてニューラルネットワークをトレーニングし、トレーニングされたニューラルネットワークを使用して画質が向上したVNC画像を生成する例示的な方法を示す。
以下は、造影スペクトルスキャン中にスペクトルCTスキャナで収集されたデータで生成されたバーチャル非造影(VNC)画像の画質を向上させるアプローチについて説明する。スペクトルCTスキャナからの低kVpデータ及び高kVpデータといった造影デュアルエネルギースペクトルデータは再構成されて、造影低kVp画像及び造影高kVp画像が生成される。これらは分解されて、VNC画像を含むスペクトル画像が生成される。2つのkVp画像は、例えば造影剤(つまり、光電吸収)によるビーム硬化アーチファクトを含む。kVp画像の再構成中にビーム硬化補償が適用されても、2つのkVp画像は、残留ビーム硬化アーチファクトを含む可能性がある。残留ビーム硬化アーチファクトは、分解中に更に低減することができるが、結果として得られるVNC画像は、依然として残留ビーム硬化アーチファクトを含み、これは、画質を低下する可能性がある。
本明細書に説明する手法では、ニューラルネットワークが、入力としてのトレーニング用スペクトル画像のセットと、参照データとしてのトレーニング用TNC画像のセットとでトレーニングされる。トレーニング中、ニューラルネットワークは、トレーニング用スペクトル画像のセットをトレーニング用TNC画像のセットにマッピングすることを学習する。一例では、トレーニング用スペクトル画像のセットは、ノイズレベルの異なるスペクトル画像を含み、参照データとしてのトレーニング用TNC画像のセットは、低ノイズTNC画像を含む。マッピングは非線形であるので、入力VNC画像における残留ビーム硬化アーチファクトは更に低減されるか、及び/又は、最終VNC画像には現れず、入力VNC画像におけるノイズは最終VNC画像では低減される。したがって、このアプローチでは、造影画像と、TNC画像と同様又は同等の画質を有するVNC画像とを生成することができる。
図1は、スペクトルCTイメージング、例えばデュアルエネルギーイメージングといったマルチエネルギーイメージング用に構成された例示的なコンピュータ断層撮影(CT)スキャナ102を概略的に示す。CTスキャナ102は、固定ガントリ104と、回転ガントリ106とを含み、回転ガントリ106は、固定ガントリ104によって回転可能に支持され、長手方向又はZ軸110を中心に検査領域108(及びその中の物体又は被験者の一部)の周りを回転する。
X線管といった放射線源112が、回転ガントリ106によって支持され、回転ガントリ106と共に検査領域108の周りを回転する。放射線源112は、コリメートされて、検査領域108を横断する概ね扇形、くさび形又は円錐形の放射X線ビームを形成する放射X線を放出する。一例では、放射線源112は、関心の単一の選択されたピーク放出電圧(kVp)について広帯域(多色)放射線を放出する単一のX線管である。
別の例では、放射線源112は、スキャン中に少なくとも2つの異なる放出電圧(例えば70keV、100keV、120keV等)間で切り替わる。更に別の例では、放射線源112は、回転ガントリ106上で角度オフセットされ、それぞれが異なる平均エネルギースペクトルを有する放射線を放出する2つ以上のX線管を含む。米国特許第8,442,184B2号は、kVpスイッチング及び複数のX線管を有するシステムを説明しており、参照によりその全体が本明細書に組み込まれる。
放射線感応検出器アレイ114が、検査領域108を挟んで放射線源112の反対側にある弧に対応する。検出器アレイ114は、Z軸方向に沿って互いに対して配置される1つ以上の検出器列を含み、検査領域108を横断する放射線を検出する。造影スキャンの場合、検出器アレイ114は、造影高エネルギー投影データ及び造影低エネルギー投影データといった造影スペクトル投影データ(線積分)を生成する。
この例では、検出器アレイ114は、多層シンチレータ/光センサ検出器(例えば参照によりその全体が本明細書に組み込まれる米国特許第7,968,853B2号)といったエネルギー分解検出器を含む。変形例では、検出器アレイ114は、光子計数(直接変換)検出器(例えば参照によりその全体が本明細書に組み込まれる国際特許公開WO2009/072056A2)を含む。これらの例では、放射線源112は、広帯域源、kVpスイッチング及び/又は複数のX線管放射線源を含む。検出器アレイ114が、非エネルギー分解検出器を含む場合、放射線源112は、kVpスイッチング及び/又は複数のX線管放射線源を含む。
再構成器116が、造影スペクトル投影データを処理し、造影高エネルギーボリュメトリック画像データ及び造影低エネルギーボリュメトリック画像データといった造影スペクトルボリュメトリック画像データを生成する。再構成器116はまた、例えば最初に高エネルギー投影データと低エネルギー投影データとを組み合わせた後に、組み合わされた投影データを再構成することによって、及び/又は、高エネルギーボリュメトリック画像と低エネルギーボリュメトリック画像データとを組み合わせることによって、造影非スペクトルボリュメトリック画像データを生成することもできる。そこからスペクトル画像及び非スペクトル画像を導き出すことができる。再構成器116は、中央処理演算ユニット(CPU)、マイクロプロセッサといったプロセッサで実装することができる。
図示される実施形態では、基礎材料分解器118が、造影投影データ(投影領域分解)、並びに/又は、造影スペクトルボリュメトリック画像データ及び/若しくは画像(画像領域分解)を基礎成分に分解する。投影領域分解用に構成される場合、再構成器116は、投影データ基礎成分を再構成して、VNCボリュメトリック画像データ及び/又は画像といった基礎成分ボリュメトリック画像データ及び/又は画像を生成する。画像領域分解用に構成される場合、基礎材料分解器118は、造影スペクトルボリュメトリック画像データ及び/又は画像を、VNCボリュメトリック画像データ及び/又は画像といった基礎成分ボリュメトリック画像データ及び/又は画像に分解する。変形例では、別々の基礎材料分解器が投影領域分解と画像領域分解とに使用される。別の変形例では、基礎材料分解器118は省略されるか又は使用されない。
1つの手法では、基礎材料分解器118は、それぞれが水-ヨウ素基礎材料ペア又はカルシウム-ヨウ素基礎材料ペアを想定する2つの材料分解アルゴリズムを使用して、ヨウ素なしのVNC画像を生成する。別のアプローチは、N次元クラスタ分析を行って、画像を、ヨウ素や、軟組織、カルシウム等を含む他の材料といった成分に分解することである。他の例は、Song他の「Virtual Non-Contrast CT Using Dual-Energy Spectral CT:Feasibility of Coronary Artery Calcium Scoring」、Korean Journal of Radiology2016、17(3)、321-329及び2014年1月13日に出願され、「Image Processing for Spectral CT」という名称の米国特許出願公開第2014/0133729Al号に説明されている。これは、参照によりその全体が本明細書に組み込まれる。分解器118は、プロセッサ等で実装することができる。
カウチといった被験者支持体120が、検査領域108内で被験者又は物体(例えばファントム)122を支持する。被検者支持体120は、被検者又は物体のローディング、スキャン及び/又はアンローディングのために被験者及び/又は物体を検査領域108に対して誘導するように、イメージング手順の実行と協調して移動可能である。
インジェクタ124が、灌流スキャンのためにスキャンされる被験者又は物体122に、1つ以上の(例えばヨウ素)造影剤といった材料を注入又は投与する。造影剤は、更に又は或いは、臨床医等によって手動で投与されてもよい。造影剤が被験者又は物体122に手動で投与される場合は、インジェクタ124は省略されてよい。
オペレータコンソール126は、ディスプレイモニタ、フィルマ等といった人間が読み取り可能な出力デバイス128と、キーボード、マウス等といった入力デバイス130とを含む。コンソール126は更に、プロセッサ132(例えばCPU、マイクロプロセッサ等)と、物理メモリといったコンピュータ可読記憶媒体134(一時的媒体を除く)を含む。図示される実施形態では、コンピュータ可読記憶媒体134は、VNC画像向上モジュール136を含む。
図示されるVNC画像向上モジュール136は、VNC画像を処理して、非造影スキャンから生成されたTNC画像と同様又は同等の向上された画質のVNC画像(画質が向上したバーチャル非造影画像)を生成するためのコンピュータ実行可能命令を含む。以下でより詳細に説明するように、命令は、トレーニング用スペクトルデータ及びTNC画像でトレーニングされたニューラルネットワークを含む。したがって、前述のように、スペクトルCTスキャナ102は、1回の造影スキャンで収集されたデータを使用して、造影画像と、非造影スキャンから生成されたTNC画像と診断上同様又は同等の画質を有するVNC画像とを生成することができる。
図2から図7は、ニューラルネットワークをトレーニングするための様々な実施形態を概略的に示す。これらの実施形態では、ニューラルネットワークは、異なる造影剤送達で収集されたデータを有するトレーニング用データセットを使用することによって、入力データにおける造影剤の送達の違いを考慮するようにトレーニングされる。一例では、トレーニング用データセットは、データ特徴(即ち、造影剤送達)の分布が、実際の臨床シナリオにおけるデータ特徴の分布をよく表すように構築される。トレーニング用データセットは、実際の臨床スキャンとシミュレーションとの両方から得ることができる。シミュレーションの利点は、すべてのデータセット(造影スペクトルCTデータ及びTNC画像を含む)における造影剤送達を操作することができるため、データ特徴の推定分布に厳密に従う多数のデータセットを生成することができる点である。別の例では、入力及びトレーニング用データは、造影剤送達がデータ全体で一貫しているように収集することができる。
図2は、材料分解を使用して生成されたトレーニング用VNC画像及び非造影CTスキャンからのトレーニング用TNC画像を用いるニューラルネットワーク202の例示的なトレーニングを概略的に示す。ニューラルネットワークのパラメータは更新され、目的関数204を使用して多数の反復を通じて収束する。目的関数204は、一例では、ニューラルネットワーク202によって出力されたVNCボリュメトリック画像データ又は画像とTNC画像との間の誤差を、例えばニューラルネットワーク202によって出力されたVNCボリュメトリック画像データ又は画像とTNC画像との間の平均二乗差及び/又は他の差に基づいて最小限にする数学関数を含む。この場合、パラメータは、誤差が所定閾値を下回るまで更新される。他の停止基準には、所定の回数の反復、所定の期間等が含まれる。
図3は、VNC画像向上モジュール136が図2のトレーニングされたニューラルネットワーク202を含む一実施形態を概略的に示す。トレーニングされたニューラルネットワーク202は、基礎材料分解器118(図1)によって生成されたVNC画像(TNC画像に比べてノイズ及び/又はアーチファクトが多い)を受け取り、入力VNC画像よりもノイズ及び/又はアーチファクトが少なく、また、診断に十分な画質を有するVNC画像を含む画質が向上したVNC画像を出力する。
図4は、ニューラルネットワーク202への入力が、トレーニング用デュアルエネルギー(高エネルギー及び低エネルギー)画像である図2の変形例を概略的に示す。図5は、VNC画像向上モジュール136が図4でトレーニングされたニューラルネットワーク202を含む一実施形態を概略的に示す。トレーニングされたニューラルネットワーク202は、再構成器116(図1)によって生成された造影デュアルエネルギー画像を受け取り、画質が向上したVNC画像を出力する。この変形例では、基礎材料分解器118は、図1のシステム102から省略されるか又は使用されない。
この変形例では、ニューラルネットワーク202は、データの非線形マッピングを行うことができるため、材料分解を担う。材料分解は、近似の後に(例えば高エネルギー及び低エネルギー用の)多項式を解くことを伴うことがよくある。材料分解にニューラルネットワーク202を使用することによって、基本材料画像におけるノイズを従来の材料分解からのノイズよりも少なくすることができる。これは、例えば従来の材料分解では、一度に1つの検出器ピクセルで測定されたデータか又は1つの画像ボクセルの再構成値しか使用されないが、ニューラルネットワークは、2つ以上のピクセル又はボクセルを含む小さいパッチを使用するからである。
一般に、投影データは、データに隠された固有のスペクトル情報を含む。例えば特定の光線経路に沿った減衰は、高Zで低密度の物体と低Zで高密度の物体との両方で同じである可能性があるが、物理的効果、即ち、コンプトン散乱及び光電吸収は異なり、低Zで高密度の物体のコンプトン散乱はより多く、これらの違いは収集された生の投影データに隠されている。別の例は、低Z材料に比べて高Z材料では、ビーム硬化アーチファクトが多いことである。これらの違いは、例えば効果的な散乱補正、ビーム硬化補償等によって低減及び/又は除去され、再構成されたボリュメトリック画像データには現れない。
変形例では、散乱補正及び/又はビーム硬化補償といったスペクトル情報を除去する補正は、再構成から省略される。例えば別の例では、図2及び図3のニューラルネットワーク202は、散乱補正及び/又はビーム硬化補償といった固有のスペクトル情報を除去する補正を使用せずに再構成されたボリュメトリック画像データから生成されたVNC画像でトレーニングされ、及び/又は、図4及び図5のニューラルネットワーク202は、散乱補正及び/又はビーム硬化補償といった固有のスペクトル情報を除去する補正を使用せずに再構成された高エネルギー画像及び低エネルギー画像でトレーニングされる。これにより、トレーニングで固有のスペクトル情報を考慮することができるようになる。
図6は、ニューラルネットワーク202への入力がトレーニング用デュアルエネルギー(高エネルギー及び低エネルギー)投影データである図2の変形例を概略的に示す。図7は、VNC画像向上モジュール136が図6でトレーニングされたニューラルネットワーク202を含む一実施形態を概略的に示す。トレーニングされたニューラルネットワーク202は、検出器アレイ114(図1)によって生成された造影デュアルエネルギー投影データを受け取り、画質が向上したVNC画像を出力する。この変形例では、基礎材料分解器118は、図1のシステム102から省略されるか又は使用されない。
一変形例は、図2から図7に説明した3つの実施形態のうちの2つを含む。例えば1つの変形例は、図2及び図4のニューラルネットワークを含む。別の変形例は、図2及び図6のニューラルネットワークを含む。別の変形例は、図4及び図6のニューラルネットワークを含む。別の変形は、図2から図7に説明した3つの実施形態の3つすべてを含む。別の変形例は、図2及び図4、図2及び図6、図2及び図6、又は、図2、図4及び図6におけるトレーニング用データといった様々なタイプのトレーニング用データでトレーニングされた1つのニューラルネットワークを含む。これらの変形例は、画質が向上したVNC画像を出力するように、図3、図5及び/又は図7の入力データを用いて使用することができる。
適切なニューラルネットワークの例は、Gouk他の「Fast Sliding Window Classification with Convolutional Neural Networks」、IVNVZ‘14 Proceedings of the 29th International Conference on Image and Vision Computing New Zealand、114~118頁、2014年11月19~21日、「Fully convolutional networks for semantic segmentation」、Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition、2015年、及び、Ronneberger他の「U-Net:Convolution Networks for Biomedical Image Segmentation」、Medical Image Computing and Computer-Assisted Intervention(MICCAI)、Springer、LNCS、第9351巻、234--241、2015年に説明されている。
図8は、本明細書における一実施形態による例示的な方法を示す。
この方法のステップの順序は限定ではないことを理解されたい。したがって、本明細書では他の順序も考えられる。更に、1つ以上のステップが省略されても及び/又は1つ以上の追加のステップが含まれてもよい。
ステップ802において、ここで及び/又は他で説明するように、トレーニング用VNC画像のセットを受け取る。
ステップ804において、ここで及び/又は他で説明するように、トレーニング用TNC画像のセットを受け取る。
ステップ806において、ここで及び/又は他で説明するように、トレーニング用VNC画像のセット及びトレーニング用TNC画像のセットを使用して、ニューラルネットワークをトレーニングする。
ステップ808において、ここで及び/又は他で説明するように、VNC画像を受け取る。
ステップ810において、受け取ったVNC画像をトレーニングされたニューラルネットワークで処理して、画質が向上したVNC画像を生成する。
図9は、本明細書における一実施形態による別の例示的な方法を示す。
この方法のステップの順序は限定ではないことを理解されたい。したがって、本明細書では他の順序も考えられる。更に、1つ以上のステップが省略されても及び/又は1つ以上の追加のステップが含まれてもよい。
ステップ902において、本明細書に説明するように及び/又は他の方法で、トレーニング用デュアルエネルギー画像のセットを受け取る。
ステップ904において、ここで及び/又は他で説明するように、トレーニング用TNC画像のセットを受け取る。
ステップ906において、ここで及び/又は他で説明するように、トレーニング用デュアルエネルギー画像のセット及びトレーニング用TNC画像のセットを使用して、ニューラルネットワークをトレーニングする。
ステップ908において、ここで及び/又は他で説明するように、造影デュアルエネルギー画像を受け取る。
ステップ910において、造影デュアルエネルギー画像をトレーニングされたニューラルネットワークで処理して、画質が向上したVNC画像を生成する。
図10は、本明細書における一実施形態による別の例示的な方法を示す。
この方法のステップの順序は限定ではないことを理解されたい。したがって、本明細書では他の順序も考えられる。更に、1つ以上のステップが省略されても及び/又は1つ以上の追加のステップが含まれてもよい。
ステップ1002において、ここで及び/又は他で説明するように、トレーニング用デュアルエネルギー投影データのセットを受け取る。
ステップ1004において、ここで及び/又は他で説明するように、トレーニング用TNC画像のセットを受け取る。
ステップ1006において、ここで及び/又は他で説明するように、トレーニング用デュアルエネルギー投影データのセット及びトレーニング用TNC画像のセットを使用して、ニューラルネットワークをトレーニングする。
ステップ1008において、ここで及び/又は他で説明するように、造影デュアルエネルギー投影データを受け取る。
ステップ1010において、造影デュアルエネルギー投影データをトレーニングされたニューラルネットワークで処理して、画質が向上したVNC画像を生成する。
上記は、コンピュータ可読記憶媒体にエンコードされ又は埋め込まれ、コンピュータプロセッサによって実行されると当該プロセッサに、説明したステップを行わせるコンピュータ可読命令によって実施することができる。更に又は或いは、コンピュータ可読命令の少なくとも1つは、信号、搬送波又はコンピュータ可読記憶媒体ではない他の一時的媒体によって運ばれる。
本発明は、図面及び上記説明において詳細に例示され、説明されたが、当該例示及び説明は、例示的に見なされるべきであり、限定的に見なされるべきではない。本発明は、開示した実施形態に限定されない。開示した実施形態の他の変形態様は、図面、開示内容及び添付の請求項の検討から、請求項に係る発明を実施する当業者によって理解され、実施される。
請求項において、「含む」との用語は、他の要素又はステップを排除するものではなく、また、単数形は、複数形を排除するものではない。単一のプロセッサ又は他のユニットが、請求項に引用される幾つかのアイテムの機能を果たしてもよい。特定の手段が相互に異なる従属請求項に記載されることだけで、これらの手段の組み合わせを有利に使用することができないことを示すものではない。
コンピュータプログラムは、他のハードウェアと共に又はその一部として供給される光学記憶媒体又はソリッドステート媒体といった適切な媒体上に記憶及び/又は分散されてもよいが、インターネット又は他の有線若しくは無線通信システムを介するといった他の形式で分散されてもよい。請求項における任意の参照符号は、範囲を限定するものと解釈されるべきではない。

Claims (13)

  1. 放射X線を放出する放射線源と、
    放射X線を検出し、造影マルチエネルギースペクトル投影データを生成する検出器アレイと、
    非造影スキャンから生成されるトレーニング用非コントラスト造影画像と同様又は同等の画質が向上したバーチャル非造影画像を生成するようにトレーニングされたニューラルネットワークを含むコンピュータ実行可能命令を含むバーチャル非造影画像向上モジュールを格納するメモリであって、前記ニューラルネットワークは、スペクトルスキャンから生成されるトレーニング用造影マルチエネルギースペクトル投影データと、前記非造影スキャンから生成される前記トレーニング用非造影画像とでトレーニングされる、メモリと、
    トレーニングされた前記ニューラルネットワークで前記造影マルチエネルギースペクトル投影データを処理して前記画質が向上したバーチャル非造影画像を生成するプロセッサと、
    を含む、スペクトルコンピュータ断層撮影イメージングシステム。
  2. 前記ニューラルネットワークは、パラメータを更新して、前記トレーニング用造影マルチエネルギースペクトル投影データから生成されるバーチャル非造影画像と、前記トレーニング用非造影画像との間の誤差を低減する、請求項1に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  3. 前記造影マルチエネルギースペクトル投影データを再構成し、造影マルチエネルギースペクトル画像を生成する再構成器を更に含み、
    前記プロセッサは更に、トレーニングされた前記ニューラルネットワークで前記造影マルチエネルギースペクトル画像を処理して前記画質が向上したバーチャル非造影画像を生成する、請求項1に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  4. 前記造影マルチエネルギースペクトル画像を処理して初期のバーチャル非造影画像を生成する基礎材料分解器を更に含み、
    前記プロセッサは更に、トレーニングされた前記ニューラルネットワークで前記初期のバーチャル非造影画像を処理して前記画質が向上したバーチャル非造影画像を生成する、請求項3に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  5. 前記造影マルチエネルギースペクトル投影データを処理してバーチャル非造影スペクトル投影データを生成する基礎材料分解器と、
    前記バーチャル非造影スペクトル投影データを再構成して初期のバーチャル非造影画像を生成する再構成器と、
    を更に含み、
    前記プロセッサは更に、トレーニングされた前記ニューラルネットワークで前記初期のバーチャル非造影画像を処理して前記バーチャル非造影画像を生成する、請求項3に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  6. 前記スペクトルスキャンから生成される前記トレーニング用造影マルチエネルギースペクトル投影データが処理されてバーチャル非造影画像を生成し、前記ニューラルネットワークは、パラメータを更新して、トレーニング用の前記バーチャル非造影画像から生成されるバーチャル非造影画像と、前記トレーニング用非造影画像との間の誤差を低減する、請求項又はに記載のスペクトルコンピュータ断層撮影イメージングシステム。
  7. 放射X線を放出する放射線源と、
    放射X線を検出し、造影マルチエネルギースペクトル投影データを生成する検出器アレイと、
    ニューラルネットワークを含むコンピュータ実行可能命令を含むバーチャル非造影画像向上モジュールを格納するメモリと、
    スペクトルスキャンから生成されるトレーニング用造影マルチエネルギースペクトル投影データと、非造影スキャンから生成されるトレーニング用非造影画像とで、前記非造影スキャンから生成される前記トレーニング用非コントラスト造影画像と同様又は同等の画質が向上したバーチャル非造影画像を生成するように前記ニューラルネットワークをトレーニングするプロセッサと、
    を含む、スペクトルコンピュータ断層撮影イメージングシステム。
  8. 前記ニューラルネットワークは、パラメータを更新して、前記トレーニング用造影マルチエネルギースペクトル投影データから生成されるバーチャル非造影画像と、前記トレーニング用非造影画像との間の誤差を低減し、オプションで、
    前記プロセッサは更に、トレーニングされた前記ニューラルネットワークで前記造影マルチエネルギースペクトル投影データを処理して前記画質が向上したバーチャル非造影画像を生成する、請求項に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  9. 前記スペクトルスキャンから生成される前記トレーニング用造影マルチエネルギースペクトル投影データが処理されてバーチャル非造影画像を生成し、
    前記ニューラルネットワークは、パラメータを更新して、トレーニング用の前記バーチャル非造影画像から生成されるバーチャル非造影画像と、前記トレーニング用非造影画像との間の誤差を低減する、請求項に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  10. 前記スペクトルコンピュータ断層撮影イメージングシステムは更に、
    前記造影マルチエネルギースペクトル投影データを再構成し、造影マルチエネルギースペクトル画像を生成する再構成器と、
    前記造影マルチエネルギースペクトル画像を処理して初期のバーチャル非造影画像を生成する基礎材料分解器と、
    を含み、
    前記プロセッサは更に、トレーニングされた前記ニューラルネットワークで前記初期のバーチャル非造影画像を処理して前記画質が向上したバーチャル非造影画像を生成する、請求項に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  11. 前記スペクトルコンピュータ断層撮影イメージングシステムは更に、
    造影マルチエネルギースペクトル画像を処理してバーチャル非造影スペクトル投影データを生成する基礎材料分解器と、
    前記バーチャル非造影スペクトル投影データを再構成して初期のバーチャル非造影画像を生成する再構成器と、
    を含み、
    前記プロセッサは更に、トレーニングされた前記ニューラルネットワークで前記初期のバーチャル非造影画像を処理して前記画質が向上したバーチャル非造影画像を生成する、請求項に記載のスペクトルコンピュータ断層撮影イメージングシステム。
  12. コンピューティングシステムのプロセッサによって実行されると、前記プロセッサに、
    放射線源で放射X線を放出させ、
    検出器アレイで、放出された放射X線を検出させ、造影スペクトル投影データを生成させ、
    非造影スキャンから生成されるトレーニング用非コントラスト造影画像と同様又は同等の画質が向上したバーチャル非造影画像を生成するようにトレーニングされたニューラルネットワークをトレーニングさせ、前記ニューラルネットワークは、スペクトルスキャンから生成されるトレーニング用造影マルチエネルギースペクトル投影データと、前記非造影スキャンから生成される前記トレーニング用非造影画像とでトレーニングされ、
    前記画質が向上したバーチャル非造影画像を生成するように、トレーニングされた前記ニューラルネットワークで前記スペクトルデータを処理させる、コンピュータ可読命令でエンコードされた、コンピュータ可読記憶媒体。
  13. 前記プロセッサは、患者スキャンからの造影剤送達の分布を表す造影剤送達の分布を有するトレーニング用データを使用して前記ニューラルネットワークをトレーニングする、請求項12に記載のコンピュータ可読記憶媒体。
JP2020541482A 2018-01-31 2019-01-30 スペクトルコンピュータ断層撮影(ct)スキャナによって生成される画質が向上したバーチャル非造影画像 Active JP7346429B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201862624415P 2018-01-31 2018-01-31
US62/624,415 2018-01-31
PCT/EP2019/052157 WO2019149711A1 (en) 2018-01-31 2019-01-30 Image quality improved virtual non-contrast images generated by a spectral computed tomography (ct) scanner

Publications (3)

Publication Number Publication Date
JP2021511874A JP2021511874A (ja) 2021-05-13
JPWO2019149711A5 JPWO2019149711A5 (ja) 2022-02-03
JP7346429B2 true JP7346429B2 (ja) 2023-09-19

Family

ID=65243573

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020541482A Active JP7346429B2 (ja) 2018-01-31 2019-01-30 スペクトルコンピュータ断層撮影(ct)スキャナによって生成される画質が向上したバーチャル非造影画像

Country Status (6)

Country Link
US (1) US11813101B2 (ja)
EP (1) EP3746979B1 (ja)
JP (1) JP7346429B2 (ja)
CN (1) CN111684486B (ja)
ES (1) ES2968760T3 (ja)
WO (1) WO2019149711A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102019200270A1 (de) * 2019-01-11 2020-07-16 Siemens Healthcare Gmbh Bereitstellen eines Differenzbilddatensatzes und Bereitstellen einer trainierten Funktion
DE102019200269A1 (de) 2019-01-11 2020-07-16 Siemens Healthcare Gmbh Bereitstellen eines Beschränkungsbilddatensatzes und/oder eines Differenzbilddatensatzes
JP7433809B2 (ja) * 2019-08-21 2024-02-20 キヤノンメディカルシステムズ株式会社 学習済みモデルの生成方法、および医用処理装置
JP7376053B2 (ja) * 2021-02-01 2023-11-08 クラリピーアイ インコーポレイテッド ディープラーニング基盤の造影増強ctイメージ対照度増幅装置及び方法
JP2022158574A (ja) * 2021-04-02 2022-10-17 富士フイルム株式会社 学習装置、方法およびプログラム、並びに放射線画像処理装置、方法およびプログラム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100034446A1 (en) 2008-08-06 2010-02-11 Siemens Corporate Research,Inc. System and Method for Coronary Digital Subtraction Angiography
US20110064292A1 (en) 2009-09-17 2011-03-17 Siemens Corporation System and method for multi-image based virtual non-contrast image enhancement for dual source ct
JP2018089301A (ja) 2016-12-07 2018-06-14 学校法人常翔学園 生体画像処理装置、出力画像製造方法、学習結果製造方法、及びプログラム

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005009206A2 (en) 2003-06-25 2005-02-03 Besson Guy M Dynamic multi-spectral imaging system
EP1876955B1 (en) 2005-04-26 2016-11-23 Koninklijke Philips N.V. Double decker detector for spectral ct
WO2009072056A2 (en) 2007-12-04 2009-06-11 Koninklijke Philips Electronics N. V. Monolithically integrated crystalline direct-conversion semiconductor detector for detecting incident x-radiation at ultra-fine pitch and method for manufacturing such an x-ray semiconductor detector
US8442184B2 (en) 2008-06-30 2013-05-14 Koninklijke Philips Electronics N.V. Spectral CT
JP6100772B2 (ja) 2011-07-15 2017-03-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 画像処理方法及びコンピューティング装置
CN107530037B (zh) * 2015-12-17 2018-11-27 皇家飞利浦有限公司 用于生成对比剂浓度图的方法
CN106097355A (zh) 2016-06-14 2016-11-09 山东大学 基于卷积神经网络的胃肠道肿瘤显微高光谱图像处理方法
WO2017223560A1 (en) * 2016-06-24 2017-12-28 Rensselaer Polytechnic Institute Tomographic image reconstruction via machine learning
CN106600568B (zh) 2017-01-19 2019-10-11 东软医疗系统股份有限公司 一种低剂量ct图像去噪方法及装置
CN107374657B (zh) * 2017-06-30 2021-05-11 上海联影医疗科技股份有限公司 对ct扫描数据进行校正的方法及ct扫描系统
CN107516330B (zh) 2017-08-31 2021-02-09 上海联影医疗科技股份有限公司 模型生成方法、图像处理方法及医学成像设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100034446A1 (en) 2008-08-06 2010-02-11 Siemens Corporate Research,Inc. System and Method for Coronary Digital Subtraction Angiography
US20110064292A1 (en) 2009-09-17 2011-03-17 Siemens Corporation System and method for multi-image based virtual non-contrast image enhancement for dual source ct
JP2018089301A (ja) 2016-12-07 2018-06-14 学校法人常翔学園 生体画像処理装置、出力画像製造方法、学習結果製造方法、及びプログラム

Also Published As

Publication number Publication date
ES2968760T3 (es) 2024-05-13
EP3746979B1 (en) 2023-11-15
CN111684486B (zh) 2024-05-21
US20200367844A1 (en) 2020-11-26
EP3746979A1 (en) 2020-12-09
US11813101B2 (en) 2023-11-14
WO2019149711A1 (en) 2019-08-08
CN111684486A (zh) 2020-09-18
JP2021511874A (ja) 2021-05-13

Similar Documents

Publication Publication Date Title
JP7346429B2 (ja) スペクトルコンピュータ断層撮影(ct)スキャナによって生成される画質が向上したバーチャル非造影画像
US10147168B2 (en) Spectral CT
Clark et al. Micro-CT of rodents: state-of-the-art and future perspectives
US9498179B1 (en) Methods and systems for metal artifact reduction in spectral CT imaging
Weidinger et al. Polychromatic iterative statistical material image reconstruction for photon‐counting computed tomography
US8615120B2 (en) Spectral imaging
CN111818851B (zh) 被配置为生成谱体积图像数据的非谱计算机断层摄影(ct)扫描器
US11419566B2 (en) Systems and methods for improving image quality with three-dimensional scout
Erath et al. Deep learning‐based forward and cross‐scatter correction in dual‐source CT
JP6945410B2 (ja) 画像再構成処理装置、x線コンピュータ断層撮像装置及び画像再構成処理方法
CN115067981A (zh) 双能ct物质分解方法、装置、系统、电子设备及存储介质
Nuyts et al. Reduction of attenuation correction artifacts in PET-CT
US11361415B2 (en) Material-selective adaptive blending of volumeiric image data
CN116172594B (zh) 用于生成患者的结果图像数据集的方法及装置
Yang et al. Sinogram-based attenuation correction in PET/CT
Maass et al. Material decomposition with inconsistent rays (MDIR) for cone-beam dual energy CT
Little Geometric modeling and sinogram restoration methods in computed tomography

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220126

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220126

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221007

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20221221

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230405

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230424

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230721

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230906

R150 Certificate of patent or registration of utility model

Ref document number: 7346429

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150