Interior region-of-interest reconstruction using a small, nearly piecewise constant subregion

Accurate image reconstruction within and beyond the field-of-view of CT system from data with truncation

Objective. Accurate image reconstruction from data with truncation in x-ray computed tomography (CT) remains a topic of research interest; and the works reported previously in the literature focus largely on reconstructing an image only within the scanning field-of-view (FOV). We develop algorithms to invert the truncated data model for numerically accurate image reconstruction within the subject support or a region slightly smaller than the subject support.Methods. We formulate image reconstruction from data with truncation as an optimization program, which includes hybrid constraints on region-based image total-variation (TV) and imageℓ1-norm (L1) for effectively suppressing truncation artifacts. An algorithm, referred to as the TV-L1 algorithm, is developed for image reconstruction (i.e. inversion of the data model) from data with truncation through solving the optimization program.Results. We perform numerical studies to evaluate accuracy and stability of the TV-L1 algorithm by using simulated and real CT data. Accurate images can be obtained stably by use of the TV-L1 algorithm within the subject support, or a region substantially larger than the FOV, from data with truncation of varying degrees.Conclusions. The TV-L1 algorithm can invert the truncated data model to accurately and stably reconstruct images within the subject support, or a region slightly smaller than the subject support but substantially larger than the FOV.Significance. Accurate image reconstruction within the subject support, or a region substantially larger than the FOV, from data with truncation can be of theoretical and practical implication. The insights and TV-L1 algorithm may also be generalized to accurate image reconstruction from data with truncation in other tomographic imaging modalities.

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.

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.

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%.

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.

Overcoming detector limitations of x-ray photon counting for preclinical microcomputed tomography

Spectral computed tomography (CT) using photon counting detectors (PCDs) can provide accurate tissue composition measurements by utilizing the energy dependence of x-ray attenuation in different materials. PCDs are especially suited for K-edge imaging, revealing the spatial distribution of select imaging probes through quantitative material decomposition. We report on a prototype spectral micro-CT system with a CZT-based PCD (DxRay, Inc.) that has 16 × 16    pixels of 0.5 × 0.5    mm 2 , a thickness of 3 mm, and four energy thresholds. Due to the PCD's limited size ( 8 × 8    mm 2 ), our system uses a translate-rotate projection acquisition strategy to cover a field of view relevant for preclinical imaging ( ∼ 4.5    cm ). Projection corrections were implemented to minimize artifacts associated with dead pixels and projection stitching. A sophisticated iterative algorithm was used to reconstruct both phantom and ex vivo mouse data. To achieve preclinically relevant spatial resolution, we trained a convolutional neural network to perform pan-sharpening between low-resolution PCD data ( 247 - μ m voxels) and high-resolution energy-integrating detector data ( 82 - μ m voxels), recovering a high-resolution estimate of the spectral contrast suitable for material decomposition. Long-term, preclinical spectral CT systems such as ours could serve in the developing field of theranostics (therapy and diagnostics) for cancer research.

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.

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.

Micro-CT of rodents: state-of-the-art and future perspectives

Micron-scale computed tomography (micro-CT) is an essential tool for phenotyping and for elucidating diseases and their therapies. This work is focused on preclinical micro-CT imaging, reviewing relevant principles, technologies, and applications. Commonly, micro-CT provides high-resolution anatomic information, either on its own or in conjunction with lower-resolution functional imaging modalities such as positron emission tomography (PET) and single-photon emission computed tomography (SPECT). More recently, however, advanced applications of micro-CT produce functional information by translating clinical applications to model systems (e.g., measuring cardiac functional metrics) and by pioneering new ones (e.g. measuring tumor vascular permeability with nanoparticle contrast agents). The primary limitations of micro-CT imaging are the associated radiation dose and relatively poor soft tissue contrast. We review several image reconstruction strategies based on iterative, statistical, and gradient sparsity regularization, demonstrating that high image quality is achievable with low radiation dose given ever more powerful computational resources. We also review two contrast mechanisms under intense development. The first is spectral contrast for quantitative material discrimination in combination with passive or actively targeted nanoparticle contrast agents. The second is phase contrast which measures refraction in biological tissues for improved contrast and potentially reduced radiation dose relative to standard absorption imaging. These technological advancements promise to develop micro-CT into a commonplace, functional and even molecular imaging modality.

Vision 20/20: Single photon counting x-ray detectors in medical imaging

Photon counting detectors (PCDs) with energy discrimination capabilities have been developed for medical x-ray computed tomography (CT) and x-ray (XR) imaging. Using detection mechanisms that are completely different from the current energy integrating detectors and measuring the material information of the object to be imaged, these PCDs have the potential not only to improve the current CT and XR images, such as dose reduction, but also to open revolutionary novel applications such as molecular CT and XR imaging. The performance of PCDs is not flawless, however, and it seems extremely challenging to develop PCDs with close to ideal characteristics. In this paper, the authors offer our vision for the future of PCD-CT and PCD-XR with the review of the current status and the prediction of (1) detector technologies, (2) imaging technologies, (3) system technologies, and (4) potential clinical benefits with PCDs.

Radial differential interior tomography and its image reconstruction with differentiated backprojection and projection onto convex sets

PURPOSE

Interior tomography has been recognized as one of the most effective approaches in computed tomography (CT) to reduce radiation dose rendered to patients. In this work, the authors propose and evaluate an imaging method of radial differential interior tomography.

METHODS

In interior tomography, an x-ray beam is collimated to only irradiate the region of interest (ROI) with suspected lesions while the surrounding area∕volume of normal tissues∕organs is spared. In the proposed imaging method of radial differential interior tomography, the outcome is a ROI image that has gone through a radial differential filtering. The image reconstruction algorithm for the radial differential interior tomography is kept in the fashion of differentiated backprojection and projection onto convex sets, but the required a priori knowledge in a small round area becomes zero and may be more readily available in practice.

RESULTS

Using the projection data simulated by computer and acquired by CT scanner, the authors evaluate and verify the performance of the proposed radial differential interior tomography method and its associated image reconstruction algorithm. The preliminary results show that the proposed imaging method can generate an image that is the radial differentiation of a conventional tomographic image and is robust over noise that inevitably exist in practice.

CONCLUSIONS

It is believed that the proposed imaging method may find its utility in advanced clinical applications wherein a ROI-based image processing and analysis is required for lesion visualization, characterization, and diagnosis.

The meaning of interior tomography

The classic imaging geometry for computed tomography is for the collection of un-truncated projections and the reconstruction of a global image, with the Fourier transform as the theoretical foundation that is intrinsically non-local. Recently, interior tomography research has led to theoretically exact relationships between localities in the projection and image spaces and practically promising reconstruction algorithms. Initially, interior tomography was developed for x-ray computed tomography. Then, it was elevated to have the status of a general imaging principle. Finally, a novel framework known as 'omni-tomography' is being developed for a grand fusion of multiple imaging modalities, allowing tomographic synchrony of diversified features.

Scout-view assisted interior micro-CT

Micro computed tomography (micro-CT) is a widely-used imaging technique. A challenge of micro-CT is to quantitatively reconstruct a sample larger than the field-of-view (FOV) of the detector. This scenario is characterized by truncated projections and associated image artifacts. However, for such truncated scans, a low resolution scout scan with an increased FOV is frequently acquired so as to position the sample properly. This study shows that the otherwise discarded scout scans can provide sufficient additional information to uniquely and stably reconstruct the interior region of interest. Two interior reconstruction methods are designed to utilize the multi-resolution data without significant computational overhead. While most previous studies used numerically truncated global projections as interior data, this study uses truly hybrid scans where global and interior scans were carried out at different resolutions. Additionally, owing to the lack of standard interior micro-CT phantoms, we designed and fabricated novel interior micro-CT phantoms for this study to provide means of validation for our algorithms. Finally, two characteristic samples from separate studies were scanned to show the effect of our reconstructions. The presented methods show significant improvements over existing reconstruction algorithms.

Volume-of-change cone-beam CT for image-guided surgery

C-arm cone-beam CT (CBCT) can provide intraoperative 3D imaging capability for surgical guidance, but workflow and radiation dose are the significant barriers to broad utilization. One main reason is that each 3D image acquisition requires a complete scan with a full radiation dose to present a completely new 3D image every time. In this paper, we propose to utilize patient-specific CT or CBCT as prior knowledge to accurately reconstruct the aspects of the region that have changed by the surgical procedure from only a sparse set of x-rays. The proposed methods consist of a 3D-2D registration between the prior volume and a sparse set of intraoperative x-rays, creating digitally reconstructed radiographs (DRRs) from the registered prior volume, computing difference images by subtracting DRRs from the intraoperative x-rays, a penalized likelihood reconstruction of the volume of change (VOC) from the difference images, and finally a fusion of VOC reconstruction with the prior volume to visualize the entire surgical field. When the surgical changes are local and relatively small, the VOC reconstruction involves only a small volume size and a small number of projections, allowing less computation and lower radiation dose than is needed to reconstruct the entire surgical field. We applied this approach to sacroplasty phantom data obtained from a CBCT test bench and vertebroplasty data with a fresh cadaver acquired from a C-arm CBCT system with a flat-panel detector. The VOCs were reconstructed from a varying number of images (10-66 images) and compared to the CBCT ground truth using four different metrics (mean squared error, correlation coefficient, structural similarity index and perceptual difference model). The results show promising reconstruction quality with structural similarity to the ground truth close to 1 even when only 15-20 images were used, allowing dose reduction by the factor of 10-20.

Time-resolved cardiac interventional cone-beam CT reconstruction from fully truncated projections using the prior image constrained compressed sensing (PICCS) algorithm

C-arm cone-beam CT could replace preoperative multi-detector CT scans in the cardiac interventional setting. However, cardiac gating results in view angle undersampling and the small size of the detector results in projection data truncation. These problems are incompatible with conventional tomographic reconstruction algorithms. In this paper, the prior image constrained compressed sensing (PICCS) reconstruction method was adapted to solve these issues. The performance of the proposed method was compared to that of FDK, FDK with extrapolated projection data (E-FDK), and total variation-based compressed sensing. A canine projection dataset acquired using a clinical C-arm imaging system supplied realistic cardiac motion and anatomy for this evaluation. Three different levels of truncation were simulated. The relative root mean squared error and the universal image quality index were used to quantify the reconstruction accuracy. Three main conclusions were reached. (1) The adapted version of the PICCS algorithm offered the highest image quality and reconstruction accuracy. (2) No meaningful variation in performance was observed when the amount of truncation was changed. (3) This study showed evidence that accurate interior tomography with an undersampled acquisition is possible for realistic objects if a prior image with minimal artifacts is available.

Practical interior tomography with radial Hilbert filtering and a priori knowledge in a small round area

PURPOSES

Interior tomography problem can be solved using the so-called differentiated backprojection-projection onto convex sets (DBP-POCS) method, which requires a priori knowledge within a small area interior to the region of interest (ROI) to be imaged. In theory, the small area wherein the a priori knowledge is required can be in any shape, but most of the existing implementations carry out the Hilbert filtering either horizontally or vertically, leading to a vertical or horizontal strip that may be across a large area in the object. In this work, we implement a practical DBP-POCS method with radial Hilbert filtering and thus the small area with the a priori knowledge can be roughly round (e.g., a sinus or ventricles among other anatomic cavities in human or animal body). We also conduct an experimental evaluation to verify the performance of this practical implementation.

METHODS

We specifically re-derive the reconstruction formula in the DBP-POCS fashion with radial Hilbert filtering to assure that only a small round area with the a priori knowledge be needed (namely radial DBP-POCS method henceforth). The performance of the practical DBP-POCS method with radial Hilbert filtering and a priori knowledge in a small round area is evaluated with projection data of the standard and modified Shepp-Logan phantoms simulated by computer, followed by a verification using real projection data acquired by a computed tomography (CT) scanner.

RESULTS

The preliminary performance study shows that, if a priori knowledge in a small round area is available, the radial DBP-POCS method can solve the interior tomography problem in a more practical way at high accuracy.

CONCLUSIONS

In comparison to the implementations of DBP-POCS method demanding the a priori knowledge in horizontal or vertical strip, the radial DBP-POCS method requires the a priori knowledge within a small round area only. Such a relaxed requirement on the availability of a priori knowledge can be readily met in practice, because a variety of small round areas (e.g., air-filled sinuses or fluid-filled ventricles among other anatomic cavities) exist in human or animal body. Therefore, the radial DBP-POCS method with a priori knowledge in a small round area is more feasible in clinical and preclinical practice.

Region-of-interest material decomposition from truncated energy-resolved CT

PURPOSE

Energy-resolved CT using photon-counting detectors has the potential to provide improved material decomposition compared to dual-kVp approaches. However, available photon-counting detectors are susceptible to pulse-pileup artifacts, especially at the periphery of the field of view (FOV) where the object attenuation is low compared to the center of the FOV. Pulse pileup may be avoided by imaging a region-of-interest (ROI) where the dynamic range is expected to be limited. This work investigated performing material decomposition and reconstructing ROI basis images from truncated energy-resolved data.

METHODS

A method is proposed to reconstruct images of basis functions primarily contained within the ROI, such as targeted or localized K-edge contrast agents. Material decomposition is performed independently for each ray in the sinogram, followed by filtered backprojection from the truncated data encompassing the ROI. A second method is proposed that uses a prior conventional energy-integrating image to estimate energy-resolved data outside the ROI. The measured and estimated energy-resolved data are decomposed into basis projections and merged into basis sinograms of the full FOV. Basis images of the ROI are then reconstructed through filtered backprojection. This method is most easily applied to objects that do not contain K-edge contrast agents outside the ROI. Simulations of a voxelized thorax phantom with iodine in the blood pool and a detector with five energy bins were performed. Full FOV, truncated, and truncated data merged with data estimated from the prior energy-integrating image were decomposed into Compton, photoelectric, and iodine basis functions. An empirical weighting factor was determined to blend the merged sinogram at the boundary of the truncated data. The effects of noise and misalignment in the prior image were also quantified. Basis images of the central 15 cm × 15 cm ROI containing the heart were reconstructed via filtered backprojection. Basis image accuracy was quantified relative to gold-standard basis images reconstructed from full FOV energy-resolved data.

RESULTS

The error in the iodine basis image reconstructed from truncated energy-resolved data without prior information was less than 1% for the central 7 cm of the 7.5-cm-radius ROI and 3% at the edge of the ROI. When the truncated and estimated basis sinograms were blended, the error was below 1% throughout the ROI for photoelectric basis images and ranged from 1% at the center of the ROI to 4% at the edge for the Compton basis image.

CONCLUSIONS

The density of localized K-edge contrast agents can be estimated to within 1% error using filtered back projection without prior information. For noncontrast and localized-contrast scans, ROI images of general basis functions can be reconstructed to within a few percent error using a prior energy-integrating image. The ability to perform material decomposition for a limited ROI may facilitate energy-resolved CT with available photon-counting detectors.