US 12,136,145 B2
Signal amplitude feature-based method for fast reconstructing a magnetic particle imaging and device
Jie Tian, Beijing (CN); Peng Zhang, Beijing (CN); Hui Hui, Beijing (CN); Yimeng Li, Beijing (CN); Lin Yin, Beijing (CN); and Xin Feng, Beijing (CN)
Assigned to Institute Of Automation, Chinese Academy of Sciences, Beijing (CN)
Filed by INSTITUTE OF AUTOMATION, CHINESE ACADEMY OF SCIENCES, Beijing (CN)
Filed on Jul. 11, 2022, as Appl. No. 17/811,738.
Claims priority of application No. 202110816843.9 (CN), filed on Jul. 20, 2021.
Prior Publication US 2023/0027988 A1, Jan. 26, 2023
Int. Cl. G06T 11/00 (2006.01); A61B 5/00 (2006.01); A61B 5/0515 (2021.01)
CPC G06T 11/006 (2013.01) [A61B 5/0515 (2013.01); A61B 5/7257 (2013.01)] 8 Claims
OG exemplary drawing
 
1. A signal amplitude feature-based method for fast reconstructing a magnetic particle imaging, comprising:
transforming a time-domain voltage signal collected by an MPI system device to a frequency domain, so as to obtain a frequency domain complex voltage signal u;
respectively calculating a square root of a square sum of a real part and an imaginary part at each frequency point of the frequency domain complex voltage signal u, so as to obtain an amplitude at the each frequency point;
arranging the amplitude at the each frequency point in a descending order so as to obtain an amplitude matrix D, and acquiring a screening threshold i0 by an amplitude ratio method;
screening an element of the amplitude matrix D based on the screening threshold i0, and constructing a frequency domain signal array u0 based on the screened elements D1˜i0;
selecting row vectors a1˜i0 of a system matrix A corresponding to each frequency point of the frequency domain signal array u0, so as to construct an update system matrix A0; and
solving, based on the frequency domain signal array u0 and the update system matrix A0, an inverse problem in a form of a least square based on an L2 constraint by an convex optimization method to obtain a three-dimensional magnetic particle concentration distribution result, so as to achieve a fast reconstruction of the MPI system.