A Feasibility Study of Low-Dose Single-Scan Dual-Energy Cone-Beam CT in Many-View Under-Sampling Framework

Real-time diagnosis of intracerebral hemorrhage by generating dual-energy CT from single-energy CT

Real-time diagnosis of intracerebral hemorrhage after thrombectomy is crucial for follow-up treatment. However, this is difficult to achieve with standard single-energy CT (SECT) due to similar CT values of blood and contrast agents under a single energy spectrum. In contrast, dual-energy CT (DECT) scanners employ two different energy spectra, which allows for real-time differentiation between hemorrhage and contrast extravasation based on energy-related attenuation characteristics. Unfortunately, DECT scanners are not as widely used as SECT scanners due to their high costs. To address this dilemma, in this paper, we generate pseudo DECT images from a SECT image for real-time diagnosis of hemorrhage. More specifically, we propose a SECT-to-DECT Transformer-based Generative Adversarial Network (SDTGAN), which is a 3D transformer-based multi-task learning framework equipped with a shared attention mechanism. In this way, SDTGAN can be guided to focus more on high-density areas (crucial for hemorrhage diagnosis) during the generation. Meanwhile, the introduced multi-task learning strategy and the shared attention mechanism also enable SDTGAN to model dependencies between interconnected generation tasks, improving generation performance while significantly reducing model parameters and computational complexity. In the experiments, we approximate real SECT images using mixed 120kV images from DECT data to address the issue of not being able to obtain the true paired DECT and SECT data. Extensive experiments demonstrate that SDTGAN can generate DECT images better than state-of-the-art methods. The code of our implementation is available at https://github.com/jiang-cw/SDTGAN.

Decomposition iteration strategy for low-dose CT denoising

In the medical field, computed tomography (CT) is a commonly used examination method, but the radiation generated increases the risk of illness in patients. Therefore, low-dose scanning schemes have attracted attention, in which noise reduction is essential. We propose a purposeful and interpretable decomposition iterative network (DISN) for low-dose CT denoising. This method aims to make the network design interpretable and improve the fidelity of details, rather than blindly designing or using deep CNN architecture. The experiment is trained and tested on multiple data sets. The results show that the DISN method can restore the low-dose CT image structure and improve the diagnostic performance when the image details are limited. Compared with other algorithms, DISN has better quantitative and visual performance, and has potential clinical application prospects.

A generalizable new figure of merit for dose optimization in dual energy cone beam CT scanning protocols

Objective. This study proposes and evaluates a new figure of merit (FOMn) for dose optimization of Dual-energy cone-beam CT (DE-CBCT) scanning protocols based on size-dependent modeling of radiation dose and multi-scale image quality.Approach. FOMn was defined using Z-score normalization and was proportional to the dose efficiency providing better multi-scale image quality, including comprehensive contrast-to-noise ratio (CCNR) and electron density (CED) for CatPhan604 inserts of various materials. Acrylic annuluses were combined with CatPhan604 to create four phantom sizes (diameters of the long axis are 200 mm, 270 mm, 350 mm, and 380 mm, respectively). DE-CBCT was decomposed using image-domain iterative methods based on Varian kV-CBCT images acquired using 25 protocols (100 kVp and 140 kVp combined with 5 tube currents).Main results. The accuracy of CED was approximately 1% for all protocols, but degraded monotonically with the increased phantom sizes. Combinations of lower voltage + higher current and higher voltage + lower current were optimal protocols balancing CCNR and dose. The most dose-efficient protocols for CED and CCNR were inconsistent, underlining the necessity of including multi-scale image quality in the evaluation and optimization of DE-CBCT. Pediatric and adult anthropomorphic phantom tests confirmed dose-efficiency of FOMn-recommended protocols.Significance. FOMn is a comprehensive metric that collectively evaluates radiation dose and multi-scale image quality for DE-CBCT. The models and data can also serve as lookup tables, suggesting personalized dose-efficient protocols for specific clinical imaging purposes.

Multi-Scale Feature Fusion Network for Low-Dose CT Denoising

Computed tomography (CT) is an imaging technique extensively used in medical treatment, but too much radiation dose in a CT scan will cause harm to the human body. Decreasing the dose of radiation will result in increased noise and artifacts in the reconstructed image, blurring the internal tissue and edge details. To get high-quality CT images, we present a multi-scale feature fusion network (MSFLNet) for low-dose CT (LDCT) denoising. In our MSFLNet, we combined multiple feature extraction modules, effective noise reduction modules, and fusion modules constructed using the attention mechanism to construct a horizontally connected multi-scale structure as the overall architecture of the network, which is used to construct different levels of feature maps at all scales. We innovatively define a composite loss function composed of pixel-level loss based on MS-SSIM-L1 and edge-based edge loss for LDCT denoising. In short, our approach learns a rich set of features that combine contextual information from multiple scales while maintaining the spatial details of denoised CT images. Our laboratory results indicate that compared with the existing methods, the peak signal-to-noise ratio (PSNR) value of CT images of the AAPM dataset processed by the new model is 33.6490, and the structural similarity (SSIM) value is 0.9174, which also achieves good results on the Piglet dataset with different doses. The results also show that the method removes noise and artifacts while effectively preserving CT images' architecture and grain information.

Dual-domain fusion deep convolutional neural network for low-dose CT denoising

BACKGROUND

In view of the underlying health risks posed by X-ray radiation, the main goal of the present research is to achieve high-quality CT images at the same time as reducing x-ray radiation. In recent years, convolutional neural network (CNN) has shown excellent performance in removing low-dose CT noise. However, previous work mainly focused on deepening and feature extraction work on CNN without considering fusion of features from frequency domain and image domain.

OBJECTIVE

To address this issue, we propose to develop and test a new LDCT image denoising method based on a dual-domain fusion deep convolutional neural network (DFCNN).

METHODS

This method deals with two domains, namely, the DCT domain and the image domain. In the DCT domain, we design a new residual CBAM network to enhance the internal and external relations of different channels while reducing noise to promote richer image structure information. For the image domain, we propose a top-down multi-scale codec network as a denoising network to obtain more acceptable edges and textures while obtaining multi-scale information. Then, the feature images of the two domains are fused by a combination network.

RESULTS

The proposed method was validated on the Mayo dataset and the Piglet dataset. The denoising algorithm is optimal in both subjective and objective evaluation indexes as compared to other state-of-the-art methods reported in previous studies.

CONCLUSIONS

The study results demonstrate that by applying the new fusion model denoising, denoising results in both image domain and DCT domain are better than other models developed using features extracted in the single image domain.

Simultaneous Image Reconstruction and Element Decomposition for Iodine Contrast Agent Visualization in Multienergy Element-Resolved Cone Beam CT

Iodine contrast agent is widely used in liver cancer radiotherapy at CT simulation stage to enhance detectability of tumor. However, its application in cone beam CT (CBCT) for image guidance before treatment delivery is still limited because of poor image quality and excessive dose of contrast agent during multiple treatment fractions. We previously developed a multienergy element-resolved (MEER) CBCT framework that included x-ray projection data acquisition on a conventional CBCT platform in a kVp-switching model and a dictionary-based image reconstruction algorithm that simultaneously reconstructed x-ray attenuation images at each kilovoltage peak (kVp), an electron density image, and elemental composition images. In this study, we investigated feasibility using MEER-CBCT for low-concentration iodine contrast agent visualization. We performed simulation and experimental studies using a phantom with inserts containing water and different concentrations of iodine solution and the MEER-CBCT scan with 600 projections in a full gantry rotation, in which the kVp level sequentially changed among 80, 100, and 120 kVps. We included iodine material in the dictionary of the reconstruction algorithm. We analyzed iodine detectability as quantified by contrast-to-noise ratio (CNR) and compared results with those of CBCT images reconstructed by the standard filter back projection (FBP) method with 600 projections. MEER-CBCT achieved similar contrast enhancement as FBP method but significantly higher CNR. At 2.5% iodine solution concentration, FBP method achieved 170 HU enhancement and CNR of 2.0, considered the standard CNR for successful tumor visualization. MEER-CBCT achieved the same CNR but at ~6.3 times lower iodine concentration of 0.4%.

Image quality of dual-energy cone-beam CT with total nuclear variation regularization

Despite the improvements in image quality of cone beam computed tomography (CBCT) scans, application remains limited to patient positioning. In this study, we propose to improve image quality by dual energy (DE) imaging and iterative reconstruction using least squares fitting with total variation (TV) regularization. The generalization of TV called total nuclear variation (TNV) was used to generate DE images. We acquired single energy (SE) and DE scans of an image quality phantom (IQP) and of an anthropomorphic human male phantom (HMP). The DE scans were dual arc acquisitions of 70kV and 130kV with a variable dose partitioning between low energy (LE) and high energy (HE) arcs. To investigate potential benefits from a larger spectral separation between LE and HE, DE scans with an additional 2 mm copper beam filtration in the HE arc were acquired for the IQP. The DE TNV scans were compared to SE scans reconstructed with FDK and iterative TV with varying parameters. The contrast-to-noise ratio (CNR), spatial frequency, and structural similarity (SSIM) were used as image quality metrics. Results showed largely improved image quality for DE TNV over FDK for both phantoms. DE TNV with the highest dose allocation in the LE arm yielded the highest CNR. Compared to SE TV, these DE TNV results had a slightly lower CNR with similar spatial resolution for the IQP. A decrease in the dose allocated to the LE arm improved the spatial resolution with a trade-off against CNR. For the HMP, DE TNV displayed a lower CNR and/or lower spatial resolution depending on the reconstruction parameters. Regarding the SSIM, DE TNV was superior to FDK and SE TV for both phantoms. The additional beam filtration for the IQP led to improved image quality in all metrics, surpassing the SE TV results in CNR and spatial resolution.

Segmental limited-angle CT reconstruction based on image structural prior

CT reconstruction from incomplete projection data is one of the key researches of X-ray CT imaging. The projection data acquired by few-view and limited-angle sampling are incomplete. In addition, few-view sampling often requires turning on and off the tube voltage, but rapid switching of tube voltage demands for high technical requirements. Limited-angle sampling is easy to realize. However, reconstructed images may encounter obvious artifacts. In this study we investigate a new segmental limited-angle (SLA) sampling strategy, which avoids rapid switching of tube voltage. Thus, the projection data has lower data correlation than limited-angle CT, which is conducive to reconstructing high-quality images. To suppress potential artifacts, we incorporate image structural prior into reconstruction model to present a reconstruction method. The limited-angle CT reconstruction experiments on digital phantoms, real carved cheese and walnut projections are used to test and verify the effectiveness of the proposed method. Several image quality evaluation indices including RMSE, PSNR, and SSIM of the reconstructions in simulation experiments are calculated and listed to show the superiority of our method. The experimental results indicate that the CT image reconstructed using the proposed new method is closer to the reference image. Images from real CT data and their residual images also show that applying the proposed new method can more effectively reduce artifacts and image structures are well preserved.

A prior image constraint robust principal component analysis reconstruction method for sparse segmental multi-energy computed tomography

Background

Multi-energy computed tomography (MECT) is a promising technique in medical imaging, especially for quantitative imaging. However, high technical requirements and system costs barrier its step into clinical practice.

Methods

We propose a novel sparse segmental MECT (SSMECT) scheme and corresponding reconstruction method, which is a cost-efficient way to realize MECT on a conventional single-source CT. For the data acquisition, the X-ray source is controlled to maintain an energy within a segmental arc, and then switch alternately to another energy level. This scan only needs to switch tube voltage a few times to acquire multi-energy data, but leads to sparse-view and limited-angle issues in image reconstruction. To solve this problem, we propose a prior image constraint robust principal component analysis (PIC-RPCA) reconstruction method, which introduces structural similarity and spectral correlation into the reconstruction.

Results

A numerical simulation and a real phantom experiment were conducted to demonstrate the efficacy and robustness of the scan scheme and reconstruction method. The results showed that our proposed reconstruction method could have achieved better multi-energy images than other competing methods both quantitatively and qualitatively.

Conclusions

Our proposed SSMECT scan with PIC-RPCA reconstruction method could lower kVp switching frequency while achieving satisfactory reconstruction accuracy and image quality.

Fast and effective single-scan dual-energy cone-beam CT reconstruction and decomposition denoising based on dual-energy vectorization

PURPOSE

Flat-panel detector (FPD) based dual-energy cone-beam computed tomography (DE-CBCT) is a promising imaging technique for dedicated clinical applications. In this paper, we proposed a fully analytical method for fast and effective single-scan DE-CBCT image reconstruction and decomposition.

METHODS

A rotatable Mo filter was inserted between an x-ray source and imaged object to alternately produce low and high-energy x-ray spectra. First, filtered-backprojection (FBP) method was applied on down-sampled projections to reconstruct low and high-energy images. Then, the two images were converted into a vectorized form represented with an amplitude and an argument image. Using amplitude image as a guide, a joint bilateral filter was applied to denoise the argument image. Then, high-quality dual-energy images were recovered from the amplitude image and the denoised argument image. Finally, the recovered dual-energy images were further used for low-noise material decomposition and electron density synthesis. Imaging was conducted on a Catphan^® 600 phantom and an anthropomorphic head phantom. The proposed method was evaluated via comparison with the traditional two-scan method and a commonly used filtering method (HYPR-LR).

RESULTS

On the Catphan^® 600 phantom, the proposed method successfully reduced streaking artifacts and preserved spatial resolution and noise-power-spectrum (NPS) pattern. In the electron density image, the proposed method increased contrast-to-noise ratio (CNR) by more than 2.5 times and achieved <1.2% error for electron density values. On the anthropomorphic head phantom, the proposed method greatly improved the soft-tissue contrast and the fine detail differentiation ability. In the selected ROIs on different human tissues, the differences between the CT number obtained by the proposed method and that by the two-scan method were less than 4 HU. In the material images, the proposed method suppressed noise by over 75.5% compared with two-scan results, and by over 40.4% compared with HYPR-LR results. Implementation of the whole algorithm took 44.5 s for volumetric imaging, including projection preprocessing, FBP reconstruction, joint bilateral filtering, and material decomposition.

CONCLUSIONS

Using down-sampled projections in single-scan DE-CBCT, the proposed method could effectively and efficiently produce high-quality DE-CBCT images and low-noise material decomposition images. This method demonstrated superior performance on spatial resolution enhancement, NPS preservation, noise reduction, and electron density accuracy, indicating better prospect in material differentiation and dose calculation.

Multi-Material Decomposition for Single Energy CT Using Material Sparsity Constraint

Multi-material decomposition (MMD) decomposes CT images into basis material images, and is a promising technique in clinical diagnostic CT to identify material compositions within the human body. MMD could be implemented on measurements obtained from spectral CT protocol, although spectral CT data acquisition is not readily available in most clinical environments. MMD methods using single energy CT (SECT), broadly applied in radiological departments of most hospitals, have been proposed in the literature while challenged by the inferior decomposition accuracy and the limited number of material bases due to the constrained material information in the SECT measurement. In this paper, we propose an image-domain SECT MMD method using material sparsity as an assistance under the condition that each voxel of the CT image contains at most two different elemental materials. L₀ norm represents the material sparsity constraint (MSC) and is integrated into the decomposition objective function with a least-square data fidelity term, total variation term, and a sum-to-one constraint of material volume fractions. An accelerated primal-dual (APD) algorithm with line-search scheme is applied to solve the problem. The pixelwise direct inversion method with the two-material assumption (TMA) is applied to estimate the initials. We validate the proposed method on phantom and patient data. Compared with the TMA method, the proposed MSC method increases the volume fraction accuracy (VFA) from 92.0% to 98.5% in the phantom study. In the patient study, the calcification area can be clearly visualized in the virtual non-contrast image generated by the proposed method, and has a similar shape to that in the ground-truth contrast-free CT image. The high decomposition image quality from the proposed method substantially facilitates the SECT-based MMD clinical applications.

A Novel Low-Dose Dual-Energy Imaging Method for a Fast-Rotating Gantry-Type CT Scanner

CT scan by use of a beam-filter placed between the x-ray source and the patient allows a single-scan low-dose dual-energy imaging with a minimal hardware modification to the existing CT systems. We have earlier demonstrated the feasibility of such imaging method with a multi-slit beam-filter reciprocating along the direction perpendicular to the CT rotation axis in a cone-beam CT system. However, such method would face mechanical challenges when the beam-filter is supposed to cooperate with a fast-rotating gantry in a diagnostic CT system. In this work, we propose a new scanning method and associated image reconstruction algorithm that can overcome these challenges. We propose to slide a beam-filter that has multi-slit structure with its slits being at a slanted angle with the CT gantry rotation axis during a scan. A streaky pattern would show up in the sinogram domain as a result. Using a notch filter in the Fourier domain of the sinogram, we removed the streaks and reconstructed an image by use of the filtered-backprojection algorithm. The remaining image artifacts were suppressed by applying l₀ norm based smoothing. Using this image as a prior, we have reconstructed low- and high-energy CT images in the iterative reconstruction framework. An image-based material decomposition then followed. We conducted a simulation study to test its feasibility using the XCAT phantom and also an experimental study using the Catphan phantom, a head phantom, an iodine-solution phantom, and a monkey in anesthesia, and showed its successful performance in image reconstruction and in material decomposition.

Estimating dual-energy CT imaging from single-energy CT data with material decomposition convolutional neural network

Dual-energy computed tomography (DECT) is of great significance for clinical practice due to its huge potential to provide material-specific information. However, DECT scanners are usually more expensive than standard single-energy CT (SECT) scanners and thus are less accessible to undeveloped regions. In this paper, we show that the energy-domain correlation and anatomical consistency between standard DECT images can be harnessed by a deep learning model to provide high-performance DECT imaging from fully-sampled low-energy data together with single-view high-energy data. We demonstrate the feasibility of the approach with two independent cohorts (the first cohort including contrast-enhanced DECT scans of 5753 image slices from 22 patients and the second cohort including spectral CT scans without contrast injection of 2463 image slices from other 22 patients) and show its superior performance on DECT applications. The deep-learning-based approach could be useful to further significantly reduce the radiation dose of current premium DECT scanners and has the potential to simplify the hardware of DECT imaging systems and to enable DECT imaging using standard SECT scanners.

Multislice spiral computed tomography imaging in evaluating hemophilic arthropathy

BACKGROUND

Certain hemophilia patients are unable to cooperate with or afford magnetic resonance imaging (MRI) examinations. The purpose of our study was to explore the value of multislice spiral computed tomography (MSCT) in evaluating hemophilic arthropathy (HA).

METHODS

Thirty-eight patients with 73 joints of HA were consecutively selected from January 2016 to May 2018 for this prospective study. All 73 joints were examined by X-ray, CT, and MRI within 2 days. The MRI scores of the joints were determined by the International Prophylaxis Study Group (IPSG) standard. The CT findings were quantified according to the IPSG standard, except for cartilage injury, which was quantified by joint space narrowing using the X-ray Pettersson score. The CT and MRI scores were compared by the paired Wilcoxon signed-rank test. The correlations between the CT score of joint space narrowing and MRI score of cartilage injury and the total CT and MRI scores were analyzed by Spearman rank correlation. The kappa test was used to compare the consistency of CT and MRI scores.

RESULTS

MRI was superior to CT based on the scores for small amount of effusion (P < 0.05), synovial hypertrophy and hemosiderin deposition in the mild groups (P < 0.05). The CT and MRI scores were not significantly different for moderate and massive effusion, synovial hypertrophy, and hemosiderin deposition in the moderate and severe groups, bone erosion or cystic changes (P > 0.05), and there was a high degree of consistency between the two scores (kappa > 0.81). The consistency between the Pettersson scores of joint space narrowing on CT and the IPSG scores of cartilage injury on MRI was high (kappa = 0. 774, P < 0.05).

CONCLUSION

The image scores of MSCT are generally consistent with MRI except for mild synovitis, which can be used as an alternative for the evaluation of HA.

Image domain multi-material decomposition using single energy CT

Multi-material decomposition (MMD) technique decomposes the CT images into basis material images and has been promising in clinical practice for material composition quantification within the human body. MMD could be implemented using the image data acquired from spectral CT or its special case, dual-energy CT (DECT) while the spectral CT data acquisition usually requires a hardware modification. In this paper, we propose an image domain MMD method using single energy CT (SECT). The proposed objective function applies a least square data fidelity term to enforce the minimization between the linear combination of decomposed material image and the measured SECT image, and an edge-preserving (EP) regularization term to meet the piecewise constant property of the material image. We apply the optimization transfer principle to form a pixel-wise separable quadratic surrogate (PWSQS) function in each iteration to decrease the objective function. The pixelwise direct inversion method assisted by the two-material assumption (TMA) is applied to obtain a good initial value. The proposed method is evaluated using a digital phantom, a Catphan phantom and the clinical data. A low-pass filtration method is implemented for a comparison purpose. In the phantom study, the proposed TMA method achieves high volume fraction accuracy (VFA) of 79.64% and the proposed EP method further increases the VFA by 15.56% and decreases the decomposition standard deviation (STD) by 81.51% compared with the TMA method. At the comparable noise level, the proposed EP method increases spatial resolution by an overall factor of 1.01 when the modulation transfer function magnitude is decreased to 50% compared with the low-pass filtration method. In clinical data study, the virtual non-contrast image generated by the proposed method achieves the root-mean-squared-relative error of 2.93% compared with the contrast-free ground-truth image.

Enhancement of soft-tissue contrast in cone-beam CT using an anti-scatter grid with a sparse sampling approach

PURPOSE

Anti-scatter grids suppress the scatter substantially thus improving image contrast in radiography. However, its active use in cone-beam CT for the purpose of improving contrast-to-noise ratio (CNR) has not been successful mainly due to the increased noise related to Poisson statistics of photons. This paper proposes a sparse-view scanning approach to address the above issue.

METHOD

Compared to the conventional cone-beam CT imaging framework, the proposed method reduces the number of projections and increases exposure in each projection to enhance image quality without an additional cost of radiation dose to patients. For image reconstruction from sparse-view data, an adaptive-steepest-descent projection-onto-convex-sets (ASD POCS) algorithm regularized by total-variation (TV) minimization was adopted. Contrast and CNR with various scattering conditions were evaluated in projection domain by a simulation study using GATE. Then we evaluated contrast, resolution, and image uniformity in CT image domain with Catphan phantom. A head phantom with soft-tissue structures was also employed for demonstrating a realistic application. A virtual grid-based estimation and reduction of scatter has also been implemented for comparison with the real anti-scatter grid.

RESULTS

In the projection domain evaluation, contrast and CNR enhancement was observed when using an anti-scatter grid compared to the virtual grid. In the CT image domain, the proposed method produced substantially higher contrast and CNR of the low-contrast structures with much improved image uniformity.

CONCLUSION

We have shown that the proposed method can provide high-quality CBCT images particularly with an increased contrast of soft-tissue at a neutral dose for image-guidance.

Image reconstruction for interrupted-beam x-ray CT on diagnostic clinical scanners

Low-dose x-ray CT is a major research area with high clinical impact. Compressed sensing using view-based sparse sampling and sparsity-promoting regularization has shown promise in simulations, but these methods can be difficult to implement on diagnostic clinical CT scanners since the x-ray beam cannot be switched on and off rapidly enough. An alternative to view-based sparse sampling is interrupted-beam sparse sampling. SparseCT is a recently-proposed interrupted-beam scheme that achieves sparse sampling by blocking a portion of the beam using a multislit collimator (MSC). The use of an MSC necessitates a number of modifications to the standard compressed sensing reconstruction pipeline. In particular, we find that SparseCT reconstruction is feasible within a model-based image reconstruction framework that incorporates data fidelity weighting to consider penumbra effects and source jittering to consider the effect of partial source obstruction. Here, we present these modifications and demonstrate their application in simulations and real-world prototype scans. In simulations compared to conventional low-dose acquisitions, SparseCT is able to achieve smaller normalized root-mean square differences and higher structural similarity measures on two reduction factors. In prototype experiments, we successfully apply our reconstruction modifications and maintain image resolution at quarter-dose reduction level. The SparseCT design requires only small hardware modifications to current diagnostic clinical scanners, opening up new possibilities for CT dose reduction.

Multiple limited-angles computed tomography reconstruction based on multi-direction total variation minimization

Accurate computed tomography (CT) reconstruction from incomplete projections is an important research topic. Sparse sampling and limited-angle sampling are two effective ways to reduce the X-ray radiation dose or the scanning time. However, it is technically complicated to realize sparse sampling in medical CT since the tube power or the pre-patient collimator is difficult to be switched frequently. Limited-angle sampling makes it difficult to reconstruct an accurate image. The developed multiple limited-angles (MLA) sampling scheme could well balance the technical implementation complexity and the CT reconstruction difficulty. It does not require frequent switching of the tube power or the pre-patient collimator. The data correlation of the acquired projections is lower than that in limited-angle sampling. Using the projections acquired by MLA sampling, CT images reconstructed by the total variation minimization (TVM) method suffer from shading artifacts. Because the artifacts are distributed in several fixed directions, the artifact-suppression reconstruction model based on multi-direction total variation was designed for MLA CT reconstruction in this work. The multi-direction total variation minimization (MDTVM) was utilized to solve the optimization model. Experiments on digital phantoms and real projections indicated that MDTVM can better suppress the shading artifacts than TVM.