CN106249183B - One kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method - Google Patents
One kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method Download PDFInfo
- Publication number
- CN106249183B CN106249183B CN201610846272.2A CN201610846272A CN106249183B CN 106249183 B CN106249183 B CN 106249183B CN 201610846272 A CN201610846272 A CN 201610846272A CN 106249183 B CN106249183 B CN 106249183B
- Authority
- CN
- China
- Prior art keywords
- steps
- xenon
- dwi
- hyperpolarization
- cssr
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Include the CSSR steps excited for hyperpolarization solubilised state xenon the present invention relates to one kind based on spectrum as integrated hyperpolarization xenon magnetic resonance method, collects magnetic resonance spectrogram, and for the DWI steps of hyperpolarization gaseous xenon excitation, collect nuclear magnetic resonance image.When the default swap time in CSSR steps being more than the repetition time of DWI steps, operation DWI steps are preset in swap time at this；The data that CSSR steps obtain are filled into the first data space, the data that DWI steps obtain are filled into the second data space；One-dimensional Fast Fourier Transform is carried out to the data in the first data space and obtains magnetic resonance spectrogram, nuclear magnetic resonance image is obtained to carrying out fast two-dimensional fourier transformation after the data permutation in the second data space successively.The present invention can speed up and synchronous acquisition obtains " gas diffusion " and " qi and blood exchange " functional information of lung's hyperpolarized xenon.
Description
Technical field
The present invention relates to mr imaging technique fields, more particularly, to one kind based on spectrum as integrated hyperpolarization xenon
Magnetic resonance method.Lung's magnetic resonance imaging suitable for hyperpolarization xenon as contrast agent.
Background technology
With the increase of number of smokers, the aggravation of air pollution (such as haze), the situation of pulmonary disease is worsening.Chest
(X-ray), computed tomography (Computer Tomography, CT) etc. are current conventional pulmonary disease clinic shadows thoroughly
As diagnostic techniques, still, they have certain radioactivity.In contrast, magnetic resonance imaging (Magnetic Resonance
Imaging, MRI) technology do not have radioactivity, and still, lung is one " blind area " of routine MRI.Use hyperpolarization xenon (Xe-
129) it is a kind of new technique to be used as contrast agent, can realize the MRI of lung mechanics.Moreover, hyperpolarization xenon can be dissolved in lung
Tissue and blood (see Fig. 1).Lung functions MRI is carried out, the important functional information of lung can be obtained.
Currently, being detected using hyperpolarization xenon MRI there are two the main methods of lung functions information.One is to be directed to super
Change gaseous xenon excitation Diffusion-Weighted MR Imaging (Diffusion Weighted Imaging, DWI) method (Ouriadov et al.,
Lung morphometry using hyperpolarized Xe-129 apparent diffusion coefficient
anisotropy in chronic obstructive pulmonary disease,Magnetic Resonance In
Medicine,2013,70(6):p.1699-1706.DOI:10.1002/mrm.24595), detection hyperpolarization xenon is in lung
Spread information, diffusion of the simulation oxygen in lung, in conjunction with corresponding lung model (the Weibel models in such as article), Jin Erneng
Estimate the lung functions parameters such as the average-size of alveolar, bronchiole.The other is for the change of hyperpolarization solubilised state xenon excitation
Displacement study saturation recovery (Chemical Shift Saturation Recovery, CSSR) method (Stewart, Neil J etc.
People, Experimental validation of the hyperpolarized Xe-129 chemical shift
saturation recovery technique in healthy volunteers and subjects with
interstitial lung disease,Magnetic Resonance In Medicine,2015,74(1):p.196-
207.DOI:10.1002/mrm.25400), detection hyperpolarization xenon, which enters from alveolar in tissue or blood, becomes hyperpolarization dissolving
The process of state xenon, simulation oxygen transports " qi and blood exchange " process into blood from alveolar, in conjunction with corresponding lung model (such as text
MOXE models in chapter), and then can estimate the lung functions parameters such as lung spacer thickness, average qi and blood swap time.
In above-mentioned two method, the detection of DWI methods is lung " gas diffusion " function, can obtain alveolar, ramuscule gas
The lung functions parameter such as average-size of pipe；The detection of CSSR methods is lung " qi and blood exchange " function, can obtain lung interval
The lung functions parameters such as thickness, average qi and blood swap time.The early detection of pulmonary disease usually requires to refer to comprehensive lung
Functional parameter, if it is possible to " gas diffusion " and " qi and blood exchange " two functions of lung are detected simultaneously, it will be to lung's disease
The early detection of disease provides help.However, the acquisition of DWI methods is lung image, the acquisition of CSSR methods is lung's spectrogram, by
In current existing hyperpolarization xenon magnetic resonance method, image and spectrogram are needed to separate in sampling at 2 times and be completed, therefore, such as
Fruit will obtain " gas diffusion " and " qi and blood exchange " two functions of patient lung simultaneously, there is several schemes below at present：
Scheme one first runs CSSR methods in same single-breath hold, and after waiting CSSR method end of runs, rerun DWI
Method；
Scheme two first runs DWI methods in same single-breath hold, after waiting DWI method end of runs, the side CSSR of reruning
Method；
Scheme three is separately operable CSSR methods and DWI methods in different hold one's breath.
In the first scheme, due to operation CSSR methods during the signal of hyperpolarization xenon have relaxation loss, because
Signal is too low when this will cause to run DWI methods.In second scheme, due to operation DWI methods during hyperpolarization xenon
Signal be excited loss, signal is too low when causing to run CSSR methods.In addition, whether above-mentioned the first or second of side
Case, 2 methods of continuous operation can all lead to single breath-hold time lengthening, cause the discomfort of patient.In the third scheme, Bing Renxu
2 hyperpolarization xenon, screen gas for sucking 2 multiple doses, increase the time for collecting hyperpolarization xenon and material cost, cause disease
People diagnoses the raising of cost and the reduction of diagnosis efficiency.
Therefore, it is necessary to a kind of magnetic resonance methods that can realize CSSR methods and the synchronous operation of DWI methods.
In order to solve the problems in above-mentioned background technology, the present invention proposes one kind based on spectrum as integrated hyperpolarization xenon
Magnetic resonance method, to collect nuclear magnetic resonance image and spectrogram simultaneously in unitary sampling so that CSSR methods and DWI methods
It can run simultaneously, to " gas diffusion " and " qi and blood exchange " functional information of synchronization gain lung.The method of the present invention can allow
As soon as patient only sucks time hyperpolarization xenon, only need unitary sampling that can be detected simultaneously by " gas diffusion " and " the qi and blood friendship of lung
Change " function, and total sampling time is shorter, shortens the breath holding time needed for patient, and ensure when two modules are run all
Possess higher magnetic resonance signal.
Invention content
It is an object of the invention in view of the deficienciess of the prior art, providing one kind based on spectrum as integrated hyperpolarization
Xenon magnetic resonance method.
The object of the invention is realized by following technical proposals：
One kind, as integrated hyperpolarization xenon magnetic resonance method, is included the following steps based on spectrum：
Step 1, the total number of repetition m for initializing CSSR steps initialize the default swap time τ of CSSR steps1~τm,
Default swap time τ1~τmBetween magnitude relationship be τ1<τ2<…<τm；
Step 2, the total excitation number Z for initializing DWI steps initialize the repetition time TR of DWI steps；
Step 3 is determined in default swap time τ according to the following formulakInterior energy run DWI steps can number of run fk
Wherein, square brackets [] expression takes the maximum positive integer less than numerical value in square brackets；
Step 4 is determined in default swap time τ according to the following formulakThe actual motion frequency n of interior operation DWI stepsk
K ∈ (1,2 ..., m)
In above formula, brace { } indicates, if the expression formula in brace is true, to take 1；If the expression formula in brace
It is not true, takes 0；
Step 5 determines the rear number of run h that DWI steps are run after the m times CSSR steps operation according to the following formula
In above formula, brace { } indicates, if the expression formula in brace is true, to take 1；If the expression formula in brace
It is not true, takes 0；
Step 6, the actual exchange time T for determining CSSR steps according to the following formulak
Tk=τk-TR×nk, k ∈ (1,2 ..., m)
Step 7, subject suck hyperpolarization Xe-129 gases and hold one's breath；
K is set as 1；
Saturation process in step 8, operation kth time CSSR steps；
Step 9, according to step 4 determine in default swap time τkThe actual motion frequency n of interior operation DWI stepsk, fortune
Row nkThe picture sampling step of secondary DWI steps obtains echo-signal Echo, and echo-signal Echo is stored in the second data space,
In nkIn the picture sampling step of secondary DWI steps, phase encoding gradient P1 and Diffusion Marker gradient DW1 are sequentially increased since 0；
The actual exchange time T of step 10, the kth time CSSR steps determined according to step 6k, in actual exchange time TkIt is interior
The exchange step in CSSR steps is run,
Exchange step in CSSR steps is specific as follows：Without any operation；
Spectrum sampling step in step 11, operation kth time CSSR steps, gains freedom and incudes deamplification (FID) and deposit
Enter in the first data space；
If step 12, k are more than or equal to m, 13 are entered step；If k is less than m, return to step 8 after k adds 1；
Step 13, the rear number of run that DWI steps are run after the m times CSSR steps operation determined according to step 5
H runs the picture sampling step of h DWI step, obtains echo-signal Echo, and echo-signal Echo is stored in the second data sky
Between, in the picture sampling step of h DWI step, phase encoding gradient P1 and Diffusion Marker gradient DW1 increase since 0 successively
Greatly；
Free induction decay signal (FID) in step 14, the first data space of deposit acquired to step 11 does one-dimensional
Fast Fourier Transform (1Dimensional Fast Fourier Transform, 1D FFT) obtains m magnetic resonance spectrogram；
Echo-signal Echo in step 15, the second data space of deposit acquired to step 9, step 13 does two-dimentional quick
Fourier transform (2Dimensional Fast Fourier Transform, 2D FFT) obtains corresponding nuclear magnetic resonance image.
Saturation process in step 8 as described above in CSSR steps is specific as follows：MR imaging apparatus emits first
One solubilised state xenon selective excitation pulse α1, solubilised state xenon selective excitation pulse α1Excitation center frequency alignment hyperpolarization
The resonant frequency of solubilised state xenon, solubilised state xenon selective excitation pulse α1Excitation angle be 90 °, and then magnetic resonance imaging is set
It is standby to apply a damage gradient G 1 on the direction perpendicular to main field direction.
The picture sampling step of DWI steps in step 9 and step 13 as described above is specific as follows, MR imaging apparatus
Emit a gaseous xenon selective excitation pulse α first3, gaseous xenon selective excitation pulse α3Excitation center frequency alignment it is super
The resonant frequency of polarized gaseous xenon, gaseous xenon selective excitation pulse α3Excitation angle be θ °, select layer terraced and at the same time applying
S1 is spent, to one imaging level of selection；MR imaging apparatus then applies a forward direction in direction initialization and one reversed
Pairs of Diffusion Marker gradient DW1 be diffused label；MR imaging apparatus then applies on phase encoding gradient direction
One phase encoding gradient P1 simultaneously applies a dephasing gradient R1 on readout gradient direction；Then it is applied on readout gradient direction
Add a readout gradient R2, and correspondingly carries out obtaining echo-signal Echo, echo-signal Echo the second numbers of deposit as sampling E2
According to space；Finally respectively apply a damage gradient G 2 in slice selective gradient direction, phase encoding gradient direction, readout gradient direction.
Excitation angle, θ ° as described above is not more than 10 °.
Spectrum sampling step in CSSR steps as described above is specific as follows：MR imaging apparatus transmitting one is non-selection
Property excitation pulse α2, non-selective excitation pulse α2Excitation center frequency alignment hyperpolarization solubilised state xenon resonant frequency, non-choosing
Selecting property excitation pulse α2The excitation of hyperpolarization solubilised state xenon and hyperpolarization gaseous xenon are excited, spectrum sampling is correspondingly carried out
E1 obtains the free induction decay signal (FID) for including hyperpolarization solubilised state xenon and hyperpolarization gaseous xenon, by free induction decay
Signal (FID) is stored in the first data space.
The present invention has the advantages that compared with the existing technology：
1, using spectrum as integrated mode, the spectrogram of CSSR and the image of DWI are collected simultaneously in unitary sampling,
So as to obtain " gas diffusion " and " qi and blood exchange " functional information of lung simultaneously.
2, the present invention can shorten and always adopt by the way of it will be run in the swap time of DWI steps insertion CSSR steps
The sample time,
3, DWI data and CSSR data can be by synchronous acquisitions in the present invention, therefore can when the operation of CSSR and DWI steps
Obtain higher magnetic resonance signal.
Description of the drawings
Fig. 1 is magnetic resonance signal distribution map of the hyperpolarization xenon in live body lung；
Fig. 2 is a kind of time sequencing figure based on spectrum as integrated hyperpolarization xenon magnetic resonance method of the present invention；
Fig. 3-1 is a kind of spectrogram signal obtained as integrated hyperpolarization xenon magnetic resonance method based on spectrum of the present invention
Figure；
Fig. 3-2 is a kind of image signal obtained as integrated hyperpolarization xenon magnetic resonance method based on spectrum of the present invention
Figure.
Specific implementation mode
Detailed description of the present invention is illustrative with reference to the accompanying drawings, non-limiting embodiment, to of the invention a kind of based on spectrum
As the technical solution of integrated hyperpolarization xenon magnetic resonance method is described in further detail, in order to fully understand this hair
It is bright.
Embodiment 1：
The present invention is based on spectrums as integrated mode, and the fortune of CSSR methods and DWI methods is completed at the same time in unitary sampling
Row.The present invention step be：
Total number of repetition m of step 1, initialization CSSR steps (chemical shift saturation recovery step), initialization CSSR steps
Suddenly the default swap time τ of (chemical shift saturation recovery step)1~τm, preset swap time τ1~τmBetween magnitude relationship
It is τ1<τ2<…<τm；
Total excitation number Z of step 2, initialization DWI steps (Diffusion-Weighted MR Imaging step), initialization DWI step (diffusion
Weighted imaging step) repetition time TR；
Step 3 is determined in default swap time τ according to the following formulakInterior energy runs DWI steps (Diffusion-Weighted MR Imaging step)
It can number of run fk
Wherein, square brackets [] expression takes the maximum positive integer less than numerical value in square brackets；
Step 4 is determined in default swap time τ according to the following formulakThe reality of interior operation DWI steps (Diffusion-Weighted MR Imaging step)
Border number of run nk
K ∈ (1,2 ..., m)
In above formula, brace { } indicates, if the expression formula in brace is true, to take 1；If the expression formula in brace
It is not true, takes 0；Z is total excitation number of DWI steps (Diffusion-Weighted MR Imaging step)；
Step 5 determines the operation after the m times CSSR steps (chemical shift saturation recovery step) operation according to the following formula
The rear number of run h of DWI steps (Diffusion-Weighted MR Imaging step)
In above formula, brace { } indicates, if the expression formula in brace is true, to take 1；If the expression formula in brace
It is not true, takes 0；
Step 6, the actual exchange time T for determining CSSR steps (chemical shift saturation recovery step) according to the following formulak
Tk=τk-TR×nk, k ∈ (1,2 ..., m)
Step 7, subject suck hyperpolarization Xe-129 gases and hold one's breath.
Subject can be patient, can also be the animal of scientific research；Air-breathing and the process of holding one's breath can actively be sucked simultaneously
It holds, can also be passively to suck and hold by lung ventilator.
K is set as 1.
Saturation process in step 8, operation kth time CSSR steps (chemical shift saturation recovery step)：
Saturation process in CSSR steps (chemical shift saturation recovery step) is specific as follows：MR imaging apparatus is first
Emit a solubilised state xenon selective excitation pulse α1, in the present embodiment, solubilised state xenon selective excitation pulse α1Select Gaussian
Shape, solubilised state xenon selective excitation pulse α1Excitation center frequency alignment hyperpolarization solubilised state xenon resonant frequency (this reality
It is 204ppm to apply in example, as shown in Figure 1), solubilised state xenon selective excitation pulse α1Excitation angle be 90 °.Solubilised state xenon selects
Selecting property excitation pulse α1Only hyperpolarization solubilised state xenon is excited, hyperpolarization gaseous xenon is not excited.To hyperpolarization solubilised state xenon
Signal by all be turned in the plane in main field direction.And then MR imaging apparatus is perpendicular to main field
Apply a damage gradient G 1 on the direction in direction；After damage gradient G 1 all damages the signal of hyperpolarization solubilised state xenon, this
When there is no any hyperpolarization solubilised state xenon signal.
Step 9, according to step 4 determine in default swap time τkInterior operation DWI steps (Diffusion-Weighted MR Imaging step)
Actual motion frequency nk, run nkThe picture sampling step of secondary DWI steps (Diffusion-Weighted MR Imaging step), in nkSecondary DWI steps
In the picture sampling step of (Diffusion-Weighted MR Imaging step), phase encoding gradient P1 and Diffusion Marker gradient DW1 since 0 successively
Increase：
The picture sampling steps of DWI steps (Diffusion-Weighted MR Imaging step) it is specific as follows (wherein, Diffusion Marker gradient DW1 and
Phase encoding gradient P1 is variable.The picture sampling step of the corresponding DWI steps (Diffusion-Weighted MR Imaging step) to rerun each time
Suddenly, phase encoding gradient P1 and Diffusion Marker gradient DW1 are sequentially increased since 0)：MR imaging apparatus emits one first
Gaseous xenon selective excitation pulse α3, in the present embodiment, gaseous xenon selective excitation pulse α3Select Gaussian in shape, gaseous xenon
Selective excitation pulse α3Excitation center frequency alignment hyperpolarization gaseous xenon resonant frequency (in the present embodiment be 0ppm, such as
Shown in Fig. 1), gaseous xenon selective excitation pulse α3Be θ ° (θ ° is the low-angle no more than 10 °), the gaseous xenon choosing of excitation angle
Selecting property excitation pulse α1Only hyperpolarization gaseous xenon is excited, hyperpolarization solubilised state xenon is not excited.Select layer terraced and at the same time applying
S1 is spent, to one imaging level of selection；MR imaging apparatus (is then phase code side in the present embodiment in direction initialization
To) apply a forward direction and a reversed pairs of Diffusion Marker gradient DW1 (i.e. gradient in Fig. 2 in dashed rectangle), it is used for
It is diffused label；MR imaging apparatus then on phase encoding gradient direction apply a phase encoding gradient P1 and
Apply a dephasing gradient R1 on readout gradient direction；Then apply a readout gradient R2, and phase on readout gradient direction
It carries out with answering obtaining echo-signal Echo as sampling E2, echo-signal Echo is stored in the second data space；Finally in slice selective gradient
Direction, phase encoding gradient direction, readout gradient direction respectively apply a damage gradient G 2, will be perpendicular to the gas in main field direction
State xenon signal is damaged.
When step 10, the actual exchange of the kth time CSSR steps (chemical shift saturation recovery step) determined according to step 6
Between Tk, in actual exchange time TkExchange step in interior operation CSSR steps (chemical shift saturation recovery step)：
Exchange step in CSSR steps (chemical shift saturation recovery step) is specific as follows：Without any operation.
Spectrum sampling step in step 11, operation kth time CSSR steps (chemical shift saturation recovery step)：
Spectrum sampling step in CSSR steps (chemical shift saturation recovery step) is specific as follows：MR imaging apparatus is sent out
Penetrate a non-selective excitation pulse α2, non-selective excitation pulse α2Excitation center frequency alignment hyperpolarization solubilised state xenon
Resonant frequency (being 204ppm in the present embodiment, as shown in Figure 1), non-selective excitation pulse α2Excitation angle be 90 °, non-choosing
Selecting property excitation pulse α2The excitation of hyperpolarization solubilised state xenon and hyperpolarization gaseous xenon are excited, spectrum sampling is correspondingly carried out
E1 just can obtain the free induction decay signal (FID) comprising hyperpolarization solubilised state xenon and hyperpolarization gaseous xenon, these are freely felt
Deamplification (FID) is answered to be stored in the first data space.
If step 12, k are more than or equal to m, 13 are entered step；If k is less than m, return to step 8 after k adds 1；
Step 13 is transported according to what step 5 determined after the m times CSSR steps (chemical shift saturation recovery step) operation
The rear number of run h of row DWI steps (Diffusion-Weighted MR Imaging step) runs the picture of h DWI step (Diffusion-Weighted MR Imaging step)
Sampling step, in the picture sampling step of h DWI step (Diffusion-Weighted MR Imaging step), phase encoding gradient P1 and diffusion are marked
Note gradient DW1 is sequentially increased since 0：
Specifically, DWI steps (Diffusion-Weighted MR Imaging step) herein as the DWI steps in sampling step and step 9
The picture sampling step of (Diffusion-Weighted MR Imaging step) is consistent.The picture sampling step of DWI steps (Diffusion-Weighted MR Imaging step) is specific
(wherein, Diffusion Marker gradient DW1 and phase encoding gradient P1 is variable as follows.The corresponding DWI steps to rerun each time
The picture sampling step of (Diffusion-Weighted MR Imaging step), phase encoding gradient P1 and Diffusion Marker gradient DW1 increase successively since 0
Greatly)：MR imaging apparatus emits a gaseous xenon selective excitation pulse α first3, in the present embodiment, gaseous xenon selectivity
Excitation pulse α3Select Gaussian in shape, gaseous xenon selective excitation pulse α3Excitation center frequency alignment hyperpolarization gaseous xenon
Resonant frequency (in the present embodiment be 0ppm, as shown in Figure 1), gaseous xenon selective excitation pulse α3Excitation angle be θ °
(θ ° is the low-angle no more than 10 °), gaseous xenon selective excitation pulse α1Only hyperpolarization gaseous xenon is excited, it is molten to hyperpolarization
Solution state xenon does not excite.And at the same time applying slice selective gradient S1, to one imaging level of selection；MR imaging apparatus then exists
Direction initialization (being phase-encoding direction in the present embodiment) applies a forward direction and a reversed pairs of Diffusion Marker gradient DW1
(i.e. gradient in Fig. 2 in dashed rectangle), for being diffused label；MR imaging apparatus is then in phase encoding gradient side
Apply a phase encoding gradient P1 upwards and applies a dephasing gradient R1 on readout gradient direction；Then in readout gradient
Apply a readout gradient R2 on direction, and correspondingly carry out obtaining echo-signal Echo as sampling E2, echo-signal Echo is deposited
Enter the second data space；Finally respectively apply a damage in slice selective gradient direction, phase encoding gradient direction, readout gradient direction
Gradient G 2, the gaseous xenon signal that will be perpendicular to main field direction are damaged.After this step end of run, subject's screen is not needed
Gas, thus subject can be with eupnea.
Free induction decay signal (FID) in step 14, the first data space of deposit acquired to step 11 does one-dimensional
Fast Fourier Transform (1 Dimensional Fast Fourier Transform, 1D FFT) obtains m magnetic resonance spectrogram；
Echo-signal Echo in step 15, the second data space of deposit acquired to step 9, step 13 does two-dimentional quick
Fourier transform (2 Dimensional Fast Fourier Transform, 2D FFT) obtains corresponding nuclear magnetic resonance image.
The nuclear magnetic resonance image that the magnetic resonance spectrogram and step 15 obtained according to step 14 obtains, respectively in connection with corresponding model meter
The functional information of corresponding " qi and blood exchange " and " gas diffusion ".
In attached embodiment shown in Fig. 2, RF, Gs、Gp、Gr, Acq indicate radio-frequency pulse, slice selective gradient, phase code respectively
Gradient, readout gradient and sample event.
Embodiment 2：
Below to some setting values in embodiment 1 to concrete values, to 3 in the method for the present invention and background technology
A scheme is compared.In the present embodiment：
The m=27 set in step 1,
[τ1,……,τm]=[0,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,
160,170,180,190,200,300,400,500,600,700,900]ms；
The Z=192 set in step 2, TR=30ms；
Then correspond to step 3
[f1,……,fm]=[0,0,0,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6,10,13,16,20,
23,30]；
Corresponding step 4
[n1,……,nm]=[0,0,0,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6,10,13,16,20,
23,30]；
The h=17 of corresponding step 5；
In corresponding step 6
[T1,……,Tm]=[0,10,20,0,10,20,0,10,20,0,10,20,0,10,20,0,10,20,0,10,
20,0,10,20,0,10,0]ms；
It is set in the present embodiment and runs saturation process in CSSR methods and compose the time of 1 needs of sampling step as Ty=
30ms。
Other are same as Example 1.
If to using scheme one, scheme two or the scheme three in background technology, total sampling time needs
But if using the scheme in the present invention, total sampling time needsAnd since CSSR methods in the present invention and DWI methods are all shifted to an earlier date,
Thus two methods can possess high RST.It is 5 ° to set the excitation angle, θ ° in DWI to gaseous xenon；Hyperpolarization xenon is in lung
Longitudinal relaxation time be 15 seconds, then the numerical value that the method for the present invention is compared with 3 schemes in background technology is shown in Table 1.It can
To see, the method for the present invention can collect CSSR signals and DWI signals in 1 time is held one's breath with shortest time synchronization, and
CSSR signals and DWI signals can keep high RST simultaneously.
Table 1 is compared using the effect of different schemes and the present invention program acquisition in background technology
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
With within principle, any modification, equivalent replacement, improvement and so on should all be included in the protection scope of the present invention god.
Claims (5)
1. one kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method, which is characterized in that include the following steps：
Step 1, the total number of repetition m for initializing CSSR steps initialize the default swap time τ of CSSR steps1~τm, preset
Swap time τ1~τmBetween magnitude relationship be τ1<τ2<…<τm；
Step 2, the total excitation number Z for initializing DWI steps initialize the repetition time TR of DWI steps；
Step 3 is determined in default swap time τ according to the following formulakInterior energy run DWI steps can number of run fk
Wherein, square brackets [] expression takes the maximum positive integer less than numerical value in square brackets；
Step 4 is determined in default swap time τ according to the following formulakThe actual motion frequency n of interior operation DWI stepsk
In above formula, brace { } indicates, if the expression formula in brace is true, to take 1；If the expression formula in brace is not
Very, 0 is taken；
Step 5 determines the rear number of run h that DWI steps are run after the m times CSSR steps operation according to the following formula
In above formula, brace { } indicates, if the expression formula in brace is true, to take 1；If the expression formula in brace is not
Very, 0 is taken；
Step 6, the actual exchange time T for determining CSSR steps according to the following formulak
Tk=τk-TR×nk, k ∈ (1,2 ..., m)
Step 7, subject suck hyperpolarization Xe-129 gases and hold one's breath；
K is set as 1；
Saturation process in step 8, operation kth time CSSR steps；
Step 9, according to step 4 determine in default swap time τkThe actual motion frequency n of interior operation DWI stepsk, run nk
The picture sampling step of secondary DWI steps obtains echo-signal Echo, and echo-signal Echo is stored in the second data space, in nk
In the picture sampling step of secondary DWI steps, phase encoding gradient P1 and Diffusion Marker gradient DW1 are sequentially increased since 0；
The actual exchange time T of step 10, the kth time CSSR steps determined according to step 6k, in actual exchange time TkInterior operation
Exchange step in CSSR steps,
Exchange step in CSSR steps is specific as follows：Without any operation；
Spectrum sampling step in step 11, operation kth time CSSR step, gaining freedom induction deamplification (FID) and is stored in the
In one data space；
If step 12, k are more than or equal to m, 13 are entered step；If k is less than m, return to step 8 after k adds 1；
Step 13, the rear number of run h that DWI steps are run after the m times CSSR steps operation determined according to step 5, fortune
The picture sampling step of h DWI step of row obtains echo-signal Echo, and echo-signal Echo is stored in the second data space,
In the picture sampling step of h DWI step, phase encoding gradient P1 and Diffusion Marker gradient DW1 are sequentially increased since 0；
Free induction decay signal (FID) in step 14, the first data space of deposit acquired to step 11 does one-dimensional quick
Fourier transform (1Dimensional Fast Fourier Transform, 1D FFT) obtains m magnetic resonance spectrogram；
Echo-signal Echo in step 15, the second data space of deposit acquired to step 9, step 13 does the quick Fourier of two dimension
Leaf transformation (2Dimensional Fast Fourier Transform, 2D FFT) obtains corresponding nuclear magnetic resonance image.
2. a kind of spectrum that is based on is as integrated hyperpolarization xenon magnetic resonance method according to claim 1, which is characterized in that institute
Saturation process in the step 8 stated in CSSR steps is specific as follows：MR imaging apparatus emits a solubilised state xenon choosing first
Selecting property excitation pulse α1, solubilised state xenon selective excitation pulse α1Excitation center frequency alignment hyperpolarization solubilised state xenon resonance
Frequency, solubilised state xenon selective excitation pulse α1Excitation angle be 90 °, and then MR imaging apparatus is perpendicular to main magnetic
Apply a damage gradient G 1 on the direction of field direction.
3. a kind of spectrum that is based on is as integrated hyperpolarization xenon magnetic resonance method according to claim 1, which is characterized in that institute
The picture sampling step of DWI steps in the step 9 and step 13 stated is specific as follows, and MR imaging apparatus emits a gas first
State xenon selective excitation pulse α3, gaseous xenon selective excitation pulse α3Excitation center frequency alignment hyperpolarization gaseous xenon be total to
Vibration frequency, gaseous xenon selective excitation pulse α3Excitation angle be θ °, and at the same time apply slice selective gradient S1, to selection one
A imaging level；MR imaging apparatus then applies a forward direction and a reversed pairs of Diffusion Marker ladder in direction initialization
Degree DW1 is diffused label；MR imaging apparatus then applies a phase encoding gradient on phase encoding gradient direction
P1 simultaneously applies a dephasing gradient R1 on readout gradient direction；Then apply a readout gradient on readout gradient direction
R2, and correspondingly carry out obtaining echo-signal Echo as sampling E2, echo-signal Echo is stored in the second data space；Finally selecting
Layer gradient direction, phase encoding gradient direction, readout gradient direction respectively apply a damage gradient G 2.
4. a kind of spectrum that is based on is as integrated hyperpolarization xenon magnetic resonance method according to claim 3, which is characterized in that institute
The excitation angle, θ ° stated is not more than 10 °.
5. a kind of spectrum that is based on is as integrated hyperpolarization xenon magnetic resonance method according to claim 1, which is characterized in that institute
Spectrum sampling step in the CSSR steps stated is specific as follows：MR imaging apparatus emits a non-selective excitation pulse α2,
Non-selective excitation pulse α2Excitation center frequency alignment hyperpolarization solubilised state xenon resonant frequency, non-selective excitation pulse
α2The excitation of hyperpolarization solubilised state xenon and hyperpolarization gaseous xenon are excited, spectrum sampling E1 is correspondingly carried out, obtains comprising super
The free induction decay signal (FID) of polarization solubilised state xenon and hyperpolarization gaseous xenon, free induction decay signal (FID) is stored in
First data space.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610846272.2A CN106249183B (en) | 2016-09-24 | 2016-09-24 | One kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610846272.2A CN106249183B (en) | 2016-09-24 | 2016-09-24 | One kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106249183A CN106249183A (en) | 2016-12-21 |
| CN106249183B true CN106249183B (en) | 2018-10-12 |
Family
ID=57611598
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610846272.2A Active CN106249183B (en) | 2016-09-24 | 2016-09-24 | One kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106249183B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107154065B (en) * | 2017-06-01 | 2019-08-16 | 中国科学院武汉物理与数学研究所 | A kind of quick function MR imaging method for human lung |
| CN108324275B (en) * | 2018-01-03 | 2021-05-28 | 上海东软医疗科技有限公司 | Method and device for acquiring magnetic resonance imaging signal and spectrum signal |
| CN109085522B (en) * | 2018-07-02 | 2021-07-09 | 上海东软医疗科技有限公司 | Method and device for acquiring magnetic resonance diffusion weighted imaging and spectrum signals |
| CN111161370B (en) * | 2019-12-30 | 2021-10-29 | 中国科学院精密测量科学与技术创新研究院 | Human body multi-core DWI joint reconstruction method based on AI |
| CN114236444B (en) * | 2021-12-03 | 2023-06-30 | 中国科学院精密测量科学与技术创新研究院 | Hyperpolarized gas lung variable sampling rate rapid magnetic resonance diffusion weighted imaging method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1224502A (en) * | 1996-03-29 | 1999-07-28 | 劳伦斯·伯克利国家实验室 | Enhancement of NMR and MRI in presence of hyperpolarized noble gases |
| CN101427918A (en) * | 2007-11-06 | 2009-05-13 | 中国人民解放军第二军医大学 | Magnetic resonance imaging method and uses thereof |
| CN103495251A (en) * | 2013-07-18 | 2014-01-08 | 中国科学院武汉物理与数学研究所 | Hyperpolarization gas breathing device for human body under non-narcosis state |
| CN105785293A (en) * | 2016-02-22 | 2016-07-20 | 中国科学院武汉物理与数学研究所 | Inert gas atomic nucleus channel apparatus and magnetic resonance imaging method |
-
2016
- 2016-09-24 CN CN201610846272.2A patent/CN106249183B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1224502A (en) * | 1996-03-29 | 1999-07-28 | 劳伦斯·伯克利国家实验室 | Enhancement of NMR and MRI in presence of hyperpolarized noble gases |
| CN101427918A (en) * | 2007-11-06 | 2009-05-13 | 中国人民解放军第二军医大学 | Magnetic resonance imaging method and uses thereof |
| CN103495251A (en) * | 2013-07-18 | 2014-01-08 | 中国科学院武汉物理与数学研究所 | Hyperpolarization gas breathing device for human body under non-narcosis state |
| CN105785293A (en) * | 2016-02-22 | 2016-07-20 | 中国科学院武汉物理与数学研究所 | Inert gas atomic nucleus channel apparatus and magnetic resonance imaging method |
Non-Patent Citations (2)
| Title |
|---|
| Regional mapping of gas uptake by blood and tissue in the human lung using hyperpolarized xenon‐129 MRI;Kun Qing et al.;《JOURNAL OF MAGNETIC RESONANCE IMAGING》;20140228;第39卷(第2期);346-359 * |
| 超级化<sup>129</sup>Xe磁共振肺成像应用研究;潘婷等;《临床放射学杂志》;20160131;第35卷(第1期);157-159 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106249183A (en) | 2016-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106249183B (en) | One kind is based on spectrum as integrated hyperpolarization xenon magnetic resonance method | |
| JP5537623B2 (en) | Magnetic resonance imaging system | |
| US9575153B2 (en) | MR imaging using a multi-point dixon technique | |
| KR101664138B1 (en) | System and method for tissue characterization using multi-slice magnetic resonance imaging | |
| US9113810B2 (en) | System and method for ungated non-contrast enhanced magnetic resonance angiography | |
| CN103800006B (en) | Motional error correction in functional magnetic resonance imaging | |
| US7885702B2 (en) | Segmentation of the airway tree using hyperpolarized noble gases and diffusion weighted magnetic resonance imaging | |
| EP3384307B1 (en) | Removal of image artifacts in sense-mri | |
| JP5269342B2 (en) | Magnetic resonance imaging apparatus and imaging condition setting method in magnetic resonance imaging apparatus | |
| US20160349339A1 (en) | System and method for magnetic resonance fingerprinting at high field strengths | |
| KR101809213B1 (en) | Method and apparatus for multi-slice imaging of t2-relaxation time | |
| JPH0871058A (en) | Magnetic resonance imaging apparatus | |
| CN102743172A (en) | Magnetic resonance method and apparatus for triggered acquisition of magnetic resonance measurement data | |
| US20140378826A1 (en) | System and Method For Selective Magnetic Resonance Imaging Angiography Of Arteries Or Veins | |
| CN109507622A (en) | MR imaging method and system | |
| CN106456046A (en) | Method of improved multiple-phase dynamic contrast-enhanced magnetic resonance imaging | |
| US10330760B2 (en) | System and method for assessing T2-relaxation times with improved accuracy | |
| US10739420B2 (en) | Motion-robust transmit radio frequency field mapping in magnetic resonance imaging using interleaved bloch-siegert shifting | |
| EP4057022A1 (en) | System and method for magnetization-prepared three-dimensional unbalanced steady-state free precession magnetic resonance imaging | |
| US20070055138A1 (en) | Accelerated whole body imaging with spatially non-selective radio frequency pulses | |
| US10459055B2 (en) | System and method for reduced field of view MR fingerprinting for parametric mapping | |
| JPH11313810A (en) | Mri device and mr image pickup method | |
| WO2020034674A1 (en) | Quantitative myocardial magnetic resonance imaging method and device, and storage medium | |
| CN106725508A (en) | Physiological movement collecting method, MR imaging method and device | |
| US11846655B2 (en) | System and method for three-dimensional (3D) unbalanced steady-state free precession magnetic resonance imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |