This research introduces Genus-Physiognomy-Ecosystem (GPE) mapping at a prefecture level through machine learning of multi-spectral and multi-temporal satellite images at 10m spatial resolution, and later integration of prefecture wise maps into country scale for dealing with 88 GPE types to be classified from a large size of training data involved in the research effectively. This research was made possible by harnessing entire archives of Level-2A product, Bottom of Atmosphere reflectance images collected by MultiSpectral Instruments onboard a constellation of two polar-orbiting Sentinel-2 mission satellites. The satellite images were pre-processed for cloud masking and monthly median composite images consisting of 10 multi-spectral bands and 7 spectral indexes were generated. The ground truth labels were extracted from extant vegetation survey maps by implementing systematic stratified sampling approach and noisy labels were dropped out for preparing a reliable ground truth database. Graphics Processing Unit (GPU) implementation of Gradient Boosting Decision Trees (GBDT) classifier was employed for classification of 88 GPE types from 204 satellite features. The classification accuracy computed with 25% test data varied from 65-81% in terms of F1-score across 48 prefectural regions. This research produced seamless maps of 88 GPE types first time at a country scale with an average 72% F1-score. In addition, mapping of solar panels and vegetation disturbance are added.

The current study aims to apply and compare the performance of six machine learning algorithms, including three basic classifiers: random forest (RF), gradient boosting decision tree (GBDT), and extreme gradient boosting (XGB), as well as their hybrid classifiers, using the logistic regression (LR) method (RF+LR, GBDT+LR, and XGB+LR), in order to map the landslide susceptibility of Zhangjiajie City, Hunan Province, China. First, a landslide inventory map was created with 206 historical landslide points and 412 non-landslide points, which was randomly divided into two datasets for model training (80%) and model testing (20%). Second, 15 landslide conditioning factors (i.e., altitude, slope, aspect, plane curvature, profile curvature, relief, roughness, rainfall, topographic wetness index (TWI), normalized difference vegetative index (NDVI), distance to roads, distance to rivers, land use/land cover (LULC), soil texture, and lithology) were initially selected to establish a landslide factor database. Thereafter, the multicollinearity test and information gain ratio (IGR) technique were applied to rank the importance of the factors. Subsequently, we used a series of metrics (e.g., accuracy, precision, recall, f-measure, area under the ROC (receiver operating characteristic) curve (AUC), kappa index, mean absolute error (MAE), and root mean square error (RMSE)) to evaluate the accuracy and performance of the six models. Based on the AUC values derived from the models, the GBDT+LR model with the highest AUC value (0.8168) was identified as the most efficient model for mapping landslide susceptibility, followed by the XGB+LR, XGB, RF+LR, GBDT, and RF models, which achieved AUC values of 0.8124, 0.8118, 0.8060, 0.7927, and 0.7883, respectively. The results from this study suggest that the stacking ensemble machine learning method is promising for use in landslide susceptibility mapping in the Zhangjiajie area and is capable of targeting the areas prone to landslides.

Remote sensing technology, especially using satellite images, has become essential support in many aspects of decision-making, particularly in disaster risk management. It requires a shorter period of data updates and less cost compared to conventional field observations and surveys. Yet, the intensive processing and high-powered computing resources are necessary to analyze satellite imagery data through Geographic Information System (GIS). In this paper, we introduce the identification and mapping of natural disaster impact in Indonesia using the open-source collaborative tool of Google Earth Engine (GEE) application which analyzes the relative temporal difference of Earth surface from three major satellite images: Sentinel-1, Sentinel-2, and Landsat-8. Taking the advantage of the geographical, geological, and demographic conditions of Indonesia's disaster-prone areas, we analyze relative difference from normalized difference vegetation index (NDVI) out of months before and after natural disaster occurrence to measure the impact of natural disaster in focus study areas. Given the high-vegetation nature of three main natural disaster impacted areas in Indonesia: Aceh, Palu, and Yogyakarta, we are able to simplify the analysis by highlighting areas with vegetative loss or gain after the event. Using an open-source GEE application, namely HazMapper, we identify and visualize the aftermath of the tsunami disaster in Aceh and Palu as well as the earthquake in Yogyakarta. Our study is potentially beneficial for government and decision-makers to utilize publicly available satellite images for disaster recovery and mitigation policy.

At present, indoor localization is one of the core technologies of location-based services (LBS), and there exist numerous scenario-oriented application solutions. Visual features, as the main semantic information to help people understand the environment and thus occupy the dominant part, many techniques about indoor scene recognition are widely adopted. However, the engineering application problem of cell phone indoor scene recognition and localization has not been well solved due to insufficient semantic constraint information of building map and the immaturity of building map location anchors (MLA) matching positioning technology. To address the above problems, this paper proposes a cell phone indoor scene recognition and localization method with building map semantic constraints. Firstly, we build a library of geocoded entities for building map location anchors (MLA), which can provide users with "immersive" real-world building maps on the one hand and semantic anchor point constraints for cell phone positioning on the other. Secondly, using the improved YOLOv5s deep learning model carried on the mobile terminal, we recognize the universal map location anchors (MLA) elements in building scenes by cell phone camera video in real-time. Lastly, the spatial location of the scene elements obtained from the cell phone video recognition is matched with the building MLA to achieve real-time positioning and navigation. The experimental results show that the model recognition accuracy of this method is above 97.2%, and the maximum localization error is within the range of 0.775 m, and minimized to 0.5 m after applying the BIMPN road network walking node constraint, which can effectively achieve high positioning accuracy in the building scenes with rich MLA element information. In addition, the building map location anchors (MLA) has universal characteristics, and the positioning algorithm based on scene element recognition is compatible with the extension of indoor map data types, so this method has good prospects for engineering applications.

Abstract: In the last five years, the urban development of Limassol City has rapidly increased in the sectors of industry, trade, real estate, and many others. This exponentially increased urban development introduces several concerns about the aggravation of the land subsidence in the Limassol coastal front. Fifty Copernicus Sentinel-1 data from 2017-2021 have been processed and analyzed using the Sentinel Application Platform (SNAP) and the Stanford Method for Persistent Scatters (StaMPS). A case study for the identification and analysis of the elements (PS) in pixels in a series of interferograms, and then, the quantity of the land displacements in the Line of Sight, in the Limassol coastal front, is presented in this research, with the subsidence rates up to about (-5 to 4 mm / year). For the validation of the detected deformation, accurate ground-based geodetic measurements along the coastal area were used. Concordantly, taking into account that there are a significant number of skyscrapers planned to be built, this study attempts a preliminary assessment of the impact these structures will pose on the coastal front of the area of Limassol.

Obtaining classification confidence at the pixel level is a challenging task for accuracy assessment in remote sensing image classification. Among the various methods for estimating classification confidence at the pixel level, interpolation-based methods have drawn special attention in the literature. Even though they have been widely recognized in the literature, their usefulness has not been rigorously evaluated. This paper conducts a comprehensive evaluation of three interpolation-based methods: local error matrix method, bootstrap method, and geostatistical method. We applied each of the three methods to three representative datasets with different spatial resolutions, spectral bands, and the number of classes. We then derive the estimated classification confidence and true classification confidence and compared the results with each other using both exploratory data analysis (bi-histogram) and statistical analysis (Willmott's d and Binned classification quality). The results indicate that the three interpolation methods provide some interesting insights on various aspects of estimating per-pixel classification confidence. Unfortunately, the interpolation assumes that classification confidence is smooth across the space, which is usually not true in practice. In other words, interpolation-based methods have limited practical use.

Satellite Crop Detection technologies are focused on detection of different types of crops on the field in the early stage before harvesting. Crop detection is usually done on a time series of satellite data by classification of the desired fields. Currently, data obtained from Remote Sensing (RS) are used to solve tasks related to the identification of the type of agricultural crops, also modern technologies using AI methods are desired in the postprocessing part. In this challenge Sentinel-1 and Sentinel-2 time series data were used due to their periodic availability. Our focus was to develop methodology for classification of time series of Sentinel 2 and Sentinel 1 data and compare how accuracy of classification can be increased, but also how to guarantee availability of data. We analyse phenology of single crops and on the basis of this analysis we started to provide crop classification. Original crop classifications were made from Enhanced Vegetation Index (EVI) layers made from Sentinel-2 time-series data and then we added also . To increase accuracy we also integrate into the process parcel borders and provide classification of fields..

Classification and mapping of plant communities is an essential step for conservation and management of ecosystems and biodiversity. We adopt the Genus-Physiognomy-Ecosystem (GPE) system developed in previous study for satellite-based classification of plant communities. This paper assesses the potential of multi-spectral and multi-temporal images collected by Sentinel-2 satellites. This research was conducted in five representative study sites in a temperate region. It consists of 44 types of plant communities including a few land cover types as well. The plant community types were enumerated in the study sites and ground truth data were prepared with reference to extant vegetation surveys, visual interpretation of high-resolution images, and onsite field observations. We acquired all Sentinel-2 Level-1C product images available for the study sites between 2017-2019 and generated monthly median composite images consisting of ten spectral and twelve spectral-indices. Gradient Boosting Decision Trees (GBDT) classifier was employed as an efficient and distributed gradient boosting technique for the supervised classification of big datasets involved in the research. The cross-validation accuracy in terms of kappa coefficient varied from 87% in Oze site with 41 land cover and plant community types to 95% in Hakkoda site with 19 land cover and plant community types; with average performance of 91% across all sites. In addition, the resulting maps demonstrated a clear distribution of plant community types involved in all sites, highlighting the potential of Sentinel-2 multi-spectral and multi-temporal images with GPE classification system for operational and broad-scale mapping of land cover and plant communities.

Deep learning techniques are increasingly being recognized as effective image classifiers. Aside from their successful performance in past studies, the accuracies have varied in complex environments in comparison with the popularly applied machine learning classifiers. This study seeks to explore the feasibility for using a U-Net deep learning architecture to classify bi-temporal high resolution county scale aerial images to determine the spatial extent and changes of land cover classes that directly or indirectly impact tidal marsh. The image set used in the analysis is a collection of a 1-m resolution collection of National Agriculture Imagery Program (NAIP) tiles from 2009 and 2019 covering Beaufort County, South Carolina. The U-net CNN classification results were compared with two machine learning classifiers, the Random Trees (RT) and the Support Vector Machine (SVM). The results revealed a significant accuracy advantage in using the U-Net classifier (92.4%) as opposed to the SVM (81.6%) and RT (75.7%) classifiers for overall accuracy. From the perspective of a GIS analyst or coastal manager, the U-Net classifier is now an easily accessible nad powerful tool for mapping large areas. Change detection analysis indicated little areal change on marsh extent, though increased land development throughout the county has the potential to negatively impact the health of the marshes. Future work should explore applying the constructed U-Net classifier to coastal environments in large geographic areas, while also implementing other data sources (e.g., LIDAR, multispectral data) to enhance classification accuracy.

The aerosol optical depth is an important indicator of aerosol particle properties and associated radiative impacts. AOD determination is therefore very important to achieve relevant climate modeling. Most remote sensing techniques to retrieve aerosol optical depth are applicable to daytime given the high level of light available. The night represents half of the time but in such conditions only a few remote sensing techniques are available. Among these techniques, the most reliable are moon photometers and star photometers. In this paper, we attempt to fill gaps in the aerosol detection performed with the aforementioned techniques using night sky brightness measurements during moonless nights with the novel CoSQM: a portable, low cost and open-source multispectral photometer. In this paper, we present an innovative method for estimating the aerosol optical depth by using an empirical relationship between the zenith night sky brightness measured at night with the CoSQM and the aerosol optical depth retrieved at daytime from the AErosol Robotic NETwork. Such a method is especially suited to light-polluted regions with light pollution sources located within a few kilometers of the observation site. A coherent day-to-night aerosol optical depth and Ångström Exponent evolution in a set of 354 days and nights from August 2019 to February 2021 was verified at the location of Santa Cruz de Tenerife on the island of Tenerife, Spain. The preliminary uncertainty of this technique was evaluated using the variance under stable day-to-night conditions, set at 0.02 for aerosol optical depth and 0.75 for Ångström Exponent. These results indicate the set of CoSQM and the proposed methodology appear to be a promising tool to add new information on the aerosol optical properties at night, which could be of key importance to improve climate predictions.

Hyper-spectral images contain a wide range of bands or wavelength due to which they are rich in information. These images are taken by specialized sensors and then investigated through various supervised or unsupervised learning algorithms. Data that is acquired by hyperspectral image contain plenty of information hence it can be used in applications where materials can be analyzed keenly, even the smallest difference can be detected on the basis of spectral signature i.e. remote sensing applications. In order to retrieve information about the concerned area, the image has to be grouped in different segments and can be analyzed conveniently. In this way, only concerned portions of the image can be studied that have relevant information and the rest that do not have any information can be discarded. Image segmentation can be done to assort all pixels in groups. Many methods can be used for this purpose but in this paper, we discussed k means clustering to assort data in AVIRIS cuprite, AVIRIS Muffet and Rosis Pavia in order to calculate the number of regions in each image and retrieved information of 1st, 10th and100th band. Clustering has been done easily and efficiently as k means algorithm is the easiest approach to retrieve information.

This paper investigates the characteristics and impact of a major Saharan dust storm during June 14^th -19^th 2020 to atmospheric radiative and thermodynamics properties over the Atlantic Ocean. The event witnessed the highest ever aerosol optical depth (close to 2 during the peak of the storm) for June since 2002. The satellites and high-resolution model reanalysis products well captured the origin, spread and the effects of the dust storm. The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) profiles, lower angstrom exponent values (~ 0.12) and higher aerosol index value (> 4) tracked the presence of elevated dust. It was found that the dust AOD was as much as 250-300% higher than their climatology resulting in an atmospheric radiative forcing ~200% larger. As a result, elevated warming ( 8-16 %) was observed, followed by a drop in relative humidity(2-4%) in the atmospheric column, as evidenced by both in-situ and satellite measurements. Quantifications such as these for extreme dust events provide significant insights that may help in understanding their climate effects, including improvements to dust simulations using chemistry-climate models

Abstract: Spatially explicit, fine-grained datasets describing historical urban extents are rarely available prior to the era of operational remote sensing. However, such data are necessary to better understand long-term urbanization and land development processes and for the assessment of coupled nature-human systems, e.g., the dynamics of the wildland-urban interface. Herein, we propose a framework that jointly uses remote sensing derived human settlement data (i.e., the Global Human Settlement Layer, GHSL) and scanned, georeferenced historical maps to automatically generate historical urban extents for the early 20th century. By applying unsupervised color segmentation to the historical maps, spatially constrained to the urban extents derived from the GHSL, our approach generates historical settlement extents for seamless integration with the multi-temporal GHSL. We apply our method to study areas in countries across four continents, and evaluate our approach for two U.S. study sites against historical settlement extents derived from the Historical Settlement Data Compilation for the US, HISDAC-US, achieving Area-under-the-Curve values >0.9. Our results are largely in agreement with model-based urban areas from the HYDE database, and demonstrate that the integration of remote sensing derived observations and historical cartographic data sources opens up new, promising avenues for assessing urbanization, and long-term land cover change in countries where historical maps are available.

Hyperspectral image (HSI) classification is a mechanism of analyzing differentiated land cover in remotely sensed hyperspectral images. In the last two decades, a number of different types of classification algorithms have been proposed for classifying hyperspectral data. These algorithms include supervised as well as unsupervised methods. Each of these algorithms has its own limitations. In this research, three different types of unsupervised classification methods are used to classify different datasets i-e Pavia Center, Pavia University, Cuprite, Moffett Field. The main objective is to assess the performance of all three classifiers K-Means, Spectral Matching, and Abundance Mapping, and observing their applicability on different datasets. This research also includes spectral feature extraction for hyperspectral datasets.

Surface temperature influences human health directly and alters the biodiversity and productivity of the environment. While previous research has identified that the composition of urban landscapes influences the physical properties of the environment such as surface temperature, a generalizable and flexible framework is needed that can be used to compare cities across time and space. This study employs the Structure of Urban Landscapes (STURLA) classification combined with remote sensing of New York City’s (NYC) surface temperature. These are then linked using machine learning and statistical modeling to identify how greenspace and the built environment influence urban surface temperature. It was observed that areas with urban units composed of largely the built environment hosted the hottest temperatures while those with vegetation and water were coolest. Likewise, this is reinforced by borough-level spatial differences in both urban structure and heat. Comparison of these relationships over the period between2008 and 2017 identified changes in surface temperature that are likely due to the changes in prevalence in water, lowrise buildings, and pavement across the city. This research reinforces how human alteration of the environment changes ecosystem function and offers units of analysis that can be used for research and urban planning.

Significant earthquakes frequently occur in Indonesia. Indonesia is situated over three active tectonic plates, resulting in the formation of faults and trenches on the land and ocean floor. For the last 120 years since 1900, there have been more than 1,250 significant earthquake events in Indonesia. In this study, we analyse Indonesia's significant earthquake events using geostatistical and geovisualisation methods to produce an appropriate geospatial analysis platform using the RShiny package to build the WebGIS application. The results show that the earthquake events were spatially distributed from the Sumatera fault in the western part of Indonesia to the southern part of Indonesia, where the Java trench was located and the eastern part of Indonesia. The WebGIS application received a positive evaluation by respondents, with a mean value of 1.617 for pragmatic quality, 1.808 for hedonic quality, and 1.713 for overall quality. This means that the WebGIS application is of good quality based on respondents' impressions. The users also more easily gained insight into information as a result of geostatistical methods. The information gained by the users during the user interaction with the WebGIS platform overlapped with the information that the researcher started with, that is, the spatial cluster of significant earthquakes in Indonesia.

Georeferencing using ground control points (GCPs) is the most common strategy in photogrammetry modeling using UAV-acquired imagery. However, with the increased availability of UAVs with onboard GNSS RTK, georeferencing without GCPs is a promising alternative. However, systematic elevation error remains a problem of this technique. We aimed to analyze the reasons for this systematic error and propose strategies for the elimination of this error. Multiple flights differing in the flight altitude and image acquisition axis were performed at two real-world sites. A flight height of 100m with vertical (nadiral) image acquisition axis was considered primary, supplemented with flight altitudes of 75 m and 125 m with vertical image acquisition axis and two flights at 100 m with oblique image acquisition axes (30° a 15°). Each of these flights was performed twice to produce a full double grid. Models were calculated from individual flights and their combinations. The elevation error from individual flights or even combinations yielded systematic elevation errors of up to several decimeters. This error was linearly dependent on the deviation of the focal length from the reference value. A combination of two flights from the same altitude (with nadiral and oblique image acquisition) was capable of reducing the systematic elevation error to less than 0.03 m. This study is the first to demonstrate the linear dependence between the systematic elevation error of the models based only on the onboard GNSS-RTK data and the deviation in the determined internal orientation parameters (focal length). Besides, we have shown that a combination of two flights with different image acquisition axis can eliminate this systematic error even in real-world conditions and that georeferencing without GCPs is, therefore, a feasible alternative to the use of GCPs.

In cotton cultivars, an insect that causes irreversible damage is the Spodoptera frugiperda, known as the fall armyworm. Since the visual detection of plants is a burdensome task for human inspection, the spectral information related to plant damage, registered on a spectral scale, can be useful. These measurements, associated with machine learning techniques, produce useful information for a rapid and non-invasive inspection method development. To contribute to this gap fulfillment, this paper proposes a machine learning framework to model the spectral response of cotton plants under the attack of the fall armyworm. Additionally, a theoretical model is presented, built from the results of the machine learning analysis, to infer this damage with up-to-date orbital sensors. The data was composed of the reflectance measurements collected at a cotton field with control plants and plants submitted to Spodoptera frugiperda damage. Their spectral response was recorded with a hand-held spectroradiometer ranging from 350 to 2,500 nm, for eight consecutive days. Different machine learning models were evaluated and the overall best model was defined by accuracies comparisons on a testing-set. A ranking approach was adopted based on the model accuracy, returning the most contributive wavelengths for the classification. Sequentially, an unsupervised neural network (Self-Organizing Map - SOM) was implemented to reduce data-dimensionality and assist in the definition of important spectral regions. The regions were associated with the spectral bands of the two sensors (OLI and MSI) and a theoretical model using a band simulation process with the overall best machine learning model was proposed. The results indicated that the Random Forest (RF) algorithm is the most suitable to predict cotton-plants damaged by insects and that the last day of analysis (8th day) was better to separate it, with F-measure equals 0.912. The ranking approach combined with the SOM method indicated the spectral regions at the red to near-infrared (650 to 1,350 nm) and shortwave infrared (1,570 to 1,640 nm) as the most important regions to the analysis. The proposed theoretical model simulated with the OLI and MSI sensor-bands returned an F-Measure of 0.865 and 0.886, respectively. In conclusion, this framework can be used to map cotton-plants under insect-attack. The theoretical model presents high accuracy to infer the insect-damaged on cotton plants based on multispectral bands from other sensors, being a useful tool for future research that intends to evaluate it in other areas and at different field scales.

This paper presents ensemble learning of multi-source satellite sensors dataset to obtain better predictive performance of the forest biomass. Spectral, spectral-indices, and spectral-textural features were generated from two optical satellite sensors, Landsat 8 Operational Land Imager (OLI) and Sentinel-2 Multispectral Instrument (MSI). In addition, two radar satellite sensors, Sentinel-1 C-band Synthetic Aperture Radar (CSAR), and Advanced Land Observing Satellite (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) were utilized to generate backscattering and backscattering-textural features. The plot-wise above ground biomass data available from five forests in New England region were utilized. Ensemble learning of multi-source satellite sensors dataset was carried out by employing four machine learning regressors namely, Support Vector Machines (SVM), Random Forests (RF), Gradient Boosting (GB), and Multilayer Perceptron (MLP). A five-fold cross-validation method was used to evaluate predictive performance of the multi-source satellite sensors. The integration of multi-source satellite features, comprising of spectral, spectral-indices, backscattering, spectral-textural, and backscattering-textural information, through ensemble learning and cross-validation approach implemented in the research showed promising results (R² = 0.81, RMSE = 46.2 Mg/ha) for the estimation of plots-level forest biomass in New England region.

This paper proposed a new evaluation model based on analytic hierarchy process to quantitatively evaluate the capability of multi-satellite cooperative remote sensing observation. The analytic hierarchical process model is a combination of qualitative and quantitative analysis of systematic decision analysis method. According to the objective of the remote sensing cooperative observation mission, we decompose the complex problem into several levels and a number of factors, compare and calculate various factors in pairs, and obtain the combination weights of different schemes. The model can be used to evaluate the observation capability of resource satellites. Taking the optical remote sensing satellites such as China’s resource satellite series and GF-4 as examples, this paper verifies and evaluates the model for three typical tasks: point target observation, regional target observation and moving target continuous observation. The results show that the model can provide quantitative reference and model support for comprehensive evaluation of the collaborative observation capability of remote sensing satellites.

With the rapid development of earth observation, satellite navigation, mobile communication and other technologies, the order of magnitude of the spatial data we acquire and accumulate is increasing, and higher requirements are put forward for the application and storage of spatial data. Under this circumstance, a new form of spatial data organization emerged-the global discrete grid. This form of data management can be used for the efficient storage and application of large-scale global spatial data, which is a digital multi-resolution the geo-reference model that helps to establish a new model of data association and fusion. It is expected to make up for the shortcomings in the organization, processing and application of current spatial data. There are different types of grid system according to the grid division form, including global discrete grids with equal latitude and longitude, global discrete grids with variable latitude and longitude, and global discrete grids based on regular polyhedrons. However, there is no accuracy evaluation index system for remote sensing images expressed on the global discrete grid to solve this problem. This paper is dedicated to finding a suitable way to express remote sensing data on discrete grids, and establishing a suitable accuracy evaluation system for modeling remote sensing data based on hexagonal grids to evaluate modeling accuracy. The results show that this accuracy evaluation method can evaluate and analyze remote sensing data based on hexagonal grids from multiple levels, and the comprehensive similarity coefficient of the images before and after conversion is greater than 98%, which further proves that the availability hexagonal grid-based remote sensing data of remote sensing images. And among the three sampling methods, the image obtained by the nearest interpolation sampling method has the highest correlation with the original image.

In this study, our goal was to research land-use change by combining spatio–temporal land use/land cover monitoring (LULC (1989–2019) and urban growth modeling (1999–2039) in Benslimane, Morocco, to determine the effect of urban growth on different groups based on cellular automata (CA) and geospatial methods. A further goal was to test the reliability of the AC algorithm for urban expansion modeling. To do this, four years of satellite data were used at the same time as population density, downtown distance, slope, and ground road distance. The LULC satellite reported a rise of 3.8 km2 (318% variation) during 1989–2019. Spatial transformation analysis reveals a good classification similarity ranging from 89% to 91% with the main component analysis (PCA) technique. The statistical accuracy between the satellite scale and the replicated built region of 2019 gave 97.23 %t of the confusion matrix overall accuracy, and the region under the receiver operational characteristics (ROC) curve to 0.94, suggesting the model's high accuracy. Although the constructed area remains low relative to the total area of the municipality's territory, the LULC project shows that the urban area will extend to 5,044 km2 in 2019, principally in the western and southwestern sections. In 2019–2039, urban development is expected to lead to a transformation of the other class (loss of 1,364 km2), followed by vegetation cover (loss of 0.345 km2). In spatial modeling and statistical calculations, the GDAL and NumPy Python 3.8 libraries were successful.

There are many new land use and land cover (LULC) products emerging yet there is still a lack of in-situ data for training, validation, and change detection purposes. The LUCAS (Land Use Cover Area frame Sample) survey is one of the few authoritative in-situ field campaigns, which takes place every three years in European Union member countries. More recently, a study has considered whether citizen science and crowdsourcing could complement LUCAS survey data, e.g., through the FotoQuest Austria mobile app and crowdsourcing campaign. Although the data obtained from the campaign were promising when compared with authoritative LUCAS survey data, there were classes that were not well classified by the citizens, and the photographs submitted through the app were not always of sufficient quality. For this reason, in the latest FotoQuest Go Europe 2018 campaign, several improvements were made to the app to facilitate interaction with the citizens contributing and to improve their accuracy in LULC identification. In addition to extending the locations from Austria to Europe, a change detection component (comparing land cover in 2018 to the 2015 LUCAS photographs) was added, as well as an improved LC decision tree and a near real-time quality assurance system to provide feedback on the distance to the target location, the LULC classes chosen and the quality of the photographs. Another modification was the implementation of a monetary incentive scheme in which users received between 1 to 3 Euros for each successfully completed quest of sufficient quality. The purpose of this paper is to present these new features and to compare the results obtained by the citizens with authoritative LUCAS data from 2018 in terms of LULC and change in LC. We also compared the results between the FotoQuest campaigns in 2015 and 2018 and found a significant improvement in 2018, i.e., a much higher match of LC between FotoQuest Go Europe and LUCAS. Finally, we present the results from a user survey to discuss challenges encountered during the campaign and what further improvements could be made in the future, including better in-app navigation and offline maps, making FotoQuest a model for enabling the collection of large amounts of land cover data at a low cost.

The recent growth of high-resolution spatial data, especially in developing urban environments, is enabling new approaches to civic activism, urban planning and the provision of services necessary for sustainable development. A special area of great potential and urgent need deals with urban expansion through informal settlements (slums). These neighborhoods are too often characterized by a lack of connections, both physical and socioeconomic, with detrimental effects to residents and their cities. Here, we show how a scalable computational approach based on the topological properties of digital maps can identify local infrastructural deficits and propose context-appropriate minimal solutions. We analyze 1 terabyte of OpenStreetMap (OSM) crowdsourced data to create worldwide indices of street block accessibility and local cadastral maps and propose infrastructure extensions with a focus on 120 Low and Middle Income Countries (LMIC) in the Global South. We illustrate how the lack of physical accessibility can be identified in detail, how the complexity and costs of solutions can be assessed and how detailed spatial proposals are generated. We discuss how these diagnostics and solutions provide a multiscalar set of new capabilities – from individual neighborhoods to global regions – that can coordinate local community knowledge with political agency, technical capability, and further research.

The quality of photogrammetric-based derived products like orthophotos, digital terrain models (DTMs) and digital line maps as well as the global digital elevation models (DEM) are rigorously dependent on the accuracy of image orientation. This paper evaluates the vertical accuracy of aerial photogrammetric Digital Terrain Model (DTM), Shuttle Radar Topography Mission (SRTM), Advanced Spaceborne Thermal Emission and Reflectance Radiometer (ASTER), and TerraSAR-X's twin satellite of TanDEM-X (TDX) datasets against in-situ orthometric heights computed from ellipsoidal heights and the 2008 Earth Gravitational Model (EGM2008) derived geoid heights in Ethiopia. The quality of the four global digital elevation models was also validated against the aerial photogrammetric DTM measurements. Besides, the accuracies of the photogrammetric DTM and the four DEM products were checked for their compliance to the American Society for Photogrammetry and Remote Sensing (ASPRS) standards as well as the Ethiopian national vertical data evaluation standards. The study showed that the photogrammetric DTM is in a good agreement with the reference orthometric heights compared to SRTM, ASTER and TDX datasets. More precisely, the result has an absolute accuracy of 1.67 m at Linear Error (LE) 95% confidence level, while the absolute accuracy of SRTM3 arc seconds (SRTM3) at LE 90% (11.91 m) is better than its product specification (16 m). The absolute accuracy of SRTM1 arc second (SRTM1) (7.70 m at LE 90%) surpasses that of SRTM3, whereas the absolute accuracy of ASTER DEM is somehow below its product specification. TDX also has the same vertical accuracy (10.29 m at LE 90%) compared to its product specification (10 m). Furthermore, the vertical accuracy of the photogrammetric DTM meets the100 cm vertical accuracy of the 2015 ASPRS standard. However, it does not meet the Ethiopian national vertical data accuracy requirement standard, i.e., RMSEz of ± 0.45 m. In general, the photogrammetric DTM, SRTM1, and TDX have been proven a superior product over the SRTM3 and ASTER DEMs, and better to use these products for high-level precision and accuracy required applications.

Remote sensing data are essential for monitoring Earth’s surface waters, especially since the number of publicly available in-situ data is declining. Satellite altimetry provides valuable information on water level and its variations for lakes, reservoirs, and rivers. In combination with satellite imagery, the derived time series allow for monitoring lake storage changes and river discharge. However, satellite altimetry is limited in spatial resolution due to its measurement geometry, only providing information in nadir direction beneath the satellite’s orbit. In a case study in the Mississippi River Basin (MRB), this study investigates the potential and the limitations of past and current satellite missions to monitor basin-wide storage changes. For that purpose an automated target detection is developed and the extracted lake surfaces are merged with the satellites’ tracks. This reveals that the current altimeter configuration misses about 80% of all lakes larger than 0.1 km2 in the MRB, and still 20% of lakes larger than 10 km2, corresponding to 30% and 7% of surface area, respectively. Past altimetry configurations perform even worse. From the water bodies represented by a global hydrology model, at least 91% of targets and 98% of storage changes are captured by the current altimeter configuration. This will significantly improve with the launch of the planned SWOT mission.

Over the past 30 years, the use of field spectroscopy has risen in importance in remote sensing studies for the characterization of the surface reflectance of materials in situ within a broad range of applications. Potential uses range from measurements of individual targets of interest (e.g. vegetation, soils, validation targets etc.), to characterizing the contributions of different materials within larger spatially-mixed areas as would be representative of the spatial resolution captured by a sensor pixel (UAV to satellite scale). As such, it is essential that a complete and rigorous assessment of both the data-acquisition procedures, and the suitability of the derived data product be carried out. The measured energy from solar-reflected range spectroradiometers is influenced by the viewing and illumination geometries and the illumination conditions which vary due to changes in solar position and atmospheric conditions. By applying corrections, the estimated absolute reflectance (Rabs) of targets can be calculated. This property is independent of illumination intensity or conditions and is the metric commonly suggested to be used to compare spectra even when data are collected by different sensors or acquired under different conditions. By standardizing the process of estimated Rabs, as is provided in the described toolkit, consistency and repeatability in processing are ensured and the otherwise labor intensive and error-prone processing steps are streamlined. The resultant end data product (Rabs) represents our best current effort to generate consistent and comparable ground spectra which have been corrected for viewing and illumination geometries as well as other factors such as the individual characteristics of the reference panel used during acquisition.

Advanced-Differential Interferometric SAR (A-DInSAR) has been used to monitor surface deformations in open pit mines and tailings dams. In this paper, ground deformations have been detected on the area of the tailings Dam-I at the Córrego do Feijão Mine (Brumadinho, Brazil) before its catastrophic failure occurred on 25 January 2019. Two techniques optimized for different scattering models, SBAS (Small BAseline Subset) and PSI (Persistent Scatterer Interferometry), were used to perform the analysis based on 26 Sentinel-1B images in IW mode, acquired on descending orbits from 03 March 2018 to 22 January 2019. A WorldDEM DSM product was used to remove the topographic phase component. The results provided by both techniques showed a synoptic and informative view of the deformation process affecting the study area, with a detection of persistent trend of deformations on the top, middle and bottom sectors of the dam face until its collapse, as well as the expected natural settlements on the tailings. It is worth noting the detection of an acceleration in the displacement time-series for a short period near the failure. The maximum accumulated displacements detected along the downstream slope face were -39 mm (SBAS) and -48 mm (PSI). It is reasonable to consider that Sentinel-1 would provide decision makers complementary motion information to the in-situ monitoring system for risk assessment and for a better understanding of on-going instability phenomena affecting the tailings dam.

NASA’s ICESat-2space-borne photon-counting lidar mission is providing global elevation measurements that will provide significant benefits to a variety of bio-geoscience research applications. Given the novelty of elevation and the derived data products from the ICESat-2 mission, the research community needs software tools that can facilitate photon-level analyses to support product validation and development new analysis methods. Here, we describe PhotonLabeler, a free graphic user interface (GUI) for manual labeling and visualization of ICESat-2 Geolocated Photon data (ATL03). Developed in MATLAB, the GUI facilitates the reading and display of ATL03 Hierarchical Data Format (HDF) files, the manual labeling of individual photons into target classes of choice using a number of point selections tools and enables eventual saving of labeled data in ASCII format. Other capabilities include saving and loading of labeling sessions to manage labeling tasks over time. We expect labeled data generated using the application to serve two main purposes. First, serve as ground truth for validating various products from ICESat-2 mission, especially for study sites around the world that do not have existing reference datasets such as airborne lidar. Second, serve as training and validation data in the development of new algorithms for generating various ICESat-2 data products. We demonstrate the first use case through a validation case study for the land and vegetation product (ATL08), which provides canopy and terrain height estimates, over two sites. For the first site, located in northwestern Zambia, we used ICESat-2 ATL03 data acquired at night and for our second site in Texas, US, we used ATL03 data acquired during the day. The PhotonLabeler application is freely available as a compiled MATLAB binary to enable free access and utilization by interested researchers.

The information on the land use and soil conservation practice based on year 2006, 2010 and 2014, hence offering an opportunity to model the impacts of land use change on erosion, deposition and surface water runoff. Limitation in the use of hydrological models had been their inability to handle the large amount of input data that describe the heterogeneity of the natural system. In this study, a procedure that takes into account soil conservation practice based on the land use change, the response of soil erosion and sediment export from the George Town Conurbation catchment area, and average annual sediment yields were estimated for each grid cell of the watershed to identify the critical erosion areas of rural and urban planning proposes. Average annual sediment yield and data on a grid basis estimated using Universal Soil Loss Equation (USLE) and an emerging technology represented by Geographic Information System (GIS) used as a tool to produce a map for erosion rate. The changing of the land use from forest to agriculture and then to an urban area is a challenging task to research on land use demand for population, and environmental impact assessment is important for the planning of natural resources management, allowing research the modification of land use properly and implement more sustainable for long term management strategies. The challenge is to formulate strategies that would promote an integrated approach to the land use planning at an appropriate level as to address the issues that arose. Modelling for creating urban growth boundary for the George Town Conurbation must have to be controlled surface runoff and soil loss and sediment export from land use of the George Town Conurbation catchment.

Failaka Island, located in the far east of Kuwait Bay, is about 20 km from the State of Kuwait’s coast, and represents a focal point for regional geography and history, with pristine beaches and archaeological sites dating to the Bronze, Iron, Hellenistic, Christian and Islamic periods. According to environmental data and in coordination with local authorities to develop an urban plan, the island is set to become the first tourist destination for the State of Kuwait. To achieve the Vision of Kuwait 2035, one of the planning objectives focuses around Urban Planning for the Establishment of Environmental Cities that Achieve (UPEECA) environmental sustainability criteria. The paper then, aims to propose the environmental urban plan for Failaka Island. Based around Environmental Analytical Hierarchical Processes (EAHP), and using the Field Calculator and ModelBuilder functions in ArcGIS, this research centres on the feasibility of carrying out an urban plan using suitability modelling that includes four factors and 13 criteria covering the island’s ecological and human composition. This study utilizes both remote sensing (UAVs for 3D imaging) and field study (ground truthing) to identify changes in land use and land cover – such as using sample analysis of the historical sites and soils for tracing evidence and creating/updating a soil map – and create the first GIS database for the island that can lead to generating a suitability model.

Sentinel-1 satellite system continuously observes European countries in a relatively high revisit frequency of 6 days per orbital track. Given the Sentinel-1 configuration, most areas in Czechia are observed every 1–2 days by different tracks in a moderate resolution. This is attractive for various types of analyses by various research groups. The starting point for processing is an original data provided by ESA, for interferometry (InSAR) this level is a Single Look Complex (SLC) data. This work represents advantages of storing data augmented to a specifically corrected level of data, SLC-C. The presented database contains Czech nation-wide Sentinel-1 data stored in burst units that have been preprocessed to the state of a consistent well-coregistered dataset of SLC-C. These are resampled SLC data with their phase values reduced by a topographic phase signature, ready for fast interferometric analyses (an interferogram is generated by a complex conjugate between two stored SLC-C files). The data can be used directly into multitemporal interferometry techniques, e.g. Persistent Scatterers (PS) or Small Baseline (SB) techniques applied here. A further development of the nation-wide system utilising SLC-C data would lead into a dynamic state where every new pre-processed burst triggers a processing update to detect unexpected changes from InSAR time series and therefore provide a signal for early warning against a potential dangerous displacement, e.g. a landslide, instability of an engineering structure or a formation of a sinkhole. An update of the processing chain would also allow use of cross-polarised Sentinel-1 data, needed for polarimetric analyses. The current system is running at a national supercomputing centre IT4Innovations in interconnection to the Czech Copernicus Collaborative Ground Segment (CESNET), providing fast on-demand InSAR results over Czech territories. A full nation-wide PS processing using data over Czechia has been performed in 2017, discovering several areas of land deformation. Its downsampled version and basic findings are demonstrated within the article.

In the last decades droughts, deforestation and wildfires have become recurring phenomena that have affected both human activities and natural ecosystems in Amazonia. The time an ecosystem requires to recover from carbon losses is a crucial metric to evaluate disturbance impacts on forests. However, the factors influencing and controlling the recovery time and its spatiotemporal patterns at the regional scale are still poorly understood. In this study, we combined forest growth model, remote sensing and field plots, to map Amazonia-wide (300-ha resolution) impact and recovery time of aboveground biomass (AGB) after drought, fire and a combination of logging and fire. Our simulated results indicate that AGB decreases by 4%, 19% and 46% in forests disturbed by drought, fire and logging + fire, respectively, with an average AGB recovery time of 27 years for drought, 44 years for burned and 63 years for logged + burned areas and with maximum values reaching 184 years in areas of high fire intensity. Our findings provide two major insights in the spatial and temporal patterns of drought and wildfire in the Amazon: 1) the recovery time of the forests takes longer in the southeastern part of the basin, and, 2) as droughts and wildfires become more frequent – since the intervals between the disturbances is getting shorter than forest regeneration – potentially causing a long-lasting damage in these fragile ecosystems and a permanent degradation.

Solar3D is an open-source software application designed to interactively calculate solar irradiation at three-dimensional (3D) surfaces in a virtual environment constructed with combinations of 3D city models, digital elevation models (DEMs), digital surface models (DSMs) and feature layers. The GRASS GIS r.sun solar radiation model computes solar irradiation based on two-dimensional (2D) raster maps for given day, latitude, surface and atmospheric conditions. With the increasing availability of 3D city models and demand for solar energy, there is an urgent need for better tools to computes solar radiation directly with 3D city models. Solar3D extends GRASS GIS r.sun from 2D to 3D by feeding the model with input, including surface slope, aspect and time-resolved shading, that is derived directly from the 3D scene using computer graphics techniques. To summarize, Solar3D offers several new features which, as a whole, distinguish itself from existing 3D solar irradiation tools: (1) the ability to consume massive heterogeneous 3D city models, including massive 3D city models such as oblique airborne photogrammetry-based 3D city models (OAP3Ds or integrated meshes); (2) the ability to perform near real-time pointwise calculation for duration from daily to annual; (3) the ability to integrate and interactively explore large-scale heterogeneous geospatial data. (4) the ability to calculate solar irradiation at arbitrary surface positions including at rooftops, facades and the ground. Solar3D is publicly available at https://github.com/jian9695/Solar3D.

Debris-covered glaciers are a notable feature in the greater Himalaya, and their ongoing mass loss under changing climate will affect the water resources of over a billion people. The current knowledge of the mass balance of Himalayan glaciers is restricted by the paucity of in-situ measurements of glaciers in both space and time, as well as the resolution of satellite remote sensing imageries. Recently, the use of Unmanned Aerial System (UAS) imagery has shown the potential to bridge this gap by enabling very detailed monitoring of inaccessible glacial areas. UAS imagery-based monitoring of Himalayan glaciers has so far been limited to a single glacier in the entire Himalaya, providing a limited understanding of spatial variability in glacier mass balance and driving factors. In the first UAS based glacial mass change estimation in the trans-Himalaya, we conducted two Unmanned Aerial System (UAS) surveys (May and November 2019) over the debris-covered Annapurna III glacier in the Himalaya. We performed Structure-from-Motion (SfM) analysis and utilized differential GPS field observations to derive geometrically accurate point clouds, ortho-mosaics and digital surface models (DSMs). The glacial volumetric loss was estimated from DSM differencing, and the magnitude and spatial variability of glacier surface change was derived from 3-D differencing of point clouds. Results revealed a heterogeneous glacial melt pattern, with an average elevation loss of 0.89 m during the monitored time period. The majority of the glacial tongue exhibited surface lowering except the area above and around the glacial snout that surprisingly exhibited significant elevation gain. Both the highest magnitude of mass loss and the highest spatial variability in mass change was observed in areas with exposed ice-cliffs and supraglacial ponds. Glacial surface velocity derived from manual feature tracking showed velocity ranging from 0-4.1 m. A detailed evaluation of specific areas allowed an improved understanding of the complex interplay of factors leading to observed surface change. Our findings expand the extent of UAS based monitoring of debris-covered glaciers in the Himalaya and conclude that UAS derived 3D topographic products will become increasingly important for monitoring of thinning debris-covered glaciers.

Remote sensing has been used as an important tool for disaster monitoring and disaster scope extraction, especially for the analysis of spatial and temporal disaster patterns of large-scale and long-duration series. Based on the Google Earth Engine cloud platform, this study used MODIS vegetation index products with 250-m spatial resolution synthesized over 16 days from the period 2005–2019 to develop a rapid and effective method for monitoring disasters across a wide spatiotemporal range. Three types of disaster monitoring and scope extraction models are proposed: the normalized difference vegetation index (NDVI) median time standardization model (R_{NDVI_TM(i)}), the NDVI median phenology standardization model (R_{NDVI_AM(i)(j)}), and the NDVI median spatiotemporal standardization model (R_{NDVI_ZM(i)(j)}). The optimal disaster extraction threshold for each model in different time phases was determined using Otsu’s method, and the extraction results were verified by medium-resolution images and ground-measured data of the same or quasi-same period. Finally, the disaster scope of cultivated land in Heilongjiang Province from 2010–2019 was extracted, and the spatial and temporal patterns of the disasters were analyzed based on meteorological data. This analysis revealed that the three aforementioned models exhibited high disaster monitoring and range extraction capabilities, with verification accuracies of 97.46%, 96.90%, and 96.67% for R_{NDVI_TM(i)}, R_{NDVI_AM(i)}, and _(j)R_{NDVI_ZM(i)(j)}, respectively. The spatial and temporal disaster distributions were found to be consistent with the disasters of the insured plots and the meteorological data across the entire province. Moreover, different monitoring and extraction methods were used for different disasters, among which wind hazard and insect disasters often required a delay of 16 days prior to observation. Each model also displayed various sensitivities and were applicable to different disasters. Compared with other techniques, the proposed method is fast and easy to implement. This new approach can be applied to numerous types of disaster monitoring as well as large-scale agricultural disaster monitoring and can easily be applied to other research areas. This study presents a novel method for large-scale agricultural disaster monitoring.

Consistent data is seldom available for whole-catchment flood modelling in many developing regions, thus this study demonstrates how the complementary strengths of open and readily available geospatial datasets and tools can be leverage to map flood risk within acceptable levels of uncertainty for flood risk management. Available fragmented remotely-sensed and in situ datasets (including hydrological data, altimetry, digital elevation model, bathymetry, aerial photos, optical and radar imageries) are systematically integrated using 2-dimensional CAESAR-LISFLOOD model to quantify and recreate the extent and impact of the historic 2012 flood in Nigeria. Experimental modelling, calibration and validation is undertaken for the whole Niger-South hydrological catchment area of Nigeria, then segmented into sub-domains for re-validation to understand how data variability and uncertainties impact on the accuracy of model outcomes. Furthermore, aerial photos are applied for the first time in the study area for flood model validation and to understand how different physio-environmental properties influence synthetic aperture radar flood delineation capacity in the Niger Delta region of Nigeria.

Remote sensing has been used as an important tool for disaster monitoring and disaster scope extraction, especially for the analysis of temporal and spatial disasters patterns of large-scale long time series. In order to find out a rapid and effective method to monitor disaster in a wide range, based on the Google Earth Engine cloud platform, this study used MODIS vegetation index products of 250 meter spatial resolution synthesized in 16 days during the year 2005-2019 and three kinds of disaster monitoring and scope extraction models are proposed: normalized vegetation index median time standardization (RNDVI_TM(i)) model, the normalized vegetation index median phenology Standardization(RNDVI_AM(i)(j)) model, normalized vegetation index median spatiotemporal Standardization (RNDVI_ZM(i)(j)) model. The optimal threshold of disaster extraction for each model in different time phases was determined by Otsu method, and the extraction results were verified by Medium resolution image and ground measured data of the same or quasi-same period. Finally, the disaster scope of cultivated land in Heilongjiang province from 2010 to 2019 was extracted and the temporal and spatial pattern of disasters was analyzed based on the meteorological data. It shows that the three above-mentioned models have high disaster monitoring and range extraction capabilities with the verification accuracy of RNDVI_TM(i) 97.46%, RNDVI_AM(i)(j) 96.90%, and RNDVI_ZM(i)(j) 96.67% respectively. The spatial and temporal distribution of disasters is consistent with the disaster of the insured plots and meteorological data in the whole province. Meanwhile, it turns out that different monitoring and extraction methods are used in different disasters, among which wind hazard and insect disasters often need to be delayed for 16 days to observe. Each model also has various sensitivity and applicability to different disasters. Compared with other methods, this method is fast, and convenient, which allows it to be used for large-scale agricultural disaster monitoring and is easy to be applied into other research areas. The research provides a new idea for large-scale agricultural disaster monitoring.

As the amount of sensor data made available online increases, it becomes more difficult for users to identify useful datasets. Semantic web technologies improve discovery with meaningful ontologies, but the decision of suitability remains with the users. The GEO label provides a visual summary of the standardised metadata to aid users in this process. This work presents novel rules for deriving the information for the GEO label's multiple facets, such as user feedback or quality information, based on the Semantic Sensor Network Ontology and related ontologies. It enhances an existing implementation of the GEO label API to generate labels for resources of the Semantic Sensor Web. The prototype is deployed to serverless cloud infrastructures. We find that serverless GEO label generation is capable of handling two evaluation scenarios for concurrent users and burst generation. More real-world semantic sensor descriptions and an integration into large scale discovery platforms are needed to develop the presented solutions further.

The reliability analysis allows to estimate the system's probability of detecting and identifying outlier. Failure to identify an outlier can jeopardise the reliability level of a system. Due to its importance, outliers must be appropriately treated to ensure the normal operation of a system. The system models are usually developed from certain constraints. Constraints play a central role in model precision and validity. In this work, we present a detailed optical investigation of the effects of the hard and soft constraints on the reliability of a measurement system model. Hard constraints represent a case in which there exist known functional relations between the unknown model parameters, whereas the soft constraints are employed for the case where such functional relations can slightly be violated depending on their uncertainty. The results highlighted that the success rate of identifying an outlier for the case of hard constraints is larger than soft constraints. This suggested that hard constraints should be used in the stage of pre-processing data for the purpose of identifying and removing possible outlying measurements. After identifying and removing possible outliers, one should set up the soft constraints to propagate the uncertainties of the constraints during the data processing. This recommendation is valid for outlier detection and identification purpose.

As people grow a custom to effortless outdoor navigation there is a rising demand for similar possibility indoors as well. Unfortunately, indoor localization, being one of the necessary requirements for navigation, continues to be problem without a clear solution. In this article we are proposing a method for an indoor positioning system using a single image. This is made possible using small preprocessed database of images with known control points as the only preprocessing needed. Using feature detection with SIFT algorithm we can look through the database and find image which is the most similar to the image taken by user. Pair of images is then used to find coordinates of database image using PnP problem. Furthermore, projection and essential matrices are determined allowing for the user image localization ~ determining the position of the user in indoor environment. Benefits of this approach lies in the single image being the only input from user and no requirements for new onsite infrastructure and thus enables a simpler realization for the building management.

In this paper we evaluate the effects of hard and soft constraints on the Iterative Data Snooping (IDS), an iterative outlier elimination procedure. Here, the measurements of a levelling geodetic network were classified according to the local redundancy and maximum absolute correlation between the outlier test statistics, referred to as clusters. We highlight that the larger the relaxation of the constraints, the higher the sensitivity indicators MDB (Minimal Detectable Bias) and MIB (Minimal Identifiable Bias) for both the clustering of measurements and the clustering of constraints. There are circumstances that increase the family-wise error rate (FWE) of the test statistics, increase the performance of the IDS. Under a scenario of soft constraints, one should set out at least three soft constraints in order to identify an outlier in the constraints. In general, hard constraints should be used in the stage of pre-processing data for the purpose of identifying and removing possible outlying measurements. In that process, one should opt to set out the redundant hard constraints. After identifying and removing possible outliers, the soft constraints should be employed to propagate their uncertainties to the model parameters during the process of least-squares estimation.

Structural diversity is a key feature of forest ecosystems that influences ecosystem functions from local to macroscales. The ability to measure structural diversity in forests with varying ecological composition and management history can improve the understanding of linkages between forest structure and ecosystem functioning. Terrestrial LiDAR has often been used to provide a detailed characterization of structural diversity at local scales, but it is largely unknown whether these same structural features are detectable using aerial LiDAR data that are available across larger spatial scales. We used univariate and multivariate analyses to quantify cross-compatibility of structural diversity metrics from terrestrial versus aerial LiDAR in seven National Ecological Observatory Network sites across the eastern USA. We found strong univariate agreement between terrestrial and aerial LiDAR metrics of canopy height, openness, internal heterogeneity, and leaf area, but found marginal agreement between metrics that describe heterogeneity of the outer most layer of the canopy. Terrestrial and aerial LiDAR both demonstrated the ability to distinguish forest sites from structural diversity metrics in multivariate space, but terrestrial LiDAR was able to resolve finer-scale detail within sites. Our findings indicate that aerial LiDAR can be of use in quantifying broad-scale variation in structural diversity across macroscales.

North and West Africa are the most vulnerable regions to drought, due to the high variation in monthly precipitation. An accurate and efficient monitoring of drought is essential. In this study, we use TRMM data with remote sensing tools for effective monitoring of drought. The Drought Severity Index (DSI), Temperature Vegetation Drought Index (TVDI), Normalized Difference Vegetation Index (NDVI), and Normalized Vegetation Supply Water Index (NVSWI) are more useful for monitoring the drought over North and West Africa. To classify the areas affected by drought, we used the TRMM spatial maps to verify the TVDI, DSI and NVSWI indexes derived from MODIS. The DSI, TVDI, NVSWI and Monthly Precipitation Anomaly (NPA) indexes with the employ of MODIS-derived ET/PET and NDVI were chosen for monitoring the drought in the study area. The seasonal spatial correlation between the DSI, NPA, NVWSI, NDVI, TVDI and TCI indicates that NVSWI, NDVI and DSI present an excellent monitor of drought indexes. The change trend of drought from 2002 to 2018 was also characterized. The frequency of drought showed a decrease during this period.

Unmanned Aerial Vehicle (UAV) based remote sensing (RS) studies in glaciology are mainly focusing on obtaining accurate high-resolution data from UAV images. Studies for identifying and minimising the challenges faced during the UAV-based RS data acquisition survey on inaccessible and harsh terrains of mountain glaciers is limited. This study aims to examine the practical challenges faced during UAV surveys of glaciers and derive strategies to minimize them. To the authors' knowledge, this is the first study that addresses such problems over the Himalayan region. Here, the UAV surveys were conducted using a fixed-wing commercial-grade off-the-shelf UAV (eBee plus, SenseFly) on three glacier sites (East Rathong, Hamtah and Panchinala-A) located in different zones and climate regimes lying within the Indian part of Himalayas. From UAV collected photos, the study was able to generate ultra-high-resolution ortho-mosaicked images and Digital Elevation Models (DEMs) at 0.1m GSD. UAV-derived DEMs was able to achieve vertical (horizontal) accuracy of 0.45 and 0.21 m (0.15 and 0.1 m) with 3 and 6 ground control points (GCPs) for an area of 0.75 km2 and 1.38 km2. Accuracy assessment of UAV DEMs generated with and without GCPs indicate that GCPs are must to obtain decimetre level accurate DEM especially on glaciers with steep-valleyed terrains. The utility of the obtained ultra-high-resolution ortho-mosaicked images was demonstrated by generating glacier surface feature maps. Based on the challenges observed during UAV surveys, the study identifies and recommends best-suited locations on a glacier and its adjacent regions for conducting UAV surveys efficiently in the glaciated terrain of Himalayas and possibly beyond. Recommendations reported in this article shall minimise the challenges faced and involved risks for data acquisition and thus enable UAVs to cover more glaciated area successfully.

Interest in drone solutions in forestry applications is growing. Using drones, datasets can be captured flexibly and at high spatial and temporal resolutions when needed. In forestry applications, fundamental tasks include the detection of individual trees, tree species classification, bio-mass estimation, etc. Deep Neural Networks (DNN) have shown superior results when comparing with conventional machine learning methods such as Multi-Layer Perceptron (MLP) in cases of huge input data. The objective of this research was to investigate 3D convolutional neural networks (3D-CNN) to classify three major tree species in a boreal forest: pine, spruce, and birch. The proposed 3D-CNN models were employed to classify tree species in a test site in Finland. The classifiers were trained with a dataset of 3039 manually labelled trees. Then the accuracies were assessed by employing independent datasets of 803 records. To find the most efficient set of feature combination, we compare the performances of 3D-CNN models trained with hyperspectral (HS) channels, RGB channels, and canopy height model (CHM), separately and combined. It is demonstrated that the proposed 3D-CNN model with RGB and HS layers produces the highest classification accuracy. The producer accuracy of the best 3D-CNN classifier on the test dataset were 99.6%, 94.8%, and 97.4% for pines, spruces, and birches, respectively. The best 3D-CNN classifier produced ~5% better classification accuracy than the MLP with all layers. Our results suggest that the proposed method provides excellent classification results with acceptable performance metrics for HS datasets. Our results show that pine class was detectable in most layers. Spruce was most detectable in RGB data, while birch was most detectable in the HS layers. Furthermore, the RGB datasets provide acceptable results for many low-accuracy applications.

Global Positioning System (GPS) stations located along coastal areas have the ability to measure tide gauge (TG) records by reflected signal reception from the sea water surface. In this study we used the GPS signal-to-noise ratio (SNR) data from the SEPT station (44.63 ⁰N, 124.05 ⁰W) located at South Beach, Oregon, USA, to estimate the TG records using only a few measurements. First, we derived the TG record from a GPS station (GPS-TG) and used traditional TG data from the National Water Level Observation Network (NWLON) sentinel station (Station ID: 9435380) located in Oregon for validation purposes because it was closest to the SEPT station. Our results show that the GPS-TG and NWLON-TG correlate well with the correlation coefficient (CC) of 0.942 and the root mean square (RMS) of their residuals was about 12.90 cm. The corresponding TG prediction by autoregressive moving average (ARMA-TG) and singular spectrum analysis (SSA-TG) models are evaluated for their effectiveness over the station. The comparative analysis demonstrates that the GPS-TG has improved correlation with ARMA-TG (CC of ~0.981 CC, RMS of ~4.80 cm), and SSA-TG (CC of ~0.998 CC, RMS of ~ 0.88 cm) compared to the NWLON-TG (CC of ~0.942 CC, RMS of ~12.90 cm) values. We believe the outcomes from this study contribute to a better understanding of the numerical modeling of TG records as well as other measurements based on reflectometry techniques.

The pile foundation in the permafrost region is in a negative temperature environment, so the concrete is affected by the negative temperature of the surrounding soil.It not only affects the formation of concrete strength, but also leads to engineering quality accidents in serious cases.Based on the actual measurement of temperature at different strata depths and the comprehensive consideration of surface temperature, terrestrial heat flux and other parameters, the law curve of temperature change along depth in Greater Khingan is established.The calculated results of the curve are consistent with the measured results of ground temperature.The results show that the variation trend of ground temperature along the strata depth at different monitoring sites is basically the same. From June to November, the ground temperature at different depths tends to be constant.From December to May, the ground temperature at any depth within the depth range of 0 to 5.5m follows the law of the cosine function.Below 5.5m, the earth temperature no longer varies with depth.The research results can be used as reference for pile foundation construction under negative temperature environment.

Archaeological topography identification from high-resolution DEMs is a current method that is used with high success in archaeological prospecting of wide areas. I present a methodology trough which burial mounds (tumuli) from LiDAR DEMS can be identified. This methodology uses geomorphometric and statistical methods to identify with high accuracy burial mound candidates. Peaks, defined as local elevation maxima are found as a first step. In the second step, local convexity watershed segments and their seeds are compared with positions of local peaks and the peaks that correspond or have in vicinity local convexity segments seeds are selected. The local convexity segments that correspond to these selected peaks are further feed to a Random Forest algorithm together with shape descriptors and descriptive statistics of geomorphometric variables in order to build a model for the classification. Multiple approaches to tune and selected the proper training dataset, settings and variables were tested. The validation of the model was performed on the full dataset where the training was performed and on an external dataset in order to test the usability of the method for other areas in a similar geomorphological and archaeological setting. The validation was performed against manually mapped and field checked burial mounds from two neighbor study areas of 100 km² each. The results show that by training the Random Forest on a dataset composed of between 75% to 100% of the segments corresponding to burial mounds and ten times more non-burial mounds segments selected using latin hypercube sampling, 93% of the burial mound segments from the external dataset are identified. There are 42 false positive cases that need to be checked, and there are two burial mound segments missed. The method shows great promise to be used for burial mound detection on wider areas by delineating a certain number of tumuli mounds for model training.

Wildfire susceptibility maps display the wildfires occurrence probability, ranked from low to high, under a given environmental context. Current studies in this field often rely on expert knowledge, including or not statistical models allowing to assess the cause-effect correlation. Machine learning (ML) algorithms can perform very well and be more generalizable thanks to their capability of learning from and make predictions on data. Italy is highly affected by wildfires due to the high heterogeneity of the territory and to the predisposing meteorological conditions. The main objective of the present study is to elaborate a wildfire susceptibility map for Liguria region (Italy) by applying Random Forest, an ensemble ML algorithm based on decision trees. Susceptibility was assessed by evaluating the probability for an area to burn in the future considering where wildfires occurred in the past and which are the geo-environmental factors that favor their spread. Different models were compared, including or not the neighboring vegetation and using an increasing number of folds for the spatial-cross validation. Susceptibility maps for the two fire seasons were finally elaborated and validated and results critically discussed highlighting the capacity of the proposed approach to identify the efficiency of fire fighting activities.

Demographic and socio-economic developments couple with other requirements to satisfy human needs result in rapid urban expansion sometimes with increasing rate of surface extent greater than the rate of growth of population, that result in continuous sealing of ground surface thereby affecting ecosystem services. This study applied remote sensing toward achieving the progress of SDGs and stage to determine the ratio of the rate of land consumption to the rate of population growth of Gombe metropolis. QGIS 2.18 was used for the image processing and classification analysis for the key Landsat ETM+ (Enhanced Thematic Mapper), Impervious Surface Indices and population data to inform on the urban trend and LCR/PGR for the periods 2000-2005, 2005-2010, and 2010-2015. The result appears that the LCR/PGR for the periods of study show split trends. During 2000-2005 the result shows that the study area expanded outward with LCR/PGR of 1.2. The result also indicate that during 2005-2010, the study area densified with little expansion with the LCR/PGR of 0.8. The result further reveals that during 2010-2015 the LCR/PGR reached 1.8. Pointing that the study area expanded outward with the rate of ground sealing getting high.

This study investigates the spatial signatures of seasonal snow in Synthetic Aperture Radar (SAR) observations at different spatial scales and for different physiographic regions. Sentinel-1 C-band (SAR) backscattering coefficients (BSC) were analyzed in the Swiss Alps (SA), in high elevation forest and grasslands in Grand Mesa (GM), Colorado, and in North Dakota (ND) croplands. GM BSC exhibit 10dB sensitivity to wetness at small scales (~100 m) over homogeneous grassland. Sensitivity decreases to 5 dB in the presence of trees, and it is demonstrated that VH BSC sensitivity enables wet snow mapping below the tree-line. Area-variance scaling relationships show minima at ~100 m and 150-250 m respectively in barren and grasslands in SA and GM, increasing up to 1 km and longer in GM forests and ND agricultural fields. The spatial organization of BSC (as described by 1D-directional BSC wavelength spectra) exhibits multi-scaling behavior in the 100 -1,000 m range with a break at (180-360 m) that is also present in UAVSAR L-band measurements in GM. Spectral slopes in GM forested areas steepen during accumulation and flatten in the melting season with mirror behavior for grasslands reflecting changes in scattering mechanisms with snow depth and wetness, and vegetation mass and structure. Overall, this study reveals persistent patterns of SAR scattering variability spatially organized by land-cover, topography and regional winds with large inter-annual variability tied to precipitation. This dynamic scaling behavior emerges as an integral physical expression of snowpack variability that can be used to model sub-km scales and for downscaling applications.

Climate change impacts and adaptations is subject to ongoing issues that attract the attention of many researchers. Insight into the wind power potential in an area and its probable variation due to climate change impacts can provide useful information for energy policymakers and strategists for sustainable development and management of the energy. In this study, spatial variation of wind power density at the turbine hub-height and its variability under future climatic scenarios are taken under consideration. An ANFIS based post-processing technique was employed to match the power outputs of the regional climate model with those obtained from the reference data. The near-surface wind data obtained from a regional climate model are employed to investigate climate change impacts on the wind power resources in the Caspian Sea. Subsequent to converting near-surface wind speed to turbine hub-height speed and computation of wind power density, the results have been investigated to reveal mean annual power, seasonal, and monthly variability for a 20-year period in the present (1981-2000) and in the future (2081-2100). The results of this study revealed that climate change does not affect the wind climate over the study area, remarkably. However, a small decrease was projected for future simulation revealing a slightly decrease in mean annual wind power in the future compared to historical simulations. Moreover, the results demonstrated strong variation in wind power in terms of temporal and spatial distribution when winter and summer have the highest values of power. The findings of this study indicated that the middle and northern parts of the Caspian Sea are placed with the highest values of wind power. However, the results of the post-processing technique using adaptive neuro-fuzzy inference system (ANFIS) model showed that the real potential of the wind power in the area is lower than those of projected from the regional climate model.

Visualization of the seabed terrain is one of the key functions of marine geographic information system. The major challenge here is that the obtained data for seabed terrain usually consists of elevation and sediment only but does not include remote sensing images, which is important for ground terrain visualization, as they can be used as textures to reveal detailed information and achieve realistic visual results. Existing seabed terrain visualizing methods (including annotations, color-mapping, and texture-mapping) have limitations in reality and intuition, as they are inadequate to express the physical characteristics of sediment accurately. This paper presents a novel and advanced 3D visualization method of seabed terrain, which introduces the Physically Based Rendering (PBR) theory into the field of geographic visualization. We analyze the main categories of seabed sediments and their optical features respectively, then develop a procedural method to generate physically based rendering materials for different sediments. Then an enhanced bidirectional reflectance distribution function is employed to achieve accurate and intuitive terrain rendering. We also refine the texture sampling method and propose a procedural seabed objects generation method to construct a more natural and realistic undersea environment. Experimental results reveal our method can make good use of the limited seabed terrain data and get significant improvements in visual effect, which can help users to cognize and analyze the seabed geographic environment more accurately and intuitively.

Inland lake variations are considered sensitive indicators of global climate change. However, human activity is playing as a more and more important role in inland lake area variations. Therefore, it is critical to identify whether anthropogenic activity or natural event is playing as the dominant factor in inland lake surface area change. In this study, we proposed a Douglas-Peucker simplification algorithm and bend simplification algorithm combined method to locate major lake surface area disturbances; these disturbances were then characterized to extract the time series change features according to documented records; and the disturbances were finally classified into anthropogenic or natural. We took the nine lakes in Yunnan Province as test sites, a 31 years long (from 1987 to 2017) time series Landsat TM/OLI images and HJ-1A/1B used as data sources, the official records was used as references to aid the feature extraction and disturbance identification accuracy. Results of our method for both disturbance location and the disturbance identification could be concluded as follows: 1) The method can accurately locate the main lake changing events based on the time series lake surface area curve. The accuracy of this model for segmenting the lake area time series curves in our study area was 95.24%. 2) Our proposed method achieved an overall accuracy of 91.67%, with F-score of 94.67 for anthropogenic disturbances and F-score of 85.71 for natural disturbances. 3) According to our results, lakes in Yunnan Provence, China, have undergone extensive disturbances, and the human-induced disturbances occurred almost twice as often as natural disturbances, indicating intensified disturbances caused by human activities. This inland lake area disturbance identification method is expected to uncover whether a disturbance to inland lake area is human activity-induced or natural event.

Rice is an essential crop for national food security in Egypt. Increasing the population calls for regular increases in rice production. At the same time, cultivated rice crop areas should be decreased because of the gradual scarcity of irrigation water. This means more rice production should be gained from less rice area. This situation calls for the annual accurate system for rice monitoring and yield estimation. Therefore, it is necessary to apply a remotely sensed based system for rice cultivation assessment using satellite imagery parallel with field measurements of some biophysical parameters. Multi-temporal normalized difference vegetation index (NDVI) extracted from twelve sentinel-2 imagery cover the whole summer season with variance and maximum value assessed by ground control points (GCPs), were used to isolate uncultivated areas, then to isolate rice areas and other vegetation covers. object-based classification methods with kappa co-efficient 0.9261 and overall accuracy 94.92% was generated to discriminate rice crop area and other summer crops on the study area. Leaf area index (LAI) for the experiment the l site was calculated using the surface energy balance algorithm for Land (SEBAL) model and then tested versus measured (LAI). NDVI and LAI were used to generate an empirical ran rice yield prediction model. Then, this model was used to produce rice to yield a map. The study was carried out in an experimental site in Kafr Elsheikh governorate with a total area of 5040 Hectare. Produced cultivated land use map showed 95% overall accuracy. High similarity was observed between measured and calculated (LAI) with high accuracy of R2 = 0.94. of Rice, yield map showed expected to yield more to than a month before harvest. The generated yield map was tested using a correlation coefficient between actual yield and estimated yield with high accuracy R2 = 0.9. This method is applicable to estimate the acreage and productivity of rice in the northern Nile delta in adequate time before harvest.

Geodetic networks are essential for most geodetic, geodynamics and civil projects, such as monitoring the position and deformation of man-made structures, monitoring the crustal deformation of the Earth, establishing and maintaining a geospatial reference frame, mapping, civil engineering projects and so on. Before the installation of geodetic marks and gathering of survey data, geodetic networks need to be designed according to the pre-established quality criteria. In this study, we present a method for designing geodetic networks based on the concept of reliability. We highlight that the method discards the use of the observation vector of Gauss-Markov model. In fact, the only needs are the geometrical network configuration and the uncertainties of the observations. The aim of the proposed method is to find the optimum configuration of the geodetic control points so that the maximum influence of an outlier on the coordinates of the network is minimum. Here, the concept of Minimal Detectable Bias defines the size of the outlier and its propagation on the parameters is used to describe the external reliability. The proposed method is demonstrated by practical application of one simulated levelling network. We highlight that the method can be applied not only for geodetic network problems, but also in any branch of modern science.

Sentinel-1 SAR data preprocessing is essential for several earth observation applications, including land cover classification, change detection, vegetation monitoring, urban growth, natural hazards, etc. The information can be extracted from the 2x2 covariance matrix [C2] of Sentinel-1 dual-pol (VV-VH) acquisitions. To generate the covariance matrix from Sentinel-1 single look complex (SLC) data, several preprocessing steps are required. The ESA SNAP S-1 toolbox can be used to preprocess the data to generate a [C2] matrix. The polarimetric analysis in respective application fields often starts with the covariance matrix. However, due to limited availability of Sentinel-1 SLC data preprocessing workflow standards for polarimetric applications in contemporary research methods, downstream applications unable to comply with these workflows directly. In this paper, we propose a couple of generic practices to preprocess Sentinel-1 SLC data in SNAP S-1 toolbox, which would be beneficial for the radar remote sensing user community.

Over the last few decades, many countries, especially Caribbean island ones, have been challenged by the devastating consequences of natural disasters, which pose a significant threat to human health and safety. Timely information related to the distribution of vulnerable population and critical infrastructure are key for an effective disaster relief. OpenStreetMap (OSM) has repeatedly been shown to be highly suitable for disaster mapping and management. However, large portions of the world, including countries exposed to natural disasters, remain unmapped. In this study, we propose a methodology that relies on remotely sensed measurements (e.g. VIIRS, Sentinel-2 and Sentinel-1) and derived classification schemes (e.g. forest and built-up land cover) to predict the completeness of OSM building footprints in three small island states (Haiti, Dominica and St. Lucia). We find that the combinatorial effects of these predictors explain up to 94% of the variation of the completeness of OSM building footprints. Our study extends the existing literature by demonstrating how remotely sensed measurements could be leveraged to evaluate the completeness of OSM database, especially in countries at high risk of natural disasters. Identifying areas that lack coverage of OSM features could help prioritize mapping efforts, especially in areas vulnerable to natural hazards and where current data gaps pose an obstacle to timely and evidence-based disaster risk management actions.

Flood is identified as one of the major disasters in the world; it destroys both human and properties across the world, where lives are lost, properties, public infrastructure, farmlands and agricultural produce with farm crops carted away as a result of flood disaster. Studies revealed that the flood in itself is not the danger, but the level of human vulnerability to flooding disaster risk, which enhances its destructive capabilities. However, based on the challenges poses by flood disaster risk, this research identifies Ala river in Akure as a potential cause of flood, considering its location and other human activities around the river. Therefore, the research used Ala-river a case study to identify and mapped out areas susceptible to flood disaster risk. The research made use of both literature review and conducted goe-data gathering with the application GIS-computer database to retrieve georeferencing relevant data from the fieldwork in the study area of Ala-river basin to mapped out locations vulnerable to achieve the research aim. The research adopted a Geo-mapping of the vulnerable area to Ala-River basin using arc-GIS tool in combination with other software such as IKONAS and OLI (Operation Land Imager) for the production of the study area imagery, ER-ITERIM was used for the collection of rainfall data and FAO was applied for digital soil mapping. These applications produced; the land use/land cover map, digital elevation map, buffer map using 30 meters setback, annual rainfall map, soil types map, vulnerability map and soil textural table for the study area. Analysis of the produced and generated maps shows 316 buildings vulnerability to flood disaster risk; the soil texture and types, and alternative use to which the soil types can be useful. The research recommends that demolition of the identified 316 buildings prone to flood disaster and compliance of building construction to 30 meters setback by developers. Others are the conversion of the future land setback for urban agricultural purposes and preservation of water retention areas for agricultural activities during the dry season among others. The study concludes that relevant government agencies in the State and in particular in Akure South Local Government should ensure prompt compliance and implementation of the recommendations to avoid potential flood disaster risks.

In the last years, the development of low-cost GNSS sensors allowed monitoring in a continuous way movement related to natural processes like landslides with increasing accuracy and limited efforts. In this work, we present the first results of an experimental low-cost GNSS continuous monitoring applied to the unstable slope affecting the Madonna del Sasso Sanctuary (NW Italy). The courtyard of Sanctuary is built of two unstable blocks delimited by high cliff. Previous studies and non-continuous monitoring showed that blocks suffer a seasonal cycle of thermal expansion and long-term trend to downslope of few millimeters per year. The presence of continuous monitoring solution, could be an essential help to better understand the kinematics of unstable slope and to recognize the beginning of a possible paroxysm phase that could end with a failure of the unstable area. We tested the accuracy of the instruments and the first year of experimental measurements are presented. We also propose a methodological approach that considers the use of automatized procedures for the identification of anomalous trends and a risk communication strategy based on monitoring data.

Proper satellite-based crop monitoring applications at the farm-level often require near-daily imagery at medium to high spatial resolution. The synthesizing of ongoing satellite missions by ESA (Sentinel 2) and NASA (Landsat7/8) provides this unprecedented opportunity at a global scale; nonetheless, this is rarely implemented because these procedures are data demanding and computationally intensive. This study developed a complete stream processing in the Google Earth Engine cloud platform to generate harmonized surface reflectance images of Landsat7,8 and Sentinel 2 missions. The harmonized images were generated for two agriculture schemes in Bekaa (Lebanon) and Ninh Thuan (Vietnam) during the period 2018-2019. We evaluated the performance of several pre-processing steps needed for the harmonization including image co-registration, brdf correction, topographic correction, and band adjustment. This study found that the miss-registration between Landsat 8 and Sentinel 2 images, varied from 10 meters in Ninh Thuan, Vietnam to 32 meters in Bekaa, Lebanon, and if not treated, posed a great impact on the quality of the harmonized dataset. Analysis of a pair overlapped L8-S2 images over the Bekaa region showed that after the harmonization, all band-to-band spatial correlations were greatly improved from (0.57, 0.64, 0.67, 0.75, 0.76, 0.75, 0.79) to (0.87, 0.91, 0.92, 0.94, 0.97, 0.97, 0.96) in bands (blue, green, red, nir,swir1,swir2, ndvi) respectively. We demonstrated that dense observation of the harmonized dataset can be very helpful for characterizing cropland in highly dynamic areas. We detected unimodal, bimodal and trimodal shapes in the temporal NDVI patterns (likely cycles of paddy rice) in Ninh Thuan province only during the year 2018. We fitted the temporal signatures of the NDVI time series using harmonic (Fourier) analysis. Derived phase (angle from the starting point to the cycle's peak) and amplitude (the cycle's height) were combined with max-NDVI to generate an R-G-B image. This image highlighted croplands as colored pixels (high phase and amplitude) and other types of land as grey/dark pixels (low phase/amplitude). Generated harmonized datasets that contain surface reflectance images (bands blue, green, red, nir, swir1, swir2, and ndvi at 30 meters) over the two studied sites are provided for public usage and testing.

Vegetation in drylands is sensitive to climatic changes and human activities. Remote sensing and spatial analyses provide us useful tools for monitoring long-term vegetation dynamics over large regional scale. In this study, we analyzed the oasis vegetation cover change of the Tarim Basin using Landsat data sets from six epochs, 1975s, 1990s, 2000s, 2005s, 2010s and 2014. The results show that vegetation cover of oases increases from 34600 km² in 1975s to 101000 km² in 2014, though there was a vegetation coverage decrease from 77600 km² in 2000s to 42680 km² in 2010s. The percentage of annul water consumption has increased from 34% in 1970s to 52% in 2010s in the upper Tarim River, and decreased from 15% in 1970s to 9% in 2010s in the lower Tarim River. The decrease of oases area from 2000s to 2010s probably resulted from the rapid urbanization and large scale land reclamation. Although there is an increasing trend for oases coverage, local degradation of oases especially in the northern part occurred. This may be caused by inadequate water supply of the Tarim River. The results of multiple regression show that human activities contribute 70% of oases area change. Human induced water resources reallocation and heat energy balance is the primary cause of total oasis change.

Mapping settlement extents at the annual time step has a wide variety of applications in demography, public health, sustainable development, and many other fields. Recently, while more multitemporal urban feature or human settlement datasets have become available, issues still exist in remotely-sensed imagery due to coverage, adverse atmospheric conditions, and expenses involved in producing such feature sets. These challenges make it difficult to increase temporal coverage while maintaining high fidelity in the spatial resolution. Here we demonstrate an interpolative and flexible modeling framework for producing annual built-settlement extents. We use a combined technique of random forest and spatio-temporal dasymetric modeling with open source subnational data to produce annual 100m x 100m resolution binary settlement maps in four test countries of varying environmental and developmental contexts for test periods of five-year gaps. We find that in the majority of years, across all study areas, the model correctly identified between 85-99% of pixels that transition to built-settlement. Additionally, with few exceptions, the model substantially out performed a model that gave every pixel equal chance of transitioning to the category “built” in each year. This modelling framework shows strong promise for filling gaps in cross-sectional urban feature datasets derived from remotely-sensed imagery, provide a base upon which to create future built/settlement extent projections, and further explore the relationships between built area and population dynamics.

This study presents an analysis of subsidence rates and their effects on Mexico City. Mexico City is well known for its subsidence as a result of excess water withdrawal for many years. This study focuses on this problem utilizing the integration of Interferometric Synthetic Aperture Radar (InSAR), Continuous Global Positioning Systems (CGPS), and optical remote sensing data. Fifty-two ENVISAT-ASAR, nine GPS stations, and one Landsat ETM+ image from Mexico City area have been analyzed to prepare a better understanding of the subsidence rates and its effects on Mexico City’s commune. This study has utilized InSAR methods. It includes differential interferometry and Persistent Scatter Interferometry (PSI) to monitor the existing subsidence in the Mexico City area. The InSAR data covers the temporal baseline between 2002 until June 2010, and the GPS data include temporal baseline from 1998 until 2012. Maximum of 352 mm annually change in Line Of Sight (LOS) direction is in agreement with the previous geodetic studies. InSAR data have been compared with CGPS data at the same time interval. The finding of this study reveals a high amount of correlation (up to 0.98) between two independent geodetic methods. We also implemented the Support Vector Machine (SVM) analysis method based on Landsat ETM+ image to classify Mexico City’s populated density area. This method performed comparing the subsidence rates with populated area buildings. This integrated study shows that the fastest subsidence zone (i.e., areas greater than 100 mm/yr) in the over mentioned temporal baseline occurs in the high and sparsely populated areas

GRACE-derived Terrestrial Water Storage Anomalies (TWSA) continue to be used in an expanding array of studies to analyze numerous processes and phenomena related to terrestrial water storage dynamics, including groundwater depletions, lake storage variations, snow, and glacial mass changes, as well as floods, droughts, among others. So far, however, few studies have investigated how the factors that affect total water storage (e.g., precipitation, runoff, soil moisture, evapotranspiration) interact and combine over space and time to produce the mass variations that GRACE detects. This paper is an attempt to fill that gap and stimulate needed research in this area. Using the Nile River Basin as case study, it explicitly analyzes nine hydroclimatic and anthropogenic processes, as well as their relationship to TWS in different climatic zones in the Nile River Basin. The analytic method employed the trends in both the dependent and independent variables applying two geographically multiple regression (GMR) approaches: (i) an unweighted or ordinary least square regression (OLS) model in which the contributions of all variables to TWS variability are deemed equal at all locations; and (ii) a geographically weighted regression (GWR) which assigns a weight to each variable at different locations based on the occurrence of trend clusters, determined by Moran’s cluster index. In both cases, model efficacy was investigated using standard goodness of fit diagnostics. The OLS showed that trends in five variables (i.e., precipitation, runoff, surface water soil moisture, and population density) significantly (p<0.0001) explain the trends in TWSA for the basin at large. However, the models R2 value is only 0.14. In contrast, the GWR produced R2 values ranging between 0.40 and 0.89, with an average of 0.86 and normally distributed standard residuals. The models retained in the GWR differ by climatic zone. The results showed that all nine variables contribute significantly to the trend in TWS in the Tropical region; population density is an important contributor to TWSA variability in all zones; ET and Population density are the only significant variables in the semiarid zone. This type of information is critical for developing robust statistical models for reconstructing time series of proxy GRACE anomalies that predate the launch of the GRACE mission and for gap-filling between GRACE and GRACE-FO.

Farm activities continued sand winning operations and the allocation of plots of land to prospective developers in Ghana pose a serious threat to the forest covers and lifespan of the Forest and game reserves. With all the positive add ups to the country from forests, Ghana has lost more than 33.7%(equivalent to 2,500,000 hectares) of its forest, since the early 1990s between 2005 and 2010, the rate of deforestation in Ghana was estimated at 2.19% per annum; the sixth highest deforestation rate globally for that period. This shows how important forest monitoring can be to the forestry commission in Ghana. Despite the frameworks which have been developed to help Ghana to protect and restore its forest resources, inadequate monitoring systems remain a barrier to effective implementation. In this study, Google earth engine was used to map and analyze the structural changes of forest cover using JavaScript to query and compute Landsat, MODIS and NOAA AVHRR satellite imageries of the study area (Ghana) with spatial resolutions 30m, 250m and 7km respectively. A supervised classification was performed on three multi-temporal satellite imageries and a total of six major land use and land cover classes were identified and mapped. By using random Forest-classification technique, from 1985 to 2018 recorded by NOAA AVHRR, forest cover has decreased by 66% and 2000 to 2018 recorded by Landsat and MODIS 61% and 47% respectively. A decrease in the forest has been as a result of anthropogenic activities in Ghana. A change detection analysis was performed on these images and it was noted that Ghana is losing forest reserves in every 5years. Overlay of the reserved forest of the 2000 and the classified map of 2018 shows vegetation changed during 2000-2018 remarkably. Therefore, forest-related institutions like the Forestry Commission can employ and use this monitoring system on Google Earth Engine for processing satellite images particularly Landsat, MODIS and NOAA AVHRR for forest cover monitoring and analysis for fast, efficient and reliable results.

Recognition of the human interaction on the unconstrained videos taken from cameras and remote sensing platforms like a drone is a challenging problem. This study presents a method to resolve issues of motion blur, poor quality of videos, occlusions, the difference in body structure or size, and high computation or memory requirement. This study contributes to the improvement of recognition of human interaction during disasters such as an earthquake and flood utilizing drone videos for rescue and emergency management. We used Support Vector Machine (SVM) to classify the high-level and stationary features obtained from Convolutional Neural Network (CNN) in key-frames from videos. We extracted conceptual features by employing CNN to recognize objects from first and last images from a video. The proposed method demonstrated the context of a scene, which is significant in determining the behaviour of human in the videos. In this method, we do not require person detection, tracking, and many instances of images. The proposed method was tested for the University of Central Florida (UCF Sports Action), Olympic Sports videos. These videos were taken from the ground platform. Besides, camera drone video was captured from Southwest Jiaotong University (SWJTU) Sports Centre and incorporated to test the developed method in this study. This study accomplished an acceptable performance with an accuracy of 90.42%, which has indicated improvement of more than 4.92% as compared to the existing methods.

The technique of obtaining information or data about any feature or object from afar, called in technical parlance as remote sensing, has proven extremely useful in diverse fields. In the ecological sphere, especially, remote sensing has enabled collection of data or information about large swaths of areas or landscapes. Even then, in remote sensing the task of identifying and monitoring of different water reservoirs has proved a tough one. This is mainly because getting correct appraisals about the spread and boundaries of the area under study and the contours of any water surfaces lodged therein becomes a factor of utmost importance. Identification of water reservoirs is rendered even tougher because of presence of cloud in satellite images, which becomes the largest source of error in identification of water surfaces. To overcome this glitch, the method of the shape matching approach for analysis of cloudy images in reference to cloud-free images of water surfaces with the help of vector data processing, is recommended. It includes the database of water bodies in vector format, which is a complex polygon structure. This analysis highlights three steps: First, the creation of vector database for the analysis; second, simplification of multi-scale vector polygon features; and third, the matching of reference and target water bodies database within defined distance tolerance. This feature matching approach provides matching of one to many and many to many features. It also gives the corrected images that are free of clouds.

The Sentinel-2 satellite mission offers high resolution multispectral time series image data, enabling the production of detailed land cover maps globally. At this scale, the trade-off between processing time and result quality is a central design decision. Currently, this machine learning task is usually performed using pixelwise classification methods. The radical shift of the computer vision field away from hand engineered image features and towards more automation by representation learning comes with many promises, including higher quality results and less engineering effort. In this paper we assess fully convolutional neural networks architectures as replacements for a Random Forest classifier in an operational context for the production of high resolution land cover maps with Sentinel-2 time series at the country scale. Our contributions include a framework for working with Sentinel-2 L2A time series image data, an adaptation of the U-Net model for dealing with sparse annotation data while maintaining high resolution output, and an analysis of those results in the context of operational production of land cover maps.

In this case study, an active runway of a civilian airport in Zonguldak, Turkey, is used to assess the suitability of spaceborne digital elevation models (DEMs) to model an anthropogenic structure. The tested DEMs include the ASTER, the AW3D30 m, the SRTM-1”, the SRTM-3”, the SRTM-X, the TanDEM-3”, and the WorldDEMTM. A photogrammetric high accuracy DEM was also available for the tests. As a reference dataset, a line leveling survey of the runway using a Leica Sprinter 150/150M instrument was performed. The selection of a runway as a testbed for this type of investigation is justified by its unique characteristics, including its flat surface, homogenous surface material, and availability for a ground survey. These characteristics are significant because DEMs over similar structures are free from environment- and target-induced error sources. The most accurate DEM is the WorldDEMTM followed by the SRTM-3” and TanDEM-3”, with vertical errors (LE90) equal to 1.291 m, 1.542 m, and 1.56 m, respectively. This investigation uses a method for identifying the vertical errors in DEMs that is known as the runway method.

The quality of digital elevation models (DEMs) is inevitably affected by the limitations of the imaging modes and the generation methods. One effective way to solve this problem is to merge the available datasets through data fusion. In this paper, a fusion-based global DEM dataset (82°S-82°N) is introduced, which we refer to as GSDEM-30. This is a 30-m DEM mainly reconstructed from the unfilled SRTM1, AW3D30, and ASTER GDEM v2 datasets combining the multi-source and multi-scale fusion techniques. A comprehensive evaluation of the GSDEM-30 data, as well as the 30-m ASTER GDEM v2 and AW3D30 DEM, was presented. Global ICESat GLAS data and the local National Elevation Dataset (NED) were used as the reference for the vertical accuracy validation, while GlobeLand30 was introduced for the landscape analysis. Furthermore, we employed the maximum slope approach to detect the potential artefacts in the DEMs. The results show that the GDEM data are seriously affected by noise and artefacts. With the advantage of the multiple datasets and the refined post-processing, the GSDEM-30 are contaminated with fewer anomalies than both ASTER GDEM and AW3D30. The fusion techniques used can also be applied to the reconstruction of other fused DEM datasets.

Several methods have been tried to estimate air temperature using satellite imagery. In this paper, the results of two machine learning algorithms, Support Vector Machine and Random Forest, are compared with Multivariate Linear Regression, TVX and Ordinary kriging. Several geographic, remote sensing and time variables are used as predictors. The validation is carried out using four different statistics on a daily basis allowing the use of ANOVA to compare the results. The main conclusion is that Random Forest with residual kriging produces the best results (R$^2$=0.612 $\pm$ 0.019, NSE=0.578 $\pm$ 0.025, RMSE=1.068 $\pm$ 0.027, PBIAS=-0.172 $\pm$ 0.046), whereas TVX produces the least accurate results. The environmental conditions in the study area are not really suited to TVX, moreover this method only takes into account satellite data. On the other hand, regression methods (Support Vector Machine, Random Forest and Multivariate Linear Regression) use several parameters that are easily calculated from a Digital Elevation Model, adding very little difficulty to the use of satellite data alone. The most important variables in the Random Forest Model were satellite temperature, potential irradiation and cdayt, a cosine transformation of the julian day.

The objective to fast-track the mapping and registration of large numbers of unrecorded land rights globally, leads to the experimental application of Artificial Intelligence (AI) in the domain of land administration, and specifically the application of automated visual cognition techniques for cadastral mapping tasks. In this research, we applied and compared the ability of rule-based systems within Object Based Image Analysis (OBIA), as opposed to human analysis, to extract visible cadastral boundaries from Very high resolution (VHR) World View-2 image, in both rural and urban settings. From our experiments, machine-based techniques were able to automatically delineate a good proportion of rural parcels with explicit polygons where the correctness of the automatically extracted boundaries was 47.4% against 74.24% for humans and the completeness of 45% for machine, as against 70.4% for humans. On the contrary, in the urban area, automatic results were counterintuitive: even though urban plots and buildings are clearly marked with visible features such as fences, roads and tacitly perceptible to eyes, automation resulted in geometrically and topologically poorly structured data, that could neither be geometrically compared with human digitised, nor actual cadastral data from the field. These results provide an updated snapshot with regards to the performance of contemporary machine-drive feature extraction techniques compared to conventional manual digitising.

Precipitation downscaling is widely employed for enhancing the resolution and accuracy of precipitation products from general circulation models (GCMs). In this study, we propose a novel statistical downscaling method to foster GCMs’ precipitation prediction resolution and accuracy for monsoon region. We develop a deep neural network composed of convolution and Long Short Term Memory (LSTM) recurrent module to estimate precipitation based on well-resolved atmospheric dynamical fields. The proposed model is compared against GCM precipitation product and classical downscaling methods in the Xiangjiang River Basin in South China. Results show considerable improvement compared to the ECMWF-Interim reanalysis precipitation. Also, the model outperforms benchmark downscaling approaches, including 1) quantile mapping, 2) support vector machine, and 3) convolutional neural network. To test the robustness of the model and its applicability in practical forecast, we apply the trained network for precipitation prediction forced by retrospective forecasts from ECMWF model. Compared to ECMWF precipitation forecast, our model makes better use of the resolved dynamical field for more accurate precipitation prediction at lead time from 1 day up to 2 weeks. This superiority decreases along forecast lead time, as GCM’s skill in predicting atmospheric dynamics being diminished by the chaotic effect. At last, we build a distributed hydrological model and force it with different sources of precipitation inputs. Hydrological simulation forced with the neural network precipitation estimation shows significant advantage over simulation forced with the original ERA-Interim precipitation (with NSE value increases from 0.06 to 0.64), and the performance is just slightly worse than the observed precipitation forced simulation (NSE=0.82). This further proves the value of the proposed downscaling method, and suggests its potential for hydrological forecasts.

The detection of temporary color steel buildings and the analysis of their spatiotemporal patterns and characteristics are of great importance in many developing cities, because areas densely distributed with color steel buildings are usually problematic regions with high population density, sustainable development and risk level, which are challenging for local governments. Using the high-resolution satellite images of the Anning District, Lanzhou City, China as the case study, this paper first extracts the information of the color steel buildings for two different periods. Then the spatiotemporal distributional characteristics of the color steel buildings are analyzed, both at small scale (temporary buildings) and large scale (planthouses, warehouses, etc.). Finally, this paper utilizes kernel density estimation, Delaunay triangulation and Voronoi diagrams to examine the spatial distribution, aggregation and proximity characteristics of the color steel buildings. Our experiments show that the temporal and spatial differences of urban development, the phased characteristics and urban space are represented by various types of color steel buildings. Thus, there are some robust coupling relationship between color steel buildings, urban spatial form and urban development.

This paper focuses on the crucial role that remote sensing plays in divining land features. Data that is collected distantly provides information in spectral, spatial, temporal and radiometric domains, with each domain having the specific resolution to information collected. Diverse sectors such as hydrology, geology, agriculture, land cover mapping, forestry, urban development and planning, oceanography and others are known to use and rely on information that is gathered remotely from different sensors. In the present study, IRS LISS IV Multi-spectral data is used for land cover mapping. It is known, however, that the task of classifying high-resolution imagery of land cover through manual digitizing consumes time and is way too costly. Therefore, this paper proposes accomplishing classifications by way of enforcing algorithms in computers. These classifications fall in three classes: supervised, unsupervised, and object-based classification. In the case of supervised classification, two approaches are relied upon for land cover classification of high-resolution LISS-IV multispectral image. These approaches are Maximum Likelihood and Support Vector Machine (SVM). Finally, the paper proposes a step-by-step procedure for optical image classification methodology. This paper concludes that in optical data classification, SVM classification gives a better result than the ML classification technique.

If given the correct remotely sensed information, machine learning can accurately describe soil moisture conditions in a heterogeneous region at the large scale based on soil moisture readings at the small scale through rule transference across scale. This paper reviews an approach to increase soil moisture resolution over a sample region over Australia using the Soil Moisture Active Passive (SMAP) sensor and Landsat 8 only and a validation experiment using Sentinal-2 and the Advanced Microwave Scanning Radiometer (AMSR-E) over Nevada. This approach uses an inductive localized approach, replacing the need to obtain a deterministic model in favor of a learning model. This model is adaptable to heterogeneous conditions within a single scene unlike traditional polynomial fitting models and has fixed variables unlike most learning models. For the purposes of this analysis, the SMAP 36 km soil moisture product is considered fully valid and accurate. Landsat bands coinciding in collection date with a SMAP capture are down sampled to match the resolution of the SMAP product. A series of indices describing the Soil-Vegetation-Atmosphere Triangle (SVAT) relationship are then produced, including two novel variables, using the down sampled Landsat bands. These indices are then related to the local coincident SMAP values to identify a series of rules or trees to identify the local rules defining the relationship between soil moisture and the indices. The defined rules are then applied to the Landsat image in the native Landsat resolution to determine local soil moisture. Ground truth comparison is done via a series of grids using point soil moisture samples and air-borne L-band Multibeam Radiometer (PLMR) observations done under the SMAPEx-5 campaign. This paper uses a random forest due to its highly accurate learning against local ground truth data yet easily understandable rules. The predictive power of the inferred learning soil moisture algorithm did well with a mean absolute error of 0.054 over an airborne L-band retrieved surface over the same region.

Data creation is often the only way for researchers to produce basic geospatial information for the pursuit of more complex tasks and procedures such as those that lead to the production of new data for studies concerning river basins, slope morphodynamics, applied geomorphology and geology, urban and territorial planning, detailed studies, for example, in architecture and civil engineering, among others. This exercise results from a reflection where specific data processing tasks executed in Google Sketchup (Pro version, 2018) can be used in a context of interoperability with Geographical Information Systems (GIS) software. The focus is based on the production of contour lines and Digital Elevation Models (DEM) using an innovative sequence of tasks and procedures in both environments (GS and GIS). It starts in Google Sketchup (GS) graphic interface, with the selection of a satellite image referring to the study area—which can be anywhere on Earth's surface; subsequent processing steps lead to the production of elevation data at the selected scale and equidistance. This new data must be exported to GIS software in vector formats such as Autodesk Design Web format—DWG or Autodesk Drawing Exchange format—DXF. In this essay the option for the use of GIS Open Source Software (gvSIG and QGIS) was made. Correcting the original SHP by removing “data noise” that resulted from DXF file conversion permits the author to create new clean vector data in SHP format and, at a later stage, generate DEM data. This means that new elevation data becomes available, using simple but intuitive and interoperable procedures and techniques which confgures a costless work flow.

This article proposes the use of a multi-scale and multi-sensor approach to collect and modelling 3D data concerning wide and complex areas in order to obtain a variety of metric information in the same 3D archive, based on a single coordinate system. The employment of these 3D georeferenced products is multifaceted and the fusion or integration among different sensors data, scales and resolutions is promising and could be useful for the generation of a model that could be defined as hybrid. The correct geometry, accuracy, radiometry and weight of the data models are hereby evaluated comparing integrated processes and results from Terrestrial Laser Scanner (TLS), Mobile Mapping System (MMS), Unmanned Aerial Vehicle (UAV), terrestrial photogrammetry, using Total Station (TS) and Global Navigation Satellite System (GNSS) as topographic survey. The entire analysis underlines the potentiality of the integration and fusion of different solutions and is a crucial part of the “Torino 1911” project whose main purpose is mapping and virtually reconstructing the 1911 Great Exhibition settled in the Valentino Park in Turin (Italy).

One of the most destructive natural disasters is the earthquake which brings enormous risks to humankind. The objective of the current study was to determine the Earthquake’s remote sensing multiparameter (i.e. land surface temperature (LST), air temperature, specific humidity, precipitation and wind speed) spatiotemporal anomaly of many earthquake samples occurred during 2018 around the world. In this research 11 earthquake (M > 6:0) studied (4 samples selected in a land with transparent sky situations, 3 samples in land within cloudy situations and 4 samples in marine earthquakes). The interquartile range (IQR) and mean ± 2σ methods utilized to improve the efficiency of anomalous differences. As a result, based on the IQR method, negative anomaly before the event detected during the daytime in Mexico and during the nighttime in Afghanistan. In addition, a negative outlier of brightness temperature (BT) detected in Alaska before, after and during the event. In contrast, based on IQR and mean ± 2σ positive anomaly detected in precipitation before and after the event in all investigated examples. According to mean ± 2σ, negative anomaly LST, specific humidity, sea surface temperature (SST_100) and wind detected in most examined earthquake samples. In contrast, positive SST_0 anomaly observed in Fiji and Honduras after the earthquake. Our results suggested in marine earthquakes, for earthquake forecasting we can merge a prior negative anomaly in the wind speed and SST_100. Regarding the in land cloudy sky earthquakes, merging anomaly parameters could be the negative prior anomaly in BT, skin temperature, in contrast, a positive anomaly in precipitation. In land transparent sky earthquake, usually negative prior anomalies in air temperature, specific humidity and LST.

Delineating the cropping area of cocoa agroforests is a major challenge for quantifying the contribution of the land use expansion to tropical deforestation. Discriminating cocoa agroforests from tropical transition forests using multi-spectral optical images is difficult due to a similarity in the spectral characteristics of their canopy; moreover, optical sensors are largely impeded by the frequent cloud cover in the tropics. This study explores multi-season Sentinel-1 C-band SAR image to discriminate cocoa agroforests from transition forests for a heterogeneous landscape in central Cameroon. We use an ensemble classifier, random forest, to average SAR image texture features of GLCM (Grey Level Co-occurrence Matrix) across seasons; next, we compare classification performance with results from RapidEye optical data. Moreover, we assess the performance of GLCM texture feature extraction at four different grey level quantization: 32bits, 8bits, 6bits, and 4bits. The classification overall accuracy (OA) of texture-based maps outperformed that from an optical image; the highest OA of 88.8% was recorded at 6bits grey level. This quantization level, in comparison to the initial 32bits in SAR images, reduced the class prediction error by 2.9%. Although this prediction gain may be large for the landscape area, the resultant thematic map reveals the decrease and fragmentation of forest cover by cocoa agroforests. According to our classification validation, the Shannon entropy (H) or uncertainty provides a reliable validation for class predictions and reveals detail inference for discriminating inherently heterogeneous vegetation categories. The texture-based classification achieved a reliable accuracy considering the heterogeneity of the landscape and vegetation classes.

Due to the increasing importance of mangroves in climate change mitigation projects, more accurate and cost-effective aboveground biomass (AGB) monitoring methods are required. However, field measurement of AGB may be a challenge because of its remote location and the difficulty to walk in these areas. This study is based on the Livelihoods Fund’ Oceanium project of 10,000 hectare mangrove plantations monitoring. In a first step, the possibility of replacing traditional field measurements of sample plots in a young mangrove plantation by a semiautomatic processing of UAV-based photogrammetric point clouds was assessed. In a second step, Sentinel-1 radar and Sentinel-2 optical imagery were used as auxiliary information to estimate AGB and its variance for the entire study area under a model-assisted framework. AGB was measured using UAV imagery in a total of 95 sample plots. UAV plot data was used in combination with non-parametric Support Vector Regression (SVR) models for the estimation of the study area AGB using model-assisted estimators. Purely UAV-based AGB estimates and their associated standard error (SE) were compared with model-assisted estimates using (1) Sentinel-1, (2) Sentinel-2 and (3) a combination of Sentinel-1 and Sentinel-2 data as auxiliary information. The validation of the UAV-based individual tree height and crown diameter measurements showed a root mean square error (RMSE) of 0.21 m and 0.32 m respectively. Relative efficiency of the three model-assisted scenarios ranged between 1.61 and 2.15. Although all SVR models improved the efficiency of the monitoring over UAV-based estimates, the best results were achieved when a combination of Sentinel-1 and Sentinel-2 data was used. Results indicated that the methodology used in this research can provide accurate and cost-effective estimates of AGB in mangrove young plantations.

Yogyakarta City is one of the big city which is located in Java Island, Indonesia. Yogyakarta City, including study area (Pleret Sub District), are very prone to earthquake hazards, because close to several active earthquake sources. For example, Sunda Megathrust which often generates a big earthquake which can affect the study area. The Sunda Megathrust extends from north to south and west to east along the Sumatra and Java Islands. Furthermore, an active normal fault called as Opak Fault pass through right in the middle of Study area and divides the study area into east and west zone. Recently, after the devastating earthquake in 2006, the population of the study area increases significantly. As a result, the housing demand is also increasing. However, due to the absence of earthquake building code in the study area, locals tend to build improper new houses. Furthermore, in some part of the mountainous area in the study area, there are some building found in unstable slopes area. Due to this condition, the multi-hazard and risk study needs to be done in Pleret. The increasing of population and improper houses in Pleret Sub-District can lead to amplify the impact. Thus, the main objective of this study is to assess the multi-hazards and risk of earthquake and other related secondary hazards such as ground amplification, liquefaction, and coseismic landslide. The method mainly utilised the geographic information system, remote sensing and was fit up by the outcrop study. The results show that the middle part of the study area has a complex geological structure. It was indicated by a lot of unchartered faults was found in the outcrops. Furthermore, the relatively prone areas to earthquake can be determined. In term of the coseismic landslide, the prone area to the coseismic landslide is located in the east part of the study area in the middle slope of Baturagung Escarpment. The highly potential area of liquefaction is dominated in the central part of the study area. In term of building collapsed probability, the result shows that the safest house based on statistical analysis is the residential house with the building attribute of wood structure, roof cast material, distance more than 15 km from the earthquake source, and located above the Nglanggran Formation. Finally, the multi-hazard and risk analysis show that the middle part of the study area is more vulnerable than the other part of Pleret Sub-District.

Since the 1970s, land subsidence has been developing rapidly in the Beijing Plain, the systematic study of its evolution mechanism is of great significance to the sustainable development of the regional economy. First, based on ENVISAT ASAT and RADARSAT2 data, the land subsidence data in Beijing Plain were obtained using permanent interferometer technology. Second, based on the GIS platform and using fishing net tools, vector data of ground settlement with different resolutions were obtained. Through a series of tests, a scale of 960 metres was selected as the research unit, and the subsidence rate of the grid was obtained from 2004 to 2015. Finally, based on the Mann-Kendall mutation test method, a trend analysis of land subsidence changes in various grids was carried out. The results showed that single-year mutation mainly distributed in the middle and lower parts of the Yongding River alluvial fan and the Chaobai River alluvial fan, mainly occurring in 2015, 2005 and 2013, respectively. The upper and middle alluvial fan of the Chaobai River, the vicinity of the emergency water source and the edge velocity of the groundwater funnel have undergone several sudden changes. Combined with hydrogeology, basic geological conditions and the impact of the South-to-North Water transfer project, we analysed the causes of the mutations in the grid. The research results can provide a basis for the study and prevention of land subsidence in this area and help to further explore the trend characteristics of land subsidence in this area.

The continuous and precise mapping of glacier calving fronts is essential for monitoring and understanding rapid glacier changes in Antarctica and Greenland, which have the potential for significant sea level rise within the current century. This effort has been mostly restricted to the slow and painstaking manual digitalization of the calving front positions in thousands of satellite imagery products. Here, we have developed a machine learning toolkit to robustly and automatically detect glacier calving front margins in satellite imagery. The toolkit is based on semantic image segmentation using Convolutional Neural Networks (CNN) with a modified U-Net architecture to isolate the calving fronts from satellite images after having been trained with a dataset of images and their corresponding manually-determined calving fronts. As a case study we train our neural network on a varied set Landsat images with lowered resolutions from Jakobshavn, Sverdrup, and Kangerlussuaq glaciers, Greenland and test the results on novel images from Helheim glacier, Greenland to evaluate the performance of the approach. The neural network is able to identify the calving front in new images with a mean deviation of 96.3 m from the true fronts, equivalent to 1.97 pixels on average, while the corresponding error for manually-determined fronts on the same resolution images is 92.5 m. We find that the trained neural network significantly outperforms common edge detection techniques, and can be used to continuously map out calving-ice fronts with a variety of data products.

Uganda is one of the poorest nations in the world. To address the developmental challenges and understand social and economic status, it is important to obtain accurate data in a timely manner. Many studies have demonstrated that nighttime lights (NTL) can be used to measure human activities. Nevertheless, methods developed from these studies (1) suffer from coarse resolutions, (2) fail to capture the nonlinearity and multi-scale variability of geospatial data, and (3) perform poorly for agriculture-dependent regions. This study proposes a new enhanced light intensity model (ELIM) to estimate the Gross Domestic Product (GDP) at sub-national scales for Uganda. This model is developed by combining the NTL data from the Suomi National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (NPP-VIIRS), the population data from the Global Human Settlement Layer (GHSL), and information on agricultural production and market prices across several commodity types. This resulted in a gridded dataset for GDP and GDP per capita for Uganda at 1 km spatial resolution and district level to capture the spatial heterogeneity in economic activity.

Remote sensing data and techniques utilized for various purposes including natural disasters such as earthquake as well as flood. The research aims to consume liberates Landsat 8 images for investigating crashed airplanes such as MH370. Overall approximately 300 Landsat images with less than 10% clouds utilized in addition processed through Google Engine Platform. Due to the materials as well as the color of airplane body different from the area which is a plane crashed there, moreover, it should be the characteristics of the plane shapefile different in terms of albedo, temperature as well as vegetation index value. The research observed Landsat 8 data as well as methods utilized in this research, especially, NDVI, albedo in addition to band 4, capable to distinguish between the plane and its surrounding green area. Therefore, our result confirms during the research period, there was no plane on the location as well as MH370 not crashed in this site.

Due to the advantages of wide coverage and continuity, remotely sensed data are widely used for large-scale drought monitoring to compensate the deficiency and discontinuity of meteorological data. However, few researches have focused on the capability of various remotely sensed drought indices (RSDIs) for representing the spatio-temporal variations of the meteorological droughts. In this study, five RSDIs, namely Vegetation Condition Index (VCI), Temperature Condition Index (TCI), Vegetation Health Index (VHI), Modified Temperature Vegetation Dryness Index (MTVDI) and Normalized Vegetation Supply Water Index (NVSWI) were calculated using Moderate Resolution Imaging Spectroradiometer (MODIS) monthly NDVI and LST. The monthly NDVI and LST data were filtered by Savitzky-Golay (S-G) filtering method. Meteorological station-based drought index represented by Standardized Precipitation Evapotranspiration Index (SPEI) was compared with RSDIs. And the dimensionless Skill Score (SS) method was adopted to identify the spatiotemporally optimal RSDIs for presenting the meteorological droughts in the Yellow River basin (YRB) from 2000 to 2015. The results indicated that (1) RSDIs revealed a decreasing trend to the overall YRB consistent with SPEI except for in winter, and different variations of seasonal trends spatially; (2) the optimal RSDIs in spring, summer, autumn and winter were VHI, TCI, MTVDI and VCI, respectively, and the average correlation coefficient between the RSDIs and SPEI was 0.577 (=0.05); (3) different RSDIs have a 0–3 months’ time-lags compared with meteorological drought index.

Satellite altimeters have been used to monitor river and reservoir water levels, from which water storage estimates can be derived. Inland water altimetry can therefore play an important role in continental water resource management. Traditionally, satellite altimeters were designed to monitor homogeneous surfaces such as oceans or ice sheets, resulting in a poor performance over small inland water bodies due to the contribution from land contamination in the returned waveforms. The advent of synthetic aperture radar (SAR) altimetry (with its improved along-track spatial resolution) has enabled the measurement of inland water levels with a better accuracy and an increased spatial resolution. This paper presents three specialized algorithms or retrackers to retrieve water levels from SAR altimeter data over inland water bodies dedicated to minimizing land contamination from the waveforms. The performances of the proposed waveform portion selection method with three retrackers, namely, the threshold retracker, Offset Centre of Gravity (OCOG) retracker and 2-step physical-based retracker, are compared. Time series of water levels are retrieved for water bodies in the Ebro River basin (Spain). The results show good agreement with in situ measurements from the Ebro Reservoir (width is approximately 1.8 km) and Ribarroja Reservoir (width is approximately 400 m) with un-biased root-mean-square errors (RMSEs) of approximately 0.28 m and 0.16 m, respectively. The performances of all three retrackers are also compared with the European Space Agency’s ocean retracker in the Sentinel-3 Level-2 product.

The study focuses on the application of R programming language towards marine geological research with a case study of Mariana Trench. Due to its logical and straightforward syntax, multi-functional standard libraries, R is especially attractive to the geologists for the scientific computing. Using R libraries, the unevenness of various factors affecting Mariana Trench geomorphic structure has been studied. These include sediment thickness, slope steepness, angle aspect, depth at the basement and magmatism of the nearby areas. Methods includes using following R libraries: {ggplot2} for regression analysis, Kernel density curves, compositional charts; {ggalt} for Dumbbell charts for data comparison by tectonic plates, ranking dot plots for correlation analysis; {vcd} for mosaic plots, silhouette plots for compositional similarities among the bathymetric profiles, association plots; {car} for ANOVA. Bathymetric GIS data processing was dome in QGIS and LaTeX. The innovativeness of the work consists in the multi-disciplinary approach combining GIS analysis and statistical methods of R which contributes towards studies of ocean trenches, aimed at geospatial analysis of big data.

Operational weather and also flood forecasting has been performed successfully for decades and is of great socioeconomic importance. Up to now, forecast products focus on atmospheric variables, such as precipitation, air temperature and, in hydrology, on river discharge. Considering the full terrestrial system from groundwater across the land surface into the atmosphere, a number of important hydrologic variables are missing especially with regard to the shallow and deeper subsurface (e.g. groundwater), which are gaining considerable attention in the context of global change. In this study, we propose a terrestrial monitoring/forecasting system using the Terrestrial Systems Modeling Platform (TSMP) that predicts all essential states and fluxes of the terrestrial hydrologic and energy cycles from groundwater into the atmosphere. Closure of the terrestrial cycles provides a physically consistent picture of the terrestrial system in TSMP. TSMP has been implemented over a regional domain over North Rhine-Westphalia and a continental domain over European in a real-time forecast/monitoring workflow. Applying a real-time forecasting/monitoring workflow over both domains, experimental forecasts are being produced with different lead times since the beginning of 2016. Real-time forecast/monitoring products encompass all compartments of the terrestrial system including additional hydrologic variables, such as plant available soil water, groundwater table depth, and groundwater recharge and storage.

An oceanic radar altimeter such as TOPEX/Poseidon (T/P) is typically for observing elevation changes over the open oceans or large inland lakes/rivers, with limited applications over solid earth due to its large footprint and susceptibility to waveform contamination and slope effect. Here we demonstrate that it is possible to construct a long-term time series of glacier elevation change from T/P-series radar altimeters over two flat surfaces near a glacier terminus and an icefield (Sites A and B, with slopes of 2° and 0.8°) in Mt. Tanggula, Tibet, at elevations over 5400 m. We retracked radar waveforms using the subwaveform threshold algorithm, selected quality altimeter data (1/4 of the original) with nearly the same slope and adjusted the original elevations by fitting with a time-varying, 2nd order surface. The glacier elevation changes at the two sites from T/P (1993–2002) show seasonal elevation oscillations with linear rates at about −3 m/year and abnormal seasonal changes around the 1997–98 El Niño. Site A is over a deep valley in southern Tanggula. Its elevation dropped about 30 m over 1993–2002 (from T/P) and the glacier almost disappeared by 2016 (from altimeters and satellite images). Despite the sporadic Jason-2 and Jason-3 altimeter data, we also derived long-term rates of glacier elevation change over 1993–2017. Landsat-derived glacier area and elevation changes near the two sites confirm the rapid glacier thinning from the altimeters. The glacier meltwater near Site A supplied increasing source water to Chibuzhang Co west of Mt. Tanggula, contributing partially to its accelerated rising lake level. The glacier thinning at Site B (icefield) was correlated with the increased discharge of the Tuotuo River in eastern Mt. Tanggula, a source region of the Yangtze River. The successful detection of glacier thinning at the two sites shows that T/P-series altimeters can serve as a virtual station at a flat glacier spot to monitor long-term glacier elevation changes in connection to climate change. This virtual station concept is particularly useful for inaccessible glaciers, but its implementation faces two challenging issues: increasing the volume of quality altimeter data and improving the ranging accuracy over a targeted mountain glacier spot.

The rapid urbanization that takes place in developing regions such as Sub-Saharan Africa is associated with a large range of environmental and social issues. In this context, remote sensing is essential to provide accurate and up-to-date spatial information to support risk assessment and decision making. However, mapping urban areas remains a challenge because of their heterogeneity, especially in developing regions where the highest rates of misclassification are observed. Nevertheless, urban areas located in arid climates --- which are among the most vulnerables to anthropogenic impacts, suffer from the spectral confusion occurring between built-up and bare soil areas when using optical imagery. Today, the increasing availability of satellite imagery from multiple sensors allow to tackle the aforementioned issues by combining optical data with Synthetic Aperture Radar (SAR). In this paper, we assess the complementarity of the Landsat 8 and Sentinel-1 sensors to map built-up areas in twelve Sub-Saharan African urban areas, using a pixel-level supervised classification based on the Random Forest classifier. We make use of textural information extracted from SAR backscattering data in order to reduce the speckle noise and to introduce contextual information at the pixel level. Results suggest that combining both optical and SAR features consistently improves classification performances, mainly by enhancing the differentiation between built-up and bare lands. However, the fusion was less beneficial in mountainous case studies, suggesting that including features derived from a Digital Elevation Model (DEM) could improve the reliability of the proposed approach. As suggested by previous studies, combining features computed from both VV and VH polarizations consistently led to better classification performances. On the contrary, introducing textures computed from different spatial scales did not improve the classification performances.

Digital Elevation Models (DEMs) are widely used in geographic and environmental studies. In the current work, the fusion of multi-source DEMs is investigated to improve the overall accuracy of public domain DEMs. Multi-scale decomposition is an important analytical method in data fusion. Three multi-scale decomposition methods – the wavelet transform (WT), bidimensional empirical mode decomposition (BEMD) and nonlinear adaptive multi-scale decomposition (N-AMD) - are applied to the 1-arc-second Shuttle Radar Topography Mission Global digital elevation model (SRTM-1 DEM) and the Advanced Land Observing Satellite World 3D – 30 m digital surface model (AW3D30 DSM) in China. Of these, the WT and BEMD are popular image fusion methods. A new approach for DEM fusion is developed using N-AMD (which is originally invented to remove the cycle from sunspots). Subsequently, a window-based rule is proposed for the fusion of corresponding frequency components obtained by these methods. Quantitative results show that N-AMD is more suitable for multi-scale fusion of multi-source DEMs, taking the ice cloud and land elevation satellite (ICESat) global land surface altimetry data as a reference. The vertical accuracy of the fused DEM shows significant improvements of 29.6% and 19.3% in a mountainous region and 27.4% and 15.5% in a low-relief region, compared to the SRTM-1 and AW3D30 respectively. Furthermore, a slope position-based linear regression method is developed to calibrate the fused DEM for different slope position classes, by investigating the distribution of the fused DEM error with topography. The results indicate that the accuracy of the DEM calibrated by this method is improved by 16% and 13.6%, compared to the fused DEM in the mountainous region and low-relief region respectively, proving that it is a practical and simple means of further increasing the accuracy of the fused DEM.

Ziyuan-3 (ZY-3) satellite is the first civilian stereo mapping satellite in China and was designed to achieve the 1: 50000 scale mapping for land and ocean. Rigorous sensor model (RSM) is required to build the relationship between the three-dimensional (3D) object space and two-dimensional (2D) image space of ZY-3 satellite imagery. However, each satellite sensor has its own imaging system with different physical sensor models, which increase the difficulty of geometric integration of multi-source images with different sensor models. Therefore, it is critical to generate generic model， especially rational polynomial coefficients (RPCs) of optical imagery. Recently, relatively a few researches have been conducted on RPCs generation to ZY-3 satellite. This paper proposes an approach to evaluate the performance of RPCs generation from RSM of ZY-3 imagery. Three scenarios experiments with different terrain features (such as ocean, city and grassland) are designed and conducted to comprehensively evaluate the replacement accuracies of this approach and analyze the RPCs fitting error. All the experimental results demonstrate that the proposed method achieved the encouraging accuracy of better than 1.946E-04 pixel in both x-axis direction and y-axis direction, and it indicates that the RPCs is suitable for ZY-3 imagery and can be used as a replacement for the RSM of ZY-3 imagery.

Digital Terrain Models (DTMs) are currently a fundamental source of information in Earth Sciences. However, DTM-based studies can contain remarkable biases if limitations and inaccuracies of these models are disregarded. In this work, four freely available datasets such as SRTM C-SAR DEM, ASTER GDEM V2 and two airborne LiDAR derived DTMs (at 5 and 1 m spatial resolution, respectively) were analysed in a comparative study in three geomorphologically contrasted catchments located in Mediterranean geoecosystems under intensive human land use influence. Vertical accuracy as well as the influence of each dataset characteristics on hydrological and geomorphological modelling applicability were assessed by using classic geometric and morphometric parameters and the more recently proposed index of sediment connectivity. Overall vertical accuracy – expressed as Root Mean Squared Error (RMSE) and Normalized Median Deviation (NMAD) – revealed the highest accuracy in the cases of the 1 m (RMSE = 1.55 m; NMAD = 0.44 m) and 5 m LiDAR DTMs (RMSE = 1.73 m; NMAD = 0.84 m). Vertical accuracy of SRTM was lower (RMSE = 6.98 m; NMAD = 5.27 m) but considerably higher than in the case of ASTER (RMSE = 16.10 m; NMAD = 11.23 m). All datasets were affected by systematic distortions. As a consequence, propagation of these errors caused negative impacts on flow routing, stream network and catchment delineation and, to a lower extent, on the distribution of slope values. These limitations should be carefully considered when applying DTMs for hydrogeomorphological modelling.

Sea ice surface roughness affects ice-atmosphere interactions, serves as an indicator of ice age, shows patterns of ice convergence and divergence, affects the spatial extent of summer melt ponds, and ice albedo. We have developed a method for mapping sea ice surface roughness using angular reflectance data from the Multi-angle Imaging SpectroRadiometer (MISR) and lidar-derived roughness measurements from the Airborne Topographic Mapper (ATM). Using an empirical data modeling approach, we derived estimates of Arctic sea ice roughness ranging from centimeters to decimeters meters within the MISR 275-m pixel size. Using independent ATM data for validation, we find that histograms of lidar and multi-angular roughness values are nearly identical for areas with roughness <20 cm but that for rougher regions, the MISR-derived roughness has a narrower range of values than the ATM data. The algorithm is able to accurately identify areas that transition between smooth and rough ice. Because of its coarser spatial scale, MISR-derived roughness data have a variance of about half that ATM roughness data.

Seaweeds are regarded as one of the valuable coastal resources because of their usage in human food, cosmetics, and other industrial items. They also play a significant role in providing nourishment, shelter, and breeding grounds for fish and many other sea species. This study introduces a newly developed seaweed enhancing index (SEI) using spectral bands of near-infrared (NIR) and shortwave infrared (SWIR) of Landsat 8 satellite data. The seaweed patches in the coastal waters of Karachi, Pakistan were mapped using SEI, and its performance was compared with other commonly used indices - Normalized Difference Vegetation Index (NDVI) and Floating Algae Index (FAI). The accuracy of the mapping results obtained from SEI, NDVI, and FAI was checked with field verified seaweed locations. The purpose of the field surveys was to validate the results of this study and to evaluate the performance of SEI with NDVI and FAI. The performance of SEI was found better than NDVI and FAI in enhancing submerged patches of the seaweed pixels what other indices failed to do.

High-resolution laboratory-based thermal infrared spectroscopy is an up-and-coming tool in the field of geological remote sensing. Its spatial resolution allows for detailed analyses at centimeter to sub-millimeter scale. However, this increase in resolution creates challenges with sample characteristics such as grain size, surface roughness and porosity that can influence the spectral signature. This research explores the effect of rock sample surface preparation on the TIR spectral signatures. We applied three surface preparation methods (split, saw and polish) to determine how the resulting differences in surface roughness affects both the spectral shape as well as the spectral contrast. The selected samples are a pure quartz sandstone, a quartz sandstone containing a small percentage of kaolinite, and an intermediate-grained gabbro. To avoid instrument or measurement type biases we conducted measurements on three TIR instruments, resulting in directional hemispherical reflectance spectra, emissivity spectra and bi-directional reflectance images. Surface imaging and analyses were performed with scanning electron microscopy and profilometer measurements. We demonstrate that surface preparation affects the TIR spectral signatures influencing both the spectral contrast as well as the spectral shape. The results show that polished surfaces predominantly display a high spectral contrast while the sawed and split surfaces display up to 25% lower reflectance values. Furthermore, the sawed and split surfaces display spectral signature shape differences at specific wavelengths, which we link to mineral transmission features, surface orientation effects and multiple reflections in fine-grained minerals. Hence, the influence of rock surface preparation should be taken in consideration to avoid an inaccurate geological interpretation.