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.

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.

Observations of middle atmospheric H2O and O3 during the 2010 major sudden stratospheric warming by a network of microwave radiometers

In this study, we present middle atmospheric water vapor (H2O) and ozone (O3) measurements obtained by ground-based microwave radiometers at three European locations in Bern (47° N), Onsala (57° N) and Sodankylä (67° N) during Northern winter 2009/2010. In January 2010, a major sudden stratospheric warming (SSW) occurred in the Northern Hemisphere whose signatures are evident in the ground-based observations of H2O and O3. The observed anomalies in H2O and O3 are mostly explained by the relative location of the polar vortex with respect to the measurement locations. The SSW started on 26 January 2010 and was most pronounced by the end of January. The zonal mean temperature in the middle stratosphere (10 hPa) increased by approximately 25 Kelvin within a few days. The stratospheric vortex weakened during the SSW and shifted towards Europe. In the mesosphere, the vortex broke down, which lead to large scale mixing of polar and midlatitudinal air. After the warming, the polar vortex in the stratosphere split into two weaker vortices and in the mesosphere, a new, pole-centered vortex formed with maximum wind speed of 70 m s−1 at approximately 40° N. The shift of the stratospheric vortex towards Europe was observed in Bern as an increase in stratospheric H2O and a decrease in O3. The breakdown of the mesospheric vortex during the SSW was observed at Onsala and Sodankylä as a sudden increase in mesospheric H2O. The following large-scale descent inside the newly formed mesospheric vortex was well captured by the H2O observations in Sodankylä. In order to combine the H2O observations from the three different locations, we applied the trajectory mapping technique on our H2O observations to derive synoptic scale maps of the H2O distribution. Based on our observations and the 3-D wind field, this method allows determining the approximate development of the stratospheric and mesospheric polar vortex and demonstrates the potential of a network of ground-based instruments.

Microwave-Assisted Phospha-Michael Reactions with Internal Alkenes

The atom efficient phospha-Michael reaction between bis 4-methylphenyl phosphine oxide and several activated internal alkenes has been shown to occur under microwave irradiation without added solvent or catalyst. The alkenes used for this study were ethyl 4-nitrocinnamate, two chalcones ((E)-3-(4-methoxy-phenyl)-1-(4- nitrophenyl)-prop-2-en-1-one and (E)-1-(4-methoxyphenyl)-3-(3-nitro-phenyl)-prop-2- en-1-one), and 2-phenylmethylene-propanedinitrile. In the case of ethyl 4-nitrocinnamate, reaction with bis 4-methylphenyl phosphine oxide for sixty minutes at 130 °C yielded the desired phospha-Michael product in a 55% yield after purification. Varying the location of the nitro group on the phenyl rings of the chalcones did not seem to have a large effect on their reactivity. By NMR, both chalcones seemed to react to the same extent when the reaction times and temperatures were held constant. Interestingly, a phospha-Michael reaction was observed at a reaction temperature of 65°C for experiments involving 2- phenyl-methylene-propanedinitrile while the other substrates required a reaction temperature of 130 °C. Similar experiments were carried out with bis mesityl phosphine oxide and two internal alkenes: 2-phenylmethylene-propanedinitrile and ethyl-2-cyano-3- methyl-2-butenoate. These experiments did not yield any of the predicted phospha- Michael products, which suggest steric limitations to the Michael donor for this reaction.

Microwave-Assisted Bond Forming Reactions

ABSTRACT FOR PART I: PHOSPHA-MICHAEL ADDITIONS TO ACTIVATED INTERNAL ALKENES: STERIC AND ELECTRONIC EFFECTS A method for the phospha-Michael addition of bis(4-methyl)phenyl phosphine oxide to activated internal alkenes has been developed. Michael acceptors including cinnamates, crotonates, chalcones, and internal alkenes containing multiple activating groups were all successfully utilized in this reaction. The reaction was fairly tolerant of electron-donating and electron-withdrawing substituents on the Michael acceptor, and moderate to excellent yields (49-96%) of the adducts were isolated. When steric bulk was increased by a second substituent on the -position of the Michael-acceptor the reaction was suppressed. This was usually overcome by adding a second activating substituent to the -position. ABSTRACT FOR PART II: MICROWAVE-ASSISTED ARYLGOLD BOND FORMATION A microwave-assisted method was developed for the formation of arylgold complexes containing (2-Biphenyl)di-tert-butylphosphine (JohnPhos) as the supporting phosphine ligand. Arylboronic acids with increasingly bulky aromatic groups were screened to determine the steric limitations of the reaction. Arylgold complexes (JohnPhos)Au(p-methoxyphenyl), (JohnPhos)Au(2,4,6-trimethylphenyl), and (JohnPhos)Au(4-bromo-10-anthracene) were all synthesized by microwave irradiation at 70ºC in the presence of Cs2CO3 in either THF or iPrOH. Reactions performed with arylboronic acids containing unhindered ortho positions were carried out in THF. Arylboronic acids with substituents on the ortho position required iPrOH as the reaction solvent. Arylboronic acids with extreme steric hindrance on the ortho position of the aryl substituent, 2,4,6-triisopropylpphenylboronic acid, were unreactive. It was determined that increasing the irradiation temperature increased the formation of side products, therefore to promote conversion to the arylgold complex the duration of the reaction time was increased while maintaining a temperature of 70ºC. Arylgold complexes (JohnPhos)Au(p-methoxyphenyl), (JohnPhos)Au(2,4,6-trimethylphenyl), and (JohnPhos)Au(4-bromo-10-anthracene) were synthesized with moderate yields (40-69%).