[发明专利]基于深度图像先验的MPI系统矩阵复原方法

说明书

本发明属于生物医学分子影像领域，具体涉及一种基于深度图像先验的MPI系统矩阵复原方法、系统、设备，旨在解决现有的系统矩阵获取方法操作繁琐，采集时间长，而且针对新的粒子，需要重新扫描，严重影响了MPI的使用效率的问题。本发明方法包括：构建欠采样的系统矩阵；将欠采样的系统矩阵按行进行RGB编码，得到RGB图像；将随机噪声作为网络输入，输入到3D‑Unet神经网络，将各RGB图像作为真值标签，每次迭代计算损失并回传，迭代结束后输出全采样的RGB图像；将全采样的RGB图像解码成复数形式，得到基于复数的系统矩阵。本发明简化了系统矩阵的获取操作，缩短了采集时间，并提升了MPI的使用效率。

技术领域

本发明属于生物医学分子影像领域，具体涉及一种基于深度图像先验的MPI系统矩阵复原方法、系统、设备。

背景技术

磁粒子成像技术(Magnetic particle imaging,MPI)是一种新兴的断层成像技术，它以高时空分辨率对磁性纳米粒子进行在体定量成像。此外，MPI只对磁性纳米粒子进行特异性成像，没有背景信号干扰，因此，采集到的信号强度与磁性纳米粒子浓度成正比，这使得MPI在生物医学应用中成为一种非常有前景的成像方式，例如心血管疾病的可视化。

为了重建磁性纳米粒子的空间浓度分布，需要借助于系统矩阵，构建测量电压信号与实际粒子浓度分布之间的线性关系。系统矩阵描述了粒子浓度与接收线圈中感应电压之间的映射关系。然而，到目前为止，确定MPI系统矩阵的最精确方法是使用一种特殊的校准程序。基本的思想是根据事先设定的三维网格数据，用一个充满磁性纳米颗粒的小样本填充网格上的一个像素或者说一个网格。这个小样本使用机器人放置在视野(FOV)的所有网格位置。在每个位置，进行MPI信号采集，采集到的信号作为系统矩阵的一列。校准方法的最大缺点是操作繁琐，采集时间长，例如中等大小34*28*20的网格图像，校准时间就大约需要6小时。一般情况下，网格的密度越高，图像的分辨率会越高，校准过程花费时间也就会越高。此外，由于磁性纳米粒子合成过程中的操作偏差，针对每一批新的粒子，系统矩阵通常需要重新扫描。而且任何采集参数的变化，如驱动场或选择场的强度变化，也需要重新测量系统矩阵。这一问题严重影响了MPI的使用效率，也会极大的限制MPI的应用前景。因此，为了进一步拓展MPI成像在生物医学研究中的应用范围，需要研究一种快速而又准确的系统矩阵复原方法以加快系统矩阵的校准过程。基于此，本发明提出了一种基于深度图像先验的MPI系统矩阵复原方法。

发明内容

为了解决现有技术中的上述问题，即为了解决现有的系统矩阵获取方法操作繁琐，采集时间长，而且针对新的粒子，需要重新扫描，严重影响了MPI的使用效率，也会极大的限制MPI的应用的问题，本发明提出了一种基于深度图像先验的MPI系统矩阵复原方法，利用欠采样的系统矩阵复原出完整的全采样的系统矩阵，进而基于复原后的系统矩阵，构建测量电压信号与实际粒子浓度分布之间的线性关系，重建磁性纳米粒子的空间浓度分布，进行磁粒子成像，该方法包括：

步骤S100，基于MPI成像时设定的欠采样网格的长宽高，构建欠采样的系统矩阵其中，M、N表示欠采样的系统矩阵的行、列数；

步骤S200，将所述欠采样的系统矩阵按行进行RGB编码，得到M个RGB图像；

步骤S300，将随机噪声作为网络输入，输入到预构建的3D-Unet神经网络，将各RGB图像作为真值标签，每次迭代计算损失并回传，迭代结束后输出全采样的RGB图像；

步骤S400，将全采样的RGB图像解码成复数形式，得到基于复数的系统矩阵，即复原后的系统矩阵。

在一些优选的实施方式中，所述欠采样的系统矩阵，其构建方法为：

根据MPI成像时设定的欠采样网格的长、宽、高，机器人将一个充满磁性纳米颗粒的小样本自动放置在视野的每个网格位置，然后进行MPI信号采集，每个位置采集到的电压信号作为系统矩阵的一列，采集完所有网格位置之后，即获得了相应的欠采样的系统矩阵。