WO2016036643A1 - Selective sampling magnetic resonance-based method for assessing structural spatial frequencies - Google Patents
Selective sampling magnetic resonance-based method for assessing structural spatial frequencies Download PDFInfo
- Publication number
- WO2016036643A1 WO2016036643A1 PCT/US2015/047676 US2015047676W WO2016036643A1 WO 2016036643 A1 WO2016036643 A1 WO 2016036643A1 US 2015047676 W US2015047676 W US 2015047676W WO 2016036643 A1 WO2016036643 A1 WO 2016036643A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- value
- gradient
- pulse
- signal
- interest
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
Definitions
- the herein claimed method relates to the field of diagnostic assessment of fine textures in biological systems for pathology assessment and disease diagnosis, and in material and structural evaluation in industry and in engineering research. More specifically, the invention employs a method for repeat measurement of k-values associated with the spatial organization of biologic tissue texture, with the MRI machine gradients turned off. This allows assessment of tissue texture on a time-scale on the order of a msec, whereby the problem of patient motion becomes negligible. The method enables in vivo assessment, towards diagnosis and monitoring, of disease and therapy- induced textural changes in tissue.
- Representative targets of the technique are: 1) for assessment of changes to trabecular architecture caused by bone disease, allowing assessment of bone health and fracture risk, 2) evaluation of fibrotic development in soft tissue diseases such as, for example, liver, lung, and heart disease, and 3) changes to fine structures in neurologic diseases, such as the various forms of dementia, or in cases of brain injury and downstream neuro -pathology as in, for example, Traumatic Brain Injury (TBI) and Chronic Traumatic Encephalopathy (CTE), or for characterization and monitoring of abnormal neurologic conditions such as autism and schizophrenia.
- Other pathology applications include assessment of vascular changes such as in in the vessel network surrounding tumors or associated with development of CVD (Cerebrovascular Disease), and of changes in mammary ducting in response to tumor growth.
- the invention also has applications in assessment of fine structures for a range of industrial purposes such as measurement of material properties in manufacturing or in geology to characterize various types of rock, as well as other uses for which measurement of fine structures/textures is needed.
- CT Tomography
- PET Positron Emission Tomography
- Bone health is compromised by ageing, by bone cancer, as a side effect of cancer treatments, diabetes, rheumatoid arthritis, and as a result of inadequate nutrition, among other causes. Bone disease affects over ten million people annually in the US alone, adversely affecting their quality of life and reducing life expectancy.
- the current diagnostic standard is Bone Mineral Density (BMD), as measured by the Dual Energy X-ray Absorptiometry (DEXA) projection technique.
- This modality yields an areal bone density integrating the attenuation from both cortical and trabecular bone, similar to the imaging mechanism of standard x-ray, but provides only limited information on trabecular architecture within the bone, which is the marker linked most closely to bone strength.
- KANIS, J. AND GLUER, C "An update on the diagnosis and assessment of osteoporosis with densitometry”; Osteoporosis International, Vol. 11, issue 3, 2000.
- LEGRAND E. et al.; "Trabecular bone microarchitecture, bone mineral density, and vertebral fractures in male osteoporosis”; JBMR, Vol. 15, issue 1, 2000.
- BMD correlates only loosely with fracture risk.
- a post-processing technique, TBS Trabecular Bone Score
- TBS Trabecular Bone Score
- a proposed MR-based technique fineSA (JAMES, T., CHASE, D.; "Magnetic field gradient structure characteristic assessment using one dimensional (ID) spatial-frequency distribution analysis” ; US patent No. 7,932,720 B2; April 26, 201 1.), attempts to circumvent the problem of patient motion by acquiring a much smaller data matrix of successive, finely-sampled, one-dimensional, frequency-encoded acquisitions which are subsequently combined to reduce noise. Imaging in this case is reduced to one dimension, reducing the size of the data matrix acquired and, hence, the acquisition time.
- SNR Signal to Noise
- Fibrotic pathology occurs in a large number of diseases, from lung and liver fibrosis, to cardiac and cystic fibrosis, pancreatic fibrosis, muscular dystrophy, bladder and heart diseases, and myelofibrosis, in which fibrotic structures replace bone marrow. Fibrotic development is attendant in several cancers, such as breast cancer. A different pathology development is seen in prostate cancer, where the disease destroys healthy organized fibrous tissue. In all cases, textural spacings highlighted in the tissue change in response to disease progression, as collagen fibers form along underlying tissue structures. In liver disease, the textural wavelength changes as the healthy tissue texture in the liver is replaced by a longer wavelength texture originating from the collagen "decoration" of the lobular structure in the organ. In other organs/diseases, textural change reflects the upset in healthy tissue with development of texture indicative of fibrotic intervention.
- tissue biopsy a highly invasive and often painful procedure with a non-negligible morbidity— and mortality— risk (patients need to stay at the hospital for post-biopsy observation for hours to overnight), and one that is prone to sampling errors and large reading variation.
- Magnetic Resonance-based Elastography which has been under development for some time for use in assessment of liver disease, is not capable of early-stage assessment— the read errors are too large prior to significant fibrotic invasion (advanced disease). Further, this technique requires expensive additional hardware, the presence of a skilled technician, and takes as much as 20 minutes total set up and scanning time, making it a very costly procedure.
- the ability to image fibrotic texture directly by MR imaging is compromised both by patient motion over the time necessary to acquire data and by lack of contrast between the fibers and the surrounding tissue. Even acquisition during a single breath hold is severely compromised by cardiac pulsatile motion and noncompliance to breath hold, which results in significant motion at many organs, such as liver and lungs.
- AD Alzheimer's Disease
- ALS Amyotrophic Lateral Sclerosis
- Parkinson's disease conditions precipitated by Traumatic Brain Injury (TBI) such as Chronic Traumatic Encephalopathy (CTE)
- TBI Traumatic Brain Injury
- CTE Chronic Traumatic Encephalopathy
- MS Multiple Sclerosis
- CVD Cerebrovascular Disease
- other neurologic diseases are often only diagnosable in advanced stages by behavioral and memory changes, precluding the ability for early stage intervention.
- conditions such as epilepsy and autism have been associated with abnormal variations in fine neuronal structures, which, if clinically diagnosable, would allow targeted selection for testing therapy response.
- PARA 14 Use of CSF biomarkers for dementia diagnosis is painful and highly invasive and cannot differentiate signal levels by anatomic position in the brain, as is possible with imaging biomarkers. As various forms of dementia are found to have different
- CVD Another disease associated with various forms of dementia is CVD
- PARA 16 Another possible neurologic application for the claimed method is to, in vivo, determine the boundaries of the various control regions of the cerebral cortex or the different Brodmann's areas of which these are comprised. Such ability would greatly aid data interpretation in brain function studies, such as those performed using, for example, FMRI (Functional Magnetic Resonance Imaging).
- FMRI Magnetic Magnetic Resonance Imaging
- PARA 17 The three classes of diseases listed above, bone disease, fibrotic diseases, and neurologic diseases are not an all-inclusive list.
- Other disease states in which pathology- induced changes of fine structures occur for instance angiogenic growth of vasculature surrounding a tumor, or fibrotic development and changes in vasculature and mammary gland ducting in response to breast tumor development, also are pathologies wherein the ability to resolve fine tissue textures would enable early detection of disease, and monitoring of response to therapy.
- a method for selective sampling to assess texture using magnetic resonance (MR) is accomplished by exciting a volume of interest. All gradients are then turned off and multiple samples of an RF signal encoded at a specific k- value are recorded in a single TR.
- exciting the volume of interest is accomplished in an exemplary embodiment by transmitting a first RF pulse with a first gradient activated for slice selection, transmitting a second RF pulse with a second gradient activated, said second gradient chosen for slice selective refocusing in a region defined by an intersection of the first slice and a second slice, apply a selected gradient pulse shaped to induce phase wrap to create a spatial encode with resulting pitch corresponding to texture wavelength of the specific k- value and then transmitting a third RF pulse with a third gradient activated.
- the third gradient is adapted to refocus excitation defining a region defined by the intersection of the first and second slices and this third slice selection.
- a selection pulse is issued on a vector combination gradient to alter to a second k- value.
- the gradient is then again turned off and multiple samples of the RF signal are recorded at the second k- value, while in the same TR.
- This embodiment may be extended by issuing additional selection gradient pulses for a predetermined plurality of pulses on the vector combination gradient with each pulse selecting a new k- value. With the vector combination gradient turned off after each pulse, multiple samples of the RF signal at the next k- value induced at each negative pulse are recorded while still in the same TR.
- An additional embodiment provides for transmitting a refocusing RF pulse with the third gradient activated and then issuing a plurality of selection gradient pulses on the vector combination gradient, each pulse selecting a new k- value. After each pulse of the vector combination gradient, multiple samples of the RF signal at the next k- value induced are recorded, again still in the same TR. This embodiment may be further refined by transmitting a second refocusing RF pulse with the third gradient activated and then issuing a plurality of selecting gradient pulses on the vector combination gradient, with each pulse selecting a new k- value. With the vector combination gradient turned off after each pulse, multiple samples of the RF signal at the next k- value induced at each pulse are recorded.
- An additional embodiment provides for recalling the signal multiple times using a GRE (Gradient Recalled Echo), thereby extending the record window to allow acquisition of additional k-value samples.
- GRE Gradient Recalled Echo
- the number of samples at each new k-value following each of the plurality of pulses on the vector combination gradient is selected based upon expected SNR, pathology, tissue contrast, texture size or texture bandwidth.
- the volume of interest in the step of exciting a volume of interest is selected using the set of intersecting slice-selective refocusing, selective excitation using phased-array transmit in combination with appropriate gradients, adiabatic pulse excitation to scramble signal from the tissue outside a region of interest, outer volume suppression sequences, and other methods of selectively exciting spins in an internal volume including physically isolating the tissue of interest.
- the initial embodiment may be employed wherein the steps of exciting a volume of interest, turning off all gradients and recording samples of an RF signal at a specific k-value are repeated multiple times within the same TR or in succeeding TRs to build a magnitude power spectrum, the specific k-values in each TR independent of intervening motion if constrained within the volume of interest for each TR.
- FIG. 1 is a simulation showing the number of data samples required for averaging to achieve an output SNR > 20dB as a function of input SNR;
- FIG. 2 is a simulation showing the number of data samples needed for averaging to achieve a SNR > 20db as a function of location in k-space;
- FIG. 3 is an example timing diagram of a pulse sequence for the claimed method showing the timing of a single TR
- FIG. 4 is a close-up of the example timing diagram of FIG. 3;
- FIG. 5 is an example of a timing diagram for the claimed method, designed to acquire multiple measures of a select set of k-values, with a different number of samples acquired at each k-value to counteract the decrease in energy density at increasing k-value;
- FIG. 6 is a simulation showing that the ability provided by the claimed method to acquire many repeats of a targeted k-value within a single TR enables robust signal averaging to boost SNR;
- FIG. 7 is a simulation showing the results of attempting to acquire 90 samples for averaging using the conventional frequency-encoded echo approach, wherein acquisition of only a small number of repeats of a particular k-value are possible in each TR due to the long record time for each echo;
- FIG. 8 is an example timing diagram for the claimed method designed to provide data acquisition over multiple refocused echoes within a single TR.
- FIGs. 9A and 9B are a depiction of two possible shapes for the acquisition volume of interest (VOI).
- 180° inversion pulse RF pulse that inverts the spins in a tissue region to allow refocusing of the MR signal.
- Adiabatic pulse excitation Adiabatic pulses are a class of amplitude- and frequency- modulated RF -pulses that are relatively insensitive to Bl inhomogeneity and frequency offset effects.
- AWGN Additive White Gaussian Noise Additive white Gaussian noise is a basic noise model used in Information theory to mimic the effect of many random processes that occur in nature.
- Biopsy A biopsy is a sample of tissue extracted from the body in order to examine it more closely.
- CJD Creutzfeld- Jakob Disease Crusher gradients Gradients applied on either side of a 180 RF refocussing slice selection pulse to reduce spurious signals generated by imperfections in the pulse.
- DEXA Dual Energy X-ray Absorptiometry is a means of measuring bone mineral density using two different energy x-ray beams.
- Gradient set the set of coils around the bore of an MR scanner used primarily to spatially encode signal or to set a particular phase wrap in a selected direction
- Isochromat A microscopic group of spins that resonate at the same frequency
- k-value coefficient The coefficient in a Fourier series or transform reflecting the relative weight of each specific k-value in the series.
- Library of k-space values the net collection of k-space coefficients acquired in a particular region of tissue for tissue characterization Machine gradients the magnetic field gradients imposable through use of the set of gradient coils in an MR scanner
- MRE Magnetic Resonance Elastography an imaging technique that measures the stiffness of soft tissues using acoustic shear waves and imaging their propagation using MRI.
- Noise floor is the measure of the signal created from the sum of all the noise sources and unwanted signals within a measurement system
- PET Positron Emission Tomography is a functional imaging technique that produces a three-dimensional image of functional processes in the body using a positron-emitting radiotracer.
- Phase coherence When referring to multiple measurements within a common VOI of a or multiple k- values indicates that the sample has the same position relative to the measurement frame of reference
- Phase encode A phase encode is used to impart a specific phase angle to a transverse magnetization vector.
- the specific phase angle depends on the location of the transverse magnetization vector within the phase encoding gradient, the magnitude of the gradient, and the duration of the gradient application.
- Phase wrap The helical precession of the phase of the transverse magnetization along a phase encoded sample
- Semi-crystalline texture a texture exhibiting regular spacing along one or more directions
- slice Used interchangeably to indicate a non-zero thickness planar section of the Slice-selective refocusing Refocussing of spins through combination of a slice selective gradient and an RF pulse such that the bandwidth of the RF pulse selects a thickness along the direction of the gradient, and the RF pulse tips the net magnetization vector away from its equilibrium position Only those spins processing at the same frequency as the RF pulse will be affected.
- T2 Defined as a time constant for the decay of transverse magnetization arising from natural interactions at the atomic or molecular levels.
- T2* In any real NMR experiment, the transverse magnetization decays much faster than would be predicted by natural atomic and molecular mechanisms; this rate is denoted T2* ("T2-star”). T2* can be considered an "observed” or “effective” T2, whereas the first T2 can be considered the “natural” or “true” T2 of the tissue being imaged. T2* is always less than or equal to T2.
- TBS Trabecular Bone Score is a technique that looks for texture patterns in the DEXA signal for correlation with bone microarchitecture for assessing bone health
- Echo Time is the time between the 90° RF pulse and MR signal sampling, corresponding to maximum of echo.
- the 180° RF pulse is applied at time TE/2.
- Repetition Time is the time between 2 excitations pulses (time between two 90° RF pulses).
- Textural frequency the number of texture wavelength repeats per unit length in a texture Texture wavelength the characteristic spacing between structural elements in a texture TR Spin Echo sequences have two parameters: Echo Time (TE) is the time between the 90° RF pulse and MR signal sampling, corresponding to maximum of echo. The 180° RF pulse is applied at time TE/2. Repetition Time is the time between 2 excitations pulses (time between two 90° RF pulses).
- TE Echo Time
- Repetition Time is the time between 2 excitations pulses (time between two 90° RF pulses).
- Vector combination gradient A magnetic gradient resulting from any vector combination of the gradient coil set
- window function In signal processing, a window function (also known as an apodization function or tapering function) is a mathematical function that is zero -valued outside of some chosen interval
- x-ray diffraction is a tool used for identifying the atomic and molecular structure of a crystal
- the embodiments disclosed herein provide an MR-based technique that enables in vivo, non- invasive, high-resolution measurement and assessment of fine biologic textures, enabling monitoring of texture formation and/or change in response to disease onset and progression in a range of pathologies.
- This same method can be applied to fine- texture characterization in other biologic and physical systems. It enables MR-based resolution of fine textures to a size scale previously unattainable in in vivo imaging.
- the method while described herein with respect to biological systems for examination of tissue, is equally applicable for assessment of fine structures in a range of industrial purposes such as measurement of material properties in manufacturing or in geology to characterize various types of rock, as well as other uses for which measurement of fine
- the acquisition TR can be repeated as many times as necessary, changing the encode as needed to span the desired extent of real and of k-space required.
- the set of one or more k-values output from each TR are now high SNR due to the ability to average repeats without motion effects, and since the measure of interest is textural spacing, and not development of an image, the lack of phase coherence between TRs is of no concern.
- the method claimed herein consists of acquiring MR signal from within an inner volume to encompass a specific tissue region of interest, such as a lesion, an organ, a location in an organ, a specific region of bone, or a number of regions in a diseased organ for sampling.
- a specific tissue region of interest such as a lesion, an organ, a location in an organ, a specific region of bone, or a number of regions in a diseased organ for sampling.
- This inner volume may be excited by one of a number of methods, including but not limited to: intersecting slice-selective refocusing, selective excitation using phased-array transmit in combination with appropriate gradients, adiabatic pulse excitation to scramble signal from the tissue outside the region of interest, outer volume suppression sequences, and other methods of selectively exciting spins in an internal volume including physically isolating the tissue of interest, to name a few, [PARA 40] After definition of a volume of interest (VOI), the gradient is turned off, and multiple samples of signal centered at a specific k-value, the spread of which is defined by receiver BW and sampling length, are acquired.
- VOI volume of interest
- This measurement is repeated only in specified directions within the VOI rather than trying to map all of k-space, as is required to generate an image.
- One or more samples of a particular k-value are acquired within an acquisition block during a single TR and the k-value subsequently incremented or decremented, allowing further multiple samples of other k-values as desired during the same TR.
- This method allows multiple sampling of each k-value of interest over a time period of milliseconds, providing immunity to subject motion. The process can then be repeated in further TRs, the requirement on motion between k-value acquisitions being only that the VOI remain within the tissue region of interest. Build up of a magnitude spectrum of spatial frequencies may be accomplished without the need to acquire it in a spatially coherent manner.
- the quantities of interest are the relative intensities of the various k-values (textural spacings) present in the sample volume, as long as the acquisition volume remains within a representative sample of tissue, any motion between the blocks does not compromise the measurement.
- use of fairly robust, real-time piloting and acquisition algorithms can be used for gross repositioning of the internal selectively excited volume and for rejecting data sets that have failed to stay in the proper tissue.
- the data acquired can, with reference to positioning images, be mapped spatially.
- Either the VOI can be moved in successive TRs or interleaved acquisition done within a single TR by exciting additional volumes during the time that the signal is recovering in advance of the next TR. The requirement is that successive VOIs be excited in new tissue, that does not overlap the previous slice selects. Spatial variation of pathology can be determined by this method. This can also be used to monitor temporal progression of a pathology through an organ if the measure is repeated longitudinally.
- ASL Arterial Spin Labeling
- DSC Dynamic Susceptibility Contrast
- PARA 46 The method claimed herein can also be used in conjunction with other MR- based measurement techniques, including DWI and DTI, to provide front end information toward parameter selection for the diffusion techniques as well as correlation with their measurements of tissue health.
- the echo record window can be designed such that recording begins with the highest k-values of interest, as signal level is highest at the echo peak. This enables recording of fine structures currently unachievable with in vivo MR imaging.
- T2* is longer with the gradient off, so SNR is improved by the longer acquisition times possible This allows acquisition of an increased number of samples, N.
- [PARA 52] Coil combination is also simplified by having higher SNR for each k-value, hence providing a significant improvement in overall SNR. This is especially beneficial as the trend in MRI is towards coil arrays composed of many small element coils. As the acquisition volumes targeted in the method claimed herein are small, correction for phase across the sample volume is not needed. Only one phase and gain value for each coil is needed for combining the multiple element channels. These can be combined using the Maximal Ratio Combining (MRC) method, which weights the coil with the highest SNR most heavily, or other multi-signal combination methods.
- MRC Maximal Ratio Combining
- Phase and gain for the elements of a given coil array can be determined once from a phantom and applied to patient data.
- PARA 53 Signal acquisition and data sampling in a standard MRI scan is done by acquiring complex- valued samples of multiple echoes, while applying a gradient sequence concurrently, as well as in sequence with the echoes. Imaging relies on frequency encode for one of the dimensions because this allows a line in k-space to be acquired with each phase encode rather than a single point. For 3 -dimensional imaging, two dimensions in k- space normally rely on phase encode to generate the targeted filling of k-space, with the third dimension frequency-encoded.
- Phase encode acquisition in imaging usually entails acquisition of on the order of 256 k-values in each of the phase-encode directions, hence is is a relatively slow process.
- Clinical MRI scans take on the order of 10-15 minutes to generate an image.
- the aim in image construction is to acquire sufficient k-space coverage to fill out all the coefficients in the 2 or 3-dimensional Fourier series, which is why in standard MR resolution is limited by subject motion.
- the method claimed herein is in direct contrast to standard MR data acquisition, with its focus on image generation. Image formation is plagued by blurring resulting from subject motion over the long time necessary to acquire the large data matrix required. Since the target of the method claimed herein is texture rather than image, the only requirement on subject motion is that the sampled volume remain within a region of similar tissue properties over the course of acquiring data. This is a much less stringent and easy to achieve target than the requirement of structural phase coherence, as the scale of the allowable motion is then large enough, and of a temporal order, to be easily correctable by real-time motion assessment and correction techniques. The speed of acquisition for the method claimed herein is such that, in most cases, real-time motion correction may not be necessary at all. While other methods have focused on post -processing of images to try to extract textural measures, the method claimed herein eliminates the need for image generation, focusing instead on directly measuring texture, hence enabling a more sensitive and robust measure.
- k-space sampling is considered synonymous with sampling of an echo in the presence of a gradient set.
- the approach to k-space filling is to acquire only the set of k-values needed for texture evaluation in the targeted pathology, with data acquired after the gradient is switched off. This method enables such rapid acquisition of single-k- value repeats for averaging for noise reduction that subject motion does not degrade the data.
- MR echo sampling provides specific samples vs. time of a time-dependent echo.
- the echo is comprised by the gradients applied concurrently (for the frequency- encode axis) and prior to (for a phase encode axis), but also contains the isochromats associated with the different chemical species of the sample, as well as the envelope (T2 & T2*) associated with spin-spin interactions.
- the gradient slope is typically made as steep as possible to minimize the apparent shift to within a narrow range (i.e. within 1 or two pixels in the image domain).
- a narrow range i.e. within 1 or two pixels in the image domain.
- this then requires a larger receiver bandwidth to accommodate the larger frequency range. This in turn increases the overall noise floor at a level proportional to the receive bandwidth.
- PARA 64 A common technique for noise reduction in signal acquisition is through repeat sampling of a signal and subsequent combination of the multiple measurements. For linear noise sources, such as Gaussian noise, this technique improves SNR through cancellation of the random noise on the signal, the cancellation effect increasing with the number of samples, N.
- [PARA 76] Within a given acquisition in standard MR practice, there are M samples which are acquired of the echo. Instead of acquiring a sample at each k-value, N ⁇ M of those samples could be used for estimation of the (complex-valued) underlying signal value at a specific k-value. Multiple samples within an acquisition can be combined with much less concern of movement than across acquisitions because they are much closer in time. [PARA 77] If the entire echo is used to measure one k- value, the receive bandwidth can be adjusted so as to pass the most abundant resonant peaks in the underlying NMR spectrum, and attenuate frequencies above them.
- the claimed method enables acquisition of values in regions of k-space which have very low signal levels, such as would be found for higher k- values (shorter textural-wavelengths)— the fine texture range that has hitherto remained elusive.
- the full NMR spectrum may be extracted (without any phase encoding gradients: just volume selection) to obtain a baseline of the underlying signal strength (and associated frequencies), which in turn will be spatially modulated, providing insight into textural wavelengths through knowledge of the chemical species expected in the textural elements under study.
- the isochromats of interest can be extracted by acquiring N samples of the echo, then taking the Fourier transform. Since the echo is being played out with no gradient, the strength of the resulting signal at the Isochromat of interest will correspond to the (complex- valued) k- value coefficient of interest.
- ⁇ represents the noise variance
- ' ' represents the squared magnitude of the isochromat(s) of interest
- 0 ⁇ £ ⁇ 1 represents the allowable error of the estimate.
- NF e ff is the effective noise figure of the receiver
- T is the temperature in Kelvin of the biological sample
- B is the receiver bandwidth.
- N can be used as a guide to the number of samples that need to be acquired within a given acquisition in order to create a reasonable estimate.
- FIG. 3 shows an example timing diagram for a pulse sequence for data acquisition using the method claimed herein.
- RF pulses included in trace 302 are employed to excite selected volumes of the tissue under investigation, as in typical MR imaging.
- a first RF pulse, 304 is transmitted coincidentally with a gradient pulse 308 on the a first magnetic field gradient, represented in trace 306 .
- This excites a single slice, or slab, of tissue the positioning of which is dependent on the orientation and magnitude of the first gradient, and the frequencies contained in the RF pulse.
- the negative gradient pulse, pulse 310 rephases the excitation within the defined thickness of the slice or slab.
- this second RF pulse 312 tips the net magnetic vector to antiparallel to Bo, it results in inversion of spins and subsequent refocusing, thus leading to a signal echo at a time after the 180 degree RF pulse equivalent to the time between the 90° and 180° RF pulses.
- An initial higher value gradient pulse, 318, at the start of gradient pulse 316 is a crusher, or "spoiler" gradient, designed to induce a large phase wrap across the tissue volume.
- the second RF pulse in combination with the applied second gradient, provides slice selective refocusing of the signal in a region defined by the intersection of the first slice and the second slice set by this second gradient.
- An encoding gradient pulse 326 on trace 314, sets an initial phase wrap, hence k-value encode, along the direction of gradient pulse 326.
- the k-value encode can be oriented in any direction, by vector combination of the machine gradients but for ease of visualization is represented as on the second gradient.
- the negative prephasing gradient pulse 326 winds up phase such that, at the signal echo following the second 180° RF pulse, signal acquisition starts at high k-value, which may then be subsequently decremented (or incremented or varied in orientation) for further acquisitions, as will be described below.
- energy density in the signal is generally proportional to k "1 , this method ensures k-values with lower SNR are acquired first, before T2 effects have caused much overall signal reduction.
- a receive gate 333 is opened to receive the RF signal, which is shown in FIG. 3 as pulse 334 on trace 336.
- the RF signal in trace 336 is a representation showing only the signal present in the receive gate window without showing the actual details of the RF signal outside the window.
- Sampling occurs as represented by trace 338 beginning with the initial k-value, 340a, seen on trace 324. Note that, at the scale of the drawing, the sampling rate is high enough that the individual triggers of the analog to digital converter (A/D) have merged together in trace 338. (The expanded time scale in FIG. 4 described below shows the individual A/D triggers.)
- a k-value selection gradient pulse 342a is then applied and the resultant k-value 340b is sampled. (Though shown in the figure as a negative pulse on the second gradient, decrementing the k-value, in practice this pulse and subsequent k-value gradient pulses can be designed through any vector combination of gradients to select any k-value or orientation.)
- the k-value selection gradient pulse 342b selects a third k-value 340c which is sampled by the A/D. Each gradient pulse changes the phase wrap, selecting a new k-value.
- Application of a k-value selection gradient pulse (342c - 342f) followed by multiple sampling of the resultant k-value coefficient is repeated as many times as desired.
- the gradient orientations for slice and k-value select may be coincident with the machine gradients, which are aligned to lie coincident or orthogonal to the Bo field.
- the acquisition directions and k-value encodes may be selected using gradients that are a vector combination of all three machine-gradient axes.
- the prewinding encoding gradient pulse can be set such that the fist k-value to be measured is the desired low SNR k-value.
- the prewinding gradient pulse can be set to zero so that the first k-value measured is kO.
- a measurement of kO may be desired for the purpose of determining the systems receiver sensitivity to the particular VOI, determining the relative prevalence of isochromats (e.g., water vs.
- a strategy for gathering a specified set of k-values for a VOI may include measuring the low SNR k-values (typically the higher k-values) in a first set of multiple TR and then measuring kO and other higher SNR k-values in other TRs while remaining in the same VOI.
- the signal magnitude reaches a maximum at the time of the spin echo. It is also shown diagrammatically that the signal magnitude is varying throughout the acquisition of the multiple RF measurements of a k- value and more so between successive blocks of measurements of k-values. Alignment in time of the measurement of the low SNR k-values with the highest echo signal enhances the SNR of the k-value measurement, alternatively alignment of higher SNR k-values with lower echo signal magnitude allows gathering additional useful k-value acquisitions during the echo.
- FIG. 4 shows a close-up of the pulse sequence of FIG. 3 during the initial portion of the RF sampling window 338 between 7.25 and 8.00 msec. Multiple samples of the same k-value, taken in rapid succession with all gradients off, provide the input for signal averaging to reduce AWGN when SNR is low.
- a first block 344a of the sampling window 338 multiple samples 346a are taken of the first k-value 340a.
- transition samples 348a are taken.
- multiple samples 346b are taken at the second k- value 340b.
- k-value selection gradient pulse 342b then occurs with associated transition samples 348b, and subsequent acquisition of samples 346c of the third k-value 340c after the gradients are switched off.
- the underlying signal is minimally impacted by motion due to the very short time window used to acquire data at each given k- value. Since the data is acquired with gradients off, there is no issue with chemical shift and the effective T2 is longer, boosting the signal value.
- a consistent number of samples at each k-value can be acquired, or an alternative sequence may be employed where, as k-values decrease, hence increasing in signal amplitude, fewer samples are acquired.
- a pulse sequence designed for this type of acquisition is illustrated in FIG. 5. Multiple samples of each k-value targeted in the acquisition are acquired in rapid succession, with all gradients off. These repeats provide the input for signal averaging in low SNR signals.
- the underlying signal is minimally impacted by motion due to the very short time window in which data is acquired for a given k-value.
- Samples within the portions of the sample window 344a - 344g outlined on FIG. 5 correspond to the number of samples acquired for a given k-value 340a - 340g each induced by an unwinding pulse 342a - 342f of the k-value selection gradient.
- Nk the number of samples associated with a given k-value, can be selected based upon expected SNR, tissue contrast, contrast to noise, pathology, texture size, and/or texture bandwidth. For the example in FIG. 5 it can be seen that a decreasing number of samples are taken for progressively smaller k-values (larger textural features).
- [PARA 100] Refocusing the echo, and/or a new TR can be used to build up a library of k- space samples. Acquisition of multiple k-values within one TR can be facilitated by application of multiple refocusing gradients and/or RF pulses, to increase the time over which the additional k-values can be sampled within a TR. These later echoes would presumably be used to acquire the coefficients of the lower k-values in the selected set, as their energy density in the continuum of values is generally higher so the effect of T2 decay on overall signal magnitude will not affect them as severely as it would the higher k-values.
- FIG. 8 shows an extension of the basic sequence of the method claimed herein, using spin-echo refocusing to extend the record time for the TR.
- Application of a refocusing RF pulse 802 with an associated gradient pulse 804 results in slice-selective refocusing.
- a second sampling window 806 is opened by the receive gate 808.
- Multiple k- value selection gradient pulses 810 are applied to increment the selected k- value and, after switching off each successive gradient pulse, multiple samples of the selected k- value are acquired in the sampling window.
- a second slice-selective refocusing RF pulse 812 with associated gradient pulse 814 again inverts the spins and, after application of each in the multiplicity of k- value selection gradient pulses 820, data is acquired in a third sampling window 816, opened by the receive gate 818.
- a third sampling window 816 opened by the receive gate 818.
- an increasing number of k-values may be sampled with each refocusing. Refocusing can be repeated until the decrease in signal level from T2 and other effects makes further signal acquisition ineffective.
- Another method to extend the record time by exciting multiple signal echoes is to use one, or a series of, gradient recalled echoes (GRE). GRE are different from the SE in that they cannot refocus the effects of stationary inhomogeneities, so T2* effects limit the number of repeats.
- GRE gradient recalled echoes
- the k-values associated with particular pathology will be part of the determination of the number of samples needed for signal averaging, Nk.
- the wavelength of pertinent textures is in the range of 400 microns, i.e. a k-value of 2.5 cycles/mm. This is similar to the textural spacings seen in fibrotic development in many other diseases, such as cardiac fibrosis.
- the spacing of elements in trabecular bone varies a lot, but the minimum spacing of interest is the width of trabecular elements, which are approximately 80 microns, setting a maximum k-value of 12.5cycles/mm.
- many of the textures of interest are very fine, on the scale of 50 microns, equivalent to a k-value of 20 cycles/mm.
- Each pathology will dictate what exactly is needed as quantitative data, i.e. what part of the continuum of k-values needs to be monitored, and with what resolution and sensitivity.
- short (long) wavelength features increase at the expense of long (short) wavelength features (e.g. liver fibrosis).
- an amplitude decrease and broadening of short wavelength features indicates disease progression— e.g. degradation of the ordered formation of cortical neuronal minicolumns (approximately 80- micron spacing) with advancing dementia.
- cortical neuronal minicolumns approximately 80- micron spacing
- the signal level obtainable will depend on anatomy to some extent. For instance, though the resolution needed is highest in brain, the proximity of the cortex to the surface of the head ensures that use of a surface coil will provide significant signal boost for cortical structures. Lower resolution is required in liver, as the structures of interest are on the order of several hundred, rather than tens of microns. But, the organ is deeper (further from the coil) reducing the measured signal. Using the in-table coil for spine data acquisition yields modest signal level and good stabilization. Also, bone is a high contrast target, so the SNR requirement is not as stringent. For all these reasons, the exact number of repeats needed for averaging depends on more than the k-value range targeted.
- FIG. 6 shows a simulation demonstrating that the ability provided by the claimed method to acquire many repeats of a targeted k-value within a single TR enables robust signal averaging to boost SNR.
- a subject displacement rate which has in practice been measured clinically over the course of several scans
- the situation is actually worse due to the fact that to acquire 90 repeats using conventional spatially-encoded echoes would require several TRs, making the acquisition time significantly longer, and the signal degradation due to motion much more severe.
- the potential SNR gain due to multiple sample combination has been nullified by the effects of motion.
- [PARA 108] This is very similar to x-ray diffraction, wherein the magnitude-only information (no phase) obtained presents the challenge of determining a best estimate of the corresponding structural profile based on this magnitude-only information. Algorithms exist towards solving the problem, the chance of success depending on the range of k-value coefficients obtained, the SNR of each averaged coefficient, and the width of values contained in a nominally single-valued acquisition of k-value. The chance of success in this effort is greatly increased using the claimed method due to its immunity to subject motion. [PARA 109] The ability to reconstruct a profile from k-value data depends somewhat on the spectral broadness of each single-k- value acquisition. While this is influenced by the VOI (Volume of Interest) size and shape, it is also influenced by k-value and pathology, as degradation of tissue often tends to lead to more textural randomness within tissue.
- VOI Volume of Interest
- VOI shape, dimensions, orientation, and positioning within an organ/anatomy affects the data measured and its interpretation.
- the VOI shape can be chosen to maximize usefulness of the acquired data.
- Data can be acquired in different directions, and at different textural wavelengths (k-values) within a VOI enabling assessment of textural anisotropy. Texture can be sampled in multiple VOIs, either interleaved within a single TR, or in successive TRs, towards assessment of pathology variation across an organ.
- Standard interleaving processes for the VOI may be used within a TR to provide additional data by applying additional encoding pulses on vector combination gradients and associated k-value selection gradient pulses for k-values in the interleaved VOI.
- additional excitation RF pulses with associated slice selection gradients may be repeated within the same TR by exciting a volume of interest with a gradient set in each repeat having at least a first gradient with an alternative orientation from the first gradient pulse 308 applied initially in the TR, to define an additional VOI for excitation in new tissue, that does not overlap any previous VOI in the TR (fourth, fifth and sixth gradients in a first repeat and succeeding incremental gradients in subsequent repeats).
- This response can be mapped, or the several measures taken and averaged, whatever is appropriate for the targeted pathology. This is similar to the multi- positioning of tissue biopsy. However, in the case of tissue biopsy, the number of repeats is limited due to the highly invasive nature of the technique. The minimum number of structural oscillations to be sampled at a specific k-value dictates a minimum VOI dimension in the direction of sampling— the length required varying inversely with targeted k-value.
- the VOI dimension in the sampled direction can be held constant for all k-values in the targeted range, with the result that the number of structural oscillations sampled will vary with k-value. This is a simple solution, requiring the sampling dimension be set by the lowest k-value (longest wavelength structure). Using this approach, the sampling dimension of the VOI is larger than required for the highest k-value in the range, thus providing less localization within the tissue than would be otherwise possible.
- the VOI may be held constant and the vector combination gradient for the encoding and k-selection pulses may be altered from TR to TR for assessing feature size.
- the VOI may be held constant and the vector combination gradient for the encoding and k-selection pulses may be altered from TR to TR for assessing feature size.
- the elliptical cross-section VOI 902 in FIG. 9B is one such possibility. Acquisition along any radial direction, as well as along the axis of the shape, would be possible within one TR.
- k-value encodes can be applied in multiple directions by changing the applied vector combination gradients for encoding and k-selection pulsing.
- the exact form of the VOI and sampling direction can be used to yield much textural information.
- the organization of cortical neuron fiber bundles is semi-crystalline, as the bundles in healthy tissue form in columns. Because of this, the measure of textural spacing perpendicular to the bundles is very sensitive to orientation. When the orientation is exactly normal to the columns, a very sharp signal maximum is expected, the signal falling off rapidly as the orientation varies on in either rotational direction away from this maximum.
- One way to measure the spacing and organizational integrity would be to "rock" the acquisition axis around this maximum looking for a resonance in signal intensity. This approach of looking for "textural resonances" by looking for signal maxima can be applied in any tissue region. As pathology degrades the organizational integrity, the sharpness of this peak will degrade and the signal maximum will be reduced.
- the randomness of the spacing in certain textures can be assessed by varying the length of tissue sampled in a specific, or in multiple directions, with subsequent change in acquisition length.
- the selected value for that length can be varied over multiple TRs to test the sensitivity of the measured coefficient to this parameter.
- the VOI can be selectively excited by a number of methods, for instance intersecting slice-selective refocusing, selective excitation using phased-array transmit in combination with appropriate gradients, adiabatic pulse excitation to scramble signal from the tissue outside the region of interest, as examples. Parameter selection for the various methods can be done with SNR optimization in mind. For instance, the VOI generated by a slice selective excitation and two additional mutually- orthogonal slice selective refocusing pulses, as by VOI 904 in FIG. 9A. Through careful RF pulse design, the shape of the VOI can be designed so that the edges are smooth and more approximate a windowing function, as shown in FIG. 9B.
- the VOI can be moved from place to place within an organ or anatomy under study to measure the variation of texture/pathology. This response can be mapped, or the several measures taken and averaged, as appropriate for the targeted pathology. This is similar to the multi-positioning of tissue biopsy. However, in the case of tissue biopsy, the number of repeats is limited due to the highly invasive nature of the technique.
- k-space is probed to reveal texture in such a way as to eliminate the loss of signal resolution that results from subject motion blurring.
- the focus here is on acquisition of a select few k-values per TR, with sufficient repeats of each to yield high SNR.
- Each of the individual acquisitions is centered on a single k- value.
- the spatial encode is, to first order, a single spatial frequency sinusoidal encode, there are a number of factors which have the effect of broadening the spatial frequency selectivity of the k-value measurement.
- One significant factor affecting the broadness, or bandwidth, of the k-value measurement is the length of the sampled tissue region.
- an aspect of the claimed method is the ability to set the bandwidth of the k-value measurements by appropriate selection of the sampling length within the VOI. Using this method, the bandwidth of the measurement can be set according to the desired k-space resolution appropriate to the tissue being evaluated.
- [PARA 129] As an example, diagnosis and monitoring of ischemic stroke has been greatly facilitated by combining DWI-MRI with standard structural MR imaging. Current studies link MR functional contrast from Blood Oxygenation Level Dependent (BOLD) imaging with both structural MR images and DWI-MRI images. Structural MR image segmentation techniques are used to assess cortical tissue degradation; this information is then combined with both glucose metabolism and amyloid plaque deposition information derived from PET studies. PET is combined with anatomical MRI and with CT in cancer diagnosis. [PARA 130] The increasing range of techniques/contrast mechanisms available in MRI enables combining different diagnostic measurements/techniques within one exam, facilitating data comparison.
- BOLD Blood Oxygenation Level Dependent
- DWI a measure of the voxel-averaged rate of diffusion of water molecules in the target tissue, has increased in use significantly in recent years, as has the derivative technique DTI which assesses the directional anisotropy in the diffusion process.
- the main clinical domain of application of diffusion imaging is neurologic disorders, especially for management of patients with acute stroke. Additionally, in conjunction with structural MRI, diffusion imaging is routinely used to diagnose Creutzfeld- Jakob Disease (CJD), diffuse axonal injury, evaluate white matter damage in MS, aid in the diagnosis of tumor metastases (breast, prostate, liver), and evaluate tissue changes associated with dementia onset and as an in vivo measure of brain circuitry.
- CJD Creutzfeld- Jakob Disease
- MS diffuse axonal injury
- evaluate white matter damage in MS aid in the diagnosis of tumor metastases (breast, prostate, liver), and evaluate tissue changes associated with dementia onset and as an in vivo measure of brain circuitry.
- [PARA 132] Diffusion imaging provides quantitative information dependent on the voxel- averaged microscopic properties of tissues, such as cell size and shape, packing, etc.
- DWI/DTI is a quantitative measure, it is not a direct measure of these microscopic properties, but is an average over a voxel of several cubic mm.
- ADC Average Diffusion Coefficient
- the diffusion signal is low and is extremely sensitive to subject motion. Long and high gradients are required to create the diffusion weighted signal, resulting in a long TE, with significant T2 and T2* decay during signal readout, hence low SNR. The strength and duration of the gradients result in high motion sensitivity, as any motion in the gradient direction during the time of gradient application will look like diffusion. (LE BIHAN, D., JOHANSEN-BERG, H.; "Diffusion MRI at 25: Exploring Brain Tissue Structure and Function"; Neuroimage, Vol. 61, 2012.) The signal is sensitive to diffusion time— a long enough diffusion time that water molecules have time to interact with the tissue structure is required.
- ALEXANDER A., et al.; "Diffusion Tensor Imaging of the Brian”; Neurotherapeutics 4(3), 2007.
- MASCALCHI M., et al.; "Diffusion-weighted MR of the brain: methodology and clinical application”; Radio logia Medica, 109(3) 2005.
- the method claimed herein offers a direct measure of very fine (near cellular) texture. As such, it provides valuable complementary diagnostic information as it measures very fine structure, and DWI/DTI measure the effects of that structure on diffusion. Along with the ability of the method claimed herein to provide intricate textural information towards pathology assessment, when used at the front- end of an MR exam that includes DTI it can provide information on the structural boundaries/hindrances in any direction, clarifying interpretation of the DTI signal.
- the method claimed herein enables robust measure of very small scale displacements— it could provide a measure of motion occurring over the course of the DTI acquisition, towards post processing reregistration.
- the different sequences could, perhaps, be intermingled in time to optimize the effect.
- Some possible targets for combining the method claimed herein with DTI would be in assessing vasculature in tumors, measuring changes in gray matter with dementia onset. (WESTON, P., et al.; "Diffusion imaging changes in grey matter in Alzheimer's disease: a potential marker of early neurodegeneration"; Alzheimer's Research and Therapy, 7: 47, 2015.), assessing hyperplasiac development of mammary ducting in breast cancer progression, and diagnosis and monitoring of MS, to name a few.
- the method claimed herein can be applied in tissue regions showing abnormalities in DWI, to measure textural variations indicative of pathology. Measurement of the variation in texture as a function of direction can be further used as an input to guide application of DTI, for correlation with and verification of the the fractional anisotropy measure of DTI.
- ASL Arterial Spin Labeling
- CBF cerebral blood flow
- CVD cerebrovascular disease
- PETCHARU PAISAN S., et al.; "Arterial spin labeling in neuroimaging”; World Journal of Radiology, 2(10), 2010.
- the method claimed herein could provide information on vasculature pathology when used as a complementary technique to ASL-acquired data or in conjunction with ASL contrast.
- the method claimed herein would use the magnetically labelled blood to heighten contrast, facilitating measurement of changes in vessel texture (vessel density and spatial randomness) as an indicator of pathology such as CVD, tumor, or the vasculature changes attendant with dementia progression.
- PARA 138 Having now described various embodiments of the invention in detail as required by the patent statutes, those skilled in the art will recognize modifications and substitutions to the specific embodiments disclosed herein. Such modifications are within the scope and intent of the present invention as defined in the following claims.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
Claims
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2015312190A AU2015312190A1 (en) | 2014-09-01 | 2015-08-31 | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies |
EP15837428.0A EP3189345A4 (en) | 2014-09-01 | 2015-08-31 | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies |
JP2017530978A JP6553723B2 (en) | 2014-09-01 | 2015-08-31 | Selective sampling magnetic resonance based method for assessing structural spatial frequencies |
CN201580053403.5A CN106716167B (en) | 2014-09-01 | 2015-08-31 | The method based on selective sampling magnetic resonance for evaluation structure spatial frequency |
CA2959621A CA2959621A1 (en) | 2014-09-01 | 2015-08-31 | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies |
AU2021201998A AU2021201998B2 (en) | 2014-09-01 | 2021-03-31 | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies |
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201462044321P | 2014-09-01 | 2014-09-01 | |
US62/044,321 | 2014-09-01 | ||
US201462064206P | 2014-10-15 | 2014-10-15 | |
US62/064,206 | 2014-10-15 | ||
US201562107465P | 2015-01-25 | 2015-01-25 | |
US62/107,465 | 2015-01-25 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2016036643A1 true WO2016036643A1 (en) | 2016-03-10 |
Family
ID=55402215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2015/047676 WO2016036643A1 (en) | 2014-09-01 | 2015-08-31 | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies |
Country Status (7)
Country | Link |
---|---|
US (1) | US9366738B2 (en) |
EP (1) | EP3189345A4 (en) |
JP (1) | JP6553723B2 (en) |
CN (1) | CN106716167B (en) |
AU (2) | AU2015312190A1 (en) |
CA (1) | CA2959621A1 (en) |
WO (1) | WO2016036643A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022155566A1 (en) * | 2021-01-18 | 2022-07-21 | Bioprotonics, Inc. | Magnetic-resonance-based method for measuring microscopic histologic soft tissue textures |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6362950B2 (en) * | 2014-07-24 | 2018-07-25 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Magnetic resonance equipment |
US10061003B2 (en) * | 2014-09-01 | 2018-08-28 | bioProtonics, L.L.C. | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms |
US9664760B2 (en) * | 2014-09-01 | 2017-05-30 | bioProtonics, L.L.C | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms |
US10955503B2 (en) * | 2014-09-01 | 2021-03-23 | Bioprotonics, Inc. | Method to measure tissue texture using NMR spectroscopy with VOI length in an analysis direction defined by receiver bandwidth |
US10215827B2 (en) * | 2014-09-01 | 2019-02-26 | Bioprotonics Inc. | Method to measure tissue texture using NMR spectroscopy to identify the chemical species of component textural elements in a targeted region of tissue |
US10085703B2 (en) * | 2015-01-27 | 2018-10-02 | Septimiu Edmund Salcudean | Dynamic computed tomography imaging of elasticity |
GB201505940D0 (en) | 2015-04-08 | 2015-05-20 | Isis Innovation | Medical imaging |
KR102449249B1 (en) * | 2015-05-27 | 2022-09-30 | 삼성전자주식회사 | Magnetic resonance imaging apparatus and method thereof |
EP3359034A4 (en) * | 2015-10-07 | 2019-08-07 | Bioprotonics, Inc. | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms |
DE102017210700A1 (en) * | 2017-06-26 | 2018-12-27 | Robert Bosch Gmbh | Method for automated quantification of an analyte and NMR measuring device for carrying out the method |
WO2019010381A1 (en) * | 2017-07-06 | 2019-01-10 | bioProtonics, LLC | A method to measure tissue texture using nmr spectroscopy to identify chemical species of component textural elements in a targeted region of tissue |
EP3517988A1 (en) * | 2018-01-29 | 2019-07-31 | Koninklijke Philips N.V. | Mr imaging using a stack-of-stars acquisition with intrinsic motion correction |
US10585154B1 (en) * | 2018-01-29 | 2020-03-10 | Quantum Valley Investment Fund LP | Nuclear magnetic resonance diffraction |
US10775466B2 (en) * | 2018-02-09 | 2020-09-15 | GE Precision Healthcare LLC | System and method for magnetic resonance imaging an object via a stochastic optimization of a sampling function |
EP3550319A1 (en) * | 2018-04-05 | 2019-10-09 | Koninklijke Philips N.V. | Emulation mode for mri |
CN108919355B (en) * | 2018-05-14 | 2020-04-07 | 中国海洋石油集团有限公司 | High-dimensional S transformation method based on structure tensor guidance |
WO2020158717A1 (en) | 2019-01-31 | 2020-08-06 | 富士フイルム株式会社 | Learned model, learning method, and program, and medical information acquisition device, method, and program |
CN112415452B (en) * | 2019-08-22 | 2024-03-19 | 西门子(深圳)磁共振有限公司 | Method and device for removing interference in signal, magnetic resonance system and storage medium |
CN111861236B (en) * | 2020-07-27 | 2023-10-27 | 东北大学 | Bar finishing process data synchronization method |
CN113156350B (en) * | 2021-03-25 | 2022-07-19 | 华南理工大学 | Non-invasive light opaque micron-sized living tissue magnetic resonance tomography method and system |
US20240219496A1 (en) * | 2021-04-30 | 2024-07-04 | Regents Of The University Of Minnesota | Systems and methods for radio frequency (rf) based spatial encoding in magnetic resonance imaging using frequency-modulated rf pulses |
CN113671427B (en) * | 2021-08-24 | 2022-07-15 | 上海交通大学 | Magnetic resonance elastography displacement acquisition method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4471306A (en) * | 1982-02-03 | 1984-09-11 | General Electric Company | Method of NMR imaging which overcomes T2 * effects in an inhomogeneous static magnetic field |
US7127092B2 (en) * | 2000-04-26 | 2006-10-24 | Mayo Foundation For Medical Education And Research | Reduction of motion artifact in NMR images using spherical navigator signals |
US20070167717A1 (en) | 2005-11-27 | 2007-07-19 | James Timothy W | Structure assessment using spatial-frequency analysis |
US7916144B2 (en) * | 2005-07-13 | 2011-03-29 | Siemens Medical Solutions Usa, Inc. | High speed image reconstruction for k-space trajectory data using graphic processing unit (GPU) |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60146138A (en) * | 1984-01-11 | 1985-08-01 | Hitachi Ltd | Examination apparatus using nuclear magnetic resonance |
IL78240A (en) * | 1986-03-24 | 1989-09-10 | Elscint Ltd | Spatially localized spectroscopy |
JPS63130056A (en) * | 1986-11-21 | 1988-06-02 | 旭化成株式会社 | Partial selective nuclear magnetic resonance data lead-out method |
JP3393895B2 (en) * | 1993-09-13 | 2003-04-07 | 株式会社東芝 | Magnetic resonance imaging |
US6134464A (en) * | 1997-09-19 | 2000-10-17 | General Electric Company | Multi-slice and multi-angle MRI using fast spin echo acquisition |
CN100451671C (en) * | 2002-11-18 | 2009-01-14 | 皇家飞利浦电子股份有限公司 | Magnetic resonance method and device |
US7218110B2 (en) * | 2005-09-02 | 2007-05-15 | Toshiba America Mri, Inc. | Method and apparatus for diffusion magnetic resonance imaging with the effects of eddy currents compensated |
US7903251B1 (en) | 2009-02-20 | 2011-03-08 | Acuitas Medical Limited | Representation of spatial-frequency data as a map |
DE102011007501B3 (en) * | 2011-04-15 | 2012-05-31 | Universitätsklinikum Freiburg | Method of magnetic resonance imaging for the selection and recording of curved layers |
WO2013040086A1 (en) * | 2011-09-13 | 2013-03-21 | Acuitas Medical Limited | Magnetic resonance based method for assessing alzheimer's disease and related pathologies |
WO2013049116A2 (en) | 2011-09-26 | 2013-04-04 | Acuitas Medical Limited | Method for detection of characteristics of organ fibrosis |
US9360542B2 (en) * | 2012-04-20 | 2016-06-07 | Wisconsin Alumni Research Foundation | System and method for spectrally-resolved three-dimensional magnetic resonance imaging without frequency-encoding gradients |
-
2015
- 2015-08-31 CN CN201580053403.5A patent/CN106716167B/en active Active
- 2015-08-31 US US14/840,327 patent/US9366738B2/en active Active
- 2015-08-31 WO PCT/US2015/047676 patent/WO2016036643A1/en active Application Filing
- 2015-08-31 JP JP2017530978A patent/JP6553723B2/en active Active
- 2015-08-31 AU AU2015312190A patent/AU2015312190A1/en not_active Abandoned
- 2015-08-31 CA CA2959621A patent/CA2959621A1/en active Pending
- 2015-08-31 EP EP15837428.0A patent/EP3189345A4/en not_active Withdrawn
-
2021
- 2021-03-31 AU AU2021201998A patent/AU2021201998B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4471306A (en) * | 1982-02-03 | 1984-09-11 | General Electric Company | Method of NMR imaging which overcomes T2 * effects in an inhomogeneous static magnetic field |
US7127092B2 (en) * | 2000-04-26 | 2006-10-24 | Mayo Foundation For Medical Education And Research | Reduction of motion artifact in NMR images using spherical navigator signals |
US7916144B2 (en) * | 2005-07-13 | 2011-03-29 | Siemens Medical Solutions Usa, Inc. | High speed image reconstruction for k-space trajectory data using graphic processing unit (GPU) |
US20070167717A1 (en) | 2005-11-27 | 2007-07-19 | James Timothy W | Structure assessment using spatial-frequency analysis |
Non-Patent Citations (3)
Title |
---|
PHILIPS HEALTHCARE: "Multiva Application Guide, Release 3.6, Volume 4 - MR Spectroscopy", TURBO SPECTROSCOPIC IMAGING (TSI, September 2012 (2012-09-01) |
See also references of EP3189345A4 |
SUSAN M. NOWOROLSKI ET AL.: "High spatial resolution 'H-MRSI and segmented MRI of cortical gray matter and subcortical white matter in three regions of the human brain", MAGNETIC RESONANCE IN MEDICINE, vol. 41, no. 1, 29 January 1999 (1999-01-29), pages 21 - 29, XP055686434, DOI: 10.1002/(SICI)1522-2594(199901)41:1<21::AID-MRM5>3.0.CO;2-V |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022155566A1 (en) * | 2021-01-18 | 2022-07-21 | Bioprotonics, Inc. | Magnetic-resonance-based method for measuring microscopic histologic soft tissue textures |
Also Published As
Publication number | Publication date |
---|---|
CN106716167A (en) | 2017-05-24 |
AU2021201998B2 (en) | 2023-01-12 |
US20160061917A1 (en) | 2016-03-03 |
CA2959621A1 (en) | 2016-03-10 |
CN106716167B (en) | 2019-03-12 |
AU2015312190A1 (en) | 2017-03-23 |
JP6553723B2 (en) | 2019-07-31 |
EP3189345A4 (en) | 2018-05-09 |
JP2017526505A (en) | 2017-09-14 |
AU2021201998A1 (en) | 2021-04-29 |
EP3189345A1 (en) | 2017-07-12 |
US9366738B2 (en) | 2016-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2021201998B2 (en) | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies | |
CA3000765C (en) | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms | |
US10330763B2 (en) | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms | |
US9664760B2 (en) | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms | |
US9664759B2 (en) | Method for assessing structural spatial frequencies using hybrid sampling with non-zero gradient for enhancement of selective sampling | |
Liney | MRI from A to Z: a definitive guide for medical professionals | |
Vargas et al. | Clinical neuroimaging using 7 T MRI: challenges and prospects | |
EP3629915A1 (en) | Selective sampling for assessing structural spatial frequencies with specific contrast mechanisms | |
WO2013116775A1 (en) | System and method for image processing with highly undersampled imaging data | |
Li et al. | Multi‐parametric MRI for radiotherapy simulation | |
Baskaran et al. | Myelographic MR imaging of the cervical spine with a 3D true fast imaging with steady-state precession technique: initial experience | |
Sammet | Magnetic resonance imaging (MRI) | |
Bihan-Poudec et al. | Cardiovascular effects on high-resolution 3D multi-shot diffusion MRI of the rhesus macaque brain | |
Stachera et al. | Heart and Brain in Magnetic Resonance Imaging | |
Jackson | Technical progress in neuroradiology and its application | |
Papanastasiou | Quantification of cardiac magnetic resonance imaging perfusion in the clinical setting at 3T | |
Coolen | Novel strategies for mouse cardiac MRI: better, faster, stronger | |
Lai | Respiratory self-gated whole-heart coronary magnetic resonance angiography | |
Koh et al. | Diffusion-weightedMRI: futuredirections |
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: 15837428 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2017530978 Country of ref document: JP Kind code of ref document: A |
|
ENP | Entry into the national phase |
Ref document number: 2959621 Country of ref document: CA |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2015312190 Country of ref document: AU Date of ref document: 20150831 Kind code of ref document: A |
|
REEP | Request for entry into the european phase |
Ref document number: 2015837428 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2015837428 Country of ref document: EP |