CN113625209A - Method and device for determining frequency drift amount of magnetic resonance system and computer equipment - Google Patents
Method and device for determining frequency drift amount of magnetic resonance system and computer equipment Download PDFInfo
- Publication number
- CN113625209A CN113625209A CN202010387615.XA CN202010387615A CN113625209A CN 113625209 A CN113625209 A CN 113625209A CN 202010387615 A CN202010387615 A CN 202010387615A CN 113625209 A CN113625209 A CN 113625209A
- Authority
- CN
- China
- Prior art keywords
- frequency
- echo signal
- navigation
- echo
- magnetic resonance
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000002592 echocardiography Methods 0.000 claims abstract description 19
- 238000004590 computer program Methods 0.000 claims description 26
- 238000006243 chemical reaction Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 abstract description 6
- 238000003384 imaging method Methods 0.000 abstract description 3
- 210000000481 breast Anatomy 0.000 description 15
- 238000010586 diagram Methods 0.000 description 7
- 238000012937 correction Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000000593 degrading effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
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/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
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The application relates to a method and a device for determining frequency drift amount of a magnetic resonance system and computer equipment. The method comprises the following steps: performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; the preset scanning sequence is a gradient double-echo sequence; after the preset duration of the first scanning, performing second scanning by adopting the preset scanning sequence to obtain a second group of echo signals; calculating according to the first group of echo signals and the second group of echoes respectively to obtain a first navigation frequency and a second navigation frequency; and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency. By adopting the method, the accurate frequency drift of the magnetic resonance system can be obtained, so that the accuracy of the system frequency is improved, and the imaging effect is improved.
Description
Technical Field
The present invention relates to the field of magnetic resonance technology, and in particular, to a method, an apparatus, and a computer device for determining a frequency drift amount of a magnetic resonance system.
Background
In the magnetic resonance imaging process, the accuracy of the system frequency directly affects the image quality. When a long-time sequence scan is performed, the magnetic resonance system generates heat due to vibration caused by continuous switching of current, and the system frequency shifts.
In the related art, in order to avoid the system frequency drift from degrading the image quality, the system frequency is usually calibrated in real time. However, when scanning the breast or neck, it is often necessary for the user to manually adjust the system frequency. At this time, if the calibration is still performed according to the previous calibration method, the manual adjustment amount of the system frequency is considered as a frequency drift, thereby causing a calibration error of the system frequency.
Disclosure of Invention
In view of the above, there is a need to provide a method, an apparatus and a computer device for determining a frequency drift amount of a magnetic resonance system, which can avoid taking a manual adjustment amount of the system frequency as a frequency drift, thereby avoiding a system frequency calibration error.
A method of determining an amount of frequency drift of a magnetic resonance system, the method comprising:
performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals;
respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
In one embodiment, each set of echo signals includes a first echo signal and a second echo signal, and the obtaining a first navigation frequency and a second navigation frequency according to the first echo signal and the second echo signal by computing includes:
respectively carrying out frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion;
and calculating the navigation frequency according to the third echo signal and the fourth echo signal.
In one embodiment, the calculating a navigation frequency according to the third echo signal and the fourth echo signal includes:
integrating the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle;
carrying out phase fitting on the system echo signal to obtain a first phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
In one embodiment, the calculating a navigation frequency according to the third echo signal and the fourth echo signal includes:
calculating according to the third echo signal, the fourth echo signal and a preset phase angle function to obtain a second phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
In one embodiment, the correspondence between the frequency and the phase angle includes that the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal.
In one embodiment, the obtaining the first navigation frequency and the second navigation frequency according to the first echo signal and the second echo signal includes:
aiming at the first group of echo signals, calculating navigation frequency according to a plurality of pairs of echo signals to obtain a plurality of candidate navigation frequencies;
an average of the plurality of candidate navigation frequencies is determined as the first navigation frequency.
In one embodiment, the determining the magnetic resonance system frequency drift amount according to the first navigation frequency and the second navigation frequency includes:
calculating a frequency difference between the first navigation frequency and the second navigation frequency;
and performing unwrapping processing on the frequency difference to obtain the frequency drift amount of the magnetic resonance system.
In one embodiment, after determining the magnetic resonance system frequency drift amount according to the first navigator frequency and the second navigator frequency, the method further comprises:
and correcting the frequency of the magnetic resonance system according to the frequency drift amount of the magnetic resonance system.
In one embodiment, after determining the magnetic resonance system frequency drift amount according to the first navigator frequency and the second navigator frequency, the method further comprises:
the frequency of the transmit coil and/or the frequency of the receive coil is modified based on the amount of magnetic resonance system frequency drift.
An apparatus for determining an amount of frequency drift of a magnetic resonance system, the apparatus comprising:
the first sequence scanning module is used for carrying out first scanning by adopting a preset scanning sequence so as to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
the second sequence scanning module is used for performing second scanning by adopting a preset scanning sequence after the preset duration of the first scanning so as to obtain a second group of echo signals;
the navigation frequency calculation module is used for calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and the drift amount determining module is used for determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
In one embodiment, each set of echo signals includes a first echo signal and a second echo signal, and the navigation frequency calculation module includes:
the frequency domain conversion sub-module is used for respectively performing frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion;
and the navigation frequency calculation sub-module is used for calculating the navigation frequency according to the third echo signal and the fourth echo signal.
In one embodiment, the navigation frequency calculation submodule is configured to integrate the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle; carrying out phase fitting on the system echo signal to obtain a first phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
In one embodiment, the navigation frequency calculation sub-module is configured to calculate according to the third echo signal, the fourth echo signal, and a preset phase angle function to obtain a second phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
In one embodiment, the correspondence between the frequency and the phase angle includes that the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal.
In one embodiment, the first set of echo signals includes a plurality of pairs of echo signals, each pair of echo signals includes two echo signals, and the first navigation frequency calculation module is specifically configured to calculate a navigation frequency according to the plurality of pairs of echo signals, so as to obtain a plurality of candidate navigation frequencies; an average of the plurality of candidate navigation frequencies is determined as the first navigation frequency.
In one embodiment, the drift amount determining module is specifically configured to calculate a frequency difference between a first navigation frequency and a second navigation frequency; and performing unwrapping processing on the frequency difference to obtain the frequency drift amount of the magnetic resonance system.
In one embodiment, the apparatus further comprises:
and the navigation frequency correction module is used for correcting the frequency of the magnetic resonance system according to the frequency drift amount of the magnetic resonance system.
In one embodiment, the apparatus further comprises:
and the coil frequency modifying module is used for modifying the frequency of the transmitting coil and/or the frequency of the receiving coil according to the frequency drift amount of the magnetic resonance system.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals;
respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals;
respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
According to the method, the device and the computer equipment for determining the frequency drift amount of the magnetic resonance system, a preset scanning sequence is adopted for carrying out first scanning so as to obtain a first group of echo signals; after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals; respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency; and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency. According to the embodiment of the application, after a user manually adjusts the system frequency, a preset scanning sequence is adopted for carrying out first scanning and obtaining a first group of echo signals, after a preset duration, a preset scanning sequence is adopted for carrying out second scanning and obtaining a second group of echo signals, a first navigation frequency and a second navigation frequency are respectively calculated according to the first group of echo signals and the second group of echo signals, and because the first navigation frequency and the second navigation frequency both contain manual adjustment quantity of the system frequency, the manual adjustment quantity of the system frequency can be offset through calculation, so that the problem that the manual adjustment quantity of the system frequency is used as frequency drift to cause system frequency calibration errors is avoided.
Drawings
FIG. 1 is a diagram of an exemplary embodiment of a method for determining an amount of frequency drift in a magnetic resonance system;
FIG. 2 is a flow chart illustrating a method for determining an amount of frequency drift of a magnetic resonance system according to one embodiment;
FIG. 3 is a waveform diagram in one embodiment;
FIG. 4 is a schematic flow chart illustrating the step of calculating a navigation frequency based on an echo signal according to an embodiment;
FIG. 5 is a flow chart illustrating a method for determining an amount of frequency drift of an MR system according to another embodiment;
FIG. 6 is a multi-phase breast magnetic resonance image obtained using a prior art method;
FIG. 7 is a multi-phase breast magnetic resonance image obtained after frequency correction of the magnetic resonance system in one embodiment;
FIG. 8 is a block diagram showing an example of an apparatus for determining the amount of frequency drift of the magnetic resonance system;
FIG. 9 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The method for determining the frequency drift amount of the magnetic resonance system can be applied to the application environment shown in fig. 1. The application environment is a magnetic resonance system comprising at least a gradient system 101, a radio frequency system 102 and a control terminal 103. Wherein the gradient system 101 comprises devices such as gradient coils for generating gradient pulses; the radio frequency system 102 includes devices such as a radio frequency coil, the radio frequency coil includes a transmitting coil and a receiving coil, and the radio frequency system 102 is configured to generate a radio frequency pulse signal and detect a magnetic resonance signal generated by nuclear spins of the subject; the control terminal 103 communicates with the gradient system 101 and the radio frequency system 102 through a network, and is used for controlling the gradient system 101 and the radio frequency system 102. The control terminal 103 may be, but is not limited to, various personal computers, notebook computers, and tablet computers.
In one embodiment, as shown in fig. 2, a method for determining a frequency drift amount of a magnetic resonance system is provided, which is described by taking the method as an example for being applied to the control terminal in fig. 1, and includes the following steps:
Since the system frequency is often manually adjusted by the user when scanning the breast or neck, the manual adjustment of the system frequency is regarded as a frequency drift in the real-time calibration, which results in an error in the calibration of the system frequency. To avoid this, the embodiment of the present application does not perform real-time frequency calibration any more, but employs dynamic frequency calibration after detecting an adjustment operation for the system frequency, that is, after a user manually adjusts the system frequency. Specifically, a preset scan sequence is first adopted to perform the first scan, wherein the preset scan sequence may be a gradient dual-echo sequence, as shown in the waveform diagram of fig. 3, RF is a radio frequency pulse, GSSFor gradient in the direction of the selected layer, GROFor frequency encoding directional gradients, the ADC represents the signal received by the analog-to-digital converter. The gradient dual echo sequence shown in fig. 3 has only frequency encoding gradients and no phase encoding gradients. In addition, in order to reduce the influence of the preset scanning sequence on the magnetic field, the excitation flip angle is as small as possible under the condition of meeting the signal-to-noise ratio, and the gradient net area of each logic axis is zero.
While performing the first scan using the predetermined scan sequence, a first set of echo signals, FIG. 3, is acquired, along GSSApplying a layer selection gradient corresponding to the radio frequency pulse in the direction and a dephasing gradient applied immediately after the layer selection gradient; along GROApplying a continuous frequency encoding gradient in a direction to avoid gradient eddy currents caused by frequent switching of the gradient; ADC samplingObtaining a first echo signal e0And a second echo signal e1G corresponding to the two echoesROThe gradient moments in the directions are the same; echo time TE (echo time) is the time from the midpoint of the RF pulse to the first echo signal e0Δ TE is the first echo signal e0And a second echo signal e1The echo time difference of (2). In this embodiment, the predetermined scan sequence does not include a gradient in the phase encoding direction, i.e. the first echo signal e is sampled in the embodiment of the present application0And a second echo signal e1Is not phase encoded. In one embodiment, the value of TE is sufficient to ensure that the water signal and the fat signal of the test subject are in phase. In one embodiment, the system frequency is calibrated once before the user manually adjusts the system frequency, so that the system frequency can be corrected more accurately after the magnetic resonance system frequency drift amount is determined. Wherein the system frequency is the larmor frequency, i.e. the precession frequency of the magnetic nuclei or the main magnetic field frequency.
In the embodiment of the present application, after the first group of echo signals are acquired, the system frequency may drift after a period of time, so that the preset scanning sequence is adopted to perform the second scanning, and the second group of echo signals are acquired at the same time.
For example, a scan is performed once using the preset scan sequence, another scan is performed after 5 minutes using the preset scan sequence, and a second set of echo signals is acquired. The preset duration is not limited in detail in the embodiment of the application, and can be set according to actual conditions.
And 203, respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency.
In the embodiment of the application, a first navigation frequency is obtained by calculation according to a first group of echo signals, and a second navigation frequency is obtained by calculation according to a second group of echo signals. Specifically, each group of echo signals is subjected to frequency domain conversion, then a phase angle is calculated according to the converted echo signals, and then navigation frequency is calculated according to the corresponding relationship between the phase angle and the frequency. The specific calculation method is not limited in detail in the embodiment of the application, and can be set according to actual conditions.
And step 204, determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
In the embodiment of the application, the first navigation frequency includes a manual adjustment amount of the system frequency, and the second navigation frequency also includes a manual adjustment amount of the system frequency, so that after the first navigation frequency and the second navigation frequency are obtained, a frequency difference between the first navigation frequency and the second navigation frequency can be calculated, so that the manual adjustment amount of the system frequency is cancelled, and further, a frequency drift amount of the magnetic resonance system is obtained.
In the method for determining the frequency drift amount of the magnetic resonance system, a preset scanning sequence is adopted to perform first scanning so as to obtain a first group of echo signals; after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals; respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency; and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency. In the embodiment of the application, the first navigation frequency and the second navigation frequency both include the manual adjustment amount of the system frequency, so that the manual adjustment amount of the system frequency can be cancelled by calculation, and the problem that the manual adjustment amount of the system frequency is taken as frequency drift to cause system frequency calibration error is avoided.
In one embodiment, as shown in FIG. 4, an optional process involving calculating a navigation frequency from an echo signal. On the basis of the above embodiment, each set of echo signals includes a first echo signal and a second echo signal, and therefore, the following steps may be adopted to calculate the first navigation frequency and the second navigation frequency:
In the embodiment of the application, a gradient double-echo sequence is adopted for scanning, and two echo signals, namely a first echo signal and a second echo signal, can be acquired. And then, performing frequency domain conversion on the first echo signal to obtain a third echo signal, and performing frequency domain conversion on the second echo signal to obtain a fourth echo signal. Specifically, the frequency domain transformation may employ a fourier transform, such as equation (1) (2):
E0＝FFT(e0)-------------------------------------------(1)
E1＝FFT(e1)--------------------------------------------(2)
wherein e is0Is a first echo signal, e1Is the second echo signal, E0Is the third echo signal, E1Is the fourth echo signal.
And step 302, calculating a navigation frequency according to the third echo signal and the fourth echo signal.
In the embodiment of the present application, after the third echo signal and the fourth echo signal in the frequency domain are obtained, the navigation frequency may be calculated in the following two ways.
The first method is as follows: integrating the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle; carrying out phase fitting on the system echo signal to obtain a first phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
Specifically, two echo signals are integrated into one system echo signal, wherein the system echo signal is represented by amplitude and phase angle. The integration process may be as in formula (3):
then, let E (i) be written as formula (4):
wherein E is0(i)*E1(i) Is the third echo signal E0And a fourth echo signal E1Conjugate multiplication; a. theiThe amplitude of the system echo signal e (i),
the phase angle of the system echo signal E (i) is shown, n represents the number of the acquired gradient double echoes, and n is more than or equal to 1.
After obtaining the system echo signal, the phase fitting can be performed by using the least square method, as shown in formula (5):
wherein, theta1Is the fitted first phase angle.
Then, calculating the navigation frequency according to the corresponding relationship between the frequency and the phase angle and the first phase angle obtained by fitting, as specifically shown in formula (6):
the correspondence between the frequency and the phase angle includes: the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal. Where f is the calculated navigation frequency and Δ TE is the time difference between the first echo signal and the second echo signal.
The second method comprises the following steps: calculating according to the third echo signal, the fourth echo signal and a preset phase angle function to obtain a second phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
Specifically, the second phase angle can be obtained by directly calculating according to the phase angle function, as shown in formula (7):
θ2＝angle(∑i(E0(i)*E1(i)))-----------------------(7)
wherein, theta1Is the second phase angle.
Then, the navigation frequency is calculated according to the corresponding relationship between the frequency and the phase angle and the calculated second phase angle, as shown in formula (8):
the correspondence between the frequency and the phase angle includes: the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal. Where f is the calculated navigation frequency and Δ TE is the time difference between the first echo signal and the second echo signal.
In the process of calculating the navigation frequency according to the echo signals, respectively performing frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion; and calculating the navigation frequency according to the third echo signal and the fourth echo signal. According to the embodiment of the application, the echo signal represented by the time domain is converted into the echo signal represented by the frequency domain, the phase angle can be calculated from the echo signal represented by the frequency domain, and the navigation frequency can be calculated according to the corresponding relation between the frequency and the phase angle, so that the first navigation frequency and the second navigation frequency are determined, and the frequency drift amount of the magnetic resonance system is obtained.
In one embodiment, as shown in fig. 5, a method for determining a frequency drift amount of a magnetic resonance system is provided, which is described by taking the method as an example for being applied to the control terminal in fig. 1, and includes the following steps:
And 402, after the preset duration of the first scanning, performing a second scanning by using a preset scanning sequence to obtain a second group of echo signals.
And 403, respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency.
In the embodiment of the application, in order to make the first navigation frequency more accurate, for the first group of echo signals, the navigation frequency can be calculated according to a plurality of pairs of echo signals to obtain a plurality of candidate navigation frequencies; an average of the plurality of candidate navigation frequencies is determined as the first navigation frequency. Specifically, a plurality of pairs of echo signals are collected, wherein each pair of echo signals comprises two echo signals; calculating navigation frequency aiming at each pair of echo signals to obtain corresponding candidate navigation frequency, so that a plurality of candidate navigation frequencies can be obtained; and finally, calculating the average value of the candidate navigation frequencies, and taking the average value as the first navigation frequency.
In one embodiment, calculating a navigation frequency for each pair of echo signals, and obtaining a corresponding candidate navigation frequency may include: respectively carrying out frequency domain conversion on the first echo signal and the second echo signal to obtain a third echo signal and a fourth echo signal after conversion; and calculating the navigation frequency according to the third echo signal and the fourth echo signal.
In one embodiment, calculating the navigation frequency from the third echo signal and the fourth echo signal comprises: integrating the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle; carrying out phase fitting on the system echo signal to obtain a first phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
In one embodiment, calculating the navigation frequency from the third echo signal and the fourth echo signal comprises: calculating according to the third echo signal, the fourth echo signal and a preset phase angle function to obtain a second phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
In one embodiment, the correspondence between frequency and phase angle includes that the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal.
The second set of echo signals includes a pair of echo signals, so that the calculation of the second navigation frequency can be referred to the above embodiments, and is not described herein again.
In the embodiment of the application, after the first navigation frequency and the second navigation frequency are obtained, the frequency difference between the first navigation frequency and the second navigation frequency is calculated to offset the manual adjustment amount of the system frequency. However, due to the phase winding problem, the frequency difference cannot directly reflect the frequency drift amount of the magnetic resonance system, and therefore, the frequency difference needs to be subjected to unwrapping processing. Specifically, the formula (9):
wherein, Deltaf is the frequency drift amount of the magnetic resonance system, fkIs the second navigation frequency, frIs a first navigation frequency. And, setting a frequency threshold value to
n is a positive integer and takes an appropriate value such that
After the amount of frequency drift of the magnetic resonance system is obtained, any one of the following steps may be employed.
In the embodiment of the present application, after the magnetic resonance system frequency drift amount is obtained, the magnetic resonance system frequency can be corrected according to the magnetic resonance system frequency drift amount.
In practical application, as shown in fig. 6, a multi-phase breast magnetic resonance image obtained by a method in the prior art is used, a breast phase 6 scan is performed, if the magnetic resonance system frequency drift amount is not calculated by adopting the above steps to perform system frequency correction, a first row in the figure sequentially corresponds to a first phase breast magnetic resonance image, a second phase breast magnetic resonance image and a third phase breast magnetic resonance image from left to right, a white artifact region (shown by an arrow in the figure) caused by uneven fat compression exists in the third phase breast magnetic resonance image, a second row sequentially corresponds to a fourth phase breast magnetic resonance image, a fifth phase breast magnetic resonance image and a sixth phase breast magnetic resonance image from left to right, and an obvious white artifact region caused by uneven fat compression exists in the three phases, that is, the fat compression effect gradually deteriorates with the increase of the phase. And as shown in fig. 7, the magnetic resonance image of the breast with multiple phases obtained after the frequency correction of the magnetic resonance system is performed by the method shown in fig. 2, which is also a 6-phase scan of the breast, and the frequency drift amount of the magnetic resonance system calculated by the above steps is used for performing the system frequency correction, so that the fat pressing effects of different phases are basically consistent. Therefore, the method for calculating the frequency drift amount of the magnetic resonance system to correct the system frequency can improve the accuracy of the system frequency and further improve the imaging effect.
The frequency of the transmit coil and/or the frequency of the receive coil is modified based on the amount of magnetic resonance system frequency drift, step 406.
In the embodiment of the application, after the frequency drift amount of the magnetic resonance system is obtained, the frequency of the transmitting coil can be corrected, or the frequency of the receiving coil can be corrected; or the frequency of the transmitting coil and the frequency of the receiving coil are simultaneously corrected, so that the effect of correcting the system frequency is achieved.
In the method for determining the frequency drift amount of the magnetic resonance system, a preset scanning sequence is adopted to perform first scanning so as to obtain a first group of echo signals; after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals; respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency; calculating a frequency difference between the first navigation frequency and the second navigation frequency; carrying out unwrapping processing on the frequency difference value to obtain the frequency drift amount of the magnetic resonance system; the navigation frequency can be corrected according to the frequency drift amount of the magnetic resonance system; the frequency of the transmit coil and/or the frequency of the receive coil may also be modified based on the amount of magnetic resonance system frequency drift. According to the embodiment of the application, the frequency drift amount of the magnetic resonance system is determined according to the first navigation frequency and the second navigation frequency, and the system frequency is corrected according to the frequency drift amount of the magnetic resonance system, so that the accuracy of the system frequency can be improved, and the imaging effect is further improved.
It should be understood that although the various steps in the flowcharts of fig. 2-5 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 2-5 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed in turn or alternately with other steps or at least some of the other steps or stages.
In one embodiment, as shown in fig. 8, there is provided an apparatus for determining a frequency drift amount of a magnetic resonance system, comprising:
a first sequence scanning module 501, configured to perform a first scanning with a preset scanning sequence to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
a second sequence scanning module 502, configured to perform a second scanning by using a preset scanning sequence after a preset duration of the first scanning to obtain a second group of echo signals;
a navigation frequency calculation module 503, configured to calculate according to the first group of echo signals and the second group of echoes respectively to obtain a first navigation frequency and a second navigation frequency;
and a drift amount determining module 504, configured to determine a magnetic resonance system frequency drift amount according to the first navigation frequency and the second navigation frequency.
In one embodiment, each set of echo signals includes a first echo signal and a second echo signal, and the navigation frequency calculation module 503 includes:
the frequency domain conversion sub-module is used for respectively performing frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion;
and the navigation frequency calculation sub-module is used for calculating the navigation frequency according to the third echo signal and the fourth echo signal.
In one embodiment, the navigation frequency calculation submodule is configured to integrate the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle; carrying out phase fitting on the system echo signal to obtain a first phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
In one embodiment, the navigation frequency calculation sub-module is configured to calculate according to the third echo signal, the fourth echo signal, and a preset phase angle function to obtain a second phase angle; and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
In one embodiment, the correspondence between the frequency and the phase angle includes that the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal.
In one embodiment, the first set of echo signals includes a plurality of pairs of echo signals, each pair of echo signals includes two echo signals, and the navigation frequency calculation module 503 is specifically configured to calculate a navigation frequency according to the plurality of pairs of echo signals, so as to obtain a plurality of candidate navigation frequencies; an average of the plurality of candidate navigation frequencies is determined as the first navigation frequency.
In one embodiment, the drift amount determining module 504 is specifically configured to calculate a frequency difference between a first navigation frequency and a second navigation frequency; and performing unwrapping processing on the frequency difference to obtain the frequency drift amount of the magnetic resonance system.
In one embodiment, the apparatus further comprises:
and the system frequency correction module is used for correcting the frequency of the magnetic resonance system according to the frequency drift amount of the magnetic resonance system.
In one embodiment, the apparatus further comprises:
and the coil frequency modifying module is used for modifying the frequency of the transmitting coil and/or the frequency of the receiving coil according to the frequency drift amount of the magnetic resonance system.
The specific definition of the determining device for the magnetic resonance system frequency drift amount can refer to the above definition of the determining method for the magnetic resonance system frequency drift amount, and is not described herein again. The modules in the device for determining the frequency drift amount of the magnetic resonance system can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 9. The computer device includes a processor, a memory, a communication interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless communication can be realized through WIFI, an operator network, NFC (near field communication) or other technologies. The computer program is executed by a processor to implement a method of determining an amount of frequency drift of a magnetic resonance system. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 9 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory and a processor, the memory having a computer program stored therein, the processor implementing the following steps when executing the computer program:
performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals;
respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
In one embodiment, each set of echo signals includes a first echo signal and a second echo signal, and the processor, when executing the computer program, implements the following steps:
respectively carrying out frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion;
and calculating the navigation frequency according to the third echo signal and the fourth echo signal.
In one embodiment, the processor, when executing the computer program, performs the steps of:
integrating the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle;
carrying out phase fitting on the system echo signal to obtain a first phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
In one embodiment, the processor, when executing the computer program, performs the steps of:
calculating according to the third echo signal, the fourth echo signal and a preset phase angle function to obtain a second phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
In one embodiment, the correspondence between the frequency and the phase angle includes that the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal.
In one embodiment, the first set of echo signals comprises a plurality of pairs of echo signals, each pair of echo signals comprising two echo signals, the processor when executing the computer program performs the steps of:
aiming at the first group of echo signals, calculating navigation frequency according to a plurality of pairs of echo signals to obtain a plurality of candidate navigation frequencies;
an average of the plurality of candidate navigation frequencies is determined as the first navigation frequency.
In one embodiment, the processor, when executing the computer program, performs the steps of:
calculating a frequency difference between the first navigation frequency and the second navigation frequency;
and performing unwrapping processing on the frequency difference to obtain the frequency drift amount of the magnetic resonance system.
In one embodiment, the processor, when executing the computer program, performs the steps of:
and correcting the frequency of the magnetic resonance system according to the frequency drift amount of the magnetic resonance system.
In one embodiment, the processor, when executing the computer program, performs the steps of:
the frequency of the transmit coil and/or the frequency of the receive coil is modified based on the amount of magnetic resonance system frequency drift.
In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of:
performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; presetting a scanning sequence as a gradient double-echo sequence;
after the preset duration of the first scanning, performing second scanning by adopting a preset scanning sequence to obtain a second group of echo signals;
respectively calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
In an embodiment, each set of echo signals comprises a first echo signal and a second echo signal, and the computer program when executed by the processor performs the steps of:
respectively carrying out frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion;
and calculating the navigation frequency according to the third echo signal and the fourth echo signal.
In one embodiment, the computer program when executed by the processor implements the steps of:
integrating the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle;
carrying out phase fitting on the system echo signal to obtain a first phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
In one embodiment, the computer program when executed by the processor implements the steps of:
calculating according to the third echo signal, the fourth echo signal and a preset phase angle function to obtain a second phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
In one embodiment, the correspondence between the frequency and the phase angle includes that the frequency is proportional to the phase angle and inversely proportional to the time difference between the first echo signal and the second echo signal.
In an embodiment, the first set of echo signals comprises a plurality of pairs of echo signals, each pair of echo signals comprising two echo signals, the computer program when executed by the processor implementing the steps of:
aiming at the first group of echo signals, calculating navigation frequency according to a plurality of pairs of echo signals to obtain a plurality of candidate navigation frequencies;
an average of the plurality of candidate navigation frequencies is determined as the first navigation frequency.
In one embodiment, the computer program when executed by the processor implements the steps of:
calculating a frequency difference between the first navigation frequency and the second navigation frequency;
and performing unwrapping processing on the frequency difference to obtain the frequency drift amount of the magnetic resonance system.
In one embodiment, the computer program when executed by the processor implements the steps of:
and correcting the frequency of the magnetic resonance system according to the frequency drift amount of the magnetic resonance system.
In one embodiment, the computer program when executed by the processor implements the steps of:
the frequency of the transmit coil and/or the frequency of the receive coil is modified based on the amount of magnetic resonance system frequency drift.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A method for determining an amount of frequency drift in a magnetic resonance system, the method comprising:
performing first scanning by adopting a preset scanning sequence to obtain a first group of echo signals; the preset scanning sequence is a gradient double-echo sequence;
after the preset duration of the first scanning, performing second scanning by adopting the preset scanning sequence to obtain a second group of echo signals;
calculating according to the first group of echo signals and the second group of echoes respectively to obtain a first navigation frequency and a second navigation frequency;
and determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
2. The method of claim 1, wherein each set of echo signals comprises a first echo signal and a second echo signal, and wherein the calculating from the first set of echo signals and the second set of echoes respectively to obtain a first navigation frequency and a second navigation frequency comprises:
respectively carrying out frequency domain conversion on the first echo signal and the second echo signal aiming at each group of echo signals to obtain a third echo signal and a fourth echo signal after conversion;
and calculating navigation frequency according to the third echo signal and the fourth echo signal.
3. The method of claim 2, wherein said calculating a navigation frequency from said third echo signal and said fourth echo signal comprises:
integrating the third echo signal and the fourth echo signal to obtain a system echo signal; the system echo signal is represented by amplitude and phase angle;
carrying out phase fitting on the system echo signal to obtain a first phase angle;
and calculating navigation frequency according to the corresponding relation between the frequency and the phase angle and the first phase angle.
4. The method of claim 2, wherein said calculating a navigation frequency from said third echo signal and said fourth echo signal comprises:
calculating according to the third echo signal, the fourth echo signal and a preset phase angle function to obtain a second phase angle;
and calculating the navigation frequency according to the corresponding relation between the frequency and the phase angle and the second phase angle.
5. The method of claim 3 or 4, wherein the correspondence between frequency and phase angle comprises a frequency proportional to phase angle and an inverse proportional to a time difference between the first echo signal and the second echo signal.
6. The method of claim 1, wherein the first set of echo signals comprises a plurality of pairs of echo signals, each pair of echo signals comprising two echo signals, and wherein the calculating from the first set of echo signals and the second set of echoes, respectively, to obtain a first navigation frequency and a second navigation frequency comprises:
aiming at the first group of echo signals, calculating navigation frequency according to the multiple pairs of echo signals to obtain multiple candidate navigation frequencies;
determining an average of a plurality of the candidate navigation frequencies as the first navigation frequency.
7. The method of claim 1, wherein determining an amount of magnetic resonance system frequency drift from the first and second navigator frequencies comprises:
calculating a frequency difference between the first navigation frequency and the second navigation frequency;
and performing unwrapping processing on the frequency difference value to obtain the frequency drift amount of the magnetic resonance system.
8. The method of claim 1, wherein after said determining an amount of magnetic resonance system frequency drift from said first and second navigator frequencies, the method further comprises:
correcting the frequency of the magnetic resonance system according to the frequency drift amount of the magnetic resonance system; or,
and correcting the frequency of the transmitting coil and/or the frequency of the receiving coil according to the frequency drift amount of the magnetic resonance system.
9. An apparatus for determining an amount of frequency drift in a magnetic resonance system, the apparatus comprising:
the first sequence scanning module is used for carrying out first scanning by adopting a preset scanning sequence so as to obtain a first group of echo signals; the preset scanning sequence is a gradient double-echo sequence;
the second sequence scanning module is used for performing second scanning by adopting the preset scanning sequence after the preset duration of the first scanning so as to obtain a second group of echo signals;
the navigation frequency calculation module is used for calculating according to the first group of echo signals and the second group of echoes to obtain a first navigation frequency and a second navigation frequency;
and the drift amount determining module is used for determining the frequency drift amount of the magnetic resonance system according to the first navigation frequency and the second navigation frequency.
10. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the method of any of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010387615.XA CN113625209B (en) | 2020-05-09 | 2020-05-09 | Method and device for determining frequency drift amount of magnetic resonance system and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010387615.XA CN113625209B (en) | 2020-05-09 | 2020-05-09 | Method and device for determining frequency drift amount of magnetic resonance system and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113625209A true CN113625209A (en) | 2021-11-09 |
| CN113625209B CN113625209B (en) | 2024-02-27 |
Family
ID=78377490
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010387615.XA Active CN113625209B (en) | 2020-05-09 | 2020-05-09 | Method and device for determining frequency drift amount of magnetic resonance system and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113625209B (en) |
Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4728906A (en) * | 1985-10-21 | 1988-03-01 | Wiltron Measurements Limited | Tuning and calibration circuits for CW and sweep frequency signal generators |
| CN1216242A (en) * | 1997-10-17 | 1999-05-12 | 日立医药株式会社 | Magnetic resonance imaging method and device therefor |
| CN1693885A (en) * | 2004-04-30 | 2005-11-09 | 西门子公司 | Method for MR dynamic frequency detection in MR spectroscopy |
| CN101034149A (en) * | 2007-04-19 | 2007-09-12 | 新奥博为技术有限公司 | Method for correcting magnetic resonance signal phase excursion based on pre-scanning survey |
| CN102078196A (en) * | 2009-11-27 | 2011-06-01 | 西门子公司 | Method to detect breathing movement of an examination subject corresponding to signal data by magnetic resonance |
| US20110267057A1 (en) * | 2010-04-30 | 2011-11-03 | Qiang He | Dynamic frequency drift correction method in magnetic resonance imaging |
| US20110304333A1 (en) * | 2009-01-13 | 2011-12-15 | Aspect Magnet Technologies Ltd. | Means And Methods For Providing High Resolution MRI |
| CN102654574A (en) * | 2012-05-21 | 2012-09-05 | 西安电子科技大学 | Doppler centroid estimation method based on field programmable gate array (FPGA) |
| US20120249138A1 (en) * | 2011-03-28 | 2012-10-04 | Josef Pfeuffer | Magnetic resonance method and system to correct phase information in mr images |
| US20120249143A1 (en) * | 2011-03-31 | 2012-10-04 | Masaaki Umeda | Magnetic resonance imaging apparatus |
| CN102866372A (en) * | 2011-07-08 | 2013-01-09 | 西门子公司 | Methods for calibrating frequency of magnetic resonance device and corresponding magnetic resonance device |
| CN102944875A (en) * | 2012-11-21 | 2013-02-27 | 西安电子科技大学 | Inverse synthetic aperture radar (ISAR) image distance unit selecting and cross-range scaling method |
| US20150035530A1 (en) * | 2013-08-02 | 2015-02-05 | Matthew Tisdall | System and method for real-time frequency correction for magnetic resonance imaging |
| CN105247382A (en) * | 2013-03-29 | 2016-01-13 | 皇家飞利浦有限公司 | Magnetic resonance imaging system with ghost artifact reduction and method of operation thereof |
| CN105997074A (en) * | 2016-04-28 | 2016-10-12 | 华东师范大学 | Magnetic resonance quantitative magnetic susceptibility imaging multi-echo phase fitting method |
| CN108375746A (en) * | 2017-01-18 | 2018-08-07 | 上海联影医疗科技有限公司 | A kind of phase warp folding method and apparatus |
| CN108627834A (en) * | 2018-06-07 | 2018-10-09 | 北京城建勘测设计研究院有限责任公司 | A kind of subway road structure monitoring method and device based on ground InSAR |
| CN108742623A (en) * | 2018-02-05 | 2018-11-06 | 苏州朗润医疗系统有限公司 | Dynamic method for shimming based on three-dimensional gradient double echo |
| CN110068780A (en) * | 2019-05-07 | 2019-07-30 | 上海东软医疗科技有限公司 | A kind of adaptive approach of magnetic resonance imaging, device and magnetic resonance imaging system |
| CN110907872A (en) * | 2018-09-18 | 2020-03-24 | 西门子医疗有限公司 | Method for generating magnetic resonance images, program, data carrier and magnetic resonance apparatus |
-
2020
- 2020-05-09 CN CN202010387615.XA patent/CN113625209B/en active Active
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4728906A (en) * | 1985-10-21 | 1988-03-01 | Wiltron Measurements Limited | Tuning and calibration circuits for CW and sweep frequency signal generators |
| CN1216242A (en) * | 1997-10-17 | 1999-05-12 | 日立医药株式会社 | Magnetic resonance imaging method and device therefor |
| CN1693885A (en) * | 2004-04-30 | 2005-11-09 | 西门子公司 | Method for MR dynamic frequency detection in MR spectroscopy |
| US20050248343A1 (en) * | 2004-04-30 | 2005-11-10 | Gunnar Kruger | Method and MR apparatus for dynamic frequency detection in MR spectroscopy |
| CN101034149A (en) * | 2007-04-19 | 2007-09-12 | 新奥博为技术有限公司 | Method for correcting magnetic resonance signal phase excursion based on pre-scanning survey |
| US20110304333A1 (en) * | 2009-01-13 | 2011-12-15 | Aspect Magnet Technologies Ltd. | Means And Methods For Providing High Resolution MRI |
| CN102078196A (en) * | 2009-11-27 | 2011-06-01 | 西门子公司 | Method to detect breathing movement of an examination subject corresponding to signal data by magnetic resonance |
| US20110267057A1 (en) * | 2010-04-30 | 2011-11-03 | Qiang He | Dynamic frequency drift correction method in magnetic resonance imaging |
| CN102232849A (en) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | Dynamic frequency drift calibrating method |
| US20120249138A1 (en) * | 2011-03-28 | 2012-10-04 | Josef Pfeuffer | Magnetic resonance method and system to correct phase information in mr images |
| US20120249143A1 (en) * | 2011-03-31 | 2012-10-04 | Masaaki Umeda | Magnetic resonance imaging apparatus |
| CN102866372A (en) * | 2011-07-08 | 2013-01-09 | 西门子公司 | Methods for calibrating frequency of magnetic resonance device and corresponding magnetic resonance device |
| CN102654574A (en) * | 2012-05-21 | 2012-09-05 | 西安电子科技大学 | Doppler centroid estimation method based on field programmable gate array (FPGA) |
| CN102944875A (en) * | 2012-11-21 | 2013-02-27 | 西安电子科技大学 | Inverse synthetic aperture radar (ISAR) image distance unit selecting and cross-range scaling method |
| CN105247382A (en) * | 2013-03-29 | 2016-01-13 | 皇家飞利浦有限公司 | Magnetic resonance imaging system with ghost artifact reduction and method of operation thereof |
| US20150035530A1 (en) * | 2013-08-02 | 2015-02-05 | Matthew Tisdall | System and method for real-time frequency correction for magnetic resonance imaging |
| CN105997074A (en) * | 2016-04-28 | 2016-10-12 | 华东师范大学 | Magnetic resonance quantitative magnetic susceptibility imaging multi-echo phase fitting method |
| CN108375746A (en) * | 2017-01-18 | 2018-08-07 | 上海联影医疗科技有限公司 | A kind of phase warp folding method and apparatus |
| CN108742623A (en) * | 2018-02-05 | 2018-11-06 | 苏州朗润医疗系统有限公司 | Dynamic method for shimming based on three-dimensional gradient double echo |
| CN108627834A (en) * | 2018-06-07 | 2018-10-09 | 北京城建勘测设计研究院有限责任公司 | A kind of subway road structure monitoring method and device based on ground InSAR |
| CN110907872A (en) * | 2018-09-18 | 2020-03-24 | 西门子医疗有限公司 | Method for generating magnetic resonance images, program, data carrier and magnetic resonance apparatus |
| CN110068780A (en) * | 2019-05-07 | 2019-07-30 | 上海东软医疗科技有限公司 | A kind of adaptive approach of magnetic resonance imaging, device and magnetic resonance imaging system |
Non-Patent Citations (2)
| Title |
|---|
| 修峰: "基于动态频率漂移校正的二项式水激发磁共振成像", 《中国优秀硕士学位论文全文数据库 (医药卫生科技辑)》, no. 5, pages 4 * |
| 赵欣欣;薄斌仕;刘田;王乙;李建奇;: "定量磁化率成像多回波相位拟合算法研究", 波谱学杂志, no. 04 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113625209B (en) | 2024-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9046590B2 (en) | Magnetic resonance imaging apparatus phase correction using one or more prescans with identical readout and slice gradients | |
| US9689945B2 (en) | Virtual coil emulation in parallel transmission MRI | |
| JP6585582B2 (en) | Magnetic resonance imaging system capable of reducing ghost artifact and method of operating the same | |
| JP2012205897A (en) | Method for correcting phase information in magnetic resonance image and magnetic resonance equipment | |
| US8854037B2 (en) | Method for compensating for eddy current fields in magnetic resonance images | |
| JPWO2008152937A1 (en) | Magnetic resonance imaging apparatus and method for correcting error caused by gradient magnetic field | |
| US9971007B2 (en) | Method and apparatus for accelerated magnetic resonance imaging | |
| US11585884B2 (en) | Continual trajectory correction in magnetic resonance imaging | |
| US20060253018A1 (en) | Calibrating method for artifact-reduced MRT imaging when there is FOV displacement | |
| US11624793B2 (en) | Method for acquiring measurement data using a magnetic resonance system with correction of k-space trajectories | |
| US10705172B2 (en) | Magnetic resonance apparatus and method for dynamic adjustment thereof with multiple adjustment parameters | |
| US10473743B2 (en) | Method and magnetic resonance apparatus for determining a scan sequence based on a representation of a pulse response in k-space of the gradient system | |
| US11721049B2 (en) | Method for generating at least one image dataset and one reference image dataset, data carrier, computer program product, and magnetic resonance system | |
| CN113625209B (en) | Method and device for determining frequency drift amount of magnetic resonance system and computer equipment | |
| US11965945B2 (en) | Magnetic resonance system and shimming method and imaging method thereof | |
| US11187760B2 (en) | Magnetic resonance imaging apparatus and quantitative-value computing program | |
| US11402454B2 (en) | Correction of distorted diffusion-weighted magnetic resonance image data | |
| CN114325531B (en) | Method, device, computer equipment and medium for correcting delay of magnetic resonance system | |
| CN115586473B (en) | Magnetic field shimming method and device, magnetic resonance imaging equipment and readable storage medium | |
| KR101778139B1 (en) | Method for correcting magnetic field of magnetic resonance imaging device and magnetic resonance imaging device employing the same | |
| CN114089242B (en) | Excitation pulse design method, excitation pulse design device, computer equipment and storage medium | |
| CN115439335A (en) | Image correction method, image correction device, computer equipment and storage medium | |
| US20220206097A1 (en) | Magnetic resonance image processing method and device, and computer readable storage medium | |
| US10823796B2 (en) | Magnetic resonance imaging apparatus and pulse setting method | |
| CN115546080A (en) | Medical image processing method and device, computer equipment and storage medium |
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 |