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.

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.

Experimental set up and characterization of a prototype for breast microwave imaging

Nella presente tesi è stato sviluppato un sistema di acquisizione automatico finalizzato allo studio del breast microwave imaging. Le misure sono state eseguite in configurazione monostatica, in cui viene acquisito un segnale da diverse posizioni lungo il perimetro dell’area di indagine. A questo scopo, è stato installato un motore ad alta precisione che permette la rotazione del fantoccio e l’esecuzione automatica delle misure da un numero di posizioni fissato. Per automatizzare il processo di acquisizione, è stato inoltre sviluppato appositamente un software in ambiente LabView. Successivamente, è stata eseguita una intensa sessione di misure finalizzate alla caratterizzazione del sistema sviluppato al variare delle condizioni di misura. Abbiamo quindi utilizzato dei fantocci di tumore di diverse dimensioni e permittività elettrica per studiare la sensibilità della strumentazione in condizione di mezzo omogeneo. Dall’analisi delle ricostruzioni multifrequenza effettuate tramite diversi algoritmi di tipo TR-MUSIC sul range di frequenze selezionato, abbiamo notato che il tumore è ricostruito correttamente in tutti gli scenari testati. Inoltre, abbiamo creato un ulteriore fantoccio per simulare la presenza di una disomogeneità nel dominio di imaging. In questo caso, abbiamo studiato le performances del sistema di acquisizione al variare della posizione del tumore, le cui caratteristiche sono state fissate, e della permittività associata al fantoccio. Dall’analisi dei risultati appare che le performances di ricostruzione sono condizionate dalla presenza della disomogeneità, in modo particolare se il tumore è posizionato all’interno di essa. Infine, abbiamo studiato delle performance di due algoritmi di ricostruzione 3D: uno di essi è basato sulla sovrappo- sizione tomografica e sfrutta metodi di interpolazione, l’altro si basa sull’utilizzo di un propagatore 3D per il dipolo Hertziano in approssimazione scalare.

Microwave imaging of the axilla to aid breast cancer diagnosis

Breast cancer is the most common type of cancer among women all over the world. An important issue that is not commonly addressed in breast cancer imaging literature is the importance of imaging the underarm region—where up to 80% of breast cancer cells can metastasise to. The first axillary lymph nodes to receive drainage from the primary tumour in the breast are called Sentinel Node. If cancer cells are found in the Sentinel Node, there is an increased risk of metastatic breast cancer which makes this evaluation crucial to decide what follow-up exams and therapy to follow. However, non-invasive detection of cancer cells in the lymph nodes is often inconclusive, leading to the surgical removal of too many nodes which causes adverse side-effects for patients. Microwave Imaging is one of the most promising non-invasive imaging modalities for breast cancer early screening and monitoring. This novel study tests the feasibility of imaging the axilla region by means of the simulation of an Ultra-Wideband Microwave Imaging system. Simulations of such system are completed in several 2D underarm models that mimic the axilla. Initial imaging results are obtained by means of processing the simulated backscattered signals by eliminating artefacts caused by the skin and beamforming the processed signals in order to time-align all the signals recorded at each antenna. In this dissertation several image formation algorithms are implemented and compared by visual inspection of the resulting images and through a range of performance metrics, such as Signal-to-Clutter Ratio and FullWidth Half Maximum calculations. The results in this study showed that Microwave Imaging is a promising technique that might allow to identify the presence and location of metastasised cancer cells in axillary lymph nodes, enabling the non-invasive evaluation of breast cancer staging.

Diffusional kurtosis imaging using a fast heuristic constrained linear least squares algorithm: a plugin for OsiriX

Diffusion Kurtosis Imaging (DKI) is a fairly new magnetic resonance imag-ing (MRI) technique that tackles the non-gaussian motion of water in biological tissues by taking into account the restrictions imposed by tissue microstructure, which are not considered in Diffusion Tensor Imaging (DTI), where the water diffusion is considered purely gaussian. As a result DKI provides more accurate information on biological structures and is able to detect important abnormalities which are not visible in standard DTI analysis. This work regards the development of a tool for DKI computation to be implemented as an OsiriX plugin. Thus, as OsiriX runs under Mac OS X, the pro-gram is written in Objective-C and also makes use of Apple’s Cocoa framework. The whole program is developed in the Xcode integrated development environ-ment (IDE). The plugin implements a fast heuristic constrained linear least squares al-gorithm (CLLS-H) for estimating the diffusion and kurtosis tensors, and offers the user the possibility to choose which maps are to be generated for not only standard DTI quantities such as Mean Diffusion (MD), Radial Diffusion (RD), Axial Diffusion (AD) and Fractional Anisotropy (FA), but also DKI metrics, Mean Kurtosis (MK), Radial Kurtosis (RK) and Axial Kurtosis (AK).The plugin was subjected to both a qualitative and a semi-quantitative analysis which yielded convincing results. A more accurate validation pro-cess is still being developed, after which, and with some few minor adjust-ments the plugin shall become a valid option for DKI computation

Active Microwave Imaging for Breast Cancer Detection

Active microwave imaging is explored as an imaging modality for early detection of breast cancer. When exposed to microwaves, breast tumor exhibits electrical properties that are significantly different from that of healthy breast tissues. The two approaches of active microwave imaging — confocal microwave technique with measured reflected signals and microwave tomographic imaging with measured scattered signals are addressed here. Normal and malignant breast tissue samples of same person are subjected to study within 30 minutes of mastectomy. Corn syrup is used as coupling medium, as its dielectric parameters show good match with that of the normal breast tissue samples. As bandwidth of the transmitter is an important aspect in the time domain confocal microwave imaging approach, wideband bowtie antenna having 2:1 VSWR bandwidth of 46% is designed for the transmission and reception of microwave signals. Same antenna is used for microwave tomographic imaging too at the frequency of 3000 MHz. Experimentally obtained time domain results are substantiated by finite difference time domain (FDTD) analysis. 2-D tomographic images are reconstructed with the collected scattered data using distorted Born iterative method. Variations of dielectric permittivity in breast samples are distinguishable from the obtained permittivity profiles.

A direct D-bar reconstruction algorithm for recovering a complex conductivity in 2D

A direct reconstruction algorithm for complex conductivities in W-2,W-infinity(Omega), where Omega is a bounded, simply connected Lipschitz domain in R-2, is presented. The framework is based on the uniqueness proof by Francini (2000 Inverse Problems 6 107-19), but equations relating the Dirichlet-to-Neumann to the scattering transform and the exponentially growing solutions are not present in that work, and are derived here. The algorithm constitutes the first D-bar method for the reconstruction of conductivities and permittivities in two dimensions. Reconstructions of numerically simulated chest phantoms with discontinuities at the organ boundaries are included.

Microwave imaging of high-contrast objects

Imaging microwave reconstruction dielectric contrast regularization iterative multiport cavity measurement

A fast Bayesian reconstruction algorithm for emission tomography with entropy prior converging to feasible images

The development and tests of an iterative reconstruction algorithm for emission tomography based on Bayesian statistical concepts are described. The algorithm uses the entropy of the generated image as a prior distribution, can be accelerated by the choice of an exponent, and converges uniformly to feasible images by the choice of one adjustable parameter. A feasible image has been defined as one that is consistent with the initial data (i.e. it is an image that, if truly a source of radiation in a patient, could have generated the initial data by the Poisson process that governs radioactive disintegration). The fundamental ideas of Bayesian reconstruction are discussed, along with the use of an entropy prior with an adjustable contrast parameter, the use of likelihood with data increment parameters as conditional probability, and the development of the new fast maximum a posteriori with entropy (FMAPE) Algorithm by the successive substitution method. It is shown that in the maximum likelihood estimator (MLE) and FMAPE algorithms, the only correct choice of initial image for the iterative procedure in the absence of a priori knowledge about the image configuration is a uniform field.

Development of active microwave imaging techniques for applications in mammography

Breast cancer is the most common non - skin malignancy in women and a leading cause of female morality. A potentially important strategy for reducing this menace is the detection at an early stage . The invention of non-invasive and non-ionizing microwave technique, to reveal the internal structure of biological objects was a break through in the field of medical diagnostics. Electrical properties of biological tissues and their interaction with electromagmetic waves have direct impact on human life. This thesis focuses on theoretical and experimental investigations of active microwave imaging techniques for breast cancer detection.

Deciduous tree reconstruction algorithm based on cylinder fitting from mobile terrestrial laser scanned point clouds

Vector reconstruction of objects from an unstructured point cloud obtained with a LiDAR-based system (light detection and ranging) is one of the most promising methods to build three dimensional models of orchards. The cylinder fitting method for woody structure reconstruction of leafless trees from point clouds obtained with a mobile terrestrial laser scanner (MTLS) has been analysed. The advantage of this method is that it performs reconstruction in a single step. The most time consuming part of the algorithm is generation of the cylinder direction, which must be recalculated at the inclusion of each point in the cylinder. The tree skeleton is obtained at the same time as the cluster of cylinders is formed. The method does not guarantee a unique convergence and the reconstruction parameter values must be carefully chosen. A balanced processing of clusters has also been defined which has proven to be very efficient in terms of processing time by following the hierarchy of branches, predecessors and successors. The algorithm was applied to simulated MTLS of virtual orchard models and to MTLS data of real orchards. The constraints applied in the method have been reviewed to ensure better convergence and simpler use of parameters. The results obtained show a correct reconstruction of the woody structure of the trees and the algorithm runs in linear logarithmic time

Immediate Locally Advanced Breast Cancer and Chest Wall Reconstruction: Surgical Planning and Reconstruction Strategies with Extended V-Y Latissimus Dorsi Myocutaneous Flap

Background: Surgical resection in locally advanced breast cancer produces large defects that may not be suitable for primary closure. Immediate reconstruction is controversial and presents a complicated scenario for breast surgeons and plastic surgeons. Methods: In this study, a different design was planned for the latissimus dorsi musculocutaneous flap with primary closure in V-Y for the correction of major lesions in the anterior chest wall. Twenty-five patients underwent immediate locally advanced breast cancer reconstruction with a V-Y latissimus dorsi musculocutaneous flap. This flap was raised from adjacent tissue located on the lateral and posterior thoracic region and presented a triangular shape whose base was the lateral aspect of the mastectomy wound. The technique was indicated in patients with large thoracic wounds. Results: Mean follow-up time was 16 months. Closure was obtained in the donor and recipient sites without the use of skin grafts or other more major procedures. Complications occurred in nine patients (36 percent), including dorsal wound dehiscence in five patients and seroma in three. All cases except one were treated by a conservative approach with a good result. No total flap loss was reported. All patients achieved a satisfactory thoracic reconstruction and adequate wound care. Conclusions: The V-Y latissimus dorsi musculocutaneous flap is a reliable technique for immediate locally advanced breast cancer reconstruction. The technique is advantageous because the V-Y design allows primary closure of the chest wound and donor defect. Success depends on patient selection, coordinated planning with the breast cancer surgeon, and careful intraoperative management. (Plast. Reconstr. Surg. 127: 2186, 2011.)