Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT

Recent advances in murine cardiac studies with three-dimensional (3D) cone beam micro-CT used a retrospective gating technique. However, this sampling technique results in a limited number of projections with an irregular angular distribution due to the temporal resolution requirements and radiation dose restrictions. Both angular irregularity and undersampling complicate the reconstruction process, since they cause significant streaking artifacts. This work provides an iterative reconstruction solution to address this particular challenge. A sparseness prior regularized weighted l2 norm optimization is proposed to mitigate streaking artifacts based on the fact that most medical images are compressible. Total variation is implemented in this work as the regularizer for its simplicity. Comparison studies are conducted on a 3D cardiac mouse phantom generated with experimental data. After optimization, the method is applied to in vivo cardiac micro-CT data.     

A registration based approach for 4D cardiac micro-CT using combined prospective and retrospective gating

Recent advances in murine cardiac studies with three-dimensional cone beam micro-computed tomography (CT) have used either prospective or retrospective gating technique. While prospective gating ensures the best image quality and the highest resolution, it involves longer sampling times and higher radiation dose. Sampling is faster and the radiation dose can be reduced with retrospective gating but the image quality is affected by the limited number of projections with an irregular angular distribution which complicate the reconstruction process, causing significant streaking artifacts. This work involves both prospective and retrospective gating in sampling. Deformable registration is used between a high quality image set acquired with prospective gating with the multiple data sets during the cardiac cycle obtained using retrospective gating. Tests were conducted on a four-dimensional (4D) cardiac mouse phantom and after optimization, the method was applied to in vivo cardiac micro-CT data. Results indicate that, by using our method, the sampling time can be reduced by a factor of 2.5 and the radiation dose can be reduced 35% compared to the prospective sampling while the image quality can be maintained. In conclusion, we proposed a novel solution to 4D cine cardiac micro-CT based on a combined prospective with retrospective gating in sampling and deformable registration post reconstruction that mixed the advantages of both strategies.     

Streaking artifacts reduction in four-dimensional cone-beam computed tomography

Cone-beam computed tomography (CBCT) using an "on-board" x-ray imaging device integrated into a radiation therapy system has recently been made available for patient positioning, target localization, and adaptive treatment planning. One of the challenges for gantry mounted image-guided radiation therapy (IGRT) systems is the slow acquisition of projections for cone-beam CT (CBCT), which makes them sensitive to any patient motion during the scans. Aiming at motion artifact reduction, four-dimensional CBCT (4D CBCT) techniques have been introduced, where a surrogate for the target's motion profile is utilized to sort the cone-beam data by respiratory phase. However, due to the limited gantry rotation speed and limited readout speed of the on-board imager, fewer than 100 projections are available for the image reconstruction at each respiratory phase. Thus, severe undersampling streaking artifacts plague 4D CBCT images. In this paper, the authors propose a simple scheme to significantly reduce the streaking artifacts. In this method, a prior image is first reconstructed using all available projections without gating, in which static structures are well reconstructed while moving objects are blurred. The undersampling streaking artifacts from static structures are estimated from this prior image volume and then can be removed from the phase images using gated reconstruction. The proposed method was validated using numerical simulations, experimental phantom data, and patient data. The fidelity of stationary and moving objects is maintained, while large gains in streak artifact reduction are observed. Using this technique one can reconstruct 4D CBCT datasets using no more projections than are acquired in a 60 s scan. At the same time, a temporal gating window as narrow as 100 ms was utilized. Compared to the conventional 4D CBCT reconstruction, streaking artifacts were reduced by 60% to 70%.     

4D micro-CT using fast prospective gating

Micro-CT is currently used in preclinical studies to provide anatomical information. But, there is also significant interest in using this technology to obtain functional information. We report here a new sampling strategy for 4D micro-CT for functional cardiac and pulmonary imaging. Rapid scanning of free-breathing mice is achieved with fast prospective gating (FPG) implemented on a field programmable gate array. The method entails on-the-fly computation of delays from the R peaks of the ECG signals or the peaks of the respiratory signals for the triggering pulses. Projection images are acquired for all cardiac or respiratory phases at each angle before rotating to the next angle. FPG can deliver the faster scan time of retrospective gating (RG) with the regular angular distribution of conventional prospective gating for cardiac or respiratory gating. Simultaneous cardio-respiratory gating is also possible with FPG in a hybrid retrospective/prospective approach. We have performed phantom experiments to validate the new sampling protocol and compared the results from FPG and RG in cardiac imaging of a mouse. Additionally, we have evaluated the utility of incorporating respiratory information in 4D cardiac micro-CT studies with FPG. A dual-source micro-CT system was used for image acquisition with pulsed x-ray exposures (80 kVp, 100 mA, 10 ms). The cardiac micro-CT protocol involves the use of a liposomal blood pool contrast agent containing 123 mg I ml(-1) delivered via a tail vein catheter in a dose of 0.01 ml g(-1) body weight. The phantom experiment demonstrates that FPG can distinguish the successive phases of phantom motion with minimal motion blur, and the animal study demonstrates that respiratory FPG can distinguish inspiration and expiration. 4D cardiac micro-CT imaging with FPG provides image quality superior to RG at an isotropic voxel size of 88 μm and 10 ms temporal resolution. The acquisition time for either sampling approach is less than 5 min. The radiation dose associated with the proposed method is in the range of a typical micro-CT dose (256 mGy for the cardiac study). Ignoring respiration does not significantly affect anatomic information in cardiac studies. FPG can deliver short scan times with low-dose 4D micro-CT imaging without sacrificing image quality. FPG can be applied in high-throughput longitudinal studies in a wide range of applications, including drug safety and cardiopulmonary phenotyping.     

High temporal resolution and streak-free four-dimensional cone-beam computed tomography

Cone-beam computed tomography (CBCT) has been clinically used to verify patient position and to localize the target of treatment in image-guided radiation therapy (IGRT). However, when the chest and the upper abdomen are scanned, respiratory-induced motion blurring limits the utility of CBCT. In order to mitigate this blurring, respiratory-gated CBCT, i.e. 4D CBCT, was introduced. In 4D CBCT, the cone-beam projection data sets acquired during a gantry rotation are sorted into several respiratory phases. In these gated reconstructions, the number of projections for each respiratory phase is significantly reduced. Consequently, undersampling streaking artifacts are present in the reconstructed images, and the image contrast resolution is also significantly compromised. In this paper, we present a new method to simultaneously achieve both high temporal resolution ( approximately 100 ms) and streaking artifact-free image volumes in 4D CBCT. The enabling technique is a newly proposed image reconstruction method, i.e. prior image constrained compressed sensing (PICCS), which enables accurate image reconstruction using vastly undersampled cone-beam projections and a fully sampled prior image. Using PICCS, a streak-free image can be reconstructed from 10-20 cone-beam projections while the signal-to-noise ratio is determined by a denoising feature of the selected objective function and by the prior image, which is reconstructed using all of the acquired cone-beam projections. This feature of PICCS breaks the connection between the temporal resolution and streaking artifacts' level in 4D CBCT. Numerical simulations and experimental phantom studies have been conducted to validate the method.     

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.     

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.     

SElf‐gated Non‐Contrast‐Enhanced FUnctional Lung imaging (SENCEFUL) using a quasi‐random fast low‐angle shot (FLASH) sequence and proton MRI

Obtaining functional information on the human lung is of tremendous interest in the characterization of lung defects and pathologies. However, pulmonary ventilation and perfusion maps usually require contrast agents and the application of electrocardiogram (ECG) triggering and breath holds to generate datasets free of motion artifacts. This work demonstrates the possibility of obtaining highly resolved perfusion‐weighted and ventilation‐weighted images of the human lung using proton MRI and the SElf‐gated Non‐Contrast‐Enhanced FUnctional Lung imaging (SENCEFUL) technique. The SENCEFUL technique utilizes a two‐dimensional fast low‐angle shot (FLASH) sequence with quasi‐random sampling of phase‐encoding (PE) steps for data acquisition. After every readout, a short additional acquisition of the non‐phase‐encoded direct current (DC) signal necessary for self‐gating was added. By sorting the quasi‐randomly acquired data according to respiratory and cardiac phase derived from the DC signal, datasets of representative respiratory and cardiac cycles could be accurately reconstructed. By application of the Fourier transform along the temporal dimension, functional maps (perfusion and ventilation) were obtained. These maps were compared with dynamic contrast‐enhanced (DCE, perfusion) as well as standard Fourier decomposition (FD, ventilation) reference datasets. All datasets were additionally scored by two experienced radiologists to quantify image quality. In addition, one initial patient examination using SENCEFUL was performed. Functional images of healthy volunteers and a patient diagnosed with hypoplasia of the left pulmonary artery and left‐sided pulmonary fibrosis were successfully obtained. Perfusion‐weighted images corresponded well to DCE‐MRI data; ventilation‐weighted images offered a significantly better depiction of the lung periphery compared with standard FD. Furthermore, the SENCEFUL technique hints at a potential clinical relevance by successfully detecting a perfusion defect in the patient scan. It can be concluded that SENCEFUL enables highly resolved ventilation‐ and perfusion‐weighted maps of the human lung to be obtained using proton MRI, and might be interesting for further clinical evaluation. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Optimization of a retrospective technique for respiratory-gated high speed micro-CT of free-breathing rodents

The objective of this study was to develop a technique for dynamic respiratory imaging using retrospectively gated high-speed micro-CT imaging of free-breathing mice. Free-breathing C57Bl6 mice were scanned using a dynamic micro-CT scanner, comprising a flat-panel detector mounted on a slip-ring gantry. Projection images were acquired over ten complete gantry rotations in 50 s, while monitoring the respiratory motion in synchrony with projection-image acquisition. Projection images belonging to a selected respiratory phase were retrospectively identified and used for 3D reconstruction. The effect of using fewer gantry rotations--which influences both image quality and the ability to quantify respiratory function--was evaluated. Images reconstructed using unique projections from six or more gantry rotations produced acceptable images for quantitative analysis of lung volume, CT density, functional residual capacity and tidal volume. The functional residual capacity (0.15 +/- 0.03 mL) and tidal volumes (0.08 +/- 0.03 mL) measured in this study agree with previously reported measurements made using prospectively gated micro-CT and at higher resolution (150 microm versus 90 microm voxel spacing). Retrospectively gated micro-CT imaging of free-breathing mice enables quantitative dynamic measurement of morphological and functional parameters in the mouse models of respiratory disease, with scan times as short as 30 s, based on the acquisition of projection images over six gantry rotations.     

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.     

Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets

When the number of projections does not satisfy the Shannon/Nyquist sampling requirement, streaking artifacts are inevitable in x-ray computed tomography (CT) images reconstructed using filtered backprojection algorithms. In this letter, the spatial-temporal correlations in dynamic CT imaging have been exploited to sparsify dynamic CT image sequences and the newly proposed compressed sensing (CS) reconstruction method is applied to reconstruct the target image sequences. A prior image reconstructed from the union of interleaved dynamical data sets is utilized to constrain the CS image reconstruction for the individual time frames. This method is referred to as prior image constrained compressed sensing (PICCS). In vivo experimental animal studies were conducted to validate the PICCS algorithm, and the results indicate that PICCS enables accurate reconstruction of dynamic CT images using about 20 view angles, which corresponds to an under-sampling factor of 32. This undersampling factor implies a potential radiation dose reduction by a factor of 32 in myocardial CT perfusion imaging.     

Fast and exact 2D image reconstruction by means of Chebyshev decomposition and backprojection

A new algorithm for the reconstruction of two-dimensional (2D) images from projections is described. The algorithm is based on the decomposition of the projections into Chebyshev polynomials of the second kind, which are the ideal basis functions for this application. The Chebyshev decomposition is done via the fast discrete sine transform. A discrete reconstruction filter is applied that corresponds to the ramp filter used in standard filtered backprojection (FBP) reconstruction. In contrast to FBP, the filter is applied to the Chebyshev coefficients and not to the Fourier coefficients of the projections. Then the reconstructed image is simply obtained by means of backprojection. Consequently, the method can be considered as a Chebyshev-domain filtered backprojection (CD-FBP). The total calculation time is dominated by the backprojection step only and is comparable to FBP. The merits of CD-FBP as compared with standard FBP are that: (a) The result is exact if the 2D function to be reconstructed can be decomposed into polynomials of finite degree, and if the sampling is adequate. Otherwise a polynomial approximation results. (b) The algorithm is inherently discrete. (c) It is particularly well suited for reconstructions from projections with non-equidistant samples that occur for instance in 2D PET (positron emission tomography) imaging and in a special form of fan beam scanning. Examples of applications comprise reconstructions of the Shepp and Logan head phantom in various sampling geometries, and a real PET test object. In the PET example an increased resolution is observed in comparison with standard FBP.     

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.     

Prospects for quantitative computed tomography imaging in the presence of foreign metal bodies using statistical image reconstruction

X-ray computed tomography (CT) images of patients bearing metal intracavitary applicators or other metal foreign objects exhibit severe artifacts including streaks and aliasing. We have systematically evaluated via computer simulations the impact of scattered radiation, the polyenergetic spectrum, and measurement noise on the performance of three reconstruction algorithms: conventional filtered backprojection (FBP), deterministic iterative deblurring, and a new iterative algorithm, alternating minimization (AM), based on a CT detector model that includes noise, scatter, and polyenergetic spectra. Contrary to the dominant view of the literature, FBP streaking artifacts are due mostly to mismatches between FBP's simplified model of CT detector response and the physical process of signal acquisition. Artifacts on AM images are significantly mitigated as this algorithm substantially reduces detector-model mismatches. However, metal artifacts are reduced to acceptable levels only when prior knowledge of the metal object in the patient, including its pose, shape, and attenuation map, are used to constrain AM's iterations. AM image reconstruction, in combination with object-constrained CT to estimate the pose of metal objects in the patient, is a promising approach for effectively mitigating metal artifacts and making quantitative estimation of tissue attenuation coefficients a clinical possibility.     

Enhanced 4D cone-beam CT with inter-phase motion model

Four-dimensional (4D) cone-beam CT (CBCT) is commonly obtained by respiratory phase binning of the projections, followed by independent reconstructions of the rebinned data in each phase bin. Due to the significantly reduced number of projections per reconstruction, the quality of the 4DCBCT images is often degraded by view-aliasing artifacts easily seen in the axial view. Acquisitions using multiple gantry rotations or slow gantry rotation can increase the number of projections and substantially improve the 4D images. However, the extra cost of the scan time may set fundamental limits to their applications in clinics. Improving the trade-off between image quality and scan time is the key to making 4D onboard imaging practical and more useful. In this article, we present a novel technique toward high-quality 4DCBCT imaging without prolonging the acquisition time, referred to as the "enhanced 4DCBCT". The method correlates the data in different phase bins and integrates the internal motion into the 4DCBCT image formulation. Several strategies of the motion derivation are discussed, and the resultant images are assessed with numerical simulations as well as a clinical case.     

Analytical reconstruction of deconvolved Fourier rebinned PET sinograms

Fully 3D PET data are often rebinned into 2D data sets in order to avoid computationally intensive fully 3D reconstruction. Then, conventional 2D reconstruction techniques are employed to obtain images from the rebinned data. In a common scenario, 2D filtered back projection (FBP) is applied to Fourier rebinned (FORE) data. This approach is suboptimal because FBP is based on an idealized mathematical model of the data and cannot account for the statistical structure of data and noise. FORE data contain some blur in all three dimensions in comparison to conventional 2D PET data. In this work, we propose methods for approximating this blur in the sinogram domain due to FORE through its point spread function (PSF). We also explore simple methods for deconvolving the rebinned data with this PSF to restore it to a more ideal state prior to FBP. Our results show that deconvolution of the approximate transaxial PSF yields no improvement. When low image noise levels are required for detection tasks, the deconvolution of the axial PSF does not provide adequate resolution or quantitative benefits to justify its application. When accurate quantitation is required and higher noise levels are acceptable, the deconvolution of the axial PSF leads to considerable gains (30%) in accuracy over conventional FORE+FBP at matched noise levels.     

Low-dose CT reconstruction via edge-preserving total variation regularization

High radiation dose in computed tomography (CT) scans increases the lifetime risk of cancer and has become a major clinical concern. Recently, iterative reconstruction algorithms with total variation (TV) regularization have been developed to reconstruct CT images from highly undersampled data acquired at low mAs levels in order to reduce the imaging dose. Nonetheless, the low-contrast structures tend to be smoothed out by the TV regularization, posing a great challenge for the TV method. To solve this problem, in this work we develop an iterative CT reconstruction algorithm with edge-preserving TV (EPTV) regularization to reconstruct CT images from highly undersampled data obtained at low mAs levels. The CT image is reconstructed by minimizing energy consisting of an EPTV norm and a data fidelity term posed by the x-ray projections. The EPTV term is proposed to preferentially perform smoothing only on the non-edge part of the image in order to better preserve the edges, which is realized by introducing a penalty weight to the original TV norm. During the reconstruction process, the pixels at the edges would be gradually identified and given low penalty weight. Our iterative algorithm is implemented on graphics processing unit to improve its speed. We test our reconstruction algorithm on a digital NURBS-based cardiac-troso phantom, a physical chest phantom and a Catphan phantom. Reconstruction results from a conventional filtered backprojection (FBP) algorithm and a TV regularization method without edge-preserving penalty are also presented for comparison purposes. The experimental results illustrate that both the TV-based algorithm and our EPTV algorithm outperform the conventional FBP algorithm in suppressing the streaking artifacts and image noise under a low-dose context. Our edge-preserving algorithm is superior to the TV-based algorithm in that it can preserve more information of low-contrast structures and therefore maintain acceptable spatial resolution.     

Dual energy CT using slow kVp switching acquisition and prior image constrained compressed sensing

Dual energy computed tomography (DECT) is currently a subject of extensive investigation. DECT is currently implemented using either a dual source scanner with high and low kVp data acquired from separate sources or a single source scanner with both high and low kVp data acquired in an alternating manner. Both methods require dedicated hardware to enable data acquisition and image reconstruction for DECT. In this paper, we present a method to enable DECT using a single x-ray source with a slow kVp switching data acquisition. The enabling reconstruction technique allowing for the reduction in slew rate is the prior image constrained compressed sensing (PICCS) algorithm. When a slow kVp switching data acquisition method is used, the projection data with high and low kVp values are undersampled and the conventional filtered backprojection (FBP) image reconstruction does not enable streaking artifact-free images for material decomposition in DECT. In this paper, all of the acquired high and low kVp projection data were used to generate a prior image using the conventional FBP method. The PICCS algorithm was then used to reconstruct both high and low kVp images to enable material decomposition in the image domain. Both numerical simulations and physical phantom experimental studies were conducted to validate the proposed DECT scheme. The results demonstrate that a slew rate corresponding to 123 views at high and low kVp (high and low kVp values used for dual energy decomposition) is sufficient for the PICCS-based DECT method. In contrast, the slew rate should be high enough to obtain over 500 projections at each kVp for artifact-free reconstruction using an FBP-based DECT method.