Iterative concurrent reconstruction algorithms for emission computed tomography

Direct reconstruction techniques, such as those based on filtered backprojection, are typically used for emission computed tomography (ECT), even though it has been argued that iterative reconstruction methods may produce better clinical images. The major disadvantage of iterative reconstruction algorithms, and a significant reason for their lack of clinical acceptance, is their computational burden. We outline a new class of 'concurrent' iterative reconstruction techniques for ECT in which the reconstruction process is reorganized such that a significant fraction of the computational processing occurs concurrently with the acquisition of ECT projection data. These new algorithms use the 10-30 min required for acquisition of a typical SPECT scan to iteratively process the available projection data, significantly reducing the requirements for post-acquisition processing. These algorithms are tested on SPECT projection data from a Hoffman brain phantom acquired with a 2 x 10(5) counts in 64 views each having 64 projections. The SPECT images are reconstructed as 64 x 64 tomograms, starting with six angular views. Other angular views are added to the reconstruction process sequentially, in a manner that reflects their availability for a typical acquisition protocol. The results suggest that if T s of concurrent processing are used, the reconstruction processing time required after completion of the data acquisition can be reduced by at least 1/3T s.     

Investigation of dynamic SPECT measurements of the arterial input function in human subjects using simulation, phantom and human studies

Computer simulations, a phantom study and a human study were performed to determine whether a slowly rotating single-photon computed emission tomography (SPECT) system could provide accurate arterial input functions for quantification of myocardial perfusion imaging using kinetic models. The errors induced by data inconsistency associated with imaging with slow camera rotation during tracer injection were evaluated with an approach called SPECT/P (dynamic SPECT from positron emission tomography (PET)) and SPECT/D (dynamic SPECT from database of SPECT phantom projections). SPECT/P simulated SPECT-like dynamic projections using reprojections of reconstructed dynamic (94)Tc-methoxyisobutylisonitrile ((94)Tc-MIBI) PET images acquired in three human subjects (1 min infusion). This approach was used to evaluate the accuracy of estimating myocardial wash-in rate parameters K(1) for rotation speeds providing 180° of projection data every 27 or 54 s. Blood input and myocardium tissue time-activity curves (TACs) were estimated using spatiotemporal splines. These were fit to a one-compartment perfusion model to obtain wash-in rate parameters K(1). For the second method (SPECT/D), an anthropomorphic cardiac torso phantom was used to create real SPECT dynamic projection data of a tracer distribution derived from (94)Tc-MIBI PET scans in the blood pool, myocardium, liver and background. This method introduced attenuation, collimation and scatter into the modeling of dynamic SPECT projections. Both approaches were used to evaluate the accuracy of estimating myocardial wash-in parameters for rotation speeds providing 180° of projection data every 27 and 54 s. Dynamic cardiac SPECT was also performed in a human subject at rest using a hybrid SPECT/CT scanner. Dynamic measurements of (99m)Tc-tetrofosmin in the myocardium were obtained using an infusion time of 2 min. Blood input, myocardium tissue and liver TACs were estimated using the same spatiotemporal splines. The spatiotemporal maximum-likelihood expectation-maximization (4D ML-EM) reconstructions gave more accurate reconstructions than did standard frame-by-frame static 3D ML-EM reconstructions. The SPECT/P results showed that 4D ML-EM reconstruction gave higher and more accurate estimates of K(1) than did 3D ML-EM, yielding anywhere from a 44% underestimation to 24% overestimation for the three patients. The SPECT/D results showed that 4D ML-EM reconstruction gave an overestimation of 28% and 3D ML-EM gave an underestimation of 1% for K(1). For the patient study the 4D ML-EM reconstruction provided continuous images as a function of time of the concentration in both ventricular cavities and myocardium during the 2 min infusion. It is demonstrated that a 2 min infusion with a two-headed SPECT system rotating 180° every 54 s can produce measurements of blood pool and myocardial TACs, though the SPECT simulation studies showed that one must sample at least every 30 s to capture a 1 min infusion input function.     

Exact inversion of the exponential x-ray transform for rotating slant-hole (RSH) SPECT

The RSH SPECT scanner provides parallel-beam attenuated projections for a fully 3D acquisition geometry. The geometry can be represented by circles on the unit sphere of projection directions, one circle for each position of the detector head. Unlike most other fully 3D geometries this one is particularly challenging because there are no 2D subsets in the data. When no attenuation is present, it is well known that an unmeasured projection can be synthesized if it lies inside one of the measured circles. The main result of this work is that under some assumptions on the attenuation distribution, attenuated projections within a circle can be synthesized from available attenuated projections. One consequence is that RSH SPECT projections can be rebinned into a conventional SPECT geometry for which analytic attenuation correction techniques are available.     

Fast generation of 4D PET-MR data from real dynamic MR acquisitions

We have implemented and evaluated a framework for simulating simultaneous dynamic PET-MR data using the anatomic and dynamic information from real MR acquisitions. PET radiotracer distribution is simulated by assigning typical FDG uptake values to segmented MR images with manually inserted additional virtual lesions. PET projection data and images are simulated using analytic forward projections (including attenuation and Poisson statistics) implemented within the image reconstruction package STIR. PET image reconstructions are also performed with STIR. The simulation is validated with numerical simulation based on Monte Carlo (GATE) which uses more accurate physical modelling, but has 150× slower computation time compared to the analytic method for ten respiratory positions and is 7000× slower when performing multiple realizations. Results are validated in terms of region of interest mean values and coefficients of variation for 65 million coincidences including scattered events. Although some discrepancy is observed, agreement between the two different simulation methods is good given the statistical noise in the data. In particular, the percentage difference of the mean values is 3.1% for tissue, 17% for the lungs and 18% for a small lesion. The utility of the procedure is demonstrated by simulating realistic PET-MR datasets from multiple volunteers with different breathing patterns. The usefulness of the toolkit will be shown for performance investigations of the reconstruction, motion correction and attenuation correction algorithms for dynamic PET-MR data.     

Cone-beam CT with megavoltage beams and an amorphous silicon electronic portal imaging device: potential for verification of radiotherapy of lung cancer

We investigate the potential of megavoltage (MV) cone-beam CT with an amorphous silicon electronic portal imaging device (EPID) as a tool for patient position verification and tumor/organ motion studies in radiation treatment of lung tumors. We acquire 25 to 200 projection images using a 22 x 29 cm EPID. The acquisition is automatic and requires 7 minutes for 100 projections; it can be synchronized with respiratory gating. From these images, volumetric reconstruction is accomplished with a filtered backprojection in the cone-beam geometry. Several important prereconstruction image corrections, such as detector sag, must be applied. Tests with a contrast phantom indicate that differences in electron density of 2% can be detected with 100 projections, 200 cGy total dose. The contrast-to-noise ratio improves as the number of projections is increased. With 50 projections (100 cGy), high contrast objects are visible, and as few as 25 projections yield images with discernible features. We identify a technique of acquiring projection images with conformal beam apertures, shaped by a multileaf collimator, to reduce the dose to surrounding normal tissue. Tests of this technique on an anthropomorphic phantom demonstrate that a gross tumor volume in the lung can be accurately localized in three dimensions with scans using 88 monitor units. As such, conformal megavoltage cone-beam CT can provide three-dimensional imaging of lung tumors and may be used, for example, in verifying respiratory gated treatments.     

A technique for on-board CT reconstruction using both kilovoltage and megavoltage beam projections for 3D treatment verification

The technologies with kilovoltage (kV) and megavoltage (MV) imaging in the treatment room are now available for image-guided radiation therapy to improve patient setup and target localization accuracy. However, development of strategies to efficiently and effectively implement these technologies for patient treatment remains challenging. This study proposed an aggregated technique for on-board CT reconstruction using combination of kV and MV beam projections to improve the data acquisition efficiency and image quality. These projections were acquired in the treatment room at the patient treatment position with a new kV imaging device installed on the accelerator gantry, orthogonal to the existing MV portal imaging device. The projection images for a head phantom and a contrast phantom were acquired using both the On-Board Imager kV imaging device and the MV portal imager mounted orthogonally on the gantry of a Varian Clinac 21EX linear accelerator. MV projections were converted into kV information prior to the aggregated CT reconstruction. The multilevel scheme algebraic-reconstruction technique was used to reconstruct CT images involving either full, truncated, or a combination of both full and truncated projections. An adaptive reconstruction method was also applied, based on the limited numbers of kV projections and truncated MV projections, to enhance the anatomical information around the treatment volume and to minimize the radiation dose. The effects of the total number of projections, the combination of kV and MV projections, and the beam truncation of MV projections on the details of reconstructed kV/MV CT images were also investigated.     

Angle-independent measure of motion for image-based gating in 3D coronary angiography

The role of three-dimensional (3D) image guidance for interventional procedures and minimally invasive surgeries is increasing for the treatment of vascular disease. Currently, most interventional procedures are guided by two-dimensional x-ray angiography, but computed rotational angiography has the potential to provide 3D geometric information about the coronary arteries. The creation of 3D angiographic images of the coronary arteries requires synchronization of data acquisition with respect to the cardiac cycle, in order to minimize motion artifacts. This can be achieved by inferring the extent of motion from a patient's electrocardiogram (ECG) signal. However, a direct measurement of motion (from the 2D angiograms) has the potential to improve the 3D angiographic images by ensuring that only projections acquired during periods of minimal motion are included in the reconstruction. This paper presents an image-based metric for measuring the extent of motion in 2D x-ray angiographic images. Adaptive histogram equalization was applied to projection images to increase the sharpness of coronary arteries and the superior-inferior component of the weighted centroid (SIC) was measured. The SIC constitutes an image-based metric that can be used to track vessel motion, independent of apparent motion induced by the rotational acquisition. To evaluate the technique, six consecutive patients scheduled for routine coronary angiography procedures were studied. We compared the end of the SIC rest period (rho) to R-waves (R) detected in the patient's ECG and found a mean difference of 14 +/- 80 ms. Two simultaneous angular positions were acquired and rho was detected for each position. There was no statistically significant difference (P = 0.79) between rho in the two simultaneously acquired angular positions. Thus we have shown the SIC to be independent of view angle, which is critical for rotational angiography. A preliminary image-based gating strategy that employed the SIC was compared to an ECG-based gating strategy in a porcine model. The image-based gating strategy selected 61 projection images, compared to 45 selected by the ECG-gating strategy. Qualitative comparison revealed that although both the SIC-based and ECG-gated reconstructions decreased motion artifact compared to reconstruction with no gating, the SIC-based gating technique increased the conspicuity of smaller vessels when compared to ECG gating in maximum intensity projections of the reconstructions and increased the sharpness of a vessel cross section in multi-planar reformats of the reconstruction.     

A pseudo-Poisson noise model for simulation of positron emission tomographic projection data.

Although radioactive decay obeys Poisson statistics, because of the corrections that are applied to the projection data prior to reconstruction, the noise in positron emission tomography (PET) projections does not follow a Poisson distribution. Use of Poisson noise when simulating PET projections in order to test the performance of reconstruction and processing techniques is therefore not appropriate. The magnitude of PET projection noise was observed to be as much as 10 to 100 times greater than Poisson noise in some instances. A quadratic function was found to fit the relationship between noise power spectral density and total projection count. The coefficients of the quadratic function were determined for projections of different tracer distributions and types. Using these observations, a method of simulating PET projections was developed based on a pseudo-Poisson noise model. Projections simulated according to this method are good approximations to real projection data and take into account the characteristics of individual PET cameras and particular tracer distributions. Such simulated projections have been valuable in predicting the performance of reconstruction algorithms. This approach can also be applied to single photon emission tomography.     

Three-dimensional SPECT reconstruction of combined cone beam and parallel beam data.

Single photon emission computed tomography (SPECT) using cone beam (CB) collimation exhibits increased sensitivity compared with acquisition geometries using parallel (P) hole collimation. However, CB collimation has a smaller field-of-view which may result in truncated projections and image artifacts. A primary objective of this work is to investigate maximum likelihood-expectation maximization (ML-EM) methods to reconstruct simultaneously acquired parallel and cone beam (P&CB) SPECT data. Simultaneous P&CB acquisition can be performed with commercially available triple camera systems by using two cone-beam collimators and a single parallel-hole collimator. The loss in overall sensitivity (relative to the use of three CB collimators) is about 15 to 20%. We have developed three methods to combine P&CB data using modified ML-EM algorithms. The first method consists of using both data sets to reconstruct a single intermediate image after each iteration using the ML-EM algorithm. The other two iterative algorithms combine intermediate parallel beam (PB) and CB source estimates to enhance image quality. For these methods, a PB estimate and a CB estimate are obtained for the first iteration. The second method consists of summing the PB and CB estimates for each subsequent iteration to obtain new PB and CB estimates. The third method is similar to the second method, with the exception that the new PB estimate simply is set equal to the PB estimate after each iteration. The combined source estimate is used in each subsequent iteration step of the EM algorithm. These algorithms are evaluated using projection data simulated using a Monte Carlo SPECT model. The P&CB SPECT images demonstrate marked improvement as compared with the CB-only reconstruction, particularly when the projections are truncated.     

Performance studies of four-dimensional cone beam computed tomography

Four-dimensional cone beam computed tomography (4DCBCT) has been proposed to characterize the breathing motion of tumors before radiotherapy treatment. However, when the acquired cone beam projection data are retrospectively gated into several respiratory phases, the available data to reconstruct each phase is under-sampled and thus causes streaking artifacts in the reconstructed images. To solve the under-sampling problem and improve image quality in 4DCBCT, various methods have been developed. This paper presents performance studies of three different 4DCBCT methods based on different reconstruction algorithms. The aims of this paper are to study (1) the relationship between the accuracy of the extracted motion trajectories and the data acquisition time of a 4DCBCT scan and (2) the relationship between the accuracy of the extracted motion trajectories and the number of phase bins used to sort projection data. These aims will be applied to three different 4DCBCT methods: conventional filtered backprojection reconstruction (FBP), FBP with McKinnon-Bates correction (MB) and prior image constrained compressed sensing (PICCS) reconstruction. A hybrid phantom consisting of realistic chest anatomy and a moving elliptical object with known 3D motion trajectories was constructed by superimposing the analytical projection data of the moving object to the simulated projection data from a chest CT volume dataset. CBCT scans with gantry rotation times from 1 to 4 min were simulated, and the generated projection data were sorted into 5, 10 and 20 phase bins before different methods were used to reconstruct 4D images. The motion trajectories of the moving object were extracted using a fast free-form deformable registration algorithm. The root mean square errors (RMSE) of the extracted motion trajectories were evaluated for all simulated cases to quantitatively study the performance. The results demonstrate (1) longer acquisition times result in more accurate motion delineation for each method; (2) ten or more phase bins are necessary in 4DCBCT to ensure sufficient temporal resolution in tumor motion and (3) to achieve the same performance as FBP-4DCBCT with a 4 min data acquisition time, MB-4DCBCT and PICCS-4DCBCT need about 2- and 1 min data acquisition times, respectively.     

Analytic image reconstruction in local phase-contrast tomography

Phase-contrast tomography is a non-interferometric imaging technique for reconstructing the refractive index distribution of a weakly absorbing object from a set of tomographic projection measurements. In many practical situations, the spatial resolution of the reconstructed image can be increased by minimizing the field of view (FOV) of the imaging system. When the object of interest is larger than the FOV, the measured projections are truncated and one is faced with a local tomography reconstruction problem. In this work, we analytically and numerically investigate the problem of reconstructing tomographic images from truncated phase-contrast projection data. A simple backprojection algorithm for reconstructing object discontinuities from truncated phase-contrast projection data is proposed and investigated that involves no explicit filtering of the projection data. We also investigate the use of the filtered backprojection algorithm and a local tomography reconstruction algorithm developed for absorption CT. These reconstruction algorithms are implemented and numerically investigated to corroborate our theoretical assertions.     

Four-dimensional cone-beam computed tomography using an on-board imager

On-board cone-beam computed tomography (CBCT) has recently become available to provide volumetric information of a patient in the treatment position, and holds promises for improved target localization and irradiation dose verification. The design of currently available on-board CBCT, however, is far from optimal. Its quality is adversely influenced by many factors, such as scatter, beam hardening, and intra-scanning organ motion. In this work we quantitatively study the influence of organ motion on CBCT imaging and investigate a strategy to acquire high quality phase-resolved [four-dimensional (4D)] CBCT images based on phase binning of the CBCT projection data. An efficient and robust method for binning CBCT data according to the patient's respiratory phase derived in the projection space was developed. The phase-binned projections were reconstructed using the conventional Feldkamp algorithm to yield 4D CBCT images. Both phantom and patient studies were carried out to validate the technique and to optimize the 4D CBCT data acquisition protocol. Several factors that are important to the clinical implementation of the technique, such as the image quality, scanning time, number of projections, and radiation dose, were analyzed for various scanning schemes. The general references drawn from this study are: (i) reliable phase binning of CBCT projections is accomplishable with the aid of external or internal marker and simple analysis of its trace in the projection space, and (ii) artifact-free 4D CBCT images can be obtained without increasing the patient radiation dose as compared to the current 3D CBCT scan.     

A novel digital tomosynthesis (DTS) reconstruction method using a deformation field map

We developed a novel digital tomosynthesis (DTS) reconstruction method using a deformation field map to optimally estimate volumetric information in DTS images. The deformation field map is solved by using prior information, a deformation model, and new projection data. Patients' previous cone-beam CT (CBCT) or planning CT data are used as the prior information, and the new patient volume to be reconstructed is considered as a deformation of the prior patient volume. The deformation field is solved by minimizing bending energy and maintaining new projection data fidelity using a nonlinear conjugate gradient method. The new patient DTS volume is then obtained by deforming the prior patient CBCT or CT volume according to the solution to the deformation field. This method is novel because it is the first method to combine deformable registration with limited angle image reconstruction. The method was tested in 2D cases using simulated projections of a Shepp-Logan phantom, liver, and head-and-neck patient data. The accuracy of the reconstruction was evaluated by comparing both organ volume and pixel value differences between DTS and CBCT images. In the Shepp-Logan phantom study, the reconstructed pixel signal-to-noise ratio (PSNR) for the 60 degrees DTS image reached 34.3 dB. In the liver patient study, the relative error of the liver volume reconstructed using 60 degrees projections was 3.4%. The reconstructed PSNR for the 60 degrees DTS image reached 23.5 dB. In the head-and-neck patient study, the new method using 60 degrees projections was able to reconstruct the 8.1 degrees rotation of the bony structure with 0.0 degrees error. The reconstructed PSNR for the 60 degrees DTS image reached 24.2 dB. In summary, the new reconstruction method can optimally estimate the volumetric information in DTS images using 60 degrees projections. Preliminary validation of the algorithm showed that it is both technically and clinically feasible for image guidance in radiation therapy.     

A preliminary investigation of local tomography for megavoltage CT imaging

We investigate the problem of reconstructing a two-dimensional (2-D) cross-sectional image of a tumor volume from a set of truncated MV projections that are produced by radiation therapy treatment beams. Our proposed approach is conceptually distinct from previously investigated approaches in that it utilizes a noniterative local tomography reconstruction algorithm. A local tomography reconstruction algorithm is implemented and systematically investigated using several sets of simulated and experimental MV projection data. We demonstrate that the conventional (nonlocal) filtered backprojection reconstruction algorithm cannot, in general, accurately reconstruct the edges and boundaries of low-contrast features from truncated MV projection data. We demonstrate that the local tomography algorithm is not adversely affected by projection truncation and can reconstruct accurately the boundaries of low-contrast structures within the region of interest from truncated MV projections.     

Projection-based motion compensation for gated coronary artery reconstruction from rotational x-ray angiograms

Three-dimensional reconstruction of coronary arteries can be performed during x-ray-guided interventions by gated reconstruction from a rotational coronary angiography sequence. Due to imperfect gating and cardiac or breathing motion, the heart's motion state might not be the same in all projections used for the reconstruction of one cardiac phase. The motion state inconsistency causes motion artefacts and degrades the reconstruction quality. These effects can be reduced by a projection-based 2D motion compensation method. Using maximum-intensity forward projections of an initial uncompensated reconstruction as reference, the projection data are transformed elastically to improve the consistency with respect to the heart's motion state. A fast iterative closest-point algorithm working on vessel centrelines is employed for estimating the optimum transformation. Motion compensation is carried out prior to and independently from a final reconstruction. The motion compensation improves the accuracy of reconstructed vessel radii and the image contrast in a software phantom study. Reconstructions of human clinical cases are presented, in which the motion compensation substantially reduces motion blur and improves contrast and visibility of the coronary arteries.     

Characterization of scattered radiation in kV CBCT images using Monte Carlo simulations

Kilovoltage (kV) cone beam computed tomography (CBCT) images suffer from a substantial scatter contribution. In this study, Monte Carlo (MC) simulations are used to evaluate the scattered radiation present in projection images. These predicted scatter distributions are also used as a scatter correction technique. Images were acquired using a kV CBCT bench top system. The EGSnrc MC code was used to model the flat panel imager, the phantoms, and the x-ray source. The x-ray source model was validated using first and second half-value layers (HVL) and profile measurements. The HVLs and the profile were found to agree within 3% and 6%, respectively. MC simulated and measured projection images for a cylindrical water phantom and for an anthropomorphic head phantom agreed within 8% and 10%. A modified version of the DOSXYZnrc MC code was used to score phase space files with identified scattered and primary particles behind the phantoms. The cone angle, the source-to-detector distance, the phantom geometry, and the energy were varied to determine their effect on the scattered radiation distribution. A scatter correction technique was developed in which the MC predicted scatter distribution is subtracted from the projections prior to reconstruction. Preliminary testing of the procedure was done with an anthropomorphic head phantom and a contrast phantom. Contrast and profile measurements were obtained for the scatter corrected and noncorrected images. An improvement of 3% for contrast between solid water and a liver insert and 11% between solid water and a Teflon insert were obtained and a significant reduction in cupping and streaking artifacts was observed.     

Exact reconstruction in 2D dynamic CT: compensation of time-dependent affine deformations

This work is dedicated to the reduction of reconstruction artefacts due to motion occurring during the acquisition of computerized tomographic projections. This problem has to be solved when imaging moving organs such as the lungs or the heart. The proposed method belongs to the class of motion compensation algorithms, where the model of motion is included in the reconstruction formula. We address two fundamental questions. First what conditions on the deformation are required for the reconstruction of the object from projections acquired sequentially during the deformation, and second how do we reconstruct the object from those projections. Here we answer these questions in the particular case of 2D general time-dependent affine deformations, assuming the motion parameters are known. We treat the problem of admissibility conditions on the deformation in the parallel-beam and fan-beam cases. Then we propose exact reconstruction methods based on rebinning or sequential FBP formulae for each of these geometries and present reconstructed images obtained with the fan-beam algorithm on simulated data.     

MR-based keyhole SPECT for small animal imaging

The rationale for multi-modality imaging is to integrate the strengths of different imaging technologies while reducing the shortcomings of an individual modality. The work presented here proposes a limited-field-of-view (LFOV) SPECT reconstruction technique that can be implemented on a multi-modality MR/SPECT system that can be used to obtain simultaneous MRI and SPECT images for small animal imaging. The reason for using a combined MR/SPECT system in this work is to eliminate any possible misregistration between the two sets of images when MR images are used as a priori information for SPECT. In nuclear imaging the target area is usually smaller than the entire object; thus, focusing the detector on the LFOV results in various advantages including the use of a smaller nuclear detector (less cost), smaller reconstruction region (faster reconstruction) and higher spatial resolution when used in conjunction with pinhole collimators with magnification. The MR/SPECT system can be used to choose a region of interest (ROI) for SPECT. A priori information obtained by the full field-of-view (FOV) MRI combined with the preliminary SPECT image can be used to reduce the dimensions of the SPECT reconstruction by limiting the computation to the smaller FOV while reducing artifacts resulting from the truncated data. Since the technique is based on SPECT imaging within the LFOV it will be called the keyhole SPECT (K-SPECT) method. At first MRI images of the entire object using a larger FOV are obtained to determine the location of the ROI covering the target organ. Once the ROI is determined, the animal is moved inside the radiofrequency (rf) coil to bring the target area inside the LFOV and then simultaneous MRI and SPECT are performed. The spatial resolution of the SPECT image is improved by employing a pinhole collimator with magnification >1 by having carefully calculated acceptance angles for each pinhole to avoid multiplexing. In our design all the pinholes are focused to the center of the LFOV. K-SPECT reconstruction is accomplished by generating an adaptive weighting matrix using a priori information obtained by simultaneously acquired MR images and the radioactivity distribution obtained from the ROI region of the SPECT image that is reconstructed without any a priori input. Preliminary results using simulations with numerical phantoms show that the image resolution of the SPECT image within the LFOV is improved while minimizing artifacts arising from parts of the object outside the LFOV due to the chosen magnification and the new reconstruction technique. The root-mean-square-error (RMSE) in the out-of-field artifacts was reduced by 60% for spherical phantoms using the K-SPECT reconstruction technique and by 48.5-52.6% for the heart in the case with the MOBY phantom. The K-SPECT reconstruction technique significantly improved the spatial resolution and quantification while reducing artifacts from the contributions outside the LFOV as well as reducing the dimension of the reconstruction matrix.     

A motion-compensated scheme for helical cone-beam reconstruction in cardiac CT angiography

Since coronary heart disease is one of the main causes of death all over the world, cardiac computed tomography (CT) imaging is an application of very high interest in order to verify indications timely. Due to the cardiac motion, electrocardiogram (ECG) gating has to be implemented into the reconstruction of the measured projection data. However, the temporal and spatial resolution is limited due to the mechanical movement of the gantry and due to the fact that a finite angular span of projections has to be acquired for the reconstruction of each voxel. In this article, a motion-compensated reconstruction method for cardiac CT is described, which can be used to increase the signal-to-noise ratio or to suppress motion blurring. Alternatively, it can be translated into an improvement of the temporal and spatial resolution. It can be applied to the entire heart in common and to high contrast objects moving with the heart in particular, such as calcified plaques or devices like stents. The method is based on three subsequent steps: As a first step, the projection data acquired in low pitch helical acquisition mode together with the ECG are reconstructed at multiple phase points. As a second step, the motion-vector field is calculated from the reconstructed images in relation to the image in a reference phase. Finally, a motion-compensated reconstruction is carried out for the reference phase using those projections, which cover the cardiac phases for which the motion-vector field has been determined.     

On-board four-dimensional digital tomosynthesis: first experimental results

The purpose of this study is to propose four-dimensional digital tomosynthesis (4D-DTS) for on-board analysis of motion information in three dimensions. Images of a dynamic motion phantom were reconstructed using acquisition scan angles ranging from 20 degrees (DTS) to full 360 degrees cone-beam computed tomography (CBCT). Projection images were acquired using an on-board imager mounted on a clinical linear accelerator. Three-dimensional (3D) images of the moving target were reconstructed for various scan angles. 3D respiratory correlated phase images were also reconstructed. For phase-based image reconstructions, the trajectory of a radiopaque marker was tracked in projection space and used to retrospectively assign respiratory phases to projections. The projections were then sorted according phase and used to reconstruct motion correlated images. By using two sets of projections centered about anterior-posterior and lateral axes, this study demonstrates how phase resolved coronal and sagittal DTS images can be used to obtain 3D motion information. Motion artifacts in 4D-DTS phase images are compared with those present in four-dimensional CT (4DCT) images. Due to the nature of data acquisition for the two modalities, superior-inferior motion artifacts are suppressed to a greater extent in 4D-DTS images compared with 4DCT. Theoretical derivations and experimental results are presented to demonstrate how optimal selection of image acquisition parameters including the frequency of projection acquisition and the phase window depend on the respiratory period. Two methods for acquiring projections are discussed. Preliminary results indicate that 4D-DTS can be used to acquire valuable kinetic information of internal anatomy just prior to radiation treatment.