Source author record

Nicolas Desassis

Nicolas Desassis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2023arXiv

Geostatistics for large datasets on Riemannian manifolds: a matrix-free approach

Large or very large spatial (and spatio-temporal) datasets have become common place in many environmental and climate studies. These data are often collected in non-Euclidean spaces (such as the planet Earth) and they often present non-stationary anisotropies. This paper proposes a generic approach to model Gaussian Random Fields (GRFs) on compact Riemannian manifolds that bridges the gap between existing works on non-stationary GRFs and random fields on manifolds. This approach can be applied to any smooth compact manifolds, and in particular to any compact surface. By defining a Riemannian metric that accounts for the preferential directions of correlation, our approach yields an interpretation of the ''local anisotropies'' as resulting from ''local'' deformations of the domain. We provide scalable algorithms for the estimation of the parameters and for optimal prediction by kriging and simulation able to tackle very large grids. Stationary and non-stationary illustrations are provided.

preprint2020arXiv

A matrix-free approach to geostatistical filtering

In this paper, we present a novel approach to geostatistical filtering which tackles two challenges encountered when applying this method to complex spatial datasets: modeling the non-stationarity of the data while still being able to work with large datasets. The approach is based on a finite element approximation of Gaussian random fields expressed as an expansion of the eigenfunctions of a Laplace--Beltrami operator defined to account for local anisotropies. The numerical approximation of the resulting random fields using a finite element approach is then leveraged to solve the scalability issue through a matrix-free approach. Finally, two cases of application of this approach, on simulated and real seismic data are presented.

preprint2018arXiv

Efficient simulation of Gaussian Markov random fields by Chebyshev polynomial approximation

This paper presents an algorithm to simulate Gaussian random vectors whose precision matrix can be expressed as a polynomial of a sparse matrix. This situation arises in particular when simulating Gaussian Markov random fields obtained by the discretization by finite elements of the solutions of some stochastic partial derivative equations. The proposed algorithm uses a Chebyshev polynomial approximation to compute simulated vectors with a linear complexity. This method is asymptotically exact as the approximation order grows. Criteria based on tests of the statistical properties of the produced vectors are derived to determine minimal orders of approximation.

preprint2015arXiv

A pairwise likelihood approach for the empirical estimation of the underlyingvariograms in the plurigaussian models

The plurigaussian model is particularly suited to describe categorical regionalized variables. Starting from a simple principle, the thresh-olding of one or several Gaussian random fields (GRFs) to obtain categories, the plurigaussian model is well adapted for a wide range ofsituations. By acting on the form of the thresholding rule and/or the threshold values (which can vary along space) and the variograms ofthe underlying GRFs, one can generate many spatial configurations for the categorical variables. One difficulty is to choose variogrammodel for the underlying GRFs. Indeed, these latter are hidden by the truncation and we only observe the simple and cross-variogramsof the category indicators. In this paper, we propose a semiparametric method based on the pairwise likelihood to estimate the empiricalvariogram of the GRFs. It provides an exploratory tool in order to choose a suitable model for each GRF and later to estimate its param-eters. We illustrate the efficiency of the method with a Monte-Carlo simulation study .The method presented in this paper is implemented in the R packageRGeostats.

preprint2014arXiv

A Generalized Convolution Model and Estimation for Non-stationary Random Functions

Standard geostatistical models assume second order stationarity of the underlying Random Function. In some instances, there is little reason to expect the spatial dependence structure to be stationary over the whole region of interest. In this paper, we introduce a new model for second order non-stationary Random Functions as a convolution of an orthogonal random measure with a spatially varying random weighting function. This new model is a generalization of the common convolution model where a non-random weighting function is used. The resulting class of non-stationary covariance functions is very general, flexible and allows to retrieve classes of closed-form non-stationary covariance functions known from the literature, for a suitable choices of the random weighting functions family. Under the framework of a single realization and local stationarity, we develop parameter inference procedure of these explicit classes of non-stationary covariance functions. From a local variogram non-parametric kernel estimator, a weighted local least-squares approach in combination with kernel smoothing method is developed to estimate the parameters. Performances are assessed on two real datasets: soil and rainfall data. It is shown in particular that the proposed approach outperforms the stationary one, according to several criteria. Beyond the spatial predictions, we also show how conditional simulations can be carried out in this non-stationary framework.

preprint2014arXiv

Estimation of Space Deformation Model for Non-stationary Random Functions

Stationary Random Functions have been successfully applied in geostatistical applications for decades. In some instances, the assumption of a homogeneous spatial dependence structure across the entire domain of interest is unrealistic. A practical approach for modelling and estimating non-stationary spatial dependence structure is considered. This consists in transforming a non-stationary Random Function into a stationary and isotropic one via a bijective continuous deformation of the index space. So far, this approach has been successfully applied in the context of data from several independent realizations of a Random Function. In this work, we propose an approach for non-stationary geostatistical modelling using space deformation in the context of a single realization with possibly irregularly spaced data. The estimation method is based on a non-stationary variogram kernel estimator which serves as a dissimilarity measure between two locations in the geographical space. The proposed procedure combines aspects of kernel smoothing, weighted non-metric multi-dimensional scaling and thin-plate spline radial basis functions. On a simulated data, the method is able to retrieve the true deformation. Performances are assessed on both synthetic and real datasets. It is shown in particular that our approach outperforms the stationary approach. Beyond the prediction, the proposed method can also serve as a tool for exploratory analysis of the non-stationarity.