Performance analysis of reconstruction algorithms for breast microwave imaging

The problem of localizing a scatterer, which represents a tumor, in a homogeneous circular domain, which represents a breast, is addressed. A breast imaging method based on microwaves is considered. The microwave imaging involves to several techniques for detecting, localizing and characterizing tumors in breast tissues. In all such methods an electromagnetic inverse scattering problem exists. For the scattering detection method, an algorithm based on a linear procedure solution, inspired by MUltiple SIgnal Classification algorithm (MUSIC) and Time Reversal method (TR), is implemented. The algorithm returns a reconstructed image of the investigation domain in which it is detected the scatterer position. This image is called pseudospectrum. A preliminary performance analysis of the algorithm vying the working frequency is performed: the resolution and the signal-to-noise ratio of the pseudospectra are improved if a multi-frequency approach is considered. The Geometrical Mean-MUSIC algorithm (GM- MUSIC) is proposed as multi-frequency method. The performance of the GMMUSIC is tested in different real life computer simulations. The performed analysis shows that the algorithm detects the scatterer until the electrical parameters of the breast are known. This is an evident limit, since, in a real life situation, the anatomy of the breast is unknown. An improvement in GM-MUSIC is proposed: the Eye-GMMUSIC algorithm. Eye-GMMUSIC algorithm needs no a priori information on the electrical parameters of the breast. It is an optimizing algorithm based on the pattern search algorithm: it searches the breast parameters which minimize the Signal-to-Clutter Mean Ratio (SCMR) in the signal. Finally, the GM-MUSIC and the Eye-GMMUSIC algorithms are tested on a microwave breast cancer detection system consisting of an dipole antenna, a Vector Network Analyzer and a novel breast phantom built at University of Bologna. The reconstruction of the experimental data confirm the GM-MUSIC ability to localize a scatterer in a homogeneous medium.

Microwave breast imaging reconstruction

In this thesis I analyzed the microwave tomography method to recognize breast can- cer. I study how identify the dielectric permittivity, the Helmoltz equation parameter used to model the real physic problem. Through a non linear least squares method I solve a problem of parameters identification; I show the theoric approach and the devel- opment to reach the results. I use the Levenberg-Marquardt algorithm, applied on COMSOL software to multiphysic models; so I do numerical proofs on semplified test problems compared to the specific real problem to solve.

Tomographic 3D reconstruction through linear inversion for breast microwave imaging

In the present thesis we address the problem of detecting and localizing a small spherical target with characteristic electrical properties inside a volume of cylindrical shape, representing female breast, with MWI. One of the main works of this project is to properly extend the existing linear inversion algorithm from planar slice to volume reconstruction; results obtained, under the same conditions and experimental setup are reported for the two different approaches. Preliminar comparison and performance analysis of the reconstruction algorithms is performed via numerical simulations in a software-created environment: a single dipole antenna is used for illuminating the virtual breast phantom from different positions and, for each position, the corresponding scattered field value is registered. Collected data are then exploited in order to reconstruct the investigation domain, along with the scatterer position, in the form of image called pseudospectrum. During this process the tumor is modeled as a dielectric sphere of small radius and, for electromagnetic scattering purposes, it's treated as a point-like source. To improve the performance of reconstruction technique, we repeat the acquisition for a number of frequencies in a given range: the different pseudospectra, reconstructed from single frequency data, are incoherently combined with MUltiple SIgnal Classification (MUSIC) method which returns an overall enhanced image. We exploit multi-frequency approach to test the performance of 3D linear inversion reconstruction algorithm while varying the source position inside the phantom and the height of antenna plane. Analysis results and reconstructed images are then reported. Finally, we perform 3D reconstruction from experimental data gathered with the acquisition system in the microwave laboratory at DIFA, University of Bologna for a recently developed breast-phantom prototype; obtained pseudospectrum and performance analysis for the real model are reported.

Compressed sensing in digital tomosynthesis reconstruction

In this work we study a model for the breast image reconstruction in Digital Tomosynthesis, that is a non-invasive and non-destructive method for the three-dimensional visualization of the inner structures of an object, in which the data acquisition includes measuring a limited number of low-dose two-dimensional projections of an object by moving a detector and an X-ray tube around the object within a limited angular range. The problem of reconstructing 3D images from the projections provided in the Digital Tomosynthesis is an ill-posed inverse problem, that leads to a minimization problem with an object function that contains a data fitting term and a regularization term. The contribution of this thesis is to use the techniques of the compressed sensing, in particular replacing the standard least squares problem of data fitting with the problem of minimizing the 1-norm of the residuals, and using as regularization term the Total Variation (TV). We tested two different algorithms: a new alternating minimization algorithm (ADM), and a version of the more standard scaled projected gradient algorithm (SGP) that involves the 1-norm. We perform some experiments and analyse the performance of the two methods comparing relative errors, iterations number, times and the qualities of the reconstructed images. In conclusion we noticed that the use of the 1-norm and the Total Variation are valid tools in the formulation of the minimization problem for the image reconstruction resulting from Digital Tomosynthesis and the new algorithm ADM has reached a relative error comparable to a version of the classic algorithm SGP and proved best in speed and in the early appearance of the structures representing the masses.

Experimental evaluation of the z-resolution in different clinical Digital Breast Tomosynthesis systems using commercial phantoms

Digital Breast Tomosynthesis (DBT) is an advanced mammography technique based on the reconstruction of a pseudo-volumetric image. To date, image quality represents the most deficient section of DBT quality control protocols. In fact, related tests are not yet characterized by either action levels or typical values. This thesis work focuses on the evaluation of one aspect of image quality: the z-resolution. The latter is studied in terms of Artifact Spread Function (ASF), a function that describes the signal spread of a detail along the reconstructed focal planes. To quantify the ASF numerically, its Full Width at Half Maximum (FWHM) is calculated and used as a representative index of z-resolution. Experimental measurements were acquired in 24 DBT systems, of 7 different models, currently in use in 20 hospital facilities in Italy. The analysis, performed on the clinical reconstructed images, of 5 different commercial phantoms, lead to the identification of characteristic FWHM values for each type of DBT system. The ASF clearly showed a dependence on the size of the detail, providing higher FWHM values for larger objects. The z-resolution was found to be positively influenced by the acquisition angle: Fujifilm sistematically showed wider ASF profiles in ST mode (15°) than in HR mode (40°). However, no clear relationship was found between angular range and ASF, among different DBT systems, due to the influence of the peculiarities of each reconstruction algorithm. The experimental approach shown in this thesis work can be proposed as a z-resolution quality control test procedure. Contextually, the values found could be used as a starting point for identifying typical values to be included in the test, in a DBT protocol. Clearly, a statistically significant number of images is needed to do this. The equipment involved in this work is located in hospitals and is not available for research purposes, so only a limited amount of data was acquired and processed.

Characterization of a prototype for cone beam breast computed tomography

The project aims to experiment the Cone Beam Breast Computed Tomography technique using a standard digital mammography system. The work is focused on the definition of a protocol of quality measurements for the pre-clinical evaluation of the machine. The paper is developed in two parts. The first is specifically concerned with the methods used to define the image quality and dosimetry aspects specific for digital mammography devices. A complete characterization of the system has been performed according to the applicable IEC standards to assure the performances of the equipment and define the quality levels. Due to the lack of a quality control protocol dedicated to CBBCT mammography scanner, a new equivalent test procedure has been proposed. The second part of the paper is focused on the evaluation, through quantitative and visual analyzes, of the CBCT exam feasibility in the hardware and software conditions currently proposed by IMS Giotto. The prototype was in fact developed differing from the technical choices of competing companies and developed for a different intended use. The main difference with respect to the existing breast CT scanners is the possibility of performing on the same system the CBBCT scanning but also all the mammographic techniques. In this thesis, we aim to assess whether, in the current setup, considering a dosimetric range very close to that used in the clinic, the tests produce results that can be considered acceptable or at least indicative of the feasibility of the entire project from a commercial point of view. For this purpose, the final reconstruction images, obtained by two previously developed software, are analyzed.