CN111460934A - Method and system for calculating center frequency of superconducting high-field magnetic resonance - Google Patents
Method and system for calculating center frequency of superconducting high-field magnetic resonance Download PDFInfo
- Publication number
- CN111460934A CN111460934A CN202010191162.3A CN202010191162A CN111460934A CN 111460934 A CN111460934 A CN 111460934A CN 202010191162 A CN202010191162 A CN 202010191162A CN 111460934 A CN111460934 A CN 111460934A
- Authority
- CN
- China
- Prior art keywords
- water
- fat
- frequency
- acquiring
- information
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/14—Classification; Matching by matching peak patterns
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
Abstract
The embodiment of the invention discloses a method and a system for calculating the center frequency of superconducting high-field magnetic resonance, wherein the method comprises the following steps: acquiring an original FID signal, and after data cleaning is carried out on the original FID signal, acquiring waveform data corresponding to human body water fat frequency; analyzing and counting the waveform data to obtain structural data of wave crest information; judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks; acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value; and acquiring the center frequency according to the water frequency position and the fat frequency position. The embodiment of the invention can automatically and quickly find out the accurate water frequency position, and can analyze and calculate the optimal water frequency position when the interference signals such as signals generated by silica gel filler, a plurality of fat frequencies exist and different channels have differences.
Description
Technical Field
The invention relates to the technical field of magnetic resonance imaging, in particular to a method and a system for calculating the center frequency of superconducting high-field magnetic resonance.
Background
In a superconducting magnetic resonance system, the center frequency is the basis for determining the physical center of a signal and the center of gradient encoding, and plays a very important role in accurate acquisition and correct reconstruction of the signal, so that the center frequency needs to be accurately searched at the signal acquisition position. And the water frequency and the fat frequency are calculated together, and an accurate central frequency position is judged. The magnetic resonance imaging system mainly collects signals of H protons in a human body, and under the condition of high-field magnetic resonance (more than or equal to 3Tesla), various complex water and fat inclusion conditions exist.
The existing method for finding the center frequency is to calculate the frequency difference between water and fat, for example, the frequency difference is about 3.4ppm under 3T condition, or add a manual judgment to the designated water frequency position to obtain the center frequency value of the water signal.
The existing center frequency algorithm searching method cannot improve the accuracy when processing complex frequency distribution or interference conditions. For example, under the 3Tesla condition, the fat frequency signal greatly exceeds the water frequency signal, and is mixed with an interference signal, and under the condition that a plurality of receiving channel signal differences exist, the accurate search of the frequency position of the water signal becomes very complicated and uncertain. The existing method mostly depends on the integration of automatic calculation and manual judgment to search and confirm, and the frequency position of the water outlet signal cannot be accurately judged by means of human eye judgment and the traditional water and fat fixed frequency calculation method.
The prior art is therefore still subject to further development.
Disclosure of Invention
In view of the above technical problems, embodiments of the present invention provide a method and a system for calculating a center frequency of a superconducting high-field magnetic resonance, which can solve the technical problem in the prior art that the frequency position of an effluent signal cannot be accurately determined by human eye determination and a conventional water-fat fixed frequency calculation method.
The first aspect of the embodiments of the present invention provides a method for calculating a center frequency of a superconducting high-field magnetic resonance, including:
acquiring an original FID signal, and after data cleaning is carried out on the original FID signal, acquiring waveform data corresponding to human body water fat frequency;
analyzing and counting the waveform data to obtain structural data of wave crest information;
judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks;
acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value;
and acquiring the center frequency according to the water frequency position and the fat frequency position.
Optionally, the obtaining an original FID signal, and after performing data cleaning on the original FID signal, obtaining waveform data corresponding to the human body water fat frequency includes:
acquiring an original FID signal, and performing smoothing processing on the FID signal to generate a smooth FID signal;
and performing band-pass filtering processing on the smooth FID signal to acquire the waveform data of the water and fat frequency of the human body data.
Optionally, the analyzing and counting the waveform data to obtain structural data of peak information includes:
obtaining amplitude of a peak position, a corresponding inflection point of a wave trough and full-height half-width information in waveform data;
and analyzing and counting the amplitude of the peak position, the inflection point of the corresponding wave trough and the full-height half-width to construct structural data of peak information.
Optionally, the determining the validity of the peak information according to the structural data of the peak information and counting the peak information of the valid peaks includes:
acquiring a water and fat frequency characteristic, and judging whether the wave crest conforms to the water and fat frequency characteristic according to the structural data of the wave crest information;
if the structural data of the current wave crest information is matched with the currently measured human body part characteristic, judging that the wave crest is effective, and counting the number of the effective wave crests; if the structural data of the current peak information is not matched with the currently measured human body part characteristics, judging that the peak is invalid;
and if the water and fat frequency characteristics are not met, judging that the wave crest is invalid.
Optionally, the analyzing and counting the waveform data to obtain structural data of peak information includes:
judging possible peak positions by using a sliding window, judging whether the peaks are correct or not through a signal amplitude threshold value, and acquiring all correct peaks;
respectively acquiring troughs on two sides of each peak by taking each peak as a starting point, and judging the waveform to determine whether the troughs are normal or not to acquire normal troughs;
and performing waveform statistics on the wave crests and the wave troughs, eliminating abnormal wave crests according to the positions and the signal sizes, and keeping normal wave crest information.
A second aspect of an embodiment of the present invention provides a center frequency calculation system of superconducting high-field magnetic resonance, including: a memory, a processor and a computer program stored on the memory and executable on the processor, the computer program when executed by the processor implementing the steps of:
acquiring an original FID signal, and after data cleaning is carried out on the original FID signal, acquiring waveform data corresponding to human body water fat frequency;
analyzing and counting the waveform data to obtain structural data of wave crest information;
judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks;
acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value;
and acquiring the center frequency according to the water frequency position and the fat frequency position.
Optionally, the computer program when executed by the processor further implements the steps of:
acquiring an original FID signal, and performing smoothing processing on the FID signal to generate a smooth FID signal;
and performing band-pass filtering processing on the smooth FID signal to acquire the waveform data of the water and fat frequency of the human body data.
Optionally, the computer program when executed by the processor further implements the steps of:
obtaining amplitude of a peak position, a corresponding inflection point of a wave trough and full-height half-width information in waveform data;
and analyzing and counting the amplitude of the peak position, the inflection point of the corresponding wave trough and the full-height half-width to construct structural data of peak information.
Optionally, the computer program when executed by the processor further implements the steps of:
acquiring a water and fat frequency characteristic, and judging whether the wave crest conforms to the water and fat frequency characteristic according to the structural data of the wave crest information;
if the structural data of the current wave crest information is matched with the currently measured human body part characteristic, judging that the wave crest is effective, and counting the number of the effective wave crests; if the structural data of the current peak information is not matched with the currently measured human body part characteristics, judging that the peak is invalid;
and if the water and fat frequency characteristics are not met, judging that the wave crest is invalid.
A third aspect of embodiments of the present invention provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer-executable instructions, and when executed by one or more processors, the computer-executable instructions cause the one or more processors to perform the above-mentioned method for calculating the center frequency of the superconducting high-field magnetic resonance.
According to the technical scheme provided by the embodiment of the invention, an original FID signal is obtained, and waveform data corresponding to the human body water fat frequency is obtained after the data of the original FID signal is cleaned; analyzing and counting the waveform data to obtain structural data of wave crest information; judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks; acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value; and acquiring the center frequency according to the water frequency position and the fat frequency position. Therefore, compared with the prior art, the embodiment of the invention can automatically and quickly find the accurate water frequency position, and can analyze and calculate the optimal water frequency position when the interference signals such as signals generated by silica gel filler and a plurality of fat frequencies exist and different channels have differences.
Drawings
FIG. 1 is a schematic flow chart illustrating an embodiment of a method for calculating a center frequency of a superconducting high field magnetic resonance according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a hardware configuration of another embodiment of a center frequency calculation system for superconducting high-field magnetic resonance according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Aiming at the problems that when a human body is scanned by high-field magnetic resonance, a plurality of signal peak values of human body signals (including a plurality of fat and water signals) are distributed in different frequency ranges and cannot be completely consistent in performance at different positions and different receiving channels, and interference signals (generated by silica gel fillers in the human body and the like) are generated, the signals need to be judged in real time and the frequency positions corresponding to the accurate water signals need to be found.
The following detailed description of embodiments of the invention refers to the accompanying drawings.
Referring to fig. 1, fig. 1 is a schematic flow chart illustrating an embodiment of a method for calculating a center frequency of a superconducting high-field magnetic resonance according to an embodiment of the present invention. As shown in fig. 1, includes:
s100, acquiring an original FID signal, and acquiring waveform data corresponding to the human body water-fat frequency after data cleaning is carried out on the original FID signal;
s200, analyzing and counting the waveform data to obtain structural data of wave crest information;
s300, judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the valid peaks;
s400, acquiring a water and fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to wave crest information of the effective wave crest and the water and fat deviation value;
and S500, acquiring the center frequency according to the water frequency position and the fat frequency position.
Specifically, cleaning original data, smoothing FID signals, and performing band-pass filtering to obtain waveform data of water and fat frequency; free Induction Decay (FID) is the simplest signal form in magnetic resonance and imaging.
Analyzing and counting peak information, namely analyzing and counting information such as amplitude of a peak position, a corresponding inflection point of a trough, full height and half width and the like to obtain structural data of the peak information;
judging and counting the effectiveness of the wave crests, judging whether the wave crests conform to the water-fat frequency characteristic, and judging the stability according to the characteristics of different parts of a human body; characteristics of different parts of a human body: the method includes that (1) the fat frequency number and the signal size are different at different parts, and (2) different types of radio frequency systems (transmitting and receiving coils) have different signals under the condition of a plurality of channels, so that the signals of a human body need to be acquired under the condition that a certain part and the radio frequency system are clear, the effective range of band-pass filtering is confirmed through the signal distribution condition, and complete parameter configuration is formed in advance to be used as a known variable when the central frequency is automatically searched subsequently.
Confirming the water and fat frequency position: according to the actual effective peak informationAdding a water and fat deviation value, judging the positions of an optimized water frequency and a corresponding fat frequency, ppm: one millionth, representing parameters such as concentration or uniformity and the like, obtaining a difference value of 3.4ppm according to a superconducting magnetic resonance experiment, wherein H protons in fat and water molecules have a difference due to precession frequency, and under the field intensity of 3Tesla, a default resonance frequency is w ═ r w according to a resonance frequency calculation formula (Larmor L amor frequency formula)) of magnetic resonance, w is an angular frequency, r gyromagnetic ratio constant 42.57748MHz/T, B is a magnetic field intensity, and when 5 resonance frequency is 127.732MHz, the frequency corresponding to 3.4ppm is △ fwf＝127.732MHz*3.4/1000000≈434Hz。
Further, acquiring an original FID signal, and after data cleaning is performed on the original FID signal, obtaining waveform data corresponding to the human body water fat frequency, wherein the waveform data comprises:
acquiring an original FID signal, and performing smoothing processing on the FID signal to generate a smooth FID signal;
and performing band-pass filtering processing on the smooth FID signal to acquire the waveform data of the water and fat frequency of the human body data.
Specifically, a multi-channel FID frequency domain signal obtained by magnetic resonance scanning of a human body is acquired. Generally, 10-channel FID frequency domain signals are used.
S1＝F1bp(F0smooth(S0(x))) (1)
In formula 1, S0(x) is a frequency domain signal, F0smoothFor smoothing filter processing function, F1bpAs a function of the bandpass filtering process, S1 is a frequency domain signal after removing the interference signal. The problem of interfering signals is solved by smoothing and band-pass filtering.
10 channels is a common example of a multi-channel application. The smoothing filter function adopts Gaussian smoothing filter function gaussian, and is mainly used for removing burrs under the condition of not influencing the signal variation trend; the band-pass filtering function takes the central position as a symmetrical point and is a finite impulse product response function type, the number of signal points in a left effective range and a right effective range is selected, tissue signals are reserved, and a background noise signal area is removed.
Further, analyzing and counting the waveform data to obtain structural data of the peak information, including:
obtaining amplitude of a peak position, a corresponding inflection point of a wave trough and full-height half-width information in waveform data;
and analyzing and counting the amplitude of the peak position, the inflection point of the corresponding wave trough and the full-height half-width to construct structural data of peak information.
Specifically, for the distribution situation of water and fat frequency peaks, all peak value information including peak position, amplitude, full-height half-width starting point and end point information is calculated first, and a peak distribution data structure is established. The water and fat frequency peak calculation functions are the same; the waveform calculation function is of the form: FWHM: full height and half width. And the distance between two points passed by the left side and the right side when the maximum height on the batch full wave chart is half is used for representing the layer thickness of the wave signal.
Further, analyzing and counting the waveform data to obtain structural data of the peak information, including:
judging possible peak positions by using a sliding window, judging whether the peaks are correct or not through a signal amplitude threshold value, and acquiring all correct peaks;
respectively acquiring troughs on two sides of each peak by taking each peak as a starting point, and judging the waveform to determine whether the troughs are normal or not to acquire normal troughs;
and performing waveform statistics on the wave crests and the wave troughs, eliminating abnormal wave crests according to the positions and the signal sizes, and keeping normal wave crest information.
During specific implementation, a sliding window is used for judging possible peak positions, whether the peaks are correct or not is judged through a signal amplitude threshold value, and all correct peaks are found; (2) respectively finding the wave troughs on two sides of each wave peak by taking each wave peak as a starting point, and judging whether the wave troughs are normal or not by waveform; and analyzing the positions of the troughs under abnormal conditions such as possible overlapping of the troughs and higher troughs when a plurality of peaks are analyzed according to FWHM calculation, (3) carrying out waveform statistics on the peaks and the troughs, eliminating abnormal peaks according to the positions and the signal sizes, and keeping normal peak information)
Pclass＝Wdetect(S1) (2)
In equation 2, the function W is the waveform calculationdetectThe structural data P containing the information listed above is obtainedclass(ii) a Typically, multi-channel signal reception is used, and thus the configuration data corresponds to data of a plurality of channels.
Further, judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the valid peaks, including:
acquiring a water and fat frequency characteristic, and judging whether the wave crest conforms to the water and fat frequency characteristic according to the structural data of the wave crest information;
if the structural data of the current wave crest information is matched with the currently measured human body part characteristic, judging that the wave crest is effective, and counting the number of the effective wave crests; if the structural data of the current peak information is not matched with the currently measured human body part characteristics, judging that the peak is invalid;
and if the water and fat frequency characteristics are not met, judging that the wave crest is invalid.
During specific implementation, the water peak position is optimally judged according to structural data of a peak, the water frequency is 450Hz faster than the adjacent fat frequency when 3Tesla is considered, the frequency difference is relatively fixed, but when different channels and other fat frequencies are relatively unfixed, the stability value differentiation configuration is carried out when different human body parts are referred, and the optimal judgment is required:
Pwater＝F2opt(Pclass) (3)
and obtaining water frequency position information in optimization judgment, wherein formula 3 is a search result of the central frequency (the best matching water frequency position according with two characteristics is calculated by using the characteristics that formula 1 only has one water peak and a plurality of fat peaks and the fixed frequency difference (the water frequency is 450Hz higher than the adjacent fat frequency in 3T magnetic resonance) of formula 2 between water and fat).
According to the method and the device, the acquisition of the FID signals is based on, the acquisition signals are filtered and cleaned, and the frequency domain signals are analyzed and judged, so that a new center frequency searching algorithm is formed, and the purpose of accurately searching the center frequency is achieved. The embodiment of the invention can automatically and quickly analyze and search to obtain an accurate result, does not depend on manual auxiliary positioning and accurate diagnosis, is not influenced by external background interference, human filler interference signals and the like, and greatly improves the accuracy and the stability.
The above description is about the method for calculating the center frequency of the superconducting high-field magnetic resonance in the embodiment of the present invention, and the following description is about the system for calculating the center frequency of the superconducting high-field magnetic resonance in the embodiment of the present invention, please refer to fig. 2, fig. 2 is a schematic hardware structure diagram of another embodiment of the system for calculating the center frequency of the superconducting high-field magnetic resonance in the embodiment of the present invention, and as shown in fig. 2, the system 10 includes: a memory 101, a processor 102 and a computer program stored on the memory and executable on the processor, the computer program realizing the following steps when executed by the processor 101:
acquiring an original FID signal, and after data cleaning is carried out on the original FID signal, acquiring waveform data corresponding to human body water fat frequency;
analyzing and counting the waveform data to obtain structural data of wave crest information;
judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks;
acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value;
and acquiring the center frequency according to the water frequency position and the fat frequency position.
The specific implementation steps are the same as those of the method embodiments, and are not described herein again.
Optionally, the computer program when executed by the processor 101 further implements the steps of:
acquiring an original FID signal, and performing smoothing processing on the FID signal to generate a smooth FID signal;
and performing band-pass filtering processing on the smooth FID signal to acquire the waveform data of the water and fat frequency of the human body data.
The specific implementation steps are the same as those of the method embodiments, and are not described herein again.
Optionally, the computer program when executed by the processor 101 further implements the steps of:
obtaining amplitude of a peak position, a corresponding inflection point of a wave trough and full-height half-width information in waveform data;
and analyzing and counting the amplitude of the peak position, the inflection point of the corresponding wave trough and the full-height half-width to construct structural data of peak information.
The specific implementation steps are the same as those of the method embodiments, and are not described herein again.
Optionally, the computer program when executed by the processor 101 further implements the steps of:
acquiring a water and fat frequency characteristic, and judging whether the wave crest conforms to the water and fat frequency characteristic according to the structural data of the wave crest information;
if the structural data of the current wave crest information is matched with the currently measured human body part characteristic, judging that the wave crest is effective, and counting the number of the effective wave crests; if the structural data of the current peak information is not matched with the currently measured human body part characteristics, judging that the peak is invalid;
and if the water and fat frequency characteristics are not met, judging that the wave crest is invalid.
The specific implementation steps are the same as those of the method embodiments, and are not described herein again.
Optionally, the computer program when executed by the processor 101 further implements the steps of:
judging possible peak positions by using a sliding window, judging whether the peaks are correct or not through a signal amplitude threshold value, and acquiring all correct peaks;
respectively acquiring troughs on two sides of each peak by taking each peak as a starting point, and judging the waveform to determine whether the troughs are normal or not to acquire normal troughs;
and performing waveform statistics on the wave crests and the wave troughs, eliminating abnormal wave crests according to the positions and the signal sizes, and keeping normal wave crest information.
The specific implementation steps are the same as those of the method embodiments, and are not described herein again.
Embodiments of the present invention provide a non-transitory computer-readable storage medium storing computer-executable instructions for execution by one or more processors, for example, to perform method steps S100-S500 of fig. 1 described above.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for calculating the center frequency of superconducting high-field magnetic resonance is characterized by comprising the following steps:
acquiring an original FID signal, and after data cleaning is carried out on the original FID signal, acquiring waveform data corresponding to human body water fat frequency;
analyzing and counting the waveform data to obtain structural data of wave crest information;
judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks;
acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value;
and acquiring the center frequency according to the water frequency position and the fat frequency position.
2. The method for calculating the center frequency of the superconducting high-field magnetic resonance according to claim 1, wherein the obtaining of the original FID signal and the data cleaning of the original FID signal to obtain the waveform data corresponding to the human body water and fat frequency comprises:
acquiring an original FID signal, and performing smoothing processing on the FID signal to generate a smooth FID signal;
and performing band-pass filtering processing on the smooth FID signal to acquire the waveform data of the water and fat frequency of the human body data.
3. The method of claim 2, wherein the analyzing and counting the waveform data to obtain the structural data of the peak information comprises:
obtaining amplitude of a peak position, a corresponding inflection point of a wave trough and full-height half-width information in waveform data;
and analyzing and counting the amplitude of the peak position, the inflection point of the corresponding wave trough and the full-height half-width to construct structural data of peak information.
4. The method according to claim 3, wherein the determining the validity of the peak information according to the structural data of the peak information and counting the peak information of the valid peaks comprises:
acquiring a water and fat frequency characteristic, and judging whether the wave crest conforms to the water and fat frequency characteristic according to the structural data of the wave crest information;
if the structural data of the current wave crest information is matched with the currently measured human body part characteristic, judging that the wave crest is effective, and counting the number of the effective wave crests; if the structural data of the current peak information is not matched with the currently measured human body part characteristics, judging that the peak is invalid;
and if the water and fat frequency characteristics are not met, judging that the wave crest is invalid.
5. The method of claim 4, wherein the analyzing and counting the waveform data to obtain the structural data of the peak information comprises:
judging possible peak positions by using a sliding window, judging whether the peaks are correct or not through a signal amplitude threshold value, and acquiring all correct peaks;
respectively acquiring troughs on two sides of each peak by taking each peak as a starting point, and judging the waveform to determine whether the troughs are normal or not to acquire normal troughs;
and performing waveform statistics on the wave crests and the wave troughs, eliminating abnormal wave crests according to the positions and the signal sizes, and keeping normal wave crest information.
6. A center frequency calculation system for superconducting high field magnetic resonance, the system comprising: a memory, a processor and a computer program stored on the memory and executable on the processor, the computer program when executed by the processor implementing the steps of:
acquiring an original FID signal, and after data cleaning is carried out on the original FID signal, acquiring waveform data corresponding to human body water fat frequency;
analyzing and counting the waveform data to obtain structural data of wave crest information;
judging the validity of the peak information according to the structural data of the peak information, and counting the peak information of the effective peaks;
acquiring a water-fat deviation value, and acquiring an optimized water frequency position and a corresponding fat frequency position according to the wave crest information of the effective wave crest and the water-fat deviation value;
and acquiring the center frequency according to the water frequency position and the fat frequency position.
7. The system of claim 6, wherein the computer program when executed by the processor further performs the steps of:
acquiring an original FID signal, and performing smoothing processing on the FID signal to generate a smooth FID signal;
and performing band-pass filtering processing on the smooth FID signal to acquire the waveform data of the water and fat frequency of the human body data.
8. The system of claim 7, wherein the computer program when executed by the processor further performs the steps of:
obtaining amplitude of a peak position, a corresponding inflection point of a wave trough and full-height half-width information in waveform data;
and analyzing and counting the amplitude of the peak position, the inflection point of the corresponding wave trough and the full-height half-width to construct structural data of peak information.
9. The system of claim 8, wherein the computer program when executed by the processor further performs the steps of:
acquiring a water and fat frequency characteristic, and judging whether the wave crest conforms to the water and fat frequency characteristic according to the structural data of the wave crest information;
if the structural data of the current wave crest information is matched with the currently measured human body part characteristic, judging that the wave crest is effective, and counting the number of the effective wave crests; if the structural data of the current peak information is not matched with the currently measured human body part characteristics, judging that the peak is invalid;
and if the water and fat frequency characteristics are not met, judging that the wave crest is invalid.
10. A non-transitory computer-readable storage medium storing computer-executable instructions that, when executed by one or more processors, cause the one or more processors to perform the method of calculating a center frequency of superconducting high field magnetic resonance of any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010191162.3A CN111460934B (en) | 2020-03-18 | 2020-03-18 | Center frequency calculation method and system for superconducting high-field magnetic resonance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010191162.3A CN111460934B (en) | 2020-03-18 | 2020-03-18 | Center frequency calculation method and system for superconducting high-field magnetic resonance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111460934A true CN111460934A (en) | 2020-07-28 |
| CN111460934B CN111460934B (en) | 2023-09-29 |
Family
ID=71682099
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010191162.3A Active CN111460934B (en) | 2020-03-18 | 2020-03-18 | Center frequency calculation method and system for superconducting high-field magnetic resonance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111460934B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112069161A (en) * | 2020-09-01 | 2020-12-11 | 上海佰贝科技发展股份有限公司 | Data cleaning method and device |
| CN112948747A (en) * | 2021-02-05 | 2021-06-11 | 深圳市贝斯达医疗股份有限公司 | Method, device, terminal equipment and storage medium for searching center frequency |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5557202A (en) * | 1993-04-27 | 1996-09-17 | Kabushiki Kaisha Toshiba | Method and system for magnetic resonance imaging |
| US5677626A (en) * | 1993-04-27 | 1997-10-14 | Kabushiki Kaisha Toshiba | System for magnetic resonance imaging |
| CN1394551A (en) * | 2001-06-28 | 2003-02-05 | Ge医疗系统环球技术有限公司 | Self-rotary exciting method and magnetic resonance imaging system |
| JP2004024902A (en) * | 1994-04-15 | 2004-01-29 | Toshiba Corp | Magnetic resonance imaging apparatus |
| CN104502873A (en) * | 2014-11-21 | 2015-04-08 | 鑫高益医疗设备股份有限公司 | Quick and active shimming method for MRI (Magnetic Resonance Imaging) system |
| CN106952261A (en) * | 2017-04-06 | 2017-07-14 | 深圳安科高技术股份有限公司 | A kind of image processing method and device for suppressing fat signal |
| CN107121653A (en) * | 2017-05-12 | 2017-09-01 | 杭州全景医学影像诊断有限公司 | The method that fat presses down water process magnetic resonance spectrum imaging data is pressed down based on frequency domain |
| CN109752682A (en) * | 2019-01-22 | 2019-05-14 | 奥泰医疗系统有限责任公司 | Centre frequency bearing calibration based on phase diagram |
-
2020
- 2020-03-18 CN CN202010191162.3A patent/CN111460934B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5557202A (en) * | 1993-04-27 | 1996-09-17 | Kabushiki Kaisha Toshiba | Method and system for magnetic resonance imaging |
| US5677626A (en) * | 1993-04-27 | 1997-10-14 | Kabushiki Kaisha Toshiba | System for magnetic resonance imaging |
| JP2004024902A (en) * | 1994-04-15 | 2004-01-29 | Toshiba Corp | Magnetic resonance imaging apparatus |
| CN1394551A (en) * | 2001-06-28 | 2003-02-05 | Ge医疗系统环球技术有限公司 | Self-rotary exciting method and magnetic resonance imaging system |
| CN104502873A (en) * | 2014-11-21 | 2015-04-08 | 鑫高益医疗设备股份有限公司 | Quick and active shimming method for MRI (Magnetic Resonance Imaging) system |
| CN106952261A (en) * | 2017-04-06 | 2017-07-14 | 深圳安科高技术股份有限公司 | A kind of image processing method and device for suppressing fat signal |
| CN107121653A (en) * | 2017-05-12 | 2017-09-01 | 杭州全景医学影像诊断有限公司 | The method that fat presses down water process magnetic resonance spectrum imaging data is pressed down based on frequency domain |
| CN109752682A (en) * | 2019-01-22 | 2019-05-14 | 奥泰医疗系统有限责任公司 | Centre frequency bearing calibration based on phase diagram |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112069161A (en) * | 2020-09-01 | 2020-12-11 | 上海佰贝科技发展股份有限公司 | Data cleaning method and device |
| CN112069161B (en) * | 2020-09-01 | 2023-11-03 | 上海佰贝科技发展股份有限公司 | Data cleaning method and device |
| CN112948747A (en) * | 2021-02-05 | 2021-06-11 | 深圳市贝斯达医疗股份有限公司 | Method, device, terminal equipment and storage medium for searching center frequency |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111460934B (en) | 2023-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111460934B (en) | Center frequency calculation method and system for superconducting high-field magnetic resonance | |
| CN114025379B (en) | Broadband multi-signal detection method, device and equipment | |
| CN112494001B (en) | PPG signal quality evaluation method and device, wearable device and storage medium | |
| CN103592609B (en) | Automated spectral fat saturation | |
| EP4343291A1 (en) | Global phase quadrature demodulation method and apparatus for optical fiber sensing data | |
| CN110313903B (en) | Pulse wave frequency domain characteristic parameter extraction method and device | |
| US20190377050A1 (en) | System and method for improved diffusion-weighted imaging | |
| US10145921B2 (en) | Method and apparatus for automated determination of the resonance frequencies of protons for magnetic resonance examinations | |
| CN104155621B (en) | Method used for accurately measuring static magnetic field B0 distribution | |
| CN116660176A (en) | Fourier spectrum automatic baseline correction method, device and storage medium | |
| US10969449B2 (en) | Local shimming system for magnetic resonance imaging and method thereof | |
| CN113375065B (en) | Method and device for eliminating trend signal in pipeline leakage monitoring | |
| CN111724386B (en) | Data preprocessing method and system for multi-frequency imaging and multi-frequency imaging system | |
| CN113269798A (en) | Medical image segmentation method and device, computer equipment and storage medium thereof | |
| WO2021203705A1 (en) | Gradient field control method and apparatus, and magnetic resonance imaging device and medium | |
| CN111239657B (en) | Spectrogram phase correction method, device and equipment | |
| US11559215B2 (en) | Method for measuring high-accuracy realtime heart rate based on continuous-wave Doppler radar and radar system therefor | |
| CN114847971A (en) | Deterministic sine intervention synchronous compression transformation spindle wave extraction method and application thereof | |
| CN109446572B (en) | Quantum measurement and control data processing method and device | |
| CN112782764B (en) | Linear interference attenuation method and system | |
| CN110742687A (en) | Method for non-invasive prediction of growth rate of pituitary adenoma tumor | |
| CN112668134B (en) | Method and device for determining harmonic impedance of traction network and terminal equipment | |
| US20230225630A1 (en) | Method for attenuating the noise in images resulting from multiple acquisitions by magnetic resonance imaging | |
| CN112948747B (en) | Method, device, terminal equipment and storage medium for searching center frequency | |
| CN116008220A (en) | Terahertz detection method, terahertz detection device and computer-readable 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 | ||
| CB02 | Change of applicant information |
Address after: 518000 workshop 1a1901, building a, beisida medical equipment building, 28 Nantong Avenue, Baolong community, Baolong street, Longgang District, Shenzhen City, Guangdong Province Applicant after: SHENZHEN BASDA MEDICAL APPARATUS Co.,Ltd. Address before: 518000, A1402, A1403, No. three, Tian Tian Digital Innovation Park, Longgang center, Longgang District, Shenzhen, Guangdong Applicant before: SHENZHEN BASDA MEDICAL APPARATUS Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |