Parallelizing remote sensing image geometric correction

Remote sensing spatial, spectral, and temporal resolutions of images, acquired over a reasonably sized image extent, result in imagery that can be processed to represent land cover over large areas with an amount of spatial detail that is very attractive for monitoring, management, and scienti c activities. With Moore's Law alive and well, more and more parallelism is introduced into all computing platforms, at all levels of integration and programming to achieve higher performance and energy e ciency. Being the geometric calibration process one of the most time consuming processes when using remote sensing images, the aim of this work is to accelerate this process by taking advantage of new computing architectures and technologies, specially focusing in exploiting computation over shared memory multi-threading hardware. A parallel implementation of the most time consuming process in the remote sensing geometric correction has been implemented using OpenMP directives. This work compares the performance of the original serial binary versus the parallelized implementation, using several multi-threaded modern CPU architectures, discussing about the approach to nd the optimum hardware for a cost-e ective execution.

A novel semi-fragile forensic watermarking scheme for remote sensing images

Peer-reviewed

Los Mapas de riesgo de inundaciones: representación de la vulnerabilidad y aportación de las innovaciones tecnológicas

Los mapas de riesgo de inundaciones deberían mostrar las inundaciones en relación con los impactos potenciales que éstas pueden llegar a producir en personas, bienes y actividades. Por ello, es preciso añadir el concepto de vulnerabilidad al mero estudio del fenómeno físico. Así pues, los mapas de riesgo de daños por inundación son los verdaderos mapas de riesgo, ya que se elaboran, por una parte, a partir de cartografía que localiza y caracteriza el fenómeno físico de las inundaciones, y, por la otra, a partir de cartografía que localiza y caracteriza los elementos expuestos. El uso de las llamadas «nuevas tecnologías», como los SIG, la percepción remota, los sensores hidrológicos o Internet, representa un potencial de gran valor para el desarrollo de los mapas de riesgo de inundaciones, que es, hoy por hoy, un campo abierto a la investigación

Linear and nonlinear terrain deformation maps from a reduced set of interferometric SAR images

In this paper, an advanced technique for the generation of deformation maps using synthetic aperture radar (SAR) data is presented. The algorithm estimates the linear and nonlinear components of the displacement, the error of the digital elevation model (DEM) used to cancel the topographic terms, and the atmospheric artifacts from a reduced set of low spatial resolution interferograms. The pixel candidates are selected from those presenting a good coherence level in the whole set of interferograms and the resulting nonuniform mesh tessellated with the Delauney triangulation to establish connections among them. The linear component of movement and DEM error are estimated adjusting a linear model to the data only on the connections. Later on, this information, once unwrapped to retrieve the absolute values, is used to calculate the nonlinear component of movement and atmospheric artifacts with alternate filtering techniques in both the temporal and spatial domains. The method presents high flexibility with respect to the required number of images and the baselines length. However, better results are obtained with large datasets of short baseline interferograms. The technique has been tested with European Remote Sensing SAR data from an area of Catalonia (Spain) and validated with on-field precise leveling measurements.

Textural classification of land cover using support vector machines: an empirical comparison with parametric, non parametric and hybrid classifiers in the Bolivian Amazon

Land cover classification is a key research field in remote sensing and land change science as thematic maps derived from remotely sensed data have become the basis for analyzing many socio-ecological issues. However, land cover classification remains a difficult task and it is especially challenging in heterogeneous tropical landscapes where nonetheless such maps are of great importance. The present study aims to establish an efficient classification approach to accurately map all broad land cover classes in a large, heterogeneous tropical area of Bolivia, as a basis for further studies (e.g., land cover-land use change). Specifically, we compare the performance of parametric (maximum likelihood), non-parametric (k-nearest neighbour and four different support vector machines - SVM), and hybrid classifiers, using both hard and soft (fuzzy) accuracy assessments. In addition, we test whether the inclusion of a textural index (homogeneity) in the classifications improves their performance. We classified Landsat imagery for two dates corresponding to dry and wet seasons and found that non-parametric, and particularly SVM classifiers, outperformed both parametric and hybrid classifiers. We also found that the use of the homogeneity index along with reflectance bands significantly increased the overall accuracy of all the classifications, but particularly of SVM algorithms. We observed that improvements in producer’s and user’s accuracies through the inclusion of the homogeneity index were different depending on land cover classes. Earlygrowth/degraded forests, pastures, grasslands and savanna were the classes most improved, especially with the SVM radial basis function and SVM sigmoid classifiers, though with both classifiers all land cover classes were mapped with producer’s and user’s accuracies of around 90%. Our approach seems very well suited to accurately map land cover in tropical regions, thus having the potential to contribute to conservation initiatives, climate change mitigation schemes such as REDD+, and rural development policies.

Optimum Redundant Array Configurations for Earth Observation Aperture Synthesis Microwave Radiometers

Two-dimensional aperture synthesis radiometry is the technologyselected for ESA's SMOS mission to provide high resolution L-bandradiometric imagery. The array topology is a Y-shaped structure. Theposition and number of redundant elements to minimise instrumentdegradation in case of element failure(s) are studied.

Enhanced radar precipitation estimates using a combined clutter and beam blockage correction technique

Weather radar observations are currently the most reliable method for remote sensing of precipitation. However, a number of factors affect the quality of radar observations and may limit seriously automated quantitative applications of radar precipitation estimates such as those required in Numerical Weather Prediction (NWP) data assimilation or in hydrological models. In this paper, a technique to correct two different problems typically present in radar data is presented and evaluated. The aspects dealt with are non-precipitating echoes - caused either by permanent ground clutter or by anomalous propagation of the radar beam (anaprop echoes) - and also topographical beam blockage. The correction technique is based in the computation of realistic beam propagation trajectories based upon recent radiosonde observations instead of assuming standard radio propagation conditions. The correction consists of three different steps: 1) calculation of a Dynamic Elevation Map which provides the minimum clutter-free antenna elevation for each pixel within the radar coverage; 2) correction for residual anaprop, checking the vertical reflectivity gradients within the radar volume; and 3) topographical beam blockage estimation and correction using a geometric optics approach. The technique is evaluated with four case studies in the region of the Po Valley (N Italy) using a C-band Doppler radar and a network of raingauges providing hourly precipitation measurements. The case studies cover different seasons, different radio propagation conditions and also stratiform and convective precipitation type events. After applying the proposed correction, a comparison of the radar precipitation estimates with raingauges indicates a general reduction in both the root mean squared error and the fractional error variance indicating the efficiency and robustness of the procedure. Moreover, the technique presented is not computationally expensive so it seems well suited to be implemented in an operational environment.

A wavelet-based method for the determination of the relative resolution between remotely sensed images

Spatial resolution is a key parameter of all remote sensing satellites and platforms. The nominal spatial resolution of satellites is a well-known characteristic because it is directly related to the area in ground that represents a pixel in the detector. Nevertheless, in practice, the actual resolution of a specific image obtained from a satellite is difficult to know precisely because it depends on many other factors such as atmospheric conditions. However, if one has two or more images of the same region, it is possible to compare their relative resolutions. In this paper, a wavelet-decomposition-based method for the determination of the relative resolution between two remotely sensed images of the same area is proposed. The method can be applied to panchromatic, multispectral, and mixed (one panchromatic and one multispectral) images. As an example, the method was applied to compute the relative resolution between SPOT-3, Landsat-5, and Landsat-7 panchromatic and multispectral images taken under similar as well as under very different conditions. On the other hand, if the true absolute resolution of one of the images of the pair is known, the resolution of the other can be computed. Thus, in the last part of this paper, a spatial calibrator that is designed and constructed to help compute the absolute resolution of a single remotely sensed image is described, and an example of its use is presented.

Detección por el satélite LANDSAT-4 de los efectos de la riada del 8-XI-82 en el río Segre

Por medio de técnicas de tratamiento de imágenes digitales se realiza un estudio de los efectos producidos por una inundación ocurrida a finales del año 1982 en el valle del río Segre, en Catalunya, a partir de la información multiespectral captada por el sensor TM del satélite LANDSAT-4. Utilizando un programa de clasificación no supervisada basado en la distancia euclídea, se diferencian cuatro tipos de suelo o de cubiertas en el rea de estudio (3.8 x 2.3 km). Se efecta un análisis cuantitativo de la calidad de los resultados, usando como referencia la información obtenida en un estudio de campo. Este análisis muestra un alto grado de correspondencia entre el mapa de campo (verdad terreno) y la cartografía realizada a partir de los datos multiespectrales.

Sea state effect on the sea surface emissivity at L-band

In May 1999, the European Space Agency (ESA) selected the Earth Explorer Opportunity Soil Moisture and Ocean Salinity (SMOS) mission to obtain global and frequent soil moisture and ocean salinity maps. SMOS' single payload is the Microwave Imaging Radiometer by Aperture Synthesis (MIRAS), an L-band two-dimensional aperture synthesis radiometer with multiangular observation capabilities. At L-band, the brightness temperature sensitivity to the sea surface salinity (SSS) is low, approximately 0.5 K/psu at 20/spl deg/C, decreasing to 0.25 K/psu at 0/spl deg/C, comparable to that to the wind speed /spl sim/0.2 K/(m/s) at nadir. However, at a given time, the sea state does not depend only on local winds, but on the local wind history and the presence of waves traveling from far distances. The Wind and Salinity Experiment (WISE) 2000 and 2001 campaigns were sponsored by ESA to determine the impact of oceanographic and atmospheric variables on the L-band brightness temperature at vertical and horizontal polarizations. This paper presents the results of the analysis of three nonstationary sea state conditions: growing and decreasing sea, and the presence of swell. Measured sea surface spectra are compared with the theoretical ones, computed using the instantaneous wind speed. Differences can be minimized using an "effective wind speed" that makes the theoretical spectrum best match the measured one. The impact on the predicted brightness temperatures is then assessed using the small slope approximation/small perturbation method (SSA/SPM).

Interferometría SAR en las islas Shetlend del Sur: modelo numérico de elevaciones de la Isla Decepción

La interferometría con imágenes de radar de apertura sintética (SAR: Synthetic Aperture Radar) desde satélite es una técnica que permite obtener información altimétrica del relieve terrestre, siendo especialmente útil en áreas remotas. Este trabajo muestra una aplicación de esta técnica en las islas Shetland del Sur (Antártida). Para ello se han utilizado imágenes SAR obtenidas por los satélites ERS (European Remote Sensing) de la Agencia Espacial Europea (ESA: European Space Agency) y un proceso interferométrico desarrollado entre el Departament de Geodinàmica i Geofísica de la Universitat de Barcelona y el Institut Cartogràfic de Catalunya.

Practices of place-making through locative media artworks

Peer-reviewed

Circumsolar energetic particle distribution on 2011 November 3

Late on 2011 November 3, STEREO-A, STEREO-B, MESSENGER, and near-Earth spacecraft observed an energetic particle flux enhancement. Based on the analysis of in situ plasma and particle observations, their correlation with remote sensing observations, and an interplanetary transport model, we conclude that the particle increases observed at multiple locations had a common single source active region and the energetic particles filled a very broad region around the Sun. The active region was located at the solar backside (as seen from Earth) and was the source of a large flare, a fast and wide coronal mass ejection, and an EIT wave, accompanied by type II and type III radio-emission. In contrast to previous solar energetic particle events showing broad longitudinal spread, this event showed clear particle anisotropies at three widely separated observation points at 1AU, suggesting direct particle injection close to the magnetic footpoint of each spacecraft, lasting for several hours.We discuss these observations and the possible scenarios explaining the extremely broad particle spread for this event.