基于信号幅值特征的磁粒子成像快速重建方法及设备与流程

1.本发明属于无磁场点的三维磁粒子成像系统中的立体重建技术领域，具体涉及了一种基于信号幅值特征的磁粒子成像快速重建方法及设备。


背景技术：

2.在临床诊断和检测中，如何准确、客观的定位肿瘤及其他病灶一直是国际上的研究热点和挑战性问题。现有的医学影像技术如ct，mri，spect等方法均存在危害大，定位差，精度低等问题。而在近些年，一种全新的基于示踪剂的成像方式——磁粒子成像(mpi)技术被提出。mpi可以通过检测对人体无害的超顺磁氧化铁纳米颗粒(spions)的空间浓度分布，对肿瘤或目标物进行精准定位，是一种不受成像深度限制的三维高时空分辨率和高灵敏度的成像。此外，mpi不显示解剖结构并且无背景信号干扰，因此信号的强度与示踪剂的浓度存在直接关系，是一种颇具医学应用潜力的新方法。
3.现今的mpi系统的重建方法基本可以分为两类：系统矩阵方法和x
‑
space方法。大量的研究表面，相比于系统矩阵重建方法，x
‑
space方法虽然具有更高的重建速度，但是其重建图像的分辨率很难提升，因此系统矩阵方法一直是图像重建的主要研究方向。目前已有研究人员通过超算平台，实现了基于系统矩阵方法的动态三维磁粒子成像技术。但由于重建计算量较大，其动态影像的帧率较低，仅能实现数帧级别的动态重建。同时由于计算过程中存在庞大的数据传输和矩阵操作，该技术严重依赖硬件设备，对包括存储，传输，计算等多方面设备提出了严苛的要求。


技术实现要素：

4.为了解决现有技术中的上述问题，即数据量、计算量庞大导致的三维磁粒子成像系统中的立体重建效率低、对设备资源要求高的问题，本发明提供了一种基于信号幅值特征的磁粒子成像快速重建方法，该重建方法包括：
5.步骤s10，将mpi系统设备采集的时域电压信号变换至频域，获得频域复数电压信号u；
6.步骤s20，分别计算所述频域复数电压信号u的每一频点的实部和虚部的平方和的平方根，获得每一频点的幅值；
7.步骤s30，将所述每一频点的幅值进行降序排列，获得幅值矩阵d，并通过幅值占比的方法，获取筛选阈值i0；
8.步骤s40，基于所述筛选阈值i0，进行所述幅值矩阵d的元素筛选，并基于筛选出的元素构建频域信号数组u0；
9.步骤s50，选取所述频域信号数组u0中每一频点对应的系统矩阵a的行向量构建更新系统矩阵a0；
10.步骤s60，基于所述频域信号数组u0以及所述更新系统矩阵a0，通过凸优化方法求解基于l2约束的最小二乘形式的逆问题，获得三维磁粒子浓度分布结果，实现mpi系统的快
速重建。
11.在一些优选的实施例中，步骤s30中通过幅值占比的方法，计算筛选阈值i0，其方法为：
12.步骤s31，计算所述幅值矩阵d的所有幅值的和n，以第1～i个降序排序后的幅值的和除以n作为幅值占比元素m
i
的值；
13.步骤s32，选取幅值占比元素m
i
中与设定第一阈值的差异值小于设定第二阈值的m个元素对应的i值，记作筛选阈值i0。
14.在一些优选的实施例中，所述设定第一阈值为0.8。
15.在一些优选的实施例中，所述幅值占比元素m
i
，其表示为：
16.其中，d
k
代表幅值矩阵d中第k个幅值，n代表幅值矩阵d的所有幅值的和，i代表1～i个降序排序后的幅值的第i个。
17.在一些优选的实施例中，所述凸优化方法为kaczmarz方法、admm方法、tv
‑
l2方法中的一种。
18.本发明的另一方面，提出了一种基于信号幅值特征的mpi快速重建系统，该重建系统包括以下模块：
19.信号转换模块，配置为将mpi系统设备采集的时域电压信号变换至频域，获得频域复数电压信号u；
20.幅值计算模块，配置为分别计算所述频域复数电压信号u的每一频点的实部和虚部的平方和的平方根，获得每一频点的幅值；
21.筛选阈值获取模块，配置为将所述每一频点的幅值进行降序排列，获得幅值矩阵d，并通过幅值占比的方法，获取筛选阈值i0；
22.幅值筛选模块，配置为基于所述筛选阈值i0，进行所述幅值矩阵d的元素筛选，并基于筛选出的元素构建频域信号数组u0；
23.系统矩阵更新模块，配置为选取所述频域信号数组u0中每一频点对应的系统矩阵a的行向量构建更新系统矩阵a0；
24.重建模块，配置为基于所述频域信号数组u0以及所述更新系统矩阵a0，通过凸优化方法求解基于l2约束的最小二乘形式的逆问题，获得三维磁粒子浓度分布结果，实现mpi系统的快速重建。
25.本发明的第三方面，提出了一种电子设备，包括：
26.至少一个处理器；以及
27.与至少一个所述处理器通信连接的存储器；其中，
28.所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于信号幅值特征的磁粒子成像快速重建方法。
29.本发明的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于信号幅值特
征的磁粒子成像快速重建方法。
30.本发明的有益效果：
31.(1)本发明基于信号幅值特征的磁粒子成像快速重建方法，可以在几乎不降低重建精度的前提下，极大降低重建所需计算量，进而提升了重建速度并降低了重建硬件要求，可以用于提升三维动态mpi成像的帧率。
32.(2)本发明基于信号幅值特征的磁粒子成像快速重建方法，由于大幅度降低计算量，使得重建可以脱离超算平台，仅使用一般工作站即可完成重建，提升了mpi设备的推广潜力。
附图说明
33.通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本技术的其它特征、目的和优点将会变得更明显：
34.图1是本发明基于信号幅值特征的磁粒子成像快速重建方法的流程示意图；
35.图2是本发明基于信号幅值特征的磁粒子成像快速重建方法一种实施例的幅值占比图以及分割阈值；
36.图3是本发明基于信号幅值特征的磁粒子成像快速重建方法一种实施例的筛选后的信号分布茎叶图；
37.图4是本发明基于信号幅值特征的磁粒子成像快速重建方法一种实施例的重建结果三维切面图。
具体实施方式
38.下面结合附图和实施例对本技术作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与有关发明相关的部分。
39.需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本技术。
40.本发明提供一种基于信号幅值特征的磁粒子成像快速重建方法，本方法利用幅值占比和阈值计算的方法，对原始信号进行筛选，使得该信号能够以较小的数据量，保持与原有信号一致的信息含量和分布特征。同时通过以频域信号为依据的系统矩阵筛选，可以将重建所需的系统矩阵大小一同降低，使得在重建过程中的计算量大幅度降低(为原数据量的1％以下)，进而大大提升了重建速度，并降低了所需的硬件设备要求。
41.本发明的一种基于信号幅值特征的磁粒子成像快速重建方法，该重建方法包括：
42.步骤s10，将mpi系统设备采集的时域电压信号变换至频域，获得频域复数电压信号u；
43.步骤s20，分别计算所述频域复数电压信号u的每一频点的实部和虚部的平方和的平方根，获得每一频点的幅值；
44.步骤s30，将所述每一频点的幅值进行降序排列，获得幅值矩阵d，并通过幅值占比的方法，获取筛选阈值i0；
45.步骤s40，基于所述筛选阈值i0，进行所述幅值矩阵d的元素筛选，并基于筛选出的
元素构建频域信号数组u0；
46.步骤s50，选取所述频域信号数组u0中每一频点对应的系统矩阵a的行向量构建更新系统矩阵a0；
47.步骤s60，基于所述频域信号数组u0以及所述更新系统矩阵a0，通过凸优化方法求解基于l2约束的最小二乘形式的逆问题，获得三维磁粒子浓度分布结果，实现mpi系统的快速重建。
48.为了更清晰地对本发明基于信号幅值特征的磁粒子成像快速重建方法进行说明，下面结合图1对本发明实施例中各步骤展开详述。
49.本发明第一实施例的基于信号幅值特征的磁粒子成像快速重建方法，包括步骤s10
‑
步骤s60，各步骤详细描述如下：
50.步骤s10，将mpi系统设备采集的时域电压信号变换至频域，获得频域复数电压信号u。
51.本发明一个实施例中，通过傅里叶变换将mpi系统设备采集的时域电压信号变换至频域，获得复数形式的频域复数电压信号u。
52.步骤s20，分别计算所述频域复数电压信号u的每一频点的实部和虚部的平方和的平方根，获得每一频点的幅值。
53.步骤s30，将所述每一频点的幅值进行降序排列，获得幅值矩阵d，并通过幅值占比的方法，获取筛选阈值i0：
54.步骤s31，计算所述幅值矩阵d的所有幅值的和n，以第1～i个降序排序后的幅值的和除以n作为幅值占比元素m
i
的值。
55.幅值占比元素m
i
，其表示如式(1)所示：
56.其中，d
k
代表幅值矩阵d中第k个幅值，n代表幅值矩阵d的所有幅值的和，i代表1～i个降序排序后的幅值的第i个。
57.步骤s32，选取幅值占比元素m
i
中与设定第一阈值的差异值小于设定第二阈值的m个元素对应的i值，记作筛选阈值i0。
58.本发明一个实施例中，设定第一阈值为0.8，获取m种最接近于0.8的元素对应的i值，记为筛选阈值i0。
59.如图2所示，为本发明基于信号幅值特征的磁粒子成像快速重建方法一种实施例的幅值占比图以及分割阈值，横坐标代表根据幅值降序排列后的幅值矩阵的位置k，纵坐标代表幅值占比m，曲线代表随着包含的幅值数量上升，幅值占比的变换，虚线代表阈值。
60.步骤s40，基于所述筛选阈值i0，进行所述幅值矩阵d的元素筛选，并基于筛选出的元素构建频域信号数组u0。
61.如图3所示，为本发明基于信号幅值特征的磁粒子成像快速重建方法一种实施例的筛选后的信号分布茎叶图，横坐标代表频点值，纵坐标代表归一化后的幅值大小，每个空心圆圈代表在该频点的幅值大小，图3为筛选出前625个较大幅值的信号分布茎叶图，可以
看到其集中在主要频点周围。
62.步骤s50，选取所述频域信号数组u0中每一频点对应的系统矩阵a的行向量构建更新系统矩阵a0。
63.可以采用一个标准浓度的较小的测量样本，移动遍历整个fov，并在每次移动后测量一次信号作为矩阵a的一行，遍历后获得整的系统矩阵a，也可以采用其他方法获取系统矩阵a，本发明在此不一一详述。
64.步骤s60，基于所述频域信号数组u0以及所述更新系统矩阵a0，通过凸优化方法求解基于l2约束的最小二乘形式的逆问题，获得三维磁粒子浓度分布结果，实现mpi系统的快速重建。
65.凸优化方法为kaczmarz方法、admm方法、tv
‑
l2方法中的一种。本发明一个实施例中，采用kaczmarz方法。
66.如图4所示，为本发明基于信号幅值特征的磁粒子成像快速重建方法一种实施例的重建结果三维切面图，横坐标为位置的x轴坐标，纵坐标为位置的y轴坐标，slice at xz&y＝10代表其对应的图为xy平面在y＝10处的切面，图4左上图、右上图、下图代表不同角度的切面
67.上述实施例中虽然将各个步骤按照上述先后次序的方式进行了描述，但是本领域技术人员可以理解，为了实现本实施例的效果，不同的步骤之间不必按照这样的次序执行，其可以同时(并行)执行或以颠倒的次序执行，这些简单的变化都在本发明的保护范围之内。
68.本发明第二实施例的基于信号幅值特征的mpi快速重建系统，该重建系统包括以下模块：
69.信号转换模块，配置为将mpi系统设备采集的时域电压信号变换至频域，获得频域复数电压信号u；
70.幅值计算模块，配置为分别计算所述频域复数电压信号u的每一频点的实部和虚部的平方和的平方根，获得每一频点的幅值；
71.筛选阈值获取模块，配置为将所述每一频点的幅值进行降序排列，获得幅值矩阵d，并通过幅值占比的方法，获取筛选阈值i0；
72.幅值筛选模块，配置为基于所述阈值i0，进行所述幅值矩阵d的元素筛选，并基于筛选出的元素构建频域信号数组u0；
73.系统矩阵更新模块，配置为选取所述频域信号数组u0中每一频点对应的系统矩阵a的行向量构建更新系统矩阵a0；
74.重建模块，配置为基于所述频域信号数组u0以及所述更新系统矩阵a0，通过凸优化方法求解基于l2约束的最小二乘形式的逆问题，获得三维磁粒子浓度分布结果，实现mpi系统的快速重建。
75.所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
76.需要说明的是，上述实施例提供的基于信号幅值特征的mpi快速重建系统，仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不
同的功能模块来完成，即将本发明实施例中的模块或者步骤再分解或者组合，例如，上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块，以完成以上描述的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称，仅仅是为了区分各个模块或者步骤，不视为对本发明的不当限定。
77.本发明第三实施例的一种电子设备，包括：
78.至少一个处理器；以及
79.与至少一个所述处理器通信连接的存储器；其中，
80.所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于信号幅值特征的磁粒子成像快速重建方法。
81.本发明第四实施例的一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于信号幅值特征的磁粒子成像快速重建方法。
82.所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的存储装置、处理装置的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
83.本领域技术人员应该能够意识到，结合本文中所公开的实施例描述的各示例的模块、方法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，软件模块、方法步骤对应的程序可以置于随机存储器(ram)、内存、只读存储器(rom)、电可编程rom、电可擦除可编程rom、寄存器、硬盘、可移动磁盘、cd
‑
rom、或技术领域内所公知的任意其它形式的存储介质中。为了清楚地说明电子硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以电子硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。本领域技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
84.术语“第一”、“第二”等是用于区别类似的对象，而不是用于描述或表示特定的顺序或先后次序。
85.术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素，而且还包括没有明确列出的其它要素，或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
86.至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征做出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。