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.

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.

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

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

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

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

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

Studying hydro-meteorological factors variations and its links to large-scale atmospheric circulation systems can facilitate the understanding of the hydrological processes and sustainable water resources management in the source region of the Yangtze River (SRYR). Currently, researches mostly focused on the temporal and spatial variation characteristics in hydro-meteorological factors; however, researches on the hydro-meteorological variations and its links to large-scale atmospheric circulation systems in the SRYR are scarce. Based on long-term hydro-meteorological and reanalysis data, this research investigated multiscale variations of hydro-meteorological factors and its links to large-scale atmospheric circulation characteristic indices during 1957~2012 in the SRYR. The results showed that the amounts of streamflows and precipitation in the SRYR declined during the 1990s. Since the 2000s, the amounts of streamflows and precipitation had increased significantly climate in the SRYR. The change trends of precipitation and streamflows in the SRYR are synergetic at annual and seasonal scales, and have three significant periods, namely 3~5 years, 15–20 years and 30–40 years. The South Asia monsoon (SAM) plays a relatively more important role in the hydro-meteorological factors changes in the SRYR. The relative contributions of SAM to streamflows and precipitation changes were 83.6% and 78%, respectively. During the driest (wettest) year, the SAM is relatively weak (strong), and brings less (more) southwest airflow into the SRYR, less (more) precipitation and streamflows will be generated in the SRYR.

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

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

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

Recent advances in open data initiatives allow us to free access to a vast amount of open LiDAR data in many cities. However, most of these open LiDAR data over cities are acquired by airborne scanning, where the points on façades are sparse or even completely missing due to the viewpoint and object occlusions in the urban environment. Integrating other sources of data, such as ground images, to complete the missing parts is an effective and practical solution. This paper presents an approach for improving open LiDAR data coverage on building façades by using point cloud generated from ground images. A coarse-to-fine strategy is proposed to fuse these two different sources of data. Firstly, the façade point cloud generated from terrestrial images is initially geolocated by matching the SFM camera positions to their GPS meta-information. Next, an improved Coherent Point Drift algorithm with normal consistency is proposed to accurately align building façades to open LiDAR data. The significance of the work resides in the use of 2D overlapping points on the outline of buildings instead of limited 3D overlap between the two point clouds and the achievement to a reliable and precise registration under possible incomplete coverage and ambiguous correspondence. Experiments show that the proposed approach can significantly improve the façades details of buildings in open LiDAR data and improving registration accuracy from up to 10 meters to less than half a meter compared to classic registration methods.

Rapid urbanization has dramatically spurred the economic development over the past three decades, especially in China, but has nevertheless had negative impacts on natural resources since it is an irreversible process. Thus, it is essential to timely monitor and quantitatively analysis the changes in land use over time and to identify the landscape pattern variation related to growth mode in different period. This study aims at inspecting spatiotemporal characteristics of landscape pattern respond to land use changes in Xuzhou city during the period from 1985 to 2015. In this connection, we proposed a new spectral index, named the Normalized Difference Enhanced Urban Index (NDEUI), which combines data from NTL (Nighttime light) from the Defense Meteorological Satellite Program/Operational Linescan System (DMSP/OLS) with annual maximum Enhanced Vegetation Index (EVI) to reduce the detection confusion between urban areas and barren land, as well as follows. NDEUI-assisted Random Forests algorithm was implemented to obtain the land use/land cover (LULC) maps of Xuzhou in 1985, 1995, 2005 and 2015, respectively. Here, four different periods viz. 1985–1995, 1995–2005, 2005–2015 and 1985–2015 are chosen for the change analysis of land use and landscape pattern. The results indicated that the urban area has increased by about 30.65%, 10.54%, 68.77%, and 143.75% during the four periods mentioned above at the main expense of agricultural land, respectively. The spatial trend maps revealed that continuous transition from other land use types into urban land has appeared a dual-core development mode throughout the urbanization process, located at the new city region and the Jiawang district, mainly affected by the construction of new city region, freeway and the high railway station. Furthermore, we quantified the patch complexity, aggregation, connectivity and diversity of landscape employing a number of landscape metrics to represent the changes of landscape pattern at both class and landscape level, affected by urbanization during the study period. The results showed that with regard to the four aspects of landscape pattern, there were considerable differences among the four years, mainly owing to the increasing dominance of urbanized land. Spatiotemporal variation of landscape pattern was also conducted on the basis of subgrids in 900 m × 900 m. Combined with the land use changes and spatiotemporal variation of landscape pattern, it can be concluded that different urbanization modes and intensity result in variously the spatiotemporal evolution of landscape patterns. For Xuzhou city, the urban growth mainly appeared a leapfrog mode alone both sides of the roads during the period of 1985 to 1995, and then shifted into edge-expansion mode during the period from 1995 to 2005, whereas the edge-expansion and leapfrog modes coexisted for the period from 2005 to 2015. The high valuable spatiotemporal information generated utilizing RS and GIS in this study may give assistance to urban planners and policymakers to well understand urban dynamics and evaluate their spatiotemporal and environmental impacts at a local level for the sake of sustainable urban planning in the future.

The Late Heavy Bombardment Period (4.1 to 3.8 billion years ago) of heightened impact cratering activity on young Earth is likely the driving force for the origin of life. During the Eoarchean, asteroids such as carbonaceous chondrites delivered the building blocks of life and water to early Earth. Asteroid collisions created innumerable hydrothermal crater lakes in the Eoarchean crust which inadvertently became the perfect cradle for prebiotic chemistry. These hydrothermal crater lakes were filled with cosmic water and the building blocks of life. forming a thick prebiotic soup. The unique combination of exogenous delivery of extraterrestrial building blocks of life, and the endogenous biosynthesis in hydrothermal impact crater lakes very likely gave rise to life. A new symbiotic model for the origin of life within the hydrothermal crater lakes is here proposed. In this scenario, life arose around four billion years ago through five hierarchical stages of increasing molecular complexity: cosmic, geologic, chemical, information, and biological. During the prebiotic synthesis, membranes first appeared in the hydrothermal crater lakes, followed by the simultaneous origin of RNA and protein molecules, creating the RNA/protein world. These proteins were noncoded protein enzymes that facilitated chemical reactions. RNA molecules formed in the hydrothermal crater basin by polymerization of the nucleotides on the montmorillonite mineral substrate. Similarly, the initial synthesis of abiotic protein enzymes was mediated by the condensation of amino acids on pyrite surfaces. The regular wet-dry cycles within the crater lakes assisted further concentration, condensation, and polymerization of the RNAs and proteins. Lipid membranes randomly encapsulated amino acids, RNA, and protein molecules from the prebiotic soup to initiate a molecular symbiosis inside the protocells, this led to the hierarchical emergence of several cell components. As the role of protein enzymes became essential for catalytic process in the RNA/protein world, Darwinian selection from noncoded to coded protein synthesis led to translation systems and the genetic code, heralding the information stage. In this stage, the biochemical pathways suggest the successive emergence of translation machineries such as tRNAs, aaRS, mRNAs, and of ribosomes for protein synthesis. The molecular attraction between tRNA and amino acid led to the emergence of translation machinery and the genetic code.  tRNA is an ancient molecule that created mRNA for the purpose of storing amino acid information like a digital strip. Each mRNA strand became the storage device for genetic information that encoded the amino acid sequences in triplet nucleotides. As information became available in the digital languages of the codon within mRNA, biosynthesis became less random and more organized and directional. The original translation machinery was simpler before the emergence of the ribosome than that of today. We review three main concepts on the origin and evolution of the genetic code: the stereochemical theory, the coevolution theory, and adaptive theory. We believe that these three theories are not mutually exclusive, but are compatible with our coevolution model of translations machines and the genetic code. We suggest biosynthetic pathways as the origin of the translation machine that provided the framework for the origin of the genetic code. During translation, the genetic code developed in three stages coincident with the refinement of the translation machinery: GNC code with four codons and four amino acids during interactions of pre-tRNA/pre-aaRS /pre-mRNA, SNS code consisting of 16 codons and 10 amino acids appeared during the tRNA/aaRS/mRNA interaction, and finally the universal genetic code evolved with the emergence of the tRNA/aaRS/mRNA/ribosome machine. The universal code consists of 64 codons and 20 amino acids, with a redundancy that minimizes errors in translation. To address the question of the origin of the biological information system in the RNA/protein world, we converted letter codons into numerical codons in the Universal Genetic Code Table. We developed a software called CATI (Codon-Amino Acid-Translator-Imitator) to translate randomly chosen numerical codons into corresponding amino acids and vice versa, gaining insight into how translation might have worked in the RNA/protein world. We simulated the likely biochemical pathways for the origin of translation and the genetic code using the Model-View-Controller (MVC) software framework, and the translation machinery step-by-step. We used AnyLogic software to simulate and visualize the evolution of the translation machines and the genetic code. We conclude that the emergence of the information age from the RNA/protein world was a watershed event in the origin of life about four billion years ago.

We present a new unsupervised method aimed to obtain a partition of a complex urban system in homogenous urban areas, called urban contexts. The area of study is initially partitioned in microzones, homogeneous portion of the urban system, that are the atomic reference elements for the census data. With the contribution of domain experts, we identify the physical, morphological, environmental and socio-economic indicators need to identify synthetic characteristics of urban contexts and create the fuzzy rule set necessary to determine the type of urban context. We implement the set of spatial analysis processes necessary to calculate the indicators for microzone and apply a Mamdani fuzzy rule system to classify the microzones. Finally, the partition of the area of study in urban contexts is obtained by dissolving continuous microzones belonging to the same type of urban context. Tests are performed on the Municipality of Pozzuoli (Naples - Italy); the reliability of out model is measured by comparing the results with the ones obtained by detailed analysis.

Satellite-based remote sensing technologies are utilized extensively to investigate urban thermal environments under rapid urban expansion. Current MODIS data is, however, unable to adequately represent the spatially detailed information because of its relatively coarser spatial resolution, while Landsat data can’t explore temporally the refined analysis due to the low temporal resolution. In order to resolve this situation, we used MODIS and Landsat data to generate “Landsat-like” data by using the flexible spatiotemporal data fusion method (FSDAF), and then studied spatiotemporal variation of land surface temperature (LST) and its driving factors. The results showed that 1) The estimated “Landsat-like” data have high precision; 2) By comparing 2013 and 2016 datasets, LST increases ranging from 1.8°C to 4°C were measurable in areas where the impervious surface area (ISA) increased, while LST decreases ranging from -3.52°C to -0.70°C were detected in areas where ISA decreased; 3) LST has a strongly negative relationship with the Normalized Difference Vegetation Index (NDVI), and a strongly positive relationship with Normalized Difference Built Index (NDBI) in summer; and 4) LST is well correlated with Building density (BD), in a complex conic mode, and LST may increase by 0.460°C to 0.786°C when BD increases by 0.1. Our findings can provide information useful for mitigating undesirable thermal conditions and for long-term urban thermal environmental management.

In this study, the fuzzy c- means classifier has been studied with nine other similarity and dissimilarity measures: Manhattan distance, chessboard distance, Bray-Curtis distance, Canberra, Cosine distance, correlation distance, mean absolute difference, median absolute difference and normalised squared Euclidean distance. Both single and composite modes were used with a varying weight constant (m) and also at different α-cuts. The two best single norms obtained were combined to study the effect of composite norms on the datasets used. An image to image accuracy check was conducted to assess the accuracy of the classified images. Fuzzy Error Matrix (FERM) was applied to measure the accuracy assessment outcomes for a Landsat-8 dataset with respect to the Formosat-2 dataset. To conclude FCM classifier with Cosine norm performed better than the conventional Euclidean norm. But, due to the incapability of the FCM classifier to handle noise properly, the classification accuracy was around 75%.

The European Space Agency (ESA) defines as Earth observation (EO) Level 2 information product a single-date multi-spectral (MS) image corrected for atmospheric, adjacency and topographic effects, stacked with its data-derived scene classification map (SCM), whose thematic map legend includes quality layers cloud and cloud-shadow. ESA EO Level 2 product generation is an inherently ill-posed computer vision (CV) problem never accomplished to date in operating mode by any EO data provider at the ground segment. Herein, it is considered: (I) necessary not sufficient pre-condition for the yet-unaccomplished dependent problems of semantic content-based image retrieval (SCBIR) and semantics-enabled information/knowledge discovery (SEIKD) in multi-source EO big data cubes. (II) Synonym of EO Analysis Ready Data (ARD) format. (III) Equivalent to a horizontal policy for background developments in Space Economy 4.0. In compliance with the GEO-CEOS Quality Assurance Framework for EO Calibration/Validation guidelines, to contribute toward filling an analytic and pragmatic information gap from multi-sensor EO big data to timely, comprehensive and operational EO value-adding information products and services, this work presents an innovative AutoCloud+ CV software toolbox for cloud and cloud-shadow quality layer detection in ESA EO Level 2 product. In vision, spatial information dominates color information. Inspired by this true-fact, the inherently ill-posed AutoCloud+ CV software was conditioned, designed and implemented to be “universal”, meaning fully automated (no human-machine interaction is required), near real-time, robust to changes in input data and scalable to changes in MS imaging sensor’s spatial and spectral resolution specifications.

In terrain analysis and hydrological modeling, surface depressions (or sinks) in a digital elevation model (DEM) are commonly treated as artifacts and thus filled and removed to create a depressionless DEM. Various algorithms have been developed to identify and fill depressions in DEMs during the past decades. However, few studies have attempted to delineate and quantify the nested hierarchy of actual depressions, which can provide crucial information for characterizing surface hydrologic connectivity and simulating the fill-merge-spill hydrological process. In this paper, we present an innovative and efficient algorithm for delineating and quantifying nested depressions in DEMs using the level-set method based on graph theory. The proposed level-set method emulates water level decreasing from the spill point along the depression boundary to the lowest point at the bottom of a depression. By tracing the dynamic topological changes (i.e., depression splitting/merging) within a compound depression, the level-set method can construct topological graphs and derive geometric properties of the nested depressions. The experimental results of two fine-resolution LiDAR-derived DEMs show that the raster-based level-set algorithm is much more efficient (~150 times faster) than the vector-based contour tree method. The proposed level-set algorithm has great potential for being applied to large-scale ecohydrological analysis and watershed modeling.

Reuse of patterns is a self-evident approach for managing interoperability concerns. Although patterns for resolving interoperability barriers exist in the literature, no study exists on adoption of interoperability patterns by Geographic Information Systems (GIS) practitioners in industry. Thus there is limited understanding of pattern re-usability, yet the advantages offered by interoperability patterns provide a reasonably sound justification for their usage. This paper examines the adoption of proven interoperability best practices in the GIS industry. An empirical study that involved the use of semi-structured interviews was employed to gather data from GIS developers on domain interoperability best practices. Results indicated that industry and communities of practice have been converging on the technical level to ensure interoperability of GIS concerns. Semantic interoperability and related patterns are least understood, yet semantic barriers still exist. This is partly due to the complexity associated with the top-down approach used to develop semantic interoperability solutions. Therefore, this study proposes research into resolving barriers in the adoption of interoperability patterns that reduce complexity while solving semantic interoperability barriers.

Salinity intrusion is one of the most serious consequences of climate change coupled with rising sea level that significantly affects agricultural activities in many parts of the world. This phenomenon has increasingly become more serious and frequently occurred in the Mekong Delta of Vietnam. As a result, Vietnam has been ranked among top five countries where have been devastatingly impacted by climate change, in particular, its Tra Vinh Province characterized by coastal plain and alluvial deposit. In addition, this area is of the tropical monsoon zone of long rainy season with source of salt brought from the sea by the tides and sea level rise. Regions that are contaminated by salt are located in lowland and often suffer from floods linking to tidal effects with salty water from river systems and channels. Soil salinity evaluation is critical for coastal protection, restoration, and agricultural planning since it can be considered as an agricultural indicator to evaluate quality of soil. Here, we attempt to estimate the soil salinity in Tra Vinh Province, in the Mekong Delta of Vietnam. Landsat 8 OLI images are utilized to derive indices for soil salinity evaluation including single bands, Vegetation Soil Salinity Index (VSSI), Soil Adjusted Vegetation Index (SAVI), Normalized Difference Vegetation Index (NDVI), and Normalized Difference Salinity Index (NDSI). Subsequently, satistical analysis between soil salinity, electrical conductivity (EC, dS/m), and environmental indices derived from Landsat 8 OLI image is performed. Results indicate that spectral value of Near Infrared (NIR) band and VSSI are highly correlated with EC (R2 = 0.7779 and R2 = 0.6957, respectively) in comparison with the other indices. Comparative results show that soil salinity derived from Landsat 8 is consistent with in situ data. Findings of this study demonstrate that Landsat 8 OLI images reveal a high potential for spatiotemporally monitoring the magnitude of soil salinity at the top soil layer. Outcomes of this study are useful for agricultural activities, planners, and farmers by providing the base map of soil salinity contamination for better selection of accomodating crop types to reduce economical lost in the context of climate change. Our proposed method that estimates soil salinity using satellite-derived variables can be applied in the other regions.

With the accelerating urbanization process, the population increasingly concentrates in urban areas. In view of the special situation in China and a series of problems in the process of rapid urbanization, there were no reasonable measures for optimizing the population pattern. This study explored the distribution pattern of the Chinese population and proposed an optimization plan for the population distribution using GIS analysis. The main findings were as follows. (1) From 2010 to 2015, the distribution of population density in China presented a pattern of high in the southeast and low in the northwest based on the county-level administrative regions. The population still showed a tendency to migrate to the southeast of the country based on the “Hu Huanyong Line”. (2) There was a great difference in the land use efficiency in terms of population and economic production in China. The economic concentration in China was higher than the population concentration. In the areas where population and economic production were aggregated, GDP per capita and land use efficiency were higher. (3) Based on the land use efficiency in terms of population and economic production, the optimized urbanization plan of “1+4+11” for China’s urbanization was put forward, namely, one national-level aggregated area of population and economic production, 4 regional-level aggregated areas of population and economic production, and 11 local regionally aggregated areas of population and economic production. This optimization plan for urbanization represents an attempt to explore the direction of China’s urbanization, and it can be used to optimize the spatial development pattern and provide scientific guidance for the new urbanization plan.

A timely inventory of agricultural areas and crop types is an essential requirement for ensuring global food security. Satellite remote sensing has proven to be an increasingly more reliable tool to identify crop types. With the Copernicus program and its Sentinel satellites, a growing source of satellite remote sensing data is publicly available at no charge. Here we use joint Sentinel-1 radar and Sentinel-2 optical imagery to create a crop map for Belgium. To ensure homogenous radar and optical input across the country, Sentinel-1 12-day backscatter composites were created after incidence angle normalization, and Sentinel-2 NDVI images were smoothed to yield dekadal cloud-free composites. An optimized random forest classifier predicted the 8 crop types with a maximum accuracy of 82% and a kappa coefficient of 0.77. We found that a combination of radar and optical imagery always outperformed a classification based on single-sensor inputs, and that classification performance increased throughout the season until July, when differences between crop types are largest. Furthermore we showed that the concept of classification confidence derived from the random forest classifier provided insight in the reliability of the predicted class for each pixel, clearly showing that parcel borders have a lower classification confidence. We concluded that the synergistic use of radar and optical data for crop classification led to richer information increasing classification accuracies compared to optical-only classification. Further work should focus on object-level classification and crop monitoring to exploit the rich potential of combined radar and optical observations.

A large variety of remote sensing-based land use/land cover (LULC) products are currently available on national and global scales. This literature review and in-situ validation study evaluates the suitability of these products for local scale applications in the complex terrain of the Ethiopian mountains. For the review, 146 research papers were analyzed. Most studies (73%) have been published since 2013 and are based on individually computed maps. Not a single study relied on readily available LULC products. Nine readily available LULC products with 20, 30, 300, 500 and 1,000 m spatial resolution have been identified at national and global scales. To complement and extend this body of research, the recent (since 2013) LULC products were validated using 185 ground truth points collected in the Bale Mountains National Park between 1,500 and 4,385 m a.s.l. The results indicate a rather poor overall accuracy (<50%).

Urban cycling is a sustainable transport mode that cities currently try to expand. However, cities still do not take advantage of geospatial technologies to understand cycling mobility based on the behavioural patterns and difficulties faced by cyclists. This study analyses a geospatial dataset crowdsourced by urban cyclists using an experimental geo-game. About 20 participants per city recorded bicycle trips during one week and, by aggregating them, we found not only the cyclists’ preferred streets but also the frictions faced during cycling. We successfully identified 284 places potentially having frictions: 71 in Münster, Germany; 70 in Castelló, Spain; and 143 in Valletta, Malta. At such places, Participants not only recorded bicycle trips but also segments at lower speeds indicating a deviation from an ideal cycling scenario. We described the potential frictions preventing from bicycle commuting based on the distance to bicycle paths, surrounding infrastructure, and location in the urban area.

Wetlands benefits can be summarized but are not limited to their ability to store floodwaters and improve water quality, providing habitats for wildlife and supporting biodiversity, as well as aesthetic values. Over the past few decades, remote sensing and geographical information technologies has proven to be a useful and frequent applications in monitoring and mapping wetlands. Combining both optical and microwave satellite data can give significant information about the biophysical characteristics of wetlands and wetlands` vegetation. Also, fusing data from different sensors, such as radar and optical remote sensing data, can increase the wetland classification accuracy. In this paper we investigate the ability of fusion two fine spatial resolution satellite data, Sentinel-2 and the Synthetic Aperture Radar Satellite, Sentinel-1, for mapping wetlands. As a study area in this paper, Balikdami wetland located in the Anatolian part of Turkey has been selected. Both Sentinel-1 and Sentinel-2 images require pre-processing before their use. After the pre-processing, several vegetation indices calculated from the Sentinel-2 bands were included in the data set. Furthermore, an object-based classification was performed. For the accuracy assessment of the obtained results, number of random points were added over the study area. In addition, the results were compared with data from Unmanned Aerial Vehicle collected on the same data of the overpass of the Sentinel-2, and three days before the overpass of Sentinel-1 satellite. The accuracy assessment showed that the results significant and satisfying in the wetland classification using both multispectral and microwave data. The statistical results of the fusion of the optical and radar data showed high wetland mapping accuracy, with an overall classification accuracy of approximately 90% in the object-based classification. Compared with the high resolution UAV data, the classification results give promising results for mapping and monitoring not just wetlands, but also the sub-classes of the study area. For future research, multi-temporal image use and terrain data collection are recommended.

Forest canopy gaps are important for the ecosystem dynamics. Depending on tree species, small canopy openings might be also associated to intra-crown porosity and to space between crowns. Yet, little is known on the relationships between the fine-scaled pattern of canopy openings and biodiversity features. This research explored the possibility of i)- mapping forest canopy gaps from a very high resolution orthomosaic (10 cm), processed from a versatile imaging platform such as unmanned aerial vehicles (UAV), ii)- to derive patch metrics that can be tested as covariates of variables of interest for forest biodiversity monitoring. This is attempted in a test area of 240 ha covered by temperate deciduous forest types in Central Italy and containing 50 forest inventory plots of about 530 m². Correlation and linear regression techniques were used to explore relationships between patch metrics and understorey (density, development and species diversity) or forest habitat biodiversity variables (density of micro-habitat bearing trees, vertical species profile, tree species diversity). The results revealed that small openings in the canopy cover (75% smaller than 7 m²) can be faithfully extracted from UAV RGB imagery, using the red band and contrast split segmentation. Highest correlations were observed in the mixed forest (beech and turkey oak), while beech forest had the poorest ones and turkey oak forest displayed intermediate results. Moderate to strong linear relationships were found between gap metrics and understorey variables in mixed forest type, with adjusted R² from linear regression ranging from 0.52 to 0.87. Equally good results, in the same forest types, were observed for forest habitat biodiversity variables (0.52<adjusted R²<0.79) with highest values found for density of trees with microhabitats and vertical species profile. In conclusion, this research highlights that UAV remote sensing can potentially provide covariate surfaces of variables of interest for forest biodiversity monitoring, conventionally collected in forest inventory plots. By integrating the two sources of data, these variables can be mapped over small forest areas with satisfactory levels of accuracy, at a much higher spatial resolution than would be possible by field-based forest inventory solely.

We studied generalization of a method for extending configurational studies to variables that are not exclusively geographical, in order to allow investigation of generic relationships in the built environment. We observed a number of limitations of the classical approach of configurational studies and we considered how a complex analytical method could be implemented in the study of non-topological variables, such as land use, noise pollution and financial or property rents. In order to do this we established a system of relationships based on a labeled primary graph. Categorization of the labels of the links was performed on a generic variable, specifically, the predominant use of roads. The graph was then contracted using an optimization algorithm, which we describe. The algorithm allows the network to be reduced to a more streamlined system that can better analyze the relationships between the different categories. Application of the method showed that it is faster to identify weaknesses in urban networks, and then take measures to resolve them. The case study concerns the pedestrianization of the Colosseum area in Rome.

This paper reports the efforts made and experiences gained in developing the Flood Prevention and Emergency Response System (FPERS) powered by Google Earth Engine, with focus on its applications at the three stages of floods. At the post-flood stage, FPERS integrates various remote sensing imageries, including Formosat-2 optical imagery, to detect and monitor barrier lakes, synthetic aperture radar imagery to derive an inundation map, and high-spatial-resolution photographs taken by unmanned aerial vehicles to evaluate damage to river channels and structures. At the pre-flood stage, a huge amount of geospatial data are integrated in FPERS and are categorized as typhoon forecast and archive, disaster prevention and warning, disaster events and analysis, or basic data and layers. At the during-flood stage, three strategies are implemented to facilitate the access of the real-time data: presenting the key information, making a sound recommendation, and supporting the decision-making. The example of Typhoon Soudelor in August of 2015 is used to demonstrate how FPERS was employed to support the work of flood prevention and emergency response from 2013 to 2016. The capability of switching among different topographic models and the flexibility of managing and searching data through a geospatial database are also explained, and suggestions are made for future works.

In this paper we will present a simplified approach for extracting the ground level – a digital terrain model (DTM) – from the surface provided in a digital surface model (DSM). Most existing algorithms try to find the ground values in a digital surface model. Our approach works the opposite direction by detecting probable above ground areas. The main advantage of our approach is the possibility to use it with incomplete DSMs containing much no data values which can be e.g. occlusions in the calculated DSM. A smoothing or filling of such original derived DSMs will destroy much information which is very useful for deriving a ground surface from the DSM. Since the presented approach needs steep edges to detect potential high objects it will fail on smoothed and filled DSMs. After presenting the algorithm it will be applied to a test area in Salzburg and compared to a terrain model freely available from the Austrian government.

It is widely acknowledged that the integration of BIM and GIS data is a crucial step forward for future 3D city modelling, but most of the research conducted so far has covered only the semantic aspects of GIS-BIM integration. We present here the results of the GeoBIM project, in which we tackled three integration problems focussing instead on aspects involving geometry processing: (i) the automated processing of complex architectural IFC models, (ii) the integration of existing GIS subsoil data in BIM, and (iii) the georeferencing of BIM models for their use in GIS software. All the problems have been studied using real world models and existing datasets made and used by practitioners in the Netherlands. For each problem, we expose in detail the issues we faced, our proposed solutions, and our recommendations for a more successful integration.

On the occasion of Kynouria and in order to achieve the protection and projection of antiquities, a web-based model is proposed for highlighting individual monuments and archaeological sites, having in mind the historical and archaeological evidence of the region, the topography, the demographic profile and the tourist infrastructure, and combining them with the development programs for cultural routes. Therefore, creating suitable databases and mapping the monuments in the area are key prerequisites of the process, as they contribute to an objective assessment of the current situation and then to make rational decisions. In this frame, modern technology provides some important planning tools (GIS, GPS, and OMS) which allow recording and mapping of data, viewing the relationships between them in the area where they appear and managing their projection. The complete study of Kynouria’s archaeological routes contains the implementation of a website using free or open-source software, which should include all the necessary procedures and the historical and archaeological information material (text, maps, and photographs).

Object. Flooding in Ukraine is a common natural phenomenon that repeats periodically and in some cases it becomes disastrous. In an average year floods on the rivers of Volyn region take place from one to three times which extend beyond the limits of the floodplain. The floodplain of Styr river is located in the historical center of Lutsk city, that`s why issues of research and forecasting of floods are very important for a given city. Methodology. Using modern technologies of geodesy and remote sensing allows to quickly determine and predict the floodplain area of settlements. Based on the statistical data of the Volyn Regional Center for Hydrometeorology during the 7 year period 2011-2017 about water levels of the river Styr. We conducted mathematical modeling of fluctuations of water levels within the territory of Lutsk, based on creating a partial Fourier series for discrete values of middle-ten-day water levels values. The post hydrological measurements of Styr river water levels in the territory of Lutsk located on the Shevchenko Street comply with an altitude 172.87 meters. Based on the data of short-term flood forecasting in February and March, and relief data from the Department of Architecture and Urban Development of Volyn State Administration, we conducted visualization of the results using geographic information system QGIS. Results. The results of mathematical processing were the basis for geoinformation simulation of flooded areas using remote sensing data that are publicly available. Use of statistical and geospatial data in this article has great potential for further application in modeling the processes of natural and technogenic origin. Scientific novelty. The mathematical model of short-term forecasting of water levels during the flood period on the river Styr with implementation of geoinformation modeling of flooded areas using remote sensing data is proposed. Practical significance. The research results of water level changes on the Styr River and flood zones within the limits of Lutsk is proposed. The spring flood in February-March 2018, with the maximum water level 5.33 m, corresponds to an absolute mark of 178.20 m, which is forecasted in this article.

A multi-decadal change analysis of the irrigation ponds in Taoyuan, Taiwan was conducted by using multi-source data including digitized ancient maps, declassified single-band CORONA satellite images, and multispectral SPOT images. Supervised LULC classifications were conducted using four textural features derived from the single-band CORONA images and spectral features derived from SPOT images. Post-classification analysis revealed that the number of irrigation ponds in the study area decreased during the post-World War II farmland consolidation period (1945 – 1965) and the subsequent industrialization period (1970 – 2000). However, efforts on restoration of irrigation ponds in recent years have resulted in gradual increases in the number (9%) and total area (12%) of irrigation ponds in the study area.

Tree species composition is an important key element for biodiversity and sustainable forest management, and hyperspectral data provide detailed spectral information, which can be used for tree species classification. There are two main challenges for using hyperspectral imagery: a) Hughes phenomena, meaning by increasing the number of bands in hyperspectral imagery, the number of required classification samples would increase exponentially, and b) in a more complex environment, such as riparian mixed forest, focusing on spectral variability per pixel may not be adequate for definability of tree species. Therefore, the focus of this study is to assess spectral-spatial dimensionality reduction of airborne hyperspectral imagery by using minim noise fraction (MNF) transformation, and object-based image analysis (OBIA). An airborne prism experiment (APEX) hyperspectral imagery was used. A study area was a riparian mixed forest located along the Salzach river, and six tree species including Picea abies, Populus (canadensis and balsamifera), Fraxinus excelsior, Alnus incana, and Salix alba were selected. A machine learning algorithm random forest (RF) was used to train and apply a prediction model for classification. Using a spectral dimensionality reduced APEX, a pixel-level classification was also done. According to a confusion matrix, the object-level classification of MNF-derived components achieved the overall accuracy of 85 %, and kappa coefficient of 0.805. The performance of classes according to producer’s accuracy varied between 80% for Fraxinus excelsior, Alnus incana, and Populus canadensis to 90% for Salix alba and Picea abies. Comparison the results to a pixel-level classification, showed a better performance of object-level classification (an overall accuracy of 63% and Kappa coefficient of 0.559 were achieved for pixel-level classification). The performance of classes using pixel-based classification varied 45 % for Alnus incana to 80% for Picea abies. In general, Spectral-spatial complexity reduction using MNF transformation and object-level classification yielded a statistically satisfactory results.

The natural and man-made landscape settings in the Khyber Pakhtunkhwa (KPK) province of Northern Pakistan have significantly changed in the last decade due to increasing demands of urbanized populace, migration, two major natural disasters, and climate change. The aim of this study is to analyze land possession, income and land-use diversification of KPK administered Charsadda district. Field data is collected through a sample survey. Furthermore, freely available Landsat 7 satellite images are used to classify land-use classes (e.g., vegetation, built-up) for two different years (e.g., 2005 and 2017) for cross-verification and comparison. The highest 45% of 80% land-owners occupy land-area between 1–10 hectares. Annually, the highest 57.5% of the total farmers / employed in-habitants of surveyed regions earn more than Rs. 100,000 or ~$970 from agriculture activities. About 41.9% land-area covered by agricultural-land is transformed into built-up area since 2007, which is attributed to the increasing demand for buildings and commercial markets. The highest diversification is reported in Naguman area of Charsadda district followed by Rajjar and Niami. Population growth and huge migration of displaced persons from neighboring tribal areas are likely to be few factors which contributed to such a drastic change in land-use pattern since 2007 to 2017. Urgent attention of the policy makers, agricultural experts and society is required to minimize land degradation and to thwart further agricultural-land loss.

The urban heat island (UHI) becomes more and more serious with the acceleration of urbanization. Many researchers have shown interest in studying the UHI by using remote sensing data. But these studies rarely examine the mountainous cities. The studies on UHI in mountainous cities often used empirical parameters to estimate the land surface temperature (LST), and lacked satellite-ground synchronous experiment to test the accuracy. This paper revised the parameters in mono-window algorithm used to retrieve the LST according to the characteristics of mountainous cities. This study examined the spatial and temporal patterns of the UHI intensity in Chongqing, a typical mountainous city, and its relationship with land cover from 2007 to 2011 based on the Landsat TM data and the improved method. The accuracy of the LST derivation increased by about 1°C compared to the traditional method. The high LST areas increased and extended from the downtown to suburban area each year, but the rate of change decreased. The UHI is dramatically impacted by the rivers. There is a good relationship between the urban sprawl and the UHI. The LST was reduced by about 1°C within a 300m distance from large urban fringe green spaces. The urban landscape parks had a strong effect relieving the UHI at a 100m distance. The LST was reduced by about 0.5°C. The study greatly improves the accuracy of LST derivation, and provides a reliable parameters for the UHI researched in mountainous city.

QPEs (Quantitative Precipitation Estimates) obtained from remote sensing or ground-based radars could complement or even be an alternative to rain gauge readings. However, to be used in operational applications, a validation process has to be carried out, usually by comparing their estimates with those of a rain gauges network. In this paper, the accuracy of two QPEs are evaluated for three extreme precipitation events in the last decade in the southeast of the Iberian Peninsula. The first QPE is PERSIANN-CCS, a satellite-based QPE. The second is a meteorological radar with Doppler capabilities that works in the C band. Pixel-to-point comparisons are made between the values offered by the QPEs and those obtained by two networks of rain gauges. The results obtained indicate that both QPEs were well below the rain gauge values, especially in extreme rainfall time slots. There seems to be a weak linear association between the value of the discrepancies and the precipitation value of the QPEs. It does not seem that radar is more accurate than PERSIANN-CCS, despite its larger spatial resolution and its commonly higher effectiveness. The main conclusion is that neither PERSIANN-CCS nor radar, without empirical calibration, are acceptable QPEs for the real-time monitoring of meteorological extremes in the southeast of the Iberian Peninsula.

This work presents datasets that can be used for getting a good understanding of an essential geoscience content knowledge that describe earth’s coordinate systems. This include coordinate system used for spherical/spheroidal earth with latitudes and longitudes and their subsequent transformations to 2d maps on a variety of media (paper as well as digital) using the process of map projections. The datasets include PDF documents that are embedded with 3d models, animations and mathematical equations. The dataset has separate PDF documents for geographic (for spherical earth) and projected (2d) coordinate systems. Additionally, the data set include individual 3d models that can be used in various digital systems (including apps) and the animations in mp4 format that can be watched on most of the modern digital devices.

Land cover change detection (LCCD) based on bi-temporal remote sensing images plays an important role in the inventory of land cover change. Due to the benefit of having spatial dependency properties within the image space while using remote sensing images for detecting land cover change, many contextual information based change detection methods have been proposed during past decades. However, there is still a space for improvement in accuracies and usability of LCCD. In this paper, a LCCD method based on adaptive contextual information is proposed. First, an adaptive region is constructed by gradually detecting the spectral similarity surrounding a central pixel. Second, the Euclidean distance between pairwise extended regions is calculated to measure the change magnitude between the pairwise central pixels of bi-temporal images. While the whole bi-temporal images are scanned pixel-by-pixel, the change magnitude image (CMI) can be generated. Then, the Otsu or a manual threshold is employed to acquire the binary change detection map (BCDM). The detection accuracies of the proposed approach are investigated by two land cover change cases with Landsat bi-temporal remote sensing images. In comparison to several widely used change detection methods, the proposed approach can achieve a land cover change inventory map with a competitive accuracy.

Crowdsourced Data (CSD) generated by citizens is becoming more popular as its potential utilisation in many applications is increasing due to its currency and availability. However, the quality of CSD, including its relevance, is often questioned as the data is not generated by professionals nor follows standard data collection procedures. The quality of CSD can be assessed according to a range of attributes including its relevance. Information relevance has been explored through using in Geographic Information Retrieval (GIR) techniques to identify relevant information. This research tested a relevance assessment approach for CSD by adapting relevance assessment techniques available in the GIR domain. The thematic and geographic relevance were assessed using the Term Frequency-Inverse Document Frequency (TF-IDF), Vector Space Model (VSM) and Natural Language Processing (NLP) techniques. The thematic and geographic specificities of the queries were calculated as 0.44 and 0.67 respectively, which indicates the queries used were more geographically specific than thematically specific. The Spearman's rho value of 0.62 indicated that the final ranked relevance lists showed reasonable agreement with a manually classified list and confirmed the potential of the approach for CSD relevance assessment for other possible crowdsourced data analysis.

Historical maps constitute unique sources of retrospective geographic information. Recently, several map archives containing map series covering large spatial and temporal extents have been systematically scanned and made available to the public. The geographic information contained in such data archives allows extending geospatial analysis retrospectively beyond the era of digital cartography. However, given the large data volumes of such archives and the low graphical quality of older map sheets, the processes to extract geographic information need to be automated to the highest degree possible. In order to understand the salient characteristics, data quality variation, and potential challenges in large-scale information extraction tasks, preparatory analytical steps are required to efficiently assess spatio-temporal coverage, approximate map content, and spatial accuracy of such georeferenced map archives across different cartographic scales. Such preparatory steps are often neglected or ignored in the map processing literature but represent highly critical phases that lay the foundation for any subsequent computational analysis and recognition. In this contribution we demonstrate how such preparatory analyses can be conducted using classical analytical and cartographic techniques as well as visual-analytical data mining tools originating from machine learning and data science, exemplified for the United States Geological Survey topographic map and Sanborn fire insurance map archives.

Laser scanning systems make use of Light Detection and Ranging (LiDAR) technology to acquire accurately georeferenced sets of dense 3D point cloud data. The information acquired using these systems produces better knowledge about the terrain objects which are inherently 3D in nature. The LiDAR data acquired from mobile, airborne or terrestrial platforms provides several benefit over conventional sources of data acquisition in terms of accuracy, resolution and attributes. However, the large volume and scale of LiDAR data have inhibited the development of automated feature extraction algorithms due to the extensive computational cost involved in it. Moreover, the heterogeneously distributed point cloud, which represents objects with varying size, point density, holes and complicated structures pose a great challenge for data processing. Currently, geospatial database systems do not provide a robust solution for efficient storage and accessibility of raw data in a way that data processing could be applied based on optimal spatial extent. In this paper, we present Global LiDAR and Imagery Mobile Processing Spatial Environment (GLIMPSE) system that provides a framework for storage, management and integration of 3D LiDAR data acquired from multiple platforms. The system facilitates an efficient accessibility to the raw dataset, which is hierarchically represented in a geographically meaningful way. We utilise the GLIMPSE system to automatically extract road median from Airborne Laser Scanning (ALS) point cloud. In the first part of this paper, we detail an approach to efficiently retrieve the point cloud data from the GLIMPSE system for a particular geographic area based on user requirements. In the second part, we present an algorithm to automatically extract road median from the retrieved LiDAR data. The developed road median extraction algorithm utilises the LiDAR elevation and intensity attributes to distinguish the median from the road surface. We successfully tested our algorithms on two road sections consisting of distinct road median types based on concrete and grass-hedge barriers. The use of GLIMPSE improved the efficiency of the road median extraction in terms of fast accessibility to ALS point cloud data for the required road sections. The developed system and its associated algorithms provide a comprehensive solution to the user's requirement for an efficient storage, integration, retrieval and processing of large volumes of LiDAR point cloud data. These findings and knowledge contribute to a more rapid, cost-effective and comprehensive approach to surveying road networks.

Degradation of environment quality is currently the prime cause of the recent occurrence of natural disasters; it also contributes in the increase of the area that is prone to natural disasters. This research is aimed to map the potential of areas around Pesanggrahan river in DKI Jakarta by segmenting the Digital Elevation Model derived from LIDAR data. The objective of this segmentation is to find the watershed lines of the DEM image. Data processing in this research is using LIDAR data which take the ground surface data, which is overlaid with Jakarta river map and subsequently, the data is then segmented the image. The expected result of the research is the flood potential area information, especially along the Pesanggrahan river in South Jakarta.

Allama Iqbal Open University (AIOU) is the largest distance learning institute of Pakistan and providing education to 1.4 million students. This is fairly a large setup across the country where students are geographically distributed. Currently the system works on a manual approach which is not efficient. Allocation of tutors and study centers to students plays a key role in distance learning for a better learning environment. Assigning tutors and study centers to distance learning students is a challenging task when there is huge geographical spread. The utilization of geospatial technologies in open and distance learning can fix allocation problems. This research analyzes the real data of twin cities Islamabad and Rawalpindi. The results show that the geospatial technologies can be used for efficient and proper resource utilization and allocation, which in turn can save the time and money. The overall idea fits into improved distance learning framework and related analytics.

Forest fire is a major ecological disaster, which has economic, social and environmental impacts on humans and also causes the loss of biodiversity. Forest officials issue the warnings to the public on the basis of fire danger index classes. There is no fire danger index for the country India due to the sparsely distributed meteorological stations. Previous studies suggest that Static Fire Danger Index (SFDI) as well as Dynamic Fire Danger Index (DFDI) have been derived from the satellite datasets. In this study, we have made an attempt to integrate both the Static and Dynamic fire danger indices and also used the Near Real Time (NRT) data sets that can be available for download through NASA FTP server after one hour of the satellite overpass. In this study, DFDI has been calculated from the MODIS Terra NRT Land Surface Temperature (MOD11_L2) and MODIS TERRA NRT surface reflectance MOD09. Finally, The Forest Fire Danger Index (FFDI) has been developed from the static and dynamic fire danger indices by the additive model and the overall accuracy was around 81.27%. Thus, the FFDI has been useful to predict the fire danger accurately and can be useful anywhere, where the meteorological stations are un-available.

 Floods in Ukraine is a common natural phenomenon that repeats periodically and in some cases it becomes disastrous signs. In an average year in the rivers of Volyn passes from one to three floods with going beyond the limits of the floodplain. Floodplain of Styr river is located in the historical center of Lutsk city, that`s why issues of research and forecasting of floods are very important for a given city. Using modern technologies of geodesy and remote sensing allows to quickly determine and predict the floodplain area of settlements. The research results of water level changes on the Styr River and flood zones within the limits of Lutsk is proposed. The mathematical model of short-term forecasting of water level in flood period on the river Styr with implementation of geoinformation modeling of flooded areas using remote sensing data is proposed.

1) Background: We analyzed the corpus of three geoscientific journals to investigate if there are enough locational references in research articles to apply a geographical search method, on the example of New Zealand. 2) Methods: Based on all available abstracts and all freely available papers of the New Zealand Journal of Geology and Geophysics, the New Zealand Journal of Marine and Freshwater Research, and the Journal of Hydrology, New Zealand, we searched title, abstracts and full texts for place name occurrences that match records from the official Land Information New Zealand (LINZ) gazetteer. We generated ISO standard compliant metadata records for each article including the spatial references and make them available in a public catalogue service. This catalogue can be queried for articles based on authors, titles, keywords, topics as well as by spatial reference. We visualize the results in a map to show which area the research articles are about. 3) Results: We outline the methodology and technical framework for the geo-referencing of the journal articles and the platform design for this knowledge inventory. The results indicate that the use of well-crafted abstracts for journal articles with carefully chosen place names of relevance for the article provides a guideline for geographically referencing unstructured information like journal articles and reports in order to make such resources discoverable through geographical queries. 4) Conclusion: This approach can actively support integrated holistic assessment of water resources and support decision making.

Semantic Web technologies are being increasingly adopted by the geospatial community during last decade through the utilization of open standards for expressing and serving geospatial data. This was also dramatically assisted by an ever increasing access and usage of geographic mapping and location-based services via smart devices in people’s daily activities. In this paper we explore the developmental framework of a pure Javascript client-side GIS platform exclusively based on invocable geospatial Web services. We also extend Javascript utilization on the server side by deploying a node server acting as a bridge between open source WPS libraries and popular geoprocessing engines. The vehicle for such an exploration is a cross platform Web browser capable of interpreting Javascript commands to achieve interaction with geospatial providers. The tool is a generic Web interface providing capabilities of acquiring spatial datasets, composing layouts and applying geospatial processes. In an ideal form the end-user will have to identify those services, which satisfy a geo-related need and put them in the appropriate row. The final output may act as a potential collector of freely available geospatial web services. Its server-side components may exploit geospatial processing suppliers composing that way a light-weight fully transparent open Web GIS platform.

Dislocation is one of the major challenges in unmanned aerial vehicle (UAV) image stitching. In this paper, we propose a new dynamic programming for seamlessly stitching UAV images using optical flow. Our solution consists of two steps: Firstly, an image-matching algorithm is used to correct the images so that they are in the same coordinate system. Secondly, a new dynamic programming algorithm is develop based on the concept of a stereo dual-channel energy accumulation using optical flow. A new energy aggregation and traversal strategy is adopted in our solution, which can find a more optimal seam line for adjacent image stitching. Our algorithm overcomes the theoretical limitation of the classical Duplaquet algorithm. Experiments show that the algorithm can effectively solve the dislocation problem in UAV image stitching. Beyond that, our solution is also direction-independent, which has more adaptability and robustness for UAV images.

The early-season area estimation of the winter wheat crop as a strategic product is important for decision makers. Classification of multi-temporal images is an approach which is affected by many factors like appropriate training sample size, proper frequency and acquisition times, vegetation indices (VIs) type, temporal gradient of spectral bands and VIs, appropriate classifier and missed values because of cloudy conditions. This paper addresses the impact of appropriate frequency and acquisition times and VIs type along with the spectral and VI gradient on random forest (RF) classifier when missed values exist in multi-temporal images. To investigate the appropriate temporal resolution for image acquisition, the study area was selected on an overlapping area between two LDCM paths. In our developed method, the miss values of cloudy bands for each pixel are retrieved by the mean of k-nearest ordinary pixels. Then the multi-temporal image analysis is performed by considering different scenarios provided by decision makers in terms of desired crop types that should be extracted at early-season in the study areas. The classification results obtained by the RF decrease by 1.6% when temporally missed values retrieved by the proposed method, which is an acceptable result. Moreover, the experimental results demonstrated that if temporal resolution of Landsat 8 increased to one week the classification task can be conducted earlier with almost better results in terms of OA and kappa. The impact of incorporating VIs along with the temporal gradients of spectral bands and VIs as new features in RF demonstrated that the OA and Kappa are improved 3.1% and 6.6%, respectively. Furthermore, the obtained result showed that if only one image from seasonal changes of crops is available, the temporal gradient of VIs and spectral bands play the main role to discriminate remarkably wheat from barley. The experiments also demonstrated that if both wheat and barley merge to a single class the crop area can be estimated two months earlier with 97.1 and 93.5 in terms of OA and kappa, respectively.

Due to their strong immersion and real-time interactivity, helmet-mounted VR devices are becoming increasingly popular. Based on these devices, an immersive virtual geographic environment (VGE) provides a promising method for research into crowd behavior in an emergency. However, the current cheaper helmet-mounted VR devices are not popular enough and will continue to coexist with PC-based systems for a long time. Therefore, a heterogeneous distributed virtual geographic environment (HDVGE) could be a feasible solution to solve the heterogeneous problems caused by various types of clients, and support implementation of virtual crowd evacuation experiments, with large numbers of concurrent participants. In this study, we developed an HDVGE framework and put forward a set of design principles to define the similarities between the real world and the VGE. We discussed the HDVGE architecture and proposed an abstract interaction layer, a protocol-based interaction algorithm and an adjusted dead reckoning algorithm to solve the heterogeneous distributed problems. We then implemented an HDVGE prototype system focusing on subway fire evacuation experiments. Two types of clients are considered in the system, PC and all-in-one VR. Finally, we evaluated the performances of the prototype system and the key algorithms. The results showed that in a low-latency LAN environment, the prototype system can smoothly support 90 concurrent users consisting of PC and all-in-one VR clients. HDVGE could serve as a new means of obtaining observational data about individual and group behavior in support of human geography research.

Cash crops are agricultural crops intended to be sold for profit as opposed to subsistence crops, meant to support the producer, or to support livestock. Since cash crops are intended for future sale, they translate into large financial value when considered on a wide geographical scale, so their production directly involves financial risk. At a national level, extreme weather events including destructive rain or hail, as well as drought, can have a significant impact on the overall economic balance. It is thus important to map such crops in order to set up insurance and mitigation strategies. Using locally generated data -such as municipality-level records of crop seeding- for mapping purposes implies facing a series of issues like data availability, quality, homogeneity etc. We thus opted for a different approach relying on global datasets. Global datasets ensure homogeneity and availability of data, although sometimes at the expense of precision and accuracy. A typical global approach makes use of spaceborne remote sensing, for which different land cover classification strategies are available in literature at different levels of cost and accuracy. We selected the optimal strategy in the perspective of a global processing chain. Thanks to a specifically developed strategy for fusing unsupervised classification results with environmental constraints and other geospatial inputs including ground-based data, we managed to obtain good classification results despite the constraints placed. The overall production process was composed using ``good-enough" algorithms at each step, ensuring that the precision, accuracy, and data-hunger of each algorithm was commensurate to the precision, accuracy, and amount of data available. This paper describes the tailored strategy developed on the occasion as a cooperation among different groups with diverse backgrounds, a strategy which is believed to be profitably reusable in other, similar contexts. The paper presents the problem, the constraints and the adopted solutions; it then summarizes the main findings including that efforts and costs can be saved on the side of Earth Observation data processing when additional ground-based data are available to support the mapping task.

Cloudburst is one of the most devastating and frequently occurring natural hazardous events in Indian Himalayan region. Localized deep cumulus convective clouds have a capability of giving enormous amount of rainfall over a limited horizontal area, within a short span of time. Whenever, such events occur, lead to flash floods causing landslides, house collapses, dislocation of traffic, and human casualties on a large scale. Therefore, it is necessary to predict the cloudburst inundation zones accurately to avoid damage associated with them. For this, high resolution Digital Elevation Model generated from CartoSat-1 (Stereo pair) were integrated in MIKE 11 Hydrodynamic 1D model to generate longitudinal profile of the study area and to find water level, peak discharge, flow velocity, flow width at different reaches along the Asi ganga and Bhagirathi river, to know the Cloudburst flood inundation scenario. On 3^rd August 2012 one of the major Cloudburst event occurred in Asi Ganga Valley in Indian Himalayan region which was considered for simulation of hydrodynamic model. For a Cloudburst event, 100 mm/hr rainfall was considered for the simulation of the hydrodynamic model. It is observed that the discharge rise from 50 m³/s to 549.164 m³/s (an abrupt increase of about 10 times) within 1 hr at Sangamchetty in Asiganga river and at Joshiyara area rise from 600 m³/s to 3378.69 m³/s (an abrupt increase of about 5 times) within 4 hr in Bhagirathi river. Similarly the water level rises around 3 m and 6m in Asi Ganga and Bhagirathi rivers respectively. Flash Flood inundation areas due to Cloudburst on 3^rd August 2012 were demarcated from the simulation results in GIS environment.

Spatio-temporal aspects of data lead to critical information. Sensors capture data at all scales continually so it is imperative that useful information be extracted ubiquitously and regularly. Location plays a vital part by helping understand relations between datasets. It is crucial to link developmental works with spatial attributes and current challenge is to create an open platform that manages real-time sensor data and provides critical spatial analytics atop expert domain knowledge provided in the system. That is a two-faced problem where the solution tackles not only data from multiple sources but also runs data management platform, a spatial data infrastructure(SDI) as backbone framework able to harness sensor web(SW). The paper proposes development of such a globally shared open spatial expert system(ES), SmaCiSENS, a first of a kind geo-enabled knowledge based(KB) ES for multiple fields, smarter cities to climate modeling. SmaCiSENS is integration of SW and SDI with domain KB on data and problems, ready to infer solutions. The paper describes an architecture for semantic enablement for SW, SDI; connect interfaces, functions of SDI and SW, and sensor data application program interfaces (APIs) to better manage climate modeling, geohazard, global changes, and other vital areas of attention and action.

The increased availability of high resolution remote sensor data for precision agriculture 1 applications permits users to aquire deeper and more relevant knowledge about crops states that lead 2 inevitably to better decisions. The algorithm libraries being developed and evolved around these 3 applications rely on multi-spectral or hyper-spectral data acquired by using manned or unmanned 4 platforms. The current state of the art makes thorough use of vegetational indicies to guide the 5 operational management of agricultural land plots. One of the most challenging sub-problems is 6 to correctly identify and separate crop from soil. Thresholding techniques based on Normalized 7 Difference Vegetation Index (NDVI) or other such similar metrics have the advantage of being simple, 8 easy to read transformations of the data packed with useful information. Obvious difficulties arise 9 when crop/tree and soil have similar spectral responses as in case of grass filled areas in vineyards. 10 In this case grass and canopy are close in terms of NDVI values and thresholding techniques will 11 generally fail. Radiometric approaches could be integrated or replaced by a geometric approach that 12 is based on terrain data like Digital Surface Models (DSMs). These models are one of the ouputs 13 of orthorectification engines usually present in data acquired by using unmanned platforms. In 14 this paper we present two approaches based on DSM that are able to segment crop/tree from soil 15 while over gradient terrain. The DSM data are processed through a two dimensional data slicing or 16 reduction technique. Each slice is separately processed as a one dimensional time series to derive the 17 terrain and tree structures separately, here interpreted as object probability densities. In particular 18 the first approach is a Cartesian grid rasterization (CARSCAN) of the terrain and the second is its 19 immediate generalisation or radial grid rasterization of the DSM model (FANSCAN). The FANSCAN 20 recovers information from the original image at greater frequencies on the Fourier plane. These 21 approaches enable the identification of crop/tree from soil in case of slopes or hilly terrain without 22 any constraint on the displacement / direction of plant/tree row. The proposed algorithm uses pure 23 DSM information even if it is possible to fuse its output with other classifiers.

This research is related to the eco-hydrological problems of herbaceous wetland drying and biodiversity loss in the floodplain lakes of the Middle Basin of the Biebrza river (Poland). An experiment was set up, whose main goals were: (i) mapping the vegetation types and the temporarily or permanently flooded areas, and (ii) comparing the usefulness of C-band Sentinel-1A (S1A) and X-band TerraSAR-X/TanDEM-X (TSX/TDX) for mapping purposes. The S1A imagery was acquired on a regular basis using the dual polarization VV/VH and the Interferometric Wide Swath Mode. The TSX/TDX data were acquired in quad-pol, a fully polarimetric mode, during the Science Phase. The paper addresses the following aspects: i) wetland mapping with S1A multi-temporal series; ii) wetland mapping with fully polarimetric TSX/TDX data; iii) comparing the wetland mapping using dual polarization TSX/TDX subsets, i.e. HH-HV, HH-VV and VV-VH; iv) comparing wetland mapping using S1A and TSX/TDX data based on the same polarization (VV-VH); v) studying the suitability of the Shannon Entropy for wetland mapping; and vi) assessing the contribution of interferometric coherence for wetland classification. The experimental results show main limitations of the S1A dataset, while they highlight the good accuracy that can be achieved using the TSX/TDX data, especially those taken in fully polarimetric mode.

Olive trees have been of economic and cultural value since pre-Roman times, and continue to dominate landscapes and agriculture in many mediterranean regions. Recent mass losses of olive trees in Southern Italy due to an exotic plant pathogen highlight the need for methods that to monitor the olive trees in a landscape or region operationally. Here, we develop a method for counting olive trees from aerial photographs and test it in areas with a high diversity of olive tree ages, sizes, and shapes. This heterogeneity complicates tree counting as centennial trees often have crowns that are split into multiple segments, resembling multiple crowns, while nearby crowns often form a semi-closed canopy comprising multiple trees. Comparisons with reference counts in two 20 ha sites and over three different years indicate the automated counts tend to be reasonably accurate (median error 13%, n = 6), but heavily influenced by a few olive orchards with particularly high planting densities and a relatively closed canopy in which distinguishing individual trees is challenging. Overall, the algorithm estimated tree densities well (counting 82 to 109 trees/ha versus 87 to 104 trees/ha in the reference counts), indicating the method is suitable to track the number of olive trees over large areas.

This study compared two object-oriented land use change detection methods—detection after classification (DAC) and classification after detection (CAD) —based on a digital elevation model, slope data, and multi-temporal Landsat images (TM image for 2000 and ETM image for 2010). We noted that the overall accuracy of the DAC (86.42%) was much higher than that of the CAD (71.71%). However, a slight difference between the accuracies of the two methods exists for deciduous broadleaf forest, evergreen coniferous forest, mixed wood, upland, paddy, reserved land, and settlement. Owing to substantial spectrum differences, these land use types can be extracted using spectral indexes. The accuracy of DAC was much higher than that of CAD for industrial land, traffic land, green shrub, reservoir, lake, river, and channel, all of which share similar spectrums. The discrepancy was mainly because DAC can completely utilize various forms of information apart from spectrum information during a two-stage classification. In addition, the change-area boundary was not limited at first, but was adjustable in the process of classification. DAC can overcome smoothing effects to a great extent using multi-scale segmentations and multi-characters in detection. Although DAC yielded better results, it was more time-consuming (28 days) because it uses a two-stage classification approach. Conversely, CAD consumed less time (15 days). Thus, a hybrid of the two methods is recommended for application in land use change detection.

The key to simulating soil erosion is to calculate the vegetation cover (C) factor. Methods that apply remote sensing to calculate C factor at regional scale cannot directly use the C factor formula. That is because the C factor formula is obtained by experiment, and needs the coverage ratio data of croplands, woodlands and grasslands at standard plot scale. In this paper, we present a C factor conversion method from a standard plot to a km-sized grid based on large sample theory and multi-scale remote sensing. Results show that: 1) Compared with the existing C factor formula, our method is based on the coverage ratio of croplands, woodlands and grasslands on a km-sized grid, takes the C factor formula obtained from the standard plot experiment and applies it to regional scale. This method improves the applicability of the C factor formula, and can satisfy the need to simulate soil erosion in large areas. 2) The vegetation coverage obtained by remote sensing interpretation is significantly consistent (paired samples t-test, t = −0.03, df = 0.12, 2-tail significance p < 0.05) and significantly correlated with the measured vegetation coverage. 3) The C factor of the study area is smaller in the middle, southern and northern regions, and larger in the eastern and western regions. The main reason for that is the distribution of woodlands, the Hunshandake and Horqin sandy lands and the valleys affected by human activities. 4) The method presented in this paper is more meticulous than the C factor method based on the vegetation index, improves the applicability of the C factor formula, and can be used to simulate soil erosion on large scale and provide strong support for regional soil and water conservation planning.

Wetlands are recognized as one of the world’s most valuable natural resources. With the increasing world population, human demands on wetland resources for agricultural expansion and urban development continue to increase. In addition, global climate change has pronounced impacts on wetland ecosystems through alterations in hydrological regimes. To better manage and conserve wetland resources, we need to know the distribution and extent of wetlands and monitor their dynamic changes. Wetland maps and inventories can provide crucial information for wetland conservation, restoration, and management. Geographic Information System (GIS) and remote sensing technologies have proven to be useful for mapping and monitoring wetland resources. Recent advances in geospatial technologies have greatly increased the availability of remotely sensed imagery with better and finer spatial, temporal, and spectral resolution. This chapter presents an introduction to the uses of GIS and remote sensing technologies for wetland mapping and monitoring. A case study is presented to demonstrate the use of high-resolution light detection and ranging (LiDAR) data and aerial photographs for mapping prairie potholes and surface hydrologic flow pathways.  

The aim of this study is to evaluate three different strategies to improve classification accuracy in a highly fragmented semiarid area. i) Using different classification algorithms: Maximum Likelihood, Random Forest, Support Vector Machines and Sequential Maximum a Posteriori, with parameter optimisation in the second and third cases; ii) using different feature sets: spectral features, spectral and textural features, and spectral, textural and terrain features; and iii) using different image-sets: winter, spring, summer, autumn, winter+summer, winter+ spring+summer; and a four seasons combination. A 3-way ANOVA is used to discern which of these approaches and their interactions significantly increases accuracy. Tukey-Kramer contrast using a heteroscedasticity-consistent estimation of the kappa covariances matrix was used to check for significant differences in accuracy. The experiment was carried out with Landsat TM, ETM, and OLI images corresponding to the period 2000-2015. A combination of four images was the best way to improve accuracy. Maximum Likelihood, Random Forest and Support Vector Machines do not significantly increase accuracy when textural information is added, but do so when terrain features are taken into account. On the other hand, Sequential Maximum a Posteriori increases accuracy when textural features are used, but reduces accuracy substantially when terrain features are included. Random Forest using the three feature subsets and Sequential Maximum a Posteriori with spectral and textural features had the largest kappa values, around 0.9.

With the rapid development of satellite remote sensing technology, the volume of image datasets in many application areas is growing exponentially and the demand for Land-Cover and Land-Use change remote sensing data is growing rapidly. It is thus becoming hard to efficiently and intelligently retrieve the change information that users need from massive image databases. In this paper, content-based image retrieval is successfully applied to change detection and a content-based remote sensing image change information retrieval model is introduced. First, the construction of a new model framework for change information retrieval in a remote sensing database is described. Then, as the target content cannot be expressed by one kind of feature alone, a multiple-feature integrated retrieval model is proposed. Thirdly, an experimental prototype system that was set up to demonstrate the validity and practicability of the model is described. The proposed model is a new method of acquiring change detection information from remote sensing imagery and so can reduce the need for image pre-processing, deal with problems related toseasonal changes as well as other problems encountered in the field of change detection. Meanwhile, the new model has important implications for improving remote sensing image management and autonomous information retrieval.

The value of landscape, as part of collective heritage, can be acquired by GIS due to the multilayer approach of the spatial configuration. Proficiency in geospatial technologies in order to collect, process, analyze, interpret, visualize and communicate geographic information is being increased by undergraduate and graduate students, but in particular by those who are training to become geography teachers at secondary education. This training can be carried out through personalized learning and distance learning methodology. Personalized GIS education aims to integrate students and enhance their understanding of landscape. Some teaching experiences are shown whereby opportunities offered by WebGIS will be described, through quantitative tools and techniques that will allow this modality of learning and improve its effectiveness. Results of this research show that students, through geospatial technologies, learn landscape as a diversity of elements but also the complexity of physical and human factors involved. Several conclusions will be highlighted: i) the contribution of geospatial training to education for sustainable development; ii) spatial analysis as a mean of skills acquisition about measures for landscape conservation; iii) expanding and applying acquired knowledge to other geographic spaces and different landscapes.

The digital elevation model (DEM) is one of the key geospatial datasets used in many fields of engineering and science for countless applications. In this contribution, we assess the vertical accuracy of the Advanced Land Observing Satellite (ALOS) World 3D-30m (AW3D30) DEM using the runway method (RWYM). The RWYM utilizes the longitudinal profiles of runways which are reliable and ubiquitous reference data. A reference dataset used in this project consists of 36 runways located at various points throughout the world. The same dataset was previously used to test the accuracy of WorldDEM^TM.  Our study indicates that AW3D30 has a remarkably high RMSE of 1.78 m (one σ). However, while analyzing the results, it has become apparent that it also contains a widespread elevation anomaly. We conclude that this anomaly is the result of uncompensated sensor noise and the data processing algorithm (downsampling of the higher resolution data). We believe that this issue should be communicated to the user community. Also, we would like to note that the traditional accuracy assessment of a DEM, e.g., statistical assessment of the elevation differences = model – reference, does not allow for identification of these type of anomalies in a DEM.

Structure from Motion with Multi-View Stereo photogrammetry (SfM) is increasingly utilised in geoscience investigations as a cost-effective method of acquiring high resolution (sub-meter) topographic data, but has not been thoroughly tested in gullied savanna systems. The aim of this study was to test the accuracy of topographic models derived from aerial (via an Unmanned Aerial Vehicle, ‘UAV’) and ground-based (via a handheld digital camera, ‘Ground’) SfM in modelling a hillslope gully system in dry-tropical savanna, and to assess the strengths and limitations of the approach at different scales. A UAV survey covered an entire hillslope gully system (0.715 km²), whereas a Ground survey covered a single gully within the broader system (650 m²). SfM topographic models, including Digital Surface Models (DSM) and dense point clouds, were compared against RTK-GPS point data and a pre-existing airborne LiDAR Digital Elevation Model (DEM). Results indicate UAV SfM can deliver topographic models with a resolution and accuracy suitable to define gully systems at a hillslope scale (e.g., 0.1 m resolution with ~ 0.5 – 1.3 m elevation error), while ground-based SfM is more capable of quantifying gully morphology (e.g., 0.01 m resolution with ~ 0.1 m elevation error). Key strengths of SfM for these applications include: the production of high resolution 3D topographic models and ortho-photo mosaics, low survey instrument costs (< $AUD 3,000); and rapid survey time (4 and 2 hours for UAV and Ground survey respectively). Current limitations of SfM include: difficulties in reconstructing vegetated surfaces; uncertainty as to optimal survey and processing designs; and high computational demands. Overall, this study has demonstrated great potential for SfM to be used as a cost-effective tool to aid in the mapping, modelling and management of hillslope gully systems at different scales, in tropical savanna landscapes and elsewhere.

Remote sensing has considerable potential to contribute to the identification and reconstruction of lost hydrological systems and networks. Remote sensing-based reconstructions of palaeo-river networks have commonly employed single or limited time-span imagery, which limits their capacity to identify features in complex and varied landscape contexts. This paper presents a seasonal multi-temporal approach to the detection of palaeo-rivers over large areas based on long-term vegetation dynamics and spectral decomposition techniques. The use of multi-temporal data has allowed the overcoming of seasonal cultivation patterns and long-term visibility issues related to crop selection, large-scale irrigation and land use patterns. The application of this approach on the Sutlej-Yamuna interfluve (northwest India), a core area for the Bronze Age Indus Civilisation, has enabled the reconstruction of an unsuspectedly complex palaeo-river network comprising more than 8000 kms of palaeo-channels. It has also enabled the definition of the morphology of these relict courses, which provides insights into the environmental conditions in which they operated. These new data will contribute to a better understanding of the settlement distribution and environmental settings in which this, often considered riverine, civilisation operated.

All empirical water column correction methods have consistently been reported to require existing depth sounding data for the purpose of calibrating a simple depth retrieval model; they yield poor results over very bright or very dark bottoms. In contrast, we set out to (i) use only the relative radiance data in the image along with published data, and several new assumptions, (ii) in order to specify and operate the simplified radiative transfer equation (RTE), (iii) for the purpose of retrieving both the satellite derived bathymetry (SDB) and the water column corrected spectral reflectance over shallow seabeds. Sea truth regressions show that SDB depths retrieved by the method only need tide correction. Therefore it shall be demonstrated that, under such new assumptions, there is no need (i) for formal atmospheric correction, (ii) nor for conversion of relative radiance into calibrated reflectance , (iii) nor for existing depth sounding data, to specify the simplified RTE and produce both SDB and spectral water column corrected radiance ready for bottom typing. Moreover, the use of the panchromatic band for that purpose is introduced. Altogether, we named this process the Self-Calibrated Supervised Spectral Shallow-sea Modeler (4SM). This approach requires a trained practitioner, though, to produce its results within hours of downloading the raw image. The ideal raw image should be a “near-nadir” view, exhibit homogeneous atmosphere and water column, include some coverage of optically deep waters and bare land, and lend itself to quality removal of haze, atmospheric adjacency effect, and sun/sky glint.

Satellite scene classification is challenging because of the high variability inherent in satellite data. Although rapid progress in remote sensing techniques has been witnessed in recent years, the resolution of the available satellite images remains limited compared with the general images acquired using a common camera. On the other hand, a satellite image usually has a greater number of spectral bands than a general image, thereby permitting the multi-spectral analysis of different land materials and promoting low-resolution satellite scene recognition. This study advocates multi-spectral analysis and explores the middle-level statistics of spectral information for satellite scene representation instead of using spatial analysis. This approach is widely utilized in general image and natural scene classification and achieved promising recognition performance for different applications. The proposed multi-spectral analysis firstly learns the multi-spectral prototypes (codebook) for representing any pixel-wise spectral data, and then based on the learned codebook, a sparse coded spectral vector can be obtained with machine learning techniques. Furthermore, in order to combine the set of coded spectral vectors in a satellite scene image, we propose a hybrid aggregation (pooling) approach, instead of conventional averaging and max pooling, which includes the benefits of the two existing methods but avoids extremely noisy coded values. Experiments on three satellite datasets validated that the performance of our proposed approach is much more accurate than even the deep learning framework for spatial analysis.

With the great development of intelligent transportation systems (ITS), travel time prediction has attracted the attentions of many researchers and a large number of prediction methods have been developed. However, as an unavoidable topic, the predictability of travel time series is the basic premise for travel time prediction has received less attention than the methodology. Based on the analysis of the complexity of travel time series, this paper defines travel time predictability to express the probability of correct travel time prediction and proposes an entropy-based method to measure the upper bound of travel time predictability. Multiscale entropy is employed to quantify the complexity of travel time series, and the relationships between entropy and the upper bound of travel time predictability are presented. Empirical studies are made with vehicle trajectory data in an express road section. The effectiveness of time scales, tolerance, and series length to entropy and travel time predictability are analysis, and some valuable suggestions about the accuracy of travel time predictability are discussed. Finally, the comparisons between travel time predictability and actual prediction results from two prediction models, ARIMA and BPNN, are conducted. Experimental results demonstrate the validity and reliability of the proposed travel time predictability.

The segmentation of urban scene mobile laser scanning (MLS) data into meaningful street objects is a great challenge due to the scene complexity of street environments, especially in the vicinity of street objects such as poles and trees. This paper proposes a three-stage method for the segmentation of urban MLS data at the object level. The original unorganized point cloud is first voxelized, and all information needed is stored in the voxels. These voxels are then classified as ground and non-ground voxels. In the second stage, the whole scene is segmented into clusters by applying a density-based clustering method based on two key parameters: local density and minimum distance. In the third stage, a merging step and a re-assignment processing step are applied to address the over-segmentation problem and noise points, respectively. We tested the effectiveness of the proposed methods on two urban MLS datasets. The overall accuracies of the segmentation results for the two test sites are 98.3% and 97%, thereby validating the effectiveness of the proposed method.

Aerial image classification has become popular and has attracted extensive research efforts in recent decades. The main challenge lies in its very high spatial resolution but relatively insufficient spectral information. To this end, spatial-spectral feature extraction is a popular strategy for classification. However, parameter determination for that feature extraction is usually time-consuming and depends excessively on experience. In this paper, an automatic spatial feature extraction approach based on image raster and segmental vector data cross-analysis is proposed for the classification of very high spatial resolution (VHSR) aerial imagery. First, multi-resolution segmentation is used to generate strongly homogeneous image objects and extract corresponding vectors. Then, to automatically explore the region of a ground target, two rules, which are derived from Tobler’s First Law of Geography (TFL) and a topological relationship of vector data, are integrated to constrain the extension of a region around a central object. Third, the shape and size of the extended region are described. A final classification map is achieved through a supervised classifier using shape, size, and spectral features. Experiments on three real aerial images of VHSR (0.1 to 0.32 m) are done to evaluate effectiveness and robustness of the proposed approach. Comparisons to state-of-the-art methods demonstrate the superiority of the proposed method in VHSR image classification.

Urban sprawl propelled by rapid population growth leads to the shrinkage of productive agricultural lands and pristine forests in the suburban areas and, in turn, substantially alters ecosystem services. Hence, the quantification of urban sprawl is crucial for effective urban planning, and environmental and ecosystem management. Like many megacities in fast growing developing countries, Chennai, the capital of Tamilnadu and one of the business hubs in India, has experienced extensive urban sprawl triggered by the doubling of total population over the past three decades. We employed the Random Forest (RF) classification on Landsat imageries from 1991, 2003, and 2016, and computed spatial metrics to quantify the extent of urban sprawl within a 10km suburban buffer of Chennai. The rate of urban sprawl was quantified using Renyi’s entropy, and the urban extent was predicted for 2027 using land-use and land-cover change modeling. A 70.35% increase in urban areas was observed for the suburban periphery of Chennai between 1991 and 2016. The Renyi’s entropy value for year 2016 was ≥ 0.9, exhibiting a two-fold rate of urban sprawl. The spatial metrics values indicate that the existing urban areas of Chennai became denser and the suburban agricultural, forests and barren lands were transformed into fragmented urban settlements. The forecasted urban growth for 2027 predicts a conversion of 13670.33ha (16.57 % of the total landscape) of existing forests and agricultural lands into urban areas with an associated increase in the entropy value of 1.7. Our findings are relevant for urban planning and environmental management in Chennai and provide quantitative measures for addressing the social-ecological consequences of urban sprawl and the protection of ecosystem services.

This paper describes the first results of an Image Recognition Based Location (IRBL) for mobile application focusing on the procedure to generate a Database of range images (RGB-D). In an indoor environment, to estimate the camera position and orientation, a prior spatial knowledge of the surrounding is needed. In order to achieve this objective a complete 3D survey of two different environment (Bangbae metro station of Seoul and E.T.R.I. building in Daejeon – Republic of Korea) was performed using LiDAR (Light Detection And Ranging) instrument and the obtained scans were processed in order to obtain a spatial model of the environments. From this, two databases of reference images were generated using a specific software realized by the Geomatics group of Politecnico di Torino (ScanToRGBDImage). This tool allow to generate synthetically different RGB-D images) centered in the each scan position in the environment. Later, the external parameters (X, Y, Z, ω, φ, κ) and the range information extracted from the DB images retrieved, are used as reference information for pose estimation of a set of acquired mobile pictures in the IRBL procedure. In this paper the survey operations, the approach for generating the RGB-D images and the IRB strategy are reported. Finally the analysis of the results and the validation test are described.

The variation of aerosols, especially dust aerosol, in time and space plays an important role in climate forcing studies. Aerosols can effectively reduce land surface longwave emission and re-emit energy at a colder temperature, making estimation of downwelling surface longwave radiation (DSLR) with satellite data difficult. Using the latest atmospheric radiative transfer code (MODTRAN 5.0), we simulate the outgoing longwave radiation (OLR) and DSLR under different land surface and atmospheric profile conditions. The results show that dust aerosol has an obvious “warming” effect to longwave radiation compared with other aerosols, that aerosol longwave radiative forcing (ALRF) increased with increasing aerosol optical depth (AOD), and that the atmospheric water vapor content (WVC) is critical to the understanding of ALRF. A method is proposed to improve the accuracy of DSLR estimation from satellite data for the skies under heavy dust aerosols. The AOD and atmospheric WVC under cloud-free conditions with a relatively simple satellite-based radiation model that yields the high accurate DSLR under heavy dust aerosol are used explicitly as model input to reduce the effects of dust aerosol on the estimation of DSLR. Validations of the proposed model with satellites data and field measurements show that it estimates the DSLR accurately under heavy dust aerosol skies. The root mean square errors (RMSEs) are 20.4 W/M² and 24.2 W/M² for Terra and Aqua satellites, respectively, at the Yingke site, and the biases are 2.7 W/M² and 9.6 W/M², respectively. For the Arvaikheer site, the RMSEs are 23.2 W/M² and 19.8 W/M² for Terra and Aqua, respectively, and the biases are 7.8 W/M² and 10.5 W/M², respectively. The proposed method is especially applicable to acquire relatively high accurate DSLR under heavy dust aerosol using MODIS data with available WVC and AOD data.

Scene classification plays an important role in the intelligent processing of high-resolution satellite (HRS) remotely sensed image. In HRS image classification, multiple features, e.g. shape, color, and texture features, are employed to represent scenes from different perspectives. Accordingly, effective integration of multiple features always results in better performance compared to methods based on a single feature in the interpretation of HRS image. In this paper, we introduce a multi-task joint sparse and low-rank representation model to combine the strength of multiple features for HRS image interpretation. Specifically, a multi-task learning formulation is applied to simultaneously consider sparse and low-rank structure across multiple tasks. The proposed model is optimized as a non-smooth convex optimization problem using an accelerated proximal gradient method. Experiments on two public scene classification datasets demonstrate that the proposed method achieves remarkable performance and improves upon the state-of-art methods in respective applications.

Supervised land-use/land-cover (LULC) classifications are typically conducted using class assignment rules derived from a set of multiclass training samples. Consequently, classification accuracy varies with the training data set and is thus associated with uncertainty. In this study, we propose a bootstrap resampling and reclassification approach that can be applied for assessing not only the uncertainty in classification results of the bootstrap-training data sets, but also the classification uncertainty of individual pixels in the study area. Two measures of pixel-specific classification uncertainty, namely the maximum class probability and Shannon entropy, were derived from the class probability vector of individual pixels and used for the identification of unclassified pixels. Unclassified pixels that are identified using the traditional chi-square threshold technique represent outliers of individual LULC classes, but they are not necessarily associated with higher classification uncertainty. By contrast, unclassified pixels identified using the equal-likelihood technique are associated with higher classification uncertainty and they mostly occur on or near the borders of different land-cover.

When the distance between an obstacle and a power line is less than the discharge distance, a discharge arc can be generated, resulting in interruption of power supplies. Therefore, regular safety inspections are necessary to ensure safe operations of power grids. Tall vegetation and buildings are the key factors threatening the safe operation of extra high voltage transmission lines within a power line corridor. Manual or LiDAR based-inspections are time consuming and expensive. To make safety inspections more efficient and flexible, a low-altitude unmanned aerial vehicle remote-sensing platform equipped with optical digital camera was used to inspect power line corridors. We propose a semi-patch matching algorithm based on epipolar constraints using both correlation coefficient and the shape of its curve to extract three dimensional (3D) point clouds for a power line corridor. Virtual photography was used to transform the power line direction from approximately parallel to the epipolar line to approximately perpendicular to epipolar line to improve power line measurement accuracy. The distance between the power lines and the 3D point cloud is taken as a criterion for locating obstacles within the power line corridor automatically. Experimental results show that our proposed method is a reliable, cost effective and applicable way for practical power line inspection, and can locate obstacles within the power line corridor with measurement accuracies better than ±0.5 m.