MRI Investigation of the Linkage Between Respiratory Motion of the Heart and Markers on Patient's Abdomen and Chest: Implications for Respiratory Amplitude Binning List-Mode PET and SPECT Studies

Toward the Use of Implanted Cardiac Leads or the Diaphragm for Active Respiratory Motion Management in Stereotactic Arrhythmia Radioablation

PURPOSE

To investigate the utility of implanted cardiac leads or the diaphragm for active respiratory motion management in stereotactic arrhythmia radioablation by quantifying the relationship between their motions.

METHODS AND MATERIALS

Seven patients treated with stereotactic arrhythmia radioablation were imaged using 5-Hz biplanar, kV x-ray fluoroscopy for 15-20 seconds under both abdominal compression (AC) and free breathing (FB) conditions. Three-dimensional motion traces for different regions of the heart were acquired by tracking and triangulating the position of all implanted cardiac leads. The heart's respiratory motion was extracted from the total motion (respiratory + cardiac) using a low-pass filter and described in optimized coordinates using principal component analysis. The existence of a relationship between the respiratory motion of different cardiac leads or the diaphragm was quantified using the Spearman rank correlation coefficient. Polynomial correlation models relating PC1 cardiac lead motion to the diaphragm were created and evaluated on the resultant errors.

RESULTS

Eighty-one respiratory motion correlations between different positions of the heart or diaphragm were calculated under both AC and FB. Consistently strong correlations between the respiratory motion of different positions in the heart and the diaphragm required accounting for phase shifts between motions. When accounting for phase shifts, the proportion of strong (>0.7) PC1 respiratory motion correlations was 100% under FB and 92.6% under AC. Linear fitting of cardiac lead motion with the diaphragm resulted in mean absolute PC1 tracking errors of (1.0 ± 0.6) mm under FB and (0.7 ± 0.4) mm under AC.

CONCLUSIONS

The respiratory motion of all combinations of implanted cardiac leads and the diaphragm are moderately to strongly correlated after accounting for phase shifts between motion traces. These phase shifts should be carefully considered to ensure patient safety during respiratory tracking or gating during stereotactic arrhythmia radioablation using cardiac leads or the diaphragm as internal surrogates.

Ultrasound-based sensors for respiratory motion assessment in multimodality PET imaging

Breathing motion can displace internal organs by up to several cm; as such, it is a primary factor limiting image quality in medical imaging. Motion can also complicate matters when trying to fuse images from different modalities, acquired at different locations and/or on different days. Currently available devices for monitoring breathing motion often do so indirectly, by detecting changes in the outline of the torso rather than the internal motion itself, and these devices are often fixed to floors, ceilings or walls, and thus cannot accompany patients from one location to another. We have developed small ultrasound-based sensors, referred to as 'organ configuration motion' (OCM) sensors, that attach to the skin and provide rich motion-sensitive information. In the present work we tested the ability of OCM sensors to enable respiratory gating during in vivo PET imaging. A motion phantom involving an FDG solution was assembled, and two cancer patients scheduled for a clinical PET/CT exam were recruited for this study. OCM signals were used to help reconstruct phantom and in vivo data into time series of motion-resolved images. As expected, the motion-resolved images captured the underlying motion. In Patient #1, a single large lesion proved to be mostly stationary through the breathing cycle. However, in Patient #2, several small lesions were mobile during breathing, and our proposed new approach captured their breathing-related displacements. In summary, a relatively inexpensive hardware solution was developed here for respiration monitoring. Because the proposed sensors attach to the skin, as opposed to walls or ceilings, they can accompany patients from one procedure to the next, potentially allowing data gathered in different places and at different times to be combined and compared in ways that account for breathing motion.

Spirometry-based reconstruction of real-time cardiac MRI: Motion control and quantification of heart-lung interactions

PURPOSE

To test the feasibility of cardiac real-time MRI in combination with retrospective gating by MR-compatible spirometry, to improve motion control, and to allow quantification of respiratory-induced changes during free-breathing.

METHODS

Cross-sectional real-time MRI (1.5T; 30 frames/s) using steady-state free precession contrast during free-breathing was combined with MR-compatible spirometry in healthy adult volunteers (n = 4). Retrospective binning assigned images to classes that were defined by electrocardiogram and spirometry. Left ventricular eccentricity index as an indicator of septal position and ventricular volumes in different respiratory phases were calculated to assess heart-lung interactions.

RESULTS

Real-time MRI with MR-compatible spirometry is feasible and well tolerated. Spirometry-based binning improved motion control significantly. The end-diastolic epicardial eccentricity index increased significantly during inspiration (1.04 ± 0.04 to 1.19 ± 0.05; P < .05). During inspiration, right ventricular end-diastolic volume (79 ± 17 mL/m² to 98 ± 18 mL/m² ), stroke volume (41 ± 8 mL/m² to 59 ± 11 mL/m² ) and ejection fraction (53 ± 3% to 60 ± 1%) increased significantly, whereas the end-systolic volume remained almost unchanged. Left ventricular end-diastolic volume, left ventricular stroke volume, and left ventricular ejection fraction decreased during inspiration, whereas the left ventricular end-systolic volume increased. The relationship between stroke volume and end-diastolic volume (Frank-Starling relationship) based on changes induced by respiration allowed for a slope estimate of the Frank-Starling curve to be 0.9 to 1.1.

CONCLUSION

Real-time MRI during free-breathing combined with MR-compatible spirometry and retrospective binning improves image stabilization, allows quantitative image analysis, and importantly, offers unique opportunities to judge heart-lung interactions.

Ultrasound-based sensors to monitor physiological motion

PURPOSE

Medical procedures can be difficult to perform on anatomy that is constantly moving. Respiration displaces internal organs by up to several centimeters with respect to the surface of the body, and patients often have limited ability to hold their breath. Strategies to compensate for motion during diagnostic and therapeutic procedures require reliable information to be available. However, current devices often monitor respiration indirectly, through changes on the outline of the body, and they may be fixed to floors or ceilings, and thus unable to follow a given patient through different locations. Here we show that small ultrasound-based sensors referred to as "organ configuration motion" (OCM) sensors can be fixed to the abdomen and/or chest and provide information-rich, breathing-related signals.

METHODS

By design, the proposed sensors are relatively inexpensive. Breathing waveforms were obtained from tissues at varying depths and/or using different sensor placements. Validation was performed against breathing waveforms derived from magnetic resonance imaging (MRI) and optical tracking signals in five and eight volunteers, respectively.

RESULTS

Breathing waveforms from different modalities were scaled so they could be directly compared. Differences between waveforms were expressed in the form of a percentage, as compared to the amplitude of a typical breath. Expressed in this manner, for shallow tissues, OCM-derived waveforms on average differed from MRI and optical tracking results by 13.1% and 15.5%, respectively.

CONCLUSION

The present results suggest that the proposed sensors provide measurements that properly characterize breathing states. While OCM-based waveforms from shallow tissues proved similar in terms of information content to those derived from MRI or optical tracking, OCM further captured depth-dependent and position-dependent (i.e., chest and abdomen) information. In time, the richer information content of OCM-based waveforms may enable better respiratory gating to be performed, to allow diagnostic and therapeutic equipment to perform at their best.

Clinical evaluation of three respiratory gating schemes for different respiratory patterns on cardiac SPECT

PURPOSE

Respiratory gating reduces respiratory blur in cardiac single photon emission computed tomography (SPECT). It can be implemented as three gating schemes: (a) equal amplitude-based gating (AG); (b) phase or time-based gating (TG); or (c) equal count-based gating (CG), that is, a variant of amplitude-based method. The goal of this study is to evaluate the effectiveness of these respiratory gating methods for patients with different respiratory patterns in myocardial perfusion SPECT.

METHODS

We reviewed 1274 anonymized patient respiratory traces obtained via the Vicon motion-tracking system during their ^99m Tc-sestamibi SPECT scans and grouped them into four breathing categories: (a) regular respiration (RR); (b) periodic respiration (PR); (c) respiration with apnea (AR); and (d) unclassified respiration (UR). For each respiratory pattern, 15 patients were randomly selected and their list-mode data were rebinned using the three gating schemes. A preliminary reconstruction was performed for each gate with the heart region segmented and registered to a reference gate to estimate the respiratory motion. A final reconstruction incorporating respiratory motion correction was done to get a final image set. The estimated respiratory motion, the full-width-at-half-maxima (FWHM) measured across the image intensity profile of the left ventricle wall, as well as the normalized standard deviation measured in a uniform cuboid region of the thorax were analyzed.

RESULTS

There are 47.1%, 24.3%, 13.5%, and 15.1% RR, PR, AR, and UR patients, respectively, among the 1274 patients in this study. The differences among the three gating schemes in RR were smaller than other respiratory patterns. The AG and CG methods showed statistically larger motion estimation than TG particularly in the AR and PR patterns. Noise of AG varied more in different gates, especially for AR and UR patterns.

CONCLUSION

More than half of the patients reviewed exhibited nonregular breathing patterns. Amplitude-based gating, that is, AG and CG, is a preferred gating method for such patterns and is a robust respiratory gating implementation method given the respiratory pattern of the patients is unknown before data acquisition. Phase gating is also a feasible option for regular respiratory pattern.

Recent advances in cardiac SPECT instrumentation and imaging methods

Cardiac SPECT continues to play a critical role in detecting and managing cardiovascular disease, in particularly coronary artery disease (CAD) (Jaarsma et al 2012 J. Am. Coll. Cardiol. 59 1719-28), (Agostini et al 2016 Eur. J. Nucl. Med. Mol. Imaging 43 2423-32). While conventional dual-head SPECT scanners using parallel-hole collimators and scintillation crystals with photomultiplier tubes are still the workhorse of cardiac SPECT, they have the limitations of low photon sensitivity (~130 count s^-1 MBq^-1), poor image resolution (~15 mm) (Imbert et al 2012 J. Nucl. Med. 53 1897-903), relatively long acquisition time, inefficient use of the detector, high radiation dose, etc. Recently our field observed an exciting growth of new developments of dedicated cardiac scanners and collimators, as well as novel imaging algorithms for quantitative cardiac SPECT. These developments have opened doors to new applications with potential clinical impact, including ultra-low-dose imaging, absolute quantification of myocardial blood flow (MBF) and coronary flow reserve (CFR), multi-radionuclide imaging, and improved image quality as a result of attenuation, scatter, motion, and partial volume corrections (PVCs). In this article, we review the recent advances in cardiac SPECT instrumentation and imaging methods. This review mainly focuses on the most recent developments published since 2012 and points to the future of cardiac SPECT from an imaging physics perspective.

Self-Gating: An Adaptive Center-of-Mass Approach for Respiratory Gating in PET

The goal is to develop an adaptive center-of-mass (COM)-based approach for device-less respiratory gating of list-mode positron emission tomography (PET) data. Our method contains two steps. The first is to automatically extract an optimized respiratory motion signal from the list-mode data during acquisition. The respiratory motion signal was calculated by tracking the location of COM within a volume of interest (VOI). The signal prominence (SP) was calculated based on Fourier analysis of the signal. The VOI was adaptively optimized to maximize SP. The second step is to automatically correct signal-flipping effects. The sign of the signal was determined based on the assumption that the average patient spends more time during expiration than inspiration. To validate our methods, thirty-one ¹⁸F-FDG patient scans were included in this paper. An external device-based signal was used as the gold standard, and the correlation coefficient of the data-driven signal with the device-based signal was measured. Our method successfully extracted respiratory signal from 30 out of 31 datasets. The failure case was due to lack of uptake in the field of view. Moreover, our sign determination method obtained correct results for all scans excluding the failure case. Quantitatively, the proposed signal extraction approach achieved a median correlation of 0.85 with the device-based signal. Gated images using optimized data-driven signal showed improved lesion contrast over static image and were comparable to those using device-based signal. We presented a new data-driven method to automatically extract respiratory motion signal from list-mode PET data by optimizing VOI for COM calculation, as well as determine motion direction from signal asymmetry. Successful application of the proposed method on most clinical datasets and comparison with device-based signal suggests its potential of serving as an alternative to external respiratory monitors.

Correction of hysteretic respiratory motion in SPECT myocardial perfusion imaging: Simulation and patient studies

PURPOSE

Amplitude-based respiratory gating is known to capture the extent of respiratory motion (RM) accurately but results in residual motion in the presence of respiratory hysteresis. In our previous study, we proposed and developed a novel approach to account for respiratory hysteresis by applying the Bouc-Wen (BW) model of hysteresis to external surrogate signals of anterior/posterior motion of the abdomen and chest with respiration. In this work, using simulated and clinical SPECT myocardial perfusion imaging (MPI) studies, we investigate the effects of respiratory hysteresis and evaluate the benefit of correcting it using the proposed BW model in comparison with the abdomen signal typically employed clinically.

METHODS

The MRI navigator data acquired in free-breathing human volunteers were used in the specially modified 4D NCAT phantoms to allow simulating three types of respiratory patterns: monotonic, mild hysteresis, and strong hysteresis with normal myocardial uptake, and perfusion defects in the anterior, lateral, inferior, and septal locations of the mid-ventricular wall. Clinical scans were performed using a Tc-99m sestamibi MPI protocol while recording respiratory signals from thoracic and abdomen regions using a visual tracking system (VTS). The performance of the correction using the respiratory signals was assessed through polar map analysis in phantom and 10 clinical studies selected on the basis of having substantial RM.

RESULTS

In phantom studies, simulations illustrating normal myocardial uptake showed significant differences (P < 0.001) in the uniformity of the polar maps between the RM uncorrected and corrected. No significant differences were seen in the polar map uniformity across the RM corrections. Studies simulating perfusion defects showed significantly decreased errors (P < 0.001) in defect severity and extent for the RM corrected compared to the uncorrected. Only for the strong hysteretic pattern, there was a significant difference (P < 0.001) among the RM corrections. The errors in defect severity and extent for the RM correction using abdomen signal were significantly higher compared to that of the BW (severity = -4.0%, P < 0.001; extent = -65.4%, P < 0.01) and chest (severity = -4.1%, P < 0.001; extent = -52.5%, P < 0.01) signals. In clinical studies, the quantitative analysis of the polar maps demonstrated qualitative and quantitative but not statistically significant differences (P = 0.73) between the correction methods that used the BW signal and the abdominal signal.

CONCLUSIONS

This study shows that hysteresis in respiration affects the extent of residual motion left in the RM-binned data, which can impact wall uniformity and the visualization of defects. Thus, there appears to be the potential for improved accuracy in reconstruction in the presence of hysteretic RM with the BW model method providing a possible step in the direction of improvement.

Data-Driven Methods for the Determination of Anterior-Posterior Motion in PET

Physiological motion combined with elongated scanning times in PET leads to image degradation and quantification errors. Correction approaches usually require 1-D signals that can be obtained with hardware-based or data-driven methods. Most of the latter are optimized or limited to capture internal motion along the superior-inferior (S-I) direction. In this work we present methods for also extracting anterior-posterior (A-P) motion from PET data and propose a set of novel weighting mechanisms that can be used to emphasize certain lines-of-response (LORs) for an increased sensitivity and better signal-to-noise ratio (SNR). The proper functioning of the methods was verified in a phantom experiment. Further, their application to clinical [18^F]-FDG-PET data of 72 patients revealed that using the weighting mechanisms leads to signals with significantly higher spectral respiratory weights, i.e. signals with higher quality. Information about multi-dimensional motion is contained in PET data and can be derived with data-driven methods. Motion models or correction techniques such as respiratory gating might benefit from the proposed methods as they allow to describe the three-dimensional movements of PET-positive structures more precisely.

Improved intrinsic motion detection using time-of-flight PET

Intrinsic or data-driven respiratory and cardiac motion detection often track the center-of-mass (COM) change in PET data to derive a motion gating signal. The effectiveness of this method depends on the contrast of the moving target to the relatively stationary background. The stationary background leads to a reduced COM displacement in PET data. Further, the COM calculated using axially truncated PET data is biased. To improve intrinsic motion detection for motion compensated image reconstruction, we use the time-of-flight (TOF) PET data of the original object f(x) to calculate the non-TOF PET data of a volume-of-interest (VOI) weighted object f(x)w(x) . The VOI-weighting w(x) can be chosen to reduce contribution from the stationary background. The reduced background in f(x)w(x) leads to an observed increase in the COM displacement. We also derive rebinning equations to obtain the exact axial COM using axially truncated PET data. To assess the quality of the motion gating signal, we analyze the variance property of the COM using different methods, including with(out) VOI weighting and with(out) compensation for axial data truncation. Analytical simulations, phantom and patient data demonstrate the effectiveness of the proposed approach in identifying the motion phase and in deriving a gating signal to be used for motion-compensated image reconstruction.

Adaptation of the modified Bouc-Wen model to compensate for hysteresis in respiratory motion for the list-mode binning of cardiac SPECT and PET acquisitions: testing using MRI

PURPOSE

Binning list-mode acquisitions as a function of a surrogate signal related to respiration has been employed to reduce the impact of respiratory motion on image quality in cardiac emission tomography (SPECT and PET). Inherent in amplitude binning is the assumption that there is a monotonic relationship between the amplitude of the surrogate signal and respiratory motion of the heart. This assumption is not valid in the presence of hysteresis when heart motion exhibits a different relationship with the surrogate during inspiration and expiration. The purpose of this study was to investigate the novel approach of using the Bouc-Wen (BW) model to provide a signal accounting for hysteresis when binning list-mode data with the goal of thereby improving motion correction. The study is based on the authors' previous observations that hysteresis between chest and abdomen markers was indicative of hysteresis between abdomen markers and the internal motion of the heart.

METHODS

In 19 healthy volunteers, they determined the internal motion of the heart and diaphragm in the superior-inferior direction during free breathing using MRI navigators. A visual tracking system (vts) synchronized with MRI acquisition tracked the anterior-posterior motions of external markers placed on the chest and abdomen. These data were employed to develop and test the Bouc-Wen model by inputting the vts derived chest and abdomen motions into it and using the resulting output signals as surrogates for cardiac motion. The data of the volunteers were divided into training and testing sets. The training set was used to obtain initial values for the model parameters for all of the volunteers in the set, and for set members based on whether they were or were not classified as exhibiting hysteresis using a metric derived from the markers. These initial parameters were then employed with the testing set to estimate output signals. Pearson's linear correlation coefficient between the abdomen, chest, average of chest and abdomen markers, and Bouc-Wen derived signals versus the true internal motion of the heart from MRI was used to judge the signals match to the heart motion.

RESULTS

The results show that the Bouc-Wen model generated signals demonstrated strong correlation with the heart motion. This correlation was slightly larger on average than that of the external surrogate signals derived from the abdomen marker, and average of the abdomen and chest markers, but was not statistically significantly different from them.

CONCLUSIONS

The results suggest that the proposed model has the potential to be a unified framework for modeling hysteresis in respiratory motion in cardiac perfusion studies and beyond.

MRI Investigation of the Linkage Between Respiratory Motion of the Heart and Markers on Patient's Abdomen and Chest: Implications for Respiratory Amplitude Binning List-Mode PET and SPECT Studies

Respiratory motion of the heart impacts the diagnostic accuracy of myocardial-perfusion emission-imaging studies. Amplitude binning has come to be the method of choice for binning list-mode based acquisitions for correction of respiratory motion in PET and SPECT. In some subjects respiratory motion exhibits hysteretic behavior similar to damped non-linear cyclic systems. The detection and correction of hysteresis between the signals from surface movement of the patient's body used in binning and the motion of the heart within the chest remains an open area for investigation. This study reports our investigation in nine volunteers of the combined MRI tracking of the internal respiratory motion of the heart using Navigators with stereo-tracking of markers on the volunteer's chest and abdomen by a visual-tracking system (VTS). The respiratory motion signals from the internal organs and the external markers were evaluated for hysteretic behavior analyzing the temporal correspondence of the signals. In general, a strong, positive correlation between the external marker motion (AP direction) and the internal heart motion (SI direction) during respiration was observed. The average ± standard deviation in the Spearman's ranked correlation coefficient (ρ) over the nine volunteer studied was 0.92 ± 0.1 between the external abdomen marker and the internal heart, and 0.87 ± 0.2 between the external chest marker and the internal heart. However despite the good correlation on average for the nine volunteers, in three studies a poor correlation was observed due to hysteretic behavior between inspiration and expiration for either the chest marker and the internal motion of the heart, or the abdominal marker and the motion of the heart. In all cases we observed a good correlation of at least either the abdomen or the chest with the heart. Based on this result, we propose the use of marker motion from both the chest and abdomen regions when estimating the internal heart motion to detect and address hysteresis when binning list-mode emission data.

Use of MRI to assess the prediction of heart motion with gross body motion in myocardial perfusion imaging by stereotracking of markers on the body surface

PURPOSE

The aim of this study is to determine using MRI in volunteers whether the rigid-body-motion (RBM) model can be approximately used to estimate the gross body-motion of the heart from that of external markers on patient's chest. Our target clinical application is to use a visual-tracking-system (VTS) which employs stereoimaging to estimate heart motion during SPECT/CT and PET∕CT myocardial perfusion imaging.

METHODS

To investigate body-motion separate from the respiration the authors had the volunteers hold their breath during the acquisition of a sequence of two sets of EKG-triggered MRI sagittal slices. The first set was acquired pre-motion, and the second postmotion. The motion of the heart within each breath-hold set of slices was estimated by registration to the semiautomatic 3D segmentation of the heart region in a baseline set acquired using the Navigator technique. The motion of the heart between the pre- and postmotion sets was then determined as the difference in the individual motions in comparison to the Navigator sets. An analysis of the combined motion of the individual markers on the chest was used to obtain an estimate of the six-degree-of-freedom RBM from the VTS system. The metric for judging agreement between the motion estimated by MRI and the VTS was the average error. This was defined as the average of the magnitudes of the differences in the vector displacements of all voxels in the heart region. Studies with the Data Spectrum Anthropomorphic Phantom and "No-Motion" studies in which the volunteer did not intentionally move were used to establish a baseline for agreement. With volunteer studies a t-test was employed to determine when statistically significant differences in Average Errors occurred compared to the No-motion studies.

RESULTS

For phantom acquisitions, the Average Error when the motion was just translation was 0.1 mm. With complex motions, which included a combination of rotations and translations, the Average Error increased to 3.6 mm. In the volunteers the Average Error averaged over all No-Motion acquisitions was 1.0 mm. For the case of translational motion, which might be expected to be RBM, the Average Error averaged over all volunteer studies increased to 2.6 mm, which was statistically different from the No-Motion studies. For the case of bends and twists of the torso, which would be expected to challenge the RBM model, the Average Error averaged over all such volunteer studies was 4.9 mm and was again statistically different. Investigations of motion of the arm including just bending at the elbow and leg motion resulted in Average Errors which were not statistically different from the No-Motion studies. However, when shoulder movement was included with arm motion the Average Error was near that of torso bends and twists, and statistically different.

CONCLUSIONS

Use of the RBM model with VTS predictions of heart motion during reconstruction should decrease the extent of artifacts for the types of patient motion studied. The impact of correction would be less for torso bends and twists, and arm motion which includes the shoulders.