CN110133559B - Magnetic resonance B0 field disturbance compensation system and method - Google Patents
Magnetic resonance B0 field disturbance compensation system and method Download PDFInfo
- Publication number
- CN110133559B CN110133559B CN201910421708.7A CN201910421708A CN110133559B CN 110133559 B CN110133559 B CN 110133559B CN 201910421708 A CN201910421708 A CN 201910421708A CN 110133559 B CN110133559 B CN 110133559B
- Authority
- CN
- China
- Prior art keywords
- disturbance
- information
- local oscillator
- calculation
- timing
- 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
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention relates to a magnetic resonance B0 field disturbance compensation system and a method, comprising the following steps: the disturbance information primary processing module is connected with the disturbance amount calculation module and used for receiving disturbance information, classifying the disturbance information, performing primary processing on classified data and sending the data after the primary processing to the disturbance amount calculation module; the disturbance quantity calculation module is connected with the local oscillator updating module and used for receiving data transmitted by the disturbance information primary processing module, performing corresponding disturbance quantity calculation according to different types of disturbance information and disturbance control parameters, and outputting the calculated disturbance quantity; the local oscillator updating module is connected with the disturbance quantity calculating module and used for receiving the local oscillator signals and the disturbance quantity transmitted by the disturbance quantity calculating module and updating different modes of the local oscillator signals according to the type of the disturbance quantity. The invention can compensate various disturbances at the same time, so that the local oscillation signal can change along with the disturbance information, and the image quality is improved.
Description
Technical Field
The invention relates to the field of magnetic resonance imaging systems, in particular to a magnetic resonance B0 field disturbance compensation system and a method.
Background
The pursuit of health has prompted the explosive growth of medical health services. Since the 21 st century, magnetic resonance imaging has been widely used due to its advantages of high resolution, large information content, and no damage.
Magnetic resonance imaging systems require the use of magnetic fields. According to lamor's equation omega0=γB0Different magnet B0 fields correspond to different local oscillator signals. Since the imaging quality of a magnetic resonance system is very dependent on the stability of the B0 field, it is very important to maintain the stability of the B0 field. But is provided withSome factors influencing the B0 field stability cannot be completely removed, such as a disturbance of 50 hz of the power supply, a disturbance caused by passing a car near the magnetic resonance imaging system, a transient disturbance caused after a gradient field is applied in the magnetic resonance imaging system, and the like. Therefore, after the B0 field is disturbed, it is very important to compensate the local oscillation for the disturbance.
Currently, many magnetic resonance imaging systems cannot compensate for disturbances, or can only compensate for a certain type of disturbance. So that the image quality deteriorates when a disturbance occurs.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a magnetic resonance B0 field disturbance compensation system and a method thereof, which can simultaneously compensate various disturbances, so that local oscillation signals can change along with disturbance information, and the image quality is improved.
The technical scheme adopted by the invention for realizing the purpose is as follows:
a magnetic resonance B0 field disturbance compensation system, comprising:
the disturbance information primary processing module is connected with the disturbance amount calculation module and used for receiving the disturbance information, classifying the disturbance information, performing primary processing on the classified data and sending the data after the primary processing to the disturbance amount calculation module; the preliminary processing comprises normalization processing and filtering processing of the disturbance information data packet;
the disturbance quantity calculation module is connected with the disturbance information primary processing module and the local oscillator updating module and used for receiving data transmitted by the disturbance information primary processing module, performing corresponding disturbance quantity calculation according to different types of disturbance information and disturbance control parameters and outputting the calculated disturbance quantity;
and the local oscillator updating module is connected with the disturbance quantity calculating module and used for receiving the local oscillator signals and the disturbance quantity transmitted by the disturbance quantity calculating module and updating different modes of the local oscillator signals according to the type of the disturbance quantity.
The classifying the disturbance information includes: dividing the disturbance information into real-time calculation and timing calculation; when the packet type of the disturbance information is 5A5A in 16 system, the disturbance information needs to be calculated in real time, and when the packet type of the disturbance information is A5A5 in 16 system, the disturbance information needs to be calculated in timing.
The normalization processing comprises the following steps: finding the maximum value in the data packet parameters, dividing each parameter in the data packet parameters by the maximum value, and updating the original parameters by the obtained numerical value.
The filtering process is as follows: and filtering the data in the data packet parameters by a high-pass FIR filter.
The disturbance amount calculation process includes: firstly, judging the type of disturbance information, and if the type of the disturbance information needs timing calculation, the formula of the timing calculation is as follows:
wherein f (x) represents the disturbance amount calculated this time, k represents the index of the parameter, akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter data packet, wherein f (x-1) represents the disturbance quantity obtained by the last time of timing calculation, and the last disturbance quantity is 0 during the first calculation;
the calculation formula of the real-time calculation is as follows:
wherein f (x) represents the disturbance amount calculated this time, k represents the index of the parameter, akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter packet.
The local oscillator updating module comprises a local oscillator timing compensation module and a local oscillator real-time compensation module:
a timer for circularly timing is arranged in the local oscillator timing compensation module, for the local oscillator signals of which the disturbance information is timing, the timer adds the timing update disturbance quantity to the local oscillator when the time of the timer is over, then the next timing period is entered again, the timer adds the timing update disturbance quantity to the local oscillator in a circulating manner, and a new local oscillator is transmitted to the local oscillator real-time compensation module;
and the local oscillator real-time compensation module is used for adding the new local oscillator sent by the local oscillator timing compensation module and the disturbance quantity updated in real time to obtain a compensated local oscillator.
A magnetic resonance B0 field disturbance compensation method, comprising the steps of:
step 1: receiving disturbance information of a magnetic resonance B0 field, classifying the disturbance information, and performing primary processing on classified data; the primary processing comprises normalization processing and filtering processing of the disturbance information data packet;
step 2: according to the disturbance information and disturbance control parameters of different types after the preliminary processing, corresponding disturbance quantity calculation is carried out to obtain the calculated disturbance quantity;
and step 3: and updating different modes of the local oscillation signal according to the type of the disturbance quantity.
The classifying the disturbance information includes: dividing the disturbance information into real-time calculation and timing calculation; when the packet type of the disturbance information is 5A5A in 16 system, the disturbance information needs to be calculated in real time, and when the packet type of the disturbance information is A5A5 in 16 system, the disturbance information needs to be calculated in timing.
The disturbance amount calculation process includes: firstly, judging the type of disturbance information, and if the type of the disturbance information needs timing calculation, the formula of the timing calculation is as follows:
where f (x) represents the amount of disturbance in the present calculation, k represents the index of the parameter, and akRepresenting the kth parameter in the packet of perturbation information, bkThe k-th parameter in the disturbance control parameter data packet is represented, f (x-1) represents the disturbance quantity obtained by the last timing calculation, and the last disturbance quantity is 0 in the first calculation. The calculation formula of the real-time calculation is as follows:
where f (x) represents the disturbance amount calculated this time, k represents the index of the parameter, ax represents the kth parameter in the disturbance information packet, bkRepresenting the kth parameter in the disturbance control parameter packet.
The updating of the local oscillator signal in different modes according to the type of the disturbance quantity comprises:
firstly, for local oscillation signals with disturbance information being timing, adding a timing update disturbance quantity to the local oscillation through a timer when the time is over, then entering the next timing period again, and circularly adding the timing update disturbance quantity to the local oscillation to generate a new local oscillation;
and then adding the new local oscillator and the disturbance quantity updated in real time to obtain a compensated local oscillator.
The invention has the following beneficial effects and advantages:
1. the disturbance source can be a multi-channel signal and can simultaneously compensate a plurality of disturbances;
2. the disturbance compensation of the invention supports various types, including linear compensation, exponential compensation and the like;
3. the invention can compensate various disturbances at the same time, so that the local oscillation signal can change along with the disturbance information, and the image quality is improved.
Drawings
FIG. 1 is a block diagram of the system of the present invention;
FIG. 2 is a functional diagram of a disturbance information primary processing module according to the present invention;
FIG. 3 is a functional diagram of a disturbance quantity calculating module according to the present invention;
FIG. 4 is a functional diagram of a local oscillation update module according to the present invention;
FIG. 5 is a flow chart of the method of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in detail below. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein, but rather should be construed as modified in the spirit and scope of the present invention as set forth in the appended claims.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention.
The structure of the system of the invention is shown in figure 1, comprising:
and a disturbance information primary processing module, configured to perform primary processing on the disturbance information, where the function of the disturbance information is implemented as shown in fig. 2. The method comprises the steps of decoding received disturbance information according to a format to obtain a data packet of the disturbance information, wherein the content of the data packet comprises the length of the data packet, the type of the data packet and the parameter of the data packet. Classifying the disturbance information according to the type of the data packet, specifically: when the data packet type is 5A5A in 16 system, the disturbance information needs to be calculated in real time, and when the data packet type is A5A5 in 16 system, the disturbance information needs to be calculated in timing. And carrying out primary treatment such as normalization, filtering and the like on the classified data packet, wherein the primary treatment of normalization is as follows: finding the maximum value in the data packet parameters, dividing each parameter in the data packet parameters by the maximum value, and updating the original parameters by using the obtained numerical values; the primary filtering treatment comprises the following steps: and filtering the data in the data packet parameters by a high-pass FIR filter. And sending the data after the preliminary processing to a disturbance quantity calculation module. The different types include, but are not limited to: disturbance information that varies linearly with time, disturbance information that varies exponentially with time, and disturbance information that varies in real time with time.
And a disturbance amount calculation module, configured to calculate the received data to generate a disturbance amount, where the function of the disturbance amount calculation module is implemented as shown in fig. 3. And calculating the corresponding disturbance amount in the disturbance amount calculation module, including but not limited to: the method comprises the following steps of calculating disturbance quantity linearly changing along with time, calculating disturbance quantity exponentially changing along with time, and calculating disturbance quantity changing along with time in real time; the time-exponentially varying disturbance quantity calculation includes one or more sets of amplitude and time term coefficients. Firstly, judging the type of disturbance information, and then carrying out corresponding timing calculation or real-time calculation according to the disturbance type and the disturbance control parameter. The calculation formula of the timing calculation is as follows:
where f (x) represents the amount of disturbance in the present calculation, k represents the index of the parameter, and akRepresenting the kth parameter in the packet of perturbation information, bkThe k-th parameter in the disturbance control parameter data packet is represented, f (x-1) represents the disturbance quantity obtained by the last timing calculation, and the last disturbance quantity is 0 in the first calculation. The calculation formula of the real-time calculation is as follows:
where f (x) represents the amount of disturbance in the present calculation, k represents the index of the parameter, and akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter packet. And sending the calculated disturbance amount to a local oscillator updating module. The disturbance amount calculation module can simultaneously calculate the multi-path interference information.
And the local oscillator updating module is used for updating the local oscillator corresponding to the magnetic resonance B0 field, and the function of the local oscillator updating module is implemented as shown in fig. 4. The updating of the different modes in the local oscillator update module includes, but is not limited to: updating the local oscillator in real time and updating the local oscillator according to a fixed time interval. The local oscillator updating module is connected with the disturbance quantity calculating module and used for receiving the disturbance quantity transmitted by the local oscillator signal and the disturbance quantity calculating module, and updating different modes of the local oscillator signal is realized according to the type of the disturbance quantity, specifically: a timer which is circularly timed is arranged in the local oscillator timing compensation module, the timing updating disturbance quantity is added to the local oscillator when the time of the timer is over, then the next timing period is entered again, the timing updating disturbance quantity is added to the local oscillator circularly, and a new local oscillator is transmitted to the local oscillator real-time compensation module; the local oscillator real-time compensation module specifically works to add the new local oscillator and the real-time updated disturbance quantity in real time to obtain a compensated local oscillator.
Fig. 5 shows a flow chart of the method of the present invention.
The method comprises the following steps:
step 1: receiving disturbance information of a magnetic resonance B0 field, classifying the disturbance information, and performing primary processing on classified data; the primary processing comprises normalization processing and filtering processing of the disturbance information data packet; the normalization processing comprises the following steps: finding the maximum value in the data packet parameters, dividing each parameter in the data packet parameters by the maximum value, and updating the original parameters by the obtained numerical value. The filtering process is as follows: and filtering the data in the data packet parameters by a high-pass FIR filter.
Step 2: according to the disturbance information and disturbance control parameters of different types after the preliminary processing, corresponding disturbance quantity calculation is carried out to obtain the calculated disturbance quantity;
and step 3: and updating different modes of the local oscillation signal according to the type of the disturbance quantity.
The classifying the disturbance information includes: dividing the disturbance information into real-time calculation and timing calculation; when the packet type of the disturbance information is 5A5A in 16 system, the disturbance information needs to be calculated in real time, and when the packet type of the disturbance information is A5A5 in 16 system, the disturbance information needs to be calculated in timing.
The disturbance amount calculation process includes: firstly, judging the type of disturbance information, and if the type of the disturbance information needs timing calculation, the formula of the timing calculation is as follows:
where f (x) the disturbance amount calculated this time, k denotes an index of a parameter, and akRepresenting the kth parameter in the packet of perturbation information, bkThe k-th parameter in the disturbance control parameter data packet is represented, f (x-1) represents the disturbance quantity obtained by the last timing calculation, and the last disturbance quantity is 0 in the first calculation. The calculation formula of the real-time calculation is as follows:
where f (x) the disturbance amount calculated this time, k denotes an index of a parameter, and akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter packet. The content of the disturbance control parameter comprises the length of a data packet, the type of the data packet and the parameter of the data packet, wherein the type of the data packet represents the type of disturbance information applicable to the disturbance control parameter, namely real-time calculation disturbance information or timing calculation disturbance information.
The updating of the local oscillator signal in different modes according to the type of the disturbance quantity comprises:
firstly, for the timing type disturbance quantity, the timing update disturbance quantity is added to a local oscillator by a timer when the time is over, then the next timing period is entered, and the timing update disturbance quantity is added to the local oscillator in a circulating manner to generate a new local oscillator;
and then adding the new local oscillator and the disturbance quantity updated in real time to obtain a compensated local oscillator.
Claims (8)
1. A magnetic resonance B0 field disturbance compensation system, comprising:
the disturbance information primary processing module is connected with the disturbance amount calculation module and used for receiving the disturbance information, classifying the disturbance information, performing primary processing on the classified data and sending the data after the primary processing to the disturbance amount calculation module; the preliminary processing comprises normalization processing and filtering processing of the disturbance information data packet;
the disturbance quantity calculation module is connected with the disturbance information primary processing module and the local oscillator updating module and used for receiving data transmitted by the disturbance information primary processing module, performing corresponding disturbance quantity calculation according to different types of disturbance information and disturbance control parameters and outputting the calculated disturbance quantity;
the disturbance amount calculation process includes: firstly, judging the type of disturbance information, and if the type of the disturbance information needs timing calculation, the formula of the timing calculation is as follows:
wherein f (x) represents the disturbance amount calculated this time, k represents the index of the parameter, akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter data packet, wherein f (x-1) represents the disturbance quantity obtained by the last time of timing calculation, and the last disturbance quantity is 0 during the first calculation;
the calculation formula of the real-time calculation is as follows:
wherein f (x) represents the disturbance amount calculated this time, k represents the index of the parameter, akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter data packet;
and the local oscillator updating module is connected with the disturbance quantity calculating module and used for receiving the local oscillator signals and the disturbance quantity transmitted by the disturbance quantity calculating module and updating different modes of the local oscillator signals according to the type of the disturbance quantity.
2. The magnetic resonance B0 field disturbance compensation system of claim 1, wherein: the classifying the disturbance information includes: dividing the disturbance information into real-time calculation and timing calculation; when the packet type of the disturbance information is 5A5A in 16 system, the disturbance information needs to be calculated in real time, and when the packet type of the disturbance information is A5A5 in 16 system, the disturbance information needs to be calculated in timing.
3. The magnetic resonance B0 field disturbance compensation system of claim 1, wherein: the normalization processing comprises the following steps: finding the maximum value in the data packet parameters, dividing each parameter in the data packet parameters by the maximum value, and updating the original parameters by the obtained numerical value.
4. The magnetic resonance B0 field disturbance compensation system of claim 1, wherein: the filtering process is as follows: and filtering the data in the data packet parameters by a high-pass FIR filter.
5. The magnetic resonance B0 field disturbance compensation system of claim 1, wherein: the local oscillator updating module comprises a local oscillator timing compensation module and a local oscillator real-time compensation module:
a timer for circularly timing is arranged in the local oscillator timing compensation module, for the local oscillator signals of which the disturbance information is timing, the timer adds the timing update disturbance quantity to the local oscillator when the time of the timer is over, then the next timing period is entered again, the timer adds the timing update disturbance quantity to the local oscillator in a circulating manner, and a new local oscillator is transmitted to the local oscillator real-time compensation module;
and the local oscillator real-time compensation module is used for adding the new local oscillator sent by the local oscillator timing compensation module and the disturbance quantity updated in real time to obtain a compensated local oscillator.
6. A magnetic resonance B0 field disturbance compensation method according to any one of claims 1-5, characterized by comprising the following steps:
step 1: receiving disturbance information of a magnetic resonance B0 field, classifying the disturbance information, and performing primary processing on classified data; the primary processing comprises normalization processing and filtering processing of the disturbance information data packet;
step 2: according to the disturbance information and disturbance control parameters of different types after the preliminary processing, corresponding disturbance quantity calculation is carried out to obtain the calculated disturbance quantity;
the disturbance amount calculation process includes: firstly, judging the type of disturbance information, and if the type of the disturbance information needs timing calculation, the formula of the timing calculation is as follows:
where f (x) represents the amount of disturbance in the present calculation, k represents the index of the parameter, and akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter data packet, wherein f (x-1) represents the disturbance quantity obtained by the last timing calculation, the last disturbance quantity is 0 in the first calculation, and the calculation formula of the real-time calculation is as follows:
wherein f (x) represents the disturbance amount calculated this time, k represents the index of the parameter, akRepresenting the kth parameter in the packet of perturbation information, bkRepresenting the kth parameter in the disturbance control parameter data packet;
and step 3: and updating different modes of the local oscillation signal according to the type of the disturbance quantity.
7. The magnetic resonance B0 field perturbation compensation method of claim 6, wherein: the classifying the disturbance information includes: dividing the disturbance information into real-time calculation and timing calculation; when the packet type of the disturbance information is 5A5A in 16 system, the disturbance information needs to be calculated in real time, and when the packet type of the disturbance information is A5A5 in 16 system, the disturbance information needs to be calculated in timing.
8. The magnetic resonance B0 field perturbation compensation method of claim 6, wherein: the updating of the local oscillator signal in different modes according to the type of the disturbance quantity comprises:
firstly, for local oscillation signals with disturbance information being timing, adding a timing update disturbance quantity to the local oscillation through a timer when the time is over, then entering the next timing period again, and circularly adding the timing update disturbance quantity to the local oscillation to generate a new local oscillation;
and then adding the new local oscillator and the disturbance quantity updated in real time to obtain a compensated local oscillator.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910421708.7A CN110133559B (en) | 2019-05-21 | 2019-05-21 | Magnetic resonance B0 field disturbance compensation system and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910421708.7A CN110133559B (en) | 2019-05-21 | 2019-05-21 | Magnetic resonance B0 field disturbance compensation system and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110133559A CN110133559A (en) | 2019-08-16 |
CN110133559B true CN110133559B (en) | 2021-05-11 |
Family
ID=67571995
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910421708.7A Active CN110133559B (en) | 2019-05-21 | 2019-05-21 | Magnetic resonance B0 field disturbance compensation system and method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110133559B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1071753A (en) * | 1991-10-16 | 1993-05-05 | 衡阳工学院 | The periodic sampling penalty method of level gauging vibration interference |
CN1328645A (en) * | 1998-09-30 | 2001-12-26 | 威尔弗里德·伯格曼 | Probe head for NMR spectrometer |
CN1367378A (en) * | 2000-10-02 | 2002-09-04 | Ge医疗系统环球技术有限公司 | Magnetic resonance imaging device with field perturbance compensation |
CN1493871A (en) * | 2002-05-15 | 2004-05-05 | 西门子公司 | Vibration compensating method of nuclear spin fault radiography |
CN103235509A (en) * | 2013-03-29 | 2013-08-07 | 北京控制工程研究所 | Rotating member disturbance compensation method based on momentum wheel |
CN104898172A (en) * | 2015-05-19 | 2015-09-09 | 吉林大学 | Cross-correlation-based nuclear magnetic resonance full wave signal noise filtering method |
CN106646637A (en) * | 2016-12-27 | 2017-05-10 | 吉林大学 | Method for removing peak noise in nuclear magnetism signal |
CN206258596U (en) * | 2016-12-09 | 2017-06-16 | 吉林大学 | A kind of ground nuclear magnetic resonance sampling apparatus |
CN107957566A (en) * | 2017-11-17 | 2018-04-24 | 吉林大学 | Magnetic resonance depth measurement method for extracting signal based on frequency selection singular spectrum analysis |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3728167B2 (en) * | 2000-02-10 | 2005-12-21 | 株式会社日立メディコ | Magnetic resonance imaging system |
-
2019
- 2019-05-21 CN CN201910421708.7A patent/CN110133559B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1071753A (en) * | 1991-10-16 | 1993-05-05 | 衡阳工学院 | The periodic sampling penalty method of level gauging vibration interference |
CN1328645A (en) * | 1998-09-30 | 2001-12-26 | 威尔弗里德·伯格曼 | Probe head for NMR spectrometer |
CN1367378A (en) * | 2000-10-02 | 2002-09-04 | Ge医疗系统环球技术有限公司 | Magnetic resonance imaging device with field perturbance compensation |
CN1493871A (en) * | 2002-05-15 | 2004-05-05 | 西门子公司 | Vibration compensating method of nuclear spin fault radiography |
CN103235509A (en) * | 2013-03-29 | 2013-08-07 | 北京控制工程研究所 | Rotating member disturbance compensation method based on momentum wheel |
CN104898172A (en) * | 2015-05-19 | 2015-09-09 | 吉林大学 | Cross-correlation-based nuclear magnetic resonance full wave signal noise filtering method |
CN206258596U (en) * | 2016-12-09 | 2017-06-16 | 吉林大学 | A kind of ground nuclear magnetic resonance sampling apparatus |
CN106646637A (en) * | 2016-12-27 | 2017-05-10 | 吉林大学 | Method for removing peak noise in nuclear magnetism signal |
CN107957566A (en) * | 2017-11-17 | 2018-04-24 | 吉林大学 | Magnetic resonance depth measurement method for extracting signal based on frequency selection singular spectrum analysis |
Also Published As
Publication number | Publication date |
---|---|
CN110133559A (en) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9570062B2 (en) | Method and device for self-adaptively eliminating noises | |
US3824501A (en) | Automatic cable equalizer | |
CN103873411B (en) | Method and device for maximum likelihood frequency offset estimation based on joint pilot frequency | |
CN110133559B (en) | Magnetic resonance B0 field disturbance compensation system and method | |
CN116405823B (en) | Intelligent audio denoising enhancement method for bone conduction earphone | |
US20190154778A1 (en) | Apparatus and Method for Analyzing Spectrum | |
CN114995638A (en) | Tactile signal generation method and device, readable medium and electronic equipment | |
WO2023273363A1 (en) | Clock frequency offset tracking method based on weighted observation fusion and timestamp-free interaction | |
CN110164461A (en) | Audio signal processing method, device, electronic equipment and storage medium | |
CN108872900A (en) | A kind of time division multiplexing nuclear magnetic resonance gradient number pre-emphasis apparatus | |
CN115905817B (en) | Multi-channel non-uniform field signal noise reduction method, device and storage medium | |
CN109039968B (en) | Data receiving method and device based on narrow-band Internet of things terminal | |
JP2001078987A (en) | Filtering method for oscillatory b0 behavior of nuclear magnetic resonance apparatus | |
US20140067739A1 (en) | Reduction or elimination of training for adaptive filters and neural networks through look-up table | |
Vucic et al. | A class of systems with symmetric impulse response | |
CN111540372B (en) | Method and device for noise reduction processing of multi-microphone array | |
Varkonyi-Koczy | Fast anytime fuzzy Fourier estimation of multisine signals | |
CN103258106A (en) | Stability analysis system and method for data processing equipment | |
CN115118296B (en) | Real-time error compensation method, system, terminal and medium for multichannel receiver | |
CN105429914A (en) | Channel correction method and apparatus based on equalizer compensation | |
CN108347394A (en) | Multi-carrier communications systems and its channel estimation method | |
CN111770413B (en) | Multi-sound-source sound mixing method and device and storage medium | |
CN101304480A (en) | Method and system for eliminating ghost of television signal based on wavelet preprocessing GCR | |
CN114567390B (en) | Method and device for eliminating cross-frequency-band phase jump points of vector network analyzer receiver | |
CN101888364B (en) | Carrier synchronization processing method and device |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211020 Address after: 100176 West, 1st floor, building 5, No.8, Hongda North Road, Beijing Economic and Technological Development Zone, Daxing District, Beijing Patentee after: Broad Teng (Beijing) medical instruments Co.,Ltd. Address before: 117004 No.9 Yaodu street, economic and Technological Development Zone, Benxi City, Liaoning Province Patentee before: LIAONING KAMPO MEDICAL SYSTEM Co.,Ltd. |