Iterative sorting for four-dimensional CT images based on internal anatomy motion

Data-driven gated positron emission tomography/computed tomography for radiotherapy

Purpose

Software-based data-driven gated (DDG) positron emission tomography/computed tomography (PET/CT) has replaced hardware-based 4D PET/CT. The purpose of this article was to review DDG PET/CT, which could improve the accuracy of treatment response assessment, tumor motion evaluation, and target tumor contouring with whole-body (WB) PET/CT for radiotherapy (RT).

Material and methods

This review covered the topics of 4D PET/CT with hardware gating, advancements in PET instrumentation, DDG PET, DDG CT, and DDG PET/CT based on a systematic literature review. It included a discussion of the large axial field-of-view (AFOV) PET detector and a review of the clinical results of DDG PET and DDG PET/CT.

Results

DDG PET matched or outperformed 4D PET with hardware gating. DDG CT was more compatible with DDG PET than 4D CT, which required hardware gating. DDG CT could replace 4D CT for RT. DDG PET and DDG CT for DDG PET/CT can be incorporated in a WB PET/CT of less than 15 min scan time on a PET/CT scanner of at least 25 cm AFOV PET detector.

Conclusions

DDG PET/CT could correct the misregistration and tumor motion artifacts in a WB PET/CT and provide the quantitative PET and tumor motion information of a registered PET/CT for RT.

Data-driven gated (DDG) CT: An automated respiratory gating method to enable DDG PET/CT

BACKGROUND

The accuracy of PET quantification and localization can be compromised if a misregistered CT is used for attenuation correction (AC) in PET/CT. As data-driven gating (DDG) continues to grow in clinical use, these issues are becoming more relevant with respect to solutions for gated CT.

PURPOSE

In this work, a new automated data-driven gated (DDG) CT method was developed to provide average CT and DDG CT for AC of PET and DDG PET, respectively.

METHODS

An automatic DDG CT was developed to provide the end-expiratory (EE) and end-inspiratory (EI) phases of images from low-dose cine CT images, with all phases being averaged to generate an average CT. The respiratory phases of EE and EI were determined according to lung region Hounsfield unit (HU) values and body outline contours. The average CT was used for AC of baseline PET and DDG CT at EE phase was used for AC of DDG PET at the quiescent or EE phase. The EI and EE phases obtained with DDG CT were used for assessing the magnitude of respiratory motion. The proposed DDG CT was compared to two commercial CT gating methods: 1) 4D CT (external device based) and 2) D4D CT (DDG based) in 38 patient data sets with respect to respiratory phase image selection, lung HU, lung volume, and image artifacts. In a separate set of twenty consecutive PET/CT studies containing a mix of ¹⁸ F-FDG, ⁶⁸ Ga-Dotatate, and ⁶⁴ Cu-Dotatate scans, the proposed DDG CT was compared with D4D CT for impacts on registration and quantification in DDG PET/CT.

RESULTS

In the EE phase, the images selected by DDG CT and 4D CT were identical 62.5±21.6% of the time, while DDG CT and D4D CT were 6.5±9.7%, and 4D CT and D4D CT were 8.6±12.2%. These differences in EE phase image selection were significant (p<0.0001). In the EI phase, the images selected by DDG CT and 4D CT were identical 68.2±18.9% of the time, DDG CT and D4D CT were 63.9±18.8%, and 4D CT and D4D CT were 61.2±19.8%. These differences were not significant. The mean lung HU and volumes were not statistically different (p > 0.1) among the three methods. In some studies, DDG CT was better than D4D or 4D CT in appropriate selection of the EE and EI phases, and D4D CT was found to reverse the EE and EI phases or not select the correct images by visual inspection. A statistically significant improvement of DDG CT over D4D CT for AC of DDG PET was also demonstrated with PET quantification analysis. When irregular breath cycles were present in the cine CT, DDG CT could be used to replace average CT for improved AC of baseline PET.

CONCLUSION

A new automatic DDG CT was developed to tackle the issues of misregistration and tumor motion in PET/CT imaging. DDG CT was significantly more consistent than D4D CT in selecting the EE phase images as the clinical standard of 4D CT. When compared to both commercial gated CT methods of 4D CT and D4D CT, DDG CT appeared to be more robust in the lower lung and upper diaphragm regions where misregistration and tumor motion often occur. DDG CT offered improved AC for DDG PET relative to D4D CT. In cases with irregular respiratory motion, DDG CT improved AC over average CT for baseline PET. The new DDG CT provides the benefits of 4D CT without the need for external device gating. This article is protected by copyright. All rights reserved.

Data-driven respiratory gating based on localized diaphragm sensing in TOF PET

It is important to measure the respiratory cycle in positron emission tomography (PET) to enhance the contrast of the tumor as well as the accuracy of its localization in organs such as the lung and liver. Several types of data-driven respiratory gating methods, such as center of mass and principal component analysis, have been developed to directly measure the breathing cycle from PET images and listmode data. However, the breathing cycle is still hard to detect in low signal-to-noise ratio (SNR) data, particularly in low dose PET/CT scans. To address this issue, a time-of-flight (TOF) PET is currently utilized for the data-driven respiratory gating because of its higher SNR and better localization of the region of interest. To further improve the accuracy of respiratory gating with TOF information, we propose an accurate data-driven respiratory gating method, which retrospectively derives the respiratory signal using a localized sensing method based on a diaphragm mask in TOF PET data. To assess the accuracy of the proposed method, the performance is evaluated with three patient datasets, and a pressure-belt signal as the ground truth is compared. In our experiments, we validate that the respiratory signal using the proposed data-driven gating method is well matched to the pressure-belt respiratory signal with less than 5% peak time errors and over 80% trace correlations. Based on gated signals, the respiratory-gated image of the proposed method provides more clear edges of organs compared to images using conventional non-TOF methods. Therefore, we demonstrate that the proposed method can achieve improvements for the accuracy of gating signals and image quality.

Unsupervised respiratory signal extraction from ungated cardiac magnetic resonance imaging at rest and during exercise

We propose and evaluate a method to estimate a respiratory signal from ungated cardiac magnetic resonance (CMR) images. Ungated CMR images were acquired in five subjects who performed exercise at different intensity levels under different physiological conditions while breathing freely. The respiratory motion was estimated by applying principal components analysis (PCA). A sign correction procedure was developed to correctly define inspiration and expiration, based on either tracking of the diaphragmatic motion or estimation of the lung volume or a combination of both. Evaluation was done using a plethysmograph signal as reference. There was a good correspondence between the plethysmograph and the estimated respiratory signals. Respiratory motion was effectively captured by one of the PCA components in 88% of the cases. Moreover, the proposed method successfully estimated the respiratory phase in 91% of the evaluated slices. The pipeline is robust, admitting a slight decline in performance with increased exercise intensity. Respiratory motion was accurately estimated by means of PCA and the application of a sign correction procedure. Our method showed promising results even for acquisitions during exercise where excessive body motion occurs. The proposed method provides a way to extract the respiratory signal from ungated CMR images, at rest as well as during exercise, in a fully unsupervised fashion, which may reduce the clinician's workload drastically.

Clinical implementations of 4D pencil beam scanned particle therapy: Report on the 4D treatment planning workshop 2016 and 2017

In 2016 and 2017, the 8th and 9th 4D treatment planning workshop took place in Groningen (the Netherlands) and Vienna (Austria), respectively. This annual workshop brings together international experts to discuss research, advances in clinical implementation as well as problems and challenges in 4D treatment planning, mainly in spot scanned proton therapy. In the last two years several aspects like treatment planning, beam delivery, Monte Carlo simulations, motion modeling and monitoring, QA phantoms as well as 4D imaging were thoroughly discussed. This report provides an overview of discussed topics, recent findings and literature review from the last two years. Its main focus is to highlight translation of 4D research into clinical practice and to discuss remaining challenges and pitfalls that still need to be addressed and to be overcome.

Retrospective data-driven respiratory gating for PET using TOF information

Traditional data-driven respiratory gating method is capable of detecting breathing cycles directly from positron emission tomography (PET) data, but usually fails at low SNR, particularly at low dose PET/CT study. Time-of-flight (TOF) PET has the potential to improve the SNR. In order for TOF information to reduce the statistical noise and boost the performance of respiratory gating, we present a robust data-driven respiratory gating method using TOF information, which retrospectively derived the respiratory signal from the acquired TOF-PET data. The PET data was acquired in list mode format and analyzed in sinogram space. The method was demonstrated with patient datasets acquired on a TOF PET/CT system. Data-driven gating methods by center of mass (COM) and principle component analysis (PCA) algorithm were successfully performed on nonTOF PET and TOF PET dataset. To assess the accuracy of the data-driven respiratory signal, a hardware-based signal was acquired for comparison. The study showed that retrospectively respiratory gating using TOF sinograms has improved the SNR, and outperforms the non-TOF gating under both COM and PCA algorithms.

Target volume and artifact evaluation of a new data-driven 4D CT

PURPOSE

Four-dimensional computed tomography (4D CT) is often used to define the internal gross target volume (IGTV) for radiation therapy of lung cancer. Traditionally, this technique requires the use of an external motion surrogate; however, a new image, data-driven 4D CT, has become available. This study aims to describe this data-driven 4D CT and compare target contours created with it to those created using standard 4D CT.

METHODS AND MATERIALS

Cine CT data of 35 patients undergoing stereotactic body radiation therapy were collected and sorted into phases using standard and data-driven 4D CT. IGTV contours were drawn using a semiautomated method on maximum intensity projection images of both 4D CT methods. Errors resulting from reproducibility of the method were characterized. A comparison of phase image artifacts was made using a normalized cross-correlation method that assigned a score from +1 (data-driven "better") to -1 (standard "better").

RESULTS

The volume difference between the data-driven and standard IGTVs was not significant (data driven was 2.1 ± 1.0% smaller, P = .08). The Dice similarity coefficient showed good similarity between the contours (0.949 ± 0.006). The mean surface separation was 0.4 ± 0.1 mm and the Hausdorff distance was 3.1 ± 0.4 mm. An average artifact score of +0.37 indicated that the data-driven method had significantly fewer and/or less severe artifacts than the standard method (P = 1.5 × 10^-5 for difference from 0).

CONCLUSIONS

On average, the difference between IGTVs derived from data-driven and standard 4D CT was not clinically relevant or statistically significant, suggesting data-driven 4D CT can be used in place of standard 4D CT without adjustments to IGTVs. The relatively large differences in some patients were usually attributed to limitations in automatic contouring or differences in artifacts. Artifact reduction and setup simplicity suggest a clinical advantage to data-driven 4D CT.

Respiratory motion correction of PET using MR-constrained PET-PET registration

Background

Respiratory motion in positron emission tomography (PET) is an unavoidable source of error in the measurement of tracer uptake, lesion position and lesion size. The introduction of PET-MR dual modality scanners opens a new avenue for addressing this issue. Motion models offer a way to estimate motion using a reduced number of parameters. This can be beneficial for estimating motion from PET, which can otherwise be difficult due to the high level of noise of the data.

Method

We propose a novel technique that makes use of a respiratory motion model, formed from initial MR scan data. The motion model is used to constrain PET-PET registrations between a reference PET gate and the gates to be corrected. For evaluation, PET with added FDG-avid lesions was simulated from real, segmented, ultrashort echo time MR data obtained from four volunteers. Respiratory motion was included in the simulations using motion fields derived from real dynamic 3D MR volumes obtained from the same volunteers.

Results

Performance was compared to an MR-derived motion model driven method (which requires constant use of the MR scanner) and to unconstrained PET-PET registration of the PET gates. Without motion correction, a median drop in uncorrected lesion \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {SUV}}_{\mathrm {peak}}$$\end{document}SUVpeak intensity to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$78.4 \pm 18.6 \,\,\%$$\end{document}78.4±18.6% and an increase in median head-foot lesion width, specified by a minimum bounding box, to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$179 \pm 63.7\,\, \%$$\end{document}179±63.7% was observed relative to the corresponding measures in motion-free simulations. The proposed method corrected these values to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$86.9 \pm 13.6\,\, \%$$\end{document}86.9±13.6% (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p<0.001$$\end{document}p<0.001) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100 \pm 29.12\,\, \%$$\end{document}100±29.12% (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p<0.001$$\end{document}p<0.001) respectively, with notably improved performance close to the diaphragm and in the liver. Median lesion displacement across all lesions was observed to be \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6.6 \pm 5.4\,\mathrm {mm}$$\end{document}6.6±5.4mm without motion correction, which was reduced to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3.5 \pm 1.8\,\mathrm {mm}$$\end{document}3.5±1.8mm (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p<0.001$$\end{document}p<0.001) with motion correction.

Discussion

This paper presents a novel technique for respiratory motion correction of PET data in PET-MR imaging. After an initial 30 second MR scan, the proposed technique does not require use of the MR scanner for motion correction purposes, making it suitable for MR-intensive studies or sequential PET-MR. The accuracy of the proposed technique was similar to both comparative methods, but robustness was improved compared to the PET-PET technique, particularly in regions with higher noise such as the liver.

Comparing approaches to correct for respiratory motion in NH3 PET-CT cardiac perfusion imaging

Aim

Respiratory motion affects cardiac PET-computed tomography (CT) imaging by reducing attenuation correction (AC) accuracy and by introducing blur. The aim of this study was to compare three approaches for reducing motion-induced AC errors and evaluate the inclusion of respiratory motion correction.

Materials and methods

AC with a helical CT was compared with averaged cine and gated cine CT, as well as with a pseudo-gated CT, which was produced by applying PET-derived motion fields to the helical CT. Data-driven gating was used to produce respiratory-gated PET and CT images, and 60 NH₃ PET scans were attenuation corrected with each of the CTs. Respiratory motion correction was applied to the gated and pseudo-gated attenuation-corrected PET images.

Results

Anterior and lateral wall intensity measured in attenuation-corrected PET images generally increased when PET-CT alignment improved and decreased when alignment degraded. On average, all methods improved PET-CT liver and cardiac alignment, and increased anterior wall intensity by more than 10% in 36, 33 and 25 cases for the averaged, gated and pseudo-gated CTAC PET images, respectively. However, cases were found where alignment worsened and severe artefacts resulted. This occurred in more cases and to a greater extent for the averaged and gated CT, where the anterior wall intensity reduced by more than 10% in 21 and 24 cases, respectively, compared with six cases for the pseudo-gated CT. Application of respiratory motion correction increased the average anterior and inferior wall intensity, but only 13% of cases increased by more than 10%.

Conclusion

All methods improved average respiratory-induced AC errors; however, some severe artefacts were produced. The pseudo-gated CT was found to be the most robust method.

Estimating the 4D respiratory lung motion by spatiotemporal registration and super-resolution image reconstruction

PURPOSE

One of the main challenges in lung cancer radiation therapy is how to reduce the treatment margin but accommodate the geometric uncertainty of moving tumor. 4D-CT is able to provide the full range of motion information for the lung and tumor. However, accurate estimation of lung motion with respect to the respiratory phase is difficult due to various challenges in image registration, e.g., motion artifacts and large interslice thickness in 4D-CT. Meanwhile, the temporal coherence across respiration phases is usually not guaranteed in the conventional registration methods which consider each phase image in 4D-CT independently. To address these challenges, the authors present a unified approach to estimate the respiratory lung motion with two iterative steps.

METHODS

First, the authors propose a novel spatiotemporal registration algorithm to align all phase images of 4D-CT (in low-resolution) to a high-resolution group-mean image in the common space. The temporal coherence of registration is maintained by a set of temporal fibers that delineate temporal correspondences across different respiratory phases. Second, a super-resolution technique is utilized to build the high-resolution group-mean image with more anatomical details than any individual phase image, thus largely alleviating the registration uncertainty especially in correspondence detection. In particular, the authors use the concept of sparse representation to keep the group-mean image as sharp as possible.

RESULTS

The performance of our 4D motion estimation method has been extensively evaluated on both the simulated datasets and real lung 4D-CT datasets. In all experiments, our method achieves more accurate and consistent results in lung motion estimation than all other state-of-the-art approaches under comparison.

CONCLUSIONS

The authors have proposed a novel spatiotemporal registration method to estimate the lung motion in 4D-CT. Promising results have been obtained, which indicates the high applicability of our method in clinical lung cancer radiation therapy.

Reduction of irregular breathing artifacts in respiration-correlated CT images using a respiratory motion model

PURPOSE

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

MRI-based motion correction of thoracic PET: initial comparison of acquisition protocols and correction strategies suitable for simultaneous PET/MRI systems

OBJECTIVES

Magnetic resonance imaging (MRI) acquired on equipment capable of simultaneous MRI and positron emission tomography (PET) could potentially provide the gold standard method for motion correction of PET. To assess the latter, in this study we compared fast 2D and 3D MRI of the torso and used deformation parameters from real MRI data to correct simulated PET data for respiratory motion.

METHODS

PET sinogram data were simulated using SimSET from a 4D pseudo-PET image series created by segmenting MR images acquired over a respiratory cycle. Motion-corrected PET images were produced using post-reconstruction registration (PRR) and motion-compensated image reconstruction (MCIR).

RESULTS

MRI-based motion correction improved PET image quality at the lung-liver and lung-spleen boundaries and in the heart but little improvement was obtained where MRI contrast was low. The root mean square error in SUV units per voxel compared to a motion-free image was reduced from 0.0271 (no motion correction) to 0.0264 (PRR) and 0.0250 (MCIR).

CONCLUSIONS

Motion correction using MRI can improve thoracic PET images but there are limitations due to the quality of fast MRI.

Characterization and identification of spatial artifacts during 4D-CT imaging

PURPOSE

The purpose of this work is twofold: First, to characterize the artifacts occurring in helical 4D-CT imaging; second, to propose a method that can automatically identify the artifacts in 4D-CT images. The authors have designed a process that can automatically identify the artifacts in 4D-CT images, which may be invaluable in quantifying the quality of 4D-CT images and reducing the artifacts from the reconstructed images on a large dataset.

METHODS

Given two adjacent stacks obtained from the same respiration phase, the authors determine if there are artifacts between them. The proposed method uses a "bridge" stack strategy to connect the two stacks. Using normalized cross correlation convolution (NCCC), the two stacks are mapped to the bridge stack and the best matching positions can be located. Using this position information, the authors can then determine if there are artifacts between the two stacks. By combining the matching positions with NCCC values, the performance can be improved.

RESULTS

To validate the method, three expert observers independently labeled over 600 stacks on five patients. The results confirmed that high performance was obtained using the proposed method. The average sensitivity was about 0.87 and the average specificity was 0.82. The proposed method also outperformed the method of using respiratory signal (sensitivity increased from 0.50 to 0.87 and specificity increased from 0.70 to 0.82).

CONCLUSIONS

This study shows that the spatial artifacts during 4D-CT imaging are characterized and can be located automatically by the proposed method. The method is relatively simple but effective. It provides a way to evaluate the artifacts more objectively and accurately. The reported approach has promising potential for automatically identifying the types and frequency of artifacts on large scale 4D-CT image dataset.

Extension of a data-driven gating technique to 3D, whole body PET studies

Respiratory gating can be used to separate a PET acquisition into a series of near motion-free bins. This is typically done using additional gating hardware; however, software-based methods can derive the respiratory signal from the acquired data itself. The aim of this work was to extend a data-driven respiratory gating method to acquire gated, 3D, whole body PET images of clinical patients. The existing method, previously demonstrated with 2D, single bed-position data, uses a spectral analysis to find regions in raw PET data which are subject to respiratory motion. The change in counts over time within these regions is then used to estimate the respiratory signal of the patient. In this work, the gating method was adapted to only accept lines of response from a reduced set of axial angles, and the respiratory frequency derived from the lung bed position was used to help identify the respiratory frequency in all other bed positions. As the respiratory signal does not identify the direction of motion, a registration-based technique was developed to align the direction for all bed positions. Data from 11 clinical FDG PET patients were acquired, and an optical respiratory monitor was used to provide a hardware-based signal for comparison. All data were gated using both the data-driven and hardware methods, and reconstructed. The centre of mass of manually defined regions on gated images was calculated, and the overall displacement was defined as the change in the centre of mass between the first and last gates. The mean displacement was 10.3 mm for the data-driven gated images and 9.1 mm for the hardware gated images. No significant difference was found between the two gating methods when comparing the displacement values. The adapted data-driven gating method was demonstrated to successfully produce respiratory gated, 3D, whole body, clinical PET acquisitions.

Estimating view parameters from random projections for Tomography using spherical MDS

Background

During the past decade, the computed tomography has been successfully applied to various fields especially in medicine. The estimation of view angles for projections is necessary in some special applications of tomography, for example, the structuring of viruses using electron microscopy and the compensation of the patient's motion over long scanning period.

Methods

This work introduces a novel approach, based on the spherical multidimensional scaling (sMDS), which transforms the problem of the angle estimation to a sphere constrained embedding problem. The proposed approach views each projection as a high dimensional vector with dimensionality equal to the number of sampling points on the projection. By using SMDS, then each projection vector is embedded onto a 1D sphere which parameterizes the projection with respect to view angles in a globally consistent manner. The parameterized projections are used for the final reconstruction of the image through the inverse radon transform. The entire reconstruction process is non-iterative and computationally efficient.

Results

The effectiveness of the sMDS is verified with various experiments, including the evaluation of the reconstruction quality from different number of projections and resistance to different noise levels. The experimental results demonstrate the efficiency of the proposed method.

Conclusion

Our study provides an effective technique for the solution of 2D tomography with unknown acquisition view angles. The proposed method will be extended to three dimensional reconstructions in our future work. All materials, including source code and demos, are available on https://engineering.purdue.edu/PRECISE/SMDS.

Retrospective data-driven respiratory gating for PET/CT

Respiratory motion can adversely affect both PET and CT acquisitions. Respiratory gating allows an acquisition to be divided into a series of motion-reduced bins according to the respiratory signal, which is typically hardware acquired. In order that the effects of motion can potentially be corrected for, we have developed a novel, automatic, data-driven gating method which retrospectively derives the respiratory signal from the acquired PET and CT data. PET data are acquired in listmode and analysed in sinogram space, and CT data are acquired in cine mode and analysed in image space. Spectral analysis is used to identify regions within the CT and PET data which are subject to respiratory motion, and the variation of counts within these regions is used to estimate the respiratory signal. Amplitude binning is then used to create motion-reduced PET and CT frames. The method was demonstrated with four patient datasets acquired on a 4-slice PET/CT system. To assess the accuracy of the data-derived respiratory signal, a hardware-based signal was acquired for comparison. Data-driven gating was successfully performed on PET and CT datasets for all four patients. Gated images demonstrated respiratory motion throughout the bin sequences for all PET and CT series, and image analysis and direct comparison of the traces derived from the data-driven method with the hardware-acquired traces indicated accurate recovery of the respiratory signal.

Quality and accuracy of cone beam computed tomography gated by active breathing control

The purpose of this study was to evaluate the quality and accuracy of cone beam computed tomography (CBCT) gated by active breathing control (ABC), which may be useful for image guidance in the presence of respiration. Comparisons were made between conventional ABC-CBCT (stop and go), fast ABC-CBCT (a method to speed up the acquisition by slowing the gantry instead of stopping during free breathing), and free breathing respiration correlated CBCT. Image quality was assessed in phantom. Accuracy of reconstructed voxel intensity, uniformity, and root mean square error were evaluated. Registration accuracy (bony and soft tissue) was quantified with both an anthropomorphic and a quality assurance phantom. Gantry angle accuracy was measured with respect to gantry speed modulation. Conventional ABC-CBCT scan time ranged from 2.3 to 5.8 min. Fast ABC-CBCT scan time ranged from 1.4 to 1.8 min, and respiratory correlated CBCT scans took 2.1 min to complete. Voxel intensity value for ABC gated scans was accurate relative to a normal clinical scan with all projections. Uniformity and root mean square error performance degraded as the number of projections used in the reconstruction of the fast ABC-CBCT scans decreased (shortest breath hold, longest free breathing segment). Registration accuracy for small, large, and rotational corrections was within 1 mm and 1 degrees. Gantry angle accuracy was within 1 degrees for all scans. For high-contrast targets, performance for image-guidance purposes was similar for fast and conventional ABC-CBCT scans and respiration correlated CBCT.