A Tutorial on Speckle Reduction in Synthetic Aperture Radar Images

Speckle is a granular disturbance, usually modeled as a multiplicative noise, that affects synthetic aperture radar (SAR) images, as well as all coherent images. Over the last three decades, several methods have been proposed for the reduction of speckle, or despeckling, in SAR images. Goal of this paper is making a comprehensive review of despeckling methods since their birth, over thirty years ago, highlighting trends and changing approaches over years. The concept of fully developed speckle is explained. Drawbacks of homomorphic filtering are pointed out. Assets of multiresolution despeckling, as opposite to spatial-domain despeckling, are highlighted. Also advantages of undecimated, or stationary, wavelet transforms over decimated ones are discussed. Bayesian estimators and probability density function (pdf) models in both spatial and multiresolution domains are reviewed. Scale-space varying pdf models, as opposite to scale varying models, are promoted. Promising methods following non-Bayesian approaches, like nonlocal (NL) filtering and total variation (TV) regularization, are reviewed and compared to spatial- and wavelet-domain Bayesian filters. Both established and new trends for assessment of despeckling are presented. A few experiments on simulated data and real COSMO-SkyMed SAR images highlight, on one side the costperformance tradeoff of the different methods, on the other side the effectiveness of solutions purposely designed for SAR heterogeneity and not fully developed speckle. Eventually, upcoming methods based on new concepts of signal processing, like compressive sensing, are foreseen as a new generation of despeckling, after spatial-domain and multiresolution-domain method

Amplitude vs intensity Bayesian despeckling in the wavelet domain for SAR images

In this paper, the problem of despeckling SAR images when the input data is either an intensity or an amplitude signal is revisited. State-of-the-art despeckling methods based on Bayesian estimators in the wavelet domain, recently proposed in the literature, are taken into consideration. First, how these methods proposed for one format (e.g., intensity) can be adapted to the other format (e.g., amplitude) is investigated. Second, the performance of such algorithms in both cases is analyzed. Experimental results carried out on simulated speckled images and on true SAR data are presented and discussed in order to assess the best strategy. From these results, it can be observed that filtering in the amplitude domain yields better performances in terms of objective quality indexes, such as preservation of structural details, as well as in terms of visual inspection of the filtered SAR dat

Blind Speckle Decorrelation for SAR Image Despeckling

In the past few decades, several methods have been developed for despeckling synthetic aperture radar (SAR) images. A considerable number of them have been derived under the assumption of a fully-developed speckle model in which the multiplicative speckle noise is supposed to be a white process. Unfortunately, the transfer function of SAR acquisition systems can introduce a statistical correlation, which decreases the despeckling efficiency of such filters. In this paper, a whitening method is proposed for processing a complex image acquired by a SAR system. We demonstrate that the proposed approach allows the successful application of classical despeckling algorithms. First, we perform an estimation of the SAR system frequency response based on some statistical properties of the acquired image and by using realistic assumptions. Then, a decorrelation process is applied on the acquired image, taking into account the presence of point targets. Finally, the image is despeckled. The experimental results show that the despeckling filters achieve better performance when they are preceded by the proposed whitening method; furthermore, the radiometric characteristics of the image are preserve

Dynamics of Vulcano Island investigated by long-term (40 years) geophysical data

Vulcano island is a composite volcanic edifice located in the south-central sector of the Aeolian Archipelago (Tyrrhenian Sea, Italy) and it is an important tourist destination. Historic activity has been characterized by frequent transitions from phereatomagmatic to minor magmatic activity. The last eruption in 1888-90 was characterized by energetic explosive pulses and defines the “vulcanian” type of activity. Since then, volcanic activity has taken the form of fumarolic emanations of variable intensity and temperature, mainly concentrated at “La Fossa” crater, with maximum temperatures ranging between 200° and 300° C; temperature increases and changes in the gas chemistry, were often observed. The most recent episode began in the 80’s when fumarole temperature progressively increased to 690°C in May 1993. Vulcano is active and this favoured monitoring and research studies, in particular focussed on the most recent structures. In the frame of DPC-INGV “V3” project, we investigate the Vulcano dynamics through ca. 40 years of ground deformation and seismicity data collected by the discrete and continuous INGV monitoring networks. We considered levelling, GPS, EDM, seismic and tilt data. EDM and levelling measurements began in the middle 1970s and since the late 1990s the EDM benchmarks have been measurered using GPS. We observed three scales of ground deformation: the first one seems to be linked to the regional tectonics, with a general transpressive kinematics; the second one affects the northern half of the island and could be related to the caldera dynamics; the third one affects only the cone of La Fossa. Regional tectonic stress seems to play an important role in the transition of the volcanic system from a phase of stability to a phase of unrest, inducing the heating and the expansion of shallow hydrothermal fluids. Ground deformation at Vulcano may be linked to the geothermal system rather than magmatic sources

Automated Heuristic Optimization of Prostate VMAT Treatment Planning

Purpose: To investigate a genetic algorithm approach to automatic treatment planning. Methods: A Python script based on genetic algorithm (GA) was implemented for VMAT treatment planning of prostate tumor. The script was implemented in RayStation treatment planning system using Python code. Two different clinical prescriptions were considered: 78 Gy prescribed to planning target volume in 39 fractions (GROUP 1) and simultaneous integrated boost (70.2 Gy to prostate bed and 61.1 Gy to seminal vesicles) in 26 fractions (GROUP 2). The script automatically optimizes doses to PTV and OARs according to GA. A comparison with corresponding plans created with Monaco TPS (M) and Auto-Planning module of Pinnacle3 (AP) was carried out. The plans were evaluated with a total score (TS) of PlanIQ software in terms of target coverage and sparing of OARs as well as clinical score (CS) performed by a Radiation Oncologist. Results: In GROUP 1, mean value of TS were 150.6 ± 30.7, 146.3 ± 36.1 and 137.4 ± 35.7 for AP, GA and M respectively. For GROUP 2, mean value for TS were 163.5 ± 16.8, 163.4 ± 24.7 and 162.9 ± 16.6 for AP, GA and M respectively with no significance differences. In terms of CS, the highest value has been attributed to GA in four patients out of five for both GROUP 1 and 2. Conclusions: Genetic approach is practicable for prostate VMAT plan generation and studies are underway in other anatomical sites such as Head and Neck and Rectum

Geophysical monitoring of Stromboli volcano: insight into recent volcanic activity

Stromboli is an open conduit strato-volcano of the Aeolian archipelago (Italy), characterized by typical Strom-bolian explosive activity, lasting for several centuries, and by the emission of huge amounts of gas. The normalactivity of Stromboli is characterized by some hundreds of moderate explosions per day. Major explosions, whichlaunch scoria up to hundreds of meters from the craters, lava flows and paroxysmal explosions, which producelarge ballistic blocks, sometimes take place. During the effusive eruption in 2002 - 2003, which caused a tsunamiwith waves of about 10 meters high along the coasts of the Island, the monitoring system was enhanced. In 2006INGV has added two Sacks-Evertson borehole volumetric dilatometers to the surveillance system, in order to mon-itor changes in the local strain field by measuring areal strain. Today we have a large amount of geophysical dataand observations that allow us to better understand how this volcano works. After a period of low explosive activitystarted in mid-2014, Stromboli has shown a more intense explosive activity in the last few months. During the re-cent phase of increased activity, the geophysical monitoring system detected four major explosions occurred on 26July, 23 October, 1 November and 1 December 2017, respectively. The current phase of reawakening of Strombolivolcano has led the Italian civil protection authorities to decree the "attention" alert level (yellow) on the Island.PublishedVienna, Austria1IT. Reti di monitoraggio e sorveglianz

MT. ETNA FEBRUARY-APRIL 2013 ERUPTION

European FP7 MED-SUV (MEditerranean SUpersite Volcanoes). Grant Number: 308665 European Research Council European FP7 (FP/2007-2013)/ERC. Grant Number: 279802 SIGMA (Sistema Integrato di sensori in ambiente cloud per la Gestione Multirischio Avanzata

Filtering of Interferometric SAR Phase Images as a Fuzzy Matching-Pursuit Blind Estimation

We present an original application of fuzzy logic to restoration of phase images from interferometric synthetic aperture radar (InSAR), which are affected by zero-mean uncorrelated noise, whose variance depends on the underlying coherence, thereby yielding a nonstationary random noise process. Spatial filtering of the phase noise is recommended, either before phase unwrapping is accomplished, or simultaneously with it. In fact, phase unwrapping basically relies on a smoothness constraint of the phase field, which is severely hampered by the noise. Space-varying linear MMSE estimation is stated as a problem of matching pursuit, in which the estimator is obtained as an expansion in series of a finite number of prototype estimators, fitting the spatial features of the different statistical classes encountered, for example, fringes and steep slope areas. Such estimators are calculated in a fuzzy fashion through an automatic training procedure. The space-varying coefficients of the expansion are stated as degrees of fuzzy membership of a pixel to each of the estimators. Neither a priori knowledge on the noise variance is required nor particular signal and noise models are assumed. Filtering performances on simulated phase images show a steady SNR improvement over conventional box filtering. Applications of the proposed filter to interferometric phase images demonstrate a superior ability of restoring fringes yet preserving their discontinuities, together with an effective noise smoothing performance, irrespective of locally varying coherence characteristics

Open-path FTIR spectroscopy of magma degassing processes during eight lava fountains on Mount Etna

International audienceIn June-July 2001 a series of 16 discrete lava fountain paroxysms occurred at the Southeast summit crater (SEC) of Mount Etna, preceding a 28-day long violent flank eruption. Each paroxysm was preceded by lava effusion, growing seismic tremor and a crescendo of Strombolian explosive activity culminating into powerful lava fountaining up to 500 m in height. During 8 of these 16 events we could measure the chemical composition of the magmatic gas phase (H2O, CO2, SO2, HCl, HF and CO), using open-path Fourier transform infrared (OP-FTIR) spectrometry at ∼1-2 km distance from SEC and absorption spectra of the radiation emitted by hot lava fragments. We show that each fountaining episode was characterized by increasingly CO2-rich gas release, with CO2/SO2 and CO2/HCl ratios peaking in coincidence with maxima in seismic tremor and fountain height, whilst the SO2/HCl ratio showed a weak inverse relationship with respect to eruption intensity. Moreover, peak values in both CO2/SO2 ratio and seismic tremor amplitude for each paroxysm were found to increase linearly in proportion with the repose interval (2-6 days) between lava fountains. These observations, together with a model of volatile degassing at Etna, support the following driving process. Prior to and during the June-July 2001 lava fountain sequence, the shallow (∼2 km) magma reservoir feeding SEC received an increasing influx of deeply derived carbon dioxide, likely promoted by the deep ascent of volatile-rich primitive basalt that produced the subsequent flank eruption. This CO2-rich gas supply led to gas accumulation and overpressure in SEC reservoir, generating a bubble foam layer whose periodical collapse powered the successive fountaining events. The anti-correlation between SO2/HCl and eruption intensity is best explained by enhanced syn-eruptive degassing of chlorine from finer particles produced during more intense magma fragmentation