WO2020190879A1 - Accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging - Google Patents

Accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging Download PDF

Info

Publication number
WO2020190879A1
WO2020190879A1 PCT/US2020/022996 US2020022996W WO2020190879A1 WO 2020190879 A1 WO2020190879 A1 WO 2020190879A1 US 2020022996 W US2020022996 W US 2020022996W WO 2020190879 A1 WO2020190879 A1 WO 2020190879A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
computer
image
contrast
degree
Prior art date
Application number
PCT/US2020/022996
Other languages
French (fr)
Inventor
JR John Thomas VAUGHAN
Sairam Geethanath
Sachin R. JAMBAWALIKAR
Pavan POOJAR
Enlin QIAN
Original Assignee
The Trustees Of Columbia University In The City Of New York
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 The Trustees Of Columbia University In The City Of New York filed Critical The Trustees Of Columbia University In The City Of New York
Priority to CA3133772A priority Critical patent/CA3133772A1/en
Priority to EP20773478.1A priority patent/EP3937775A4/en
Publication of WO2020190879A1 publication Critical patent/WO2020190879A1/en
Priority to US17/475,687 priority patent/US11864863B2/en

Links

Classifications

    • 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/004Features 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 image acquisition of a particular organ or body part
    • A61B5/0042Features 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 image acquisition of a particular organ or body part for the brain
    • 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/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • A61B5/7425Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image
    • 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/5602Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • 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/561Image 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • 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/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging

Definitions

  • the present disclosure relates generally to magnetic resonance imaging, and more specifically, to exemplary embodiments of exemplary system, method and computer- accessible medium for accelerated, simultaneous quantitative and non-synthetic multi- contrast imaging.
  • qMRI quantitative magnetic resonance imaging
  • References 1-5 See, e.g., References 1-5.
  • this involves acquiring data from multiple acquisitions to facilitate a regression analysis. This results in prolonged acquisition times, specifically in the case of multi-parametric magnetic resonance imaging (“MRI”) to obtain maps of relaxation (see, e.g., Reference 6), diffusion (see, e.g., References 7 and 8), pharmacokinetic parameters (see, e.g., Reference 8), etc.
  • MRI multi-parametric magnetic resonance imaging
  • Synthetic MRI in general see, e.g., References 9 and 10
  • MRF Magnetic Fingerprinting
  • Synthetically generated contrast images can be derived from MRF reconstructed parametric maps.
  • it can be challenging to estimate the multitude of phase terms involved in the MR signal equation see, e.g., Reference 13), resulting from diffusion, flow, susceptibility, off-resonance, etc.
  • the exemplary system, method, and computer- accessible medium can be used for simultaneous, natural (e.g non-synthetic), multi-contrast and quantitative MR imaging through tailoring of the MRF acquisition schedule
  • An exemplary system, method, and computer-accessible medium for generating a particular image which can be a quantitative image(s) of at least one section(s) of a patient(s) or (ii) a non-synthetic contrast image(s) of the section(s) of the patient(s), can include, for example, generating a first magnetic resonance (MR) signal and directing the first MR signal to patient(s), receiving a second MR signal from the patient(s) that can be based on the first MR signal, and generating the particular image(s) based on the second MR signal.
  • the first MR signal can be a configured MR signal.
  • the configured MR signal can be configured for a particular contrast.
  • the first MR signal can have a constant signal intensity.
  • the first MR signal can be generated based on a degree of a plurality of flip angles that maintains the constant signal intensity. A degree of flip angles can be selected for the first MR signal based on the particular contrast.
  • the degree of the flip angles can vary within a particular range.
  • the degree of the flip angles can vary about a mean value.
  • the degree of the flip angles can vary monotonously about the mean value.
  • the degree of the flip angles can vary pseudo randomly within the particular range.
  • the particular range can be about - 5 +/- 4 degrees, about 45 +/- 5 degrees, about 75+/- 5 degrees, or about 75 +/- 5 degrees.
  • the particular contrast can include T1, T2, proton density, water, fat, off resonance, diffusion, perfusion, or flow.
  • the non-synthetic contrast image(s) can be a non-synthetic multi-contrast image(s).
  • the particular image(s) can be generated using a reconstruction procedure.
  • the reconstruction procedure can be a sliding window
  • the reconstruction procedure can include converting the second MR signal to an image using a Non-Uniform Fast Fourier Transform.
  • MR information can be generated based on the second MR signal by pre-processing the second MR sign by compensating for a calibrated gradient delay, scaling k-space of the second MR signal with a ratio of a field of view to a matrix size, removing spikes in the second MR signal, and weighting k-space data in the second MR signal with a predetermined density compensation factor, where the particular image(s) can be generated based on the MR information.
  • the particular image(s) can be generated by vector-dot product matching L2-norm normalized dictionary entries with voxel signal evolutions in the second MR signal.
  • Figure 1 A is a set of exemplary graphs illustrating a tailored magnetic resonance fingerprint design for TR, FA and TE values according to an exemplary embodiment of the present disclosure
  • Figure IB is an exemplary graph illustrating Extended Phase Graph simulations of the magnetic resonance fingerprint acquisition schedules from white matter, gray matter, and cerebrospinal fluid according to an exemplary embodiment of the present disclosure
  • Figure 1C is an exemplary graph illustrating Extended Phase Graph simulations for tailored magnetic resonance fingerprint with targeted windows according to an exemplary embodiment of the present disclosure
  • Figure ID is an exemplary graph illustrating three voxels for white matter, gray matter and cerebrospinal fluid from a representative magnetic resonance fingerprint acquired in vivo according to an exemplary embodiment of the present disclosure
  • Figure IE is an exemplary graph illustrating data voxel plots for tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure
  • Figure 2A is a set of exemplary images of a mid-volume slice of four volunteers magnetic resonance fingerprint acquired over 1000 frames according to an exemplary embodiment of the present disclosure
  • Figure 2B is a set of exemplary images of tailored magnetic resonance fingerprint depicting proton density contrast and T1 contrast according to an exemplary embodiment of the present disclosure
  • Figure 2C is a set of exemplary images of contrast windows according to an exemplary embodiment of the present disclosure.
  • Figure 2D is a set of exemplary images of plots for tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure
  • Figure 3A is a set of exemplary images of mid-volume slices of four volunteers magnetic resonance fingerprint acquired data for PDw, T1w, T2w and flow Maximum Intensity Projection (“MIP”) at various contrast window frames according to an exemplary embodiment of the present disclosure
  • Figure 3B is a set of exemplary images of contrast volumes according to an exemplary embodiment of the present disclosure.
  • Figure 3C is a set of exemplary images of tailored magnetic resonance fingerprint contrast windows according to an exemplary embodiment of the present disclosure
  • Figure 3D is a set of exemplary images of contrast volumes from tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure
  • Figure 4A is a set of exemplary images of mid-volume slices of four volunteers magnetic resonance fingerprint acquired data for proton density, T1, T2 maps according to an exemplary embodiment of the present disclosure
  • Figure 4B is a set of volumes depicted by five slices according to an exemplary embodiment of the present disclosure
  • Figure 4C is a set of exemplary tailored magnetic resonance fingerprint parametric maps according to an exemplary embodiment of the present disclosure
  • Figure 4D is a set of exemplary volumes from tailored magnetic resonance fingerprint acquisition according to an exemplary embodiment of the present disclosure
  • Figure 5A is a set of exemplary images of mid-volume slices of four volunteers depicting three regions of interest each for white matter, gray matter and cerebrospinal fluid according to an exemplary embodiment of the present disclosure
  • Figure 5B is an exemplary graph illustrating average signal intensity plots for the regions of interest from Figure 5A according to an exemplary embodiment of the present disclosure
  • Figure 5C is an exemplary graph illustrating tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure.
  • Figure 5D is an exemplary graph illustrating mean and standard deviation for the three regions of interest from Figure 5A for T1 maps according to an exemplary embodiment of the present disclosure
  • Figure 5E is an exemplary graph illustrating mean and standard deviation for the three regions of interest from Figure 5A for T2 maps according to an exemplary embodiment of the present disclosure
  • Figures 6A-6C are exemplary images of slice(s) of a representative brain data set for three contrasts according to an exemplary embodiment of the present disclosure.
  • Figures 7A and 7B are exemplary water-fat images using the exemplary TMRF according to an exemplary embodiment of the present disclosure
  • Figures 8 A and 8B are exemplary quantitative images according to an exemplary embodiment of the present disclosure
  • Figure 8C is an exemplary graph illustrating the comparison of the mean to standard deviation in each of the contrast blocks according to an exemplary embodiment of the present disclosure
  • Figure 9 is a set of exemplary T 1 and T 2 graphs according to an exemplary embodiment of the present disclosure.
  • Figure 10 is a set of exemplary reconstructed images according to an exemplary embodiment of the present disclosure.
  • Figures 11 A-l II are exemplary images comparing image quality according to an exemplary embodiment of the present disclosure.
  • Figure 12 is a flow diagram illustrating an exemplary method for generating a quantitative image of a patient or a non-synthetic contrast image of the patient according to an exemplary embodiment of the present disclosure.
  • Figure 13 is an illustration of an exemplary block diagram of an exemplary system in accordance with certain exemplary embodiments of the present disclosure.
  • a steady state free precession (“SSFP”) sequence can be provided with a thousand spiral readouts designed to derive unique signal evolutions for different tissue matters.
  • MR image contrast derived from such a sequence can be modulated more by the flip angle than the repetition time.
  • Figure 1A shows a set of exemplary graphs illustrating a tailored magnetic resonance fingerprint design for TR, FA and TE values according to an exemplary embodiment of the present disclosure.
  • Figure IB shows an exemplary graph illustrating Extended Phase Graph simulations of the magnetic resonance fingerprint acquisition schedules from white matter, gray matter, and cerebrospinal fluid according to an exemplary embodiment of the present disclosure.
  • MR signal simulation tools such as Extended Phase Graph (“EPG”) (see. e.g., Reference 17)) can be leveraged to vary signal intensities of tissue matters based on their properties (e.g., T 1 , T 2 , etc.) in“blocks” of time. Also, the signal intensities of these tissue matters can remain relatively constant, or slowly varying, if the flip angles were varied smoothly, monotonously, and in a small range within one block.
  • Magnetization preparation in the form of an inversion pulse can be efficiently utilized to include suppression of short and long relaxation components like fat and liquids at different temporal points.
  • This tissue matter dependent tailored design choice can facilitate different relaxation contrast“windows”, signal constancy in a given contrast window, and a meaningful sliding window reconstructed signal intensity image.
  • This exemplary approach referred to as the Tailored MRF (“TMRF”) can retain the MRF derived qMRI benefits of generating multiple parametric maps simultaneously.
  • the MRF sequence includes oscillations in flip angle resulting in corresponding fluctuating signal evolutions. Sliding window reconstruction of non-synthetic signal intensity images can be challenging due to the highly undersampled acquisition and the superposition of highly varying signal intensities.
  • Exemplary Methods [0044] Exemplary Tailored acquisition design and dictionary generation.
  • a one thousand time point (e.g., frames) MRF schedule was used to tailor the magnetization evolution of three tissue types: White Matter (“WM”), Gray Matter (“GM”) and
  • Cerebrospinal Fluid This was performed by designing acquisition blocks targeting one or more contrast windows in each block.
  • three acquisition blocks for (i) PD and T 1 contrast windows; (ii) T 2 plus flow contrast window; and (iii) T 2 contrast window were generated.
  • Each block facilitated a total of 250 frames (e.g. , block size).
  • Each block was generated by choosing different Flip Angles (“FA”) with the minimum TR (e.g., 12.8ms) possible: 5° for PD and T 1 , 45° for T 2 and flow, 70° for T 2 weighting. These values were determined based on SSFP literature (see, e.g., References 13, 15 and 16), and EPG simulations.
  • FA Flip Angles
  • FIG. 1C shows an exemplary graph illustrating Extended Phase Graph simulations for TMRP with targeted windows according to an exemplary embodiment of the present disclosure.
  • a 90° pulse was introduced at the 250 th frame (e.g., end of PD and T 1 acquisition block) to facilitate constant signal intensities during the T 2 plus flow contrast window (e.g., shown in boxes 920 in Figure 1C between 250 th to 500 th frame).
  • a fourth block included repeated values from the first block. This was performed to match the length of the MRF schedule.
  • the values of TE were chosen to be the minimum (e.g., 1.908ms) except in frames between 500 and 625 (e.g., Dixon contrast window). In this contrast window, the TE was increased by 2.2ms to facilitate an out of phase acquisition.
  • the resulting TR, FA and TE vectors for the three blocks were concatenated to form the TR/FA/TE schedule.
  • FIG. 1A illustrates TR, FA and TE values for the 1000 frame acquisition for MRF (e.g., lines 925) and TMRF (e.g., lines 930), the TR and FA plot depict the acquisition“blocks” of length 250 frames.
  • EPG was used to simulate the magnetization evolution dictionary for a range of T 1 (e.g., 0 to 4000ms in steps of 20ms) and T2 (e.g., 0 to 400ms in steps of 20ms; 450 to 600ms in steps of 50ms; 700ms to 2000ms in steps of 500ms) values.
  • T 1 e.g., 0 to 4000ms in steps of 20ms
  • T2 e.g., 0 to 400ms in steps of 20ms; 450 to 600ms in steps of 50ms; 700ms to 2000ms in steps of 500ms
  • the dictionary was then sliding window reconstructed (see, e.g., Reference 18), with a window length of
  • Exemplary MR acquisition In vivo brain imaging of four healthy volunteers was performed. Each of the volunteers was scanned with the MRF and TMRF schedules with 20 slices (e.g., for whole brain coverage), slice thickness of 5mm, minimum TR of 12.6 - 12.8ms, field of view 225 - 240mm in each direction with a final matrix size of 256 x 256. The slice planning for both sequences per volunteer were maintained the same to facilitate spatial comparisons. Both sequences leveraged an 89-shot spiral with 608-point readout, maximum gradient strength of 33mT/m and a maximum slew rate of 120T/m/s on a 3T GE 750w scanner using an eight-channel head coil. The resulting acquisition times for the MRF and TMRF cases were 5: 11 and 4:41 (e.g., minute: seconds) respectively.
  • Exemplary Reconstruction - quantitative maps L2-norm normalized dictionary entries were vector-dot product matched with the voxel signal evolutions of the contrast images. The index of maximum match was chosen as the key to the dictionary entries for the values of T 1 , T 2 . PD was determined as the ratio of the resulting match (e.g., maximum value) to the norm of the chosen entry. Repeating this process for all voxels provided the spatial maps of the quantitative parameters. The maps from both schedules were compared with the MRF maps as the reference. In addition, Regions-Of-Interest (“ROIs”) were drawn on the central slice of each of the four volunteers. Three ROIs each for WM, GM and CSF were drawn and their mean and standard deviation computed.
  • ROIs Regions-Of-Interest
  • ROI analysis Three ROIs each corresponding to WM, GM and CSF were manually drawn on the central slices (#9) of the four volunteers.
  • the ROI masks were multiplied with the thousand-point time series of contrast images to obtain the mean ROI signal evolution of the three tissue matters for each of the four volunteers for both schedules.
  • These ROI masks were also multiplied with the parametric maps to obtain mean T 1 , T 2 and PD values for each volunteer.
  • a mean of means was then computed for the T 1 , T 2 values for MRF and the TMRF schedules for the three tissue matters.
  • TMRF signal evolutions have blocks of slowly varying magnitude compared to MRF, shown by dashed boxes 920.
  • Figure ID illustrates representative sliding window reconstructed signal evolutions from three voxels belonging to three tissue types from a volunteer brain data for the MRF case. Exemplary plots from the corresponding voxel locations for TMRF are shown in Figure IE. The representative voxel data follow the simulation predictions.
  • FIG. 2A The sliding window reconstructed central slice for each of the four volunteers for every 100 frames starting with the 10 th frame is shown in Figure 2A.
  • Figure 2A shows the mid-volume slice of the four volunteers MRF acquired data over 1000 frames with an interval of 100 frames depicting the changes in contrast over time.
  • Figure 2A shows the different, non-synthetic, contrasts that can be derived from TMRF (e.g., frame# 10 - PD; frame # 110 - T 1 weighting with CSF suppression, frame# 610 for T 2 weighting). It also indicates that the signal strengths varied relatively less across the four volunteers in TMRF as compared to MRF.
  • TMRF e.g., frame# 10 - PD; frame # 110 - T 1 weighting with CSF suppression, frame# 610 for T 2 weighting
  • Figure 2B shows corresponding illustrations for TMRF depicting PD contrast at the 10th frame, T1 contrast at the 110th frame, T2 contrast around the 610th frame (e.g., chosen frame for T2 was 575th) and the effect of the 900 pulse on the sliding window reconstruction shown by arrow 205 on frame 210.
  • the exemplary TMRF provides improved contrast for T 1 , PD and T 2 with suppressed flow artifacts (e.g., shown by arrows 210 for the MRF case).
  • the CSF suppression and increased contrast can be particularly noticeable in the case of T 1 Fluid Attenuated Inversion Recovery (“FLAIR”)-like weighting (e.g., the first two frames in Figures 2C and 2D).
  • FLAIR Fluid Attenuated Inversion Recovery
  • Figure 2C shows various contrast windows: five frames of the central slice of the representative MRF data set to depict PD contrast with an interval of ten frames between adjacent images, T1w contrast images of the central slice with an interval of 20 frames, T2w contrast images with an interval of 50 frames, with the arrows 210 depicting unsuppressed CSF and blood flow.
  • the relative signal strength of TMRF based T 2 weighting can be higher than T 1 which in turn can be higher than PD (e.g., shown on the same window level of 0 - 250 au. in Figure 2C).
  • This trend can be validated by the simulated and observed voxel signal evolutions for the three tissue types in the top and bottom rows of Figure 1C.
  • Figure 2D shows corresponding exemplary plots for TMRF data with identical frame intervals within the contrast windows all shown for a signal intensity range between 0 and 250 a.u.
  • the SNR of PD weighted (“PDw”) images in TMRF can be lesser than the MRF case (e.g., top rows in Figures 2C and 2D). This can be in line with the simulated response of the TMRF design.
  • the PDw images in both cases can be normalized to unit amplitude for better visualization, but was not performed to demonstrate validation of simulation results and comparison of SNR for all three contrasts.
  • the T 1 and T 2 weighted images have increased signal strength compared to PD as can be validated by Figures IB and 1C.
  • the relative signal strengths between T 1 and T 2 weighted images do not follow a trend as the signal evolution oscillates by design. This also does not facilitate the simple and efficient use of the sliding window reconstruction.
  • Figure 3A shows a set of exemplary images of mid-volume slices of four volunteers magnetic resonance fingerprint acquired data for PDw, T1w, T2w and flow Maximum Intensity Projection (“MIP”) at various contrast window frames according to an exemplary embodiment of the present disclosure.
  • Figure 3A shows Mid- volume slice of the four volunteers for MRF acquired data for PDw, T1w, T2w and flow MIP at the chosen contrast window frames of 10, 120, 575 and 410.
  • Figure 3B depicts the multiple contrasts obtained from specific temporal points for the central slice of the four volunteers and five equally spaces slices (#1, 5, 9, 13, and 17) out of the twenty slices, with an interval of 4 slices between adjacent images for the
  • FIG. 3C shows TMRF contrast windows in comparison to MRF data shown in Figure 3A with the additional information of water and uncorrected fat images from frames 450 (in phase) and 575 (out of phase). Arrows 305 show the effect of off-resonance on the uncorrected fat images.
  • Figure 3D shows corresponding contrast volumes from TMRF acquisitions for the representative data set.
  • the improved CSF suppression with the exemplary TMRF can be attributed to lower flip angles at the beginning of the schedule following the inversion pulse.
  • the higher SNR for the MIPs can be attributed to the higher flip angles in the second acquisition block leveraging the recovery of long T 2 components like CSF.
  • the MIPs include artifacts, such as capturing information related to eyes and any flow component including CSF rather than blood (e.g., vasculature) alone.
  • TMRF panels also show the results of Dixon imaging.
  • the TMRF fat images suffer from off-resonance artifacts shown arrows 305 in Figures 3C and 3D.
  • FIG. 4A shows mid-volume slices of the four volunteers MRF acquired data for PD, T1, T2 maps.
  • Figure 4B shows corresponding volumes depicted by five slices with an interval of 4 slices between adjacent images for the representative data set.
  • Figure 4C shows TMRF parametric maps in comparison to MRF data shown in Figure 4A with the additional information of uncorrected Proton Density Fat Fraction (“PDFF”) map including with arrow 405 pointing to the off-resonance artifact.
  • PDFF Proton Density Fat Fraction
  • Figure 4D shows corresponding exemplary volumes from TMRF acquisitions for the representative data set.
  • FIG. 5A The ROIs drawn on the four volunteers for WM (e.g., areas 505), GM (e.g., areas 510) and CSF (e.g., areas 515) are shown in Figure 5A.
  • the mean ROI intensity plot over time validated the simulation results shown in Figures 1 A-1E. These curves can also serve as inputs to choose appropriate contrast images in the exemplary TMRF.
  • Figure 5B shows an exemplary graph illustrating average signal intensity plots for the regions of interest (e.g., WM 520, GM 525, and CSF 530) from Figure 5A according to an exemplary embodiment of the present disclosure.
  • the increased CSF suppression in the first and second acquisition blocks can be seen in Figure 5C.
  • the mean ROI values of T 1 for each volunteer is shown as one datum point in Figure 5D, resulting in four data points each for the two schedules.
  • the TMRF T 1 values for GM and WM can be higher compared to corresponding MRF values. However, these values can be within the range of T 1 values reported in literature (e.g., Table 1, (see, e.g.. Reference 21)).
  • the mean of means and standard deviation are also shown to demonstrate the close grouping of these observations for each schedule.
  • Figure 5E shows the corresponding T 2 plots for both schedules.
  • the mean ROI values can be closely grouped across volunteers and corresponding values for the two schedules can be similar to one another.
  • the TMRF schedule can simultaneously produce five different, natural contrast images and four quantitative maps. This can facilitate collapsing a traditional protocol into a sequence, reduce acquisition time for multiple parametric maps, and provide radiologists contrast windows to pick images from.
  • the exemplary TMRF provides a simpler and robust alternative to synthetic MRI approaches that can be used for multiple desired contrasts, measurement based and straightforward to reconstruct and visualize. These acquisitions were included to facilitate a comparison with the thousand frame MRF and subsequent dictionary matching operations. Also, this remainder block can be utilized for other contrast, such as diffusion.
  • the inclusion of a 90° pulse at the beginning of the second acquisition block introduced artifacts in sliding window reconstruction, but was utilized to achieve contrast windows in block 2.
  • TMRF and MRF schedules can benefit from smoothly varying magnetizations or by eliminating reconstructed images at such transients.
  • the relatively slowly changing magnetization evolutions can be leveraged to assign an effective TR, TE and FA to map them to“contrast images” routinely viewed by the radiologist.
  • the TRs and FAs can be generated based on the mean/base value in each acquisition block as shown in Figure 1A.
  • the TMRF schedule can facilitate inclusion of Dixon imaging. These corrections can be leveraged from demonstrated methods. (See, e.g., Reference 24).
  • the choice of FA, TR and TE was based on the knowledge of the tissue matters and contrasts for short TRs.
  • An exemplary EPG look ahead procedure (see, e.g., Reference 17), can be utilized to jointly optimize for TMRF contrast between tissue types over time, rather than echo intensities alone.
  • Figures 6A-6C show exemplary images of slices of a representative brain data set for three contrasts (e.g., MRF illustrated in Figure 6 A, synthetic images from MRF computed maps shown in Figure 6B, and the exemplary TMRF reconstructions illustrated in Figure 6C) according to an exemplary embodiment of the present disclosure.
  • Figure 6C illustrates the choice of the exemplary flip angle in the T1 block that can facilitate suppression of signals from liquids compared to MRF.
  • the growth of CSF and other liquid-like matters can be facilitated in the exemplary TMRF providing a maximum intensity projection image depicting flow.
  • Figure 6A-6C show that the flow effects can be captured well by the exemplary TMRF.
  • TMRF provides for better T1 contrast without flow artifacts (e.g., arrows 605) as compared to the other two contrasts.
  • the PD and T2 images from the exemplary TMRF can also suppress flow artifacts compared to MRF.
  • the window for all three contrasts are maintained identically in order to compare the relative signal strength for each contrast.
  • Figures 7A and 7B show exemplary water-fat images using the exemplary TMRF according to an exemplary embodiment of the present disclosure.
  • the top two rows provided in Figure 7A depict the in and out of phase images for one slice over 5 adjacent time points.
  • Figure 7B shows these exemplary images for three continuous slices to demonstrate uniform separation across slices.
  • the third and fourth rows in Figures 7A and 7B illustrate five frames and the three slices of water and fat images, respectively. These images show the separated fat component along with flow artifacts.
  • the bottom row shows the generated temporal and spatial profiles of the proton density fat fraction map. This information related to water-fat imaging demonstrates the flexibility of the exemplary TMRF to include a desired contrast.
  • Figures 8A and 8B show exemplary quantitative images according to an exemplary embodiment of the present disclosure.
  • Figure 8A shows top map T1 and bottom map T2 over three slices for a representative data set computed using MRF.
  • Figure 8B shows the corresponding maps for the exemplary TMRF. It can be observed that both methods can yield similar maps for T1 and T2, which conforms to previously published literature values.
  • Figure 8C shows an exemplary graph illustrating the comparison of the mean to standard deviation in each of the contrast blocks according to an exemplary embodiment of the present disclosure.
  • the exemplary TMRF provides superior mean to standard deviation for PD and T2. This can also observed in Figures 2A-2D and Figures 6A-6C.
  • T1 contrast from the exemplary TMRF has increased contrast but can have lower mean values, which can be attributed to the decreased signal strength seen from TMRF T1 block.
  • TMRF generated data was validated and optimized in following two steps during the first six months.
  • T 1 and T 2 estimation of the exemplary TMRF was compared to MRF and the Gold Standard (“GS”) measurements at two sites.
  • the GS included Spin Echo based measurements of T 1 and T 2 . All data was acquired on ISMRM/NIST (QalibreMD Inc., CO, USA) which has three layers of spheres with a large range of T 1 , T 2 and Proton Density values. For in vitro studies, data for 10 days was acquired on phantom using three methods (e.g., GS, MRF and TMRF) at two sites with the same vendor and field strength (e.g., on GE 3T Discovery 750w).
  • GS T 1 and T 2 measurements were conducted using inversion recovery spin echo (“IRSE”) and spin echo (“SE”) sequences. The data were reconstructed and matched with an EPG simulated dictionary. ROI analysis was performed to get T 1 and T 2 estimation for 14 spheres. From the exemplary results, the underestimation of TMRF of long T 1 values can be corrected by a constant bias term.
  • site 1 can be close to literature reported values by NIST while site 2 can be close to GS values. The higher order of standard deviation in T 2 values can be due to Bi variations.
  • Figure 9 shows the comparison of T 1 map and T 2 map obtained from MRF and TMRF across two sites. MRF and TMRF data were averaged from 10-day repeatability scans.
  • the data were obtained from the distinct 14 spheres present in the T 1 and T 2 plates of ISMRM/NIST phantom (QalibreMD Inc., CO, USA).
  • a manual ROI analysis was performed on eight datasets. 100 pixels were selected for each sphere in T 1 maps while 70 pixels were selected for T 2 maps due to differences in field of view; MRF and TMRF.
  • Figure 10 shows a set of exemplary reconstructed images according to an exemplary embodiment of the present disclosure for the GS, MRF synthetic contrast, TMRF synthetic contrast, and TMRF natural contrasts.
  • Each row includes four different contrasts - T 1 weighted, T 1 FLAIR, T 2 weighted and Short T1 inversion Recovery (“STIR”).
  • Arrows 1005 and ellipses 1010 in the second and third columns indicate the flow artifacts (e.g ., hyper- intensities or holes) which are not present in the GS and TMRF natural contrast images.
  • the gold standard acquired was magnetization prepared T 1 weighted instead of T 1 weighted. All images were acquired on a 3T GE 750w scanner.
  • MRF and TMRF studies were performed on 4 normal volunteers. Both methods involved an axial, 20 slice brain coverage with a slice thickness of 5mm, and lmm in plane resolution.
  • GS sequences - T 1 weighted, T 1 FLAIR, T 2 weighted and STIR were acquired on the human brain.
  • the spatio-temporal profiles of T 1 , T 2 , PD, water-fat and flow contrasts were reconstructed block-wise along with relaxometric maps.
  • These MRF and TMRF maps were used to generate the synthetic contrast images using the signal equation and were compared.
  • the TMRF natural contrast images were compared with the GS.
  • Figure 10 shows the representative images obtained from GS (e.g., natural contrast), MRF (e.g, synthetic contrast) and TMRF (e.g, synthetic and natural contrasts) for the following sequences - magnetization prepared T 1 weighted, T 1 FLAIR, T 2 weighted and STIR.
  • GS e.g., natural contrast
  • MRF e.g, synthetic contrast
  • TMRF synthetic and natural contrasts
  • Exemplary EOF Analysis Three ROIs each corresponding to WM, GM and CSF were manually drawn on the central slices (#9). The ROI masks were multiplied with the thousand-point time series of contrast images to obtain the mean ROI signal evolution of the three types of tissue matter for each of the volunteers for both schedules. These ROI masks were also multiplied with the parametric maps to obtain mean T 1 , T 2 and PD values for each volunteer. A mean of means was computed for T 1 , T 2 values for MRF and the TMRF schedules for the three types of tissue matter. The TMRF T 1 values for GM and WM can be higher compared to corresponding MRF values. However, these values can be within the range of T 1 values as reported in literature. The mean ROI values can be closely grouped across volunteers and corresponding values for the two schedules can be similar to one another.
  • Figures 1 lA-1 II illustrated exemplary images comparing image quality according to an exemplary embodiment of the present disclosure.
  • the synthetically generated images e.g., MRF and TMRF
  • the synthetically generated images can be compared with the natural contrast images obtained from GS (see. e.g., Figures 1 lA-1 ID) and the TMRF images produced after sliding window reconstruction as shown in Figures 1 IB-1 ID.
  • the magnified version of the edges (e.g., between white matter and gray matter) and the presence/absence of flow artifact are shown in Figures 11E-1I.
  • the exemplary GS images can be similar to that of natural contrast obtained from TMRF.
  • the data shown here can be representative T 1 FLAIR contrast and the results can be similar to other contrast.
  • the exemplary DIXON images obtained from the natural TMRF data do not suppress the water/fat accurately.
  • the exemplary system, method and computer-accessible medium can optimize the DIXON method for better suppression of the fat and water. These optimized images can be compared with the LAVA Flex sequence.
  • FIG. 11 A The comparison of the image quality between the GS is shown in Figure 11 A, between exemplary synthetic images generated from MRF is shown in Figure G1B, between synthetic images from TMRF is shown in Figure 1 1C, and between natural images obtained after sliding window reconstruction from TMRF is shown in Figure 1 1C.
  • Ail illustrated images are T 1 FLAIR contrast.
  • a part of the image is magnified (e.g., as shown in the boxes) to see the effect of flow as shown in Figures 11E-1 1G and boundaries between gray matter and white matter, as shown in Figures 11H and 11I.
  • Rapid comprehensive MRI exams for example, pediatric neuroimaging, multi- parametric prostate imaging, whole body imaging oncology, diabetes studies inclusive of NASH, study of fat types such as brown, white and brite fat, etc.
  • MR value driven protocols for example, 5 -minute stroke protocol, as an alternate to EPImix, MAGIC, etc.
  • Multi-scale, multi-modality image fusion e.g., rapid MR-PET exams for
  • Figure 12 shows a flow diagram illustrating an exemplary method for generating a quantitative image of a patient or a non-synthetic contrast image of the patient according to an exemplary embodiment of the present disclosure.
  • a degree of flip angles for the first MR signal can be selected based on the particular contrast.
  • a first magnetic resonance (MR) signal can be generated, which can be directed to the patient.
  • a second MR signal can be received from the patient that can be based on the first MR signal.
  • the second MR signal can be pre-processed by compensating for a calibrated gradient delay.
  • k- space of the second MR signal can be scaled with a ratio of a field of view to a matrix size.
  • spikes in the second MR signal can be removed.
  • k- space data in the second MR signal can be weighted with a predetermined density compensation factor.
  • the particular image can be generated based on the second MR signal, for example, using a reconstruction procedure.
  • Figure 13 shows a block diagram of an exemplary embodiment of a system according to the present disclosure.
  • exemplary procedures in accordance with the present disclosure described herein can be performed by a processing arrangement and/or a computing arrangement (e.g., computer hardware arrangement) 1305.
  • a processing arrangement and/or a computing arrangement e.g., computer hardware arrangement
  • processing/computing arrangement 1305 can be, for example entirely or a part of, or include, but not limited to, a computer/processor 1310 that can include, for example one or more microprocessors, and use instructions stored on a computer-accessible medium (e.g., RAM, ROM, hard drive, or other storage device).
  • a computer-accessible medium e.g., RAM, ROM, hard drive, or other storage device.
  • a computer-accessible medium 1315 e.g., as described herein above, a storage device such as a hard disk, floppy disk, memory stick, CD- ROM, RAM, ROM, etc., or a collection thereol
  • the computer-accessible medium 1315 can contain executable instructions 1320 thereon.
  • a storage arrangement 1325 can be provided separately from the computer-accessible medium 1315, which can provide the instructions to the processing arrangement 1305 so as to configure the processing arrangement to execute certain exemplary procedures, processes, and methods, as described herein above, for example.
  • the exemplary processing arrangement 1305 can be provided with or include an input/output ports 1335, which can include, for example a wired network, a wireless network, the internet, an intranet, a data collection probe, a sensor, etc.
  • the exemplary processing arrangement 1305 can be in communication with an exemplary display arrangement 1330, which, according to certain exemplary embodiments of the present disclosure, can be a touch-screen configured for inputting information to the processing arrangement in addition to outputting information from the processing arrangement, for example.
  • the exemplary display arrangement 1330 and/or a storage arrangement 1325 can be used to display and/or store data in a user-accessible format and/or user-readable format.
  • Weigel M Extended phase graphs: Dephasing, RF pulses, and echoes - pure and simple.

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Signal Processing (AREA)
  • Neurology (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

An exemplary system, method, and computer-accessible medium for generating a particular image which can be a quantitative image(s) of at least one section(s) of a patient(s) or (ii) a non-synthetic contrast image(s) of the section(s) of the patient(s), can include, for example, generating a first magnetic resonance (MR) signal and directing the first MR signal to patient(s), receiving a second MR signal from the patient(s) that can be based on the first MR signal, and generating the particular image(s) based on the second MR signal. The first MR signal can be a configured MR signal. The configured MR signal can be configured for a particular contrast. The first MR signal can have a constant signal intensity. The first MR signal can be generated based on a degree of a plurality of flip angles that maintains the constant signal intensity. A degree of flip angles can be selected for the first MR signal based on the particular contrast

Description

ACCELERATED, SIMULTANEOUS QUANTITATIVE AND NON-SYNTHETIC
MULTI-CONTRAST IMAGING
CROSS-REFERENCE TO RELATED APPLICATION(S)
[0001] This application relates to and claims priority from U.S. Patent Application No. 62/819,159, filed on March 15, 2019, the entire disclosure of which is incorporated herein by reference.
FIELD OF THE DISCLOSURE
[0002] The present disclosure relates generally to magnetic resonance imaging, and more specifically, to exemplary embodiments of exemplary system, method and computer- accessible medium for accelerated, simultaneous quantitative and non-synthetic multi- contrast imaging.
BACKGROUND INFORMATION
[0003] The role of quantitative magnetic resonance imaging (“qMRI”) is well established. (See, e.g., References 1-5). Typically, this involves acquiring data from multiple acquisitions to facilitate a regression analysis. This results in prolonged acquisition times, specifically in the case of multi-parametric magnetic resonance imaging (“MRI”) to obtain maps of relaxation (see, e.g., Reference 6), diffusion (see, e.g., References 7 and 8), pharmacokinetic parameters (see, e.g., Reference 8), etc. Synthetic MRI in general (see, e.g., References 9 and 10), and MR Fingerprinting (“MRF”) (see, e.g., References 11 and 12), in particular mitigate this challenge by simultaneously acquiring data for multiple parametric maps (e.g., T1, T2, etc.). However, one issue can be the unavailability of contrast images routinely obtained from conventional imaging protocols. Synthetically generated contrast images can be derived from MRF reconstructed parametric maps. However, it can be challenging to estimate the multitude of phase terms involved in the MR signal equation (see, e.g., Reference 13), resulting from diffusion, flow, susceptibility, off-resonance, etc. Recently, a deep learning based MRF reconstruction approach directly reconstructing contrast images from MRF data has been reported. (See, e.g., Reference 14). This work overcomes challenges with simulation model limitations and associated artifacts. It assumes that the routinely obtained training data can be the ground truth. However, these images can also be“synthetically” generated based on MRF data and training data, rather than from raw data directly, and have not been explored in the cases of pathology. Thus, the performance of such methods can significantly depend on the variety, volume and veracity of training data.
[0004] Thus, it may be beneficial to provide the exemplary system, method and computer- accessible medium for an accelerated, simultaneous quantitative and non-synthetic multi- contrast imaging which can overcome at least some of the deficiencies described herein above.
SUMMARY OF EXEMPLARY EMBODIMENTS
[0005] To overcome these limitations, the exemplary system, method, and computer- accessible medium, according to an exemplary embodiment of the present disclosure, can be used for simultaneous, natural ( e.g non-synthetic), multi-contrast and quantitative MR imaging through tailoring of the MRF acquisition schedule
[0006] An exemplary system, method, and computer-accessible medium for generating a particular image which can be a quantitative image(s) of at least one section(s) of a patient(s) or (ii) a non-synthetic contrast image(s) of the section(s) of the patient(s), can include, for example, generating a first magnetic resonance (MR) signal and directing the first MR signal to patient(s), receiving a second MR signal from the patient(s) that can be based on the first MR signal, and generating the particular image(s) based on the second MR signal. The first MR signal can be a configured MR signal. The configured MR signal can be configured for a particular contrast. The first MR signal can have a constant signal intensity. The first MR signal can be generated based on a degree of a plurality of flip angles that maintains the constant signal intensity. A degree of flip angles can be selected for the first MR signal based on the particular contrast.
[0007] In some exemplary embodiments of the present disclosure, the degree of the flip angles can vary within a particular range. The degree of the flip angles can vary about a mean value. The degree of the flip angles can vary monotonously about the mean value. The degree of the flip angles can vary pseudo randomly within the particular range. The particular range can be about - 5 +/- 4 degrees, about 45 +/- 5 degrees, about 75+/- 5 degrees, or about 75 +/- 5 degrees. The particular contrast can include T1, T2, proton density, water, fat, off resonance, diffusion, perfusion, or flow. The non-synthetic contrast image(s) can be a non-synthetic multi-contrast image(s). The particular image(s) can be generated using a reconstruction procedure. The reconstruction procedure can be a sliding window
reconstruction procedure. The reconstruction procedure can include converting the second MR signal to an image using a Non-Uniform Fast Fourier Transform.
[0008] In certain exemplary embodiments of the present disclosure, MR information can be generated based on the second MR signal by pre-processing the second MR sign by compensating for a calibrated gradient delay, scaling k-space of the second MR signal with a ratio of a field of view to a matrix size, removing spikes in the second MR signal, and weighting k-space data in the second MR signal with a predetermined density compensation factor, where the particular image(s) can be generated based on the MR information. The particular image(s) can be generated by vector-dot product matching L2-norm normalized dictionary entries with voxel signal evolutions in the second MR signal. [0009] These and other objects, features and advantages of the exemplary embodiments of the present disclosure will become apparent upon reading the following detailed description of the exemplary embodiments of the present disclosure, when taken in conjunction with the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
[0010] Further objects, features and advantages of the present disclosure will become apparent from the following detailed description taken in conjunction with the accompanying Figures showing illustrative embodiments of the present disclosure, in which:
[0011] Figure 1 A is a set of exemplary graphs illustrating a tailored magnetic resonance fingerprint design for TR, FA and TE values according to an exemplary embodiment of the present disclosure;
[0012] Figure IB is an exemplary graph illustrating Extended Phase Graph simulations of the magnetic resonance fingerprint acquisition schedules from white matter, gray matter, and cerebrospinal fluid according to an exemplary embodiment of the present disclosure;
[0013] Figure 1C is an exemplary graph illustrating Extended Phase Graph simulations for tailored magnetic resonance fingerprint with targeted windows according to an exemplary embodiment of the present disclosure;
[0014] Figure ID is an exemplary graph illustrating three voxels for white matter, gray matter and cerebrospinal fluid from a representative magnetic resonance fingerprint acquired in vivo according to an exemplary embodiment of the present disclosure;
[0015] Figure IE is an exemplary graph illustrating data voxel plots for tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure; [0016] Figure 2A is a set of exemplary images of a mid-volume slice of four volunteers magnetic resonance fingerprint acquired over 1000 frames according to an exemplary embodiment of the present disclosure;
[0017] Figure 2B is a set of exemplary images of tailored magnetic resonance fingerprint depicting proton density contrast and T1 contrast according to an exemplary embodiment of the present disclosure;
[0018] Figure 2C is a set of exemplary images of contrast windows according to an exemplary embodiment of the present disclosure;
[0019] Figure 2D is a set of exemplary images of plots for tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure;
[0020] Figure 3A is a set of exemplary images of mid-volume slices of four volunteers magnetic resonance fingerprint acquired data for PDw, T1w, T2w and flow Maximum Intensity Projection (“MIP”) at various contrast window frames according to an exemplary embodiment of the present disclosure;
[0021] Figure 3B is a set of exemplary images of contrast volumes according to an exemplary embodiment of the present disclosure;
[0022] Figure 3C is a set of exemplary images of tailored magnetic resonance fingerprint contrast windows according to an exemplary embodiment of the present disclosure;
[0023] Figure 3D is a set of exemplary images of contrast volumes from tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure;
[0024] Figure 4A is a set of exemplary images of mid-volume slices of four volunteers magnetic resonance fingerprint acquired data for proton density, T1, T2 maps according to an exemplary embodiment of the present disclosure;
[0025] Figure 4B is a set of volumes depicted by five slices according to an exemplary embodiment of the present disclosure; [0026] Figure 4C is a set of exemplary tailored magnetic resonance fingerprint parametric maps according to an exemplary embodiment of the present disclosure;
[0027] Figure 4D is a set of exemplary volumes from tailored magnetic resonance fingerprint acquisition according to an exemplary embodiment of the present disclosure;
[0028] Figure 5A is a set of exemplary images of mid-volume slices of four volunteers depicting three regions of interest each for white matter, gray matter and cerebrospinal fluid according to an exemplary embodiment of the present disclosure;
[0029] Figure 5B is an exemplary graph illustrating average signal intensity plots for the regions of interest from Figure 5A according to an exemplary embodiment of the present disclosure;
[0030] Figure 5C is an exemplary graph illustrating tailored magnetic resonance fingerprint according to an exemplary embodiment of the present disclosure;
[0031] Figure 5D is an exemplary graph illustrating mean and standard deviation for the three regions of interest from Figure 5A for T1 maps according to an exemplary embodiment of the present disclosure;
[0032] Figure 5E is an exemplary graph illustrating mean and standard deviation for the three regions of interest from Figure 5A for T2 maps according to an exemplary embodiment of the present disclosure;
[0033] Figures 6A-6C are exemplary images of slice(s) of a representative brain data set for three contrasts according to an exemplary embodiment of the present disclosure.
[0034] Figures 7A and 7B are exemplary water-fat images using the exemplary TMRF according to an exemplary embodiment of the present disclosure;
[0035] Figures 8 A and 8B are exemplary quantitative images according to an exemplary embodiment of the present disclosure; [0036] Figure 8C is an exemplary graph illustrating the comparison of the mean to standard deviation in each of the contrast blocks according to an exemplary embodiment of the present disclosure;
[0037] Figure 9 is a set of exemplary T1 and T2 graphs according to an exemplary embodiment of the present disclosure;
[0038] Figure 10 is a set of exemplary reconstructed images according to an exemplary embodiment of the present disclosure;
[0039] Figures 11 A-l II are exemplary images comparing image quality according to an exemplary embodiment of the present disclosure;
[0040] Figure 12 is a flow diagram illustrating an exemplary method for generating a quantitative image of a patient or a non-synthetic contrast image of the patient according to an exemplary embodiment of the present disclosure; and
[0041] Figure 13 is an illustration of an exemplary block diagram of an exemplary system in accordance with certain exemplary embodiments of the present disclosure.
[0042] Throughout the drawings, the same reference numerals and characters, unless otherwise stated, are used to denote like features, elements, components or portions of the illustrated embodiments. Moreover, while the present disclosure will now be described in detail with reference to the figures, it is done so in connection with the illustrative embodiments and is not limited by the particular embodiments illustrated in the figures and the appended claims.
DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS
[0043] A steady state free precession (“SSFP”) sequence can be provided with a thousand spiral readouts designed to derive unique signal evolutions for different tissue matters. MR image contrast derived from such a sequence can be modulated more by the flip angle than the repetition time. (See, e.g.. References 13, 15 and 16, and Figures 1A and IB). Figure 1A shows a set of exemplary graphs illustrating a tailored magnetic resonance fingerprint design for TR, FA and TE values according to an exemplary embodiment of the present disclosure. Figure IB shows an exemplary graph illustrating Extended Phase Graph simulations of the magnetic resonance fingerprint acquisition schedules from white matter, gray matter, and cerebrospinal fluid according to an exemplary embodiment of the present disclosure. This has been shown by MRF to significantly vary flip angles while restricting the TR to a much smaller range above the minimum TR achievable. Further, MR signal simulation tools such as Extended Phase Graph (“EPG”) (see. e.g., Reference 17)) can be leveraged to vary signal intensities of tissue matters based on their properties (e.g., T1, T2, etc.) in“blocks” of time. Also, the signal intensities of these tissue matters can remain relatively constant, or slowly varying, if the flip angles were varied smoothly, monotonously, and in a small range within one block. Magnetization preparation in the form of an inversion pulse can be efficiently utilized to include suppression of short and long relaxation components like fat and liquids at different temporal points. This tissue matter dependent tailored design choice can facilitate different relaxation contrast“windows”, signal constancy in a given contrast window, and a meaningful sliding window reconstructed signal intensity image. This exemplary approach, referred to as the Tailored MRF (“TMRF”) can retain the MRF derived qMRI benefits of generating multiple parametric maps simultaneously. In comparison, the MRF sequence includes oscillations in flip angle resulting in corresponding fluctuating signal evolutions. Sliding window reconstruction of non-synthetic signal intensity images can be challenging due to the highly undersampled acquisition and the superposition of highly varying signal intensities.
Exemplary Methods [0044] Exemplary Tailored acquisition design and dictionary generation. A one thousand time point (e.g., frames) MRF schedule was used to tailor the magnetization evolution of three tissue types: White Matter (“WM”), Gray Matter (“GM”) and
Cerebrospinal Fluid (“CSF”). This was performed by designing acquisition blocks targeting one or more contrast windows in each block. In this exemplary implementation, three acquisition blocks for (i) PD and T1 contrast windows; (ii) T2 plus flow contrast window; and (iii) T2 contrast window were generated. Each block facilitated a total of 250 frames (e.g. , block size). Each block was generated by choosing different Flip Angles (“FA”) with the minimum TR (e.g., 12.8ms) possible: 5° for PD and T1, 45° for T2 and flow, 70° for T2 weighting. These values were determined based on SSFP literature (see, e.g., References 13, 15 and 16), and EPG simulations. These FA values can be referred to as base FAs for each of the contrast blocks. In each of the three cases, a normally distributed noise with a zero mean and a standard deviation of 0.5 was added to the minimum TR and the base FAs. These two vectors were then sorted in ascending order to avoid spike like transients in the magnetization evolution. Figure 1C shows an exemplary graph illustrating Extended Phase Graph simulations for TMRP with targeted windows according to an exemplary embodiment of the present disclosure. A 90° pulse was introduced at the 250th frame (e.g., end of PD and T1 acquisition block) to facilitate constant signal intensities during the T2 plus flow contrast window (e.g., shown in boxes 920 in Figure 1C between 250th to 500th frame). A fourth block (e.g., Frame# 751 to 1000) included repeated values from the first block. This was performed to match the length of the MRF schedule. The values of TE were chosen to be the minimum (e.g., 1.908ms) except in frames between 500 and 625 (e.g., Dixon contrast window). In this contrast window, the TE was increased by 2.2ms to facilitate an out of phase acquisition. The resulting TR, FA and TE vectors for the three blocks were concatenated to form the TR/FA/TE schedule. These TMRF parameters in comparison to the corresponding MRF schedule (see, e.g.. Reference 12), are shown in Figure 1A, which illustrates TR, FA and TE values for the 1000 frame acquisition for MRF (e.g., lines 925) and TMRF (e.g., lines 930), the TR and FA plot depict the acquisition“blocks” of length 250 frames. EPG was used to simulate the magnetization evolution dictionary for a range of T1 (e.g., 0 to 4000ms in steps of 20ms) and T2 (e.g., 0 to 400ms in steps of 20ms; 450 to 600ms in steps of 50ms; 700ms to 2000ms in steps of 500ms) values. The dictionary was then sliding window reconstructed (see, e.g., Reference 18), with a window length of 89. The simulation was repeated to account for any changes in the minimum TR resulting from the acquisition. Examples of flip angles can include, but is not limited to, to the following: T1 - 5 +/- 4 degrees (plus or minus 15%), T2 = 45 +/- 5 degrees (plus or minus 15%), T1 Flair = 75+/- 5 degrees (plus or minus 15%), and Dixon = 75 +/- 5 degrees (plus or minus 15%).
[0045] Exemplary MR acquisition . In vivo brain imaging of four healthy volunteers was performed. Each of the volunteers was scanned with the MRF and TMRF schedules with 20 slices (e.g., for whole brain coverage), slice thickness of 5mm, minimum TR of 12.6 - 12.8ms, field of view 225 - 240mm in each direction with a final matrix size of 256 x 256. The slice planning for both sequences per volunteer were maintained the same to facilitate spatial comparisons. Both sequences leveraged an 89-shot spiral with 608-point readout, maximum gradient strength of 33mT/m and a maximum slew rate of 120T/m/s on a 3T GE 750w scanner using an eight-channel head coil. The resulting acquisition times for the MRF and TMRF cases were 5: 11 and 4:41 (e.g., minute: seconds) respectively.
[0046] Exemplary Reconstruction - non-synthetic contrast images. Raw data for both cases were pre-processed by compensating for a calibrated gradient delay (e.g.,
3.5 microseconds in the exemplary case), scaling k-space with the ratio of FOV to matrix size, removal of spikes (e.g., threshold of twice the standard deviation of the FID) and weighting the k-space data with the pre-computed density compensation factor for the 89 shot spiral trajectory. This data was then converted to images using the Non-Uniform Fast Fourier Transform using the image reconstruction toolbox (see, e.g.. Reference 19), and complex coil combined and sliding window reconstructed to provide 2D multi-slice images over time.
[0047] Exemplary Reconstruction - quantitative maps. L2-norm normalized dictionary entries were vector-dot product matched with the voxel signal evolutions of the contrast images. The index of maximum match was chosen as the key to the dictionary entries for the values of T1, T2. PD was determined as the ratio of the resulting match (e.g., maximum value) to the norm of the chosen entry. Repeating this process for all voxels provided the spatial maps of the quantitative parameters. The maps from both schedules were compared with the MRF maps as the reference. In addition, Regions-Of-Interest (“ROIs”) were drawn on the central slice of each of the four volunteers. Three ROIs each for WM, GM and CSF were drawn and their mean and standard deviation computed. For the PDFF maps in the case of TMRF, the In Phase (“IP”) and Out of Phase (“OP”) images were added and subtracted to yield water-only (“W”) and fat-only (“F”) images as typically performed in Dixon imaging. (See, e.g., Reference 20). A PDFF map was computed as F/(F + W).
[0048] Exemplary ROI analysis. Three ROIs each corresponding to WM, GM and CSF were manually drawn on the central slices (#9) of the four volunteers. The ROI masks were multiplied with the thousand-point time series of contrast images to obtain the mean ROI signal evolution of the three tissue matters for each of the four volunteers for both schedules. These ROI masks were also multiplied with the parametric maps to obtain mean T1, T2 and PD values for each volunteer. A mean of means was then computed for the T1, T2 values for MRF and the TMRF schedules for the three tissue matters.
Exemplary Results [0049] Exemplary Acquisition schedule design and simulation. Figures IB and 1C show representative MRF and TMRF sliding window simulated signal evolutions for the three tissue types of WM shown by line 905 (T1=850ms, T2=80ms for simulation), GM shown by line 910 (T1=1330ms, T2=110ms) and CSF shown by line 915 (T1=4500ms, T2=1700ms). TMRF signal evolutions have blocks of slowly varying magnitude compared to MRF, shown by dashed boxes 920. Figure ID illustrates representative sliding window reconstructed signal evolutions from three voxels belonging to three tissue types from a volunteer brain data for the MRF case. Exemplary plots from the corresponding voxel locations for TMRF are shown in Figure IE. The representative voxel data follow the simulation predictions.
The reduction in variation of the TMRF signal evolutions can be easily observed.
[0050] Exemplary Non-synthetic MRI reconstruction - temporal and spatial profiles : The sliding window reconstructed central slice for each of the four volunteers for every 100 frames starting with the 10th frame is shown in Figure 2A. In particular, Figure 2A shows the mid-volume slice of the four volunteers MRF acquired data over 1000 frames with an interval of 100 frames depicting the changes in contrast over time. Figure 2A shows the different, non-synthetic, contrasts that can be derived from TMRF (e.g., frame# 10 - PD; frame # 110 - T1 weighting with CSF suppression, frame# 610 for T2 weighting). It also indicates that the signal strengths varied relatively less across the four volunteers in TMRF as compared to MRF. This can be attributed to the slowly and monotonously increasing values of the TMRF schedule. In the case of the sudden variation at the end of block 1 (e.g., due to the 90° pulse at the 250th frame), an aliasing-like artifact was caused shown arrows 210. Given that this frame was not part of the targeted contrast blocks, this artifact was ignored. Furthermore, the illustrations of Figures 2B shows five adjacent frames for the PD, T1 and T2 contrasts for both acquisition procedures for a representative data set. In particular, Figure 2B shows corresponding illustrations for TMRF depicting PD contrast at the 10th frame, T1 contrast at the 110th frame, T2 contrast around the 610th frame (e.g., chosen frame for T2 was 575th) and the effect of the 900 pulse on the sliding window reconstruction shown by arrow 205 on frame 210.
[0051] Qualitatively, the exemplary TMRF provides improved contrast for T1, PD and T2 with suppressed flow artifacts (e.g., shown by arrows 210 for the MRF case). The CSF suppression and increased contrast can be particularly noticeable in the case of T1 Fluid Attenuated Inversion Recovery (“FLAIR”)-like weighting (e.g., the first two frames in Figures 2C and 2D). For example, Figure 2C shows various contrast windows: five frames of the central slice of the representative MRF data set to depict PD contrast with an interval of ten frames between adjacent images, T1w contrast images of the central slice with an interval of 20 frames, T2w contrast images with an interval of 50 frames, with the arrows 210 depicting unsuppressed CSF and blood flow. The relative signal strength of TMRF based T2 weighting can be higher than T1 which in turn can be higher than PD (e.g., shown on the same window level of 0 - 250 au. in Figure 2C). This trend can be validated by the simulated and observed voxel signal evolutions for the three tissue types in the top and bottom rows of Figure 1C. Figure 2D shows corresponding exemplary plots for TMRF data with identical frame intervals within the contrast windows all shown for a signal intensity range between 0 and 250 a.u.
The SNR of PD weighted (“PDw”) images in TMRF can be lesser than the MRF case (e.g., top rows in Figures 2C and 2D). This can be in line with the simulated response of the TMRF design. The PDw images in both cases can be normalized to unit amplitude for better visualization, but was not performed to demonstrate validation of simulation results and comparison of SNR for all three contrasts. In the case of MRF, the T1 and T2 weighted images have increased signal strength compared to PD as can be validated by Figures IB and 1C. However, the relative signal strengths between T1 and T2 weighted images do not follow a trend as the signal evolution oscillates by design. This also does not facilitate the simple and efficient use of the sliding window reconstruction.
[0052] Figure 3A shows a set of exemplary images of mid-volume slices of four volunteers magnetic resonance fingerprint acquired data for PDw, T1w, T2w and flow Maximum Intensity Projection (“MIP”) at various contrast window frames according to an exemplary embodiment of the present disclosure. In particular, Figure 3A shows Mid- volume slice of the four volunteers for MRF acquired data for PDw, T1w, T2w and flow MIP at the chosen contrast window frames of 10, 120, 575 and 410.
[0053] Figure 3B depicts the multiple contrasts obtained from specific temporal points for the central slice of the four volunteers and five equally spaces slices (#1, 5, 9, 13, and 17) out of the twenty slices, with an interval of 4 slices between adjacent images for the
representative data set, resulting in full brain coverage for the representative data set. The corresponding TMRF panels shown in Figures 3C and 3D qualitatively show improved T1 ( e.g ., WM > GM) and T2 (e.g., GM > WM) contrast, better CSF suppression for T1 weighted images and higher SNR for flow MIP. For example, Figure 3C shows TMRF contrast windows in comparison to MRF data shown in Figure 3A with the additional information of water and uncorrected fat images from frames 450 (in phase) and 575 (out of phase). Arrows 305 show the effect of off-resonance on the uncorrected fat images. Figure 3D shows corresponding contrast volumes from TMRF acquisitions for the representative data set.
[0054] The improved CSF suppression with the exemplary TMRF can be attributed to lower flip angles at the beginning of the schedule following the inversion pulse. The higher SNR for the MIPs can be attributed to the higher flip angles in the second acquisition block leveraging the recovery of long T2 components like CSF. The MIPs include artifacts, such as capturing information related to eyes and any flow component including CSF rather than blood (e.g., vasculature) alone. TMRF panels also show the results of Dixon imaging. The TMRF fat images suffer from off-resonance artifacts shown arrows 305 in Figures 3C and 3D.
[0055] Exemplary Reconstruction - quantitative maps. Dictionary matched T1, T2 and PD maps for MRF and TMRF schedules are shown in Figures 4A and 4B for the central slice for the four volunteers. For example, Figure 4A shows mid-volume slices of the four volunteers MRF acquired data for PD, T1, T2 maps. Figure 4B shows corresponding volumes depicted by five slices with an interval of 4 slices between adjacent images for the representative data set. Figure 4C shows TMRF parametric maps in comparison to MRF data shown in Figure 4A with the additional information of uncorrected Proton Density Fat Fraction (“PDFF”) map including with arrow 405 pointing to the off-resonance artifact.
Figure 4D shows corresponding exemplary volumes from TMRF acquisitions for the representative data set.
[0056] Parametric maps of the five slices for the representative data set shown in Figures 2B and 2D are shown in Figures 3B and 3D. It can be observed that both approaches result in similar maps and the range of values (e.g., for WM, GM) conform to the literature. (See. e.g.. Reference 21). This can be validated qualitatively by the window range and intensities shown. The value of CSF T2 can be much shorter than the literature value and follows values reported in MRF literature. (See, e.g., References 11 and 12). These values can be attributed to flow artifacts and the absence of a longer TR.
[0057] Exemplary ROI analysis. The ROIs drawn on the four volunteers for WM (e.g., areas 505), GM (e.g., areas 510) and CSF (e.g., areas 515) are shown in Figure 5A. The mean ROI intensity plot over time validated the simulation results shown in Figures 1 A-1E. These curves can also serve as inputs to choose appropriate contrast images in the exemplary TMRF. Figure 5B shows an exemplary graph illustrating average signal intensity plots for the regions of interest (e.g., WM 520, GM 525, and CSF 530) from Figure 5A according to an exemplary embodiment of the present disclosure. The increased CSF suppression in the first and second acquisition blocks can be seen in Figure 5C. The mean ROI values of T1 for each volunteer is shown as one datum point in Figure 5D, resulting in four data points each for the two schedules. The TMRF T1 values for GM and WM can be higher compared to corresponding MRF values. However, these values can be within the range of T1 values reported in literature (e.g., Table 1, (see, e.g.. Reference 21)). The mean of means and standard deviation are also shown to demonstrate the close grouping of these observations for each schedule. Figure 5E shows the corresponding T2 plots for both schedules. The mean ROI values can be closely grouped across volunteers and corresponding values for the two schedules can be similar to one another.
[0058] The TMRF schedule can simultaneously produce five different, natural contrast images and four quantitative maps. This can facilitate collapsing a traditional protocol into a sequence, reduce acquisition time for multiple parametric maps, and provide radiologists contrast windows to pick images from. The exemplary TMRF provides a simpler and robust alternative to synthetic MRI approaches that can be used for multiple desired contrasts, measurement based and straightforward to reconstruct and visualize. These acquisitions were included to facilitate a comparison with the thousand frame MRF and subsequent dictionary matching operations. Also, this remainder block can be utilized for other contrast, such as diffusion. The inclusion of a 90° pulse at the beginning of the second acquisition block introduced artifacts in sliding window reconstruction, but was utilized to achieve contrast windows in block 2. Therefore, TMRF and MRF schedules can benefit from smoothly varying magnetizations or by eliminating reconstructed images at such transients. The relatively slowly changing magnetization evolutions can be leveraged to assign an effective TR, TE and FA to map them to“contrast images” routinely viewed by the radiologist. For TMRF, the TRs and FAs can be generated based on the mean/base value in each acquisition block as shown in Figure 1A. The TMRF schedule can facilitate inclusion of Dixon imaging. These corrections can be leveraged from demonstrated methods. (See, e.g., Reference 24). The choice of FA, TR and TE was based on the knowledge of the tissue matters and contrasts for short TRs. An exemplary EPG look ahead procedure (see, e.g., Reference 17), can be utilized to jointly optimize for TMRF contrast between tissue types over time, rather than echo intensities alone.
[0059] Figures 6A-6C show exemplary images of slices of a representative brain data set for three contrasts (e.g., MRF illustrated in Figure 6 A, synthetic images from MRF computed maps shown in Figure 6B, and the exemplary TMRF reconstructions illustrated in Figure 6C) according to an exemplary embodiment of the present disclosure. Figure 6C illustrates the choice of the exemplary flip angle in the T1 block that can facilitate suppression of signals from liquids compared to MRF. The growth of CSF and other liquid-like matters can be facilitated in the exemplary TMRF providing a maximum intensity projection image depicting flow. Additionally, Figure 6A-6C show that the flow effects can be captured well by the exemplary TMRF. It can be observed that TMRF provides for better T1 contrast without flow artifacts (e.g., arrows 605) as compared to the other two contrasts. The PD and T2 images from the exemplary TMRF can also suppress flow artifacts compared to MRF.
The window for all three contrasts are maintained identically in order to compare the relative signal strength for each contrast.
Exemplary Water-Fat Imaging Using TMRF
[0060] Figures 7A and 7B show exemplary water-fat images using the exemplary TMRF according to an exemplary embodiment of the present disclosure. The top two rows provided in Figure 7A depict the in and out of phase images for one slice over 5 adjacent time points. Figure 7B shows these exemplary images for three continuous slices to demonstrate uniform separation across slices. The third and fourth rows in Figures 7A and 7B illustrate five frames and the three slices of water and fat images, respectively. These images show the separated fat component along with flow artifacts. The bottom row shows the generated temporal and spatial profiles of the proton density fat fraction map. This information related to water-fat imaging demonstrates the flexibility of the exemplary TMRF to include a desired contrast.
[0061] Figures 8A and 8B show exemplary quantitative images according to an exemplary embodiment of the present disclosure. In particular, Figure 8A shows top map T1 and bottom map T2 over three slices for a representative data set computed using MRF. Figure 8B shows the corresponding maps for the exemplary TMRF. It can be observed that both methods can yield similar maps for T1 and T2, which conforms to previously published literature values. Figure 8C shows an exemplary graph illustrating the comparison of the mean to standard deviation in each of the contrast blocks according to an exemplary embodiment of the present disclosure. As shown in Figures 8A-8C, the exemplary TMRF provides superior mean to standard deviation for PD and T2. This can also observed in Figures 2A-2D and Figures 6A-6C. T1 contrast from the exemplary TMRF has increased contrast but can have lower mean values, which can be attributed to the decreased signal strength seen from TMRF T1 block.
Exemplary Utility of TMRF Generated Data for Diagnostic Viability
[0062] TMRF generated data was validated and optimized in following two steps during the first six months.
[0063] Exemplary In Vitro Studies - A two-site repeatability study was performed using the National Institute of Standards and Technology (“NIST”) quantitative MRI phantom.
The repeatability of T1 and T2 estimation of the exemplary TMRF was compared to MRF and the Gold Standard (“GS”) measurements at two sites. The GS included Spin Echo based measurements of T1 and T2. All data was acquired on ISMRM/NIST (QalibreMD Inc., CO, USA) which has three layers of spheres with a large range of T1, T2 and Proton Density values. For in vitro studies, data for 10 days was acquired on phantom using three methods (e.g., GS, MRF and TMRF) at two sites with the same vendor and field strength (e.g., on GE 3T Discovery 750w). GS T1 and T2 measurements were conducted using inversion recovery spin echo (“IRSE”) and spin echo (“SE”) sequences. The data were reconstructed and matched with an EPG simulated dictionary. ROI analysis was performed to get T1 and T2 estimation for 14 spheres. From the exemplary results, the underestimation of TMRF of long T1 values can be corrected by a constant bias term. For T2 MRF values, site 1 can be close to literature reported values by NIST while site 2 can be close to GS values. The higher order of standard deviation in T2 values can be due to Bi variations. Figure 9 shows the comparison of T1 map and T2 map obtained from MRF and TMRF across two sites. MRF and TMRF data were averaged from 10-day repeatability scans. The data were obtained from the distinct 14 spheres present in the T1 and T2 plates of ISMRM/NIST phantom (QalibreMD Inc., CO, USA). A manual ROI analysis was performed on eight datasets. 100 pixels were selected for each sphere in T1 maps while 70 pixels were selected for T2 maps due to differences in field of view; MRF and TMRF.
[0064] Further Exemplary In Vivo Studies : Multi-contrast studies on healthy human brain. Simultaneous non-synthetic multi-contrast and quantitative images were rapidly acquired by tailoring the MR fingerprinting acquisition schedule in contrast blocks. TMRF for five contrasts was designed, simulated and demonstrated on five healthy volunteers.
Figure 10 shows a set of exemplary reconstructed images according to an exemplary embodiment of the present disclosure for the GS, MRF synthetic contrast, TMRF synthetic contrast, and TMRF natural contrasts. Each row includes four different contrasts - T1 weighted, T1 FLAIR, T2 weighted and Short T1 inversion Recovery (“STIR”). Arrows 1005 and ellipses 1010 in the second and third columns indicate the flow artifacts ( e.g ., hyper- intensities or holes) which are not present in the GS and TMRF natural contrast images. The gold standard acquired was magnetization prepared T1 weighted instead of T1 weighted. All images were acquired on a 3T GE 750w scanner.
[0065] MRF and TMRF studies were performed on 4 normal volunteers. Both methods involved an axial, 20 slice brain coverage with a slice thickness of 5mm, and lmm in plane resolution. GS sequences - T1 weighted, T1 FLAIR, T2 weighted and STIR were acquired on the human brain. The spatio-temporal profiles of T1, T2, PD, water-fat and flow contrasts were reconstructed block-wise along with relaxometric maps. These MRF and TMRF maps were used to generate the synthetic contrast images using the signal equation and were compared. The TMRF natural contrast images were compared with the GS. Figure 10 shows the representative images obtained from GS (e.g., natural contrast), MRF (e.g, synthetic contrast) and TMRF (e.g, synthetic and natural contrasts) for the following sequences - magnetization prepared T1 weighted, T1 FLAIR, T2 weighted and STIR. The acquisition times for MRF and TMRF were 5:57 and 5:27 (min:see) respectively.
Exemplary Image Quality
[0066] Exemplary EOF Analysis: Three ROIs each corresponding to WM, GM and CSF were manually drawn on the central slices (#9). The ROI masks were multiplied with the thousand-point time series of contrast images to obtain the mean ROI signal evolution of the three types of tissue matter for each of the volunteers for both schedules. These ROI masks were also multiplied with the parametric maps to obtain mean T1, T2 and PD values for each volunteer. A mean of means was computed for T1, T2 values for MRF and the TMRF schedules for the three types of tissue matter. The TMRF T1 values for GM and WM can be higher compared to corresponding MRF values. However, these values can be within the range of T1 values as reported in literature. The mean ROI values can be closely grouped across volunteers and corresponding values for the two schedules can be similar to one another.
[0067] Exemplary Image Analysis: Figures 1 lA-1 II illustrated exemplary images comparing image quality according to an exemplary embodiment of the present disclosure. The synthetically generated images (e.g., MRF and TMRF) can be compared with the natural contrast images obtained from GS (see. e.g., Figures 1 lA-1 ID) and the TMRF images produced after sliding window reconstruction as shown in Figures 1 IB-1 ID. The magnified version of the edges (e.g., between white matter and gray matter) and the presence/absence of flow artifact are shown in Figures 11E-1I. It can be observed that natural contrast from GS and TMRF can be resistant to flow artifact and also the boundaries can be smooth as compared to synthetic contrast. Additionally, the exemplary GS images can be similar to that of natural contrast obtained from TMRF. The data shown here can be representative T1 FLAIR contrast and the results can be similar to other contrast.
[0068] The exemplary DIXON images obtained from the natural TMRF data do not suppress the water/fat accurately. The exemplary system, method and computer-accessible medium can optimize the DIXON method for better suppression of the fat and water. These optimized images can be compared with the LAVA Flex sequence.
[0069] The comparison of the image quality between the GS is shown in Figure 11 A, between exemplary synthetic images generated from MRF is shown in Figure G1B, between synthetic images from TMRF is shown in Figure 1 1C, and between natural images obtained after sliding window reconstruction from TMRF is shown in Figure 1 1C. Ail illustrated images are T1 FLAIR contrast. A part of the image is magnified (e.g., as shown in the boxes) to see the effect of flow as shown in Figures 11E-1 1G and boundaries between gray matter and white matter, as shown in Figures 11H and 11I. The synthetic images shown in Figures 1 1H and HI show' pixel ated and patchy behavior due to the fitting of a single tissue type per voxel as compared to those shown in Figures 11E- 11G, which are smoother due to the absence of regression fits. Natural contrast obtained from GS and TMRF may not have any flow-' artifacts and can be smoother as compared to synthetic contrast data.
Exemplary Benefits of TMRF
1) Simultaneous non-synthetic multi-contrast and quantitative imaging.
2) Simple, quick and robust MRF protocol implementation that is flexible to include new contrasts such as water-fat imaging.
3) Challenges involving sliding window reconstruction along the temporal dimension due to constancy can be overcome.
4) Delivery of high insensitivity to flow artifacts caused by acquisition.
5) Scalable fingerprinting framework for tissue parameters measurable directly and indirectly with MRI - conductivity, temperature, etc.
6) Increased degrees of freedom in acquisition - randomization of trajectories, combination of sequence parameters for diverse contrasts such as, but not limited to, perfusion (e.g., contrast and non-contrast methods), diffusion, blood flow, etc.
7) Significant reduction in reconstruction computation times - as compared to analytical (including multiple variants of Fourier transform) methods relying on gridding or iterative reconstruction for non-Cartesian and/or under-sampled acquisitions.
8) Use of histopathological data as reference for fingerprints of pathology - training MRF sequences on stack of histopathological slides to understand the MRF signatures of such data. 9) Potential avoidance of biopsies in such anatomies.
10) Rapid comprehensive MRI exams, for example, pediatric neuroimaging, multi- parametric prostate imaging, whole body imaging oncology, diabetes studies inclusive of NASH, study of fat types such as brown, white and brite fat, etc. 11) MR value driven protocols, for example, 5 -minute stroke protocol, as an alternate to EPImix, MAGIC, etc.
12) Multi-scale, multi-modality image fusion e.g., rapid MR-PET exams for
oncological applications, whole body metabolic disorders and neuro-psychiatric diseases such as AD, PD, MS and SZ.
13) Atlas creation at higher field strengths to deliver increased information content at lower fields - synthesis of tissue parametric maps at higher fields could be utilized to train FCNs and employed for data generated from lower field strengths with appropriate correction factors that are field dependent.
[0070] Figure 12 shows a flow diagram illustrating an exemplary method for generating a quantitative image of a patient or a non-synthetic contrast image of the patient according to an exemplary embodiment of the present disclosure. For example, at procedure 1205, a degree of flip angles for the first MR signal can be selected based on the particular contrast. At procedure 1210, a first magnetic resonance (MR) signal can be generated, which can be directed to the patient. At procedure 1215, a second MR signal can be received from the patient that can be based on the first MR signal. At procedure 1220, the second MR signal can be pre-processed by compensating for a calibrated gradient delay. At procedure 1225, k- space of the second MR signal can be scaled with a ratio of a field of view to a matrix size. At procedure 1230, spikes in the second MR signal can be removed. At procedure 1235, k- space data in the second MR signal can be weighted with a predetermined density compensation factor. At procedure 1240, the particular image can be generated based on the second MR signal, for example, using a reconstruction procedure.
[0071] Figure 13 shows a block diagram of an exemplary embodiment of a system according to the present disclosure. For example, exemplary procedures in accordance with the present disclosure described herein can be performed by a processing arrangement and/or a computing arrangement (e.g., computer hardware arrangement) 1305. Such
processing/computing arrangement 1305 can be, for example entirely or a part of, or include, but not limited to, a computer/processor 1310 that can include, for example one or more microprocessors, and use instructions stored on a computer-accessible medium (e.g., RAM, ROM, hard drive, or other storage device).
[0072] As shown in Figure 13, for example a computer-accessible medium 1315 (e.g., as described herein above, a storage device such as a hard disk, floppy disk, memory stick, CD- ROM, RAM, ROM, etc., or a collection thereol) can be provided (e.g., in communication with the processing arrangement 1305). The computer-accessible medium 1315 can contain executable instructions 1320 thereon. In addition or alternatively, a storage arrangement 1325 can be provided separately from the computer-accessible medium 1315, which can provide the instructions to the processing arrangement 1305 so as to configure the processing arrangement to execute certain exemplary procedures, processes, and methods, as described herein above, for example.
[0073] Further, the exemplary processing arrangement 1305 can be provided with or include an input/output ports 1335, which can include, for example a wired network, a wireless network, the internet, an intranet, a data collection probe, a sensor, etc. As shown in Figure 13, the exemplary processing arrangement 1305 can be in communication with an exemplary display arrangement 1330, which, according to certain exemplary embodiments of the present disclosure, can be a touch-screen configured for inputting information to the processing arrangement in addition to outputting information from the processing arrangement, for example. Further, the exemplary display arrangement 1330 and/or a storage arrangement 1325 can be used to display and/or store data in a user-accessible format and/or user-readable format.
[0074] The foregoing merely illustrates the principles of the disclosure. Various modifications and alterations to the described embodiments will be apparent to those skilled in the art in view of the teachings herein. It will thus be appreciated that those skilled in the art will be able to devise numerous systems, arrangements, and procedures which, although not explicitly shown or described herein, embody the principles of the disclosure and can be thus within the spirit and scope of the disclosure. Various different exemplary embodiments can be used together with one another, as well as interchangeably therewith, as should be understood by those having ordinary skill in the art. In addition, certain terms used in the present disclosure, including the specification, drawings and claims thereof, can be used synonymously in certain instances, including, but not limited to, for example, data and information. It should be understood that, while these words, and/or other words that can be synonymous to one another, can be used synonymously herein, that there can be instances when such words can be intended to not be used synonymously. Further, to the extent that the prior art knowledge has not been explicitly incorporated by reference herein above, it is explicitly incorporated herein in its entirety. All publications referenced are incorporated herein by reference in their entireties. EXEMPLARY REFERENCES
[0075] The following references are hereby incorporated by reference in their entireties:
1. Yankeelov TE, Abramson RG, Quarles CC. Quantitative multimodality imaging in cancer research and therapy. Nat. Rev. Clin. Oncol. 2014;11 :670-680.
2. Abramson R, Burton K, John-Paul J, radiology ES-A, 2015 undefined. Methods and challenges in quantitative imaging biomarker development. Elsevier.
3. Doran SJ. Chapter 21 Quantitative magnetic resonance imaging: applications and
estimation of errors. Data Handl Sci. Technol. 1996;18:452-488 doi: 10.1016/S0922- 3487(96)80058-X.
4. Jara H. Theory of quantitative magnetic resonance imaging.
5. Seow P, Wong JHD, Ahmad-Annuar A, Mahajan A, Abdullah NA, Ramli N. Quantitative magnetic resonance imaging and radiogenomic biomarkers for glioma characterisation: a systematic review. Br. J. Radiol. 2018;91 :20170930 doi: 10.1259/bjr.20170930.
6. Deoni SCL, Rutt BK, Peters TM. Rapid combinedT1 andT2 mapping using gradient recalled acquisition in the steady state. Magn. Reson. Med. 2003;49:515-526 doi:
10.1002/mrm.10407.
7. Frost R, Jezzard P, Douaud G, Clare S, Porter DA, Miller KL. Scan time reduction for readout-segmented EPI using simultaneous multislice acceleration: Diffusion-weighted imaging at 3 and 7 Tesla. Magn. Reson. Med. 2015;74: 136-149 doi:
10.1002/mrm.25391.
8. Shukla-Dave A, Obuchowski NA, Chenevert TL, et al. Quantitative imaging biomarkers alliance (QIBA) recommendations for improved precision of DWI and DCE-MRI derived biomarkers in multicenter oncology trials. J. Magn. Reson. Imaging 2018 doi:
10.1002/jmri.26518. 9. Tanenbaum LN, Tsiouris AJ, Johnson AN, et al. Synthetic MRI for clinical
neuroimaging: results of the Magnetic Resonance Image Compilation (MAGiC) prospective, multicenter, multireader trial. Am. J. Neuroradiol. 2017.
10. Gonçalves F, Serai S, ... GZMR, 2018 undefined. Synthetic Brain MRI: Review of
Current Concepts and Future Directions joumals.lww.com.
11. Ma D, Gulani V, Seiberlich N, et al. Magnetic resonance fingerprinting. Nature 2013 doi:
10.1038/nature 11971.
12. Jiang Y, Ma D, Seiberlich N, Gulani V, Griswold MA. MR fingerprinting using fast imaging with steady state precession (FISP) with spiral readout. Magn. Reson. Med. 2015;74: 1621- 1631 doi: 10.1002/mrm.25559.
13. Brown R, Haacke E, Cheng Y, Thompson M. Magnetic resonance imaging: physical principles and sequence design.; 2014.
14. Virtue P, Tamir JI, Doneva M, Yu SX, Lustig M. Learning Contrast Synthesis from MR Fingerprinting.
15. Hanicke W, Vogel HU. An analytical solution for the SSFP signal in MRI. Magn.
Reson. Med. 2003;49:771-775 doi: 10.1002/mrm. l0410.
16. Bernstein M a., King KF, Zhou XJ. Handbook of MRI Pulse Sequences, Appendix II.;
2004. doi: 10.1016/B978-012092861-3/50026-1.
17. Weigel M. Extended phase graphs: Dephasing, RF pulses, and echoes - pure and simple.
J. Magn. Reson. Imaging 2015;41 :266-295 doi: 10.1002/jmri.24619.
18. Cao X, Liao C, Wang Z, et al. Robust sliding- window reconstruction for Accelerating the acquisition of MR fingerprinting. Magn. Reson. Med. 2017;78: 1579-1588 doi:
10.1002/mrm.26521.
19. Fessler J. Image reconstruction toolbox. Available at website http://www.eecs.umich.edu
/fessler/code 2012 doi: 10.1007/978-0-387-85922-4. 20. Dixon WT. Simple proton spectroscopic imaging. Radiology 1984;153: 189-194 doi:
10.1148/radiology.153.1.6089263.
21. Wansapura JP, Holland SK, Dunn RS, Ball WS. NMR relaxation times in the human brain at 3.0 tesla. J. Magn. Reson. Imaging 1999;9:531-538 doi: 10.1002/(SICI)1522- 2586(199904)9:4<531 : : AID-JMRI4>3.0.CO;2-L.
22. Wansapura JP, Holland SK, Dunn RS, Ball WS. NMR relaxation times in the human brain at 3.0 Tesla. J. Magn. Reson. Imaging 1999;9:531-538 doi: 10.1002/(SICI)1522- 2586(199904)9:4<531 : : AID-JMRI4>3.0.CO;2-L.
23. Smith SM, Jenkinson M, Woolrich MW, et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 2004;23:S208-S219 doi:
10.1016/J.NEUROIMAGE.2004.07.051.
24. Reeder SB, McKenzie CA, Pineda AR, et al. Water-fat separation with IDEAL gradient- echo imaging. J. Magn. Reson. Imaging 2007;25:644-652 doi: 10.1002/jmri.20831.).
25. Ravi, K. S., Geethanath, S., Jochen, W., & Vaughan, J. T. (2018a). MR value driven autonomous mri using imr-framework. In ISMRM workshop on machine learning part ii, Washington, D.C.
26. Ravi, K. S., Geethanath, S., & Vaughan, J. T. (2018b). Imr-framework for rapid design and deployment of non-cartesian sequences. In I2i workshop. Retrieved from
cai2r.net/i2i.
27. Ravi, K. S., Geethanath, S., & Vaughan, J. T. (2019a). Autonomous scanning using imr framework to improve mr accessibility. In ISMRM workshop on accessible mri for the world. New Delhi, India.
28. Ravi, K. S., Geethanath, S., & Vaughan, J. T. (2019b). Self-administered exam using autonomous magnetic resonance imaging (amri). In ISMRM 27th annual meeting and exhibition. Montreal, Canada. 29. Ravi, K. S., Potdar, S., Poojar, P., Reddy, A. K., Kroboth, S., Nielsen, J.-F., Zaitsev, M., et al. (2018). Pulseq graphical programming interface: Open source visual environment for prototyping pulse sequences and integrated magnetic resonance imaging algorithm development. Magnetic resonance imaging, 52, 9-15. doi: 10.1016/j.mri.2018.03.008. 30. Keerthi Sravan Ravi , Sairam Geethanath , Girish Srinivasan , Rahul Sharma , Sachin R
Jambawalikar , Angela Lignelli-Dipple , and John Thomas Vaughan Jr; Deep learning Assisted Radiological reporT (DART). In ISMRM 28th Annual Meeting & Exhibition in Sydney, Australia .
31. Keerthi Sravan Ravi, Sairam Geethanath, Patrick Quarterman, Maggie Fung, John
Thomas Vaughan Jr.,“Intelligent Protocolling for Autonomous MRI". ISMRM Annual Meeting & Exhibition. ISMRM, 2020.

Claims

WHAT IS CLAIMED IS:
1. A non-transitory computer-accessible medium having stored thereon computer-executable instructions for generating at least one particular image which is at least one of (i) at least one quantitative image of at least one section of at least one patient or (ii) at least one non- synthetic contrast image of the at least one section of the at least one patient, wherein, when a computer arrangement executes the instructions, the computer arrangement is configured to perform procedures comprising:
generating a first magnetic resonance (MR) signal and directing the first MR signal to the at least one patient;
receiving a second MR signal from the at least one patient that is based on the first MR signal; and
generating the at least one particular image based on the second MR signal.
2. The computer-accessible medium of claim 1, wherein the first MR signal is a configured MR signal.
3. The computer-accessible medium of claim 2, wherein the tailored MR signal is configured for a particular contrast.
4. The computer-accessible medium of claim 3, wherein the first MR signal has a constant signal intensity.
5. The computer-accessible medium of claim 4, wherein the computer arrangement is configured to generate the first MR signal based on a degree of a plurality of flip angles that maintains the constant signal intensity.
6. The computer-accessible medium of claim 3, wherein the computer arrangement is configured to select a degree of flip angles for the first MR signal based on the particular contrast.
7. The computer-accessible medium of claim 6, wherein the degree of the flip angles varies within a particular range.
8. The computer-accessible medium of claim 7, wherein the degree of the flip angles varies about a mean value.
9. The computer-accessible medium of claim 8, wherein the degree of the flip angles varies monotonously about the mean value.
10. The computer-accessible medium of claim 7, wherein the degree of the flip angles varies pseudo randomly within the particular range.
11. The computer-accessible medium of claim 7, wherein the particular range is about - 5 +/- 4 degrees, about 45 +/- 5 degrees, about 75+/- 5 degrees, or about 75 +/- 5 degrees.
12. The computer-accessible medium of claim 3, wherein the particular contrast includes at least one of T1, T2, proton density, water, fat, off resonance, diffusion, perfusion, or flow.
13. The computer-accessible medium of claim 1, wherein the at least one non-synthetic contrast image is at least one non-synthetic multi-contrast image.
14. The computer-accessible medium of claim 1, wherein the computer arrangement is configured to generate the at least one particular image using a reconstruction procedure.
15. The computer-accessible medium of claim 14, wherein the reconstruction procedure is a sliding window reconstruction procedure.
16. The computer-accessible medium of claim 15, wherein the reconstruction procedure includes converting the second MR signal to an image using a Non-Uniform Fast Fourier Transform.
17. The computer-accessible medium of claim 1, wherein the computer arrangement is further configured to generate MR information based on the second MR signal by:
pre-processing the second MR signal by compensating for a calibrated gradient delay; scaling k-space of the second MR signal with a ratio of a field of view to a matrix size;
removing spikes in the second MR signal; and
weighting k-space data in the second MR signal with a predetermined density compensation factor, and
wherein the at least one particular image is generated based on the MR information.
18. The computer-accessible medium of claim 1, wherein the computer arrangement is configured to generate the at least one quantitative image by vector-dot product matching L2- norm normalized dictionary entries with voxel signal evolutions in the second MR signal.
19. A method for generating at least one particular image which is at least one of (i) at least one quantitative image of at least one section of at least one patient or (ii) at least one non- synthetic contrast image of the at least one section of the at least one patient, comprising: generating a first magnetic resonance (MR) signal and directing the first MR signal to the at least one patient;
receiving a second MR signal from the at least one patient that is based on the first MR signal; and
generating the at least one particular image based on the second MR signal.
20. The method of claim 19, wherein the first MR signal is a configured MR signal.
21. The method of claim 20, wherein the tailored MR signal is configured for a particular contrast.
22. The method of claim 21, wherein the first MR signal has a constant signal intensity.
23. The method of claim 22, further comprising generating the first MR signal based on a degree of a plurality of flip angles that maintains the constant signal intensity.
24. The method of claim 21, further comprising selecting a degree of flip angles for the first MR signal based on the particular contrast.
25. The method of claim 24, wherein the degree of the flip angles varies within a particular range.
26. The method of claim 25, wherein the degree of the flip angles varies about a mean value.
27. The method of claim 26, wherein the degree of the flip angles varies monotonously about the mean value.
28. The method of claim 25, wherein the degree of the flip angles varies pseudo randomly within the particular range.
29. The method of claim 25, wherein the particular range is about - 5 +/- 4 degrees, about 45 +/- 5 degrees, about 75+/- 5 degrees, or about 75 +/- 5 degrees.
30. The method of claim 21, wherein the particular contrast includes at least one of T1, T2, proton density, water, fat, off resonance, diffusion, perfusion, or flow.
31. The method of claim 19, wherein the at least one non-synthetic contrast image is at least one non-synthetic multi-contrast image.
32. The method of claim 19, further comprising generating the at least one particular image using a reconstruction procedure.
33. The method of claim 32, wherein the reconstruction procedure is a sliding window reconstruction procedure.
34. The method of claim 33, wherein the reconstruction procedure includes converting the second MR signal to an image using a Non-Uniform Fast Fourier Transform.
35. The method of claim 19, further comprising generating MR information based on the second MR signal by:
pre-processing the second MR signal by compensating for a calibrated gradient delay; scaling k-space of the second MR signal with a ratio of a field of view to a matrix size;
removing spikes in the second MR signal; and
weighting k-space data in the second MR signal with a predetermined density compensation factor, and
wherein the at least one particular image is generated based on the MR information.
36. The method of claim 19, further comprising generating the at least one quantitative image by vector-dot product matching L2-norm normalized dictionary entries with voxel signal evolutions in the second MR signal.
37. A system for generating at least one particular image which is at least one of (i) at least one quantitative image of at least one section of at least one patient or (ii) at least one non- synthetic contrast image of the at least one section of the at least one patient, comprising: a computer hardware arrangement configured to:
generate a first magnetic resonance (MR) signal and directing the first MR signal to the at least one patient;
receive a second MR signal from the at least one patient that is based on the first MR signal; and
generate the at least one particular image based on the second MR signal.
38. The system of claim 37, wherein the first MR signal is a configured MR signal.
39. The system of claim 38, wherein the tailored MR signal is configured for a particular contrast.
40. The system of claim 39, wherein the first MR signal has a constant signal intensity.
41. The system of claim 40, wherein the computer hardware arrangement is configured to generate the first MR signal based on a degree of a plurality of flip angles that maintains the constant signal intensity.
42. The system of claim 39, wherein the computer hardware arrangement is configured to select a degree of flip angles for the first MR signal based on the particular contrast.
43. The system of claim 42, wherein the degree of the flip angles varies within a particular range.
44. The system of claim 43, wherein the degree of the flip angles varies about a mean value.
45. The system of claim 44, wherein the degree of the flip angles varies monotonously about the mean value.
46. The system of claim 43, wherein the degree of the flip angles varies pseudo randomly within the particular range.
47. The system of claim 43, wherein the particular range is about - 5 +/- 4 degrees, about 45 +/- 5 degrees, about 75+/- 5 degrees, or about 75 +/- 5 degrees.
48. The system of claim 39, wherein the particular contrast includes at least one of T1, T2, proton density, water, fat, off resonance, diffusion, perfusion, or flow.
49. The system of claim 37, wherein the at least one non-synthetic contrast image is at least one non-synthetic multi-contrast image.
50. The system of claim 37, wherein the computer hardware arrangement is configured to generate the at least one particular image using a reconstruction procedure.
51. The system of claim 50, wherein the reconstruction procedure is a sliding window reconstruction procedure.
52. The system of claim 51, wherein the reconstruction procedure includes converting the second MR signal to an image using a Non-Uniform Fast Fourier Transform.
53. The system of claim 37, wherein the computer hardware arrangement is further configured to generate MR information based on the second MR signal by:
pre-processing the second MR signal by compensating for a calibrated gradient delay; scaling k-space of the second MR signal with a ratio of a field of view to a matrix size;
removing spikes in the second MR signal; and weighting k-space data in the second MR signal with a predetermined density compensation factor, and
wherein the at least one particular image is generated based on the MR information.
54. The system of claim 37, wherein the computer hardware arrangement is configured to generate the at least one quantitative image by vector-dot product matching L2-norm normalized dictionary entries with voxel signal evolutions in the second MR signal.
PCT/US2020/022996 2019-03-15 2020-03-16 Accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging WO2020190879A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CA3133772A CA3133772A1 (en) 2019-03-15 2020-03-16 Accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging
EP20773478.1A EP3937775A4 (en) 2019-03-15 2020-03-16 Accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging
US17/475,687 US11864863B2 (en) 2019-03-15 2021-09-15 System, method and computer-accessible medium for accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201962819159P 2019-03-15 2019-03-15
US62/819,159 2019-03-15

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/475,687 Continuation US11864863B2 (en) 2019-03-15 2021-09-15 System, method and computer-accessible medium for accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging

Publications (1)

Publication Number Publication Date
WO2020190879A1 true WO2020190879A1 (en) 2020-09-24

Family

ID=72521260

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2020/022996 WO2020190879A1 (en) 2019-03-15 2020-03-16 Accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging

Country Status (4)

Country Link
US (1) US11864863B2 (en)
EP (1) EP3937775A4 (en)
CA (1) CA3133772A1 (en)
WO (1) WO2020190879A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT202000023257A1 (en) * 2020-10-02 2022-04-02 Esaote Spa SYSTEM AND METHOD FOR DIAGNOSTIC IMAGING OF WATER AND FAT

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140085318A1 (en) * 2012-09-26 2014-03-27 Siemens Corporation Multi-GPU FISTA Implementation for MR Reconstruction with Non-Uniform K-Space Sampling
US20170007148A1 (en) * 2015-07-07 2017-01-12 Tesla Health, Inc Fast scaning based on magnetic resonance history
US20170082715A1 (en) * 2015-09-18 2017-03-23 Samsung Electronics Co., Ltd. Magnetic resonance imaging apparatus and method of scanning blood vessel
WO2018108643A1 (en) * 2016-12-15 2018-06-21 Koninklijke Philips N.V. Multi-state magnetic resonance fingerprinting

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10761171B2 (en) * 2015-06-22 2020-09-01 Case Western Reserve University Systems and methods for free-breathing three-dimensional magnetic resonance fingerprinting

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140085318A1 (en) * 2012-09-26 2014-03-27 Siemens Corporation Multi-GPU FISTA Implementation for MR Reconstruction with Non-Uniform K-Space Sampling
US20170007148A1 (en) * 2015-07-07 2017-01-12 Tesla Health, Inc Fast scaning based on magnetic resonance history
US20170082715A1 (en) * 2015-09-18 2017-03-23 Samsung Electronics Co., Ltd. Magnetic resonance imaging apparatus and method of scanning blood vessel
WO2018108643A1 (en) * 2016-12-15 2018-06-21 Koninklijke Philips N.V. Multi-state magnetic resonance fingerprinting

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3937775A4 *

Also Published As

Publication number Publication date
EP3937775A4 (en) 2022-12-14
US20220071490A1 (en) 2022-03-10
EP3937775A1 (en) 2022-01-19
CA3133772A1 (en) 2020-09-24
US11864863B2 (en) 2024-01-09

Similar Documents

Publication Publication Date Title
US10345408B2 (en) Low field magnetic resonance imaging (MRI) scanner for cardiac imaging
CN109477877B (en) Magnetic resonance imaging system and method
US10561337B2 (en) Rapid 3D dynamic arterial spin labeling with a sparse model-based image reconstruction
CN102920455B (en) For operating the method for Fat and moisture amount in nuclear magnetic resonance independently
US11796618B2 (en) Systems and methods for magnetic resonance imaging
US10845443B2 (en) System and method for magnetic resonance imaging reconstruction using novel k-space sampling sequences
CN1973211A (en) Contrast prepared MRI involving non-cartesian trajectories with oversampling of the center of K-space
US10267883B2 (en) System and method for motion resolved MRI
Dietrich et al. 3D Free‐breathing multichannel absolute Mapping in the human body at 7T
Chen et al. Self‐calibrating wave‐encoded variable‐density single‐shot fast spin echo imaging
US8704518B2 (en) System and method of high signal-to-noise ratio magnetic resonance imaging screening
US20220381855A1 (en) Systems and methods for pulse sequence configuration in magnetic resonance imaging
CN111012349B (en) System and method for signal characterization determination in magnetic resonance imaging
US11864863B2 (en) System, method and computer-accessible medium for accelerated, simultaneous quantitative and non-synthetic multi-contrast imaging
EP4071494A1 (en) Method for acquiring a three-dimensional magnetic resonance image dataset and for generating a motion-corrected image dataset
US11287504B2 (en) Systems and methods for magnetic resonance imaging
US9709651B2 (en) Compensated magnetic resonance imaging system and method for improved magnetic resonance imaging and diffusion imaging
US11069063B2 (en) Systems and methods for noise analysis
Ma et al. Mitigating transmit‐B1 artifacts by predicting parallel transmission images with deep learning: A feasibility study using high‐resolution whole‐brain diffusion at 7 Tesla
US11935211B2 (en) Systems and methods for image processing
US11249159B2 (en) Systems and methods for enhancement of resolution for strategically acquired gradient echo (stage) imaging
Poojar et al. Tailored Magnetic Resonance Fingerprinting
Marchetto Improved motion-correction for MRI with markerless face-tracking
Ge et al. Three-dimensional high resolution venography using susceptibility weighted imaging at 7T
Qiu et al. Direct synthesis of multi‐contrast brain MR images from MR multitasking spatial factors using deep learning

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

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3133772

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2020773478

Country of ref document: EP