Background technique
Magnetic resonance temperature real time monitoring is to realize the important prerequisite of guide of magnetic resonant image.The magnetic that can be used for temperature imaging is total
Vibration parameter has proton density, spin-lattice relaxation time, spin spin relaxation time, hydrone apparent diffusion coefficient, magnetization
Transfer, Proton Resonance Frequency displacement (Proton Resonance Frequency Shift, PRFS) and intermolecular quantum coherent
Deng.Temperature imaging based on Proton Resonance Frequency transfer (Proton Resonance Frequency Shift, PRFS) principle
Method is one of the method being most widely used in magnetic resonance temperature imaging, this is because this method has good time-space resolution
Rate has the fabulous linearity to temperature change, and does not change with the variation of tissue and environment with the linear relationship of temperature.
Hydrogen Proton variation of resonant frequency in hydrone is utilized in PRFS temperature imaging method can variation phenomenon, resonance with temperature change
The relationship of frequency variation and temperature change, which is about that temperature is every, increases 1 DEG C, and Hydrogen Proton variation of resonant frequency reduces 0.01ppm.When
The preceding magnetic resonance temperature measurement precision based on PRFS principle is mainly by temperature-insensitive ingredient (i.e. fatty), magnetic susceptibility artifact and movement
The influence of artifact.
It is existing to be had based on PRFS temp measuring method: have reference (reference) subtraction (Poorter, J.D., et al.,
Noninvasive MRI Thermometry with the Proton Resonance Frequency Method:In
Vivo Results in Human Muscle.Magnetic Resonance in Medicine,1995.33(1):p.74-
81), without with reference to (reference-less) thermometry (Rieke, V., et al., Referenceless PRF shift
Thermometry.Magn Reson Med, 2004.51 (6): p.1223-31.), fat with reference to (fat referenced) survey
Warm therapy (Hofstetter, L.W., et al., Fat-referenced MR thermometry in the breast and
Prostate using IDEAL.J Magn Reson Imaging, 2012.36 (3): p.722-32.) etc..
There is the method with reference to subtraction, subtracted each other using current phase diagram and reference phase.This method is mainly to histokinesis's ratio
It is more sensitive.During heating, the reference phase mismatch due to caused by histokinesis or deformation can generate biggish thermometric and miss
Difference.To overcome this disadvantage, there is researcher to propose without reference temp measuring method, wherein assuming phase diagram by the sky of one group of continuous and derivable
Between Polynomial approximation of functions, the phase diagram in unmanaged region is used for fit phase spatial model, and the background phase of heat affected zone
Position is obtained by the model extrapolation.Acquisition of this method independent of benchmark image, benchmark image pass through phase by present image
Spatial model is fitted and obtains.In practical application, in order to improve temperature measurement accuracy, whether has with reference to subtraction or surveyed without reference
Warm therapy requires to carry out fat suppression.And under the conditions of High-Field or super high field, the existing fat suppression in practical thermometric tissue
Technology may have in the signal temperature-resistant ingredient in local error, directly affect thermometric accuracy.In order to
The defect for overcoming Fat-suppression technique has researcher to propose the fat based on PRFS principle with reference to method.This method utilizes water rouge point
From technology, water and fat phase map are obtained respectively.Phase caused by main field is drifted about is corrected using fat phase variation to become
Change.The phase change and fat phase difference in change of front and back water are heated by calculating, and obtain temperature-induced phase change.The party
The major defect of method be it be a kind of modification reference subtraction, still can not eliminate influence of the histokinesis to temperature measurement accuracy.
Summary of the invention
Based on this, it is necessary to it needs to carry out fat suppression for the existing magnetic resonance temperature measurement method based on PRFS principle,
Under the conditions of High-Field or super high field, thermometric accuracy and existing base can be influenced in local error in practical thermometric tissue
Influence problem of the histokinesis to temperature measurement accuracy can not be eliminated with reference to method in the fat of PRFS principle, is provided a kind of based on water rouge point
From without reference temperature measurement method and system, not only without carrying out fat suppression but also insensitive to histokinesis.
It is a kind of based on water rouge separation without reference temperature measurement method comprising:
Obtain the more echo sequence images of two-dimensional gradient;
Area-of-interest and reference zone will be divided into as tissue in described image;
Based on the method for water rouge separation, the water of the reference zone and the phase diagram and the reference zone of fat are obtained
Field figure;
Water and the phase diagram of fat and the field figure of the reference zone to the reference zone carry out multinomial respectively
Models fitting acquires three groups of multinomial coefficients respectively;
Using three groups of multinomial coefficients, interpolation is carried out, obtains the water of the area-of-interest and the baseline phase of fat
The field figure of bitmap and the area-of-interest;
Based on the method for water rouge separation, using the more echo sequence image signal values of the two-dimensional gradient of the area-of-interest and
The field figure of the area-of-interest carries out the separation of water rouge to the area-of-interest in described image, obtains the area-of-interest
The phase diagram of water and fat;
According to the baseline phase figure and phase diagram of the water of the area-of-interest and fat, the area-of-interest is obtained
The phase change of water and fat;
The temperature change of area-of-interest is obtained according to the phase change of the water of the area-of-interest and fat.
The more echo sequences of the two-dimensional gradient generally acquire 3 to 5 echoes in one of the embodiments,.
The reference zone surrounds the area-of-interest in one of the embodiments,.
In one of the embodiments, the area-of-interest is circle.
The reference zone is annular shape in one of the embodiments,.
The field of the water of the reference zone, the phase diagram of fat and the reference zone in one of the embodiments,
Figure is separated from the more echo sequence images of two-dimensional gradient using IDEAL algorithm.
The method of the water rouge separation uses IDEAL algorithm in one of the embodiments,.
The phase change of the water according to the area-of-interest and fat obtains feeling emerging in one of the embodiments,
The process of the temperature change in interesting region the following steps are included:
Calculate the area-of-interest by temperature-induced phase change;
According to no reference temperature measurement method, the linear relationship of phase change and temperature change based on water proton in tissue
Calculate the temperature change of the area-of-interest.
The phase change of the water of the area-of-interest is by area-of-interest water phase position in one of the embodiments,
The baseline phase figure that figure subtracts the area-of-interest water acquires.
The phase change of the fat of the area-of-interest is by the area-of-interest fat in one of the embodiments,
The baseline phase figure that phase diagram subtracts the area-of-interest fat acquires.
In one of the embodiments, the area-of-interest by temperature-induced phase change by described interested
The difference that the phase change of regional water subtracts the phase change of the area-of-interest fat is calculated.
The temperature change of the area-of-interest and the area-of-interest are lured by temperature in one of the embodiments,
The phase change led is in a linear relationship.
It is a kind of based on water rouge separation without reference temperature measuring system, comprising:
Acquisition module, for obtaining the more echo sequence images of two-dimensional gradient;
Discriminating module, for area-of-interest and reference zone will to be divided into as tissue in described image;
Separation module, the method for being separated based on water rouge, obtain the reference zone water and fat phase diagram, with
And the field figure of the reference zone;
Fitting module, for the phase diagram of water and fat to the reference zone and the field figure point of the reference zone
Not carry out multinomial model fitting, acquire three groups of multinomial coefficients;
Extraction module, for utilize three groups of multinomial coefficients, carry out interpolation, obtain the area-of-interest water and
The baseline phase figure of fat and the field figure of the area-of-interest;
Processing module, the method for being separated based on water rouge utilize the more echo sequences of the two-dimensional gradient of the area-of-interest
The field figure of column image signal value and the area-of-interest carries out the separation of water rouge to the area-of-interest in described image, obtains institute
State the water of area-of-interest and the phase diagram of fat;
Module is obtained, for obtaining institute according to the water of the area-of-interest and the baseline phase figure and phase diagram of fat
State the water of area-of-interest and the phase change of fat;And
Computing module, for obtaining area-of-interest according to the water of the area-of-interest and the phase change of fat
Temperature change.
It is above-mentioned based on water rouge separation without reference temperature measurement method and system, water and rouge are obtained by water rouge isolation technics
The phase diagram of fat reduces pressure rouge technology bring error, and combines with no reference thermometric algorithm, and the method for the present invention was both right
Histokinesis is insensitive and without carrying out fat suppression, improves the precision of thermometric.
Specific embodiment
As depicted in figs. 1 and 2, it is a kind of based on water rouge separation without reference temperature measurement method, comprising the following steps:
In the step s 100, the more echo sequence images of two-dimensional gradient are obtained.
Preferably, the more echo sequences of above-mentioned two-dimensional gradient generally acquire 3 to 5 echoes.
It in step s 200, will be in above-mentioned image (referring to the more echo sequence images of two-dimensional gradient obtained in step S100)
It is divided into area-of-interest (ROI) and reference zone (ROR) to as tissue, as shown in Figure 2.Preferably, above-mentioned reference zone surrounds
The area-of-interest.Preferably, above-mentioned area-of-interest is circle.Preferably, above-mentioned reference zone is not comprising any temperature
Spend the region of variation.Preferably, above-mentioned reference zone is annular shape.
In step S300, based on the method for water rouge separation, the phase diagram of the water of above-mentioned reference zone, the phase of fat are obtained
The field figure of bitmap and above-mentioned reference zone.Here the preferred IDEAL algorithm of method of water rouge separation.Specifically, above-mentioned reference
The field figure of the phase diagram of the water in region, fatty phase diagram and above-mentioned reference zone is more from two-dimensional gradient using IDEAL algorithm
It is separated in echo sequence image, it is preferable that above-mentioned IDEAL algorithm uses T2* correction algorithm simultaneously.
In step S400, to the phase diagram of the water of above-mentioned reference zone, fatty phase diagram and above-mentioned reference zone
Field figure carry out respectively multinomial model (referring to following formula (1-1)) fitting, acquire three groups of multinomial coefficients.Specifically, more
The item for example following formula (1-1) of formula model, utilizing least square method to calculate least residual can get multinomial coefficient.
φ (x, y)=a0+a1x+a2y+a3xy+... (1-1)
In one particular embodiment of the present invention, above-mentioned steps 400 the following steps are included:
Multinomial model fitting is carried out to the phase diagram of water in above-mentioned reference zone, acquires first group of multinomial coefficient;
Multinomial model fitting is carried out to phase diagram fatty in above-mentioned reference zone, acquires second group of multinomial coefficient;
Multinomial model fitting is carried out to the field figure in above-mentioned reference zone, acquires third group multinomial coefficient.
In step S500, using three groups of above-mentioned multinomial coefficients, interpolation is carried out, the water of above-mentioned area-of-interest is obtained
Baseline phase and fat baseline phase and above-mentioned area-of-interest field distribution.Specifically, at of the invention one
In specific embodiment, above-mentioned steps 500 the following steps are included:
Using above-mentioned first group of multinomial coefficient, interpolation is carried out, the baseline phase of water in above-mentioned area-of-interest is obtained;
Using above-mentioned second group of multinomial coefficient, interpolation is carried out, baseline phase fatty in above-mentioned area-of-interest is obtained;
Using above-mentioned third group multinomial coefficient, interpolation is carried out, the field figure of above-mentioned area-of-interest is obtained.
In step S600, based on the method for water rouge separation, the more echo sequences of the two-dimensional gradient of above-mentioned area-of-interest are utilized
The field figure of column image signal value and above-mentioned area-of-interest carries out the separation of water rouge to the area-of-interest in above-mentioned image, in acquisition
State the water phase bitmap of area-of-interestAnd fat phase mapHere the preferred IDEAL algorithm of method of water rouge separation.
Specifically, the phase diagram of the water of above-mentioned area-of-interestThe phase diagram of fatAnd the field figure of above-mentioned area-of-interest
It is separated from the more echo sequence images of two-dimensional gradient using IDEAL algorithm, it is preferable that above-mentioned IDEAL algorithm is rectified using T2* simultaneously
Normal operation method.
Step S700 obtains above-mentioned sense according to the baseline phase figure and phase diagram of the water of above-mentioned area-of-interest and fat
The water in interest region and the phase change of fat.
Preferably, the phase change of the water of above-mentioned area-of-interestBy the phase diagram of the water of above-mentioned area-of-interestSubtract the baseline phase figure of the water of above-mentioned area-of-interestIt acquires, i.e., are as follows:
Preferably, the phase change of the fat of above-mentioned area-of-interestBy the phase of the fat of above-mentioned area-of-interest
FigureSubtract the baseline phase figure of the fat of above-mentioned area-of-interestIt acquires, i.e., are as follows:
Step S800 obtains the temperature of area-of-interest according to the phase change of the water of above-mentioned area-of-interest and fat
Variation.
In one embodiment of the invention, above-mentioned steps S800 the following steps are included:
Firstly, calculate above-mentioned area-of-interest by temperature-induced phase changeIn one of them of the invention
In embodiment, above-mentioned area-of-interest by temperature-induced phase changeBecome by the phase of above-mentioned area-of-interest water
The difference for changing the phase change for subtracting above-mentioned area-of-interest fat is calculated, i.e., are as follows:
Wherein, phase change caused by main field is drifted about is corrected with the phase change of fatIt is i.e. above-mentioned interested
The phase change of the fat in region
Then, according to no reference temperature measurement method, the line of phase change and temperature change based on water proton in tissue
Sexual intercourse calculates the change in temperature Δ T of area-of-interest, referring specifically to following formula (5).
Wherein,Above-mentioned area-of-interest by temperature-induced phase change, γ indicates gyromagnetic ratio, α=
0.01ppm/ DEG C, B0 represents main field strength, and TE represents the echo time.
Based on the above method, in one embodiment of the invention, as shown in figure 3, it is a kind of based on water rouge separation without ginseng
Examine temperature measurement system, comprising: acquisition module 10, discriminating module 20, separation module 30, fitting module 40, extraction module 50, place
It manages module 60, obtain module 70 and computing module 80.Wherein,
Acquisition module 10, for obtaining the more echo sequence images of two-dimensional gradient;
Discriminating module 20, for area-of-interest and reference zone will to be divided into as tissue in described image;
Separation module 30, the method for being separated based on water rouge, obtain the reference zone water and fat phase diagram,
And the field figure of the reference zone;
Fitting module 40, for the phase diagram of water and fat to the reference zone and the field figure of the reference zone
Multinomial model fitting is carried out respectively, acquires three groups of multinomial coefficients;
Extraction module 50 carries out interpolation, obtains the water of the area-of-interest for utilizing three groups of multinomial coefficients
With the baseline phase figure of fat and the field figure of the area-of-interest;
Processing module 60, the method for being separated based on water rouge utilize the more echoes of the two-dimensional gradient of the area-of-interest
The field figure of sequence image signal value and the area-of-interest carries out the separation of water rouge to the area-of-interest in described image, obtains
The water of the area-of-interest and the phase diagram of fat;
Module 70 is obtained, for obtaining according to the water of the area-of-interest and the baseline phase figure and phase diagram of fat
The water of the area-of-interest and the phase change of fat;And
Computing module 80, for obtaining area-of-interest according to the water of the area-of-interest and the phase change of fat
Temperature change.
In one embodiment of the invention, above-mentioned computing module 80 includes with lower unit
For calculating the unit by temperature-induced phase change of the area-of-interest;With
For according to no reference temperature measurement method, phase change and temperature change based on water proton in tissue it is linear
Relationship calculates the unit of the temperature change of the area-of-interest.
It is above-mentioned based on water rouge separation without reference temperature measurement method and system, water and rouge are obtained by water rouge isolation technics
The phase diagram of fat reduces pressure rouge technology bring error, and combines with no reference thermometric algorithm, and this method is both to tissue
Motion-insensitive is not necessarily to carry out fat suppression again, improves the precision of thermometric.
The embodiments described above only express several embodiments of the present invention, and the description thereof is more specific and detailed, but simultaneously
Limitations on the scope of the patent of the present invention therefore cannot be interpreted as.It should be pointed out that for those of ordinary skill in the art
For, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to guarantor of the invention
Protect range.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.