CN111354056B - Method for accelerating diffusion magnetic resonance imaging acquisition - Google Patents
Method for accelerating diffusion magnetic resonance imaging acquisition Download PDFInfo
- Publication number
- CN111354056B CN111354056B CN202010445650.2A CN202010445650A CN111354056B CN 111354056 B CN111354056 B CN 111354056B CN 202010445650 A CN202010445650 A CN 202010445650A CN 111354056 B CN111354056 B CN 111354056B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- diffusion
- machine learning
- diffusion magnetic
- resonance image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a method for accelerating diffusion magnetic resonance imaging acquisition. The method comprises the following steps: 1) sorting a data set prepared for machine learning, wherein the data set comprises diffusion magnetic resonance images and corresponding diffusion parameter maps; 2) constructing and training a machine learning model to enable the machine learning model to learn the corresponding relation between the diffusion magnetic resonance image and the diffusion parameter map; 3) calculating the importance of each diffusion magnetic resonance image in the training set to the trained machine learning model; 4) preserving the diffusion magnetic resonance image with the importance meeting the set condition as a new data set; 5) retraining the machine learning model by using the new data set to enable the machine learning model to learn the corresponding relation between the diffusion magnetic resonance image and the diffusion parameter map; 6) calculating a diffusion parameter map by using the machine learning model trained in the step 5), verifying the calculation effect of the diffusion parameter map, if the verification is passed, keeping the corresponding scanning condition of the diffusion magnetic resonance image used in the step 5), and then collecting the diffusion magnetic resonance image by using the scanning condition.
Description
Technical Field
The invention belongs to the technical field of medical image and Magnetic Resonance Image (MRI) image processing, and particularly relates to a diffusion-accelerated MRI acquisition method based on a machine learning feature selection technology.
Background
Diffusion magnetic resonance imaging (Diffusion MRI) in Magnetic Resonance Imaging (MRI) is also commonly referred to clinically as Diffusion magnetic resonance imaging. Diffusion magnetic resonance imaging indirectly observes microstructural properties of brain tissue, such as white matter nerve fiber tract orientation, by measuring the degree and direction of water molecule diffusion within the tissue. It is the only technology that can show the trend of the white matter nerve fiber bundle of the living body without damage at present, therefore is widely used in clinical and scientific research. Diffusion magnetic resonance imaging generally requires acquiring images of a plurality of scanning conditions (for example, setting different diffusion sensitive gradient strengths and directions), and sufficient images acquired under different conditions can be brought into model calculation to obtain a diffusion parameter map reflecting microstructure information.
There are currently many models for calculating Diffusion magnetic resonance Imaging, including Diffusion Tensor Imaging (DTI), Diffusion Kurtosis Imaging (DKI), Neurite direction Dispersion and Density Imaging (NODDI), Constrained Spherical Deconvolution (CSD), Fractional Motion (FM), etc. In order to calculate the diffusion parameters more stably and accurately, these models require the acquisition of images of many different scan conditions. This results in an excessively long image acquisition time, which severely limits the practical application of the model.
In recent years, there have been many researchers attempting to calculate diffusion parameters using machine learning techniques. The result shows that compared with the prior model fitting method, the calculation method of machine learning has higher precision and higher speed. Especially, the requirement on the original image is lower, and the result which can be calculated by using a large number of images by using the model fitting method can be obtained only by a small number of images under different conditions. Therefore, the total number of images required to be acquired can be reduced and the image acquisition process can be accelerated by estimating the diffusion parameter map by using a machine learning method.
Although there have been studies showing that machine learning methods can reduce the number of images that need to be acquired, they all choose randomly which images to reject and which images to retain. This is likely to eliminate images containing critical information, resulting in loss of computational accuracy. At present, a method which can reduce the images required to be acquired and can ensure the calculation accuracy to the maximum extent does not exist.
Disclosure of Invention
In view of the above, the main objective of the present invention is to provide a method for accelerating diffusion mri acquisition based on a machine learning feature selection technique. After the machine model is trained, the importance of each input image to the machine learning model is calculated. Based on the method, the images with low importance on the model are removed during acquisition, the number of the images required to be acquired is reduced, and the image acquisition process is accelerated.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a method for accelerating diffusion magnetic resonance imaging acquisition based on a machine learning feature selection technology comprises the following steps:
A. sorting a data set prepared for machine learning, wherein the data set needs to comprise a diffusion magnetic resonance image and a corresponding diffusion parameter map thereof;
B. constructing and training a machine learning model to enable the machine learning model to learn the corresponding relation between the diffusion magnetic resonance image and the diffusion parameter map;
C. calculating the importance of each diffusion magnetic resonance image in the data set to the trained machine learning model;
D. sorting the importance calculated in the step C, excluding images with low importance, and keeping diffusion magnetic resonance images with the importance meeting set conditions (for example, if the scanning time needs to be changed to half of the original scanning time, namely the aim is to reduce the image acquisition time by 50%, the images with the lowest importance need to be excluded, and for example, if the final calculation precision deviation is required to be accelerated within 20%, the three steps of the step C, D, E need to be iterated, and the least importance is eliminated at the end of each time until the precision requirement is exceeded at the end of the elimination) as a new data set;
E. retraining the machine learning model by using a new data set to enable the machine learning model to learn the corresponding relation between the diffusion magnetic resonance images and the diffusion parameter map;
F. and E, verifying whether the calculation precision meets the preset requirement or not, using the machine learning model trained in the step E after verification for calculating a diffusion parameter map later, and using the scanning condition corresponding to the image still reserved at the moment as a scheme for acquiring the diffusion magnetic resonance image later.
Wherein: the preparation of the data set in step a aims to obtain a set of diffusion magnetic resonance images and corresponding diffusion parameter maps thereof, and can be divided into the following two methods:
a1, establishing a data set by using the simulation data, specifically: presetting a diffusion parameter map in a reasonable range, designing scanning conditions, calculating corresponding diffusion magnetic resonance signals from the diffusion parameters and the scanning conditions according to a model, synthesizing a diffusion magnetic resonance image, and adding noise interference of a reasonable degree to the synthesized diffusion magnetic resonance image;
a2, establishing a data set by using the real data, specifically: and collecting the originally collected diffusion magnetic resonance image, and calculating a diffusion parameter map by using a model fitting method.
And B, constructing and training the machine learning model, namely taking the diffusion magnetic resonance image as the input of the machine learning model, and taking the diffusion parameter map as a target.
The calculation of the importance of each input image in step C may be implemented in various ways. For some machine learning models with simple structures, weighting coefficients after model training can be used for measurement; for machine learning models such as decision trees, the division standard difference value before and after node division can be used for calculation; for more complex machine learning models, the calculation can be done by using random ordering tests.
The method for accelerating diffusion magnetic resonance imaging acquisition based on the machine learning feature selection technology has the following beneficial effects:
the method is suitable for various diffusion magnetic resonance models, can be directly used for calculating parameters, and can be used for follow-up white matter fiber tract tracking and the like. Many diffusion magnetic resonance models need to acquire many images to perform fitting calculation, so that the image acquisition time is too long, and the diffusion magnetic resonance models are difficult to apply in actual situations. By adopting the method, the number of images required to be acquired can be reduced under the condition of ensuring the calculation accuracy of the diffusion parameter map, so that a large amount of image acquisition time is saved, and the popularization and application of each diffusion magnetic resonance model are facilitated.
Drawings
Fig. 1 is a schematic flow chart of a method for accelerating diffusion mri acquisition based on a machine learning feature selection technique according to the present invention.
Fig. 2 is a comparison graph of the importance of 18 diffuse magnetic resonance images in the absence of noise to the trained model.
Figure 3 is a graph of the significance of 18 diffusion magnetic resonance images for a signal-to-noise ratio of 50 for the trained model.
Figure 4 is a graph of the significance of 18 diffusion magnetic resonance images for a signal-to-noise ratio of 40 for the trained model.
FIG. 5 is a graph of the significance of 18 diffusion MR images for a signal-to-noise ratio of 30 for the trained model.
FIG. 6 is a graph of the significance of 18 diffusion MR images for a signal-to-noise ratio of 20 for the trained model.
Fig. 7 shows the accuracy of calculating the noria index by the machine learning model under different signal-to-noise ratios after images with low importance are sequentially excluded.
Fig. 8 shows the accuracy of the machine learning model in computing the hurst exponent under different snr conditions after sequentially excluding the less important images.
Detailed Description
The present invention will be described in further detail below with reference to the accompanying drawings and embodiments thereof. It should be noted that the present invention is not limited to the specific diffusion models, image acquisition parameters, calculation methods, etc. described, and other embodiments, or combinations of other embodiments, are possible. Some of the processing steps in the present invention may be provided in plural.
The present example uses a fractional motion model as the diffusion magnetic resonance parameter model to be calculated, which includes two diffusion parameters: noah indexAnd Hers' specific number. When acquiring images using a Pulse Gradient Spin Echo (PGSE) sequence, these two parameters are related to the diffusion magnetic resonance acquisition parameters and the diffusion magnetic resonance signals by:
whereinS 0 Is the signal acquired without diffusion sensitive gradients;Sis the signal acquired in the presence of a diffusion sensitive gradient;is the generalized diffusion coefficient;is the gyromagnetic ratio constant;Gis the pulse gradient strength of the PGSE sequence;is the time interval of two pulse gradients of the PGSE sequence;is a dimensionless parameter:
Parameters of PGSE sequencesG、、Are all parameters of image acquisition. It can be seen from formula (1) that a plurality of groups of PGSE parameters need to be set during diffusion magnetic resonance image acquisition to obtain corresponding magnetic resonance signals, and then the diffusion parameters can be calculated in a fitting mannerAnd. In a previous study in which a fractional motion model was used for tumor grading, 18 sets of PGSE parameters were set at image acquisitionEach group contains different pulse gradient strength, two pulse gradient time intervals and pulse gradient duration. Plus the image without diffusion sensitive gradients, the total time for image acquisition is about 8 minutes 42 seconds.
Fig. 1 is a schematic flow chart of a method for accelerating diffusion mri acquisition based on a machine learning feature selection technique according to the present invention. Reducing the image acquisition time of this example in this way would include the steps of:
And 2, constructing and training a machine learning model. This example uses a Random Forest (RF) regression model to learn 18 signal attenuation and diffusion parametersAndthe corresponding relationship of (1). And respectively training a random forest regression model for each data set of the signal-to-noise ratio grade. Each random forest contains 200 decision trees, and the depth of each decision tree is automatically determined in the training stage. The trained loss function selects the Mean Squared Error function (MSE).
And 3, calculating the importance of each diffusion magnetic resonance image to the trained model. In this example, there are 18 signal attenuations that are characteristic of the input random forest model, and this step will calculate the importance of each signal attenuation to the model result. The random forest model may calculate the importance of a feature from the sum of the reduction in the degree of invisibility (Gini Impurity) of the feature split nodes.
And 4, selecting the images with high importance. The results of this example are shown in fig. 2-6, where the importance of 18 signal attenuations to random forests is given for each data set of signal-to-noise ratio levels, and the values with importance exceeding the average value are shown by the diagonal lines in fig. 2-6. Based on this, the low importance can be excluded from the 18 signal attenuations.
And 5, retraining the machine learning model by using the retained data. In this example, each time an input signal of lowest current importance is excluded, a random forest model is retrained with the remaining data. The random forest model may adopt the same setup as step 2, i.e. comprising 200 decision trees, the depth of each decision tree is automatically determined in the training stage, and the mean square error function is selected as the loss function.
And 6, verifying the calculation effect and confirming the scanning condition of the reserved image. Fig. 7 and 8 show the calculation effect of the machine learning model after images with low importance are gradually excluded. The accuracy of the results is obtained in this example using a weighted machine learning model with a quantitative determination Coefficient (coeffient determination). The closer the decision coefficient is to 1, the more accurate the result. As can be seen from fig. 7 and 8, even if only 6 signal attenuations are used for the calculation of the machine learning model, the accuracy is not greatly degraded. Therefore, only the acquisition conditions corresponding to the 6 input images which are most important to the machine learning model can be acquired during the subsequent actual image acquisition, so that the image acquisition time can be reduced to about one third of the original time.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention.
Claims (6)
1. A method of accelerating a diffusion magnetic resonance imaging acquisition, comprising the steps of:
1) sorting a data set prepared for machine learning, wherein the data set comprises diffusion magnetic resonance images and corresponding diffusion parameter maps;
2) constructing and training a machine learning model to enable the machine learning model to learn the corresponding relation between the diffusion magnetic resonance image and the diffusion parameter map;
3) calculating the importance of each diffusion magnetic resonance image in the training set to the trained machine learning model;
4) sorting according to importance, and keeping the diffusion magnetic resonance images with the importance meeting set conditions as a new data set;
5) retraining the machine learning model by using the new data set to enable the machine learning model to learn the corresponding relation between the diffusion magnetic resonance image and the diffusion parameter map;
6) calculating a diffusion parameter map by using the machine learning model trained in the step 5), verifying the calculation effect of the diffusion parameter map, if the verification is passed, keeping the corresponding scanning condition of the diffusion magnetic resonance image used in the step 5), and then collecting the diffusion magnetic resonance image by using the scanning condition.
2. The method of claim 1, wherein in step 1), the method of collating data sets prepared for machine learning is: firstly, setting a diffusion parameter graph and scanning conditions; and then calculating corresponding diffusion magnetic resonance signals according to the diffusion parameter map and the scanning conditions, synthesizing a diffusion magnetic resonance image, and adding noise interference to the synthesized diffusion magnetic resonance image.
3. The method of claim 1, wherein in step 1), the method of collating data sets prepared for machine learning is: and collecting the originally collected diffusion magnetic resonance image, and calculating a diffusion parameter image corresponding to the diffusion magnetic resonance image by using a model fitting method.
4. A method as claimed in claim 1, characterized in that the diffusion magnetic resonance image of which the importance exceeds the mean value of the importance is retained as a new data set.
5. The method of claim 1, wherein the diffusion magnetic resonance images with significance exceeding a set threshold are retained as a new data set or the K diffusion magnetic resonance images with significance ranking top are retained as a new data set.
6. The method of claim 1, wherein the machine learning model is a random forest model; each diffusion magnetic resonance image comprises a plurality of signal features; and 5) in the step of retraining the machine learning model by using new data, eliminating a signal feature with the lowest current importance every time, and retraining the random forest model by using the rest data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010445650.2A CN111354056B (en) | 2020-05-25 | 2020-05-25 | Method for accelerating diffusion magnetic resonance imaging acquisition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010445650.2A CN111354056B (en) | 2020-05-25 | 2020-05-25 | Method for accelerating diffusion magnetic resonance imaging acquisition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111354056A CN111354056A (en) | 2020-06-30 |
CN111354056B true CN111354056B (en) | 2020-10-27 |
Family
ID=71197729
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010445650.2A Active CN111354056B (en) | 2020-05-25 | 2020-05-25 | Method for accelerating diffusion magnetic resonance imaging acquisition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111354056B (en) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014162300A1 (en) * | 2013-04-05 | 2014-10-09 | Isis Innovation Ltd. | Acceleration of low-rank mri data acquisition |
CN104424386A (en) * | 2013-08-23 | 2015-03-18 | 北京大学 | Multi-parameter magnetic resonance image based prostate cancer computer auxiliary identification system |
CN105395198B (en) * | 2015-06-23 | 2018-10-26 | 高家红 | A kind of method and its application obtaining completely new diffusion magnetic resonance image contrast |
-
2020
- 2020-05-25 CN CN202010445650.2A patent/CN111354056B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111354056A (en) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Karayumak et al. | Retrospective harmonization of multi-site diffusion MRI data acquired with different acquisition parameters | |
CN104574298B (en) | A kind of noise-reduction method of more b values diffusion weightings images based on mutual information | |
US20200300957A1 (en) | Method for Estimating Systematic Imperfections in Medical Imaging Systems with Deep Learning | |
CN112614126A (en) | Magnetic resonance image brain region dividing method, system and device based on machine learning | |
CN106548213A (en) | Blood vessel recognition methodss and device | |
CN110969614A (en) | Brain age prediction method and system based on three-dimensional convolutional neural network | |
CN111445546A (en) | Image reconstruction method and device, electronic equipment and storage medium | |
CN111353935A (en) | Magnetic resonance imaging optimization method and device based on deep learning | |
CN110458842B (en) | Brain tumor segmentation method based on two-channel three-dimensional dense connection network | |
CN111583356B (en) | Magnetic resonance image synthesis method and device based on convolutional neural network | |
CN109799472B (en) | Magnetic resonance eddy current compensation method based on deep learning | |
CN109544652B (en) | Nuclear magnetic resonance multi-weighted imaging method based on depth generation antagonistic neural network | |
CN111354056B (en) | Method for accelerating diffusion magnetic resonance imaging acquisition | |
CN114983389B (en) | Quantitative evaluation method for human brain axon density based on magnetic resonance diffusion tensor imaging | |
CN105759233B (en) | A kind of rapid chemical exchanges saturation transfer imaging method and system | |
CN108596900B (en) | Thyroid-associated ophthalmopathy medical image data processing device and method, computer-readable storage medium and terminal equipment | |
CN111640097A (en) | Skin mirror image identification method and equipment | |
CN115861228A (en) | Brain glioma segmentation method based on UNet network and attention mechanism | |
CN114723937A (en) | Method and system for classifying blood vessel surrounding gaps based on nuclear magnetic resonance image | |
CN114782371A (en) | Brain area automatic segmentation method based on artificial intelligence | |
CN111723521B (en) | Method for recovering limited caliber data into full caliber data based on GAN network | |
CN115186564A (en) | Analog circuit fault diagnosis method based on feature fusion and improved particle swarm algorithm | |
CN113553991A (en) | Tongue picture recognition method based on deep learning | |
CN114519367A (en) | Motor imagery electroencephalogram frequency characteristic analysis method and system based on sample learning | |
CN106611107A (en) | Method for removing noise of sequencing data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |