JP6802177B2 - Mrcat画像の不確定領域に対するフォールバックソリューション - Google Patents

Mrcat画像の不確定領域に対するフォールバックソリューション Download PDF

Info

Publication number
JP6802177B2
JP6802177B2 JP2017548434A JP2017548434A JP6802177B2 JP 6802177 B2 JP6802177 B2 JP 6802177B2 JP 2017548434 A JP2017548434 A JP 2017548434A JP 2017548434 A JP2017548434 A JP 2017548434A JP 6802177 B2 JP6802177 B2 JP 6802177B2
Authority
JP
Japan
Prior art keywords
electron density
density map
map
alternative
calculating
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
JP2017548434A
Other languages
English (en)
Other versions
JP2018510705A5 (ja
JP2018510705A (ja
Inventor
ステフェン レニスヒ
ステフェン レニスヒ
ニコル スハーデワルツ
ニコル スハーデワルツ
スヴェン プレフルハル
スヴェン プレフルハル
ハインリッヒ スフルツ
ハインリッヒ スフルツ
トマス ブラフェルト
トマス ブラフェルト
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2018510705A publication Critical patent/JP2018510705A/ja
Publication of JP2018510705A5 publication Critical patent/JP2018510705A5/ja
Application granted granted Critical
Publication of JP6802177B2 publication Critical patent/JP6802177B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • 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
    • 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/45For evaluating or diagnosing the musculoskeletal system or teeth
    • 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/1031Treatment planning systems using a specific method of dose optimization
    • 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/4808Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]
    • G01R33/4812MR combined with X-ray or computed tomography [CT]
    • 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
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • 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/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475User input or interface means, e.g. keyboard, pointing device, joystick
    • A61B5/748Selection of a region of interest, e.g. using a graphics tablet
    • 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/543Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Rheumatology (AREA)
  • Pulmonology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Radiation-Therapy Devices (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Image Processing (AREA)

Description

本発明は、放射線治療計画のための磁気共鳴(MR)イメージングの分野に関する。より具体的には、本発明は、検査ボリュームの代わりの電子密度マップを計算するシステム及び方法に関する。
コンピュータ断層撮影法(CT)は、放射線治療計画の分野で使用され、CTイメージングによって収集される情報は電子密度(ED)に直接関連し、治療シミュレーション中の放射線量の計算の基礎として使用されることができる。
放射線被ばくを低減するために、磁気共鳴(MR)イメージングによって検査されるべき患者の電子密度(ED)マップを生成することが提案されている。磁気共鳴イメージング(MRI)は、強力な磁場及び電波を用いて、体内の器官及び内部構造の断面画像を生成する。 MRI装置によって検出される信号は、体内の特定の領域の含水率及び局所的な磁気特性に依存して変化するため、調査画像において異なる組織又は物質が互いに区別されることができる。
電子密度マップは磁気共鳴信号から導出され、磁気共鳴イメージングから導出されるこのような電子密度マップはMR-CAT画像(「MR-計算減衰」)と呼ばれる。このような画像は、コンピュータ断層撮影画像と同様の情報を表すが、X線による有害な照射は必要とされないため、「合成CT画像」と命名されてもよい。磁気共鳴ベースの放射線治療プランニングでは、MR-CAT画像から収集される情報に基づいて線量分布が計算される。
画像収集の物理学のために、磁気共鳴強度は電子密度に一意に対応しない。したがって、CT画像からこれらの電子密度マップを推定するときに一般的に行われるように、単純なルックアップ操作によって電子密度マップはMR画像から導出されることはできない。 MR画像からCTのような画像を作成するための第1のアプローチは、「MRのみに基づく放射線治療計画のための骨盤におけるディクソンに基づく軟組織及び骨の分類の評価」(Michael Helle、Nicole Schadewaldt他、ISMRM 2014要約)に記載される。この第1のアプローチは、多数の所定の組織クラスへのバルク密度の割り当てに依存している。他のアプローチは、「超短エコー時間を有するMRIシーケンスから得られるCT代替」と題される科学論文(Adam Johansson、Mikael Karlsson及びTufve Nyholm、Medical Physics 38(5) 2011年5月)に記載のように、多次元MR情報からCT値依存をモデリングすることにある。
しかしながら、これらのアプローチによって得られるMR-CAT画像は、例えば、 MR収集又は画像処理からの画像アーチファクトによる不正確な電子密度マップ等価物である。これは、その後、計算される線量分布の差につながる可能性がある。線量分布に関してMR-CATアーチファクトの実際の影響を推定することは困難である。したがって、MR画像から推定電子密度マップを作成することは問題である。
上記を考慮して、本発明の目的は、磁気共鳴イメージングから、改善された電子密度マップのためのシステム及び方法を提供することにある。
本発明によれば、目的は、請求項1の特徴を有する処理システム、請求項5の特徴を有する計画システム、請求項6の特徴を有する方法、請求項10の特徴を有する方法によって達成される。好ましい実施例はそれぞれの従属請求項に規定されている。
本開示の一態様では、検査ボリュームの代わりの電子密度マップを計算するための処理システムが提供される。処理システムは、複数のイメージングデータを使用して第1の電子密度マップを計算し、第1の電子密度マップの簡略化バージョンである第2の電子密度マップを計算し、第1の電子密度マップと第2の電子密度マップとを用いて代わりの電子密度マップを計算するように構成される。したがって、本開示は、未知のエラーとともに得られる電子密度マップを代わりの電子密度マップで置き換えることを教示する。本発明は、特に磁気共鳴イメージングによる電子密度(ED)マップの計算に関する。本発明の計算には2つの段階がある。まず、第1のEDマップが計算される。第1のEDマップの簡略化バージョンである第2のEDマップが計算される。特に、第2のEMマップは、第1のEDマップと同じ解剖学的構造の部分を表すが、細部を含まないという点で簡略化されている。一般に、第2のEDマップはアーチファクトの影響を受けにくい。次に、本発明によれば、代わりのEDマップを計算するために第2のEDマップが第1のEDマップとともに使用される。本発明の好ましい実施例では、代わりのEDマップは、第1のEDマップのアーチファクト領域を第2の(より単純な)EDマップ(のデータ)で置き換えることによって形成される。本発明の洞察は、例えば放射線治療計画のために、代わりのEDマップが第一のEDマップより信頼性が高いことにある。アーチファクト領域のより少ない細部で置き換えることは有利であり、第一のEDマップでアーチファクトを完全に除去又は回避する労力を必要としない。
本開示の他の態様では、処理システムは、代わりの電子密度マップを得るために、一つ又はそれより多くのアーチファクト領域において、第1の電子密度マップを第2の電子密度マップで置き換えるように構成される。したがって、本発明は、アーチファクトの影響を受けにくい簡略化された電子密度マップを使用することを提案する。代わりの電子密度マップは線量分布のエラーにつながる可能性があるが、これらのエラーは臨床的許容限界内にあることが知られている。言い換えれば、線量分布に及ぼす未知の影響を伴う一つ又はそれより多くのアーチファクト領域は、線量分布に対する許容可能な影響を有するアーチファクトに置き換えられてもよい。
本システムのさらなる態様では、ユーザインターフェースシステムが提供され、第1の電子密度マップを表示するように構成され、ユーザインターフェースシステムは、ユーザが第1の電子密度マップ上に一つ又はそれより多くのアーチファクト領域を識別し、画定するユーザデータを入力することを可能にするユーザ入力システムを有する。したがって、ユーザは、第1の電子密度マップをインタラクティブに分析して、どのゾーンがアーチファクト領域とみなされるべきかを決定することができる。
一態様では、第2の電子密度マップは、第1の電子密度マップの輪郭に対応し、均一な値で満たされ、特に均一な値は等価水密度である。均一な値で満たされる輪郭は、計算が容易で合理的なエラーをもたらす簡略化される第1の電子密度マップの例である。
本開示は、線量分布を計算するための計画システムを提案する。計画システムは、上述したような処理システムを有する。計画システムは、複数の画像データ、第2の電子密度マップ及び代わりの電子密度マップを使用して第1の電子密度マップを計算するように構成され、代わりの電子密度マップは前記第1及び第2の電子密度マップを使用して計算される。計画システムは、第1の電子密度マップを使用して第1の線量分布を計算し、代わりの電子密度マップを使用して第2の線量分布を再計算するように構成される放射線治療計画システムをさらに有する。
本発明は、検査ボリュームの代わりの電子密度マップを計算する方法であって、複数の画像データを用いて第1の電子密度マップを計算するステップと、第2の電子密度マップを計算するステップであって、第2の電子密度マップは第1の電子密度マップの簡略化されたバージョンである、ステップと、前記第1の電子密度マップ(80)及び前記第2の電子密度マップ(82)を用いて前記代わりの電子密度マップ(84)を計算するステップとを含む方法を教示する。
本開示の一態様では、本方法は、一つ又はそれより多くの識別されるアーチファクト領域において、第1の電子密度マップを第2の電子密度マップで置き換えて、代わりの電子密度マップを得るステップを有する。したがって、本方法は、未知のエラーを許容可能であると知られているエラーで置き換えることを提案する。
本開示のさらなる態様では、この方法は、第1の電子密度マップを表示し、一つ又はそれより多くのアーチファクト領域を識別するユーザ入力データを得るステップを有する。これにより、ユーザは、彼の判断に従って、第1の電子密度マップの一つ又はそれより多くのゾーンがアーチファクトとみなされるべきかを評価し、これらの一つ又はそれより多くのアーチファクト領域を境界付けすることができる。
本開示の一態様では、本方法は、第2の電子密度マップに、均一な密度値、特に水等価密度値を割り当てるステップを有する。水等価密度値のような均一な密度値は、許容可能なエラーにつながることが知られている。均一な値は、代わりの電子密度マップの計算のために、さらに容易に処理される。
本開示は、線量分布を計算する方法であって、第1の電子密度マップを計算するステップと、第1の電子密度マップを使用して第1の線量分布を計算するステップと、上記の方法及び/又は上記のシステムを用いて代わりの密度マップを計算するステップと、代わりの電子密度マップを使用して第2の線量分布を計算するステップとを有する、方法も提案する。したがって、この方法は、ユーザが線量分布上の一つ又はそれより多くのアーチファクト領域の潜在的影響を評価することを可能にする。
本開示はさらに、プロセッサに、提案される方法のステップを実行させる命令をコンピュータに実行させるためにここに記憶される制御論理を有するコンピュータプログラムプロダクトを提案する。本開示は、複数のイメージングデータを使用して第1の電子密度マップを計算し、第1の電子密度マップの簡略化バージョンである第2の電子密度マップを計算し、第1の電子密度マップと第2の電子密度マップとを用いて、代わりの電子密度マップを計算するように構成されるコンピュータ手段を有するコンピュータプログラムプロダクトを提案する。
本発明のこれら及び他の態様は、添付の図面を参照して以下に記載される好ましい実施形態を参照して明らかになるであろう。
本発明の好ましい実施形態によるシステムのブロック図である。 図1のシステムで得られた画像密度マップを示す。 図1のシステムで得られた他の画像密度マップを示す。 第1の画像密度マップ及び対応する線量分布を示す。 図1のシステム及び図2の方法で得られる代わりの画像密度マップ及び対応する線量分布を示す。 本開示の一態様による代わりの電子密度マップを計算する方法を示す。 本開示の一態様による放射線計画方法を示す。
図1は、検査ボリュームVの電子密度マップを推定するためのシステム1を示す。システム1は、推定電子密度マップに基づいて線量分布を計算するようにさらに構成される。
この目的のために、システム1は、検査ボリューム内の物理的特性の特定分布を収集するための複数の画像データ20を収集するためのイメージングシステム10を有する。イメージングシステム10は、記載される実施形態において、当業者に知られているように、磁気共鳴(MR)イメージングシステムである。この場合、物理的特性は磁気共鳴強度である。しかしながら、イメージングシステムは、イメージングデータが電子密度マップを推定するために使用され得る他のタイプのイメージングシステムになり得る。
システム1は、複数の画像データ20を処理するように構成される処理システム30をさらに有する。処理システム30は、収集される複数の画像データ20に基づいて第1の電子密度マップ80を生成するように構成される。第1の電子密度マップ80は骨及び軟部組織の内部構造を示す。第1の電子密度マップ80は、当業者に公知の計算技術を使用して計算されることができる。
第1電子密度マップ80の一例が図3Aに示される。当業者は、第1の電子密度マップ80が、イメージングシステム10及び当業者に知られている計算イメージング技術を使用して可能な限り高い詳細なレベルを有するように計算されることができることを理解するであろう。冒頭で説明されたように、第1の電子密度マップ80は、アーチファクト領域90と識別される一つ又はそれより多くのゾーンを有する、MR-CAT画像(「MR-計算減衰」)に対応する。図3Aに示されるように、第1の電子密度マップ80は、アーチファクト値91を備えるアーチファクトを有する1つのアーチファクト領域90と、図3Aの輝点として示される近傍値92を備えるその近傍とを有する。このアーチファクト領域90は、磁気共鳴収集からもたらされ、又は複数のイメージングデータの処理中及び第1の電子密度マップ80の計算中の割り当てエラーからもたらされる。
処理システム30は、第2の電子密度マップ82を生成するようにさらに構成される。第2の電子密度マップ82は、第1の電子密度マップ80と比較すると、細部のより低いレベルだがアーチファクトのより少ないリスクを有する、簡略化電子密度マップである。記載される実施形態では、第2の電子密度マップ82は、第1の電子密度マップ80の輪郭92を表すことができる。第2の電子密度マップ82は、均一な密度値93で満たされることができる。第2の電子密度マップ82の例は図3Bに示される。この例において、第2の電子密度マップ82は、輪郭が描かれ、水分等価密度値で満たされるセグメント化体を表す。これは非限定的な例であり、他の密度等価値が使用されてもよい。
処理システム30はさらに、第1の電子密度マップ80及び第2の電子密度マップ82を使用して第3の電子密度マップ84を導出するように構成される。より正確には、処理システム30は、相関される第3の電子密度マップ84を得るために、第1の電子密度マップ80において識別される一つ又はそれより多くのアーチファクト90において、第2の電子密度マップ82によって第1の電子密度マップ80を置換する。すなわち、アーチファクト値91を備えるアーチファクト及び第1の電子密度マップ80における近傍値92を備えるその近傍は、識別されるアーチファクトの第2の電子密度マップ82からの第2の電子密度マップ密度値93によって置換される。第2電子密度マップ密度値93は、置換値93である。
ユーザインタフェースシステム40は、システム1に設けられ、ユーザへのインタラクティブな表示のために構成される。ユーザインターフェースシステム40は、処理システム30内に統合されてもよく、又は2次計算デバイスであってもよい。
ユーザインターフェースシステム40は、第1電子密度マップ80、第2電子密度マップ82、及び第3電子密度マップ84のうちの少なくとも1つを表示するように構成される。
ユーザインタフェースシステム40は、ユーザが第1の電子密度マップ80上の一つ又はそれより多くのアーチファクト領域90を識別する識別データ44を入力することを可能にする入力システム42を備えるように構成される。例えば、ユーザは、ペン型デバイスを使用して、識別データ44として、一つ又はそれより多くのアーチファクト領域90のアーチファクト輪郭をディスプレイ上に描く。他の入力システム42は、例えば、シード配置及び領域成長、ペイントボール、領域を囲むライブワイヤ、領域内にボールを配置することを含むことができる。
ユーザは、アーチファクトであると考えられる第1の電子密度マップ80上の領域を識別し、画定する専門的な経験又はアーチファクト領域を識別し、画定するように画像処理システムへの入力又は簡便性のために提供される閾値に頼る。
処理システム30は、アーチファクト領域であるとみなされる領域を識別して画定するように構成することができる。例えば、一つ又はそれより多くのアーチファクト領域90は、体内の皮質骨セグメントにおいてサイズ閾値によって検出されてもよい。処理システム30は、画像内の検出される組織タイプの間の関係を考慮に入れて、この関係をアーチファクトを示すための有効性指標として使用することができる。そのような関係の一例は、骨における皮質骨であり得る。最小及び最大閾値を使用する方法は、一つ又はそれより多くのアーチファクト領域を境界付けするために使用されることができる。
処理システム30は、第1の電子密度マップ80内の識別される一つ又はそれより多くのアーチファクト領域90を使用し、識別される一つ又はそれより多くのアーチファクト領域90内の第2の電子密度マップ82によって第1の電子密度マップ80を置き換え、第3の電子密度マップ84をもたらすように構成される。
したがって、第3の電子密度マップ84は、アーチファクトフリーであるか、又は識別されるアーチファクト領域90内の検査ボリュームの実際の密度値に対する少なくともずれを有し、それはアーチファクト値91からもたらされる実際の密度値に対するずれよりも小さくなる。言い換えれば、処理システム30は、未知のエラーを有する識別されるアーチファクト領域90を、既知で許容可能なエラーを示す領域で置き換える。第3の電子密度マップ84は、例えば、放射線治療計画に使用されることができる代わりの電子密度マップである。
放射線治療計画システム45は、第1の電子密度マップ80又は第3の電子密度マップ84のうちの1つを使用して、撮像ボリュームのための線量分布を計算するように構成される。特に、放射線治療計画システム45は、第1の電子密度マップ80を用いる撮像ボリュームの第一の線量分布6と、第3の電子密度マップ84を用いる撮像ボリュームの第2の線量分布7とを計算するように構成される。
第1の線量分布6は、図4Aに示されるように、第1の電子密度マップ80上に表示され重ね合わされてもよい。第2の線量分布7は、第3の電子密度マップ84(図4B)上に表示され、重ね合わされてもよい。
第1の電子密度マップ80で計算される第1の線量分布35又は第3の電子密度マップ84上で計算される第2の線量分布37を示す画像は、インタラクティブディスプレイ40又は二次ディスプレイ44によって表示されることができる。
第1の電子密度マップ80で計算される第1の線量分布35又は第3の電子密度マップ84上で計算される第2の線量分布37を示す画像は、インタラクティブディスプレイ40の別々の表示ゾーン、又は二次ディスプレイ44によって同時に表示されてもよい。
当業者であれば、本発明は、特定の領域、例えばありそうもない領域としてユーザによってマークされる領域におけるMR-CAT画像を、一般的に真のEDマップに対してより大きなずれを示すが、治療線量計算にとって許容可能な結果を与えることが知られている他のEDマップからのパッチによって置換することを提案する。換言すれば、未知のエラーを有する領域は、線量分布における既知で許容エラーを示す領域で置き換えられる。
本システムの他の利点は、ユーザが、識別される一つ又はそれより多くのアーチファクト領域を変化させることによって、線量分布に対するアーチファクトの潜在的影響を評価することができることにある。
処理システム30及び放射線治療計画システム45は、一つ又はそれより多くのコンピュータプロセッサを有する何れかのタイプの計算デバイスを使用して実装されることができる。
図4は、本開示の一態様による代わりの電子密度マップ又はMR-CAT画像を計算する方法を示す。この方法は図1のシステムを参照して説明される。
第1のステップS1において、処理システム30は、複数の画像データを使用して、第1の電子密度マップ80を計算する。複数の画像データは、磁気共鳴イメージングシステムから得られる収集データであってもよい。第1の電子密度マップ80は、可能な最良の詳細レベルを示す元のMR-CAT画像に対応する。第1の電子密度マップ80は、妥当ではないと考えられ、したがってアーチファクト領域90であると考えられる一つ又はそれより多くのゾーンを有する。
アーチファクト領域は、ユーザインターフェースシステム40上で第1の電子密度マップ80を視覚化するユーザによって識別されマークアップされてもよい。例えば、ユーザは、第1の電子密度マップ80上に一つ又はそれより多くのアーチファクト領域90を画定する輪郭をペンによって入力することができる。識別及び画定ステップは、適切な画像処理システム及びソフトウェアによって自動的に実行されることもできる。処理システム30は、例えば、アーチファクトを示すための有効性指標として画像内の検出組織タイプ間の関係を使用することができる。そのような関係の一例は、骨の皮質骨分画になり得る。
処理システム30は、磁気共鳴測定によって収集されるイメージングデータに基づいてMR-CAT画像を計算するために当業者に知られている何れかの計算方法を使用することができる。
第2のステップS2において、処理システムは、第1の電子密度マップ80よりも少ない細部を有する、第2の電子密度マップ82を計算する。換言すれば、第2の電子密度マップ82は、第1の電子密度マップ80の簡略化バージョンである。例は図2Bに示されており、第2の電子密度マップ82は、簡略化される第1の電子密度マップ80に対応し、本体輪郭のみがセグメント化され、均一な値で満たされる。例えば、水等価密度値が用いられることができる。
このステップS2は、画像アーチファクトに対してよりロバストである一方、真の密度マップからのずれは、元のMR-CAT画像と比較して一般に大きくなる。
ステップS1及びS2は、並列に又は順次に実行されてもよい。第2の電子密度マップ82は、複数のイメージングデータを用いて計算されることができる。本発明の代わりの態様では、第2の電子密度マップ82は、第1の電子密度マップ80から導き出されることができる。例えば、処理システム30は、第1の電子密度マップの輪郭を導出し、第1の電子密度マップを他の密度値で満たして、第2の電子密度マップ82を提供する。第2の電子密度マップ82は、第1の電子密度マップ80に最大閾値を設定し、この最大閾値を超える何れかの密度値を他の密度で設定することによって生成される。好ましい他の密度値は、水の密度値である。
第3のステップS3において、処理システム30は、第1の電子密度マップ80及び第2の電子密度マップ82を用いて代わりの電子密度マップ84を計算する。処理システム30は、一つ又はそれより多くのアーチファクト90において、第2の電子密度マップ82によって第1の電子密度マップ80を置換して、代わりの電子密度マップ84を得る。
図5は、本開示の一態様における線量分布を計算する方法のワークフローを示す。この方法は、図1に記載されるシステムを参照して説明される。
第1のステップS21において、第1の電子密度マップ80が、磁気共鳴技術又は同様の技術によって得られるイメージングデータを用いて計算される。第1の線量分布6は、第1の電子密度マップを用いて計算され、第1の電子密度マップと共に表示される。
第2のステップS22では、代わりの密度マップ84が、図4を参照して説明される方法を使用して計算される。代わりの密度マップ84を使用して、第2の線量分布7が計算される。第2の線量分布7は、代わりの電子密度マップ84と共に表示される。
第2の線量分布7は、第1の電子密度マップ80のアーチファクト領域90に起因する未知のエラーで生成される第1の線量分布6の代わりに、許容され得るエラーで生成されることが明らかであろう。本開示は、ユーザが、線量分布に対するアーチファクトの潜在的影響を評価することを可能にする。
本発明は、MRベースの密度マップ、例えばPinnacle又はMRベースのRT計画の潜在的なソリューションを使用するRT計画システムに適用可能である。更に、PET / MRシステム上のPET再構成にも適用されることができる。
本発明は、勿論説明される又は図示される実施形態に限定されるものではなく、前述の説明及び図面に照らして見られるように概して添付の特許請求の範囲内に入る何れの実施形態にも及ぶ。本発明の特定の特徴は、例示される実施形態のうちの1つのみに関して上述されているが、何れかの所与の特定の用途に所望され、有利であるように、他の実施形態の一つ又はそれより多くの特徴と組み合わされることができる。本発明の上記説明から、当業者は、改良、変更及び修正を認識するであろう。そのような改良、変更、及び修正は、添付の特許請求の範囲によってカバーされることが意図される。特許請求の範囲内のいかなる参照符号も、本発明の範囲を限定するものではない。用語「有する(comprising)」は、他の要素又はステップを排除するものではなく、用語「a」又は「an」は複数を除外しないことを理解される。
1:システム
6:第一の線量分布
7:第二の線量分布
10:イメージングシステム
20:複数の画像データ
30:処理システム
40:ユーザーインターフェースシステム
42:ユーザー入力部
44:二次ディスプレイ
80:第一の電子密度マップ
82:第二の電子密度マップ
84:第三の電子密度マップ
90:アーチファクトゾーン
91:密度値アーチファクト
92:密度値の近似
93:置換密度値
95:輪郭

Claims (11)

  1. 検査ボリュームの代わりの電子密度マップを計算するための処理システムであって、
    -複数の磁気共鳴イメージングデータを使用して第1の電子密度マップを計算し、
    -第2の電子密度マップを計算し、前記第2の電子密度マップは前記第1の電子密度マップの簡略化バージョンであり、前記第2の電子密度マップは、前記第1の電子密度マップの輪郭に対応し、均一な値で満たされ、
    -前記第1の電子密度マップと前記第2の電子密度マップとを用いて前記代わりの電子密度マップを計算する
    ように構成され
    前記処理システムは、前記代わりの電子密度マップを得るために、一つ又はそれより多くのアーチファクト領域において、前記第1の電子密度マップを前記第2の電子密度マップによって置き換えるように構成される、
    処理システム。
  2. ユーザインターフェースシステムを有し、前記ユーザインターフェースシステムは、前記第1の電子密度マップを表示するように構成され、前記ユーザインターフェースシステムは、ユーザが前記第1の電子密度マップ上で前記一つ又はそれより多くのアーチファクト領域を識別し、画定するユーザデータを入力することを可能にするユーザ入力システムを有する、請求項1に記載の処理システム。
  3. 前記均一な値は等価水密度である、請求項1乃至2の何れか一項に記載の処理システム。
  4. 線量分布を計算するための計画システムであって、
    -請求項1乃至3の何れか一項に記載の処理システムであって、複数のイメージングデータ、第2の電子密度マップ、及び代わりの電子密度マップを使用して第1の電子密度マップを計算するように構成され、前記代わりの電子密度マップは前記第1及び第2の電子密度マップを使用して計算される、処理システムと、
    -前記第1の電子密度マップを使用して第1の線量分布を計算し、前記代わりの電子密度マップを使用して第2の線量分布を再計算するように構成される放射線治療計画システムと
    を有する、計画システム。
  5. 検査ボリュームの代わりの電子密度マップを計算する方法であって、
    -複数の磁気共鳴イメージングデータを用いて第1の電子密度マップを計算するステップと、
    -第2の電子密度マップを計算するステップであって、前記第2の電子密度マップは前記第1の電子密度マップの簡略化バージョンであり、前記第2の電子密度マップは、前記第1の電子密度マップの輪郭に対応し、均一な値で満たされる、ステップと、
    -前記第1の電子密度マップ及び前記第2の電子密度マップを用いて前記代わりの電子密度マップを計算するステップと、
    一つ又はそれより多くの識別されるアーチファクト領域において、前記第1の電子密度マップを前記第2の電子密度マップで置き換えて、前記代わりの電子密度マップを得るステップと
    を有する、方法。
  6. 前記第1の電子密度マップを表示するステップと、前記一つ又はそれより多くのアーチファクト領域を識別するユーザ入力データを得るステップと
    を有する、請求項5に記載の方法。
  7. 前記第2の電子密度マップに、水等価密度値を割り当てるステップ
    を有する、請求項5乃至6の何れか一項に記載の方法。
  8. 線量分布を計算する方法であって、
    -第1の電子密度マップを計算するステップと、
    -前記第1の電子密度マップを使用して第1の線量分布を計算するステップと、
    -請求項5乃至7の何れか一項に記載の方法及び/又は請求項1乃至4の何れか一項に記載のシステムを用いて代わりの密度マップを計算するステップと、
    -前記代わりの電子密度マップを使用して第2の線量分布を計算するステップと
    を有する、方法。
  9. プロセッサが請求項5乃至7の何れか一項に記載の方法のステップを実行することを可能にする命令を、コンピュータに実行させるために記憶される制御論理を有する、コンピュータプログラム。
  10. コンピュータプログラムであって、
    -複数の磁気共鳴イメージングデータを使用して第1の電子密度マップを計算し、
    -前記第1の電子密度マップの簡略化バージョンであり、前記第1の電子密度マップの輪郭に対応し、均一な値で満たされる第2の電子密度マップを計算し、
    -前記第1の電子密度マップと前記第2の電子密度マップとを用いて、前記代わりの電子密度マップを計算し、
    一つ又はそれより多くの識別されるアーチファクト領域において、前記第1の電子密度マップを前記第2の電子密度マップで置き換えて、前記代わりの電子密度マップを得る
    ように構成されるコンピュータ手段を有する、コンピュータプログラム。
  11. プロセッサが請求項8に記載の方法のステップを実行することを可能にする命令を、コンピュータに実行させるために記憶される制御論理を有する、コンピュータプログラム。
JP2017548434A 2015-03-20 2016-03-18 Mrcat画像の不確定領域に対するフォールバックソリューション Active JP6802177B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP15160158.0 2015-03-20
EP15160158 2015-03-20
PCT/EP2016/056072 WO2016150888A1 (en) 2015-03-20 2016-03-18 Fall-back solution for uncertain regions in mrcat images

Publications (3)

Publication Number Publication Date
JP2018510705A JP2018510705A (ja) 2018-04-19
JP2018510705A5 JP2018510705A5 (ja) 2020-06-25
JP6802177B2 true JP6802177B2 (ja) 2020-12-16

Family

ID=52824019

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017548434A Active JP6802177B2 (ja) 2015-03-20 2016-03-18 Mrcat画像の不確定領域に対するフォールバックソリューション

Country Status (6)

Country Link
US (1) US10870017B2 (ja)
EP (1) EP3270779B1 (ja)
JP (1) JP6802177B2 (ja)
CN (1) CN107427690B (ja)
RU (1) RU2721932C2 (ja)
WO (1) WO2016150888A1 (ja)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1057371A (ja) * 1996-08-27 1998-03-03 Ge Yokogawa Medical Syst Ltd 医用装置および標的領域設定方法
US8718232B2 (en) * 2009-04-30 2014-05-06 Siemens Aktiengesellschaft Reconstruction of an emission tomography image
RU2571374C2 (ru) * 2009-08-31 2015-12-20 Конинклейке Филипс Электроникс, Н.В. Интерактивный компьютеризованный редактор для компенсаторов, используемых в планировании лучевой терапии
EP2724172A2 (en) * 2011-06-27 2014-04-30 Koninklijke Philips N.V. Bone mri using an ultrashort echo time pulse sequence with fid and multiple gradient echo acquisition and water-fat separation processing
US9245208B2 (en) 2011-08-11 2016-01-26 The Regents Of The University Of Michigan Patient modeling from multispectral input image volumes
DE102012211146A1 (de) 2012-06-28 2014-01-02 Siemens Aktiengesellschaft Verfahren und Röntgensystem zur Erzeugung einer Phasenkontrastdarstellung
JP5966112B1 (ja) * 2013-04-02 2016-08-10 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 磁気共鳴画像法を使用した骨組織の検出
EP3039441A1 (en) * 2013-08-30 2016-07-06 Koninklijke Philips N.V. Dixon magnetic resonance imaging
JP6378778B2 (ja) * 2014-03-28 2018-08-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 組織種類分離の助けを借りて磁気共鳴画像に基づき1つ以上のコンピュータ断層撮影画像を生成する方法及び装置
JP6927880B2 (ja) 2015-02-17 2021-09-01 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 医療撮像検出器
DE102016218359B3 (de) * 2016-09-23 2017-11-09 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Korrektur einer synthetischen Elektronendichtekarte

Also Published As

Publication number Publication date
RU2017135135A (ru) 2019-04-05
CN107427690A (zh) 2017-12-01
RU2721932C2 (ru) 2020-05-25
US10870017B2 (en) 2020-12-22
WO2016150888A1 (en) 2016-09-29
RU2017135135A3 (ja) 2019-08-01
CN107427690B (zh) 2020-09-01
US20180071550A1 (en) 2018-03-15
EP3270779B1 (en) 2019-11-06
EP3270779A1 (en) 2018-01-24
JP2018510705A (ja) 2018-04-19

Similar Documents

Publication Publication Date Title
Wagenknecht et al. MRI for attenuation correction in PET: methods and challenges
CN106133790B (zh) 用于在组织种类分离的帮助下基于磁共振图像生成一幅或多幅计算机断层摄影图像的方法和设备
US10557904B2 (en) Detection of bone tissue using magnetic resonance imaging
JP6506495B2 (ja) 磁気共鳴撮像法を使用した陽電子放出断層撮像法における減衰補正
JP7019568B2 (ja) 白質脳病変の分析のための医療機器
US10114098B2 (en) Method estimating a pseudo Hounsfield Unit value
US20110007959A1 (en) Ct surrogate by auto-segmentation of magnetic resonance images
RU2608975C2 (ru) Способ магниторезонансной томографии (mri) для назначения индивидуальным пикселам или вокселам специфических для ткани значений ослабления позитронно-эмиссионной томографии (рет)
US10052032B2 (en) Stenosis therapy planning
Abadi et al. Patient‐specific quantification of image quality: an automated technique for measuring the distribution of organ Hounsfield units in clinical chest CT images
Malkus et al. A method to extract image noise level from patient images in CT
US20180253838A1 (en) Systems and methods for medical imaging of patients with medical implants for use in revision surgery planning
US8831323B2 (en) Method and apparatus for measuring activity of a tracer
US9760991B2 (en) System and method for image intensity bias estimation and tissue segmentation
US20180085080A1 (en) Method and apparatus for correction of a synthetic electron density map
US20110148861A1 (en) Pet data processing system, an arrangement, a method and a computer program product for determining a distribution of a tracer uptake
Boydev et al. Zero echo time MRI-only treatment planning for radiation therapy of brain tumors after resection
CN108475419A (zh) 用于计算机断层摄影的数据处理的方法
Solomon et al. New image analysis technique for quantitative longitudinal assessment of lung pathology on CT in infected rhesus macaques
JP6802177B2 (ja) Mrcat画像の不確定領域に対するフォールバックソリューション
Chen et al. Low dose cone-beam computed tomography reconstruction via hybrid prior contour based total variation regularization (hybrid-PCTV)
JP2017072482A (ja) 散乱補正方法
JP2010507438A (ja) セグメント化の向上
Anderla et al. Suppression of metal artefacts in CT using virtual singorams and corresponding MR images
Debus Medical Image Processing

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190306

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191010

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20191010

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200123

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20200210

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200218

A524 Written submission of copy of amendment under article 19 pct

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20200515

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200618

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200916

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201126

R150 Certificate of patent or registration of utility model

Ref document number: 6802177

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250