Image reconstruction and performance evaluation for ECG-gated spiral scanning with a 16-slice CT system

We present an image reconstruction approach and a performance evaluation for ECG-gate cardiac spiral scanning with recently introduced 16-slice CT equipment. We present an extension of the Adaptive Cardio Volume (ACV) reconstruction approach for ECG-gated multislice spiral scanning. We discuss the image z reformation introduced to control the spiral slice width of the final images and give an overview of the reformation functions chosen. We investigate image quality and discuss the maximum number of slices that can be reconstructed without severe cone-beam artifacts. Slice sensitivity profiles (SSPs) and transverse resolution are evaluated as a function of the patient's heart rate. We demonstrate the influence of slice width on the visualization of stents and plaques and show the impact of reduced gantry rotation time (0.42 s) on temporal resolution. Deviating from general purpose spiral scanning cone-beam reconstruction is not required for ECG-gated cardiac CT with up to 16 slices. Using the ACV approach with image reformation, SSPs are well defined and independent of the patient's heart rate. With 0.75 mm collimated slice width, the measured full width at half-maximum (FWHM) of the smallest reconstructed slice is about 0.83 mm. Using this slice width and overlapping image reconstruction, cylindrical holes 0.6-0.7 mm in diameter can be resolved in a z-resolution phantom. Adequate visualization of the coronary arteries requires reconstruction slice widths not larger than 1.5 mm. Visualization of stents and severe calcifications is significantly improved with sub-mm slice width. Experimental evidence for the theoretically predicted temporal resolution and for the variation of temporal resolution depending on the position in the field of measurement (FOM) is presented. With 0.42 s gantry rotation temporal resolution reaches its optimum of 105 ms in the center of the FOM at 81 bpm. First scans on human subjects demonstrate the potential to expand the range of heart rates accessible to routine clinical examinations. A 16-slice platform can cover the heart with sub-mm slices within short breath-hold times, allowing for improved cardiac imaging due to isotropic sub-mm spatial resolution.     

ECG-correlated image reconstruction from subsecond multi-slice spiral CT scans of the heart

Subsecond spiral computed tomography (CT) offers great potential for improving heart imaging. The new multi-row detector technology adds significantly to this potential. We therefore developed and validated dedicated cardiac reconstruction algorithms for imaging the heart with subsecond multi-slice spiral CT utilizing electrocardiogram (ECG) information. The single-slice cardiac z-interpolation algorithms 180 degrees CI and 180 degrees CD [Med. Phys. 25, 2417-2431 (1998)] were generalized to allow imaging of the heart for M-slice scanners. Two classes of algorithms were investigated: 180 degrees MCD (multi-slice cardio delta), a partial scan reconstruction of 180 degrees + delta data with a < phi (fan angle) resulting in effective scan times of 250 ms (central ray) when a 0.5 s rotation mode is available, and 180 degrees MCI (multi-slice cardio interpolation), a piecewise weighted interpolation between successive spiral data segments belonging to the same heart phase, potentially providing a relative temporal resolution of 12.5% of the heart cycle when a four-slice scanner is used and the table increment is chosen to be greater than or equal to the collimated slice thickness. Data segments are selected by correlation with the simultaneously recorded ECG signal. Theoretical studies, computer simulations, as well as patient measurements were carried out for a multi-slice scanner providing M = 4 slices to evaluate these new approaches and determine the optimal scan protocol. Both algorithms, 180 degrees MCD and 180 degrees MCI, provide significant improvements in image quality, including extremely arythmic cases. Artifacts in the reconstructed images as well as in 3D displays such as multiplanar reformations were largely reduced as compared to the standard z-interpolation algorithm 180 degrees MLI (multi-slice linear interpolation). Image quality appears adequate for precise calcium scoring and CT angiography of the coronary arteries with conventional subsecond multislice spiral CT. It turned out that for heart rates fH > or = 70 min(-1) the partial scan approach 180 degrees MCD yields unsatisfactory results as compared to 180 degrees MCI. Our theoretical considerations show that a freely selectable scanner rotation time chosen as a function of the patient's heart rate, would further improve the relative temporal resolution and thus further reduce motion artifacts. In our case an additional 0.6 s mode besides the available 0.5 s mode would be very helpful. Moreover, if technically feasible, lower rotation times such as 0.3 s or even less would result in improved image quality. The use of multi-slice techniques for cardiac CT together with the new z-interpolation methods improves the quality of heart imaging significantly. The high temporal resolution of 180 degrees MCI is adequate for spatial and temporal tracking of anatomic structures of the heart (4D reconstruction).     

High temporal resolution for multislice helical computed tomography

Multislice helical computed tomography (CT) substantially reduces scanning time. However, the temporal resolution of individual images is still insufficient for imaging rapidly moving organs such as the heart and adjacent pulmonary vessels. It may, in some cases, be worse than with current single-slice helical CT. The purpose of this study is to describe a novel image reconstruction algorithm to improve temporal resolution in multislice helical CT, and to evaluate its performance against existing algorithms. The proposed image reconstruction algorithm uses helical interpolation followed by data weighting based on the acquisition time. The temporal resolution, the longitudinal (z-axis) spatial resolution, the image noise, and the in-plane image artifacts created by a moving phantom were compared with those from the basic multislice helical reconstruction (helical filter interpolation, HFI) algorithm and the basic single-slice helical reconstruction algorithm (180 degrees linear interpolation, 180LI) using computer simulations. Computer simulation results were verified with CT examinations of the heart and lung vasculature using a 0.5 second multislice scanner. The temporal resolution of HFI algorithm varies from 0.28 and 0.86 s, depending on helical pitch. The proposed method improves the resolution to a constant value of 0.29 s, independent of pitch, allowing moving objects to be imaged with reduced blurring or motion artifacts. The spatial (z) resolution was slightly worse than with the HFI algorithm; the image noise was worse than with the HFI algorithm but was comparable to axial (step-and-shoot) CT. The proposed method provided sharp images of the moving objects, portraying the anatomy accurately. The proposed algorithm for multislice helical CT allowed us to obtain CT images with high temporal resolution. It may improve the image quality of clinical cardiac, lung, and vascular CT imaging.     

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.     

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.     

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.     

A technical solution to avoid partial scan artifacts in cardiac MDCT

Quantitative evaluation of cardiac image data obtained using multidetector row computed tomography (CT) is compromised by partial scan reconstructions, which improve the temporal resolution but significantly increase image-to-image CT number variations for a fixed region of interest compared to full reconstruction images. The feasibility of a new approach to solve this problem is assessed. An anthropomorphic cardiac phantom and an anesthetized pig were scanned on a dual-source CT scanner using both full and partial scan acquisition modes under different conditions. Additional scans were conducted with the electrocardiogram (ECG) signal being in synchrony with the gantry rotation. In the animal study, a simple x-ray detector was used to generate a signal once per gantry rotation. This signal was then used to pace the pig's heart. Phantom studies demonstrated that partial scan artifacts are strongly dependent on the rotational symmetry of angular projections, which is determined by the object shape and composition and its position with respect to the isocenter. The degree of partial scan artifacts also depends on the location of the region of interest with respect to highly attenuating materials (bones, iodine, etc.) within the object. Single-source partial scan images (165 ms temporal resolution) were significantly less affected by partial scan artifacts compared to dual-source partial scan images (82 ms temporal resolution). When the ECG signal was in synchrony with the gantry rotation, the same cardiac phase always corresponded to the same positions of the x-ray tube(s) and, hence, the same scattering and beam hardening geometry. As a result, the range of image-to-image CT number variations for partial scan reconstruction images acquired in synchronized mode was decreased to that achieved using full reconstruction image data. The success of the new approach, which synchronizes the ECG signal with the position of the x-ray tube(s), was demonstrated both in the phantom and animal experiments.     

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.     

Adaptive temporal resolution optimization in helical cardiac cone beam CT reconstruction

Cone beam computed tomography scanners in combination with heart rate adaptive reconstruction schemes have the potential to enable cardiac volumetric computed tomography (CT) imaging for a larger number of patients and applications. In this publication, an adaptive scheme for the automatic and patient-specific reconstruction optimization is introduced to improve the temporal resolution and image quality. The optimization method permits the automatic determination of the required amount of gated helical cone beam projection data for the reconstruction volume. It furthermore allows one to optimize subvolume reconstruction yielding an increased temporal resolution. In addition, methods for the assessment of the temporal resolution are given which enable a quantitative documentation of the reconstruction improvements. Results are presented for patient data sets acquired in low pitch helical mode using a 16-slice cone beam CT system with parallel ECG recording.     

A knowledge-based cone-beam x-ray CT algorithm for dynamic volumetric cardiac imaging

With the introduction of spiral/helical multislice CT, medical x-ray CT began a transition into cone-beam geometry. The higher speed, thinner slice, and wider coverage with multislice/cone-beam CT indicate a great potential for dynamic volumetric imaging, with cardiac CT studies being the primary example. Existing ECG-gated cardiac CT algorithms have achieved encouraging results, but they do not utilize any time-varying anatomical information of the heart, and need major improvements to meet critical clinical needs. In this paper, we develop a knowledge-based spiral/helical multislice/cone-beam CT approach for dynamic volumetric cardiac imaging. This approach assumes the relationship between the cardiac status and the ECG signal, such as the volume of the left ventricle as a function of the cardiac phase. Our knowledge-based cardiac CT algorithm is evaluated in numerical simulation and patient studies. In the patient studies, the cardiac status is estimated initially from ECG data and subsequently refined with reconstructed images. Our results demonstrate significant image quality improvements in cardiac CT studies, giving clearly better clarity of the chamber boundaries and vascular structures. In conclusion, this approach seems promising for practical cardiac CT screening and diagnosis.     

4D micro-CT for cardiac and perfusion applications with view under sampling

Micro-CT is commonly used in preclinical studies to provide anatomical information. There is growing interest in obtaining functional measurements from 4D micro-CT. We report here strategies for 4D micro-CT with a focus on two applications: (i) cardiac imaging based on retrospective gating and (ii) pulmonary perfusion using multiple contrast injections/rotations paradigm. A dual source micro-CT system is used for image acquisition with a sampling rate of 20 projections per second. The cardiac micro-CT protocol involves the use of a liposomal blood pool contrast agent. Fast scanning of free breathing mice is achieved using retrospective gating. The ECG and respiratory signals are used to sort projections into ten cardiac phases. The pulmonary perfusion protocol uses a conventional contrast agent (Isovue 370) delivered by a micro-injector in four injections separated by 2 min intervals to allow for clearance. Each injection is synchronized with the rotation of the animal, and each of the four rotations is started with an angular offset of 22.5 from the starting angle of the previous rotation. Both cardiac and perfusion protocols result in an irregular angular distribution of projections that causes significant streaking artifacts in reconstructions when using traditional filtered backprojection (FBP) algorithms. The reconstruction involves the use of the point spread function of the micro-CT system for each time point, and the analysis of the distribution of the reconstructed data in the Fourier domain. This enables us to correct for angular inconsistencies via deconvolution and identify regions where data is missing. The missing regions are filled with data from a high quality but temporally averaged prior image reconstructed with all available projections. Simulations indicate that deconvolution successfully removes the streaking artifacts while preserving temporal information. 4D cardiac micro-CT in a mouse was performed with adequate image quality at isotropic voxel size of 88 μm and 10 ms temporal resolution. 4D pulmonary perfusion images were obtained in a mouse at 176 μm and 687 ms temporal resolution. Compared with FBP reconstruction, the streak reduction ratio is 70% and the contrast to noise ratio is 2.5 times greater in the deconvolved images. The radiation dose associated with the proposed methods is in the range of a typical micro-CT dose (0.17 Gy for the cardiac study and 0.21 Gy for the perfusion study). The low dose 4D micro-CT imaging presented here can be applied in high-throughput longitudinal studies in a wide range of applications, including drug safety and cardiopulmonary phenotyping.     

Image reconstruction with shift-variant filtration and its implication for noise and resolution properties in fan-beam computed tomography

In computed tomography (CT), the fan-beam filtered backprojection (FFBP) algorithm is used widely for image reconstruction. It is known that the FFBP algorithm can significantly amplify data noise and aliasing artifacts in situations where the focal lengths are comparable to or smaller than the size of the field of measurement (FOM). In this work, we propose an algorithm that is less susceptible to data noise, aliasing, and other data inconsistencies than is the FFBP algorithm while retaining the favorable resolution properties of the FFBP algorithm. In an attempt to evaluate the noise properties in reconstructed images, we derive analytic expressions for image variances obtained by use of the FFBP algorithm and the proposed algorithm. Computer simulation studies are conducted for quantitative evaluation of the spatial resolution and noise properties of images reconstructed by use of the algorithms. Numerical results of these studies confirm the favorable spatial resolution and noise properties of the proposed algorithm and verify the validity of the theoretically predicted image variances. The proposed algorithm and the derived analytic expressions for image variances can have practical implications for both estimation and detection/classification tasks making use of CT images, and they can readily be generalized to other fan-beam geometries.     

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.     

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.     

Fully 4D motion-compensated reconstruction of cardiac SPECT images

In this paper, we investigate the benefits of a spatiotemporal approach for reconstruction of image sequences. In the proposed approach, we introduce a temporal prior in the form of motion compensation to account for the statistical correlations among the frames in a sequence, and reconstruct all the frames collectively as a single function of space and time. The reconstruction algorithm is derived based on the maximum a posteriori estimate, for which the one-step late expectation-maximization algorithm is used. We demonstrated the method in our experiments using simulated single photon emission computed tomography (SPECT) cardiac perfusion images. The four-dimensional (4D) gated mathematical cardiac-torso phantom was used for simulation of gated SPECT perfusion imaging with Tc-99m-sestamibi. In addition to bias-variance analysis and time activity curves, we also used a channelized Hotelling observer to evaluate the detectability of perfusion defects in the reconstructed images. Our experimental results demonstrated that the incorporation of temporal regularization into image reconstruction could significantly improve the accuracy of cardiac images without causing any significant cross-frame blurring that may arise from the cardiac motion. This could lead to not only improved detection of perfusion defects, but also improved reconstruction of the heart wall which is important for functional assessment of the myocardium.     

16层螺旋CT冠状动脉图像质量评估

目的量化评估16层螺旋CT冠状动脉各分支图像质量,探讨冠状动脉CT成像的临床应用效果。方法采用16层螺旋CT对102例患者行回顾性心电门控冠状动脉成像,男60例,女42例,平均年龄（57．8±9．4）岁,平均心率（62．8±10．2）次／min。心率〈60次／min（n=40）为I组,60～70次／min（n=35）为Ⅱ组,〉70次／min（n=27）为Ⅲ组。4条冠状动脉分支（左主干、左前降支、左回旋支、右冠状动脉）分别用于图像质量分析。扫描原始数据以间隔10％在20％～80％时相分别回顾性重建冠状动脉图像,用横断面、曲面重建、容积再现等方法对冠状动脉显示率评估。用曲面重建方法测量冠状动脉各分支长度及近段和远段对比噪声比。结果①冠状动脉各分支平均显示长度：左主干（10．9±3．0）mm,左前降支（130．3±29．7）mm,左回旋支（82．8±19．8）mm,右冠状动脉（150．5±27．5）mm;②Ⅰ组60％和70％为最佳时相图像质量优,Ⅱ组60％时相为最佳时相图像质量优,Ⅲ组RCA较优的时相为40％,LM、LAD、LCX较优时相为60％。③所有冠状动脉分支平均对比噪声比10．9±3．2。结论16层螺旋CT有较好质量的冠状动脉图像,心率和时相影响冠状动脉的成像质量。     

Scanning-beam digital x-ray (SBDX) technology for interventional and diagnostic cardiac angiography

The scanning-beam digital x-ray (SBDX) system is designed for x-ray dose reduction in cardiac angiographic applications. Scatter reduction, efficient detection of primary x-rays, and an inverse beam geometry are the main components of the entrance dose reduction strategy. This paper reports the construction of an SBDX prototype, image reconstruction techniques, and measurements of spatial resolution and x-ray output. The x-ray source has a focal spot that is electronically scanned across a large-area transmission target. A multihole collimator beyond the target defines a series of x-ray beams directed at a distant small-area detector array. The prototype has a 23 cm X 23 cm target, 100 X 100 focal spot positions, and a 5 cm X 5 cm CdTe detector positioned 150 cm from the target. With this nonmechanical method of beam scanning, patient images with low detected scatter are generated at up to 30 frame/s. SBDX data acquisition is tomosynthetic. The prototype simultaneously reconstructs 16 planes spaced throughout the cardiac volume using shift-and-add backprojection. Image frames analogous to conventional projection images are generated with a multiplane compositing algorithm. Single-plane versus multiplane reconstruction of contrast-filled coronary arteries is demonstrated with images of the porcine heart. Phantom and porcine imaging studies show multiplane reconstruction is practicable under clinically realistic levels of patient attenuation and cardiac motion. The modulation transfer function for an in-plane slit at mechanical isocenter measured 0.41-0.56 at 1 cycle/mm, depending on the detector element to image pixel interpolation technique. Modeling indicates that desired gains in spatial resolution are achievable by halving the detector element width. The x-ray exposure rate 15 cm below isocenter, without table or patient in the beam, measured 11.5 R/min at 120 kVp, 24.3 kWp and 3.42 R/min at 70 kVp, 14.2 kWp.     

Iterative 4D cardiac micro-CT image reconstruction using an adaptive spatio-temporal sparsity prior

Temporal-correlated image reconstruction, also known as 4D CT image reconstruction, is a big challenge in computed tomography. The reasons for incorporating the temporal domain into the reconstruction are motions of the scanned object, which would otherwise lead to motion artifacts. The standard method for 4D CT image reconstruction is extracting single motion phases and reconstructing them separately. These reconstructions can suffer from undersampling artifacts due to the low number of used projections in each phase. There are different iterative methods which try to incorporate some a priori knowledge to compensate for these artifacts. In this paper we want to follow this strategy. The cost function we use is a higher dimensional cost function which accounts for the sparseness of the measured signal in the spatial and temporal directions. This leads to the definition of a higher dimensional total variation. The method is validated using in vivo cardiac micro-CT mouse data. Additionally, we compare the results to phase-correlated reconstructions using the FDK algorithm and a total variation constrained reconstruction, where the total variation term is only defined in the spatial domain. The reconstructed datasets show strong improvements in terms of artifact reduction and low-contrast resolution compared to other methods. Thereby the temporal resolution of the reconstructed signal is not affected.     

Half-scan fan-beam computed tomography with improved noise and resolution properties

Half-scan strategy can be used for reducing scanning time and radiation dose delivered to the patient in fan-beam computed tomography (CT). In helical CT, the data weighting/interpolation functions are often devised based upon half-scan configurations. The half-scan fan-beam filtered backprojection (FFBP) algorithm is generally used for image reconstruction from half-scan data. It can, however, be susceptible to sample aliasing and data noise for configurations with short focal lengths and/or large fan-angles, leading to nonuniform resolution and noise properties in reconstructed images. Uniform resolution and noise properties are generally desired because they may lead to an increased utility of reconstructed images in estimation and/or detection/classification tasks. In this work, we propose an algorithm for reconstruction of images with uniform noise and resolution properties in half-scan CT. In an attempt to evaluate the image-noise properties, we derive analytic expressions for image variances obtained by use of the half-scan algorithms. We also perform numerical studies to assess quantitatively the resolution and noise properties of the algorithms. The results in these studies confirm that the proposed algorithm yields images with more uniform spatial resolution and with lower and more uniform noise levels than does the half-scan FFBP algorithm. Empirical results obtained in noise studies also verify the validity of the derived expressions for image variances. The proposed algorithm would be particularly useful for image reconstruction from data acquired by use of configurations with short focal lengths and large field of measurement, which may be encountered in compact micro-CT and radiation therapeutic CT applications. The analytic results of the image-noise properties can be used for image-quality assessment in detection/classification tasks by use of model-observers.