ISPRS Int. J. Geo-Inf., Volume 8, Issue 10 (October 2019) – 41 articles

The past several years have witnessed much progress in landscape ecology and fractal analysis. In landscape ecology, the gradient model of landscape patterns (i.e., landscape gradient) has emerged as a new operating paradigm, where most landscape metrics do not apply because they were developed for the patch mosaic model. In the fractal analysis, a new definition of fractal has been proposed, and various new fractal metrics have been developed. This technical note aims to provide an intersection of these two lines of advance, which will further present an opportunity to advance geo-informatics by considering the latest progress in both landscape ecology and fractal analysis. We first present an overview of the new definition of fractal and all the fractal metrics developed under this new definition. Since the chief obstacle to geographers and landscape ecologists in applying these metrics is the lack of readily accessible methods for their easy computation, we then develop an integrated tool to compute them on landscape gradients. The developed tool facilitates the computation of these new fractal metrics. A case study was carried out with real-life landscape gradients, namely a digital terrain model. These new fractal metrics and the developed tool can be expected to facilitate the fractal characterization of the patterns of gradient landscapes and the understanding of landscape dynamics from a new perspective. Full article

Triangulated irregular networks (TINs) are widely used in terrain visualization due to their accuracy and efficiency. However, the conventional algorithm for multi-scale terrain rendering, based on TIN, has many problems, such as data redundancy and discontinuities in scale transition. To solve these issues, a method based on a detail-increment model for the construction of a continuous-scale hierarchical terrain model is proposed. First, using the algorithm of edge collapse, based on a quadric error metric (QEM), a complex terrain base model is processed to a most simplified model version. Edge collapse records at different scales are stored as compressed incremental information in order to make the rendering as simple as possible. Then, the detail-increment hierarchical terrain model is built using the incremental information and the most simplified model version. Finally, the square root of the mean minimum quadric error (MMQE), calculated by the points at each scale, is considered the smallest visible object (SVO) threshold that allows for the scale transition with the required scale or the visual range. A point cloud from Yanzhi island is converted into a hierarchical TIN model to verify the effectiveness of the proposed method. The results show that the method has low data redundancy, and no error existed in the topology. It can therefore meet the basic requirements of hierarchical visualization. Full article

Considering the multiscale characteristics of the human visual system and any natural scene, the spatial autocorrelation of remotely sensed imagery, and the multilevel spatial structure of ground targets in remote sensing images, an information-measurement approach based on a single-level geometrical mapping model can only reflect partial feature information at a single level (e.g., global statistical information and local spatial distribution information). The single mapping model cannot validly characterize the information of the multilevel and multiscale features of the spatial structures inherent in remotely sensed images. Additionally, the validity, practicability, and application range of the results of single-level mapping models are greatly limited in practical applications. In this paper, we present the multilevel geometrical mapping entropy (MGME) model to evaluate the information content of related attribute characteristics contained in remotely sensed images. Subsequently, experimental images with different types of objects, including reservoir area, farmland, water area (i.e., water and trees), and mountain area, were used to validate the performance of the proposed method. Experimental results show that the proposed method can not only reflect the difference in the information of images in terms of spectrum features, spatial structural features, and visual perception but also eliminates the inadequacy of a single-level mapping model. That is, the multilevel mapping strategy is feasible and valid. Additionally, the vector set of the MGME method and its standard deviation (Std) value can be used to further explore and study the spatial dependence of ground scenes and the difference in the spatial structural characteristics of different objects. Full article

Rapid urbanization in cities can result in a decrease in green urban areas. Reductions in green urban infrastructure pose a threat to the sustainability of cities. Up-to-date maps are important for the effective planning of urban development and the maintenance of green urban infrastructure. There are many possible ways to map vegetation; however, the most effective way is to apply machine learning methods to satellite imagery. In this study, we analyze four machine learning methods (support vector machine, random forest, artificial neural network, and the naïve Bayes classifier) for mapping green urban areas using satellite imagery from the Sentinel-2 multispectral instrument. The methods are tested on two cities in Croatia (Varaždin and Osijek). Support vector machines outperform random forest, artificial neural networks, and the naïve Bayes classifier in terms of classification accuracy (a Kappa value of 0.87 for Varaždin and 0.89 for Osijek) and performance time. Full article

Small earthquakes following a large event in the same area are typically aftershocks, which are usually less destructive than mainshocks. These aftershocks are considered mainshocks if they are larger than the previous mainshock. In this study, records of aftershocks (M > 2.5) of the Kermanshah Earthquake (M 7.3) in Iran were collected from the first second following the event to the end of September 2018. Different machine learning (ML) algorithms, including naive Bayes, k-nearest neighbors, a support vector machine, and random forests were used in conjunction with the slip distribution, Coulomb stress change on the source fault (deduced from synthetic aperture radar imagery), and orientations of neighboring active faults to predict the aftershock patterns. Seventy percent of the aftershocks were used for training based on a binary (“yes” or “no”) logic to predict locations of all aftershocks. While untested on independent datasets, receiver operating characteristic results of the same dataset indicate ML methods outperform routine Coulomb maps regarding the spatial prediction of aftershock patterns, especially when details of neighboring active faults are available. Logistic regression results, however, do not show significant differences with ML methods, as hidden information is likely better discovered using logistic regression analysis. Full article

A new Python package, pyjeo, that deals with the analysis of geospatial data has been created by the Joint Research Centre (JRC). Adopting the principles of open science, the JRC strives for transparency and reproducibility of results. In this view, it has been decided to release pyjeo as free and open software. This paper describes the design of pyjeo and how its underlying C/C++ library was ported to Python. Strengths and limitations of the design choices are discussed. In particular, the data model that allows the generation of on-the-fly data cubes is of importance. Two uses cases illustrate how pyjeo can contribute to open science. The first is an example of large-scale processing, where pyjeo was used to create a global composite of Sentinel-2 data. The second shows how pyjeo can be imported within an interactive platform for image analysis and visualization. Using an innovative mechanism that interprets Python code within a C++ library on-the-fly, users can benefit from all functions in the pyjeo package. Images are processed in deferred mode, which is ideal for prototyping new algorithms on geospatial data, and assess the suitability of the results created on the fly at any scale and location. Full article

In this paper, a MATLAB tool for the automatic detection of the best locations to install a wireless sensor network (WSN) is presented. The implemented code works directly on high-resolution 3D point clouds and aims to help in positioning sensors that are part of a network requiring inter-visibility, namely, a clear line of sight (LOS). Indeed, with the development of LiDAR and Structure from Motion technologies, there is an opportunity to directly use 3D point cloud data to perform visibility analyses. By doing so, many disadvantages of traditional modelling and analysis methods can be bypassed. The algorithm points out the optimal deployment of devices following mainly two criteria: inter-visibility (using a modified version of the Hidden Point Removal operator) and inter-distance. Furthermore, an option to prioritize significant areas is provided. The proposed method was first validated on an artificial 3D model, and then on a landslide 3D point cloud acquired from terrestrial laser scanning for the real positioning of an ultrawide-band WSN already installed in 2016. The comparison between collected data and data acquired by the WSN installed following traditional patterns has demonstrated its ability for the optimal deployment of a WSN requiring inter-visibility. Full article

Making transport systems sustainable is a topic that has attracted the attention of many researchers and urban planners. The Association of Southeast Asian Nations (ASEAN) Smart Cities Network (ASCN) was initiated to develop a sustainable transport system in the ASEAN countries. A comprehensive understanding of street networks (SNs) can contribute significantly to the achievement of this initiative. Therefore, this paper measured and compared characteristics of drivable street networks (DSNs) and walkable street networks (WSNs) of the 26 ASCN pilot cities by applying multiple network metrics. The OSMnx tool was used to download and analyse WSNs and DSNs from the OpenStreetMap. The findings present the topological and geometric characteristics of WSNs and DSNs that are diverse and characterized by different factors. The cities with orthogonal street grids, high street density, intersection density, and fewer cul-de-sacs have good accessibility to reach destinations. In contrast, some other cities have more curvilinear and circuitous SNs with many missing links to other streets, which in turn are prone to traffic disruption. The study highlights the important features of SNs that have significant implications for future designs of SNs in the ASCN whose goal is to make cities smart and liveable for ASEAN members. Full article

Understanding the spatial variability of soil health and identifying areas that share similar soil properties can help nations transition to sustainable agricultural practices. This information is particularly applicable to management decisions such as tillage, nutrient application, and soil and water conservation. This study evaluated the spatial variability and derived the optimal number of spatially contiguous regions of Nigeria’s 774 Local Government Areas (LGAs) using three soil health indicators, organic carbon (OC), bulk density (BD) and total nitrogen (TN) extracted from the Africa Soil Information Service database. Missing data were imputed using the random forest imputation method with topography and normalized difference vegetation index (NDVI) as auxiliary variables. Using an exponential covariance function, the spatial ranges for BD, SN, and OC were calculated as 18, 42, and 60 km, respectively. These were the maximum distances at which there was no correlation between the sample data points. This finding suggests that OC has high variability across Nigeria as compared with other tested indicators. The ordinary kriging (OK) technique revealed spatial dependency (positive correlation) among TN and OC on interpolated surfaces, with high values in the southern part of the county and low values in the north. The BD values were also high in the northern regions where the soils are sandy; correspondingly, TN and OC had low values. The “regionalization with dynamically constrained agglomerative clustering and partitioning” (REDCAP) technique was used to divide LGAs into a possible number of regions while optimizing a sum of squares deviation (SSD). Optimal division was not observed in the resulting number of regional partitions. Conducting the Markov Chain Monte Carlo (MCMC) method on within-zone heterogeneity (WZH) revealed three partitions (two, five, and 15 regions) as optimal, in other words, there would be no significant change in WZH after three partitions. Ensuring a proper understanding of soil spatial variability and heterogeneities (or homogeneities) could facilitate agricultural planning that combines or merges state and local governments that share the same soil health properties, rather than basing decisions on geopolitical, racial, or ethnoreligious factors. The findings of this study could be applied to understand the importance of soil heterogeneities in hydrologic modeling applications. In addition, the findings may aid decision-making bodies such as the United Nations’ Food and Agricultural Organization, the International Fund for Agricultural Development, or the World Bank in their efforts to alleviate poverty, meet future food needs, mitigate the impacts of climate change, and provide financial funding through sustainable agriculture and intervention in developing countries such as Nigeria. Full article

Land prices are among the most important parameters of urbanization and have been an important subject of urban geography studies for many years. The relationship between urban geography and land prices was examined in the first established models, which had linear and static structures. In these models, which have a radial form, cities are considered to be commercial centers. However, since the 20th century, it has been accepted that cities have structures without obvious order, consisting of many subsystems related to political, social, and economic life and space. This irregular structure that repeats itself independent of scale has a fractal geometry. Developments in the field of geographic information systems in the last 30 years have provided great convenience in analyzing the structure of cities with fractal dimensions. The geometric shapes of buildings, streets, and blocks that create the physical city form at the same time constitute the urban geometry. This study, which aims to investigate the spatial relationship between urban geometry and land prices, examines the relationship between the fractal dimension values of buildings, streets, blocks, and land prices and whether the factors of population and distance to the center have an impact on this relationship by using geostatistical methods. In this context, the fractal dimension values of urban geometry components were calculated separately in the study area, consisting of 65 neighborhoods. A two-step cluster analysis was used to determine how these obtained fractal values are dispersed geographically within the study area. By measuring the success of clustering through the independent samples t-test, it was decided which data would be used in the regression model in which the relationship between urban geometry and land prices would be established. By using exploratory factor analysis, intercorrelated data to be used in the regression model were eliminated. According to the results of the multivariate regression model, it was revealed that there was a directly proportional relationship between the fractal dimension values of building-block geometry and land prices, and an inversely proportional relationship between the fractal dimension values of street geometry and land prices. Full article

China’s rapid urbanization and industrialization have continually placed massive pressure on the country’s natural resources. The fragmented departmental administration of natural resources also intensifies the problem of sustainable use. Accordingly, China’s central government has launched natural resource administration reform from decentralization to unification. This study systematically analyzes the reform requirements from legal, organizational, and technical aspects. The right structure of China’s natural resource assets for fulfilling such requirements is examined in this work through a review of relevant legal text, and such a right structure is converted into a draft national technical standard of China’s natural resource administration on the basis of the land administration domain model (LADM). Results show that China’s natural resource administration covers lands, buildings, structures, forests, grasslands, waters, beaches, sea areas, minerals, and other fields. The types of private rights over natural resources include ownerships, land-contracted management rights (cultivated land, forest land, grassland, and water area), rights to use construction land (state-owned and collective-owned), rights to use agricultural land, rights to use homestead land, breeding rights on water areas and beaches, rights to use sea areas, rights to use uninhabited islands, and mining rights. The types of public rights over natural resources include comprehensive land use, urban and rural, sea use, and territory space planning. Furthermore, various types of these property rights can be converted into corresponding classes in LADM on the basis of the analysis of the property subject, object, and rights. Full article

Historical maps constitute an essential information for investigating the ecological and landscape features of a region over time. The integration of heritage maps in GIS models requires their digitalization and classification. This paper presents a semi-automatic procedure for the digitalization of heritage maps and the successive filtering of undesirable features such as text, symbols and boundary lines. The digitalization step is carried out using Object-based Image Analysis (OBIA) in GRASS GIS and R, combining image segmentation and machine-learning classification. The filtering step is performed by two GRASS GIS modules developed during this study and made available as GRASS GIS add-ons. The first module evaluates the size of the filter window needed for the removal of text, symbols and lines; the second module replaces the values of pixels of the category to be removed with values of the surrounding pixels. The procedure has been tested on three maps with different characteristics, the “Historical Cadaster Map for the Province of Trento” (1859), the “Italian Kingdom Forest Map” (1926) and the “Map of the potential limit of the forest in Trentino” (1992), with an average classification accuracy of 97%. These results improve the performance of classification of heritage maps compared to more classical methods, making the proposed procedure that can be applied to heterogeneous sets of maps, a viable approach. Full article

The Cellular Automata Markov model combines the cellular automata (CA) model’s ability to simulate the spatial variation of complex systems and the long-term prediction of the Markov model. In this research, we designed a parallel CA-Markov model based on the MapReduce framework. The model was divided into two main parts: A parallel Markov model based on MapReduce (Cloud-Markov), and comprehensive evaluation method of land-use changes based on cellular automata and MapReduce (Cloud-CELUC). Choosing Hangzhou as the study area and using Landsat remote-sensing images from 2006 and 2013 as the experiment data, we conducted three experiments to evaluate the parallel CA-Markov model on the Hadoop environment. Efficiency evaluations were conducted to compare Cloud-Markov and Cloud-CELUC with different numbers of data. The results showed that the accelerated ratios of Cloud-Markov and Cloud-CELUC were 3.43 and 1.86, respectively, compared with their serial algorithms. The validity test of the prediction algorithm was performed using the parallel CA-Markov model to simulate land-use changes in Hangzhou in 2013 and to analyze the relationship between the simulation results and the interpretation results of the remote-sensing images. The Kappa coefficients of construction land, natural-reserve land, and agricultural land were 0.86, 0.68, and 0.66, respectively, which demonstrates the validity of the parallel model. Hangzhou land-use changes in 2020 were predicted and analyzed. The results show that the central area of construction land is rapidly increasing due to a developed transportation system and is mainly transferred from agricultural land. Full article

Prime numbers are routinely used in a variety of applications, e.g., cryptography and hashing. A prime number can only be divided by one and the number itself. A semi-prime number is a product of two or more prime numbers (e.g., 5 × 3 = 15) and can only be formed by these numbers (e.g., 3 and 5). Exploiting this mathematical property allows schema-free encoding of geographical data in nominal or ordinal measurement scales for thematic maps. Schema-free encoding becomes increasingly important in the context of data variety. In this paper, I investigate the encoding of categorical thematic map data using prime numbers instead of a sequence of all natural numbers (1, 2, 3, 4, ..., n) as the category identifier. When prime numbers are multiplied, the result as a single value contains the information of more than one location category. I demonstrate how this encoding can be used on three use-cases, (1) a hierarchical legend for one theme (CORINE land use/land cover), (2) a combination of multiple topics in one theme (Köppen-Geiger climate classification), and (3) spatially overlapping regions (tree species distribution). Other applications in the field of geocomputation in general can also benefit from schema-free approaches with dynamic instead of handcrafted encoding of geodata. Full article

The history of modern maps in Japan begins with the Japan maps (called INŌ’s maps) prepared by Tadataka Inō after he thoroughly surveyed the whole of Japan around 200 years ago. The purpose of this study was to investigate the precision degree of INŌ’s Tokyo map by overlaying it with present maps and analyzing the map style (map projection, map scale, etc.). Specifically, we quantitatively examined the spatial distortion of INŌ’s maps through comparisons with the present map using GIS (geographic information system), a spatial analysis tool. Furthermore, by examining various factors that caused the positional gap and distortion of features, we explored the actual situation of surveying in that age from a geographical viewpoint. As a result of the analysis, a particular spatial regularity was confirmed in the positional gaps with the present map. We found that INŌ’s Tokyo map had considerably high precision. The causes of positional gaps from the present map were related not only to natural conditions, such as areas and land but also to social and cultural phenomena. Full article

The number of tidal flood events has been increasing in Indonesia in the last decade, especially along the north coast of Java. Hydrodynamic models in combination with Geographic Information System applications are used to assess the impact of high tide events upon the salt production in Cirebon, West Java. Two major flood events in June 2016 and May 2018 were selected for the simulation within inputs of tidal height records, national seamless digital elevation dataset of Indonesia (DEMNAS), Indonesian gridded national bathymetry (BATNAS), and wind data from OGIMET. We used a finite method on MIKE 21 to determine peak water levels, and validation for the velocity component using TPXO9 and Tidal Model Driver (TMD). The benchmark of the inundation is taken from the maximum water level of the simulation. This study utilized ArcGIS for the spatial analysis of tidal flood distribution upon solar salt production area, particularly where the tides are dominated by local factors. The results indicated that during the peak events in June 2016 and May 2018, about 83% to 84% of salt ponds were being inundated, respectively. The accurate identification of flooded areas also provided valuable information for tidal flood assessment of marginal agriculture in data-scarce region. Full article

Open areas, along with their non-forest vegetation, are often threatened by secondary succession, which causes deterioration of biodiversity and the habitat’s conservation status. The knowledge about characteristics and dynamics of the secondary succession process is very important in the context of management and proper planning of active protection of the Natura 2000 habitats. This paper presents research on the evaluation of the possibility of using selected methods of textural analysis to determine the spatial extent of trees and shrubs based on archival aerial photographs, and consequently on the investigation of the secondary succession process. The research was carried out on imagery from six different dates, from 1971 to 2015. The images differed from each other in spectral resolution (panchromatic, in natural colors, color infrared), in original spatial resolution, as well as in radiometric quality. Two methods of textural analysis were chosen for the analysis: Gray level co-occurrence matrix (GLCM) and granulometric analysis, in a number of variants, depending on the selected parameters of these transformations. The choice of methods has been challenged by their reliability and ease of implementation in practice. The accuracy assessment was carried out using the results of visual photo interpretation of orthophotomaps from particular years as reference data. As a result of the conducted analyses, significant efficacy of the analyzed methods has been proved, with granulometric analysis as the method of generally better suitability and greater stability. The obtained results show the impact of individual image features on the classification efficiency. They also show greater stability and reliability of texture analysis based on granulometric/morphological operations. Full article

With the rapid development of big data, numerous industries have turned their focus from information research and construction to big data technologies. Earth science and geographic information systems industries are highly information-intensive, and thus there is an urgent need to study and integrate big data technologies to improve their level of information. However, there is a large gap between existing big data and traditional geographic information technologies. Owing to certain characteristics, it is difficult to quickly and easily apply big data to geographic information technologies. Through the research, development, and application practices achieved in recent years, we have gradually developed a common geospatial big data solution. Based on the formation of a set of geospatial big data frameworks, a complete geospatial big data platform system called BiGeo was developed. Through the management and analysis of massive amounts of spatial data from Sichuan Province, China, the basic framework of this platform can be better utilized to meet our needs. This paper summarizes the design, implementation, and experimental experience of BiGeo, which provides a new type of solution to the research and construction of geospatial big data. Full article

The spatio-temporal relationship between tourism product similarity and spatial proximity has not been adequately studied empirically because of data and methodological limitations. New forms of data available at high temporal frequencies and low levels of spatial aggregation, together with large commercial data and expanding computational ability allow a variety of theories, old and new to be explored and evaluated more meticulously and systemically than has been possible hitherto. This study uses spatial visualization and data harvesting to synthesize a variety of data for exploring the evolution of hotel clusters and co-location synergies in US cities. The findings question the reliability of the current data to be used for identifying and analyzing the formation of tourist destination clusters and their dynamics. We conclude that synthesizing social media and large commercial data can generate a more robust database for research on tourism development and planning and improving opportunities for the examining spatial patterns of tourism activities. We also devise a protocol to combine ‘social media’ sources with big commercial sources for tourism development and planning, and eventually other sectors. Full article

Due to the high demand for cheap construction materials, clay-made brick manufacturing has become a thriving industry in Bangladesh, with manufacturing kilns heavily concentrated in the peripheries of larger cities and towns. These manufacturing sites, known as brickfields, operate using centuries-old technologies which expel dust, ash, black smoke and other pollutants into the atmosphere. This in turn impacts the air quality of cities and their surroundings and may also have broader impacts on health, the environment, and potentially contribute to global climate change. Using remotely sensed Landsat imagery, this study identifies brickfield locations and areal expansion between 1990 and 2015 in Dhaka, and employs spatial statistics methods including quadrat analysis and Ripley’s K-function to analyze the spatial variation of brickfield locations. Finally, using nearest neighbor distance as density functions, the distance between brickfield locations and six major geographical features (i.e., urban, rural settlement, wetland, river, highway, and local road) were estimated to investigate the threat posed by the presence of such polluting brickfields nearby urban, infrastructures and other natural areas. Results show significant expansion of brickfields both in number and clusters between 1990 and 2015 with brickfields increasing in number from 247 to 917 (total growth rate 271%) across the Dhaka urban center. The results also reveal that brickfield locations are spatially clustered: 78% of brickfields are located on major riverbanks and 40% of the total are located in ecologically sensitive wetlands surrounding Dhaka. Additionally, the average distance from the brick manufacturing plant to the nearest urban area decreased from 1500 m to 500 m over the study period. This research highlights the increasing threats to the environment, human health, and the sustainability of the megacity Dhaka from brickfield expansion in the immediate peripheral areas of its urban center. Findings and methods presented in this study can facilitate data-driven decision making by government officials and city planners to formulate strategies for improved brick production technologies and decreased environmental impacts for this urban region in Bangladesh. Full article

This paper presents a novel approach to the modelling and simulation of real estate transactions. The main purpose of the study was to develop the theoretical foundations for building simulation models of transaction locations and real estate prices. Pursuing this objective involved a spatial market analysis based on geostatistics to develop maps of the dynamics and spatial activity of the real estate market. The research was conducted by presenting the issue against the background of the literature of the subject and by conducting an experiment, which involved developing an original procedure of providing simulated market data. The study deals with the market for non-built-up land real estate with a residential function in the city of Olsztyn (Poland). The time range concerned the years 2004–2015. Information on 932 real estate transactions was adopted for the study. A set of additional information on virtual transactions was generated during the study; this information can supplement market data for markets of low activity or if there are information gaps. Geoinformation analyses were performed in order to determine new trends in price levels and spatial activity of a real estate market. Overall, this resulted in generating maps of simulated transaction densities, a map of simulated prices and a map of the probability of a specific price occurring. Full article

Public transport plays an important role in developing sustainable cities. A better understanding of how different public transit modes (bus, metro, and taxi) interact with each other will provide better sustainable strategies to transport and urban planners. However, most existing studies are either limited to small-scale surveys or focused on the identification of general interaction patterns during times of regular traffic. Transient demographic changes in a city (i.e., many people moving out and in) can lead to significant changes in such interaction patterns and provide a useful context for better investigating the changes in these patterns. Despite that, little has been done to explore how such interaction patterns change and how they are linked to the built environment from the perspective of transient demographic changes using urban big data. In this paper, the tap-in-tap-out smart card data of bus/metro and taxi GPS trajectory data before and after the Chinese Spring Festival in Shenzhen, China, are used to explore such interaction patterns. A time-series clustering method and an elasticity change index (ECI) are adopted to detect the changing transit mode patterns and the underlying dynamics. The findings indicate that the interactions between different transit modes vary over space and time and are competitive or complementary in different parts of the city. Both ordinary least-squares (OLS) and geographically weighted regression (GWR) models with built environment variables are used to reveal the impact of changes in different transit modes on ECIs and their linkage with the built environment. The results of this study will contribute to the planning and design of multi-modal transport services. Full article

Digital elevation models (DEMs) have become ubiquitous and remarkably effective in the field of earth sciences as a tool to characterize surface topography. All DEMs have a degree of inherent error and uncertainty that is propagated to subsequent models and analyses, which can lead to misinterpretation and inaccurate estimates. A new method was developed to estimate local DEM errors and implement corrections while quantifying the uncertainties of the implemented corrections. The method is based on the flexibility and ability to model complex problems with ensemble neural networks (ENNs). The method was developed to be applied to any DEM created from a corresponding set of elevation points (point cloud) and a set of ground truth measurements. The method was developed and tested using hyperspatial resolution terrestrial laser scanning (TLS) data (sub-centimeter point spacing) collected from a marsh site located along the southern portion of the Texas Gulf Coast, USA. ENNs improve the overall DEM accuracy in the study area by 68% for six model inputs and by 75% for 12 model inputs corresponding to root mean square errors (RMSEs) of 0.056 and 0.045 m, respectively. The 12-input model provides more accurate tolerance interval estimates, particularly for vegetated areas. The accuracy of the method is confirmed based on an independent data set. Although the method still underestimates the 95% tolerance interval, 8% below the 95% target, results show that it is able to quantify the spatial variability in uncertainties due to a relationship between vegetation/land cover and accuracy of the DEM for the study area. There are still opportunities and challenges in improving and confirming the applicability of this method for different study sites and data sets. Full article

The China–Pakistan Economic Corridor (CPEC) is the flagship project of the Belt and Road Initiative. At the end of the CPEC, the Gwadar port on the Arabian Sea is being built quickly, providing an important economical route for the flow of Central Asia’s natural resources to the world. Gwadar city is in a rapid urbanization process and will be developed as a modern, world-class port city in the near future. Therefore, monitoring the urbanization process of Gwadar at both high spatial and temporal resolution is vital for its urban planning, city ecosystem management, and the sustainable development of CPEC. The impervious surface percentage (ISP) is an essential quantitative indicator for the assessment of urban development. Through the integration of remote sensing images and ISP estimation models, ISP can be routinely and periodically estimated. However, due to clouds’ influence and spatial–temporal resolution trade-offs in sensor design, it is difficult to estimate the ISP with both high spatial resolution and dense temporal frequency from only one satellite sensor. In recent years, China has launched a series of Earth resource satellites, such as the HJ (Huangjing, which means environment in Chinese)-1A/B constellation, showing great application potential for rapid Earth surface mapping. This study employs the Random Forest (RF) method for a long-term and fine-scale ISP estimation and analysis of the city of Gwadar, based on the density in temporal and multi-source Chinese satellite images. In the method, high spatial resolution ISP reference data partially covering Gwadar city was first extracted from the 1–2 meter (m) GF (GaoFen, which means high spatial resolution in Chinese)-1/2 fused images. An RF retrieval model was then built based on the training samples extracted from ISP reference data and multi-temporal 30-m HJ-1A/B satellite images. Lastly, the model was used to generate the 30-m time series ISP from 2009 to 2017 for the whole city area based on the HJ-1A/B images. Results showed that the mean absolute error of the estimated ISP was 6.1–8.1% and that the root mean square error (RMSE) of the estimation results was 12.82–15.03%, indicating the consistently high performance of the model. This study highlights the feasibility and potential of using multi-source Chinese satellite images and an RF model to generate long-term ISP estimations for monitoring the urbanization process of the key node city in the CPEC. Full article

In the domain of computer vision, object recognition aims at detecting and classifying objects in data sets. Model-driven approaches are typically constrained through their focus on either a specific type of data, a context (indoor, outdoor) or a set of objects. Machine learning-based approaches are more flexible but also constrained as they need annotated data sets to train the learning process. That leads to problems when this data is not available through the specialty of the application field, like archaeology, for example. In order to overcome such constraints, we present a fully semantic-guided approach. The role of semantics is to express all relevant knowledge of the representation of the objects inside the data sets and of the algorithms which address this representation. In addition, the approach contains a learning stage since it adapts the processing according to the diversity of the objects and data characteristics. The semantic is expressed via an ontological model and uses standard web technology like SPARQL queries, providing great flexibility. The ontological model describes the object, the data and the algorithms. It allows the selection and execution of algorithms adapted to the data and objects dynamically. Similarly, processing results are dynamically classified and allow for enriching the ontological model using SPARQL construct queries. The semantic formulated through SPARQL also acts as a bridge between the knowledge contained within the ontological model and the processing branch, which executes algorithms. It provides the capability to adapt the sequence of algorithms to an individual state of the processing chain and makes the solution robust and flexible. The comparison of this approach with others on the same use case shows the efficiency and improvement this approach brings. Full article

Three-dimensional (3D) pipe network modeling plays an essential part in high performance-based smart city applications. Given that massive 3D pipe networks tend to be difficult to manage and to visualize, we propose in this study a hybrid framework for high-performance modeling of a 3D pipe network, including pipe network data model and high-performance modeling. The pipe network data model is devoted to three-dimensional pipe network construction based on network topology and building information models (BIMs). According to the topological relationships of the pipe point pipelines, the pipe network is decomposed into multiple pipe segment units. The high-performance modeling of 3D pipe network contains a spatial 3D model, the instantiation, adaptive rendering, and combination parallel computing. Spatial 3D model (S3M) is proposed for spatial data transmission, exchange, and visualization of massive and multi-source 3D spatial data. The combination parallel computing framework with GPU and OpenMP was developed to reduce the processing time for pipe networks. The results of the experiments showed that the hybrid framework achieves a high efficiency and the hardware resource occupation is reduced. Full article

Extracting features from crowd flow analysis has become an important research challenge due to its social cost and the impact of inadequate planning of high-quality services and security monitoring on the lives of citizens. This paper descriptively reviews and compares existing crowd analysis approaches based on different data sources. This survey provides the fundamentals of crowd analysis and considers three main approaches: crowd video analysis, crowd spatio-temporal analysis, and crowd social media analysis. The key research contributions in each approach are presented, and the most significant techniques and algorithms used to improve the precision of results that could be integrated into solutions to enhance the quality of services in a smart city are analyzed. Full article

The relationship between crime and income inequality is a complex and controversial issue. While there is some consensus that a relationship exists, the nature of it is still the subject of much debate. In this paper, this relationship is investigated in the context of urban geography and whether income inequality can explain the geography of crime within cities. This question is examined for the specific case of residential burglaries in the city of Belo Horizonte, Brazil, where I tested how much burglary rates are affected by local average household income and by local exposure to poverty, while I controlled for other variables relevant to criminological theory, such as land-use type, density and accessibility. Different scales were considered for testing the effect of exposure to poverty. This study reveals that, in Belo Horizonte, the rate of burglaries per single family house is significantly and positively related to income level, but a higher exposure to poverty has no significant independent effect on these rates at any scale tested. The rate of burglaries per apartment, on the other hand, is not significantly affected by either average household income or exposure to poverty. These results seem consistent with a description where burglaries follow a geographical distribution based on opportunity, rather than being a product of localized income disparity and higher exposure between different economic groups. Full article

Soil erosion processes are a type of geological hazard. They cause soil loss and sediment production, landscape dissection, and economic damage, which can, in the long term, result in land abandonment. Thus, identification of soil erosion processes is necessary for sustainable land management in an area. This study presents the potential of visual interpretation of high resolution LiDAR (light detection and ranging) imagery for direct and unambiguous identification and mapping of soil erosion processes, which was tested in the study area of the Vinodol Valley (64.57 km²), in Croatia. Eight LiDAR images were derived from the 1 m airborne LiDAR DTM (Digital Terrain Model) and were used to identify and map gully erosion, sheet erosion, and the combined effect of rill and sheet erosion, with the ultimate purpose to create a historical erosion inventory. The two-step procedure in a visual interpretation of LiDAR imagery was performed: preliminary and detailed. In the preliminary step, possibilities and limitations for unambiguous identification of the soil erosion processes were determined for representative portions of the study area, and the exclusive criteria for the accurate and precise manual delineation of different types of erosion phenomena were established. In the detailed step, the findings from the preliminary step were used to map the soil erosion phenomena in the entire studied area. Results determined the highest potential for direct identification and mapping of the gully erosion phenomena. A total of 236 gullies were identified and precisely delineated, although most of them were previously unknown, due to the lack of previous investigations on soil erosion processes in the study area. On the other hand, the used method was proven to be inapplicable for direct identification and accurate mapping of the sheet erosion. Sheet erosion, however, could have been indirectly identified on certain LiDAR imagery, based on recognition of colluvial deposits accumulated at the foot of the eroded slopes. Furthermore, the findings of this study present which of the used LiDAR imagery, and what features of the imagery used, are most effective for identification and mapping of different types of erosion processes. Full article

Quickly and efficiently monitoring soil nutrient contents using remote sensing technology is of great significance for farmland soil productivity, food security and sustainable agricultural development. Current research has been conducted to estimate and map soil nutrient contents in large areas using hyper-spectral techniques, however, it is difficult to obtain accurate estimates. In order to improve the estimation accuracy of soil nutrient contents, we introduced a GA-BPNN method, which combined a back propagation neural network (BPNN) with the genetic algorithm optimization (GA). This study was conducted in Guangdong, China, based on soil nutrient contents and hyperspectral data. The prediction accuracies from a partial least squares regression (PLSR), BPNN and GA-BPNN were compared using field observations. The results showed that (1) Among three methods, the GA-BPNN provided the most accurate estimates of soil total nitrogen (TN), total phosphorus (TP) and total potassium (TK) contents; (2) Compared with the BPNN models, the GA-BPNN models significantly improved the estimation accuracies of the soil nutrient contents by decreasing the relative root mean square error (RRMSE) values by 15.9%, 5.6% and 20.2% at the sample point level, and 20.1%, 16.5% and 47.1% at the regional scale for TN, TP and TK, respectively. This indicated that by optimizing the parameters of BPNN, the GA-BPNN provided greater potential to improving the estimation; and (3) Soil TK content could be more accurately mapped by the GA-BPNN method using HuanJing-1A Hyperspectral Imager (HJ-1A HSI) (manufacturer: China Aerospace Science and Technology Corporation; Beijing, China) data with a RRMSE value of 20.37% than the soil TN and TP with the RRMSE values of 40.41% and 34.71%, respectively. This implied that the GA-BPNN model provided the potential to map the soil TK content for the large area. The research results provided an important reference for high-accuracy prediction of soil nutrient contents. Full article