EPR oximetry in three spatial dimensions using sparse spin distribution

Implementation of GPU-accelerated back projection for EPR imaging

Electron paramagnetic resonance (EPR) Imaging (EPRI) is a robust method for measuring in vivo oxygen concentration (pO2). For 3D pulse EPRI, a commonly used reconstruction algorithm is the filtered backprojection (FBP) algorithm, in which the backprojection process is computationally intensive and may be time consuming when implemented on a CPU. A multistage implementation of the backprojection can be used for acceleration, however it is not flexible (requires equal linear angle projection distribution) and may still be time consuming. In this work, single-stage backprojection is implemented on a GPU (Graphics Processing Units) having 1152 cores to accelerate the process. The GPU implementation results in acceleration by over a factor of 200 overall and by over a factor of 3500 if only the computing time is considered. Some important experiences regarding the implementation of GPU-accelerated backprojection for EPRI are summarized. The resulting accelerated image reconstruction is useful for real-time image reconstruction monitoring and other time sensitive applications.

Comparison of parabolic filtration methods for 3D filtered back projection in pulsed EPR imaging

Pulse electron paramagnetic resonance imaging (Pulse EPRI) is a robust method for noninvasively measuring local oxygen concentrations in vivo. For 3D tomographic EPRI, the most commonly used reconstruction algorithm is filtered back projection (FBP), in which the parabolic filtration process strongly influences image quality. In this work, we designed and compared 7 parabolic filtration methods to reconstruct both simulated and real phantoms. To evaluate these methods, we designed 3 error criteria and 1 spatial resolution criterion. It was determined that the 2 point derivative filtration method and the two-ramp-filter method have unavoidable negative effects resulting in diminished spatial resolution and increased artifacts respectively. For the noiseless phantom the rectangular-window parabolic filtration method and sinc-window parabolic filtration method were found to be optimal, providing high spatial resolution and small errors. In the presence of noise, the 3 point derivative method and Hamming-window parabolic filtration method resulted in the best compromise between low image noise and high spatial resolution. The 3 point derivative method is faster than Hamming-window parabolic filtration method, so we conclude that the 3 point derivative method is optimal for 3D FBP.

Clinical EPR: unique opportunities and some challenges

Electron paramagnetic resonance (EPR) spectroscopy has been well established as a viable technique for measurement of free radicals and oxygen in biological systems, from in vitro cellular systems to in vivo small animal models of disease. However, the use of EPR in human subjects in the clinical setting, although attractive for a variety of important applications such as oxygen measurement, is challenged with several factors including the need for instrumentation customized for human subjects, probe, and regulatory constraints. This article describes the rationale and development of the first clinical EPR systems for two important clinical applications, namely, measurement of tissue oxygen (oximetry) and radiation dose (dosimetry) in humans. The clinical spectrometers operate at 1.2 GHz frequency and use surface-loop resonators capable of providing topical measurements up to 1 cm depth in tissues. Tissue pO2 measurements can be carried out noninvasively and repeatedly after placement of an oxygen-sensitive paramagnetic material (currently India ink) at the site of interest. Our EPR dosimetry system is capable of measuring radiation-induced free radicals in the tooth of irradiated human subjects to determine the exposure dose. These developments offer potential opportunities for clinical dosimetry and oximetry, which include guiding therapy for individual patients with tumors or vascular disease by monitoring of tissue oxygenation. Further work is in progress to translate this unique technology to routine clinical practice.

Compressed sensing of spatial electron paramagnetic resonance imaging

PURPOSE

To improve image quality and reduce data requirements for spatial electron paramagnetic resonance imaging (EPRI) by developing a novel reconstruction approach using compressed sensing (CS).

METHODS

EPRI is posed as an optimization problem, which is solved using regularized least-squares with sparsity promoting penalty terms, consisting of the l1 norms of the image itself and the total variation of the image. Pseudo-random sampling was employed to facilitate recovery of the sparse signal. The reconstruction was compared with the traditional filtered back-projection reconstruction for simulations, phantoms, isolated rat hearts, and mouse gastrointestinal (GI) tracts labeled with paramagnetic probes.

RESULTS

A combination of pseudo-random sampling and CS was able to generate high-fidelity EPR images at high acceleration rates. For three-dimensional (3D) phantom imaging, CS-based EPRI showed little visual degradation at nine-fold acceleration. In rat heart datasets, CS-based EPRI produced high quality images with eight-fold acceleration. A high resolution mouse GI tract reconstruction demonstrated a visual improvement in spatial resolution and a doubling in signal-to-noise ratio (SNR).

CONCLUSION

A novel 3D EPRI reconstruction using compressed sensing was developed and offers superior SNR and reduced artifacts from highly undersampled data.

Electron paramagnetic resonance oxygen imaging in vivo

This review covers the last 15 years of the development of EPR in vivo oxygen imaging. During this time, a number of major technological and methodological advances have taken place. Narrow line width, long relaxation time, and non-toxic triaryl methyl radicals were introduced in the late 1990s. These not only improved continuous wave (CW) imaging, but also enabled the application of pulse EPR imaging to animals. Recent developments in pulse technology have brought an order of magnitude increase in image acquisition speed, enhancement of sensitivity, and considerable improvement in the precision and accuracy of oxygen measurements. Consequently, pulse methods take up a significant part of this review.

Dual-scan acquisition for accelerated continuous-wave EPR oximetry

Statistical analysis reveals that, given a fixed acquisition time, linewidth (and thus pO(2)) can be more precisely determined from multiple scans with different modulation amplitudes and sweep widths than from a single-scan. For a Lorentzian lineshape and an unknown but spatially uniform modulation amplitude, the analysis suggests the use of two scans, each occupying half of the total acquisition time. We term this mode of scanning as dual-scan acquisition. For unknown linewidths in a range [Γ(min), Γ(max)], practical guidelines are provided for selecting the modulation amplitude and sweep width for each dual-scan component. Following these guidelines can allow for a 3-4 times reduction in spectroscopic acquisition time versus an optimized single-scan, without requiring hardware modifications. Findings are experimentally verified using L-band spectroscopy with an oxygen-sensitive particulate probe.

Multisite EPR oximetry from multiple quadrature harmonics

Multisite continuous wave (CW) electron paramagnetic resonance (EPR) oximetry using multiple quadrature field modulation harmonics is presented. First, a recently developed digital receiver is used to extract multiple harmonics of field modulated projection data. Second, a forward model is presented that relates the projection data to unknown parameters, including linewidth at each site. Third, a maximum likelihood estimator of unknown parameters is reported using an iterative algorithm capable of jointly processing multiple quadrature harmonics. The data modeling and processing are applicable for parametric lineshapes under nonsaturating conditions. Joint processing of multiple harmonics leads to 2-3-fold acceleration of EPR data acquisition. For demonstration in two spatial dimensions, both simulations and phantom studies on an L-band system are reported.

Digital detection and processing of multiple quadrature harmonics for EPR spectroscopy

A quadrature digital receiver and associated signal estimation procedure are reported for L-band electron paramagnetic resonance (EPR) spectroscopy. The approach provides simultaneous acquisition and joint processing of multiple harmonics in both in-phase and out-of-phase channels. The digital receiver, based on a high-speed dual-channel analog-to-digital converter, allows direct digital down-conversion with heterodyne processing using digital capture of the microwave reference signal. Thus, the receiver avoids noise and nonlinearity associated with analog mixers. Also, the architecture allows for low-Q anti-alias filtering and does not require the sampling frequency to be time-locked to the microwave reference. A noise model applicable for arbitrary contributions of oscillator phase noise is presented, and a corresponding maximum-likelihood estimator of unknown parameters is also reported. The signal processing is applicable for Lorentzian lineshape under nonsaturating conditions. The estimation is carried out using a convergent iterative algorithm capable of jointly processing the in-phase and out-of-phase data in the presence of phase noise and unknown microwave phase. Cramér-Rao bound analysis and simulation results demonstrate a significant reduction in linewidth estimation error using quadrature detection, for both low and high values of phase noise. EPR spectroscopic data are also reported for illustration.

In vivo multisite oximetry using EPR-NMR coimaging

Coimaging employing electron paramagnetic resonance (EPR) imaging and MRI is used for rapid in vivo oximetry conducted simultaneously across multiple organs of a mouse. A recently developed hybrid EPR-NMR coimaging instrument is used for both EPR and NMR measurements. Oxygen sensitive particulate EPR probe is implanted in small localized pockets, called sites, across multiple regions of a live mouse. Three dimensional MRI is used to generate anatomic visualization, providing precise locations of implant sites. The pO₂ values, one for every site, are then estimated from EPR measurements. To account for radio frequency (RF) phase inhomogeneities inside a large resonator carrying a lossy sample, a generalization of an existing EPR data model is proposed. Utilization of known spectral lineshape, sparse distribution, and known site locations reduce the EPR data collection by more than an order of magnitude over a conventional spectral-spatial imaging, enhancing the feasibility of in vivo EPR oximetry for clinically relevant models.

Implantable resonators--a technique for repeated measurement of oxygen at multiple deep sites with in vivo EPR

EPR oximetry using implantable resonators allows measurements at much deeper sites than are possible with surface resonators (> 80 vs. 10 mm) and achieves greater sensitivity at any depth. We report here the development of an improved technique that enables us to obtain the information from multiple sites and at a variety of depths. The measurements from the various sites are resolved using a simple magnetic field gradient. In the rat brain multi-probe implanted resonators measured pO(2) at several sites simultaneously for over 6 months under normoxic, hypoxic, and hyperoxic conditions. This technique also facilitates measurements in moving parts of the animal such as the heart, because the orientation of the paramagnetic material relative to the sensing loop is not altered by the motion. The measured response is fast, enabling measurements in real time of physiological and pathological changes such as experimental cardiac ischemia in the mouse heart. The technique also is quite useful for following changes in tumor pO(2), including applications with simultaneous measurements in tumors and adjacent normal tissues.