Xavier Descombes

A Gibbs Point Process for Road Extraction from Remotely Sensed Images

In this paper we propose a new method for the extraction of roads from remotely sensed images. Under the assumption that roads form a thin network in the image, we approximate such a network by connected line segments.To perform this task, we construct a point process able to simulate and detect thin networks. The segments have to be connected, in order to form a line-network. Aligned segments are favored whereas superposition is penalized. These constraints are enforced by the interaction model (called the Candy model). The specific properties of the road network in the image are described by the data term. This term is based on statistical hypothesis tests.The proposed probabilistic model can be written within a Gibbs point process framework. The estimate for the network is found by minimizing an energy function. In order to avoid local minima, we use a simulated annealing algorithm, based on a Monte Carlo dynamics (RJMCMC) for finite point processes. Results are shown on SPOT, ERS and aerial images.

Structural Approach for Building Reconstruction from a Single DSM

We present a new approach for building reconstruction from a single Digital Surface Model (DSM). It treats buildings as an assemblage of simple urban structures extracted from a library of 3D parametric blocks (like a LEGO set). First, the 2D-supports of the urban structures are extracted either interactively or automatically. Then, 3D-blocks are placed on the 2D-supports using a Gibbs model which controls both the block assemblage and the fitting to data. A Bayesian decision finds the optimal configuration of 3D--blocks using a Markov Chain Monte Carlo sampler associated with original proposition kernels. This method has been validated on multiple data set in a wide-resolution interval such as 0.7 m satellite and 0.1 m aerial DSMs, and provides 3D representations on complex buildings and dense urban areas with various levels of detail.

Spatio-temporal fMRI analysis using Markov random fields

Functional magnetic resonance images (fMRIs) provide high-resolution datasets which allow researchers to obtain accurate delineation and sensitive detection of activation areas involved in cognitive processes. To preserve the resolution of this noninvasive technique, refined methods are required in the analysis of the data. In this paper, the authors first discuss the widely used methods based on a statistical parameter map (SPM) analysis exposing the different shortcomings of this approach when considering high-resolution data. First, the often used Gaussian filtering results in a blurring effect and in delocalization of the activated area. Secondly, the SPM approach only considers false alarms due to noise but not rejections of activated voxels. The authors propose to embed the fMRI analysis problem into a Bayesian framework consisting of two steps: (i) data restoration and (ii) data analysis. They, therefore, propose two Markov random fields (MRFs) to solve these two problems. Results on three protocols (visual, motor and word recognition) are shown for two SPM approaches and compared with the proposed MRF-approach.

Estimation of Markov random field prior parameters using Markov chain Monte Carlo maximum likelihood

Recent developments in statistics now allow maximum likelihood estimators for the parameters of Markov random fields (MRFs) to be constructed. We detail the theory required, and present an algorithm that is easily implemented and practical in terms of computation time. We demonstrate this algorithm on three MRF models--the standard Potts model, an inhomogeneous variation of the Potts model, and a long-range interaction model, better adapted to modeling real-world images. We estimate the parameters from a synthetic and a real image, and then resynthesize the models to demonstrate which features of the image have been captured by the model. Segmentations are computed based on the estimated parameters and conclusions drawn.

Building Outline Extraction from Digital Elevation Models Using Marked Point Processes

This work presents an automatic algorithm for extracting vectorial land registers from altimetric data in dense urban areas. We focus on elementary shape extraction and propose a method that extracts rectangular buildings. The result is a vectorial land register that can be used, for instance, to perform precise roof shape estimation. Using a spatial point process framework, we model towns as configurations of and unknown number of rectangles. An energy is defined, which takes into account both low level information provided by the altimetry of the scene, and geometric knowledge about the disposition of buildings in towns. Estimation is done by minimizing the energy using simulated annealing. We use an MCMC sampler that is a combination of general Metropolis Hastings Green techniques and the Geyer and Møller algorithm for point process sampling. We define some original proposition kernels, such as birth or death in a neighborhood and define the energy with respect to an inhomogeneous Poisson point process. We present results on real data provided by the IGN (French National Geographic Institute). Results were obtained automatically. These results consist of configurations of rectangles describing a dense urban area.

Texture feature analysis using a gauss-Markov model in hyperspectral image classification

Texture analysis has been widely investigated in the monospectral and multispectral imagery domains. At the same time, new image sensors with a large number of bands (more than ten) have been designed. They are able to provide images with both fine spectral and spatial sampling, and are called hyperspectral images. The aim of this work is to perform a joint texture analysis in both discrete spaces. To achieve this goal, we propose a probabilistic vector texture model, using a Gauss-Markov random field (MRF). The MRF parameters allow the characterization of different hyperspectral textures. A possible application of this work is the classification of urban areas. These areas are not well characterized by radiometry alone, and so we use the MRF parameters as new features in a maximum-likelihood classification algorithm. The results obtained on Airborne Visible/Infrared Imaging Spectrometer hyperspectral images demonstrate that a better classification is achieved when texture information is included in the analysis.

Object Extraction Using a Stochastic Birth-and-Death Dynamics in Continuum

We define a new birth and death dynamics dealing with configurations of disks in the plane. We prove the convergence of the continuous process and propose a discrete scheme converging to the continuous case. This framework is developed to address image processing problems consisting in detecting a configuration of objects from a digital image. The derived algorithm is applied for tree crown extraction and bird detection from aerial images. The performance of this approach is shown on real data.

Building Development Monitoring in Multitemporal Remotely Sensed Image Pairs with Stochastic Birth-Death Dynamics

In this paper, we introduce a new probabilistic method which integrates building extraction with change detection in remotely sensed image pairs. A global optimization process attempts to find the optimal configuration of buildings, considering the observed data, prior knowledge, and interactions between the neighboring building parts. We present methodological contributions in three key issues: 1) We implement a novel object-change modeling approach based on Multitemporal Marked Point Processes, which simultaneously exploits low-level change information between the time layers and object-level building description to recognize and separate changed and unaltered buildings. 2) To answer the challenges of data heterogeneity in aerial and satellite image repositories, we construct a flexible hierarchical framework which can create various building appearance models from different elementary feature-based modules. 3) To simultaneously ensure the convergence, optimality, and computation complexity constraints raised by the increased data quantity, we adopt the quick Multiple Birth and Death optimization technique for change detection purposes, and propose a novel nonuniform stochastic object birth process which generates relevant objects with higher probability based on low-level image features.

Geometric Feature Extraction by a Multimarked Point Process

This paper presents a new stochastic marked point process for describing images in terms of a finite library of geometric objects. Image analysis based on conventional marked point processes has already produced convincing results but at the expense of parameter tuning, computing time, and model specificity. Our more general multimarked point process has simpler parametric setting, yields notably shorter computing times, and can be applied to a variety of applications. Both linear and areal primitives extracted from a library of geometric objects are matched to a given image using a probabilistic Gibbs model, and a Jump-Diffusion process is performed to search for the optimal object configuration. Experiments with remotely sensed images and natural textures show that the proposed approach has good potential. We conclude with a discussion about the insertion of more complex object interactions in the model by studying the compromise between model complexity and efficiency.

A Marked Point Process of Rectangles and Segments for Automatic Analysis of Digital Elevation Models

This work presents a framework for automatic feature extraction from images using stochastic geometry. Features in images are modeled as realizations of a spatial point process of geometrical shapes. This framework allows the incorporation of a priori knowledge on the spatial repartition of features. More specifically, we present a model based on the superposition of a process of segments and a process of rectangles. The former is dedicated to the detection of linear networks of discontinuities, whereas the latter aims at segmenting homogeneous areas. An energy is defined, favoring connections of segments, alignments of rectangles, and a relevant interaction between both types of objects. The estimation is performed by minimizing the energy using a simulated annealing algorithm. The proposed model is applied to the analysis of digital elevation models (DEMs). These images are raster data representing the altimetry of a dense urban area. We present results on real data provided by the French National Geographic Institute (IGN) consisting in low-quality DEMs of various types.