CN100346171C - Receiving method of nuclear magnetic resonance imaging signal - Google Patents
Receiving method of nuclear magnetic resonance imaging signal Download PDFInfo
- Publication number
- CN100346171C CN100346171C CNB2004100531543A CN200410053154A CN100346171C CN 100346171 C CN100346171 C CN 100346171C CN B2004100531543 A CNB2004100531543 A CN B2004100531543A CN 200410053154 A CN200410053154 A CN 200410053154A CN 100346171 C CN100346171 C CN 100346171C
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- signal
- receiver
- nmr signal
- 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.)
- Expired - Fee Related
Links
Images
Abstract
The present invention relates to a nuclear magnetic resonance imaging method and a nuclear magnetic resonance imaging instrument class, more specifically a receiving method of a nuclear magnetic resonance imaging signal. In the receiving method, a counter which is used for measuring time of errors each time is arranged on a digital receiver. The counter transfers the error time obtained through measurement to a host computer. The linear phase of a spectrum line of the nuclear magnetic resonance imaging signal is corrected in the spectrum width of the receiver. The present invention has the advantages that under the premise that hardware is changed to the minimal extent, the phase jitter of the nuclear magnetic resonance receiver, which is caused by the digital receiver is fundamentally eliminated. The method is very compact and is easily enforced. Simultaneously, because rapid Fourier conversion and inverse Fourier conversion are used, the calculation burden of the host computer is greatly relieved.
Description
Technical field
The present invention relates to magnetic resonance imaging method employing and NMR imaging instrument device class, more particularly relate to a kind of method of reseptance of Magnetic resonance imaging signal.
Background technology
Along with the development of magnetic resonance digital quadrature demodulation method, the application of digital receiver in imaging system more and more widely.Fig. 7 is typical nuclear magnetic resonance digital receiver system synoptic diagram.By the NMR signal that collects of probe (1), after low-noise preamplifier (2) amplifies, enter frequency mixer (3) and carry out detection, carry out digitized sampling through entering analog to digital converter ADC (5) after intermediate-freuqncy signal conditioner (4) amplification and the filtering high fdrequency component.According to being Qwest's sampling thheorem, the sample frequency fs of ADC is usually than the resonant frequency f of magnetic resonance signal
FidAt least the above undistorted collection that could realize of high twice to magnetic resonance signal.If in order further to reduce the quantizing noise of ADC, will realize the over-sampling to NMR signal, the sample frequency of ADC is often high tens to tens times than intermediate frequency.For now, from the difficult real-time processing of the high-speed digital signal of ADC output, there is the people once to solve this problem with the little processing concurrent working of a plurality of high-performance with microprocessor realization digital signal.Obviously, this can increase the complexity and the cost of system to a great extent.More easy way is to adopt special IC (ASIC) to realize.As Fig. 7, ADC sampling gained data are given digital down converter (6), wherein (7A), (7B) be digital quadrature mixing device with 90 degree fixed skew, magnetic resonance signal is down-converted near the zero-frequency by intermediate frequency through it, enter cascade comb filter (CIC) then (8A), (8B) carry out the first order of signal is extracted and filtering, obtain the slower code stream of speed, finite impulse response decimation filter (DEC1 able to programme, DEC2) (9A), (9B) then magnetic resonance signal is carried out some grades of filtering and extraction again, (I Q) deposits main frame (10) in the complex signal that obtains the nuclear magnetic resonance base band at last.
Use the magnetic resonance spectrogram that digital receiver can obtain not having the zero-frequency peak, the mirror image peak disturbs, owing to used digital filter in digital receiver, linear phase shift and out-of-band noise that it also has in the extraordinary passband suppress ability in addition.As can be seen from Figure 7, the digital signal (bandwidth is generally tens megahertzes) of high-speed ADC sampling output, after handling by digital receiver, at a slow speed the narrow-band digital signal (bandwidth is generally tens KHz) often of output.Fig. 2 is the sequential chart of digital receiver output signal, and wherein, A is the nuclear magnetic resonance complex digital signal of digital receiver output.A ' is the digital receiver data frame signal, and each signal has been represented the data sampling point one time.B is the sampling enabling signal, and it is risen along beginning sampling thereon by host computer control.C is actual nuclear magnetic resonance simulating signal.D represents N-1 sampled point of digital receiver and the position on actual NMR signal thereof.D ' then represents N sampled point of digital receiver and the position on actual NMR signal thereof.It is E that sampling is engraved in position corresponding on the actual NMR signal when starting.As shown in Figure 2, E and digital received receiver output sampled point D has t
pTime error.And following relation is arranged:
0≤t
p<t
s
t
s=1/SW
Wherein, t
sIt is the time interval of each sampling; SW is a receiver sampling spectrum width.
This time error shows as the randomized jitter of NMR signal phase place on the unitary sampling of NMR signal, thereby has influence on adding up or offsetting of NMR signal.Then showing as the pseudo-shadow that produces on the phase-encoding direction on the Magnetic resonance imaging.
In order to reduce the randomized jitter of this phase place, a direct way is the bandwidth that increases receiver, makes SW increase, and reduces t
pTime error.But, can make main frame that the signal of receiver output is carried out filtering extraction once more to obtain the suitable visual field (FOV) like this, export the increase of spectrum width in addition along with receiver, need transmit data with computer data bus faster, with bigger memory headroom store data, this meeting is soft to system, and hardware proposes very high requirement.
Summary of the invention
The objective of the invention is provides a kind of method of reseptance of Magnetic resonance imaging signal according to above-mentioned the deficiencies in the prior art part, and this method is determined the deviate of phase place by setting up counter, signal is proofreaied and correct eliminated phase jitter again.
The object of the invention realizes being finished by following technical scheme:
A kind of method of reseptance of Magnetic resonance imaging signal, it is characterized in that the nuclear magnetic resonance digital receiver is provided with the counter that is used to measure each error time, in each sampling, this counter is measuring after the error time value sends to main frame, spectral line to NMR signal carries out the linear phase correction in the receiver spectrum width again, thereby eliminates the phase jitter of NMR signal.
Before carrying out the linear phase correction, need the NMR signal that collects is done fast Fourier transform (FFT), NMR signal is transformed from the time domain to frequency domain.
Spectral line is done linear phase in the described receiver spectrum width and proofread and correct, refer in frequency domain, be multiplied by phase correction factor α for the data point of each Frequency point correspondence of spectrum width scope.
Proofread and correct the nuclear magnetic resonance spectral line through linear phase and need carry out quick inverse-Fourier transform, with NMR signal from the frequency domain transform to the time domain, thereby obtain NMR signal through overcorrection.
Advantage of the present invention is, and compares by increasing the method that the receiver spectrum width reduces NMR signal phase jitter, fundamentally eliminated the NMR signal phase jitter that is caused by digital receiver; And under the prerequisite that minimal hardware is changed, eliminated the phase jitter of nuclear magnetic resonance, very succinct, the easy row of its implementation; Simultaneously, alleviated the computational burden of main frame greatly owing to used Fourier transform and inverse-Fourier transform fast to realize correction to NMR signal phase jitter.
Summary of drawings
Accompanying drawing 1 is digital receiver structured flowchart among the present invention;
Accompanying drawing 3 is NMR signal phase correction synoptic diagram among the present invention;
Accompanying drawing 4 is self-rotary echo-pulse series synoptic diagram among the present invention;
Accompanying drawing 5A is original nmr echo signal time domain synoptic diagram;
Accompanying drawing 5B is a nmr echo signal frequency domain synoptic diagram;
Accompanying drawing 5C is the echoed signal time domain synoptic diagram behind phase correction;
Accompanying drawing 6A is SE sequence (not adding a phase encoding) NMR signal raw data;
Accompanying drawing 6B is the data of the SE sequence (not adding phase encoding) behind phase correction;
Accompanying drawing 7 is a digital receiver structured flowchart in the prior art;
Concrete technical scheme
Feature of the present invention and other correlated characteristic are described in further detail by embodiment below in conjunction with accompanying drawing, so that technician's of the same trade understanding:
Shown in Fig. 1-7, label 1-11 represents respectively: probe (1), low-noise preamplifier (2), frequency mixer (3), intermediate-freuqncy signal conditioner (4), ADC (5), digital down converter (6), digital quadrature mixing device (7A) (7B), cascade comb filter (CIC) (8A) (8B), decimation filter (9A) (9B), main frame (10), counter (11), register (12).
On the nuclear magnetic resonance digital receiver, connect the counter (11) that is used to measure each error time, in each sampling, this counter (11) is measuring after the error time value sends to main frame (10), spectral line to NMR signal carries out the linear phase correction in the receiver spectrum width again, thereby eliminates the phase jitter of NMR signal.
Specific embodiments is general divided for three steps carried out:
The first step is done fast Fourier transform (FFT) to the NMR signal that collects, and NMR signal is transformed from the time domain to frequency domain.
Second step, in frequency domain, spectral line is done the one-level phase correction in the receiver spectrum width, promptly the data point for each Frequency point correspondence of spectrum width scope is multiplied by phase correction factor α.
The 3rd step, the nuclear magnetic resonance spectral line is carried out quick inverse-Fourier transform, obtain NMR signal through overcorrection.
From Fig. 2 as seen, between the data point D of the data point E of sampling startup moment correspondence and digital receiver output t is arranged
pTime error.If the frequency of NMR signal is f, then Shi Ji NMR signal can be expressed as:
A(t)=A(0)e
i2πf×t (1)
If digital receiver sampling spectrum width is SW, the NMR signal that then obtains at last is
Relatively 1 formula and 2 formulas can obtain, and for each signal sampling, the error of phase jitter is
0<Δφ<2πf/SW (3)
For each sampling, if error time t
pValue can measure, we just can do linear phase to NMR signal and correct by introducing the phase place modifying factor, the expression formula of phase place modifying factor α is as follows:
2 formulas be multiply by 4 formulas, can get
A′(t)×α=A(t) (5)
This shows, on the basis of raw data, multiply by phase correction factor and just can eliminate the phase error that causes by digital receiver.
Fig. 3 is a NMR signal phase correction synoptic diagram among the present invention.The counting clock welding system clock (frequency is generally tens megahertzes) of counter (11) wherein.The sampling enabling signal is used the position of rest of data frame signal as counter as the count enable signal of counter.When sampling started, at the rising edge of sampling enabling signal, counter began counting, stops counting when running into frame signal for high level, and the value of counter is deposited in register (12).Among Fig. 3, starting to counter from counter, to stop the used time be t
c, the count frequency of establishing counter is f
Clk, then the count value N of counter and error time t
pFollowing relation is arranged:
t
s-N/f
clk=t
p (6)
T wherein
sIn the time interval at expression digital receiver consecutive number strong point,, can measure t by reading the count value N of counter
pValue.
In order to eliminate error time t
pTo the influence that NMR signal causes, also can utilize phase correction factor to carry out that in time domain NMR signal is carried out convolution algorithm and eliminate phase jitter.
With spin echo (Spin Echo) sequence of Magnetic resonance imaging for for example shown in Figure 4, Gs wherein, Gp, Gr are gradient fields, their effect is that sample is carried out space encoding.90 ° and 180 ° is that radio-frequency pulse and sample magnetization vector are pulled the angle down, and its effect is the nuclear spin system of excited sample, produces NMR signal.After sampling started, digital receiver was gathered NMR signal.For the imaging sequence of Fig. 4, the frequency of its radio-frequency pulse is 10.7MHz, and the spectrum width SW of digital receiver is 20KHz, and the bandwidth of NMR signal is 2KHz, and sampling number TD is 256 points.Then according to 3 formulas, when receiving data, digital receiver produces the maximum phase error and is 36 ° of 2 π * 2KHz/20KHz=0.2 π.Among Fig. 4, be t=5ms the interval time of 90 ° and 180 ° radio-frequency pulses, and according to the nuclear magnetic resonance theory, the time t ' on the echo summit that its signal produces also should be 5ms.Fig. 5 A is the time domain parent magnetic resonance signal that digital receiver collects.Because the generation randomized jitter of NMR signal phase place, thus from the figure echo summit moment of 5ms not as can be seen.This shows that if NMR signal is not carried out phase correction, then can introduce serious phase place randomized jitter at the phase-encoding direction of nuclear magnetic resonance data and disturb, this must cause very serious image artifacts.For NMR signal is carried out phase correction, adopt following steps:
The first step is carried out obtaining behind the fast Fourier transform (FFT) nuclear magnetic resoance spectrum such as Fig. 5 B of frequency domain to the time domain NMR signal.
Second step, in frequency domain, nuclear magnetic resoance spectrum is carried out phase correction, its method is as follows:
Because the sampling bandwidth of receiver is 20KHz, sampling number is 256 points, and according to the character of Fourier transform, the frequency values f that obtains n data point on the spectrogram is
f
n=n×78.125Hz(0≤n<256) (7)
Then the phase correction factor of each data point can be expressed as in the spectrum width:
For this embodiment, phase correction factor can be expressed as:
α
n=e
i2π×n×78.125×tp (9)
Wherein, t
pIt is the time error value of each sampling.For this embodiment, the spectrum width of establishing digital receiver is 20KHz, then t
pMaximal value be 50us.When the system counts clock was 50MHz, count value N was 50us*50MHz=2500 to the maximum.For this sampling, the count value of establishing counter is 1000, then according to 6 formulas, obtains time error value t
p=50us-20us=30us.Again according to 9 formulas, each data point in the original spectrogram that digital receiver is obtained multiply by phase correction factor α
n
Calculate α
nAfter the factor, begin from 0Hz in the nmr spectrum of gained that the data of each Frequency point are multiplied by the respective phase correction factor, the nuclear magnetic resoance spectrum after so just obtaining proofreading and correct in the SW spectrum width scope.
The 3rd step, the nuclear magnetic resoance spectrum through overcorrect is done inverse-Fourier transform, obtain not have the time domain NMR signal of phase jitter, as Fig. 5 C, echo summit A has eliminated the phase jitter that each sampling causes through after the phase correction.
Fig. 6 A is a parent magnetic resonance K spatial data, each line data in the K space is corresponding to the sampling of a NMR signal, its signal amplitude is represented with gray scale, the high expression of brightness signal amplitude is big, the low expression of brightness signal amplitude is little, can see clearly because the each sampling of digital receiver causes the phase jitter of NMR signal from the K space.Fig. 6 B is the NMR signal K spatial data that obtains after this method for correcting phase of process of present embodiment is handled, can it is evident that,, eliminated to a large extent after the NMR signal correction by the method because the signal phase shake of each sampling that digital receiver causes.
Claims (1)
1, a kind of method of reseptance of Magnetic resonance imaging signal is characterized in that the nuclear magnetic resonance digital receiver is provided with timer conter, in each sampling, and the each sampling time error value of this counter measures; By the sampling time error value that usage counter obtains, in the receiver spectrum width, the nuclear magnetic resonance spectral line is carried out linear phase and proofread and correct; The time error value of utilizing counter to obtain calculates phase correction factor α; Before carrying out the linear phase correction, need the NMR signal that collects is done fast fourier transform, NMR signal is transformed from the time domain to frequency domain; In the receiver spectrum width, spectral line is done linear phase and proofread and correct, promptly in frequency domain, be multiplied by phase correction factor α for the data point of each Frequency point correspondence of spectrum width scope; The nuclear magnetic resonance spectral line of proofreading and correct through linear phase need carry out quick inverse-Fourier transform, with NMR signal from the frequency domain transform to the time domain.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100531543A CN100346171C (en) | 2004-07-23 | 2004-07-23 | Receiving method of nuclear magnetic resonance imaging signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100531543A CN100346171C (en) | 2004-07-23 | 2004-07-23 | Receiving method of nuclear magnetic resonance imaging signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1588113A CN1588113A (en) | 2005-03-02 |
CN100346171C true CN100346171C (en) | 2007-10-31 |
Family
ID=34602753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100531543A Expired - Fee Related CN100346171C (en) | 2004-07-23 | 2004-07-23 | Receiving method of nuclear magnetic resonance imaging signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100346171C (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103195B (en) * | 2009-12-18 | 2014-01-15 | 东软飞利浦医疗设备系统有限责任公司 | Device and method for realizing broadband digital magnetic resonance radio frequency receiving |
US10302716B2 (en) * | 2014-09-25 | 2019-05-28 | Koninklijke Philips N.V. | Digital receiver coil with built-in received phase noise indicator |
CN104836547B (en) * | 2015-06-05 | 2017-09-19 | 中国科学院武汉物理与数学研究所 | A kind of short group delay digital filtering method |
CN113466280B (en) * | 2018-02-27 | 2022-07-26 | 华东师范大学 | Simulated nuclear magnetic resonance spectrum analysis method and system convenient for expanding molecular information base and application thereof |
CN112014781B (en) * | 2020-09-02 | 2021-04-20 | 无锡鸣石峻致医疗科技有限公司 | Phase correction method and device for magnetic resonance echo signals, computer equipment and computer readable storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4673880A (en) * | 1985-08-16 | 1987-06-16 | Technicare Corporation | Phase sensitive detection in multislice magnetic resonance imaging systems |
US4857849A (en) * | 1987-03-25 | 1989-08-15 | Mitsubishi Denki Kabushiki Kaisha | High frequency magnetic field generator for nuclear magnetic resonance |
CN1289920A (en) * | 1999-09-28 | 2001-04-04 | 通用电器横河医疗系统株式会社 | NMR imaging device |
CN1433737A (en) * | 2002-01-23 | 2003-08-06 | Ge医疗系统环球技术有限公司 | Magnetic resonant imaging device |
-
2004
- 2004-07-23 CN CNB2004100531543A patent/CN100346171C/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4673880A (en) * | 1985-08-16 | 1987-06-16 | Technicare Corporation | Phase sensitive detection in multislice magnetic resonance imaging systems |
US4857849A (en) * | 1987-03-25 | 1989-08-15 | Mitsubishi Denki Kabushiki Kaisha | High frequency magnetic field generator for nuclear magnetic resonance |
CN1289920A (en) * | 1999-09-28 | 2001-04-04 | 通用电器横河医疗系统株式会社 | NMR imaging device |
CN1433737A (en) * | 2002-01-23 | 2003-08-06 | Ge医疗系统环球技术有限公司 | Magnetic resonant imaging device |
Also Published As
Publication number | Publication date |
---|---|
CN1588113A (en) | 2005-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu et al. | An enhanced empirical wavelet transform for noisy and non-stationary signal processing | |
CN104897962B (en) | The short sample high-precision frequency measurement method of simple signal and its device based on coprime perception | |
CN100504400C (en) | Oscilloscope high speed signal reconstruction method | |
US20140247176A1 (en) | Extension of adc dynamic range using post-processing logic | |
CN110837002B (en) | Spectrum scanning measuring device and time domain waveform obtaining method | |
CN100346171C (en) | Receiving method of nuclear magnetic resonance imaging signal | |
CN1956339A (en) | Device for testing an analog-/digital converter | |
CN1372718A (en) | Digit filter design | |
CN1645163A (en) | Generating method for linear digital frequency modulation signal | |
CN1480736A (en) | Method for testing electronic component and its instrument | |
CN111224672A (en) | Multi-harmonic signal undersampling method based on multi-channel time delay | |
CN108181486B (en) | The processing method and processing device of acceleration signal | |
CN109308453A (en) | Undersampled signal frequency estimating methods and device based on pattern clustering and spectrum correction | |
CN1276262C (en) | Current sensor | |
CN111812404B (en) | Signal processing method and processing device | |
CN109541309B (en) | Spectrum analyzer and signal processing method thereof | |
CN111521110B (en) | Rotary transformer signal envelope detection method | |
CN108535542B (en) | Peak-seeking phase discrimination method | |
Akujuobi et al. | Wavelet-based differential nonlinearity testing of mixed signal system ADCs | |
Hongwei | Fft basics and case study using multi-instrument | |
US10848168B1 (en) | Real-time digital waveform averaging with sub-sampling resolution | |
CN1761882A (en) | Wave detection device, method, program, and recording medium | |
CN210469275U (en) | Analog-to-digital converter | |
CN111865312B (en) | Analog-digital local oscillator synchronization method for digital bandwidth alternating system | |
CN103251428B (en) | Ultrasonic scanning system and blocking filter module and method for same |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20071031 Termination date: 20100723 |