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

Unconstrained simultaneous scheme to fully couple reconstruction and registration for digital breast tomosynthesis: a feasible study

Digital breast tomosynthesis (DBT) provides a pseudo-3D reconstruction which addresses the limitation of superimposition of dense fibro-glandular tissue associated with conventional mammography. Reg- istration of temporal DBT volumes searches for the optimum deforma- tion to transform two observed images of the same object into a common reference frame. This aligns the two images via minimising an objective function that calculates the similarity between the two datasets. In this paper, we present a novel algorithm which combines recon- struction of a pair of temporal DBT acquisitions with their simultaneous registration. We approach this nonlinear inverse problem using a generic unconstrained optimisation scheme. To evaluate the performance of our method we use 2D and 3D software phantoms and demonstrate that this simultaneous approach has comparable results to performing these tasks sequentially or iteratively w.r.t both the reconstruction fidelity and the registration accuracy

A Digital X-Ray Tomosynthesis Coupled Near Infrared Spectral Tomography System for Dual-Modality Breast Imaging

A Near Infrared Spectral Tomography (NIRST) system has been developed and integrated into a commercial Digital Breast Tomosynthesis (DBT) scanner to allow structural and functional imaging of breast in vivo. The NIRST instrument uses an 8-wavelength continuous wave (CW) laser-based scanning source assembly and a 75-element silicon photodiode solid-state detector panel to produce dense spectral and spatial projection data from which spectrally constrained 3D tomographic images of tissue chromophores are produced. Integration of the optical imaging system into the DBT scanner allows direct co-registration of the optical and DBT images, while also facilitating the synergistic use of x-ray contrast as anatomical priors in optical image reconstruction. Currently, the total scan time for a combined NIRST-DBT exam is ~50s with data collection from 8 wavelengths in the optical scan requiring ~42s to complete. The system was tested in breast simulating phantoms constructed using intralipid and blood in an agarose matrix with a 3 cm x 2 cm cylindrical inclusion at 1 cm depth from the surface. Diffuse image reconstruction of total hemoglobin (HbT) concentration resulted in accurate recovery of the lateral size and position of the inclusion to within 6% and 8%, respectively. Use of DBT structural priors in the NIRST reconstruction process improved the quantitative accuracy of the HbT recovery, and led to linear changes in imaged versus actual contrast, underscoring the advantages of dual-modality optical imaging approaches. The quantitative accuracy of the system can be further improved with independent measurements of scattering properties through integration of frequency or time domain data

Breast Mass Characterization Using 3‐Dimensional Automated Ultrasound as an Adjunct to Digital Breast Tomosynthesis

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/135628/1/jum201332193.pd

Quantifying the Tibiofemoral Joint Space Using X-ray Tomosynthesis

Purpose: Digital x-ray tomosynthesis (DTS) has the potential to provide 3D information about the knee joint in a load-bearing posture, which may improve diagnosis and monitoring of knee osteoarthritis compared with projection radiography, the current standard of care. Manually quantifying and visualizing the joint space width (JSW) from 3D tomosynthesis datasets may be challenging. This work developed a semiautomated algorithm for quantifying the 3D tibiofemoral JSW from reconstructed DTS images. The algorithm was validated through anthropomorphic phantom experiments and applied to three clinical datasets. Methods: A user-selected volume of interest within the reconstructed DTS volume was enhanced with 1D multiscale gradient kernels. The edge-enhanced volumes were divided by polarity into tibial and femoral edge maps and combined across kernel scales. A 2D connected components algorithm was performed to determine candidate tibial and femoral edges. A 2D joint space width map (JSW) was constructed to represent the 3D tibiofemoral joint space. To quantify the algorithm accuracy, an adjustable knee phantom was constructed, and eleven posterior–anterior (PA) and lateral DTS scans were acquired with the medial minimum JSW of the phantom set to 0–5 mm in 0.5 mm increments (VolumeRadTM, GE Healthcare, Chalfont St. Giles, United Kingdom). The accuracy of the algorithm was quantified by comparing the minimum JSW in a region of interest in the medial compartment of the JSW map to the measured phantom setting for each trial. In addition, the algorithm was applied to DTS scans of a static knee phantom and the JSW map compared to values estimated from a manually segmented computed tomography (CT) dataset. The algorithm was also applied to three clinical DTS datasets of osteoarthritic patients. Results: The algorithm segmented the JSW and generated a JSW map for all phantom and clinical datasets. For the adjustable phantom, the estimated minimum JSW values were plotted against the measured values for all trials. A linear fit estimated a slope of 0.887 (R2¼0.962) and a mean error across all trials of 0.34 mm for the PA phantom data. The estimated minimum JSW values for the lateral adjustable phantom acquisitions were found to have low correlation to the measured values (R2¼0.377), with a mean error of 2.13 mm. The error in the lateral adjustable-phantom datasets appeared to be caused by artifacts due to unrealistic features in the phantom bones. JSW maps generated by DTS and CT varied by a mean of 0.6 mm and 0.8 mm across the knee joint, for PA and lateral scans. The tibial and femoral edges were successfully segmented and JSW maps determined for PA and lateral clinical DTS datasets. Conclusions: A semiautomated method is presented for quantifying the 3D joint space in a 2D JSW map using tomosynthesis images. The proposed algorithm quantified the JSW across the knee joint to sub-millimeter accuracy for PA tomosynthesis acquisitions. Overall, the results suggest that x-ray tomosynthesis may be beneficial for diagnosing and monitoring disease progression or treatment of osteoarthritis by providing quantitative images of JSW in the load-bearing knee

Calibration and Optimization of 3D Digital Breast Tomosynthesis Guided Near Infrared Spectral Tomography

Calibration of a three-dimensional multimodal digital breast tomosynthesis (DBT) x-ray and non-fiber based near infrared spectral tomography (NIRST) system is challenging but essential for clinical studies. Phantom imaging results yielded linear contrast recovery of total hemoglobin (HbT) concentration for cylindrical inclusions of 15 mm, 10 mm and 7 mm with a 3.5% decrease in the HbT estimate for each 1 cm increase in inclusion depth. A clinical exam of a patient\u27s breast containing both benign and malignant lesions was successfully imaged, with greater HbT was found in the malignancy relative to the benign abnormality and fibroglandular regions (11 μM vs. 9.5 μM). Tools developed improved imaging system characterization and optimization of signal quality, which will ultimately improve patient selection and subsequent clinical trial results

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

X-ray absorption tomography employing a conical shell beam

We demonstrate depth-resolved absorption imaging by scanning an object through a conical shell of X-rays. We measure ring shaped projections and apply tomosynthesis to extract optical sections at different axial focal plane positions. Three-dimensional objects have been imaged to validate our theoretical treatment. The novel principle of our method is scalable with respect to both scan size and X-ray energy. A driver for this work is to complement previously reported methods concerning the measurement of diffracted X-rays for structural analysis. The prospect of employing conical shell beams to combine both absorption and diffraction modalities would provide enhanced analytical utility and has many potential applications in security screening, process control and diagnostic imaging