Helical cardiac cone beam reconstruction using retrospective ECG gating

Cardiac Motion Correction for Helical CT Scan With an Ordinary Pitch

Cardiac X-ray computed tomography (CT) imaging is still challenging due to the cardiac motion during CT scanning, which leads to the presence of motion artifacts in the reconstructed image. In response, many cardiac X-ray CT imaging algorithms have been proposed, based on motion estimation (ME) and motion compensation (MC), to improve the image quality by alleviating the motion artifacts in the reconstructed image. However, these ME/MC algorithms are mainly based on an axial scan or a low-pitch helical scan. In this paper, we propose a ME/MC-based cardiac imaging algorithm for the data set acquired from a helical scan with an ordinary pitch of around 1.0 so as to obtain the whole cardiac image within a single scan of short time without ECG gating. In the proposed algorithm, a sequence of partial angle reconstructed (PAR) images is generated by using consecutive parts of the sinogram, each of which has a small angular span. Subsequently, an initial 4-D motion vector field (MVF) is obtained using multiple pairs of conjugate PAR images. The 4-D MVF is then refined based on an image quality metric so as to improve the quality of the motion-compensated image. Finally, a time-resolved cardiac image is obtained by performing motion-compensated image reconstruction by using the refined 4-D MVF. Using digital XCAT phantom data sets and a human data set commonly obtained via a helical scan with a pitch of 1.0, we demonstrate that the proposed algorithm significantly improves the image quality by alleviating motion artifacts.

Cardiac Image Reconstruction via Nonlinear Motion Correction Based on Partial Angle Reconstructed Images

Even though the X-ray Computed Tomography (CT) scan is considered suitable for fast imaging, motion-artifact-free cardiac imaging is still an important issue, because the gantry rotation speed is not fast enough compared with the heart motion. To obtain a heart image with less motion artifacts, a motion estimation (ME) and motion compensation (MC) approach is usually adopted. In this paper, we propose an ME/MC algorithm that can estimate a nonlinear heart motion model from a sinogram with a rotation angle of less than 360°. In this algorithm, we first assume the heart motion to be nonrigid but linear, and thereby estimate an initial 4-D motion vector field (MVF) during a half rotation by using conjugate partial angle reconstructed images, as in our previous ME/MC algorithm. We then refine the MVF to determine a more accurate nonlinear MVF by maximizing the information potential of a motion-compensated image. Finally, MC is performed by incorporating the determined MVF into the image reconstruction process, and a time-resolved heart image is obtained. By using a numerical phantom, a physical cardiac phantom, and an animal data set, we demonstrate that the proposed algorithm can noticeably improve the image quality by reducing motion artifacts throughout the image.

Data-Driven Affine Deformation Estimation and Correction in Cone Beam Computed Tomography

In computed tomography (CT), motion and deformation during the acquisition lead to streak artefacts and blurring in the reconstructed images. To remedy these artefacts, we introduce an efficient algorithm to estimate and correct for global affine deformations directly on the cone beam projections. The proposed technique is data driven and thus removes the need for markers and/or a tracking system. A relationship between affine transformations and the cone beam transform is proved and used to correct the projections. The deformation parameters that describe deformation perpendicular to the projection direction are estimated for each projection by minimizing a plane-based inconsistency criterion. The criterion compares each projection of the main scan with all projections of a fast reference scan, which is acquired prior or posterior to the main scan. Experiments with simulated and experimental data show that the proposed affine deformation estimation method is able to substantially reduce motion artefacts in cone beam CT images.

Completeness map evaluation demonstrated with candidate next-generation cardiac CT architectures

PURPOSE

In this report, the authors introduce the general concept of the completeness map, as a means to evaluate the completeness of data acquired by a given CT system design (architecture and scan mode). They illustrate the utility of completeness map by applying the completeness map concept to a number of candidate CT system designs, as part of a study to advance the state-of-the-art in cardiac CT.

METHODS

In order to optimally reconstruct a point within a volume of interest (VOI), the Radon transform on all possible planes through that point should be measured. The authors quantified the extent to which this ideal condition is satisfied for the entire image volume. They first determined a Radon completeness number for each point in the VOI, as the percentage of possible planes that is actually measured. A completeness map is then defined as a 3D matrix of the completeness numbers for the entire VOI. The authors proposed algorithms to analyze the projection datasets in Radon space and compute the completeness number for a fixed point and apply these algorithms to various architectures and scan modes that they are evaluating. In this report, the authors consider four selected candidate architectures, operating with different scan modes, for a total of five system design alternatives. Each of these alternatives is evaluated using completeness map.

RESULTS

If the detector size and cone angle are large enough to cover the entire cardiac VOI, a single-source circular scan can have ≥99% completeness over the entire VOI. However, only the central z-slice can be exactly reconstructed, which corresponds to 100% completeness. For a typical single-source architecture, if the detector is limited to an axial dimension of 40 mm, a helical scan needs about five rotations to form an exact reconstruction region covering the cardiac VOI, while a triple-source helical scan only requires two rotations, leading to a 2.5x improvement in temporal resolution. If the source and detector of an inverse-geometry (IGCT) system have the same axial extent, and the spacing of source points in the axial and transaxial directions is sufficiently small, the IGCT can also form an exact reconstruction region for the cardiac VOI. If the VOI can be covered by the x-ray beam in any view, a composite-circling scan can generate an exact reconstruction region covering the VOI.

CONCLUSIONS

The completeness map evaluation provides useful information for selecting the next-generation cardiac CT system design. The proposed completeness map method provides a practical tool for analyzing complex scanning trajectories, where the theoretical image quality for some complex system designs is impossible to predict, without yet-undeveloped reconstruction algorithms.

Sequential Turning Acquisition and Reconstruction (STAR) method for four-dimensional imaging of cyclically moving structures

Optical coherence tomography allows for dynamic, three-dimensional (3D+T) imaging of the heart within animal embryos. However, direct 3D+T imaging frame rates remain insufficient for cardiodynamic analysis. Previously, this limitation has been addressed by reconstructing 3D+T representations of the beating heart based on sets of two-dimensional image sequences (2D+T) acquired sequentially at high frame rate and in fixed (and parallel) planes throughout the heart. These methods either require additional hardware to trigger the acquisition of each 2D+T series to the same phase of the cardiac cycle or accumulate registration errors as the slices are synchronized retrospectively by pairs, without a gating signal. Here, we present a sequential turning acquisition and reconstruction (STAR) method for 3D+T imaging of periodically moving structures, which does not require any additional gating signal and is not prone to registration error accumulation. Similarly to other sequential cardiac imaging methods, multiple fast image series are consecutively acquired for different sections but in between acquisitions, the imaging plane is rotated around the center line instead of shifted along the direction perpendicular to the slices. As the central lines of all image-sequences coincide and represent measurements of the same spatial position, they can be used to accurately synchronize all the slices to a single inherent reference signal. We characterized the accuracy of our method on a simulated dynamic phantom and successfully imaged a beating embryonic rat heart. Potentially, this method can be applied for structural or Doppler imaging approaches with any direct space imaging modality such as computed tomography, ultrasound, or light microscopy.

Adaptive two-pass cone-beam artifact correction using a FOV-preserving two-source geometry: a simulation study

The evolution to ever wider detector arrays that are able to cover whole organs with a single circular gantry sweep has revitalized the research efforts toward finding improved axial scanning algorithms and protocols. The authors propose a computed tomography scan and reconstruction concept using two sources, a single detector and a two-pass cone-beam correction method, as an integral part of the reconstruction. Compared with standard circular acquisition and reconstruction methods, the new concept excels with improved coverage and very low cone-beam artifact level also for short scan acquisitions, which makes it especially attractive for cardiac applications.

Improvement of cardiac CT reconstruction using local motion vector fields

The motion of the heart is a major challenge for cardiac imaging using CT. A novel approach to decrease motion blur and to improve the signal to noise ratio is motion compensated reconstruction which takes motion vector fields into account in order to correct motion. The presented work deals with the determination of local motion vector fields from high contrast objects and their utilization within motion compensated filtered back projection reconstruction. Image registration is applied during the quiescent cardiac phases. Temporal interpolation in parameter space is used in order to estimate motion during strong motion phases. The resulting motion vector fields are during image reconstruction. The method is assessed using a software phantom and several clinical cases for calcium scoring. As a criterion for reconstruction quality, calcium volume scores were derived from both, gated cardiac reconstruction and motion compensated reconstruction throughout the cardiac phases using low pitch helical cone beam CT acquisitions. The presented technique is a robust method to determine and utilize local motion vector fields. Motion compensated reconstruction using the derived motion vector fields leads to superior image quality compared to gated reconstruction. As a result, the gating window can be enlarged significantly, resulting in increased SNR, while reliable Hounsfield units are achieved due to the reduced level of motion artefacts. The enlargement of the gating window can be translated into reduced dose requirements.

Motion-compensated iterative cone-beam CT image reconstruction with adapted blobs as basis functions

This paper presents a three-dimensional method to reconstruct moving objects from cone-beam X-ray projections using an iterative reconstruction algorithm and a given motion vector field. For the image representation, adapted blobs are used, which can be implemented efficiently as basis functions. Iterative reconstruction requires the calculation of line integrals (forward projections) through the image volume, which are compared with the actual measurements to update the image volume. In the existence of a divergent motion vector field, a change in the volumes of the blobs has to be taken into account in the forward and backprojections. An efficient method to calculate the line integral through the adapted blobs is proposed. It solves the problem, how to compensate for the divergence in the motion vector field on a grid of basis functions. The method is evaluated on two phantoms, which are subject to three different known motions. Moreover, a motion-compensated filtered back-projection reconstruction method is used, and the reconstructed images are compared. Using the correct motion vector field with the iterative motion-compensated reconstruction, sharp images are obtained, with a quality that is significantly better than gated reconstructions.

Effect of novel amplitude/phase binning algorithm on commercial four-dimensional computed tomography quality

PURPOSE

Respiratory motion is a significant source of anatomic uncertainty in radiotherapy planning and can result in errors of portal size and the subsequent radiation dose. Although four-dimensional computed tomography allows for more accurate analysis of the respiratory cycle, breathing irregularities during data acquisition can cause considerable image distortions. The aim of this study was to examine the effect of respiratory irregularities on four-dimensional computed tomography, and to evaluate a novel image reconstruction algorithm using percentile-based tagging of the respiratory cycle.

METHODS AND MATERIALS

Respiratory-correlated helical computed tomography scans were acquired for 11 consecutive patients. The inspiration and expiration data sets were reconstructed using the default phase-based method, as well as a novel respiration percentile-based method with patient-specific metrics to define the ranges of the reconstruction. The image output was analyzed in a blinded fashion for the phase- and percentile-based reconstructions to determine the prevalence and severity of the image artifacts.

RESULTS

The percentile-based algorithm resulted in a significant reduction in artifact severity compared with the phase-based algorithm, although the overall artifact prevalence did not differ between the two algorithms. The magnitude of differences in respiratory tag placement between the phase- and percentile-based algorithms correlated with the presence of image artifacts.

CONCLUSION

The results of our study have indicated that our novel four-dimensional computed tomography reconstruction method could be useful in detecting clinically relevant image distortions that might otherwise go unnoticed and to reduce the image distortion associated with some respiratory irregularities. Additional work is necessary to assess the clinical impact on areas of possible irregular breathing.

Diagnostic accuracy of myocardial hypoenhancement on multidetector computed tomography in identifying myocardial infarction in patients admitted with acute chest pain syndrome

OBJECTIVE

To evaluate prevalence and diagnostic accuracy of myocardial hypoenhancement (MH) using multidetector computed tomography (MDCT) in patients admitted for acute chest pain syndromes.

METHODS

Sixty-nine patients underwent first-pass MDCT, coronary angiography, and echocardiography. Using a standardized analysis protocol, left ventricular short-axis reformations were evaluated for presence, size, and density of MH in 16 myocardial segments. These were correlated with the presence and location of myocardial infarction (MI), regional myocardial dysfunction, and coronary artery disease.

RESULTS

Myocardial hypoenhancement was found in acute MI (27/35), healed MI (6/14), unstable angina (3/9), and atypical chest pain (0/11). Sensitivity, specificity, and positive and negative predictive values of MH for diagnosing any MI were 67%, 85%, 92% and 52%, respectively.

CONCLUSIONS

The presence of MH on MDCT in acute chest pain patients has high positive predictive value and specificity but only moderate sensitivity for presence of acute or healed MI using the strict criteria proposed in this study.

Is functional improvement after myocardial infarction predicted with myocardial enhancement patterns at multidetector CT?

PURPOSE

To prospectively evaluate the sensitivity of myocardial early perfusion defects (EDs) and late enhancement (LE) at multidetector computed tomography (CT) following acute myocardial infarction (AMI) to predict segment myocardial dysfunction and myocardial functional recovery (MFR), by using echocardiography as the reference standard.

MATERIALS AND METHODS

Institutional review board approval and informed consent were obtained. Twenty-six patients (25 men, one woman; mean age, 53 years+/-9 [standard deviation]), underwent baseline multidetector CT, coronary angiography, and echocardiography within a week of AMI and a follow-up echocardiography at 3 months. ED, LE, and late hypoattenuation were compared with regional left ventricular function and MFR. A logistic regression model and generalized estimating equation analysis were applied to estimate the predictive effect of ED and LE. Differences between groups were evaluated by using nonpaired Student t tests.

RESULTS

All EDs and LE corresponded with AMI location determined by using angiography and echocardiography. For occluded arteries (n=5), no relationship was found between the presence of ED or LE and MFR. For patent arteries (n=21), presence of LE had a respective sensitivity and specificity of 73% and 85% for predicting follow-up segment dysfunction, compared with 57% and 90% for ED. In abnormal baseline segments, nonrecovery was clearly related to the presence and size of segment defect area for both ED (odds ratio: 1.95 [95% confidence interval: 0.9, 4.1] per square centimeter) and LE (odds ratio: 1.85 [95% confidence interval: 1.2, 2.9] per square centimeter). Segments that recovered had significantly lower prevalence of ED and LE, and if present, were significantly smaller than in segments remaining abnormal (P<.05).

CONCLUSION

The presence and size of ED and LE at multidetector CT is closely related to follow-up segment myocardial dysfunction and MFR.

Effective doses in subjects undergoing computed tomography cardiac imaging with the 256-multislice CT scanner

BACKGROUND

The 256-multislice CT (256MSCT) obtains volumetric data with 128-mm coverage in a single rotation. This coverage allows satisfactory visualization of the whole heart, allowing the 256MSCT to visualize the cardiac chambers and coronary arteries by cine scan without ECG gating. These characteristics provide a solution to the problems of MSCT. Although a wider beam width provides more efficient imaging over a wider coverage area, patient doses with the 256MSCT are of considerable concern.

OBJECTIVE

We assessed potential radiation exposure with the 256MSCT in a cardiac CT protocol and compared the results to those with 16- and 64MSCT (collimated 64x0.5mm using 256MSCT).

METHODS

Organ or tissue doses were measured in an anthropomorphic phantom under a coronary artery imaging protocol with the 256MSCT in cine scan mode without ECG gating, and with the 16- and 64MSCT in helical scan mode with ECG gating.

RESULTS

Average effective doses were 22.8mSv for the 16MSCT, 27.8mSv for the 64MSCT and 14.1mSv for the 256MSCT. The 16- and 64MSCT doses were thus approximately 1.6- and 2.0-fold higher than those of the 256MSCT.

CONCLUSIONS

Use of the 256MSCT in cardiac volumetric cine imaging offers lower radiation exposure than 16- and 64MSCT, and suggests the potential of this equipment in single-beat cardiac imaging without ECG gating. This effective dose is acceptable for routine cardiac imaging.

Quantitative assessment of four-dimensional computed tomography image acquisition quality

The purpose of the present work was to describe the development and validation of a series of tests to assess the quality of four-dimensional (4D) computed-tomography (CT) imaging as it is applied to radiation treatment planning. Using a commercial respiratory motion phantom and a programmable moving platform with a CT phantom, we acquired 4D CT datasets on two commercial multislice helical CT scanners that use different approaches to 4D CT image reconstruction. Datasets were obtained as the platform moved in various patterns designed to simulate breathing. Known inserts in the phantom were contoured, and statistics were generated to evaluate properties important to radiation therapy--namely, accuracy of phase-binning, shape, volume, and CT number. Phase-binning accuracy varied by as much as 5% for a 4D procedure in which images were reconstructed and then binned, but exhibited no variation for a 4D procedure in which projections were binned before reconstruction. The magnitude of geometric distortion was found to be small for both approaches, as was the magnitude of volume error. Partial-volume effects in the direction perpendicular to the transverse planes of reconstruction affected volume accuracy, however. Computed tomography numbers were reproduced accurately, but 4D images exhibited more variation in CT number than static CT images did. Characterization of such properties can be used to better understand and optimize the various parameters that affect 4D CT image quality.

Evaluation of the influence of acquisition and reconstruction parameters for 16-row multidetector CT on coronary calcium scoring using a stationary and dynamic cardiac phantom

A calcium-scoring phantom with hydroxyapatite-filled cylindrical holes (0.5 to 4 mm) was used. High-resolution scans were performed for an accuracy baseline. The phantom was mounted to a moving heart phantom. Non-moving data with the implementation of an ECG-signal were acquired for different pitches (0.2/0.3), heart rates (60/80/95 bpm) and collimations (16 x 0.75/16 x 1.5 mm). Images were reconstructed with a cone-beam multi-cycle algorithm at a standard thickness/increment of 3 mm/1.5 mm and the thinnest possible thickness (0.8/0.4 and 2/1). Subsequently, ECG-gated moving calcium-scoring phantom data were acquired. The calcium volume and Agatston score were measured. The temporal resolution and reconstruction cycles were calculated. High-resolution scans determine the calcium volume with a high accuracy (mean overestimation, 0.8%). In the non-moving measurements, the volume underestimation ranged from about 6% (16 x 0.75 mm; 0.8/0.4 mm) to nearly 25% (16 x 1.5 mm; 3/1.5 mm). Moving scans showed increased measurement errors depending on the reconstructed RR interval, collimation, pitch, heart rate and gantry rotation time. Also, a correlation with the temporal resolution could be found. The reliability of calcium-scoring results can be improved with the use of a narrower collimation, a lower pitch and the reconstruction of thinner images, resulting in higher patient doses. The choice of the correct cardiac phase within the RR interval is essential to minimize measurement errors.

Selectable Source Rotational Velocity for Cardiac Computed Tomography

The article analyzes the existing algorithms of cardiac imaging reconstructions using multidata segments and derives the ideal data segment conditions. A scheme of selectable velocity scanning is then proposed for cardiac computed tomography with a circular trajectory, which can cover any cardiac phase exactly or approximately to achieve a desirable temporal resolution. Numerical tests show that the selectable source rotational velocity method outperforms the traditional multidata segment algorithms. The proposed method can be also extended into the cases of helical scanning and variable heart rates.

Algebraic and analytic reconstruction methods for dynamic tomography

In this work, we discuss algebraic and analytic approaches for dynamic tomography. We present a framework of dynamic tomography for both algebraic and analytic approaches. We finally present numerical experiments.

Improvement in banding artefacts in four-dimensional computed tomography for radiotherapy planning

Respiratory-gated CT (RGCT) and four-dimensional CT (4DCT) scan techniques cover consecutive segments of the respiratory cycle. However, motion artefacts may occur in fast respiratory phases such as mid-inhalation and -exhalation. CT imaging involves the use of a number of x-ray tube positions for each couch position. We investigated the fundamental nature of motion artefacts using a constant-velocity moving phantom in motion in the CT plane or perpendicular to the CT plane, and in pigs to simulate a human model. Artefacts and movement distance were evaluated in a moving phantom and artificially ventilated pigs with a 256-multi-detector row CT (256MDCT). The phantom moved in the CT plane or perpendicular to the CT plane with a constant velocity. Backprojection used variable initial backprojection angles (IBAs). The phantom length for motion perpendicular to the CT plane was independent of IBA but was represented by phantom diameter plus the distance of movement per gantry rotation. In contrast, that for the motion in the CT plane was dependent on IBA, as represented by phantom diameter plus the distance of movement per rotation for IBA perpendicular to the phantom movement direction, and phantom diameter plus half the distance of movement per gantry rotation for other IBAs. Results for volumetric CT images with different IBAs showed the presence of banding artefacts. Similar findings were seen in artificially ventilated pigs. Motion artefacts are unavoidable in both conventional CT and 256MDCT. Banding artefacts will be improved if the same IBAs at each couch position are accounted for during image reconstruction. This improvement will be beneficial in respiratory gated and 4D radiation therapies.

CEnPiT: Helical cardiac CT reconstruction

Computer tomography (CT) scanners with an increasing number of detector rows offer the potential of shorter scanning times. Nevertheless, the reconstruction problem becomes more challenging, since cone beam artifacts are likely to enter. Here, we consider helical cardiac CT. We analyze how a relationship can be established between exact reconstruction algorithms and the demand to perform a cardiac gating. Utilizing the redundancies requires the consideration of all kinds of Radon planes. For the reconstruction algorithm proposed here, we separate the data into two parts. The first part contains contributions of Radon planes, which are measured with a large number of redundancies. The second part contains the remaining contributions. As it turns out, the second part contributes rather to the low-frequency contents of trans-axial slices. Therefore, we propose to perform a gated back-projection only for the first part, while the second part is back-projected in an ungated way. Data from the complete source trajectory are employed in the reconstruction process in contrary to conventional helical cardiac reconstruction methods. Moreover, all different types of Radon planes are taken into account in the reconstruction, though an ECG-dependent cardiac gating is applied. The reconstruction results, which we present for clinical and simulated data, demonstrate the high potential of CEnPiT for helical cardiac CT with large cone angle systems.

Aperture weighted cardiac reconstruction for cone-beam CT

Multi-row detectors together with fast rotating gantries made cardiac imaging possible for CT. Due to the cardiac motion, ECG gating has to be integrated into the reconstruction of the data measured on a low pitch helical trajectory. Since the first multi-row scanners were introduced, it has been shown that approximative true cone-beam reconstruction methods are most suitable for the task of retrospectively gated cardiac volume CT. In this paper, we present the aperture weighted cardiac reconstruction (AWCR), which is a three-dimensional reconstruction algorithm of the filtered back-projection type. It is capable of handling all illumination intervals of an object point, which occur as a consequence of a low pitch helical cone-beam acquisition. Therefore, this method is able to use as much redundant data as possible, resulting in an improvement of the image homogeneity, the signal to noise ratio and the temporal resolution. Different optimization techniques like the heart rate adaptive cardiac weighting or the automatic phase determination can be adopted to AWCR. The excellent image quality achieved by AWCR is presented for medical datasets acquired with both a 40-slice and a 64-slice cone-beam CT scanner.

Multithreaded cardiac CT

Phase-correlated CT, as it is used for cardiac imaging, is the most popular and the most important but also the most demanding special CT application in the clinical routine, today. Basically, it fulfills the four-dimensional imaging task of depicting a quasiperiodically moving object at any desired motion phase with significantly reduced motion artifacts. Although image quality with phase-correlated reconstruction is far better than with standard reconstruction, there are motion artifacts remaining and improvements of temporal resolution are required. As a well-known alternative to simply decreasing rotation time, we consider a spiral cone-beam CT scanner that has G x-ray guns and detectors mounted. We call this a multisource or a multithreaded CT scanner. Aiming for improved temporal resolution the relative temporal resolution tau, which measures the fraction of a motion period that enters the image, is studied as a function of the motion rate (heart rate) and the degree of scan overlap (pitch value) for various configurations. The parameters to optimize are the number of threads G and the interthread parameters delta alpha and delta z, which are the angular and the longitudinal separation between adjacent threads, respectively. To demonstrate the improvements approximate image reconstruction of multithreaded raw data is performed by using a generalization of the extended parallel back projection cone-beam reconstruction algorithm [Med. Phys. 31(6), 1623-1641 (2004)] to the case of multithreaded CT. Reconstructions of a simulated cardiac motion phantom and of simulated semi-antropomorphic phantoms are presented for two and three threads and compared to the single-threaded case to demonstrate the potential of multithreaded cardiac CT. Patient data were acquired using a clinical double-threaded CT scanner to validate the theoretical results. The optimum angle delta alpha between the tubes is 90 degrees for a double-threaded system, and for triple-threaded scanners it is 60 degrees or 120 degrees. In all cases, delta z = 0 results as an optimum, which means that the threads should be mounted in the same transversal plane. However, the dependency of the temporal resolution on delta z is very weak and a longitudinal separation delta z not = 0 would not deteriorate image quality. The mean temporal resolution achievable with an optimized multithreaded CT scanner is a factor of G better than the mean temporal resolution obtained with a single-threaded scanner. The standard reconstructions showed decreased cone-beam artifacts with multithreaded CT compared to the single-threaded case. Our phase-correlated reconstructions demonstrate that temporal resolution is significantly improved with multithreaded CT. The clinical patient data confirm our results.

Improved motion compensation in 3D-CT using respiratory-correlated segment reconstruction: diagnostic and radiotherapy applications

Conventional respiratory-gated CT and four-dimensional CT (4DCT) are disadvantaged by their low temporal resolution, which results in the inclusion of anatomic motion-induced artefacts. These represent a significant source of error both in radiotherapy treatment planning for the thorax and upper abdomen and in diagnostic procedures. In particular, temporal resolution and image quality are vitally important to accurate diagnosis and the minimization of planning target volume margin due to respiratory motion. To improve both temporal resolution and signal-to-noise ratio (SNR), we developed a respiratory-correlated segment reconstruction method (RS) and adapted it to the Feldkamp-Davis-Kress algorithm (FDK) with a 256 multidetector row CT (256MDCT). The 256MDCT scans approximately 100 mm in the craniocaudal direction with a 0.5 mm slice thickness in one rotation. Data acquisition for the RS-FDK relies on the assistance of a respiratory sensing system operating in cine scan mode (continuous axial scan with the table stationary). We evaluated the RS-FDK for volume accuracy and image noise in a phantom study with the 256MDCT and compared results with those for a full scan (FS-FDK), which is usually employed in conventional 4DCT and in half scan (HS-FDK). Results showed that the RS-FDK gave a more accurate volume than the others and had the same SNR as the FS-FDK. In a subsequent animal study, we demonstrated a practical sorting process for projection data which was unaffected by variations in respiratory period, and found that the RS-FDK gave the clearest visualization among the three algorithms of the margins of the liver and pulmonary vessels. In summary, the RS-FDK algorithm provides multi-phase images with higher temporal resolution and better SNR. This method should prove useful when combined with new radiotherapeutic and diagnostic techniques.

Direct cone-beam cardiac reconstruction algorithm with cardiac banding artifact correction

Multislice helical computed tomography (CT) is a promising noninvasive technique for coronary artery imaging. Various factors can cause inconsistencies in cardiac CT data, which can result in degraded image quality. These inconsistencies may be the result of the patient physiology (e.g., heart rate variations), the nature of the data (e.g., cone-angle), or the reconstruction algorithm itself. An algorithm which provides the best temporal resolution for each slice, for example, often provides suboptimal image quality for the entire volume since the cardiac temporal resolution (TRc) changes from slice to slice. Such variations in TRc can generate strong banding artifacts in multiplanar reconstruction images or three-dimensional images. Discontinuous heart walls and coronary arteries may compromise the accuracy of the diagnosis. A beta-blocker is often used to reduce and stabilize patients' heart rate but cannot eliminate the variation. In order to obtain robust and optimal image quality, a software solution that increases the temporal resolution and decreases the effect of heart rate is highly desirable. This paper proposes an ECG-correlated direct cone-beam reconstruction algorithm (TCOT-EGR) with cardiac banding artifact correction (CBC) and disconnected projections redundancy compensation technique (DIRECT). First the theory and analytical model of the cardiac temporal resolution is outlined. Next, the performance of the proposed algorithms is evaluated by using computer simulations as well as patient data. It will be shown that the proposed algorithms enhance the robustness of the image quality against inconsistencies by guaranteeing smooth transition of heart cycles used in reconstruction.

Displacement-based binning of time-dependent computed tomography image data sets

Respiration can cause tumors in the thorax or abdomen to move by as much as 3 cm; this movement can adversely affect the planning and delivery of radiation treatment. Several techniques have been used to compensate for respiratory motion, but all have shortcomings. Manufacturers of computed tomography (CT) equipment have recently used a technique developed for cardiac CT imaging to track respiratory-induced anatomical motion and to sort images according to the phase of the respiratory cycle they represent. Here we propose a method of generating CT images that accounts for respiratory-induced anatomical motion on the basis of displacement, i.e., displacement-binned CT image sets. This technique has shown great promise, however, it is not fully supported by currently used CT image reconstruction software. As an interim solution, we have developed a method for extracting displacement-binned CT image data sets from data sets assembled on the basis of a prospectively determined breathing phase acquired on a multislice helical CT scanner. First, the projection data set acquired from the CT scanner was binned at small phase intervals before reconstruction. The manufacturer's software then generated image sets identified as belonging to particular phases of the respiratory cycle. All images were then individually correlated to the displacement of an external fiducial marker. Next, CT image data sets were resorted on the basis of the displacement and assigned an appropriate phase. Finally, displacement-binned image data sets were transferred to a treatment-planning system for analysis. Although the technique is currently limited by the phase intervals allowed by the CT software, some improvement in image reconstruction was seen, indicating that this technique is useful at least as an interim measure.

Accurate estimation of global and regional cardiac function by retrospectively gated multidetector row computed tomography

Retrospective reconstruction of ECG-gated images at different parts of the cardiac cycle allows the assessment of cardiac function by multi-detector row CT (MDCT) at the time of non-invasive coronary imaging. We compared the accuracy of such measurements by MDCT to cine magnetic resonance (MR). Forty patients underwent the assessment of global and regional cardiac function by 16-slice MDCT and cine MR. Left ventricular (LV) end-diastolic and end-systolic volumes estimated by MDCT (134+/-51 and 67+/-56 ml) were similar to those by MR (137+/-57 and 70+/-60 ml, respectively; both P=NS) and strongly correlated (r=0.92 and r=0.95, respectively; both P<0.001). Consequently, LV ejection fractions by MDCT and MR were also similar (55+/-21 vs. 56+/-21%; P=NS) and highly correlated (r=0.95; P<0.001). Regional end-diastolic and end-systolic wall thicknesses by MDCT were highly correlated (r=0.84 and r=0.92, respectively; both P<0.001), but significantly lower than by MR (8.3+/-1.8 vs. 8.8+/-1.9 mm and 12.7+/-3.4 vs. 13.3+/-3.5 mm, respectively; both P<0.001). Values of regional wall thickening by MDCT and MR were similar (54+/-30 vs. 51+/-31%; P=NS) and also correlated well (r=0.91; P<0.001). Retrospectively gated MDCT can accurately estimate LV volumes, EF and regional LV wall thickening compared to cine MR.

Cardiac imaging using multislice computed tomography scanners: technical considerations

Conventional coronary angiography is currently the gold standard in the detection and diagnosis of coronary artery disease. This modality, however, is invasive in nature. Hence, there is a need for noninvasive imaging techniques to provide comprehensive assessment of coronary artery disease, especially in stable patients at low to moderate risk of disease. In recent years, a number of noninvasive modalities have found wide applications in cardiac imaging. Most recent investigations have used magnetic resonance imaging, multislice computed tomography and electron-beam computed tomography scanners. This review discusses the clinical challenges existing in the field of cardiac imaging and focuses on the technical advancements of multislice computed tomography scanners that have made them a very attractive noninvasive option for the detection and diagnosis of coronary artery disease.

Tilted plane Feldkamp type reconstruction algorithm for spiral cone beam CT

An approximate image reconstruction method for spiral cone beam computed tomography (CT), called tilted plane Feldkamp type reconstruction algorithm (TPFR), is presented in this paper, which extends Feldkamp cone beam reconstruction algorithm to deal with its inaccuracy and artifact problems caused by large cone angle. This is done by tilting the reconstructing planes to minimize the cone angle and optimally fit the spiral segment of the source. The tilted plane image reconstruction requires reforming the three-dimensional projection data set for the tilted plane and application of Feldkamp algorithm to the reformed data set. Analytical and computational results can show that the image reconstruction performance of the proposed TPFR algorithm is superior to that of the Feldkamp reconstruction algorithm in the image quality, volume coverage speed, maximum achievable pitch value, and slice sensitivity profiles. Moreover, it provides more accurate image reconstruction than the existing two-dimensional reconstruction algorithms.

Four-dimensional cardiac imaging in living embryos via postacquisition synchronization of nongated slice sequences

Being able to acquire, visualize, and analyze 3D time series (4D data) from living embryos makes it possible to understand complex dynamic movements at early stages of embryonic development. Despite recent technological breakthroughs in 2D dynamic imaging, confocal microscopes remain quite slow at capturing optical sections at successive depths. However, when the studied motion is periodic--such as for a beating heart--a way to circumvent this problem is to acquire, successively, sets of 2D+time slice sequences at increasing depths over at least one time period and later rearrange them to recover a 3D+time sequence. In other imaging modalities at macroscopic scales, external gating signals, e.g., an electro-cardiogram, have been used to achieve proper synchronization. Since gating signals are either unavailable or cumbersome to acquire in microscopic organisms, we have developed a procedure to reconstruct volumes based solely on the information contained in the image sequences. The central part of the algorithm is a least-squares minimization of an objective criterion that depends on the similarity between the data from neighboring depths. Owing to a wavelet-based multiresolution approach, our method is robust to common confocal microscopy artifacts. We validate the procedure on both simulated data and in vivo measurements from living zebrafish embryos.

Cardiac cone-beam CT volume reconstruction using ART

Modern computed tomography systems allow volume imaging of the heart. Up to now, approximately two-dimensional (2D) and 3D algorithms based on filtered backprojection are used for the reconstruction. These algorithms become more sensitive to artifacts when the cone angle of the x-ray beam increases as it is the current trend of computed tomography (CT) technology. In this paper, we investigate the potential of iterative reconstruction based on the algebraic reconstruction technique (ART) for helical cardiac cone-beam CT. Iterative reconstruction has the advantages that it takes the cone angle into account exactly and that it can be combined with retrospective cardiac gating fairly easily. We introduce a modified ART algorithm for cardiac CT reconstruction. We apply it to clinical cardiac data from a 16-slice CT scanner and compare the images to those obtained with a current analytical reconstruction method. In a second part, we investigate the potential of iterative reconstruction for a large area detector with 256 slices. For the clinical cases, iterative reconstruction produces excellent images of diagnostic quality. For the large area detector, iterative reconstruction produces images superior to analytical reconstruction in terms of cone-beam artifacts.

EnPiT: filtered back-projection algorithm for helical CT using an n-Pi acquisition

In this paper, we formulate a reconstruction algorithm for an n-Pi acquisition, where n can be any positive odd integer. The algorithm is a generalization of the method presented in (Bontus et al. 2003). It is based on the results obtained by Katsevich (2004). For the algorithm, different sets of filter-lines have to be defined. We describe the variation of these lines along the detector in some detail, before we discuss, how the method gives all Radon-plane contributions the correct weighting. The different sets of filter-lines are all contained within the n-Pi window, such that a practical realization is possible. Reconstruction results, which we present in the final section, show convincing image quality.

Automatic determination of minimal cardiac motion phases for computed tomography imaging: initial experience

Low motion phases for cardiac computed tomography reconstructions are currently detected manually in a user-dependent selection process which is often time consuming and suboptimal. The concept of motion maps was recently introduced to achieve automatic phase selection. This pilot study compared the accuracy of motion-map phase selection to that with manual iterative selection. The study included 20 patients, consisting of one group with low and one with high heart rate. The technique automatically derives a motion strength function between multiple low-resolution reconstructions through the cardiac cycle, with periods of lowest difference between neighboring phases indicating minimal cardiac motion. A high level of agreement was found for phase selection achieved with the motion map approach compared with the manual iterative selection process. The motion maps allowed automated quiescent phase detection of the cardiac cycle in 85% of cases, with best results at low heart rates and for the left coronary artery. They can also provide additional information such as the presence of breathing artifacts. Motion maps show promise as a rapid off-line tool to automatically detect quiescent cardiac phases in a variety of patients.

Semiautomated four-dimensional computed tomography segmentation using deformable models

The purpose of this work is to demonstrate a proof of feasibility of the application of a commercial prototype deformable model algorithm to the problem of delineation of anatomic structures on four-dimensional (4D) computed tomography (CT) image data sets. We acquired a 4D CT image data set of a patient's thorax that consisted of three-dimensional (3D) image data sets from eight phases in the respiratory cycle. The contours of the right and left lungs, cord, heart, and esophagus were manually delineated on the end inspiration data set. An interactive deformable model algorithm, originally intended for deforming an atlas-based model surface to a 3D CT image data set, was applied in an automated fashion. Triangulations based on the contours generated on each phase were deformed to the CT data set on the succeeding phase to generate the contours on that phase. Deformation was propagated through the eight phases, and the contours obtained on the end inspiration data set were compared with the original manually delineated contours. Structures defined by high-density gradients, such as lungs, cord, and heart, were accurately reproduced, except in regions where other gradient boundaries may have confused the algorithm, such as near bronchi. The algorithm failed to accurately contour the esophagus, a soft-tissue structure completely surrounded by tissue of similar density, without manual interaction. This technique has the potential to facilitate contour delineation in 4D CT image data sets; and future evolution of the software is expected to improve the process.

Motion compensated coronary interventional navigation by means of diaphragm tracking and elastic motion models

Current catheter tracking in the x-ray catheter laboratory during coronary interventions is performed using 2D fluoroscopy. Although this features real-time navigation on high-resolution images, drawbacks such as overlap and foreshortening exist and hamper the diagnosis and treatment process. An alternative to fluoroscopy-based tracking is device tracking by means of a magnetic tracking system (MTS). Having measured the 3D location of the interventional device, its position can be reconstructed on 3D images or virtual roadmaps of the organ or vessel structure under examination. In this paper, a method is presented which compensates the interventional device location measured by the MTS for organ motion and thus registers it dynamically to a 3D virtual roadmap. The motion compensation is accomplished by using an elastic motion model which is driven by the ECG signal and a respiratory sensor signal derived from ultrasonic diaphragm tracking. The model is updated during the intervention itself, thus allowing for a local refinement in regions which bear a complex geometric structure, such as stenoses and bifurcations. The evaluation is done by means of a phantom-based study using a dynamic heart-phantom. The mean displacement caused by the overall motion of the heart is improved from 10.4+/-4.8 mm in the uncompensated case to 2.1+/-1.2 mm in the motion compensated case.

Automatic phase determination for retrospectively gated cardiac CT

The recent improvements in CT detector and gantry technology in combination with new heart rate adaptive cone beam reconstruction algorithms enable the visualization of the heart in three dimensions at high spatial resolution. However, the finite temporal resolution still impedes the artifact-free reconstruction of the heart at any arbitrary phase of the cardiac cycle. Cardiac phases must be found during which the heart is quasistationary to obtain outmost image quality. It is challenging to find these phases due to intercycle and patient-to-patient variability. Electrocardiogram (ECG) information does not always represent the heart motion with an adequate accuracy. In this publication, a simple and efficient image-based technique is introduced which is able to deliver stable cardiac phases in an automatic and patient-specific way. From low-resolution four-dimensional data sets, the most stable phases are derived by calculating the object similarity between subsequent phases in the cardiac cycle. Patient-specific information about the object motion can be determined and resolved spatially. This information is used to perform optimized high-resolution reconstructions at phases of little motion. Results based on a simulation study and three real patient data sets are presented. The projection data were generated using a 16-slice cone beam CT system in low-pitch helical mode with parallel ECG recording.

Helical cardiac cone beam CT reconstruction with large area detectors: a simulation study

Retrospectively gated cardiac volume CT imaging has become feasible with the introduction of heart rate adaptive cardiac CT reconstruction algorithms. The development in detector technology and the rapid introduction of multi-row detectors has demanded reconstruction schemes which account for the cone geometry. With the extended cardiac reconstruction (ECR) framework, the idea of approximate helical cone beam CT has been extended to be used with retrospective gating, enabling heart rate adaptive cardiac cone beam reconstruction. In this contribution, the ECR technique is evaluated for systems with an increased number of detector rows, which leads to larger cone angles. A simulation study has been carried out based on a 4D cardiac phantom consisting of a thorax model and a dynamic heart insert. Images have been reconstructed for different detector set-ups. Reconstruction assessment functions have been calculated for the detector set-ups employing different rotation times, relative pitches and heart rates. With the increased volume coverage of large area detector systems, low-pitch scans become feasible without resulting in extensive scan times, inhibiting single breath hold acquisitions. ECR delivers promising image results when being applied to systems with larger cone angles.

Evaluation of biventricular ejection fraction with ECG-gated 16-slice CT: preliminary findings in acute pulmonary embolism in comparison with radionuclide ventriculography

This study aimed to assess the feasibility of cardiac global function evaluation during a whole-chest multi-slice CT (MSCT) acquisition in patients referred for suspicion of pulmonary embolism (PE), and to compare the results with planar equilibrium radionuclide ventriculography (ERNA). Ten consecutive haemodynamically stable patients (six female, four male; mean age 69.7 years; heart rate 65-99 bpm) with suspicion of PE underwent an MSCT and ERNA within a 6 h period. CT acquisition was performed after contrast medium injection by using 16x1.5 mm collimation and retrospective ECG gating. Left ventricular (LVEF) and right ventricular (RVEF) ejection fractions were calculated using dedicated three-dimensional software. Relationships between measurements obtained with MSCT and ERNA were assessed using linear regression analysis and reliability of MSCT was assessed with intra-class correlation coefficient. Bland-Altman analysis was performed to calculate limits of agreement between MSCT and ERNA. MSCT was performed successfully in ten patients with a mean acquisition time of 16.5+/-2.8 s. Functional cardiac evaluation was possible on CT for all patients except for one due to poor opacification of right ventricle. Linear regression analysis showed a good correlation between MSCT and ERNA for the LVEF (R=0.91) and the RVEF (R=0.89) measurements. Intra-class correlation was superior for LVEF (0.92) than for the RVEF (0.68). Bland-Altman plots demonstrated that MSCT substantially overestimated the ERNA RVEF. Morphological CT data demonstrated PE in four of ten of patients and alternative diagnoses in five of ten patients. Our study reveals that MSCT with retrospective ECG gating may provide in one modality a morphological and a functional cardiopulmonary evaluation. Comparison with ERNA demonstrated a good correlation for both ventricular ejection fractions.

Evaluation of spatial and temporal resolution for ECG-gated 16-row multidetector CT using a dynamic cardiac phantom

Measurements of spatial and temporal resolution for ECG-gated scanning of a stationary and moving heart phantom with a 16-row MDCT were performed. A resolution phantom with cylindrical holes from 0.4 to 3.0 mm diameter was mounted to a cardiac phantom, which simulates the motion of a beating heart. Data acquisition was performed with 16x0.75 mm at various heart rates (HR, 60-120 bpm), pitches (0.15-0.30) and scanner rotation times (RT, 0.42 and 0.50 s). Raw data were reconstructed using a multi-cycle real cone-beam reconstruction algorithm at multiple phases of the RR interval. Multi-planar reformations (MPR) were generated and analyzed. Temporal resolution and cardiac cycles used for image reconstruction were calculated. In 97.2% (243/250) of data obtained with the stationary phantom, the complete row of holes with 0.6 mm was visible. These results were independent of heart rate, pitch, scanner rotation time and phase point of reconstruction. For the dynamic phantom, spatial resolution was determined during phases of minimal motion (116/250). In 40.5% (47/116), the resolution was 0.6 mm and in 37.1% (43/116) 0.7 mm. Temporal resolution varied between 63 and 205 ms, using 1.5-4.37 cardiac cycles for image reconstruction.

A novel platform simulating irregular motion to enhance assessment of respiration-correlated radiation therapy procedures

Respiratory motion continues to present challenges in the delivery of radiation therapy to tumors in the thorax and abdomen by causing movement of structures within those areas. Several approaches to account for this movement in the planning and delivery of treatment have been developed over the past several years. To assist in the development and assessment of various techniques for respiration‐correlated radiation therapy, a platform capable of programmable irregular longitudinal motion has been designed and fabricated to simulate intrafractional respiratory motion. A sliding platform and the base on which it was mounted were constructed from polycarbonate plastic, and a stepper motor provided platform motion. Respiratory motion data, either artificially generated on a spreadsheet or extracted from respiratory monitoring files, were converted to a format appropriate for driving the stepper motor. Various phantoms were placed on top of the platform and used in studies related to respiration‐correlated radiation therapy. Several applications of the platform were demonstrated, such as improving the quality of acquisition of time‐dependent computed tomography image datasets, comparing various methods of acquiring such datasets, and implementing feedback‐guided breath hold treatment delivery procedures. This study showed that a platform capable of programmable irregular motion is a useful tool for the development and assessment of procedures related to the effects of respiratory motion in radiation therapy.

PACS number: 87.66.Xa

Noninvasive coronary angiography with 16-detector row CT: effect of heart rate

PURPOSE

To evaluate the effect of heart rate on the quality of coronary angiograms obtained with 16-detector row computed tomography (CT) by using temporally enhanced three-dimensional (3D) approaches.

MATERIALS AND METHODS

The local ethics committee approved the study, and informed consent was obtained from all patients. Fifty patients underwent coronary CT angiography (heart rate range, 45-103 beats per minute). Raw data from helical CT and electrocardiography (ECG) were saved in a combined data set. Retrospectively ECG-gated images were reconstructed at preselected phases (50% and 80%) of the cardiac cycle. A 3D voxel-based approach with cardiac phase weighting was used for reconstruction. Testing for correlation between heart rate, cardiac phase reconstruction window, and image quality was performed with Kruskal-Wallis analysis. Image quality (freedom from cardiac motion-related artifacts) was referenced against findings at conventional angiography in a secondary evaluation step. Regression analysis was performed to calculate heart rate thresholds for future beta-blocker application.

RESULTS

A significant negative correlation was observed between heart rate and image quality (r = 0.80, P < .001). Motion artifact-free images were available for 44 (88%) patients and were achieved consistently at a heart rate of 80 or fewer beats per minute (n = 39). Best image quality was achieved at 75 or fewer beats per minute. Segmental analysis revealed that 97% of arterial segments (diameter > or = 1.5 mm according to conventional angiography) were assessable at 80 or fewer beats per minute. Premature ventricular contractions and rate-contained arrhythmia did not impede diagnostic assessment of the coronary arteries in 10 (83%) of the 12 patients affected.

CONCLUSION

Motion-free coronary angiograms can be obtained consistently with 16-detector row CT scanners and adaptive multicyclic reconstruction algorithms in patients with heart rates of less than 80 beats per minute.

Artifact analysis and reconstruction improvement in helical cardiac cone beam CT

With the introduction of cone beam (CB) scanners, cardiac volumetric computed tomography (CT) imaging has the potential to become a noninvasive imaging tool in clinical routine for the diagnosis of various heart diseases. Heart rate adaptive reconstruction schemes enable the reconstruction of high-resolution volumetric data sets of the heart. Artifacts, caused by strong heart rate variations, high heart rates and obesity, decrease the image quality and the diagnostic value of the images. The image quality suffers from streak artifacts if suboptimal scan and reconstruction parameters are chosen, demanding improved gating techniques. In this paper, an artifact analysis is carried out which addresses the artifacts due to the gating when using a three-dimensional CB cardiac reconstruction technique. An automatic and patient specific cardiac weighting technique is presented in order to improve the image quality. Based on the properties of the reconstruction algorithm, several assessment techniques are introduced which enable the quantitative determination of the cycle-to-cycle transition smoothness and phase homogeneity of the image reconstruction. Projection data of four patients were acquired using a 16-slice CBCT system in low pitch helical mode with parallel electrocardiogram recording. For each patient, image results are presented and discussed in combination with the assessment criteria.

Extended parallel backprojection for standard three-dimensional and phase-correlated four-dimensional axial and spiral cone-beam CT with arbitrary pitch, arbitrary cone-angle, and 100% dose usage

We have developed a new approximate Feldkamp-type algorithm that we call the extended parallel backprojection (EPBP). Its main features are a phase-weighted backprojection and a voxel-by-voxel 180 degrees normalization. The first feature ensures three-dimensional (3-D) and 4-D capabilities with one and the same algorithm; the second ensures 100% detector usage (each ray is accounted for). The algorithm was evaluated using simulated data of a thorax phantom and a cardiac motion phantom for scanners with up to 256 slices. Axial (circle and sequence) and spiral scan trajectories were investigated. The standard reconstructions (EPBPStd) are of high quality, even for as many as 256 slices. The cardiac reconstructions (EPBPCI) are of high quality as well and show no significant deterioration of objects even far off the center of rotation. Since EPBPCI uses the cardio interpolation (CI) phase weighting the temporal resolution is equivalent to that of the well-established single-slice and multislice cardiac approaches 180 degrees CI, 180 degrees MCI, and ASSRCI, respectively, and lies in the order of 50 to 100 ms for rotation times between 0.4 and 0.5 s. EPBP appears to fulfill all required demands. Especially the phase-correlated EPBP reconstruction of cardiac multiple circle scan data is of high interest, e.g., for dynamic perfusion studies of the heart.

Acquiring 4D thoracic CT scans using a multislice helical method

Respiratory motion degrades anatomic position reproducibility during imaging, necessitates larger margins during radiotherapy planning and causes errors during radiation delivery. Computed tomography (CT) scans acquired synchronously with the respiratory signal can be used to reconstruct 4D CT scans, which can be employed for 4D treatment planning to explicitly account for respiratory motion. The aim of this research was to develop, test and clinically implement a method to acquire 4D thoracic CT scans using a multislice helical method. A commercial position-monitoring system used for respiratory-gated radiotherapy was interfaced with a third generation multislice scanner. 4D cardiac reconstruction methods were modified to allow 4D thoracic CT acquisition. The technique was tested on a phantom under different conditions: stationary, periodic motion and non-periodic motion. 4D CT was also implemented for a lung cancer patient with audio-visual breathing coaching. For all cases, 4D CT images were successfully acquired from eight discrete breathing phases, however, some limitations of the system in terms of respiration reproducibility and breathing period relative to scanner settings were evident. Lung mass for the 4D CT patient scan was reproducible to within 2.1% over the eight phases, though the lung volume changed by 20% between end inspiration and end expiration (870 cm3). 4D CT can be used for 4D radiotherapy, respiration-gated radiotherapy, 'slow' CT acquisition and tumour motion studies.