Volumetric Parameterization of the Placenta to a Flattened Template.

We present a volumetric mesh-based algorithm for parameterizing the placenta to a flattened template to enable effective visualization of local anatomy and function. MRI shows potential as a research tool as it provides signals directly related to placental function. However, due to the curved and highly variable in vivo shape of the placenta, interpreting and visualizing these images is difficult. We address interpretation challenges by mapping the placenta so that it resembles the familiar ex vivo shape. We formulate the parameterization as an optimization problem for mapping the placental shape represented by a volumetric mesh to a flattened template. We employ the symmetric Dirichlet energy to control local distortion throughout the volume. Local injectivity in the mapping is enforced by a constrained line search during the gradient descent optimization. We validate our method using a research study of 111 placental shapes extracted from BOLD MRI images. Our mapping achieves sub-voxel accuracy in matching the template while maintaining low distortion throughout the volume. We demonstrate how the resulting flattening of the placenta improves visualization of anatomy and function. Our code is freely available at https://github.com/mabulnaga/placenta-flattening.

Quantitative T1 and T2 mapping by magnetic resonance fingerprinting (MRF) of the placenta before and after maternal hyperoxia.

INTRODUCTION: MR relaxometry has been used to assess placental exchange function, but methods to date are not sufficiently fast to be robust to placental motion. Magnetic resonance fingerprinting (MRF) permits rapid, voxel-wise, intrinsically co-registered T1 and T2 mapping. After characterizing measurement error, we scanned pregnant women during air and oxygen breathing to demonstrate MRF's ability to detect placental oxygenation changes. METHODS: The accuracy of FISP-based, sliding-window reconstructed MRF was tested on phantoms. MRF scans in 9-s breath holds were acquired at 3T in 31 pregnant women during air and oxygen breathing. A mixed effects model was used to test for changes in placenta relaxation times between physiological states, to assess the dependency on gestational age (GA), and the impact of placental motion. RESULTS: MRF estimates of known phantom relaxation times resulted in mean absolute errors for T1 of 92 ms (4.8%), but T2 was less accurate at 16 ms (13.6%). During normoxia, placental T1 = 1825 ± 141 ms (avg ± standard deviation) and T2 = 60 ± 16 ms (gestational age range 24.3-36.7, median 32.6 weeks). In the statistical model, placental T2 rose and T1 remained contant after hyperoxia, and no GA dependency was observed for T1 or T2. DISCUSSION: Well-characterized, motion-robust MRF was used to acquire T1 and T2 maps of the placenta. Changes with hyperoxia are consistent with a net increase in oxygen saturation. Toward the goal of whole-placenta quantitative oxygenation imaging over time, we aim to implement 3D MRF with integrated motion correction to improve T2 accuracy.

Rapid head-pose detection for automated slice prescription of fetal-brain MRI.

In fetal-brain MRI, head-pose changes between prescription and acquisition present a challenge to obtaining the standard sagittal, coronal and axial views essential to clinical assessment. As motion limits acquisitions to thick slices that preclude retrospective resampling, technologists repeat ~55-second stack-of-slices scans (HASTE) with incrementally reoriented field of view numerous times, deducing the head pose from previous stacks. To address this inefficient workflow, we propose a robust head-pose detection algorithm using full-uterus scout scans (EPI) which take ~5 seconds to acquire. Our ~2-second procedure automatically locates the fetal brain and eyes, which we derive from maximally stable extremal regions (MSERs). The success rate of the method exceeds 94% in the third trimester, outperforming a trained technologist by up to 20%. The pipeline may be used to automatically orient the anatomical sequence, removing the need to estimate the head pose from 2D views and reducing delays during which motion can occur.

STRESS: Super-Resolution for Dynamic Fetal MRI using Self-Supervised Learning.

Fetal motion is unpredictable and rapid on the scale of conventional MR scan times. Therefore, dynamic fetal MRI, which aims at capturing fetal motion and dynamics of fetal function, is limited to fast imaging techniques with compromises in image quality and resolution. Super-resolution for dynamic fetal MRI is still a challenge, especially when multi-oriented stacks of image slices for oversampling are not available and high temporal resolution for recording the dynamics of the fetus or placenta is desired. Further, fetal motion makes it difficult to acquire high-resolution images for supervised learning methods. To address this problem, in this work, we propose STRESS (Spatio-Temporal Resolution Enhancement with Simulated Scans), a self-supervised super-resolution framework for dynamic fetal MRI with interleaved slice acquisitions. Our proposed method simulates an interleaved slice acquisition along the high-resolution axis on the originally acquired data to generate pairs of low- and high-resolution images. Then, it trains a super-resolution network by exploiting both spatial and temporal correlations in the MR time series, which is used to enhance the resolution of the original data. Evaluations on both simulated and in utero data show that our proposed method outperforms other self-supervised super-resolution methods and improves image quality, which is beneficial to other downstream tasks and evaluations.

Placental MRI: Development of an MRI compatible ex vivo system for whole placenta dual perfusion.

PURPOSE: Placental dysfunction plays a key role in diseases that affect the fetus in utero and after birth. Aiming to develop a platform for validating in vivo placental MRI and investigations into placental physiology, we designed and built a prototype MRI-compatible perfusion chamber with an integrated MRI receive coil for high SNR ex vivo placental imaging. PRINCIPAL RESULTS: After optimizing placenta vascular clearing and perfusion protocols, we performed contrast enhanced MR angiography and MR relaxometry on eight carefully selected placentas while they were perfused via the umbilical arteries (UAs). Additionally, two of these placentas underwent maternal perfusion via the intervillous space (IVS). Despite striving for homogenous perfusion across the whole placenta, imaging results were highly heterogeneous for both UA and IVS perfused placentas. By histology, we observed blood congestion in the villi in regions that showed low UA perfusion during MRI. In two placentas prominent chorionic arteries followed by adjacent veins underwent contrast enhancement in the absence of villous capillary blush. The single placenta from a pregnancy affected by IUGR had the most homogeneous villous capillary perfusion. MAJOR CONCLUSIONS: A dual perfusion system for ex vivo placentas compatible with MRI permitted assessment of UA and IVS placental perfusion. We observed spatial UA perfusion heterogeneity and evidence for arteriovenous shunting in placentas from normal pregnancies and deliveries, but relative villous capillary perfusion homogeneity in a single IUGR placenta. Future work will focus on system optimization, followed by physiological manipulation and validation of in vivo placental MRI.

Semi-Supervised Learning for Fetal Brain MRI Quality Assessment with ROI consistency.

Fetal brain MRI is useful for diagnosing brain abnormalities but is challenged by fetal motion. The current protocol for T2-weighted fetal brain MRI is not robust to motion so image volumes are degraded by inter- and intra-slice motion artifacts. Besides, manual annotation for fetal MR image quality assessment are usually time-consuming. Therefore, in this work, a semi-supervised deep learning method that detects slices with artifacts during the brain volume scan is proposed. Our method is based on the mean teacher model, where we not only enforce consistency between student and teacher models on the whole image, but also adopt an ROI consistency loss to guide the network to focus on the brain region. The proposed method is evaluated on a fetal brain MR dataset with 11,223 labeled images and more than 200,000 unlabeled images. Results show that compared with supervised learning, the proposed method can improve model accuracy by about 6% and outperform other state-of-the-art semi-supervised learning methods. The proposed method is also implemented and evaluated on an MR scanner, which demonstrates the feasibility of online image quality assessment and image reacquisition during fetal MR scans.

Enhanced detection of fetal pose in 3D MRI by Deep Reinforcement Learning with physical structure priors on anatomy.

Fetal MRI is heavily constrained by unpredictable and substantial fetal motion that causes image artifacts and limits the set of viable diagnostic image contrasts. Current mitigation of motion artifacts is predominantly performed by fast, single-shot MRI and retrospective motion correction. Estimation of fetal pose in real time during MRI stands to benefit prospective methods to detect and mitigate fetal motion artifacts where inferred fetal motion is combined with online slice prescription with low-latency decision making. Current developments of deep reinforcement learning (DRL), offer a novel approach for fetal landmarks detection. In this task 15 agents are deployed to detect 15 landmarks simultaneously by DRL. The optimization is challenging, and here we propose an improved DRL that incorporates priors on physical structure of the fetal body. First, we use graph communication layers to improve the communication among agents based on a graph where each node represents a fetal-body landmark. Further, additional reward based on the distance between agents and physical structures such as the fetal limbs is used to fully exploit physical structure. Evaluation of this method on a repository of 3-mm resolution in vivo data demonstrates a mean accuracy of landmark estimation within 10 mm of ground truth as 87.3%, and a mean error of 6.9 mm. The proposed DRL for fetal pose landmark search demonstrates a potential clinical utility for online detection of fetal motion that guides real-time mitigation of motion artifacts as well as health diagnosis during MRI of the pregnant mother.

Computer-Vision Techniques for Water-Fat Separation in Ultra High-Field MRI Local Specific Absorption Rate Estimation.

OBJECTIVE: The purpose of this paper is to prove that computer-vision techniques allow synthesizing water-fat separation maps for local specific absorption rate (SAR) estimation, when patient-specific water-fat images are not available. METHODS: We obtained ground truth head models by using patient-specific water-fat images. We obtained two different label-fusion water-fat models generating a water-fat multiatlas and applying the STAPLE and local-MAP-STAPLE label-fusion methods. We also obtained patch-based water-fat models applying a local group-wise weighted combination of the multiatlas. Electromagnetic (EM) simulations were performed, and B1+ magnitude and 10 g averaged SAR maps were generated. RESULTS: We found local approaches provide a high DICE overlap (72.6 ± 10.2% fat and 91.6 ± 1.5% water in local-MAP-STAPLE, and 68.8 ± 8.2% fat and 91.1 ± 1.0% water in patch-based), low Hausdorff distances (18.6 ± 7.7 mm fat and 7.4 ± 11.2 mm water in local-MAP-STAPLE, and 16.4 ± 8.5 mm fat and 7.2 ± 11.8 mm water in patch-based) and a low error in volume estimation (15.6 ± 34.4% fat and 5.6 ± 4.1% water in the local-MAP-STAPLE, and 14.0 ± 17.7% fat and 4.7 ± 2.8% water in patch-based). The positions of the peak 10 g-averaged local SAR hotspots were the same for every model. CONCLUSION: We have created patient-specific head models using three different computer-vision-based water-fat separation approaches and compared the predictions of B1+ field and SAR distributions generated by simulating these models. Our results prove that a computer-vision approach can be used for patient-specific water-fat separation, and utilized for local SAR estimation in high-field MRI. SIGNIFICANCE: Computer-vision approaches can be used for patient-specific water-fat separation and for patient specific local SAR estimation, when water-fat images of the patient are not available.

Placental Flattening via Volumetric Parameterization.

We present a volumetric mesh-based algorithm for flattening the placenta to a canonical template to enable effective visualization of local anatomy and function. Monitoring placental function in vivo promises to support pregnancy assessment and to improve care outcomes. We aim to alleviate visualization and interpretation challenges presented by the shape of the placenta when it is attached to the curved uterine wall. To do so, we flatten the volumetric mesh that captures placental shape to resemble the well-studied ex vivo shape. We formulate our method as a map from the in vivo shape to a flattened template that minimizes the symmetric Dirichlet energy to control distortion throughout the volume. Local injectivity is enforced via constrained line search during gradient descent. We evaluate the proposed method on 28 placenta shapes extracted from MRI images in a clinical study of placental function. We achieve sub-voxel accuracy in mapping the boundary of the placenta to the template while successfully controlling distortion throughout the volume. We illustrate how the resulting mapping of the placenta enhances visualization of placental anatomy and function. Our implementation is freely available at https://github.com/mabulnaga/placenta-flattening.

Fetal Pose Estimation in Volumetric MRI using a 3D Convolution Neural Network.

The performance and diagnostic utility of magnetic resonance imaging (MRI) in pregnancy is fundamentally constrained by fetal motion. Motion of the fetus, which is unpredictable and rapid on the scale of conventional imaging times, limits the set of viable acquisition techniques to single-shot imaging with severe compromises in signal-to-noise ratio and diagnostic contrast, and frequently results in unacceptable image quality. Surprisingly little is known about the characteristics of fetal motion during MRI and here we propose and demonstrate methods that exploit a growing repository of MRI observations of the gravid abdomen that are acquired at low spatial resolution but relatively high temporal resolution and over long durations (10-30 minutes). We estimate fetal pose per frame in MRI volumes of the pregnant abdomen via deep learning algorithms that detect key fetal landmarks. Evaluation of the proposed method shows that our framework achieves quantitatively an average error of 4.47 mm and 96.4% accuracy (with error less than 10 mm). Fetal pose estimation in MRI time series yields novel means of quantifying fetal movements in health and disease, and enables the learning of kinematic models that may enhance prospective mitigation of fetal motion artifacts during MRI acquisition.

Radiofrequency heating studies on anesthetized swine using fractionated dipole antennas at 10.5 T.

PURPOSE: To validate electromagnetic and thermal simulations with in vivo temperature measurements, and to demonstrate a framework that can be used to predict temperature increase caused by radiofrequency (RF) excitation with dipole transmitter arrays. METHODS: Dipole arrays were used to deliver RF energy in the back/neck region of the swine using different RF excitation patterns (n = 2-4 per swine) for heating. The temperature in anesthetized swine (n = 3) was measured using fluoroscopic probes (n = 12) and compared against thermal modeling from animal-specific electromagnetic simulations. RESULTS: Simulated temperature curves were in agreement with the measured data. The root mean square error between simulated and measured temperature rise at all locations (at the end of each RF excitation) is calculated as 0.37°C. The mean experimental temperature rise at the maximum temperature rise locations (averaged over all experiments) is calculated as 2.89°C. The root mean square error between simulated and measured temperature at the maximum temperature rise location is calculated as 0.57°C. (Error values are averaged over all experiments.) CONCLUSIONS: Electromagnetic and thermal simulations were validated with experiments. Thermal effects of RF excitation at 10.5 Tesla with dipoles were investigated. Magn Reson Med 79:479-488, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Spatiotemporal alignment of in utero BOLD-MRI series.

PURPOSE: To present a method for spatiotemporal alignment of in-utero magnetic resonance imaging (MRI) time series acquired during maternal hyperoxia for enabling improved quantitative tracking of blood oxygen level-dependent (BOLD) signal changes that characterize oxygen transport through the placenta to fetal organs. MATERIALS AND METHODS: The proposed pipeline for spatiotemporal alignment of images acquired with a single-shot gradient echo echo-planar imaging includes 1) signal nonuniformity correction, 2) intravolume motion correction based on nonrigid registration, 3) correction of motion and nonrigid deformations across volumes, and 4) detection of the outlier volumes to be discarded from subsequent analysis. BOLD MRI time series collected from 10 pregnant women during 3T scans were analyzed using this pipeline. To assess pipeline performance, signal fluctuations between consecutive timepoints were examined. In addition, volume overlap and distance between manual region of interest (ROI) delineations in a subset of frames and the delineations obtained through propagation of the ROIs from the reference frame were used to quantify alignment accuracy. A previously demonstrated rigid registration approach was used for comparison. RESULTS: The proposed pipeline improved anatomical alignment of placenta and fetal organs over the state-of-the-art rigid motion correction methods. In particular, unexpected temporal signal fluctuations during the first normoxia period were significantly decreased (P < 0.01) and volume overlap and distance between region boundaries measures were significantly improved (P < 0.01). CONCLUSION: The proposed approach to align MRI time series enables more accurate quantitative studies of placental function by improving spatiotemporal alignment across placenta and fetal organs. LEVEL OF EVIDENCE: 1 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:403-412.

Approximate Fourier domain expression for Bloch-Siegert shift.

PURPOSE: In this study, a new simple Fourier domain-based analytical expression for the Bloch-Siegert (BS) shift-based B1 mapping method is proposed to obtain |B1+| more accurately while using short BS pulse durations and small off-resonance frequencies. THEORY AND METHODS: A new simple analytical expression for the BS shift is derived by simplifying the Bloch equations. In this expression, the phase is calculated in terms of the Fourier transform of the radiofrequency pulse envelope, and thus making the off- and on-resonance effects more easily understandable. To verify the accuracy of the proposed expression, Bloch simulations and MR experiments are performed for the hard, Fermi, and Shinner-Le Roux pulse shapes. RESULTS: Analyses of the BS phase shift-based B1 mapping method in terms of radiofrequency pulse shape, pulse duration, and off-resonance frequency show that |B1+| can be obtained more accurately with the aid of this new expression. CONCLUSIONS: In this study, a new simple frequency domain analytical expression is proposed for the BS shift. Using this expression, |B1+| values can be predicted from the phase data using the frequency spectrum of the radiofrequency pulse. This method works well even for short pulse durations and small offset frequencies.

Reduction of the radiofrequency heating of metallic devices using a dual-drive birdcage coil.

In this work, it is demonstrated that a dual-drive birdcage coil can be used to reduce the radiofrequency heating of metallic devices during magnetic resonance imaging. By controlling the excitation currents of the two channels of a birdcage coil, the radiofrequency current that is induced near the lead tip could be set to zero. To monitor the current, the image artifacts near the lead tips were measured. The electric field distribution was controlled using a dual-drive birdcage coil. With this method, the lead currents and the lead tip temperatures were reduced substantially [<0.3 °C for an applied 4.4 W/kg SAR compared to >4.9 °C using quadrature excitation], as demonstrated by phantom and animal experiments. The homogeneity of the flip angle distribution was preserved, as shown by volunteer experiments. The normalized root-mean-square error of the flip angle distribution was less than 10% for all excitations. The average specific absorption rate increased as a trade-off for using different excitation patterns.