Monte Carlo Analysis and Physics Characterization of a Novel Nanoparticle Detector for Medical and Micro-dosimetry Applications

The outcomes for both (i) radiation therapy and (ii) preclinical small animal radio- biology studies are dependent on the delivery of a known quantity of radiation to a specific and intentional location. Adverse effects can result from these procedures if the dose to the target is too high or low, and can also result from an incorrect spatial distribution in which nearby normal healthy tissue can be undesirably damaged by poor radiation delivery techniques. Thus, in mice and humans alike, the spatial dose distributions from radiation sources should be well characterized in terms of the absolute dose quantity, and with pin-point accuracy. When dealing with the steep spatial dose gradients consequential to either (i) high dose rate (HDR) brachytherapy or (ii) within the small organs and tissue inhomogeneities of mice, obtaining accurate and highly precise dose results can be very challenging, considering commercially available radiation detection tools, such as ion chambers, are often too large for in-vivo use.

In this dissertation two tools are developed and applied for both clinical and preclinical radiation measurement. The first tool is a novel radiation detector for acquiring physical measurements, fabricated from an inorganic nano-crystalline scintillator that has been fixed on an optical fiber terminus. This dosimeter allows for the measurement of point doses to sub-millimeter resolution, and has the ability to be placed in-vivo in humans and small animals. Real-time data is displayed to the user to provide instant quality assurance and dose-rate information. The second tool utilizes an open source Monte Carlo particle transport code, and was applied for small animal dosimetry studies to calculate organ doses and recommend new techniques of dose prescription in mice, as well as to characterize dose to the murine bone marrow compartment with micron-scale resolution.

Hardware design changes were implemented to reduce the overall fiber diameter to <0.9 mm for the nano-crystalline scintillator based fiber optic detector (NanoFOD) system. Lower limits of device sensitivity were found to be approximately 0.05 cGy/s. Herein, this detector was demonstrated to perform quality assurance of clinical 192Ir HDR brachytherapy procedures, providing comparable dose measurements as thermo-luminescent dosimeters and accuracy within 20% of the treatment planning software (TPS) for 27 treatments conducted, with an inter-quartile range ratio to the TPS dose value of (1.02-0.94=0.08). After removing contaminant signals (Cerenkov and diode background), calibration of the detector enabled accurate dose measurements for vaginal applicator brachytherapy procedures. For 192Ir use, energy response changed by a factor of 2.25 over the SDD values of 3 to 9 cm; however a cap made of 0.2 mm thickness silver reduced energy dependence to a factor of 1.25 over the same SDD range, but had the consequence of reducing overall sensitivity by 33%.

For preclinical measurements, dose accuracy of the NanoFOD was within 1.3% of MOSFET measured dose values in a cylindrical mouse phantom at 225 kV for x-ray irradiation at angles of 0, 90, 180, and 270˝. The NanoFOD exhibited small changes in angular sensitivity, with a coefficient of variation (COV) of 3.6% at 120 kV and 1% at 225 kV. When the NanoFOD was placed alongside a MOSFET in the liver of a sacrificed mouse and treatment was delivered at 225 kV with 0.3 mm Cu filter, the dose difference was only 1.09% with use of the 4x4 cm collimator, and -0.03% with no collimation. Additionally, the NanoFOD utilized a scintillator of 11 µm thickness to measure small x-ray fields for microbeam radiation therapy (MRT) applications, and achieved 2.7% dose accuracy of the microbeam peak in comparison to radiochromic film. Modest differences between the full-width at half maximum measured lateral dimension of the MRT system were observed between the NanoFOD (420 µm) and radiochromic film (320 µm), but these differences have been explained mostly as an artifact due to the geometry used and volumetric effects in the scintillator material. Characterization of the energy dependence for the yttrium-oxide based scintillator material was performed in the range of 40-320 kV (2 mm Al filtration), and the maximum device sensitivity was achieved at 100 kV. Tissue maximum ratio data measurements were carried out on a small animal x-ray irradiator system at 320 kV and demonstrated an average difference of 0.9% as compared to a MOSFET dosimeter in the range of 2.5 to 33 cm depth in tissue equivalent plastic blocks. Irradiation of the NanoFOD fiber and scintillator material on a 137Cs gamma irradiator to 1600 Gy did not produce any measurable change in light output, suggesting that the NanoFOD system may be re-used without the need for replacement or recalibration over its lifetime.

For small animal irradiator systems, researchers can deliver a given dose to a target organ by controlling exposure time. Currently, researchers calculate this exposure time by dividing the total dose that they wish to deliver by a single provided dose rate value. This method is independent of the target organ. Studies conducted here used Monte Carlo particle transport codes to justify a new method of dose prescription in mice, that considers organ specific doses. Monte Carlo simulations were performed in the Geant4 Application for Tomographic Emission (GATE) toolkit using a MOBY mouse whole-body phantom. The non-homogeneous phantom was comprised of 256x256x800 voxels of size 0.145x0.145x0.145 mm3. Differences of up to 20-30% in dose to soft-tissue target organs was demonstrated, and methods for alleviating these errors were suggested during whole body radiation of mice by utilizing organ specific and x-ray tube filter specific dose rates for all irradiations.

Monte Carlo analysis was used on 1 µm resolution CT images of a mouse femur and a mouse vertebra to calculate the dose gradients within the bone marrow (BM) compartment of mice based on different radiation beam qualities relevant to x-ray and isotope type irradiators. Results and findings indicated that soft x-ray beams (160 kV at 0.62 mm Cu HVL and 320 kV at 1 mm Cu HVL) lead to substantially higher dose to BM within close proximity to mineral bone (within about 60 µm) as compared to hard x-ray beams (320 kV at 4 mm Cu HVL) and isotope based gamma irradiators (137Cs). The average dose increases to the BM in the vertebra for these four aforementioned radiation beam qualities were found to be 31%, 17%, 8%, and 1%, respectively. Both in-vitro and in-vivo experimental studies confirmed these simulation results, demonstrating that the 320 kV, 1 mm Cu HVL beam caused statistically significant increased killing to the BM cells at 6 Gy dose levels in comparison to both the 320 kV, 4 mm Cu HVL and the 662 keV, 137Cs beams.

Carbon fluxes and pelagic ecosystem dynamics near two western Antarctic Peninsula Adélie penguin colonies: an inverse model approach

An inverse food-web model for the western Antarctic Peninsula (WAP) pelagic food web was constrained with data from Palmer Long Term Ecological Research (PAL-LTER) project annual austral summer sampling cruises. Model solutions were generated for 2 regions with Adelie penguin Pygoscelis adeliae colonies presenting different population trends (a northern and a southern colony) for a 12 yr period (1995-2006). Counter to the standard paradigm, comparisons of carbon flow through bacteria, microzooplankton, and krill showed that the diatom-krill-top predator food chain is not the dominant pathway for organic carbon exchanges. The food web is more complex, including significant contributions by microzooplankton and the microbial loop. Using both inverse model results and network indices, it appears that in the northern WAP the food web is dominated by the microbial food web, with a temporal trend toward its increasing importance. The dominant pathway for the southern WAP food web varies from year to year, with no detectable temporal trend toward dominance of microzooplankton versus krill. In addition, sensitivity analyses indicated that the northern colony of Adelie penguins, whose population size has been declining over the past 35 yr, appears to have sufficient krill during summer to sustain its basic metabolic needs and rear chicks, suggesting the importance of other processes in regulating the Adelie population decline.

Zooplankton associations in a Mediterranean long-term time-series

Coastal zooplankton have been investigated since 1984 at a Long Term Ecological Research station MC (LTER-MC) in the inner Gulf of Naples (Tyrrhenian Sea, Western Mediterranean). The sampling site, located between the littoral and the open sea systems, has very active hydrography that affects plankton communities. The present work was aimed at establishing whether, in such a dynamic and variable environment, species associations and homogeneous periods could be identified as characteristic and stable features of the mesozooplankton over the period 1984–2006. Hierarchical clustering was applied to assess species associations based on a matrix of similarities between species (R-mode), and homogeneous periods based on a matrix of similarities between observations (Q-mode). The Indicator Value index [IndVal, Dufrene and Legendre (1997) Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol. Monogr., 67, 345–366] was calculated to identify species characterizing each period. Five taxonomic groups with well-defined composition and abundance were identified as robust associations that likely reflect different modes of community functioning. The temporal course of these associations was largely shaped by strong seasonal forcing comprising both physical and biological (e.g. trophic) signals. These associations persisted over the long term, thus indicating some stable characters in the Naples zooplankton time-series, providing evidence of resilience in communities in highly variable coastal conditions.

Oithona similislikes it cool: evidence from two long-term time series

We compare the long-term and seasonal patterns of abundance and phenology of the cyclopoid copepod Oithona similis at the L4 site (1988–2013) in the North Atlantic and at the LTER-MC site (1984–2013) in the Mediterranean Sea to investigate whether high temperature limits the occurrence of this species with latitudinal cline. The two sites are well suited to testing this hypothesis as they are characterized by similar chlorophyll a concentration (Chl a) but different temperature [sea surface temperature (SST)]. The abundance of O. similis at L4 was ∼10 times higher than at LTER-MC. Moreover, this species had several peaks of abundance during the year at L4 but a single peak in spring at LTER-MC. The main mode of temporal variability in abundance was seasonal at both sites. The abundance of O. similis was negatively correlated with SST only at LTER-MC, whereas it was positively correlated with Chl a at both sites. Oithona similis had a temperature optimum between 15 and 20°C reaching maximum abundance at ∼16.5°C at LTER-MC, but showed no Chl a optimum at either site. We conclude that the abundance of O. similis increases with prey availability up to 16.5°C and that temperature >20°C represents the main limiting factor for population persistence.

Metronidazole Population Pharmacokinetics in Preterm Neonates Using Dried Blood-Spot Sampling

OBJECTIVES: To characterize the population pharmacokinetics of metronidazole in preterm neonates.
PATIENTS AND METHODS: Data were collected prospectively from 32 preterm neonates who received intravenous metronidazole for the treatment of or prophylaxis against necrotizing enterocolitis. Dried
blood spots (n 203) on ?lter paper were analyzed by highperformance liquid chromatography, and the data were subjected to pharmacokinetic analysis performed by using nonlinear mixed-effect modeling.
RESULTS: A 1-compartment model best described the data. Signi?cant covariates were weight (WT) and postmenstrual age (PMA). The ?nal population models for metronidazole clearance (CL) and volume of distribution (V) were: CL 0.0247 (WT/1.00)0.75 (1 0.107 [PMA 30]) and V 0.726 WT, where CL is in liters per hour, WT is in kilograms, PMA is in weeks, and V is in liters. This model predicts that the half-life of metronidazole decreases rapidly from 40 hours at 25 weeks’ PMA to 19 hours at 32 weeks’ PMA, after which it starts to plateau. This decrease in half-life is the result of a 5-fold increase in CL compared with only a 2.5-fold increase in V during the same period.
CONCLUSIONS: Currently, there are no speci?c dose recommendations for metronidazole in preterm neonates. However, a dosing scheme for preterm neonates that takes into consideration both the weight and PMA has been suggested and should avoid administration of doses that are excessive or more frequent than necessary.

A Frequency Domain Algorithm for Wave Separation

We propose a frequency domain adaptive algorithm for
wave separation in wind instruments. Forward and backward travelling waves are obtained from the signals acquired by two microphones placed along the tube, while the
separation ?lter is adapted from the information given by a
third microphone. Working in the frequency domain has a
series of advantages, among which are the ease of design of
the propagation ?lter and its differentiation with respect to
its parameters.
Although the adaptive algorithm was developed as a ?rst
step for the estimation of playing parameters in wind instruments it can also be used, without any modi?cations, for
other applications such as in-air direction of arrival (DOA)
estimation. Preliminary results on these applications will
also be presented.

Web Data Extraction from Query Result Pages Based on Visual and Content Features

A rapidly increasing number of Web databases are now become accessible via
their HTML form-based query interfaces. Query result pages are dynamically generated
in response to user queries, which encode structured data and are displayed for human
use. Query result pages usually contain other types of information in addition to query
results, e.g., advertisements, navigation bar etc. The problem of extracting structured data
from query result pages is critical for web data integration applications, such as comparison
shopping, meta-search engines etc, and has been intensively studied. A number of approaches
have been proposed. As the structures of Web pages become more and more complex, the
existing approaches start to fail, and most of them do not remove irrelevant contents which
may a®ect the accuracy of data record extraction. We propose an automated approach for
Web data extraction. First, it makes use of visual features and query terms to identify data
sections and extracts data records in these sections. We also represent several content and
visual features of visual blocks in a data section, and use them to ¯lter out noisy blocks.
Second, it measures similarity between data items in di®erent data records based on their
visual and content features, and aligns them into di®erent groups so that the data in the
same group have the same semantics. The results of our experiments with a large set of
Web query result pages in di®erent domains show that our proposed approaches are highly
e®ective.

Ecological quality assessment in transitional systems

Estuaries are poles of attraction for human settlement which is a source of pressures to surface water bodies. The implementation of the European Water Framework Directive (WDF, 2000/60/EC) has increased the investigation in order to develop methodologies to assess the Ecological Quality Status (EQS) of aquatic ecosystems. Transitional systems are naturally stressed and characterized by highly dynamic physical, chemical and hydro-morphologic conditions and by species with a higher level of tolerance to change, being more difficult to develop suitable quality indicators for these systems. The general purpose of this study is to test the ability of synthesis descriptors, including primary (S, taxa richness) and derived biological variable (H’, Shannon-Wiener diversity), biotic indices (AMBI and M-AMBI), body size properties (abundance distribution by body size classes, length, weight and length-weight relationships) and non-taxonomic indices (ISS), as well as functional indicators related to the decomposition rates of various experimental substrates, a macrophyte (Phragmites australis) and an alga (Fucus vesiculosus), to evaluate the environmental quality in transitional systems. This study was carried out in one of the most pristine channels of the Ria the Aveiro, Mira Channel, along a full salinity gradient and in a metals and metalloid sediment contamination area, the Estarreja Channel, and two reference channels (Canelas and Salreu). In this study were used different sampling techniques, the leaf-bag technique and a hand-held corer. In Mira Channel, the alga and the macrophyte presented an opposite trend in the decomposition rate along the salinity gradient, with the decomposition rates of the alga always higher than those of the macrophyte. The decomposition rates of the macrophyte and the alga were higher in the mid estuary and in higher salinity areas, respectively, corresponding to the preferencial distribution areas of each species. The macrobenthic fauna associated with the decaying and an artificial substrate (control) showed equally well the benthic succession from the marine to the freshwater areas and, despite the strong differences in the decay rates, no significant differences were found between the benthic communities associated with the alga and the macrophyte. The body size properties of the macrobenthic fauna associated with the P. australis leaf-bag (1mm and 5mm) and corer samples were studied along the full salinity gradient. The dominant species of the sub-set of measured specimens were not the same of the original macrobenthic fauna sampled but, despite that, the sub-set of measured specimens was also able to show the benthic succession from the marine to the freshwater areas. The body size abundance distribution of the benthic macroinvertebrates according to the ISS size classes did not show a particular trend in any sampler along the salinity gradient. Significant differences were found in the length, weight and length-weight relationships of Annelids, , Molluscs and even some species along the salinity gradient. No significant differences were found in the AMBI, M-AMBI and ISS values along the salinity gradient for all the samplers. The EQS of the corer samples obtained using the M-AMBI was lower than that of the leaf-bags. The EQS obtained with the ISS was higher than that obtained with the M-AMBI in the leaf-bags but not in the corer samples. The ecological effects of contaminated sediments associated with the industrial chemical effluents discharged in the Estarreja Channel were studied a decade after ceasing the emissions, using the Sediment Quality Triad approach and two reference channels. The results showed that despite the emissions ceased in 2004, the sediment remains polluted with high levels of metals and metalloid, available to bioaccumulation and with severe consequences at the community level. The sediment contamination problem was also studied using the leaf-bag technique with a macrophyte, an alga and a control substrate. The results showed that the decay rates, the associated macrofauna and the application of the AMBI, M-AMBI and ISS indices to the mesh-bag samples were not able to identify the sediment contamination. Contrarily to the AMBI, the M-AMBI and the ISS showed significant differences between the contaminated and the reference channels for the corer samples. Although such statistical significance, the interest of using these complex biotic indices could be questioned, when much simple ones, like the S and H’ allow to reach the same conclusions.

A time frequency approach to blind deconvolution in multipath underwater channels

Blind deconvolution is studied in the underwater acoustic channel context, by time-frequency (TF) processing. The acoustic propagation environment is modelled by ray tracing and mathematically described by a multipath propagation channel. Representation of the received signal by means of a signal-dependent TF distribution (radially Gaussian kernel distribution) allowed to visualize the resolved replicas of the emitted signal, while signi cantly attenuating the inherent interferences of classic quadratic TF distributions. The source signal instantaneous frequency estimation was the starting point for both source and channel estimation. Source signature estimation was performed by either TF inversion, based on the Wigner-Ville distribution of the received signal, or a subspace- -based method. The channel estimate was obtained either via a TF formulation of the conventional matched- lter, or via matched- - ltering with the previously obtained source estimate. A shallow water realistic scenario is considered, comprising a 135-m depth water column and an acoustic source located at 90-m depth and 5.6-km range from the receiver. For the corresponding noiseless simulated data, the quality of the best estimates was 0.856 for the source signal, and 0.9664 and 0.9996 for the amplitudes and time-delays of the impulse response, respectively. Application of the proposed deconvolution method to real data of the INTIMATE '96 sea trial conduced to source and channel estimates with the quality of 0.530 and 0.843, respectively. TF processing has proved to remove the typical ill-conditioning of single sensor deterministic deconvolution techniques.

Soll man in der klinischen Praxis Frailty abschätzen?

Frailty, ein geriatrisches Syndrom mit altersbedingter Gebrechlichkeit assoziiert, heisst Verlust physiologischer Reserven verschiedener Organe. Die Folge ist eine erhöhte Verletzlichkeit durch Stress. Eine allseits anerkannte Definition fehlt, und die Abschätzung erfolgt hauptsächlich aufgrund zweier Modelle: des Phänotyps nach Fried und des Index nach Rockwood. Frailty bringt ein deutlich erhöhtes Risiko für funktionelle Hilfebedürftigkeit, Hospitalisation oder Einweisung in eine Institution für Langzeitpflege sowie Tod mit sich. Der Spontanverlauf ist progressiv, kann aber auch reversibel sein, und daher könnte Prävention grundsätzlich in Frage kommen. Unsere Kenntnisse sind derzeit aber immer noch zu lückenhaft, als dass wir gezielt gegen Frailty angehen könnten. Klinische Forschung über Frailty wird erst seit kurzem betrieben. Die Empfehlung, den Frailtygrad durch Leistungsmessung (wie Ganggeschwindigkeit, Griffkraft etc.) zu evaluieren, stützt sich auf die Beobachtung eines Zusammenhangs mit einer späteren ungünstigen Entwicklung, und nicht darauf, dass wir wirksam eingreifen könnten. Klinische und epidemiologische Studien über Frailty sind wichtig. Da die Babyboomer nun älter werden, ist die Gefahr einer Epidemie funktioneller Hilfebedürftigkeit Mitte des Jahrhunderts absehbar, und somit stellt das Problem der Frailty eine grosse Herausforderung dar.

High Selectivity Filter Employing Stepped Impedance Resonators, Series Capacitors and Defected Ground Structures for Ultra Wide Band Applications

The paper presents a compact planar Ultra Wide Band ¯lter employing folded stepped impedance resonators with series capacitors and dumb bell shaped defected ground structures. An interdigital quarter wavelength coupled line is used for achieving the band pass characteristics. The transmission zeros are produced by stepped impedance resonators. The ¯lter has steep roll o® rate and good attenuation in its lower and upper stop bands, contributed by the series capacitor and defected ground structures respectively.

Characterization of Novel Fabry Pérot Filter Arrays for Nanospectrometers in Medical Applications

Optische Spektroskopie ist eine sehr wichtige Messtechnik mit einem hohen Potential für zahlreiche Anwendungen in der Industrie und Wissenschaft. Kostengünstige und miniaturisierte Spektrometer z.B. werden besonders für moderne Sensorsysteme “smart personal environments” benötigt, die vor allem in der Energietechnik, Messtechnik, Sicherheitstechnik (safety and security), IT und Medizintechnik verwendet werden. Unter allen miniaturisierten Spektrometern ist eines der attraktivsten Miniaturisierungsverfahren das Fabry Pérot Filter. Bei diesem Verfahren kann die Kombination von einem Fabry Pérot (FP) Filterarray und einem Detektorarray als Mikrospektrometer funktionieren. Jeder Detektor entspricht einem einzelnen Filter, um ein sehr schmales Band von Wellenlängen, die durch das Filter durchgelassen werden, zu detektieren. Ein Array von FP-Filter wird eingesetzt, bei dem jeder Filter eine unterschiedliche spektrale Filterlinie auswählt. Die spektrale Position jedes Bandes der Wellenlänge wird durch die einzelnen Kavitätshöhe des Filters definiert. Die Arrays wurden mit Filtergrößen, die nur durch die Array-Dimension der einzelnen Detektoren begrenzt werden, entwickelt. Allerdings erfordern die bestehenden Fabry Pérot Filter-Mikrospektrometer komplizierte Fertigungsschritte für die Strukturierung der 3D-Filter-Kavitäten mit unterschiedlichen Höhen, die nicht kosteneffizient für eine industrielle Fertigung sind. Um die Kosten bei Aufrechterhaltung der herausragenden Vorteile der FP-Filter-Struktur zu reduzieren, wird eine neue Methode zur Herstellung der miniaturisierten FP-Filtern mittels NanoImprint Technologie entwickelt und präsentiert. In diesem Fall werden die mehreren Kavitäten-Herstellungsschritte durch einen einzigen Schritt ersetzt, die hohe vertikale Auflösung der 3D NanoImprint Technologie verwendet. Seit dem die NanoImprint Technologie verwendet wird, wird das auf FP Filters basierende miniaturisierte Spectrometer nanospectrometer genannt. Ein statischer Nano-Spektrometer besteht aus einem statischen FP-Filterarray auf einem Detektorarray (siehe Abb. 1). Jeder FP-Filter im Array besteht aus dem unteren Distributed Bragg Reflector (DBR), einer Resonanz-Kavität und einen oberen DBR. Der obere und untere DBR sind identisch und bestehen aus periodisch abwechselnden dünnen dielektrischen Schichten von Materialien mit hohem und niedrigem Brechungsindex. Die optischen Schichten jeder dielektrischen Dünnfilmschicht, die in dem DBR enthalten sind, entsprechen einen Viertel der Design-Wellenlänge. Jeder FP-Filter wird einer definierten Fläche des Detektorarrays zugeordnet. Dieser Bereich kann aus einzelnen Detektorelementen oder deren Gruppen enthalten. Daher werden die Seitenkanal-Geometrien der Kavität aufgebaut, die dem Detektor entsprechen. Die seitlichen und vertikalen Dimensionen der Kavität werden genau durch 3D NanoImprint Technologie aufgebaut. Die Kavitäten haben Unterschiede von wenigem Nanometer in der vertikalen Richtung. Die Präzision der Kavität in der vertikalen Richtung ist ein wichtiger Faktor, der die Genauigkeit der spektralen Position und Durchlässigkeit des Filters Transmissionslinie beeinflusst.

Efficient topology estimation for large scale optical mapping

Large scale image mosaicing methods are in great demand among scientists who study different aspects of the seabed, and have been fostered by impressive advances in the capabilities of underwater robots in gathering optical data from the seafloor. Cost and weight constraints mean that lowcost Remotely operated vehicles (ROVs) usually have a very limited number of sensors. When a low-cost robot carries out a seafloor survey using a down-looking camera, it usually follows a predetermined trajectory that provides several non time-consecutive overlapping image pairs. Finding these pairs (a process known as topology estimation) is indispensable to obtaining globally consistent mosaics and accurate trajectory estimates, which are necessary for a global view of the surveyed area, especially when optical sensors are the only data source. This thesis presents a set of consistent methods aimed at creating large area image mosaics from optical data obtained during surveys with low-cost underwater vehicles. First, a global alignment method developed within a Feature-based image mosaicing (FIM) framework, where nonlinear minimisation is substituted by two linear steps, is discussed. Then, a simple four-point mosaic rectifying method is proposed to reduce distortions that might occur due to lens distortions, error accumulation and the difficulties of optical imaging in an underwater medium. The topology estimation problem is addressed by means of an augmented state and extended Kalman filter combined framework, aimed at minimising the total number of matching attempts and simultaneously obtaining the best possible trajectory. Potential image pairs are predicted by taking into account the uncertainty in the trajectory. The contribution of matching an image pair is investigated using information theory principles. Lastly, a different solution to the topology estimation problem is proposed in a bundle adjustment framework. Innovative aspects include the use of fast image similarity criterion combined with a Minimum spanning tree (MST) solution, to obtain a tentative topology. This topology is improved by attempting image matching with the pairs for which there is the most overlap evidence. Unlike previous approaches for large-area mosaicing, our framework is able to deal naturally with cases where time-consecutive images cannot be matched successfully, such as completely unordered sets. Finally, the efficiency of the proposed methods is discussed and a comparison made with other state-of-the-art approaches, using a series of challenging datasets in underwater scenarios