JP6567179B2 - 特徴回帰モデルを用いたmrデータからの疑似ct生成 - Google Patents

特徴回帰モデルを用いたmrデータからの疑似ct生成 Download PDF

Info

Publication number
JP6567179B2
JP6567179B2 JP2018519442A JP2018519442A JP6567179B2 JP 6567179 B2 JP6567179 B2 JP 6567179B2 JP 2018519442 A JP2018519442 A JP 2018519442A JP 2018519442 A JP2018519442 A JP 2018519442A JP 6567179 B2 JP6567179 B2 JP 6567179B2
Authority
JP
Japan
Prior art keywords
image
pseudo
training
generating
parameter
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
JP2018519442A
Other languages
English (en)
Other versions
JP2018530401A (ja
JP2018530401A5 (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.)
Elekta Inc
Original Assignee
Elekta Inc
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 Elekta Inc filed Critical Elekta Inc
Publication of JP2018530401A publication Critical patent/JP2018530401A/ja
Publication of JP2018530401A5 publication Critical patent/JP2018530401A5/ja
Application granted granted Critical
Publication of JP6567179B2 publication Critical patent/JP6567179B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • 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 classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0035Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
    • 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/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7278Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1039Treatment planning systems using functional images, e.g. PET or MRI
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/60Image enhancement or restoration using machine learning, e.g. neural networks
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B90/00Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
    • A61B90/36Image-producing devices or illumination devices not otherwise provided for
    • A61B90/37Surgical systems with images on a monitor during operation
    • A61B2090/376Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy
    • A61B2090/3762Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy using computed tomography systems [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • 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/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Artificial Intelligence (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Theoretical Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Radiation-Therapy Devices (AREA)

Description

(関連出願との相互参照)
[001]
本出願は、2015年10月13日に提出され、「組織パラメータ推定を用いたMRデータからの疑似CT生成」と題された代理人整理番号12475.0050−00000に関連し、その全内容は参照により本明細書に組み込まれる。
(技術分野)
[002]
本開示は、一般に、放射線治療(radiation therapy)又は放射線療法(radiotherapy)に関する。より詳しくは、本開示は、放射線治療中に使用される放射線治療計画を開発する際に使用するために、MRデータから疑似CT画像を生成するためのシステム及び方法に関する。
[003]
放射線療法は、哺乳動物(例えば、ヒト及び動物)組織における癌及び他の病気を治療するために使用される。そのような放射線治療技術のひとつはガンマナイフであり、それによって患者は標的(例えば、腫瘍)において高強度かつ高精度に収束する多数の低強度ガンマ線によって照射される。別の実施形態では、線形加速器を使用して放射線療法が提供され、それによって腫瘍が高エネルギー粒子(例えば、電子、プロトン、イオンなど)によって照射される。腫瘍が所定の放射線を確実に受けるように、放射線ビームの配置及び線量を正確に制御しなければならず、ビームの配置は周囲の健康な組織(しばしばリスク臓器(OARs)と呼ばれる)への損傷を最小限に抑えるようなものでなければならない。医学の処方箋と同様に、医師が、腫瘍及び周囲の器官に所定量の放射線を指示するので、放射線は「処方された」と呼ばれる。
[004]
従来、臨床及び線量測定の目的及び制約(例えば、腫瘍及び重要な器官への放射線の最大、最小及び平均線量)に基づく最適化技術を使用して、各患者について、放射線療法治療計画(「治療計画」)を作成することができる。治療計画手順は、標的領域(例えば、腫瘍)を同定し、腫瘍の近くの重要な器官を同定するために、患者の三次元画像を使用することを含むことができる。治療計画の作成は、計画作成者が個々の重要性(例えば、重み付け)を考慮して、臨床的に許容される治療計画を生成するために、様々な治療目的又は制約(例えば、線量量ヒストグラム(DVH)目的)に従うことを試みる時間のかかるプロセスとなり得る。この作業は、リスク臓器(OARs)の数が増える(例えば、頭頸部治療では13個まで)につれてプロセスの複雑さが増すので、様々なリスク臓器(OARs)によって複雑になり時間のかかる試行錯誤のプロセスとなり得る。腫瘍から離れたリスク臓器(OARs)は、放射線から容易に逃げることができ、一方、標的腫瘍に近いか又は重複するリスク臓器(OARs)は、逃げることが困難である。
[005]
コンピュータ断層撮影(CT)画像は、伝統的に、放射線治療の治療計画のための画像データの主要なソースとして役立つ。CT画像は患者の幾何形状の正確な表現を提供し、CT値は放射線量計算のために電子密度(例えば、ハウンスフィールド単位(Hounsfield units))に直接変換することができる。しかし、CTを使用すると、患者は追加の放射線量にさらされる。CT画像に加えて、磁気共鳴画像(MRI)スキャンは、CT画像と比較してより優れた軟組織コントラストのために、放射線療法で使用することができる。MRIは電離放射線を含まず、組織代謝及び機能性などの人体の機能情報を捕捉するために使用することができる。
[006]
したがって、より正確な構造輪郭形成のためにCTを補完するMRIを使用することができる。しかしながら、MRIの強度値は電子密度に直接関係せず、線量計算に直接使用することはできない。したがって、MR画像を対応する導出画像、通常はCT画像(しばしば「疑似CT画像(pseudo-CT image)」と呼ばれる)に変換することが望ましい。実際のCT画像のような疑似CT画像は、線量計算のために電子密度に直接変換可能なCT値を示すデータ点の集合を有する。したがって、MR画像から得られた疑似CT画像を用いて、放射線療法治療計画における患者線量計算を容易にすることができる。従って、CT撮像から生じる追加の放射線被曝から患者を守るためには、MR画像データを用いて疑似CT画像を正確に生成することが望ましい。必要なことは、疑似CT画像が「実際の」CT画像に置き換えることができることである。
[007]
典型的には、疑似CT画像を生成するために、アトラス画像が使用される。アトラス画像は、派生画像を生成するために新しい画像をどのように変換するかを容易にするために参照として使用される既存の画像である。例えば、疑似CT画像生成コンテキストでは、アトラスMR画像及びアトラスCT画像を、新しいMR画像から誘導CT画像を生成するための基準として使用することができる。新しいMR画像の対象となる同じ患者の同じ関心領域のアトラス画像を事前に生成することができ、そこではこれらのアトラス画像を解析して関心のある構造を識別する。例えば、多くの治療又は診断状況では、患者は、治療又は診断の過程にわたって異なる時間に撮像に供される必要がある。しかしながら、これは必ずしも真実である必要はなく、例えば、アトラス画像は同じ人物の画像である必要はない。
[008]
アトラスMR画像とアトラスCT画像は、位置合わせ技術によって互いに位置合わせされる(アトラスMR画像とアトラスCT画像とが互いに「位置合わせ」されるか、アトラスMR画像とアトラスCT画像双方が「位置合わせ」するか)ことが好ましい。そのような位置合わせでは、被検体の特定の位置についてのアトラスMR画像内のポイントを、同じ特定の位置についてのアトラスCT画像内の所定のポイントにマッピングすることができる(逆もまた同様である)。しかしながら、この位置合わせにはある程度の誤差が存在する可能性がある。それ自体、アトラスMR画像とアトラスCT画像との間の位置合わせは完全ではない可能性がある。
[009]
実際のCT画像を置換するためには、疑似CT画像は、放射線療法治療計画における線量計算の目的で、又は画像のデジタル再構成放射線写真(DRR)を生成するために、患者の実際のCT画像に可能な限り近接しているべきである。しかしながら、CT画像強度値(CT値)とMR強度値との間には単純な数学的関係はない。MR強度値は標準化されておらず、またMR強度値は異なるMRスキャナ設定又は異なるMR画像シーケンスパラメータに応じて大きく変化する可能性があるので、困難が生じる。したがって、MR画像の組織セグメンテーションに基づくCT値を割り当てる、又は点比較及び重み付けされた組み合わせに基づくCT値を割り当てるなどの既存の技法は非常に粗い割り当てのみを提供するだけで、結果的に真の解剖学的詳細を欠く疑似CT画像をもたらすことなる。
[010]
したがって、治療計画の線量計算のために実際のCT画像を置き換えることが可能な改良された品質の疑似CT画像を生成し、画像誘導などのためにデジタル再構成放射線写真(DRRs:Digitally Reconstructed Radiographs)を生成する必要がある。
[011]
一態様において、本開示は、疑似CT予測モデルを生成するシステムを含む。システムは、複数の訓練対象のMRデータ及びCTデータを含む訓練データを格納するように構成されたデータベースを含むことができる。各訓練対象は、少なくともひとつのMR画像及び少なくともひとつのCT画像を有することができる。システムはまた、データベースに格納された情報にアクセスするためにデータベースに通信可能に結合されたプロセッサを含むことができる。システムは、プロセッサに通信可能に結合されたメモリをさらに含むことができる。メモリは、プロセッサによって実行されると、プロセッサが様々な操作を実行するように構成する命令を格納することができる。この操作は、データベースにアクセスして、複数の訓練対象のそれぞれについて少なくともひとつのMR画像及び少なくともひとつのCT画像を含む訓練データを取り出すことを含むことができる。この操作は、各訓練対象について、少なくともひとつのMR画像の各画像点から複数の特徴を抽出すること、抽出された特徴に基づいて各画像点の特徴ベクトルを作成すること、及び少なくともひとつのCT画像の各画像点からCT値を抽出することを含むことができる。この操作は、また、複数の訓練対象の特徴ベクトル及びCT値に基づいて疑似CT予測モデルを生成することを含むことができる。
[012]
他の一態様において、本開示は、疑似CT画像を生成するためのシステムに関する。システムは、プロセッサと、プロセッサに通信可能に結合されたメモリとを含むことができる。メモリは、プロセッサによって実行されると、プロセッサが様々な操作を実行するように構成する命令を格納することができる。この操作は、患者のMR画像を受信し、MR画像の各画像点から複数の特徴を抽出することを含むことができる。この操作は、また、抽出された特徴に基づいて各画像点の特徴ベクトルを生成することを含むことができる。この操作は、予測モデルを使用してその画像ポイントに対して作成された特徴ベクトルに基づいて各画像ポイントのCT値を決定することをさらに含むことができる。この操作は、さらに、全ての画像点に対して決定されたCT値に基づいて疑似CT画像を生成することを含むことができる。
[013]
更なる一態様において、本開示は、患者の疑似CT予測画像を生成するシステムを含む。システムは、プロセッサと、プロセッサに通信可能に結合されたメモリとを含むことができる。メモリは、プロセッサによって実行されると、プロセッサが様々な操作を実行するように構成する命令を格納することができる。この操作は、患者のMR画像を受信し、MR画像から複数の特徴を抽出することを含むことができる。この操作は、また、抽出された特徴に基づいて予測モデルを使用して中間画像を生成することを含むことができる。この操作は、中間画像からひとつ又はそれ以上の特徴を抽出することをさらに含むことができる。操作は、さらに、MR画像から抽出された複数の特徴及び中間画像から抽出されたひとつ又はそれ以上の特徴に基づいて患者の疑似CT画像を生成することを含むことができる。
[014]
更なる一態様において、本開示は、疑似CT予測モデルを生成するためのシステムを含む。システムは、複数の訓練対象のマルチチャネルMRデータ及びCTデータを含む訓練データを格納するように構成されたデータベースを含むことができる。各訓練対象は、複数のMR画像及び少なくともひとつのCT画像を有することができる。システムは、また、データベースに格納された情報にアクセスするためにデータベースに通信可能に結合されたプロセッサを含むことができる。システムは、プロセッサに通信可能に結合されたメモリをさらに含むことができる。メモリは、プロセッサによって実行されると、プロセッサが様々な操作を実行するように構成する命令を格納することができる。この操作は、データベースにアクセスして、複数のMR画像を含む訓練データ及び複数の訓練対象のそれぞれのための少なくともひとつのCT画像を取り出すことを含むことができる。この操作は、各訓練対象について、複数のMR画像に基づいて少なくともひとつの組織パラメータマップを決定し、少なくともひとつのCT画像に基づいてCT値を取得することを含むことができる。この操作は、また、組織パラメータマップ及び複数の訓練対象のCT値に基づいて疑似CT予測モデルを生成することを含むことができる。
[015]
更なる一態様において、本開示は、疑似CT画像を生成するためのシステムに関する。システムは、プロセッサと、プロセッサに通信可能に結合されたメモリとを含むことができる。メモリは、プロセッサによって実行されると、プロセッサが様々な操作を実行するように構成する命令を格納することができる。この操作は、患者の複数のマルチチャネルMR画像を受信し、複数のマルチチャネルMR画像を少なくともひとつの組織パラメータマップに変換することを含むことができる。この操作は、また、少なくともひとつの組織パラメータマップに予測モデルを適用することによってCT値を生成することを含むことができる。この操作は、予測モデルによって生成されたCT値に基づいて疑似CT画像を生成することをさらに含むことができる。
[016]
更なる一態様において、本開示は、疑似CT予測モデルを生成するためのコンピュータ実装方法を含む。この方法は、複数のマルチチャネルMR画像を含む訓練データと、複数の訓練対象のそれぞれについての少なくともひとつのCT画像とを取り出すことを含むことができる。この方法は、各訓練対象について、複数のマルチチャネルMR画像に基づいて少なくともひとつの組織パラメータマップを決定し、少なくともひとつのCT画像に基づいてCT値を取得することを含むことができる。この方法は、また、組織パラメータマップ及び複数の訓練対象のCT値に基づいて疑似CT予測モデルを生成することを含むことができる。
[017]
前述の一般的な説明及び以下の詳細な説明は、例示的及び説明的なものに過ぎず、特許請求されるような開示された実施形態を限定するものではないことを理解されたい。本開示のこれら及び他の特徴及び利点は、以下の説明、図面及び特許請求の範囲に記載される教示を検討することにより、当業者には明らかであろう。
[018]
添付の図面は、本明細書に組み込まれ、本明細書の一部を構成し、開示された実施形態を説明し、説明及び請求項とともに開示された実施形態を説明する役割を果たす。そのような実施形態は例証的であり、本装置、システム、又は方法の網羅的又は排他的な実施形態であることを意図するものではない。必ずしも一定の縮尺で描かれていない図面では、同様の参照番号は異なる図で類似の構成要素を示すことがある。文字の接尾辞又は異なる文字の接尾辞を有する同様の数字は、同様の構成要素の異なる例を表すことがある。
[019]
図1は、疑似CT予測モデルを構築するための例示的なプロセスの図である。
[020]
図2は、MR画像の各ボクセルから特徴を抽出する例示的なプロセスの図である。
[021]
図3は、図1の予測モジュールを使用して患者の疑似CT画像を生成する例示的なプロセスの図である。
[022]
図4Aは、例示的な放射線治療システムを示す。
[023]
図4Bは、例示的な放射線治療装置であるガンマナイフを示す。
[024]
図4Cは、線形加速器である例示的な放射線治療装置を示す。
[025]
図5は、疑似CT予測モデルを構築し、疑似CT画像を生成するための例示的なシステムを示す。
[026]
図6は、疑似CT予測モデルを訓練及び構築するための例示的なプロセスのフローチャートである。
[027]
図7は、疑似CT画像を生成するために疑似CT予測モデルを使用する例示的なプロセスのフローチャートである。
[028]
図8は、多段式疑似CT予測モデルを訓練及び構築するための例示的プロセスの図である。
[029]
図9は、疑似CT画像を生成するために多段式疑似CT予測モデルを使用する例示的プロセスの図である。
[030]
図10は、マルチチャネルMRスキャンから推定された組織パラメータを用いて疑似CT予測モデルを訓練するための例示的プロセスの図である。
[031]
図11は、マルチチャネルMRデータを用いて疑似CT予測モデルを構築する例示的な方法のフローチャートである。
[032]
図12は、マルチチャネルMR画像を用いて患者の疑似CT画像を生成する例示的な方法のフローチャートである。
発明の詳細な説明
[033]
ここで、開示された実施形態を詳細に参照し、その例を添付図面に示す。便宜上、図面全体を通して同じ又は同様の部分を指すために同じ参照番号が使用される。
[034]
一実施形態では、MR画像から疑似CT画像(合成CT画像又は誘導CT画像とも称される)を作成するために、訓練モジュール及び予測モジュールを含む学習ベースの手法が提供される。訓練モジュールは、選択された位置についてひとつ又はそれ以上のMR画像から抽出された特徴に基づいて任意の所与のボクセルのCT値を予測するのに使用できる予測モデル(回帰モデルとも呼ばれる)を構築する。訓練中、複数の既存の患者からMRスキャン及びCTスキャンを収集して訓練データを形成する。訓練データは、既存の患者からの事前に整列されたCT画像及びMR画像の対を含む。CT画像及びMR画像の各対について、対応するMR値及びCT値は既知であり、各ピクセル又はボクセル(2D及び3Dシナリオの両方を含む画像ポイントとも呼ばれる)について登録される。
[035]
予測モデルは、訓練データを使用して訓練することができる。訓練段階の間、収集された訓練データに回帰法(例えば、統計学習、回帰分析、機械学習など)を使用してモデルを訓練することができる。予測モデルが訓練された後、そのモデルは、予測モジュールによって使用されて、患者画像の各画像ポイントのCT値を予測することができる。したがって、訓練されたモデルを使用して、同じ又は異なる患者に対して、任意の将来のスキャンのMRデータからの疑似CT画像を作成することができる。
[036]
図1は、開示された実施形態と一致する、疑似CT予測モデル150を構築するための例示的なプロセスのフローチャートを示す。図に示すように、一実施形態は、訓練モジュールと予測モジュールとを含む学習ベースのアプローチである。訓練モジュールは、ひとつ又はそれ以上の新しいMRスキャンに基づいて疑似CT値を予測するために予測モジュールによって使用されることができる回帰モデル(例えば、疑似CTモデル150)を作成する。
[037]
一実施形態では、訓練データ110は、既存の患者又は対象(総称して「訓練対象」と呼ばれる)から収集することができる。訓練対象は、疑似CT予測モデル150を構築するために使用される、以前に撮影されたMRスキャン及び対応するCTスキャンの両方を有することができる。訓練データ110は、少なくともひとつのMRスキャン及び少なくともひとつのCTスキャン(例えば、訓練対象データ110a−110N)を有する複数の訓練対象のデータを含むことができる。訓練データを提供する訓練対象の数が多いほど(例えば、データセットが大きいほど)、
より小さいデータセットからなるモデルと比較して、より良好な疑似CT予測モデルが生成される。訓練対象データ110a−110Nは、予め整列されたMR画像及びCT画像の対を含む。MR画像とCT画像は別々に取得することができる。したがって、画像が互いに重ね合わされている場合、それらは通常一致しない。したがって、当技術分野で知られているように、MR画像及びCT画像を事前に整列させるために画像位置合わせが使用される。いくつかの実施形態によれば、訓練対象に関連するMRスキャンは、疑似CT画像が所望される新しい患者のMRスキャンのものと同じMRスキャナによって生成することができる。他の実施形態では、訓練対象に関連するMRスキャンは、異なるMRスキャナによって生成されてもよい。さらに、単一の訓練対象に関連する複数のMRスキャンは、より正確な疑似CT生成結果を提供するために、異なるコントラスト特性(例えば、T1加重、T2加重など)のMRスキャンを含むことができる。
[038]
画像特徴抽出モジュール111を使用して、訓練データ110に関連するMR画像から画像特徴を抽出することができる。画像特徴は、MRボクセルの数値(例えば、強度値、特徴の座標位置など)又はカテゴリ特性(例えば、組織タイプ、構造ラベルなど)を参照することができる。例えば、「濃淡特徴」は、MRボクセルの強度値を参照することができる。しかしながら、単一のMR特徴は、疑似CT画像を生成するためにMRボクセルを適切に表現するには不十分である可能性がある。例えば、単独で取られたMRボクセルの強度値は、CT推定のためのあいまいな表現を提供する。単一強度値は、とりわけ、同じ強度レベルの2つのMRボクセルが、異なるCT値を有する異なる組織(例えば、骨及び空気)に属し得るため、曖昧である。本明細書で使用される用語「組織」は分類を参照し、単に特定のタイプの組織を示唆するものではない。例えば、空気は組織ではない。したがって、MRボクセルのより明確な記述を提供するために、MRスキャンの各MRボクセルの複数の特徴タイプが抽出される。
[039]
複数のMR画像又はマルチチャネルMR画像では、豊富な一連の画像ベースの特徴を抽出することができ、より多くの情報を提供し、より正確な疑似CT予測に導かれる。画像特徴抽出モジュール111を使用して、各画像又は各チャネルから別々に特徴(例えば、MR特徴ベクトル120)を抽出することができる。
[040]
結果として得られるMR特徴ベクトル120は、各々が訓練対象のMRスキャン(例えば、画像ベクトルセット120a−120N)に関連する収集された画像特徴ベクトルの複数のセットを含むことができる。各画像特徴ベクトルセット120a−120Nは、複数の特徴ベクトルを含むことができる。例えば、所与の特徴ベクトルセット(例えば、120a)の各列は、ベクトルのエレメントとして複数の特徴を含む特徴ベクトルを表すことができる。特徴ベクトルの複数の特徴は、例えば、訓練対象のMR走査/画像の画像ポイント(例えば、ボクセル)に関連する異なるタイプの画像特徴を表す。特徴ベクトル内の特徴要素の数(例えば、列内のエレメントの数)は、特徴ベクトルの次元に参照される。いくつかの実施形態では、特徴ベクトルを、行又は他の適切な形態で配置することもできる。開示された実施形態による特徴抽出は、図2を参照して、以下でさらに詳細に説明される。
[041]
図2は、開示された実施形態と一致する特徴抽出プロセスを示す。患者(例えば、MR対象210)のMRスキャンの各MRボクセルについて複数の画像特徴を抽出することができる。図2に示すように、抽出された画像特徴は、局所パターン特徴212、ランドマーク特徴214、コンテキスト特徴216、及び様々な他のタイプの特徴218を含むことができる。これらの特徴の各々は、集合的に図2の特徴ベクトル列を形成する小さなブロックとして示されるひとつ又はそれ以上の特徴要素によって表されてもよい。特徴は、画像ピクセル(2D)、画像ボクセル(3D)、又は画像ポイントの集合(例えば、2D又は3Dの画像パッチ)に関連付けることができる。例えば、図2は、ボクセルi、ボクセルi+1、...、ボクセルMに関連する特徴ベクトル(例えば、列)を示す。複数の特徴ベクトル、例えば、MR画像の複数のボクセル(例えば、ボクセルi〜M)に関連する特徴ベクトルの集合は、訓練対象aに対応する特徴ベクトル集合120aを形成することができる。
[042]
潜在的な画像特徴の非限定的なリストは次のものを含む。
強度の特徴:複数のスケールでのMR画像強度値−生の強度値、又はMR強度バイアス補正及び/又はMR強度標準化・正規化のような前処理の後のいずれか。
ランドマークベースの特徴:ひとつ又はそれ以上のランドマーク点に関して所与のボクセルについて計算される相対位置、距離または他の幾何学的特徴(例えば、脳の前交連−後交連(anterior commissure-posterior commissure)(AC−PC)点、各眼球の中心など)。
コンテキストの特徴:所与の点の特定の近傍位置で計算される任意の他の画像特徴。
位置の特徴:各ボクセルの正規化された座標。正規化は、例えば、線形又は非線形の画像位置合わせのいずれかを使用して、各画像を共通の基準フレームに位置合わせすることによって達成することができる。
パッチの特徴:パッチは、いくつかの実施形態では、特徴が計算される画像ボクセルを取り囲む画像のサブ領域またはサブセットを参照することができる。例えば、パッチは5×5×5ボクセルのサイズを含み、125ボクセル位置の画像強度値はパッチの中心の点の125個の特徴要素に関連させることができる。
高レベルの特徴はひとつ又はそれ以上のパッチから導き出すことができる:これらのタイプの特徴は、SIFT(スケール不変特徴変換:Scale-invariant feature transform)、SURF(スピードアップロバスト特徴:Speeded Up Robust Features)、GLOH(勾配位置及びオリエンテーションヒストグラム:Gradient Location and Orientation Histogram)、LBP(ローカルバイナリパターン:local binary patterns)、又はHOG(指向勾配のヒストグラム:Histogram of Oriented Gradients)などのような当該技術分野で知られている様々な特徴記述子を含むことができる。そのような特徴は、考慮中のボクセルを含む各2D画像スライスについて計算されてもよく、さらに、そのような特徴は、一実施形態では3D画像に拡張されてもよい、
テクスチャの特徴:エネルギー、エントロピー、コントラスト、均質性、及び局所画像のグレースケール共出現行列の相関性、及びガボールフィルタ(Gabor filters)などで画像をフィルタリングすることによって計算されたものなど。
ジョイントの特徴:複数のMR画像(例えば、T1強調、T2強調など)が所定の訓練対象に関連付けられている場合など。そのような場合、強度、パッチ、テクスチャなどの特徴は、その後の組み合わせのために各MRスキャンから独立して抽出することができる。さらに、複数のMRスキャン間の相関を特徴付ける特徴は、各ボクセル位置、例えば局所関節ヒストグラム及び/又は局所相互相関、又は複数のMRチャネルの共分散で計算することができる。
・少なくともひとつの線形又は非線形フィルタ(例えば、局部位相、勾配、曲率、エッジ検出器、又はコーナー検出器など)を有する画像の畳み込みから得られる特徴。
・画像の変換(例えば、フーリエ変換、ヒルベルト変換、ラドン変換、距離変換、離散コサイン変換、ウェーブレット変換など)によって得られる特徴。
・領域共分散特徴:局所的な小領域内の上記の点ごとの特徴のいずれかの共分散。
・分類ベースの特徴。これについては、以下で詳しく説明する。
[043]
図2に示すように、MR画像ボクセルに関連する特徴の集合は、単一のベクトル(例えば、ボクセルi、ボクセルi+1、...、ボクセルMに関連するベクトル)で表すことができる。
[044]
図1に戻ると、画像特徴の抽出に際して、MR特徴ベクトル120は、複数の次元を有することができる(例えば、図2に示すように、特徴ベクトル内の各特徴要素は次元と見なすことができる)。しかしながら、MR画像からの抽出された特徴の数が増加すると、予測モデルを作成するタスクを達成することがより困難になる。これは、各患者の画像が通常何百万ものボクセルを含み、各ボクセルが多数の特徴と関連している可能性があるからである。したがって、複数の既存の全ての患者からの全ての画像の全てのボクセルから抽出された特徴が予測モデルを構築するために使用される場合、そのような莫大な量のデータを処理するための計算コストは非常に高価になる。その結果、実用的な次元数は、計算コストと比較して利用可能なコンピュータの処理能力に依存する。さらに、抽出された特徴を処理することから生じる予測モデルの性能は、特徴の次元の数に比例しないことがある。場合によっては、特徴の次元の数が増えるにつれて、両方の特徴が処理に含まれると、ある特徴の影響が他の特徴によりキャンセルされたり、弱められたりすることがあるので、予測モデルの性能が低下することがある。多数の特徴は、また、新しいMRデータに基づいて疑似CT画像を決定するために予測モデルを使用する際に容認できない計算コストを引き起こす可能性がある。したがって、一実施形態では、次元低減モジュール132を使用して、MR特徴によって提供される識別情報を実質的に失うことなく、低減された次元特徴ベクトル140を生成することができる。次元低減モジュール132は、元の次元数を減らすときに、元の特徴ベクトルから関連情報の大部分を取得するために使用する。例えば、MR特徴ベクトル120のいくつかの次元は、除去可能な疑似CT画像を生成することと無関係なノイズ又は他の情報を含む。他の次元には、特徴によって提供される特有の情報のよりコンパクトな表現のために組み合わせる、又は合理化できる冗長な情報が含まれている場合がある。例えば、元のデータがガウス分布に適合する場合、元のデータの平均及び標準偏差を用いて元のデータを表すことによって、元のデータの全体的な寸法を縮小することができる。このような次元削減方法により、元のデータが変換される。いくつかの実施形態では、次元数の削減のレベルは、元の特徴ベクトルを使用すること(すなわち削減なし)から任意の所定のレベルの次元を使用すること(例えば、特徴ベクトルの削減されたセットを使用すること)までに及ぶことができる。したがって、一実施形態では、次元削減モジュール132はオプションであり、元の特徴ベクトルを使用して疑似CTモデル150を生成することができる。
[045]
次元削減モジュール132が利用される場合、次元削減モジュール132によって使用される次元削減技法には、(1)教師なし次元削減と(2)教師あり次元削減の少なくとも2つのタイプの技法を含むことができる。典型的には、教師あり次元削減は、後述するように、教師なし次元削減よりも優れている。
[046]
教師なし次元削減は、重要でない雑音及びデータ冗長性を除去し、入力としてMR特徴ベクトル120のみを必要とすることができる。一般的な教師なし次元削減技法には、例えば、主成分分析(PCA)及びその非線形バージョン、カーネル主成分分析(KPCA)が含まれる。
[047]
教師なし次元削減は、他の関心あるデータを利用して、疑似CT画像を生成するのに無関係な次元をさらに除くことができる。例えば、CT値130は、次元削減のために使用されてもよい。訓練データ110のCTスキャンデータからCT値130(例えば、元のCT値又はCT数)を得ることができる。教師あり次元削減は、MR特徴ベクトル120とCT値130の両方を入力とすることができる。教師あり次元削減技法には、例えば、正準成分分析(CCA)、メトリック学習(ML)、主成分分析(SPCA)、局所感性ハッシング(LSH)、局所感性弁別分析(LSDA)などがある。関心のあるデータを離散クラスラベルに関連付けることを必要とする次元削減技法では、訓練データ110のCTスキャン又はMRスキャンに画像セグメンテーションを適用して、クラスラベルとして使用できるセグメンテーションクラスを得ることができる。
[048]
CT値130は、次元削減モジュール132によって利用されて、訓練データ110内のどの信号が基礎となるCT値に関連するかを決定することができる。元のCT値を使用すると、無関係な信号を抑制し、関連する信号を維持することができる。一般に、各訓練対象について少なくともひとつのCT画像が利用可能であるべきである。いくつかの実施形態では、複数のCT画像が利用可能であってもよい。画像ノイズを低減するために、より多くのCTスキャンを平均化することができ、それによって、次元削減モジュール132の有効性が改善される。次元削減モジュール132の出力は、縮小次元特徴ベクトル140である。
[049]
訓練データが収集され、処理されると(例えば、画像特徴抽出、次元削減技術などに供されると)、疑似CT予測モデル150は、統計学習又は機械学習技術のいずれかを使用して構築することができる。ひとつの実施形態では、回帰分析を使用して、疑似CT予測モデル150を構築することができる。回帰分析は、変数間の関係を推定するための統計的プロセスである。回帰分析を行うための多くの既知の方法がある。例えば、線形回帰又は一般最小二乗回帰は、回帰関数が訓練データから推定できる有限数の未知モデルパラメータに関して定義される点で「パラメトリック」である。疑似CT画像生成の場合、回帰モデル(例えば、式1)は、例えば、以下のように定義することができる。
ここで、「H」はCT値を示し、「X」は入力変数のベクトル(例えば、MR特徴ベクトル120又は削減された次元特徴ベクトル140のいずれか)を示し、「β」は回帰モデルに対する決定又は訓練された未知のパラメータのベクトルを示す。ある実施形態では、CT値は、CTスキャンのハウンスフィールド値( Hounsfield values)であってもよい。
[050]
MRスキャン及びCTスキャンを含む訓練データ110は、対応するX値(例えば、同じ訓練対象のMRスキャンから抽出された特徴ベクトル)を有する既知のH値のセット(例えば、訓練対象のCTスキャンに関連するCT値)を提供する。これらのデータを使用して、モデルパラメータβは、最小2乗、最大尤度などのデータフィッティング技法を使用して計算できる。βが推定されると、モデルは新しいX値の集合(例えば、新しいMRスキャンから抽出された特徴ベクトル)についてH(例えば疑似CT値)を計算することができる。
[051]
別の実施形態では、疑似CT予測モデル150を構築するために、機械学習及び教師あり学習を使用することができる。教師あり学習は、学習データのセットを与えられた予測モデルを推定する機械学習の分派である。訓練データの各サンプルは、入力データ(例えば、測定値又は特徴のベクトル)と所望の出力値(教室信号とも呼ばれる)からなるペアである。教師あり学習アルゴリズムは、訓練データを分析し、出力変数が数値又は連続であるときの回帰関数である予測関数を生成する。開示された実施形態として、多くの異なるアルゴリズムを適用することができる。アルゴリズムは、kNN(k最近傍)回帰(kNN (k-nearest neighbors) regression)、サポートベクトルマシン(support vector machines)、ニューラルネットワーク(neural networks)、決定木(decision trees)、ランダムフォレスト(random forests)、勾配ブースティングマシン(gradient boosting machines)を含むがそれらに限定されない。
[052]
図3は、開示された実施形態として、疑似CTモデル150を使用することができる予測モジュールの例示的なプロセスのフローチャートを示す。疑似CTモデル150が作成されて訓練されると、その疑似CTモデル150は、同じ患者又は新しい患者のいずれかについて、新しいMRスキャンから疑似CT画像を生成するために、適用段階の予測モジュール301によって使用される。図3に示すように、疑似CT画像350を生成するプロセスは、先に生成され訓練された疑似CT予測モデル150が適用段階で利用されることを除いて、図1について上述したプロセスと同様である。このプロセスでは、新しいMRスキャン310が予測モジュール301に入力される。一実施形態では、予測モジュール301は、画像特徴抽出モジュール311及び疑似CT予測モデル150を含むことができる。この実施形態では、MRスキャン301は対応するCTスキャンを有さない。特徴は、MRスキャン301から抽出され、MR特徴ベクトル120の生成について上述と同様の方法で、患者MR特徴ベクトル320を生成することができる。次元削減モジュール321が含まれて、患者MR特徴ベクトル320の大きさを縮小することができる。あるいは、患者MR特徴ベクトル320は、破線331により示されるように、次元の削減なしに疑似CT予測モデル150によって使用されてもよい。
[053]
このように、予測モジュール301は、訓練段階中に作成された疑似CT予測モデル150を使用して、新しいMRスキャンに対応してもともとCTスキャンが提供されなかったので、患者MR画像310の各位置で疑似CT値を予測する 。疑似CT予測モデル150は、例えば各画像位置で「点ごと」に動作することができるので、疑似CT値は、MRスキャン310内の特定の位置における特定のボクセルの特徴ベクトルに基づいて導出された値を表す。したがって、予測モデル301は、疑似CT値340を生成することができる。疑似CT値340は、疑似CT画像350の複数の強度値を表す。疑似CT画像350を生成するために、疑似CT値340は、典型的には、ボクセルのグリッド上のそれらの適切な位置に配置される。一実施形態では、予測モデル301は、画像がボクセルのグリッドであるため(例えば、すべての画像ボクセルが予測されるわけではない)、ボクセルグリッドのいくつかの値(例えば、疑似CT値340)を予測することができ、患者の解剖学的詳細の正確な視覚的表現を示すために疑似CT画像350を生成するために補間を使用することができる。
[054]
疑似CT予測モデル150は、全ての利用可能な患者の訓練データ110を使用して一度訓練され、そして、疑似CT予測モデル150は、全ての将来の新しい患者に使用され得る。あるいは、同じ疑似CT予測モデル150をすべての患者に使用することはできない。疑似CT予測モデル150は、特定の患者のためにカスタマイズすることができる。例えば、訓練データは、新しい患者に最も類似又は関連するデータを含むように選択されてもよく、モデルは、新しい患者に対して特異的に構築されてもよい。
[055]
図4Aは、本開示のいくつかの実施形態による、例示的な放射線治療システム400を示す。放射線療法システム400は、訓練モジュール412、予測モジュール414、訓練データベース422、検査データベース424、放射線治療装置430、及び画像取得装置440を含むことができる。放射線療法システム400は、患者情報を提供することができる治療計画システム(TPS)442及び腫瘍学情報システム(OIS)444に接続することもできる。さらに、放射線治療システム400は、表示装置及びユーザインターフェース(図示せず)を含むことができる。
[056]
図4Bは、本開示のいくつかの実施形態による、放射線治療装置430ひとつのタイプの例(例えば、スウエーデン、ストックホルムのエレクタ アクチボラゲット(パブル)(Elekta、AB、Stockholm、Sweden)によって製造されたLeksell Gamma Knife)を示す。図4Bに示すように、放射線治療セッションでは、患者452は、手術又は放射線療法を受けている患者の体部分(例えば、頭部)を安定に保つために座標フレーム454を着用してもよい。座標フレーム454及び患者位置決めシステム456は、患者を撮像している間、又は放射線手術中に使用することができる空間座標系を確立することができる。放射線治療装置430は、複数の放射線源462を囲む保護ハウジング464を含むことができる。放射源462は、ビームチャネル466を介して複数の放射ビーム(例えば、ビームレット)を生成することができる。複数の放射ビームは、異なる方向からアイソセンタ458に集束するように構成することができる。個々の放射ビームの強度は比較的低いが、アイソセンタ458は、異なる放射ビームからの複数の線量がアイソセンタ458に蓄積するとき、相対的に高いレベルの放射を受け取ることができる。特定の実施形態では、アイソセンタ458は、腫瘍などの外科手術又は治療下の標的に対応し得る。放射線治療装置430(例えば、スウエーデン、ストックホルムのエレクタ アクチボラゲット(パブル)(Elekta、AB、Stockholm、Sweden)によって製造されたLeksell Gamma Knife)は、一実施形態では、割り当てられたバルク密度を有するMR画像、又はMR画像と融合されたCT画像を利用し、開示された実施形態として生成される。
[057]
図4Cは、本開示のいくつかの実施形態による、放射線治療装置430(例えば、線形加速器470)の別の例を示す。線形加速器470を使用して、患者472を患者テーブル473上に配置して、治療計画によって決定された放射線量を受けることができる。線形加速器470は、放射ビーム476を生成する放射ヘッド475を含むことができる。放射ヘッド475全体は、水平軸477の周りを回転可能であってもよい。さらに、患者テーブル473の下には、アイソセンタ471の周りで放射線ヘッド475と同期して回転するフラットパネルシンチレータ検出器474が設けられてもよい。放射ヘッド475によって生成された軸477とビーム476の中心との交差点は、通常アイソセンタと呼ばれる。患者テーブル473は、患者472がアイソセンタ471又はその近くに腫瘍部位とともに配置されるように電動化されてもよい。放射線ヘッド475は、ガントリ478の周りを回転して、治療計画に従って複数の変化する放射線量を患者472に提供することができる。別の実施形態では、線形加速器470は、MR線形加速器(「MR−LINAC」)であってもよい。線形加速器10及びMR−LINACは、一実施形態では、MR画像、CT画像を利用することができ、開示された実施形態により生成された疑似CT画像を使用することができる。
[058]
図5は、開示された実施形態と一致して、疑似CT予測モデルを構築し、疑似CT画像を生成するための例示的なシステム500である。いくつかの実施形態によれば、システム500は、開示された実施形態と一致する大量のデータを識別、分析、維持、生成、及び/又は提供することができるひとつ又はそれ以上の高性能コンピューティングデバイスであってもよい。システム500は、スタンドアローンでもよいし、サブシステムの一部であってもよく、サブシステムは、より大きなシステムの一部であってもよい。例えば、システム500は、遠隔地に位置し、インターネットなどのネットワーク又はLAN又はWANなどの専用ネットワークを介して通信する分散型高性能サーバを表すことができる。いくつかの実施形態では、システム500は、ひとつ又はそれ以上の遠隔配置された高性能コンピューティングデバイスと通信する組み込みシステム、MRスキャナ、及び/又はタッチスクリーンディスプレイデバイスを含むことができる。
[059]
一実施形態では、システム500は、ひとつ又はそれ以上のプロセッサ514、ひとつ又はそれ以上のメモリ510、及びひとつ又はそれ以上の通信インターフェース515を含むことができる。プロセッサ514は、処理装置であってもよく、マイクロプロセッサ、中央処理装置(CPU)、グラフィック処理装置(GPU)などのようなひとつ以上の汎用処理装置を含む。より詳細には、プロセッサ514は、CISCマイクロプロセッサ、RISCマイクロプロセッサ、VLIWマイクロプロセッサ、他の命令セットを実装するプロセッサ、又は命令セットの組み合わせを実装するプロセッサであってもよい。プロセッサ514は、特定用途向け集積回路(ASIC)、フィールドプログラマブルゲートアレイ(FPGA)、デジタル信号プロセッサ(DSP)、システムオンチップ(SoC)などのひとつ以上の専用処理装置であってもよい。当業者に理解されるように、いくつかの実施形態では、プロセッサ514は、汎用プロセッサではなく、専用プロセッサであってもよい。プロセッサ514は、Intel(登録商標)によって製造された、Pentium(登録商標)又はXeon(登録商標)ファミリーや、AMD(登録商標)によって製造されたTurion(登録商標)ファミリーからのマイクロプロセッサ、Sun Microsystemsによって製造された様々なプロセッサのような、ひとつ又はそれ以上の既知の処理デバイスを含むことができる。プロセッサ514はまたNvidia(登録商標)により製造されたグラフィカル処理ユニットを含むことができる。開示された実施形態は、大量のイメージングデータ、又は開示された実施形態と一致する任意の他のタイプのデータの識別、分析、維持、生成、及び/又は提供の計算命令を満たすように構成された任意のタイプのプロセッサに限定されない。
[060]
メモリ510は、開示された実施形態に関連する機能を実行するためにプロセッサ514によって使用されるコンピュータ実行可能命令を格納するように構成されたひとつ又はそれ以上の記憶装置を含むことができる。
例えば、メモリ510は、治療計画ソフトウェア511、オペレーティングシステムソフトウェア512、及び訓練/予測ソフトウェア513のためのコンピュータ実行可能ソフトウェア命令を格納することができる。プロセッサ514は、メモリ/記憶装置510に通信可能に結合されてもよく、プロセッサ514は、開示された実施形態によるひとつ又はそれ以上の動作を実行するために記憶されたコンピュータ実行可能命令を実行するように構成されてもよい。例えば、プロセッサ514は、訓練/予測ソフトウェア513を実行して、訓練モジュール412及び予測モジュール414の機能を実装することができる。さらに、プロセッサ装置514は、訓練/予測ソフトウェア513にインターフェースされる治療計画ソフトウェア511(例えば、Elektaによって製造されたMonaco(登録商標)ソフトウェアなど)を実行することができる。
[061]
開示された実施形態は、専用のタスクを実行するように構成された別々のプログラム又はコンピュータに限定されない。例えば、メモリ510は、システム500又は複数のプログラム(例えば、治療計画ソフトウェア511及び/又は訓練/予測ソフトウェア513)の機能を実行する単一のプログラムを含むことができる。さらに、プロセッサ514は、データベース520に格納されたプログラムなど、システム500から遠隔に位置するひとつ又はそれ以上のプログラムを実行してもよく、このような遠隔プログラムは、腫瘍情報システムソフトウェア又は治療計画ソフトウェアを含むことができる。メモリ510はまた開示された実施形態による動作を実行するためにシステムが使用することができる任意のフォーマットの画像データ又は任意の他のタイプのデータ/情報を格納することができる。
[062]
通信インターフェース515は、データがシステム500によって受信及び/又は送信されるように構成されたひとつ又はそれ以上のデバイスであってもよい。通信インターフェース515は、システム500がシステム500の遠隔に位置するコンポーネント、データベース520、又は病院データベース530などの他の機械及びデバイスと通信することを可能にするひとつ又はそれ以上のデジタル及び/又はアナログ通信デバイスを含むことができる。例えば、プロセッサ514は、通信インターフェース515を介してデータベース520又は病院データベース530に通信可能に接続されてもよい。例えば、通信インターフェース515は、インターネットなどのコンピュータネットワーク、又はLAN又はWANなどの専用ネットワークとすることができる。あるいは、通信インターフェース515は、プロセッサ514がデータベース520、530のいずれかとの間でデータを送受信することを可能にする衛星通信リンク又は任意の形式のデジタル又はアナログ通信リンクであってもよい。
[063]
データベース520及び病院データベース530は、情報を記憶し、システム500を介してアクセス及び管理されるひとつ又はそれ以上のメモリデバイスを含むことができる。例えば、データベース520、病院データベース530、又はその両方は、Oracle(登録商標)データベース、Sybase(登録商標)データベースなどのリレーショナルデータベース、又は、Hadoopシーケンスファイル、HBase、Cassandraなどの非リレーショナルデータベースを含むことができる。データベース又は他のファイルは、例えば、訓練対象に関連するMRスキャン又はCTスキャンからの生データ、MR特徴ベクトル120、CT値130、縮小次元特徴ベクトル140、疑似CT予測モデル150、疑似CT値340、疑似CT画像350、DICOMデータ等を含む。しかしながら、開示された実施形態のシステム及び方法は、別々のデータベースに限定されない。一態様では、システム500は、データベース520又は病院データベース530を含むことができる。あるいは、データベース520及び/又は病院データベース530は、システム500から遠隔に位置してもよい。データベース520及び病院データベース530は、データベース520又は病院データベースのメモリデバイスに格納されたデータの要求を受信及び処理し、データベース530又は病院データベース530からデータを提供するように構成された計算コンポーネント(例えば、データベース管理システム、データベースサーバなど)を含むことができる。
[064]
システム500は、ネットワーク(図示せず)を介してシステム500の他のデバイス及びコンポーネントと通信することができる。ネットワークは、通信を提供し、情報を交換し、又は情報の交換を容易にし、ネットワーク(図示せず)を介して他のデバイス及び/又はシステム500のコンポーネント間で情報の送受信を可能にする任意のタイプのネットワーク(インフラストラクチャを含む)であることができる。他の実施形態では、システム500のひとつ又はそれ以上のコンポーネントは、システム500とデータベース520と病院データベース530との間のリンク(例えば、有線リンク、無線リンク、又は衛星リンク、又は他の通信リンク)のような専用通信リンクを介して直接的に通信することができる。
[065]
システム500の機能構築ブロックの構成及び境界は、ここでは説明の便宜のために定義されている。指定された機能及びそれらの関係が適切に実行される限り、代替の境界を定義することができる。本明細書に含まれる教示に基づいて、当業者には代替物(本明細書に記載されているものの均等物、拡張物、変形物、偏差などを含む)は明らかであろう。そのような代替案は、開示された実施形態の範囲及び精神の範囲内に入る。
[066]
図6は、開示された実施形態による、疑似CT予測モデルを訓練及び構築するための例示的なプロセス600のフローチャートである。プロセス600は、複数のステップを含み、そのいくつかは任意であってもよい。ステップ610において、システム500は、例えば、データベース520及び/又は病院データベース530から、複数の訓練対象に関連する訓練データ110にアクセスする。訓練データ110は、各訓練対象(例えば、図1に示すような、訓練対象データ110a−110N)について、少なくともひとつのMRスキャン及び少なくともひとつのCTスキャンを含むことができる。いくつかの実施形態では、訓練データ110は、同じ患者に対する少なくともひとつのMRスキャン及び複数のCTスキャンを含むことができる。
[067]
いくつかの実施形態によれば、システム500は、疑似CT予測モデル150を訓練及び構築するために使用される前に、訓練データ110のいくつか又はすべてが前処理を必要とするか否かを決定することができる。ステップ620において、プロセッサ514は、訓練データ内のひとつ又はそれ以上の訓練対象について、MRスキャン及び対応するCTスキャンが整列されているかどうかを判定する(例えば、各MRボクセルについて、対応するCTボクセルからのCT値は、知られている)。MR画像及びCT画像が整列していない場合、プロセス600は、ステップ624で走査を整列させるために分岐621(例えば「NO」)に進む。システム600は、必要に応じて、当業者に公知の方法に従ってMRスキャン及び対応するCTスキャンを整列させることができる。あるいは、MR画像及びCT画像が整列されている場合、プロセス600は分岐622(例えば「YES」)に進み、ステップ630に続く。
[068]
任意に、ステップ630において、プロセッサ514は、訓練データ110が同じ訓練対象について複数のCTスキャンを含むかどうかを検証する。複数のCTスキャンがある場合、プロセッサ514は、同じ患者の画像ノイズを低減するために、複数のCTスキャン間の対応するCTボクセルの平均CT値を決定する。そうでない場合、プロセス600は、ステップ620からステップ640に直接進む。
[069]
ステップ640において、前処理の一部として、プロセッサ514は、例えば、治療計画ソフトウェア511又は訓練/予測ソフトウェア513に反映されたシステム設定に基づいて、MRスキャンから画像アーチファクトを低減するか又は消滅させるかを決定する。画像アーチファクトの低減が所望される場合、プロセス600は、分岐642(「YES」)からステップ644に進む。ステップ644において、プロセッサ514は、前処理の一部として、画像アーチファクト低減技術を適用する。MR走査を前処理することにより、プロセッサ514は、強度不均一性(MR画像バイアスフィールドとしても知られている)及び画像ノイズなどの画像アーチファクトを除去又は低減することができる。さらに、前処理は、異なるMRスキャナタイプ(例えば、GE、Siemensなどによって製造されたもの、又は0.5テスラ、1.5テスラなどの様々な磁場強度)によってMR画像強度値を正規化/標準化することができる。前処理技術は、新しい特許MRスキャン310(図3に示す)の画像アーチファクトを除去又は低減するために使用することもできる。画像アーチファクト低減が行われない場合、プロセス600は、特徴を抽出するためにステップ650に進む。いくつかの実施形態では、前処理のひとつ又はそれ以上のステップ(例えば、図6の破線で囲まれた)を省略することができる。
[070]
ステップ650では、訓練データ110から特徴を抽出する。いくつかの実施形態では、システム500は、訓練データ110内の各MRスキャンの各ボクセルから特徴を抽出することができる。例えば、MR画像自体を使用することができ、MR画像からの各ボクセル又は選択されたボクセルを使用して特徴を抽出することができる。あるいは、プロセッサ514は、MR画像を異なる組織タイプにセグメント化し、組織タイプに基づいて各MRスキャンの画像ボクセルをセグメント化することができる。これは、例えば、他の抽出された特徴に加えて、追加の特徴として組織タイプを使用することができるので、場合によっては有利であり得る。
[071]
ステップ660において、システム500は、MRスキャンの各ボクセルについて抽出された画像特徴に基づいてMR特徴ベクトルを作成する。したがって、MRスキャンの各ボクセルについて複数の特徴を含むベクトルをプロセッサ514によって生成することができる。複数のMR特徴ベクトル120は、MR走査の複数のボクセルに対してプロセッサ514によって生成され得る。
[072]
ステップ670において、システム500は、訓練データ内の各CTスキャンの各ボクセルからCT値を抽出する。
[073]
ステップ680において、システム500は、MR特徴ベクトル120に関連する次元の数を減らすかどうかを決定する。
[074]
例えば、システム500のプロセッサ514は、疑似CT予測モデル150によって処理されたときに、MR特徴ベクトル120に関連する次元の数が高い計算コストをもたらすか、又は潜在的に性能問題を引き起こすと判断することができる。別の例では、システム500は、MR特徴ベクトル120が、疑似CT予測モデル150の精度に影響を及ぼすと考えられる閾値を超えるノイズ又は複製データを含むと判断することができる。別の実施形態では、システム500は、性能及び/又は出力品質に影響を及ぼす要因に基づいて次元削減を行うかどうかを決定することができる。従って、プロセッサ514が、次元削減が必要であると判定した場合、プロセス600は、分岐682(例えば「YES」)からステップ686に進み、プロセッサ514は、MR特徴ベクトル120に関連する次元を縮小することができる。代替的に、いくつかの実施形態では、システム500は、MR特徴ベクトル120の次元削減を行わないように(例えば、ユーザからの)入力を受け取ることができる。
[075]
次元削減が必要ない場合、プロセス600は分岐684から直接ステップ690に進むことができる。ステップ690において、システム500は、MR特徴ベクトル(例えば、120)及び抽出されたCT値に基づいて疑似CT予測モデル150を生成するために、統計又は機械学習技術を利用することができる。いくつかの実施形態では、次元削減された特徴ベクトル140を利用することができる。
[076]
いくつかの実施形態によれば、訓練データ110のサブセットは、疑似CT予測モデル150を訓練して構築するための基礎として使用することができる。したがって、システム500は、疑似CT予測モデル150を訓練及び構築するために、訓練データ110のサブセットを(例えば、治療計画ソフトウェア511及び/又は訓練/予測ソフトウェア513に反映されたユーザ入力及び/又はシステム設定に基づいて)決定することができる。別の実施形態では、訓練データ110のサブセットは、特定の画像領域に基づいて分類することができる。例えば、訓練データ110のサブセットは、1)特定の解剖学的領域に関して、2)様々な組織分類に関して、又は3)訓練対象の特徴に関して、であることができる。
[077]
例えば、ひとつ又はそれ以上の特徴は、所与の患者のMRスキャンのための基礎となる解剖学的構造の優れた解釈を提供することができる。したがって、疑似CT予測モデル150を訓練する際に使用するために、訓練データ110から上位の特徴のサブセットのみを抽出することができる。優れた特徴を使用して、所与の患者のMRスキャンに対する対応する疑似CT値340を推定する疑似CT予測モデル150の予測力を改善することができる。特徴のサブセットは、ひとつ又はそれ以上の疑似CTモデルを生成及び訓練するために使用することができる。
[078]
一実施形態では、特定の解剖学的領域に関して疑似CT予測モデルを構築する場合、訓練データ110のサブセット(例えば、関心のある身体領域に関連する訓練データ110のみ)を使用して、疑似−CT予測モデル150を訓練し、構築することができる。ひとつの疑似CT予測モデルの代わりに、プロセッサ514は、身体の特定の解剖学的領域(例えば、頭部、上半身、下半身など)に関して複数の疑似CT予測モデルを生成することができる。したがって、プロセッサ514は、MR特徴ベクトル120(又は削減次元特徴ベクトル140)及び所定の解剖学的位置のCT値130を利用して、MRスキャン内に表れた所定の解剖学的関心位置の疑似CT予測モデル150を生成することができる。
[079]
例えば、システム500は、患者MRスキャン310が、患者の前立腺のMR画像を含むと判断することができる。したがって、開示された実施形態により、システム500は、訓練データとして前立腺のひとつ又はそれ以上のMRスキャン及びCTスキャンを利用して訓練データ110に基づいて構築及び訓練された疑似CT予測モデルを特定することができる。一実施形態では、複数の疑似CT予測モデルが利用可能であり、各モデルは、例えば、前立腺の様々な解剖学的側面を表すことができる。したがって、複数の疑似CT予測モデルを生成することができ、各疑似CT予測モデルは、特定の解剖学的領域(例えば、前立腺の疑似CT予測モデル、右肺の疑似CT予測モデル、左肺の疑似CT予測モデル、脳の疑似CT予測モデルなど)に対するものである。
[080]
別の実施形態では、疑似CT予測モデル150は、組織分類のような分類ベースの特徴に基づいて生成されてもよい。例えば、システム500は、画像特徴抽出モジュール111を使用して、組織クラス(例えば、骨、脂肪、筋肉、水分、空気、及び、心臓組織、肺組織、肝臓組織、脳組織などのような構造クラス)に従って、訓練データ110内の各MRスキャンの画像ボクセルをセグメント化することができる。セグメント化された画像ボクセルに基づいて、各MRスキャンに対する複数のセグメント化マップを生成することができる。画像特徴は、セグメンテーションマップから抽出することができる。セグメンテーションマップ抽出画像特徴は、各ボクセルのMRスキャン抽出画像特徴と組み合わせることができる。MR特徴ベクトルは、結合された画像特徴に基づいて各訓練対象について決定されてもよい。合成されたMR特徴ベクトルと抽出されたCT値に基づく疑似CT予測モデルを生成することができる。上述のように、用語「組織」は、単に特定のタイプの組織(例えば、空気は組織ではない)を示唆するためではなく、分類として使用されている。
[081]
さらに別の実施形態では、疑似CT画像を生成するプロセスは、訓練対象の特性に従って選択された訓練データを使用することに基づくことができる。例えば、システム500は、訓練対象のサブセットのうちのひとつ又はそれ以上の共通特性を識別することができる。例えば、システム500は、各訓練対象に関連する年齢、性別、体重クラス等を特定し、ひとつ以上の共通特性を有する訓練対象を選択してもよい。他の例では、システム500は、訓練データ110内のMR及びCTスキャンに基づいて訓練対象のひとつ又はそれ以上の特性を識別することができる。さらに、システム500は、訓練対象のサブセットと共通の患者(例えば、新しい患者)のひとつ又はそれ以上の特性を識別することができる。例えば、システム500は、患者のひとつ又はそれ以上の特性を識別し、患者の特性を、共通特性を識別するために訓練対象について識別された特性と比較することができる。次いで、システム500は、疑似CT予測モデル150を訓練し構築するために、訓練データとしてひとつ又はそれ以上の共通の特性を有する訓練対象を選択することができる。
[082]
画像特徴は、訓練対象の特徴に関連するCTスキャンからのMRスキャン及びCT番号から抽出することができる。例えば、システム500は、新しい患者と共通の特性を有する訓練対象のサブセットに関連する訓練データ110のCTスキャンからMRスキャン及びCT値130から画像特徴を抽出することができる。次に、抽出された画像特徴に基づいて、サブセットの各訓練対象についてMR特徴ベクトルを決定することができる。これらのMR特徴ベクトルと抽出されたCT値とに基づいて疑似CT予測モデルを生成することができる。
[083]
疑似CT予測モデル150は、すべての訓練データ110を使用して訓練され、次いで、新しい患者のための新しいMRスキャンのために利用され得る。疑似CT予測モデル150は、すべての将来の新しい患者にも使用することができる。いくつかの実施形態では、同じ疑似CT予測モデル150をすべての患者に使用することはできない。疑似CT予測モデル150は、特定の患者についてカスタム生成されてもよい。例えば、訓練データは、新しい患者に類似又は関連する訓練対象に基づいて選択されてもよく、モデルは新しい患者のために特別に構築されてもよい。
[084]
医療従事者は、最適な治療又は診断を決定するために、患者の関心領域のMR特性及びCT特性の両方を評価することが有用であることが分かる。さらに、疑似CTモデルを使用して、MR画像からCT画像を導出して、放射線療法治療計画における患者線量計算を容易にすることができる。これは、CT画像から生じる追加の放射線被曝から患者を守るために、MR画像から疑似CT画像を正確に生成するために望ましい。実際のCT画像を置換するためには、疑似CT画像は、放射線療法治療計画における線量計算の目的で、又は画像ガイダンスに対するデジタル再構成放射線写真(DRR)を生成するために、患者の実際のCT画像にできるだけ近づけるべきである。しかし、CT画像強度値(CT値)とMR強度値との間には単純な数学的関係はない。MR強度値は標準化されておらず、異なるMRスキャナ設定又は異なるMR画像シーケンスパラメータに応じて大きく変化し得るので、困難が生じる。
[085]
図7は、開示された実施形態のように、モデルが(図3に記載されているように)疑似CT値および疑似CT画像を生成するように訓練された後、(図1及び図6により説明されるように)疑似CT予測モデルを使用するための例示的なプロセス700のフローチャートである。ステップ710において、システム500は、患者(例えば、新しい患者)に関連する少なくともひとつのMRスキャン(例えば、MRスキャン310)を受信することができる。少なくともひとつのMRスキャンは、対応するCTスキャンを有していなくてもよい。MRスキャンは、疑似CT画像を生成するために使用される。
[086]
ステップ720において、プロセッサ514は、MR画像ボクセルがセグメント化されるべきかどうかを決定することができる。MRスキャンをセグメント化することは任意である。プロセッサ514は、治療計画ソフトウェア511(図5に示す)又はユーザインターフェース(図示せず)から命令を受信して、MRスキャンを分割すべきかどうかを示すことができる。そうである場合、プロセス700は、分岐722(例えば「YES」)に進み、MRスキャンを分割する。ステップ730では、MRスキャンの画像ボクセルが、例えば、組織分類に従ってセグメント化される。セグメント化は、当業者に知られているセグメント化技術に従って実施することができる。例えば、プロセッサ514は、ひとつ又はそれ以上のセグメンテーションマップを作成するために、k−meansクラスタリングセグメンテーション法(k-means clustering segmentation method)、ファジーc−meansセグメンテーション法(fuzzy C-means segmentation method)などを採用することができる。
[087]
プロセッサ514はさらに進んだセグメンテーション方法を使用することができる。例えば、プロセッサ514は、セグメンテーションを実行するための学習ベースまたは特徴ベースのアプローチを採用することができ、それは、例えば、アルゴリズム(例えば、局所パターン特徴、ランドマーク特徴、コンテキスト特徴など)を使用して分類予測モデルを構築することを含み、画像ボクセルの特徴に基づいて各画像ボクセルの組織ラベルを予測することができる。
[088]
ステップ732において、プロセッサ514は、分類された画像ボクセルに基づいて各MRスキャンについて複数のセグメント化マップを生成して、分類予測モデルを作成することができる。例えば、バイナリボーンセグメンテーションマップ(binary bone segmentation map)は、ボーンとしてラベル付けされたボクセルでは「1」に等しい値を有し、他のすべてのボクセルでは「0」である画像であってもよい。プロセッサ514は、セグメンテーションマップを使用して元のMR画像から追加の特徴を抽出することができる。開示された実施形態による、上記に開示された学習ベースの方法を用いて、ひとつ又はそれ以上のセグメンテーションマップを生成するための予測モデルを訓練し構築することができる。あるいは、MRスキャンのセグメンテーションが必要でない場合、プロセス700は、分岐724(例えば「NO」)に進む。
[089]
ステップ740において、プロセッサ514は、患者のMRスキャンから画像特徴を抽出することができる。上記で説明したMRスキャンをセグメント化する任意の経路が実行された場合、抽出された画像特徴を経路722に沿って提供することができる。ステップ734において、プロセッサ514は、MR画像からセグメンテーションマップを使用して抽出された追加の特徴と、MR画像から直接抽出された特徴とを組み合わせて、各データポイント(例えば、ボクセル)の特徴の組合せセットを形成する。
[090]
MRスキャンが分割されているかどうかにかかわらず、画像特徴が抽出された後、プロセス700は(例えば、経路744又は736に沿って)ステップ750に進む。ステップ750において、プロセッサ514は、抽出された画像特徴から各訓練対象についての特徴ベクトル120を決定することができる。
[091]
ステップ760において、プロセッサ514は、MR特徴ベクトル120を疑似CT予測モデル150に入力することができる。ステップ770において、プロセッサ514は、疑似CT予測モデル150を入力MR特徴ベクトル120に適用して、患者MR画像310の各ボクセルのCT数(例えば、疑似CT値340)を決定することができる。
[092]
ステップ780において、プロセッサ514は、疑似CT値340に基づいて、患者の疑似CT画像350を生成することができる。得られた疑似CT画像350は、治療計画における線量計算、画像誘導のためのDRRの生成などのために使用されてもよい。
[093]
図8は、開示された実施形態と一致する、疑似CT予測モデルを構築するための訓練データを増強するための例示的な方法の図を示す。図8に示す方法は、カスケード訓練技法又は多段階訓練技法とも呼ばれ、最初に訓練された予測モデルが中間予測結果を生成するために使用され、中間予測結果は、データを訓練して、述語モデルをさらに洗練させる。カスケード訓練技術は、複数の訓練段階を含むことができる。各段階において、予測モデルは、前の段階で生成された予測モデルによって生成された予測結果と組み合わされた初期データ(例えば、訓練データ110)を用いて訓練される。
[094]
一実施形態では、図1に関して上述したように、疑似CT予測モデル150は、当初構築され、訓練データの初期セットを使用して訓練されてもよい。最初の訓練プロセスの一部として、複数の画像走査から画像特徴が抽出され、抽出された画像特徴から特徴ベクトルが決定される。さらに、CTスキャンから対応するCT値を決定することができる。例えば、各訓練主題のための少なくともひとつのCT画像が利用可能であるべきである。一実施形態では、複数のCT画像が利用可能であってもよい。複数のCT画像が利用可能であれば、CT画像を平均して画像ノイズを低減することができる。疑似CTモデルは、カスケード訓練プロセスの任意の段階で訓練され得る。
[095]
別の実施形態では、分類予測モデルは、初期モデル又は任意の中間モデルとして訓練されてもよい。上述のように、分類予測モデルは、分類ベースの特徴を抽出して次の訓練段階で使用することができるひとつ以上の区分マップを予測するために使用されてもよい。
[096]
例示的なカスケード訓練プロセスでは、疑似CT予測モデルが最終段階で訓練され構築されている限り、疑似CT予測モデル及び分類予測モデルを複数のステージ間の任意の組み合わせで使用することができる。
[097]
図8に示すように、最初に構築された予測モデルはモデル#1として示されている。モデル#1は、訓練段階#1の間、各MRスキャン、各CTスキャン、又はMR及びCTスキャンの各ペア(例えば、110a、110b)のような元の訓練データを処理することによって生成される。上述したように、モデル#1は、疑似CT予測モデル又は分類予測モデルであってもよい。モデル#1は、次の段階の訓練プロセスで使用できる。例えば、モデル#1は、複数の予測結果(例えば、予測結果1、ステージ1;予測結果2、ステージ1;...;予測結果N、ステージ1)を生成するために使用することができる。例えば、モデル#1が疑似CT予測モデルである場合、訓練対象aのMRスキャンをモデル#1への入力として使用して(例えば、訓練対象aが新しい患者であり、訓練対象aのMRスキャンが新しいMRスキャンであるかのように)、疑似CT予測結果1120a:予測結果1、ステージ1を生成することができる。別の例では、モデル#1が分類予測モデルである場合、モデル#1の予測結果はセグメントマップであってもよい。他の予測結果1120b...1120Nも同様に生成することができる。次いで、これらの予測結果は、初期訓練データに増強されて、拡張訓練データ1120を形成することができる。例えば、モデル#1から生成された対応する予測結果(例えば、1120a、1120b、...1120N)とMR及びCTスキャンのペア(例えば、110a、110b、...110N)を関連付けることによって、拡張された訓練データ1120を作成することができる。拡張訓練データ1120は、別のモデル(例えば、モデル#2)を生成するために、次の訓練段階で訓練モジュール(例えば、訓練モジュール412)によって使用され得る。
[098]
このようにして、各段階で、新しく洗練された訓練されたモデル(例えば、モデル#2、モデル#3、...モデル#M)を生成することができる。前段階で開発されたモデルは、拡張訓練データ(例えば、拡張訓練データ1120,1130など)を訓練モジュールに適用することによって改良することができる。例えば、モデル#2は、モデル#1によって生成された予測結果を含む拡張訓練データ1120を使用して生成することができる。モデル#2を生成するために、拡張された訓練データ1120は、訓練段階#2において訓練モジュール(例えば、訓練モジュール412)に入力され得る。いくつかの実施形態では、モデル#1によって生成される予測結果は、各訓練対象についての画像(例えば、疑似CT画像)又はマップ(例えば、セグメンテーションマップ)であってもよい。次いで、画像特徴は、例えば、画像特徴抽出モジュール111を使用して、予測結果(例えば、画像又はマップ)から抽出することができる。予測結果から抽出することができる画像特徴は、輝度特徴、文脈特徴、パッチ特徴、局所パターン特徴、ランドマーク特徴など、従来の文章で論じた任意の特徴を含むことができる。予測結果から抽出された特徴は、元のMR画像から抽出された特徴と組み合わされて、各画像点のための新たな拡張特徴ベクトルを形成することができる。拡張特徴ベクトルは、例えばモデル#2を訓練するために次の訓練段階で使用されてもよい。各後続の予測モデル(例えば、モデル#2)は、その以前の予測モデル(例えば、モデル#1)の予測結果を使用して構築され、訓練されるので、予測結果から明らかにされた新しい情報を訓練プロセスに加え、後続の予測モデルを改善することができる。各逐次モデル(例えば、モデル#1、モデル#2、...など)を訓練するために拡張訓練データ(例えば、拡張されたトレーニングデータ1120、1130など)を使用するこのプロセスは、最終予測モデルモデル#Mが訓練され構築されるまで続く。カスケード訓練プロセスの目標が疑似CT予測モデルを構築することである場合、最後のモデルであるモデル#Mは疑似CT予測モデルであり、中間ステージの他のモデルはいずれのモデルでもよい。利用される段階の数は、CT値を正確に予測するためのモデル#Mの検証に依存し得る。例えば、最新のモデルによって生成された疑似CT予測値と元のCT値との間の差が所定の閾値未満である場合、この反復プロセスは停止することができる。別の例では、反復プロセスは、連続する疑似CT予測モデルの予測結果の間の差が所定の閾値未満であるときに停止することができる。
[099]
上述したように、カスケード訓練技術は、疑似CT予測モデル及び/又は組織分類モデルの訓練及び構築に適用可能である。組織分類モデルを訓練及び構築する場合、各中間予測モデル(例えば、モデル#1、モデル#2、...モデル#M−1)は、例えば組織ラベルを反映するセグメントマップを提供することができる組織分類モデルであり得る疑似CT値を生成する。上述の多段階訓練プロセスを使用して、各組織分類モデルは、以前のモデルの予測結果を含む拡張データを使用して構築され、訓練され、各段階でモデルを絶えず改良することができる。さらに、組織分類及び疑似CT予測は、多段階プロセスに混在させることができる。例えば、最初のK段階では、訓練されたモデルは組織分類モデルであってもよく、予測結果は組織分類結果であってもよい。組織分類モデルが訓練されると、K+1段階で組織分類結果を使用して疑似CT値を予測することができ、ここで組織分類結果を使用して、抽出された他の特徴と共に特徴の集合をMRスキャンから抽出することができる。M段階に達するまで(例えば、M>K+1の場合)、単一の余分な段階を実行することができる(例えば、M=K+1である場合)又は追加の段階を実行することができる。プロセスの最後に、最終的な予測モデルであるモデル#Mが疑似CT値を生成するように訓練され構築される。
[0100]
図9は、開示された実施形態により、疑似CT値を予測するために多段階モデルを適用する例示的プロセスの図を示す。
[0101]
図9に示すように、プロセッサ(例えば、図5に示すプロセッサ514)は、画像取得装置440又はデータベース520,530(図5に示す)からひとつ又はそれ以上の患者のMRスキャン310を取得することができる。
[0102]
患者のMRスキャンが取得されると、複数の画像特徴をMRスキャンの各MRボクセルについて抽出することができる。図2に関して上述したように、抽出された画像特徴は、局所パターン特徴212、ランドマーク特徴214、コンテキスト特徴216、及び様々な他のタイプの特徴218を含むことができる。特徴は、画像ポイント、画像ボクセル、又は画像サンプルに関連付けることができる。図9に示すように、画像特徴抽出モジュール311を使用して、MR走査から画像特徴を抽出することができる。画像特徴は、MRボクセルの数値(例えば、強度値、特徴の座標位置など)又はカテゴリ特性(例えば、組織タイプ、構造ラベルなど)を参照することができる。上述したように、各画像点について抽出された画像特徴は、特徴ベクトルを形成することができる。例えば、MRスキャン310のすべての画像ポイントに対する特徴ベクトルのセット(例えば320)は、モデル#1(図8の訓練段階1で生成されたモデル)に入力されてもよい。言い換えれば、プロセッサ514は、図8の多段階訓練プロセスにおける訓練段階1で生成されたモデル#1をMRスキャン310から抽出された一連の特徴ベクトルに適用することができる。モデル#1は、予測結果#1 940(例えば、モデル#1のタイプに応じて、疑似CT値、セグメンテーションマップなど)を出力することができる。その後、予測結果#1 940とMRスキャン310とを組み合わせて別の画像特徴抽出311を行うことができる。より多くの情報が予測結果#1 940によって提供されるので、より多くの画像特徴が第2の抽出から生じ得るか、又は第2の抽出からの画像特徴結果が第1の抽出からの結果よりも優れた品質を有し得る。第2の抽出からの画像特徴結果は、図8の訓練段階2で生成されたモデル#2に入力されて予測結果#2 950を生成することができるMR特徴ベクトルのセットを形成することができ、MRスキャン310と組み合わせることができる。画像特徴抽出モジュール311は、組み合わされた予測結果#2 950及びMRスキャン310からMR特徴ベクトルのセットを再び抽出することができる。このプロセスは、最終予測モデルモデル#Mが疑似CT値340を生成するために適用されるまで繰り返される。いくつかの実施形態によれば、モデル#1からモデル#Mの予測モデルは、訓練プロセスで生成されたモデルと同じ順序で適用されるべきである。疑似CT値340は、患者の解剖学的形状の正確な視覚的表現を示す疑似CT画像350を生成するために使用することができる。
[0103]
いくつかの実施形態では、予測モデルのモデル#1、モデル#2、モデル#M−1のうちの任意のひとつは、組織分類モデルであってもよい。例えば、モデル#1は分類モデルであってもよい。MRスキャン310から抽出された特徴ベクトルにモデル#1を適用することにより、分類マップ(例えば、組織分類マップ)を生成することができる。モデル分類マップは、MRスキャン310と共に、画像特徴抽出モジュール311への入力として使用されて、モデル#2が適用され得るより多くの情報又は改善された品質を有する特徴ベクトルを生成することができる。特徴ベクトルをさらに精細化するために、同様のステップを繰り返すことができる。一連の分類モデルは、多段階訓練プロセスから提供されてもよい(図8)。例えば、モデル#1は、分類モデルAであってもよく、モデル#2は、分類モデルBであってもよい。分類モデルA及び分類モデルBは、同一の組織クラスであるが異なる精緻に関連してもよく、又は異なる組織クラスに関連してもよい。上述のように、最終予測モデル#Mは、開示された実施形態により、疑似CT値340を生成する疑似CT予測モデルである。
[0104]
いくつかの実施形態では、予測モデルは、マルチチャネルMRスキャン及び対応するCT値を含む訓練データ110(図1に示す)で構築及び訓練することができる。マルチチャネルMRスキャンは、シングルチャネルMRスキャンより多くの情報を提供する。マルチチャネルMR画像で利用可能な情報の増加は、疑似CT画像を生成するためにCT値のより正確でより堅牢な予測を可能にする。例えば、マルチチャネルMR画像は、MR強度値から固有の組織パラメータ値への変換を可能にする。
[0105]
MRIは、磁気及び無線周波数(RF)フィールドの操作を通じて、人体の様々な特性(例えば、構造的及び機能的の両方)の研究を可能にする、非常に汎用性の高いイメージング技術である。標準的な構造(又は解剖学的)イメージングのために、測定されたMR信号(例えば、MR画像強度)は、いくつかの固有の組織パラメータ、すなわちプロトン密度(P)、縦緩和時間(T)及び横緩和時間(T又はT (磁場の不均一性の影響を考慮した場合))の関数である。例えば、FLASH及びSPGRイメージングプロトコル(例えばイメージングシーケンスとしても知られている)の両方について、MR信号強度(S)は、式2に従って固有の組織パラメータ(P、T、及びT )の関数として表すことができる。
ここで、TR、TE、及びαは、ユーザが自由に変更できるMR取得パラメータである。異なるパラメータを使用して、異なる画像コントラストを生成することができる。
[0106]
MR画像からCT数を予測するために、予測モデルは、シーケンスパラメータ(TR、TE、及びα)に依存するMR信号強度Sの代わりに、主に固有組織パラメータ(P、T、及びT )に依存してもよい。何故ならシーケンスパラメータ(TR、TE、及びα)は患者の解剖学的特性を少なくとも直接には表していないからである。
[0107]
マルチチャネルMR画像を使用すると、シーケンスパラメータ(TR、TE、及びα)の異なる設定を有する各画像を有する複数の画像を提供することができるため、これらの固有の組織パラメータの推定を可能にする。したがって、マルチチャネルMR画像は、Sの複数の値とシーケンスパラメータ(TR、TE、及びα)が既知である式2を解くことによって、固有の組織パラメータ(P、T、及びT )の推定を可能にする。例えば、3つの未知のパラメータ(P、T、及びT )すべて推定するために、3つのMR画像(例えば、3つのチャネル)が必要とされる。追加のチャネルを使用することは、画像ノイズを低減することによってパラメータ推定のロバスト性(robustness)を改善する。
[0108]
いくつかの実施形態では、開示された方法を実施するためのワークフローは、訓練(例えば、モデル構築)ステージと適用(例えば、疑似CT生成)ステージの2つのステージを含む。いくつかの実施形態では、訓練ステージは、訓練データが収集された後に一度計算されるだけでよい。(c)
予測モデルを訓練した後、適用ステージで、訓練されたモデルを新しいマルチチャネルMRスキャンに適用して、マルチチャネルMRスキャンのみを有する新しい患者の疑似CT画像を生成することができる。以下の説明では、図10及び図11は訓練ステージについてでありし、図12は予測ステージについてである。
[0109]
図10は、マルチチャネルMRスキャンから推定された組織パラメータを用いて疑似CT予測モデルを訓練するための例示的プロセスの図を示す。図1に示すプロセスと同様に、CTスキャンとマルチチャネルMRスキャンの両方を有する複数の訓練対象からの訓練データを収集することができる。例えば、図10に示すように、被検体aに対して、マルチチャンネルMR画像データ1010a及びCT画像データ1014aを収集することができる。マルチチャネルMRデータ1010aは、異なるシーケンスパラメータセットを使用して得られた複数のMR画像を含むことができる。同様に、被検体bについて、マルチチャネルMR画像データ1010b、CT画像データ1014bを取得することができる。マルチチャネルMR画像データ及びCT画像データは、画像取得装置(例えば、440)又は画像データベース(例えば、520,530)から収集することができる。データ取得/収集のプロセスは、n番目のデータセット(例えば1010nと1014n)が訓練データに含まれるまで続きます。一実施形態では、図10に示す訓練データは、図10のMR画像データがマルチチャネルMRデータであることを除いて、図1の訓練データ110と同様であるが、一方、図1のMRデータは、チャネル又はマルチチャネルのMRデータのいずれかである。
[0110]
いくつかの実施形態では、CTスキャン(例えば、1014a)及びマルチチャネルMRスキャン(例えば、1010a)が整列される。そうでない場合には、自動又は半自動画像登録又は整列手順を適用してそれらを整列させることができる。上述のように、CT画像とMR画像の位置合わせされたペアは、各ボクセル(例えば、画像の空間的位置を示す)について、対応するCT及びMR画像値が既知であるか、又は対応関係が分かっていることを意味する。加えて、MR画像は、幾何学的歪みを補正するために何らかの手順を経ることもある。
[0111]
本出願は、訓練データを使用して予測モデルを構築して訓練する学習ベースのアプローチを開示する。一実施形態では、予測モデルは、予測モデルの出力が連続変数(例えば、CT値)であり得るように、回帰モデル又は回帰関数であり得る。
[0112]
MR画像から得られた特徴に基づいてCT画像強度(CT値又はCT数としても知られている)を予測することができる予測モデルを構築及び/又は訓練するために、多くの統計学的又は機械学習法を使用することができる。例えば、教師付き学習は、学習データの集合に基づいて予測モデルを決定するために使用することができる機械学習のブランチである。訓練データの各サンプルは、入力データ(例えば、測定値又は特徴のベクトル)及び所望の出力値(例えば、監視信号)を含むペアである。教師付き学習アルゴリズムは、疑似CT画像を生成するアプリケーションにおいて、通常、出力変数が数値又は連続である場合、訓練データを分析して予測関数(例えば、回帰関数)を生成することができる。様々なアルゴリズムを適用して、サポートベクトルマシン、ニューラルネットワーク、決定木、及びランダムフォレストを含むが、これらに限定されない予測モデルを決定することができる。
[0113]
訓練データを用いて訓練された予測モデルは、同一又は異なる患者の任意の新しいマルチチャネルMRスキャンセットに対して疑似CT画像を生成するために使用することができる。
[0114]
本願の実施形態は、MR強度値(例えば、S)を固有の組織パラメータ(例えば、P、T、T 、又はT )に変換し、固有の組織パラメータに基づいて予測モデルを構築することができる。上述のように、MRイメージングを使用する能力は、CTイメージングよりも大きなフレキシビリティとより少ない放射線被曝を提供するが、MR強度値は、シーケンス依存性であるため、投与量計算に直接使用することはできない。生のMR強度値を使用する代わりに、固有の組織パラメータを使用して予測モデルを訓練することは、シーケンスに依存しない予測モデルを提供することができる。配列に依存しないことは、画像化配列又は配列パラメータが容易に変更可能であり、しばしば異なる診療所間で著しく異なるため、有利であり得る。予測モデルを配列に依存しないように設計することにより、MRシーケンスを用いて固有の組織パラメータを推定することができるならば、異なるMRスキャナ、異なるMR画像シーケンス、又は異なる診療所から取得したデータを一緒に使用することができる。さらに、新しい患者のMRイメージングシーケンスは、訓練データによって使用されるMRイメージングシーケンスと同じである必要はない。したがって、ユーザは、予測モデルを訓練するために新たな訓練データを取得する必要なく、将来の患者のために新しいMR画像シーケンスを自由に設計することができる。
[0115]
組織パラメータに基づいて予測モデルを構築するには、MR画像強度を各患者の組織パラメータ値に変換する必要がある。これは、患者のMR画像の各画像点(例えば、ボクセル)においてMR画像化方程式(例えば、式2)を解くことによって達成することができる。組織パラメータ画像(組織パラメータマップとも呼ばれる)のセットを生成することができる。このセットは、各組織パラメータのひとつの組織パラメータマップを含むことができる。例えば、セットは、Pのマップと、Tのマップと、T又はT のマップを含むことができる。組織パラメータ値は、患者の体の下にある組織又は器官の特性を反映する固有の値である。さらに、CT画像は各訓練対象についてMR画像と整列されるので、CT画像はMR画像から生成された組織パラメータマップとさらに整列する。
[0116]
図10に示すように、各訓練対象の組織パラメータマップは、マルチチャネルMR画像データに基づいて生成されてもよい。例えば、被験体aの組織パラメータマップ1012aは、3つの組織パラメータマップのすべてのセットを含むことができる。いくつかの実施形態では、マルチチャネルMR画像からのすべての組織パラメータの推定が困難な場合、モデルは、組織パラメータのサブセットのみを使用して構築することもできる。各訓練対象の組織パラメータマップ(例えば、1012b、...、1012n)が得られると、訓練モジュール1045を使用して疑似CT予測モデル1050を構築して訓練するために、組織パラメータマップを対応するCT画像と共に使用することができる(例えば、統計的又は機械的学習技術を使用して)。いくつかの実施形態では、各画像点(例えば、ボクセル)について、組織パラメータのセットは、図1で説明したように特徴ベクトルに含まれる特徴として扱うことができる。例えば、3つの組織パラメータのすべてが使用される場合、[P、T、T ]の特徴ベクトルを構成することができる。この組織パラメータ特徴ベクトルは、予測モデルを構築し訓練するために、図2に関連して説明したように、単独で、又は他の特徴と組み合わせて使用することができる。上述の予測モデル150を構築し訓練する開示された技術は、予測モデル1050を構築し訓練するプロセスにも適用可能である。
[0117]
図11は、マルチチャネルMRデータを使用して疑似CT予測モデルを構築する例示的な方法1100のフローチャートである。方法1100は、システム500によって実装することができる。ステップ1110において、プロセッサ514は、マルチチャネルMRデータ(例えば、1010a、1010bなど)及びCTデータ(例えば、1014a、1014bなど)を含む訓練データを受信することができる。マルチチャネルMRデータは、1人又はそれ以上の患者のマルチチャネルMR画像を含むことができる。いくつかの実施形態では、マルチチャネルMRデータは、各患者について少なくとも2つのマルチチャネルMR画像を含むことができる。同じ患者のマルチチャネルMR画像は、異なる撮像シーケンスパラメータを使用して得ることができる。CTデータは、1人以上の患者のCT画像を含むことができる。いくつかの実施形態では、CTデータは、各患者について少なくともひとつのCT画像を含むことができる。複数のCT画像が利用可能であれば、それらを平均して画像ノイズを低減することができる。所与の患者に対して、CT画像及びマルチチャネルMR画像は整列され得る。いくつかの実施形態では、CT画像及びマルチチャネルMR画像が位置合わせされていない場合、画像位置合わせ技術を使用してそれらを位置合わせすることができる。
[0118]
ステップ1120において、プロセッサ514は、マルチチャネルMRデータからのMR強度に基づいて、少なくともひとつの組織パラメータマップ(例えば、1012a、1012bなど)を決定することができる。P、T、T の少なくともひとつの組織パラメータマップは、数式2を解くことによって患者の複数のMR画像のMR強度値を用いて推定することができる。数式2は非線形方程式であるので、フィッティング技法を使用して、マルチチャネルMRデータに基づいて既知である複数のMR強度(S)及びシーケンスパラメータ(TR、TE、及びα)に基づいて組織パラメータ値を推定することができる。上述のように、3つの組織パラメータマップのすべてが好ましいが、組織パラメータマップのサブセットが使用されてもよい。
[0119]
いくつかの実施形態では、組織パラメータマップは、個々の画像点を推定することによって生成することができる。例えば、いくつかの種類の組織パラメータ(例えば、P、T、T 、又はT )を含む組織パラメータ値のセットは、あらゆる画像ポイントで推定されてもよい。対応する組織パラメータマップは、特定の種類の全ての組織パラメータ値の集合として形成することができる。
[0120]
ステップ1130において、プロセッサ514は、ステップ1120で生成された組織パラメータマップに対応するCT値を取得することができる。いくつかの実施形態では、CT値は、組織パラメータマップと同じ患者のCT画像のCT強度値であってもよい。上述したように、CT画像はMR画像と位置合わせされるので、CT画像はまた、MR画像から変換された組織パラメータマップと位置合わせされる。
[0121]
ステップ1140において、プロセッサ514は、CT値及び組織パラメータマップに基づいて疑似CT予測モデル(例えば、1050)を生成することができる。いくつかの実施形態では、疑似CTモデル1050を訓練するために、CT値及び組織パラメータマップを訓練モジュール(例えば、1045)に入力することができる。統計的学習又は機械学習技術のような回帰方法は、予測モデルを訓練するために訓練モジュール1045によって使用されてもよい。訓練された予測モデル1050は、ひとつ又はそれ以上の組織パラメータ値(例えば、P、T、T 、T )に基づいてCT値を予測するために使用され得る数学的又は統計的モデルであり得る。上述のように、全ての組織パラメータを使用することが好ましいが、組織パラメータのサブセットのみを使用してモデルを構築することもできる。
[0122]
図12は、マルチチャネルMR画像を用いて患者の疑似CT画像を生成するための方法例1200のフローチャートである。方法1200は、システム500によって実施され、方法1100によって構築された予測モデル(例えば、1050)を使用して新しい患者の疑似CT画像を生成するために使用されてもよい。
[0123]
ステップ1210において、プロセッサ514は、例えば画像取得装置440又はデータベース520,530から患者(例えば、新しい患者)のマルチチャネルMR画像を受信することができる。受信されたマルチチャネルMR画像は、対応するCT画像を有していなくてもよい。ステップ1220において、プロセッサ514は、マルチチャネルMR画像を少なくともひとつの組織パラメータマップに変換することができる。例えば、プロセッサ514は、フィッティング技術などを用いて数式2を解くことによって、マルチチャネルMR画像を組織パラメータマップに変換することができる。P、T、T 、及び/又はT のマップのような、ひとつ又はそれ以上の組織パラメータマップが、変換プロセス1220によって生成されてもよい。上述のように、すべての組織パラメータマップを有することが好ましいが、いくつかの実施形態では、マップのサブセットもまた使用され得る。
[0124]
ステップ1230において、プロセッサ514は、組織パラメータマップに予測モデル(例えば、1050)を適用することができる。例えば、変換された組織パラメータマップは、予測モデルへの入力として使用されてもよい。いくつかの実施形態では、予測モデルは、ポイントモデル方式で動作することができ、ここで、予測モデルは、その場所で計算された組織パラメータ値に基づいて患者画像の各位置(例えば、疑似CT画像)でのCT値を予測する。より複雑な予測モデルを使用することもできる。例えば、モデルは、各点について近傍の組織パラメータ値を考慮に入れることができ、データノイズに関する頑健性を改善することができる。モデルはまた、組織パラメータ値と、画像ポイント位置のような他の情報、又は、例えば、テクスチャ、グラジエントなどの組織パラメータマップから導き出され得る他の特徴との組み合わせに基づいて予測を行うことができる。
[0125]
ステップ1240において、患者の疑似CT画像は、ステップ1230から得られた疑似CT値を組み立てることによって生成することができる。
[0126]
予測モデルを訓練し、モデルを使用して疑似CT値を予測するには、線形回帰又は一般最小二乗回帰などのパラメトリック手法を使用することができる。これらのパラメトリック法では、回帰関数は、入力データから推定できる有限数の未知モデルパラメータによって定義される。例えば、回帰モデルは次のように定義できる。
ここで、Hは予測されるCT値を示し、Xは入力変数のベクトル、例えば組織パラメータ値(P、T、及びT (又はT))を示し、βは回帰モデルの未知モデルパラメータのベクトルを示す。ひとつの例示的なモデルは次のように定義される線形回帰モデルである。
[0127]
訓練ステージでは、訓練データは多数の観測、例えば、(例えば、ステップ1130でのCTスキャンにより提供される)既知のH値と、(例えば、ステップ1120でのMRスキャンを変換することにより提供される)対応するP、T、及びT 値、を提供することができる。これらの観測データを使用して、モデルパラメータを(例えば、最小二乗フィッティングを使用して)計算又は訓練することができる。訓練後にβが得られると、モデルは、(例えば、方法1200を使用して予測ステージにおいて)P、T、及びT 値の新しいセットに対するHを計算するために使用することができる。
[0128]
いくつかの実施形態では、予測モデルの入力データとして組織パラメータを使用する代わりに、又はそれに加えて、他の情報も収集することができる。たとえば、計算又は収集できる余分な特徴には、以下のものが含まれるが、これらに限定されない。
・イメージポイントの座標、又は外部参照空間又はひとつ又はいくつかの内部ランドマークポイントに対する正規化座標。
・サンプルポイント位置で計算された組織パラメータマップの曲率。
・サンプルポイント位置で計算された組織パラメータマップのテクスチャ測定値。
・局所組織パラメータ値のパッチ、例えば、サンプル点の5×5×5近傍内の組織パラメータ値。
[0129]
本明細書の記述及び方法に基づいたコンピュータプログラムは、ソフトウェア開発者の技術の範囲内である。様々なプログラム又はプログラムモジュールは、様々なソフトウェアプログラミング技術を使用して作成することができる。例えば、プログラムセクション又はプログラムモジュールは、Java、Python、C、C++、アセンブリ言語、又は既知のプログラミング言語で設計することができる。このようなソフトウェアのセクション又はモジュールのひとつ又はそれ以上は、コンピュータシステム及び/又はコンピュータ読み土地可能媒体に統合することができる。
[0130]
さらに、例示的な実施形態が本明細書において説明したが、本開示に基づく等価な要素、修正、省略、(様々な実施形態にわたる態様の)組み合わせ、適合又は変更を有する任意の及びすべての実施形態を含む。特許請求の範囲の要素は、特許請求の範囲で使用されている言語に基づいて広く解釈されるべきであり、本明細書に記載された実施例に限定されるものではなく、アプリケーションの実行中に限定されるものではない。したがって、明細書及び実施例は、単なる例示としてみなされ、真の範囲及び精神は、添付の特許請求の範囲及びそれらの等価物の全範囲によって示されるものである。

Claims (20)

  1. 疑似CT画像を生成するシステムであって、
    複数の訓練対象の磁気共鳴(MR)データ及びコンピュータ断層撮影(CT)データを含む訓練データを格納するように構成されたデータベースであって、各訓練対象が少なくともひとつのMR画像及び少なくともひとつのCT画像を有するデータベースと、
    前記データベースに格納された情報にアクセスするために前記データベースに通信可能に結合されたコンピュータプロセッサ回路装置と、
    前記コンピュータプロセッサ回路装置に通信可能に結合された記憶装置とを有し、
    前記記憶装置は、前記コンピュータプロセッサ回路装置によって実行されると、前記コンピュータプロセッサ回路装置が操作を実行するように構成された命令を格納しており、
    前記操作は、
    前記データベースにアクセスして、前記複数の訓練対象のそれぞれについて、前記少なくともひとつのMR画像及び前記少なくともひとつのCT画像を含む訓練データを取り出すステップと、
    各訓練対象について、前記少なくともひとつのMR画像の各画像点から複数の特徴を抽出し、前記抽出された特徴のそれぞれを前記特徴ベクトルの複数の次元のうちのそれぞれの次元で格納することにより、前記抽出された特徴に基づいて各画像の特徴ベクトルを生成し、前記少なくともひとつのCT画像の各画像点からCT値を抽出し、前記複数の訓練対象のそれぞれに対して前記少なくとも1つのMR画像と前記少なくとも1つのCT画像とを位置合わせするステップと、
    CT値パラメータH、組織値パラメータ、およびβパラメータの間の関係をモデル化するために機械学習技術を訓練するステップであって、前記訓練することは、前記CT値パラメータHを各画像点の前記抽出されたCT値と表示し、前記組織値パラメータを各画像点の前記特徴ベクトルと表示し、前記Hが前記βパラメータおよび前記組織値パラメータの数学的演算にほぼ等しくなるように、前記表示されたCT値パラメータHおよび前記組織値パラメータに対してデータフィッティング技法を実行することにより前記βパラメータを計算することを含むステップと、
    前記訓練された機械学習技術に基づいて、新しいMR画像の1つまたは複数の特徴から疑似CT画像を生成するステップであって、前記機械学習技術によってモデル化された前記組織値パラメータを前記新しいMR画像の前記1つまたは複数の特徴であると表示し、前記機械学習技術によってモデル化された前記計算されたβパラメータを前記新しいMR画像の前記1つまたは複数の特徴であると表示された前記組織値パラメータに適用することにより前記機械学習技術によりモデル化された前記CT値パラメータHを計算することにより、前記疑似CT画像を生成するステップとを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  2. 請求項1記載のシステムにおいて、
    前記訓練データは、前記複数の訓練対象のうちの少なくともひとつに関する複数のCT画像を含み、
    前記操作は、前記複数のCT画像の対応する画像点のCT値を平均化するステップを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  3. 請求項1記載のシステムにおいて、
    前記操作は、前記複数の訓練対象のそれぞれに対する前記少なくともひとつのMR画像と前記少なくともひとつのCT画像とを整列させるステップを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  4. 請求項1記載のシステムにおいて、
    前記機械学習技術によりモデル化された関係は、次式
    ただし、Hは前記CT値パラメータHに対応し、β 、β 、β は前記βパラメータに対応し、P、T 、T は前記組織パラメータに対応する
    に対応する
    ことを特徴とする疑似CT画像を生成するシステム。
  5. 請求項1記載のシステムにおいて、
    前記操作は、
    前記特徴ベクトルの次元を削減するステップと、
    削減された次元数を有する前記特徴ベクトルと前記複数の訓練対象の前記CT値とに基づいて前記疑似CT画像を生成するステップとを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  6. 請求項5記載のシステムにおいて、
    前記特徴ベクトルの次元を削減するステップは、
    前記特徴ベクトルの雑音又は冗長性を低減することによって、前記複数の訓練対象の前記MR画像に基づいて前記特徴ベクトルの次元を削減するステップを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  7. 請求項5記載のシステムにおいて、
    前記特徴ベクトルの次元を削減するステップは、
    前記特徴ベクトル内の前記CT画像に無関係な信号を除去することによって、前記複数の訓練対象の前記MR画像と前記CT画像に基づいて前記特徴ベクトルの次元を削減するステップを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  8. 請求項1記載のシステムにおいて、
    前記抽出された特徴が、局所パターン特徴、ランドマーク特徴、又は文脈特徴のうちの少なくともひとつを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  9. 疑似CT画像を生成するシステムであって、
    コンピュータプロセッサ回路装置と、
    前記コンピュータプロセッサ回路装置に通信可能に結合された記憶装置を有し、
    前記記憶装置は、前記コンピュータプロセッサ回路装置によって実行されると、前記コンピュータプロセッサ回路装置が操作を実行するように構成された命令を格納しており、
    前記操作は、
    患者の磁気共鳴(MR)画像を受信するステップと、
    機械学習技術は、コンピュータ断層撮影(CT)値パラメータH、組織値パラメータ、βパラメータの間の関係をモデル化するために訓練されており、
    1つ又は複数のMR画像の各画像点から複数の特徴を抽出するステップであって、各抽出された特徴は異なる種類のMR画像情報に対応するステップと、
    前記抽出された特徴に基づいて各画像に対する特徴ベクトルを作成するステップと、
    1つまたは複数のMR画像の各画像点についてCT値を抽出するステップと、
    複数の訓練対象のそれぞれについて前記1つまたは複数のMR画像と前記抽出されたCT値に対応するCT画像を位置合わせするステップと、
    前記CT値パラメータHを前記抽出されたCT値と表示し、前記組織値パラメータを各画像点の前記特徴ベクトルと表示するステップと、
    前記Hが前記βパラメータおよび前記組織値パラメータの数学的演算にほぼ等しくなるように、前記表示されたCT値パラメータHおよび前記組織値パラメータに対してデータフィッティング技法を実行することにより前記βパラメータを計算するステップと、
    すべての画像点に対する前記機械学習技術に基づいて前記疑似CT画像を生成するステップであって、前記機械学習技術によってモデル化された前記組織値パラメータを前記新しいMR画像の前記1つまたは複数の特徴であると表示し、前記機械学習技術によってモデル化された前記計算されたβパラメータを前記新しいMR画像の前記1つまたは複数の特徴であると表示された前記組織値パラメータに適用することにより前記機械学習技術によりモデル化された前記CT値パラメータHを計算することにより、前記疑似CT画像を生成するステップとを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  10. 請求項9記載のシステムにおいて、
    前記異なる種類のMR画像情報は、少なくとも1つの数値と少なくとも1つのカテゴリカル特性とを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  11. 請求項10記載のシステムにおいて、
    前記少なくとも1つの数値は、強度値、エネルギー、コントラスト、均質性、相関関係、関節特徴、または特徴の座標位置のうちの少なくとも1つを含み、前記少なくとも1つのカテゴリカル特性は、組織タイプ、ローカルパターン特徴、ランドマーク特徴、コンテキスト特徴、または構造ラベルの少なくとも1つを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  12. 請求項9記載のシステムにおいて、
    前記操作は、
    組織クラスに基づいて前記MR画像をセグメント化するステップと、
    前記セグメント化に基づいて分類マップを生成するステップと、
    前記分類マップからひとつ又はそれ以上の画像特徴を抽出するステップと、
    前記分類マップから抽出されたひとつ又はそれ以上の画像特徴を、前記MR画像から抽出された前記特徴と組み合わせるステップとを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  13. 患者の疑似CT画像を生成するシステムであって、
    コンピュータプロセッサ回路装置と、
    前記コンピュータプロセッサ回路装置に通信可能に結合された記憶装置を有し、
    前記記憶装置は、前記コンピュータプロセッサ回路装置によって実行されると、前記コンピュータプロセッサ回路装置が操作を実行するように構成された命令を格納しており、
    前記操作は、
    前記患者の磁気共鳴(MR)画像を受信するステップと、
    機械学習技術は、コンピュータ断層撮影(CT)値パラメータH、組織値パラメータ、βパラメータの間の関係をモデル化するために訓練され、
    少なくとも1つのMR画像の各画像点から複数の特徴を抽出するステップと、
    前記抽出された特徴のそれぞれを前記特徴ベクトルの複数の次元のうちのそれぞれの次元に格納することにより、前記抽出された特徴に基づいて各画像点に対する特徴ベクトルを作成するステップと、
    少なくとも1つのCT画像の各画像点からCT値を抽出するステップと、
    複数の訓練対象のそれぞれについて前記少なくとも1つのMR画像と前記少なくとも1つのCT画像を位置合わせするステップと、
    前記CT値パラメータHを各画像点の前記抽出されたCT値と表示し、前記組織値パラメータを各画像点の前記特徴ベクトルと表示するステップと、
    前記Hが前記βパラメータおよび前記組織値パラメータの数学的演算にほぼ等しくなるように、前記表示されたCT値パラメータHおよび前記組織値パラメータに対してデータフィッティング技法を実行することにより前記βパラメータを計算するステップと、
    前記訓練された機械学習技術に基づいて新しいMR画像の1つまたは複数の特徴から前記疑似CT画像を生成するステップであって、前記機械学習技術によってモデル化された前記組織値パラメータを前記新しいMR画像の前記1つまたは複数の特徴であると表示し、前記機械学習技術によってモデル化された前記計算されたβパラメータを前記新しいMR画像の前記1つまたは複数の特徴であると表示された前記組織値パラメータに適用することにより前記機械学習技術によりモデル化された前記CT値パラメータHを計算することにより、前記疑似CT画像を生成するステップとを含む
    ことを特徴とする疑似CT画像を生成するシステム。
  14. 請求項13記載のシステムにおいて、
    前記操作は、
    画像取得装置から前記MR画像を受信するステップを含む
    ことを特徴とする患者の疑似CT画像を生成するシステム。
  15. 請求項13記載のシステムにおいて、
    前記操作は、
    MR画像データベースから前記MR画像を受信するステップを含む
    ことを特徴とする患者の疑似CT画像を生成するシステム。
  16. 請求項13記載のシステムにおいて、
    前記新しい画像は分類画像である
    ことを特徴とする患者の疑似CT画像を生成するシステム。
  17. 請求項13記載のシステムにおいて、
    前記新しい画像は中間疑似CT画像である
    ことを特徴とする患者の疑似CT画像を生成するシステム。
  18. 請求項1記載のシステムにおいて、
    前記操作は、
    前記疑似CT画像に基づいて対象の解剖学的特徴を視覚的に表すステップを含む
    ことを特徴とする患者の疑似CT画像を生成するシステム。
  19. 請求項9記載のシステムにおいて、
    前記操作は、
    前記疑似CT画像に基づいて対象の解剖学的特徴を視覚的に表すステップを含む
    ことを特徴とする患者の疑似CT画像を生成するシステム。
  20. 請求項13記載のシステムにおいて、
    前記機械学習技術によりモデル化された関係は、次式
    ただし、Hは前記CT値パラメータHに対応し、β 、β 、β は前記βパラメータに対応し、P、T 、T は前記組織パラメータに対応する
    に対応する
    ことを特徴とする疑似CT画像を生成するシステム。
JP2018519442A 2015-10-13 2016-10-12 特徴回帰モデルを用いたmrデータからの疑似ct生成 Active JP6567179B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US14/881,903 2015-10-13
US14/881,903 US10307108B2 (en) 2015-10-13 2015-10-13 Pseudo-CT generation from MR data using a feature regression model
PCT/US2016/056538 WO2017066247A1 (en) 2015-10-13 2016-10-12 Pseudo-ct generation from mr data using a feature regression model

Publications (3)

Publication Number Publication Date
JP2018530401A JP2018530401A (ja) 2018-10-18
JP2018530401A5 JP2018530401A5 (ja) 2019-06-27
JP6567179B2 true JP6567179B2 (ja) 2019-08-28

Family

ID=57227090

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018519442A Active JP6567179B2 (ja) 2015-10-13 2016-10-12 特徴回帰モデルを用いたmrデータからの疑似ct生成

Country Status (7)

Country Link
US (3) US10307108B2 (ja)
EP (2) EP3362984B1 (ja)
JP (1) JP6567179B2 (ja)
CN (1) CN108770373B (ja)
AU (1) AU2016338923B2 (ja)
RU (1) RU2703344C1 (ja)
WO (1) WO2017066247A1 (ja)

Families Citing this family (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10102451B2 (en) * 2015-10-13 2018-10-16 Elekta, Inc. Pseudo-CT generation from MR data using tissue parameter estimation
US10307108B2 (en) 2015-10-13 2019-06-04 Elekta, Inc. Pseudo-CT generation from MR data using a feature regression model
WO2017103237A1 (en) * 2015-12-18 2017-06-22 Koninklijke Philips N.V. Method for determining a patient specific locally varying margin
DE112016006218B4 (de) * 2016-02-15 2022-02-10 Mitsubishi Electric Corporation Schallsignal-Verbesserungsvorrichtung
US10794977B2 (en) * 2016-06-23 2020-10-06 Siemens Healthcare Gmbh System and method for normalized reference database for MR images via autoencoders
US11478212B2 (en) * 2017-02-16 2022-10-25 Siemens Healthcare Gmbh Method for controlling scanner by estimating patient internal anatomical structures from surface data using body-surface and organ-surface latent variables
CN110692107B (zh) * 2017-05-31 2023-11-14 皇家飞利浦有限公司 用于临床决策支持的对原始医学成像数据的机器学习
US11276208B2 (en) * 2017-06-12 2022-03-15 The Research Foundation For The State University Of New York System, method and computer-accessible medium for ultralow dose computed tomography image reconstruction
US10753997B2 (en) * 2017-08-10 2020-08-25 Siemens Healthcare Gmbh Image standardization using generative adversarial networks
CN108416846A (zh) * 2018-03-16 2018-08-17 北京邮电大学 一种无标识三维注册算法
EP3543911A1 (en) * 2018-03-22 2019-09-25 Koninklijke Philips N.V. Anomaly detection using magnetic resonance fingerprinting
EP3583982A1 (en) * 2018-06-18 2019-12-25 Koninklijke Philips N.V. Therapy planning device
EP3617733A1 (en) * 2018-08-30 2020-03-04 Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. Method and apparatus for processing magnetic resonance data using machine learning
EP3633622A1 (en) * 2018-10-02 2020-04-08 Koninklijke Philips N.V. Generation of pseudo radiographic images from optical images
CN109710793B (zh) * 2018-12-25 2021-08-17 科大讯飞股份有限公司 一种哈希参数确定方法、装置、设备及存储介质
CN109785405A (zh) * 2019-01-29 2019-05-21 子威·林 利用多组磁共振图像的多变量回归生成准ct图像的方法
US10467504B1 (en) * 2019-02-08 2019-11-05 Adhark, Inc. Systems, methods, and storage media for evaluating digital images
US11531840B2 (en) 2019-02-08 2022-12-20 Vizit Labs, Inc. Systems, methods, and storage media for training a model for image evaluation
CN109918519A (zh) * 2019-02-26 2019-06-21 重庆善功科技有限公司 一种面向海量延时摄影胚胎图像存储与查询的方法及***
CN109893784A (zh) * 2019-04-19 2019-06-18 深圳先进技术研究院 一种实现超声穿颅聚焦的方法以及电子设备
EP3958965A4 (en) 2019-04-26 2023-01-11 Spectronic AB GENERATION OF SYNTHETIC ELECTRON DENSITY IMAGES FROM MAGNETIC RESONANCE IMAGES
CN110084318B (zh) * 2019-05-07 2020-10-02 哈尔滨理工大学 一种结合卷积神经网络和梯度提升树的图像识别方法
CN110270015B (zh) * 2019-05-08 2021-03-09 中国科学技术大学 一种基于多序列MRI的sCT生成方法
JP7451293B2 (ja) 2019-06-13 2024-03-18 キヤノンメディカルシステムズ株式会社 放射線治療システム
US11024034B2 (en) * 2019-07-02 2021-06-01 Acist Medical Systems, Inc. Image segmentation confidence determination
GB2586791B (en) * 2019-08-30 2022-11-16 Elekta ltd Pseudo-CT image generation
US11562482B2 (en) 2020-03-30 2023-01-24 Varian Medical Systems International Ag Systems and methods for pseudo image data augmentation for training machine learning models
US12008695B2 (en) * 2020-09-25 2024-06-11 GE Precision Healthcare LLC Methods and systems for translating magnetic resonance images to pseudo computed tomography images
CN112435554B (zh) * 2020-11-25 2022-04-12 皖南医学院 一种ct教学模拟***及其控制方法
CN112509022A (zh) * 2020-12-17 2021-03-16 安徽埃克索医疗机器人有限公司 一种术前三维影像与术中透视图像的无标定物配准方法
CN114692665B (zh) * 2020-12-25 2024-05-24 西南电子技术研究所(中国电子科技集团公司第十研究所) 基于度量学习的辐射源开集个体识别方法
US11966454B2 (en) * 2021-10-28 2024-04-23 Shanghai United Imaging Intelligence Co., Ltd. Self-contrastive learning for image processing
CN114190889A (zh) * 2021-11-19 2022-03-18 上海联影智能医疗科技有限公司 心电信号的分类方法、***、电子设备及可读存储介质
CN114638745B (zh) * 2022-03-16 2023-08-18 江南大学 一种基于多借鉴信息的医学影像智能转换方法
WO2023215726A2 (en) * 2022-05-01 2023-11-09 The General Hospital Corporation System for and method of planning and real-time navigation for transcranial focused ultrasound stimulation

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1271570C (zh) * 2004-06-17 2006-08-23 上海交通大学 核磁共振多谱图像分割方法
CN101061490A (zh) * 2004-11-19 2007-10-24 皇家飞利浦电子股份有限公司 利用支持向量机(svm)在计算机辅助检测(cad)中进行假阳性降低的***和方法
EP2053972B1 (en) * 2006-08-17 2013-09-11 Koninklijke Philips Electronics N.V. Computed tomography image acquisition
US8068654B2 (en) * 2007-02-02 2011-11-29 Siemens Akteingesellschaft Method and system for detection and registration of 3D objects using incremental parameter learning
CN101959454B (zh) * 2008-03-07 2013-11-20 皇家飞利浦电子股份有限公司 通过磁共振图像的自动分割进行ct替代
WO2010035163A1 (en) * 2008-09-29 2010-04-01 Koninklijke Philips Electronics, N.V. Method for increasing the robustness of computer-aided diagnosis to image processing uncertainties
DE102008058488B4 (de) * 2008-11-21 2018-09-20 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Aufbereitung von kombinierten MR-Emissionstomographieaufnahmen
WO2011015957A1 (en) * 2009-08-06 2011-02-10 Koninklijke Philips Electronics N.V. Method and apparatus for generating computed tomography images with offset detector geometries
US8774482B2 (en) 2010-05-20 2014-07-08 Siemens Aktiengesellschaft Generating pseudo-CT image volumes from ultra-short echo time MR
US8706668B2 (en) * 2010-06-02 2014-04-22 Nec Laboratories America, Inc. Feature set embedding for incomplete data
JP2014528284A (ja) * 2011-09-30 2014-10-27 チルドレンズ ホスピタル メディカル センター コンピュータ断層撮影(ct)放射線量の一貫性のある検証可能な最適化の方法
US9256965B2 (en) * 2013-01-30 2016-02-09 Impac Medical Systems, Inc. Method and apparatus for generating a derived image using images of different types
CN105263407B (zh) * 2013-04-02 2018-09-14 皇家飞利浦有限公司 使用磁共振成像的对骨组织的探测
US9235781B2 (en) * 2013-08-09 2016-01-12 Kabushiki Kaisha Toshiba Method of, and apparatus for, landmark location
US10307108B2 (en) 2015-10-13 2019-06-04 Elekta, Inc. Pseudo-CT generation from MR data using a feature regression model
KR102426894B1 (ko) 2016-10-26 2022-08-01 나이키 이노베이트 씨.브이. 자동화된 신발류 플랫폼을 위한 변형 가능한 레이스 가이드
KR102187717B1 (ko) 2017-03-15 2020-12-07 나이키 이노베이트 씨.브이. 케이블 및 갑피 텐셔너를 구비하는 자동화된 신발

Also Published As

Publication number Publication date
AU2016338923A1 (en) 2018-04-26
US20220079529A1 (en) 2022-03-17
US11944463B2 (en) 2024-04-02
EP3362984A1 (en) 2018-08-22
JP2018530401A (ja) 2018-10-18
US11234654B2 (en) 2022-02-01
US20170100078A1 (en) 2017-04-13
AU2016338923B2 (en) 2019-06-27
RU2703344C1 (ru) 2019-10-16
CN108770373B (zh) 2022-08-19
WO2017066247A1 (en) 2017-04-20
CN108770373A (zh) 2018-11-06
US20190209099A1 (en) 2019-07-11
US10307108B2 (en) 2019-06-04
EP3940624A1 (en) 2022-01-19
EP3362984B1 (en) 2023-07-12

Similar Documents

Publication Publication Date Title
JP6567179B2 (ja) 特徴回帰モデルを用いたmrデータからの疑似ct生成
US10664723B2 (en) Pseudo-CT generation from MR data using tissue parameter estimation
US10149987B2 (en) Method and system for generating synthetic electron density information for dose calculations based on MRI
EP3123443B1 (en) Method and device for generating one or more computer tomography images based on magnetic resonance images with the help of tissue class separation
US8787648B2 (en) CT surrogate by auto-segmentation of magnetic resonance images
Nasief et al. Predicting necessity of daily online adaptive replanning based on wavelet image features for MRI guided adaptive radiation therapy
Opposits et al. Population based ranking of frameless CT-MRI registration methods
Wodzinski et al. Application of B-splines FFD image registration in breast cancer radiotherapy planning
Wodzinski et al. Usage of ICP algorithm for initial alignment in B-splines FFD image registration in breast cancer radiotherapy planning

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190521

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190521

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20190521

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20190618

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190730

R150 Certificate of patent or registration of utility model

Ref document number: 6567179

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250