JP7209003B2 - 定量分析が向上された低放射線量コンピュータ断層撮影灌流(ctp) - Google Patents

定量分析が向上された低放射線量コンピュータ断層撮影灌流(ctp) Download PDF

Info

Publication number
JP7209003B2
JP7209003B2 JP2020541568A JP2020541568A JP7209003B2 JP 7209003 B2 JP7209003 B2 JP 7209003B2 JP 2020541568 A JP2020541568 A JP 2020541568A JP 2020541568 A JP2020541568 A JP 2020541568A JP 7209003 B2 JP7209003 B2 JP 7209003B2
Authority
JP
Japan
Prior art keywords
perfusion
timeline
image data
sequential
volumetric image
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
JP2020541568A
Other languages
English (en)
Other versions
JPWO2019149718A5 (ja
JP2021511881A (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 JP2021511881A publication Critical patent/JP2021511881A/ja
Publication of JPWO2019149718A5 publication Critical patent/JPWO2019149718A5/ja
Application granted granted Critical
Publication of JP7209003B2 publication Critical patent/JP7209003B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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/481Diagnostic techniques involving the use of contrast agents
    • 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
    • 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/30196Human being; Person
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/444Low dose acquisition or reduction of radiation dose

Landscapes

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

Description

以下は、一般にイメージングに関連し、より具体的には、定量分析が向上された低放射線量コンピュータ断層撮影灌流(CTP)に関する。
コンピュータ断層撮影灌流(CTP)イメージングは、血管及び/又は臓器といった関心脈管組織を通る投与された造影剤(例えばヨウ素)の通過を捕捉するイメージング技術である。一般にCTPイメージングでは、造影剤ボーラスが患者に投与され、関心脈管組織を含む患者の関心領域がスキャンされる。造影剤は、造影剤が脈管組織を通って流れるときに、関心脈管組織におけるX線密度を一時的に増加させる。典型的な灌流スキャンには、例えば造影剤の到着、取り込み及びウォッシュアウトを含む複数の時間間隔にわたって同じ領域のデータを収集することが含まれる。
収集されたデータを分析して、例えば造影剤の動態に基づいて、関心脈管組織の灌流状態を決定することができる。心臓用途では、これには、時間の経過に伴う心筋内の造影剤分布の定量化が含まれうる。このような分析には、時間減衰曲線、血流、血液量、平均通過時間、最大上昇勾配、ピークまでの時間等といった関心脈管組織の様々な灌流関連情報を決定することが含まれうる。この情報を使用して、虚血組織を特定したり、不可逆的に損傷した(又は壊死した)組織と潜在的に可逆的に損傷した(又はリスクのある)組織とを区別したりすることができる。
しかし、放射X線は、電離放射線であり、細胞を殺したり、損傷したり、がんのリスクを高めたりすることがある。CTPイメージングでは、X線が30秒以上継続してオンにされ、これは患者X線量が高いと見なされる。X線量を減らす1つの方法は、患者が連続的にスキャンされない時間領域でスパーススキャンを行うことである。残念ながら、スキャン間の時間間隔が長いと、定量測定の精度が低下する可能性がある。別のアプローチは、高度反復再構成(IR)を行うことである。しかし、IRは集中的な計算を必要とし、全体的な再構成時間を増加させる。更に、X線量を減らすとノイズが増加し、画質が低下する。
本明細書に説明する態様は、上記問題及びその他に対処する。
一態様では、コンピュータ断層撮影スキャナは、放射X線を放出する放射線源と、放射X線を検出し、それを示す信号を生成する検出器アレイと、信号を再構成して、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを生成する再構成器とを含む。コンピュータ断層撮影スキャナは更に、潅流データ強化モジュールのトレーニングされたニューラルネットワークを使用してシーケンシャルスパースタイムライン灌流ボリュメトリック画像データを処理してシーケンシャルデンスタイムライン灌流ボリュメトリック画像データを生成するプロセッサを含む。
別の態様では、コンピュータ断層撮影スキャナは、放射X線を放出する放射線源と、放射X線を検出し、それを示す信号を生成する検出器アレイと、信号を再構成して、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを生成する再構成器とを含む。コンピュータ断層撮影スキャナは更に、ニューラルネットワークをトレーニングして、トレーニング用シーケンシャルスパースタイムライン灌流ボリュメトリック画像データ及びトレーニング用シーケンシャルデンスタイムライン灌流ボリュメトリック画像データからシーケンシャルデンスタイムライン灌流ボリュメトリック画像データを予測するプロセッサを含む。
別の態様では、コンピュータ可読記憶媒体は、コンピューティングシステムのコンピュータプロセッサによって実行されると、コンピュータプロセッサに、トレーニングされたニューラルネットワークを使用して、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを処理して、シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを生成させるコンピュータ可読命令でエンコードされる。
当業者であれば、添付の説明を読んで理解すると、本願の更に他の態様を認識するであろう。
本発明は、様々なコンポーネント及びコンポーネントの構成、また、様々なステップ及びステップの構成の形を取ってよい。図面は、好適な実施形態を例示するに過ぎず、本発明を限定するものと解釈されるべきではない。
図1は、シーケンシャルスパースタイムラインCTPボリュメトリック画像データからシーケンシャルデンスタイムラインCTPボリュメトリック画像データを予測するためのニューラルネットワークを有するCTPデータ強化モジュールを含む例示的なCTスキャナを概略的に示す。 図2は、トレーニング用シーケンシャルスパースタイムラインCTPボリュメトリック画像データのセット及びトレーニング用シーケンシャルデンスタイムラインCTPボリュメトリック画像データのセットを用いたニューラルネットワークのトレーニングの一例を概略的に示す。 図3は、図2に示すようにトレーニングすることができる例示的なニューラルネットワークを概略的に示す。 図4は、図2に示すようにトレーニングされた図3のニューラルネットワークを使用するCTスキャナを概略的に示す。 図5は、トレーニング用シーケンシャルデンスタイムラインCTPボリュメトリック画像データのセットを用いて生成されるCTPパラメータのセットを用いたニューラルネットワークの更なるトレーニングの一例を概略的に示す。 図6は、図5に示すようにトレーニングすることができる例示的なニューラルネットワークを概略的に示す。 図7は、図5に示すようにトレーニングされた図6のニューラルネットワークを使用するCTスキャナを概略的に示す。 図8は、本明細書の一実施形態による例示的な方法を示す。 図9は、本明細書の一実施形態による別の例示的な方法を示す。
以下は、定量分析が向上される低放射線量CTPのアプローチについて説明する。CTPイメージングは、高X線量の手順と見なされている。X線量は、例えばX線が連続的にオンではない時間領域でスパーススキャンを行うことによって減らすことができる。しかし、これにより時間分解能が低下するため、定量測定の精度も低下する。線形補間は、時間スパース画像から時間デンス画像を生成することができる。しかし、造影剤は血管内を一定の速度で流れないので、造影剤の分布は時間に対して線形ではない場合がある。結果として、線形補間された画像は恣意的であり、情報を失う傾向がある。
本明細書に説明する実施形態では、ニューラルネットワークをトレーニングして、血管/組織の構造、血管/組織内の空間強度分布及び造影剤が異なる構造及び異なる強度分布を有する血管/組織内をどのように流れるかを学習する。次に、トレーニングされたニューラルネットワークは、スパースタイムラインCTP画像のセットを処理し、時間分解能が増加したことで定量的測定の精度も増加したデンスタイムラインCTP画像のセットを生成することができる。
図1は、例示的なコンピュータ断層撮影(CT)スキャナ102を概略的に示す。CTスキャナ102は、固定ガントリ104と、回転ガントリ106とを含み、回転ガントリ106は、固定ガントリ104によって回転可能に支持され、長手方向又はz軸110を中心に検査領域108及びその中の物体又は被験者の一部の周りを回転する。
X線管といった放射線源112が、回転ガントリ106によって支持され、検査領域104の周りを回転ガントリ106と共に回転する。放射線源112は、コリメートされて、検査領域108を横断する概ね扇形、くさび形又は円錐形の放射X線ビームを形成する放射X線を放出する。
線源コントローラ(SC)114が、放射線源112を制御する。これには、放射線源112を作動させて、時系列灌流イメージング中に放射線を放出させて、時間フレームデータを収集することが含まれる。このような作動は、放射線源112の各回転の全体又は一部を含む放射線源112の1回以上の回転についてのイメージング検査中に時間的に連続的又は時間的に断続的であってよい。
放射線感応検出器アレイ116が、検査領域108を挟んで放射線源112の反対側にある弧に対応する。検出器アレイ116は、z軸110に沿って互いに対して配置される1つ以上の検出器列を含み、検査領域108を横断する放射線を検出する。検出器アレイ116は、検出された放射線を示す投影データ(線積分)を生成する。
再構成器118が、投影データを再構成して、ボリュメトリック画像データを生成する。CTP手順では、これには、以下で更に詳細に説明するように、シーケンシャル収集の投影データを再構成することと、時系列CTPボリュメトリック画像データを生成することとが含まれる。適切な再構成アルゴリズムの例には、フィルタ補正逆投影、統計的、反復的なスパースサンプリング及び/又は他の再構成アルゴリズムが含まれる。
カウチといった被験者支持体120が、検査領域108内の物体(例えばファントム)又は被験者122を支持する。被験者支持体120は、物体又は被験者122のローディング、スキャン及び/又はアンローディングのために物体又は被験者122を検査領域108に対して誘導するように、イメージング手順の実行と連携して移動可能である。
インジェクタ124が、灌流スキャン中にスキャンされる物体又は被験者122に、1つ以上の造影剤といった材料を注入又は投与する。造影剤(例えばヨウ素ベース)は、更に又は或いは、臨床医等によって手動で投与されてもよい。造影剤が物体又は被験者122に手動で投与される場合、インジェクタ124を省略することができる。
オペレーターコンソール126が、ディスプレイモニタ、フィルマ等といった人間が読み取り可能な出力デバイス128と、キーボード、マウス等といった入力デバイス130とを含む。コンソール126は更に、プロセッサ132(例えば中央処理演算ユニット(CPU)、マイクロプロセッサ等)及び物理メモリといったコンピュータ可読記憶媒体134(一時的な媒体を除く)を含む。図示される実施形態では、コンピュータ可読記憶媒体134は、CTPデータ強化モジュール136を含む。変形例では、CTPデータ強化モジュール136は、スキャナ102から離れており、スキャナ102の一部ではない別のコンピューティングシステムの一部である。
図示されるCTPデータ強化モジュール136は、CTPイメージングスキャン中にCTスキャナ102(又は他のスキャナ)によって生成されるシーケンシャルスパースタイムラインCTPボリュメトリック画像データからシーケンシャルデンスタイムラインCTPボリュメトリック画像データを予測するためのコンピュータ実行可能命令を含む。以下でより詳細に説明するように、コンピュータ実行可能命令は、シーケンシャルデンスタイムラインCTPボリュメトリック画像データを予測するためにトレーニングされたニューラルネットワークを含む。一例では、これにより、例えば時間的に連続する(デンス(dense))CTPイメージングスキャンからシーケンシャルデンスタイムラインCTPボリュメトリック画像データが生成されたかのように、時間的に断続的な(スパース(sparse))CTPスキャンを行い、シーケンシャルデンスタイムラインCTPボリュメトリック画像データを生成することができる。したがって、患者のX線量は(時間的に連続するCTPイメージングスキャンと比較して)減少するが、正確なCTPパラメータ測定のためのシーケンシャルデンスタイムラインCTPボリュメトリック画像データが生成される。
コンピュータ可読記憶媒体134は更に、CTPデータ分析モジュール138を含み、これは、CTPデータ強化モジュール136によって生成されたシーケンシャルデンスタイムラインCTPボリュメトリック画像データを処理して、時間減衰曲線、血流、血液量、平均通過時間、最大上昇勾配、ピークまでの時間及び/又は他の灌流パラメータといった1つ以上のCTPパラメータを決定するためのコンピュータ実行可能命令を含む。変形例では、CTPデータ強化モジュール136のニューラルネットワークは更に、トレーニング用シーケンシャルデンスタイムラインCTPボリュメトリック画像データから生成された1つ以上のトレーニング用CTPパラメータでトレーニングされ、これらのCTPパラメータのうちの1つ又は複数を、予測されたシーケンシャルデンスタイムラインCTPボリュメトリック画像データから予測する。この変形例では、CTPデータ分析モジュール138を省略することができる。
図2は、シーケンシャルスパースタイムラインCTPボリュメトリック画像データのトレーニングセット及びシーケンシャルデンスタイムラインCTPボリュメトリック画像データのトレーニングセットを用いるニューラルネットワーク202の例示的なトレーニングを概略的に示す。一例では、トレーニング用スパースボリュメトリック画像及びデンスボリュメトリック画像は、同じデンスタイムラインCTPスキャンから取得される。例えばシーケンシャルスパースタイムラインCTPボリュメトリック画像データのトレーニングセットは、シーケンシャルデンスタイムラインCTPボリュメトリック画像データのトレーニングセットから、例えばそのサブセットとして(例えば時間的に1つおきのセット等)導出することができる。或いは、トレーニング用スパースCTPボリュメトリック画像データ及びデンスCTPボリュメトリック画像データは、別々のCTPイメージングスキャンを通じて取得される。
ニューラルネットワーク202の出力は、目的関数204を用いて評価される。この例では、目的関数204は、シーケンシャルスパースタイムラインCTPボリュメトリック画像データのトレーニングセットとシーケンシャルデンスタイムラインCTPボリュメトリック画像データのトレーニングセットとの間のエラーを、例えばそれらの間の平均平方差に基づいて及び/又は別の方法で最小にする数学関数を含む。ニューラルネットワークのトレーニング可能なパラメータが更新され、幾つかの反復を通じて収束する。この例では、パラメータは、エラーが所定の閾値を下回るまで更新される。他の停止基準には、所定の反復回数、所定の継続時間等が含まれる。
一例では、目的関数204は、例えば式1:
Figure 0007209003000001
に示すように、費用関数L(w)として実装される。ここで、S1は、トレーニング用スパースタイムラインCTPボリュメトリック画像データのj番目のセットであり、D1は、トレーニング用デンスタイムラインCTPボリュメトリック画像データのj番目のセットであり、Tは、トレーニング可能なパラメータwを有する変換を表す。正則化項
Figure 0007209003000002
は、小さいパラメータwに向けて式を駆動し、λは、ハイパーパラメータを表し、kは、ニューラルネットワークの層を表す。
図3は、図2に適したニューラルネットワーク202の一例を概略的に示す。入力にはシーケンシャルスパースタイムラインCTPボリュメトリック画像データ(スパースデータ)が含まれ、出力にはシーケンシャルデンスタイムラインCTPボリュメトリック画像データ(デンスデータ)が含まれる。
この例では、ニューラルネットワーク202は、3次元(3D)畳み込み層302、正規化線形ユニット304及びバッチ正規化306を有する3D畳み込みニューラルネットワークであり、3D畳み込み層/正規化線形ユニット/バッチ正規化中間層のN個(Nは1以上の正数)のセット308、…、308を有する。一例では、各畳み込み層で使用される畳み込みカーネルは3Dカーネルである。適切なサイズのカーネルの例には、3×3×3、5×5×5、7×7×7及び/又は他のサイズのカーネルが含まれる。層数、カーネルサイズ、チャネル数、学習率等を調整することができる。
一般に、各畳み込み層は、入力に畳み込み演算を適用する畳み込みニューロンを含む。一例では、これには、各ボリュームの幅及び高さにわたってカーネルのセットを畳み込み、フィルタのエントリとボリュームとのドット積を計算することが含まれる。各正規化線形ユニットは、その入力に非線形活性化関数を適用する。これにより、ネットワーク全体の非線形特性が増加する。バッチ正規化は、出力が正規化線形ユニットによって処理される前に、3D畳み込み層の出力を正規化する。これにより、内部共変量シフトが減少し、より高い学習率の使用につながる可能性がある。
畳み込みネットワークの例については、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年に更に説明されている。
図4は、CTPデータ強化モジュール136が図2でトレーニングされた図3のニューラルネットワーク202を含む一実施形態を概略的に示す。トレーニングされたニューラルネットワーク202は、CTスキャナ102(又は他のスキャナ)によって生成されるシーケンシャルスパースタイムラインCTPボリュメトリック画像データを受け取り、シーケンシャルデンスタイムラインCTPボリュメトリック画像データを出力する。CTPデータ分析モジュール138は、この出力を処理し、時間減衰曲線、血流、血液量、平均通過時間、最大上昇勾配、ピークまでの時間及び/又は他の灌流パラメータといった1つ以上のCTPパラメータを生成する。
図5は、ニューラルネットワーク202が、トレーニング用シーケンシャルデンスタイムラインCTPボリュメトリック画像データから生成されるトレーニング用灌流パラメータで更にトレーニングされている図2の変形例を概略的に示す。この例では、目的関数204は、式2:
Figure 0007209003000003
に示す費用関数といった費用関数L(w,wl)として実装される。ここで、T1は、トレーニング可能なパラメータw及び灌流パラメータNNトレーニング可能なパラメータw1を有する変換を表し、Bfmは、デンスタイムラインのCTPボリュメトリック画像データから生成されるトレーニング用灌流パラメータを表す。この例では、ニューラルネットワーク202は、真のデンスタイムラインCTPボリュメトリック画像データとできるだけ近いデンスタイムラインCTPボリュメトリック画像データを、真のデンスタイムラインCTPボリュメトリック画像データから計算される真の灌流パラメータとできるだけ近い灌流パラメータと同時に予測するようにトレーニングされる。
図6は、図5に示す構成のニューラルネットワーク202に適したニューラルネットワークの一例を概略的に示す。この例では、ニューラルネットワーク202は更に、CTPパラメータ層602を含む。CTPパラメータ層602は、出力されたシーケンシャルデンスタイムラインCTPボリュメトリック画像データを同時に受け取り、そこから1つ以上のCTPパラメータを生成する。
図7は、CTPデータ強化モジュール136が図5でトレーニングされた図6のニューラルネットワーク202を含む一実施形態を概略的に示す。トレーニングされたニューラルネットワーク202は、CTスキャナ102(又は他のスキャナ)によって生成されるシーケンシャルスパースタイムラインCTPボリュメトリック画像データを受け取り、シーケンシャルデンスタイムラインCTPボリュメトリック画像データ及び1つ以上のCTPパラメータの両方を出力する。
図8は、図2~図4に説明する実施形態による例示的な方法を示す。
この方法のステップの順序は限定ではないことを理解されたい。したがって、本明細書では他の順序も考えられる。更に、1つ以上のステップが省略されても及び/又は1つ以上の追加のステップが含まれてもよい。
ステップ802において、本明細書及び/又は別の方法で説明するように、トレーニング用スパース灌流画像のセットを受け取る。
ステップ804において、本明細書及び/又は別の方法で説明するように、トレーニング用デンス灌流画像のセットを受け取る。
ステップ806において、本明細書及び/又は別の方法で説明するように、トレーニング用スパース灌流画像のセット及びトレーニング用デンス灌流画像のセットを使用して、ニューラルネットワークをトレーニングする。
ステップ808において、本明細書及び/又は別の方法で説明するように、スパース灌流画像を受け取る。
ステップ810において、スパース灌流画像を、トレーニングされたニューラルネットワークで処理してデンス灌流画像を生成する。
ステップ812において、デンス灌流画像を処理してCTPパラメータを生成する。
図9は、図5~図7に説明する実施形態による例示的な方法を示す。
この方法のステップの順序は限定ではないことを理解されたい。したがって、本明細書では他の順序も考えられる。更に、1つ以上のステップが省略されても及び/又は1つ以上の追加のステップが含まれてもよい。
ステップ902において、本明細書及び/又は別の方法で説明するように、トレーニング用スパース灌流画像のセットを受け取る。
ステップ904において、本明細書及び/又は別の方法で説明するように、トレーニング用デンス灌流画像のセット及びそこから生成されるCTPパラメータのセットを受け取る。
ステップ906において、本明細書及び/又は別の方法で説明するように、トレーニング用スパース灌流画像のセット、トレーニング用デンス灌流画像のセット及びCTPパラメータのセットを使用してニューラルネットワークをトレーニングする。
ステップ908において、本明細書及び/又は別の方法で説明するように、スパース灌流画像を受け取る。
ステップ910において、スパース灌流画像を、トレーニングされたニューラルネットワークで処理してデンス灌流画像及びそこからCTPパラメータを生成する。
上記は、コンピュータ可読記憶媒体にエンコードされ又は埋め込まれ、コンピュータプロセッサによって実行されると当該プロセッサに、説明したステップを行わせるコンピュータ可読命令によって実施することができる。更に又は或いは、コンピュータ可読命令の少なくとも1つは、信号、搬送波又はコンピュータ可読記憶媒体ではない他の一時的媒体によって運ばれる。
本発明は、図面及び上記説明において詳細に例示され、説明されたが、当該例示及び説明は、例示的に見なされるべきであり、限定的に見なされるべきではない。本発明は、開示した実施形態に限定されない。開示した実施形態の他の変形態様は、図面、開示内容及び添付の請求項の検討から、請求項に係る発明を実施する当業者によって理解され、実施される。
請求項において、「含む」との用語は、他の要素又はステップを排除するものではなく、また、単数形も、複数形を排除するものではない。単一のプロセッサ又は他のユニットが、請求項に引用される幾つかのアイテムの機能を果たしてもよい。特定の手段が相互に異なる従属請求項に記載されることだけで、これらの手段の組み合わせを有利に使用することができないことを示すものではない。
コンピュータプログラムは、他のハードウェアと共に又はその一部として供給される光学記憶媒体又はソリッドステート媒体といった適切な媒体上に記憶及び/又は分散されてもよいが、インターネット又は他の有線若しくは無線通信システムを介するといった他の形式で分散されてもよい。請求項における任意の参照符号は、範囲を限定するものと解釈されるべきではない。

Claims (15)

  1. 放射X線を放出する放射線源と、
    放射X線を検出し、それを示す信号を生成する検出器アレイと、
    前記信号を再構成して、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを生成する再構成器と、
    流データ強化モジュールのトレーニングされたニューラルネットワークを使用して前記シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを処理してシーケンシャルデンスタイムライン灌流ボリュメトリック画像データを生成するプロセッサと、
    を含み、
    前記プロセッサは、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データのトレーニングセット、シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセット、及び、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットから生成される灌流パラメータのトレーニングセットを用いて前記ニューラルネットワークをトレーニングする、コンピュータ断層撮影スキャナ。
  2. 前記プロセッサは更に、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを処理して1つ以上の灌流パラメータを生成する、請求項1に記載のコンピュータ断層撮影スキャナ。
  3. 前記プロセッサは、前記トレーニングされたニューラルネットワークを使用して、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データ及び前記1つ以上の灌流パラメータの両方を生成する、請求項2に記載のコンピュータ断層撮影スキャナ。
  4. 前記プロセッサは、異なるニューラルネットワークを使用して、前記1つ以上の灌流パラメータを生成する、請求項2に記載のコンピュータ断層撮影スキャナ。
  5. 前記プロセッサは、灌流データ分析モジュールを使用して前記1つ以上の灌流パラメータを生成する、請求項2に記載のコンピュータ断層撮影スキャナ。
  6. 前記シーケンシャルスパースタイムライン灌流ボリュメトリック画像データのトレーニングセットは、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットのサブセットである、請求項1から5のいずれか一項に記載のコンピュータ断層撮影スキャナ。
  7. 放射X線を放出する放射線源と、
    放射X線を検出し、それを示す信号を生成する検出器アレイと、
    前記信号を再構成して、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを生成する再構成器と、
    シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを予測するニューラルネットワークを、トレーニング用シーケンシャルスパースタイムライン灌流ボリュメトリック画像データ、トレーニング用シーケンシャルデンスタイムライン灌流ボリュメトリック画像データ、及び、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットから生成される灌流パラメータのトレーニングセットを用いてトレーニングするプロセッサと、
    を含む、コンピュータ断層撮影スキャナ。
  8. 前記シーケンシャルスパースタイムライン灌流ボリュメトリック画像データのトレーニングセットは、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットのサブセットである、請求項7に記載のコンピュータ断層撮影スキャナ。
  9. 前記プロセッサは、費用関数の出力が、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データの生成されたセットと前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットとの間の差を最小にするまで前記ニューラルネットワークをトレーニングする、請求項7又は8に記載のコンピュータ断層撮影スキャナ。
  10. 前記プロセッサは、費用関数の出力が、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データの生成されたセットと前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットとの間の差、及び、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データの生成されたセットから生成される灌流パラメータのセットと前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットから生成される灌流パラメータのセットとの間の差を最小にするまで前記ニューラルネットワークをトレーニングする、請求項7に記載のコンピュータ断層撮影スキャナ。
  11. 前記プロセッサは、前記トレーニングされたニューラルネットワークを用いて前記シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを処理して、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを生成する、請求項7から10のいずれか一項に記載のコンピュータ断層撮影スキャナ。
  12. 前記プロセッサは更に、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを処理して1つ以上の灌流パラメータを生成する、請求項11に記載のコンピュータ断層撮影スキャナ。
  13. 前記プロセッサは、前記トレーニングされたニューラルネットワークを用いて前記シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを処理して、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データ及び1つ以上の灌流パラメータを生成する、請求項7から12のいずれか一項に記載のコンピュータ断層撮影スキャナ。
  14. コンピューティングシステムのプロセッサによって実行されると、前記プロセッサに、
    トレーニングされたニューラルネットワークを使用して、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データを処理して、シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを生成させ、前記ニューラルネットワークは、シーケンシャルスパースタイムライン灌流ボリュメトリック画像データのトレーニングセット、シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセット、及び、前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データのトレーニングセットから生成される灌流パラメータのトレーニングセットを用いてトレーニングされる、コンピュータ可読命令でエンコードされた、コンピュータ可読記憶媒体。
  15. 前記コンピュータ可読命令を実行することにより、前記プロセッサに更に、
    前記シーケンシャルデンスタイムライン灌流ボリュメトリック画像データを処理して、1つ以上の灌流パラメータを生成させる、請求項14に記載のコンピュータ可読記憶媒体。
JP2020541568A 2018-02-01 2019-01-30 定量分析が向上された低放射線量コンピュータ断層撮影灌流(ctp) Active JP7209003B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201862624871P 2018-02-01 2018-02-01
US62/624,871 2018-02-01
PCT/EP2019/052168 WO2019149718A1 (en) 2018-02-01 2019-01-30 Low radiation dose computed tomography perfusion (ctp) with improved quantitative analysis

Publications (3)

Publication Number Publication Date
JP2021511881A JP2021511881A (ja) 2021-05-13
JPWO2019149718A5 JPWO2019149718A5 (ja) 2022-02-03
JP7209003B2 true JP7209003B2 (ja) 2023-01-19

Family

ID=65243577

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020541568A Active JP7209003B2 (ja) 2018-02-01 2019-01-30 定量分析が向上された低放射線量コンピュータ断層撮影灌流(ctp)

Country Status (5)

Country Link
US (1) US20210085273A1 (ja)
EP (1) EP3746986B1 (ja)
JP (1) JP7209003B2 (ja)
CN (1) CN111670461B (ja)
WO (1) WO2019149718A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112419303B (zh) * 2020-12-09 2023-08-15 上海联影医疗科技股份有限公司 神经网络训练方法、系统、可读存储介质和设备
CN113034438B (zh) * 2021-03-01 2021-11-23 北京安德医智科技有限公司 一种ct灌注参数预测方法及装置、电子设备和存储介质
CN114332043B (zh) * 2021-12-31 2022-12-13 深圳市铱硙医疗科技有限公司 基于ct灌注影像的侧枝循环测量方法、装置、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070092055A1 (en) 2005-10-14 2007-04-26 Vives Pau M Dynamic computed tomography method and apparatus with temporal interpolation of data in perfusion studies
JP2008136800A (ja) 2006-11-08 2008-06-19 Toshiba Corp X線診断装置及び画像処理装置
JP2014144372A (ja) 2006-11-08 2014-08-14 Toshiba Corp X線診断装置及び画像処理装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6888914B2 (en) * 2002-11-26 2005-05-03 General Electric Company Methods and apparatus for computing volumetric perfusion
US9107581B2 (en) * 2003-04-15 2015-08-18 Koninklijke Philips N.V. Elastography device and method for determining and imaging of mechanical and elastic parameters of an examination object
US7218702B2 (en) * 2004-05-10 2007-05-15 Wisconsin Alumni Research Foundation X-ray system for use in image guided procedures
JP5911610B2 (ja) * 2012-03-06 2016-04-27 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 介入的x線かん流画像化のための定期的造影剤注射とハーモニクスの分析
WO2016161308A1 (en) * 2015-04-02 2016-10-06 Heartflow, Inc. Systems and methods for predicting perfusion deficits from physiological, anatomical, and patient characteristics
WO2017192629A1 (en) * 2016-05-02 2017-11-09 The Regents Of The University Of California System and method for estimating perfusion parameters using medical imaging
EP3412208B1 (de) * 2018-01-29 2021-05-26 Siemens Healthcare GmbH Bereitstellen eines medizinischen bildes

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070092055A1 (en) 2005-10-14 2007-04-26 Vives Pau M Dynamic computed tomography method and apparatus with temporal interpolation of data in perfusion studies
JP2008136800A (ja) 2006-11-08 2008-06-19 Toshiba Corp X線診断装置及び画像処理装置
JP2014144372A (ja) 2006-11-08 2014-08-14 Toshiba Corp X線診断装置及び画像処理装置

Also Published As

Publication number Publication date
EP3746986B1 (en) 2023-07-05
US20210085273A1 (en) 2021-03-25
JP2021511881A (ja) 2021-05-13
CN111670461B (zh) 2023-10-17
EP3746986A1 (en) 2020-12-09
WO2019149718A1 (en) 2019-08-08
CN111670461A (zh) 2020-09-15

Similar Documents

Publication Publication Date Title
US9959640B2 (en) Iterative image reconstruction with a sharpness driven regularization parameter
US7526060B2 (en) Artifact correction
US20180330233A1 (en) Machine learning based scatter correction
US8055050B2 (en) Motion compensation in energy-sensitive computed tomography
JP2020168353A (ja) 医用装置及びプログラム
EP2476099B1 (en) Apparatus and method for processing projection data
JP7209003B2 (ja) 定量分析が向上された低放射線量コンピュータ断層撮影灌流(ctp)
US20080267480A1 (en) Iterative Image Reconstruction of a Moving Object From Projection Data
US10722178B2 (en) Method and apparatus for motion correction in CT imaging
US10420519B2 (en) Computed tomography system
US20140126685A1 (en) Angiographic examination method
CN111684486B (zh) 谱计算机断层摄影成像系统和计算机可读存储介质
US11986336B2 (en) Non-spectral computed tomography (CT) scanner configured to generate spectral volumetric image data
US7961839B2 (en) Advanced convergence for multiple iterative algorithm
US9858688B2 (en) Methods and systems for computed tomography motion compensation
Cant et al. Modeling blurring effects due to continuous gantry rotation: Application to region of interest tomography
EP3123446B1 (en) Image generation apparatus
US10984564B2 (en) Image noise estimation using alternating negation
US9996951B2 (en) Computerized tomographic image exposure and reconstruction method
TW201417769A (zh) 一種影像品質的改善處理方法及其造影系統
Pock Generative Regularizers in Computed Tomography
JP2024515698A (ja) コンピュータ断層撮影スキャナによって生成された投影領域データの処理

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

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221124

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230106

R150 Certificate of patent or registration of utility model

Ref document number: 7209003

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150