Source author record

Peder E. Z. Larson

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

9works
6topics
4close collaborators

Published work

9 publications

preprint2022arXiv

Physics-driven Deep Learning for PET/MRI

In this paper, we review physics- and data-driven reconstruction techniques for simultaneous positron emission tomography (PET) / magnetic resonance imaging (MRI) systems, which have significant advantages for clinical imaging of cancer, neurological disorders, and heart disease. These reconstruction approaches utilize priors, either structural or statistical, together with a physics-based description of the PET system response. However, due to the nested representation of the forward problem, direct PET/MRI reconstruction is a nonlinear problem. We elucidate how a multi-faceted approach accommodates hybrid data- and physics-driven machine learning for reconstruction of 3D PET/MRI, summarizing important deep learning developments made in the last 5 years to address attenuation correction, scattering, low photon counts, and data consistency. We also describe how applications of these multi-modality approaches extend beyond PET/MRI to improving accuracy in radiation therapy planning. We conclude by discussing opportunities for extending the current state-of-the-art following the latest trends in physics- and deep learning-based computational imaging and next-generation detector hardware.

preprint2021arXiv

Anisotropic field-of-views in radial imaging

Radial imaging techniques, such as projection-reconstruction (PR), are used in magnetic resonance imaging (MRI) for dynamic imaging, angiography, and short-imaging. They are robust to flow and motion, have diffuse aliasing patterns, and support short readouts and echo times. One drawback is that standard implementations do not support anisotropic field-of-view (FOV) shapes, which are used to match the imaging parameters to the object or region-of-interest. A set of fast, simple algorithms for 2-D and 3-D PR, and 3-D cones acquisitions are introduced that match the sampling density in frequency space to the desired FOV shape. Tailoring the acquisitions allows for reduction of aliasing artifacts in undersampled applications or scan time reductions without introducing aliasing in fully-sampled applications. It also makes possible new radial imaging applications that were previously unsuitable, such as imaging elongated regions or thin slabs. 2-D PR longitudinal leg images and thin-slab, single breath-hold 3-D PR abdomen images, both with isotropic resolution, demonstrate these new possibilities. No scan time to volume efficiency is lost by using anisotropic FOVs. The acquisition trajectories can be computed on a scan by scan basis.

preprint2021arXiv

Di-chromatic Interpolation of Magnetic Resonance Metabolic Imagery

Magnetic resonance imaging with hyperpolarized contrast agents can provide unprecedented \textit{in-vivo} measurements of metabolism, but yields images that are lower resolution than that achieved with proton anatomical imaging. In order to spatially localize the metabolic activity, the metabolic image must be interpolated to the size of the proton image. The most common methods for choosing the unknown values rely exclusively on values of the original un-interpolated image. In this work, we present an alternative method that uses the higher-resolution proton image to provide additional spatial structure. The interpolated image is the result of a convex optimization algorithm which is solved with the Fast Iterative Shrinkage Threshold Algorithm (FISTA). Results are shown with images of hyperpolarized pyruvate, lactate, and bicarbonate using data of the heart and brain from healthy human volunteers, a healthy porcine heart, and a human with prostate cancer.

preprint2020arXiv

A Metabolite Specific 3D Stack-of-Spiral bSSFP Sequence for Improved Lactate Imaging in Hyperpolarized [1-13^{13}C]Pyruvate Studies on a 3T Clinical Scanner

Purpose: The balanced steady-state free precession sequence has been previously explored to improve the efficient use of non-recoverable hyperpolarized 13^{13}C magnetization, but suffers from poor spectral selectivity and long acquisition time. The purpose of this study was to develop a novel metabolite-specific 3D bSSFP ("MS-3DSSFP") sequence with stack-of-spiral readouts for improved lactate imaging in hyperpolarized [1-13^{13}C]pyruvate studies on a clinical 3T scanner. Methods: Simulations were performed to evaluate the spectral response of the MS-3DSSFP sequence. Thermal 13^{13}C phantom experiments were performed to validate the MS-3DSSFP sequence. In vivo hyperpolarized [1-13^{13}C]pyruvate studies were performed to compare the MS-3DSSFP sequence with metabolite specific gradient echo ("MS-GRE") sequences for lactate imaging. Results: Simulations, phantom and in vivo studies demonstrate that the MS-3DSSFP sequence achieved spectrally selective excitation on lactate while minimally perturbing other metabolites. Compared with MS-GRE sequences, the MS-3DSSFP sequence showed approximately a 2.5-fold SNR improvement for lactate imaging in rat kidneys, prostate tumors in a mouse model and human kidneys. Conclusions: Improved lactate imaging using the MS-3DSSFP sequence in hyperpolarized [1-13^{13}C]pyruvate studies was demonstrated in animals and humans. The MS-3DSSFP sequence could be applied for other clinical applications such as in the brain or adapted for imaging other metabolites such as pyruvate and bicarbonate.

preprint2020arXiv

A Regional Bolus Tracking and Real-time B1_1 Calibration Method for Hyperpolarized 13^{13}C MRI

Purpose: Acquisition timing and B1_1 calibration are two key factors that affect the quality and accuracy of hyperpolarized 13^{13}C MRI. The goal of this project was to develop a new approach using regional bolus tracking to trigger Bloch-Siegert B1_1 mapping and real-time B1_1 calibration based on regional B1_1 measurements, followed by dynamic imaging of hyperpolarized 13C^{13}C metabolites in vivo. Methods: The proposed approach was implemented on a system which allows real-time data processing and real-time control on the sequence. Real-time center frequency calibration upon the bolus arrival was also added. The feasibility of applying the proposed framework for in vivo hyperpolarized 13^{13}C imaging was tested on healthy rats, tumor-bearing mice and a healthy volunteer on a clinical 3T scanner following hyperpolarized [1-13^{13}C]pyruvate injection. Multichannel receive coils were used in the human study. Results: Automatic acquisition timing based on either regional bolus peak or bolus arrival was achieved with the proposed framework. Reduced blurring artifacts in real-time reconstructed images were observed with real-time center frequency calibration. Real-time computed B1_1 scaling factors agreed with real-time acquired B1_1 maps. Flip angle correction using B1_1 maps results in a more consistent quantification of metabolic activity (i.e, pyruvate-to-lactate conversion, kPL_{PL}). Experiment recordings are provided to demonstrate the real-time actions during the experiment. Conclusion: The proposed method was successfully demonstrated on animals and a human volunteer, and is anticipated to improve the efficient use of the hyperpolarized signal as well as the accuracy and robustness of hyperpolarized 13^{13}C imaging.

preprint2020arXiv

Bone Material Analogues for PET/MRI Phantoms

Purpose: To develop bone material analogues that can be used in construction of phantoms for simultaneous PET/MRI systems. Methods: Plaster was used as the basis for the bone material analogues tested in this study. It was mixed with varying concentrations of an iodinated CT contrast, a gadolinium-based MR contrast agent, and copper sulfate to modulate the attenuation properties and MRI properties (T1 and T2*). Attenuation was measured with CT and 68Ge transmission scans, and MRI properties were measured with quantitative ultrashort echo time pulse sequences. A proof-of-concept skull was created by plaster casting. Results: Undoped plaster has a 511 keV attenuation coefficient (~0.14 cm-1) similar to cortical bone (0.10-0.15 cm-1), but slightly longer T1 (~500 ms) and T2* (~1.2 ms) MR parameters compared to bone (T1 ~ 300 ms, T2* ~ 0.4 ms). Doping with the iodinated agent resulted in increased attenuation with minimal perturbation to the MR parameters. Doping with a gadolinium chelate greatly reduced T1 and T2*, resulting in extremely short T1 values when the target T2* values were reached, while the attenuation coefficient was unchanged. Doping with copper sulfate was more selective for T2* shortening and achieved comparable T1 and T2* values to bone (after 1 week of drying), while the attenuation coefficient was unchanged. Conclusions: Plaster doped with copper sulfate is a promising bone material analogue for a PET/MRI phantom, mimicking the MR properties (T1 and T2*) and 511 keV attenuation coefficient of human cortical bone.

preprint2020arXiv

Extreme MRI: Large-Scale Volumetric Dynamic Imaging from Continuous Non-Gated Acquisitions

Purpose: To develop a framework to reconstruct large-scale volumetric dynamic MRI from rapid continuous and non-gated acquisitions, with applications to pulmonary and dynamic contrast enhanced (DCE) imaging. Theory and Methods: The problem considered here requires recovering hundred-gigabytes of dynamic volumetric image data from a few gigabytes of k-space data, acquired continuously over several minutes. This reconstruction is vastly under-determined, heavily stressing computing resources as well as memory management and storage. To overcome these challenges, we leverage intrinsic three dimensional (3D) trajectories, such as 3D radial and 3D cones, with ordering that incoherently cover time and k-space over the entire acquisition. We then propose two innovations: (1) A compressed representation using multi-scale low rank matrix factorization that constrains the reconstruction problem, and reduces its memory footprint. (2) Stochastic optimization to reduce computation, improve memory locality, and minimize communications between threads and processors. We demonstrate the feasibility of the proposed method on DCE imaging acquired with a golden-angle ordered 3D cones trajectory and pulmonary imaging acquired with a bit-reversed ordered 3D radial trajectory. We compare it with "soft-gated" dynamic reconstruction for DCE and respiratory resolved reconstruction for pulmonary imaging. Results: The proposed technique shows transient dynamics that are not seen in gating based methods. When applied to datasets with irregular, or non-repetitive motions, the proposed method displays sharper image features. Conclusion: We demonstrated a method that can reconstruct massive 3D dynamic image series in the extreme undersampling and extreme computation setting.

preprint2020arXiv

multiMap: A Gradient Spoiled Sequence for Simultaneously Measuring B1+, B0, T1/M0, T2, T2*, and Fat Fraction of a Slice

We propose multiMap, a single scan that can generate several quantitative maps simultaneously. The sequence acquires multiple images in a time-efficient manner, which can be modeled for T_2, T2*, main- and transmit-field inhomogeneity, T_1:equilibrium magnetization, and water and fat content. The sequence is constructed so that cross-dependencies between parameters are isolated. Thus, each parameter can be estimated independently. Estimates of all parameters are shown on bottle phantoms, the brain, and the knee. The results are compared to estimates from established techniques.

preprint2015arXiv

Imaging Renal Urea Handling in Rats at Millimeter Resolution using Hyperpolarized Magnetic Resonance Relaxometry

\textit{In vivo} spin spin relaxation time (T2T_2) heterogeneity of hyperpolarized \textsuperscript{13}C urea in the rat kidney was investigated. Selective quenching of the vascular hyperpolarized \textsuperscript{13}C signal with a macromolecular relaxation agent revealed that a long-T2T_2 component of the \textsuperscript{13}C urea signal originated from the renal extravascular space, thus allowing the vascular and renal filtrate contrast agent pools of the \textsuperscript{13}C urea to be distinguished via multi-exponential analysis. The T2T_2 response to induced diuresis and antidiuresis was performed with two imaging agents: hyperpolarized \textsuperscript{13}C urea and a control agent hyperpolarized bis-1,1-(hydroxymethyl)-1-\textsuperscript{13}C-cyclopropane-2H8^2\textrm{H}_8. Large T2T_2 increases in the inner-medullar and papilla were observed with the former agent and not the latter during antidiuresis suggesting that T2T_2 relaxometry may be used to monitor the inner-medullary urea transporter (UT)-A1 and UT-A3 mediated urea concentrating process. Two high resolution imaging techniques - multiple echo time averaging and ultra-long echo time sub-2 mm3^3 resolution 3D imaging - were developed to exploit the particularly long relaxation times observed.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.