Validation of surface velocity estimated from satellite images

This report concerns the validation of surface velocity estimated from satellite images. The estimation is obtained with a dynamic model based on shallow-water equations. We first compare the stationary assumption to the shallow-water heuristics to justify our choice. Second, we quantify the quality of the estimation by measuring the misfit between the model output and the altimetry measures. Experiments are achieved on Sea Surface Temperature data acquired by the NOAA/AVHRR satellites over the Black Sea. The altimetry measures are obtained by two radar sensors: Envisat and GFO. The good adequacy between the shallow-water output and the altimetry data validates our motion estimation approach.Ce rapport de recherche concerne la validation de l'estimation de la vitesse de surface à partir d'images satellite. Cette estimation est effectuée avec un modèle de la dynamique, basé sur les équations shallow-water. Nous comparons d'abord l'hypothèse de stationnarité aux équations shallow-water afin de justifier notre choix. Puis, nous quantifions la qualité des estimations en mesurant l'écart entre la sortie du modèle et les mesures d'altimétrie. Les expérimentations sont effectuées en utilisant des données de température de surface, acquises au-dessus de la Mer Noire avec les satellites NOAA/AVHRR. Les mesures altimétriques proviennent de deux capteurs radar : Envisat et GFO. La bonne adéquation entre la sortie du modèle shallow-water et les données altimétriques valide notre approche d'estimation du mouvement

Vedel-objektiiv abil salvestatud kaugseire piltide analüüs kasutades super-resolutsiooni meetodeid

Väitekirja elektrooniline versioon ei sisalda publikatsiooneKäesolevas doktoritöös uuriti nii riist- kui ka tarkvaralisi lahendusi piltide töötlemiseks. Riist¬varalise poole pealt pakuti lahenduseks uudset vedelläätse, milles on dielekt¬rilisest elastomeerist kihilise täituriga membraan otse optilisel teljel. Doktoritöö käigus arendati välja kaks prototüüpi kahe erineva dielektrilisest elastomeerist ki¬hilise täituriga, mille aktiivne ala oli ühel juhul 40 ja teisel 20 mm. Läätse töö vas¬tas elastomeeri deformatsiooni mehaanikale ja suhtelistele muutustele fookuskau¬guses. Muutuste demonstreerimiseks meniskis ja läätse fookuskauguse mõõtmiseks kasutati laserkiirt. Katseandmetest selgub, et muutuste tekitamiseks on vajalik pinge vahemikus 50 kuni 750 volti. Tarkvaralise poole pealt pakuti uut satelliitpiltide parandamise süsteemi. Paku¬tud süsteem jagas mürase sisendpildi DT-CWT laineteisenduse abil mitmeteks sagedusalamribadeks. Pärast müra eemaldamist LA-BSF funktsiooni abil suu¬rendati pildi resolutsiooni DWT-ga ja kõrgsagedusliku alamriba piltide interpo¬leerimisega. Interpoleerimise faktor algsele pildile oli pool sellest, mida kasutati kõrgsagedusliku alamriba piltide interpoleerimisel ning superresolutsiooniga pilt rekonst¬rueeriti IDWT abil. Käesolevas doktoritöös pakuti tarkvaraliseks lahenduseks uudset sõnastiku baasil töötavat super-resolutsiooni (SR) meetodit, milles luuakse paarid suure resolutsiooniga (HR) ja madala resolut-siooniga (LR) piltidest. Kõigepealt jagati vastava sõnastiku loomiseks HR ja LR paarid omakorda osadeks. Esialgse HR kujutise saamiseks LR sisendpildist kombineeriti HR osi. HR osad valiti sõnastikust nii, et neile vastavad LR osad oleksid võimalikult lähedased sisendiks olevale LR pil¬dile. Iga valitud HR osa heledust korrigeeriti, et vähendada kõrvuti asuvate osade heleduse erine¬vusi superresolutsiooniga pildil. Plokkide efekti vähendamiseks ar¬vutati saadud SR pildi keskmine ning bikuupinterpolatsiooni pilt. Lisaks pakuti käesolevas doktoritöös välja kernelid, mille tulemusel on võimalik saadud SR pilte teravamaks muuta. Pakutud kernelite tõhususe tõestamiseks kasutati [83] ja [50] poolt pakutud resolutsiooni parandamise meetodeid. Superreso¬lutsiooniga pilt saadi iga kerneli tehtud HR pildi kombineerimise teel alpha blen¬dingu meetodit kasutades. Pakutud meetodeid ja kerneleid võrreldi erinevate tavaliste ja kaasaegsete meetoditega. Kvantita-tiivsetest katseandmetest ja saadud piltide kvaliteedi visuaal¬sest hindamisest selgus, et pakutud meetodid on tavaliste kaasaegsete meetoditega võrreldes paremad.In this thesis, a study of both hardware and software solutions for image enhance¬ment has been done. On the hardware side, a new liquid lens design with a DESA membrane located directly in the optical path has been demonstrated. Two pro¬totypes with two different DESA, which have a 40 and 20 mm active area in diameter, were developed. The lens performance was consistent with the mechan¬ics of elastomer deformation and relative focal length changes. A laser beam was used to show the change in the meniscus and to measure the focal length of the lens. The experimental results demonstrate that voltage in the range of 50 to 750 V is required to create change in the meniscus. On the software side, a new satellite image enhancement system was proposed. The proposed technique decomposed the noisy input image into various frequency subbands by using DT-CWT. After removing the noise by applying the LA-BSF technique, its resolution was enhanced by employing DWT and interpolating the high-frequency subband images. An original image was interpolated with half of the interpolation factor used for interpolating the high-frequency subband images, and the super-resolved image was reconstructed by using IDWT. A novel single-image SR method based on a generating dictionary from pairs of HR and their corresponding LR images was proposed. Firstly, HR and LR pairs were divided into patches in order to make HR and LR dictionaries respectively. The initial HR representation of an input LR image was calculated by combining the HR patches. These HR patches are chosen from the HR dictionary corre-sponding to the LR patches that have the closest distance to the patches of the in¬put LR image. Each selected HR patch was processed further by passing through an illumination enhancement processing order to reduce the noticeable change of illumination between neighbor patches in the super-resolved image. In order to reduce the blocking effect, the average of the obtained SR image and the bicubic interpolated image was calculated. The new kernels for sampling have also been proposed. The kernels can improve the SR by resulting in a sharper image. In order to demonstrate the effectiveness of the proposed kernels, the techniques from [83] and [50] for resolution enhance¬ment were adopted. The super-resolved image was achieved by combining the HR images produced by each of the proposed kernels using the alpha blending tech-nique. The proposed techniques and kernels are compared with various conventional and state-of-the-art techniques, and the quantitative test results and visual results on the final image quality show the superiority of the proposed techniques and ker¬nels over conventional and state-of-art technique

Divergence-Free Motion Estimation

International audienceThis paper describes an innovative approach to estimate motion from image observations of divergence-free flows. Unlike most state-of-the-art methods, which only minimize the divergence of the motion field, our approach utilizes the vorticity-velocity formalism in order to construct a motion field in the subspace of divergence free functions. A 4DVAR-like image assimilation method is used to generate an estimate of the vorticity field given image observations. Given that vorticity estimate, the motion is obtained solving the Poisson equation. Results are illustrated on synthetic image observations and compared to those obtained with state-of-the-art methods, in order to quantify the improvements brought by the presented approach. The method is then applied to ocean satellite data to demonstrate its performance on the real images

Self-similar prior and wavelet bases for hidden incompressible turbulent motion

This work is concerned with the ill-posed inverse problem of estimating turbulent flows from the observation of an image sequence. From a Bayesian perspective, a divergence-free isotropic fractional Brownian motion (fBm) is chosen as a prior model for instantaneous turbulent velocity fields. This self-similar prior characterizes accurately second-order statistics of velocity fields in incompressible isotropic turbulence. Nevertheless, the associated maximum a posteriori involves a fractional Laplacian operator which is delicate to implement in practice. To deal with this issue, we propose to decompose the divergent-free fBm on well-chosen wavelet bases. As a first alternative, we propose to design wavelets as whitening filters. We show that these filters are fractional Laplacian wavelets composed with the Leray projector. As a second alternative, we use a divergence-free wavelet basis, which takes implicitly into account the incompressibility constraint arising from physics. Although the latter decomposition involves correlated wavelet coefficients, we are able to handle this dependence in practice. Based on these two wavelet decompositions, we finally provide effective and efficient algorithms to approach the maximum a posteriori. An intensive numerical evaluation proves the relevance of the proposed wavelet-based self-similar priors.Comment: SIAM Journal on Imaging Sciences, 201

Variational fluid flow measurements from image sequences: synopsis and perspectives

[Departement_IRSTEA]Ecotechnologies [TR1_IRSTEA]SPEEVariational approaches to image motion segmentation has been an active field of study in image processing and computer vision for two decades. We present a short overview over basic estimation schemes and report in more detail recent modifications and applications to fluid flow estimation. Key properties of these approaches are illustrated by numerical examples. We outline promising research directions and point out the potential of variational techniques in combination with correlation-based PIV methods, for improving the consistency of fluid flow estimation and simulation

Automatic Decomposition of Geodetic Time Series for Studies of Ground Deformation

Geodetic measurements of surface deformation have been used for several decades to study how the Earth's surface responds to a wide range of geophysical processes. Geodetic time series acquired over a finite spatial extent can be used to quantify the time dependence of surface strain for a wide range of spatial and temporal scales. In this thesis, we present a new method for automatically decomposing geodetic time series into temporal components corresponding to different geophysical processes. This method relies on constructing an overcomplete temporal dictionary of reference functions such that any geodetic signal can be described by a linear combination of the functions in the dictionary. By solving a linear least squares problem with sparsity-inducing regularization, we can limit the total number of dictionary elements needed to reconstruct a signal. In Chapter 2, we present the development of this method in the context of transient detection, where we define transient deformation as nonperiodic, nonsecular accumulation of strain in the crust. The sparsity regularization term automatically localizes the dominant timescales and onset times of any transient signals. We apply this method to Global Positioning System (GPS) data for a slow slip event in the Cascadia subduction zone while incorporating a spatial weighting scheme that filters for spatially coherent signals. In Chapter 3, we use a combination of unique space geodetic measurements and seismic observations to study the 2014 collapse of Bárðarbunga Caldera in Iceland associated with a major eruption event. The eruption sequence, which involved deflation of a magma chamber underneath the caldera and emplacement of a dike leading to lava flow, resulted in rapid subsidence of the glacier surface overlying the caldera and wide-scale ground deformation encompassing the rift zone associated with the dike emplacement. We present a model of the collapse that suggests that the majority of the observed subsidence occurs aseismically via a deflating sill-like magma chamber. In Chapter 4, we extend upon the transient detection framework presented in Chapter 2 to study complex surface deformation over groundwater basins near Los Angeles, California. We develop a distributed time series analysis framework based on the sparse estimation techniques of Chapter 2 and apply it to an 18-year interferometric synthetic aperture radar (InSAR) time series covering the Los Angeles area. We compare long- and short-term ground deformation signals to hydraulic head data from monitoring wells to understand the mechanical link between pressure variations in subsurface aquifers and observed ground deformation

Assessment and application of wavelet-based optical flow velocimetry (wOFV) to wall-bounded turbulent flows

The performance of a wavelet-based optical flow velocimetry (wOFV) algorithm to extract high accuracy and high resolution velocity fields from particle images in wall-bounded turbulent flows is assessed. wOFV is first evaluated using synthetic particle images generated from a channel flow DNS of a turbulent boundary layer. The sensitivity of wOFV to the regularization parameter (lambda) is quantified and results are compared to PIV. Results on synthetic particle images indicated different sensitivity to under-regularization or over-regularization depending on which region of the boundary layer is analyzed. Synthetic data revealed that wOFV can modestly outperform PIV in vector accuracy across a broad lambda range. wOFV showed clear advantages over PIV in resolving the viscous sublayer and obtaining highly accurate estimates of the wall shear stress. wOFV was also applied to experimental data of a developing turbulent boundary layer. Overall, wOFV revealed good agreement with both PIV and PIV + PTV. However, wOFV was able to successfully resolve the wall shear stress and correctly normalize the boundary layer streamwise velocity to wall units where PIV and PIV + PTV showed larger deviations. Analysis of the turbulent velocity fluctuations revealed spurious results for PIV in close proximity to the wall, leading to significantly exaggerated and non-physical turbulence intensity. PIV + PTV showed a minor improvement in this aspect. wOFV did not exhibit this same effect, revealing that it is more accurate in capturing small-scale turbulent motion in the vicinity of boundaries. The enhanced vector resolution of wOFV enabled improved estimation of instantaneous derivative quantities and intricate flow structure both closer to the wall. These aspects show that, within a reasonable lambda range, wOFV can improve resolving the turbulent motion occurring in the vicinity of physical boundaries

Cassini observations reveal a regime of zonostrophic macroturbulence on Jupiter

In December 2000, the Cassini fly-by near Jupiter delivered high-resolution images of Jupiter’s clouds over the entire planet in a band between 50°N and 50°S. Three daily-averaged two-dimensional velocity snapshots extracted from these images are used to perform spectral analysis of jovian atmospheric macroturbulence. A similar analysis is also performed on alternative data documented by Choi and Showman (Choi, D., Showman, A. [2011]. Icarus 216, 597–609), based on a different method of image processing. The inter-comparison of the products of both analyses ensures a better constraint of the spectral estimates. Both analyses reveal strong anisotropy of the kinetic energy spectrum. The zonal spectrum is very steep and most of the kinetic energy resides in slowly evolving, alternating zonal (west–east) jets, while the non-zonal, or residual spectrum obeys the Kolmogorov–Kraichnan law specific to two-dimensional turbulence in the range of the inverse energy cascade. The spectral data is used to estimate the inverse cascade rate ∊ and the zonostrophy index Rβ for the first time. Although both datasets yield somewhat different values of ∊, it is estimated to be in the range 0.5–1.0 × 10−5 m2 s−3. The ensuing values of Rβ ≳ 5 belong well in the range of zonostrophic turbulence whose threshold corresponds to Rβ ≃ 2.5. We infer that the large-scale circulation is maintained by an anisotropic inverse energy cascade. The removal of the Great Red Spot from both datasets has no significant effect upon either the spectra or the inverse cascade rate. The spectral data are used to compute the rate of the energy exchange, W, between the non-zonal structures and the large-scale zonal flow. It is found that instantaneous values of W may exceed ∊ by an order of magnitude. Previous numerical simulations with a barotropic model suggest that W and ∊ attain comparable values only after averaging of W over a sufficiently long time. Near-instantaneous values of W that have been routinely used to infer the rate of the kinetic energy supply to Jupiter’s zonal flow may therefore significantly overestimate ∊. This disparity between W and ∊ may resolve the long-standing conundrum of an unrealistically high rate of energy transfer to the zonal flow. The meridional diffusivity Kϕ in the regime of zonostrophic turbulence is given by an expression that depends on ∊. The value of Kϕ estimated from the spectra is compared against data from the dispersion of stratospheric gases and debris resulting from the Shoemaker-Levy 9 comet and Wesley asteroid impacts in 1994 and 2009 respectively. Not only is Kϕ found to be consistent with estimates for both impacts, but the eddy diffusivity found from observations appears to be scale-independent. This behaviour could be a consequence of the interaction between anisotropic turbulence and Rossby waves specific to the regime of zonostrophic macroturbulence

The development of virtual leaf surface models for interactive agrichemical spray applications

This project constructed virtual plant leaf surfaces from digitised data sets for use in droplet spray models. Digitisation techniques for obtaining data sets for cotton, chenopodium and wheat leaves are discussed and novel algorithms for the reconstruction of the leaves from these three plant species are developed. The reconstructed leaf surfaces are included into agricultural droplet spray models to investigate the effect of the nozzle and spray formulation combination on the proportion of spray retained by the plant. A numerical study of the post-impaction motion of large droplets that have formed on the leaf surface is also considered