Spectrum reconstruction from a scattering measurement using the adjoint Boltzmann transport equation for photons

Quality control of medical radiological systems is of fundamental importance, and requires efficient methods for accurately determine the X-ray source spectrum. Straightforward measurements of X-ray spectra in standard operating require the limitation of the high photon flux, and therefore the measure has to be performed in a laboratory. However, the optimal quality control requires frequent in situ measurements which can be only performed using a portable system. To reduce the photon flux by 3 magnitude orders an indirect technique based on the scattering of the X-ray source beam by a solid target is used. The measured spectrum presents a lack of information because of transport and detection effects. The solution is then unfolded by solving the matrix equation that represents formally the scattering problem. However, the algebraic system is ill-conditioned and, therefore, it is not possible to obtain a satisfactory solution. Special strategies are necessary to circumvent the ill-conditioning. Numerous attempts have been done to solve this problem by using purely mathematical methods. In this thesis, a more physical point of view is adopted. The proposed method uses both the forward and the adjoint solutions of the Boltzmann transport equation to generate a better conditioned linear algebraic system. The procedure has been tested first on numerical experiments, giving excellent results. Then, the method has been verified with experimental measurements performed at the Operational Unit of Health Physics of the University of Bologna. The reconstructed spectra have been compared with the ones obtained with straightforward measurements, showing very good agreement.

High-resolution reconstruction of Holocene climate variability and environmental changes in the North Pacific using bivalve shells

Bivalve mollusk shells are useful tools for multi-species and multi-proxy paleoenvironmental reconstructions with a high temporal and spatial resolution. Past environmental conditions can be reconstructed from shell growth and stable oxygen and carbon isotope ratios, which present an archive for temperature, freshwater fluxes and primary productivity. The purpose of this thesis is the reconstruction of Holocene climate and environmental variations in the North Pacific with a high spatial and temporal resolution using marine bivalve shells. This thesis focuses on several different Holocene time periods and multiple regions in the North Pacific, including: Japan, Alaska (AK), British Columbia (BC) and Washington State, which are affected by the monsoon, Pacific Decadal Oscillation (PDO) and El Niño/Southern Oscillation (ENSO). Such high-resolution proxy data from the marine realm of mid- and high-latitudes are still rare. Therefore, this study contributes to the optimization and verification of climate models. However, before using bivalves for environmental reconstructions and seasonality studies, life history traits must be well studied to temporally align and interpret the geochemical record. These calibration studies are essential to ascertain the usefulness of selected bivalve species as paleoclimate proxy archives. This work focuses on two bivalve species, the short-lived Saxidomus gigantea and the long-lived Panopea abrupta. Sclerochronology and oxygen isotope ratios of different shell layers of P. abrupta were studied in order to test the reliability of this species as a climate archive. The annual increments are clearly discernable in umbonal shell portions and the increments widths should be measured in these shell portions. A reliable reconstruction of paleotemperatures may only be achieved by exclusively sampling the outer shell layer of multiple contemporaneous specimens. Life history traits (e.g., timing of growth line formation, duration of the growing season and growth rates) and stable isotope ratios of recent S. gigantea from AK and BC were analyzed in detail. Furthermore, a growth-temperature model based on S. gigantea shells from Alaska was established, which provides a better understanding of the hydrological changes related to the Alaska Coastal Current (ACC). This approach allows the independent measurement of water temperature and salinity from variations in the width of lunar daily growth increments of S. gigantea. Temperature explains 70% of the variability in shell growth. The model was calibrated and tested with modern shells and then applied to archaeological specimens. The time period between 988 and 1447 cal yrs BP was characterized by colder (~1-2°C) and much drier (2-5 PSU) summers, and a likely much slower flowing ACC than at present. In contrast, the summers during the time interval of 599-1014 cal yrs BP were colder (up to 3°C) and fresher (1-2 PSU) than today. The Aleutian Low may have been stronger and the ACC was probably flowing faster during this time.

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn

Oblique ionograms automatic scaling and eikonal based ray tracing

A method for automatic scaling of oblique ionograms has been introduced. This method also provides a rejection procedure for ionograms that are considered to lack sufficient information, depicting a very good success rate. Observing the Kp index of each autoscaled ionogram, can be noticed that the behavior of the autoscaling program does not depend on geomagnetic conditions. The comparison between the values of the MUF provided by the presented software and those obtained by an experienced operator indicate that the procedure developed for detecting the nose of oblique ionogram traces is sufficiently efficient and becomes much more efficient as the quality of the ionograms improves. These results demonstrate the program allows the real-time evaluation of MUF values associated with a particular radio link through an oblique radio sounding. The automatic recognition of a part of the trace allows determine for certain frequencies, the time taken by the radio wave to travel the path between the transmitter and receiver. The reconstruction of the ionogram traces, suggests the possibility of estimating the electron density between the transmitter and the receiver, from an oblique ionogram. The showed results have been obtained with a ray-tracing procedure based on the integration of the eikonal equation and using an analytical ionospheric model with free parameters. This indicates the possibility of applying an adaptive model and a ray-tracing algorithm to estimate the electron density in the ionosphere between the transmitter and the receiver An additional study has been conducted on a high quality ionospheric soundings data set and another algorithm has been designed for the conversion of an oblique ionogram into a vertical one, using Martyn's theorem. This allows a further analysis of oblique soundings, throw the use of the INGV Autoscala program for the automatic scaling of vertical ionograms.

The reconstruction of skeletal movement: the soft tissue artefact issue

In 3D human movement analysis performed using stereophotogrammetric systems and skin markers, bone pose can only be estimated in an indirect fashion. During a movement, soft tissue deformations make the markers move with respect to the underlying bone generating soft tissue artefact (STA). STA has devastating effects on bone pose estimation and its compensation remains an open question. The aim of this PhD thesis was to contribute to the solution of this crucial issue. Modelling STA using measurable trial-specific variables is a fundamental prerequisite for its removal from marker trajectories. Two STA model architectures are proposed. Initially, a thigh marker-level artefact model is presented. STA was modelled as a linear combination of joint angles involved in the movement. This model was calibrated using ex-vivo and in-vivo STA invasive measures. The considerable number of model parameters led to defining STA approximations. Three definitions were proposed to represent STA as a series of modes: individual marker displacements, marker-cluster geometrical transformations (MCGT), and skin envelope shape variations. Modes were selected using two criteria: one based on modal energy and another on the selection of modes chosen a priori. The MCGT allows to select either rigid or non-rigid STA components. It was also empirically demonstrated that only the rigid component affects joint kinematics, regardless of the non-rigid amplitude. Therefore, a model of thigh and shank STA rigid component at cluster-level was then defined. An acceptable trade-off between STA compensation effectiveness and number of parameters can be obtained, improving joint kinematics accuracy. The obtained results lead to two main potential applications: the proposed models can generate realistic STAs for simulation purposes to compare different skeletal kinematics estimators; and, more importantly, focusing only on the STA rigid component, the model attains a satisfactory STA reconstruction with less parameters, facilitating its incorporation in an pose estimator.

Scene Reconstruction And Understanding By RGB-D Sensors

This thesis investigates interactive scene reconstruction and understanding using RGB-D data only. Indeed, we believe that depth cameras will still be in the near future a cheap and low-power 3D sensing alternative suitable for mobile devices too. Therefore, our contributions build on top of state-of-the-art approaches to achieve advances in three main challenging scenarios, namely mobile mapping, large scale surface reconstruction and semantic modeling. First, we will describe an effective approach dealing with Simultaneous Localization And Mapping (SLAM) on platforms with limited resources, such as a tablet device. Unlike previous methods, dense reconstruction is achieved by reprojection of RGB-D frames, while local consistency is maintained by deploying relative bundle adjustment principles. We will show quantitative results comparing our technique to the state-of-the-art as well as detailed reconstruction of various environments ranging from rooms to small apartments. Then, we will address large scale surface modeling from depth maps exploiting parallel GPU computing. We will develop a real-time camera tracking method based on the popular KinectFusion system and an online surface alignment technique capable of counteracting drift errors and closing small loops. We will show very high quality meshes outperforming existing methods on publicly available datasets as well as on data recorded with our RGB-D camera even in complete darkness. Finally, we will move to our Semantic Bundle Adjustment framework to effectively combine object detection and SLAM in a unified system. Though the mathematical framework we will describe does not restrict to a particular sensing technology, in the experimental section we will refer, again, only to RGB-D sensing. We will discuss successful implementations of our algorithm showing the benefit of a joint object detection, camera tracking and environment mapping.

The reconstruction of paleoclimate variability in Northern Germany, during marine isotope stage 5

In order to obtain a better understanding about the influence of post-depositional diagenesis on speleothem 230Th/U-ages and paleoclimate variability during Marine Isotope Stage (MIS) 5 in northern Germany, four stalagmites from the Riesenberghöhle (RBH) were investigated by thin section analysis, 230Th/U-dating as well as stable oxygen and carbon isotope and laser ablation inductively coupled mass spectrometry (LA-ICPMS) trace element analysis. The RBH is located in the Weser Hills and is one of the northernmost limestone caves in Germany.rnMulti collector (MC) ICPMS 230Th/U-ages and thin section analysis of the RBH stalagmites shows that some growth phases of the stalagmites were diagenetically altered after their deposition. The impact of post-depositional diagenesis (PDD) on the 230Th/U-ages is modeled, and potential processes leading to PDD are discussed. In this context, it is suggested that PDD may be induced by rapid climate change at the inception of the GIS.rnDespite of the dating uncertainties resulting from PDD, 230Th/U-dating shows that the RBH stalagmites grew during the Eemian and most of the Greenland Interstadials (GIS) during MIS 5. Thus, the growth phases of the RBH stalagmites might be related to a reorganization of the Atlantic Meridional Overturning Circulation (AMOC). The stable isotope (δ13C and δ18O) and the trace element variability of the stalagmites reflects rapid changes of past temperature and precipitation on millennial and sub-millennial timescales. These past climate changes can be amplified by orbitally forced variations of the July solar insolation at 65°N.

Calibration study of bivalve shells : implications for paleoenvironmental reconstruction

The accretionary shells of bivalve mollusks can provide environmental information, such as water temperature, precipitation, freshwater fluxes, primary productivity and anthropogenic activities in the form of variable growth rates and variable geochemical properties, such as stable oxygen and carbon isotopes. However, paleoenvironmental reconstructions are constrained by uncertainties about isotopic equilibrium fractionation during shell formation, which is generally acknowledged as a reasonable assumption for bivalves, but it has been disputed in several species. Furthermore, the variation in shell growth rates is accepted to rely on multiple environmental variables, such as temperature, food availability and salinity, but can differ from species to species. Therefore, it is necessary to perform species-specific calibration studies for both isotope proxies and shell growth rates before they can be used with confidence for environmental interpretations of the past. Accordingly, the principal objective of this Ph.D research is to examine the reliability of selected bivalve species, the long-lived Eurhomalea exalbida (Dillwyn), the short-lived and fast growing species Paphia undulata (Born 1778), and the freshwater mussel Margaritifera falcata (Gould 1850), as paleoenvironmental proxy archives.rnThe first part is focused on δ18Oshell and shell growth history of live-collected E. exalbida from the Falkland Islands. The most remarkable finding, however, is that E. exalbida formed its shell with an offset of -0.48‰ to -1.91‰ from the expected oxygen isotopic equilibrium with the ambient water. If this remained unnoticed, paleotemperature estimates would overestimate actual water temperatures by 2.1-8.3°C. With increasing ontogenetic age, the discrepancy between measured and reconstructed temperatures increased exponentially, irrespective of the seasonally varying shell growth rates. This study clearly demonstrates that, when the disequilibrium fractionation effect is taken into account, E. exalbida can serve as a high-resolution paleoclimate archive for the southern South America. The species therefore provides quantifiable temperature estimates, which yields new insights into long-term paleoclimate dynamics for mid to high latitudes on the southern hemisphere.rnThe stable carbon isotope of biogenic carbonates is generally considered to be useful for reconstruction of seawater dissolved inorganic carbon. The δ13Cshell composition of E. exalbida was therefore, investigated in the second part of this study. This chapter focuses on inter-annual and intra-annual variations in δ13Cshell. Environmental records in δ13Cshell are found to be strongly obscured by changes in shell growth rates, even if removing the ontogenetic decreasing trend. This suggests that δ13Cshell in E. exalbida may not be useful as an environmental proxy, but a potential tool for ecological investigations. rnIn addition to long-lived bivalve species, short-lived species that secrete their shells extremely fast, can also be useful for environmental reconstructions, especially as a high-resolution recorder. Therefore, P. undulata from Daya Bay, South China Sea was utilized in Chapter 4 to evaluate and establish a potential proxy archive for past variations of the East Asian monsoon on shorter time-scales. The δ18Oshell can provide qualitative estimates of the amount of monsoonal rain and terrestrial runoff and the δ13Cshell likely reflect the relative amount of isotopically light terrestrial carbon that reaches the ocean during the summer monsoon season. Therefore, shells of P. undulata can provide serviceable proxy archives to reconstruct the frequency of exceptional summer monsoons in the past. The relative strength of monsoon-related precipitation and associated changes in ocean salinity and the δ13C ratios of the dissolved inorganic carbon signature (δ13CDIC) can be estimated from the δ18Oshell and δ13Cshell values as well as shell growth patterns. rnIn the final part, the freshwater pearl shell M. falcata from four rivers in British Columbia, Canada was preliminarily studied concerning the lifespans and the shell growth rates. Two groups separated by the Georgia Strait can be clearly distinguished. Specimens from the western group exhibit a shorter lifespan, while the eastern group live longer. Moreover, the average lifespan seems to decrease from south to north. The computed growth equations from the eastern and western groups differ as well. The western group exhibits a lower growth rate, while bivalves from the eastern group grow faster. The land use history seems to be responsible for the differences in lifespans of the specimens from the two groups. Differences in growth rate may be induced by differences in water temperature or nutrient input also related to the land use activities.

Paleoclimate reconstruction using LA-ICP-MS analyses of ostracods and stalagmite

Als paläoklimatische Archive können unter anderem Stalagmiten und Klappen von Ostrakoden herangezogen werden. Bisher gab es noch keine veröffentlichten Ergebnisse räumlich hochaufgelöster Spurenelementverteilungen in einzelnen Ostrakodenklappen. Das gleiche gilt für die Bestimmung radiogener Blei-Isotope in Stalagmiten. Um klimatische Prozesse vergangener Zeiten zu untersuchen, wurde eine neue LA-ICP-MS (Laserablations – Massenspektrometrie mit induktiv gekoppeltem Plasma) Technik für hochaufgelöste in-situ Messungen von Spurenelementen und Bleiisotopen entwickelt. rnrnZunächst wurden geeignete Materialien für die Kalibrierung der Technik untersucht; als Proben dienten die Silikatreferenzgläser BAM-S005-A und BAM-S005B der Bundesanstalt für Materialforschung und –prüfung (BAM). Die Homogenität dieser Referenzgläser wurde mit LA-ICP-MS und den anderen mikroanalytischen Methoden EPMA (Elektronenmikrosonde) und SIMS (Ionensonde) überprüft. Die Ergebnisse zeigten, dass alle Haupt- und die meisten Spurenelemente in beiden Gläsern selbst im Mikrometerbereich homogen verteilt sind. Ausnahmen sind einige Spurenelemente, wie Cs, Cl, Cr, Mo und Ni, die heterogen in den Gläsern verteilt sind. Die Hauptelementzusammensetzung von BAM-S005-A und BAM-S005-B wurde mit Hilfe der EPMA bestimmt, wobei die Ergebnisse die Referenzwerte des BAM-Zertifikats bestätigten. Mit Ausnahme von Sr, Ba, Ce und Pb, stimmten die LA-ICP-MS-Spurenelementwerte mit den zertifizierten Werten innerhalb der angegebenen Fehlergrenzen überein. Gründe für die Diskrepanz der vier oben erwähnten Elemente sind noch unklar, aber sind möglicherweise durch fehlerhafte Referenzwerte zu erklären. Zusätzlich wurden 22 Spurenelemente gemessen, deren Gehalte von BAM nicht zertifiziert wurden. Aufgrund dieser Untersuchungen konnte festgestellt werden, dass beide BAM-Gläser für mikroanalytische Anwendungen geeignet sind.rnrnUm neuartige paläoklimatische Proxies im Calcit von Ostrakoden aus tibetanischen Seesedimenten zu untersuchen, wurde die Spurenelementvariabilität in einzelnen Ostrakodenklappen durch eine neue LA-ICP-MS-Technik bestimmt. Klappen von Ostrakoden der drei Arten (Leucocytherella sinensis Huang, 1982, ?Leucocythere dorsotuberosa Huang, 1982 und ?L. dorsotuberosa f. postilirata sensu Pang, 1985) wurden aus zwei Sedimentkernen des Nam Co Sees auf dem Hochplateau von Tibet gewonnen. Zwei LA-ICP-MS-Varianten, Spot- bzw. Linienanalyse, wurden verwendet, um die Elementkonzentrationen der Spurenelemente Mg, Sr, Ba, U und die der Seltenen Erdelemente (SEE) in den Klappen einzelner Ostrakoden zu bestimmen. Die Ergebnisse zeigten, dass die Linienanalyse präzisere Daten als die Spotanalyse liefert und sie wurde daher vorgezogen. Signifikante Unterschiede in der Spurenelementzusammensetzung zwischen den verschiedenen Arten der Ostrakoden wurden nicht gefunden. Variationen der Elementverhältnisse Mg/Ca und Sr/Ca in den Klappen stimmen mit veröffentlichten Seenspiegelschwankungen während des Holozäns überein, was zeigt, dass Mg- und Sr-Messungen in den Ostrakoden zur Untersuchung paläohydrochemischer Prozesse in diesem Gebiet herangezogen werden kann. Die gute Korrelation, die in dieser Arbeit zwischen Ba/Ca und Sr/Ca gefunden wurde, ist ein Hinweis darauf, dass der Einbau von Ba und Sr in die Klappen durch den gleichen Mechanismus erfolgte. Eine mögliche Beziehung zwischen dem U/Ca-Verhältnis in den Ostrakoden und den Redoxbedingungen auf dem Boden des Sees in der Vergangenheit wird diskutiert. Relativ geringe und konstante La/Ca-Verhältnisse wurden festgestellt, deren Ursache möglicherweise entweder auf der SEE-Charakteristik des Seewassers, auf biologischen Prozessen in den Ostrakoden oder auf Kontamination von Fe-Mn und/oder organischen Substanzen beruhen. Weitere Untersuchungen an Proben aus diesem Gebiet, speziell Klappen von lebenden Ostrakoden, sind notwendig, um den Gehalt von Ba, U und den SEE in Ostrakoden als paläoklimatische Proxies von Umweltbedingungen zu verwenden.

Image quality and dose evaluation of filtered back projection versus iterative reconstruction algorithm in multislice computed tomography.

Il presente lavoro di tesi è stato svolto presso il servizio di Fisica Sanitaria del Policlinico Sant'Orsola-Malpighi di Bologna. Lo studio si è concentrato sul confronto tra le tecniche di ricostruzione standard (Filtered Back Projection, FBP) e quelle iterative in Tomografia Computerizzata. Il lavoro è stato diviso in due parti: nella prima è stata analizzata la qualità delle immagini acquisite con una CT multislice (iCT 128, sistema Philips) utilizzando sia l'algoritmo FBP sia quello iterativo (nel nostro caso iDose4). Per valutare la qualità delle immagini sono stati analizzati i seguenti parametri: il Noise Power Spectrum (NPS), la Modulation Transfer Function (MTF) e il rapporto contrasto-rumore (CNR). Le prime due grandezze sono state studiate effettuando misure su un fantoccio fornito dalla ditta costruttrice, che simulava la parte body e la parte head, con due cilindri di 32 e 20 cm rispettivamente. Le misure confermano la riduzione del rumore ma in maniera differente per i diversi filtri di convoluzione utilizzati. Lo studio dell'MTF invece ha rivelato che l'utilizzo delle tecniche standard e iterative non cambia la risoluzione spaziale; infatti gli andamenti ottenuti sono perfettamente identici (a parte le differenze intrinseche nei filtri di convoluzione), a differenza di quanto dichiarato dalla ditta. Per l'analisi del CNR sono stati utilizzati due fantocci; il primo, chiamato Catphan 600 è il fantoccio utilizzato per caratterizzare i sistemi CT. Il secondo, chiamato Cirs 061 ha al suo interno degli inserti che simulano la presenza di lesioni con densità tipiche del distretto addominale. Lo studio effettuato ha evidenziato che, per entrambi i fantocci, il rapporto contrasto-rumore aumenta se si utilizza la tecnica di ricostruzione iterativa. La seconda parte del lavoro di tesi è stata quella di effettuare una valutazione della riduzione della dose prendendo in considerazione diversi protocolli utilizzati nella pratica clinica, si sono analizzati un alto numero di esami e si sono calcolati i valori medi di CTDI e DLP su un campione di esame con FBP e con iDose4. I risultati mostrano che i valori ricavati con l'utilizzo dell'algoritmo iterativo sono al di sotto dei valori DLR nazionali di riferimento e di quelli che non usano i sistemi iterativi.

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.

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.

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.

Vascularised fibula grafts for early tibia reconstruction in infants with congenital pseudarthrosis

Congenital pseudarthrosis of the tibia (CPT) is caused by an ill-defined, segmental disturbance of periosteal bone formation leading to spontaneous bowing, followed by fracture and subsequent pseudarthrosis in the first 2 years of life. The results of conventional treatment modalities (e.g., bracing, internal and external fixation and bone grafting) are associated with high failure rates in terms of persisting pseudarthrosis, malunion and impaired growth. As a more promising alternative, a more aggressive approach, including wide resection of the affected bone, reconstruction with free vascularised fibula grafts from the healthy contralateral leg and stable external fixation at a very early stage has been suggested. Between 1995 and 2007, 10 children (age 12-31 months, median 20 months) suffering from CPT were treated at our institutions according to this principle. Two patients were treated before a fracture had occurred. The length of the fibula graft was 7-9cm. End-to-end anastomoses were performed at the level of the distal tibia stump. The follow-up was 80 months (median, range 12 months to 12 years). Radiologic examination at 6 weeks postoperatively showed normal bone density and structure of the transplanted fibula in all cases and osseous consolidation at 19 of the 20 graft/tibia junctions. One nonunion was sucessfully treated with bone grafting and plate osteosynthesis. Pin-tract infection occurred in three patients. Five children sustained graft fractures that were successfully treated with internal or external fixation. Two patients developed diminished growth of the affected limb or foot; all others had equal limb length and shoe size. At long-term follow-up, tibialisation of the transplant had occurred, and normal gait and physical activities were possible in all children. We conclude that in spite of a relatively high complication rate and the reluctance to perform free flap surgery in infants at this young age, the present concept may successfully prevent the imminent severe sequelae associated with CPT.