Edge-preserving reconstruction from sparse projections of limited-angle computed tomography using ℓ0-regularized gradient prior

Image reconstruction method for incomplete CT projection based on self-guided image filtering

In some fields of medical diagnosis or industrial nondestructive testing, it is difficult to obtain complete computed tomography (CT) data due to the limitation of radiation dose or other factors. Therefore, image reconstruction of incomplete projection data is the focus of this paper. In this paper, a new image reconstruction model based on self-guided image filtering (SGIF) term is proposed for few-view and segmental limited-angle (SLA) CT reconstruction. Then the alternating direction method (ADM) is used to solve this model. For simplicity, we call it ADM-SGIF method. The key idea of ADM-SGIF method is to use the reconstructed image itself as a reference and utilize its structural features to guide CT reconstruction. This method can effectively preserve image structures and remove shading artifacts. To validate the effectiveness of the proposed reconstruction method, we conduct digital phantom and real CT data experiments. The results indicate that ADM-SGIF method outperforms competing methods, including total variation (TV), relative total variation (RTV), and L0-norm minimization solved by ADM (ADM-L0) methods, in both subjective and objective evaluations.

Parallel processing model for low-dose computed tomography image denoising

Low-dose computed tomography (LDCT) has gained increasing attention owing to its crucial role in reducing radiation exposure in patients. However, LDCT-reconstructed images often suffer from significant noise and artifacts, negatively impacting the radiologists' ability to accurately diagnose. To address this issue, many studies have focused on denoising LDCT images using deep learning (DL) methods. However, these DL-based denoising methods have been hindered by the highly variable feature distribution of LDCT data from different imaging sources, which adversely affects the performance of current denoising models. In this study, we propose a parallel processing model, the multi-encoder deep feature transformation network (MDFTN), which is designed to enhance the performance of LDCT imaging for multisource data. Unlike traditional network structures, which rely on continual learning to process multitask data, the approach can simultaneously handle LDCT images within a unified framework from various imaging sources. The proposed MDFTN consists of multiple encoders and decoders along with a deep feature transformation module (DFTM). During forward propagation in network training, each encoder extracts diverse features from its respective data source in parallel and the DFTM compresses these features into a shared feature space. Subsequently, each decoder performs an inverse operation for multisource loss estimation. Through collaborative training, the proposed MDFTN leverages the complementary advantages of multisource data distribution to enhance its adaptability and generalization. Numerous experiments were conducted on two public datasets and one local dataset, which demonstrated that the proposed network model can simultaneously process multisource data while effectively suppressing noise and preserving fine structures. The source code is available at https://github.com/123456789ey/MDFTN .

An interactive method based on multi-objective optimization for limited-angle CT reconstruction

Objective. Limited-angle x-ray computed tomography (CT) is a typical ill-posed inverse problem, leading to artifacts in the reconstructed image due to the incomplete projection data. Most iteration CT reconstruction methods involve optimization for a single object. This paper explores a multi-objective optimization model and an interactive method based on multi-objective optimization to suppress the artifacts of limited-angle CT.Approach. The model includes two objective functions on the dual domain within the data consistency constraint. In the interactive method, the structural similarity index measure (SSIM) is regarded as the value function of the decision maker (DM) firstly. Secondly, the DM arranges the objective functions of the multi-objective optimization model to be optimized according to their absolute importance. Finally, the SSIM and the simulated annealing (SA) method help the DM choose the desirable reconstruction image by improving the SSIM value during the iteration process.Main results. Simulation and real data experiments demonstrate that the artifacts can be suppressed by the proposed method, and the results were superior to those reconstructed by the other three reconstruction methods in preserving the edge structure of the image.Significance. The proposed interactive method based on multi-objective optimization shows some potential advantages over classical single object optimization methods.

Sparse-view cone-beam computed tomography iterative reconstruction based on new multi-gradient direction total variation

AIM

The accurate reconstruction of cone-beam computed tomography (CBCT) from sparse projections is one of the most important areas for study. The compressed sensing theory has been widely employed in the sparse reconstruction of CBCT. However, the total variation (TV) approach solely uses information from the i-coordinate, j-coordinate, and k-coordinate gradients to reconstruct the CBCT image.

MATERIALS AND METHODS

It is well recognized that the CBCT image can be reconstructed more accurately with more gradient information from different directions. Thus, this study introduces a novel approach, named the new multi-gradient direction total variation minimization method. The method uses gradient information from the ij-coordinate, ik-coordinate, and jk-coordinate directions to reconstruct CBCT images, which incorporates nine different types of gradient information from nine directions.

RESULTS

This study assessed the efficacy of the proposed methodology using under-sampled projections from four different experiments, including two digital phantoms, one patient's head dataset, and one physical phantom dataset. The results indicated that the proposed method achieved the lowest RMSE index and the highest SSIM index. Meanwhile, we compared the voxel intensity curves of the reconstructed images to assess the edge structure preservation. Among the various methods compared, the curves generated by the proposed method exhibited the highest level of consistency with the gold standard image curves.

CONCLUSION

In summary, the proposed method showed significant potential in enhancing the quality and accuracy of CBCT image reconstruction.

Bilateral Weighted Relative Total Variation for Low-Dose CT Reconstruction

Low-dose computed tomography (LDCT) has been widely used for various clinic applications to reduce the X-ray dose absorbed by patients. However, LDCT is usually degraded by severe noise over the image space. The image quality of LDCT has attracted aroused attentions of scholars. In this study, we propose the bilateral weighted relative total variation (BRTV) used for image restoration to simultaneously maintain edges and further reduce noise, then propose the BRTV-regularized projections onto convex sets (POCS-BRTV) model for LDCT reconstruction. Referring to the spacial closeness and the similarity of gray value between two pixels in a local rectangle, POCS-BRTV can adaptively extract sharp edges and minor details during the iterative reconstruction process. Evaluation indexes and visual effects are used to measure the performances among different algorithms. Experimental results indicate that the proposed POCS-BRTV model can achieve superior image quality than the compared algorithms in terms of the structure and texture preservation.

Structure-guided computed tomography reconstruction from limited-angle projections

Limited-angle computed tomography (CT) imaging is one of the common imaging problems. The reconstructed images often encounter obvious artifacts and structure degradation. In recent years, the recoverability prior of image structure has been widely explored in limited-angle CT reconstruction, and the image quality has been greatly improved. However, the artifacts and structure degradation still exist. In this study, we establish a new reconstruction model based on weighted relative structure (wRS) determined by image gradients, which serves as weights to guide image reconstruction in order to reduce artifacts and preserve structures. Then, we develop an efficient algorithm using a surrogate function to solve this model. Moreover, this method is compared with some of other popular reconstruction methods, such as anisotropic total variation method and image gradient L0 norm minimization method and so on. Experiments on digital phantoms, real carved cheese and walnut projection are reported to demonstrate its superiority. Several quantitative indices including RMSE, PSNR, and SSIM of the reconstruction images from 90°-data of FORBILD head phantom are 0.0120, 43.52, and 0.9961. The experimental results indicate that the image obtained by our method is the closest to reference image. By comparing reconstruction images or their residual images, images reconstructed from real CT data, the experimental results of the residual images and the respective quantitative data analysis also demonstrate that the images reconstructed using our new method suffer from less artifacts and structure degradation.

Superiorization-inspired unrolled SART algorithm with U-Net generated perturbations for sparse-view and limited-angle CT reconstruction

Objective.Unrolled algorithms are a promising approach for reconstruction of CT images in challenging scenarios, such as low-dose, sparse-view and limited-angle imaging. In an unrolled algorithm, a fixed number of iterations of a reconstruction method are unrolled into multiple layers of a neural network, and interspersed with trainable layers. The entire network is then trained end-to-end in a supervised fashion, to learn an appropriate regularizer from training data. In this paper we propose a novel unrolled algorithm, and compare its performance with several other approaches on sparse-view and limited-angle CT.Approach.The proposed algorithm is inspired by the superiorization methodology, an optimization heuristic in which iterates of a feasibility-seeking method are perturbed between iterations, typically using descent directions of a model-based penalty function. Our algorithm instead uses a modified U-net architecture to introduce the perturbations, allowing a network to learn beneficial perturbations to the image at various stages of the reconstruction, based on the training data.Main Results.In several numerical experiments modeling sparse-view and limited angle CT scenarios, the algorithm provides excellent results. In particular, it outperforms several competing unrolled methods in limited-angle scenarios, while providing comparable or better performance on sparse-view scenarios.Significance.This work represents a first step towards exploiting the power of deep learning within the superiorization methodology. Additionally, it studies the effect of network architecture on the performance of unrolled methods, as well as the effectiveness of the unrolled approach on both limited-angle CT, where previous studies have primarily focused on the sparse-view and low-dose cases.

DIOR: Deep Iterative Optimization-Based Residual-Learning for Limited-Angle CT Reconstruction

Limited-angle CT is a challenging problem in real applications. Incomplete projection data will lead to severe artifacts and distortions in reconstruction images. To tackle this problem, we propose a novel reconstruction framework termed Deep Iterative Optimization-based Residual-learning (DIOR) for limited-angle CT. Instead of directly deploying the regularization term on image space, the DIOR combines iterative optimization and deep learning based on the residual domain, significantly improving the convergence property and generalization ability. Specifically, the asymmetric convolutional modules are adopted to strengthen the feature extraction capacity in smooth regions for deep priors. Besides, in our DIOR method, the information contained in low-frequency and high-frequency components is also evaluated by perceptual loss to improve the performance in tissue preservation. Both simulated and clinical datasets are performed to validate the performance of DIOR. Compared with existing competitive algorithms, quantitative and qualitative results show that the proposed method brings a promising improvement in artifact removal, detail restoration and edge preservation.

Multi-domain integrative Swin transformer network for sparse-view tomographic reconstruction

Decreasing projection views to a lower X-ray radiation dose usually leads to severe streak artifacts. To improve image quality from sparse-view data, a multi-domain integrative Swin transformer network (MIST-net) was developed and is reported in this article. First, MIST-net incorporated lavish domain features from data, residual data, image, and residual image using flexible network architectures, where a residual data and residual image sub-network was considered as a data consistency module to eliminate interpolation and reconstruction errors. Second, a trainable edge enhancement filter was incorporated to detect and protect image edges. Third, a high-quality reconstruction Swin transformer (i.e., Recformer) was designed to capture image global features. The experimental results on numerical and real cardiac clinical datasets with 48 views demonstrated that our proposed MIST-net provided better image quality with more small features and sharp edges than other competitors.

Tensor Dictionary Learning with an Enhanced Sparsity Constraint for Sparse-View Spectral CT Reconstruction

Spectral computed tomography (CT) can divide collected photons into multi-energy channels and gain multi-channel projections synchronously by using photon-counting detectors. However, reconstructed images usually contain severe noise due to the limited number of photons in the corresponding energy channel. Tensor dictionary learning (TDL)-based methods have achieved better performance, but usually lose image edge information and details, especially from an under-sampling dataset. To address this problem, this paper proposes a method termed TDL with an enhanced sparsity constraint for spectral CT reconstruction. The proposed algorithm inherits the superiority of TDL by exploring the correlation of spectral CT images. Moreover, the method designs a regularization using the L0-norm of the image gradient to constrain images and the difference between images and a prior image in each energy channel simultaneously, further improving the ability to preserve edge information and subtle image details. The split-Bregman algorithm has been applied to address the proposed objective minimization model. Several numerical simulations and realistic preclinical mice are studied to assess the effectiveness of the proposed algorithm. The results demonstrate that the proposed method improves the quality of spectral CT images in terms of noise elimination, edge preservation, and image detail recovery compared to the several existing better methods.

Image-domain Material Decomposition for Spectral CT using a Generalized Dictionary Learning

The spectral computed tomography (CT) has huge advantages by providing accurate material information. Unfortunately, due to the instability or overdetermination of material decomposition model, the accuracy of material decomposition can be compromised in practice. Very recently, the dictionary learning based image-domain material decomposition (DLIMD) can obtain high accuracy for material decompositions from reconstructed spectral CT images. This method can explore the correlation of material components to some extent by training a unified dictionary from all material images. In addition, the dictionary learning based prior as a penalty is applied on material components independently, and many parameters would be carefully elaborated in practice. Because the concentration of contrast agent in clinical applications is low, it can result in data inconsistency for dictionary based representation during the iteration process. To avoid the aforementioned limitations and further improve the accuracy of materials, we first construct a generalized dictionary learning based image-domain material decomposition (GDLIMD) model. Then, the material tensor image is unfolded along the mode-1 to enhance the correlation of different materials. Finally, to avoid the data inconsistency of low iodine contrast, a normalization strategy is employed. Both physical phantom and tissue-synthetic phantom experiments demonstrate the proposed GDLIMD method outperforms the DLIMD and direct inversion (DI) methods.

Anisotropic structure property based image reconstruction method for limited-angle computed tomography

Limited-angle computed tomography (CT) may appear in restricted CT scans. Since the available projection data is incomplete, the images reconstructed by filtered back-projection (FBP) or algebraic reconstruction technique (ART) often encounter shading artifacts. However, using the anisotropy property of the shading artifacts that coincide with the characteristic of limited-angle CT images can reduce the shading artifacts. Considering this concept, we combine the anisotropy property of the shading artifacts with the anisotropic structure property of an image to develop a new algorithm for image reconstruction. Specifically, we propose an image reconstruction method based on adaptive weighted anisotropic total variation (AwATV). This method, termed as AwATV method for short, is designed to preserve image structures and then remove the shading artifacts. It characterizes both of above properties. The anisotropy property of the shading artifacts accounts for reducing artifacts, and the anisotropic structure property of an image accounts for preserving structures. In order to evaluate the performance of AwATV, we use the simulation projection data of FORBILD head phantom and real CT data for image reconstruction. Experimental results show that AwATV can always reconstruct images with higher SSIM and PSNR, and smaller RMSE, which means that AwATV enables to reconstruct images with higher quality in term of artifact reduction and structure preservation.

Low-dose spectral CT reconstruction based on image-gradient L0-norm and adaptive spectral PICCS

The photon-counting detector based spectral computed tomography (CT) is promising for lesion detection, tissue characterization, and material decomposition. However, the lower signal-to-noise ratio within multi-energy projection dataset can result in poorly reconstructed image quality. Recently, as prior information, a high-quality spectral mean image was introduced into the prior image constrained compressed sensing (PICCS) framework to suppress noise, leading to spectral PICCS (SPICCS). In the original SPICCS model, the image gradient L₁-norm is employed, and it can cause blurred edge structures in the reconstructed images. Encouraged by the advantages in edge preservation and finer structure recovering, the image gradient L₀-norm was incorporated into the PICCS model. Furthermore, due to the difference of energy spectrum in different channels, a weighting factor is introduced and adaptively adjusted for different channel-wise images, leading to an L₀-norm based adaptive SPICCS (L₀-ASPICCS) algorithm for low-dose spectral CT reconstruction. The split-Bregman method is employed to minimize the objective function. Extensive numerical simulations and physical phantom experiments are performed to evaluate the proposed method. By comparing with the state-of-the-art algorithms, such as the simultaneous algebraic reconstruction technique, total variation minimization, and SPICCS, the advantages of our proposed method are demonstrated in terms of both qualitative and quantitative evaluation results.

An Algorithm of l 1-Norm and l 0-Norm Regularization Algorithm for CT Image Reconstruction from Limited Projection

The l ₁-norm regularization has attracted attention for image reconstruction in computed tomography. The l ₀-norm of the gradients of an image provides a measure of the sparsity of gradients of the image. In this paper, we present a new combined l ₁-norm and l ₀-norm regularization model for image reconstruction from limited projection data in computed tomography. We also propose an algorithm in the algebraic framework to solve the optimization effectively using the nonmonotone alternating direction algorithm with hard thresholding method. Numerical experiments indicate that this new algorithm makes much improvement by involving l ₀-norm regularization.

Deep learning based image reconstruction algorithm for limited-angle translational computed tomography

As a low-end computed tomography (CT) system, translational CT (TCT) is in urgent demand in developing countries. Under some circumstances, in order to reduce the scan time, decrease the X-ray radiation or scan long objects, furthermore, to avoid the inconsistency of the detector for the large angle scanning, we use the limited-angle TCT scanning mode to scan an object within a limited angular range. However, this scanning mode introduces some additional noise and limited-angle artifacts that seriously degrade the imaging quality and affect the diagnosis accuracy. To reconstruct a high-quality image for the limited-angle TCT scanning mode, we develop a limited-angle TCT image reconstruction algorithm based on a U-net convolutional neural network (CNN). First, we use the SART method to the limited-angle TCT projection data, then we import the image reconstructed by SART method to a well-trained CNN which can suppress the artifacts and preserve the structures to obtain a better reconstructed image. Some simulation experiments are implemented to demonstrate the performance of the developed algorithm for the limited-angle TCT scanning mode. Compared with some state-of-the-art methods, the developed algorithm can effectively suppress the noise and the limited-angle artifacts while preserving the image structures.

Image gradient L0-norm based PICCS for swinging multi-source CT reconstruction

Dynamic computed tomography (CT) is usually employed to image motion objects, such as beating heart, coronary artery and cerebral perfusion, etc. Recently, to further improve the temporal resolution for aperiodic industrial process imaging, the swinging multi-source CT (SMCT) systems and the corresponding swinging multi-source prior image constrained compressed sensing (SM-PICCS) method were developed. Since the SM-PICCS uses the L1-norm of image gradient, the edge structures in the reconstructed images are blurred and motion artifacts are still present. Inspired by the advantages in terms of image edge preservation and fine structure recovering, the L0-norm of image gradient is incorporated into the prior image constrained compressed sensing, leading to an L0-PICCS Algorithm 1Table 1The parameters of L0-PICCS (δ1,δ2,λ1*,λ2*) for numerical simulation.Sourceswδ1(10-2)δ2(10-2)λ1*(10-2)λ2*(10-8)Noise-free510522.001.525522.001.55035002.00471014.33332.00500025522.00500050222.005000Noise51062002.505002554502.501.55054502.901.571027.385.91.5810000258.285.91.5850050522.001.5. The experimental results confirm that the L0-PICCS outperforms the SM-PICCS in both visual inspection and quantitative analysis.

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.

Impact of statistical reconstruction and compressed sensing algorithms on projection data elimination during X-ray CT image reconstruction

OBJECTIVES

To examine the effect of incomplete, or total elimination of, projection data on computed tomography (CT) images subjected to statistical reconstruction and/or compressed sensing algorithms.

METHODS

Multidetector row CT images were used. The algebraic reconstruction technique (ART) and the maximum likelihood-expectation maximization (ML-EM) method were compared with filtered back-projection (FBP). Effects on reconstructed images were studied when the projection data of 360° (360 projections) were decreased to 180 or 90 projections by reducing the collection angle or thinning the image data. The total variation (TV) regularization method using compressed sensing was applied to images processed by the ART. Image noise was subjectively evaluated using the root-mean-square error and signal-to-noise ratio.

RESULTS

When projection data were reduced by one-half or three-quarters, ART and ML-EM produced better image quality than FBP. Both ART and ML-EM resulted in high quality at a spread of 90 projections over 180° rotation. Computational loading was high for statistical reconstruction, but not for ML-EM, compared with the ART. TV regularization made it possible to use only 36 projections while still achieving acceptable image quality.

CONCLUSIONS

Incomplete projection data-accomplished by reducing the angle to collect image data or thinning the projection data without reducing the angle of rotation over which it is collected-made it possible to reduce the radiation dose while retaining image quality with statistical reconstruction algorithms and/or compressed sensing. Despite heavier computational calculation loading, these methods should be considered for reducing radiation doses.

Low-dose spectral CT reconstruction using image gradient ℓ 0-norm and tensor dictionary

Spectral computed tomography (CT) has a great superiority in lesion detection, tissue characterization and material decomposition. To further extend its potential clinical applications, in this work, we propose an improved tensor dictionary learning method for low-dose spectral CT reconstruction with a constraint of image gradient ℓ ₀-norm, which is named as ℓ ₀TDL. The ℓ ₀TDL method inherits the advantages of tensor dictionary learning (TDL) by employing the similarity of spectral CT images. On the other hand, by introducing the ℓ ₀-norm constraint in gradient image domain, the proposed method emphasizes the spatial sparsity to overcome the weakness of TDL on preserving edge information. The split-bregman method is employed to solve the proposed method. Both numerical simulations and real mouse studies are perform to evaluate the proposed method. The results show that the proposed ℓ ₀TDL method outperforms other competing methods, such as total variation (TV) minimization, TV with low rank (TV+LR), and TDL methods.

Spatial-Spectral Cube Matching Frame for Spectral CT Reconstruction

Spectral computed tomography (CT) reconstructs the same scanned object from projections of multiple narrow energy windows, and it can be used for material identification and decomposition. However, the multi-energy projection dataset has a lower signal-noise-ratio (SNR), resulting in poor reconstructed image quality. To address this thorny problem, we develop a spectral CT reconstruction method, namely spatial-spectral cube matching frame (SSCMF). This method is inspired by the following three facts: i) human body usually consists of two or three basic materials implying that the reconstructed spectral images have a strong sparsity; ii) the same basic material component in a single channel image has similar intensity and structures in local regions. Different material components within the same energy channel share similar structural information; iii) multi-energy projection datasets are collected from the subject by using different narrow energy windows, which means images reconstructed from different energy-channels share similar structures. To explore those information, we first establish a tensor cube matching frame (CMF) for a BM4D denoising procedure. Then, as a new regularizer, the CMF is introduced into a basic spectral CT reconstruction model, generating the SSCMF method. Because the SSCMF model contains an L0-norm minimization of 4D transform coefficients, an effective strategy is employed for optimization. Both numerical simulations and realistic preclinical mouse studies are performed. The results show that the SSCMF method outperforms the state-of-the-art algorithms, including the simultaneous algebraic reconstruction technique, total variation minimization, total variation plus low rank, and tensor dictionary learning.

Edge information diffusion based reconstruction (EIDR) for cone beam computed laminography

Computed laminography (CL) is a prospective nondestructive testing technique for flat object inspection in industrial applications. However, CL image reconstruction is a challenging task because incomplete projection data are acquired from the CL scan. When a conventional computed tomography (CT) reconstruction method is applied to cone beam CL data, the vertical edges (singularities in the z-direction) in the reconstructed image would be blurred. On the contrary, the horizontal edges (singularities within slices) can be quite accurately reconstructed. Based on this key observation, an edge information diffusion method is developed, which fixes the horizontal edges and propagates their values within the slices. An effective CL reconstruction method is then proposed for flat object inspection by combining the edge information diffusion procedure, which plays the role of regularization, with conventional CT image reconstruction algorithms. Experiments on both simulated data and real data are performed to verify the effectiveness of the proposed method. The results show that the proposed method can effectively suppress the inter-slice aliasing and blurring caused by incompleteness of the CL scan data, and that it outperforms other state-of-the-art methods.

Geometric artifacts reduction for cone-beam CT via L0-norm minimization without dedicated phantoms

For cone-beam computed tomography (CBCT), transversal shifts of the rotation center exist inevitably, which will result in geometric artifacts in CT images. In this work, we propose a novel geometric calibration method for CBCT, which can also be used in micro-CT. The symmetry property of the sinogram is used for the first calibration, and then L0-norm of the gradient image from the reconstructed image is used as the cost function to be minimized for the second calibration. An iterative search method is adopted to pursue the local minimum of the L0-norm minimization problem. The transversal shift value is updated with affirmatory step size within a search range determined by the first calibration. In addition, graphic processing unit (GPU)-based FDK algorithm and acceleration techniques are designed to accelerate the calibration process of the presented new method. In simulation experiments, the mean absolute difference (MAD) and the standard deviation (SD) of the transversal shift value were less than 0.2 pixels between the noise-free and noisy projection images, which indicated highly accurate calibration applying the new calibration method. In real data experiments, the smaller entropies of the corrected images also indicated that higher resolution image was acquired using the corrected projection data and the textures were well protected. Study results also support the feasibility of applying the proposed method to other imaging modalities.