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

Improving Image Reconstruction for Digital Breast Tomosynthesis

Digital breast tomosynthesis (DBT) has been developed to reduce the issue of overlapping tissue in conventional 2-D mammography for breast cancer screening and diagnosis. In the DBT procedure, the patient’s breast is compressed with a paddle and a sequence of x-ray projections is taken within a small angular range. Tomographic reconstruction algorithms are then applied to these projections, generating tomosynthesized image slices of the breast, such that radiologists can read the breast slice by slice. Studies have shown that DBT can reduce both false-negative diagnoses of breast cancer and false-positive recalls compared to mammography alone. This dissertation focuses on improving image quality for DBT reconstruction. Chapter I briefly introduces the concept of DBT and the inspiration of my study. Chapter II covers the background of my research including the concept of image reconstruction, the geometry of our experimental DBT system and figures of merit for image quality. Chapter III introduces our study of the segmented separable footprint (SG) projector. By taking into account the finite size of detector element, the SG projector improves the accuracy of forward projections in iterative image reconstruction. Due to the more efficient access to memory, the SG projector is also faster than the traditional ray-tracing (RT) projector. We applied the SG projector to regular and subpixel reconstructions and demonstrated its effectiveness. Chapter IV introduces a new DBT reconstruction method with detector blur and correlated noise modeling, called the SQS-DBCN algorithm. The SQS-DBCN algorithm is able to significantly enhance microcalcifications (MC) in DBT while preserving the appearance of the soft tissue and mass margin. Comparisons between the SQS-DBCN algorithm and several modified versions of the SQS-DBCN algorithm indicate the importance of modeling different components of the system physics at the same time. Chapter V investigates truncated projection artifact (TPA) removal algorithms. Among the three algorithms we proposed, the pre-reconstruction-based projection view (PV) extrapolation method provides the best performance. Possible improvements of the other two TPA removal algorithms have been discussed. Chapter VI of this dissertation examines the effect of source blur on DBT reconstruction. Our analytical calculation demonstrates that the point spread function (PSF) of source blur is highly shift-variant. We used CatSim to simulate digital phantoms. Analysis on the reconstructed images demonstrates that a typical finite-sized focal spot (~ 0.3 mm) will not affect the image quality if the x-ray tube is stationary during the data acquisition. For DBT systems with continuous-motion data acquisition, the motion of the x-ray tube is the main cause of the effective source blur and will cause loss in the contrast of objects. Therefore modeling the source blur for these DBT systems could potentially improve the reconstructed image quality. The final chapter of this dissertation discusses a few future studies that are inspired by my PhD research.PHDElectrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/144059/1/jiabei_1.pd

An Enhanced Visualization of DBT Imaging Using Blind Deconvolution and Total Variation Minimization Regularization

Digital Breast Tomosynthesis (DBT) presents out-of-plane artifacts caused by features of high intensity. Given observed data and knowledge about the point spread function (PSF), deconvolution techniques recover data from a blurred version. However, a correct PSF is difficult to achieve and these methods amplify noise. When no information is available about the PSF, blind deconvolution can be used. Additionally, Total Variation (TV) minimization algorithms have achieved great success due to its virtue of preserving edges while reducing image noise. This work presents a novel approach in DBT through the study of out-of-plane artifacts using blind deconvolution and noise regularization based on TV minimization. Gradient information was also included. The methodology was tested using real phantom data and one clinical data set. The results were investigated using conventional 2D slice-by-slice visualization and 3D volume rendering. For the 2D analysis, the artifact spread function (ASF) and Full Width at Half Maximum (FWHMMASF) of the ASF were considered. The 3D quantitative analysis was based on the FWHM of disks profiles at 90°, noise and signal to noise ratio (SNR) at 0° and 90°. A marked visual decrease of the artifact with reductions of FWHMASF (2D) and FWHM90° (volume rendering) of 23.8% and 23.6%, respectively, was observed. Although there was an expected increase in noise level, SNR values were preserved after deconvolution. Regardless of the methodology and visualization approach, the objective of reducing the out-of-plane artifact was accomplished. Both for the phantom and clinical case, the artifact reduction in the z was markedly visible

A comparative study of limitedâ angle coneâ beam reconstruction methods for breast tomosynthesis

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134781/1/mp7543.pd

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

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

A fast Total Variation-based iterative algorithm for digital breast tomosynthesis image reconstruction:

In this work, we propose a fast iterative algorithm for the reconstruction of digital breast tomosynthesis images. The algorithm solves a regularization problem, expressed as the minimization of the sum of a least-squares term and a weighted smoothed version of the Total Variation regularization function. We use a Fixed Point method for the solution of the minimization problem, requiring the solution of a linear system at each iteration, whose coefficient matrix is a positive definite approximation of the Hessian of the objective function. We propose an efficient implementation of the algorithm, where the linear system is solved by a truncated Conjugate Gradient method. We compare the Fixed Point implementation with a fast first order method such as the Scaled Gradient Projection method, that does not require any linear system solution. Numerical experiments on a breast phantom widely used in tomographic simulations show that both the methods recover microcalcifications very fast while the Fixed Point is more efficient in detecting masses, when more time is available for the algorithm execution

QUANTITATIVE THREE DIMENSIONAL ELASTICITY IMAGING

Neoplastic tissue is typically highly vascularized, contains abnormal concentrations of extracellular proteins (e.g. collagen, proteoglycans) and has a high interstitial fluid pres- sure compared to most normal tissues. These changes result in an overall stiffening typical of most solid tumors. Elasticity Imaging (EI) is a technique which uses imaging systems to measure relative tissue deformation and thus noninvasively infer its mechanical stiffness. Stiffness is recovered from measured deformation by using an appropriate mathematical model and solving an inverse problem. The integration of EI with existing imaging modal- ities can improve their diagnostic and research capabilities. The aim of this work is to develop and evaluate techniques to image and quantify the mechanical properties of soft tissues in three dimensions (3D). To that end, this thesis presents and validates a method by which three dimensional ultrasound images can be used to image and quantify the shear modulus distribution of tissue mimicking phantoms. This work is presented to motivate and justify the use of this elasticity imaging technique in a clinical breast cancer screening study. The imaging methodologies discussed are intended to improve the specificity of mammography practices in general. During the development of these techniques, several issues concerning the accuracy and uniqueness of the result were elucidated. Two new algorithms for 3D EI are designed and characterized in this thesis. The first provides three dimensional motion estimates from ultrasound images of the deforming ma- terial. The novel features include finite element interpolation of the displacement field, inclusion of prior information and the ability to enforce physical constraints. The roles of regularization, mesh resolution and an incompressibility constraint on the accuracy of the measured deformation is quantified. The estimated signal to noise ratio of the measured displacement fields are approximately 1800, 21 and 41 for the axial, lateral and eleva- tional components, respectively. The second algorithm recovers the shear elastic modulus distribution of the deforming material by efficiently solving the three dimensional inverse problem as an optimization problem. This method utilizes finite element interpolations, the adjoint method to evaluate the gradient and a quasi-Newton BFGS method for optimiza- tion. Its novel features include the use of the adjoint method and TVD regularization with piece-wise constant interpolation. A source of non-uniqueness in this inverse problem is identified theoretically, demonstrated computationally, explained physically and overcome practically. Both algorithms were test on ultrasound data of independently characterized tissue mimicking phantoms. The recovered elastic modulus was in all cases within 35% of the reference elastic contrast. Finally, the preliminary application of these techniques to tomosynthesis images showed the feasiblity of imaging an elastic inclusion.CenSSIS, the Center for Subsurface Sensing and Imaging Systems, under the Engineering Research Centers Program of the National Science Foundation (award number EEC-9986821

Image reconstruction and imaging configuration optimization with a novel nanotechnology enabled breast tomosynthesis multi-beam X-ray system

Digital breast tomosynthesis is a new technology that provides three-dimensional information of the breast and makes it possible to distinguish the cancer from overlying breast tissues. We are dedicated to optimizing image reconstruction and imaging configuration for a new multi-beam parallel digital breast tomosynthesis prototype system. Several commonly used algorithms from the typical image reconstruction models which were used for iso-centric tomosynthesis systems were investigated for our multi-beam parallel tomosynthesis imaging system. The representative algorithms, including back-projection (BP), filtered back-projection (FBP), matrix inversion tomosynthesis reconstruction (MITS), maximum likelihood expectation maximization (MLEM), ordered-subset maximum likelihood expectation maximization (OS-MLEM), simultaneous algebraic reconstruction technique (SART), were implemented to fit our system design. An accelerated MLEM algorithm was proposed, which significantly reduced the running time but had the same image quality. Furthermore, two statistical variants of BP reconstruction were validated for our tomosynthesis prototype system. Experiments based on phantoms and computer simulations show that the prototype system combined with our algorithms is capable of providing three-dimensional information of the objects with good image quality and has great potentials to improve digital breast tomosynthesis technology. Four methodologies were employed to optimize the reconstruction algorithms and different imaging configurations for the prototype system. A linear tomosynthesis imaging analysis tool was used to investigate blurring-out reconstruction algorithms. Computer simulations of sphere and wire objects aimed at the performance of out-of-plane artifact removal. A frequency-domain-based methodology, relative NEQ(f) analysis, was investigated to evaluate the overall system performance based on the propagation of signal and noise. Conclusions were made to determine the optimal image reconstruction algorithm and imaging configuration of this new multi-beam parallel digital breast tomosynthesis prototype system for better image quality and system performance