CN107582057A - MR imaging method and device - Google Patents
MR imaging method and device Download PDFInfo
- Publication number
- CN107582057A CN107582057A CN201710915393.2A CN201710915393A CN107582057A CN 107582057 A CN107582057 A CN 107582057A CN 201710915393 A CN201710915393 A CN 201710915393A CN 107582057 A CN107582057 A CN 107582057A
- Authority
- CN
- China
- Prior art keywords
- data
- multilayer
- imaging
- coil
- coefficient
- 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
Abstract
The invention provides a kind of MR imaging method and device.The MR imaging method includes：To the individual layer calibration data in K spaces and multilayer, imaging data carries out the first conversion simultaneously, and to obtain the first individual layer calibration data and the first multilayer while imaging data in intermediate field, the intermediate field is between the K spaces and image area；Coil is calculated according to the first individual layer calibration data and merges coefficient；Merging coefficient and first multilayer according to the coil, imaging data carries out data hierarchy simultaneously, to obtain multiple monolithic layer datas in intermediate field；And multiple monolithic layer datas in intermediate field are carried out with the second conversion, to obtain MRI.The MR imaging method and device have the characteristics that reconstruction amount of calculation is small, it is fast to rebuild speed.
Description
Technical field
The invention mainly relates to magnetic resonance imaging arts, more particularly to the image reconstruction that a kind of magnetic resonance multilayer is imaged simultaneously
Method and device.
Background technology
Magnetic resonance imaging (Magnetic Resonance Imaging, MRI) technology is in clinical diagnosis and scientific research
More and more extensive application has been obtained, there is safety, more contrasts and have good resolution capability to soft tissue.So
And compared with other medical imaging technologies (such as ultrasonic imaging, CT imagings), imaging time needed for mr imaging technique compared with
It is long, not only reduce the comfort level of patient, add sensitiveness of the picture quality to motion, and it is imaged in dynamic process
In application bring certain challenge.
Multilayer as described by document 1 is imaged (Simultaneous MultiSlice, SMS) technology simultaneously can be effective
Ground shortens the MRI scan time, has obtained very big development in recent years.The basic process of imaging technique is multilayer simultaneously
Excite multiple lamellas of space diverse location simultaneously first, be spatially encoded using gradient fields, pass through coil array again afterwards
Signal is received.Piece number of layers if necessary to collection is N, and the time that every layer of collection needs is t, then conventional magnetic resonance
Total time needed for imaging method scanning is N × t.On the other hand, if multilayer is imaged the lamella for exciting and gathering simultaneously simultaneously
Number is m, then scans and required total time shorten to N/m × t.Fig. 1 shows two lamellas while the piece for exciting while gathering
Layer schematic diagram.In scheme illustrated in fig. 1, lamella 1 and lamella N/2+1 are excited while gathered simultaneously, lamella 2 and lamella N/2
+ 2 are excited while gather simultaneously, and lamella 3 and lamella N/2+3 are excited while gathered simultaneously, and the rest may be inferred.
Different from conventional imaging method, the multilayer information that the data that directly gather of imaging technique contain multiple lamellas simultaneously,
It is simple to rebuild the aliasing artefacts that brought between multiple lamellas.Document 2 describes a kind of more common multilayer imaging reconstruction simultaneously
Method-slice-GRAPPA (GeneRalized Autocalibrating Partial Parallel Acquisition).
This process employs the difference of each coil sensitivities in receiving array, the calibration data being collected into is excited by advance individual layer,
Calculated in K spaces and the coil (passage) of aliased data layering is merged into coefficient (as shown in Figure 2 a), this is merged again afterwards
Layering is realized in the K space data that coefficient acting is imaged to multilayer simultaneously (as shown in Figure 2 b).Data after layering can then lead to
Cross traditional reconstructing method (such as Fourier transformation) and obtain image.
Imaging technique is to calculate coil in K spaces to merge coefficient to existing multilayer as described above simultaneously, and in K
Multilayer aliased image is layered in space, this image rebuilding method is computationally intensive, low to the read-write efficiency of internal memory, and then
It is slow to result in image reconstruction speed.
Document 1：Barth M,Breuer F,Koopmans PJ,Norris DG,Poser BA.Simultaneous
multislice(SMS)imaging techniques.Magn Reson Med.2016；75(1):63-81.
Document 2：Setsompop K,Gagoski BA,Polimeni JR,Witzel T,Wedeen VJ,Wald
LL.Blipped-controlled aliasing in parallel imaging for simultaneous
multislice echo planar imaging with reduced g-factor penalty.Magn Reson
Med.2012；67(5):1210-24.
The content of the invention
The technical problem to be solved in the present invention is to provide the image rebuilding method and dress that a kind of magnetic resonance multilayer is imaged simultaneously
Put, it has the characteristics that, and reconstruction amount of calculation is small, it is fast to rebuild speed.
In order to solve the above technical problems, the invention provides a kind of MR imaging method, including：To the list in K spaces
Imaging data carries out the first conversion simultaneously for layer calibration data and multilayer, to obtain the first individual layer calibration data in intermediate field and the
One multilayer while imaging data, the intermediate field is between the K spaces and image area；Number is calibrated according to first individual layer
Merge coefficient according to coil is calculated；Merging coefficient and first multilayer according to the coil, imaging data carries out data point simultaneously
Layer, to obtain the monolithic layer data in intermediate field；And multiple monolithic layer datas in intermediate field are carried out with the second conversion, to obtain
Obtain MRI.
In one embodiment of this invention, it is described that coil merging coefficient bag is calculated according to the first individual layer calibration data
Include：Apply first phase modulation to the first individual layer calibration data, to construct the first coefficient matrix；Adjusted according to phase is not applied
The first individual layer calibration data construction first object vector of system；And according to first coefficient matrix and first mesh
Mark coil described in Vector operation and merge coefficient.
In one embodiment of this invention, it is described that coefficient and first multilayer are merged according to the coil while are imaged number
Include according to data hierarchy is carried out：To first multilayer, imaging data applies second phase modulation simultaneously, to construct the second coefficient
Matrix；And the monolithic layer of coefficient and second coefficient matrix calculating acquisition in intermediate field is merged according to the coil
Data.
In one embodiment of this invention, described first it is transformed to the Fourier transformation along frequency coding direction.
In one embodiment of this invention, imaging data is by repeatedly exciting target area to obtain simultaneously for the multilayer, institute
State target area and include multiple scanning slices, and the quantity of the scanning slice is more than and excites number.
In one embodiment of this invention, the target area is at least performed and excites and excite for the second time for the first time, institute
State and excite corresponding first group of scanning slice for the first time, it is described to excite corresponding second group of scanning slice for the second time.
In one embodiment of this invention, first group of scanning slice or second group of scanning slice include spaced
Multiple scanning slices, and during second group of scanning slice is excited, imaging data corresponding to first group of scanning slice is entered
Row first converts.
In one embodiment of this invention, described second it is transformed to the Fourier transformation along phase-encoding direction.
A kind of MR imaging apparatus is provided in another aspect of this invention, including：First conversion module, for K
Imaging data carries out the first conversion simultaneously for individual layer calibration data and multilayer in space, to obtain the first individual layer school in intermediate field
Quasi- data and the first multilayer while imaging data, the intermediate field is between the K spaces and image area；Coil merges coefficient
Computing module, merge coefficient for calculating coil according to the first individual layer calibration data；Data hierarchy module, for according to institute
State coil and merge coefficient and first multilayer data hierarchy of imaging data progress simultaneously, to obtain the monolithic number of plies in intermediate field
According to；And
Second conversion module, for multiple monolithic layer datas in intermediate field to be carried out with the second conversion, to obtain magnetic resonance
Image.
In one embodiment of this invention, the coil merges coefficients calculation block and is configured to perform following steps：To institute
State the first individual layer calibration data and apply first phase modulation, to construct the first coefficient matrix；According to the institute for not applying phase-modulation
State the first individual layer calibration data construction first object vector；And according to first coefficient matrix and the first object vector
Calculate the coil and merge coefficient.
In one embodiment of this invention, the data hierarchy module is configured to perform following steps：More than described first
Imaging data applies second phase modulation to layer simultaneously, to construct the second coefficient matrix；And according to the coil merge coefficient and
Second coefficient matrix calculates the monolithic layer data obtained in intermediate field.
In one embodiment of this invention, described first it is transformed to the Fourier transformation along frequency coding direction.
In one embodiment of this invention, imaging data is by repeatedly exciting target area to obtain simultaneously for the multilayer, institute
State target area and include multiple scanning slices, and the quantity of the scanning slice is more than and excites number.
In one embodiment of this invention, the target area is at least performed and excites and excite for the second time for the first time, institute
State and excite corresponding first group of scanning slice for the first time, it is described to excite corresponding second group of scanning slice for the second time.
In one embodiment of this invention, first group of scanning slice or second group of scanning slice include spaced
Multiple scanning slices, and during second group of scanning slice is excited, imaging data corresponding to first group of scanning slice is entered
Row first converts.
In one embodiment of this invention, described second it is transformed to the Fourier transformation along phase-encoding direction.
A kind of MR imaging apparatus, including memory, processor and storage are provided in another aspect of this invention
On a memory and the computer instruction that can run on a processor, wherein, described in the computing device during computer instruction
Realize method as described above.
A kind of computer-readable medium is provided in another aspect of this invention, is stored thereon with computer instruction, wherein
When computer instruction is executed by processor, method as described above is performed.
A kind of magnetic resonance imaging system, including magnetic resonance imaging as described above dress are provided in another aspect of this invention
Put.
Compared with prior art, the present invention has advantages below：
Merge coefficient compared to traditional coil that calculated in K spaces, and multilayer aliased image is carried out in K spaces
The mode of layering, the present invention are that the calculating of coil merging coefficient is carried out in intermediate field, and to multilayer aliasing figure in intermediate field
As being layered.The read-write of data can be carried out in units of whole piece frequency coding line in intermediate field, improves the reading of memory
Write efficiency.It is also possible to the first conversion is carried out when gathering multilayer imaging data simultaneously, calculating coil merges coefficient
One or more of with data hierarchy, shorten from the time required for gathered data to output final image.
Brief description of the drawings
Fig. 1 is multilayer two lamellas of imaging technique while the lamella schematic diagram for exciting while gathering simultaneously.
Fig. 2 a, 2b are the basic step schematic diagrames of slice-GRAPPA multilayers while imaging reconstruction method.
Fig. 3 is the structural representation of magnetic resonance imaging system.
Fig. 4 is the schematic diagram of the basic process of magnetic resonance imaging.
Fig. 5 is the sub-step schematic diagram that image area is transformed into from K spaces.
Fig. 6 is the basic flow sheet of the magnetic resonance multilayer of one embodiment of the invention while the image rebuilding method of imaging.
Fig. 7 is the basic flow sheet for calculating coil and merging coefficient of one embodiment of the invention.
Fig. 8 is the basic flow sheet being layered to multilayer aliased image of one embodiment of the invention.
Fig. 9 is the magnetic resonance multilayer of one embodiment of the invention while the fundamental block diagram of the equipment for reconstructing image of imaging.
Figure 10 is the result schematic diagram being layered to multilayer aliased image of one embodiment of the invention.
Figure 11 is the structural representation of the magnetic resonance multilayer of one embodiment of the invention while the equipment for reconstructing image of imaging.
Figure 12 is the schematic diagram of the computer-readable medium of one embodiment of the invention.
Embodiment
For the above objects, features and advantages of the present invention can be become apparent, the tool below in conjunction with accompanying drawing to the present invention
Body embodiment elaborates.
Many details are elaborated in the following description to facilitate a thorough understanding of the present invention, still the present invention can be with
It is different from other manner described here using other to implement, therefore the present invention is not limited by following public specific embodiment
System.
As shown in the application and claims, unless context clearly prompts exceptional situation, " one ", "one", " one
The word such as kind " and/or "the" not refers in particular to odd number, may also comprise plural number.It is, in general, that term " comprising " is only prompted to wrap with "comprising"
Include clearly identify the step of and element, and these steps and element do not form one it is exclusive enumerate, method or equipment
The step of may also including other or element.
Magnetic resonance imaging system generally include to have the magnet of certain pore size, the transmitting coil for launching radiofrequency signal and
For receive the receiving coil of magnetic resonance signal, for carrying out sterically defined gradient coil to magnetic resonance signal, for generating
The impulse generator and control system of scanning sequence.Magnetic resonance imaging system is controlled and controlled by operator (clinician)
The console work of system connection, console may include keyboard or other input equipments, control panel and display, to input life
Order and the image of display generation.
Fig. 3 is the structural representation of magnetic resonance imaging system, and when carrying out magnetic resonance examination, clinician will be examined first
Person 3 is placed on scanning bed 1, and (is not shown in figure in the local coil that the body surface of person under inspection 3 is placed for receiving magnetic resonance signal
Go out)；Then the console that clinician is connected by operation and control system 5 is controlled in the scanning bed aperture formed toward magnet 2
It is mobile, after the instruction that magnetic resonance imaging system monitors that clinician sends scanning bed 1 movement, pass through control system 5 immediately
Scanning bed moving range is monitored, when the edge that scanning bed 1 enters scanning imagery region 4, pulse sequence is controlled by control system 5
Row generator generates corresponding sequence and is scanned.In scanning bed 1 moving process, connecing for person under inspection's body surface is positioned over
Take-up circle can move with scanning bed 1 in the inner space of magnet spaces, and the receiving coil in diverse location is in control system 5
Effect is lower in opening or closed mode, to receive corresponding magnetic resonance signal.
Information of each signal containing holostrome of magnetic resonance, it is therefore desirable to which space orientation volume is carried out to magnetic resonance signal
Code, i.e. frequency coding and phase code.The MR signals that magnetic resonance reception coil collects are really with space encoding information
Radio wave, belong to analog signal and nonnumerical information, it is necessary to by analog-to-digital conversion (analog-digital conversion,
ADC digital information) is become, the latter is filled to K spaces, turns into numerical data dot matrix.K spaces are with space orientation coding letter
The packing space of the MR signal initial numberical datas of breath, each width MR images have its corresponding K space data dot matrix.Magnetic resonance
The basic process of imaging mainly includes step as shown in Figure 4：First, pulse signal caused by radio-frequency drive, the pulse signal
Person under inspection's body region can be excited to produce precession nuclear spin；Secondly, by space encoding, receiving coil can collect simultaneously to be had
The response signal of space orientation coding information, the response signal represent MR letters caused by person under inspection's body region precession nuclear spin
Number；Then, these MR signals are filled to K spaces through analog-to-digital conversion；Finally, Fourier transformation is carried out to the data in K spaces, with regard to energy
Space orientation coding information in initial data is decoded, decomposites different frequency, the MR signals of phase and amplitude, it is different
Frequency and phase represent different locus, and amplitude then represents MR signal intensities, different frequencies, phase and signal
The MR digital informations of intensity are assigned in corresponding pixel, have just obtained MR view data, that is, reconstruct MR images.
It would be appreciated that Fourier transformation to K space data dot matrix, being assigned to the medium processing of each pixel and can be
It is performed in the von Neumann architecture computer such as laptop computer, desktop computer, server, tablet personal computer.Ability
Field technique personnel are further appreciated that these processing can also be by general processor, digital signal processor (DSP), special collection
Into circuit (ASIC), field programmable gate array (FPGA) or other PLDs, discrete door or transistor logic,
Discrete nextport hardware component NextPort or its any combinations are realized or performed.
Fig. 5 shows the sub-step schematic diagram that image area is transformed into from K spaces.The processing of image area is transformed into from K spaces
It can be generally divided into two sub-steps：(a) the first conversion is carried out to the data in K spaces, so as to obtain the number of corresponding intermediate field
According to；(b) the second conversion is carried out to the data of intermediate field, so as to obtain the image in image area.In practical operation, it can typically incite somebody to action
Linear transformation operation is merged into first conversion and the second conversion.Exemplified by performing two-dimensional Fourier transform to K space data, this two
Dimension Fourier transformation, which can resolve into, first have been carried out one-dimensional Fourier transform (corresponding first conversion) to K space data and obtains centre
The data in domain, the one-dimensional Fourier transform may include along phase-encoding direction or along frequency coding direction；Again to the number of intermediate field
The image in image area is obtained according to one-dimensional Fourier transform (corresponding second conversion) has been carried out again, the one-dimensional Fourier transform can wrap
Include along frequency coding direction or along phase-encoding direction.
As introduced in the introduction, imaging technique is to calculate coil in K spaces to merge system to conventional multilayer simultaneously
Number, and multilayer aliased image is layered in K spaces, this image rebuilding method need to be in raw k-space Data Collection
After start data processing, it is and low to the read-write efficiency of internal memory, and then it is slow to result in image reconstruction speed.It is existing in order to overcome
Disadvantages mentioned above present in technology, the present inventor proposes first carries out first to multilayer while the K space data of collection
Conversion, K space data is transformed in intermediate field, and the calculating that coil merges coefficient is then carried out in intermediate field, and in
Between multilayer aliased image is layered in domain, to obtain monolithic layer data, can finally carry out second to monolithic layer data and become
Change, obtain final image.The monolithic layer data that this kind of processing mode obtains is in intermediate field, can reduce subsequent reconstruction
Amount of calculation, and can cause simultaneously in calculating process more than read-write to internal memory in units of whole piece frequency coding line, it is blocking
Data processing improves the efficiency of memory read-write.Therefore the present invention can accelerate multilayer while the image reconstruction speed of imaging.
The feasibility of the present invention is proved by taking slice-GRAPPA methods as an example below.
The slice-GRAPPA methods that coil merging coefficient calculates and lamella separates are carried out in K spaces to briefly explain such as
Under.Before magnetic resonance imaging, sweep parameter is set for control system 5, the sweep parameter may include sweep parameter corresponding to individual layer
Corresponding sweep parameter is excited simultaneously with multilayer, further, corresponding to sweep parameter corresponding to individual layer or multilayer excite simultaneously
Sweep parameter may include gradient parameter, radio frequency parameter etc.；Control system 5 sends sweep parameter to spectrometer system, corresponding to generation
Gradient scanning sequence or rf pulse sequence；Gradient scanning sequence driving gradient coil produces corresponding gradient fields, radio-frequency pulse
Sequence driving radio-frequency coil produces corresponding radiofrequency field；The each layer in region to be detected of person under inspection produced under radio frequency field excitation into
Kinetonucleus spins, and carries out gradient coding to precession nuclear spin corresponding to each layer under gradient field action, produces corresponding to every layer
Magnetic resonance signal, the magnetic resonance signal are filled into K spaces, can obtain individual layer collection K space data.In one embodiment, it is single
Layer collection K space data can be corresponding every layer K space center region.Similarly, the region multilayer to be detected of person under inspection is (such as
Two layersth, three layers etc.) in radiofrequency field lower generation precession nuclear spin is excited simultaneously, and to precession corresponding to multilayer under gradient field action
Nuclear spin carries out gradient coding simultaneously, produces magnetic resonance signal corresponding to multilayer, and the magnetic resonance signal is filled into K spaces, can obtained
Obtain multilayer and gather K space data simultaneously.
In one embodiment, multichannel K space data individual layer excited is as calibration data.K spaces are being determined
After convolution kernel size, solving the coefficient matrix B that coil merges needed for coefficient can be generated by the calibration data in K spaces, and be contained
Multi-channel data.If the piece number of layers excited simultaneously every time is p, then coefficient matrix B can be expressed as (A1+A2+…Ap),
Wherein AiFor the monolithic layer gathered data of calibration data, i is the numbering of each lamella, 1≤i≤p.
The coil for needing to solve merges coefficient X and meets following relation：
BX=Atar (1)
Wherein X is that multiple coils (passage) merge coefficient, AtarIt is the plan extracted from the monolithic layer calibration data of K spaces
Close object vector.Here B, AtarFor, it is known that coil merging coefficient X can be calculated by following formula：
X=(BTB)-1BTAtar (2)
In formula, the T in the upper right corner represents to ask the transposition computing of matrix；- the 1 of the upper left corner represents to ask inverse of a matrix computing.
Multilayer is imaged simultaneously can be used to construct coefficient matrix C containing aspect aliased data.Isolated in note K spaces
Monolithic layer data is D, and the coil above calculated can be utilized to merge coefficient X and isolate monolithic layer data D：
D=CX (3)
The demonstration present invention is set up below.Distinguish premultiplication at equation (1) both ends to carry out in (segmentation) Fu along frequency coding direction
The matrix of leaf transformation, is obtained：
FroBX=FroAtar (4)
IfAF_tar=FroAtar, coil, which merges coefficient X, to be calculated as follows：
It can be seen that it is equal that formula (5) merges coefficient X with the coil that formula (3) calculates, wherein, ro represents frequency coding
Direction or readout direction；FroRepresent the Fourier transformation along frequency coding direction；AF_tarRepresent along frequency coding direction Fourier
Fit object vector after conversion.
When being layered to aspect aliased data, it is necessary to first build along the coefficient after the Fourier transformation of frequency coding direction
MatrixMonolithic layer is along the data D after frequency direction Fourier transformationF_tar(the monolithic number of plies i.e. in intermediate field
According to) can be calculated by following formula：
DF_tar=CFX (6)
For coefficient matrix B, C, its column direction data can be raw k-space data using segmented construction along frequency coding direction
Certain translation is carried out, maximal translation amount is determined by convolution kernel size.AtarIt is the raw k-space data without translation along frequency
Encode direction arrangement.According to the property of Fourier transformation, the result that the translation before conversion in K spaces can be given after conversion brings one
Linear phase-modulation.The power of linear phase modulation is relevant with K spatial translation sizes, in the translation of each unit of K spaces
2 π of adjacent pixel/N phase place change can be brought to the result of Fourier transformation, wherein N is the total pixel in Fourier transformation direction
Number, N is positive integer and 1≤N≤128.
From above-mentioned demonstration, the coil that the present invention calculates merges coefficient and the phase calculated using traditional approach
Together, the present invention calculates the monolithic layer data D in the intermediate field obtainedF_tarWith the monolithic in the K spaces that calculate traditional approach
Layer data D is transformed to obtained by the identical therefore proposed by the invention processing mode of the monolithic layer data obtained after intermediate field
Result is correct, and the correctness for the image finally rebuild will not be had an impact.
In addition, those skilled in the art it would be appreciated that, the above-mentioned processing side to K space data dot matrix proposed by the present invention
Formula be equally applicable to such as CAIPIRINHA (Breuer FA, Blaimer M, Heidemann RM, Mueller MF,
Griswold MA,Jakob PM.Controlled Aliasing in Parallel Imaging Results in
Higher Acceleration(CAIPIRINHA)for Multi-Slice Imaging.Magn Reson Med.2005；53
(3):684-91.)、SP-SG(Stephen F.Cauley,Jonathan R.Polimeni,Himanshu Bhat,Dingxin
Wang,Lawrence L.Wald,and Kawin Setsompop.Inter-slice Leakage Artifact
Reduction Technique for Simultaneous Multi-Slice Acquisitions.Magn Reson
Med.2014；72(1):93-102.) in the methods of.As space is limited, herein not to proposed by the present invention above-mentioned to K space data
The feasibility for the methods of processing mode of dot matrix is applied to CAIPIRINHA, SP-SG is proved.
Fig. 6 is the basic flow sheet of the magnetic resonance multilayer of one embodiment of the invention while the image rebuilding method of imaging.Magnetic
The image rebuilding method 100 that resonance multilayer is imaged simultaneously may include steps of：
Step 110：To the individual layer calibration data in K spaces and multilayer, imaging data carries out the first conversion simultaneously, to obtain
The first individual layer calibration data and the first multilayer while imaging data in intermediate field.Wherein, intermediate field is located at K spaces and image area
Between.
Step 120：Coil is calculated according to the first individual layer calibration data and merges coefficient.
Step 130：Merging coefficient and the first multilayer according to coil, imaging data carries out data hierarchy simultaneously, to obtain centre
Monolithic layer data in domain.
In one embodiment, the image rebuilding method 100 that magnetic resonance multilayer is imaged simultaneously can also include step 140：It is right
Monolithic layer data in intermediate field performs the second conversion, to obtain MRI.
In step 110, the first conversion can be the Fourier transformation along frequency coding direction.
In the step 120, sub-step as shown in Figure 7 can be included：
Sub-step 121：Apply first phase modulation to the first individual layer calibration data, to construct the first coefficient matrix；
Sub-step 122：The first individual layer calibration data according to phase-modulation is not applied constructs first object vector；
Sub-step 123：Coefficient is merged according to the first coefficient matrix and first object Vector operation coil.
It is appreciated that sub-step 121 and sub-step 122 can perform in any order, you can to first carry out sub-step
121, it is rear to perform sub-step 122；Can be rear to perform sub-step 121 with restricted sub-step 122；Sub-step 121 and sub-step
122 can also perform simultaneously.
In step 130, sub-step as shown in Figure 8 can be included：
Sub-step 131：To the first multilayer, imaging data applies second phase modulation simultaneously, to construct the second coefficient matrix；
Sub-step 132：Coefficient is merged according to coil and the second coefficient matrix calculates the monolithic number of plies obtained in intermediate field
According to.
It would be appreciated that sub-step 131 can after step 110, any time before sub-step 132 is performed.
For example, in one embodiment, sub-step 131 can be held simultaneously with any one or more in sub-step 121,122 and 123
OK；In another embodiment, sub-step 131 can be performed between sub-step 121 and 122.
In step 140, the second conversion can be the Fourier transformation along phase-encoding direction.
In one embodiment, the image rebuilding method 100 that magnetic resonance multilayer is imaged is performed simultaneously when, mediant therein
According to can be stored in memory, such as imaging data, coil merge coefficient simultaneously for the first individual layer calibration data, the first multilayer
It can be stored in one or more of monolithic layer data in memory.Because these data are the numbers in intermediate field
According to when being written and read to these data, being carried out in units of whole piece frequency coding line, in this way, data can be blocking
Ground processing, improve the read-write efficiency of memory.It would be appreciated that memory can be RAM memory, flash memory, ROM storages
Device, eprom memory, eeprom memory, register, hard disk, removable disk, CD-ROM or known in the art are any
The storage medium of other forms.
In a non-limiting embodiment, the part for the image rebuilding method 100 that magnetic resonance multilayer is imaged simultaneously calculates
It can be performed when gathering multilayer imaging data simultaneously.That is, multilayer imaging data simultaneously can be gathered with side, side is carried out
The data processing of the present invention.Relative to traditional method for needing to have gathered all data and just proceeding by processing, the present invention
It can shorten from the time required for gathered data to output final image.
In addition, it would be appreciated that, the method pre-processed in K spaces or the data in intermediate field can also be applied
Into the present invention, such as Bruder H, Fischer H, Reinfelder HE, Schmitt F.Image
reconstruction for echo planar imaging with nonequidistant k-space
sampling.Magn Reson Med.1992；23(2):311-23. disclosed Echo-Planar Imaging data calibrations
Deng.But this part is not the emphasis of the present invention, therefore this is not deployed to describe.
Fig. 9 is the magnetic resonance multilayer of one embodiment of the invention while the fundamental block diagram of the equipment for reconstructing image of imaging.Magnetic is total to
Multilayer of shaking at the same imaging equipment for reconstructing image 200 can include：
First conversion module 210, for the individual layer calibration data in K spaces and multilayer imaging data progress simultaneously first
Conversion, to obtain the first individual layer calibration data and the first multilayer while imaging data in intermediate field.Wherein, intermediate field is located at K
Between space and image area, the first conversion can be the Fourier transformation along frequency coding direction.
Coil merges coefficients calculation block 220, merges coefficient for calculating coil according to the first individual layer calibration data.
Data hierarchy module 230, for merging coefficient and the first multilayer according to coil, imaging data carries out data point simultaneously
Layer, to obtain the monolithic layer data in intermediate field.
In one embodiment, the equipment for reconstructing image 200 that magnetic resonance multilayer is imaged simultaneously can also include second and become mold changing
Block 240, for performing the second conversion to monolithic layer data, to obtain MRI.
In one embodiment, coil merges coefficients calculation block 220 and can closed by performing following steps to calculate coil
And coefficient：Apply first phase modulation to the first individual layer calibration data, to construct the first coefficient matrix；Adjusted according to phase is not applied
The first individual layer calibration data construction first object vector of system；And according to the first coefficient matrix and first object Vector operation line
Circle merges coefficient.
In one embodiment, data hierarchy module 230 can be carried out by performing following steps to aspect aliased data
Layering：To the first multilayer, imaging data applies second phase modulation simultaneously, to construct the second coefficient matrix；And closed according to coil
And coefficient and the second coefficient matrix calculate the monolithic layer data obtained in intermediate field.
In one embodiment, the equipment for reconstructing image 200 that magnetic resonance multilayer is imaged simultaneously can also include memory 250,
For storing intermediate data.Intermediate data can be that imaging data, coil merge simultaneously for the first individual layer calibration data, the first multilayer
One or more of coefficient and monolithic layer data.Because these data are the data in intermediate field, enter to these data
During row read-write, it can be carried out in units of whole piece frequency coding line, in this way, data can by blocking be handled, improve storage
The read-write efficiency of device.It would be appreciated that memory 250 can be RAM memory, flash memory, ROM memory, eprom memory,
The storage of eeprom memory, register, hard disk, removable disk, CD-ROM or any other form known in the art is situated between
Matter.
In a non-limiting embodiment, the equipment for reconstructing image 200 that magnetic resonance multilayer is imaged simultaneously can gather
The first conversion is carried out when multilayer while imaging data, calculates coil merging one or more of coefficient and data hierarchy.Also
It is to say, multilayer imaging data simultaneously can be gathered with side, side carries out the data processing of the present invention.Need to have gathered relative to traditional
All data can just proceed by the device of processing, and the present invention can shorten required for gathered data to output final image
Time.
Figure 10 is the result schematic diagram being layered to multilayer aliased image of one embodiment of the invention.It can be seen by Figure 10
Go out, original aliased image has obtained distinct layer 1 and layer 2 after being handled via the method, apparatus of the present invention.
Figure 11 is the structural representation of the magnetic resonance multilayer of one embodiment of the invention while the equipment for reconstructing image of imaging.Please
With reference to figure 11, the equipment for reconstructing image 300 that magnetic resonance multilayer is imaged simultaneously includes memory 310 and processor 320.Memory
Computer code is stored with 310, promotes device 300 when the computer code is configured as running on the processor 320
At least perform the image rebuilding method of magnetic resonance multilayer as described above while imaging.
Figure 12 is the schematic diagram of the computer-readable medium of one embodiment of the invention.Stored on computer-readable medium 400
There is computer code, magnetic resonance multilayer as described above is performed when the computer code is configured as running on a processor
The image rebuilding method of imaging simultaneously.
Those skilled in the art will further appreciate that, with reference to the embodiments described herein come the various illustratives that describe
Logic plate, module, circuit and algorithm steps can be realized as electronic hardware, computer software or combination of the two.To be clear
Explain to Chu this interchangeability of hardware and software, various illustrative components, frame, module, circuit and step be above with
Its functional form makees vague generalization description.Such feature be implemented as hardware or software depend on concrete application and
Put on the design constraint of total system.Technical staff can be realized described for every kind of application-specific with different modes
Feature, but such realize that decision-making should not be interpreted to cause departing from the scope of the present invention.
With reference to presently disclosed embodiment describe various illustrative logic modules and circuit can use general processor,
Digital signal processor (DSP), application specific integrated circuit (ASIC), field programmable gate array (FPGA) or other FPGAs
Device, discrete door or transistor logic, discrete nextport hardware component NextPort or its be designed to carry out any group of function described herein
Close to realize or perform.General processor can be microprocessor, but in alternative, the processor can be any routine
Processor, controller, microcontroller or state machine.Processor is also implemented as the combination of computing device, such as DSP
Combination, multi-microprocessor with microprocessor, one or more microprocessors to be cooperated with DSP core or any other this
Class configures.
It can be embodied directly in hardware, in by processor with reference to the step of method or algorithm that embodiment disclosed herein describes
Embodied in the software module of execution or in combination of the two.Software module can reside in RAM memory, flash memory, ROM and deposit
Reservoir, eprom memory, eeprom memory, register, hard disk, removable disk, CD-ROM or known in the art appoint
In the storage medium of what other forms.Exemplary storage medium is coupled to processor to enable the processor from/to the storage
Medium is read and write-in information.In alternative, storage medium can be integrated into processor.Processor and storage medium can
Reside in ASIC.ASIC can reside in user terminal.In alternative, processor and storage medium can be used as discrete sets
Part is resident in the user terminal.
In one or more exemplary embodiments, described function can be in hardware, software, firmware or its any combinations
Middle realization.If being embodied as computer program product in software, each function can be used as the instruction of one or more bars or generation
Code storage is transmitted on a computer-readable medium or by it.Computer-readable medium includes computer-readable storage medium and communication
Both media, it includes any medium for facilitating computer program to shift from one place to another.Storage medium can be can quilt
Any usable medium that computer accesses.It is non-limiting as example, such computer-readable medium may include RAM, ROM,
EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage apparatus can be used to carrying or store instruction
Or desirable program code and any other medium that can be accessed by a computer of data structure form.Any connection is also by by rights
Referred to as computer-readable medium.For example, if software is using coaxial cable, fiber optic cables, twisted-pair feeder, digital subscriber line
(DSL) or the wireless technology of such as infrared, radio and microwave etc passes from web site, server or other remote sources
Send, then the coaxial cable, fiber optic cables, twisted-pair feeder, DSL or such as infrared, radio and microwave etc is wireless
Technology is just included among the definition of medium.Disk (disk) and dish (disc) as used herein include compact disc
(CD), laser disc, laser disc, digital versatile disc (DVD), floppy disk and blu-ray disc, which disk (disk) are often reproduced in a manner of magnetic
Data, and dish (disc) laser reproduce data optically.Combinations of the above should also be included in computer-readable medium
In the range of.
Although the present invention describes with reference to current specific embodiment, those of ordinary skill in the art
It should be appreciated that the embodiment of the above is intended merely to the explanation present invention, can also make in the case of without departing from spirit of the invention
Go out various equivalent change or replacement, therefore, as long as change, change in the spirit of the present invention to above-described embodiment
Type will all fall in the range of following claims.
Claims (10)
1. a kind of MR imaging method, including：
To the individual layer calibration data in K spaces and multilayer, imaging data carries out the first conversion simultaneously, to obtain the in intermediate field
One individual layer calibration data and the first multilayer while imaging data, the intermediate field is between the K spaces and image area；
Coil is calculated according to the first individual layer calibration data and merges coefficient；
Merging coefficient and first multilayer according to the coil, imaging data carries out data hierarchy simultaneously, to obtain in intermediate field
Multiple monolithic layer datas；And
Multiple monolithic layer datas in intermediate field are carried out with the second conversion, to obtain MRI.
2. according to the method for claim 1, it is characterised in that coil is calculated according to the first individual layer calibration data and merged
Coefficient includes：
Apply first phase modulation to the first individual layer calibration data, to construct the first coefficient matrix；
The first individual layer calibration data according to phase-modulation is not applied constructs first object vector；And
Coefficient is merged according to coil described in first coefficient matrix and the first object Vector operation.
3. according to the method for claim 1, it is characterised in that described that coefficient and described more than first is merged according to the coil
The layer data hierarchy of imaging data progress simultaneously includes：
To first multilayer, imaging data applies second phase modulation simultaneously, to construct the second coefficient matrix；And
Coefficient is merged according to the coil and second coefficient matrix calculates the monolithic layer data obtained in intermediate field.
4. according to the method for claim 1, it is characterised in that described first is transformed to the Fourier along frequency coding direction
Conversion.
5. according to the method for claim 1, it is characterised in that the multilayer simultaneously imaging data by repeatedly exciting target
Region obtains, and the target area include multiple scanning slices, and the quantity of the scanning slice is more than exciting number.
6. according to the method for claim 5, it is characterised in that the target area is at least performed and excited for the first time and the
Secondary excitation, the first time excites corresponding first group of scanning slice, described to excite corresponding second group of scanning slice for the second time.
7. according to the method for claim 6, it is characterised in that first group of scanning slice or second group of scanning slice bag
Spaced multiple scanning slices are included, and during second group of scanning slice is excited, it is corresponding to first group of scanning slice
Imaging data carry out first conversion.
8. according to the method for claim 1, it is characterised in that described second is transformed to the Fourier along phase-encoding direction
Conversion.
9. a kind of MR imaging apparatus, including：
First conversion module, for being converted to the individual layer calibration data in K spaces and multilayer imaging data progress simultaneously first, with
Obtaining the first individual layer calibration data in intermediate field and the first multilayer, imaging data, the intermediate field are located at the K spaces simultaneously
Between image area；
Coil merges coefficients calculation block, merges coefficient for calculating coil according to the first individual layer calibration data；
Data hierarchy module, for merging coefficient and first multilayer according to the coil, imaging data carries out data point simultaneously
Layer, to obtain the monolithic layer data in intermediate field；And
Second conversion module, for multiple monolithic layer datas in intermediate field to be carried out with the second conversion, to obtain MRI.
10. a kind of MR imaging apparatus, including memory, processor and storage can transport on a memory and on a processor
Capable computer instruction, it is characterised in that realized described in the computing device during computer instruction as in claim 1 to 8
Method described in any one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710915393.2A CN107582057B (en) | 2017-09-30 | 2017-09-30 | Magnetic resonance imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710915393.2A CN107582057B (en) | 2017-09-30 | 2017-09-30 | Magnetic resonance imaging method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107582057A true CN107582057A (en) | 2018-01-16 |
| CN107582057B CN107582057B (en) | 2020-10-27 |
Family
ID=61052488
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710915393.2A Active CN107582057B (en) | 2017-09-30 | 2017-09-30 | Magnetic resonance imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107582057B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111103562A (en) * | 2019-12-09 | 2020-05-05 | 中国科学院深圳先进技术研究院 | Reconstruction method and device for simultaneously imaging multiple slice layers |
| CN114624640B (en) * | 2022-02-24 | 2024-05-17 | 深圳技术大学 | Magnetic resonance three-dimensional motion correction method and system based on simultaneous multi-layer acquisition |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100244825A1 (en) * | 2009-03-31 | 2010-09-30 | General Electric Company | System And Method Of Parallel Imaging With Calibration To A Virtual Coil |
| CN103901376A (en) * | 2012-12-30 | 2014-07-02 | 上海联影医疗科技有限公司 | Magnetic resonance imaging method and device |
| CN104248435A (en) * | 2013-06-27 | 2014-12-31 | 上海联影医疗科技有限公司 | Magnetic resonance imaging method and device |
| CN104635188A (en) * | 2013-11-12 | 2015-05-20 | 上海联影医疗科技有限公司 | K-space reconstruction method and magnetic resonance imaging method |
| CN104778341A (en) * | 2014-01-09 | 2015-07-15 | 上海联影医疗科技有限公司 | Magnetic resonance coil combination coefficient calculation method, magnetic resonance imaging method and device thereof |
| CN106597337A (en) * | 2016-12-09 | 2017-04-26 | 深圳先进技术研究院 | Magnetic resonance T2 * weighted rapid imaging method and device |
-
2017
- 2017-09-30 CN CN201710915393.2A patent/CN107582057B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100244825A1 (en) * | 2009-03-31 | 2010-09-30 | General Electric Company | System And Method Of Parallel Imaging With Calibration To A Virtual Coil |
| CN103901376A (en) * | 2012-12-30 | 2014-07-02 | 上海联影医疗科技有限公司 | Magnetic resonance imaging method and device |
| CN104248435A (en) * | 2013-06-27 | 2014-12-31 | 上海联影医疗科技有限公司 | Magnetic resonance imaging method and device |
| CN104635188A (en) * | 2013-11-12 | 2015-05-20 | 上海联影医疗科技有限公司 | K-space reconstruction method and magnetic resonance imaging method |
| CN104778341A (en) * | 2014-01-09 | 2015-07-15 | 上海联影医疗科技有限公司 | Magnetic resonance coil combination coefficient calculation method, magnetic resonance imaging method and device thereof |
| CN107092805A (en) * | 2014-01-09 | 2017-08-25 | 上海联影医疗科技有限公司 | magnetic resonance parallel imaging device |
| CN106597337A (en) * | 2016-12-09 | 2017-04-26 | 深圳先进技术研究院 | Magnetic resonance T2 * weighted rapid imaging method and device |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111103562A (en) * | 2019-12-09 | 2020-05-05 | 中国科学院深圳先进技术研究院 | Reconstruction method and device for simultaneously imaging multiple slice layers |
| CN111103562B (en) * | 2019-12-09 | 2022-03-04 | 中国科学院深圳先进技术研究院 | Reconstruction method and device for simultaneously imaging multiple slice layers |
| CN114624640B (en) * | 2022-02-24 | 2024-05-17 | 深圳技术大学 | Magnetic resonance three-dimensional motion correction method and system based on simultaneous multi-layer acquisition |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107582057B (en) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111513716B (en) | Method and system for magnetic resonance image reconstruction using an extended sensitivity model and a deep neural network | |
| CN111656392B (en) | System and method for synthesizing magnetic resonance images | |
| JP6270709B2 (en) | Method and apparatus for reconstructing raw magnetic resonance data | |
| CN104204838B (en) | Have Dixon type water/fat separate and with regard to main field inhomogeneities priori MRI | |
| Kasten et al. | Data-driven MRSI spectral localization via low-rank component analysis | |
| RU2702859C2 (en) | Parallel multi-slice mr-imaging with suppression of side band artefacts | |
| Chen et al. | Self‐calibrating wave‐encoded variable‐density single‐shot fast spin echo imaging | |
| JP6762284B2 (en) | Magnetic resonance imaging device and noise removal method | |
| US11860257B2 (en) | Method for acquiring a three-dimensional magnetic resonance image dataset and for generating a motion-corrected image dataset | |
| CN115639510A (en) | Magnetic resonance imaging method, spectroscopic imaging method, apparatus, device, and storage medium | |
| US10690741B2 (en) | Method and systems for reducing artifacts in magnetic resonance imaging | |
| US11307278B2 (en) | Reconstruction of MR image data | |
| US8823372B2 (en) | Method for homogenizing resolution in magnet resonance tomography measurements using non-linear encoding fields | |
| JP7350568B2 (en) | Magnetic resonance data acquisition method and magnetic resonance imaging apparatus | |
| CN107582057A (en) | MR imaging method and device | |
| JP5636058B2 (en) | Magnetic resonance imaging device | |
| US20190170842A1 (en) | Determining slice leakage in accelerated magnetic resonance imaging | |
| JP6618786B2 (en) | Magnetic resonance imaging apparatus and image processing apparatus | |
| CN110907871B (en) | Method and system for reducing shadows and blurring in magnetic resonance imaging | |
| CN115113114A (en) | System and method for deep learning based accelerated magnetic resonance imaging with extended field of view coil sensitivity calibration | |
| CN114062988B (en) | Magnetic resonance spectrum imaging method, apparatus, computer device and storage medium | |
| Chen | Fast Motion-robust Magnetic Resonance Imaging | |
| Daval‐Frérot et al. | Deep learning‐assisted model‐based off‐resonance correction for non‐Cartesian SWI | |
| Menon et al. | Data‐driven optimization of sampling patterns for MR brain T1ρ mapping | |
| Chappell et al. | Fourier Reconstruction—MRI |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |