WO2016151232A1 - Système et procédé irm pour estimer un paramètre physiologique à partir de deux autres paramètres physiologiques estimés - Google Patents
Système et procédé irm pour estimer un paramètre physiologique à partir de deux autres paramètres physiologiques estimés Download PDFInfo
- Publication number
- WO2016151232A1 WO2016151232A1 PCT/FR2016/050626 FR2016050626W WO2016151232A1 WO 2016151232 A1 WO2016151232 A1 WO 2016151232A1 FR 2016050626 W FR2016050626 W FR 2016050626W WO 2016151232 A1 WO2016151232 A1 WO 2016151232A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- physiological parameter
- estimating
- magnetic resonance
- parameters
- physiological
- 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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
-
- 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
-
- 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/4816—NMR imaging of samples with ultrashort relaxation times such as solid samples, e.g. MRI using ultrashort TE [UTE], single point imaging, constant time imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/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
-
- 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
-
- 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
Definitions
- the present invention provides a system and method for estimating a physiological parameter from data resulting from the acquisition of medical images.
- the invention differs in particular methods known for its accuracy and speed of execution.
- the invention is based in particular on Magnetic Resonance Imaging techniques (also known by the abbreviation "MRI” or “Magnetic Resonance Imaging” - MRI - according to an Anglo-Saxon terminology). These techniques make it possible to quickly obtain valuable information on the organs of human beings or animals. This information is particularly crucial for a practitioner seeking to establish a diagnosis and make a therapeutic decision in the treatment of pathologies.
- Magnetic Resonance Imaging techniques also known by the abbreviation "MRI” or “Magnetic Resonance Imaging” - MRI - according to an Anglo-Saxon terminology.
- a Nuclear Magnetic Resonance imaging apparatus As illustrated by way of non-limiting example in FIGS. 1 and 2, is generally used.
- the latter can deliver a plurality of digital image sequences 12 of one or more parts of the body of a patient, such as, by way of non-limiting examples, the brain, the heart, the lungs.
- said apparatus applies a combination of high-frequency electromagnetic waves to the part of the body in question and measures the signal re-emitted by certain atoms, such as, by way of non-limiting example, hydrogen for magnetic resonance imaging. Nuclear.
- the apparatus thus makes it possible to determine the chemical composition and hence the nature of the biological tissues in each elementary volume, which is commonly called a voxel, of the imaged volume.
- the nuclear magnetic resonance imaging apparatus 1 is controlled by means of a console 2. A user can thus choose parameters 11 for controlling the apparatus 1. From information produced by said apparatus 1, a plurality of digital image sequences 12 are obtained from a part of a body of a human or an animal.
- the image sequences 12 may optionally be stored within a server 3 and constitute a medical file 13 of a patient. Such a file 13 may include images of different types, such as perfusion or diffusion images.
- the image sequences 12 are analyzed by means of a dedicated processing unit 4.
- Said processing unit 4 comprises means for communicating with the outside world to collect the images.
- Said means for communicating further allow the processing unit 4 to deliver in fine, by means of rendering means 5 proposing a graphic, sound or other rendering, to a user 6 of the analysis system, in particular a practitioner or a researcher, an estimate of one or more physiological parameters, possibly formatted in the form of a content, from the images obtained by Magnetic Resonance Imaging by means of a suitable human-machine interface.
- rendering means any device alone or in combination for rendering a representation, for example graphic, sound or otherwise, of an estimated physiological parameter, to the user 6 of a system imaging analysis by Magnetic resonance.
- Such rendering means 5 may consist, in a non-exhaustive manner, of one or more screens, loudspeakers or other man-machine interfaces.
- Said user 6, advantageously a practitioner, of the analysis system can thus confirm or invalidate a diagnosis, decide on a therapeutic action that he deems appropriate, to deepen research work ...
- this user 6 can parameterize the operation of the processing unit 4 or the restitution means 5 by means of parameters 16. For example, it can thus define display thresholds or choose the estimated parameters for which it wishes to have a representation, for example graphic.
- the processing unit 4 generally comprises processing means, such as a computer, for implementing an estimation method in the form of a program previously loaded in storage means cooperating with said processing means. More generally, the processing unit may consist of one or more microprocessors or microcontrollers and / or internal memories cooperating with said microprocessors or microcontrollers.
- processing unit can also be extended to any software resource of the operating system type, implemented by said hardware elements, which provides services to facilitate the management of the hardware resources of said processing unit for any application method implemented by it.
- the acquisition of data advantageously signals, by Magnetic Resonance Imaging, hereinafter called MRI
- MRI Magnetic Resonance Imaging
- the two-dimensional images obtained consist of pixels with a thickness, corresponding to the thickness of the section and called voxel.
- the signal S obtained using such an MRI acquisition depends on two types of parameters.
- such a signal S depends on physiological parameters, that is to say the magnetic properties of the tissue which are, for example:
- Tl longitudinal relaxation time (spin-network) has the longitudinal relaxation is the process returning the magnetization to equilibrium along the direction of the Bo magnetic field.
- Tl is the characteristic time for the establishment of magnetization when the sample is placed in the magnetic field or that which characterizes the return to equilibrium after an inversion.
- T1 is also the time interval corresponding to the recovery of 63% of the initial longitudinal magnetization;
- T2 spin-spin transverse relaxation time
- the transverse relaxation is the process of return to equilibrium, that is to say to zero, of a magnetization brought into the plane perpendicular to the magnetic field Bo. This magnetization decreases with a characteristic time T2.
- T2 is the time interval corresponding to the loss of 63% of the initial transverse magnetization since the cessation of the application of a radiofrequency;
- said signal S depends on acquisition parameters directly related to the imaging apparatus 1 by Nuclear Magnetic Resonance, said parameters therefore being applicable to all the voxels.
- acquisition parameters are for example:
- the echo time TE time interval between an excitation by means of a pulse and the occurrence of the MRI signal in response to said excitation
- the inversion time 77 time interval between two characteristic pulses of a sequence for specific acquisitions in the context of inversion-recovery MRI.
- the signal S can be defined according to the following proportionality relationship: PD [le- TR, T1 ⁇ - ⁇ / ⁇ 2
- T2 weighted image By cleverly and manually combining the acquisition parameters by means of a prior configuration of a Nuclear Magnetic Resonance Imaging apparatus, a user of the apparatus is able to obtain images or sequences of Tl weighted images. , T2, PD, or obscure certain types of tissue. Thus, the user 6 can influence the generation of images.
- the user chooses, for example, a small value of TR, the term dependent on T1 can then be neglected and finally the signal S is only and substantially dependent on the physiological parameter T2.
- T2 weighted image also known as "T2 weighted image”
- the signal S can be defined according to the following proportionality relation:
- the user when the user judiciously chooses the acquisition parameter T1, the user can thus generate images missing or omitting a certain type of fabric, as non-limiting examples. fat or tissue.
- the transverse relaxation time T2 is modified by the effects of magnetic field heterogeneity. Indeed, the magnetic field applied within the imaging apparatus is not perfect, the magnet inducing the magnetic field is not uniform. The transverse relaxation time is then called T2 *.
- the signal S can then be defined according to the following proportionality relation:
- FIGS. 3A, 3B and 3C present three examples of Tl-weighted T2-weighted and inversion-recovery maps or images obtained by a Nuclear Magnetic Resonance Imaging apparatus according to a user-defined choice of parameters. Such figures have different contrasts given, highlighting certain parts of the brain.
- the user has chosen acquisition parameter settings TR, TE and 77 so as to suppress the water-induced signal.
- the user using a suitable imaging analysis system, can generate different types of weighted images which, as mentioned above, make it possible to highlight different organs of interest.
- a user advantageously a practitioner, can then use these weighted images to establish a diagnosis, for example the location and / or characterization of a tumor.
- the user or practitioner must perform several acquisition sequences to obtain different types of weighted images and consequently different contrasts.
- the acquisition sequences are generally quite long, of the order of several minutes. The multiplication of the sequences considerably increases the duration of the examination and generates several negative consequences, such as, without limitation:
- methods advantageously implemented by a processing unit of a magnetic resonance imaging analysis system, have been developed to first proceed to the acquisition of images. from multi-echo spin echo sequences, then calculate and / or estimate the physiological parameter T2 (also known as "T2 mapping sequences") and then generate weighted images regardless of the value of the TE acquisition parameter.
- T2 physiological parameter
- the creation of weighted images is, however, limited to the creation of T2 weighted synthetic maps, since the acquisition parameter TR can not vary.
- analogous methods can be applied to estimate and / or calculate the physiological parameter Tl, in advantageously varying the parameter TR or the flip-flop angle ⁇ (also known as "Tl mapping sequences").
- Tl mapping sequences also known as "Tl mapping sequences”
- T2 * the physiological parameter
- Said method 100 comprises a series of three successive steps:
- the most efficient is the method based on the QRAPMASTER sequence. Indeed, this method makes it possible to evaluate in a single acquisition all the relevant physiological parameters necessary for the generation of images, in particular of spin echo or inversion-recovery.
- this method has a number of disadvantages. Indeed, such a method related to this type of sequence is very specific and can not be applied to any magnetic resonance imaging analysis system. Consequently, such a method has very high implementation and maintenance costs for the establishment wishing to use it, of the order of hundreds of thousands of dollars. Moreover, such a method is applicable to only one type of organ, the brain. Other organs can not be analyzed.
- the step of estimating the parameters T1, T2 or PD generally consists of a step of calculation by linearization of the equations connecting the physiological parameters to the signal S. Admittedly, the calculation times are reduced. Nevertheless, said estimates are very sensitive to noise.
- These methods thus require the use of sequences for which the signal-to-noise ratio (also known under the English terminology "SNR for Signal to Noise Ratio") is high or risk of obtaining unreliable or erroneous weighted maps or images. .
- SNR Signal to Noise Ratio
- the increase in the acquisition time considerably increasing the duration of the examination, generates the same negative consequences, as previously explained, such as the discomfort for the patient, a low frequency of examination of patients due to exams of relatively long duration, and a high examination cost.
- the invention makes it possible to meet the great majority of the disadvantages raised by the known solutions.
- a method for producing an estimate of a physiological parameter of an elementary volume - called voxel - of an organ, said method being implemented by means of processing a processing unit. of a Magnetic Resonance imaging analysis system, and comprising a step for estimating said physiological parameter.
- a method comprises a step for estimating a first physiological parameter from first experimental signals resulting from a first acquisition of signals, as well as a step for estimating a second physiological parameter from second experimental signals resulting from a second acquisition of signals.
- said step for estimating the physiological parameter consists in producing the estimated value of said physiological parameter from the respective estimates of the first and second physiological parameters.
- a method according to the invention may further comprise steps for producing first and second experimental signals respectively from first and second signal acquisitions.
- the invention provides that the step for estimating a first physiological parameter may consist of a Bayesian estimation step of estimating said first physiological parameter by calculating its marginalized posterior distribution. .
- the invention provides that the step for estimating a second physiological parameter consists of a step of estimation by a Bayesian method, of estimating said second physiological parameter by the calculation of its marginalized posterior distribution.
- a method according to the invention may comprise a subsequent step for triggering a restitution of the estimated physiological parameter and / or the first, second physiological parameters.
- a method according to the invention may further comprise a preliminary step of pretreatment of the first and / or second experimental signals obtained respectively from the first and / or second acquisitions by Magnetic Resonance, said step being arranged to correct said first and / or second experimental signals .
- the invention relates to a method for producing an estimate of a physiological parameter of a region of interest, said region comprising at least one voxel.
- said physiological parameter is estimated by voxel by means of a method according to the first subject of the invention.
- a method according to the invention may further include a step subsequent to triggering the return of the estimate of the physiological parameter, the first and / or second physiological parameters for each voxel of the region of interest in the form of a map describing a physiological parameter.
- a method according to the invention may comprise in addition a subsequent step for generating a weighted image from the values produced of the estimated physiological parameter, first and second physiological parameters for a predetermined acquisition sequence.
- the invention relates to a processing unit comprising means for communicating with the outside world and processing means cooperating with storage means.
- the means for communicating are arranged to receive from the outside world first and second experimental signals from the first and / or second acquisitions of Magnetic Resonance signals and the storage means comprise instructions that can be executed or interpreted by the means of communication. processing whose interpretation or execution of said instructions by said processing means causes the implementation of a method according to the first subject of the invention.
- the means for communicating a treatment unit according to the invention can deliver an estimated physiological parameter in a format appropriate to restitution means able to restore it to an user.
- the invention relates to a Magnetic Resonance imaging analysis system comprising a processing unit according to the invention and restitution means able to restore to a user a physiological parameter according to a method according to the first object of the invention and implemented by said processing unit.
- the invention relates to a computer program product comprising one or more instructions interpretable or executable by the processing means of a processing unit according to the invention.
- Said processing unit further comprises storage means or cooperating with such storage means, said program being loadable in said storage means.
- Said instructions by said processing means are such that their interpretation or execution causes the implementation of a method according to the first subject of the invention.
- FIGS. 1 and 2 previously described, show two alternative embodiments of a medical imaging analysis system, possibly by Magnetic Resonance;
- FIGS. 3A, 3B and 3C previously described, show three examples of cards or weighted images obtained by a Nuclear Magnetic Resonance Imaging apparatus according to the state of the
- FIG. 4 previously described, presents a simplified flowchart of a method according to the state of the art
- FIG. 5 schematically depicts a simplified flowchart of a method according to the invention
- FIGS. 6A, 6B and 6C present three examples of physiological parameter maps, estimated according to a method according to the invention
- FIGS. 7A, 7B and 7C present three examples of weighted images generated and restored according to a method according to the invention.
- FIG. 5 schematizes a method 200 for estimating a physiological parameter of an elementary volume - called voxel - of an organ.
- voxel means any pixel with a thickness.
- a method 200 according to the invention advantageously comprises a step 230 for estimating said physiological parameter.
- a method 200 according to the invention also comprises a step 221 for estimating a first physiological parameter from first experimental signals resulting from a first acquisition of signals.
- the method also includes a step 222 for estimating a second physiological parameter from second experimental signals resulting from a second signal acquisition.
- first and second experimental signals can be directly downloaded from a server, advantageously arranged to store said first and second signals.
- the method 200 according to the invention may comprise steps 211, 212 for producing first and second experimental signals respectively from first and second signal acquisitions.
- a step 211 may advantageously consist in the implementation of a first acquisition of signals from a first acquisition sequence determined to estimate a first physiological parameter.
- a step 212 may advantageously consist in the implementation of a second acquisition of signals from a second acquisition sequence determined to estimate a second physiological parameter.
- the selection of the first and / or second acquisition sequences can be performed automatically or manually, during a preliminary step of configuring the implementation of a method 200 according to the invention, for example via the parameters 16 described. previously in connection with Figures 1 and 2.
- such first and second physiological parameters to be estimated may be respectively the physiological parameters T1 or T2.
- the sequences of the first and second signal acquisitions may advantageously consist of two respective sequences of type T1 and T2 mapping, or even T2 * mapping in the case of a sequence echo gradient type acquisition. Therefore, by way of nonlimiting example, said step 211 for producing first experimental signals from a first signal acquisition may comprise the use of a gradient echo type sequence with different flip-flop angles. to estimate the first physiological parameter T1.
- a sequence is particularly appreciated since it is very fast and available for any magnetic resonance imaging analysis system.
- step 221 for estimating the first physiological parameter may consist of a linearization calculation step of an equation connecting said first physiological parameter with the first experimental signals S.
- step 221 for estimating the first physiological parameter T1 can consist of an estimation step by a Bayesian method.
- said Bayesian method is described in the document WO2012049421 filed by the company OLEA MEDICAL or in the document WO2010139895A1 also filed by the company OLEA MEDICAL.
- Said Bayesian methods may consist in estimating the first physiological parameter by calculating its marginalized posterior distribution.
- Bayesian methods make it possible in particular to increase the accuracy of the estimates and to reduce the sensitivity to noise.
- step 212 for producing second experimental signals from a second signal acquisition may comprise the use of a multi-valued spin echo sequence. echo to estimate the second physiological parameter ⁇ 2.
- a second sequence is particularly preferred since it is very fast and available for any magnetic resonance imaging analysis system.
- step 222 for estimating the second physiological parameter may consist of a step of calculation by linearization of an equation connecting said second physiological parameter to the second experimental signals S.
- step 222 for estimating the second physiological parameter T2 can consist of a Bayesian estimation step. Said Bayesian methods may consist in estimating the second physiological parameter by calculating its marginalized posterior distribution.
- Such a Bayesian method makes it possible in particular to increase the precision of the estimates and to reduce the sensitivity to noise.
- the second acquisition sequence can then be shorter in time, obtaining qualitatively identical results to those obtained with a longer sequence.
- the spatial resolution of the images, and therefore the noise level can be significantly improved without degrading the estimates.
- step 230 for estimating the physiological parameter of a method 200 according to the invention consists in producing the estimated value of said physiological parameter from the respective estimates of the first and second physiological parameters.
- the physiological parameter to be estimated may be the physiological parameter PD.
- the Magnetic Resonance imaging analysis system may comprise means for restoring to a user 6, said reproduction means 5 cooperating advantageously with the processing unit 4.
- Such means of rendering make it possible to have an advantageously graphic, sound or other rendering and may comprise, for example, a screen or loudspeakers.
- a method 200 according to the invention may also comprise a subsequent step for triggering a restitution of the estimated physiological parameter and / or the first, second physiological parameter in an appropriate format.
- the first and second physiological parameters are respectively the physiological parameters T1 and T2
- the physiological parameter estimated is the physiological parameter PD
- such a reproduction can consist of a graphical representation in the form of maps of the physiological parameters.
- a method 200 according to the invention may also comprise one or more steps of pretreatment of the first and / or second experimental signals respectively obtained from the first and / or second acquisitions of signals by magnetic resonance, said step or steps of correcting said first and / or second experimental signals, in particular by the correction of artifacts or the application of any other hotfixes.
- steps may consist of steps of:
- a rigid or non-rigid registration algorithm may be chosen; - Co-registration or registration between the first and second acquisition sequences if the field of view of the latter is changed, or if the patient has moved between the first and second sequences.
- Such a co-registation can advantageously be translated in the form of a rigid or non-rigid co-registration algorithm;
- noise reduction step in the acquisitions of the two sequences.
- a noise reduction step may advantageously be embodied in the form of an algorithm for convolution smoothing of images with a Gaussian nucleus;
- the invention further relates to a method 300 for producing an estimate of a physiological parameter of a region of interest.
- "Region of interest” means any area with at least one voxel. Nevertheless, a region of interest can not be limited to a single voxel, but may comprise a plurality of voxels, advantageously selected manually or automatically.
- said physiological parameter can be estimated for each voxel by means of a method 200 according to the invention, as described above, in particular in connection with FIG. 5, said method being implemented iteratively for each voxel by the processing means of the processing unit 4.
- a method 300 in accordance with the invention may also comprise a subsequent step 350 for triggering a restitution of the estimated physiological parameter and / or the first, second parameter.
- the Magnetic Resonance imaging analysis system comprises means for rendering 5 to a user 6, said reproduction means 5 cooperating advantageously with the processing unit 4.
- the physiological parameter PD such a reproduction can consist in the display or printing of a graphical representation in the form of maps of the first, second physiological parameter T1 and T2 and / or the physiological parameter PD or one or several estimated values of the first, second physiological parameters T1 and T2 and / or the physiological parameter PD. Examples of such parameter cards will be described later in connection with FIGS. 6A, 6B and 6C.
- a method 300 for estimating a hemodynamic parameter of a region of interest may further comprise a subsequent step 340 for generating a weighted image from the values produced of the estimated physiological parameter, first and second physiological parameters for a predetermined acquisition sequence, when the Magnetic Resonance imaging analysis system comprises means for restoring said system, said restitution means co-operating with the processing unit 4.
- a step 340 allows in particular to obtain valuable information on the physiological parameters and to generate one or more weighted images from any type of acquisition sequence chosen, without requiring the implementation of a new examination, therefore a new acquisition, extremely expensive in time and money.
- the first, second physiological parameters and the estimated physiological parameter may be, advantageously and non-limitatively, the physiological parameters T1, T2 and PD, respectively.
- the method 300 may comprise a configuration step (not shown in FIG. 5), prior to step 340 for generating a weighted image, for selecting an acquisition sequence and the associated acquisition parameters, such as as examples not limiting parameters TR, TE and 77.
- Such a selection of sequences and parameters can be performed manually by a user or can be implemented automatically. Examples of such weighted images will be described later in connection with FIGS. 7A, 7B and 7C.
- a method 300 in accordance with the invention may advantageously comprise a step (not shown in FIG. 5) to selectively segment tissue from known theoretical values of said tissue. For example, suppose the values of T1 and T2 of the white matter are known and are respectively 560130 ms and 77 ⁇ 5 ms. A segmentation based on thresholds of the estimated values of T1 and T2 makes it possible to extract the white matter voxels according to the following equation:
- weighted images for segmentation purposes would be to use the estimated values of T1, T2 and PD as input data of a partitioning algorithm such as the k-algorithm. means ("k-means" according to an English terminology).
- step 212 may advantageously consist in the implementation of a second acquisition of signals from a second acquisition sequence determined to estimate a second physiological parameter.
- a second acquisition sequence may advantageously be a T2 mapping sequence.
- said T2 mapping sequence implemented by processing means of a processing unit 4 of a Magnetic Resonance imaging analysis system, can advantageously be an echo sequence.
- multi-echo spin When using such a multi echo spin echo sequence, the experimental signal in each voxel can be calculated by means of a decreasing exponential function, such that:
- step 222 for estimating the second physiological parameter T2 consists of a sub-step of linearization computation of the preceding equation by taking the logarithm of the experimental signal combined with a linear regression substep. Nevertheless, the use of such sub-steps is not satisfactory, since said substeps have large computational uncertainties.
- step 222 for estimating the second physiological parameter T2 may consist of a Bayesian estimation, such as, as specified above, that described in document WO2012049421 or that described in document WO2010139895A1.
- a template is predefined manually or automatically.
- the Bayes theorem can then be applied, producing an equation linking the posterior distribution of the parameters of said predefined model to the prior distributions of these same parameters P (T 2 ), P (S 0 ), P (o) and to the likelihood function , the likelihood function being defined as the probability distribution of the data knowing the parameters, such that:
- ⁇ J is the standard deviation of the noise affecting the data D in a voxel of interest.
- the data D correspond to the second experimental signals obtained by the acquisition of a second sequence.
- the estimation of any parameter of interest is made using the marginalized posterior distribution estimate of said parameter of interest.
- the estimation of the marginalized posterior distribution of the second physiological parameter T2 can be calculated for a voxel of interest by the evaluation of the relation:
- an estimate of the second physiological parameter T2 can finally be calculated as the maximum a posteriori
- the method includes sub-steps for calculating, estimating and / or selecting the prior distributions of these same parameters P (T 2 ),
- N is the number of echo times used to make the acquisition.
- Prior distributions of the parameters can be chosen, manually or automatically, also during a preliminary step of configuring the implementation of a method 200 according to the invention, for example via the parameters 16. described above with reference to FIGS. 1 and 2, as non-limiting examples:
- the marginalized posterior distribution of the second physiological parameter T2 for a given voxel can then be produced, such that:
- the parameters S 0 and ⁇ J at the voxel of interest can be produced analytically, such as:
- FIG. 6A shows an estimation map of the second physiological parameter T2, resulting from an iterative implementation of the method 200 for a plurality of voxels. Let's describe in a second step steps 211 and
- step 211 may advantageously consist in the implementation of a first acquisition of signals from a first acquisition sequence determined to estimate a first physiological parameter.
- a first acquisition sequence can advantageously be a sequence of Tl mapping.
- said Tl mapping sequence implemented by processing means of a processing unit 4 of a Magnetic Resonance imaging analysis system, may advantageously be a sequence of FIG. inversion-recovery, a look-locker sequence or a sequence with variable flip-flops.
- variable flip-flops also known by the abbreviation "VFA for Variable Flip Angle”
- VFA Variable Flip Angle
- This sequence is indeed the fastest of the sequences compared to previous ones.
- the experimental signal in each voxel of interest can be expressed as a proportionality relation, such that:
- a conventional approach for estimating the parameters T1 and 0 is to perform a linearization calculation of the preceding proportionality relation.
- step 221 for estimating the first physiological parameter T1 may advantageously consist of a Bayesian estimation, such as, as previously stated, that described in document WO2012049421 or that described in document WO2010139895A1.
- a template is predefined manually or automatically.
- the Bayes theorem can then be applied producing an equation linking the posterior distribution of the parameters PT ⁇ M Q ⁇ D) of the predefined model to the prior distributions of these same parameters ⁇ (), P (M 0 ), P (o) and the likelihood function , the likelihood function being defined as the probability distribution of the data knowing the parameters, such that:
- the data D correspond to the first experimental signals obtained by the acquisition of said first sequence.
- the estimation of any parameter of interest is made using the marginalized posterior distribution estimation of said parameter of interest.
- the estimation of the marginalized posterior distribution of the first physiological parameter T1 can be calculated advantageously for a voxel of interest such as:
- an estimate of the first physiological parameter T1 can be calculated for example in the form of the maximum a posteriori
- the method advantageously comprises sub-steps for calculating, estimating and / or selecting the prior distributions of these same parameters ⁇ ( ⁇ ⁇ ,
- N is the number of flip-flops used to acquire.
- the prior distributions of the parameters of the model, for their part, can be chosen, manually or automatically, also during a preliminary step of configuring the implementation of a method 200 according to the invention, for example via the parameters 16 previously described in connection with FIGS. 1 and 2, as non-limiting examples:
- the marginalized posterior distribution the first physiological parameter T1 for a given voxel can then be produced, such that:
- parameters 0 and O at the voxel of interest can be computed anal tically, such as:
- FIG. 6B shows an estimation map of the first physiological parameter T1, resulting from an iterative implementation of the method 200 for a plurality of voxels.
- step 230 for estimating the physiological parameter PD.
- the parameter S 0 has been estimated thanks to step 222 to estimate the second physiological parameter T2 and to the sequence of T2 mapping.
- a parameter S 0 depends on the parameter physiological PD, and is weighted 7 ⁇ .
- the parameter can therefore be calculated according to the following proportionality relation, such that:
- the physiological parameter PD can be estimated at 230 to the voxel of interest as:
- acquisition parameter TR corresponds to the repetition time of the T2 mapping sequence.
- Said calculations, used in step 230 to estimate the physiological parameter PD are advantageously implemented by the processing means of a processing unit 4 of a magnetic resonance analysis system according to the invention.
- Such an estimate of the physiological parameter PD is relative and proportional to the actual value of the estimated physiological parameter PD.
- the proportionality factor between the estimate and the true value of PD depends solely on the properties of the Magnetic Resonance Imaging device.
- indetermination, in the form of a relative value does not pose a problem to the weighted image generation or synthetic MRI card.
- FIG. 6C shows an estimation map of the physiological parameter PD, resulting from an iterative implementation of the method 200 for a plurality of voxels.
- a user of said apparatus can instruct the apparatus to generate images, maps, or image sequences weighted in T1, T2, PD, or occult and / or hide certain types of tissue from a chosen acquisition sequence, as well as associated acquisition parameters.
- the invention provides that such card images or Tl, T2, PD weighted image sequences can be generated automatically and reproduced by a magnetic resonance imaging analysis system advantageously comprising a processing unit 4 and restitution means 5 cooperating advantageously with said processing unit 4.
- the acquisition sequence and associated acquisition parameters can thus be automatically chosen.
- the method for estimating and generating may comprise a step for calculating in each vox 1 i:
- the Magnetic Resonance imaging analysis system can generate and restore, using its processing unit and its restitution means, weighted synthetic images T1 or T2.
- the method for estimating and generating may comprise a step for calculating in each voxel i:
- the method may then comprise a step for generating an inversion-recovery image.
- This type of sequence makes it possible to suppress certain types of tissue, for example liquids.
- the method for estimating and generating may include a step for calculating in each voxel i:
- FIGS. 7A, 7B and 7C show three examples of weighted images generated according to a method according to the invention. Respectively, said FIGS. 7A, 7B and 7C respectively show T2-weighted synthetic images, weighted in T1 and in inversion-recovery, said inversion-recovery image highlighting the suppression of water.
- Figure 7A shows a T2-weighted image from a selected spin echo sequence with an echo time set to one hundred and twenty milliseconds and a repetition time set to one thousand five hundred milliseconds.
- 7B shows a T1-weighted image from a selected spin echo sequence with an echo time set to thirty milliseconds and a repetition time set to five thousand milliseconds.
- i has 7C presents a weighted image from a sequence of recovery inversion- selected with an echo time set to fifty milliseconds, a repetition time set to twenty thousand milliseconds and an inversion time defined one thousand seven hundred milliseconds.
- the invention makes it possible to make available to a practitioner a whole set of relevant and coherent information, information that is rapidly available thanks to the use of a method according to the invention.
- This provision is made possible by an adaptation of the processing unit 4 according to FIG. 1 or 2, in that the processing means implement a method for estimating a physiological parameter of a voxel or a region of interest comprising producing the estimated value of said physiological parameter from the respective estimates of the first and second physiological parameters.
- Such an implementation is advantageously made possible by the loading or the recording, within memory means cooperating with said processing means, of a computer program product. The latter has indeed interpretable instructions and / or executable by said processing means.
- the interpretation or execution of said instructions triggers the implementation of a method 200 or 300 according to the invention.
- the means for communicating with the outside world of said processing unit can deliver a physiological parameter, namely the estimated parameters 14, in a format appropriate to restitution means able to restore it to a user 6, said estimated physiological parameter being advantageously restored in the form, for example, of cards or weighted images such as those illustrated in Figures 6A to 6C and 7A to 7C. Thanks to the invention, the information delivered is more numerous, consistent, reproducible and accurate. The information available to the practitioner is thus likely to increase the confidence and speed of the practitioner in determining a diagnosis and making decisions.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Signal Processing (AREA)
- Heart & Thoracic Surgery (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
Claims
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/560,647 US20180095152A1 (en) | 2015-03-23 | 2016-03-22 | System and method for estimating a physiological parameter of an elementary volume |
EP16723794.0A EP3274734A1 (fr) | 2015-03-23 | 2016-03-22 | Système et procédé irm pour estimer un paramètre physiologique à partir de deux autres paramètres physiologiques estimés |
AU2016238675A AU2016238675A1 (en) | 2015-03-23 | 2016-03-22 | MRI system and method for estimating a physiological parameter on the basis of two other estimated pysiological parameters |
JP2017549700A JP2018509240A (ja) | 2015-03-23 | 2016-03-22 | 基本ボリュームの生理的パラメータを推定するシステムおよび方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562136931P | 2015-03-23 | 2015-03-23 | |
US62/136,931 | 2015-03-23 | ||
FR1555656 | 2015-06-19 | ||
FR1555656A FR3037496A1 (fr) | 2015-06-19 | 2015-06-19 | Systeme et procede pour estimer un parametre physiologique d'un volume elementaire |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2016151232A1 true WO2016151232A1 (fr) | 2016-09-29 |
Family
ID=55345887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/FR2016/050626 WO2016151232A1 (fr) | 2015-03-23 | 2016-03-22 | Système et procédé irm pour estimer un paramètre physiologique à partir de deux autres paramètres physiologiques estimés |
Country Status (6)
Country | Link |
---|---|
US (1) | US20180095152A1 (fr) |
EP (1) | EP3274734A1 (fr) |
JP (1) | JP2018509240A (fr) |
AU (1) | AU2016238675A1 (fr) |
FR (1) | FR3037496A1 (fr) |
WO (1) | WO2016151232A1 (fr) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3060803A1 (fr) * | 2016-12-20 | 2018-06-22 | Ondilo | Procede pour maintenir un equilibre d'un parametre physico-chimique d'un milieu, produit programme d'ordinateur et moyen electronique associes |
JP2019005557A (ja) * | 2017-06-22 | 2019-01-17 | キヤノンメディカルシステムズ株式会社 | 画像処理装置、磁気共鳴イメージング装置及び画像処理プログラム |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6979151B2 (ja) * | 2017-05-26 | 2021-12-08 | 富士フイルムヘルスケア株式会社 | 磁気共鳴イメージング装置及び磁気共鳴画像処理方法 |
CN109938733B (zh) * | 2019-05-14 | 2020-06-30 | 浙江大学 | 一种高质量磁共振图像合成方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010139895A1 (fr) | 2009-06-05 | 2010-12-09 | Olea Medical | Procédé pour estimer des paramètres hémodynamiques par estimation conjointe des paramètres d'un modèle global de perfusion |
WO2012049421A2 (fr) | 2010-10-11 | 2012-04-19 | Olea Medical | Système et procédé pour estimer une quantité d'intérêt d'un système dynamique artère/tissu/veine |
US20140330106A1 (en) * | 2011-12-13 | 2014-11-06 | Isis Innovation Limited | Multi-parametric magnetic resonance diagnosis & staging of liver disease |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6012931B2 (ja) * | 2010-03-30 | 2016-10-25 | 東芝メディカルシステムズ株式会社 | 画像処理装置及び画像処理装置の制御プログラム |
FR2996667B1 (fr) * | 2012-10-05 | 2015-12-11 | Olea Medical | Systeme et procede pour estimer une quantite d'interet dans un systeme cinematique par tomographie par agent de contraste |
US20140213912A1 (en) * | 2013-01-29 | 2014-07-31 | Covidien Lp | Low power monitoring systems and method |
US20140275886A1 (en) * | 2013-03-14 | 2014-09-18 | Streamline Automation, Llc | Sensor fusion and probabilistic parameter estimation method and apparatus |
-
2015
- 2015-06-19 FR FR1555656A patent/FR3037496A1/fr active Pending
-
2016
- 2016-03-22 AU AU2016238675A patent/AU2016238675A1/en not_active Abandoned
- 2016-03-22 EP EP16723794.0A patent/EP3274734A1/fr not_active Withdrawn
- 2016-03-22 US US15/560,647 patent/US20180095152A1/en not_active Abandoned
- 2016-03-22 JP JP2017549700A patent/JP2018509240A/ja active Pending
- 2016-03-22 WO PCT/FR2016/050626 patent/WO2016151232A1/fr active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010139895A1 (fr) | 2009-06-05 | 2010-12-09 | Olea Medical | Procédé pour estimer des paramètres hémodynamiques par estimation conjointe des paramètres d'un modèle global de perfusion |
WO2012049421A2 (fr) | 2010-10-11 | 2012-04-19 | Olea Medical | Système et procédé pour estimer une quantité d'intérêt d'un système dynamique artère/tissu/veine |
US20140330106A1 (en) * | 2011-12-13 | 2014-11-06 | Isis Innovation Limited | Multi-parametric magnetic resonance diagnosis & staging of liver disease |
Non-Patent Citations (6)
Title |
---|
FRANCESCO DE PASQUALE ET AL: "Bayesian estimation of relaxation times T1 in MR images of irradiated Fricke-agarose gels", MAGNETIC RESONANCE IMAGING, vol. 18, no. 6, 1 July 2000 (2000-07-01), TARRYTOWN, NY, US, pages 721 - 731, XP055272491, ISSN: 0730-725X, DOI: 10.1016/S0730-725X(00)00149-1 * |
IVES LEVESQUE ET AL: "The role of edema and demyelination in chronic T1 black holes: A quantitative magnetization transfer study", JOURNAL OF MAGNETIC RESONANCE IMAGING, vol. 21, no. 2, 21 January 2005 (2005-01-21), US, pages 103 - 110, XP055272444, ISSN: 1053-1807, DOI: 10.1002/jmri.20231 * |
KELVIN LAYTON ET AL: "Improved quantification of MRI relaxation rates using Bayesian estimation", ACOUSTICS SPEECH AND SIGNAL PROCESSING (ICASSP), 2010 IEEE INTERNATIONAL CONFERENCE ON, IEEE, PISCATAWAY, NJ, USA, 14 March 2010 (2010-03-14), pages 481 - 484, XP031697621, ISBN: 978-1-4244-4295-9 * |
NEEB H ET AL: "A new method for fast quantitative mapping of absolute water content in vivo", NEUROIMAGE, ACADEMIC PRESS, ORLANDO, FL, US, vol. 31, no. 3, 1 July 2006 (2006-07-01), pages 1156 - 1168, XP024906405, ISSN: 1053-8119, [retrieved on 20060701], DOI: 10.1016/J.NEUROIMAGE.2005.12.063 * |
NEEB H ET AL: "Fast quantitative mapping of absolute water content with full brain coverage", NEUROIMAGE, ACADEMIC PRESS, ORLANDO, FL, US, vol. 42, no. 3, 11 April 2008 (2008-04-11), pages 1094 - 1109, XP026446471, ISSN: 1053-8119, [retrieved on 20080411] * |
WARNTJES J B M ET AL: "Novel method for rapid, simultaneous T1, T2* and proton density quantification", MAGNETIC RESONANCE IN MEDICINE, JOHN WILEY & SONS, INC, US, vol. 57, 26 February 2007 (2007-02-26), pages 528 - 537, XP002639788, ISSN: 0740-3194, [retrieved on 20070226], DOI: 10.1002/MRM.21165 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3060803A1 (fr) * | 2016-12-20 | 2018-06-22 | Ondilo | Procede pour maintenir un equilibre d'un parametre physico-chimique d'un milieu, produit programme d'ordinateur et moyen electronique associes |
WO2018115751A1 (fr) * | 2016-12-20 | 2018-06-28 | Ondilo | Procede pour maintenir un equilibre d'un parametre physico-chimique d'un milieu, produit programme d'ordinateur et module electronique associes |
US11543839B2 (en) | 2016-12-20 | 2023-01-03 | Ondilo | Method for maintaining an equilibrium of a physico-chemical parameter of a medium, associated computer program product and electronic module |
JP2019005557A (ja) * | 2017-06-22 | 2019-01-17 | キヤノンメディカルシステムズ株式会社 | 画像処理装置、磁気共鳴イメージング装置及び画像処理プログラム |
JP7246864B2 (ja) | 2017-06-22 | 2023-03-28 | キヤノンメディカルシステムズ株式会社 | 画像処理装置、磁気共鳴イメージング装置及び画像処理プログラム |
Also Published As
Publication number | Publication date |
---|---|
JP2018509240A (ja) | 2018-04-05 |
FR3037496A1 (fr) | 2016-12-23 |
US20180095152A1 (en) | 2018-04-05 |
AU2016238675A1 (en) | 2017-10-26 |
EP3274734A1 (fr) | 2018-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2628018B1 (fr) | Système et procédé pour estimer une quantité d'intérêt d'un système dynamique artère/tissu/veine | |
EP2437664B1 (fr) | Procédé pour estimer des paramètres hémodynamiques par estimation conjointe des paramètres d'un modèle global de perfusion | |
EP3274734A1 (fr) | Système et procédé irm pour estimer un paramètre physiologique à partir de deux autres paramètres physiologiques estimés | |
EP2904582A1 (fr) | Système et procédé pour estimer une quantité d'intérêt dans un système cinématique par tomographie par agent de contraste | |
EP2435985B1 (fr) | Procede de quantification de l'evolution de pathologies impliquant des changements de volumes de corps, notamment de tumeurs | |
EP2924586A1 (fr) | Système et procédé pour estimer une quantité d'intérêt d'un système dynamique artère/tissu/veine | |
EP3207397B1 (fr) | Système et procédé pour estimer une quantité d'intérêt d'un système dynamique artère/tissu/veine | |
EP4147062B1 (fr) | Procédé de post-traitement d'un échantillonnage d'une séquence d'acquisition wasab1 | |
WO2023275478A1 (fr) | Procédé d'atténuation du bruit d'images résultant d'acquisitions multiples par résonance magnétique | |
WO2011064512A1 (fr) | Procédé pour estimer l'adéquation d'un processus d'estimation de paramètres hémodynamiques de perfusion | |
EP4128164A1 (fr) | Procédé de reconstruction par tomographie rayon x et dispositif associé | |
EP2345007A1 (fr) | Procédé d'estimation de la concentration d'un traceur dans un ensemble de structures tissulaires, support d'information et dispositif correspondants | |
EP3509483B1 (fr) | Systeme et procede pour reconstruire un signal physiologique d'un systeme dynamique artere/tissu/veine d'un organe dans un espace surfacique | |
FR3137827A1 (fr) | Procédé de post-traitement d’une séquence d’acquisition de perfusion par un appareil d’imagerie médicale | |
WO2021234278A1 (fr) | Système et procédé pour estimer un indicateur d'activité tissulaire d'un organe |
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: 16723794 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2017549700 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 15560647 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REEP | Request for entry into the european phase |
Ref document number: 2016723794 Country of ref document: EP |
|
ENP | Entry into the national phase |
Ref document number: 2016238675 Country of ref document: AU Date of ref document: 20160322 Kind code of ref document: A |