The invention relates to a nano test technology field, and more particularly to a method and a system for in-vivo temperature measurement based on AC magnetization of magnetic nanoparticles.
In-vivo temperature refers to temperature of tissues in complete and survival individuals. In the biological and medical field, in-vivo temperature measurement methods are divided into invasive measurement and noninvasive measurement. The invasive measurement method is simple, real-time and accurate, and it is convenient to monitor measurement positions by videos. However, problems with the method are that, wound caused by the method is comparatively large, insertion of probes often leads to metastasis of pathological cells, interaction between radiation fields of heating sources and the probes reduces measurement accuracy, and measured temperature is temperature at a certain point, not temperature distribution over a whole coagulation area. The noninvasive temperature measurement method is now widely accepted and used in the medical field as being capable of effectively preventing infection of wound or expansion of cancer cell, and of providing real-time information of in-vivo temperature and temperature (field) distribution images deeply in tissues with comparatively high accuracy.
Presently, the noninvasive temperature measurement method comprises a ultrasound temperature measurement method, a microwave temperature measurement method, a nuclear magnetic resonance (NMR) temperature measurement method, a remote magnetic nanoparticle measurement method and so on.
Limited by size of organisms, the ultrasound temperature measurement method employs reflection echo, and a disadvantages thereof is that it must measure sound and temperature characteristics of different tissues in advance, but temperature characteristics of the tissues is significantly different and unstable. The microwave temperature measurement method can only be used to tumor at superficial layers, and resolution is to be significantly reduced at deep layers. However, this method must predetermine temperature distribution and measure structure and electrical parameters of organisms, this is because that thermal noise microwave measured at body surface is related to temperature distribution, structure of tissues and electrical properties within a measurement range. Problems with the NMR temperature measurement method are that it is too expensive to be widely used, spatial and temperature resolution thereof are limited, and more importantly, it can only obtain variation in temperature, not absolute temperature of tissues. A remote noninvasive temperature measurement method using magnetic nanoparticles is a high-accuracy temperature (field) measurement technique deeply in tissues being capable of addressing the above-mentioned problems. This method implements measurement and control of in-vivo temperature, and thus being capable of assessing and adjusting hyperthermia treatment of tumor in real time at a molecular level. In addition, a magnetic particle imaging (MPI) based on magnetic nanoparticles features an extremely high real-time rate in remote noninvasive in-vivo temperature measurement. Superior to a traditional noninvasive method for in-vivo temperature measurement—NMR, the MPI system is more simple, cost-effective and accurate, and generates pixels by changing magnetization of the magnetic nanoparticles whereby enabling the system to have high temperature sensibility. Therefore, the MPI is capable of measuring and imaging the in-vivo temperature in theory.
Temperature-related attributes of the magnetic nanoparticle comprise a particle size, saturated magnetic moment, concentration distribution thereof in the organism, and so on. The particle size and the saturated magnetic moment can be determined via repeated test outside the organism, but it is difficult to measure concentration of the magnetic nanoparticle in the organism, and no measurement method is available at present. Meanwhile, uncertainty of concentration and spatial distribution of the magnetic nanoparticle in the organism is to cause a significant error in the in-vivo temperature measurement. To summarize, to facilitate remote temperature measurement without knowing concentration of the magnetic nanoparticle has become an unsolved problem in magnetic nano hyperthermia treatment of tumor.
In view of the above-mentioned problems, it is an objective of the invention to provide a method for in-vivo temperature measurement based on AC magnetization of magnetic nanoparticles that is capable of accurately measuring in-vivo temperature without obtaining concentration of the magnetic nanoparticles.
A method for in-vivo temperature measurement, comprising steps of:
(1) positioning magnetic nano agent at an area of a measured object;
(2) applying an AC excitation magnetic field to the area of the magnetic nano agent;
(3) collecting an AC magnetization signal of the magnetic nano agent under the AC excitation magnetic field;
(4) detecting amplitudes C1, C3, . . . , C2n-1 of odd harmonics of the AC magnetization signal, the number of odd harmonics n≧1; and
(5) calculating in-vivo temperature T according to a relationship between the odd harmonics and the in-vivo temperature X=AY, where
the amplitudes of the odd harmonics form a column vector
a coefficient matrix A is defined as: the AC magnetization signal is expressed via the Langevin's function, and expression of the amplitudes of the odd harmonics is deduced from a Taylor series expansion of the Langevin's function, and thus the coefficient matrix
where Ms is effective magnetic moment of a magnetic nanoparticle, k is a Boltzmann's constant, H0 is an amplitude of the AC excitation magnetic field, au,v is a coefficient of an element at the uth row and the vth column of the coefficient matrix A, u=1, 2, . . . , n, v=1, 2 . . . , m, m is the number of Taylor's expansion terms, and m≧n;
a related column vector Y of the in-vivo temperature is expressed as
where T is the in-vivo temperature of the measured object, and N is magnetic nano concentration at the measured object.
Advantageously, as the number of the Taylor's expansion terms m is equal to the number of the odd harmonics n, step (5) comprises: calculating the related column vector of the in-vivo temperature Y=A−1x, and further the in-vivo temperature T=2(q-p)√{square root over (Yp/Yq)}, 1≦p<q≦2n−1, where Yp and Yq respectively represents the pth and the qth element of the related column vector Y of the in-vivo temperature; and as the number of the Taylor's expansion terms m is greater than the number of the odd harmonics n, step (5) comprises: performing least square fitting via B=[N,T]T as a variable and the relationship between the odd harmonics and the in-vivo temperature X=AY as an objective function whereby determining the in-vivo temperature T.
Advantageously, step (4) uses a digital phase-sensitive detection method to detect the amplitudes of the odd harmonics, and comprises sub-steps of:
(a) sampling the AC magnetization signal M of the magnetic nano agent to obtain a discrete intensity signal M (i), where i=1, 2 . . . , L, and L is the number of sampling points;
(b) setting j=1;
(c) calculating f′=(2j−1)f, where f is a frequency of the AC excitation magnetic field;
(d) generating two orthogonal reference signals, and sampling the orthogonal reference signals to obtain two discrete sampling signals R1(i) and R2(i);
(e) determining whether j is equal to 1, proceeding to step (g) if yes, and proceeding to step (f) otherwise;
(f) filtering the 2j−3th odd harmonic in the discrete intensity signal M(i);
(g) performing mutual correlation calculation and mean filtering on the discrete sampling signal R1(i) and the discrete intensity signal M(i) to obtain a first component S2j-1, performing mutual correlation calculation and mean filtering on the discrete sampling signal R2(i) and the discrete intensity signal M(i) to obtain a second component O2j-1, and calculating an amplitude of a 2j−1th odd harmonic C2j-1=2√{square root over (S2j-12+O2j-12)}; and
(h) determining whether j is less than n, setting j=j+1 and retuning step (c) if yes, and ending otherwise.
Advantageously, step (4) uses a least square system parameter identification method to detect the amplitudes of the odd harmonics, and comprises sub-steps of:
sampling the AC magnetization signal M of the magnetic nano agent to obtain a discrete intensity signal M(i), where i=1, 2 . . . , L, and L is the number of sampling points;
assuming
finding matrices φ and Z so that J=[Z−φX]T[Z−φX] is the minimum by adjusting a sampling frequency, and determining the column vector X=(φTφ)−1φTZ formed by the amplitudes Cj of the odd harmonics, where a digital angular frequency
fs represents a sampling frequency, π represents a circumference ratio, and a superscript T represents transpose.
Advantageously, intensity of the AC excitation magnetic field ranges from 20 Oe to 1000 Oe, and a frequency of the AC excitation magnetic field ranges from 20 Hz to 20 kHz.
Advantageously, the number of the odd harmonics ranges from 2 to 4, and that of the Taylor's expansion terms m ranges from 2 to 6.
It is another objective of the invention to provide a method for in-vivo temperature measurement based on AC magnetization of magnetic nanoparticles that is capable of accurately measuring in-vivo temperature without obtaining concentration of the magnetic nanoparticles.
An in-vivo temperature measurement system, comprising:
a first module configured to position magnetic nano agent at an area of a measured object;
a second module configured to apply an AC excitation magnetic field to the area of the magnetic nano agent;
a third module configured to collect an AC magnetization signal of the magnetic nano agent under the AC excitation magnetic field;
a fourth module configured to detect amplitudes C1, C3, . . . , C2n-1 of odd harmonics of the AC magnetization signal, the number of odd harmonics n 1; and
a fifth module configured to calculate in-vivo temperature T according to a relationship between the odd harmonics and the in-vivo temperature X=AY, where
the amplitudes of the odd harmonics form a column vector
a coefficient matrix A is defined as: the AC magnetization signal is expressed via the Langevin's function, and expression of the amplitudes of the odd harmonics is deduced from a Taylor series expansion of the Langevin's function, and thus the coefficient matrix
where Ms is effective magnetic moment of a magnetic nanoparticle, k is a Boltzmann's constant, H0 is an amplitude of the AC excitation magnetic field, au,v is a coefficient of an element at the uth row and the vth column of the coefficient matrix A, u=1, 2, . . . , n, v=1, 2 . . . , m, m is the number of Taylor's expansion terms, and m≧n;
a related column vector Y of the in-vivo temperature is expressed as
where T is the in-vivo temperature of the measured object, and N is magnetic nano concentration at the measured object.
Advantageously, the fifth module comprises:
a fifty-first sub-module configured to calculate the related column vector of the in-vivo temperature Y=A−1 X, and further the in-vivo temperature T=2(q-p)√{square root over (Yp/Yq)}, 1≦p<q≦2n−1 as the number of the Taylor's expansion terms m is equal to the number of the odd harmonics n, where Yp and Yq respectively represents the pth and the qth element of the related column vector Y of the in-vivo temperature; and
a fifty-second sub-module configured to perform least square fitting via B=[N,T]T as a variable and the relationship between the odd harmonics and the in-vivo temperature X=AY as an objective function whereby determining the in-vivo temperature T as the number of the Taylor's expansion terms m is greater than the number of the odd harmonics n.
Advantageously, the fourth module uses a digital phase-sensitive detection method to detect the amplitudes of the odd harmonics, and comprises:
a forty-first sub-module configured to sample the AC magnetization signal M of the magnetic nano agent to obtain a discrete intensity signal M(i), where i=1, 2 . . . , L, and L is the number of sampling points;
a forty-second sub-module configured to set j=1;
a forty-third sub-module configured to calculate f′=(2j−1)f, where f is a frequency of the AC excitation magnetic field;
a forty-fourth sub-module configured to generate two orthogonal reference signals, and to sample the orthogonal reference signals to obtain two discrete sampling signals R1(i) and R2(i);
a forty-fifth sub-module configured to determine whether j is equal to 1, to proceed to step (g) if yes, and to proceed to step (f) otherwise;
a forty-sixth sub-module configured to filter the 2j−3th odd harmonic in the discrete intensity signal M(i);
a forty-seventh sub-module configured to perform mutual correlation calculation and mean filtering on the discrete sampling signal R1(i) and the discrete intensity signal M (i) to obtain a first component S2j-1, to perform mutual correlation calculation and mean filtering on the discrete sampling signal R2(i) and the discrete intensity signal M (i) to obtain a second component O2j-1, and to calculate an amplitude of a 2j−1th odd harmonic C2j-1=2√{square root over (S2j-12+O2j-12)}; and
a forty-eighth sub-module configured to determine whether j is less than n, to set j=j+1 and return step (c) if yes, and to end otherwise.
Advantageously, the fourth module uses a least square system parameter identification method to detect the amplitudes of the odd harmonics,
the fourth module is configured to sample the AC magnetization signal M of the magnetic nano agent to obtain a discrete intensity signal M(i), where i=1, 2 . . . , L, and L is the number of sampling points;
assuming
to find matrices φ and Z so that J=[Z−φX]T[Z−φX] is the minimum by adjusting a sampling frequency, and to determine the column vector X=(φTφ)−1φTZ formed by the amplitudes Cj of the odd harmonics, where a digital angular frequency
fs represents a sampling frequency, π represents a circumference ratio, and a superscript T represents transpose.
Advantages of the Invention Comprise:
The invention applies the AC excitation magnetic field to the field of the magnetic nano agent, collects the AC magnetization of the magnetic nano agent under the AC excitation magnetic field, detects the odd harmonics of the AC magnetization signal of the magnetic nanoparticles, analyzes spectrum information of a magnetic intensity model of the magnetic nanoparticles via the discrete Langevin's function and the Fourier transformation, establishes a matrix relationship between the odd harmonics and temperature and concentration of the magnetic nanoparticles, and solves the in-vivo temperature according to the relationship without considering concentration of the magnetic nanoparticles, and effective moment as temperature varies, and thus accurately detecting the in-vivo temperature. In addition, the invention is capable of accurately and rapidly detecting the in-vivo temperature, and is especially applicable to detect thermal motion of biological molecules, and an average error thereof after ten times of measurement experiments at a SNR of 80 dB is less than 0.2K.
For clear understanding of the objectives, features and advantages of the invention, detailed description of the invention will be given below in conjunction with accompanying drawings and specific embodiments. It should be noted that the embodiments are only meant to explain the invention, and not to limit the scope of the invention.
As shown in
(1) positioning magnetic nano agent at an area of a measured object;
in details, firstly surface of magnetic nanoparticles is covered with modification agent, and then the magnetic nanoparticles are redirected to the measured objects, such as various biological organs, along with blood circulation systems via injection.
(2) applying an AC excitation magnetic field to the area of the magnetic nano agent;
Using a Taylor series expansion of the Langevin's function of magnetization to calculate amplitudes of odd harmonics in the following steps is based on a weak excitation magnetic field, but a too weak excitation magnetic field is to lead to a small SNR of the magnetization, and make it difficult to extract high harmonics. Therefore, it is very important to select an appropriate strength for the AC excitation magnetic field, normally 20˜1000 Oe, and advantageously 100 Oe. In addition, a frequency of the AC excitation magnetic field is 20 Hz to 20 kHz, and advantageously 1 kHz.
(3) collecting an AC magnetization signal of the magnetic nano agent under the AC excitation magnetic field;
a detection coil as a sensor is placed in the vicinity of the measured object containing the magnetic nano agent, magnetization analog signal induced thereby is pre-processed by a conditioning circuit, and then sampled and transmitted to a computer for further procession via a data acquisition card.
(4) detecting amplitudes of odd harmonics of the AC magnetization signal;
firstly, the AC magnetization signal is modeled: AC magnetization of magnetic nanoparticles with paramagnetic characteristics follows the Langevin's function as follows:
where
is the Langevin's function, N is concentration of the magnetic nano agent, the excitation magnetic field H=H0 sin (2π ft), H0 is a amplitude of the excitation magnetic field, f is a frequency of the excitation magnetic field, Ms is effective magnetic moment of the magnetic nanoparticle, k is the Boltzmann's constant, T is temperature of the magnetic nanoparticle, namely in-vivo temperature of the measured object.
Then the above equation is divided into multiple frequency mixing signals containing multiple odd harmonics, namely
where C2j-1 is an amplitude of the 2j−1th odd harmonic, an angular frequency ω=2πf, t is time, n is the preset number of polynominal expansion terms that determines the number of all measured odd harmonics, if n is too small, an approximation error of finite terms of the AC magnetization and thus a temperature estimation error become large, if n is too large, it is difficult to accurately detect weak amplitude signals of high harmonics, and thus the temperature estimation error is to be increased. It is suggested that n range from 2 to 4.
By replacing the Fourier equation
with the Taylor series expansions
it is possible to obtain an expression for an amplitude C2j-1 of each harmonic, bm is a constant, the number of Taylor's expansion terms m ranges from 2 to 6, and m≧n
An example in which the number of odd harmonics n=4 and the number of Taylor's expansion terms m=4 will be given below:
the odd harmonics are proportional to concentration of the particles per unit volume, and harmonic ratios are functionally related to temperature and irrelevant to the concentration, and can prove approximation of the temperature without knowing particle concentration.
Based on the above model, this step employs a digital phase-sensitive detection (DPSD) method and a least square system parameter identification method to measure the amplitudes of the odd harmonics of different AC magnetization.
Measurement of an amplitude of a fundamental harmonic C1, an amplitude of a third harmonic C3, an amplitude of a fifth harmonic C5, and an amplitude of a seventh harmonic C7 of different AC magnetization using the DPSD and the least square system parameter identification method is described below with reference to an example in which the number of the harmonics n=4.
A. Using the DPSD to Measure the Amplitudes of the Odd Harmonics
where A is a amplitude of the particular frequency, θ is a phase angle of the particular frequency, and finally the amplitude and the phase angle can be obtained according to equations A=2√{square root over (S2+O2)} and
A1. Signal Sampling
The AC magnetization signal M of the magnetic nanoparticles is M(i)=C1 sin (Ωi)+C3 sin (3 Ωi)+C5 sin (5 Ωi)+C7 sin (7 Ωi)+ε(i) after analog-digital conversion, and the two reference signals are R1(i)=cos(Ω′i), R2(i)=sin(Ω′i) after analog-digital conversion, where f is an excitation frequency, fs is a sampling frequency, a digital angular frequency
sampling points i=1, 2 . . . , L, L is the length of discrete sampling series, ε(.) is an error function containing high harmonics.
A2. Extraction of Amplitudes of Harmonics
(1) To obtain an amplitude of the fundamental harmonic, assume f′=f and use the DPSD method to get:
then the amplitude of the fundamental harmonic is
C1=2√{square root over (S12+O12)}
(2) after the amplitude C1 of the fundamental harmonic is obtained, to obtain the amplitude C3 of the third harmonic with high accuracy, assume f′=3 f and use the DPSD method to process a magnetization signal M′ after filtering the fundamental harmonic and get:
then the amplitude of the third harmonic is
C3=2√{square root over (S32+O32)}
(3) after the amplitude C3 of the third harmonic is obtained, for the magnetization M″ after the third harmonic of M′ is filtered, assume f′=5 f and use the DPSD method to get:
then the amplitude of the fifth harmonic is
C5=2√{square root over (S52+O52)}
(4) after the amplitude C5 of the fifth harmonic is obtained, for the magnetization M′″ after the third harmonic of M″ is filtered, assume f′=7 f and use the DPSD method to get:
then the amplitude of the seventh harmonic is
C7=2√{square root over (S72+O72)}
By performing the above steps, it is possible to extract the amplitudes of the harmonics in the AC magnetization signal of the magnetic nanoparticles via the DPSD method.
B. using the least square system parameter identification method to measure harmonics
Signal sampling: The AC magnetization signal M of the magnetic nanoparticles is M(i)=C1 sin (Ωi)+C3 sin (3 Ωi)+C5 sin (5 Ωi)+C7 sin (7 Ωi)+ε(i) after analog-digital conversion, where f is an excitation frequency, fs is a sampling frequency, a digital angular frequency
sampling points i=1, 2 . . . , L, ε(.) is an error function containing high harmonics, and L is the length of discrete sampling series.
Then the above equation is further converted to a matrix relationship:
In the model, the least square system parameter identification method assesses fitting of measurement data based on estimation of a unknown parameter vector X being capable of minimizing square sum of errors, namely
To minimize J=[Z−φX]T[Z−φX], it is required to find appropriate matrices φ and Z by adjusting a sampling frequency. Assume an optimum estimation column vector {circumflex over (X)}=(φTφ)−1φTZ, let
namely J is the minimum. Thus, the amplitude of the fundamental harmonic C1, the amplitude of the third harmonic C3, the amplitude of the fifth harmonic C5, and the amplitude of the seventh harmonic C7 can be obtained from the optimum estimation column vector {circumflex over (X)}.
(5) calculating in-vivo temperature according to the matrix relationship between the amplitudes C2j-1 of the odd harmonics and the in-vivo temperature T;
The amplitudes C2j-1 of the odd harmonics and the in-vivo temperature T follow a matrix relationship X=AY, where
X is a column vector formed by the amplitudes C2j-1 of the odd harmonics, and
Y is a related column vector of the in-vivo temperature, and
the coefficient matrix A is defined as: the AC magnetization signal is expressed via the Langevin's function, and expression of the amplitudes of the odd harmonics is deduced from a Taylor series expansion of the Langevin's function, and thus the coefficient matrix
where au,v is a coefficient of an element at the uth row and the vth column of the coefficient matrix A, u=1, 2, . . . , n, v=1, 2 . . . , m, m≧n.
Using the above relationship between the harmonics and the temperature X=AY to solve the in-vivo temperature comprises the following two scenarios:
(a) as the number of the Taylor's expansion terms m is equal to the number of the odd harmonics n, step (5) comprises: calculating the related column vector of the in-vivo temperature Y=A−1 X, and further the in-vivo temperature T=2(q-p)√{square root over (Yp/Yq)}, 1≦p<q≦2n−1, where Yp and Yq respectively represents the pth and the qth element of the related column vector Y of the in-vivo temperature;
(b) as the number of the Taylor's expansion terms m is greater than the number of the odd harmonics n, step (5) comprises: performing least square fitting via B=[N,T]T as a variable and the relationship between the odd harmonics and the in-vivo temperature X=AY as an objective function whereby determining the in-vivo temperature T.
An example in which m=n=4 is described as follows:
is converted to a matrix relationship:
assume a related column vector of the temperature
the amplitudes of the odd harmonics form a column vector
coefficient matrix A is
then the related column vector of the in-vivo temperature is calculated Y=A−1 X, thus the temperature T can be measured by any two dimensional data (such as the pth dimension and the qth dimension) without knowing concentration of the magnetic nanoparticles, namely
where Yp and Yq respectively represents the pth and the qth element of the related column vector Y of the in-vivo temperature. Advantageously, as P=1, q=2, highly accurate results can be obtained.
To study effectiveness of the temperature measurement method and to optimize design, the example uses simulation data with noise to test the algorithm of the invention. During simulation and testing, effective moment Ms of the particles in the agent is 4×10−19 (the value is determined by the type of the agent). Since a quadrinomial approximation model of the Fourier series of the AC magnetization causes error effect, an AC excitation magnetic field with a small amplitude must be used. In this example, the amplitude of the AC excitation magnetic field H0=100 Oe. A noise model uses a function in the MATLAB software named ‘awgn’ to apply noise with a preset SNR to the AC magnetization of the sampling signal. A SNR during temperature measurement is set to 80 dB, a frequency of the excitation magnetic field is set to 1 kHz (low frequency point), 5 kHz (middle frequency point) and 10 kHz (high frequency point) for different purposes, and a temperature test range is between 300 K to 340 K.
To test accuracy of the temperature measurement method of the invention, the simulation uses the DPSD method and the least square system parameter identification method to test temperature estimation errors obtained from different harmonics in the AC magnetization of the magnetic nanoparticles. Meanwhile, to test repeatability of the method, an average temperature estimation error after ten times of repeated measurement is used.
While preferred embodiments of the invention have been described above, the invention is not limited to disclosure in the embodiments and the accompanying drawings. Any changes or modifications without departing from the spirit of the invention fall within the scope of the invention.
Number | Date | Country | Kind |
---|---|---|---|
2013100658967 | Mar 2013 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2013/077917 | 6/25/2013 | WO | 00 |