Combined Reconstruction and Registration of Digital Breast Tomosynthesis

Digital breast tomosynthesis (DBT) has the potential to en- hance breast cancer detection by reducing the confounding e ect of su- perimposed tissue associated with conventional mammography. In addi- tion the increased volumetric information should enable temporal datasets to be more accurately compared, a task that radiologists routinely apply to conventional mammograms to detect the changes associated with ma- lignancy. In this paper we address the problem of comparing DBT data by combining reconstruction of a pair of temporal volumes with their reg- istration. Using a simple test object, and DBT simulations from in vivo breast compressions imaged using MRI, we demonstrate that this com- bined reconstruction and registration approach produces improvements in both the reconstructed volumes and the estimated transformation pa- rameters when compared to performing the tasks sequentially

Numerical methods for coupled reconstruction and registration in digital breast tomosynthesis.

Digital Breast Tomosynthesis (DBT) provides an insight into the fine details of normal fibroglandular tissues and abnormal lesions by reconstructing a pseudo-3D image of the breast. In this respect, DBT overcomes a major limitation of conventional X-ray mam- mography by reducing the confounding effects caused by the superposition of breast tissue. In a breast cancer screening or diagnostic context, a radiologist is interested in detecting change, which might be indicative of malignant disease. To help automate this task image registration is required to establish spatial correspondence between time points. Typically, images, such as MRI or CT, are first reconstructed and then registered. This approach can be effective if reconstructing using a complete set of data. However, for ill-posed, limited-angle problems such as DBT, estimating the deformation is com- plicated by the significant artefacts associated with the reconstruction, leading to severe inaccuracies in the registration. This paper presents a mathematical framework, which couples the two tasks and jointly estimates both image intensities and the parameters of a transformation. Under this framework, we compare an iterative method and a simultaneous method, both of which tackle the problem of comparing DBT data by combining reconstruction of a pair of temporal volumes with their registration. We evaluate our methods using various computational digital phantoms, uncom- pressed breast MR images, and in-vivo DBT simulations. Firstly, we compare both iter- ative and simultaneous methods to the conventional, sequential method using an affine transformation model. We show that jointly estimating image intensities and parametric transformations gives superior results with respect to reconstruction fidelity and regis- tration accuracy. Also, we incorporate a non-rigid B-spline transformation model into our simultaneous method. The results demonstrate a visually plausible recovery of the deformation with preservation of the reconstruction fidelity

An attempt to estimate out-of-plane lung nodule elongation in tomosynthesis images

In chest tomosynthesis (TS) the most commonly used reconstruction methods are based on Filtered Back Projection (FBP) algorithms. Due to the limited angular range of x-ray projections, FBP reconstructed data is typically associated with a low spatial resolution in the out-of-plane dimension. Lung nodule measures that depend on depth information such as 3D shape and volume are therefore difficult to estimate. In this paper the relation between features from FBP reconstructed lung nodules and the true out-of-plane nodule elongation is investigated and a method for estimating the out-of-plane nodule elongation is proposed. In order to study these relations a number of steps that include simulation of spheroidal-shaped nodules, insertion into synthetic data volumes, construction of TS-projections and FBP-reconstruction were performed. In addition, the same procedure was used to simulate nodules and insert them into clinical chest TS projection data. The reconstructed nodule data was then investigated with respect to in-plane diameter, out-of-plane elongation, and attenuation coefficient. It was found that the voxel value in each nodule increased linearly with nodule elongation, for nodules with a constant attenuation coefficient. Similarly, the voxel value increased linearly with in-plane diameter. These observations indicate the possibility to predict the nodule elongation from the reconstructed voxel intensity values. Such a method would represent a quantitative approach to chest tomosynthesis that may be useful in future work on volume and growth rate estimation of lung nodules

Modeling the Anisotropic Resolution and Noise Properties of Digital Breast Tomosynthesis

Digital breast tomosynthesis (DBT) is a 3D imaging modality in which a reconstruction of the breast is generated from various x-ray projections. Due to the newness of this technology, the development of an analytical model of image quality has been on-going. In this thesis, a more complete model is developed by addressing the limitations found in the previous linear systems (LS) model [Zhao, Med. Phys. 2008, 35(12): 5219-32]. A central assumption of the LS model is that the angle of x-ray incidence is approximately normal to the detector in each projection. To model the effect of oblique x-ray incidence, this thesis generalizes Swank\u27s calculations of the transfer functions of x-ray fluorescent screens to arbitrary incident angles. In the LS model, it is also assumed that the pixelation in the reconstruction grid is the same as the detector; hence, the highest frequency that can be resolved is the detector alias frequency. This thesis considers reconstruction grids with smaller pixelation to investigate super-resolution, or visibility of higher frequencies. A sine plate is introduced as a conceptual test object to analyze super-resolution. By orienting the long axis of the sine plate at various angles, the feasibility of oblique reconstruction planes is also investigated. This formulation differs from the LS model in which reconstruction planes are parallel to the breast support. It is shown that the transfer functions for arbitrary angles of x-ray incidence can be modeled in closed form. The high frequency modulation transfer function (MTF) and detective quantum efficiency (DQE) are degraded due to oblique x-ray incidence. In addition, using the sine plate, it is demonstrated that a reconstruction can resolve frequencies exceeding the detector alias frequency. Experimental images of bar patterns verified the existence of super-resolution. Anecdotal clinical examples showed that super-resolution improves the visibility of microcalcifications. The feasibility of oblique reconstructions was established theoretically with the sine plate and was validated experimentally with bar patterns. This thesis develops a more complete model of image quality in DBT by addressing the limitations of the LS model. In future studies, this model can be used as a tool for optimizing DBT

Comparison of different image reconstruction algorithms for Digital Breast Tomosynthesis and assessment of their potential to reduce radiation dose

Tese de mestrado, Engenharia Física, 2022, Universidade de Lisboa, Faculdade de CiênciasDigital Breast Tomosynthesis is a three-dimensional medical imaging technique that allows the view of sectional parts of the breast. Obtaining multiple slices of the breast constitutes an advantage in contrast to conventional mammography examination in view of the increased potential in breast cancer detectability. Conventional mammography, despite being a screening success, has undesirable specificity, sensitivity, and high recall rates owing to the overlapping of tissues. Although this new technique promises better diagnostic results, the acquisition methods and image reconstruction algorithms are still under research. Several articles suggest the use of analytic algorithms. However, more recent articles highlight the iterative algorithm’s potential for increasing image quality when compared to the former. The scope of this dissertation was to test the hypothesis of achieving higher quality images using iterative algorithms acquired with lower doses than those using analytic algorithms. In a first stage, the open-source Tomographic Iterative GPU-based Reconstruction (TIGRE) Toolbox for fast and accurate 3D x-ray image reconstruction was used to reconstruct the images acquired using an acrylic phantom. The algorithms used from the toolbox were the Feldkamp, Davis, and Kress, the Simultaneous Algebraic Reconstruction Technique, and the Maximum Likelihood Expectation Maximization algorithm. In a second and final state, the possibility of further reducing the radiation dose using image postprocessing tools was evaluated. A Total Variation Minimization filter was applied to the images reconstructed with the TIGRE toolbox algorithm that provided the best image quality. These were then compared to the images of the commercial unit used for the image acquisitions. With the use of image quality parameters, it was found that the Maximum Likelihood Expectation Maximization algorithm performance was the best of the three for lower radiation doses, especially with the filter. In sum, the result showed the potential of the algorithm in obtaining images with quality for low doses

Recent advances in x-ray cone-beam computed laminography

X-ray computed tomography is a well established volume imaging technique used routinely in medical diagnosis, industrial non-destructive testing, and a wide range of scientific fields. Traditionally, computed tomography uses scanning geometries with a single axis of rotation together with reconstruction algorithms specifically designed for this setup. Recently there has however been increasing interest in more complex scanning geometries. These include so called X-ray computed laminography systems capable of imaging specimens with large lateral dimensions, or large aspect ratios, neither of which are well suited to conventional CT scanning procedures. Developments throughout this field have thus been rapid, including the introduction of novel system trajectories, the application and refinement of various reconstruction methods, and the use of recently developed computational hardware and software techniques to accelerate reconstruction times. Here we examine the advances made in the last several years and consider their impact on the state of the art

Iterative CT reconstruction from few projections for the nondestructive post irradiation examination of nuclear fuel assemblies

The core components (e.g. fuel assemblies, spacer grids, control rods) of the nuclear reactors encounter harsh environment due to high temperature, physical stress, and a tremendous level of radiation. The integrity of these elements is crucial for safe operation of the nuclear power plants. The Post Irradiation Examination (PIE) can reveal information about the integrity of the elements during normal operations and off‐normal events. Computed tomography (CT) is a tool for evaluating the structural integrity of elements non-destructively. CT requires many projections to be acquired from different view angles after which a mathematical algorithm is adopted for reconstruction. Obtaining many projections is laborious and expensive in nuclear industries. Reconstructions from a small number of projections are explored to achieve faster and cost-efficient PIE. Classical reconstruction algorithms (e.g. filtered back projection) cannot offer stable reconstructions from few projections and create severe streaking artifacts. In this thesis, conventional algorithms are reviewed, and new algorithms are developed for reconstructions of the nuclear fuel assemblies using few projections. CT reconstruction from few projections falls into two categories: the sparse-view CT and the limited-angle CT or tomosynthesis. Iterative reconstruction algorithms are developed for both cases in the field of compressed sensing (CS). The performance of the algorithms is assessed using simulated projections and validated through real projections. The thesis also describes the systematic strategy towards establishing the conditions of reconstructions and finds the optimal imaging parameters for reconstructions of the fuel assemblies from few projections. --Abstract, page iii

System Characterizations and Optimized Reconstruction Methods for Novel X-ray Imaging

In the past decade there have been many new emerging X-ray based imaging technologies developed for different diagnostic purposes or imaging tasks. However, there exist one or more specific problems that prevent them from being effectively or efficiently employed. In this dissertation, four different novel X-ray based imaging technologies are discussed, including propagation-based phase-contrast (PB-XPC) tomosynthesis, differential X-ray phase-contrast tomography (D-XPCT), projection-based dual-energy computed radiography (DECR), and tetrahedron beam computed tomography (TBCT). System characteristics are analyzed or optimized reconstruction methods are proposed for these imaging modalities. In the first part, we investigated the unique properties of propagation-based phase-contrast imaging technique when combined with the X-ray tomosynthesis. Fourier slice theorem implies that the high frequency components collected in the tomosynthesis data can be more reliably reconstructed. It is observed that the fringes or boundary enhancement introduced by the phase-contrast effects can serve as an accurate indicator of the true depth position in the tomosynthesis in-plane image. In the second part, we derived a sub-space framework to reconstruct images from few-view D-XPCT data set. By introducing a proper mask, the high frequency contents of the image can be theoretically preserved in a certain region of interest. A two-step reconstruction strategy is developed to mitigate the risk of subtle structures being oversmoothed when the commonly used total-variation regularization is employed in the conventional iterative framework. In the thirt part, we proposed a practical method to improve the quantitative accuracy of the projection-based dual-energy material decomposition. It is demonstrated that applying a total-projection-length constraint along with the dual-energy measurements can achieve a stabilized numerical solution of the decomposition problem, thus overcoming the disadvantages of the conventional approach that was extremely sensitive to noise corruption. In the final part, we described the modified filtered backprojection and iterative image reconstruction algorithms specifically developed for TBCT. Special parallelization strategies are designed to facilitate the use of GPU computing, showing demonstrated capability of producing high quality reconstructed volumetric images with a super fast computational speed. For all the investigations mentioned above, both simulation and experimental studies have been conducted to demonstrate the feasibility and effectiveness of the proposed methodologies

Application of boundary detection information in breast tomosynthesis reconstruction

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/135006/1/mp1968.pd