JP7284876B2 - 乳房磁気共鳴ctにおける勧告領域の特定 - Google Patents

乳房磁気共鳴ctにおける勧告領域の特定 Download PDF

Info

Publication number
JP7284876B2
JP7284876B2 JP2022555738A JP2022555738A JP7284876B2 JP 7284876 B2 JP7284876 B2 JP 7284876B2 JP 2022555738 A JP2022555738 A JP 2022555738A JP 2022555738 A JP2022555738 A JP 2022555738A JP 7284876 B2 JP7284876 B2 JP 7284876B2
Authority
JP
Japan
Prior art keywords
region
regions
magnetic resonance
tissue
weighted magnetic
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
JP2022555738A
Other languages
English (en)
Other versions
JP2023509237A (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 JP2023509237A publication Critical patent/JP2023509237A/ja
Application granted granted Critical
Publication of JP7284876B2 publication Critical patent/JP7284876B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/5607Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reducing the NMR signal of a particular spin species, e.g. of a chemical species for fat suppression, or of a moving spin species for black-blood imaging
    • 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/43Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4312Breast evaluation or disorder diagnosis
    • 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/50NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341Diffusion imaging
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56554Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/174Segmentation; Edge detection involving the use of two or more images
    • 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/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30068Mammography; Breast

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Medical Informatics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Vascular Medicine (AREA)
  • Quality & Reliability (AREA)
  • Reproductive Health (AREA)
  • Gynecology & Obstetrics (AREA)
  • Artificial Intelligence (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は磁気共鳴CTに関し、特に、磁気共鳴CTを用いた乳房のスクリーニングに関する。
大きな静磁場は患者の身体に画像を生成するためのプロシージャの一部として原子の核スピンを位置合わせさせるために、磁気共鳴撮像(MRI)スキャナによって使用される。この大きな静磁場は、B0場又は主磁場と呼ばれる。対象の様々な量又は特性は、MRIを使用して空間的に測定することができる。
国際特許出願WO2012070951A1は、動的造影磁気共鳴画像(DCE-MRI)及び動的感度増強磁気共鳴画像(DSC)からの測定データの位置合わせ及び比較を通して、軟質組織、特に母乳組織などの組織における、ex vivoと悪性腫瘍を区別するためのプロセスを開示している。プロセスはMR造影物質の非経口投与中にインターリーブモードで2つの動的MRIパルスシーケンスを最初に実行するように調整され、ここで、前記パルスシーケンスの1つは高い空間分解能のために最適化され、他のパルスシーケンスは高時間分解能のために調整され、高時間溶解シーケンスは薬物動態バイオマーカーデータ、記述的DCEバイオマーカー及び記述的DSCバイオマーカーなどのいくつかの異なるバイオマーカーデータを生成するために、DCE及びDSCの両方に対して感受性である二重エコー収集をさらに含み、その後、対応する良性腫瘍及び悪性腫瘍からの対応するデータとそれぞれ正規化及び比較する。
各腫瘍体積から、(関連するバイオマーカーに属する)異常な動力学的特性の大部分を同定する95パーセンタイルについて収集されるバイオマーカーがあることが好ましい。
本発明は、独立請求項における医用システム、コンピュータプログラム、及び方法を提供する。実施形態は、従属請求項において与えられる。
異常組織についての乳房組織のスクリーニングは、典型的にはX線に基づくマンモグラフィを用いて行われる。実施形態は、磁気共鳴画像を使用して乳房組織内の勧告領域を位置特定する改善される手段を提供することができる。これを達成するために、2つのタイプの画像が使用される。エコー平面拡散強調磁気共鳴画像は、高い拡散率領域を位置付けるために使用される。高拡散率領域は例えば、所定の拡散率閾値を超える拡散率を有する領域のラベルであってもよい。この技術は、異常な組織領域を同定するのに非常に有効である。しかしながら、エコー平面拡散強調磁気共鳴画像は、空間歪みを含むことがある。異常な組織を同定することができるが、その真の位置は不明であり得る。
高拡散率領域は例えば、交互に、識別される拡散率領域と呼ばれてもよい。
適切なローカライゼーションを提供するために、脂肪抑制T2強調磁気共鳴画像が使用される。脂肪抑制T2強調磁気共鳴画像は、組織領域を識別するためにセグメント化される。脂肪又は脂質抑制が使用されるので、これらのセグメント化は、非脂肪組織領域又は非脂質組織領域であると考えられ得る。次いで、得られた組織領域は、乳房内の正常組織構造ならびに任意の異常構造のためのセグメント化を含む。次いで、画像処理モジュールを使用して、組織領域の一部を勧告領域として識別する。
撮像処理モジュールは、高拡散率領域と組織領域との間の空間的、体積的、又は形状に基づく関係を使用して、高拡散率領域を引き起こす原因となる組織領域の空間的位置を識別するように構成され得る。
一態様では、本発明が機械実行可能命令を記憶するメモリを備える医用システムを提供する。医用システムは、医用システムを制御するように構成される計算システムをさらに備える。医用システムは、異なる実施例において異なる形態をとることができる。一例では、医用システムはワークステーションとすることができる。別の例では、医用システムが磁気共鳴CTのための画像処理及び再構成サービスを提供するリモートサーバなどであり得る。別の例では、医用システムが磁気共鳴撮像システムなどの医用撮像システムを含むシステムである。
機械実行可能命令の実行はコンピュータシステムに、乳房組織を記述する関心領域のエコー平面拡散強調磁気共鳴画像を受信させる。
機械実行可能命令の実行はさらに、計算システムに、関心領域を記述する脂肪抑制T2強調磁気共鳴画像を受信させる。関心領域が同一又は同等であるように、エコー平面拡散強調磁気共鳴画像と脂肪抑制T2強調磁気共鳴画像との間に位置合わせがあってもよい。機械実行可能命令の実行はさらに、計算システムに、エコー平面拡散強調磁気共鳴画像をセグメント化させて、高拡散率領域を識別させる。
高拡散率領域を識別する様々な手段があり得るが、一方法はエコー平面拡散強調磁気共鳴画像をしきい値にすることである。高拡散率領域の識別は、悪性腫瘍であるなどの病理を有する組織を識別することができる。エコー平面拡散強調磁気共鳴画像はそのような病状又は組織タイプを検出するのに良好であり得るが、エコー平面拡散強調磁気共鳴画像の歪みがしばしば存在する。したがって、そのような組織タイプを識別することは可能であり得るが、正確な位置は不明であり得る。
機械実行可能命令の実行はさらに、計算システムに、脂肪抑制T2強調磁気共鳴画像をセグメント化して、組織領域を識別させる。これらの組織領域は例えば、非脂肪組織又は非脂質組織であってもよい。画像は脂肪抑制されているので、識別される組織領域は、脂肪又は脂質組織を除外することができる。機械実行可能命令の実行はさらに、計算システムに、高拡散率領域及び組織領域を画像処理モジュールに入力することによって、組織領域の一部を勧告領域として識別させる。
勧告領域という用語の使用は、さらなるレビュー又は検査のために医師又は他の医療提供者に提供又は強調され得る領域を示す。画像処理モジュールは例えば、組織領域を高拡散率領域と相関させるために使用されてもよい。脂肪抑制T2強調磁気共鳴画像において、同定され得る組織の多くの部分が存在する。例えば、乳房内の腺組織、ならびに線維性、支持性又は結合性組織が存在し得る。さらに、靭帯及び瘢痕組織があってもよい。これは、存在し得る任意の癌又は腫瘍に加えてである。
磁気共鳴画像では、これらの組織型を区別することはしばしば実用的ではない。一方、エコー平面拡散強調磁気共鳴画像は、高い拡散率領域を有する乳房の領域に関する情報を提供する。すなわち、拡散率領域は、所定の閾値を超える拡散を有する。高拡散率領域を組織領域と関連付けることは、勧告領域の識別を可能にし、医師を支援することができる。機械実行可能命令の実行はさらに、計算システムに、脂肪抑制T2強調磁気共鳴画像のセグメンテーションとして勧告領域を提供させる。したがって、T2強調磁気共鳴画像は、マーキング又は強調表示される勧告領域とともに表示され得る。これは、例えば、医師がさらなる医学的調査又はプロシージャを計画するのを助けることができる。
上述のことは、2次元及び3次元の両方の画像において実行され得ることが言及され得る。脂肪抑制T2強調磁気共鳴画像及びエコー平面拡散強調磁気共鳴画像は2次元画像、2次元スライスの集まり、又は完全な3次元データセットであり得る。
別の実施形態では、エコー平面拡散強調磁気共鳴画像が高b値エコー平面拡散強調磁気共鳴画像である。高b値エコー平面拡散強調磁気共鳴画像の使用は、腫瘍組織又は癌組織を含む可能性がある高い拡散率領域を識別するのに有用であり得るので、有益であり得る。しかしながら、高b値エコー平面拡散強調磁気共鳴画像の使用は、画像のより大きな歪みを引き起こす可能性がある。脂肪抑制T2強調磁気共鳴画像からの組織領域を用いた勧告領域の同定は、位置及び高拡散率領域の両方の識別を可能にする。ラベル「高b値」の使用は、画像の幾何学的歪みを引き起こし得るb値を使用して画像が取得されることを示すものとして解釈され得るラベルである。
別の実施形態では、機械実行可能命令の実行がさらに、計算システムに、関心領域を記述する導電率マップを受信させる。機械実行可能命令の実行はさらに、計算システムに、導電率マップを使用して、勧告領域の各々の導電率測定を計算させる。導電率測定は、勧告領域として識別されるボクセルについて計算される統計値又は測定であってもよい。これは、平均導電率ならびに導電率の最小値及び最大値及び平均値などの様々な統計量を含むことができる。
別の実施形態では、勧告領域が導電率マップの計算中に境界を定義するために使用される。
機械実行可能命令の実行はさらに、計算システムに、導電率測定に従って、勧告領域の各々の分類を割り当てさせる。導電率測定値は、ナトリウムなどの様々なイオンタイプに非常に依存し得る。勧告領域が特定の高い拡散レート値を超える高い拡散率、ならびに所定の閾値を超える平均導電率又は他の導電率測定の両方を有する場合、それは、特定の勧告領域が癌性であるか、又は医師によってさらに調査されるべきであることを示し得る。
別の実施形態では、機械実行可能命令の実行がさらに、計算システムに、脂肪抑制T2強調磁気共鳴画像上のオーバーレイとして勧告領域をレンダリングさせる。
別の実施形態では、機械実行可能命令の実行がさらに、計算システムに、導電率測定に従って、又は導電率測定を使用して、勧告領域の各々の分類をレンダリングさせる。
例えば、所定の導電率測定閾値を超える導電率測定を有する勧告領域は、ディスプレイ又はレンダリング上で強調することができる。別の例では、レンダリングが勧告領域が調査されることが推奨される順序をランク付け又は割り当てることができる。ユーザインターフェースはレンダリング上にランクを表示することができ、又はユーザインターフェースは、オペレータがランクに従って勧告領域を順次強調することができる制御を有することができる。
別の実施形態では、医用システムが撮像ゾーンからk空間データを取得するように構成される磁気共鳴撮像システムをさらに備える。メモリは、第1のパルスシーケンスコマンド及び第2のパルスシーケンスコマンドをさらに含む。第1のパルスシーケンスコマンドは、エコー平面拡散強調磁気共鳴撮像プロトコルに従って第1のk空間データを取得するように構成される。第2のパルスシーケンスコマンドは、T2強調磁気共鳴撮像プロトコルに従って第2のk空間データを取得するように構成される。第2のパルスシーケンスコマンドは、脂肪又は脂質抑制である。
機械実行可能命令の実行はさらに、計算システムに、第1のパルスシーケンスコマンドで磁気共鳴撮像システムを制御することによって、乳房組織を含む関心領域から第1のk空間データを取得させる。機械実行可能命令の実行はさらに、計算システムに、第2のパルスシーケンスコマンドを用いて磁気共鳴撮像システムを制御することによって、乳房組織を含む関心領域から第2のk空間データを取得させる。第1及び第2のパルスシーケンスコマンドは、両方の取得のための関心領域が同じ空間体積からであるように構成され得る。機械実行可能命令の実行はさらに、計算システムに、第1のk空間データからエコー平面拡散加重磁気共鳴画像を再構成させる。機械実行可能命令の実行はさらに、計算システムに、第2のk空間データから脂肪抑制T2強調磁気共鳴画像を再構成させる。
別の実施形態では、第1のパルスシーケンスコマンドが0s/mm2乃至3000s/mm2のb値を有する第1のk空間データを取得するように構成される。
別の実施形態では、第1のパルスシーケンスコマンドが800s/mm2と1600s/mm2 との間のb値を有する第1のk空間データを取得するように構成される。この範囲のいくつかの例では、b値が高b値であると見なされ得る。この実施形態は腫瘍などの組織型のDWI検出がこれらのb値が使用されるときに良好に機能するが、有意な幾何学的歪みがあるので、有益であり得る。
別の実施形態では、第1のパルスシーケンスコマンドが1200s/mm2乃至1400s/mm2のb値を有する第1のk空間データを取得するように構成される。この範囲のいくつかの例では、b値が高b値であると見なされ得る。この実施形態は腫瘍などの組織型のDWI検出がこれらのb値が使用されるときに非常に効果的に機能するが、有意な幾何学的歪みがあるので、有益であり得る。
別の実施形態では、第2のパルスシーケンスコマンドが電気特性トモグラフィ磁気共鳴CTプロトコルに従う。電気特性トモグラフィ磁気共鳴CTプロトコルは例えば、B1磁場の歪みを測定するために使用される磁気共鳴CTプロトコルであってもよい。
多くの脂肪抑制パルスは、電気的特性トモグラフィに適合する。様々なアプローチがあるが、いくつかの主要なアプローチが挙げられる。第1に、脂肪抑制RFパルスの周波数を介して脂肪プロトンを選択的に飽和させるRFパルスを使用する脂肪抑制があり得る。別のアプローチは、化学シフト撮像に類似するDixon型技術を使用することである。異なるエコー時間を有する2つの画像があり、これらは、脂肪のみを得るために取得され、処理され、水のみの画像を得る。この例の水のみの画像は、脂肪抑制T2強調磁気共鳴画像に使用することができる。別のアプローチは脂肪が典型的には水よりも短いT1値を有することを利用する反転回収方法である。スピンエコーパルスシーケンスで使用される上述の3つの方法は全て、EPT計算に使用することができる位相測定値をもたらす。
別の実施形態では、第2のパルスシーケンスコマンドはスピンエコーパルスシーケンスコマンドである。
別の実施形態において、第2のパルスシーケンスコマンドは、超短エコー時間パルスシーケンスコマンドである。
別の実施形態では、第2のパルスシーケンスコマンドがゼロエコー時間パルスシーケンスコマンドである。
別の実施形態では、第2のパルスシーケンスコマンドがマルチエコー勾配エコーパルスシーケンスコマンドである。
別の実施形態では、第2のパルスシーケンスコマンドが平衡勾配エコーパルスシーケンスコマンドである。
別の実施形態では、第2のパルスシーケンスコマンドが定常状態プレセッションパルスシーケンスコマンドである。上述のタイプのパルスシーケンスコマンドの全ては、上述のように脂肪抑制を実行するだけでなく、電気的特性トモグラフィに使用され得る位相の測定に適合する。
別の実施形態では、機械実行可能命令の実行はさらに、計算システムに、第2のk空間データから関心領域を記述する導電率マップを再構成させる。この実施形態は、空間的に依存する導電率が様々なイオン濃度に依存することが多いため、有益であり得る。これは、高い拡散率領域である領域、又は所定の基準を超える拡散率を有する領域を識別することに関連して、組織を分類する際に有用であり得る。例えば、分類は、プライマリケア医又は他の医療専門家によって調査するために勧告領域を調査又は提案すべき順序をランク付けするために使用され得る。
別の実施形態では、高拡散率領域がエコー平面拡散強調磁気共鳴画像を閾値処理することによって識別される。この例では、高い拡散率領域がその領域を高い拡散率を有すると定義する閾値又は所定の拡散率を提供することによって定義される。
別の実施形態では、画像処理モジュールが勧告領域をアルゴリズム的に識別するように構成される。一例では組織領域の1つが、所定のオーバラップを超える高拡散率領域の何れかのオーバラップを有する場合、2つの組織領域の1つを勧告領域の1つとして識別することによって、アルゴリズム的に識別される。
別の例では画像処理モジュールが組織領域の1つが、所定の距離未満の他の高拡散率領域の何れかとの中心距離を有する場合、勧告領域の1つとして組織領域の1つを識別する。オーバラップはなくてもよいが、中心が十分に近い場合、それらは同じ領域であると分類され得る。中心は例えば、重心であってもよい。別の例では画像処理モジュールが組織領域の1つが、所定の体積差内の高拡散率領域の何れかに一致する体積を有する場合、勧告領域の1つとして組織領域の1つを識別する。エコー平面拡散強調磁気共鳴画像は、高拡散率領域の形状及び位置の歪みを引き起こし得る。体積を比較することは、高拡散率領域を適切な組織領域に割り当てる1つの方法であり得る。
別の例における画像処理モジュールは組織領域の1つが、所定の歪み内の高拡散率領域の何れかに一致する形状を有する場合、組織領域の1つを勧告領域の1つとして識別することができる。
別の実施形態では、画像処理モジュールが組織領域及び高拡散率領域の入力に応答して、組織領域の一部を勧告領域としてラベル付けするように構成される、訓練されるニューラルネットワークである。ニューラルネットワークは例えば、最初に、人間によって分類される多数の組織領域及び高拡散率領域を有することによって訓練され得る。これは、トレーニングデータを作成するために使用され得る。例えば、次いで、ニューラルネットワークは、深層学習アルゴリズムを使用して訓練され得る。
別の実施形態では、組織領域が非脂肪組織領域又は非脂質組織領域である。
別の態様では、本発明が医用システムを制御する計算システムによる実行のための機械実行可能命令を含むコンピュータプログラムを提供する。例えば、コンピュータプログラムは、非一時的記憶媒体に記憶されるコンピュータプログラム製品であってもよい。機械実行可能命令の実行はコンピュータシステムに、乳房組織を記述する関心領域のエコー平面拡散強調磁気共鳴画像を受信させる。機械実行可能命令の実行はさらに、計算システムに、関心領域を記述する脂肪抑制T2強調磁気共鳴画像を受信させる。機械実行可能命令の実行はさらに、計算システムに、エコー平面拡散強調磁気共鳴画像をセグメント化させて、高拡散率領域を識別させる。機械実行可能命令の実行はさらに、計算システムに、脂肪抑制T2強調磁気共鳴画像をセグメント化して、組織領域を識別させる。機械実行可能命令の実行はさらに、計算システムに、高拡散率領域及び組織領域を画像処理モジュールに入力することによって、組織領域の一部を勧告領域として識別させる。機械実行可能命令の実行はさらに、計算システムに、脂肪抑制T2強調磁気共鳴画像のセグメンテーションとして勧告領域を提供させる。脂肪抑制T2強調磁気共鳴画像は例えば、セグメンテーションをオーバーレイとして表示することができる。
別の態様では、本方法が医療又は磁気共鳴CTの方法を提供する。この方法は、乳房組織を記述する関心領域のエコー平面拡散強調磁気共鳴画像を受信することを含む。方法は、関心領域を記述する脂肪抑制T2強調磁気共鳴画像を受信することをさらに含む。本方法は高拡散率領域を識別するために、エコー平面拡散強調磁気共鳴画像をセグメント化することをさらに含む。本方法は、脂肪抑制T2強調磁気共鳴画像をセグメント化して組織領域を識別することをさらに含む。本方法は、高拡散率領域及び組織領域を画像処理モジュールに入力することによって、組織領域の一部を勧告領域として識別することをさらに含む。本方法は、脂肪抑制T2強調磁気共鳴画像のセグメント化として勧告領域を提供することをさらに含む。
本発明の前述の実施形態の1つ又は複数は組み合わされる実施形態が相互に排他的でない限り、組み合わされ得ることが理解される。
当業者によって理解されるように、本発明の態様は、装置、方法、又はコンピュータプログラム製品として具現化され得る。したがって、本発明の態様は完全にハードウェアの実施形態、完全にソフトウェアの実施形態(ファームウェア、常駐ソフトウェア、マイクロコードなどを含む)、又は本明細書ではすべて一般に「回路」、「モジュール」、又は「システム」と呼ぶことができるソフトウェア態様とハードウェア態様とを組み合わせた実施形態の形態をとることができ、さらに、本発明の態様は、その上に具現化されるコンピュータ実行可能コードを有する1つ又は複数のコンピュータ可読媒体内に具現化されるコンピュータプログラム製品の形態をとることができる。
1つ又は複数のコンピュータ可読媒体の任意の組合せが利用され得る。コンピュータ可読媒体は、コンピュータ可読信号媒体又はコンピュータ可読記憶媒体であり得る。本明細書で使用される「コンピュータ可読記憶媒体」は、コンピューティングデバイスのプロセッサ又は計算システムによって実行可能な命令を記憶することができる任意の有形記憶媒体を包含する。コンピュータ可読記憶媒体は、コンピュータ可読非一時的記憶媒体と呼ばれることがある。コンピュータ可読記憶媒体は、有形のコンピュータ可読媒体と呼ばれることもある。いくつかの実施形態では、コンピュータ可読記憶媒体はまた、コンピューティングデバイスの計算システムによってアクセスされることが可能であるデータを記憶することが可能であり得る。コンピュータ可読記憶媒体の例としてはフロッピー(登録商標)ディスク、磁気ハードディスクドライブ、ソリッドステートハードディスク、フラッシュメモリ、USBサムドライブ、ランダムアクセスメモリ(RAM)、リードオンリーメモリ(ROM)、光ディスク、光磁気ディスク、及び計算システムのレジスタファイルが挙げられるが、これらに限定されない。光ディスクの例としては、コンパクトディスク(CD)及びデジタル多用途ディスク(DVD)、例えば、CDROM、CDRW、CDR、DVDROM、DVDRW、又はDVDRディスクが挙げられる。コンピュータ可読記憶媒体という用語はまた、ネットワーク又は通信リンクを介してコンピュータ装置によってアクセスされることが可能な様々なタイプの記録媒体を指す。たとえば、データは、モデムを介して、インターネットを介して、又はローカルエリアネットワークを介して取り出され得る。コンピュータ可読媒体上に具現化されるコンピュータ実行可能コードはワイヤレス、ワイヤライン、光ファイバケーブル、RFなど、又は前述の任意の適切な組合せを含むが、これらに限定されない、任意の適切な媒体を使用して送信され得る。
コンピュータ可読信号媒体は例えば、ベースバンドにおいて、又は搬送波の一部として、その中に具現化されるコンピュータ実行可能コードを有する伝搬データ信号を含み得る。そのような伝搬信号は電磁気、光学、又はそれらの任意の適切な組合せを含むが、それらに限定されない、様々な形態の何れかをとり得る。コンピュータ可読信号媒体は、コンピュータ可読記憶媒体ではなく、命令実行システム、装置、又はデバイスによって、又はそれに関連して使用するためのプログラムを通信、伝播、又は移送することができる任意のコンピュータ可読媒体であり得る。
「コンピュータメモリ」又は「メモリ」は、コンピュータ可読記憶媒体の一例である。コンピュータメモリは、計算システムに直接アクセス可能な任意のメモリである。「コンピュータ記憶装置」又は「ストレージ」は、コンピュータ可読記憶媒体のさらなる例である。コンピュータ記憶装置は、任意の不揮発性メモリコンピュータ可読記憶媒体である。いくつかの実施形態では、コンピュータ記憶装置がコンピュータメモリであってもよく、逆もまた同様である。本明細書で使用される「計算システム」は、プログラム、機械実行可能命令、又はコンピュータ実行可能コードを実行することができる電子コンポーネントを包含する。
「計算システム」の例を含む計算システムへの言及は、場合によっては2つ以上の計算システム又は処理コアを含むものとして解釈されるべきである。計算システムは例えば、マルチコアプロセッサであってもよい。計算システムはまた、単一のコンピュータシステム内の、又は複数のコンピュータシステム間で分散される計算システムの集合を指し得る。計算システムという用語はまた、おそらく、各々がプロセッサ又は計算システムを備えるコンピューティングデバイスの集合又はネットワークを指すと解釈されるべきである。機械実行可能コード又は命令は、同じコンピューティングデバイス内にあってもよく、又は複数のコンピューティングデバイスにわたって分散されてもよい複数の計算システム又はプロセッサによって実行されてもよい。
機械実行可能命令又はコンピュータ実行可能コードは、プロセッサ又は他の計算システムに本発明の態様を実行させる命令又はプログラムを備え得る。本発明の態様のための動作を実行するためのコンピュータ実行可能コードは、Java、Smalltalk、C++などのオブジェクト指向プログラミング言語、及び「C」プログラミング言語又は同様のプログラミング言語などの従来の手続き型プログラミング言語を含む、1つ又は複数のプログラミング言語の任意の組合せで書かれ、機械実行可能命令にコンパイルされ得る。いくつかの事例では、コンピュータ実行可能コードが高級言語の形態又は事前コンパイルされる形態であり得、オンザフライで機械実行可能命令を生成するインタプリタと併せて使用され得る。他の事例では、機械実行可能命令又はコンピュータ実行可能コードがプログラマブル論理ゲートアレイのためのプログラムの形成であり得る。
コンピュータ実行可能コードは、完全にユーザのコンピュータ上で、部分的にユーザのコンピュータ上で、スタンドアロンソフトウェアパッケージとして、部分的にユーザのコンピュータ上で、部分的にリモートコンピュータ上で、又は完全にリモートコンピュータもしくはサーバ上で実行することができる。後者のシナリオでは、リモートコンピュータがローカルエリアネットワーク(LAN)又はワイドエリアネットワーク(WAN)を含む任意のタイプのネットワークを介してユーザのコンピュータに接続されてもよく、又は接続は外部コンピュータに(たとえば、インターネットサービスプロバイダを使用してインターネットを介して)行われてもよい。
本発明の態様は、本発明の実施形態による方法、装置(システム)、及びコンピュータプログラム製品のフローチャート図及び/又はブロック図を参照して説明される。フローチャート、図、及び/又はブロック図の各ブロック又はブロックの部分は、適用可能な場合、コンピュータ実行可能コードの形態のコンピュータプログラム命令によって実装され得ることを理解される。さらに、相互に排他的ではない場合、異なるフローチャート、図、及び/又はブロック図におけるブロックの組合せが組み合わされ得ることに留意される。
これらのコンピュータプログラム命令は、汎用コンピュータ、専用コンピュータ、又は他のプログラマブルデータ処理装置の計算システムに提供されて、コンピュータ又は他のプログラマブルデータ処理装置の計算システムを介して実行される命令がフローチャート及び/又はブロックダイアグラムブロック又はブロックで指定される機能/動作を実施するための手段を作成するように、機械を生成することができる。
これらの機械実行可能命令又はコンピュータプログラム命令はコンピュータ、他のプログラマブルデータ処理装置、又は他のデバイスに特定の方法で機能するように指示することができるコンピュータ可読媒体に記憶されてもよく、その結果、コンピュータ可読媒体に記憶される命令はフローチャート及び/又はブロック図の一つ又は複数のブロックに特定される機能/動作を実施する命令を含む製造品を生成する。
機械実行可能命令又はコンピュータプログラム命令はまた、コンピュータ、他のプログラマブルデータ処理装置、又は他のデバイス上にロードされて、コンピュータ、他のプログラマブル装置、又は他のデバイス上で実行される一連の動作ステップをコンピュータ上で実行させて、コンピュータ又は他のプログラマブル装置上で実行される命令がフローチャート及び/又はブロック図のブロック又はブロック内で指定される関数/動作を実装するためのプロセスを提供するように、コンピュータ実装プロセスを生成することができる。本明細書で使用される「ユーザインターフェース」は、ユーザ又はオペレータがコンピュータ又はコンピュータシステムと対話することを可能にするインターフェースである。
「ユーザインターフェース」は「ヒューマンインターフェースデバイス」と呼ばれることもあり、ユーザインターフェースはオペレータに情報又はデータを提供し、及び/又はオペレータから情報又はデータを受信することができる。ユーザインターフェースはオペレータからの入力がコンピュータによって受信されることを可能にすることができ、コンピュータからユーザに出力を提供することができる。言い換えれば、ユーザインターフェースはオペレータがコンピュータを制御又は操作することを可能にすることができ、インターフェースは、コンピュータがオペレータの制御又は操作の効果を示すことを可能にすることができる。ディスプレイ又はグラフィカルユーザインタフェース上のデータ又は情報の表示は、オペレータに情報を提供する一例である。キーボード、マウス、トラックボール、タッチパッド、ポインティングスティック、グラフィックスタブレット、ジョイスティック、ゲームパッド、ウェブカメラ、ヘッドセット、ペダル、ワイヤードグローブ、リモートコントロール、及び加速度計を介したデータの受信は、オペレータからの情報又はデータの受信を可能にするユーザインターフェース構成要素のすべての例である。
本明細書で使用される「ハードウェアインターフェース」はコンピュータシステムの計算システムが外部コンピューティングデバイス及び/又は装置と対話し、及び/又はそれを制御することを可能にするインターフェースを包含する。ハードウェアインターフェースは、コンピュータシステムが制御信号又は命令を外部コンピューティングデバイス及び/又は装置に送信することを可能にし得る。ハードウェアインターフェースはまた、計算システムが外部コンピューティングデバイス及び/又は装置とデータを交換することを可能にし得る。ハードウェアインターフェースの例としてはユニバーサルシリアルバス、IEEE 1394ポート、パラレルポート、IEEE 1284ポート、シリアルポート、RS232ポート、IEEE488ポート、Bluetooth接続、無線ローカルエリアネットワーク接続、TCP/IP接続、イーサネット接続、制御電圧インターフェース、MIDIインターフェース、アナログ入力インターフェース、及びデジタル入力インターフェースが挙げられるが、これらに限定されない。
本明細書で使用される「ディスプレイ」又は「ディスプレイデバイス」は、画像又はデータを表示するように適合される出力デバイス又はユーザインターフェースを包含する。ディスプレイは、視覚的データ、音声データ、及び/又は触覚データを出力することができる。ディスプレイの例としてはコンピュータモニタ、テレビスクリーン、タッチスクリーン、触覚電子ディスプレイ、点字スクリーン、陰極線管(CRT)、記憶管、双安定ディスプレイ、電子ペーパー、ベクトルディスプレイ、フラットパネルディスプレイ、真空蛍光ディスプレイ(VF)、発光ダイオード(LED)ディスプレイ、エレクトロルミネセンスディスプレイ(ELD)、プラズマディスプレイパネル(PDP)、液晶ディスプレイ(LCD)、有機発光ダイオードディスプレイ(OLED)、プロジェクタ、及びヘッドマウントディスプレイが挙げられるが、これらに限定されない。
K空間データは、本明細書では磁気共鳴CTスキャン中に磁気共鳴装置のアンテナを使用して原子スピンによって放出される高周波信号の記録される測定値であると定義される。
磁気共鳴データは、断層画像データの一例である。磁気共鳴撮像(MRI)画像又はMR画像は、本明細書では磁気共鳴撮像データ内に含まれる解剖学的データの再構成される2次元又は3次元視覚化として定義される。この視覚化は、コンピュータを使用して実行することができる。
以下において、本発明の好ましい実施形態が、単なる例として、図面を参照して説明される。
医用システムの一例を示す。 図1の医用システムを動作させる方法を示すフローチャートを示す。 医用システムのさらなる例を示す。 助言領域を選択する方法を示す。
これらの図における同様の番号を付される要素は、同等の要素であるか、又は同じ機能を実行する。前述の要素は、機能が等価である場合、必ずしも後の図で論じられるとは限らない。
図1は、医用システム100の一例を示す。この例では、医用システム100が計算システム106を有するコンピュータを備える。コンピュータ102は、共にネットワーク化される1つ又は複数のコンピュータシステムを表すことができる。コンピュータシステム106は例えば、1つ以上の処理装置、マイクロ制御装置、又はフィールドプログラマブルゲートアレイを表すことができる。計算システム106は、オプションのハードウェアインターフェース104、オプションのユーザインターフェース108、及びメモリ110に接続されているものとして示されている。ハードウェアインターフェース104は例えば、医用システム100の他の構成要素が存在する場合、それらを制御するために計算システム106によって使用され得る。ユーザインターフェース108は例えば、オペレータが医用装置100を制御し、磁気共鳴画像又は他の情報をディスプレイするための手段を提供することができる。メモリ110は、計算システム106がアクセス可能な任意のタイプのメモリを表す。メモリ110の内容は例えば、複数の記憶装置又はメモリデバイス間で分割されてもよい。
メモリ110は、機械実行可能命令120を含むものとして示されている。機械実行可能命令120は、計算システム106が医用システム100の動作及び機能を制御すること、ならびに任意の他の追加成分を制御することを可能にする命令を含む。機械実行可能命令120はまた、計算システム106が様々な数学的演算、データ処理タスク、及び画像処理タスクを実行することを可能にし得る。メモリ110はさらに、エコー平面拡散強調磁気共鳴画像122を含むものとして示されている。メモリ110は、脂肪抑制T2強調磁気共鳴画像124を含むものとしてさらに示されている。
メモリ110は、エコー平面拡散強調磁気共鳴画像122と脂肪抑制T2強調磁気共鳴画像124との間の任意の位置合わせ126を含むものとして示されている。いくつかの事例では、2つの磁気共鳴画像122及び124がそれらのボクセル間に1対1の関係があるように取得され得る。他の事例では被験者が場合によっては移動していてもよく、位置合わせ126は2つの間のデータを使用することができるようにマップを含んでもよい。エコー平面拡散強調磁気共鳴画像122はその中に歪みを有することがあり、その場合、位置合わせからのわずかな逸脱は、医用システム100の動作に悪影響を及ぼさないことがあることにも留意される。これはまた、位置合わせを不要にし得る。
メモリ110はさらに、エコー平面拡散強調磁気共鳴画像122において識別される多数の高拡散率領域128を含むものとして示されている。これは、例えば、拡散のための所定の閾値に従って拡散率を閾値化することによって実行されてもよい。メモリ110は、脂肪抑制T2強調磁気共鳴画像124において識別される組織領域130を含むものとしてさらに示されている。組織領域130は例えば、対象の組織の位置を識別する非脂肪又は非脂質組織領域であってもよい。メモリ110は、画像処理モジュール132を含むものとしてさらに示されている。画像処理モジュール132は、組織領域130から選択されるいくつかの勧告領域134を出力する。それらは、高拡散率領域128を組織領域130に相関させるか又はマッチングさせることによって選択される。メモリ110は、脂肪抑制T2強調磁気共鳴画像124上に重ねられた又は強調される勧告領域134を示す任意選択の増強T2強調磁気共鳴画像136を含むものとして更に示されている。
図2は、図1の医用システム100を動作させる方法を示すフローチャートを示す。まず、ステップ200において、エコー平面拡散強調磁気共鳴画像122が受信される。次に、ステップ202において、脂肪抑制T2強調磁気共鳴画像124が受信される。次に、ステップ204において、エコー平面拡散強調磁気共鳴画像122がセグメント化されて、高拡散率領域128が識別される。次に、ステップ206において、脂肪抑制T2強調磁気共鳴画像124は、組織領域130を識別するためにセグメント化される。ステップ208では、高拡散率領域128及び組織領域130が画像処理モジュール132に入力されて、組織領域130の一部を勧告領域134として提供又は識別する。最後に、ステップ210において、勧告領域134は、脂肪抑制T2強調磁気共鳴画像のセグメント化として提供される。
図3は、医用システム300のさらなる例を示す。図3の医用システム300は図1の医用システム100と同様であるが、計算システム106によって制御することができる磁気共鳴撮像システム302をさらに備える点が異なる。
磁気共鳴撮像システム302は、磁石304を備える。磁石304は、それをタイプボア306を有する超伝導円筒形磁石である。異なるタイプの磁石の使用も可能であり、例えば、分割円筒形磁石といわゆる開放磁石の両方を使用することも可能である。分割円筒形磁石はクライオスタットが磁石の等平面へのアクセスを可能にするために2つのセクションに分割されていることを除いて、規格円筒形磁石に類似しており、そのような磁石は例えば、荷電粒子ビーム治療と併せて使用され得る。開放磁石は2つの磁石セクションを有し、一方が他方の上に、被験者を受け入れるのに十分な大きさの間の空間、すなわちヘルムホルツコイルの領域と同様の2つのセクション領域の配置を有する。開放磁石は、被検体の閉じ込めがより少ないため、人気がある。円筒形磁石のクライオスタットの内部には、超伝導コイルの集合体がある。
円筒形磁石304のボア306内には、磁場が強く、磁気共鳴CTを実行するのに十分に均一である撮像ゾーン308がある。関心領域309が撮像ゾーン308内に示されている。磁気共鳴データは、関心領域から取得される。被験者318は被験者318の少なくとも一部が撮像ゾーン308及び関心領域309内にあるように、被験者サポート320によって支持されているように示されている。関心領域309内には、被験者318の乳房組織が存在することが分かる。
磁石のボア306内には、磁石304の撮像ゾーン308内の磁気スピンを空間的に符号化するための予備的な磁気共鳴データの取得のために使用される一組の磁場勾配コイル310もある。傾斜磁場コイル310は、傾斜磁場コイル電源312に接続される。磁場勾配コイル310は、代表的であることが意図される。典型的には、磁場勾配コイル310が3つの直交する空間方向に空間的に符号化するための3つの別個のコイルセットを含む。傾斜磁場電源は、傾斜磁場コイルに電流を供給する。磁場勾配コイル310に供給される電流は、時間の関数として制御され、ランプ又はパルス化され得る。
撮像ゾーン308に隣接して、撮像ゾーン308内の磁気スピンの向きを操作し、また撮像ゾーン308内のスピンからの無線送信を受信するための無線周波数コイル314がある。無線周波数アンテナは、複数のコイル要素を含むことができる。無線周波数アンテナは、チャネル又はアンテナと呼ばれることもある。高周波コイル314は、無線周波数トランシーバ316に接続される。無線周波数コイル314及び無線周波数トランシーバ316は、別個の送信コイル及び受信コイル、ならびに別個の送信機及び受信機によって置き換えられ得る。無線周波数コイル314及び無線周波数トランシーバ316が代表的であることが理解される。無線周波数コイル314はまた、専用送信アンテナ及び専用受信アンテナを表すことを意図している。同様に、トランシーバ316は、別個の送信機及び受信機を表すこともできる。無線周波数コイル314はまた、複数の受信/送信素子を有することができ、無線周波数トランシーバ316は、複数の受信/送信チャネルを有することができる。例えば、検知などの並列撮像技術が実行される場合、無線周波数314は、複数のコイル要素を有することができる。
トランシーバ316及び勾配コントローラ312は、コンピュータシステム102のハードウェアインターフェース106に接続されているものとして示されている。
メモリ110はさらに、エコー平面拡散強調磁気共鳴撮像プロトコルに従って第1のk空間データを取得するように構成される第1のパルスシーケンスコマンド330を含むものとして示されている。メモリ110は、第1のパルスシーケンスコマンド330で磁気共鳴撮像システム302を制御することによって取得される第1のk空間データ332を含むものとしてさらに示されている。エコー平面拡散強調磁気共鳴画像122は、第1のk空間データ332から再構成される。メモリ110は、第2のパルスシーケンスコマンド334を含むものとしてさらに示されている。第2のパルスシーケンスコマンド334は、T2強調磁気共鳴撮像プロトコルに従って第2のk空間データを取得するように構成される。
図3に示す例では、第2のパルスシーケンスコマンド334が電気的特性トモグラフィに適した位相測定値を取得するように構成される。メモリ110は、第2のパルスシーケンスコマンド334で磁気共鳴撮像システム302を制御することによって取得される第2のk空間データ336を含むものとして示されている。第2のk空間データ336は、脂肪抑制T2強調磁気共鳴画像124及び導電率マップ338を再構成するために使用される。導電率マップ338は、メモリ110に記憶されているものとして示されている。
乳癌研究のための磁気共鳴(MR)方法は通常、造影剤の投与(例えば、動的造影(DCE)撮像のため)に基づいており、これは、一般的なMRベースの乳癌スクリーニングのための主要な欠点と見なされる。
「電気的特性断層撮影」(EPT)は一般に、DCEと同じ乳房病変特徴付けのポテンシャルを有し得る。しかしながら、最新技術における導電率の決定は信頼できる腫瘍境界検出のための造影剤の支持を必要とし、したがって、乳癌EPTの造影剤を含まないバージョンは、依然として欠けている。
例は(a)腫瘍組織のみを示すが、幾何学的に歪んだ(エコー平面拡散強調磁気共鳴画像122)第1の高感度MRシーケンスA 330(例えば、高b値EPIDWI)を実行し、(b)多くの潜在的病変を示すが、幾何学的歪みを伴わない、第2のEPTに適したMRシーケンスB 334(例えば、脂肪抑制T2強調画像)を実行し(脂肪抑制T2強調磁気共鳴画像)、(c)シーケンスAの歪んだ病変とシーケンスBの潜在的病変との間の最大幾何学的オーバラップを探索することによって、シーケンスBの正しい腫瘍病変を識別し(助言領域134を識別し)、任意選択的に(d)シーケンスBから再構成されるその導電率によってこの組織のタイプを特徴付けることができる。
EPTの主な可能性は、腫瘍の同定(ローカライゼーション)ではなく、疑わしい勧告領域又は腫瘍の可能な特徴付けである。腫瘍の同定/ローカライゼーションは、通常、DCEによって行われる。例は、DCEを第1の配列(配列A)で置き換えることを示唆する。シーケンスAは例えば、高b値を有するEPIベースの拡散強調画像化(DWI)によって与えられるように、DCEと同様に感受性である。第2のステップにおいて、EPI-DWIによって同定される潜在的病変は、EPTによって特徴付けられる。EPTは後処理EPI-DWIによって実行することができないので、第2の専用EPTスキャン(シーケンスB、例えば、脂肪抑制T2強調(T2w)撮像)が必要である。EPTを可能にするために、シーケンスBは(i)純粋にB1関連であり、(ii)腫瘍対脂肪及び腫瘍対腺組織の間の十分な内因性造影を有する位相を有する。(i) 純粋なB1関連相は、スピンエコーにベースシーケンス、平衡勾配を有する配列(定常状態自由歳差運動)、又は超短/ゼロTE配列によって与えることができる。マルチエコー勾配エコーシーケンスは原理的に、B1関連位相と等価なTE=0への位相発展のボクセル毎の外挿を可能にするが、潜在的に不安定であり、したがって誤差が起こりやすい外挿アルゴリズムを必要とする。(ii) 腫瘍対脂肪と腫瘍対腺組織との間の十分な固有のコントラストは有益であり得、なぜなら、位相の数値的区別は、好ましくは同様の導電率を有するボクセルのみを含む、標的ボクセル(いわゆる「カーネル」)の周りの隣接ボクセルの集合体の使用を必要とするからである。したがって、腫瘍境界についての先験的な知識は、腫瘍体積内に留まるようにカーネルを局所的に形作るために必要とされる。これらの腫瘍境界は、異なるスキャンから取得することができる。しかしながら、別個スキャンから腫瘍境界及び位相を取得することは、不完全な画像位置合わせのリスクのために常に問題である。特に、小さな腫瘍の場合、誤って数個の誤って登録される非腫瘍ボクセルのみを含むと、導電率の結果を著しく改ざんする可能性がある。
図4は、勧告領域134を識別する一手段を示す。図4は、脂肪抑制T2強調磁気共鳴画像124を表すスケッチを示す。130及び134とラベル付けされる領域は、組織領域130を表す。領域128は、エコー平面拡散強調磁気共鳴画像122において識別される高拡散率領域を表す。エコー平面拡散強調磁気共鳴画像122は歪みを含み、勧告領域134の位置を確実に識別することができない。これを調べると、曲線128と曲線134との間にオーバラップ領域400があることが分かる。次いで、領域134は、勧告領域として識別される。2つの比較から、領域134は高い拡散率を有する領域であり、したがって、医師によってさらに調査されるべきであることが明らかである。
オーバラップ400の領域は、曲線134の面積と比較することができる。これが一定の割合を上回る場合、これは、勧告領域134の識別をトリガするために使用され得る所定のオーバラップであり得る。同様に、領域128及び134の中心又は重心も識別することができ、それらの間の距離を比較することができる。また、図4を調べると、領域128の体積と領域134の体積とが最も類似していることが分かる。これは、勧告領域134を識別する代替手段であってもよい。領域128の形状及び領域134の形状は、他の組織領域130と比較して最も類似していることも分かる。したがって、領域128乃至134を一致させるために必要とされる形状又は変形も、勧告領域134を識別するために使用され得る。
本発明は図面及び前述の説明において詳細に図示及び説明されているが、そのような図示及び説明は例示的又は例示的であり、限定的ではないと見なされるべきであり、本発明は開示される実施形態に限定されない。
開示される実施形態に対する他の変形は図面、開示、及び添付の特許請求の範囲の研究から、特許請求される発明を実施する際に当業者によって理解され、達成され得る。請求項において、単語「有する(comprising)」は他の要素又はステップを除外せず、不定冠詞「a」又は「an」は複数を除外しない。単一のプロセッサ又は他のユニットは、特許請求の範囲に列挙されるいくつかのアイテムの機能を満たすことができる。特定の手段が相互に異なる従属請求項に記載されているという単なる事実は、これらの手段の組み合わせが有利に使用され得ないことを示すものではない。コンピュータプログラムは他のハードウェアとともに、又は他のハードウェアの一部として供給される光記憶媒体又はソリッドステート媒体などの適切な媒体上に記憶/配布され得るが、インターネット又は他の有線もしくはワイヤレス電気通信システムなどの他の形態で配布されてもよい。請求項におけるいかなる参照符号も、範囲を限定するものとして解釈されるべきではない。
100 医療システム
102 コンピューター
104 ハードウェアインターフェイス
106 計算システム
108 ユーザー インターフェイス
110 メモリ
120 機械実行可能命令
122 エコー平面拡散強調磁気共鳴画像
124 脂肪抑制 T2 強調磁気共鳴画像
126 任意レジストレーション
128 高拡散率領域
130 組織領域
132 画像処理モジュール
134 勧告領域
136 拡張 T2 強調磁気共鳴画像
200 乳房組織を表す関心領域のエコー平面拡散強調磁気共鳴画像を受信
202 関心領域を表す脂肪抑制T2強調磁気共鳴画像を受信
204 エコー平面拡散強調磁気共鳴画像をセグメント化して、拡散率の高い領域を特定
206 脂肪抑制T2強調磁気共鳴画像をセグメント化して、組織領域を特定
208 高拡散率領域及び組織領域を画像処理モジュールに入力することにより、組織領域の一部を勧告領域として識別
210 脂肪抑制T2強調磁気共鳴画像のセグメンテーションとして勧告領域を提供
300 医療システム
302 磁気共鳴画像システム
304 磁石
306 ボアの磁石
308 イメージングゾーン
309 関心領域
310 磁場傾斜コイル
312 磁場傾斜コイル電源
314 高周波コイル
316 トランシーバー
318 被検体
320 被検体サポート
330 第一のパルス シーケンス コマンド
332 第一のk空間データ
334 第2のパルスシーケンスコマンド
336 第2の k 空間データ
338 導電率マップ
400 オーバーラップ領域

Claims (14)

  1. 医用システムであって、
    機械実行可能命令を記憶するメモリと、
    前記医用システムを制御するように構成される計算システムであって、前記機械実行可能命令は前記計算システムに、
    乳房組織を記述する関心領域のエコー平面拡散強調磁気共鳴画像を受信するステップと、
    前記関心領域を記述する脂肪抑制T2強調磁気共鳴画像を受信するステップと、
    前記エコー平面拡散強調磁気共鳴画像をセグメント化して高拡散率領域を識別するステップであって、前記高拡散率領域は、所定の拡散率閾値を上回る拡散率を有する領域である、ステップと、
    前記脂肪抑制T2強調磁気共鳴画像をセグメント化して組織領域を識別するステップと、
    画像処理モジュールに前記高拡散率領域及び前記組織領域を入力することによって前記組織領域の一部を勧告領域として識別するステップと、
    前記関心領域を記述する導電率マップを受信するステップと、
    前記導電率マップを用いて前記勧告領域の各々の導電率測定値を計算するステップと、
    前記導電率測定値に従って前記勧告領域の各々の分類を割り当てるステップと、
    前記脂肪抑制T2強調磁気共鳴画像のセグメント化画像として各々に前記分類が割り当てられた前記勧告領域を提供するステップと
    を実行させる、医用システム。
  2. 前記エコー平面拡散強調磁気共鳴画像は、高b値エコー平面拡散強調磁気共鳴画像である、請求項1に記載の医用システム。
  3. 前記医用システムは、撮像ゾーンからk空間データを取得するように構成される磁気共鳴撮像システムをさらに有し、前記メモリは、第1のパルスシーケンスコマンド及び第2のパルスシーケンスコマンドをさらに含み、前記第1のパルスシーケンスコマンドはエコー平面拡散強調磁気共鳴撮像プロトコルに従って第1のk空間データを取得するように構成され、前記第2のパルスシーケンスコマンドはT2強調磁気共鳴撮像プロトコルに従って第2のk空間データを取得するように構成され、前記第2のパルスシーケンスコマンドは電気的特性トモグラフィ磁気共鳴撮像プロトコルに従っており、前記第2のパルスシーケンスコマンドは脂肪抑制であり、前記機械実行可能命令は前記計算システムに、
    前記第1のパルスシーケンスコマンドを用いて前記磁気共鳴撮像システムを制御することによって、前記乳房組織を含む前記関心領域から前記第1のk空間データを取得するステップと、
    前記第2のパルスシーケンスコマンドを用いて前記磁気共鳴撮像システムを制御することによって、前記乳房組織を含む前記関心領域から前記第2のk空間データを取得するステップと、
    前記第1のk空間データから前記エコー平面拡散強調磁気共鳴画像を再構成するステップと、
    前記第2のk空間データから前記脂肪抑制T2強調磁気共鳴画像を再構成するステップと、
    前記第2のk空間データから前記関心領域を記述する導電率マップを再構成するステップと
    を実行させる、請求項1又は2に記載の医用システム。
  4. 前記第1のパルスシーケンスコマンドは、0s/mm2及び3000s/mm2の間のb値を有する前記第1のk空間データを取得するように構成される、請求項3に記載の医用システム。
  5. 前記第2のパルスシーケンスコマンドは、
    スピンエコーパルスシーケンスコマンドと、
    超短エコー時間パルスシーケンスコマンドと、
    ゼロエコー時間パルスシーケンスコマンドと、
    マルチエコーグラジエントエコーパルスシーケンスコマンドと、
    バランスド・グラジエントエコーパルスシーケンスコマンドと、
    定常状態自由歳差パルスシーケンスコマンドと
    の何れか1つである、請求項3又は4に記載の医用システム。
  6. 前記勧告領域は、前記導電率マップの計算中に境界を定義するために使用される、請求項3、4、又は5に記載の医用システム。
  7. 前記機械実行可能命令の実行は前記計算システムにさらに、前記勧告領域を、前記脂肪抑制T2強調磁気共鳴画像上のオーバーレイとしてレンダリングさせる、請求項1乃至6の何れか一項に記載の医用システム。
  8. 前記機械実行可能命令の実行は前記計算システムに、前記導電率測定値を使用して前記勧告領域のレンダリングをランク付けさせる、請求項7に記載の医用システム。
  9. 前記高拡散率領域は、前記エコー平面拡散強調磁気共鳴画像を閾値処理することによって識別される、請求項1乃至8の何れか一項に記載の医用システム。
  10. 前記画像処理モジュールは、
    前記組織領域の1つが、高拡散率領域の何れかと、所定の割合を上回るオーバラップを有する場合、前記組織領域の1つを前記勧告領域の1つとして識別するステップと、
    前記組織領域の1つが、所定の距離を下回る高い拡散率領域の何れかを有する中心距離を有する場合、前記組織領域の1つを前記勧告領域の1つとして識別するステップと、
    前記組織領域の1つが、所定の体積差内の高拡散率領域の何れかと一致する体積を有する場合、前記組織領域の1つを前記勧告領域の1つとして識別するステップと、
    前記組織領域の1つが、所定の歪み内の高拡散率領域の何れかに一致する形状を有する場合、前記組織領域の1つを前記勧告領域の1つとして識別するステップと、
    それらの組み合わせと
    の何れか1つを実行することによって、前記勧告領域をアルゴリズム的に識別するように構成される、請求項1乃至9の何れか1項に記載の医用システム。
  11. 前記画像処理モジュールは、前記組織領域及び前記高拡散率領域を入力することに応答して、前記組織領域の一部を勧告領域としてラベル付けするように構成される、訓練されるニューラルネットワークである、請求項1乃至10の何れか一項に記載の医用システム。
  12. 前記組織領域は非脂肪組織領域である、請求項1乃至11の何れか一項に記載の医用システム。
  13. 計算システムによる実行のための機械実行可能命令を有するコンピュータプログラムであって、前記機械実行可能命令は前記計算システムに、
    乳房組織を記述する関心領域のエコー平面拡散強調磁気共鳴画像を受信するステップと、
    前記関心領域を記述する脂肪抑制T2強調磁気共鳴画像を受信するステップと、
    前記エコー平面拡散強調磁気共鳴画像をセグメント化して高拡散率領域を識別するステップであって、前記高拡散率領域は、所定の拡散率閾値を上回る拡散率を有する領域である、ステップと、
    前記脂肪抑制T2強調磁気共鳴画像をセグメント化して組織領域を識別するステップと、
    前記高拡散率領域及び前記組織領域を画像処理モジュールに入力することによって、前記組織領域の一部を勧告領域として識別するステップと、
    前記関心領域を記述する導電率マップを受信するステップと、
    前記導電率マップを用いて前記勧告領域の各々の導電率測定値を計算するステップと、
    前記導電率測定値に従って前記勧告領域の各々の分類を割り当てるステップと、
    前記脂肪抑制T2強調磁気共鳴画像のセグメンテーションとして前記勧告領域を提供するステップと
    を実行させる、コンピュータプログラム。
  14. 医用撮像の方法であって、前記方法は、
    乳房組織を記述する関心領域のエコー平面拡散強調磁気共鳴画像を受信するステップと、
    前記関心領域を記述する脂肪抑制T2強調磁気共鳴画像を受信するステップと、
    前記エコー平面拡散強調磁気共鳴画像をセグメント化して高拡散率領域を識別するステップであって、前記高拡散率領域は、所定の拡散率閾値を上回る拡散率を有する領域である、ステップと、
    前記脂肪抑制T2強調磁気共鳴画像をセグメント化して組織領域を識別するステップと、
    前記高拡散率領域及び前記組織領域を画像処理モジュールに入力することによって前記組織領域の一部を勧告領域として識別するステップと、
    前記関心領域を記述する導電率マップを受信するステップと、
    前記導電率マップを用いて前記勧告領域の各々の導電率測定値を計算するステップと、
    前記導電率測定値に従って前記勧告領域の各々の分類を割り当てるステップと、
    前記脂肪抑制T2強調磁気共鳴画像のセグメントテーションとして前記勧告領域を提供するステップと、
    を有する、方法。
JP2022555738A 2020-03-26 2021-03-17 乳房磁気共鳴ctにおける勧告領域の特定 Active JP7284876B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP20165864.8 2020-03-26
EP20165864.8A EP3886039A1 (en) 2020-03-26 2020-03-26 Identification of advisory regions in breast magnetic resonance imaging
PCT/EP2021/056746 WO2021191007A1 (en) 2020-03-26 2021-03-17 Identification of advisory regions in breast magnetic resonance imaging

Publications (2)

Publication Number Publication Date
JP2023509237A JP2023509237A (ja) 2023-03-07
JP7284876B2 true JP7284876B2 (ja) 2023-05-31

Family

ID=70049857

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022555738A Active JP7284876B2 (ja) 2020-03-26 2021-03-17 乳房磁気共鳴ctにおける勧告領域の特定

Country Status (5)

Country Link
US (1) US20230152404A1 (ja)
EP (2) EP3886039A1 (ja)
JP (1) JP7284876B2 (ja)
CN (1) CN114173650A (ja)
WO (1) WO2021191007A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2617371A (en) * 2022-04-06 2023-10-11 Perspectum Ltd Method and apparatus for characterization of breast tissue using multiparametric MRI

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100277174A1 (en) 2007-12-21 2010-11-04 Koninklijke Philips Electronics N.V. PASSIVE SHIMS TO INCREASE THE EFFECTIVE B0 and B1 UNIFORMITY IN A BODY COIL
US20150153431A1 (en) 2013-11-27 2015-06-04 General Electric Company Systems and methods for determining electrical properties using magnetic resonance imaging

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9411033B2 (en) * 2005-05-11 2016-08-09 Regents Of The University Of Minnesota Methods and apparatus for imaging with magnetic induction
NO20101638A1 (no) 2010-11-22 2012-05-23 Sunnmore Mr Klinikk As Fremgangsmate ved ex vivo distingvering mellom maligne og benigne tumorer ved anvendelse av kontrastmiddelbasert MR-skanning
WO2015162694A1 (ja) * 2014-04-22 2015-10-29 株式会社日立製作所 医用画像診断支援装置、磁気共鳴イメージング装置および医用画像診断支援方法
EP3833246A4 (en) * 2018-08-12 2022-05-11 The Trustees of Columbia University in the City of New York SYSTEM, METHOD AND COMPUTER-ACCESSIBLE MEDIUM FOR NON-INVASIVE TEMPERATURE ESTIMATION

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100277174A1 (en) 2007-12-21 2010-11-04 Koninklijke Philips Electronics N.V. PASSIVE SHIMS TO INCREASE THE EFFECTIVE B0 and B1 UNIFORMITY IN A BODY COIL
US20150153431A1 (en) 2013-11-27 2015-06-04 General Electric Company Systems and methods for determining electrical properties using magnetic resonance imaging

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Jaewook SHIN et al.,"Initial study on in vivo conductivity mapping of breast cancer using MRI",Journal of Magnetic Resonance Imaging,2014年11月21日,Vol. 42, No. 2,p.371-378,DOI: 10.1002/jmri.24803

Also Published As

Publication number Publication date
EP3886039A1 (en) 2021-09-29
US20230152404A1 (en) 2023-05-18
JP2023509237A (ja) 2023-03-07
WO2021191007A1 (en) 2021-09-30
EP3977404B1 (en) 2023-05-10
CN114173650A (zh) 2022-03-11
EP3977404A1 (en) 2022-04-06

Similar Documents

Publication Publication Date Title
Yang et al. Evaluation of radiomic texture feature error due to MRI acquisition and reconstruction: a simulation study utilizing ground truth
US10330757B2 (en) MRI method for calculating derived values from B0 and B1 maps
CN110476075B (zh) 针对解剖区域的磁共振指纹识别词典的选择
US10823798B2 (en) Virtual CT images from magnetic resonance images
US10591562B2 (en) Bone MRI using B0 inhomogeneity map and a subject magnetic susceptibility map
US11092659B2 (en) Sub voxel resolution magnetic resonance fingerprinting imaging
JP7284876B2 (ja) 乳房磁気共鳴ctにおける勧告領域の特定
US11353532B2 (en) (3-n) dimensional determination of electric conductivity
US20230377112A1 (en) Magnetic Resonance Imaging of Breast Micro-Calcifications
EP3892191A1 (en) Measurement of intra-cellular conductivity using magnetic resonance imaging
US20210259569A1 (en) Electical properties tomography mapping of conductivity changes

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220915

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20220915

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230214

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230424

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230519

R150 Certificate of patent or registration of utility model

Ref document number: 7284876

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150