Information-based Unified Design Framework for NMR Magnetic Tracking System
Authors: Seçkin Eren Yetim, Mehmet Efe Tiryaki
Organizations: METU-ROMER, Center for Robotics and AI, Middle East Technical University, Ankara, Türkiye. · Department of Mechanical Engineering, Middle East Technical University, Ankara, Türkiye.
Nuclear magnetic resonance (NMR) magnetic tracking offers a promising alternative to magnetic resonance imaging (MRI)-based localization for medical robotic guidance. However, scalar field measurements and field-aligned marker moments constrain position and orientation tracking, requiring careful configuration design. We present a unified information-based framework that uses local covariance analysis to evaluate receiver geometry, marker separation, and external source configuration in on/off-board NMR tracking. Numerical evaluation shows that planar-array precision depends jointly on receiver count and geometry, with triangular lattice arrays achieving the lowest best-case position uncertainty among the evaluated layouts for 9--20 receivers. Two spatially separated, field-aligned markers enable device-axis tracking, with predicted angular uncertainty below 2∘ at the evaluated workspace center for suitable marker separations. We further establish point-dipole source--receiver duality, unifying the analysis of on/off-board NMR tracking. These
Figures & tables
Fig. 1: Magnetic-tracking architectures considered in this study.MRI setting and representative distance-dependent measurements δBi(ri) from seven receivers. (A) Off-board sensing of a device-mounted, field-aligned marker m by an external receiver array. (B) On-board sensing of fields δBk from independently activated external sources with moments mk .
Fig. 2: (A) Single marker–receiver geometry with depth z , lateral separation ρ , and local field δB .(B) Two identical field-aligned markers separated by L define the unoriented device axis u .
Fig. 3: Local 3-D position-uncertainty score for the three-receiver planar array. (A) Measurement geometry. (B) Top view of the equilateral reference array. (C) Ap map in the x=0 plane. (D–F) Ap maps at z=5 , 20, and 60 mm. All maps use a logarithmic scale with Ap reported in millimeters.
Fig. 4: Planar receiver-array comparison with 20-mm minimum pairwise spacing. (A) Representative layouts from the four evaluated families. (B) Minimum sampled Ap versus receiver count N at z=60 mm over the 100×100 mm 2 region; the inset enlarges N=18 –20.
Fig. 5: Worst-over-direction two-marker uncertainty. (A) Au workspace map at z=60 mm for L=20 mm; contours at 1∘ and 2∘ . (B) Corresponding Ap map; contours at 0.1 and 1 mm. (C) Center-point Au versus L . (D) Center-point Ap versus L . Open circles mark Lref=20 mm; (A–C) use logarithmic scales.
Fig. 6: Sequential scalar on-board sensing with three independently activated sources. (A) Distributed–parallel configuration. (B) Co-located–orthogonal configuration. (C) Distributed–parallel Ap map in the y=0 plane. (D) Co-located–orthogonal Ap map. Contours mark 0.01 , 0.1 , and 1 mm; both maps share a 0–2 mm scale, with larger values saturated.
Wireless localization of permanent magnets enables occlusion-free guidance for medical interventions, yet its practical accuracy is fundamentally limited by two coupled challenges: the poor observability of conventional planar sensor arrays and the simulation-to-reality (Sim-to-Real) gap of learning-based estimators. To address these issues, this article presents a unified framework that combines information-theoretic sensor geometry optimization with physics-aware deep learning. First, a rigorous Fisher Information Matrix (FIM)-based evaluation framework is established to quantify geometry-induced observability limitations. The results show that a staggered split-array topology provides a substantially stronger observability foundation for localization while remaining compatible with practical external deployment. Second, building on this optimized sensing configuration, we propose Phy-GAANet, a calibration-free estimator trained entirely on hardware-aware synthetic data. By incorporating Physics-Informed Features (PIF) for saturation modeling and Geometry-Aware Attention (GAA) for preserving cross-layer vector structure, the network effectively bridges the Sim-to-Real gap. Extensive real-world experiments demonstrate state-of-the-art performance, achieving a position error of 1.84 mm and an orientation error of 3.18 degrees at a refresh rate exceeding 270 Hz. The proposed method consistently outperforms classical Levenberg--Marquardt solvers and generic convolutional baselines, particularly in suppressing catastrophic outliers and maintaining robustness in challenging near-field boundary regions. Beyond the proposed network, the FIM-guided analysis also provides a framework for sensor geometry design in magnetic localization systems under practical deployment constraints.
Wenxuan Xie, Yuelin Zhang, Qingpeng Ding +4
Department of Mechanical and Automation Engineering and T Stone Robotics Institute, The Chinese University of Hong Kong, Hong Kong.
Embedded magnetic tracking holds highly attractive prospects for remote navigation of endoluminal medical devices. However, existing six-degree-of-freedom pose recovery approaches often require pre-calibrated workspace field maps or iterative nonlinear optimization. This letter presents a Gradiometer-Based Electromagnetic Localization System (GELS), a closed-form tracking framework that uses a compact magnetometer array as an embedded quasi-gradiometer to estimate local magnetic fields and gradient tensors. These quantities are mapped by the Euler homogeneous relation to displacements between source and array, from which multi-source Procrustes registration recovers the array orientation and position using at least three non-collinear sources. The algorithm requires known source positions and array geometry, but no pre-calibrated workspace field maps, initial pose guesses, or calibrated excitation-source moments. The recovered pose also enables a proof-of-concept sub-level dipole localization task by serving as a mobile magnetic reference frame. Benchtop experiments across sensor-array configurations and excitation modes demonstrate sequence-averaged position errors of \SI{10.80}{\milli\meter}--\SI{15.57}{\milli\meter}, a fastest update rate of \SI{14.49}{\hertz}, and a median solver runtime of \SI{172.00}{\micro\second}. A perturbation-based error propagation analysis further identifies inter-sensor inconsistency and dipole-model mismatch as the dominant accuracy limits, thereby informing future sensor array and magnetic source design for further reducing pose-estimation error.
Zhiwei Wu, Jiahao Luo, Yubo Pu +3
School of Automation, Beijing Institute of Technology, Beijing, 100081, China
Robots that localize a radio transmitter need more than a point estimate: in cluttered rooms, one measurement is often consistent with several transmitter locations and, because upper-mid-band antennas are directional, several headings. We present MAGNETAR, which infers a joint posterior over planar transmitter position and heading from a single asynchronous radio-frequency (RF) multipath snapshot, represented by angle-of-arrival and signal-to-noise-ratio estimates, given the room layout and receiver pose. Among our five neural scorers, MAGNETAR adopts a shared 2D U-Net conditioned on each candidate heading, jointly normalizing scores over a discretized position-heading grid. Training uses real-to-sim-calibrated 10 GHz simulations and a small measured subset. Grid-based joint posteriors outperform parametric ones on held-out simulations, the heading-conditioned scorer transfers best to robotic measurements, and fusing joint posteriors improves on fusing position-only marginals.
Haozhe Lei, Ruibin Chen, Yuhan Jiang +3
NYU WIRELESS, Tandon School of Engineering, New York University, Brooklyn, NY 11201, USA