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.

Spectral unmixing represents both an application per se and a pre-processing step for several applications involving data acquired by imaging spectrometers. However, there is still a lack of publicly available reference data sets suitable for the validation and comparison of different spectral unmixing methods. In this paper we introduce the DLR HyperSpectral Unmixing (DLR HySU) benchmark dataset, acquired over German Aerospace Center (DLR) premises in Oberpfaffenhofen. The dataset includes airborne hyperspectral and RGB imagery of targets of different materials and sizes, complemented by simultaneous ground-based reflectance measurements. The DLR HySU benchmark allows a separate assessment of all spectral unmixing main steps: dimensionality estimation, endmember extraction (with and without pure pixe assumption), and abundance estimation. Results obtained with traditional algorithms for each of these steps are reported. To the best of our knowledge, this is the first time that real imaging spectrometer data with accurately measured targets are made available for hyperspectral unmixing experiments. The DLR HySU benchmark dataset is openly available online and the community is welcome to use it for spectral unmixing and other applications.

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.

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.

In daily travel and activities, pathfinding is a significant process. They are often used in transportation routes calculation. They have now evolved to be able to solve most situations of the pathfinding and its related problems. This review describes previous and recent studies on the pathfinding algorithms. It reviews the development of pathfinding algorithms in a classification base on their usage. The aim is to summarize the application of the pathfinding algorithms for the readers interested in the subject that can be used as a supplement.

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.

COVID-19 pandemic event requires a rapid response from various organizations at the international and national levels. One important response is the provision of information sharing facilities and monitoring of the spread of cases around the world, JHU CSSE developed the Dashboard in January 2020 and followed by WHO the same month for the WHO COVID-19 dashboard. Both dashboards have distributed information as expected by the user with their respective pros and cons. JHU CSSE Dashboard provides faster information with good access to mobile device users even though the display and color selection are less attractive. Information on the WHO COVID-19 Dashboard is often late but more data appearances and variations and comparisons between countries can be made. In the Indonesian context, ESSC for COVID-19 Geoportal as Esri Indonesia initiative has been developed with the support of data and information from various parties and developed with the principles of big data management which are fully supported by adequate spatial portal developer software from Esri. Particularly in Indonesia, there is not yet an adequate system to support spatial based decision making at the local level, therefore the development of a GIS dashboard to support provincial and district governments is highly recommended.

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.

We have an unprecedented ability to analyze and map the Earth’s surface, as deep learning technologies are applied to an abundance of Earth observation systems collecting images of the planet daily. In order to realize the potential of these data to improve conservation outcomes, simple, free, and effective methods are needed to enable a wide variety of stakeholders to derive actionable insights from these tools. In this paper we demonstrate simple methods and workflows using free, open computing resources to train well-studied convolutional neural networks and use these to delineate objects of interest in publicly available Earth observation images. With limited training datasets (<1000 observations), we used Google Earth Engine and Tensorflow to process Sentinel-2 and National Agricultural Imaging Program data, and use these to train U-Net and DeepLab models that delineate ground mounted solar arrays and parking lots in satellite imagery. The trained models achieved 81.5% intersection over union between predictions and ground-truth observations in validation images. These images were generated at different times and from different places from those upon which they were trained, indicating the ability of models to generalize outside of data on which they were trained. The two case studies we present illustrate how these methods can be used to inform and improve the development of renewable energy in a manner that is consistent with wildlife conservation.

The article describes the results of comparison of occurence of buildings (and address points) in Poland with delimitations of land use belonging to particular classes in the CORINE Land Cover (CLC) 2018 dataset. Large discrepancies have been identified, which reach on average approx. 34% of addresses and 35% of buildings located outside class 1 (artificial surfaces), mainly on terrains of class 2 (agricultural areas). Among single-family buildings it was 37% and among new addresses (forecasted or ‘under construction’ buildings) – as much as 50%. This puts a question mark over the possibility of using CLC data with resolution of 25 ha for monitoring of spatial planning and development in Poland for purposes of the diagnosis and assessment of the scale of dispersion of built-up areas. It is worth carrying out similar analyses in other countries, known for the deconcentration processes and a relatively large share of dispersed settlement e.g. other CEE countries, Spain, Portugal, Italy.

The world tsunami hazard map is developed to understand the pattern of occurrence of tsunami across the globe and for using this database for preparing effective evacuation plan and rescue operations. It has been developed by using Web GIS tools. HTML, CSS, Java Script, JQuery, IIS and Maps API are utilized to develop hazard map. These data can be even used by government authorities as well as by other nodal agencies.

India is one of the large sources of the anthropogenic pollutants and their increasing emission due to the recent economic growth in India. In this study we analyzed the annual and seasonal behaviors of ozone (O₃) gas using satellite remote sensing dataset from the sources Ozone Monitoring Instrument (OMI) over India region from 2006-2015. The study focuses on the seasonal behaviors of O₃ gas i.e., monthly, seasonal, annual mean variations of trace gas and also trend analysis of O₃ gas and comparison of the seasonal behavior of the ozone gas by trend analysis were assessed. In this study we also taken eleven cities to show the increment and decrement in four seasons of O₃ gas by taking 2006 as a base year and investigate the behaviors of gases during (2007-2015) years. Higher concentrations of O₃ south-to-north gradient, indicating the variations due to the impact of emissions and local meteorology. Ozone concentrations were higher during the warmer months. However, in winter season lowest concentration of O₃ seen due to the less amount of heat and due to cold days and ozone holes in the stratosphere. Instead, total O₃ concentrations rises over Delhi, Lucknow and Kolkata due to large population density, high traffic emission, highly polluted air and larger industrial activities.

This study evaluates spatial analyses of groundwater quality and environmental changes to obtain information on the groundwater contamination in the Permian Basin, Texas. Coupled with the U.S. government’s open data, these analyses can identify regions where environmental change could have potentially effected groundwater quality. A total of thirty-six wells were selected within the six counties: Andrews, Martin, Ector, Midland, Crane, and Upton. Spatial distribution maps were created for six different parameters: pH, total dissolved solids (TDS), chloride, fluoride, nitrate, and arsenic. Total groundwater quality maps incorporate all the contaminants and denote regions of poor, medium, and optimum conditions. To identify spatial changes in groundwater quality, maps were separated into two different time intervals, 1992-2005 and 2006-2019. We found that groundwater contamination resulted primarily from the mobilization of the contaminant from natural sources or anthropogenic activities such as chemical fertilizers. Overall, groundwater quality decreased during the study period from 1992 to 2019 as population and urban growth began to develop in the Permian Basin. This study contributes on understanding of the response of groundwater quality associated with environmental change in the Permian Basin. Therefore, this research provides important information for groundwater managements in developing plans for the use of water resource in the future for Texas.

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.

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.

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.

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.

Topographic mapping using stereo plotting is not effective because it takes much time and labor-intensive. Thus, this research was conducted to find the effective way to extract building footprint for mapping acceleration. Building extraction method in this process comprises four steps: ground / non-ground filtering, building classification, segmentation, and building extraction. Non-ground points from filtering process were classified as building with the algorithm based on multi-scale local dimensionality to separate points at the maximum separability plane. Segmentation using segment growing was used to separate each building, so edge detection could be conducted for each segment to create boundary of each building. Lastly, building extraction was conducted through three steps: edge points detection, building delineation, and building regularization. With 10 samples and step 0.5, classification resulted quality and miss factor of 0.597 and 0.524, respectively. The quality was improved by segmentation process to 0.604, while miss factor was getting worse to 0.561. Meanwhile, on average shape index value from extracted building had 0.02 difference and the number of errors was 30% for line segment comparison. Regarding positional accuracy using centroid accuracy assessment, this method could produce RMSE of 1.169 meters.

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.

Spatial planning processes generally consider three levels of planning, which are applied to three types of territory: state, county and city. As the coastal areas are of a significant natural, cultural, economic and social value, as well as are characterized by a diverse range of involved society with specific interests and needs, there is a necessity for an innovative and new approach to sustainable development planning in accordance with the modern age of growth, as well as to work with local communities in specific areas. Planning of a small populated area like village territory is more diverse and subject to the wishes and needs of the population. Small territory planning involves a very narrow circle of individuals or communities that identify spatial development needs for the future, including socio-economic, cultural, and environmental and climate change scenarios. In order to assess the development opportunities and needs of the area, it is necessary to monitor the area by regularly updating data. As it is well known, methodically derived data (facts) provide objectivity and transparency. Nowadays, when information about the present and the past is circulating very fast, it is possible to analyze the current situation, to forecast the future using databases, and to show several constructed realities (scenarios) using the geographic information system (GIS). Therefore, it is crucial to explore and find out the local needs-based planning approach to the development of village in coastal areas.

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.

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.

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.

Spatial Data Infrastructures (SDI) support the harvesting, curating, storage, and sharing of data along with providing access to development, analytic, and visualization tools that enable the building of innovative applications to address broad or specific challenges. SDIs can be especially powerful in bringing together data and tools supporting a particular theme – and this paper discusses and demonstrates the value of an SDI focused on Health. Many potential benefits of a Health SDI are proposed, and the case of supporting emergency response efforts is developed in detail. Leveraging a Health SDI, a Health Risk Index was created that provides emergency response personnel (both Emergency Operations Managers and Emergency Medical Responders) key insights into the unique health risks the impacted population faces due to the disaster. In order to establish the Health Risk Index, datasets from multiple national and global sources representing health data and social data that influences health outcomes – typically called social determinants of health – are harvested, merged, and republished to support further efforts at advancing the Health Risk Index. Visualizations of the Health Risk Index at the global, national, and sub-national levels down to the address level are presented along with demonstrations of its use.

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.

In this paper we analyse spatial variation in Japanese dialectal lexicon by assembling a set of methodologies using theories in variationist linguistics and GIScience, and tools used in historical GIS. Based on historical dialect atlas data, we calculate a linguistic distance matrix across survey localities. The linguistic variation expressed through this distance is contrasted with several measurements, based on spatial distance, utilised to estimate language contact potential across Japan, historically and at present. Further, administrative boundaries are tested for their separation effect. Measuring aggregate association within linguistic variation can contrast previous notions of dialect area formation by detecting continua. Depending on local geographies in spatial subsets, great circle distance, travel distance and travel times explain a similar proportion of the variance in linguistic distance despite the limitations of the latter two. While they explain the majority, two further measurements estimating contact have lower explanatory power: least cost paths modelling contact before the industrial revolution, based on DEM and seafaring, and a linguistic influence index based on settlement hierarchy. Historical domain boundaries and present day prefecture boundaries are found to have a statistically significant effect on dialectal variation. However, the interplay of boundaries and distance is yet to be identified. We claim that a similar methodology can address spatial variation in other digital humanities, given a similar spatial and attribute granularity.

We assemble a set of methodologies using theories in variationist linguistic and GIScience, and tools used in historical GIS to analyse spatial variation in Japanese dialectal lexicon. Based on historical dialect atlas data, we calculate a linguistic distance matrix across survey localities. The linguistic variation expressed through this distance is contrasted with several distance based measurements utilised to estimate the potential of language contact across Japan historically and at present. Besides, administrative boundaries are tested for their separation effect. Aggregate association measures within linguistic variation can challenge previous notions of dialect area formation. Depending on local geographies in spatial subsets, great circle distance, travel distance and travel times explain a similar proportion of the variance in dialectal variation despite the limitations of the latter two. While they explain the majority, two further contact estimations have lower explanatory power: least cost paths implemented to model contact before the industrial revolution, based on DEM and seafaring, and a linguistic influence index based on the law of gravity. Historical domain boundaries and present day prefecture boundaries are found to have a substantial effect on dialectal variation. However, the interplay of boundaries and distance is yet to be identified. We claim that a similar methodology can address spatial variation in other digital humanities, given a similar spatial and attribute granularity.

With the pressure of the increasing density of urban areas, some public infrastructures are moving to the underground to free up space above, such as utility lines, rail lines and roads. In the big data era, the three dimensional (3D) data can be beneficial to understand the complex urban area. Comparing to spatial data and information of the above ground, we lack of the precise and detailed information about underground infrastructures, such as the spatial information of underground infrastructure, the ownership of underground objects and the interdependence of infrastructures in the above and below ground. How to map reliable 3D underground utility networks and use it in the land administration? First, to explain the importance of this work and find a possible solution, this paper observes the current issues of the existing underground utility database in Singapore. A framework for utility data governance is proposed to manage the work process from the underground utility data capture to data usage. This is the backbone to support the coordination of different roles in the utility data governance and usage. Then, an initial design of the 3D underground utility data model is introduced to describe the 3D geometric and spatial information about underground utility data and connect it to the cadastral parcel for land administration. In the case study, the newly collected data from mobile Ground Penetrating Radar is integrated with the existing utility data for 3D modelling. It is expected to explore the integration of new collected 3D data, the existing 2D data and cadastral information for land administration of underground utilities.

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.

Conceived and developed by Christopher Alexander through his life’s work: The Nature of Order, wholeness is defined as a mathematical structure of physical space in our surroundings. Yet, there was no mathematics, as Alexander admitted then, that was powerful enough to capture his notion of wholeness. Recently, a mathematical model of wholeness, together with its topological representation, has been developed that is capable of addressing not only why a space is good, but also how much goodness the space has. This paper develops a structural perspective on goodness of space – both large- and small-scale – in order to bridge two basic concepts of space and place through the very concept of wholeness. The wholeness provides a de facto recursive definition of goodness of space from a holistic and organic point of view. A space is good, genuinely and objectively, if its adjacent spaces are good, the larger space to which it belongs is good, and what is contained in the space is also good. Eventually, goodness of space – sustainability of space – is considered a matter of fact rather than of opinion under the new view of space: space is neither lifeless nor neutral, but a living structure capable of being more living or less living, or more sustainable or less sustainable. Under the new view of space, geography or architecture will become part of complexity science, not only for understanding complexity, but also for making and remaking complex or living structures.

With the growth of air traffic demand in busy airspace, there is an urgent need for airspace sectorization to increase air traffic throughput and ease the pressure on controllers. The purpose of this paper is to develop a method framework that can perform airspace sectorization automatically, reasonably, which can be used as an advisory tool for controllers as an automatic system, especially for eliminating irregular sector shapes generated by simulated annealing algorithm (SAA) based on region growth method. Two graph cutting method, dynamic Monte Carlo method by changing location of flexible vertices (MC-CLFV) and Monte Carlo method by radius changing (MC-RC) were developed to eliminating irregular sector shapes generated by SAA in post-processing. The experimental results show that the proposed method framework of AS can automatically and reasonably generate sector design schemes that meet the design criteria. Our methodology framework and software can provide assistant design and analysis tools for airspace planners to design airspace, improve the reliability and efficiency of airspace design, and reduce the burden of airspace planners. In addition, this lays the foundation for reconstructing airspace with more intelligent method.

The problem of identifying functional regions in an urban setting has been approached in literature using two general methodologies: top-down, encoding expert knowledge on urban planning and design (e.g. into patterns) and using that knowledge for identification, and bottom-up, relying on crowdsourcing and Volunteered Geographic Information (VGI) to train learning models, using techniques such as Latent Dirichlet Allocation (LDA) topic modeling. Both approaches have their advantages but also face important limitations, with knowledge-based approaches being criticized for scalability and transferability issues and data-driven approaches for lacking interpretability and depending heavily on data quality. To mitigate these disadvantages, we propose a novel framework that fuses data and knowledge in three different ways: functional regions identified from individual approaches are evaluated against each other, knowledge from patterns is used to adjust learning model results and topic models are used to adjust pattern-based results. The proposed methodologies are demonstrated through the use case of identifying shopping-related functional regions in the Los Angeles metropolitan area. Results show that the combination of results from knowledge-based and data-driven techniques can help uncover discrepancies between the two different approaches and smoothen inaccuracies caused by the limitations of each approach.

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.

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.

Spatio-temporal Land-Use and Land-Cover (LULC) changes have been affecting geo-environmental and climate change globally. This study aims to analyze LULC changes in Bahir Dar city and its surrounds. Landsat 5 TM (1987), Landsat 7 ETM+ (2002) and Landsat 8 OLI (2017) and SPOT images, and aerial photographs, master plan map and Google Earth Landsat images were used to analyze changes. In Bahir Dar city and its surrounds, LULC has been changing in space and time. During 1987-2017, more than 50% of the study area was covered with cropland. Settlement areas have increased from 3.3% in 1987 to 9.13% in 2017. However, wetland vegetation, shrubland, grassland, forest, and waterbodies have degraded. These changes are mainly attributed to population growth and its effect on the environment. Land-use and land-cover is a serious problem and it causes land and environmental degradation, climate change and loss of the biological environment.

Authorities define cities – or human settlements in general – through imposing top-down rules in terms of whether buildings belong to cities. Emerging geospatial big data makes it possible to define cities from the bottom up, i.e., buildings determine themselves whether they belong to a city based on the notion of natural cities that is defined based on head/tail breaks, a classification and visualization tool for data with a heavy-tailed distribution. In this paper, we used 125 million building locations – all building footprints of America (mainland) or their centroids more precisely – to derive 2.1 million natural cities in the country (http://lifegis.hig.se/uscities/). These natural cities – in contrast to government defined city boundaries – constitute a valuable data source for city-related research.

Cities around the world face an increasing need for land as density in urban areas increases rapidly. The pressure to expand a city's space is especially acute for a city-state like Singapore. In the big data era, a data-driven approach of underground spaces is necessary for the sustainable development of a city along with rapid urbanization. A reliable three dimensional (3D) digital map of utility networks is crucial for urban planners to understand one of the most impactful aspects of underground space planning. How to map reliable 3D underground utility networks and use it in the land administration? This is a challenging issue, especially for cities with limited land resources, congested underground spaces, and a lack of uniform existing practices. First, this paper proposes a framework for utility data governance from the underground utility data survey to data usage. This is the backbone to support coordination of different roles in the utility data management and usage. Then, an initial design of the 3D utility cadastral data model is introduced, which aims to support the 3D modelling of utility networks and connect it to the cadastral parcel. It is expected that reliable and accurate information on underground utility networks can lead to a better understanding and management of underground space, which eventually contributes to better city planning, making the unseen structures visible.

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.

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.

The authors introduce the GRASS GIS add-on module g.citation as an initial implementation of a fine-grained software citation concept. The module extends the existing citation capabilities of GRASS GIS, which until now only provide for automated citation of the software project as a whole, authored by the GRASS Development Team, without reference to individual persons. The functionalities of the new module enable individual code citation for each of the over 500 implemented functionalities, including add-on modules. Three different classes of citation output are provided in a variety human- and machine-readable formats. The implications of this reference implementation of scientific software citation for both for the GRASS GIS project and the OSGeo foundation are outlined.

Indoor scene recognition has great significance for intelligent applications such as mobile robots, location-based services (LBS) and so on. Wherever we are or whatever we do, we are under a specific scene. The human brain can easily discern a scene with a quick glance. However, for a machine to achieve this purpose, on one hand, it often requires plenty of well-annotated data which is time-consuming and labor-intensive. On the other hand, it is hard to learn effective visual representations due to large intra-category variation and inter-categories similarity of indoor scenes. To solve these problems, in this paper, we adopted an unsupervised visual representation learning method which can learn from unlabeled data with a Siamese Convolutional Neural Network (Siamese ConvNet) and graph-based constraints. Specifically, we first mined relationships between unlabeled samples with a graph structure. And then, these relationships can be used as supervision for representation learning with a Siamese network. In this method, firstly, a k-NN graph would be constructed by taking each image as a node in the graph and its k nearest neighbors are linked to form the edges. Then, with this graph, cycle consistency and geodesic distance would be considered as criteria for positive and negative pairs mining respectively. In other words, by detecting cycles in the graph, images with large differences but in the same cycle can be considered as same category (positive pairs). By computing geodesic distance instead of Euclidean distance from one node to another, two nodes with large geodesic distance can be regarded as in different categories (negative pairs). After that, visual representations of indoor scenes can be learned by a Siamese network in an unsupervised manner with the mined pairs as inputs. In order to evaluate the proposed method, we tested it on two scene-centric datasets, MIT67 and Places365. Experiments with different number of categories have been conducted to excavate the potential of proposed method. The results demonstrated that semantic visual representations for indoor scenes can be learned in this unsupervised manner. In addition, with the learned visual representations, indoor scene recognition models trained with the learned representations and a few of labeled samples can achieve competitive performance compared to the state-of-the-art approaches.

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.

The General Data Protection Regulation (GDPR) protects the personal data of natural persons and at the same time allows the free movement of such data within the European Union (EU). Hailed as majestic by admirers and dismissed as protectionist by critics, the Regulation is expected to have a profound impact around the world, including in the African Union (AU). For European–African consortia conducting research that may affect the privacy of African citizens, the question is ‘how to protect personal data of data subjects while at the same time ensuring a just distribution of the benefits of a global digital ecosystem?’ We use location privacy as a point of departure, because information about an individual’s location is different from other kinds of personally identifiable information. We analyse privacy at two levels, individual and cultural. Our perspective is interdisciplinary: we draw from computer science to describe three scenarios of transformation of volunteered/observed information to inferred information about a natural person and from cultural theory to distinguish four privacy cultures emerging within the EU in the wake of GDPR. We highlight recent data protection legislation in the AU and discuss factors that may accelerate or inhibit the alignment of data protection legislation in the AU with the GDPR.

Spatial modeling in Geographic Information Systems (GIS) always involves choices. The existence of constraints, either of a financial nature or related to the specifics of the software itself, to the algorithms, the uncertainty and even the reliability of the data, the purposes and the applications of the studies, make this a kind of guiding compass for GIS analysts. Building on a previous exercise of data acquisition (check-ins) based on two Digital Social Networks (DSN — Facebook and Foursquare) and on the awareness of the use of Volunteered Geographic Information (VGI) generated by tourists sharing their topophilic ties through DSN, the present analysis aims to evaluate the contribution of modern techniques of spatial analysis applied to tourism in the “Alta and University of Coimbra” area. Concepts and procedural tasks related to density determination, cluster analysis and identification of patterns associated with regionalized variables have thus been implemented with the purpose of evaluating and comparing the results obtained through the application of two techniques of spatial analysis, Kernel Density Estimation (KDE) and Optimized Hot-Spot Analysis (OHSA) & Inverse Distance Weighting (IDW) Interpolation.

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).

Location-Based Services (LBS) have been widely deployed for the connected vehicle (CV) applications such as vehicle navigation,vehicle tracking and location-based augmented reality. The current LBS deployments have limitations in supporting time-critical CV use cases, including vehicle to vehicle (V2V), vehicle to infrastructure (V2I) and vehicle-to-people (V2P) safety applications. The paper presents the new LBS framework based on the publish-subscribe communication paradigm, to enable device-to-device (D2D) connections through use of selected application protocols in the application layer of the TCP/IP layered protocol model. Two publish-subscribe application protocols, Distributed Data Service (DDS) real-time publish and subscribe (DDS-RTPS) and Message Queue Telemetry Transport (MQTT), are introduced to support the LBS D2D applications. A number of test scenarios with Mosquitto MQTT and OpenDDS under 4G-mobile broadband (MBB) services are designed to assess the transmit/receive round-trip time (RTT) and packet-loss rate (PLR) with settings of a publisher to multiple subscribers, to simulate the connections to multiple vehicles. The transmission frequency is set for 10 Hz and the message sizes vary from 100 to 2000 Bytes. The PLRs are defined as the percentages of the delayed messages beyond a delay limit. Static test results with OpenDDS show that for the RTT delay beyond the limit of 100 ms, the total PLRs range between 5.25% and 8.76% for the message size of 50 to 2000 Bytes. Vehicle testing results with Mosquitto show that PLRs for the RTT delays between 200 ms and 1000 ms are 0.63%, 3.58% and 5.77%, for connections with 1, 4 and 10 vehicles, respectively. The results demonstrate the potential of the D2D LBS framework for medium-demanding CV safety applications such as V2P and V2I use cases, taking advantages of the 4G-MBB services and 5G extreme mobile broadband (eMBB) services and mobile devices generally available with all road users.

As noted in the introductory quotation, an ideal map was long ago seen as the map of the map, the map of the map, of the map, and so on endlessly. This recursive perspective on maps, however, has received little attention in cartography. Cartography, as a scientific discipline, is essentially founded on Euclidean geometry and Gaussian statistics, which deal with respectively regular shapes, and more or less similar things. It is commonly accepted that geographic features – such as rivers, cities, streets and building – are not regular and that the Earth’s surface is full of fractal or scaling or living phenomena with far more small things than large ones at different levels of scale. This paper argues for a new paradigm in mapping, based on fractal or living geometry and Paretian statistics, and – more critically – on the new conception of space, conceived and developed by Christopher Alexander, that space is neither lifeless nor neutral, but a living structure capable of being more living or less living. The fractal geometry is not limited to Benoit Mandelbrot’s framework, but is extended towards Christopher Alexander’s living geometry and based upon the third definition of fractal: A set or pattern is fractal if the scaling of far more small things than large ones recurs multiple times. Paretian statistics deals with far more small things than large ones, so it differs fundamentally from Gaussian statistics, which deals with more or less similar things. Under the new paradigm, I make several claims about maps and mapping: (1) Topology of geometrically coherent things – in addition to that of geometric primitives – enables us to see a scaling or fractal or living structure; (2) Under the third definition, all geographic features are fractal or living, given the right perspective and scope; (3) Exactitude is not truth – to paraphrase Henri Matisse – but the living structure is; and (4) Töpfer’s law is not universal, but scaling law is. All these assertions are supported by evidence, drawn from a series of previous studies. This paper demands a monumental shift in perspective and thinking from what we are used to on the legacy of cartography and GIS.

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.

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.

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.

Coal combustion in thermal power plants releases ash. Ash is reported to cause different adverse health hazards in humans and other organisms. Owing to the presence of radionuclides, it is also considered as a potential radiation hazard. In this study, based on the surface radiation measurements and relevant ancillary data, expected radiation risk zones were identified with regard to the human population residing near the Thermal Power Plant. With population density as the risk determining criteria, about 20% of the study area was at ‘High’ risk and another 20% of the study area was at ‘Low’ risk zone. The remaining 60% was under medium risk zone. Based on the findings remedial measures which may be adopted have been suggested.

There is much discussion regarding the Sustainable Development Goals’ (SDGs) capacity to promote inclusive development. While some argue that they represent an opportunity for goal-led alignment of stakeholders and evidence-based decision-making, other voices express concerns as they perceive them as a techno-managerial framework that measures development according to quantitatively defined parameters and does not allow for local variation. We argue that the extent to which the positive or negative aspects of the SDGs prevail depends on the monitoring system’s ability to account for multiple and intersecting inequalities. The need for sub-nationally (urban) representative indicators poses an additional methodological challenge – especially in cities with intra-urban inequalities related to socio-spatial variations across neighbourhoods. This paper investigates the extent to which the SDG indicators’ representativeness could be affected by inequalities. It does so by proposing a conceptual framing for understanding the relation between inequalities and SDG monitoring, which is then applied to analyse the current methodological proposals for the indicator framework of the “urban SDG”, Goal 11. The outcome is a call for 1) a more explicit attention to intra-urban inequalities, 2) the development of a methodological approach to “recalibrate” the city-level indicators to account for the degree of intra-urban inequalities, and 3) an alignment between methodologies and data practices applied for monitoring SDG 11 and the extent of the underlying inequalities within the city. This would enable an informed decision regarding the trade-off in indicator representativeness between conventional data sources, such as censuses and household surveys, and emerging methods, such as participatory geospatial methods and citizen-generated data practices.

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.