CN103892809B - 一种磁纳米温度成像方法及系统 - Google Patents

Abstract

本发明公开一种磁纳米温度成像方法，首先，对磁纳米粒子样品所在区域同时施加恒定直流磁场和交流磁场，采集磁纳米粒子的交流磁化强度信号，检测出各奇次谐波幅值；然后，将恒定直流梯度场替换为含梯度磁场的组合直流磁场，采集磁纳米粒子的交流磁化强度信号，检测出各奇次谐波幅值；计算两次谐波幅值差值；利用朗之万函数的泰勒级数展开建立奇次谐波差值与温度的关系式，求解关系式获得在体温度；最后，改变直流梯度场至下一位置，直到完成整个一维空间的温度测量。本发明对磁纳米粒子施加不同的激励磁场，从而一维空间的温度成像转变成了对每一个小区间的点温度求解，从而在不知磁纳米粒子浓度的情况下精密、快速地获得一维空间温度场。

Description

一种磁纳米温度成像方法及系统

技术领域

本发明涉及纳米测试技术领域，具体涉及一种磁纳米温度成像方法，更具体地说，是一种基于磁纳米粒子顺磁特性的一维在体温度成像方法。

背景技术

在体(invivo)温度成像是指进行于完整且存活的个体内的组织的温度成像。在肿瘤热疗等生物医疗领域内，由于活体内部的温度场分布信息难以准确获取，导致很多医疗手段不能有效地使用。目前，在体温度测量方式分为侵入式测量和非侵入式测量。侵入式测温方法简单，便于直接实时高精度地监控病灶温度。但创伤性较大，插针容易引起病变细胞的转移，加热源的辐射场直接与探针作用引起测量精度下降，测得的温度数据为点温度，而非整个病灶的温度场分布。而非侵入式温度测量能够有效避免创口感染或癌细胞扩散，同时可以提供较高精度的在体温度场实时成像，正是由于这些优点，该方式在生物医疗领域具有广阔的应用潜力。

非侵入式温度测量主要有红外测温、超声测温、核磁共振测温和磁纳米远程测温等。红外测温是根据被测物的红外辐射强度确定其温度，用于对不同温度物体的表面温度测量，不能测量组织深处温度场，且容易受到物体发射率、气雾的影响。超声测温的关键是超声波传播时间的精确测量，但必须预先测出各种组织的声特性及其温度特性，而各组织的温度特性存在较大差异且不稳定给温度测量带来较大影响。核磁共振测温的缺陷在于价格昂贵，不利于普及应用，且空间分辨率及温度分辨率有限。利用磁纳米粒子进行非侵入式温度场成像可克服上述缺点。此方法实现在体温度成像，可对肿瘤热疗过程进行实时监测以便及时做出有效调整。

此外，目前基于磁纳米粒子的非侵入式在体温度测量方法，只能实现点温度的测量，不能得到组织深处的温度场分布图，且温度测量精度受磁纳米粒子在组织深处的浓度分布影响。因此探索一种在不知磁纳米粒子浓度分布的情况下实现在体温度场成像的方法成为磁纳米肿瘤热疗领域亟待解决的问题。

发明内容

针对现有技术的缺陷，本发明的目的在于提供一种基于磁纳米粒子顺磁特性的在体温度成像方法，旨在实现不知磁纳米粒子浓度分布的情况下准确地检测到一维在体温度场。

一种磁纳米温度成像方法，包括如下步骤：

（1）将磁纳米试剂置放于一维待测空间；

（2）向磁纳米试剂所在一维空间同时施加恒定直流磁场H_dc＝b和交流激励磁场，采集待测空间磁纳米试剂的交流磁化强度信号，获取该交流磁化强度信号的各奇次谐波幅值A₁,A₃…A_2j-1，谐波个数j≥2，b为恒定直流磁场的幅值；

（3）保持步骤（2）中的交流激励磁场不变，将步骤（2）中的恒定直流磁场替换为组合直流磁场 ${H_{\mathrm{dc}}}^{\&prime;}=\left\{\begin{array}{cc}b& x<x_1\\ f\left(x\right)& x_1\&le;x\&le;x_1+\mathrm{\&Delta;x}\\ -b& x>x_1+\mathrm{\&Delta;x}\end{array}\right.,$ 其中x表示位置变量，x₁表示直流梯度磁场f(x)相对一维待测空间的起始位置，Δx为直流梯度磁场f(x)的宽度，再次采集待测空间磁纳米试剂的交流磁化强度信号，获取该交流磁化强度信号的各奇次谐波幅值B₁,B₃…B_2j-1；

（4）计算步骤（3）中各奇次谐波幅值B₁,B₃…B_2j-1与步骤（2）中各奇次谐波幅值A₁,A₃…A_2j-1之间的差值S₁,S₃…S_2j-1；

（5）根据各奇次谐波幅值差与在体温度的关系式 $\begin{array}{c}{S}_1=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_1(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_1(y,z){|}_{z=b}\\ S_3=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_3(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_3(y,z){|}_{\text{z=b}}\\ \begin{array}{ccc}.& .& .\\ .& .& .\\ .& .& .\end{array}\\ S_{2j-1}=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_{2j-1}(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_{2j-1}(y,z){|}_{z=b}\end{array},$ 求解[x₁,x₁+Δx]区间的y，进而计算得到[x₁,x₁+Δx]区间的在体温度

其中，c＝NM_s，N为磁纳米粒子浓度，M_s为磁纳米样品原子有效磁矩，k₀为波尔兹曼常数，磁纳米粒子的磁化强度由朗之万函数描述，c·f_2j-1(y,z(r_k))为朗之万函数的有限项泰勒级数展开得到的第^2j-1次谐波幅值表达式，z(r_k)为[x₁,x₁+Δx]区间内的第k个离散点r_k对应的直流磁场大小，m为[x₁,x₁+Δx]区间被离散化的点数；

（6）改变直流梯度磁场f(x)的起始位置使得宽度为Δx的直流梯度磁场f(x)扫描至下一个区间，重复步骤（3）～（5）得到下一个区间的在体温度，按照此方式重复测试直到一维待测空间温度检测完毕。

进一步地，通过移动直流磁场产生装置改变直流梯度磁场相对一维待测空间的起始位置或者通过改变激励线圈的电流改变直流梯度磁场相对一维待测空间的位置。

进一步地，所述步骤（2）和（3）采用数字相敏检波方法或最小二乘系统参数辨识方法检测磁纳米粒子交流磁化强度信号的各奇次谐波幅值。

一种磁纳米温度成像系统，包括：

恒定直流磁场产生装置，用于向磁纳米试剂所在一维空间施加恒定直流磁场H_dc＝b，其中，b为恒定直流磁场的幅值；

组合直流磁场产生装置，用于向磁纳米试剂所在一维空间施加组合直流磁场 ${H_{\mathrm{dc}}}^{\&prime;}=\left\{\begin{array}{cc}b& x<x_1\\ f\left(x\right)& x_1\&le;x\&le;x_1+\mathrm{\&Delta;x}\\ -b& x>x_1+\mathrm{\&Delta;x}\end{array}\right.,$ 其中，x表示位置变量，x₁表示直流梯度磁场f(x)相对一维待测空间的起始位置，Δx为直流梯度磁场f(x)的宽度；

交流激励磁场产生装置，用于向磁纳米试剂所在一维空间施加交流激励磁场；

磁化强度采集装置，用于采集待测空间磁纳米试剂的交流磁化强度信号；

处理器，用于对同时施加恒定直流磁场和交流激励磁场时采集的交流磁化强度信号进行处理得到该交流磁化强度信号的各奇次谐波幅值A₁,A₃…A_2j-1，对在同时施加组合直流磁场和交流激励磁场时采集的交流磁化强度信号进行处理得到该交流磁化强度信号的各奇次谐波幅值B₁,B₃…B_2j-1，计算各奇次谐波幅值B₁,B₃…B_2j-1与A₁,A₃…A_2j-1之间的差值S₁,S₃…S_2j-1；根据各奇次谐波幅值差与在体温度的关系式 $\begin{array}{c}{S}_1=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_1(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_1(y,z){|}_{z=b}\\ S_3=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_3(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_3(y,z){|}_{\text{z=b}}\\ \begin{array}{ccc}.& .& .\\ .& .& .\\ .& .& .\end{array}\\ S_{2j-1}=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_{2j-1}(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_{2j-1}(y,z){|}_{z=b}\end{array},$ 求解[x₁,x₁+Δx]区间的y，进而计算得到[x₁,x₁+Δx]区间的在体温度

其中，c＝NM_s，N为磁纳米粒子浓度，M_s为磁纳米样品原子有效磁矩，k₀为波尔兹曼常数，磁纳米粒子的磁化强度由朗之万函数描述，c·f_2j-1(y,z(r_k))为朗之万函数的有限项泰勒级数展开得到的第^2j-1次谐波幅值表达式，z(r_k)为[x₁,x₁+Δx]区间内的第k个离散点r_k对应的直流磁场大小，m为[x₁,x₁+Δx]区间被离散化的点数，谐波个数j≥2。

本发明的技术效果体现在：

本发明首先对磁纳米粒子样品所在区域同时施加一个恒定直流磁场和一个交流磁场，采集磁纳米粒子的交流磁化强度信号，检测出交流磁化强度信号中各次谐波幅值；然后，对磁纳米粒子样品所在区域同时施加一个组合直流磁场和一个交流磁场，采集磁纳米粒子的交流磁化强度信号，检测出交流磁化强度信号中各次谐波幅值；接着，用第二次检测出的各次谐波幅值中的奇次谐波幅值减去第一次检测出的各次谐波幅值中的奇次谐波幅值得到谐波幅值差值。在两次不同磁场的激励下，整个待测区域内的奇次谐波幅值的差值只包含了所处位置Δx区域内磁纳米粒子的信号，即只与Δx区域的磁纳米粒子的温度和浓度有关。因此，整个一维待测区域内的磁纳米粒子的浓度分布对温度求解不产生影响，只有在Δx范围内的磁纳米粒子的浓度分布会影响温度求解精度。由于Δx很小，在Δx内的磁纳米粒子的浓度可以看成一个常数。这样一来，对磁纳米粒子施加不同的激励磁场，在温度成像的过程中可把整个一维空间分割成一个一个的小区间Δx，从而一维空间的温度成像转变成了对每一个小区间的点温度求解。本发明能够在不知磁纳米粒子浓度的情况下精密、快速地获得一维空间温度场，特别适用于生物分子层面热运动的温度成像，试验表明在信噪比为80dB的噪声环境下测量误差小于0.79K。

附图说明：

图1为本发明方法流程图；

图2为同时施加恒定直流磁场和交流激励磁场的示意图；

图3为同时施加组合直流磁场和交流激励磁场的示意图；

图4为一次、二次和三次谐波随直流磁场变化的示意图；

图5为总磁场为H(t)＝H₀sin(2πft)+H_dc时，待测区域内各处磁纳米粒子一次谐波大小的示意图；

图6为总磁场为H(t)＝H₀sin(2πft)+H_dc′时，待测区域内各处磁纳米粒子一次谐波大小的示意图；

图7为激励频率160Hz，信噪比90dB，温度测试范围308K至318K时，用DPSD算法测得谐波幅值的仿真结果，其中，图7（a）为一维温度成像，图7（b）为温度成像误差；

图8为激励频率160Hz，信噪比80dB，温度测试范围308K至318K时，用DPSD算法测得谐波幅值的仿真结果，其中，图8（a）为一维温度成像，图8（b）为温度成像误差；

图9为激励频率160Hz，信噪比90dB，温度测试范围300K至310K时，用DPSD算法测得谐波幅值的仿真结果，其中，图9（a）为一维温度成像，图9（b）为温度成像误差；

图10为激励频率160Hz，信噪比80dB，温度测试范围300K至310K时，用DPSD算法测得谐波幅值的仿真结果，其中，图10（a）为一维温度成像，图10（b）为温度成像误差；

图11为激励频率160Hz，信噪比80dB，磁场强度存在2%的偏差，温度测试范围300K至310K时，用DPSD算法测得谐波幅值的仿真结果，其中，图11（a）为一维温度成像，图11（b）为温度成像误差；

图12为激励频率160Hz，信噪比80dB，磁纳米粒子浓度分布存在10%的偏差，温度测试范围300K至310K时，用DPSD算法测得谐波幅值的仿真结果，其中，图12（a）为一维温度场成像，图12（b）为温度成像误差；

图13为直流梯度磁场为非线性时的组合直流磁场示意图；

图14为激励频率1.6kHz，信噪比90dB，温度测试范围300K至310K时，用DPSD算法测得谐波幅值的仿真结果，其中，图14（a）为一维温度成像，图14（b）为温度成像误差；

图15为激励频率1.6kHz，信噪比80dB，温度测试范围300K至310K时，用DPSD算法测得谐波幅值的仿真结果，其中，图15（a）为一维温度成像，图15（b）为温度成像误差。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

参见图1，本发明一种磁纳米温度成像方法，包括如下步骤：

（1）将磁纳米试剂置放于待测对象处。

将已包裹表面生物分子修饰剂的磁纳米粒子注射到人或动物体内，使其跟随血液循环系统靶向至待测对象处，例如人或动物体内的癌变组织。

（2）首次奇次谐波幅值采集

（21）向磁纳米试剂所在区域同时施加恒定直流磁场H_dc＝b和交流激励磁场H(t)＝H₀sin(2πft)，则总磁场为H(t)＝H₀sin(2πft)+H_dc，H₀表示交流磁场幅值，f表示频率，t表示时间，如图2所示。

由于后面步骤在求解在体温度时，只用到了朗之万函数的有限项泰勒级数展开式，因此考虑到模型的截断误差，交流磁场幅值H₀和直流磁场幅值b都应该尽量小一些。但是若施加的磁场过于微弱，背景噪声会增大，信号信噪比会降低，不便于有用信号的提取。因此，合理地选择交流磁场和直流磁场的强度至关重要，幅值具体大小可根据实验结果进行调整。

（22）采集待测区域磁纳米试剂的交流磁化强度信号。

利用螺线管或者小线圈作为传感器，探测待测区域内的磁纳米粒子的交流磁化强度，经过放大等调理电路后被数据采集卡采集并存储于计算机以便后续的数据处理。

（23）检测交流磁化强度信号的各奇次谐波幅值A₁,A₃…A_2j-1，谐波个数j≥2。

利用数字相敏检波方法（DPSD）或最小二乘系统参数辨识方法检测由数据采集卡采集的数据，并得到交流磁化强度信号的各次谐波幅值。

（3）再次奇次谐波幅值采集

（31）去除恒定直流磁场，向磁纳米试剂所在区域施加一个组合直流磁场 ${H_{\mathrm{dc}}}^{\&prime;}=\left\{\begin{array}{cc}b& x<x_1\\ f\left(x\right)& x_1\&le;x\&le;x_1+\mathrm{\&Delta;x}\\ -b& x>x_1+\mathrm{\&Delta;x}\end{array}\right.,$ 其中x表示位置变量，x₁表示直流梯度磁场f(x)相对一维待测空间的起始位置，Δx为直流梯度磁场f(x)的宽度，在Δx内的直流磁场z形成一个线性或非线性的梯度磁场，梯度磁场两侧分别为z＝b和z＝-b，同时保持交流激励磁场H(t)＝H₀sin(2πft)不变，则总磁场为H(t)＝H₀sin(2πft)+H_dc′，如图3所示。

必须说明此时施加的直流梯度磁场的两侧磁场强度必须与步骤（2）中的强度相同。同时，交流磁场的幅值和频率也应保持不变。此外，若磁场强度b过大，会导致位于梯度磁场内的磁纳米粒子的奇次谐波的正负符号发生变化，给在体温度的计算带来误差。

（32）采集待测区域磁纳米试剂的交流磁化强度信号。

与步骤（22）一致，利用螺线管或者小线圈作为传感器，探测待测区域内的磁纳米粒子的交流磁化强度，经过放大等调理电路后被数据采集卡采集并存储于计算机以便后续的数据处理。

（33）检测交流磁化强度信号的各奇次谐波幅值B₁,B₃…B_2j-1，谐波个数j≥2。

与步骤（23）一致，利用数字相敏检波方法（DPSD）或最小二乘系统参数辨识方法检测由数据采集卡采集的数据，并得到交流磁化强度信号的各次谐波幅值。

（4）计算步骤（3）中各奇次谐波幅值B₁,B₃…B_2j-1与步骤（2）中各奇次谐波幅值A₁,A₃…A_2j-1之间的差值S₁,S₃…S_2j-1；

（5）根据各奇次谐波幅值的差值与在体温度的关系式，计算在体温度。

磁纳米粒子的磁化强度可以由朗之万函数描述为：

$M={\mathrm{NM}}_s(\coth \left(\mathrm{\&alpha;}\right)-\frac{1}{\mathrm{\&alpha;}})$

其中，H为对磁纳米粒子施加的磁场，N为磁纳米粒子浓度，M_s为磁纳米粒子原子有效磁矩，k₀为波尔兹曼常数，T为待测对象的绝对温度。

当同时对磁纳米粒子施加交流磁场和直流磁场H(t)＝H₀sin(2πft)+z时，各次谐波幅值的表达式可由朗之万函数的有限项泰勒级数展开式近似获得，一次谐波幅值表达式为：

${\mathrm{Amp}}_1=c[\frac{1}{3}{H}_{0}y-(\frac{1}{60}{H}_0^3+\frac{1}{15}{H}_0{z}^2)y^3+(\frac{1}{756}{H_0}^5+\frac{1}{63}{H}_0^3{z}^2+H_0{z}^4)y^5+...],$

简记Amp₁＝c·f₁(y,z)；

二次谐波幅值表达式为：

${\mathrm{Amp}}_2=c[\frac{1}{30}{H}_0^2{\mathrm{zy}}^3-(\frac{1}{189}{H}_0^{4}z+\frac{1}{189}{H}_0^2{z}^3+H_0^2{z}^3)y^5+...],$ 简记Amp₂＝c·f₂(y,z)；

三次谐波幅值表达式为：

${\mathrm{Amp}}_3=c[\frac{1}{180}{H}_0^3{y}^3-(\frac{1}{1512}{H}_0^5+\frac{1}{189}{H}_0^3{z}^2)y^5+...],$ 简记Amp₃＝c·f₃(y,z)；

n次谐波幅值表达式简记Amp_n＝c·f_n(y,z)。

其中，c＝NM_s，N为磁纳米粒子浓度，M_s为样品原子有效磁矩，k₀为波尔兹曼常数，T为待测对象的绝对温度，z为直流磁场大小，H₀为交流磁场幅值。

通过数学归纳可知各奇次谐波是关于直流磁场z的偶函数，各偶次谐波是关于直流磁场z的奇函数，即f_2j-1(z)＝f_2j-1(-z)，f_2j(-z)＝-f_2j(z)，j≥1。如图4所示，可看出一次、三次谐波是直流磁场z的偶函数，二次谐波是直流磁场z的奇函数。因此，当施加的磁场从H(t)＝H₀sin(2πft)+H_dc变为H(t)＝H₀sin(2πft)+H_dc′时，在直流磁场z＝±b处的磁纳米粒子的奇次谐波的大小和符号保持不变，偶次谐波则大小不变符号相反。奇次谐波唯一发生变化的磁纳米粒子在区域Δx内，即直流磁场z形成的一个线性或非线性的梯度磁场内。

为了更清楚明白地表达本发明的思想，下面结合图5和图6详细分析谐波幅值大小随施加的激励磁场的变化情况，并推导出奇次谐波幅值的差值与在体温度的关系式。由于其他较高次的奇次谐波与一次谐波的分析方法相同，因此本发明只对一次谐波进行详细分析。

将待测区域看作有T₁,T₂…T_i…T_n个温度点，且每个温度点的空间范围与组合直流磁场中梯度磁场宽度Δx相同，由于Δx很小，在Δx内的磁纳米粒子的浓度可以认为是一个常数。当在待测区域施加磁场H(t)＝H₀sin(2πft)+H_dc时，由朗之万函数可知在同一温度T_i的Δx内的磁纳米粒子的一次谐波幅值相同，如图5所示。因此，整个一维空间内磁纳米粒子的一次谐波幅值为各个温度点处的磁纳米粒子一次谐波幅值的总和，记为A₁。

当在待测区域施加磁场H(t)＝H₀sin(2πft)+H_dc′，宽度为Δx的梯度磁场位于温度点T_i处，由于一次谐波的大小是关于直流磁场z的偶函数，即f₁(z)＝f₁(-z)，则在直流磁场z＝±b处的磁纳米粒子的一次谐波的大小和符号保持不变，唯一发生变化的磁纳米粒子位于区域为直流磁场z形成一个线性或非线性的梯度磁场内，此时各温度点处磁纳米粒子的一次谐波幅值如图6所示。因此，整个一维空间内磁纳米粒子的一次谐波幅值的总和为A₁与梯度磁场内的变化量之和，记为B₁。

在两次不同磁场的激励下，整个待测区域内的一次谐波幅值的差值可以表示为S_1,i＝B₁-A₁，结合图5和图6，不难看出S_1,i实际上表示在第i个温度点所处位置Δx区域内的磁纳米粒子在两次不同磁场的激励下的一次谐波幅值的差值，其只包含了第i个温度点所处位置Δx区域内磁纳米粒子的信号，即只与第i个温度点处的磁纳米粒子的温度和浓度有关。因此，不难发现整个一维待测区域内的磁纳米粒子的浓度分布对温度求解不产生影响，只有在Δx范围内的磁纳米粒子的浓度分布会影响温度的求解精度。此外，由于Δx很小，在Δx内的磁纳米粒子的浓度可以看成一个常数。

基于上述分析，可知，[x₁,x₁+Δx]区间的一次谐波幅值的差值可表示为

$S_1=\underset{\mathrm{\&Omega;}}{\&Integral;}c\&CenterDot;f_1(y,z\left(r\right))\mathrm{dr}-\underset{\mathrm{\&Omega;}}{\&Integral;}c\&CenterDot;f_1(y,z){|}_{z=b}\mathrm{dr}$

其中，c＝NM_s，N为磁纳米粒子浓度，M_s为磁纳米样品原子有效磁矩，c·f₁(y,z(r))表示Ω＝[x₁,x₁+Δx]区间内r∈Ω处的一次谐波大小。

同理，[x₁,x₁+Δx]区间的三次谐波幅值的差值可表示为

$S_3=\underset{\mathrm{\&Omega;}}{\&Integral;}c\&CenterDot;f_3(y,z\left(r\right))\mathrm{dr}-\underset{\mathrm{\&Omega;}}{\&Integral;}c\&CenterDot;f_3(y,z){|}_{z=b}\mathrm{dr}$

其中，c·f₃(y,z(r))表示Ω＝[x₁,x₁+Δx]区间内r∈Ω处的三次谐波大小。

同理，可得整个待测区域内的2j-1次谐波幅值的差值可表示为

$S_{2j-1}=\underset{\mathrm{\&Omega;}}{\&Integral;}c\&CenterDot;f_{2j-1}(y,z\left(r\right))\mathrm{dr}-\underset{\mathrm{\&Omega;}}{\&Integral;}c\&CenterDot;f_{2j-1}(y,z){|}_{z=b}\mathrm{dr}$

其中，c·f_2j-1(y,z(r))表示Ω＝[x₁,x₁+Δx]区间内r∈Ω处的2j-1次谐波大小，j≥2。

综上所述，磁纳米粒子各奇次谐波幅值的差值与在体温度之间的关系式得以建立。

将第i个温度点在Δx区域内r∈Ω离散化为m段，则各奇次谐波幅值的差值S_1,i…S_2j-1,i与在体温度之间的关系式构成如下方程组为

$\begin{array}{c}{S}_1=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_1(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_1(y,z){|}_{z=b}\\ S_3=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_3(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_3(y,z){|}_{\text{z=b}}\\ \begin{array}{ccc}.& .& .\\ .& .& .\\ .& .& .\end{array}\\ S_{2j-1}=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_{2j-1}(y,z\left(r_k\right))-\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}c\&CenterDot;f_{2j-1}(y,z){|}_{z=b}\end{array}$

对离散化的方程组求解即可得到第i个温度点处的c_i、y_i，从而获得第i个温度点的在体温度磁纳米粒子浓度

当检测的奇次谐波个数j＝2，可直接求解方程组获得c_i、y_i，从而得到第i个温度点的在体温度

当检测的奇次谐波个数j≥3，可采用最小二乘法等算法求解超定方程组获得c_i、y_i，从而得到在体温度

（6）采用机械扫描或激励磁场扫描的方式，改变直流梯度磁场相对一维待测空间的位置x₁，使宽度为Δx的梯度场扫描下一个位置，并记录空间坐标。若待测区域未全部检测完，则跳回步骤（2）。若待测区域已全部检测，则结束。

当使用机械扫描方式时，可利用电机来移动磁场装置使梯度磁场扫描整个待测区域或直接移动磁纳米粒子样品，并记录坐标位置信息。

当使用激励磁场扫描方式时，可通过改变激励线圈的电流等方式来控制梯度磁场的位置，使其扫描整个待测区域，并记录坐标位置信息。

当扫描速度足够快或温度变化缓慢时，第一次测得的A₁基本上不变，则可直接跳回步骤（3）完成整个成像区域的扫描，不必跳回步骤（2）。将每次扫描得到的信号直接与第一次测得的A₁相减即可，这样不仅不会影响温度求解精度还可以缩短温度成像时间；

当一维空间中磁纳米粒子分布较均匀时，可利用二次谐波幅值的差值来定位各个温度点的坐标位置。画出二次谐波幅值差值的变化趋势，可发现每个温度点的二次谐波幅值差值在坐标轴上的投影便是该温度点的坐标位置。

仿真实例1：

1.仿真模型与测试说明：

为了研究磁纳米温度成像方法的有效性和可行性，本仿真实例将采用含有噪声的仿真数据来测试该方法。仿真测试过程中采用的试剂粒子有效磁矩M_s测定为8.5×10^-19(补充说明有效磁矩的测定数值由试剂类型参数决定)，直流梯度磁场的宽度Δx＝5mm，模拟机械移动的方式每次移动1mm。考虑到采用朗之万函数的前八项泰勒级数展开式获得的近似模型的截断误差，本仿真实例采用交流磁场幅值H₀＝50Oe，频率160Hz，直流磁场强度为30Oe，并采用MATLAB中awgn函数给交流磁化强度信号添加噪声。为了达到不同的测试目的：

I磁纳米粒子浓度分布均匀时，在不同温度场分布的在体温度场成像过程中添加信噪比为90dB、80dB的噪声，温度测试范围分别是308K至318K和300K至310K，仿真结果如图7至图10所示。

II直流磁场幅值存在2%的偏差时，在温度成像过程中添加信噪比为80dB的噪声，温度测试范围300K至310K，仿真结果如图11所示。

III磁纳米粒子浓度存在10%的偏差时，在温度成像过程中添加信噪比为80dB的噪声，温度测试范围300K至310K，仿真结果如图12所示。

2、仿真试验结果：

图7至图10反映了在不同的温度场分布下，磁纳米温度成像方法的温度测量误差小于0.63K。由于测试的信噪比分别为80dB、90dB，可见该方法具有一定的抗噪能力。

图11反映了磁场幅值出现偏差时，一维在体温度场成像方法的精度和适用性。在信噪比为80dB的条件下，温度测量误差小于0.49K。表明在磁场幅值有所变化时，该方法仍可以得到较好温度求解精度。

图12反映了磁纳米粒子浓度分布不均匀时，一维在体温度场成像方法的精度和适用性。在信噪比为80dB的条件下，温度测量误差小于0.79K。表明在不知浓度且浓度分布不均匀的条件下，该方法仍可以得到较好的温度求解精度。

仿真实例2：

1.仿真模型与测试说明：

为了研究直流梯度磁场为非线性时磁纳米温度成像方法的有效性，仿真测试过程中采用的试剂粒子有效磁矩M_s测定为8.5×10^-19(补充说明有效磁矩的测定数值由试剂类型参数决定)，成像范围40mm，直流梯度磁场的宽度Δx＝20mm，模拟机械移动的方式每次移动1mm。本仿真实例采用非线性的直流梯度磁场，交流磁场幅值H₀＝60Oe，频率1.6kHz，直流磁场强度为30Oe，并采用MATLAB中awgn函数给交流磁化强度信号添加噪声。当磁纳米粒子浓度分布均匀时，在温度成像过程中添加信噪比为90dB、80dB的噪声，温度测试范围为300K至310K，仿真结果如图14和图15所示。

2.仿真试验结果：

图14反映了在直流梯度磁场为非线性时，磁纳米温度成像方法的精度和适用性。在信噪比为90dB的条件下，温度测量误差小于0.36K。表明在直流梯度磁场为非线性时，该方法仍可以得到较好温度求解精度。

图15反映了在直流梯度磁场为非线性时，磁纳米温度成像方法的精度和适用性。在信噪比为80dB的条件下，温度测量误差小于0.54K。表明在直流梯度磁场为非线性时，该方法仍可以得到较好温度求解精度，且具备一定的抗噪声能力。

因此，这种磁纳米温度成像方法的精度、稳定性以及重复性是有保证的。为在复杂环境下完成精密快速的生物体温度场非侵入式成像提供了可靠的方法。

Claims (1)

1.一种磁纳米温度成像系统，其特征在于，包括：

恒定直流磁场产生装置，用于向磁纳米试剂所在一维空间施加恒定直流磁场H_dc＝b，其中，b为恒定直流磁场的幅值；

组合直流磁场产生装置，用于向磁纳米试剂所在一维空间施加组合直流磁场 ${H_{dc}}^{\&prime;}=\left\{\begin{array}{cc}b& x<x_1\\ f\left(x\right)& x_1\&le;x\&le;x_1+\mathrm{\&Delta;}x\\ -b& x>x_1+\mathrm{\&Delta;}x\end{array}\right.,$ 其中，x表示位置变量，x₁表示直流梯度磁场f(x)相对一维待测空间的起始位置，Δx为直流梯度磁场f(x)的宽度；

交流激励磁场产生装置，用于向磁纳米试剂所在一维空间施加交流激励磁场；

磁化强度采集装置，用于采集待测空间磁纳米试剂的交流磁化强度信号；

处理器，用于对同时施加恒定直流磁场和交流激励磁场时采集的交流磁化强度信号进行处理得到该交流磁化强度信号的各奇次谐波幅值A₁,A₃…A_2j-1，对在同时施加组合直流磁场和交流激励磁场时采集的交流磁化强度信号进行处理得到该交流磁化强度信号的各奇次谐波幅值B₁,B₃…B_2j-1，计算各奇次谐波幅值B₁,B₃…B_2j-1与A₁,A₃…A_2j-1之间的差值S₁,S₃…S_2j-1；根据各奇次谐波幅值差与在体温度的关系式求解[x₁,x₁+Δx]区间的y，进而计算得到[x₁,x₁+Δx]区间的在体温度

其中，c＝NM_s，N为磁纳米粒子浓度，M_s为磁纳米样品原子有效磁矩，k₀为波尔兹曼常数，磁纳米粒子的磁化强度由朗之万函数描述，c·f_2j-1(y,z(r_k))为朗之万函数的有限项泰勒级数展开得到的第2j-1次谐波幅值表达式，z(r_k)为[x₁,x₁+Δx]区间内的第k个离散点r_k对应的直流磁场大小，m为[x₁,x₁+Δx]区间被离散化的点数，谐波个数j≥2。