Cardiac motion and spillover correction for quantitative PET imaging using dynamic MRI.

Accurate fusion of late gadolinium enhancement magnetic resonance imaging (MRI) and electro-anatomical voltage mapping (EAM) is required to evaluate the potential of MRI to identify the substrate of ventricular tachycardia. However, both datasets are not acquired at the same cardiac phase and EAM data is corrupted with respiratory motion limiting the accuracy of current rigid fusion techniques. Knowledge of cardiac and respiratory motion during EAM is thus required to enhance the fusion process. In this study, we propose a novel approach to characterize both cardiac and respiratory motion from EAM data using the temporal evolution of the 3D catheter location recorded from clinical EAM systems. Cardiac and respiratory motion components are extracted from the recorded catheter location using multi-band filters. Filters are calibrated for each EAM point using estimates of heart rate and respiratory rate. The method was first evaluated in numerical simulations using 3D models of cardiac and respiratory motions of the heart generated from real time MRI data acquired in 5 healthy subjects. An accuracy of 0.6–0.7 mm was found for both cardiac and respiratory motion estimates in numerical simulations. Cardiac and respiratory motions were then characterized in 27 patients who underwent LV mapping for treatment of ventricular tachycardia. Mean maximum amplitude of cardiac and respiratory motion was 10.2±2.7 mm (min = 5.5, max = 16.9) and 8.8±2.3 mm (min = 4.3, max = 14.8), respectively. 3D Cardiac and respiratory motions could be estimated from the recorded catheter location and the method does not rely on additional imaging modality such as X-ray fluoroscopy and can be used in conventional electrophysiology laboratory setting.

Cardiac stereotactic body radiation therapy for refractory ventricular tachycardia: Impact of cardiac and respiratory movement on the target volume definition

Cardiac applications in radiation therapy are rapidly expanding including magnetic resonance guided radiation therapy (MRgRT) for real-time gating for targeting and avoidance near the heart or treating ventricular tachycardia (VT). This work describes the development and implementation of a novel multi-modality and magnetic resonance (MR)-compatible cardiac phantom. The patient-informed 3D model was derived from manual contouring of a contrast-enhanced Coronary Computed Tomography Angiography scan, exported as a Stereolithography model, then post-processed to simulate female heart with an average volume. The model was 3D-printed using Elastic50A to provide MR contrast to water background. Two rigid acrylic modules containing cardiac structures were designed and assembled, retrofitting to an MR-safe programmable motor to supply cardiac and respiratory motion in superior-inferior directions. One module contained a cavity for an ion chamber (IC), and the other was equipped with multiple interchangeable cavities for plastic scintillation detectors (PSDs). Images were acquired on a 0.35 T MR-linac for validation of phantom geometry, motion, and simulated online treatment planning and delivery. Three motion profiles were prescribed: patient-derived cardiac (sine waveform, 4.3mm peak-to-peak, 60 beats/min), respiratory (cos4 waveform, 30mm peak-to-peak, 12 breaths/min), and a superposition of cardiac (sine waveform, 4mm peak-to-peak, 70 beats/min) and respiratory (cos4 waveform, 24mm peak-to-peak, 12 breaths/min). The amplitude of the motion profiles was evaluated from sagittal cine images at eight frames/s with a resolution of 2.4mm × 2.4mm. Gated dosimetry experiments were performed using the two module configurations for calculating dose relative to stationary. A CT-based VT treatment plan was delivered twice under cone-beam CT guidance and cumulative stationary doses to multi-point PSDs were evaluated. No artifacts were observed on any images acquired during phantom operation. Phantom excursions measured 49.3±25.8%/66.9±14.0%, 97.0±2.2%/96.4±1.7%, and 90.4±4.8%/89.3±3.5% of prescription for cardiac, respiratory, and cardio-respiratory motion profiles for the 2-chamber (PSD) and 12-substructure (IC) phantom modules respectively. In the gated experiments, the cumulative dose was<2% from expected using the IC module. Real-time dose measured for the PSDs at 10Hz acquisition rate demonstrated the ability to detect the dosimetric consequences of cardiac, respiratory, and cardio-respiratory motion when sampling of different locations during a single delivery, and the stability of our phantom dosimetric results over repeated cycles for the high dose and high gradient regions. For the VT delivery, high dose PSD was<1% from expected (5-6 cGy deviation of 5.9Gy/fraction) and high gradient/low dose regions had deviations<3.6% (6.3 cGy less than expected 1.73Gy/fraction). A novel multi-modality modular heart phantom was designed, constructed, and used for gated radiotherapy experiments on a 0.35 T MR-linac. Our phantom was capable of mimicking cardiac, cardio-respiratory, and respiratory motion while performing dosimetric evaluations of gated procedures using IC and PSD configurations. Time-resolved PSDs with small sensitive volumes appear promising for low-amplitude/high-frequency motion and multi-point data acquisition for advanced dosimetric capabilities. Illustrating VT planning and delivery further expands our phantom to address the unmet needs of cardiac applications in radiotherapy.

Catheter ablation is widely accepted as the best remaining option for the treatment of atrial fibrillation if drug therapy fails. Ablation procedures can be guided by 3-D overlay images projected onto live fluoroscopic X-ray images. These overlay images are generated from either MR, CT or C-Arm CT volumes. As the alignment of the overlay is often compromised by cardiac and respiratory motion, motion compensation methods are desirable. The most recent and promising approaches use either a catheter in the coronary sinus vein, or a circumferential mapping catheter placed at the ostium of one of the pulmonary veins. As both methods suffer from different problems, we propose a novel method to achieve motion compensation for fluoroscopy guided cardiac ablation procedures. Our new method localizes the coronary sinus catheter. Based on this information, we estimate the position of the circumferential mapping catheter. As the mapping catheter is placed at the site of ablation, it provides a good surrogate for respiratory and cardiac motion. To correlate the motion of both catheters, our method includes a training phase in which both catheters are tracked together. The training information is then used to estimate the cardiac and respiratory motion of the left atrium by observing the coronary sinus catheter only. The approach yields an average 2-D estimation error of 1.99 ± 1.20 mm.KeywordsPulmonary VeinLeft AtriumRespiratory MotionAblation ProcedureMotion CompensationThese keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

To obtain multiphase cardiac cine images with high resolution, a novel self-gating method for both cardiac and respiratory motions is proposed. The proposed method uses the phase of projection data obtained from a separate axial slice to measure cardiac and respiratory motion, after the acquisition of every k-space line in the image plane. Cardiac motion is estimated from the phase of the projection data passing through the aorta, which is amplified by superior-inferior directional bipolar gradients, whereas respiratory motion is estimated from the phase of the left-right directional projection data of the abdomen. To verify the proposed self-gating method, a simulation and in vivo steady state free precession cardiac imaging were performed. The proposed method provides high resolution multiphase cardiac cine images. Using the proposed self-gating method, the phase variation of the projection data offers information about cardiac and respiratory motions that is equivalent to external gating devices. The proposed method can capture time-resolved cardiac and respiratory motion from the phase information of the projection data. Because the projection data is obtained from a separate gating slice, the self-gating signals are not affected by imaging planes. Magn Reson Med 77:1216-1222, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Respiratory motion (RM) and cardiac motion (CM) degrade the quality and resolution in cardiac PET scans. We have developed non-rigid motion estimation methods to estimate both RM and CM based on 4D cardiac gated PET data alone, and compensate the dual respiratory and cardiac (R&C) motions after (MCAR), during (MCDR), and before (MCBR) image reconstruction. In all three R&C motion correction methods, attenuation-activity mismatch effect was modeled by using transformed attenuation maps using the estimated RM. The difference of using activity preserving and non-activity preserving models in R&C correction was also studied. Realistic Monte Carlo simulated 4D cardiac PET data using the 4D XCAT phantom and accurate models of the scanner design parameters and performance characteristics at different noise levels were employed as the known truth and for method development and evaluation. Results from the simulation study suggested that all three dual R&C motion correction methods provide substantial improvement in the quality of 4D cardiac gated PET images as compared with no motion correction. Specifically, the MCDR method yields the best performance for all different noise levels compared with the MCAR and MCBR methods. While MCBR reduces computational time dramatically but the resultant 4D cardiac gated PET images has overall inferior image quality when compared to that from the MCAR and MCDR approaches in the ‘almost’ noise free case. Also, the MCBR method has better noise handling properties when compared with MCAR and provides better quantitative results in high noise cases. When the goal is to reduce scan time or patient radiation dose, MCDR and MCBR provide a good compromise between image quality and computational times.

Gastrointestinal (GI) tract motility is one of the main sources for intra/inter-fraction variability and uncertainty in radiation therapy for abdominal targets. Models for GI motility can improve the assessment of delivered dose and contribute to the development, testing, and validation of deformable image registration (DIR) and dose-accumulation algorithms. To implement GI tract motion in the 4D extended cardiac-torso (XCAT) digital phantom of human anatomy. Motility modes that exhibit large amplitude changes in the diameter of the GI tract and may persist over timescales comparable to online adaptive planning and radiotherapy delivery were identified based on literature research. Search criteria included amplitude changes larger than planning risk volume expansions and durations of the order of tens of minutes. The following modes were identified: peristalsis, rhythmic segmentation, high amplitude propagating contractions (HAPCs), and tonic contractions. Peristalsis and rhythmic segmentations were modeled by traveling and standing sinusoidal waves. HAPCs and tonic contractions were modeled by traveling and stationary Gaussian waves. Wave dispersion in the temporal and spatial domain was implemented by linear, exponential, and inverse power law functions. Modeling functions were applied to the control points of the nonuniform rational B-spline surfaces defined in the reference XCAT library. GI motility was combined with the cardiac and respiratory motions available in the standard 4D-XCAT phantom. Default model parameters were estimated based on the analysis of cine MRI acquisitions in 10 patients treated in a 1.5T MR-linac. We demonstrate the ability to generate realistic 4D multimodal images that simulate GI motility combined with respiratory and cardiac motion. All modes of motility, except tonic contractions, were observed in the analysis of our cine MRI acquisitions. Peristalsis was the most common. Default parameters estimated from cine MRI were used as initial values for simulation experiments. It is shown that in patients undergoing stereotactic body radiotherapy for abdominal targets, the effects of GI motility can be comparable or larger than the effects of respiratory motion. The digital phantom provides realistic models to aid in medical imaging and radiation therapy research. The addition of GI motility will further contribute to the development, testing, and validation of DIR and dose accumulation algorithms for MR-guided radiotherapy.

Rotational coronary angiography enables 3D reconstruction but suffers from intra-scan cardiac and respiratory motion. While gating handles cardiac motion, respiratory motion requires compensation. State-of-the-art algorithms rely on 3D-2D registration that depends on initial reconstructions of sufficient quality. We propose a compensation method that is applied directly in projection domain. It overcomes the need for reconstruction and thus complements the state-of-the-art. Virtual single-frame background subtraction based on vessel segmentation and spectral deconvolution yields non-truncated images of the contrasted lumen. This allows motion compensation based on data consistency conditions. We compensate craniocaudal shifts by optimizing epipolar consistency to (a) devise an image-based surrogate for cardiac motion and (b) compensate for respiratory motion. We validate our approach in two numerical phantom studies and three clinical cases. Correlation of the image-based surrogate for cardiac motion with the ECG-based ground truth was excellent yielding a Pearson correlation of 0.93±0.04. Considering motion compensation, the target error measure decreased by 98% and 69%, respectively, for the phantom experiments while for the clinical cases the same figure of merit improved by 46±21%. The proposed method is entirely image-based and accurately estimates craniocaudal shifts due to respiration and cardiac contraction. Future work will investigate experimental trajectories and possibilities for simplification of the single-frame subtraction pipeline.

Background Electro-anatomical voltage mapping (EAM) is an invasive technique used for the identification of ventricular tachycardia (VT) substrate and subsequent guidance of VT ablation [1]. The mapping of VT substrate is very timeconsuming procedure, requires highly skilled electrophysiologist, is associated with patient risk and is an invasive procedure. Late gadolinium enhancement (LGE) MRI allows non-invasive evaluation of 3D structure of scar. Therefore, LGE has the potential to identify the VT substrate and can now be integrated in the current clinical platform for guidance of VT ablation as a roadmap. However, fusion of the two imaging modality is very challenging due to respiratory and cardiac motion during the mapping which results in large errors in data fusion. Our aim in this study is to develop a novel algorithm to detect the respiratory and cardiac-induced motion of the mapping catheter during the VT ablation to facilitate integration of LGE MRI to EAM data. Methods

Effects of Intrafractional Motion on Water Equivalent Pathlength in Respiratory-Gated Heavy Charged Particle Beam Radiotherapy

Subject motion including cardiac and respiration motion during data acquisition can cause degraded image reconstruction. Time-of-flight (TOF) in PET provides additional useful information, which can be utilized for motion estimation. We aim to develop accurate motion field estimation and motion compensation from TOF PET data fully utilizing TOF information for improved quantitative PET imaging. In this work, we develop a new motion displacement field estimation method using TOF PET histoimages. The TOF list-mode data can be efficiently deposited into a TOF histoimage, and radioactive tracer distribution can be efficiently reconstructed using the DIRECT (Direct Image REConstruction for TOF) approaches since TOF histoimage has the same spatial coordinates and geometry as the reconstructed image. First we derive the optical flow equation for TOF histoimage, from which we formulate the global energy function including the total squared error of data consistency based on the optical flow equation and a displacement regularizer. We then obtain Euler–Lagrange equation of TOF histoimage by minimizing the global energy function based on the principle of least action. Inspired by Horn–Schunck method, we obtain an iterative algorithm for motion displacement field estimation, and displacement smoothness regularization is naturally included. We simulate a generic 2D TOF PET scanner using XCAT phantom with cardiac motion, and show accurate motion displacement field can be obtained from both noise-free and noisy TOF histoimages using the proposed method. In summary, we present a motion displacement field estimation method directly using TOF histoimages, which can be important for motion compensation and improved dynamic image reconstructions.

This paper describes an acquisition and reconstruction strategy for cardiac cine MRI that does not require the use of electrocardiogram or breath holding. The method has similarities with self-gated techniques as information about cardiac and respiratory motion is derived from the imaging sequence itself; here, by acquiring the center k-space line at the beginning of each segment of a balanced steady-state free precession sequence. However, the reconstruction step is fundamentally different: a generalized reconstruction by inversion of coupled systems is used instead of conventional gating. By correcting for nonrigid cardiac and respiratory motion, generalized reconstruction by inversion of coupled systems (GRICS) uses all acquired data, whereas gating rejects data acquired in certain motion states. The method relies on the processing and analysis of the k-space central line data: local information from a 32-channel cardiac coil is used in order to automatically extract eigenmodes of both cardiac and respiratory motion. In the GRICS framework, these eigenmodes are used as driving signals of a motion model. The motion model is defined piecewise, so that each cardiac phase is reconstructed independently. Results from six healthy volunteers, with various slice orientations, show improved image quality compared to combined respiratory and cardiac gating.

Image Guided, Motion-free Patient Body Setup using 3D Volumetric Image Registration of Classified Stable Bony Landmarks

To develop a novel framework for free-breathing MRI called XD-GRASP, which sorts dynamic data into extra motion-state dimensions using the self-navigation properties of radial imaging and reconstructs the multidimensional dataset using compressed sensing. Radial k-space data are continuously acquired using the golden-angle sampling scheme and sorted into multiple motion-states based on respiratory and/or cardiac motion signals derived directly from the data. The resulting undersampled multidimensional dataset is reconstructed using a compressed sensing approach that exploits sparsity along the new dynamic dimensions. The performance of XD-GRASP is demonstrated for free-breathing three-dimensional (3D) abdominal imaging, two-dimensional (2D) cardiac cine imaging and 3D dynamic contrast-enhanced (DCE) MRI of the liver, comparing against reconstructions without motion sorting in both healthy volunteers and patients. XD-GRASP separates respiratory motion from cardiac motion in cardiac imaging, and respiratory motion from contrast enhancement in liver DCE-MRI, which improves image quality and reduces motion-blurring artifacts. XD-GRASP represents a new use of sparsity for motion compensation and a novel way to handle motions in the context of a continuous acquisition paradigm. Instead of removing or correcting motion, extra motion-state dimensions are reconstructed, which improves image quality and also offers new physiological information of potential clinical value.

Data-driven respiratory gating techniques in cardiac PET extract respiratory motion signal from the PET data to guide respiratory gating and motion estimation. However, the influence of the quality of the list-mode data including the uptake ratio between the myocardium and background and the detected count level on its performance has never been studied thoroughly. In this project, list-mode PET data derived from realistic Monte Carlo simulation of different uptake ratios and count levels was employed to quantitatively evaluate the respiratory motion signal extraction and estimation techniques. Simulated projection data were generated from the 4D XCAT phantom with known respiratory motion. The list-mode data were rebinned into projection data at 200 msec intervals and the centers-of-mass of a region-of-interest over the heart were calculated. Contributions from the background were subtracted in the centroid calculation. The peak in the Fourier transform of the time sequence of the centroid locations, or the frequency spectrum, was identified as the respiratory motion signal and its SNR was measured. The high frequency components of the frequency spectrum were removed for noise smoothing before extracting the respiratory motion signal for use in respiratory gating estimation. The results show that with background correction, the amplitude of the estimated respiratory motion signal is increased and is closer to the truth. Both lower myocardium to background uptake ratio and lower count level reduce the SNR of the respiratory motion signal. When they become too low, the extraction of the respiratory motion signal becomes difficult or fails completely.