Tiny a priori knowledge solves the interior problem in computed tomography

Deep-Interior: A new pathway to interior tomographic image reconstruction via a weighted backprojection and deep learning

BACKGROUND

In recent years, deep learning strategies have been combined with either the filtered backprojection or iterative methods or the direct projection-to-image by deep learning only to reconstruct images. Some of these methods can be applied to address the interior reconstruction problems for centered regions of interest (ROIs) with fixed sizes. Developing a method to enable interior tomography with arbitrarily located ROIs with nearly arbitrary ROI sizes inside a scanning field of view (FOV) remains an open question.

PURPOSE

To develop a new pathway to enable interior tomographic reconstruction for arbitrarily located ROIs with arbitrary sizes using a single trained deep neural network model.

METHODS

The method consists of two steps. First, an analytical weighted backprojection reconstruction algorithm was developed to perform domain transform from divergent fan-beam projection data to an intermediate image feature space,B ( x ⃗ ) $B(\vec{x})$ , for an arbitrary size ROI at an arbitrary location inside the FOV. Second, a supervised learning technique was developed to train a deep neural network architecture to perform deconvolution to obtain the true imagef ( x ⃗ ) $f(\vec{x})$ from the new feature spaceB ( x ⃗ ) $B(\vec{x})$ . This two-step method is referred to as Deep-Interior for convenience. Both numerical simulations and experimental studies were performed to validate the proposed Deep-Interior method.

RESULTS

The results showed that ROIs as small as a diameter of 5 cm could be accurately reconstructed (similarity index 0.985 ± 0.018 on internal testing data and 0.940 ± 0.025 on external testing data) at arbitrary locations within an imaging object covering a wide variety of anatomical structures of different body parts. Besides, ROIs of arbitrary size can be reconstructed by stitching small ROIs without additional training.

CONCLUSION

The developed Deep-Interior framework can enable interior tomographic reconstruction from divergent fan-beam projections for short-scan and super-short-scan acquisitions for small ROIs (with a diameter larger than 5 cm) at an arbitrary location inside the scanning FOV with high quantitative reconstruction accuracy.

Hybrid source translation scanning mode for interior tomography

Interior tomography is a promising technique that can be used to image large objects with high acquisition efficiency. However, it suffers from truncation artifacts and attenuation value bias due to the contribution from the parts of the object outside the ROI, which compromises its ability of quantitative evaluation in material or biological studies. In this paper, we present a hybrid source translation scanning mode for interior tomography, called hySTCT-where the projections inside the ROI and outside the ROI are finely sampled and coarsely sampled respectively to mitigate truncation artifacts and value bias within the ROI. Inspired by our previous work-virtual projection-based filtered backprojection (V-FBP) algorithm, we develop two reconstruction methods-interpolation V-FBP (iV-FBP) and two-step V-FBP (tV-FBP)-based on the linearity property of the inverse Radon transform for hySTCT reconstruction. The experiments demonstrate that the proposed strategy can effectively suppress truncated artifacts and improve the reconstruction accuracy within the ROI.

Interior Reconstruction from Truncated Projection Data in Cone-beam Computed Tomography

The interior reconstruction of completely truncated projection data is a frontier research hotspot in cone-beam computed tomography (CBCT) application. It is difficult to find a method with acceptable accuracy and high efficiency to solve it. Based on the simplified algebraic reconstruction technique (S-ART) algorithm and the filtered back projection (FBP) algorithm with the new filter, an efficient and feasible interior reconstruction algorithm is proposed in this paper. The algorithm uses the S-ART algorithm to quickly recover the complete projection data and then uses the new ramp filter which can suppress the high-frequency noise in the projection data to filter the recovered complete projection data. Finally, the interior reconstructed images are obtained by back projection. The computational complexity of the proposed algorithm is close to that of the FBP algorithm for the reconstruction of the whole object, and the reconstructed image quality is acceptable, which provides an effective method for interior reconstruction in CBCT. Simulation results show the effectiveness of the method.

Statistical Interior Tomography via L1 Norm Dictionary Learning without Assuming an Object Support

Interior tomography of X-ray computed tomography (CT) has many advantages, such as a lower radiation dose and lower detector hardware cost compared to traditional CT. However, this imaging technique only uses the projection data passing through the region of interest (ROI) for imaging; accordingly, the projection data are truncated at both ends of the detector, so the traditional analytical reconstruction algorithm cannot satisfy the demand of clinical diagnosis. To solve the above limitations, in this paper we propose a high-quality statistical iterative reconstruction algorithm that uses the zeroth-order image moment as novel prior knowledge; the zeroth-order image moment can be estimated in the projection domain using the Helgason–Ludwig consistency condition. Then, the L1norm of sparse representation, in terms of dictionary learning, and the zeroth-order image moment constraints are incorporated into the statistical iterative reconstruction framework to construct an objective function. Finally, the objective function is minimized using an alternating minimization iterative algorithm. The chest CT image simulated and CT real data experimental results demonstrate that the proposed approach can remove shift artifacts effectively and has superior performance in removing noise and persevering fine structures than the total variation (TV)-based approach.

Impedance Imaging of Cells and Tissues: Design and Applications

Due to their label-free and noninvasive nature, impedance measurements have attracted increasing interest in biological research. Advances in microfabrication and integrated-circuit technology have opened a route to using large-scale microelectrode arrays for real-time, high-spatiotemporal-resolution impedance measurements of biological samples. In this review, we discuss different methods and applications of measuring impedance for cell and tissue analysis with a focus on impedance imaging with microelectrode arrays in in vitro applications. We first introduce how electrode configurations and the frequency range of the impedance analysis determine the information that can be extracted. We then delve into relevant circuit topologies that can be used to implement impedance measurements and their characteristic features, such as resolution and data-acquisition time. Afterwards, we detail design considerations for the implementation of new impedance-imaging devices. We conclude by discussing future fields of application of impedance imaging in biomedical research, in particular applications where optical imaging is not possible, such as monitoring of ex vivo tissue slices or microelectrode-based brain implants.

Multi source translation based projection completion for interior region of interest tomography with CBCT

Interior tomography by rotary computed tomography (RCT) is an effective method to improve the detection efficiency and achieve high-resolution imaging for the region of interest (ROI) within a large-scale object. However, because only the X-rays through the ROI can be received by detector, the projection data is inevitably truncated, resulting in truncation artifacts in the reconstructed image. When the ROI is totally within the object, the solution of the problem is not unique, which is named interior problem. Fortunately, projection completion (PC) is an effective technique to solve the interior problem. In this study, we proposed a multi source translation CT based PC method (mSTCT-PC) to cope with the interior problem. Firstly, mSTCT-PC employs multi-source translation to sparsely obtain the global projection which covered the whole object. Secondly, the sparse global projection is utilized to fill up the truncated projection of ROI. The global projection and truncated projection are obtained under the same geometric parameters. Therefore, it omits the registration of projection. To verify the feasibility of this method, simulation and practical experiments were implemented. Compared with the results of ROI reconstructed by filtered back-projection (FBP), simultaneous iterative reconstruction technique-total variation (SIRT-TV) and the multi-resolution based method (mR-PC), the proposed mSTCT-PC is good at mitigating truncation artifacts, preserving details and improving the accuracy of ROI images.

Projection-to-image transform frame: a lightweight block reconstruction network (LBRN) for computed tomography

To solve the problem of learning based computed tomography (CT) reconstruction, several reconstruction networks were invented. However, applying neural network to tomographic reconstruction still remains challenging due to unacceptable memory space requirement. In this study, we presents a novel lightweight block reconstruction network (LBRN), which transforms the reconstruction operator into a deep neural network by unrolling the filter back-projection (FBP) method. Specifically, the proposed network contains two main modules, which, respectively, correspond to the filter and back-projection of FBP method. The first module of LBRN decouples the relationship of Radon transform between the reconstructed image and the projection data. Therefore, the following module, block back-projection module, can use the block reconstruction strategy. Due to each image block is only connected with part filtered projection data, the network structure is greatly simplified and the parameters of the whole network is dramatically reduced. Moreover, this approach is trained end-to-end, working directly from raw projection data and does not depend on any initial images. Five reconstruction experiments are conducted to evaluate the performance of the proposed LBRN: full angle, low-dose CT, region of interest (ROI), metal artifacts reduction and real data experiment. The results of the experiments show that the LBRN can be effectively introduced into the reconstruction process and has outstanding advantages in terms of different reconstruction problems.

Data Extrapolation From Learned Prior Images for Truncation Correction in Computed Tomography

Data truncation is a common problem in computed tomography (CT). Truncation causes cupping artifacts inside the field-of-view (FOV) and anatomical structures missing outside the FOV. Deep learning has achieved impressive results in CT reconstruction from limited data. However, its robustness is still a concern for clinical applications. Although the image quality of learning-based compensation schemes may be inadequate for clinical diagnosis, they can provide prior information for more accurate extrapolation than conventional heuristic extrapolation methods. With extrapolated projection, a conventional image reconstruction algorithm can be applied to obtain a final reconstruction. In this work, a general plug-and-play (PnP) method for truncation correction is proposed based on this idea, where various deep learning methods and conventional reconstruction algorithms can be plugged in. Such a PnP method integrates data consistency for measured data and learned prior image information for truncated data. This shows to have better robustness and interpretability than deep learning only. To demonstrate the efficacy of the proposed PnP method, two state-of-the-art deep learning methods, FBPConvNet and Pix2pixGAN, are investigated for truncation correction in cone-beam CT in noise-free and noisy cases. Their robustness is evaluated by showing false negative and false positive lesion cases. With our proposed PnP method, false lesion structures are corrected for both deep learning methods. For FBPConvNet, the root-mean-square error (RMSE) inside the FOV can be improved from 92HU to around 30HU by PnP in the noisy case. Pix2pixGAN solely achieves better image quality than FBPConvNet solely for truncation correction in general. PnP further improves the RMSE inside the FOV from 42HU to around 27HU for Pix2pixGAN. The efficacy of PnP is also demonstrated on real clinical head data.

Generative adversarial networks improve interior computed tomography angiography reconstruction

In interior computed tomography (CT), the x-ray beam is collimated to a limited field-of-view (FOV) (e.g. the volume of the heart) to decrease exposure to adjacent organs, but the resulting image has a severe truncation artifact when reconstructed with traditional filtered back-projection (FBP) type algorithms. In some examinations, such as cardiac or dentomaxillofacial imaging, interior CT could be used to achieve further dose reductions. In this work, we describe a deep learning (DL) method to obtain artifact-free images from interior CT angiography. Our method employs the Pix2Pix generative adversarial network (GAN) in a two-stage process: (1) An extended sinogram is computed from a truncated sinogram with one GAN model, and (2) the FBP reconstruction obtained from that extended sinogram is used as an input to another GAN model that improves the quality of the interior reconstruction. Our double GAN (DGAN) model was trained with 10 000 truncated sinograms simulated from real computed tomography angiography slice images. Truncated sinograms (input) were used with original slice images (target) in training to yield an improved reconstruction (output). DGAN performance was compared with the adaptive de-truncation method, total variation regularization, and two reference DL methods: FBPConvNet, and U-Net-based sinogram extension (ES-UNet). Our DGAN method and ES-UNet yielded the best root-mean-squared error (RMSE) (0.03 ± 0.01), and structural similarity index (SSIM) (0.92 ± 0.02) values, and reference DL methods also yielded good results. Furthermore, we performed an extended FOV analysis by increasing the reconstruction area by 10% and 20%. In both cases, the DGAN approach yielded best results at RMSE (0.03 ± 0.01 and 0.04 ± 0.01 for the 10% and 20% cases, respectively), peak signal-to-noise ratio (PSNR) (30.5 ± 2.6 dB and 28.6 ± 2.6 dB), and SSIM (0.90 ± 0.02 and 0.87 ± 0.02). In conclusion, our method was able to not only reconstruct the interior region with improved image quality, but also extend the reconstructed FOV by 20%.

Reducing axial truncation artifacts in iterative cone-beam CT for radiation therapy using a priori preconditioned information

PURPOSE

Cone-beam computed tomography (CBCT) is increasingly utilized in radiation therapy for image guidance and adaptive applications. While iterative reconstruction algorithms have been shown to outperform traditional filtered back-projection methods in improving image quality and reducing imaging dose, they cannot handle data truncation in the axial view, which frequently occurs in the full-fan partial-trajectory acquisition mode. This proof-of-concept study presents a novel approach on truncation artifact reduction by utilizing a priori preconditioned information as the initial input for the iterative algorithm.

METHODS

Projections containing axial truncation were used for image reconstruction in extended axial field-of-view (AFOV) using the conjugate gradient least-squares (CGLS) algorithm. A priori information in the form of a planning fan-beam CT (FBCT) was repositioned in the expected CBCT imaging geometry, then further processed to dampen high-density features and convolved with a cubic Gaussian kernel to ensure differentiability for the gradient descent method. Anatomical and positional differences between the estimated and the actual imaging object were introduced to verify the efficacy of the proposed method.

RESULTS

Extending the reconstruction AFOV alone could partially reduce truncation artifact. Using a priori information directly resulted in ghosting artifact when there were anatomical and positional differences between the estimated and the actual imaging object. Using a priori preconditioned information was shown to effectively reduce truncation artifact and recover peripheral information.

CONCLUSIONS

Using a priori preconditioned information can effectively alleviate truncation artifact and assist recovery of peripheral information in iterative CBCT reconstruction.

Development of X-ray phase CT with a hybrid configuration of Lau and Talbot-Lau interferometers

X-ray phase computed tomography (CT) is used to observe the inside of light materials. In this paper, we report a new study to develop and test a laboratory assembled X-ray phase CT system that comprises an X-ray Lau interferometer, a rotating Mo anode X-ray tube, and a detector with high spatial resolution. The system has a high spatial resolution lower than 10 μm, which is evaluated by differentiating neighbouring carbon fibres in a polymer composite material. The density resolution is approximately 0.035 g/cm3, which enables to successfully distinguish the high-density polyethylene (HDPE, 0.93 g/cm3) from the ultra-low-density polyethylene (ULDPE, 0.88 g/cm3) in the sample. Moreover, the system can be switched to operate on another mode based on a Talbot-Lau interferometer that provides a wider field of view with a moderate spatial resolution (approximately 100 μm). By analyzing sample images of the biological, this study demonstrates the feasibility and advantages of using hybrid configuration of this X-ray phase CT system.

Interior photon counting computed tomography for quantification of coronary artery calcium: pre-clinical phantom study

Computed tomography (CT) is the reference method for cardiac imaging, but concerns have been raised regarding the radiation dose of CT examinations. Recently, photon counting detectors (PCDs) and interior tomography, in which the radiation beam is limited to the organ-of-interest, have been suggested for patient dose reduction. In this study, we investigated interior PCD-CT (iPCD-CT) for non-enhanced quantification of coronary artery calcium (CAC) using an anthropomorphic torso phantom and ex vivo coronary artery samples. We reconstructed the iPCD-CT measurements with filtered back projection (FBP), iterative total variation (TV) regularization, padded FBP, and adaptively detruncated FBP and adaptively detruncated TV. We compared the organ doses between conventional CT and iPCD-CT geometries, assessed the truncation and cupping artifacts with iPCD-CT, and evaluated the CAC quantification performance of iPCD-CT. With approximately the same effective dose between conventional CT geometry (0.30 mSv) and interior PCD-CT with 10.2 cm field-of-view (0.27 mSv), the organ dose of the heart was increased by 52.3% with interior PCD-CT when compared to CT. Conversely, the organ doses to peripheral and radiosensitive organs, such as the stomach (55.0% reduction), were often reduced with interior PCD-CT. FBP and TV did not sufficiently reduce the truncation artifact, whereas padded FBP and adaptively detruncated FBP and TV yielded satisfactory truncation artifact reduction. Notably, the adaptive detruncation algorithm reduced truncation artifacts effectively when it was combined with reconstruction detrending. With this approach, the CAC quantification accuracy was good, and the coronary artery disease grade reclassification rate was particularly low (5.6%). Thus, our results confirm that CAC quantification can be performed with the interior CT geometry, that the artifacts are effectively reduced with suitable interior reconstruction methods, and that interior tomography provides efficient patient dose reduction.

Diagonalization of the finite Hilbert transform on two adjacent intervals: the Riemann-Hilbert approach

In this paper we study the spectra of bounded self-adjoint linear operators that are related to finite Hilbert transformsH L : L 2 ( [ b L , 0 ] ) → L 2 ( [ 0 , b R ] ) andH R : L 2 ( [ 0 , b R ] ) → L 2 ( [ b L , 0 ] ) . These operators arise when one studies the interior problem of tomography. The diagonalization ofH R , H L has been previously obtained, but only asymptotically whenb L ≠ - b R . We implement a novel approach based on the method of matrix Riemann-Hilbert problems (RHP) which diagonalizesH R , H L explicitly. We also find the asymptotics of the solution to a related RHP and obtain error estimates.

Learning to Reconstruct Computed Tomography Images Directly From Sinogram Data Under A Variety of Data Acquisition Conditions

Computed tomography (CT) is widely used in medical diagnosis and non-destructive detection. Image reconstruction in CT aims to accurately recover pixel values from measured line integrals, i.e., the summed pixel values along straight lines. Provided that the acquired data satisfy the data sufficiency condition as well as other conditions regarding the view angle sampling interval and the severity of transverse data truncation, researchers have discovered many solutions to accurately reconstruct the image. However, if these conditions are violated, accurate image reconstruction from line integrals remains an intellectual challenge. In this paper, a deep learning method with a common network architecture, termed iCT-Net, was developed and trained to accurately reconstruct images for previously solved and unsolved CT reconstruction problems with high quantitative accuracy. Particularly, accurate reconstructions were achieved for the case when the sparse view reconstruction problem (i.e., compressed sensing problem) is entangled with the classical interior tomographic problems.

Efficient solving algorithm for determining the exact sampling condition of limited-angle computed tomography reconstruction

Total variation (TV) regularization-based iterative reconstruction algorithms have an impressive potential to solve limited-angle computed tomography with insufficient sampling projections. The analysis of exact reconstruction sampling conditions for a TV-minimization reconstruction model can determine the minimum number of scanning angle and minimize the scanning range. However, the large-scale matrix operations caused by increased testing phantom size are the computation bottleneck in determining the exact reconstruction sampling conditions in practice. When the size of the testing phantom increases to a certain scale, it is very difficult to analyze quantitatively the exact reconstruction sampling condition using existing methods. In this paper, we propose a fast and efficient algorithm to determine the exact reconstruction sampling condition for large phantoms. Specifically, the sampling condition of a TV minimization model is modeled as a convex optimization problem, which is derived from the sufficient and necessary condition of solution uniqueness for the L1 minimization model. An effective alternating direction minimization algorithm is developed to optimize the objective function by alternatively solving two sub-problems split from the convex problem. The Cholesky decomposition method is used in solving the first sub-problem to reduce computational complexity. Experimental results show that the proposed method can efficiently solve the verification problem of the accurate reconstruction sampling condition. Furthermore, we obtain the lower bounds of scanning angle range for the exact reconstruction of a specific phantom with the larger size.

K-edge-based interior tomography

Interior tomography reconstructs a region of interest using truncated projection data, but it is subject to the ill-posedness of interior tomography. With the photon-counting detector, K-edge imaging uses data in the low and high energy bins around the K-edge of a contrast agent, and can faithfully recover true image contrast for improved diagnosis. The purpose of this paper is to reconstruct a region of interest inside a patient assuming the existence of a known K-edge material. In this case, there is a significant difference in x-ray attenuation around the K-edge, but these attenuation coefficients are inter-related to guide updating an intermediate reconstruction until a stopping criterion is satisfied. In our study, new interior tomography algorithms were developed without any major computational overhead, and several phantoms were used to validate the algorithms. The proposed methods are advantageous relative to the existing interior tomography algorithms, because of the available spectral information in the form of a known K-edge material.

Quantitative region-of-interest tomography using variable field of view

In X-ray computed tomography, the task of imaging only a local region of interest (ROI) inside a larger sample is very important. However, without a priori information, this ROI cannot be exactly reconstructed using only the image data limited to the ROI. We propose here an approach of region-of-interest tomography, which reconstructs a ROI within an object from projections of different fields of view acquired on a specific angular sampling scheme in the same tomographic experiment. We present a stable procedure that not only yields high-quality images of the ROI but keeps as well the quantitative contrast on the reconstructed images. In addition, we analyze the minimum number of projections required for ROI tomography from the point of view of the band region of the Radon transform, which confirms this number must be estimated based on the size of the entire object and not only on the size of the ROI.

Optimization-based region-of-interest reconstruction for X-ray computed tomography based on total variation and data derivative

Region-of-interest (ROI) and interior reconstructions for computed tomography (CT) have drawn much attention and can be of practical value for potential applications in reducing radiation dose and hardware cost. The conventional wisdom is that the exact reconstruction of an interior ROI is very difficult to be obtained by only using data associated with lines through the ROI. In this study, we propose and investigate optimization-based methods for ROI and interior reconstructions based on total variation (TV) and data derivative. Objective functions are built by the image TV term plus the data finite difference term. Different data terms in the forms of L1-norm, L2-norm, and Kullback-Leibler divergence are incorporated and investigated in the optimizations. Efficient algorithms are developed using the proximal alternating direction method of multipliers (ADMM) for each program. All sub-problems of ADMM are solved by using closed-form solutions with high efficiency. The customized optimizations and algorithms based on the TV and derivative-based data terms can serve as a powerful tool for interior reconstructions. Simulations and real-data experiments indicate that the proposed methods can be of practical value for CT imaging applications.

Dual-rotation C-arm cone-beam computed tomography to increase low-contrast detection

PURPOSE

This paper investigates the capabilities of a dual-rotation C-arm cone-beam computed tomography (CBCT) framework to improve non-contrast-enhanced low-contrast detection for full volume or volume-of-interest (VOI) brain imaging.

METHOD

The idea is to associate two C-arm short-scan rotational acquisitions (spins): one over the full detector field of view (FOV) at low dose, and one collimated to deliver a higher dose to the central densest parts of the head. The angular sampling performed by each spin is allowed to vary in terms of number of views and angular positions. Collimated data is truncated and does not contain measurement of the incoming X-ray intensities in air (air calibration). When targeting full volume reconstruction, the method is intended to act as a virtual bow-tie. When targeting VOI imaging, the method is intended to provide the minimum full detector FOV data that sufficiently corrects for truncation artifacts. A single dedicated iterative algorithm is described that handles all proposed sampling configurations despite truncation and absence of air calibration.

RESULTS

Full volume reconstruction of dual-rotation simulations and phantom acquisitions are shown to have increased low-contrast detection for less dose, with respect to a single-rotation acquisition. High CNR values were obtained on 1% inserts of the Catphan® 515 module in 0.94 mm thick slices. Image quality for VOI imaging was preserved from truncation artifacts even with less than 10 non-truncated views, without using the sparsity a priori common to such context.

CONCLUSION

A flexible dual-rotation acquisition and reconstruction framework is proposed that has the potential to improve low-contrast detection in clinical C-arm brain soft-tissue imaging.

Quest for the ultimate cardiac CT scanner

PURPOSE

To quantitatively evaluate and compare six proposed system architectures for cardiac CT scanning.

METHODS

Starting from the clinical requirements for cardiac CT, we defined six dedicated cardiac CT architectures. We selected these architectures based on a previous screening study and defined them in sufficient detail to comprehensively analyze their cost and performance. We developed rigorous comparative evaluation methods for the most important aspects of performance and cost, and we applied these evaluation criteria to the defined cardiac CT architectures.

RESULTS

We found that CT system architectures based on the third-generation geometry provide nearly linear performance improvement versus the increased cost of additional beam lines (i.e., source-detector pairs), although similar performance improvement could be achieved with advanced motion-correction algorithms. The third-generation architectures outperform even the most promising of the proposed architectures that deviate substantially from the traditional CT system architectures.

CONCLUSION

This work confirms the validity of the current trend in commercial CT scanner design. However, we anticipate that over time, CT hardware and software technologies will evolve, the relative importance of the performance criteria will change, the relative costs of components will vary, some of the remaining challenges will be addressed, and perhaps new candidate architectures will be identified; therefore, the conclusion of a comparative analysis like this may change. The evaluation methods that we used can provide a framework for other researchers to analyze their own proposed CT architectures.

Initial analysis of the middle problem in CT image reconstruction

The interior and exterior problems have been extensively studied in the field of reconstruction of computed tomography (CT) images, which lead to important theoretical and practical results. In this study, we formulate a middle problem of CT image reconstruction, which is more challenging than either the interior or exterior problems. In the middle problem of CT image reconstruction, projection data are measured through and only through the middle dough-like region, so that each projection profile misses data not only internally but also on both sides. For an object with a radially symmetric exterior, we proved that the middle problem could be uniquely solved if the middle ring-shaped zone is piecewise constant or there is a known sub-region inside this middle region. Then, we designed and evaluated a POCS-based algorithm for middle tomography, which is to reconstruct a middle image only from the available data. Finally, the remaining issues are also discussed for further research.

Efficient implementation of a local tomography reconstruction algorithm

We propose an efficient implementation of an interior tomography reconstruction method based on a known subregion. This method iteratively refines a reconstruction, aiming at reducing the local tomography artifacts. To cope with the ever increasing data volumes, this method is highly optimized on two aspects: firstly, the problem is reformulated to reduce the number of variables, and secondly, the operators involved in the optimization algorithms are efficiently implemented. Results show that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4096^2$$\end{document}40962 slices can be processed in tens of seconds, while being beyond the reach of equivalent exact local tomography method.

Singular value decomposition-based 2D image reconstruction for computed tomography

Singular value decomposition (SVD)-based 2D image reconstruction methods are developed and evaluated for a broad class of inverse problems for which there are no analytical solutions. The proposed methods are fast and accurate for reconstructing images in a non-iterative fashion. The multi-resolution strategy is adopted to reduce the size of the system matrix to reconstruct large images using limited memory capacity. A modified high-contrast Shepp-Logan phantom, a low-contrast FORBILD head phantom, and a physical phantom are employed to evaluate the proposed methods with different system configurations. The results show that the SVD methods can accurately reconstruct images from standard scan and interior scan projections and that they outperform other benchmark methods. The general SVD method outperforms the other SVD methods. The truncated SVD and Tikhonov regularized SVD methods accurately reconstruct a region-of-interest (ROI) from an internal scan with a known sub-region inside the ROI. Furthermore, the SVD methods are much faster and more flexible than the benchmark algorithms, especially in the ROI reconstructions in our experiments.

Interior tomography with curvelet-based regularization

The interior problem, i.e. reconstruction from local truncated projections in computed tomography (CT), is common in practical applications. However, its solution is non-unique in a general unconstrained setting. To solve the interior problem uniquely and stably, in recent years both the prior knowledge- and compressive sensing (CS)-based methods have been developed. Those theoretically exact solutions for the interior problem are called interior tomography. Along this direction, we propose here a new CS-based method for the interior problem based on the curvelet transform. A curvelet is localized in both radial and angular directions in the frequency domain. A two-dimensional (2D) image can be represented in a curvelet frame. We employ the curvelet transform coefficients to regularize the interior problem and obtain a curvelet frame based regularization method (CFRM) for interior tomography. The curvelet coefficients of the reconstructed image are split into two sets according to their visibility from the interior data, and different regularization parameters are used for these two sets. We also presents the results of numerical experiments, which demonstrate the feasibility of the proposed approach.

A practical local tomography reconstruction algorithm based on a known sub-region

A new method to reconstruct data acquired in a local tomography setup is proposed. This method uses an initial reconstruction and refines it by correcting the low-frequency artifacts, known as the cupping effect. A basis of Gaussian functions is used to correct the initial reconstruction. The coefficients of this basis are found by optimizing iteratively a fidelity term under the constraint of a known sub-region. Using a coarse basis reduces the degrees of freedom of the problem while actually correcting the cupping effect. Simulations show that the known region constraint yields an unbiased reconstruction, in accordance with uniqueness theorems stated in local tomography.

Fast iterative reconstruction of data in full interior tomography

This paper introduces two novel strategies for iterative reconstruction of full interior tomography (FINT) data, i.e. when the field of view is entirely inside the object support and knowledge of the object support itself or the attenuation coefficients inside specific regions of interest are not available. The first approach is based on data edge-padding. The second technique creates an intermediate virtual sinogram, which is, then, reconstructed by a standard iterative algorithm. Both strategies are validated in the framework of the alternate direction method of multipliers plug-and-play with gridding projectors that provide a speed-up of three orders of magnitude with respect to standard operators implemented in real space. The proposed methods are benchmarked on synchrotron-based X-ray tomographic microscopy datasets of mouse lung alveoli. Compared with analytical techniques, the proposed methods substantially improve the reconstruction quality for FINT underconstrained datasets, facilitating subsequent post-processing steps.

Multi-resolution statistical image reconstruction for mitigation of truncation effects: application to cone-beam CT of the head

A prototype cone-beam CT (CBCT) head scanner featuring model-based iterative reconstruction (MBIR) has been recently developed and demonstrated the potential for reliable detection of acute intracranial hemorrhage (ICH), which is vital to diagnosis of traumatic brain injury and hemorrhagic stroke. However, data truncation (e.g. due to the head holder) can result in artifacts that reduce image uniformity and challenge ICH detection. We propose a multi-resolution MBIR method with an extended reconstruction field of view (RFOV) to mitigate truncation effects in CBCT of the head. The image volume includes a fine voxel size in the (inner) nontruncated region and a coarse voxel size in the (outer) truncated region. This multi-resolution scheme allows extension of the RFOV to mitigate truncation effects while introducing minimal increase in computational complexity. The multi-resolution method was incorporated in a penalized weighted least-squares (PWLS) reconstruction framework previously developed for CBCT of the head. Experiments involving an anthropomorphic head phantom with truncation due to a carbon-fiber holder were shown to result in severe artifacts in conventional single-resolution PWLS, whereas extending the RFOV within the multi-resolution framework strongly reduced truncation artifacts. For the same extended RFOV, the multi-resolution approach reduced computation time compared to the single-resolution approach (viz. time reduced by 40.7%, 83.0%, and over 95% for an image volume of 600³, 800³, 1000³ voxels). Algorithm parameters (e.g. regularization strength, the ratio of the fine and coarse voxel size, and RFOV size) were investigated to guide reliable parameter selection. The findings provide a promising method for truncation artifact reduction in CBCT and may be useful for other MBIR methods and applications for which truncation is a challenge.

An adaptive reconstruction algorithm for spectral CT regularized by a reference image

The photon counting detector based spectral CT system is attracting increasing attention in the CT field. However, the spectral CT is still premature in terms of both hardware and software. To reconstruct high quality spectral images from low-dose projections, an adaptive image reconstruction algorithm is proposed that assumes a known reference image (RI). The idea is motivated by the fact that the reconstructed images from different spectral channels are highly correlated. If a high quality image of the same object is known, it can be used to improve the low-dose reconstruction of each individual channel. This is implemented by maximizing the patch-wise correlation between the object image and the RI. Extensive numerical simulations and preclinical mouse study demonstrate the feasibility and merits of the proposed algorithm. It also performs well for truncated local projections, and the surrounding area of the region- of-interest (ROI) can be more accurately reconstructed. Furthermore, a method is introduced to adaptively choose the step length, making the algorithm more feasible and easier for applications.

A spectral interior CT by a framelet-based reconstruction algorithm

Reducing radiation dose is an important goal in medical computed tomography (CT), for which interior tomography is an effective approach. There have been interior reconstruction algorithms for monochromatic CT, but in reality, X-ray sources are polychromatic. Using a polychromatic acquisition model and motivated by framelet-based image processing algorithms, in this paper, we propose an interior reconstruction algorithm to obtain an image with spectral information assuming only one scan with a current energy-integrating detector. This algorithm is a new nonlinear iterative method by minimizing a special functional under a polychromatic acquisition model for X-ray CT, where the attenuation coefficients are energy-dependent. Experimental results validate that our algorithm can effectively reduce the beam-hardening artifacts and metal artifacts. It also produces color overlays which are useful in tumor identification and quantification.

A two-step filtering-based iterative image reconstruction method for interior tomography

The optimization-based method that utilizes the additional sparse prior of region-of-interest (ROI) image, such as total variation, has been the subject of considerable research in problems of interior tomography reconstruction. One challenge for optimization-based iterative ROI image reconstruction is to build the relationship between ROI image and truncated projection data. When the reconstruction support region is smaller than the original object, an unsuitable representation of data fidelity may lead to bright truncation artifacts in the boundary region of field of view. In this work, we aim to develop an iterative reconstruction method to suppress the truncation artifacts and improve the image quality for direct ROI image reconstruction. A novel reconstruction approach is proposed based on an optimization problem involving a two-step filtering-based data fidelity. Data filtering is achieved in two steps: the first takes the derivative of projection data; in the second step, Hilbert filtering is applied in the differentiated data. Numerical simulations and real data reconstructions have been conducted to validate the new reconstruction method. Both qualitative and quantitative results indicate that, as theoretically expected, the proposed method brings reasonable performance in suppressing truncation artifacts and preserving detailed features. The presented local reconstruction method based on the two-step filtering strategy provides a simple and efficient approach for the iterative reconstruction from truncated projections.

Vision 20/20: Simultaneous CT-MRI--Next chapter of multimodality imaging

Multimodality imaging systems such as positron emission tomography-computed tomography (PET-CT) and MRI-PET are widely available, but a simultaneous CT-MRI instrument has not been developed. Synergies between independent modalities, e.g., CT, MRI, and PET/SPECT can be realized with image registration, but such postprocessing suffers from registration errors that can be avoided with synchronized data acquisition. The clinical potential of simultaneous CT-MRI is significant, especially in cardiovascular and oncologic applications where studies of the vulnerable plaque, response to cancer therapy, and kinetic and dynamic mechanisms of targeted agents are limited by current imaging technologies. The rationale, feasibility, and realization of simultaneous CT-MRI are described in this perspective paper. The enabling technologies include interior tomography, unique gantry designs, open magnet and RF sequences, and source and detector adaptation. Based on the experience with PET-CT, PET-MRI, and MRI-LINAC instrumentation where hardware innovation and performance optimization were instrumental to construct commercial systems, the authors provide top-level concepts for simultaneous CT-MRI to meet clinical requirements and new challenges. Simultaneous CT-MRI fills a major gap of modality coupling and represents a key step toward the so-called "omnitomography" defined as the integration of all relevant imaging modalities for systems biology and precision medicine.

Interior tomographic imaging of mouse heart in a carbon nanotube micro-CT

BACKGROUND

The relatively high radiation dose from micro-CT is a cause for concern in preclinical research involving animal subjects. Interior region-of-interest (ROI) imaging was proposed for dose reduction, but has not been experimentally applied in micro-CT.

OBJECTIVE

Our aim is to implement interior ROI imaging in a carbon nanotube (CNT) x-ray source based micro-CT, and present the ROI image quality and radiation dose reduction for interior cardiac micro-CT imaging of a mouse heart in situ.

METHODS

An aperture collimator was mounted at the source-side to induce a small-sized cone beam (10 mm width) at the isocenter. Interior in situ micro-CT scans were conducted on a mouse carcass and several micro-CT phantoms. A GPU-accelerated hybrid iterative reconstruction algorithm was employed for volumetric image reconstruction. Radiation dose was measured for the same system operated at the interior and global micro-CT modes.

RESULTS

Visual inspection demonstrated comparable image quality between two scan modes. Quantitative evaluation demonstrated high structural similarity index (up to 0.9614) with improved contrast-noise-ratio (CNR) on interior micro-CT mode. Interior micro-CT mode yielded significant reduction (up to 83.9%) for dose length product (DLP).

CONCLUSIONS

This work demonstrates the applicability of using CNT x-ray source based interior micro-CT for preclinical imaging with significantly reduced radiation dose.

Accelerated statistical reconstruction for C-arm cone-beam CT using Nesterov's method

PURPOSE

To accelerate model-based iterative reconstruction (IR) methods for C-arm cone-beam CT (CBCT), thereby combining the benefits of improved image quality and/or reduced radiation dose with reconstruction times on the order of minutes rather than hours.

METHODS

The ordered-subsets, separable quadratic surrogates (OS-SQS) algorithm for solving the penalized-likelihood (PL) objective was modified to include Nesterov's method, which utilizes "momentum" from image updates of previous iterations to better inform the current iteration and provide significantly faster convergence. Reconstruction performance of an anthropomorphic head phantom was assessed on a benchtop CBCT system, followed by CBCT on a mobile C-arm, which provided typical levels of incomplete data, including lateral truncation. Additionally, a cadaveric torso that presented realistic soft-tissue and bony anatomy was imaged on the C-arm, and different projectors were assessed for reconstruction speed.

RESULTS

Nesterov's method provided equivalent image quality to OS-SQS while reducing the reconstruction time by an order of magnitude (10.0 ×) by reducing the number of iterations required for convergence. The faster projectors were shown to produce similar levels of convergence as more accurate projectors and reduced the reconstruction time by another 5.3 ×. Despite the slower convergence of IR with truncated C-arm CBCT, comparison of PL reconstruction methods implemented on graphics processing units showed that reconstruction time was reduced from 106 min for the conventional OS-SQS method to as little as 2.0 min with Nesterov's method for a volumetric reconstruction of the head. In body imaging, reconstruction of the larger cadaveric torso was reduced from 159 min down to 3.3 min with Nesterov's method.

CONCLUSIONS

The acceleration achieved through Nesterov's method combined with ordered subsets reduced IR times down to a few minutes. This improved compatibility with clinical workflow better enables broader adoption of IR in CBCT-guided procedures, with corresponding benefits in overcoming conventional limits of image quality at lower dose.

Extended ellipse-line-ellipse trajectory for long-object cone-beam imaging with a mounted C-arm system

Recent reports show that three-dimensional cone-beam (CB) imaging with a floor-mounted (or ceiling-mounted) C-arm system has become a valuable tool in interventional radiology. Currently, a circular short scan is used for data acquisition, which inevitably yields CB artifacts and a short coverage in the direction of the patient table. To overcome these two limitations, a more sophisticated data acquisition geometry is needed. This geometry should be complete in terms of Tuy's condition and should allow continuous scanning, while being compatible with the mechanical constraints of mounted C-arm systems. Additionally, the geometry should allow accurate image reconstruction from truncated data. One way to ensure such a feature is to adopt a trajectory that provides full R-line coverage within the field-of-view (FOV). An R-line is any segment of line that connects two points on a source trajectory, and the R-line coverage is the set of points that belong to an R-line. In this work, we propose a novel geometry called the extended ellipse-line-ellipse (ELE) for long-object imaging with a mounted C-arm system. This trajectory is built from modules consisting of two elliptical arcs connected by a line. We demonstrate that the extended ELE can be configured in many ways so that full R-line coverage is guaranteed. Both tight and relaxed parametric settings are presented. All results are supported by extensive mathematical proofs provided in appendices. Our findings make the extended ELE trajectory attractive for axially-extended FOV imaging in interventional radiology.

Patient-bounded extrapolation using low-dose priors for volume-of-interest imaging in C-arm CT

PURPOSE

Three-dimensional (3D) volume-of-interest (VOI) imaging with C-arm systems provides anatomical information in a predefined 3D target region at a considerably low x-ray dose. However, VOI imaging involves laterally truncated projections from which conventional reconstruction algorithms generally yield images with severe truncation artifacts. Heuristic based extrapolation methods, e.g., water cylinder extrapolation, typically rely on techniques that complete the truncated data by means of a continuity assumption and thus appear to be ad-hoc. It is our goal to improve the image quality of VOI imaging by exploiting existing patient-specific prior information in the workflow.

METHODS

A necessary initial step prior to a 3D acquisition is to isocenter the patient with respect to the target to be scanned. To this end, low-dose fluoroscopic x-ray acquisitions are usually applied from anterior-posterior (AP) and medio-lateral (ML) views. Based on this, the patient is isocentered by repositioning the table. In this work, we present a patient-bounded extrapolation method that makes use of these noncollimated fluoroscopic images to improve image quality in 3D VOI reconstruction. The algorithm first extracts the 2D patient contours from the noncollimated AP and ML fluoroscopic images. These 2D contours are then combined to estimate a volumetric model of the patient. Forward-projecting the shape of the model at the eventually acquired C-arm rotation views gives the patient boundary information in the projection domain. In this manner, we are in the position to substantially improve image quality by enforcing the extrapolated line profiles to end at the known patient boundaries, derived from the 3D shape model estimate.

RESULTS

The proposed method was evaluated on eight clinical datasets with different degrees of truncation. The proposed algorithm achieved a relative root mean square error (rRMSE) of about 1.0% with respect to the reference reconstruction on nontruncated data, even in the presence of severe truncation, compared to a rRMSE of 8.0% when applying a state-of-the-art heuristic extrapolation technique.

CONCLUSIONS

The method we proposed in this paper leads to a major improvement in image quality for 3D C-arm based VOI imaging. It involves no additional radiation when using fluoroscopic images that are acquired during the patient isocentering process. The model estimation can be readily integrated into the existing interventional workflow without additional hardware.

An iterative deconvolution algorithm for image recovery in clinical CT: A phantom study

PURPOSE

To study the feasibility of using an iterative reconstruction algorithm to improve previously reconstructed CT images which are judged to be non-diagnostic on clinical review. A novel rapidly converging, iterative algorithm (RSEMD) to reduce noise as compared with standard filtered back-projection algorithm has been developed.

MATERIALS AND METHODS

The RSEMD method was tested on in-silico, Catphan(®)500, and anthropomorphic 4D XCAT phantoms. The method was applied to noisy CT images previously reconstructed with FBP to determine improvements in SNR and CNR. To test the potential improvement in clinically relevant CT images, 4D XCAT phantom images were used to simulate a small, low contrast lesion placed in the liver.

RESULTS

In all of the phantom studies the images proved to have higher resolution and lower noise as compared with images reconstructed by conventional FBP. In general, the values of SNR and CNR reached a plateau at around 20 iterations with an improvement factor of about 1.5 for in noisy CT images. Improvements in lesion conspicuity after the application of RSEMD have also been demonstrated. The results obtained with the RSEMD method are in agreement with other iterative algorithms employed either in image space or with hybrid reconstruction algorithms.

CONCLUSIONS

In this proof of concept work, a rapidly converging, iterative deconvolution algorithm with a novel resolution subsets-based approach that operates on DICOM CT images has been demonstrated. The RSEMD method can be applied to sub-optimal routine-dose clinical CT images to improve image quality to potentially diagnostically acceptable levels.

A dynamic attenuator improves spectral imaging with energy-discriminating, photon counting detectors

Energy-discriminating, photon counting (EDPC) detectors have high potential in spectral imaging applications but exhibit degraded performance when the incident count rate approaches or exceeds the characteristic count rate of the detector. In order to reduce the requirements on the detector, we explore the strategy of modulating the X-ray flux field using a recently proposed dynamic, piecewise-linear attenuator. A previous paper studied this modulation for photon counting detectors but did not explore the impact on spectral applications. In this work, we modeled detection with a bipolar triangular pulse shape (Taguchi et al., 2011) and estimated the Cramer-Rao lower bound (CRLB) of the variance of material selective and equivalent monoenergetic images, assuming deterministic errors at high flux could be corrected. We compared different materials for the dynamic attenuator and found that rare earth elements, such as erbium, outperformed previously proposed materials such as iron in spectral imaging. The redistribution of flux reduces the variance or dose, consistent with previous studies on benefits with conventional detectors. Numerical simulations based on DICOM datasets were used to assess the impact of the dynamic attenuator for detectors with several different characteristic count rates. The dynamic attenuator reduced the peak incident count rate by a factor of 4 in the thorax and 44 in the pelvis, and a 10 Mcps/mm (2) EDPC detector with dynamic attenuator provided generally superior image quality to a 100 Mcps/mm (2) detector with reference bowtie filter for the same dose. The improvement is more pronounced in the material images.

An algorithm to estimate the object support in truncated images

PURPOSE

Truncation artifacts in CT occur if the object to be imaged extends past the scanner field of view (SFOV). These artifacts impede diagnosis and could possibly introduce errors in dose plans for radiation therapy. Several approaches exist for correcting truncation artifacts, but existing correction algorithms do not accurately recover the skin line (or support) of the patient, which is important in some dose planning methods. The purpose of this paper was to develop an iterative algorithm that recovers the support of the object.

METHODS

The authors assume that the truncated portion of the image is made up of soft tissue of uniform CT number and attempt to find a shape consistent with the measured data. Each known measurement in the sinogram is interpreted as an estimate of missing mass along a line. An initial estimate of the object support is generated by thresholding a reconstruction made using a previous truncation artifact correction algorithm (e.g., water cylinder extrapolation). This object support is iteratively deformed to reduce the inconsistency with the measured data. The missing data are estimated using this object support to complete the dataset. The method was tested on simulated and experimentally truncated CT data.

RESULTS

The proposed algorithm produces a better defined skin line than water cylinder extrapolation. On the experimental data, the RMS error of the skin line is reduced by about 60%. For moderately truncated images, some soft tissue contrast is retained near the SFOV. As the extent of truncation increases, the soft tissue contrast outside the SFOV becomes unusable although the skin line remains clearly defined, and in reformatted images it varies smoothly from slice to slice as expected.

CONCLUSIONS

The support recovery algorithm provides a more accurate estimate of the patient outline than thresholded, basic water cylinder extrapolation, and may be preferred in some radiation therapy applications.

Interior micro-CT with an offset detector

PURPOSE

The size of field-of-view (FOV) of a microcomputed tomography (CT) system can be increased by offsetting the detector. The increased FOV is beneficial in many applications. All prior investigations, however, have been focused to the case in which the increased FOV after offset-detector acquisition can cover the transaxial extent of an object fully. Here, the authors studied a new problem where the FOV of a micro-CT system, although increased after offset-detector acquisition, still covers an interior region-of-interest (ROI) within the object.

METHODS

An interior-ROI-oriented micro-CT scan with an offset detector poses a difficult reconstruction problem, which is caused by both detector offset and projection truncation. Using the projection completion techniques, the authors first extended three previous reconstruction methods from offset-detector micro-CT to offset-detector interior micro-CT. The authors then proposed a novel method which combines two of the extended methods using a frequency split technique. The authors tested the four methods with phantom simulations at 9.4%, 18.8%, 28.2%, and 37.6% detector offset. The authors also applied these methods to physical phantom datasets acquired at the same amounts of detector offset from a customized micro-CT system.

RESULTS

When the detector offset was small, all reconstruction methods showed good image quality. At large detector offset, the three extended methods gave either visible shading artifacts or high deviation of pixel value, while the authors' proposed method demonstrated no visible artifacts and minimal deviation of pixel value in both the numerical simulations and physical experiments.

CONCLUSIONS

For an interior micro-CT with an offset detector, the three extended reconstruction methods can perform well at a small detector offset but show strong artifacts at a large detector offset. When the detector offset is large, the authors' proposed reconstruction method can outperform the three extended reconstruction methods by suppressing artifacts and maintaining pixel values.

Reconstructing cone-beam CT with spatially varying qualities for adaptive radiotherapy: a proof-of-principle study

With the aim of maximally reducing imaging dose while meeting requirements for adaptive radiation therapy (ART), we propose in this paper a new cone beam CT (CBCT) acquisition and reconstruction method that delivers images with a low noise level inside a region of interest (ROI) and a relatively high noise level outside the ROI. The acquired projection images include two groups: densely sampled projections at a low exposure with a large field of view (FOV) and sparsely sampled projections at a high exposure with a small FOV corresponding to the ROI. A new algorithm combining the conventional filtered back-projection algorithm and the tight-frame iterative reconstruction algorithm is also designed to reconstruct the CBCT based on these projection data. We have validated our method on a simulated head-and-neck (HN) patient case, a semi-real experiment conducted on a HN cancer patient under a full-fan scan mode, as well as a Catphan phantom under a half-fan scan mode. Relative root-mean-square errors (RRMSEs) of less than 3% for the entire image and ~1% within the ROI compared to the ground truth have been observed. These numbers demonstrate the ability of our proposed method to reconstruct high-quality images inside the ROI. As for the part outside ROI, although the images are relatively noisy, it can still provide sufficient information for radiation dose calculations in ART. Dose distributions calculated on our CBCT image and on a standard CBCT image are in agreement, with a mean relative difference of 0.082% inside the ROI and 0.038% outside the ROI. Compared with the standard clinical CBCT scheme, an imaging dose reduction of approximately 3-6 times inside the ROI was achieved, as well as an 8 times outside the ROI. Regarding computational efficiency, it takes 1-3 min to reconstruct a CBCT image depending on the number of projections used. These results indicate that the proposed method has the potential for application in ART.

Prior-based artifact correction (PBAC) in computed tomography

PURPOSE

Image quality in computed tomography (CT) often suffers from artifacts which may reduce the diagnostic value of the image. In many cases, these artifacts result from missing or corrupt regions in the projection data, e.g., in the case of metal, truncation, and limited angle artifacts. The authors propose a generalized correction method for different kinds of artifacts resulting from missing or corrupt data by making use of available prior knowledge to perform data completion.

METHODS

The proposed prior-based artifact correction (PBAC) method requires prior knowledge in form of a planning CT of the same patient or in form of a CT scan of a different patient showing the same body region. In both cases, the prior image is registered to the patient image using a deformable transformation. The registered prior is forward projected and data completion of the patient projections is performed using smooth sinogram inpainting. The obtained projection data are used to reconstruct the corrected image.

RESULTS

The authors investigate metal and truncation artifacts in patient data sets acquired with a clinical CT and limited angle artifacts in an anthropomorphic head phantom data set acquired with a gantry-based flat detector CT device. In all cases, the corrected images obtained by PBAC are nearly artifact-free. Compared to conventional correction methods, PBAC achieves better artifact suppression while preserving the patient-specific anatomy at the same time. Further, the authors show that prominent anatomical details in the prior image seem to have only minor impact on the correction result.

CONCLUSIONS

The results show that PBAC has the potential to effectively correct for metal, truncation, and limited angle artifacts if adequate prior data are available. Since the proposed method makes use of a generalized algorithm, PBAC may also be applicable to other artifacts resulting from missing or corrupt data.

An OpenPET scanner with bridged detectors to compensate for incomplete data

We are developing an open-type PET 'OpenPET' geometry. One possible geometry is a dual-ring OpenPET, which consists of two detector rings separated by a gap for entrance of a radiotherapy beam or for inserting other modalities. In our previous simulations and experiments the OpenPET imaging geometry was shown to be feasible by applying iterative reconstruction methods. However, the gap violates Orlov's completeness condition for accurate tomographic reconstruction. In this study, we propose a solution for the incompleteness problem by adding bridge detectors to fill in parts of the gaps of the OpenPET geometry; we call this bridged OpenPET. Although this geometry was considered previously, its analytical property was not discussed. Therefore, we applied the direct Fourier method as an analytical reconstruction method to the bridged OpenPET, dual-ring OpenPET and conventional cylindrical PET for comparison. Numerical simulations showed that the additional bridge detectors compensate for the incompleteness of the OpenPET by covering one direction perpendicular to the transaxial slices of the imaging subjects.

Impact of extraneous mispositioned events on motion-corrected brain SPECT images of freely moving animals

PURPOSE

Single photon emission computed tomography (SPECT) brain imaging of freely moving small animals would allow a wide range of important neurological processes and behaviors to be studied, which are normally inhibited by anesthetic drugs or precluded due to the animal being restrained. While rigid body motion of the head can be tracked and accounted for in the reconstruction, activity in the torso may confound brain measurements, especially since motion of the torso is more complex (i.e., nonrigid) and not well correlated with that of the head. The authors investigated the impact of mispositioned events and attenuation due to the torso on the accuracy of motion corrected brain images of freely moving mice.

METHODS

Monte Carlo simulations of a realistic voxelized mouse phantom and a dual compartment phantom were performed. Each phantom comprised a target and an extraneous compartment which were able to move independently of each other. Motion correction was performed based on the known motion of the target compartment only. Two SPECT camera geometries were investigated: a rotating single head detector and a stationary full ring detector. The effects of motion, detector geometry, and energy of the emitted photons (hence, attenuation) on bias and noise in reconstructed brain regions were evaluated.

RESULTS

The authors observed two main sources of bias: (a) motion-related inconsistencies in the projection data and (b) the mismatch between attenuation and emission. Both effects are caused by the assumption that the orientation of the torso is difficult to track and model, and therefore cannot be conveniently corrected for. The motion induced bias in some regions was up to 12% when no attenuation effects were considered, while it reached 40% when also combined with attenuation related inconsistencies. The detector geometry (i.e., rotating vs full ring) has a big impact on the accuracy of the reconstructed images, with the full ring detector being more advantageous.

CONCLUSIONS

Motion-induced inconsistencies in the projection data and attenuation/emission mismatch are the two main causes of bias in reconstructed brain images when there is complex motion. It appears that these two factors have a synergistic effect on the qualitative and quantitative accuracy of the reconstructed images.

The piecewise-linear dynamic attenuator reduces the impact of count rate loss with photon-counting detectors

Photon counting x-ray detectors (PCXDs) offer several advantages compared to standard energy-integrating x-ray detectors, but also face significant challenges. One key challenge is the high count rates required in CT. At high count rates, PCXDs exhibit count rate loss and show reduced detective quantum efficiency in signal-rich (or high flux) measurements. In order to reduce count rate requirements, a dynamic beam-shaping filter can be used to redistribute flux incident on the patient. We study the piecewise-linear attenuator in conjunction with PCXDs without energy discrimination capabilities. We examined three detector models: the classic nonparalyzable and paralyzable detector models, and a 'hybrid' detector model which is a weighted average of the two which approximates an existing, real detector (Taguchi et al 2011 Med. Phys. 38 1089-102). We derive analytic expressions for the variance of the CT measurements for these detectors. These expressions are used with raw data estimated from DICOM image files of an abdomen and a thorax to estimate variance in reconstructed images for both the dynamic attenuator and a static beam-shaping ('bowtie') filter. By redistributing flux, the dynamic attenuator reduces dose by 40% without increasing peak variance for the ideal detector. For non-ideal PCXDs, the impact of count rate loss is also reduced. The nonparalyzable detector shows little impact from count rate loss, but with the paralyzable model, count rate loss leads to noise streaks that can be controlled with the dynamic attenuator. With the hybrid model, the characteristic count rates required before noise streaks dominate the reconstruction are reduced by a factor of 2 to 3. We conclude that the piecewise-linear attenuator can reduce the count rate requirements of the PCXD in addition to improving dose efficiency. The magnitude of this reduction depends on the detector, with paralyzable detectors showing much greater benefit than nonparalyzable detectors.