WO2021191748A1 - Method for extracting multiple t2* values from single set of cpmg data - Google Patents

Method for extracting multiple t2* values from single set of cpmg data Download PDF

Info

Publication number
WO2021191748A1
WO2021191748A1 PCT/IB2021/052234 IB2021052234W WO2021191748A1 WO 2021191748 A1 WO2021191748 A1 WO 2021191748A1 IB 2021052234 W IB2021052234 W IB 2021052234W WO 2021191748 A1 WO2021191748 A1 WO 2021191748A1
Authority
WO
WIPO (PCT)
Prior art keywords
values
cpmg
signals
time domain
segments
Prior art date
Application number
PCT/IB2021/052234
Other languages
French (fr)
Inventor
Saul Stokar
Tal Bareket
Zachi Peles
Original Assignee
Clear-Cut Medical Ltd.
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 Clear-Cut Medical Ltd. filed Critical Clear-Cut Medical Ltd.
Publication of WO2021191748A1 publication Critical patent/WO2021191748A1/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/445—MR involving a non-standard magnetic field B0, e.g. of low magnitude as in the earth's magnetic field or in nanoTesla spectroscopy, comprising a polarizing magnetic field for pre-polarisation, B0 with a temporal variation of its magnitude or direction such as field cycling of B0 or rotation of the direction of B0, or spatially inhomogeneous B0 like in fringe-field MR or in stray-field imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/448—Relaxometry, i.e. quantification of relaxation times or spin density
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE

Definitions

  • the present invention relates generally to NMR or MRI methods and systems, and particularly to a method for extracting multiple T2* values from a single set of CPMG data.
  • NMR nuclear magnetic resonance
  • MRI magnetic resonance imaging
  • a weighted average of all the points in the spectrum is calculated.
  • Zero weights are given to points outside the spectrum from the tissue (i.e. noise points), in order to improve the signal-to-noise ratio (SNR), and weights that depend on the amplitude are given to points within the spectrum.
  • SNR signal-to-noise ratio
  • the non-zero weights are given to the points whose frequencies lie approximately in the range [-75, +50] KHz, with the largest weights given to the points approximately in the range [-25, 25] KHz.
  • the prior art version of the CLEAR SIGHT system uses such weights.
  • the signal of weighted average vs. echo number is fitted to an exponential, known as the T2 value; high T2 values are an indication that the voxel contains cancerous tissue.
  • T2 value an indication that the voxel contains cancerous tissue.
  • the presence of a Bo gradient during signal sampling implies that the different points in the spectrum (i.e., the different frequencies) contain the contributions of different positions in space. Indeed, if the Bo field gradient were purely in the Z direction (that is, perpendicular to the tissue surface), the different points in the spectrum would correspond to different depths into the tissue. In practice, the gradient may have a X-Y component as well (i.e. parallel to the tissue surface), complicating the simple geometrical interpretation, but the essential concept remains the same; namely, different points in the spectrum (i.e., the different frequencies) contain the contributions of different positions in space.
  • the present invention seeks to provide improved methods for extracting multiple T2* values from a single set of CPMG data, as is describe in detail hereinbelow.
  • a method for extracting multiple T2* values from a single set of CPMG data including using an MRI system that includes an RF coil to image tissue and generate NMR time domain signals via a CPMG sequence, converting the signals to a set of spectra, one spectrum per echo, by performing a Fourier transform of complex time domain signals associated with the signals, and dividing each of the spectra into a number of segments and calculating a T2* value for each segment, wherein the T2* values represent decay rates for a subset of frequencies in the spectrum of the signal.
  • the method may use equal frequency limits for all echoes in the CPMG sequence.
  • the segments may include overlapping segments.
  • Fig. 1A is a graphical illustration of a simulated signal in the time domain of a prior art MRI system, the signal being generated by the MRI system imaging tissue via a CPMG sequence in the presence of a strong magnetic field with a strong gradient;
  • Fig. IB is a graphical illustration of the Fourier transform of the time-domain signal of Fig. 1A.
  • Fig. 2 is a simplified flow chart of a method for extracting multiple T2* values from a single set of CPMG data, in accordance with a non-limiting embodiment of the present invention.
  • Fig. 2 illustrates a method for extracting multiple T2* values from a single set of CPMG data in accordance with a non-limiting embodiment of the present invention.
  • An MRI system is used to image tissue and generate NMR time domain signals via a CPMG sequence in the presence of a strong magnetic field with a strong gradient (step 1).
  • the signals are converted to a set of spectra (one spectrum per echo) by performing a Fourier transform of the (complex) time domain signals (step 2).
  • the spectrum is divided into a number of segments, and a T2 * value is calculated for each segment (step 3).
  • the same frequency limits may be used for all the echoes in the CPMG sequence (step 4). Accordingly, multiple T2 values are calculated for each voxel. For example, if the spectrum is split into two segments, two T2 * values, T2 * (l) and T 2 * (2) are calculated. If the gradient is purely in the Z direction, then T2 (1) represents the decay rate for tissue closer to the RF coil, and
  • T2 (2) represents the decay rate for tissue further from the RF coil (i.e., deeper into the tissue).
  • the gradient is in an arbitrary direction
  • T2 (1) and T2 (2) represent the signals from different volumes in the tissue.
  • More than two segments may be chosen (as many as there are points in the spectrum), including overlapping segments (step 5); e.g., the first 50% of the frequencies, the middle 50% of the frequencies and the last 50% of the frequencies and many other variations. If any one of the segments exhibits a high T2 value, it can be displayed to the user as a possible indication that cancer is present in that segment of this voxel.
  • the advantages of splitting the spectrum into multiple segments is that (a) this provides knowledge of where the high T2 value comes from (i.e., how deep into the tissue) and (b) this effectively narrows the size of the effective voxel, narrowing the diluting effect of healthy tissue (i.e., the partial volume effect that occurs if the voxel contains a mixture of tumor and healthy tissue).
  • the disadvantage of splitting the spectrum is that the signal-to-noise ratio (SNR) of each segment is lower than the full spectrum.

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

A method for extracting multiple T2* values from a single set of CPMG data includes using an MRI system that includes an RF coil to image tissue and generate NMR time domain signals via a CPMG sequence, converting the signals to a set of spectra, one spectrum per echo, by performing a Fourier transform of complex time domain signals associated with the signals, and dividing each of the spectra into a number of segments and calculating a T2* value for each segment. The T2* values represent decay rates as a function of distance from the RF coil.

Description

METHOD FOR EXTRACTING MULTIPLE T2* VALUES FROM SINGLE SET OF
CPMG DATA
FIELD OF THE INVENTION
The present invention relates generally to NMR or MRI methods and systems, and particularly to a method for extracting multiple T2* values from a single set of CPMG data.
BACKGROUND OF THE INVENTION
It is known in the art of NMR (nuclear magnetic resonance) or MRI (magnetic resonance imaging) that the NMR signal in a so-called CPMG (Carr-Purcell-Meiboom- Gill) sequence decays at a rate that is strongly affected by the self-diffusion of the water in the tissue under examination. It is also well established that in the presence of a substantially uniform constant gradient in the Bo field (static magnetic field), the decay is exponential for a single value of D and T2 in the voxel under study. (As is known, D is the self-diffusion constant and T2 is the spin-spin decay constant.)
In the CLEAR SIGHT system of Clear Cut Medical Ltd., excised tissue is examined via a CPMG sequence in the presence of a strong magnetic field with a strong gradient. The RF coil and the hardware are described in US Patent 9310450. The sampled signal is converted to a set of spectra (one spectrum per echo) by performing the Fourier transform of the (complex) time domain signal. An example of a prior-art simulated signal and its spectrum are shown in Fig. 1A (which illustrates the simulated signal in the time domain) and Fig. IB (which illustrates the Fourier transform of the time-domain signal).
For each echo, a weighted average of all the points in the spectrum is calculated. Zero weights are given to points outside the spectrum from the tissue (i.e. noise points), in order to improve the signal-to-noise ratio (SNR), and weights that depend on the amplitude are given to points within the spectrum. For example, in the spectrum shown in Fig. IB, the non-zero weights are given to the points whose frequencies lie approximately in the range [-75, +50] KHz, with the largest weights given to the points approximately in the range [-25, 25] KHz. Using methods well known to those skilled in the art of mathematics, it is possible to calculate the weights that optimize the final SNR. The prior art version of the CLEAR SIGHT system uses such weights. The signal of weighted average vs. echo number is fitted to an exponential, known as the T2 value; high T2 values are an indication that the voxel contains cancerous tissue. The presence of a Bo gradient during signal sampling implies that the different points in the spectrum (i.e., the different frequencies) contain the contributions of different positions in space. Indeed, if the Bo field gradient were purely in the Z direction (that is, perpendicular to the tissue surface), the different points in the spectrum would correspond to different depths into the tissue. In practice, the gradient may have a X-Y component as well (i.e. parallel to the tissue surface), complicating the simple geometrical interpretation, but the essential concept remains the same; namely, different points in the spectrum (i.e., the different frequencies) contain the contributions of different positions in space.
SUMMARY OF THE INVENTION
The present invention seeks to provide improved methods for extracting multiple T2* values from a single set of CPMG data, as is describe in detail hereinbelow.
There is provided in accordance with a non-limiting embodiment of the invention a method for extracting multiple T2* values from a single set of CPMG data including using an MRI system that includes an RF coil to image tissue and generate NMR time domain signals via a CPMG sequence, converting the signals to a set of spectra, one spectrum per echo, by performing a Fourier transform of complex time domain signals associated with the signals, and dividing each of the spectra into a number of segments and calculating a T2* value for each segment, wherein the T2* values represent decay rates for a subset of frequencies in the spectrum of the signal.
The method may use equal frequency limits for all echoes in the CPMG sequence.
The segments may include overlapping segments.
BRIEF DESCRIPTION OF THE DRAWINGS
The present invention will be understood and appreciated more fully from the following detailed description, taken in conjunction with the drawings in which:
Fig. 1A is a graphical illustration of a simulated signal in the time domain of a prior art MRI system, the signal being generated by the MRI system imaging tissue via a CPMG sequence in the presence of a strong magnetic field with a strong gradient;
Fig. IB is a graphical illustration of the Fourier transform of the time-domain signal of Fig. 1A; and
Fig. 2 is a simplified flow chart of a method for extracting multiple T2* values from a single set of CPMG data, in accordance with a non-limiting embodiment of the present invention.
DETAILED DESCRIPTION OF EMBODIMENTS Reference is now made to Fig. 2, which illustrates a method for extracting multiple T2* values from a single set of CPMG data in accordance with a non-limiting embodiment of the present invention.
An MRI system is used to image tissue and generate NMR time domain signals via a CPMG sequence in the presence of a strong magnetic field with a strong gradient (step 1). The signals are converted to a set of spectra (one spectrum per echo) by performing a Fourier transform of the (complex) time domain signals (step 2). In accordance with a non-limiting embodiment of the present invention, instead of taking the weighted sum of all the (non-zero) points in the spectrum, the spectrum is divided into a number of segments, and a T2* value is calculated for each segment (step 3). The same frequency limits may be used for all the echoes in the CPMG sequence (step 4). Accordingly, multiple T2 values are calculated for each voxel. For example, if the spectrum is split into two segments, two T2* values, T2*(l) and T2 *(2) are calculated. If the gradient is purely in the Z direction, then T2 (1) represents the decay rate for tissue closer to the RF coil, and
T2 (2) represents the decay rate for tissue further from the RF coil (i.e., deeper into the tissue). However, even in the general case where the gradient is in an arbitrary direction,
T2 (1) and T2 (2) represent the signals from different volumes in the tissue.
More than two segments may be chosen (as many as there are points in the spectrum), including overlapping segments (step 5); e.g., the first 50% of the frequencies, the middle 50% of the frequencies and the last 50% of the frequencies and many other variations. If any one of the segments exhibits a high T2 value, it can be displayed to the user as a possible indication that cancer is present in that segment of this voxel.
The advantages of splitting the spectrum into multiple segments is that (a) this provides knowledge of where the high T2 value comes from (i.e., how deep into the tissue) and (b) this effectively narrows the size of the effective voxel, narrowing the diluting effect of healthy tissue (i.e., the partial volume effect that occurs if the voxel contains a mixture of tumor and healthy tissue). The disadvantage of splitting the spectrum is that the signal-to-noise ratio (SNR) of each segment is lower than the full spectrum.

Claims

CLAIMS What is claimed is:
1. A method for extracting multiple T2* values from a single set of Carr-Purcell- Meiboom-Gill (CPMG) data comprising: using a magnetic resonance imaging (MRI) system that includes a radio-frequency (RF) coil to image tissue and generate nuclear magnetic resonance (NMR) time domain signals via a CPMG sequence; converting said signals to a set of spectra, one spectrum per echo, by performing a Fourier transform of complex time domain signals associated with said signals; and dividing each of said spectra into a number of segments and calculating a T2* value for each segment, wherein said T2* values represent decay rates as a function of distance from said RF coil.
2. The method according to claim 1, comprising using equal frequency limits for all echoes in the CPMG sequence.
3. The method according to claim 1, wherein said segments comprise overlapping segments.
PCT/IB2021/052234 2020-03-22 2021-03-17 Method for extracting multiple t2* values from single set of cpmg data WO2021191748A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US202062992994P 2020-03-22 2020-03-22
US62/992,994 2020-03-22

Publications (1)

Publication Number Publication Date
WO2021191748A1 true WO2021191748A1 (en) 2021-09-30

Family

ID=75625610

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IB2021/052234 WO2021191748A1 (en) 2020-03-22 2021-03-17 Method for extracting multiple t2* values from single set of cpmg data

Country Status (1)

Country Link
WO (1) WO2021191748A1 (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2565663A1 (en) * 2010-02-01 2013-03-06 Clear-Cut Medical Ltd Tumor margin assessment of ex-vivo sample
WO2019030620A1 (en) * 2017-08-06 2019-02-14 Clear-Cut Medical Ltd. Hybrid nmr and oct system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2565663A1 (en) * 2010-02-01 2013-03-06 Clear-Cut Medical Ltd Tumor margin assessment of ex-vivo sample
US9310450B2 (en) 2010-02-01 2016-04-12 Clear-Cut Medical Ltd. Margin assessment of ex-vivo sample
WO2019030620A1 (en) * 2017-08-06 2019-02-14 Clear-Cut Medical Ltd. Hybrid nmr and oct system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DANIELI ERNESTO ET AL: "Single-sided magnetic resonance profiling in biological and materials science", JOURNAL OF MAGNETIC RESONANCE, ACADEMIC PRESS, ORLANDO, FL, US, vol. 229, 8 December 2012 (2012-12-08), pages 142 - 154, XP028996800, ISSN: 1090-7807, DOI: 10.1016/J.JMR.2012.11.023 *
LANDEGHEM MAXIME VAN ET AL: "Low-gradient single-sided NMR sensor for one-shot profiling of human skin", JOURNAL OF MAGNETIC RESONANCE, vol. 215, 28 December 2011 (2011-12-28), pages 74 - 84, XP028891294, ISSN: 1090-7807, DOI: 10.1016/J.JMR.2011.12.010 *

Similar Documents

Publication Publication Date Title
US20190046052A1 (en) Distinguishing Diseased Tissue From Healthy Tissue Based On Tissue Component Fractions Using Magnetic Resonance Fingerprinting (MRF)
CN108027409B (en) Time domain MRI
US8605967B2 (en) Magnetic resonance imaging relating to correction of chemical shift artifact and intensity inhomogeneity
US8749236B2 (en) MRI with fat suppression using fat decoupling RF during pre-sequence shimming
US9297873B2 (en) Method to select an undersampling scheme for magnetic resonance imaging, and magnetic resonance imaging method and system using such a selected undersampling scheme
JPH074350B2 (en) NMR signal receiving method
EP3198291B1 (en) Dixon mr imaging with suppression of flow artifacts
EP3081955A1 (en) Mri method for determining signature indices of an observed tissue from signal patterns obtained by motion-probing pulsed gradient mri
JP2015512317A (en) MR imaging using APT contrast enhancement and sampling at multiple echo times
US20080211499A1 (en) Low Power Decoupling for Multi-Nuclear Spectroscopy
US20160313423A1 (en) Mri with dixon-type water/fat separation with estimation of the main magnetic field variations
US10295613B2 (en) Method and apparatus for updating multiple magnetic resonance datasets acquired using a parallel acquisition technique
JP6424170B2 (en) Method for nuclear magnetic resonance analysis
US8224609B2 (en) Excitation uniformity
JP6420348B2 (en) Nuclear magnetic resonance analysis method
WO2018160839A1 (en) System and method for improved diffusion-weighted imaging
Peng et al. Studies of the interactions of an MRI system with the shielding in a combined PET/MRI scanner
US7956611B2 (en) Magnetic resonance imaging apparatus and method
CN111638479B (en) One-dimensional scout image acquisition method and magnetic resonance system
WO2021191748A1 (en) Method for extracting multiple t2* values from single set of cpmg data
EP2388611B1 (en) Production of spectra with improved signal-to-noise ratio
CN111512173B (en) Apparatus and method for diffusion MR imaging with fat suppression
US10302733B2 (en) NMR spin-echo amplitude estimation
Renou et al. Radio-frequency pulse calibration using the MISSTEC sequence
Wörtge et al. Quantitative stray-field T1 relaxometry with the matrix pencil method

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21720558

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21720558

Country of ref document: EP

Kind code of ref document: A1