JP7535575B2 - 組織特性を予測、予想、および/または査定するためのシステム、方法、およびコンピュータプログラム製品 - Google Patents

組織特性を予測、予想、および/または査定するためのシステム、方法、およびコンピュータプログラム製品 Download PDF

Info

Publication number
JP7535575B2
JP7535575B2 JP2022517507A JP2022517507A JP7535575B2 JP 7535575 B2 JP7535575 B2 JP 7535575B2 JP 2022517507 A JP2022517507 A JP 2022517507A JP 2022517507 A JP2022517507 A JP 2022517507A JP 7535575 B2 JP7535575 B2 JP 7535575B2
Authority
JP
Japan
Prior art keywords
tissue
voxels
image
parameters
time points
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
JP2022517507A
Other languages
English (en)
Other versions
JP2022548930A (ja
JPWO2021053585A5 (ja
JP2022548930A5 (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.)
Bayer AG
Original Assignee
Bayer AG
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 Bayer AG filed Critical Bayer AG
Publication of JP2022548930A publication Critical patent/JP2022548930A/ja
Publication of JPWO2021053585A5 publication Critical patent/JPWO2021053585A5/ja
Publication of JP2022548930A5 publication Critical patent/JP2022548930A5/ja
Application granted granted Critical
Publication of JP7535575B2 publication Critical patent/JP7535575B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00Three-dimensional [3D] image rendering
    • G06T15/08Volume rendering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/42Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
    • A61B5/4222Evaluating particular parts, e.g. particular organs
    • A61B5/4244Evaluating particular parts, e.g. particular organs liver
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5601Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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 OR CALCULATING; 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/0464Convolutional networks [CNN, ConvNet]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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/0499Feedforward networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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 OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/09Supervised learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • G06T7/0016Biomedical image inspection using an image reference approach involving temporal comparison
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/62Extraction of image or video features relating to a temporal dimension, e.g. time-based feature extraction; Pattern tracking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56366Perfusion imaging
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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/10088Magnetic resonance imaging [MRI]
    • G06T2207/10096Dynamic contrast-enhanced magnetic resonance imaging [DCE-MRI]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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 OR CALCULATING; 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 OR CALCULATING; 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
    • G06T2207/30024Cell structures in vitro; Tissue sections in vitro
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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
    • G06T2207/30056Liver; Hepatic
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Quality & Reliability (AREA)
  • Gastroenterology & Hepatology (AREA)
  • Signal Processing (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Computer Graphics (AREA)
  • Multimedia (AREA)
  • Endocrinology (AREA)
  • Physiology (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本開示は、概して、画像診断、ならびに、いくつかの非限定的な実施形態または態様において、医用画像の後期相および/または組織特性を、選択された相から予測、予想、および/または査定することに関する。
磁気共鳴画像法、略してMRIは、ヒトまたは動物の身体内の組織および器官の構造および機能を描写するための、特に医用診断において使用される撮像方法である。
MRIにおいて、検査対象における陽子の磁気モーメントは、基本磁場内に整列され、その結果、長手方向に沿って巨視的磁化が存在する。これは、続いて、高周波(HF)パルスの入射放射線によって静止位置から偏向される(励起)。励起状態から静止位置への戻り(緩和)または磁化ダイナミクスは、続いて、1つまたは複数のHF受信器コイルによって緩和信号として検出される。
空間符号化の場合、迅速に切り替えられた勾配磁場が、基本磁場に重ね合わせられる。捕捉された緩和信号、または検出され空間的に分解されたMRIデータは、最初、空間周波数空間内に生データとして存在し、その後のフーリエ変換により実空間(画像空間)へ変換され得る。
ネイティブMRIの場合、組織コントラストは、異なる緩和時間(T1およびT2)ならびに陽子密度によって生成される。
T1緩和は、縦磁化のその平衡状態への遷移を説明し、T1は、共鳴励起の前に63.21%の平衡磁化に到達するために必要とされる時間である。それは、縦緩和時間またはスピン-格子緩和時間とも呼ばれる。
類似して、T2緩和は、横方向磁化のその平衡状態への遷移を説明する。
MRI造影剤は、造影剤を取り込む構造体の緩和時間を変更することによって、それらの作用を進展させる。常磁性体および超常磁性体という2つの物質グループの区別がなされ得る。両方の物質グループは、個々の原子または分子の周りに磁場を誘起する不対電子を有する。
超常磁性造影剤は、T2の顕著な短縮をもたらす一方、常磁性造影剤は主に、T1の短縮をもたらす。T1時間の短縮は、T1重み付けMRI画像内の信号強度の増加をもたらし、T2時間の短縮は、T2重み付けMRI画像内の信号強度の減少をもたらす。
造影剤自体が信号を出すことはなく、代わりに、その周辺の水素陽子の信号強度にだけ影響を及ぼすため、上記造影剤の作用は間接的である。
細胞外、細胞内、および血管内造影剤は、組織内のそれらの広がりのパターンに従って区別され得る。
超常磁性造影剤の例は、酸化鉄ナノ粒子(SPIO:超常磁性酸化鉄)である。超常磁性造影剤は、一般的に、血管内剤として使用される。
常磁性造影剤の例は、ガドペンテト酸ジメグルミン(商品名:Magnevist(登録商標)他)、ガドベン酸メグルミン(商品名:Multihance(登録商標))、ガドテル酸(Dotarem(登録商標)、Dotagita(登録商標)、Cyclolux(登録商標))、ガドジアミド(Omniscan(登録商標))、ガドテリドール(ProHance(登録商標))、およびガドブトロール(Gadovist(登録商標))などのカドリニウムキレートである。
細胞外、細胞内、および血管内造影剤は、組織内のそれらの広がりのパターンに従って区別され得る。
CTでは、造影剤は、CT機によって作成されるX線放射線を、X線放射線が患者を通過するときに吸収する。分子撮像の一形態である核医学では、造影剤は、撮像機器によって検出される信号を生成するために減衰する放射性原子を含む。過分極MRIでは、大半の核は、磁場と整列され、撮像のために患者内へ注入される。
ガドキセト酸に基づいた造影剤は、肝細胞、肝実質細胞による明確な取り込み、機能性組織(実質組織)における濃縮、および健康な肝臓組織内の造影の強化を特徴とする。嚢胞、転移、および大部分の肝細胞がん腫の細胞は、もはや正常な肝細胞のようには機能せず、造影剤を取り込まないか、またはそれをほとんど取り込まず、強調して描写されず、結果として識別可能および位置特定可能である。
ガドキセト酸に基づいた造影剤の例は、米国特許第6,039,931A号に説明されており、それらは、例えば、商品名Primovist(登録商標)またはEovist(登録商標)の下で市販されている。
Primovist(登録商標)/Eovist(登録商標)の造影強化効果は、安定したガドリニウム錯体Gd-EOB-DTPA(ガドリニウムエトキシベンジルジエチレントリアミンペンタ酢酸)によって仲介される。DTPAは、常磁性ガドリニウムイオンと共に、極端に高い熱力学的安定性を有する錯体を形成する。エトキシベンジルラジカル(EOB)は、造影剤の肝胆道取り込みの仲介物である。
Primovist(登録商標)は、肝臓内の腫瘍の検出のために使用され得る。健康な肝臓組織への血液供給は、主に、門脈(vena portae)により達成される一方、肝臓動脈(arteria hepatica)は、大部分の原発腫瘍を供給する。造影剤のボーラスの静脈内注射後、相応して、健康な肝実質の信号立ち上がりと腫瘍の信号立ち上がりとの間の時間遅延を観察することが可能である。悪性腫瘍のほかにも、肝臓内で頻繁に見られるのは、嚢胞、血管腫、および限局性結節性過形成(FNH)などの良性病変である。治療の適切な計画は、これらが悪性腫瘍から区別されることが必要とされる。Primovist(登録商標)は、良性および悪性限局性肝臓病変の識別のために使用され得る。T1重み付けMRIによって、それは、上記病変の性質に関する情報を提供する。区別は、肝臓および腫瘍への異なる血液供給、ならびに造影強化の経時的プロファイルを使用することによって達成される。
Primovist(登録商標)によって達成される造影強化は、動的相(いわゆる動脈相、門脈相、および後期相を含む)、および肝実質細胞へのPrimovist(登録商標)の著しい取り込みがすでに発生している肝胆道相という少なくとも2つの相に分割され得る。
分布相の間にPrimovist(登録商標)によって達成される造影強化の場合、観察されるのは、病変の特徴付けのための情報を提供する典型的な灌流パターンである。血管新生を描写することにより、病変タイプを特徴付けること、および腫瘍と血管との空間的関係を決定することを助ける。
T1重み付けMRI画像の場合、Primovist(登録商標)は、注入(肝胆道相において)から10~20分後、健康な肝実質においては明確な信号強調をもたらすが、肝実質細胞を含まない、またはほんの少しだけ肝実質細胞を含む病変、例えば、転移、または中~低分化肝細胞がん(HCC)は、より暗い領域として現れる。
動的相および肝胆道相にわたる経時的な造影剤の広がりを追跡することは、限局性肝臓病変の検出および鑑別診断の可能性を十分に提供するが、検査は比較的長い時間スパンに及ぶ。上記時間スパンにわたって、患者による運動は、MRI画像内の運動アーチファクトを最小限にするために回避されなければならない。運動の長きにわたる制限は、患者にとっては不快となり得、実際に一貫して達成するのが困難であり得る。
本開示の非限定的な実施形態または態様は、例えば、走査時間の加速または短縮における、画像診断検査の改善を対象とする(例えば、動的造影強化磁気共鳴画像法(MRI)などを用いた限局性肝臓病変の検出および鑑別診断のため)。本開示の非限定的な実施形態または態様は、選択した相(例えば、初期相など)、例えば、肝胆道相の間の肝臓の相から、医用画像の後期相を予測、予想、および/または査定し、以て、データ画像取得時間を低減するための方法、システム、およびコンピュータプログラム製品を提供する。追加的に、または代替的に、パラメータは、後期画像を合成するために、取って代わり得る、および/または使用され得る撮像シーケンスからのデータ、例えば、組織の体積(例えば、ボクセルなど)の薬物動力学的パラメータを使用して決定され得る。
造影剤を用いた肝臓組織のMRIに関して主に説明されるが、本開示の非限定的な実施形態または態様は、それに限定されず、非限定的な実施形態または態様に従う方法、システム、およびコンピュータプログラム製品は、造影剤ありまたはなしで任意の画像診断法を使用して取得される任意の器官の医用画像の選択した相(例えば、初期相など)から医用画像の後期相を予測、予想、および/または査定し得る。
いくつかの非限定的な実施形態または態様によると、複数のMRI画像を受信するステップであって、MRI画像が、第1の時間スパンの間の検査領域を示す、ステップと、複数のMRI画像を予測モデルに供給するステップであって、予測モデルが、第1の時間スパンの間の検査領域を示すMRI画像に基づいて、第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を予測するように、教師あり学習によって訓練されている、ステップと、予測によって第2の時間スパンの間の検査領域を示す1つまたは複数の予測MRI画像を生成するステップと、1つもしくは複数の予測MRI画像を表示および/もしくは出力する、ならびに/または1つもしくは複数の予測MRI画像をデータ記憶媒体に記憶するステップと、を含む方法が提供される。
いくつかの非限定的な実施形態または態様によると、受信装置、制御および計算装置、ならびに出力装置を備えるシステムであって、制御および計算装置は、受信装置に複数のMRI画像を受信するように促すように構成され、受信したMRI画像は、第1の時間スパンの間の検査領域を示し、制御および計算装置は、受信したMRI画像に基づいて1つまたは複数のMRI画像を予測するように構成され、1つまたは複数の予測MRI画像は、第2の時間スパンの間の検査領域を示し、制御および計算装置は、出力装置に、1つまたは複数の予測MRI画像を表示するように、それらを出力するように、またはそれらをデータ記憶媒体に記憶するように促すように構成される、システムが提供される。
いくつかの非限定的な実施形態または態様によると、コンピュータのメモリに読み込まれ得るコンピュータプログラムを備えるコンピュータプログラム製品であって、コンピュータプログラムは、コンピュータに、以下のステップ:複数のMRI画像を受信するステップであって、MRI画像が、第1の時間スパンの間の検査領域を示す、ステップと、受信したMRI画像を予測モデルに供給するステップであって、予測モデルが、第1の時間スパンの間の検査領域を示すMRI画像に基づいて、第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を予測するように、教師あり学習によって訓練されている、ステップと、予測モデルからの出力として、第2の時間スパンの間の検査領域を示す1つまたは複数の予測MRI画像を受信するステップと、1つもしくは複数の予測MRI画像を表示および/もしくは出力する、ならびに/または1つもしくは複数の予測MRI画像をデータ記憶媒体に記憶するステップとを実行することを促す、コンピュータプログラム製品が提供される。
いくつかの非限定的な実施形態または態様によると、MRI法における造影剤の使用が提供され、本MRI法は、以下のステップ:造影剤を投与するステップであって、造影剤が検査領域に広がる、ステップと、第1の時間スパンの間の検査領域の複数のMRI画像を生成するステップと、生成したMRI画像を予測モデルに供給するステップであって、予測モデルが、第1の時間スパンの間の検査領域を示すMRI画像に基づいて、第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を予測するように、教師あり学習によって訓練されている、ステップと、予測モデルからの出力として、第2の時間スパンの間の検査領域を示す1つまたは複数の予測MRI画像を受信するステップと、1つもしくは複数の予測MRI画像を表示および/もしくは出力する、ならびに/または1つもしくは複数の予測MRI画像をデータ記憶媒体に記憶するステップとを含む。
いくつかの非限定的な実施形態または態様によると、MRI法における使用のための造影剤が提供され、本MRI法は、以下のステップ:造影剤を投与するステップであって、造影剤が検査領域に広がる、ステップと、第1の時間スパンの間の検査領域の複数のMRI画像を生成するステップと、生成したMRI画像を予測モデルに供給するステップであって、予測モデルが、第1の時間スパンの間の検査領域を示すMRI画像に基づいて、第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を予測するように、教師あり学習によって訓練されている、ステップと、予測モデルからの出力として、第2の時間スパンの間の検査領域を示す1つまたは複数の予測MRI画像を受信するステップと、1つもしくは複数の予測MRI画像を表示および/もしくは出力する、ならびに/または1つもしくは複数の予測MRI画像をデータ記憶媒体に記憶するステップとを含む。
本開示の実施形態または態様に従う造影剤およびコンピュータプログラム製品を備えるキットがさらに提供される。
本開示の非限定的な実施形態または態様は、実施形態の主題(方法、システム、コンピュータプログラム製品、使用、使用のための造影剤、キット)を区別することなく以下により具体的に解明される。逆に、以下の解明は、それらが発生する文脈(方法、システム、コンピュータプログラム製品、使用、使用のための造影剤、キット)とは関係なく、すべての実施形態のすべての主題に類似的に当てはまることが意図される。
ステップが、本説明または特許請求項においてある順序で記載される場合、これは、実施形態または態様が記載された順序に制限されることは必ずしも意味しない。逆に、ステップはまた、1つのステップが別のステップの上に積み上げられる(これは積み上げステップが続いて実行されることを絶対的に必要とする(しかしながら、これは個々の場合に明白である))ということがない限り、異なる順序で、さもなければ互いに並列して実行されることが想定される。故に、記載された順序は、好ましい実施形態であり得る。
本開示の非限定的な実施形態または態様は、MRI画像の生成における検査対象の検査の時間スパンを短縮する。いくつかの非限定的な実施形態または態様において、これは、検査対象の検査領域のMRI画像が第1の時間スパンにおいて測定されること(磁気共鳴測定)、および測定されたMRI画像が次いで、自己学習アルゴリズムの助けを借りて、第2の時間スパンにおける検査領域を示す1つまたは複数のMRI画像を予測するために使用されること、によって達成される。検査対象に対する実際の磁気共鳴測定は、故に、本開示に制限され、第2の時間スパンを包含しない。第1の時間スパンの間の検査領域を示すMRI画像は、第2の時間スパンの予測を可能にする情報を含む。
“検査対象”は、通常、生物、好ましくは哺乳動物、非常に具体的には、好ましくはヒトであり得る。検査領域は、検査対象の一部分、例えば、器官または器官の一部分であり得る。非限定的な実施形態または態様において、検査領域は、哺乳動物(例えば、ヒト)の肝臓または肝臓の一部分である。
画像体積とも呼ばれる“検査領域”(視野、FOV)は、特に、磁気共鳴画像において撮像される体積であり得る。一般的に、3次元視野は、1つまたは複数の体積要素またはボクセルからなると言うことができる。2次元FOVが検討されている場合、1つもしくは複数のボクセルまたは1つもしくは複数のピクセル(画素)からなると言うことができる。検査領域は、典型的には、例えば概観画像に対して、放射線科医によって規定され得る(ローカライザ)。検査領域が、代替的または追加的に、例えば、選択されたプロトコルに基づいて、自動的に規定され得ることは自明である。
検査領域は、基本磁場内へ導入される。検査領域は、MRI法の対象となり、これが、第1の時間スパンの間の検査領域を示す複数のMRI画像を生成する。
複数という用語は、少なくとも2つのMRI画像または測定値、好ましくは少なくとも3つ、非常に具体的には、好ましくは少なくとも4つのMRI画像または測定値が生成されることを意味する。
検査領域に広がる造影剤は、検査対象に投与される。造影剤は、好ましくは、体重適合様式で、ボーラスとして静脈内に投与される。
“造影剤”は、磁気共鳴測定におけるその存在が改変信号をもたらす物質または物質混合物を意味すると理解される。好ましくは、造影剤は、T1緩和時間および/またはT2緩和時間の短縮をもたらす。
好ましくは、造影剤は、例えば、Gd-EOB-DTPAまたはGd-BOPTAなどの肝胆道造影剤である。
特に好ましい実施形態において、造影剤は、造影強化活性物質として、ガドキセト酸またはガドキセト酸塩を有する物質または物質混合物である。非常に具体的に好ましいのは、ガドキセト酸のジナトリウム塩(Gd-EOB-DTPAジナトリウム)である。
好ましくは、第1の時間スパンは、造影剤の投与の前、または造影剤の投与と共に、開始する。造影剤なしに検査領域を示す1つまたは複数のMRI画像が生成されるときが有利であるが、それは、放射線科医がそのような画像内で検査対象の健康の状態に関する重要な情報を既に獲得することができるためである。例えば、放射線科医は、そのようなネイティブMRI画像内で出血を識別することができる。
第1の時間スパンは、好ましくは、検査領域内に分布する造影剤を包含する。好ましくは、第1の時間スパンは、検査対象の肝臓または肝臓の一部分の動的造影強化磁気共鳴断層撮影における動脈相および/または門脈相および/または後期相を包含する。記載の相は、例えば、以下の刊行物:J.Magn.Reson.Imaging,2012,35(3):492-511,doi:10.1002/jmri.22833; Clujul Medical,2015,Vol.88 no.4:438-448,DOI:10.15386/cjmed-414; Journal of Hepatology,2019,Vol.71:534-542, http://dx.doi.org/10.1016/j.jhep.2019.05.005)に規定および説明される。
好ましい実施形態において、第1の時間スパンは、検査対象の肝臓または肝臓の一部分のそのようなMRI画像が生成されるように選択され、
造影剤なしの検査領域を示し、
造影剤が動脈を介して検査領域内に広がる動脈相の間の検査領域を示し、
造影剤が門脈を介して検査領域に到達する門脈相の間の検査領域を示し、
動脈および静脈内の造影剤の濃度が低下し、血管外組織および/または肝細胞内の造影剤の濃度が上昇する後期相の間の検査領域を示す。
好ましくは、第1の時間スパンは、造影剤の投影前1分~1秒の時間スパン内で、または造影剤の投与と共に、開始し、造影剤の投与から2分~15分、好ましくは2分~13分、さらにより好ましくは3分~10分の時間スパンにわたって継続する。造影剤は非常に緩徐に腎臓および/または胆道で排出されるため、広範な時間スパンは、造影剤の投与後、最大で2時間以上におよび得る。
造影剤は異なる検査対象においては様々な速さで広がり得るため、第1の時間スパンはまた、検査領域の異なるエリアにおける造影剤の濃度により規定され得る。1つの可能性が図1に描写される。
図1は、肝臓動脈(A)、肝臓静脈(V)、および健康な肝細胞(P)における造影剤の濃度の経時的プロファイルを概略的に示す。濃度は、時間tの関数としての磁気共鳴測定における記載のエリア(肝臓動脈、肝臓静脈、肝細胞)内の信号強度Iの形態で描写される。静脈内ボーラス注射の際、造影剤の濃度は、第一に肝臓動脈(A)内で上昇する(破線曲線)。濃度は、最大を通過し、次いで降下する。肝臓静脈(V)内の濃度は、肝臓動脈内よりも緩徐に上昇し、後にその最大に到達する(点線曲線)。健康な肝細胞(P)内の造影剤の濃度は、緩徐に上昇し(連続曲線)、はるかに後の時間点(図1には描写されない)においてのみその最大に到達する。いくつかの特徴的な時間点が規定され得る。時間点TP0において、造影剤は、ボーラスとして静脈内に投与される。時間点TP1において、肝臓動脈内の造影剤の濃度(信号強度)は、その最大に到達する。時間点TP2において、肝臓動脈および肝臓静脈の信号強度の曲線は交差する。時間点TP3において、肝臓静脈内の造影剤の濃度(信号強度)は、その最大を通過する。時間点TP4において、肝臓動脈および肝細胞の信号強度の曲線は交差する。時間点T5において、肝臓動脈および肝臓静脈内の濃度は、それらがもはや測定可能な造影強化を引き起こさないレベルまで降下している。
好ましい実施形態において、第1の時間スパンは、少なくとも時間点TP0、TP1、TP2、TP3、およびTP4を包含する。
好ましい実施形態において、TP0~TP1の時間スパン、TP1~TP2の時間スパン、TP2~TP3の時間スパン、およびTP3~TP4の時間スパンにおけるすべての相の少なくともMRI画像が生成される(測定によって)。
TP0~TP1、TP1~TP2、TP2~TP3、TP3~TP4の時間スパンにおいて、1つまたは複数のMRI画像がその都度(測定によって)生成されることが想定される。1つまたは複数の時間スパンの間、MRI画像のシーケンスが(測定によって)生成されることも想定される。
シーケンスという用語は、時系列、すなわち、生成されるものが連続した時間点における検査領域を示す複数のMRI画像であることを意味する。
時間点は、各MRI画像に割り当てられるか、または時間点は、各MRI画像に割り当てられ得る。通常、この時間点は、MRI画像が生成されている時間点(絶対時間)である。当業者は、MRI画像の生成が特定の時間スパンを使用することを認識している。MRI画像に割り当てられ得るものは、例えば、取得の開始の時間点または取得の完了の時間点である。しかしながら、任意の時間点がMRI画像に割り当てられる(例えば、相対時間点)ことも想定される。
時間点に基づいて、MRI画像は、別のMRI画像に対して経時的に配置され得、このMRI画像の時間点に基づいて、MRI画像内に示される瞬間が、別のMRI画像内に示される瞬間より経時的に前に、または経時的に後に起きたかを確立することが可能である。
好ましくは、MRI画像は、一連および複数のMRI画像において経時的に順序付けされ、その結果として、検査領域の初期状態を示すMRI画像は、一連および複数のMRI画像において、検査領域の後期状態を示すMRI画像の前に配置される。
一連および複数のMRI画像内の互いの直後の2つのMRI画像の間の時間スパンは、好ましくは、一連および複数のMRI画像内の互いの直後のMRI画像のすべての対について同一であり、すなわち、MRI画像は、好ましくは、一定の取得速度で生成されたものである。
第1の時間スパンの間に(測定によって)生成されるMRI画像に基づいて、1つのMRI画像が予測されるか、または第2の時間スパンの間の検査領域を示す複数のMRI画像が予測される。
いくつかの非限定的な実施形態または態様において、第2の時間スパンは、第1の時間スパンの後に続く。
第2の時間スパンは、好ましくは、肝胆道相内の時間スパン、好ましくは造影剤の投与後少なくとも10分、好ましくは造影剤の投与後少なくとも20分に開始する時間スパンである。
第1の時間スパンの間の検査領域を示す複数の測定MRI画像は、予測モデルに供給される。予測モデルは、第1の時間スパンの間の検査領域を示す複数のMRI画像に基づいて、第2の時間スパンを示す検査領域を示す1つまたは複数のMRI画像を予測するように構成されるモデルである。
これに関連して、用語“予測”は、第2の時間スパンの間の検査領域を示すMRI画像が、第1の時間スパンの間の検査領域を示すMRI画像を使用して計算されることを意味する。
予測モデルは、好ましくは、教師あり機械学習プロセスにおいて自己学習アルゴリズムの助けを借りて作成されたものである。学習は、第1および第2の時間スパンの多様なMRI画像を含む訓練データを使用することによって達成される。
自己学習アルゴリズムは、機械学習の間、訓練データに基づく統計的モデルを生成する。これは、例が単に暗記されるのではなく、アルゴリズムが訓練データ内のパターンおよび規則性を“認識する”ことを意味する。故に、予測モデルは、未知のデータを査定することができる。検証データは、未知のデータの査定の質を試験するために使用され得る。
予測モデルは、教師あり学習によって訓練され、すなわち、第1の時間スパンからの複数のMRI画像が、アルゴリズムに連続して提示され、第2の時間スパン内のMRI画像のうちのどれがこれらの複数のMRI画像と関連付けられるかが通知される。次いで、アルゴリズムは、第1の時間スパンの未知の複数のMRI画像ついて第2の時間スパン内の1つまたは複数のMRI画像を予測するために、第1の時間スパンの複数のMRI画像と第2の時間スパンのMRI画像との関係性を学習する。
教師あり学習によって訓練される自己学習システムは、例えば、C. Perez: Machine Learning Techniques: Supervised Learning and Classification, Amazon Digital Services LLC - Kdp Print Us,2019,ISBN 1096996545, 9781096996545に広く説明される。
好ましくは、予測モデルは、人工ニューラルネットワークである。
そのような人工ニューラルネットワークは、処理要素の少なくとも3つの層:入力ニューロン(ノード)を有する第1の層、少なくとも1つの出力ニューロン(ノード)を有する第Nの層、および第N-2の内層を備え、Nは、自然数であり、2よりも大きい。
入力ニューロンは、入力値としてデジタルMRI画像を受信する役割を果たす。通常、デジタルMRI画像のピクセルまたはボクセルごとに1つの入力ニューロンが存在する。追加の入力値(例えば、検査領域に関する情報、検査対象に関する情報、および/またはMRI画像を生成するときに優先される健康状態に関する情報)のための追加の入力ニューロンが存在し得る。
そのようなネットワークにおいて、出力ニューロンは、第1の時間スパンの複数のMRI画像について第2の時間スパンの1つまたは複数のMRI画像を予測する役割を果たす。
入力ニューロンと出力ニューロンとの間の層の処理要素は、既定の結合重みを伴う既定のパターンで互いに結合される。
好ましくは、人工ニューラルネットワークは、いわゆる畳み込みニューラルネットワーク(略してCNN)である。
畳み込みニューラルネットワークは、行列の形態にある入力データを処理することができる。これは、行列(例えば、幅×高さ×色チャネル)として描写されるデジタルMRI画像を入力データとして使用することを可能にする。対照的に、例えば多層パーセプトロン(MLP)の形態にある、通常のニューラルネットワークは、入力としてベクトルを必要とし、すなわち、MRI画像を入力として使用するため、MRI画像のピクセルまたはボクセルは、長鎖で連続的にロールアウトされる必要がある。その結果、通常のニューラルネットワークは、例えば、MRI画像内の物体の位置から独立してMRI画像内の物体を認識することはできない。MRI画像内の異なる位置における同じ物体は、全く異なる入力ベクトルを有する。
CNNは、交互に繰り返されるフィルタ(畳み込み層)および集約層(プーリング層)、ならびに最後に、“通常の”完全に結合したニューロン(高密度層/全結合層)の1つの層または複数の層から本質的になる。
シーケンス(MRI画像のシーケンス)を分析するとき、空間および時間は、同等の次元として取り扱われ、例えば、3D折り畳みにより処理され得る。これは、Baccoucheら(Sequential Deep Learning for Human Action Recognition; International Workshop on Human Behavior Understanding, Springer 2011,pages 29-39)、およびJiら(3D Convolutional Neural Networks for Human Action Recognition, IEEE Transactions on Pattern Analysis and Machine Intelligence,35(1),221-231)による論文に示されている。さらには、Karpathyら(Large-scale Video Classification with Convolutional Neural Networks; Proceedings of the IEEE conference on Computer Vision and Pattern Recognition,2014,pages 1725-1732)、およびSimonyan&Zisserman(Two-stream Convolutional Networks for Action Recognition in Videos; Advances in Neural Information Processing Systems,2014,pages 568-576)による刊行物において説明されるように、時間および空間に関与している異なるネットワークを訓練し、最後に特徴を併合することが可能である。
リカレントニューラルネットワーク(RNN)は、層同士のフィードバック結合を含むいわゆるフィードフォワードニューラルネットワークのファミリーである。RNNは、ニューラルネットワークの異なる部分を介したパラメータデータの共通利用によって、順次データのモデリングを可能にする。RNNのためのアーキテクチャは、サイクルを含む。サイクルは、RNNからの出力データの少なくとも一部分が、シーケンス内の後続入力を処理するためのフィードバックとして使用されるため、変数の現在の値の、将来の時間点における独自の値に対する影響を表す。
詳細は、例えば、S.Khanら:A Guide to Convolutional Neural Networks for Computer Vision, Morgan & Claypool Publishers 2018,ISBN 1681730227, 9781681730226から収集され得る。
ニューラルネットワークの訓練は、例えば、バックプロパゲーション法によって実行され得る。これに関連して、ネットワークに対して求められるのは、できる限り信頼性の高い、所与の出力ベクトルへの所与の入力ベクトルのマッピングである。マッピング品質は、誤差関数により説明される。目標は、誤差関数を最小限にすることである。バックプロパゲーション法の場合、人工ニューラルネットワークは、結合重みを変更することによって教示される。
訓練された状態において、処理要素間の結合重みは、第1の時間スパンの間の検査領域を示す新規の複数のMRI画像について第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を予測するために使用され得る、第1の時間スパンの複数のMRI画像と第2の時間スパンのMRI画像との関係性に関する情報を含む。
交差検証法は、データを訓練データセットおよび検証データセットに分割するために使用され得る。訓練データセットは、ネットワーク重みのバックプロパゲーション訓練に使用される。検証データセットは、訓練されたネットワークが未知の複数のMRI画像に適用され得る予測の正確性をチェックするために使用される。
既に示されるように、検査対象、検査領域、および/または検査条件に関する更なる情報はまた、訓練、検証、および予測のために使用され得る。
検査対象に関する情報の例は、性別、年齢、体重、身長、既往歴、既に摂取した薬剤の働きおよび持続期間および量、血圧、中心静脈圧、呼吸数、血清アルブミン、総ビリルビン量、血糖、鉄含量、肺活量等である。これらは、例えば、データベースまたは電子患者ファイルから収集され得る。
検査領域に関する情報の例は、以前から存在する健康状態、手術、肺部分切除術、肝移植、肝臓の鉄、脂肪肝等である。
第1の時間スパンの間の検査領域を示す複数のMRI画像は、それらが予測モデルに供給される前に、運動補正の対象となることが想定される。そのような運動補正は、第1のMRI画像のピクセルまたはボクセルが、経時的に下流にある第2のMRI画像の対応するピクセルまたはボクセルと同じ検査領域を示すことを確実にする。運動補正法は、欧州特許第3118644号、欧州特許第3322997号、米国特許第20080317315号、米国特許第20170269182号、米国特許第20140062481号、欧州特許第2626718号に説明される。
非限定的な実施形態または態様に従うシステムは、非限定的な実施形態または態様に従う方法を実行し得る。
本システムは、受信装置、制御および計算装置、ならびに出力装置を備える。
記載された装置は、単一のコンピュータシステムの構成要素であることが想定されるが、記載された装置は、1つの装置から別の装置へデータおよび/または制御信号を送信するため、ネットワークを介して互いに結合される、複数の別個のコンピュータシステムの構成要素であることも想定される。
“コンピュータシステム”は、プログラム可能な計算規則によってデータを処理する電子データ処理のためのシステムである。そのようなシステムは、通常、“コンピュータ”を備え、そのような装置は、論理演算を実施するためのプロセッサ、およびさらに周辺装置を備える。
コンピュータ技術において、“周辺装置”とは、コンピュータに結合され、コンピュータの制御に役立つ、ならびに/または入力および出力デバイスとしての役割を果たす、すべてのデバイスを指す。それらの例は、モニタ(スクリーン)、プリンタ、スキャナ、マウス、キーボード、ドライブ、カメラ、マイク、ラウドスピーカなどである。内部ポートおよび拡張カードもまた、コンピュータ技術においては周辺装置であると見なされる。
今日のコンピュータシステムは、しばしば、デスクトップPC、ポータブルPC、ラップトップ、ノートブック、ネットブック、およびタブレットPC、ならびにいわゆるハンドヘルド(例えば、スマートフォン)に分割され、すべてのこれらのシステムは、本開示の実施形態または態様を実行するために利用され得る。
コンピュータシステムへの入力は、例えば、キーボード、マウス、マイク、および/または同様のものなど、入力手段により達成される。
本システムは、第1の時間スパンの間の検査領域を示す複数のMRI画像を受信し、これらのデータおよび任意選択的に更なるデータに基づいて、第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を生成する(予測する、計算する)ように構成され得る。
制御および計算装置は、受信装置の制御、様々な装置間のデータおよび信号フローの整合、ならびにMRI画像の計算に役立つ。複数の制御および計算装置が存在することが想定される。
受信装置は、複数のMRI画像の受信に役立つ。複数のMRI画像は、例えば、磁気共鳴システムから送信され得るか、またはデータ記憶媒体から読み出され得る。磁気共鳴システムは、システムの構成要素であり得る。しかしながら、システムは、磁気共鳴システムの構成要素であることも想定される。
MRI画像のシーケンスおよび任意選択的に更なるデータは、受信装置から制御および計算装置へ送信される。
制御および計算装置は、第1の時間スパンの間の検査領域を示す複数のMRI画像に基づいて、1つまたは複数のMRI画像を予測するように構成され、予測MRI画像は、第2の時間スパンの間の検査領域を示す。好ましくは、制御および計算装置のメモリ内に読み込まれ得るのは、第2の時間スパンのMRI画像を計算するために使用される予測モデルである。予測モデルは、好ましくは、教師あり学習によって、自己学習アルゴリズムの助けを借りて生成(訓練)されたものである。
出力装置を介して、予測MRI画像は、(例えば、スクリーン上に)表示されるか、(例えば、プリンタを介して)出力されるか、またはデータ記憶媒体に記憶され得る。
更なる非限定的な実施形態または態様は、以下の番号付きの条項に明記される。
条項1.2つ以上の時間点において測定される患者の組織のボクセルのパラメータと関連付けられた測定情報を獲得するステップであって、2つ以上の時間点は、組織のボクセルの1つまたは複数の特性が、2つ以上の時間点のうちの単一の時間点において測定されるボクセルのパラメータに基づいて生成される画像において分離可能である前に発生する、ステップと、2つ以上の時間点におけるボクセルのパラメータに基づいて、組織のボクセルの1つまたは複数の特性を決定するステップとを含む、コンピュータ実施方法。
条項2.組織のボクセルの1つまたは複数の特性は、患者および患者の健康状態のうちの少なくとも1つと関連付けられた情報に基づいてさらに決定される、条項1に記載のコンピュータ実施方法。
条項3.組織のボクセルの1つまたは複数の特性は、少なくとも1つの2つ以上の時間点に対応する時間点について決定される、条項1および2に記載のコンピュータ実施方法。
条項4.組織のボクセルの1つまたは複数の特性は、2つ以上の時間点の後の時間点について決定される、条項1~3のいずれか一項に記載のコンピュータ実施方法。
条項5.1つまたは複数の特性に基づいて、2つ以上の時間点の後の時間点における組織のボクセルの1つまたは複数の特性を含む1つまたは複数の画像を生成するステップをさらに含む、条項1~4のいずれか一項に記載のコンピュータ実施方法。
条項6.患者の組織のボクセルのパラメータと関連付けられた測定情報が、組織のボクセルの1つまたは複数の特性を決定することと関連付けられた測定情報のしきい値量を含むことを決定するステップと、測定情報が測定情報のしきい値量を含むことを決定することに応答して、測定情報の取得を自動的に停止するように撮像システムを制御するステップとをさらに含む、条項1~5のいずれか一項に記載のコンピュータ実施方法。
条項7.1つまたは複数の特性を決定するステップは、患者の組織のボクセルのパラメータと関連付けられた測定情報を予測モデルに供給することを含み、予測モデルは、2つ以上の時間点におけるパラメータと関連付けられた測定情報に基づいて、組織のボクセルの1つまたは複数の特性を予測するように、教師あり学習によって訓練されている、条項1~6のいずれか一項に記載のコンピュータ実施方法。
条項8.1つまたは複数の特性を決定するステップは、組織のボクセルの薬物動力学的/薬力学的(PK/PD)モデルを、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされたPK/PDモデルに基づいて、組織のボクセルの1つまたは複数の特性を決定することを含む、条項1~7のいずれか一項に記載のコンピュータ実施方法。
条項9.1つまたは複数の特性を決定するステップは、パラメータについて事前計算される複数の複数のPK/PD曲線のうちの1つのPK/PD曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点におけるパラメータにフィッティングされたPK/PD曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することを含む、条項1~8のいずれか一項に記載のコンピュータ実施方法。
条項10.1つまたは複数の特性を決定するステップは、基底関数のセットを用いて、組織のボクセルの1つまたは複数の特性を表す曲線を近似すること、近似された曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた近似された曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することを含む、条項1~9のいずれか一項に記載のコンピュータ実施方法。
条項11.1つまたは複数の特性を決定するステップは、基底関数のセットを用いてパラメータについて事前計算される複数の曲線のうちの1つの曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することを含む、条項1~10のいずれか一項に記載のコンピュータ実施方法。
条項12.2つ以上の時間点において測定される患者の組織のボクセルのパラメータと関連付けられた測定情報を獲得することであって、2つ以上の時間点は、組織のボクセルの1つまたは複数の特性が、2つ以上の時間点のうちの単一の時間点において測定されるボクセルのパラメータに基づいて生成される画像において分離可能である前に発生する、獲得すること、および2つ以上の時間点におけるボクセルのパラメータに基づいて、組織のボクセルの1つまたは複数の特性を決定すること、を行うようにプログラムおよび/または構成される1つまたは複数のプロセッサを備えるシステム。
条項13.組織のボクセルの1つまたは複数の特性は、患者および患者の健康状態のうちの少なくとも1つと関連付けられた情報に基づいてさらに決定される、条項12に記載のシステム。
条項14.組織のボクセルの1つまたは複数の特性は、少なくとも1つの2つ以上の時間点に対応する時間点について決定される、条項12および13のいずれか一項に記載のシステム。
条項15.組織のボクセルの1つまたは複数の特性は、2つ以上の時間点の後の時間点について決定される、条項12~14のいずれか一項に記載のシステム。
条項16.1つまたは複数のプロセッサは、1つまたは複数の特性に基づいて、2つ以上の時間点の後の時間点における組織のボクセルの1つまたは複数の特性を含む1つまたは複数の画像を生成するようにさらにプログラムおよび/または構成される、条項12~15のいずれか一項に記載のシステム。
条項17.1つまたは複数のプロセッサは、患者の組織のボクセルのパラメータと関連付けられた測定情報が、組織のボクセルの1つまたは複数の特性を決定することと関連付けられた測定情報のしきい値量を含むことを決定し、測定情報が測定情報のしきい値量を含むことを決定することに応答して、測定情報の取得を自動的に停止するように撮像システムを制御するようにさらにプログラムおよび/または構成される、条項12~16のいずれか一項に記載のシステム。
条項18.1つまたは複数のプロセッサは、患者の組織のボクセルのパラメータと関連付けられた測定情報を予測モデルに供給することによって1つまたは複数の特性を決定するようにさらにプログラムおよび/または構成され、予測モデルは、2つ以上の時間点におけるパラメータと関連付けられた測定情報に基づいて、組織のボクセルの1つまたは複数の特性を予測するように、教師あり学習によって訓練されている、条項12~17のいずれか一項に記載のシステム。
条項19.1つまたは複数のプロセッサは、組織のボクセルの薬物動力学的/薬力学的(PK/PD)モデルを、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされたPK/PDモデルに基づいて、組織のボクセルの1つまたは複数の特性を決定することによって1つまたは複数の特性を決定するようにさらにプログラムおよび/または構成される、条項12~18のいずれか一項に記載のシステム。
条項20.1つまたは複数のプロセッサは、パラメータについて事前計算される複数の複数のPK/PD曲線のうちの1つのPK/PD曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点におけるパラメータにフィッティングされたPK/PD曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することによって1つまたは複数の特性を決定するようにさらにプログラムおよび/または構成される、条項12~19のいずれか一項に記載のシステム。
条項21.1つまたは複数のプロセッサは、基底関数のセットを用いて、組織のボクセルの1つまたは複数の特性を表す曲線を近似すること、近似された曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた近似された曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することによって1つまたは複数の特性を決定するようにさらにプログラムおよび/または構成される、条項12~20のいずれか一項に記載のシステム。
条項22.1つまたは複数のプロセッサは、基底関数のセットを用いてパラメータについて事前計算される複数の曲線のうちの1つの曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することによって1つまたは複数の特性を決定するようにさらにプログラムおよび/または構成される、条項12~21のいずれか一項に記載のシステム。
条項23.プログラム命令を含む少なくとも1つの非一時的なコンピュータ可読媒体を備えるコンピュータプログラム製品であって、プログラム命令は、少なくとも1つのプロセッサによって実行されるとき、少なくとも1つのプロセッサに、2つ以上の時間点において測定される患者の組織のボクセルのパラメータと関連付けられた測定情報を獲得することであって、2つ以上の時間点は、組織のボクセルの1つまたは複数の特性が、2つ以上の時間点のうちの単一の時間点において測定されるボクセルのパラメータに基づいて生成される画像において分離可能である前に発生する、獲得すること、および2つ以上の時間点におけるボクセルのパラメータに基づいて、組織のボクセルの1つまたは複数の特性を決定すること、を行わせる、コンピュータプログラム製品。
条項24.組織のボクセルの1つまたは複数の特性は、患者および患者の健康状態のうちの少なくとも1つと関連付けられた情報に基づいてさらに決定される、条項23に記載のコンピュータプログラム製品。
条項25.組織のボクセルの1つまたは複数の特性は、少なくとも1つの2つ以上の時間点に対応する時間点について決定される、条項23および24のいずれか一項に記載のコンピュータプログラム製品。
条項26.組織のボクセルの1つまたは複数の特性は、2つ以上の時間点の後の時間点について決定される、条項23~25のいずれか一項に記載のコンピュータプログラム製品。
条項27.命令は、少なくとも1つのプロセッサに、1つまたは複数の特性に基づいて、2つ以上の時間点の後の時間点における組織のボクセルの1つまたは複数の特性を含む1つまたは複数の画像を生成することをさらに行わせる、条項23~26のいずれか一項に記載のコンピュータプログラム製品。
条項28.命令は、少なくとも1つのプロセッサに、患者の組織のボクセルのパラメータと関連付けられた測定情報が、組織のボクセルの1つまたは複数の特性を決定することと関連付けられた測定情報のしきい値量を含むことを決定し、測定情報が測定情報のしきい値量を含むことを決定することに応答して、測定情報の取得を自動的に停止するように撮像システムを制御することをさらに行わせる、条項23~27のいずれか一項に記載のコンピュータプログラム製品。
条項29.命令は、少なくとも1つのプロセッサに、患者の組織のボクセルのパラメータと関連付けられた測定情報を予測モデルに供給することによって1つまたは複数の特性を決定することを行わせ、予測モデルは、2つ以上の時間点におけるパラメータと関連付けられた測定情報に基づいて、組織のボクセルの1つまたは複数の特性を予測するように、教師あり学習によって訓練されている、条項23~28のいずれか一項に記載のコンピュータプログラム製品。
条項30.命令は、少なくとも1つのプロセッサに、組織のボクセルの薬物動力学的/薬力学的(PK/PD)モデルを、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされたPK/PDモデルに基づいて、組織のボクセルの1つまたは複数の特性を決定することによって、1つまたは複数の特性を決定することを行わせる、条項23~29のいずれか一項に記載のコンピュータプログラム製品。
条項31.命令は、少なくとも1つのプロセッサに、パラメータについて事前計算される複数の複数のPK/PD曲線のうちの1つのPK/PD曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点におけるパラメータにフィッティングされたPK/PD曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することによって、1つまたは複数の特性を決定することを行わせる、条項23~30のいずれか一項に記載のコンピュータプログラム製品。
条項32.命令は、少なくとも1つのプロセッサに、基底関数のセットを用いて、組織のボクセルの1つまたは複数の特性を表す曲線を近似すること、近似された曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた近似された曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することによって、1つまたは複数の特性を決定することを行わせる、条項23~31のいずれか一項に記載のコンピュータプログラム製品。
条項33.命令は、少なくとも1つのプロセッサに、基底関数のセットを用いてパラメータについて事前計算される複数の曲線のうちの1つの曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングすること、および2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定することによって、1つまたは複数の特性を決定することを行わせる、条項23~32のいずれか一項に記載のコンピュータプログラム製品。
本開示の実施形態または態様の追加の利点および詳細は、添付の概略図に例証される例示的な実施形態を参照して、以下により詳細に説明される。
肝臓動脈(A)、肝臓静脈(V)、および肝細胞(P)内の造影剤の濃度の経時的プロファイル、ならびに肝臓動脈組織、肝臓静脈組織、および肝細胞組織の各々を含むボクセルについての総計(summed enhancement)(S)のグラフである。
非限定的な実施形態または態様に従うシステムを示す図である。
組織特性を予測、予想、および/または査定するためのプロセスの非限定的な実施形態または態様のフローチャートである。
動的相および肝胆道相の間の肝臓のMRI画像を示す図である。
第1の時間スパン内の肝臓を示す3つのMRI画像(1)、(2)、および(3)、ならびに第2の時間スパン内の肝臓を示すMRI画像(4)を示す図である。
本明細書に説明されるシステム、デバイス、製品、装置、および/または方法が実装され得る環境の非限定的な実施形態または態様を示す図である。
組織特性を予測、予想、および/または査定するためのプロセスの非限定的な実施形態または態様のフローチャートである。
肝臓組織の薬物動力学的モデルの非限定的な例を示す図である。
“無地の白い肝臓”画像を生成するために、肝特異性造影剤についてHBPから知られる造影情報を肝臓の初期相のMRI画像(1)、(2)と組み合わることによって合成されるMRI画像(3)を示す図である。
測定された造影強化画像から合成される、非強化T1w画像(1)、(2)、(3)を示す図である。
1つまたは複数の第2の注入後の画像から1つまたは複数の平衡後の更新画像を減ずることによって導出される平衡取り込みとは別の血管強化を示す画像を示す図である。
本開示は、様々な代替的な変異形およびステップシーケンスを仮定し得るが、逆のことが明示的に記される場合はこの限りではない、ということを理解されたい。添付の図面に例証される、および以下の説明に記載される特定のデバイスおよびプロセスは、単に例示的および非限定的な実施形態または態様すぎないということも理解されたい。故に、本明細書に開示される実施形態または態様に関連した特定の寸法および他の物理的特性は、限定と見なされるべきではない。
以後説明の目的のため、用語“端部”、“上方”、“下方”、“右”、“左”、“垂直”、“水平”、“上部”、“下部”、“横”、および“縦”、ならびにそれらの派生語は、それらが描写図において配向される際の実施形態または態様に関連するものとする。しかしながら、実施形態または態様は、様々な代替的な変異形およびステップシーケンスを仮定し得るが、逆のことが明示的に記される場合はこの限りではない、ということを理解されたい。添付の図面に例証される、および以下の説明に記載される特定のデバイスおよびプロセスは、単に非限定的な例示的な実施形態または態様にすぎないということも理解されたい。故に、本明細書に開示される実施形態または態様の実施形態または態様に関連した特定の寸法および他の物理的特性は、別途示されない限り、限定と見なされるべきではない。
本明細書内で使用される態様、構成要素、要素、構造、作用、ステップ、機能、命令、および/または同様のものは、重要または必須と解釈されるべきではないが、そのようなものとして明示的に説明される場合はその限りではない。また、本明細書で使用される場合、冠詞“a(1つの)”および“an(1つの)”は、1つまたは複数の事項を含むことが意図され、“1つまたは複数”および“少なくとも1つ”と同義で使用され得る。さらには、本明細書で使用される場合、用語“セット”は、1つまたは複数の事項(例えば、関連事項、非関連事項、関連事項および非関連事項の組み合わせなど)を含むことが意図され、“1つまたは複数”および“少なくとも1つ”と同義で使用され得る。1つの事項のみが意図される場合、用語“1つ”または同様の言葉が使用される。また、本明細書で使用される場合、用語“has(有する)”、“have(有する)”、“having(有すること)”、または同様のものは、オープンエンドの用語であることが意図される。さらに、語句“基づく”は、別途明示的に記されない限り、“少なくとも部分的に基づく”を意味することが意図される。
本明細書で使用される場合、用語“通信”および“通信する”は、情報(例えば、データ、信号、メッセージ、命令、コマンド、および/または同様のもの)の受信、受領、送信、転送、提供、および/または同様のものを指し得る。1つの装置(例えば、デバイス、システム、デバイスもしくはシステムの構成要素、それらの組み合わせ、および/または同様のもの)が別の装置と通信状態にあるとは、1つの装置が、直接的または間接的に、別の装置から情報を受信すること、および/または別の装置へ情報を送信することができることを意味する。これは、本質的に有線および/またはワイヤレスである、直接または間接接続を指し得る。追加的に、2つの装置は、送信される情報が、第1の装置と第2の装置との間で改変、処理、中継、および/またはルーティングされ得るとしても、互いと通信状態にあり得る。例えば、第1の装置が情報を受動的に受信し、その情報を第2の装置に能動的に送信しないとしても、第1の装置は、第2の装置と通信状態にあり得る。別の例として、少なくとも1つの中継装置(例えば、第1の装置と第2の装置との間に位置する第3の装置)が、第1の装置から受信される情報を処理し、処理した情報を第2の装置に通信するとしても、第1の装置は、第2の装置と通信状態にあり得る。いくつかの非限定的な実施形態または態様において、メッセージは、データを含むネットワークパケット(例えば、データパケットおよび/または同様のもの)を指し得る。多数の他の配置が可能であるということを理解されたい。
本明細書で使用される場合、用語“コンピューティングデバイス”は、1つもしくは複数のネットワークと、またはこれを介して、直接的または間接的に通信するように構成される1つまたは複数の電子デバイスを指し得る。コンピューティングデバイスは、モバイルもしくはポータブルコンピューティングデバイス、デスクトップコンピュータ、サーバ、および/または同様のものであり得る。さらには、用語“コンピュータ”は、データを受信、処理、および出力するために必要な構成要素を含み、通常、ディスプレイ、プロセッサ、メモリ、入力デバイス、およびネットワークインターフェースを含む、任意のコンピューティングデバイスを指し得る。“コンピューティングシステム”は、1つまたは複数のコンピューティングデバイスまたはコンピュータを含み得る。“アプリケーション”または“アプリケーションプログラムインターフェース”(API)は、クライアント側フロントエンドおよび/またはクライアントからデータを受信するためのサーバ側バックエンドなど、ソフトウェア構成要素間の対話を促進するためにプロセッサによって実行され得る、コンピュータ可読媒体に記憶されたコンピュータコードまたは他のデータを指す。“インターフェース”は、ユーザが、直接的または間接的に(例えば、キーボード、マウス、タッチスクリーンなどを通じて)対話し得る1つまたは複数のグラフィカルユーザインターフェース(GUI)などの生成ディスプレイを指す。さらに、ネットワーク環境内で直接的または間接的に通信する、複数のコンピュータ、例えば、サーバ、または、車両コンピューティングシステムを含む自律車両などの他のコンピュータ化デバイスが、“システム”または“コンピューティングシステム”を構成し得る。
本明細書に説明されるシステムおよび/または方法は、ハードウェア、ソフトウェア、またはハードウェアおよびソフトウェアの組み合わせの形態で実装され得ることは明白であるものとする。これらのシステムおよび/または方法を実装するために使用される実際の専用制御ハードウェアまたはソフトウェアコードは、実装形態の限定ではない。故に、本システムおよび/または方法の動作および挙動は、特定のソフトウェアコードを参照せずに本明細書に説明され、ソフトウェアおよびハードウェアは、本明細書内の説明に基づいてシステムおよび/または方法を実装するために設計され得ることを理解されたい。
いくつかの非限定的な実施形態または態様は、しきい値と関連して本明細書に説明される。本明細書で使用される場合、しきい値を満足することは、値が、しきい値よりも大きいこと、しきい値超であること、しきい値よりも高いこと、しきい値以上であること、しきい値未満であること、しきい値よりも小さいこと、しきい値よりも低いこと、しきい値以下であること、しきい値に等しいことなどを指し得る。
これより図1を参照すると、図1は、肝臓動脈(A)、肝臓静脈(V)、肝細胞(P)内の造影剤の濃度の経時的プロファイル、および3つのタイプの組織すべて(例えば、肝臓動脈組織、肝臓静脈組織、および肝細胞組織の各々)を含むボクセルについての総計(S)を概略的に示し、上で本明細書に既に詳細に説明されている。
図2も参照すると、図2は、非限定的な実施形態または態様に従うシステムを概略的に示す。システム(10)は、受信装置(11)、制御および計算装置(12)、ならびに出力装置(13)を備える。
図3も参照すると、図3は、非限定的な実施形態または態様に従う方法のフローチャートである。方法(100)は、
複数のMRI画像を受信するステップ(110)であって、MRI画像が第1の時間スパンの間の検査領域を示す、ステップ(110)と、
複数のMRI画像を予測モデルに供給するステップ(120)であって、予測モデルが、第1の時間スパンの間の検査領域を示すMRI画像に基づいて、第2の時間スパンの間の検査領域を示す1つまたは複数のMRI画像を予測するように、教師あり学習によって訓練されている、ステップ(120)と、
第2の時間スパンの間の検査領域を示す1つまたは複数の予測MRI画像を、予測モデルによって生成するステップ(130)と、
1つもしくは複数の予測MRI画像を表示および/もしくは出力する、ならびに/または1つもしくは複数の予測MRI画像をデータ記憶媒体に記憶するステップ(140)と、を含む。
図4も参照すると、図4は、動的相および肝胆道相の間の肝臓のMRI画像を示す。図4(a)、図4(b)、図4(c)、図4(d)、図4(e)、および図4(f)では、異なる時間点における肝臓の同じ横断面が常に描写される。図4(a)、図4(b)、図4(d)、および図4(f)に入力された参照記号は、図4(a)、図4(b)、図4(c)、図4(d)、図4(e)、および図4(f)のすべてに当てはまり、それらは各々、単に明白性の目的のために一度のみ入力される。
図4(a)は、肝胆道造影剤の静脈内投与前の肝臓(L)の横断面を示す。図4(a)および図4(b)によって描写される時間点の間の時間点において、肝胆道造影剤は、ボーラスとして静脈内に投与された。これは、図4(b)において、肝臓動脈(A)を介して肝臓に到達する。したがって、肝臓動脈は、信号強化を伴って描写される(動脈相)。主に動脈を介して血液が供給される腫瘍(T)は、同様に、より明るい(信号強化された)領域として肝細胞組織よりも目立つ。図4(c)に描写される時間点において、造影剤は、静脈を介して肝臓に到達する。図4(d)では、静脈血管(V)は、明るい(信号強化された)領域(静脈相)として肝臓組織よりも目立つ。同時に、主に静脈を介して造影剤が供給される健康な肝細胞内の信号強度は、継続的に上昇する(図4(c)→4(d)→4(e)→4(f))。図4(f)に描写される肝胆道相では、肝細胞(P)は、信号強化を伴って描写され、血管および腫瘍は、もはや造影剤を有さず、したがって暗く描写される。
図5も参照すると、図5は、第1の時間スパン内の肝臓を示す3つのMRI画像(1)、(2)、および(3)がどのように予測モデル(PM)に供給されるかを例示的および概略的に示す。予測モデルは、3つのMRI画像(1)、(2)、および(3)から、第2の時間スパン内の肝臓を示すMRI画像(4)を計算する。MRI画像(1)、(2)、および(3)は、例えば、図4(b)、図4(c)、および図4(d)に示されるMRI画像を示すことができ、MRI画像(4)は、例えば、図4(f)に示されるMRI画像であり得る。
既存の撮像システムおよび/または画像分析システムによって満たされていない医学的必要性としては、少なくとも以下のことが挙げられる:健康な組織および/または異なるグレードの病理組織の識別および査定、異なる疾病および/または病理ステージの区別、改善された撮像/走査ワークフロー、細胞特異的および/または細胞機能的な情報の生成、ならびに中空または管状の器官または器官系に関する機能情報。
本開示の非限定的な実施形態または態様は、少なくとも以下のことを提供および/または改善する:健康な組織および/もしくは異なるグレードの病理組織(例えば、組織線維化/肝硬変、炎症、脂肪浸潤、機能障害/死亡、および/または同様のものなどの広範性病理、良性もしくは悪性腫瘍、および/または同様のものなどの限局性病理など)の識別および査定、異なる疾病および/もしくは病理ステージ(例えば、広範性病理、限局性病理など)の区別、撮像/走査ワークフロー(例えば、より高速の画像生成および/または取得など)、細胞特異的および/もしくは細胞機能的な情報の生成(例えば、特定の薬物もしくは代謝物を取り込む、および分泌する細胞能力および速度、ならびに/または同様のものなどの細胞機能、細胞酸化の量、特定の細胞抗原の発現、チャネル/膜タンパク質、ならびに/または同様のものなどの分子情報など)、ならびに/または中空もしくは管状の器官もしくは器官系に関する機能情報(例えば、門脈系などの、血管系、胆管系、血管圧など)。本開示の非限定的な実施形態または態様は、動脈循環体積、門脈体積、静脈体積、細胞外体積、胆管体積、正常肝細胞体積、脂肪細胞体積、線維化体積、クッパー細胞体積、幹細胞体積、他の肝細胞体積、および/または転移性もしくは他の病変もしくは非肝細胞体積のうちの1つまたは複数のボクセル固有の推測を提供する。
これより図6を参照すると、図6は、本明細書に説明されるデバイス、システム、方法、および/または製品が実装され得る、例となる環境600の図である。図6に示されるように、環境600は、造影剤注入器602、注入器制御および演算システム604、注入器ユーザインターフェース606、撮像装置608、撮像装置制御および演算システム610、撮像装置ユーザインターフェース612、画像分析および演算システム614、画像分析インターフェース616、画像およびデータレビューインターフェース618、病院情報システム620、ならびに/またはクラウドコンピューティングおよびオフサイトリソース622を含む。環境600のシステムおよび/またはデバイスは、有線接続、ワイヤレス接続、または有線およびワイヤレス接続の組み合わせにより相互接続し得る。例えば、環境600のシステムおよび/またはデバイスは、セルラネットワーク(例えば、ロングタームエボリューション(LTE)ネットワーク、第三世代(3G)ネットワーク、第四世代(4G)ネットワーク、第五世代ネットワーク(5G)、符号分割多重アクセス(CDMA)ネットワークなど)、公衆陸上移動ネットワーク(PLMN)、ローカルエリアネットワーク(LAN)、広域ネットワーク(WAN)、メトロポリタンエリアネットワーク(MAN)、電話ネットワーク(例えば、公衆交換電話ネットワーク(PSTN))、プライベートネットワーク、アドホックネットワーク、イントラネット、インターネット、光ファイバベースのネットワーク、クラウドコンピューティングネットワーク、近距離無線通信ネットワーク、および/もしくは同様のもの、ならびに/またはこれらもしくは他のタイプのネットワークの組み合わせを介して、相互接続(例えば、情報および/またはデータなどを通信)し得る。造影剤注入器602、注入器制御および演算システム604、および/または注入器ユーザインターフェース606は、1つまたは複数の注入プロトコルを設定し、その1つまたは複数の注入プロトコルに従って1つまたは複数の造影剤を患者に送達するために、それぞれのソフトウェアおよび/またはハードウェアを有する造影剤注入システムを備え得る。いくつかの非限定的な実施形態または態様において、造影剤注入器602、注入器制御および演算システム604、ならびに/または注入器ユーザインターフェース606は、米国特許第6,643,537号および同第7,937,134号に説明され、国際出願第2019046299A1号に公開されるような造影剤注入システムを含み得、これらの文献の各々の内容全体が参照により本明細書に組み込まれる。いくつかの非限定的な実施形態または態様において、造影剤注入器602、注入器制御および演算システム604、ならびに/または注入器ユーザインターフェース606は、MEDRAD(登録商標)Stellant FLEX CT注入システム、MEDRAD(登録商標)MRXperion MR注入システム、MEDRAD(登録商標)Mark 7 Arterion 注入システム、MEDRAD(登録商標)Intego PET注入システム、MEDRAD(登録商標)Spectris Solaris EP MR注入システム、Certegra(登録商標)Workstationを用いたMEDRAD(登録商標)Stellant CT注入システム、および/または同様のものを備え得る。
撮像装置608、撮像装置制御および演算システム610、ならびに/または撮像装置ユーザインターフェース612は、撮像プロトコルを設定し、患者の非造影および造影強化スキャンを取得するために、それぞれのソフトウェアおよび/またはハードウェアを有する撮像システムを備え得る。いくつかの非限定的な実施形態または態様において、撮像装置608、撮像装置制御および演算システム610、ならびに/または撮像装置ユーザインターフェース612は、MRIシステム(例えば、T1、T2、TWI、PD、マッピング(脂肪、鉄)、多パラメータ手法、過分極MRI、MRフィンガープリンティング、エラストグラフィなどに基づいたMRI)、コンピュータ断層撮影(CT)システム、超音波システム、単光子放射型コンピュータ断層撮影(SPECT)システム、ポジトロン放出断層撮影-磁気共鳴(PET/MRI)システム、ポジトロン放出断層撮影-コンピュータ断層撮影(PET/CT)システム、および/または他の画像診断システムを含み得る。いくつかの非限定的な実施形態または態様において、撮像装置608、撮像装置制御および演算システム610、ならびに/または撮像装置ユーザインターフェース612は、2019年12月11日に提出された米国特許出願第16/710,118号に説明されるような撮像システムを含み得、この文献の内容全体が参照により本明細書に組み込まれる。いくつかの非限定的な実施形態または態様において、撮像装置608、撮像装置制御および演算システム610、ならびに/または撮像装置ユーザインターフェース612は、シーメンスヘルシニアーズのSomatom Go CTシステム、ゼネラルエレクトリック社のSigna MRシステム、および/または同様のものを含み得る。
いくつかの非限定的な実施形態または態様において、造影剤注入器602、注入器制御および演算システム604、注入器ユーザインターフェース606、撮像装置608、撮像装置制御および演算システム610、撮像装置ユーザインターフェース612、画像分析および演算システム614、画像分析インターフェース616、画像およびデータレビューインターフェース618、病院情報システム620、ならびに/またはクラウドコンピューティングおよびオフサイトリソース622のうちの1つまたは複数は、本明細書に説明されるようなコンピュータシステム、ならびに/またはコンピュータシステムの1つもしくは複数の構成要素および/もしくは周辺機器を含み得る。
患者は、異なるタイプの細胞を有する複数の組織および/もしくは器官(例えば、肝実質細胞、クッパー細胞、免疫細胞、幹細胞などを含む肝臓)、求心性および/または遠心性供給/循環系(例えば、動脈、静脈、リンパ、または胆道系など)、ならびに/または異なる区画/空間(例えば、血管空間、間質空間、細胞内空間など)を含む生存生物(例えば、哺乳動物、ヒトなど)を含み得る。
造影剤注入システムによって患者に送達される造影剤は、患者を走査するために使用される撮像システムのタイプに従って選択または構成され得る。造影剤は、ガドリニウムベースの造影剤(GBCA)(例えば、MRIにおける使用のためなど)、ヨウ素ベースの造影剤(例えば、CTにおける使用のためなど)、超音波造影媒体(例えば、マイクロバブルなど)、および/または他のより非一般的な造影剤、例えば、鉄ベースの薬剤(例えば、小型もしくは超小型の超常磁性酸化鉄、マンガンベース、CO2、または他の薬剤)、血液プール剤(例えば、血管内の長い血液半減期を有する)、腎臓優性もしくは肝胆道優性排泄作用を有する薬剤、細胞内取り込み、および器官特異性もしくは細胞特異性取り込みを伴う薬剤、器官もしくは細胞特異性結合を伴う薬剤(例えば、細胞内取り込みなし)、長期“保留”を伴う薬剤(例えば、FDG)、ならびに/または同様のものなどを含み得る。造影剤は、放射性であってもよい。造影剤は、細胞マーカ特異性であってもよく、これは、造影剤が特定の細胞表面または細胞内部マーカと結合または相互作用することを意味する。関与する造影剤は、自然造影、例えば、ヘモグロビン内の酸素レベルであってもよい。造影剤は、例えば、ヘマトクリット値を減少させる陰性造影、および故に、血液をガスで置き換える、血液またはマイクロバブル内の赤血球信号であってもよい。通常のMRIでは、ガスは、信号を提供せず、CTでは、ガスは、増大した送信および故に減少したハウンスフィールド単位をもたらす。
造影剤注入システムは、単一の薬剤(例えば、それ自体によって送達される単一の造影剤など)、または同時間における組み合わせた複数の造影剤(例えば、複数の並行注入、一緒に混合された2つの流体の注入、または一方が他方の後に(例えば、複数の連続的注入))を送達し得る。撮像システムは、1つもしくは複数の時間点(相)における単一の画像取得または走査、1つもしくは複数の時間点における複数の取得(例えば、MRマッピング技術などを使用して)、および/または連続した撮像期間にわたる取得を実施することができる。
いくつかの非限定的な実施形態または態様において、造影剤注入は、静脈血管区画、動脈血管区画、リンパ管区画、および/または同様のものなどの1つまたは複数の異なる場所または区画に送達され得る。
これより図7を参照すると、図7は、組織特性を予測、予想、および/または査定するためのプロセス700の非限定的な実施形態または態様のフローチャートである。いくつかの非限定的な実施形態または態様において、プロセス700のステップのうちの1つまたは複数は、画像分析および演算システム614(例えば、画像分析および演算システム614の1つまたは複数のデバイス)によって、(例えば、完全に、部分的になど)実施され得る。いくつかの非限定的な実施形態または態様において、プロセス700のステップのうちの1つまたは複数は、撮像装置制御および演算システム610(例えば、撮像装置制御および演算システム610の1つまたは複数のデバイス)、病院情報システム620(例えば、病院情報システムの1つまたは複数のデバイス)、クラウドコンピューティングおよびオフサイトリソースシステム622(例えば、クラウドコンピューティングおよびオフサイトリソースシステム622の1つまたは複数のデバイス)、ならびに/または同様のものなど、画像分析および演算システム614とは別個の、またはこれを含む、別のデバイスまたはデバイスのグループによって、(例えば、完全に、部分的になど)実施され得る。
図7に示されるように、ステップ702において、プロセス700は、造影剤を患者に送達することを含む。例えば、造影剤注入システムは、造影剤を患者に送達(例えば、注入など)し得る。例として、造影剤注入器602、注入器制御および演算システム604、ならびに/または注入器ユーザインターフェース606は、1つまたは複数の注入プロトコルを設定し、その1つまたは複数の注入プロトコルに従って1つまたは複数の造影剤を患者に送達し得る。
いくつかの非限定的な実施形態または態様において、造影剤は、個々の相(例えば、動脈相、門脈相、静脈相など)が区別されることを可能にするために、小型で、短く、および比較的素早いボーラスで患者に送達され得る。例えば、造影剤濃度は、造影剤が身体内に分散する前の短期間の“安定状態”の間は比較的均質であり、短期の安定状態撮像を可能にし得る。いくつかの非限定的な実施形態または態様において、造影剤は、米国特許第9,436,989号に説明されるように、造影剤の注入時間を遅らせる、または引き延ばすように患者に送達され得、この文献の内容全体が参照により本明細書に組み込まれる。いくつかの非限定的な実施形態または態様において、造影剤は、血液および/または組織内で所望のMR造影剤濃度を達成するように患者に送達され得、その結果として、MRフィンガープリントは、2017年11月21日提出された米国特許出願第16/462,410号に説明されるような1つまたは複数の造影剤関連パラメータを含み得、この文献の内容全体が参照により本明細書に組み込まれる。いくつかの非限定的な実施形態または態様において、造影剤の注入速度/送達は、単一の高速ボーラスとして(例えば、相を区別するため)、時間Xから時間Yまでの実際の時間内の撮像を可能にするために単一のボーラスとして、2つの別々の画像期間(例えば、0.5~2分および次いで10または20分など)として、単一の緩徐なボーラス(例えば、視覚的に区別可能でない正常相を伴うなど)として、または二重ボーラス(例えば、第1のボーラスの送達、待機、撮像、第2のボーラスの送達、および撮像などを含むシーケンスで)として実施され得る。
図7に示されるように、ステップ704において、プロセス700は、患者の組織と関連付けられた測定情報を獲得することを含む。例えば、撮像システムは、1つまたは複数の撮像プロトコルを設定し、その1つまたは複数の撮像プロトコルに従って患者の組織と関連付けられた測定情報(例えば、患者の組織の1つまたは複数の画像またはスキャンなど)を獲得するために、患者の組織の非造影および/または造影強化スキャンを取得し得る。例として、撮像装置608、撮像装置制御および演算システム610、ならびに/または撮像装置ユーザインターフェース612(ならびに/または画像分析および演算システム614)は、2つ以上の時間点において測定される患者の組織のボクセルのパラメータと関連付けられた測定情報(例えば、2つ以上の時間点における組織のボクセルのパラメータの値、2つ以上の時間点におけるパラメータを含むか、これを示す組織の画像またはスキャンなど)を獲得し得る。いくつかの非限定的な実施形態または態様において、2つ以上の時間点は、組織のボクセルの1つまたは複数の特性が、2つ以上の時間点のうちの単一の時間点におけるボクセルのパラメータに基づいて生成される画像内で、分離可能である(例えば、測定データおよび/または測定データから作成される画像内の1つまたは複数の他の特性および/または雑音から分離可能であるなど)、および/または識別可能である(例えば、人間の目によって識別可能であるなど)前に発生する。
いくつかの非限定的な実施形態または態様において、測定情報は、以下のパラメータ/画像取得相:自然相(例えば、造影前)、動脈相、門脈相、静脈相、平衡相のうちの1つまたは複数、肝胆道取り込みおよび排泄の場合、肝胆道相(HBP)を結果としてもたらす、造影剤注入の前、間、および/もしくは後に、異なる(例えば、別々の)時間に、もしくは連続して、および/または1つもしくは複数の相にわたる連続的な画像取得として、獲得(例えば、走査、撮像など)され得る。
いくつかの非限定的な実施形態または態様において、測定データを獲得するための画像取得または走査は、欧州特許出願第20161359.3号に説明されるような患者の循環時間に基づいて決定される時間点(例えば、最適な時間点など)において開始され得、この文献の内容全体が参照により本明細書に組み込まれる。
いくつかの非限定的な実施形態または態様において、組織のボクセルの1つまたは複数のパラメータを含む測定情報は、以下の技術のうちの少なくとも1つ:MRI取得技術(例えば、T1、T2、TWI、PD、マッピング(脂肪、鉄)、多パラメータ手法、過分極MRI、MRフィンガープリンティング、エラストグラフィなどに基づくMRI)、CT取得技術、超音波技術、SPECT技術、PET/MRI技術、PET/CT技術、別の画像診断技術、またはそれらの任意の組み合わせを使用して獲得され得る。
いくつかの非限定的な実施形態または態様において、撮像システムによって測定される組織のボクセルのパラメータは、以下のうちの少なくとも1つ:T1重み付き(T1w)、T2重み付き(T2w)、陽子密度重み付き(PDw)、拡散重み付き(DWI)、x線吸収量、T1および/またはT2緩和時間の短縮量、x線吸収量の変化、トレーサ取り込み量/代謝および放射登録、1つもしくは複数の薬物動力学的モデルパラメータ、またはそれらの任意の組み合わせを含み得る。例えば、撮像システムによって取得される画像は、パラメータが撮像システムによって測定される時間における組織のボクセルのパラメータを含み得るか、またはこれを示し得る。
図7に示されるように、ステップ706において、プロセス700は、患者の組織と関連付けられた1つまたは複数の特性を決定することを含む。例えば、画像分析および演算システム614は、2つ以上の時間点におけるボクセルのパラメータに基づいて、1つまたは複数の特性を決定し得る。例として、画像分析および演算システム614は、2つ以上の時間点におけるボクセルの1つまたは複数のパラメータを含む測定情報(および/または、2つ以上の時間点におけるボクセルのパラメータを含むか、これを示す画像)に基づいて、患者の組織のボクセルの1つまたは複数の特性を含むか、これを示す1つまたは複数の画像を生成し得る。そのような例において、画像分析および演算システム614は、2つ以上の時間点におけるボクセルのパラメータを含む測定情報を取得または測定した撮像システムによって取得されるおよび/または取得されない1つまたは複数の画像を合成する、強化する、組み合わせる、および/または置き換える/削除することができる。例えば、単一の特性が、グレースケール画像として表示され得る。別の例として、2つ以上の特性が、グレースケール画像上のカラーオーバーレイとして表示され得る。
いくつかの非限定的な実施形態または態様において、1つまたは複数の特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像)は、2つ以上の時間点のうちの1つまたは複数に対応する(例えば、これと同じ時間などの)時間点および/または時間期間について決定され得る。いくつかの非限定的な実施形態または態様において、1つまたは複数の特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像)は、2つ以上の時間点の後の(例えば、2つ以上の時間点の後に続くなどの)時間点および/または時間期間について決定され得る。
いくつかの非限定的な実施形態または態様において、患者の組織と関連付けられた特性は、以下のうちの少なくとも1つ:ボクセルの動脈内(例えば、肝臓動脈(A)内など)の造影剤の濃度、ボクセルの静脈内(例えば、肝臓静脈(V)内など)の造影剤の濃度、ボクセルの細胞内(例えば、肝細胞(P)内など)の造影剤の濃度、ボクセルの動脈、静脈、および細胞内の造影剤の濃度の総計(例えば、肝臓組織のボクセルの肝臓動脈(A)、肝臓静脈(V)、肝細胞(P)内の造影剤の濃度の総計)、ボクセルの組織空間を通る造影剤の動きと関連付けられた薬物動力学的パラメータのうちの1つもしくは複数、またはそれらの任意の組み合わせを含み得る。いくつかの非限定的な実施形態または態様において、患者の組織と関連付けられた特性は、例えば、電子密度、水素密度、T1、T2、見かけの拡散係数(ADC)、またはそれらの任意の組み合わせを含む、注入された造影剤と関連付けられない少なくとも1つの特性を含み得る。
医用撮像が始まった当初、データの取得および画像の作成は、X線が、患者を横断し、画像を作成するために現像される写真フィルムと対になった発光スクリーンによって吸収されることによって実施された。次いでそのフィルムは、診断を行う放射線科医によって読み取られた。現在の電子およびコンピュータ技術では、データの2Dまたは3Dアレイを(時間(4D)にわたって取得し、このデータを、人間が理解可能な2D、3D、または4Dデータ表現を意味する画像へと変換する異なる方式を有する複数の画像診断法が存在する。目的の患者領域内のボクセルごとに各時間点においてデータの複数のパラメータを収集する能力である第5の次元についても述べることができる。図7にさらに示されるように、例となるが非限定的なデータ分析ステップは、取得したデータまたは画像を人間が読むことができる画像へ変えるデータ分析711、例えば、2つ以上の画像を使用して以前に存在していない画像を人間が読むことができる形式で作成するデータ分析712、例えば、関与する医療専門家による確認のため、読み手に通知し得るか、推奨し得るか、または診断を行い得るデータ分析713を含み得る。コンピュータ支援診断アルゴリズムは、このデータ分析の例である。例えば、画像分析および演算システム614は、以下の技術のうちの少なくとも1つ:ウィンドウ&レベル技術、ガンマ補正技術、フィルタ逆投影技術、反復再構成技術、モデルフィッティング技術、曲線フィッティング技術の辞書モード(例えば、MRフィンガープリンティングなど)、空間的、時間的フィルタリングおよび強化技術、CADもしくはAIを使用した特徴抽出技術、フラクタル解析技術、動き補正技術、柔軟な登録技術、パラメータ定量化技術、単一画像の2D AIもしくは機械学習ベースの分析、ボクセルの時間シーケンスの2D AIもしくは機械学習ベースの分析技術、画像技術の時間シーケンスの3D AIもしくは機械学習ベースの分析技術、AIもしくは機械学習ベースの雑音除去および鮮明化技術、AIもしくは機械学習ベースのコントラスト増幅技術、またはそれらの任意の組み合わせを使用して、2つ以上の時間点におけるボクセルのパラメータを含む測定情報に基づいて、1つまたは複数の特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像、および/またはそれと関連付けられた診断情報)を決定し得る。
いくつかの非限定的な実施形態または態様において、画像分析および演算システム614(ならびに/または撮像装置制御および演算システム610)は、患者の組織のボクセルのパラメータと関連付けられた測定情報が、組織のボクセルの1つまたは複数の特性を決定することと関連付けられた測定情報のしきい値量(例えば、1つもしくは複数の特性および/または1つまたは複数の特性を含む1つもしくは複数の画像を決定するのに十分な量の情報および/またはデータ)を含むことを決定し、測定情報が測定情報のしきい値量を含むことを決定することに応答して、測定情報の取得を自動的に停止する(例えば、患者を走査または撮像するのを停止するなど)ように撮像システム(例えば、撮像装置608など)を制御し得る。いくつかの非限定的な実施形態または態様において、組織のボクセルの1つまたは複数の特性を決定することと関連付けられた測定情報のしきい値量は、任意選択的に、既定のレベルの、パラメータ微分、予測の確信、誤差の余裕を伴って、患者についての1つまたは複数の選択した、または予め決められた診断を提供することと関連付けられた情報のしきい値量を含み得る。
いくつかの非限定的な実施形態または態様において、画像分析および演算システム614は、2つ以上の時間点におけるボクセルのパラメータを含む測定情報(および/または、2つ以上の時間点におけるボクセルのパラメータを含むか、これを示す1つもしくは複数の画像)に基づいて、患者の組織のボクセルの1つまたは複数の画像を合成(例えば、1つまたは複数の合成または強化画像、1つまたは複数の非強化/造影強化/トレーサ画像、1つまたは複数のT1w/T2w/PDw画像(および/または、例えば、fat-satあり/なしで、逆パルスで、など、それらの異形)、動脈/門脈/静脈/平衡/HB相の1つまたは複数の画像、MRI連続取得の1つまたは複数の画像、1つまたは複数の非強化/造影強化/トレーサx線吸収画像、1つまたは複数の動脈/門脈/静脈/平衡/または連続取得x線吸収画像などを生成または作成する)ことができる。例として、画像分析および演算システム614は、HBPの前に発生する1つまたは複数の相において測定される2つ以上の時間点におけるボクセルのパラメータを含む測定情報から、HBPにおける1つまたは複数の画像を合成し得る。例として、画像分析および演算システム614は、肝臓特異性造影剤についてのHBPから知られる造影剤情報を、肝臓の初期相の1つまたは複数の画像(例えば、動脈、門脈、および/または静脈相などにおいて測定または取得される1つまたは複数の画像など)と比較する(例えば、追加する、減じる、乗じるなど)ことによって1つまたは複数の画像を合成して、2つもしくは複数の画像からの情報全体、または2つもしくは複数の画像からの情報の一部のみを組み合わせて新規合成画像にし得る、図9に示され、2019年9月18日に提出された欧州特許出願第19197986.3号、および2020年9月10日に提出された国際特許出願第PCT/EP2020/075288号に詳細に示されるような、“無地の白い肝臓”を生成し得、これらの文献の各々の内容全体が参照により本明細書に組み込まれる。例として、画像分析および演算システム614は、図10に示され、2019年10月8日に提出された欧州特許出願第19201919.8号にさらに詳細に説明されるような1つまたは複数の測定された造影強化画像(例えば、肝臓内など)から、非強化T1w画像を合成し得、これらの文献の内容全体が参照により本明細書に組み込まれる。例として、画像分析および演算システム614は、209年10月8日に提出された欧州特許出願第19201937.0号、および2019年12月5日に提出された米国仮特許出願第62/943,980号に説明されるように、低減された量の造影剤で組織および腫瘍の高コントラスト高精度の撮像を可能にするために特定の画像造影(例えば、造影強化されたT1w MR画像内の動脈相強化の造影を強化し得、これらの文献の内容全体が参照により本明細書に組み込まれる。例として、画像分析および演算システム614は、2020年4月3日に提出された欧州特許出願第20167879.4号に説明されるように、平衡相の造影強化されたT1w画像(例えば、仮想血液プール造影画像など)を合成し得、この文献の内容全体が参照により本明細書に組み込まれる。例として、画像分析および演算システム614は、従来のガドキセト酸強化画像から、遷移相または後期相の間の肝実質細胞内のガドキセト酸誘発型造影変化なしに、画像を合成し得る。例として、画像分析および演算システム614は、物質(例えば、造影剤)が組織を通過する間の薬物動力学的時間フレームにわたる組織の完全な緩和曲線に関する情報を含む、従来通りに取得されたマッピングシーケンスのセットに基づいて、“疑似画像”を合成し得る。
例として、造影剤の二重注入プロトコルは、図11に示され、2019年10月8日に提出された欧州特許出願第19201919.8号、2019年10月11日に提出された欧州特許出願第19202711.8号、および2020年5月6日に提出された欧州特許出願第20173089.2号にさらに詳細に説明されるように、使用され得、これらの文献の各々の内容全体が参照により本明細書に組み込まれる。二重注入の第1の注入は、患者がスキャナ内にいない間に実施され得る。第1の注入は、試験がMRI試験となる予定の場合、例えば、ガドキセト酸であり得る。患者は、造影剤が体内で平衡し、標的組織によって取り込まれる間、スキャナの外に留まる。患者は、スキャナ内に置かれ、1つまたは複数の画像が、その平衡取り込みについて撮られ得、ならびに取り込みによる影響を受けない画像も撮られ得る。二重注入の第2の造影剤注入が次いで付与され、1つまたは複数の追加の画像が撮られる。これらの第2の注入後の画像は、例えば、平衡取り込み撮像に加えて、血管強化の可視化および/または測定を可能にするようにタイミングが計られ得る。例として、画像分析および演算システム614は、1つまたは複数の第2の注入後の画像から1つまたは複数の平衡取り込み後の画像を減じて、対応する時間点における平衡取り込みとは別個の血管強化を示す画像を導出し得る。第2の造影剤注入は、第1の造影剤注入と同じタイプの造影剤であり得るか、または既に取り込みが発生しており、撮像されていることから、それは、異なるタイプの造影剤、例えば、ガドキセト酸に対して異なる、およびいくつかの点で有利な、性質を有するガドブトロールであり得る。いくつかの非限定的な実施形態または態様において、2つの異なる造影剤は、混合されるか、または同時に注入され得る。例えば、造影剤を血管造影剤として好ましいものにする性質を有する造影剤は、血管相造影剤としてはそれほど有用ではないが、細胞または組織特異性マーカまたは造影剤としての性質を有する造影剤と混合され得る。
図7に示されるように、ステップ708において、プロセス700は、患者の組織と関連付けられた1つまたは複数の特性を分析することを含む。いくつかの非限定的な実施形態または態様において、放射線科医および/または他の訓練を受けた専門家は、患者の組織と関連付けられた1つまたは複数の特性を含む1つまたは複数の画像を読み取って、1つまたは複数の画像から診断情報を導出し得る。いくつかの非限定的な実施形態または態様において、画像分析および演算システム614は、患者の組織のボクセルと関連付けられた1つもしくは複数の特性、および/または1つもしくは複数の画像、患者の組織と関連付けられた1つもしくは複数の特性を分析して、1つまたは複数の画像から診断情報を導出し得る。例えば、画像分析および演算システム614は、もう1つの特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像)を、分類モデルに供給し得、この分類モデルは、もう1つの特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像に基づいて、1つまたは複数の診断と関連付けられた組織を分類するように、教師あり学習によって訓練されている。例として、画像分析および演算システム614は、組織の1つまたは複数の特性を決定するために使用されるステップ706に関して上に本明細書で説明される技術のうちの1つまたは複数を使用して、組織のための1つまたは複数の診断を決定し得る。
いくつかの非限定的な実施形態または態様において、診断情報は、患者の組織のボクセルの状態の識別(例えば、健康な組織、不健康な組織、不健康な組織のタイプ、組織と関連付けられた病気のタイプなど)、患者の組織のボクセルの状態の識別についての1つもしくは複数の理由、またはそれらの任意の組み合わせを含み得る。
図7に示されるように、ステップ710において、プロセス700は、患者の組織と関連付けられた診断を提供することを含む。例えば、放射線科医および/または他の訓練された専門家は、患者の組織と関連付けられた診断情報を提供し得る。例として、画像分析および演算システム614は、患者の組織と関連付けられた診断情報を提供し得る。
依然として図7を参照すると、いくつかの非限定的な実施形態または態様において、ステップ720で、患者および/または患者の健康状態に関する追加情報が、ステージ702~710のうちの1つまたは複数において入力および使用されて、各ステージにおけるオプション間の選択ならびに/または各ステージで決定および/もしくは生成される情報および/もしくはデータの決定に、注視し得るか、またはこれらに影響を及ぼし得る。例えば、患者および/または患者の健康状態に関する追加情報は、注入プロトコル、撮像プロトコル、測定情報、1つまたは複数の特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像)、診断情報、および/または診断を決定するために使用される、プロセス、アルゴリズム、機械学習モデルもしくはニューラルネットワーク、モデルフィッティング分析、および/または同様のものに入力され得る。例として、患者および/もしくは患者の健康状態ならびに/または試験のための所望の条件に関する追加情報は、以下のうちの少なくとも1つ:患者の身長、患者の体重、患者の年齢、患者の性別、患者の心拍数、患者の心拍出量、患者の臨床症状、患者のビリルビンレベル、患者に送達される造影剤の濃度の所望の変化率(rate)、上昇、および/もしくはプラトーレベル、またはそれらの任意の組み合わせを含み得る。患者固有のデータの追加の例は、米国特許第5,840,026号および同第9,616,166号に見ることができ、これらの文献の各々の内容全体が参照により本明細書に組み込まれる。
データ分析ステップ711、712、713、ならびに関連したコンピュータハードウェアおよび/またはソフトウェア、例えば、画像分析および演算システム614はまた、例えば、病院情報システム620および/または他のクラウドコンピューティングデータストア622から獲得され得る患者に関する追加データ720を分析および使用し得る。
これより図8を参照すると、図8は、撮像されている患者の器官のボクセル800(例えば、体積要素など)の簡略概略図を示す。この概略図は、肝臓を表しているが、一般に、適切な適合または簡略化を伴って他の器官に当てはまる。肝臓は、医学および生理学においては良く知られているように、比較的複雑な器官である。
例となる要素の各々は、ボクセルのいくらかの断片的な体積を占有し得る。肝臓は、肝門脈304および肝動脈302からの二重血液供給を受ける。肝門脈304は、肝臓の血液供給の約75%を送達し、脾臓、消化管、およびその関連器官から排出される静脈血を運ぶ。肝静脈306は、血液を心臓へ戻す。各ボクセルには、細胞が存在する。いくつかの細胞311は、撮像に使用されている造影剤を処理するか、または取り込み得る。これらの細胞311は、造影剤を胆管310内へ移動させ得る。ボクセルは、造影剤を効果的に取り込まない他の細胞313を含み得る。すべてのこれらの他の構成要素を適所に保持する流体および結合分子を表す細胞外空間301も存在する。当然ながら、器官のいくつかの部分には、完全に1つのタイプの構成要素、例えば、動脈、静脈、または胆管内にあるボクセルが存在し得る。他のボクセルは、ボクセル要素の他の断片を有する。
分子は、血液から細胞外体積内へ、およびそこから細胞または胆管内へ移動し得、その逆もまた然りである。異なる分子は、それらの特性および関与する構造体または細胞の特性に応じて異なる速度で、拡散するか、または輸送される。PK/PDモデリング(薬物動力学的/薬力学的モデリング)は、薬物動力学および薬力学の2つの古典的な薬理学的規律を組み合わせた技術である。それは、薬物動力学的および薬力学的モデル構成要素を、ある投薬量の投与に応答した効果強度の時間的経過の説明を可能にする数学的表現の1つのセットへと統合する。単純なPK/PDモデルにおいては、K1AおよびK2Aは、それぞれ、肝動脈302および関連毛細血管から出る、またはそこへ入る輸送定数である。K1VおよびK2Vは、門脈304の定数である。K3およびK4は、造影剤を取り込む細胞の定数であり、K5およびK6は、胆管内への、および胆管から戻る輸送の定数である。PK/PDモデルはまた、パラメータとして、様々な区画:動脈血液、門脈血液、静脈血液、細胞外流体、細胞外基質、ならびに様々なタイプの細胞、例えば、肝実質細胞および非肝実質細胞、によって占有される各ボクセルの断片的な体積のうちの1つまたは複数を決定する。
典型的なMR撮像において、取得される画像は、放射線科医がそれらを快適に見ることができ、それらを確実に読み取って自らの診断に達することができる十分に高い品質のものである。これは一般的に、数十秒あるいは数十分の走査時間が単一の画像を作成するために必要とされることを意味する。
表示条件を満たしたMRI画像の作成を高速化するために技術が開発されている。これらの技術は、例えば、パラレル撮像、圧縮センシング、スパース撮像、および多くの他の技術を含む。これらの技術を使用して作成される画像は、例えば、本明細書に説明されるように患者の組織のボクセルの1つまたは複数の特性を決定するために、本開示の様々な非限定的な実施形態または態様において使用され得る。
再び図7を参照すると、いくつかの非限定的な実施形態または態様において、ステップ704で、より高速であるがより雑音の多い画像が、ほんの数秒または数十秒の取得時間を使用して作成され得る。故に、ボクセルごとに捕捉されるデータ(例えば、測定情報、2つ以上の時間点において測定される組織のボクセルのパラメータなど)は、人間の視覚観察者にとって、依然として不快であるか、または混乱を招き得る測定値内の何らかの雑音があるとしても、図1の曲線Sの形状を、より良好に測定または近似することができる。
いくつかの非限定的な実施形態または態様において、図7のステップ706で、より雑音の多い画像(例えば、測定情報、2つ以上の時間点において測定される組織のボクセルのパラメータなど)の時間シーケンスは、予測モデル(例えば、人工ニューラルネットワーク、本明細書に説明されるような予測モデルなど)に供給され得、予測モデルは、2つ以上の時間点において測定される組織のボクセルのパラメータを含む測定情報に基づいて、2つ以上の時間点のうちの1つの時間点に対応する時間点および/もしくは時間期間における、ならびに/または2つ以上の時間点の後の(例えば、2つ以上の時間点の後に続くなど)時間点および/もしくは時間期間における1つまたは複数の特性(および/または、1つもしくは複数の特性を含むか、これを示す1つもしくは複数の画像)を予測するように、教師あり学習によって訓練されている。
いくつかの非限定的な実施形態または態様において、図7のステップ706で、より雑音の多い画像の時間シーケンスは、様々なモデルパラメータをデータ(例えば、2つ以上の時間点において測定される組織のボクセルのパラメータなど)にフィッティングするために使用され得る。いくつかの非限定的な実施形態または態様において、図8のPK/PDモデルは、2つ以上の時間点において測定される組織のボクセルのパラメータ(例えば、ボクセルごとに測定される曲線など)にフィッティングされ得る。例えば、画像分析および演算システム614は、組織のボクセルのPK/PDモデルを、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングし、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされたPK/PDモデルに基づいて、組織のボクセルの1つまたは複数の特性を決定し得る。例として、パラメータまたは必要な入力関数は、完全に肝動脈および門脈であるボクセルを使用して測定され得る。文献内で知られている多くの曲線フィッティング法が存在する。最小二乗曲線フィッティングは、単純であるが、必ずしも最も効率的または効果的な方法ではない。本開示の非限定的な実施形態または態様における使用に有利であり得るものは、肝臓灌流を測定する方式として単純な細胞外造影(K3~K6はゼロであることを意味する)について、RSNA 2018、コースRC629Cの間にDr.Nicole Seiberlichによって議論されたように、パラメータの離散集合の各々について事前計算され、データベースまたは辞書内に置かれる理想のPK/PD曲線のセットにおいて最良適合を見つけ出すために、比較、マッチング、またはフィンガープリンティングを用いる。ボクセルごとのPK/PDパラメータは、辞書検索により見つけ出され、対応する時間および任意の将来の時間点における画像は、ボクセルごとにPK/PD方程式を進展させ、濃度を信号強度へと変換することによって推測され得、例えば、腎臓、肝臓、および血液脳関門の細胞機能に対するある特定の結論を可能にする、薬物取り込み、分散、および排泄など、経時的なPK/PD尺度は、経時的に相対的または絶対的造影剤誘導組織強化から決定され得る。例えば、画像分析および演算システム614は、パラメータについて事前計算される複数の複数のPK/PD曲線のうちの1つのPK/PD曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングし、2つ以上の時間点におけるパラメータにフィッティングされたPK/PD曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定し得る。
PK/PD分析においては、入力関数が一般的に使用される。例えば、入力関数は、画像内、例えば、肝臓に隣接する大動脈の画像内で測定され得る。米国特許出願第16/346,219号に説明されるように、造影剤の、数秒の長さの、迅速で短い点滴またはボーラスに続いて、ボーラスを腕を通じて中枢循環へ移動させるのに十分な生理食塩水を使用する、比較的狭い入力関数が好ましいと一般的に考えられており、この文献の内容全体が参照により本明細書に組み込まれる。この短い造影剤ボーラスの欠点は、造影剤ボーラスが、患者の中枢循環を通って移動する際に広がることであり、幅10~15秒のボーラスに広がる。更なる欠点は、ボーラス形状が主に患者に依存することである。代替的に、任意選択的に>15秒の、より長いボーラスが、米国特許第9,867,589号に説明されるように使用され得、この文献の内容全体が参照により本明細書に組み込まれる。より長いボーラスは、入力関数が、患者にあまり可変的に基づかず、むしろ注入持続時間によって決定されることを引き起こす。これは、モデル内で発生することが予期され得るパラメータの範囲を制限することにより、本明細書に説明されるいくつかのタイプの分析、モデリング、または曲線フィッティングにおいて好ましい場合がある。より長いボーラスの1つの欠点は、正常な動脈、門脈、または他の相画像が、時間の重複が理由で、利用可能ではないことであり得る。放射線科医に対して“標準”画像を提示することが望ましい場合、画像分析および演算システム614は、ボクセルごとにモデルの1つまたは複数の特性を分析し得、“標準画像”は、そのボクセルについてのモデルパラメータから構築され得る。このより長いボーラス手法は、検出器飽和またはパルスパイルアップおよびデッドスペース補正効果を回避するためにPET撮像において特に有益であり得る。造影レベルがより遅い速度で変化しており、タイミングがより良好に予測され得、患者の生理機能に応じてそれほど可変ではないことから、より少ないスキャンが撮られることを可能にすることによって、CTにとっても有益であり得る。問われている診断的疑問および査定されている特性に相応しい場合、Patlak分析などの簡略化されたPK/PDモデル、および/または他のそのようなモデルが使用される。
別の非限定的な実施形態または態様において、曲線Sは、選択された基底関数のセット(例えば、多項式関数のセット、ラプラス関数のセット、フーリエ関数のセットなど)によって近似されるか、または曲線Sは、選択された基底関数のセットへと分解される。例えば、基底関数は、関数空間内のある特定の基底の要素である。ベクトル空間内の全ベクトルが基底ベクトルの線形結合として表され得るのと同じように、関数空間内の全連続関数は、基底関数の線形結合として表され得る。例として、画像分析および演算システム614は、基底関数のセットを用いて、組織のボクセルの1つまたは複数の特性を表す曲線を近似し、近似された曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングし、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた近似された曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定し得る。参照によりその内容全体が本明細書に組み込まれる、H K.Thompsonらによる、“Indicator Transit Time Considered as a Gamma Variate”, Circulation Research, Volume XIV, June 1964,pp502-515の記事において説明されるように、時間曲線にわたる初回通過造影強化は、以下の方程式(1)に従ってガンマ変量曲線としてモデリングまたは近似され得、
C(t)=K(t-AT)exp(-(t-AT)/b (1)
式中、t=注入後の時間、C(t)=時間における濃度、K=定数スケール係数、AT=出現時間、およびa,b=t>ATについての任意パラメータである。
曲線Sは、A曲線の初回通過部分、V曲線の初回通過部分を表すために2つ以上のガンマ変量曲線によって、定常再循環および再分布を表すために任意選択的に第3のガンマ変量曲線、および1つの線形に増加する曲線C(t)=m(t-AT’)+nを用いて、近似され得、式中、mおよびnは、任意パラメータであり、AT’は、線形上昇の開始であり、tは、t>AT’についての注入後の時間である。しかしながら、非限定的な実施形態または態様は、モデリングのためにガンマ変量曲線を使用することに限定されず、基底関数の任意の他のセット(例えば、多項式関数のセット、ラプラス関数のセット、フーリエ関数のセットなど)がモデリングのために使用され得、任意の計算効率の良い曲線フィッティングプログラムが、測定された曲線 Sに対する最良適合パラメータを決定するために使用され得る。ボクセルごとの基底関数パラメータが、曲線フィッティングプロセスを通じて見つけ出されると、対応する時間および任意の時間点における画像は、ボクセルごとに方程式を進展または後退させて、所望される時間tについて濃度を信号強度へと変換することによって推測され得る。例えば、画像分析および演算システム614は、基底関数のセットを用いて1つまたは複数のパラメータについて事前計算される複数の曲線のうちの1つの曲線を、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングし、2つ以上の時間点において測定される組織のボクセルのパラメータにフィッティングされた曲線に基づいて、組織のボクセルの1つまたは複数の特性を決定し得る。
本明細書内の他の場所に説明されるように、より長い造影剤ボーラスまたは注入は、曲線Sがより遅い上昇および降下を有することを引き起こし得、このことが、モデル内で予期される基底関数パラメータに影響を及ぼすか、またはこれを制約し得るため、曲線フィッティング活動および/または結果として生じるモデルを、簡略化する、高速化する、またはその正確性を増大し得る。
モデルまたはアルゴリズム、例えば、PK/PD、基底関数、AIまたは当該技術分野において当業者に知られる他のモデルが、ボクセルごとの2つ以上の時間点にフィッティングされると、画像が作成され得るか、または物理的状況において決して独立して測定されることがない情報が導出され得る。図9の“白い肝臓”は、これの非限定的な例である。別の非限定的な実施形態または態様は、造影剤が門脈だけを通って肝臓に入ったのであり、門動脈を通って入らなかった、またはその逆であるかのように、経時的に肝臓を通って流れる造影剤を示すために作成される画像または画像のシーケンスである。同様に、より長いボーラスが実際には撮像プロトコル全体において使用されていたとしても、非常に高速または短いボーラスの使用により発生した画像が作成され得る。使用される実際の撮像プロトコルを考慮すると測定することが物理的に不可能であった画像を作成するこのような能力は、診断、微妙な差異もしくは健康状態、または疾病の治療に対する処置および応答における新たな理解をもたらし得る。これは、PK/PDパラメータ自体の現在行われる画像に加えて、例えば血液量通過または様々なKもしくは複合定数マップである。
依然として図7を参照すると、いくつかの非限定的な実施形態または態様において、ステップ708で、画像分析および演算システム614は、患者の組織のボクセルと関連付けられた1つもしくは複数の特性、および/または1つもしくは複数の画像患者の組織と関連付けられた1つもしくは複数の特性を分析して、動脈相強化(例えば、肝臓)と比較して、経時的に病変のウォッシュアウト(washout)についての絶対値を計算し得る。現在、病変のウォッシュアウトは通常、動脈相の後に続く造影相の間、周辺組織の強化を病変の強化と比較することによって、視覚的に査定される。しかしながら、これは、いくつかの造影剤(肝臓特異性薬剤Primovist(登録商標))が、病変(例えば、潜在的に転移または他の腫瘍)の部分を囲む、または病変の部分である肝細胞によって取り込まれ得ることから、誤った解釈のリスクを明るみにする。Primovist(登録商標)のこのような初期取り込みは、病変と周囲組織との間のより強いコントラストの錯覚を作り出し得、実際には何も存在しないときにウォッシュアウトとして誤って解釈され得る。これは、本開示の非限定的な実施形態または態様の、これらの態様の別個または別々の分解および/または可視化を可能にする分析を伴わない、誤診断のリスクを明るみにする。
いくつかの非限定的な実施形態または態様において、検査領域は、哺乳動物(好ましくは、ヒト)の肝臓または肝臓の一部分を含む。いくつかの非限定的な実施形態または態様において、検査領域は、肺またはその一部分を含み得る。肺は、右心から非酸素化血液、また左心から酸素化血液の循環を受容する。加えて、肺は、気道を通じてガスを受容する。故に、肺は、静脈内、または気道を通じたガス状造影剤のいずれかで、造影剤を受容し得る。したがって、本開示の非限定的な実施形態または態様は、吸入造影剤、ならびにIV注入造影剤に当てはまり得る。加えて、大半の組織は、組織の細胞外空間を通る流体の循環のためにリンパ系(例えば、脳内のグリンパティック系)として作用する。本開示の非限定的な実施形態または態様は、このリンパ系およびグリンパティック系を通る造影剤の流れ、またはその欠如を含み得る。
本明細書に説明される様々な非限定的な実施形態または態様において、様々な組織特性およびボクセルパラメータは、理解および開示のために列挙されている。これらは、例示的であり、限定的または徹底的なリストではないということを理解されたい。医術において知られる他の特性および/またはパラメータが使用され得る。加えて、研究中または今後発見される特性および/またはパラメータもまた、撮像研究および/または一連の研究からのデータの処理、撮像、および/または分析において、本開示の非限定的な実施形態または態様の適用から恩恵を受け得る。

Claims (8)

  1. 者の組織の画像のボクセルのパラメータと関連付けられた測定情報を獲得するステップであって、前記測定情報は、2つ以上の時間点のうちの第1時間点における前記ボクセルの前記パラメータと関連付けられた第1測定情報と、前記2つ以上の時間点のうちの第2時間点における前記ボクセルの前記パラメータと関連付けられた第2測定情報とを提供するように、前記2つ以上の時間点において測定され、前記2つ以上の時間点は、記ボクセルの1つまたは複数の組織特性が、1つの時間点において測定される前記ボクセルの前記パラメータに基づいて生成される別の画像内で分離可能または識別可能である前に発生する、ステップを含み
    前記1つまたは複数の組織特性は、
    動脈内の造影剤の濃度、
    静脈内の造影剤の濃度、
    細胞内の造影剤の濃度、
    動脈、静脈、および細胞内の造影剤の濃度の総計、
    組織空間を通る造影剤の動きと関連付けられた1つまたは複数の薬物動力学的モデルパラメータ、または
    これらの任意の組み合わせを含み、
    前記第1時間点における前記ボクセルの前記パラメータと関連付けられた前記第1測定情報と、前記第2時間点における前記ボクセルの前記パラメータと関連付けられた前記第2測定情報と、前記患者に送達される造影剤の濃度の所望の変化率および/またはプラトーレベルとに基づいて、前記組織の画像の前記ボクセルの前記1つまたは複数の特性を決定するステップをさらに含む
    コンピュータ実施方法。
  2. 前記組織の画像の前記ボクセルの前記1つまたは複数の特性は、前記2つ以上の時間点の後の時間点について決定される、請求項1に記載のコンピュータ実施方法。
  3. 前記1つまたは複数の組織特性に基づいて、前記2つ以上の時間点の後の前記時間点における前記組織の画像の前記ボクセルの1つまたは複数の特性を含む1つまたは複数の画像を生成するステップをさらに含む、請求項2に記載のコンピュータ実施方法。
  4. 前記組織の画像の前記ボクセルの前記1つまたは複数の特性を決定するステップは、
    前記患者の前記組織の画像の前記ボクセルの前記パラメータと関連付けられた前記測定情報を予測モデルに供給することを含み、前記予測モデルは、前記2つ以上の時間点における前記パラメータと関連付けられた前記測定情報に基づいて、前記組織の画像の前記ボクセルの前記1つまたは複数の特性を予測するように、教師あり学習によって訓練されている、請求項2に記載のコンピュータ実施方法。
  5. プログラム命令を含む少なくとも1つの非一時的なコンピュータ可読媒体を備えるコンピュータプログラム製品であって
    前記プログラム命令は、少なくとも1つのプロセッサによって実行されるとき、前記少なくとも1つのプロセッサに、
    者の組織の画像のボクセルのパラメータと関連付けられた測定情報を獲得することであって、前記測定情報は、2つ以上の時間点のうちの第1時間点における前記ボクセルの前記パラメータと関連付けられた第1測定情報と、前記2つ以上の時間点のうちの第2時間点における前記ボクセルの前記パラメータと関連付けられた第2測定情報とを提供するように、前記2つ以上の時間点において測定され、前記2つ以上の時間点は、記ボクセルの1つまたは複数の組織特性が、1つの時間点において測定される前記ボクセルの前記パラメータに基づいて生成される別の画像内で分離可能または識別可能である前に発生する、獲得することを行わせ
    前記1つまたは複数の組織特性は、
    動脈内の造影剤の濃度、
    静脈内の造影剤の濃度、
    細胞内の造影剤の濃度、
    動脈、静脈、および細胞内の造影剤の濃度の総計、
    組織空間を通る造影剤の動きと関連付けられた1つまたは複数の薬物動力学的モデルパラメータ、または
    これらの任意の組み合わせを含み、
    前記プログラム命令は、前記少なくとも1つのプロセッサによって実行されるとき、前記少なくとも1つのプロセッサに、
    前記第1時間点における前記ボクセルの前記パラメータと関連付けられた前記第1測定情報と、前記第2時間点における前記ボクセルの前記パラメータと関連付けられた前記第2測定情報と、前記患者に送達される造影剤の濃度の所望の変化率および/またはプラトーレベルとに基づいて、前記組織の画像の前記ボクセルの前記1つまたは複数の特性を決定することをさらに行わせる、
    コンピュータプログラム製品。
  6. 前記組織の画像の前記ボクセルの前記1つまたは複数の特性は、前記2つ以上の時間点の後の時間点について決定される、請求項5に記載のコンピュータプログラム製品。
  7. 前記命令は、前記少なくとも1つのプロセッサに、
    前記1つまたは複数の組織特性に基づいて、前記2つ以上の時間点の後の前記時間点における前記組織の画像の前記ボクセルの1つまたは複数の特性を含む1つまたは複数の画像を生成することをさらに行わせる、請求項6に記載のコンピュータプログラム製品。
  8. 前記命令は、前記少なくとも1つのプロセッサに、
    前記患者の前記組織の画像の前記ボクセルの前記パラメータと関連付けられた前記測定情報を予測モデルに供給することによって、前記組織の画像の前記ボクセルの前記1つまたは複数の特性を決定することを行わせ、前記予測モデルは、前記2つ以上の時間点における前記パラメータと関連付けられた前記測定情報に基づいて、前記組織の画像の前記ボクセルの前記1つまたは複数の特性を予測するように、教師あり学習によって訓練されている、請求項6に記載のコンピュータプログラム製品。
JP2022517507A 2019-09-18 2020-09-17 組織特性を予測、予想、および/または査定するためのシステム、方法、およびコンピュータプログラム製品 Active JP7535575B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
EP19197989.7 2019-09-18
EP19197989 2019-09-18
US201962943969P 2019-12-05 2019-12-05
US62/943,969 2019-12-05
PCT/IB2020/058688 WO2021053585A1 (en) 2019-09-18 2020-09-17 System, method, and computer program product for predicting, anticipating, and/or assessing tissue characteristics

Publications (4)

Publication Number Publication Date
JP2022548930A JP2022548930A (ja) 2022-11-22
JPWO2021053585A5 JPWO2021053585A5 (ja) 2024-06-18
JP2022548930A5 JP2022548930A5 (ja) 2024-06-18
JP7535575B2 true JP7535575B2 (ja) 2024-08-16

Family

ID=72521684

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022517507A Active JP7535575B2 (ja) 2019-09-18 2020-09-17 組織特性を予測、予想、および/または査定するためのシステム、方法、およびコンピュータプログラム製品

Country Status (7)

Country Link
US (1) US11915361B2 (ja)
EP (1) EP4031895B1 (ja)
JP (1) JP7535575B2 (ja)
CN (2) CN117582208A (ja)
AU (1) AU2020349121B2 (ja)
CA (1) CA3154689A1 (ja)
WO (1) WO2021053585A1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11922314B1 (en) * 2018-11-30 2024-03-05 Ansys, Inc. Systems and methods for building dynamic reduced order physical models
WO2022076775A1 (en) * 2020-10-09 2022-04-14 Massachusetts Institute Of Technology Methods and systems for assessment of liver fibrosis and steatosis
US20230213604A9 (en) * 2020-12-15 2023-07-06 Washington University Diffusion dictionary imaging (ddi) of microstructure and inflammation
CN112541876B (zh) * 2020-12-15 2023-08-04 北京百度网讯科技有限公司 卫星图像处理方法、网络训练方法、相关装置及电子设备
EP4156100A1 (en) * 2021-09-23 2023-03-29 Siemens Healthcare GmbH Providing result image data
US20230104028A1 (en) * 2021-10-05 2023-04-06 Hitachi, Ltd. System for failure prediction for industrial systems with scarce failures and sensor time series of arbitrary granularity using functional generative adversarial networks
EP4343786A1 (en) * 2022-09-20 2024-03-27 Siemens Healthineers AG Building a machine-learning model to predict sematic context information for contrast-enhanced medical imaging measurements
WO2024148577A1 (zh) * 2023-01-12 2024-07-18 深圳先进技术研究院 动态pet成像方法、装置、电子设备及存储介质
WO2025053243A1 (ja) * 2023-09-07 2025-03-13 武田薬品工業株式会社 3次元画像予測装置、3次元画像学習装置、3次元画像予測方法、3次元画像学習方法、3次元画像予測プログラム及び3次元画像学習プログラム
CN119090844A (zh) * 2024-08-28 2024-12-06 江苏省人民医院(南京医科大学第一附属医院) 一种基于ct和超声影像构建纳米颗粒肿瘤积聚预测模型的系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002165775A (ja) 2000-12-01 2002-06-11 Hitachi Medical Corp 磁気共鳴イメージング装置
US20100198054A1 (en) 2007-05-18 2010-08-05 Ewing James R Mri estimation of contrast agent concentration using a neural network approach
US20190108634A1 (en) 2017-10-09 2019-04-11 The Board Of Trustees Of The Leland Stanford Junior University Contrast Dose Reduction for Medical Imaging Using Deep Learning

Family Cites Families (112)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US732697A (en) 1903-02-19 1903-07-07 John Davis Bates Lumber-surface meter.
US6039931A (en) 1989-06-30 2000-03-21 Schering Aktiengesellschaft Derivatized DTPA complexes, pharmaceutical agents containing these compounds, their use, and processes for their production
US5840026A (en) 1994-09-21 1998-11-24 Medrad, Inc. Patient specific dosing contrast delivery systems and methods
US6397098B1 (en) 1994-09-21 2002-05-28 Medrad, Inc. Data communication and control for medical imaging systems
US5732697A (en) 1995-11-22 1998-03-31 Arch Development Corporation Shift-invariant artificial neural network for computerized detection of clustered microcalcifications in mammography
US6339718B1 (en) 1999-07-30 2002-01-15 Medrad, Inc. Programmable injector control
US6754376B1 (en) 2000-11-22 2004-06-22 General Electric Company Method for automatic segmentation of medical images
US6819790B2 (en) 2002-04-12 2004-11-16 The University Of Chicago Massive training artificial neural network (MTANN) for detecting abnormalities in medical images
US7458936B2 (en) 2003-03-12 2008-12-02 Siemens Medical Solutions Usa, Inc. System and method for performing probabilistic classification and decision support using multidimensional medical image databases
US7545965B2 (en) 2003-11-10 2009-06-09 The University Of Chicago Image modification and detection using massive training artificial neural networks (MTANN)
US20060018524A1 (en) 2004-07-15 2006-01-26 Uc Tech Computerized scheme for distinction between benign and malignant nodules in thoracic low-dose CT
CN102488527B (zh) 2004-11-16 2015-06-10 拜耳医疗保健公司 从成像过程确定病人生理参数的系统和方法
US7738683B2 (en) 2005-07-22 2010-06-15 Carestream Health, Inc. Abnormality detection in medical images
US7564990B2 (en) 2005-08-18 2009-07-21 Nu Skin International, Inc. Imaging system and method for physical feature analysis
US7343032B2 (en) 2005-09-01 2008-03-11 Fujifilm Corporation Method and apparatus for automatic and dynamic vessel detection
EP1780677A1 (en) 2005-10-25 2007-05-02 BRACCO IMAGING S.p.A. Image processing system, particularly for use with diagnostics images
EP1952347B1 (en) 2005-11-01 2019-03-06 Edda Technology, Inc. Method and system for liver lobe segmentation and pre-operative surgical planning
US20220031270A1 (en) * 2007-03-08 2022-02-03 Sync-Rx, Ltd Identification an dpresentation of device-tovessel relative motion
DE102007028660B3 (de) 2007-06-21 2009-01-29 Siemens Ag Verfahren zur Korrektur von Bewegungsartefakten bei der Aufnahme von MR-Bildern
CN101938939B (zh) 2008-02-07 2014-11-26 皇家飞利浦电子股份有限公司 用于测量和预测患者呼吸稳定性的装置
WO2009135923A1 (en) 2008-05-08 2009-11-12 Oslo Universitetssykehus Hf Vessel segmentation in dce mr imaging
WO2009146388A1 (en) * 2008-05-28 2009-12-03 The Trustees Of Columbia University In The City Of New York Voxel-based methods for assessing subjects using positron emission tomography
US7949167B2 (en) 2008-06-12 2011-05-24 Siemens Medical Solutions Usa, Inc. Automatic learning of image features to predict disease
GB0919269D0 (en) * 2009-11-03 2009-12-16 Univ Manchester Lung function analysis method and apparatus
EP2528493B1 (en) 2010-01-28 2017-04-19 Brainlab AG Method and apparatus for estimating temperature in a body
US8478698B1 (en) 2010-03-17 2013-07-02 James Mah Methods and systems for employing artificial intelligence in automated orthodontic diagnosis and treatment planning
CN101912262B (zh) * 2010-07-22 2012-07-04 中国科学院深圳先进技术研究院 磁共振成像定量参数计算方法及系统
DK3281586T3 (da) 2010-12-08 2019-06-24 Bayer Healthcare Llc Generering af et estimat for patientstrålingsdosis som følge af medicinske billeddannelsesscanninger
US9436989B2 (en) 2011-06-03 2016-09-06 Bayer Healthcare Llc System and method for rapid quantitative dynamic molecular imaging scans
FR2979453B1 (fr) 2011-08-26 2016-01-08 Olea Medical Systeme et procede pour estimer une quantite d'interet d'un systeme dynamique artere/tissu/veine
WO2013036842A2 (en) 2011-09-08 2013-03-14 Radlogics, Inc. Methods and systems for analyzing and reporting medical images
JP5878009B2 (ja) 2011-12-12 2016-03-08 株式会社根本杏林堂 医療用画像処理ネットワークシステム
WO2016007734A1 (en) 2014-07-09 2016-01-14 The Trustees Of Dartmouth College Methods for quantitative and enhanced-contrast molecular medical imaging using cross-modality correction for differing tracer kinetics
EP2626718A1 (en) 2012-02-09 2013-08-14 Koninklijke Philips Electronics N.V. MRI with motion correction using navigators acquired using a Dixon technique
WO2013121374A2 (en) 2012-02-17 2013-08-22 Koninklijke Philips N.V. Acute lung injury (ali)/acute respiratory distress syndrome (ards) assessment and monitoring
US20150125398A1 (en) 2012-05-11 2015-05-07 University Of Iowa Research Foundation Multimodal imaging methods using mesoporous silica nanoparticles
TW201408262A (en) * 2012-08-30 2014-03-01 Univ Chang Gung Method of extracting arterial input function and its application for dynamic contrast enhanced magnetic resonance imaging (DCE-MRI)
DE102012215718B4 (de) 2012-09-05 2022-05-12 Siemens Healthcare Gmbh Verfahren und Magnetresonanzanlage zur MR-Bildgebung eines vorbestimmten Volumenabschnitts eines lebenden Untersuchungsobjekts mittels Stimulieren des Untersuchungsobjekts
US9449381B2 (en) 2012-09-10 2016-09-20 Arizona Board Of Regents, A Body Corporate Of The State Of Arizona, Acting For And On Behalf Of Arizona State University Methods, systems, and media for generating and analyzing medical images having elongated structures
WO2014042902A1 (en) * 2012-09-13 2014-03-20 The Regents Of The University Of California Lung, lobe, and fissure imaging systems and methods
US20140094679A1 (en) * 2012-10-03 2014-04-03 General Electric Company Systems and methods for performing organ detection
RU2015120325A (ru) * 2012-10-31 2016-12-20 Конинклейке Филипс Н.В. Визуализация перфузии
EP2750102B2 (en) 2012-12-27 2025-12-10 General Electric Company Method, system and computer readable medium for liver analysis
EP2964270A4 (en) 2013-02-21 2017-03-22 University Of Rochester Methods for evaluating brain-wide paravascular pathway for waste clearance function and methods for treating neurodegenerative disorders based thereon
JP6530743B2 (ja) 2013-04-03 2019-06-12 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 血管セグメント化
EP2988668B1 (en) 2013-04-24 2019-07-24 Tel HaShomer Medical Research Infrastructure and Services Ltd. Magnetic resonance maps for analyzing tissue
EP3052013A4 (en) * 2013-10-04 2017-07-05 Howard University Voxel-resolution myocardial blood flow analysis
US9754371B2 (en) 2014-07-31 2017-09-05 California Institute Of Technology Multi modality brain mapping system (MBMS) using artificial intelligence and pattern recognition
US20160038092A1 (en) 2014-08-11 2016-02-11 Douglas A. Golay Applying non-real time and non-user attended algorithms to stored non-imaging data and existing imaging data for obtaining a dental diagnosis
US10061003B2 (en) 2014-09-01 2018-08-28 bioProtonics, L.L.C. Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms
US10176573B2 (en) * 2014-10-13 2019-01-08 Agency For Science, Technology And Research Automatic region-of-interest segmentation and registration of dynamic contrast-enhanced images of colorectal tumors
JP6740136B2 (ja) * 2014-11-27 2020-08-12 国立大学法人広島大学 シミュレータ、該シミュレータを備える注入装置又は撮像システム、及びシミュレーションプログラム
JP2018511791A (ja) 2015-02-27 2018-04-26 バイエル・ヘルスケア・エルエルシーBayer HealthCare LLC マルチイメージングモダリティ用定量的ファントム
US9959615B2 (en) 2015-07-01 2018-05-01 Arizona Board Of Regents On Behalf Of Arizona State University System and method for automatic pulmonary embolism detection
US20170016972A1 (en) 2015-07-13 2017-01-19 Siemens Medical Solutions Usa, Inc. Fast Prospective Motion Correction For MR Imaging
CN107850652B (zh) 2015-07-15 2021-10-15 皇家飞利浦有限公司 具有运动检测的mr成像
US10157467B2 (en) 2015-08-07 2018-12-18 Arizona Board Of Regents On Behalf Of Arizona State University System and method for detecting central pulmonary embolism in CT pulmonary angiography images
US10176408B2 (en) 2015-08-14 2019-01-08 Elucid Bioimaging Inc. Systems and methods for analyzing pathologies utilizing quantitative imaging
US20180242917A1 (en) * 2015-08-21 2018-08-30 Duke University Systems and methods for determining image quality and models based on contrast injection protocols
KR101728046B1 (ko) * 2015-08-27 2017-04-18 삼성전자주식회사 단층 영상 복원 장치 및 그에 따른 단층 영상 복원 방법
ES2952520T3 (es) 2015-08-28 2023-10-31 Bayer Healthcare Llc Sistema y procedimiento para la verificación del llenado del fluido de la jeringa y el reconocimiento de imágenes de las características del sistema inyector de energía
US10645359B2 (en) 2015-12-16 2020-05-05 B<>Com Method for processing a digital image, device, terminal equipment and associated computer program
US9786058B2 (en) 2016-02-08 2017-10-10 Sony Corporation Method and system for segmentation of vascular structure in a volumetric image dataset
WO2017139653A1 (en) 2016-02-10 2017-08-17 The Charles Stark Draper Laboratory, Inc. Systems and methods for imaging
DE102016204198B4 (de) 2016-03-15 2018-06-07 Siemens Healthcare Gmbh Verfahren zur Erzeugung von MR-Bildern mit prospektiver Bewegungskorrektur und teilvolumenspezifischer Gewichtung der Bildinformation
US10282588B2 (en) 2016-06-09 2019-05-07 Siemens Healthcare Gmbh Image-based tumor phenotyping with machine learning from synthetic data
US10650512B2 (en) * 2016-06-14 2020-05-12 The Regents Of The University Of Michigan Systems and methods for topographical characterization of medical image data
EP3490454A4 (en) * 2016-07-29 2020-08-05 Novadaq Technologies ULC METHODS AND SYSTEMS FOR CHARACTERIZING A SUBJECT FABRIC USING MACHINE LEARNING
WO2018048507A1 (en) 2016-09-06 2018-03-15 Han Xiao Neural network for generating synthetic medical images
CN109690577A (zh) 2016-09-07 2019-04-26 皇家飞利浦有限公司 利用堆叠式自动编码器进行的半监督式分类
WO2018075386A1 (en) 2016-10-17 2018-04-26 Bayer Healthcare Llc Fluid injector with syringe engagement mechanism
JP2019535424A (ja) 2016-11-22 2019-12-12 ハイパーファイン リサーチ,インコーポレイテッド 磁気共鳴画像における自動検出のためのシステムおよび方法
EP3544655B1 (en) 2016-11-22 2025-01-15 Bayer Healthcare LLC System and method for delivering a fluid with a consistent total volumetric flowrate
US10096109B1 (en) 2017-03-31 2018-10-09 The Board Of Trustees Of The Leland Stanford Junior University Quality of medical images using multi-contrast and deep learning
CN110753935A (zh) 2017-04-25 2020-02-04 小利兰·斯坦福大学托管委员会 用于医学成像的使用深度卷积神经网络的剂量减少
US20180315183A1 (en) * 2017-04-28 2018-11-01 General Electric Company System and method for monitoring an amount of a contrast agent within an object
EP3619631A1 (de) 2017-05-02 2020-03-11 Bayer Aktiengesellschaft Verbesserungen bei der radiologischen erkennung chronisch thromboembolischer pulmonaler hypertonie
JP6920904B2 (ja) 2017-07-04 2021-08-18 国立大学法人広島大学 シミュレータ、該シミュレータを備える注入装置又は撮像システム、及びシミュレーションプログラム
CN107316334B (zh) * 2017-07-31 2020-02-14 华东师范大学 个性化精准磁共振影像方法
EP3676854A1 (en) 2017-08-31 2020-07-08 Bayer Healthcare LLC Fluid path impedance assessment for improving fluid delivery performance
WO2019063520A1 (en) 2017-09-26 2019-04-04 Koninklijke Philips N.V. SYSTEM AND METHOD FOR JOINT CLINICAL DECISION FOR PHARMACEUTICAL PRODUCTS
US10335106B2 (en) 2017-09-29 2019-07-02 Infinitt Healthcare Co., Ltd. Computing system and method for identifying and visualizing cerebral thrombosis based on medical images
US11100621B2 (en) 2017-10-20 2021-08-24 Imaging Biometrics, Llc Simulated post-contrast T1-weighted magnetic resonance imaging
CA3078313A1 (en) 2017-10-31 2019-05-09 Ge Healthcare Limited Medical system for diagnosing cognitive disease pathology and/or outcome
CN116741352A (zh) 2017-11-24 2023-09-12 佳能医疗系统株式会社 医用数据处理装置、磁共振成像装置及学习完毕模型生成方法
CN108324244A (zh) 2018-01-03 2018-07-27 华东师范大学 用于ai+mri影像辅助诊断的自动增广训练样本的构建方法及系统
KR102001398B1 (ko) 2018-01-25 2019-07-18 재단법인 아산사회복지재단 기계학습을 통한 뇌질환 상태변화 예측방법, 장치 및 프로그램
US11880962B2 (en) 2018-02-15 2024-01-23 General Electric Company System and method for synthesizing magnetic resonance images
WO2019183555A1 (en) * 2018-03-22 2019-09-26 The University Of North Carolina At Chapel Hill Methods, systems and computer readable media for processing digital subtraction angiography (dsa) and computed tomography (ct) images for reducing radiation exposure in dsa and ct subjects
US11100632B2 (en) 2018-04-13 2021-08-24 Elekta, Inc. Image synthesis using adversarial networks such as for radiation therapy
US11982726B2 (en) * 2018-04-13 2024-05-14 University Of Southern California Method for joint arterial input function and tracer kinetic parameter estimation in accelerated DCE-MRI using a model consistency constraint
CN111989039B (zh) 2018-04-19 2024-11-22 长沙微妙医疗科技有限公司 使用深度学习改进磁共振成像的系统和方法
CN108703764B (zh) * 2018-05-29 2021-11-05 北京东软医疗设备有限公司 血管造影方法、装置、系统、设备及存储介质
EP3807814A4 (en) 2018-06-15 2022-03-16 Subtle Medical, Inc. MAGNETIC RESONANCE IMAGING STANDARDIZATION SYSTEMS AND METHODS USING DEEP LEARNING
EP3649955A1 (en) 2018-11-07 2020-05-13 Koninklijke Philips N.V. Deep learning based spectral bolus tracking
EP3657433B1 (en) 2018-11-22 2022-10-19 Siemens Healthcare GmbH Medical image data
DE102018222606A1 (de) 2018-12-20 2020-06-25 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Detektion eines anatomischen Merkmals eines Blutgefäßabschnittes
CN109726752A (zh) * 2018-12-25 2019-05-07 脑玺(上海)智能科技有限公司 基于时间信号曲线的灌注动态影像的分割方法及系统
US10867375B2 (en) * 2019-01-30 2020-12-15 Siemens Healthcare Gmbh Forecasting images for image processing
US11398026B2 (en) 2019-03-28 2022-07-26 The Board Of Trustees Of The Leland Stanford Junior University Systems and methods for synthetic medical image generation
DE102019207639A1 (de) * 2019-05-24 2020-11-26 Siemens Healthcare Gmbh Verbesserte Berechnung der Kontrastmittelkonzentration
EP3756547A1 (en) * 2019-06-28 2020-12-30 Koninklijke Philips N.V. Automated coronary angiography analysis
US20210012486A1 (en) 2019-07-09 2021-01-14 Shenzhen Malong Technologies Co., Ltd. Image synthesis with generative adversarial network
CN110391016B (zh) * 2019-07-10 2021-11-09 浙江大学 一种动态对比增强磁共振图像的分析方法
JP2022548716A (ja) 2019-09-18 2022-11-21 バイエル、アクチエンゲゼルシャフト 肝臓のmri画像の生成
CA3156921A1 (en) 2019-10-08 2021-04-15 Bayer Aktiengesellschaft Generation of mri images of the liver without contrast enhancement
EP3804615A1 (de) 2019-10-08 2021-04-14 Bayer AG Erzeugung von mrt-aufnahmen der leber
CA3157469A1 (en) 2019-10-11 2021-04-15 Bayer Aktiengesellschaft Acceleration of mri examinations
CN110853738B (zh) * 2019-10-12 2023-08-18 拜耳股份有限公司 一种造影剂作用下的成像方法与设备
EP3875979A1 (de) 2020-03-06 2021-09-08 Bayer Aktiengesellschaft Optimiertes verfahren für dynamische kontrastverstärkende magnetresonanztomographie
US12394058B2 (en) 2020-04-03 2025-08-19 Bayer Aktiengesellschaft Generation of radiological images
US11308613B2 (en) 2020-06-09 2022-04-19 Siemens Healthcare Gmbh Synthesis of contrast enhanced medical images

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002165775A (ja) 2000-12-01 2002-06-11 Hitachi Medical Corp 磁気共鳴イメージング装置
US20100198054A1 (en) 2007-05-18 2010-08-05 Ewing James R Mri estimation of contrast agent concentration using a neural network approach
US20190108634A1 (en) 2017-10-09 2019-04-11 The Board Of Trustees Of The Leland Stanford Junior University Contrast Dose Reduction for Medical Imaging Using Deep Learning

Also Published As

Publication number Publication date
WO2021053585A1 (en) 2021-03-25
CN113329688A (zh) 2021-08-31
JP2022548930A (ja) 2022-11-22
EP4031895B1 (en) 2026-04-22
CN113329688B (zh) 2023-11-28
AU2020349121A1 (en) 2022-03-31
US11915361B2 (en) 2024-02-27
US20220414972A1 (en) 2022-12-29
CA3154689A1 (en) 2021-03-25
CN117582208A (zh) 2024-02-23
AU2020349121B2 (en) 2026-02-05
EP4031895A1 (en) 2022-07-27

Similar Documents

Publication Publication Date Title
JP7535575B2 (ja) 組織特性を予測、予想、および/または査定するためのシステム、方法、およびコンピュータプログラム製品
JP7587575B2 (ja) コントラスト強調のない肝臓のmri画像の生成
CN113330483B (zh) 由监督学习训练的预测模型预测磁共振成像图像
JP7729836B2 (ja) 放射線画像の生成
CN104217398B (zh) 图像处理装置、图像处理方法和医学图像设备
US8774898B2 (en) Oxygen-enhanced MR imaging technique
CN114502069A (zh) 加速mrt检查
CN108283017B (zh) 在使用磁共振系统的诊断成像中优化造影剂浓度的时间分布的预先确定的方法和磁共振系统
US20240153163A1 (en) Machine learning in the field of contrast-enhanced radiology
CA3154652C (en) Forecast of mri images by means of a forecast model trained by supervised learning
Salm et al. Validation of a high-resolution, phase contrast cardiovascular magnetic resonance sequence for evaluation of flow in coronary artery bypass grafts
Seemann et al. Characterizing extravascular lung water—A dual-contrast agent extracellular volume approach by cardiovascular magnetic resonance
CN117083629A (zh) 对比度增强放射学领域中的机器学习
HK1255646B (en) Method for optimizing the predetermination of the time profile of a contrast agent concentration in diagnostic imaging using a magnetic resonance system and magnetic resonance system

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230626

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230626

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240126

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240301

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240521

A524 Written submission of copy of amendment under article 19 pct

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20240521

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240805

R150 Certificate of patent or registration of utility model

Ref document number: 7535575

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150