Reduction of breathing irregularity-related motion artifacts in low-pitch spiral 4D CT by optimized projection binning.
Ontology highlight
ABSTRACT: Background
Respiration-correlated CT (4D CT) is the basis of radiotherapy treatment planning of thoracic and abdominal tumors. Current clinical 4D CT images suffer, however, from artifacts due to unfulfilled assumptions concerning breathing pattern regularity. We propose and evaluate modifications to existing low-pitch spiral 4D CT reconstruction protocols to counteract respective artifacts.Methods
The proposed advanced reconstruction (AR) approach consists of two steps that build on each other: (1) statistical analysis of the breathing signal recorded during CT data acquisition and extraction of a patient-specific reference breathing cycle for projection binning; (2) incorporation of an artifact measure into the reconstruction. 4D CT data of 30 patients were reconstructed by standard phase- and local amplitude-based reconstruction (PB, LAB) and compared with images obtained by AR. The number of artifacts was evaluated and artifact statistics correlated to breathing curve characteristics.Results
AR reduced the number of 4D CT artifacts by 31% and 27% compared to PB and LAB; the reduction was most pronounced for irregular breathing curves.Conclusions
We described a two-step optimization of low-pitch spiral 4D CT reconstruction to reduce artifacts in the presence of breathing irregularity and illustrated that the modifications to existing reconstruction solutions are effective in terms of artifact reduction.
SUBMITTER: Werner R
PROVIDER: S-EPMC5477247 | biostudies-literature | 2017 Jun
REPOSITORIES: biostudies-literature
Publications
Reduction of breathing irregularity-related motion artifacts in low-pitch spiral 4D CT by optimized projection binning.
Werner René R Hofmann Christian C Mücke Eike E Gauer Tobias T
Radiation oncology (London, England) 20170619 1
<h4>Background</h4>Respiration-correlated CT (4D CT) is the basis of radiotherapy treatment planning of thoracic and abdominal tumors. Current clinical 4D CT images suffer, however, from artifacts due to unfulfilled assumptions concerning breathing pattern regularity. We propose and evaluate modifications to existing low-pitch spiral 4D CT reconstruction protocols to counteract respective artifacts.<h4>Methods</h4>The proposed advanced reconstruction (AR) approach consists of two steps that buil ...[more]
Similar Datasets
Project description:Respiration-correlated CT (RCCT) images produced with commonly used phase-based sorting of CT slices often exhibit discontinuity artifacts between CT slices, caused by cycle-to-cycle amplitude variations in respiration. Sorting based on the displacement of the respiratory signal yields slices at more consistent respiratory motion states and hence reduces artifacts, but missing image data (gaps) may occur. The authors report on the application of a respiratory motion model to produce an RCCT image set with reduced artifacts and without missing data.Input data consist of CT slices from a cine CT scan acquired while recording respiration by monitoring abdominal displacement. The model-based generation of RCCT images consists of four processing steps: (1) displacement-based sorting of CT slices to form volume images at 10 motion states over the cycle; (2) selection of a reference image without gaps and deformable registration between the reference image and each of the remaining images; (3) generation of the motion model by applying a principal component analysis to establish a relationship between displacement field and respiration signal at each motion state; (4) application of the motion model to deform the reference image into images at the 9 other motion states. Deformable image registration uses a modified fast free-form algorithm that excludes zero-intensity voxels, caused by missing data, from the image similarity term in the minimization function. In each iteration of the minimization, the displacement field in the gap regions is linearly interpolated from nearest neighbor nonzero intensity slices. Evaluation of the model-based RCCT examines three types of image sets: cine scans of a physical phantom programmed to move according to a patient respiratory signal, NURBS-based cardiac torso (NCAT) software phantom, and patient thoracic scans.Comparison in physical motion phantom shows that object distortion caused by variable motion amplitude in phase-based sorting is visibly reduced with model-based RCCT. Comparison of model-based RCCT to original NCAT images as ground truth shows best agreement at motion states whose displacement-sorted images have no missing slices, with mean and maximum discrepancies in lung of 1 and 3 mm, respectively. Larger discrepancies correlate with motion states having a larger number of missing slices in the displacement-sorted images. Artifacts in patient images at different motion states are also reduced. Comparison with displacement-sorted patient images as a ground truth shows that the model-based images closely reproduce the ground truth geometry at different motion states.Results in phantom and patient images indicate that the proposed method can produce RCCT image sets with reduced artifacts relative to phase-sorted images, without the gaps inherent in displacement-sorted images. The method requires a reference image at one motion state that has no missing data. Highly irregular breathing patterns can affect the method's performance, by introducing artifacts in the reference image (although reduced relative to phase-sorted images), or in decreased accuracy in the image prediction of motion states containing large regions of missing data.
| S-EPMC3360690 | biostudies-other
Project description:We have previously developed a retrospective 4D-MRI technique using body area as the respiratory surrogate, but generally, the reconstructed 4D MR images suffer from severe or mild artifacts mainly caused by irregular motion during image acquisition. Those image artifacts may potentially affect the accuracy of tumor target delineation or the shape representation of surrounding nontarget tissues and organs. So the purpose of this study is to propose an approach employing principal component analysis (PCA), combined with a linear polynomial fitting model, to remodel the displacement vector fields (DVFs) obtained from deformable image registration (DIR), with the main goal of reducing the motion artifacts in 4D MR images. Seven patients with hepatocellular carcinoma (2/7) or liver metastases (5/7) in the liver, as well as a patient with non-small cell lung cancer (NSCLC), were enrolled in an IRB-approved prospective study. Both CT and MR simulations were performed for each patient for treatment planning. Multiple-slice, multiple-phase, cine-MRI images were acquired in the axial plane for 4D-MRI reconstruction. Single-slice 2D cine-MR images were acquired across the center of the tumor in axial, coronal, and sagittal planes. For a 4D MR image dataset, the DVFs in three orthogonal direction (inferior–superior (SI), anterior–posterior (AP), and medial–lateral (ML)) relative to a specific reference phase were calculated using an in-house DIR algorithm. The DVFs were preprocessed in three temporal and spatial dimensions using a polynomial fitting model, with the goal of correcting the potential registration errors introduced by three-dimensional DIR. Then PCA was used to decompose each fitted DVF into a linear combination of three principal motion bases whose spanned subspaces combined with their projections had been validated to be sufficient to represent the regular respiratory motion. By wrapping the reference MR image using the remodeled DVFs, 'synthetic' MR images with reduced motion artifacts were generated at selected phase. Tumor motion trajectories derived from cine-MRI, 4D CT, original 4D MRI, and 'synthetic' 4D MRI were analyzed in the SI, AP, and ML directions, respectively. Their correlation coefficient (CC) and difference (D) in motion amplitude were calculated for comparison. Of all the patients, the means and standard deviations (SDs) of CC comparing 'synthetic' 4D MRI and cine-MRI were 0.98 ± 0.01, 0.98 ± 0.01, and 0.99 ± 0.01 in SI, AP, and ML directions, respectively. The mean ± SD Ds were 0.59 ± 0.09 mm, 0.29± 0.10 mm, and 0.15 ± 0.05 mm in SI, AP and ML directions, respectively. The means and SDs of CC comparing 'synthetic' 4D MRI and 4D CT were 0.96 ± 0.01, 0.95± 0.01, and 0.95 ± 0.01 in SI, AP, and ML directions, respectively. The mean ± SD Ds were 0.76 ± 0.20 mm, 0.33 ± 0.14 mm, and 0.19± 0.07 mm in SI, AP, and ML directions, respectively. The means and SDs of CC comparing 'synthetic' 4D MRI and original 4D MRI were 0.98 ± 0.01, 0.98± 0.01, and 0.97± 0.01 in SI, AP, and ML directions, respectively. The mean ± SD Ds were 0.58 ± 0.10 mm, 0.30 ± 0.09mm, and 0.17 ± 0.04 mm in SI, AP, and ML directions, respectively. In this study we have proposed an approach employing PCA combined with a linear polynomial fitting model to capture the regular respiratory motion from a 4D MR image dataset. And its potential usefulness in reducing motion artifacts and improving image quality has been demonstrated by the preliminary results in oncological patients.
| S-EPMC5690092 | biostudies-literature
Project description:Liver dynamic contrast-enhanced MRI (DCE-MRI) requires high spatiotemporal resolution and large field of view to clearly visualize all relevant enhancement phases and detect early-stage liver lesions. The low-rank plus sparse (L + S) reconstruction outperforms standard sparsity-only-based reconstruction through separation of low-rank background component (L) and sparse dynamic components (S). However, the L + S decomposition is sensitive to respiratory motion so that image quality is compromised when breathing occurs during long time data acquisition. To enable high quality reconstruction for free-breathing liver 4D DCE-MRI, this paper presents a novel method called SMC-LS, which incorporates Sliding Motion Compensation into the standard L + S reconstruction. The global superior-inferior displacement of the internal abdominal organs is inferred directly from the undersampled raw data and then used to correct the breathing induced sliding motion which is the dominant component of respiratory motion. With sliding motion compensation, the reconstructed temporal frames are roughly registered before applying the standard L + S decomposition. The proposed method has been validated using free-breathing liver 4D MRI phantom data, free-breathing liver 4D DCE-MRI phantom data, and in vivo free breathing liver 4D MRI dataset. Results demonstrated that SMC-LS reconstruction can effectively reduce motion blurring artefacts and preserve both spatial structures and temporal variations at a sub-second temporal frame rate for free-breathing whole-liver 4D DCE-MRI.
| S-EPMC6411681 | biostudies-literature
Project description:PurposeTo develop a continuous-acquisition cardiac self-gated spiral pulse sequence and a respiratory motion-compensated reconstruction strategy for free-breathing cine imaging.MethodsCine data were acquired continuously on a 3T scanner for 8 seconds per slice without ECG gating or breath-holding, using a golden-angle gradient echo spiral pulse sequence. Cardiac motion information was extracted by applying principal component analysis on the gridded 8 × 8 k-space center data. Respiratory motion was corrected by rigid registration on each heartbeat. Images were reconstructed using a low-rank and sparse (L+S) technique. This strategy was evaluated in 37 healthy subjects and 8 subjects undergoing clinical cardiac MR studies. Image quality was scored (1-5 scale) in a blinded fashion by 2 experienced cardiologists. In 13 subjects with whole-heart coverage, left ventricular ejection fraction (LVEF) from SPiral Acquisition with Respiratory correction and Cardiac Self-gating (SPARCS) was compared to that from a standard ECG-gated breath-hold balanced steady-state free precession (bSSFP) cine sequence.ResultsThe self-gated signal was successfully extracted in all cases and demonstrated close agreement with the acquired ECG signal (mean bias, -0.22 ms). The mean image score across all subjects was 4.0 for reconstruction using the L+S model. There was good agreement between the LVEF derived from SPARCS and the gold-standard bSSFP technique.ConclusionSPARCS successfully images cardiac function without the need for ECG gating or breath-holding. With an 8-second data acquisition per slice, whole-heart cine images with clinically acceptable spatial and temporal resolution and image quality can be acquired in <90 seconds of free-breathing acquisition.
| S-EPMC6510595 | biostudies-literature
Project description:Head motion represents one of the greatest technical obstacles in magnetic resonance imaging (MRI) of the human brain. Accurate detection of artifacts induced by head motion requires precise estimation of movement. However, head motion estimates may be corrupted by artifacts due to magnetic main field fluctuations generated by body motion. In the current report, we examine head motion estimation in multiband resting state functional connectivity MRI (rs-fcMRI) data from the Adolescent Brain and Cognitive Development (ABCD) Study and comparison 'single-shot' datasets. We show that respirations contaminate movement estimates in functional MRI and that respiration generates apparent head motion not associated with functional MRI quality reductions. We have developed a novel approach using a band-stop filter that accurately removes these respiratory effects from motion estimates. Subsequently, we demonstrate that utilizing a band-stop filter improves post-processing fMRI data quality. Lastly, we demonstrate the real-time implementation of motion estimate filtering in our FIRMM (Framewise Integrated Real-Time MRI Monitoring) software package.
| S-EPMC7307712 | biostudies-literature
Project description:We introduce a kernel low-rank algorithm to recover free-breathing and ungated dynamic MRI from spiral acquisitions without explicit k-space navigators. It is often challenging for low-rank methods to recover free-breathing and ungated images from undersampled measurements; extensive cardiac and respiratory motion often results in the Casorati matrix not being sufficiently low-rank. Therefore, we exploit the non-linear structure of the dynamic data, which gives the low-rank kernel matrix. Unlike prior work that rely on navigators to estimate the manifold structure, we propose a kernel low-rank matrix completion method to directly fill in the missing k-space data from variable density spiral acquisitions. We validate the proposed scheme using simulated data and in-vivo data. Our results show that the proposed scheme provides improved reconstructions compared to the classical methods such as low-rank and XD-GRASP. The comparison with breath-held cine data shows that the quantitative metrics agree, whereas the image quality is marginally lower.
| S-EPMC7806246 | biostudies-literature
Project description:Nuclear magnetic resonance (NMR) is a mainstay of biophysical studies that provides atomic level readouts to formulate molecular mechanisms. Side chains are particularly important to derive mechanisms involving proteins as they carry functional groups, but NMR studies of side chains are often limited by challenges in assigning their signals. Here, we designed a novel computational method that combines spectral derivatives and matrix square-rooting to produce reliable 4D covariance maps from routinely acquired 3D spectra and facilitates side chain resonance assignments. Thus, we generate two 4D maps from 3D-HcccoNH and 3D-HCcH-TOCSY spectra that each help overcome signal overlap or sensitivity losses. These 4D maps feature HC-HSQCs of individual side chains that can be paired to assigned backbone amide resonances of individual aliphatic signals, and both are obtained from a single modified covariance calculation. Further, we present 4D maps produced using conventional triple resonance experiments to easily assign asparagine side chain amide resonances. The 4D covariance maps encapsulate the lengthy manual pattern recognition used in traditional assignment methods and distill the information as correlations that can be easily visualized. We showcase the utility of the 4D covariance maps with a 10 kDa peptidyl carrier protein and a 52 kDa cyclization domain from a nonribosomal peptide synthetase.
| S-EPMC8480538 | biostudies-literature
Project description:Recent mobile brain/body imaging (MoBI) techniques based on active electrode scalp electroencephalogram (EEG) allow the acquisition and real-time analysis of brain dynamics during active unrestrained motor behavior involving whole body movements such as treadmill walking, over-ground walking and other locomotive and non-locomotive tasks. Unfortunately, MoBI protocols are prone to physiological and non-physiological artifacts, including motion artifacts that may contaminate the EEG recordings. A few attempts have been made to quantify these artifacts during locomotion tasks but with inconclusive results due in part to methodological pitfalls. In this paper, we investigate the potential contributions of motion artifacts in scalp EEG during treadmill walking at three different speeds (1.5, 3.0, and 4.5 km/h) using a wireless 64 channel active EEG system and a wireless inertial sensor attached to the subject's head. The experimental setup was designed according to good measurement practices using state-of-the-art commercially available instruments, and the measurements were analyzed using Fourier analysis and wavelet coherence approaches. Contrary to prior claims, the subjects' motion did not significantly affect their EEG during treadmill walking although precaution should be taken when gait speeds approach 4.5 km/h. Overall, these findings suggest how MoBI methods may be safely deployed in neural, cognitive, and rehabilitation engineering applications.
| S-EPMC4710850 | biostudies-literature
Project description:PURPOSE:The utility of four-dimensional (4D) spiral flow in imaging of stenotic flows in both phantoms and human subjects with aortic stenosis is investigated. METHODS:The method performs 4D flow acquisitions through a stack of interleaved spiral k-space readouts. Relative to conventional 4D flow, which performs Cartesian readout, the method has reduced echo time. Thus, reduced flow artifacts are observed when imaging high-speed stenotic flows. Four-dimensional spiral flow also provides significant savings in scan times relative to conventional 4D flow. RESULTS:In vitro experiments were performed under both steady and pulsatile flows in a phantom model of severe stenosis (one inch diameter at the inlet, with 87% area reduction at the throat of the stenosis) while imaging a 6-cm axial extent of the phantom, which included the Gaussian-shaped stenotic narrowing. In all cases, gradient strength and slew rate for standard clinical acquisitions, and identical field of view and resolution were used. For low steady flow rates, quantitative and qualitative results showed a similar level of accuracy between 4D spiral flow (echo time [TE] = 2 ms, scan time = 40 s) and conventional 4D flow (TE = 3.6 ms, scan time = 1:01 min). However, in the case of high steady flow rates, 4D spiral flow (TE = 1.57 ms, scan time = 38 s) showed better visualization and accuracy as compared to conventional 4D flow (TE = 3.2 ms, scan time = 51 s). At low pulsatile flow rates, a good agreement was observed between 4D spiral flow (TE = 2 ms, scan time = 10:26 min) and conventional 4D flow (TE = 3.6 ms, scan time = 14:20 min). However, in the case of high flow-rate pulsatile flows, 4D spiral flow (TE = 1.57 ms, scan time = 10:26 min) demonstrated better visualization as compared to conventional 4D flow (TE = 3.2 ms, scan time = 14:20 min). The feasibility of 4D spiral flow was also investigated in five normal volunteers and four subjects with mild-to-moderate aortic stenosis. The approach achieved TE = 1.68 ms and scan time = 3:44 min. The conventional sequence achieved TE = 2.9 ms and scan time = 5:23 min. In subjects with aortic stenosis, we also compared both MRI methods with Doppler ultrasound (US) in the measurement of peak velocity, time to peak systolic velocity, and eject time. Bland-Altman analysis revealed that, when comparing peak velocities, the discrepancy between Doppler US and 4D spiral flow was significantly less than the discrepancy between Doppler and 4D Cartesian flow (2.75 cm/s vs. 10.25 cm/s), whereas the two MR methods were comparable (-5.75 s vs. -6 s) for time to peak. However, for the estimation of eject time, relative to Doppler US, the discrepancy for 4D conventional flow was smaller than that of 4D spiral flow (-16.25 s vs. -20 s). CONCLUSION:Relative to conventional 4D flow, 4D spiral flow achieves substantial reductions in both the TE and scan times; therefore, utility for it should be sought in a variety of in vivo and complex flow imaging applications.
| S-EPMC7264044 | biostudies-literature
Project description:The purpose of this study was to develop a methodology to quantitatively measure the thorax-to-abdomen breathing ratio from a 4DCT dataset for breathing motion modeling and breathing motion studies.The thorax-to-abdomen breathing ratio was quantified by measuring the rate of cross-sectional volume increase throughout the thorax and abdomen as a function of tidal volume. Twenty-six 16-slice 4DCT patient datasets were acquired during quiet respiration using a protocol that acquired 25 ciné scans at each couch position. Fifteen datasets included data from the neck through the pelvis. Tidal volume, measured using a spirometer and abdominal pneumatic bellows, was used as breathing-cycle surrogates. The cross-sectional volume encompassed by the skin contour when compared for each CT slice against the tidal volume exhibited a nearly linear relationship. A robust iteratively reweighted least squares regression analysis was used to determine η(i), defined as the amount of cross-sectional volume expansion at each slice i per unit tidal volume. The sum Ση(i) throughout all slices was predicted to be the ratio of the geometric expansion of the lung and the tidal volume; 1.11. The Xiphoid process was selected as the boundary between the thorax and abdomen. The Xiphoid process slice was identified in a scan acquired at mid-inhalation. The imaging protocol had not originally been designed for purposes of measuring the thorax-to-abdomen breathing ratio so the scans did not extend to the anatomy with η(i) = 0. Extrapolation of η(i)-η(i) = 0 was used to include the entire breathing volume. The thorax and abdomen regions were individually analyzed to determine the thorax-to-abdomen breathing ratios. There were 11 image datasets that had been scanned only through the thorax. For these cases, the abdomen breathing component was equal to 1.11 - Ση(i) where the sum was taken throughout the thorax.The average Ση(i) for thorax and abdomen image datasets was found to be 1.20 ± 0.17, close to the expected value of 1.11. The thorax-to-abdomen breathing ratio was 0.32 ± 0.24. The average Ση(i) was 0.26 ± 0.14 in the thorax and 0.93 ± 0.22 in the abdomen. In the scan datasets that encompassed only the thorax, the average Ση(i) was 0.21 ± 0.11.A method to quantify the relationship between abdomen and thoracic breathing was developed and characterized.
| S-EPMC3669136 | biostudies-other