A New Combined Assessment of Mixed Uncertainty in Spatial Models: Conceptualization and Implementation

Uncertainty quantification is not often performed in spatial modeling applications, especially when there is a mixture of probabilistic and non‐probabilistic uncertainties. Furthermore, the effect of positional uncertainty is often not assessed, despite its relevance to geographical applications. Although there has been much work in investigating the aforementioned types of uncertainty in isolation, combined approaches have not been much researched. This has resulted in a lack of tools for conducting mixed uncertainty analyses that include positional uncertainty. This research addresses the issue by first presenting a new, flexible, simulation‐oriented conceptualization of positional uncertainty in geographic objects called F‐Objects. F‐Objects accommodates various representations of uncertainty, while remaining conceptually simple. Second, a new Python‐based framework is introduced, termed Wiggly and capable of conducting mixed uncertainty propagation using fuzzy Monte Carlo simulation (FMCS). FMCS combines both traditional Monte Carlo with fuzzy analysis in a so‐called hybrid approach. F‐Objects is implemented within the Wiggly framework, resulting in a tool capable of considering any combination of: (1) probabilistic variables; (2) fuzzy variables; and (3) positional uncertainty of objects (probabilistic/fuzzy). Finally, a realistic GIS‐based groundwater contamination problem demonstrates how F‐Objects and Wiggly can be used to assess the effect of positional uncertainty.     

A Novel Approach of Despeckling SAR Images Using Nonlocal Means Filtering

Synthetic aperture radar (SAR) is a day and night, all weather satellite imaging technology. Inherent property of SAR image is speckle noise which produces granular patterns in the image. Speckle noise occurs due to the interference of backscattered echo from earth’s rough surface. There are various speckle reduction techniques in spatial domain and transform domain. Non local means filtering (NLMF) is the technique used for denoising which uses Gaussian weights. In NLMF algorithm, the filtering is performed by taking the weighted mean of all the pixels in a selected search area. The weight given to the pixel is based on the similarity measure calculated as the weighted Euclidean distance over the two windows. Non local means filtering smoothes out homogeneous areas but edges are not preserved. So a discontinuity adaptive weight is used in order to preserve heterogeneous areas like edges. This technique is called as discontinuity adaptive non local means filtering and is well-adapted and robust in the case of Additive White Gaussian Noise (AWGN) model. But speckle is a multiplicative random noise and hence Euclidean distance is not a good choice. This paper presents evaluation results of using different distance measures for improving the accuracy of the Non local means filtering technique. The results are verified using real and synthetic images and from the results it can be concluded that the usage of Manhattan distance improves the accuracy of NLMF technique. Non local approach is used as a preprocessing or post processing technique for many denoising algorithms. So improving NLMF technique would help improving many of the existing denoising techniques.     

Contextual Spatiospectral Postreconstruction of Cloud-Contaminated Images

A general method has been proposed recently for the contextual reconstruction of cloud-contaminated areas in multitemporal multispectral images. It is based on the idea of making the prediction process learn from information available in the cloud-free neighborhood of contaminated areas. Though promising, this method does not fully exploit all available information, thus leaving room for further methodological enhancements. This letter presents a postreconstruction methodology for improving the contextual reconstruction process by opportunely capturing spatial and spectral correlations characterizing the considered image. In addition, we propose a solution to a problem that has not yet been addressed in the remote sensing literature, i.e., the generation of an error map beside the reconstructed images to provide end-users with helpful indications about reconstruction reliability. Thorough experiments conducted on a multitemporal sequence of Landsat-7 ETM+ images are reported and discussed.     

Unsupervised Change Detection From Multichannel SAR Images

Multichannel synthetic aperture radar (SAR) data present a good potential for environmental monitoring and disaster management, owing both to their insensitivity to atmospheric and sun-illumination conditions, and to the improved discrimination capability they may provide as compared to single-channel SAR. However, this requires accurate and possibly automatic techniques to generate change maps from multichannel SAR images acquired from the same geographic area at different times. In this letter, an automatic unsupervised contextual change-detection method is proposed for two-date multichannel SAR images, by integrating a SAR-specific extension of the Fisher transform with a variant of the expectation-maximization algorithm and with Markov random fields. The method is validated by experiments on SIR-C/XSAR data     

Testing satellite and ground thermal imaging of low-temperature fumarolic fields: The dormant Nisyros Volcano (Greece)

The Nisyros Volcano (Greece) was monitored by satellite and ground thermal imaging during the period 2000–2002. Three night-scheduled Landsat-7 ETM⁺ thermal (band 6) images of Nisyros Island were processed to obtain land surface temperature. Ground temperature data were also collected during one of the satellite overpasses. Processed results involving orthorectification and 3-D atmospheric correction clearly show the existence of a thermal anomaly inside the Nisyros Caldera. This anomaly is associated mainly with the largest hydrothermal craters and has land surface temperatures 5–10 °C warmer than its surroundings. The ground temperature generally increased by about 4 °C inside the main crater over the period 2000–2002. Ground thermal images of the hydrothermal Stephanos Crater were also collected in 2002 using a portable thermal infrared camera. These images were calibrated to ground temperature data and orthorectified. A difference of about 0–2 °C was observed between the ground thermal images and the ground temperature data. The overall study demonstrates that satellite remote sensing of low-temperature fumarolic fields within calderas can provide a reliable long-term monitoring tool of dormant volcanoes that have the potential to reactivate. Similarly, a portable thermo-imager can easily be deployed for real-time monitoring using telemetric data transfer. The operational costs for both systems are relatively low for an early warning system.     

Evaluating a thermal image sharpening model over a mixed agricultural landscape in India

Fine spatial resolution (e.g., <300 m) thermal data are needed regularly to characterise the temporal pattern of surface moisture status, water stress, and to forecast agriculture drought and famine. However, current optical sensors do not provide frequent thermal data at a fine spatial resolution. The TsHARP model provides a possibility to generate fine spatial resolution thermal data from coarse spatial resolution (≥1 km) data on the basis of an anticipated inverse linear relationship between the normalised difference vegetation index (NDVI) at fine spatial resolution and land surface temperature at coarse spatial resolution. The current study utilised the TsHARP model over a mixed agricultural landscape in the northern part of India. Five variants of the model were analysed, including the original model, for their efficiency. Those five variants were the global model (original); the resolution-adjusted global model; the piecewise regression model; the stratified model; and the local model. The models were first evaluated using Advanced Space-borne Thermal Emission Reflection Radiometer (ASTER) thermal data (90 m) aggregated to the following spatial resolutions: 180 m, 270 m, 450 m, 630 m, 810 m and 990 m. Although sharpening was undertaken for spatial resolutions from 990 m to 90 m, root mean square error (RMSE) of <2 K could, on average, be achieved only for 990–270 m in the ASTER data. The RMSE of the sharpened images at 270 m, using ASTER data, from the global, resolution-adjusted global, piecewise regression, stratification and local models were 1.91, 1.89, 1.96, 1.91, 1.70 K, respectively. The global model, resolution-adjusted global model and local model yielded higher accuracy, and were applied to sharpen MODIS thermal data (1 km) to the target spatial resolutions. Aggregated ASTER thermal data were considered as a reference at the respective target spatial resolutions to assess the prediction results from MODIS data. The RMSE of the predicted sharpened image from MODIS using the global, resolution-adjusted global and local models at 250 m were 3.08, 2.92 and 1.98 K, respectively. The local model consistently led to more accurate sharpened predictions by comparison to other variants.     

Simulation of the Ku-band Radar altimeter sea ice effective scattering surface

A radiative transfer model is used to simulate the sea ice radar altimeter effective scattering surface variability as a function of snow depth and density. Under dry snow conditions without layering these are the primary snow parameters affecting the scattering surface variability. The model is initialized with in situ data collected during the May 2004 GreenIce ice camp in the Lincoln Sea (73/spl deg/W; 85/spl deg/N). Our results show that the snow cover is important for the effective scattering surface depth in sea ice and thus for the range measurement, ice freeboard, and ice thickness estimation.     

The Global Impact of Clouds on the Production of MODIS Bidirectional Reflectance Model-Based Composites for Terrestrial Monitoring

A global data set of cloud occurrence probability derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua gridded daily data is analyzed to investigate the probability of obtaining at least a minimum number of cloud-free observations within various compositing periods. The probabilities derived from Terra and Aqua, with morning and afternoon overpass times, respectively, are similar and increase with compositing period. Compositing both Terra and Aqua observations results in considerably higher probabilities of obtaining a sufficient number of observations for bidirectional reflectance model-based compositing. Given that the only alternative to obtaining sufficient samples is to extend the observation period, which can cause significant problems when the surface state changes, it is concluded that using data from the two MODIS sensors provides the most effective way of generating composited products. Findings with respect to the availability of cloud-free composites when n-day composites are generated on a temporally overlapping daily rolling basis, i.e., every day, rather than every n-days, are also discussed for regional and global applications