Weighted iterative reconstruction for magnetic particle imaging

被引：153

作者：

Knopp, T. ^{[1
]}

Rahmer, J. ^{[2
]}

Sattel, T. F. ^{[1
]}

Biederer, S. ^{[1
]}

Weizenecker, J. ^{[2
]}

Gleich, B. ^{[2
]}

Borgert, J. ^{[2
]}

Buzug, T. M. ^{[1
]}

机构：

[1] Med Univ Lubeck, Inst Med Engn, D-23538 Lubeck, Germany

[2] Philips Technol GmbH Forschungslaboratorien, Hamburg, Germany

来源：

PHYSICS IN MEDICINE AND BIOLOGY | 2010年 / 55卷 / 06期

关键词：

ART ALGORITHM; CONVERGENCE; PROJECTIONS; PICTURES;

D O I：

10.1088/0031-9155/55/6/003

中图分类号：

R318 [生物医学工程];

学科分类号：

0831 ;

摘要：

Magnetic particle imaging (MPI) is a new imaging technique capable of imaging the distribution of superparamagnetic particles at high spatial and temporal resolution. For the reconstruction of the particle distribution, a system of linear equations has to be solved. The mathematical solution to this linear system can be obtained using a least-squares approach. In this paper, it is shown that the quality of the least-squares solution can be improved by incorporating a weighting matrix using the reciprocal of the matrix-row energy as weights. A further benefit of this weighting is that iterative algorithms, such as the conjugate gradient method, converge rapidly yielding the same image quality as obtained by singular value decomposition in only a few iterations. Thus, the weighting strategy in combination with the conjugate gradient method improves the image quality and substantially shortens the reconstruction time. The performance of weighting strategy and reconstruction algorithms is assessed with experimental data of a 2D MPI scanner.

引用

页码：1577 / 1589

页数：13

共 22 条

[1]

ANDERSON E., 1999, LAPACK USERSGUIDE, V3rd

[2]

[Anonymous], 1993, MATRIX COMPUTATIONS

[3] SMOOTHING NOISY DATA WITH SPLINE FUNCTIONS [J].