[go: up one dir, main page]

Next Issue
Volume 9, November
Previous Issue
Volume 9, September
 
 

ISPRS Int. J. Geo-Inf., Volume 9, Issue 10 (October 2020) – 55 articles

Cover Story (view full-size image): This work implements an approach for geospatial time-series processing in the analysis of deforestation. The method includes machine learning satellite image classification deployed on the cloud-based processing service Google Earth Engine. Twenty years of forest dynamics in the Amazonas were analysed, where classified maps were validated against high-resolution imagery through open-source software CollectEarth. An artificial neural network-based simulation, built thanks to the derived historical forest trends, was implemented to obtain a simulation of future forest loss. The processing approach is adaptable to different regions and periods of time, allowing also the monitoring of short/long-term implementation of forest preserving policies. View this paper
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Select all
Export citation of selected articles as:
19 pages, 3301 KiB  
Article
A Spatially Explicit Approach for Targeting Resource-Poor Smallholders to Improve Their Participation in Agribusiness: A Case of Nyando and Vihiga County in Western Kenya
by Mwehe Mathenge, Ben G. J. S. Sonneveld and Jacqueline E. W. Broerse
ISPRS Int. J. Geo-Inf. 2020, 9(10), 612; https://doi.org/10.3390/ijgi9100612 - 21 Oct 2020
Cited by 3 | Viewed by 2619
Abstract
The majority of smallholder farmers in Sub-Saharan Africa face myriad challenges to participating in agribusiness markets. However, how the spatially explicit factors interact to influence household decision choices at the local level is not well understood. This paper’s objective is to identify, map, [...] Read more.
The majority of smallholder farmers in Sub-Saharan Africa face myriad challenges to participating in agribusiness markets. However, how the spatially explicit factors interact to influence household decision choices at the local level is not well understood. This paper’s objective is to identify, map, and analyze spatial dependency and heterogeneity in factors that impede poor smallholders from participating in agribusiness markets. Using the researcher-administered survey questionnaires, we collected geo-referenced data from 392 households in Western Kenya. We used three spatial geostatistics methods in Geographic Information System to analyze data—Global Moran’s I, Cluster and Outliers Analysis, and geographically weighted regression. Results show that factors impeding smallholder farmers exhibited local spatial autocorrelation that was linked to the local context. We identified distinct local spatial clusters (hot spots and cold spots clusters) that were spatially and statistically significant. Results affirm that spatially explicit factors play a crucial role in influencing the farming decisions of smallholder households. The paper has demonstrated that geospatial analysis using geographically disaggregated data and methods could help in the identification of resource-poor households and neighborhoods. To improve poor smallholders’ participation in agribusiness, we recommend policymakers to design spatially targeted interventions that are embedded in the local context and informed by locally expressed needs. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Geographical location of Nyando and Vihiga study areas.</p>
Full article ">Figure 2
<p>Steps applied in geocoded household survey design; (<b>a</b>) distribution of randomized GIS points; (<b>b</b>) uploaded KML layers on ‘GPS Essential App’ in Android phone, and (<b>c</b>) Actual surveyed household GPS points.</p>
Full article ">Figure 3
<p>Graph showing the most optimal statistically significant peak z-scores of spatial clustering.</p>
Full article ">Figure 4
<p>Map of Nyando showing local spatial clusters with a higher concentration of poorer households (hot spots) and richer households (cold spots).</p>
Full article ">Figure 5
<p>Map of Vihiga showing local spatial clusters with a higher concentration of poorer households (hot spots) and richer households (cold spots).</p>
Full article ">Figure 6
<p>The two maps reveals that, in both Nyando and Vihiga areas, there were more poorer households (red dots) than richer households (blue dots) located within a I Kilometer buffer (crosshatched areas) from the tarmac road, main town, and rivers.</p>
Full article ">Figure 7
<p>Survey results showing households head education level vs agribusiness skills possession.</p>
Full article ">
19 pages, 4761 KiB  
Article
Evolution and Optimization of Urban Network Spatial Structure: A Case Study of Financial Enterprise Network in Yangtze River Delta, China
by Yizhen Zhang, Tao Wang, Agus Supriyadi, Kun Zhang and Zhi Tang
ISPRS Int. J. Geo-Inf. 2020, 9(10), 611; https://doi.org/10.3390/ijgi9100611 - 21 Oct 2020
Cited by 23 | Viewed by 3700
Abstract
The urban network is an important method of spatial optimization, and measuring the development level of the urban network is a prerequisite for spatial optimization. Combining geographic information system (GIS) spatial analysis, social network analysis, and multidimensional scaling models, we explored the evolution [...] Read more.
The urban network is an important method of spatial optimization, and measuring the development level of the urban network is a prerequisite for spatial optimization. Combining geographic information system (GIS) spatial analysis, social network analysis, and multidimensional scaling models, we explored the evolution of the urban network spatial structure in the Yangtze River Delta from 1990 to 2017 and proposed corresponding optimization measures. The results showed that the urban network spatial structure of the Yangtze River Delta has evolved from a single-center cluster with Shanghai as its core to a multi-center network with Shanghai as its core and Nanjing, Hangzhou, and Hefei as secondary cores. The density of the urban network has gradually expanded, but the strength of the connection between edge cities such as Chizhou, Suqian, and Quzhou and the core cities needs to be further improved. We found that the evolution of the urban network spatial structure has been driven by preferential attachment, path dependence, and differences in economic and industrial development. Finally, we propose optimizing the urban network spatial structure by strengthening the driving ability of the core cities, clarifying urban functions and development directions, and establishing a unified coordination mechanism. This paper enriches and deepens our understanding of the characteristics of the city network in the Yangtze River Delta, and provides a reference for the optimization of the urban network spatial structure. Full article
Show Figures

Figure 1

Figure 1
<p>Study area and location of the Yangtze River Delta in China.</p>
Full article ">Figure 2
<p>Location and number of financial companies in the Yangtze River Delta.</p>
Full article ">Figure 3
<p>Operation principle of the modularity model. We divided closely connected cities into different communities through a modularity model. (<b>a</b>) The blue nodes represent all the different cities. (<b>b</b>) After modularity optimization, some closely connected cities are selected, and cities of the same color are in the same community. (<b>c</b>) The cities of the same color are grouped into the same community after modularity division.</p>
Full article ">Figure 4
<p>Intercity network for the Yangtze River Delta from 1990 to 2017. The weight of the edge represents the R value (the total connection strength between cities) calculated above.</p>
Full article ">Figure 4 Cont.
<p>Intercity network for the Yangtze River Delta from 1990 to 2017. The weight of the edge represents the R value (the total connection strength between cities) calculated above.</p>
Full article ">Figure 5
<p>Network community division in the Yangtze River Delta from 1990 to 2017. The weight of the edge represents the R value (the total connection strength between cities) as calculated above. The same color is used to represent the cities in the same community, and the size of the node represents the degree centrality. (<b>a</b>–<b>d</b>) show the evolution of network community in 1990, 2000, 2008, and 2017, respectively.</p>
Full article ">Figure 6
<p>Evolution Mode of Urban Network Spatial Structure in Yangtze River Delta from 1990 to 2017. <b>Note:</b> <span class="html-fig-inline" id="ijgi-09-00611-i005"> <img alt="Ijgi 09 00611 i005" src="/ijgi/ijgi-09-00611/article_deploy/html/images/ijgi-09-00611-i005.png"/></span> represents a type I city, <span class="html-fig-inline" id="ijgi-09-00611-i006"> <img alt="Ijgi 09 00611 i006" src="/ijgi/ijgi-09-00611/article_deploy/html/images/ijgi-09-00611-i006.png"/></span> represents a type II city, <span class="html-fig-inline" id="ijgi-09-00611-i007"> <img alt="Ijgi 09 00611 i007" src="/ijgi/ijgi-09-00611/article_deploy/html/images/ijgi-09-00611-i007.png"/></span> represents a type III city, <span class="html-fig-inline" id="ijgi-09-00611-i008"> <img alt="Ijgi 09 00611 i008" src="/ijgi/ijgi-09-00611/article_deploy/html/images/ijgi-09-00611-i008.png"/></span> represents a type IV city. The green surrounding the red node refers to the research of Li et al., and to a certain extent represents the urban hinterland. [<a href="#B14-ijgi-09-00611" class="html-bibr">14</a>].</p>
Full article ">
17 pages, 10168 KiB  
Article
Articulated Trajectory Mapping for Reviewing Walking Tours
by Iori Sasaki, Masatoshi Arikawa and Akinori Takahashi
ISPRS Int. J. Geo-Inf. 2020, 9(10), 610; https://doi.org/10.3390/ijgi9100610 - 21 Oct 2020
Cited by 5 | Viewed by 2752
Abstract
This paper addresses how to enrich a map-based representation for reviewing walking tours with the features of trajectory mapping and tracing animation. Generally, a trajectory generated by raw GPS data can often be difficult to browse through on a map. To resolve this [...] Read more.
This paper addresses how to enrich a map-based representation for reviewing walking tours with the features of trajectory mapping and tracing animation. Generally, a trajectory generated by raw GPS data can often be difficult to browse through on a map. To resolve this issue, we first illustrated tangled trajectory lines, inaccurate indoor positioning, and unstable trajectory lines as problems encountered when mapping raw trajectory data. Then, we proposed a new framework that focuses on GPS horizontal accuracy to locate indoor location points and find stopping points on an accelerometer. We also applied a conventional line simplification algorithm to make the trajectory cleaner and then integrated the extracted points with the clean trajectory line. Furthermore, our experiments with some actual logs of walking tours demonstrated that articulated trajectory mapping, which comprises simplification and characterization methods, sufficiently reliable and effective for better reviewing experiences. The paper contributes to the research on cleaning up map-based displays and tracing animations of raw trajectory GPS data by using not only location data but also sensor data that smartphones can collect. Full article
(This article belongs to the Special Issue Recent Trends in Location Based Services and Science)
Show Figures

Figure 1

Figure 1
<p>Example displays of our mobile application during walking tours: (<b>a</b>) guide content is automatically displayed based on the user’s current location, and (<b>b</b>) users can find where they are on the illustrated maps.</p>
Full article ">Figure 2
<p>Example display of trajectory log data in the map-based reviewing mode of our application (without any data processing), which illustrates three kinds of problems (A, B, and C). The blue line is a GPS trajectory.</p>
Full article ">Figure 3
<p>Flow of generating articulated trajectory data from raw trajectory data as human-centered readable spatial content.</p>
Full article ">Figure 4
<p>Algorithm for extracting a set of indoor location points from the trajectory data using GPS horizontal accuracy values.</p>
Full article ">Figure 5
<p>Generating a center point representing an indoor location from multiple indoor points using a minimum bounding box.</p>
Full article ">Figure 6
<p>Examples of horizontal GPS errors while walking. The duration from the 10th point to the 58th point shows that a user had remained indoors.</p>
Full article ">Figure 7
<p>Improved algorithm for extracting a set of points of staying from trajectory data using GPS horizontal accuracy values (including tolerant distance buffering).</p>
Full article ">Figure 8
<p>Algorithm for extracting a set of stopping points from trajectory data using an acceleration sensor.</p>
Full article ">Figure 9
<p>Example of changes in acceleration values with a low-pass filter (a tester repeated ten walking steps and then stopped for five seconds three times. This test suggested that <math display="inline"><semantics> <mrow> <mi>T</mi> <msub> <mi>h</mi> <mn>2</mn> </msub> </mrow> </semantics></math> can be set to around 0.15 [G]).</p>
Full article ">Figure 10
<p>The process of the Douglas–Peucker algorithm.</p>
Full article ">Figure 11
<p>Illustrated map that examinees used in this experiment (Tomachi Area (traditional downtown area), Akita City, Japan).</p>
Full article ">Figure 12
<p>(<b>a</b>) Raw trajectory mapping, (<b>b</b>) articulated trajectory mapping (blue line: trajectory; blue round label with emoji: point of UGC, e.g., photos and textual notes; star label: point of staying).</p>
Full article ">Figure 13
<p>Results at which points of staying indoors were recognized with two different threshold settings, i.e., <math display="inline"><semantics> <mrow> <mi>T</mi> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10.0</mn> <mo> </mo> <mi>m</mi> </mrow> </semantics></math> (<b>left</b>) and <math display="inline"><semantics> <mrow> <mi>T</mi> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>20.0</mn> <mo> </mo> <mi>m</mi> </mrow> </semantics></math> (<b>right</b>).</p>
Full article ">Figure 14
<p>Examples of variations in the radius of certainty for the location as the horizontal GPS error while a user is staying indoors. <math display="inline"><semantics> <mrow> <mi>T</mi> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>10.0</mn> <mo> </mo> <mi>m</mi> </mrow> </semantics></math> is set as the threshold, and the graph does not display values higher than <math display="inline"><semantics> <mrow> <mn>40</mn> <mo> </mo> <mi>m</mi> </mrow> </semantics></math>. The duration from the 151st to the 313rd point actually represents the situation of a user being indoors. However, several parts shift below the threshold <math display="inline"><semantics> <mrow> <mi>T</mi> <msub> <mi>h</mi> <mn>1</mn> </msub> </mrow> </semantics></math> even if a user is moving indoors.</p>
Full article ">Figure 15
<p>Examples of effectiveness in articulated trajectory mapping. User experiments for testing articulated trajectory mapping reveal its ability to be reviewed and played back after walking tours.</p>
Full article ">Figure 16
<p>Experimental data of three cases with variations of acceleration sensors’ values. Light green areas indicate the durations corresponding to tangled trajectory lines when users were stopping outdoors, e.g., waiting for traffic signals.</p>
Full article ">Figure 17
<p>Examples of recognizing walking and stopping with the acceleration sensors’ threshold value set at <math display="inline"><semantics> <mrow> <mi>T</mi> <msub> <mi>h</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.1</mn> <mo> </mo> <mi>G</mi> </mrow> </semantics></math>. Blue dots show every 15-second point of a user’s trajectory data. Simplified trajectory data are effective in reviewing and playing back a user’s trajectory on the screen.</p>
Full article ">Figure 18
<p>Pedestrian subway.</p>
Full article ">Figure 19
<p>Typical architectural style of a Japanese sweets shop.</p>
Full article ">Figure 20
<p>Locations of the points of interest (POIs) that appears in this paper.</p>
Full article ">
23 pages, 7669 KiB  
Article
A 3D Geodatabase for Urban Underground Infrastructures: Implementation and Application to Groundwater Management in Milan Metropolitan Area
by Davide Sartirana, Marco Rotiroti, Chiara Zanotti, Tullia Bonomi, Letizia Fumagalli and Mattia De Amicis
ISPRS Int. J. Geo-Inf. 2020, 9(10), 609; https://doi.org/10.3390/ijgi9100609 - 21 Oct 2020
Cited by 8 | Viewed by 4415
Abstract
The recent rapid increase in urbanization has led to the inclusion of underground spaces in urban planning policies. Among the main subsurface resources, a strong interaction between underground infrastructures and groundwater has emerged in many urban areas in the last few decades. Thus, [...] Read more.
The recent rapid increase in urbanization has led to the inclusion of underground spaces in urban planning policies. Among the main subsurface resources, a strong interaction between underground infrastructures and groundwater has emerged in many urban areas in the last few decades. Thus, listing the underground infrastructures is necessary to structure an urban conceptual model for groundwater management needs. Starting from a municipal cartography (Open Data), thus making the procedure replicable, a GIS methodology was proposed to gather all the underground infrastructures into an updatable 3D geodatabase (GDB) for the metropolitan city of Milan (Northern Italy). The underground volumes occupied by three categories of infrastructures were included in the GDB: (a) private car parks, (b) public car parks and (c) subway lines and stations. The application of the GDB allowed estimating the volumes lying below groundwater table in four periods, detected as groundwater minimums or maximums from the piezometric trend reconstructions. Due to groundwater rising or local hydrogeological conditions, the shallowest, non-waterproofed underground infrastructures were flooded in some periods considered. This was evaluated in a specific pilot area and qualitatively confirmed by local press and photographic documentation reviews. The methodology emerged as efficient for urban planning, particularly for urban conceptual models and groundwater management plans definition. Full article
(This article belongs to the Special Issue Measuring, Mapping, Modeling, and Visualization of Cities)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Geographical setting of the study area; (<b>b</b>) the subway network of Milan; the location of the terminal stations of each line is provided. Line AA’ points to the location of the cross section represented in <a href="#ijgi-09-00609-f002" class="html-fig">Figure 2</a>. Satellite image Source: Geoportale Regione Lombardia.</p>
Full article ">Figure 2
<p>Hydrogeological schematic N–S cross section of the study area, showing the location of some subway stations.</p>
Full article ">Figure 3
<p>Flowchart of the proposed methodology.</p>
Full article ">Figure 4
<p>Schematization of the procedure to calculate the underground volume occupied by private buildings. (<b>a</b>) Identification of the dressing lines. (<b>b</b>) Digitization of the ramp polygon associated to the dressing lines. (<b>c</b>) Spatial analysis to associate the building with the ramp. Satellite image source: Geoportale Regione Lombardia.</p>
Full article ">Figure 5
<p>Location and volumes of all the infrastructural elements contained in the 3D geodatabase (GDB) (private and public car parks and subway lines).</p>
Full article ">Figure 6
<p>(<b>a</b>) Monitoring wells (MW) time series for the considered period (January 1990–December 2019). (<b>b</b>) Location of the MWs in the study area.</p>
Full article ">Figure 7
<p>Volumes lying below the groundwater (GW) table for the local minimum of Sep07 (<b>a</b>) and for the global maximum of Dec14 (<b>b</b>). Colour coding indicates percentages of the volumes below the groundwater table.</p>
Full article ">Figure 8
<p>Volumes and portions below the groundwater table (cyan expressed as percentages of the total volume) over time for: (<b>a</b>) Private Car Parks, (<b>b</b>) Public Car Parks, (<b>c</b>) Subway Line M1, (<b>d</b>) Subway Line M2, (<b>e</b>) Subway Line M3, (<b>f</b>) Underground Railway. <span class="html-italic">Y</span>-axis scale is the same for all the categories, except for Private Car Parks, where 10⁷ has been kept as an order of magnitude.</p>
Full article ">Figure 9
<p>(<b>a</b>) Geographical setting of the pilot area. Locations of volumes lying below the GW table for the maximum groundwater level condition of Dec14 are shown. (<b>b</b>) Underground infrastructures within the pilot area. (P) means public car park; (S) means station. Satellite image Source: Geoportale Regione Lombardia.</p>
Full article ">Figure 10
<p>Volumes lying below the GW table in (<b>a</b>) Jan90, (<b>b</b>) Dec02, (<b>c</b>) Sep07, (<b>d</b>) Dec14. Colour coding indicates percentages of the volumes below the groundwater table. Satellite image Source: Geoportale Regione Lombardia.</p>
Full article ">Figure 11
<p>(<b>a</b>) Three-dimensional surface reconstruction close to Sant’Agostino station. Sant’Agostino station is visible below the road network. (<b>b</b>) Three-dimensional underground reconstruction of Sant’Agostino station. Volumes below the GW table of Sant’Agostino station in (<b>c</b>) Jan90, (<b>d</b>) Dec02, (<b>e</b>) Sep07, (<b>f</b>) Dec14. Images were realized with ArcGIS Pro.</p>
Full article ">Figure 12
<p>(<b>a</b>) Newspaper article of “La Repubblica” (2 July 2013), dealing with flooding episodes in Sant’Agostino station. (<b>b</b>) Flooding evidence in Sant’Agostino station (8 September 2020). (Image credits to the authors). (<b>c</b>) Flooding evidence in Numa Pompilio public car park. (<b>d</b>) Absence of flooding evidence after waterproofing works in Numa Pompilio public car park. (<b>c</b>,<b>d</b>): images were provided by Rete Irene.</p>
Full article ">
14 pages, 4883 KiB  
Article
Accurate Road Marking Detection from Noisy Point Clouds Acquired by Low-Cost Mobile LiDAR Systems
by Ronghao Yang, Qitao Li, Junxiang Tan, Shaoda Li and Xinyu Chen
ISPRS Int. J. Geo-Inf. 2020, 9(10), 608; https://doi.org/10.3390/ijgi9100608 - 20 Oct 2020
Cited by 22 | Viewed by 4185
Abstract
Road markings that provide instructions for unmanned driving are important elements in high-precision maps. In road information collection technology, multi-beam mobile LiDAR scanning (MLS) is currently adopted instead of traditional mono-beam LiDAR scanning because of the advantages of low cost and multiple fields [...] Read more.
Road markings that provide instructions for unmanned driving are important elements in high-precision maps. In road information collection technology, multi-beam mobile LiDAR scanning (MLS) is currently adopted instead of traditional mono-beam LiDAR scanning because of the advantages of low cost and multiple fields of view for multi-beam laser scanners; however, the intensity information scanned by multi-beam systems is noisy and current methods designed for road marking detection from mono-beam point clouds are of low accuracy. This paper presents an accurate algorithm for detecting road markings from noisy point clouds, where most nonroad points are removed and the remaining points are organized into a set of consecutive pseudo-scan lines for parallel and/or online processing. The road surface is precisely extracted by a moving fitting window filter from each pseudo-scan line, and a marker edge detector combining an intensity gradient with an intensity statistics histogram is presented for road marking detection. Quantitative results indicate that the proposed method achieves average recall, precision, and Matthews correlation coefficient (MCC) levels of 90%, 95%, and 92%, respectively, showing excellent performance for road marking detection from multi-beam scanning point clouds. Full article
(This article belongs to the Special Issue Measuring, Mapping, Modeling, and Visualization of Cities)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>An overview of the proposed detection method.</p>
Full article ">Figure 2
<p>Pre-processing for road surface extraction: (<b>a</b>) The POS (positioning and orientation system) height, <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">H</mi> <mrow> <mi>POS</mi> </mrow> </msub> </mrow> </semantics></math>, from the POS centre to the road surface; (<b>b</b>) pseudo-scan lines rendered by random colours; (<b>c</b>) local coordinate system of points inside a pseudo-scan line.</p>
Full article ">Figure 3
<p>Two typical types of road surfaces: (<b>a</b>) Road surfaces I; (<b>b</b>) road surfaces II. The crosses are drainage channel points that have not been scanned.</p>
Full article ">Figure 4
<p>Intensity information of road markings: (<b>a</b>) Intensity values of a road marking in a pseudo-scan line before smoothing; (<b>b</b>) intensity values of a road marking in a pseudo-scan line after smoothing using a density-based adaptive window median filter; (<b>c</b>) intensity gradients of a road marking in a pseudo-scan line; (<b>d</b>) intensity statistics histogram of the smoothed road surface. Point #1 is the point that is not fully smoothed, point #2 is the left-hand edge point, and point #3 is the right-hand edge point. The intensity value at the red vertical line was chosen as the intensity threshold <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">I</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Overviews of the two datasets: (<b>a</b>) Dataset I rendered by the intensity attribute; (<b>b</b>) a detailed view of dataset I; (<b>c</b>) dataset II rendered by the intensity attribute; (<b>d</b>) a detailed view of dataset II.</p>
Full article ">Figure 6
<p>Road markings extracted by the marker edge constraint detector (MECD) method with different thresholds: (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">P</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = 2; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">P</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = 4; (<b>c</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">P</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = 6; (<b>d</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">N</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = −6; (<b>e</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">N</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = −4; (<b>f</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">N</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = −2; (<b>g</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">I</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = 8; (<b>h</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">I</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math>= 10; (<b>i</b>) <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">I</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> = 12. <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">P</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">N</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">I</mi> <mrow> <mi>th</mi> </mrow> </msub> </mrow> </semantics></math> denote the positive gradient threshold, the negative gradient threshold, and the intensity threshold.</p>
Full article ">Figure 7
<p>Extraction results of the two datasets: (<b>a</b>,<b>b</b>) Extracted road surfaces; (<b>c</b>,<b>d</b>) extracted road markings; (<b>e</b>,<b>f</b>) refined road markings.</p>
Full article ">Figure 8
<p>Road markings extracted by Yu’s method [<a href="#B16-ijgi-09-00608" class="html-bibr">16</a>]: (<b>a</b>) From dataset I; (<b>b</b>) from dataset II.</p>
Full article ">
20 pages, 2913 KiB  
Article
Privacy-Aware Visualization of Volunteered Geographic Information (VGI) to Analyze Spatial Activity: A Benchmark Implementation
by Alexander Dunkel, Marc Löchner and Dirk Burghardt
ISPRS Int. J. Geo-Inf. 2020, 9(10), 607; https://doi.org/10.3390/ijgi9100607 - 20 Oct 2020
Cited by 17 | Viewed by 3775
Abstract
Through volunteering data, people can help assess information on various aspects of their surrounding environment. Particularly in natural resource management, Volunteered Geographic Information (VGI) is increasingly recognized as a significant resource, for example, supporting visitation pattern analysis to evaluate collective values and improve [...] Read more.
Through volunteering data, people can help assess information on various aspects of their surrounding environment. Particularly in natural resource management, Volunteered Geographic Information (VGI) is increasingly recognized as a significant resource, for example, supporting visitation pattern analysis to evaluate collective values and improve natural well-being. In recent years, however, user privacy has become an increasingly important consideration. Potential conflicts often emerge from the fact that VGI can be re-used in contexts not originally considered by volunteers. Addressing these privacy conflicts is particularly problematic in natural resource management, where visualizations are often explorative, with multifaceted and sometimes initially unknown sets of analysis outcomes. In this paper, we present an integrated and component-based approach to privacy-aware visualization of VGI, specifically suited for application to natural resource management. As a key component, HyperLogLog (HLL)—a data abstraction format—is used to allow estimation of results, instead of more accurate measurements. While HLL alone cannot preserve privacy, it can be combined with existing approaches to improve privacy while, at the same time, maintaining some flexibility of analysis. Together, these components make it possible to gradually reduce privacy risks for volunteers at various steps of the analytical process. A specific use case demonstration is provided, based on a global, publicly-available dataset that contains 100 million photos shared by 581,099 users under Creative Commons licenses. Both the data processing pipeline and resulting dataset are made available, allowing transparent benchmarking of the privacy–utility tradeoffs. Full article
(This article belongs to the Special Issue Volunteered Geographic Information and Citizen Science)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Illustration of the system model and the two cases of possible adversaries discussed in this work.</p>
Full article ">Figure 2
<p>Transformation steps applied to a single character string, such as a user ID, for generating a HyperLogLog (HLL) set, and the final estimation of cardinality (Example values were generated with real data, but different values may be produced based on various parameter settings).</p>
Full article ">Figure 3
<p>Percentage of global spatial outlier volume (k = 1) in the Yahoo Flickr Creative Commons 100 Million (YFCC100M) dataset, for decreasing precision levels (GeoHash) and different metrics used in this paper (to reproduce this graphic, see <a href="#app1-ijgi-09-00607" class="html-app">Supplementary Materials, S5</a>).</p>
Full article ">Figure 4
<p>Comparison of automatic classification of raw and HLL user days for Europe (100 km grid).</p>
Full article ">Figure 5
<p>Screenshot of map for user counts per 100 km grid bin, allowing interactive comparison of estimated values (HLL) and exact counts (raw) (see <a href="#ijgi-09-00607-f0A1" class="html-fig">Figure A1</a> for a static, worldwide view of the map, and <a href="#app1-ijgi-09-00607" class="html-app">Supplementary Materials S8</a> for the interactive version).</p>
Full article ">Figure 6
<p>Estimated post count with a reduced grid size of 50 km for Europe.</p>
Full article ">Figure 7
<p>Analyzing spatial relationships with HLL intersection, based on incremental union of user sets from benchmark data (100 km-grid) for France, Germany and the United Kingdom (<b>left</b>). The Venn Diagram (<b>right</b>) shows estimation of common user counts for different groups, and the percentage of error compared to raw data. The same graphic, generated for 50 km grid size, is available in <a href="#ijgi-09-00607-f0A2" class="html-fig">Figure A2</a>.</p>
Full article ">Figure 8
<p>Alex case study, evaluation of scenario “Sandy”.</p>
Full article ">Figure 9
<p>Alex case study, evaluation of scenario “Robert”.</p>
Full article ">Figure A1
<p>Worldwide map of estimated user counts (YFCC) per 100 km grid bin.</p>
Full article ">Figure A2
<p><a href="#ijgi-09-00607-f007" class="html-fig">Figure 7</a> generated with 50 km grid size parameter and corresponding error rates.</p>
Full article ">
26 pages, 3605 KiB  
Article
Towards Deriving Freight Traffic Measures from Truck Movement Data for State Road Planning: A Proposed System Framework
by Ahmed Karam, Thorbjørn M. Illemann, Kristian Hegner Reinau, Goran Vuk and Christian O. Hansen
ISPRS Int. J. Geo-Inf. 2020, 9(10), 606; https://doi.org/10.3390/ijgi9100606 - 14 Oct 2020
Cited by 6 | Viewed by 3533
Abstract
To make the right decisions on investments, operations, and policies in the public road sector, decision makers need knowledge about traffic measures of trucks, such as average travel time and the frequency of trips among geographical zones. Private logistics companies daily collect a [...] Read more.
To make the right decisions on investments, operations, and policies in the public road sector, decision makers need knowledge about traffic measures of trucks, such as average travel time and the frequency of trips among geographical zones. Private logistics companies daily collect a large amount of freight global positioning system (GPS) and shipment data. Processing such data can provide public decision makers with detailed freight traffic measures, which are necessary for making different planning decisions. The present paper proposes a system framework to be used in the research project “A new system for sharing data between logistics companies and public infrastructure authorities: improving infrastructure while maintaining competitive advantage”. Previous studies ignored the fact that the primary step for delivering valuable and usable data processing systems is to consider the final user’s needs when developing the system framework. Unlike existing studies, this paper develops the system framework through applying a user-centred design approach combining three main steps. The first step is to identify the specific traffic measures that satisfy the public decision makers’ planning needs. The second step aims to identify the different types of freight data required as inputs to the data processing system, while the third step illustrates the procedures needed to process the shared freight data. To do so, the current work employs methods of literature review and users’ need identification in applying a user-centralized approach. In addition, we develop a systematic assessment of the coverage and sufficiency of the currently acquired data. Finally, we illustrate the detailed functionality of the data processing system and provide an application case to illustrate its procedures. Full article
(This article belongs to the Special Issue Recent Trends in Location Based Services and Science)
Show Figures

Figure 1

Figure 1
<p>Number of research articles per year.</p>
Full article ">Figure 2
<p>The proposed requirements gathering method.</p>
Full article ">Figure 3
<p>Distributions of the identified articles in relation to the three categories (<b>a</b>), and time analysis of the articles in each category (<b>b</b>).</p>
Full article ">Figure 4
<p>Distributions of the identified articles in categories of infrastructure-planning (<b>a</b>), freight movement regulations (<b>b</b>), and freight data analyses (<b>c</b>).</p>
Full article ">Figure 5
<p>The percentage of the GPS data currently available at the main regions of Denmark (Adapted from reference [<a href="#B92-ijgi-09-00606" class="html-bibr">92</a>]).</p>
Full article ">Figure 6
<p>Zoom in on raw GPS-data from trucks, a snapshot from the developed system.</p>
Full article ">Figure 7
<p>Overall structure of the developed data processing system.</p>
Full article ">Figure 8
<p>Average values of indices, data availability index and data reliability index for selected road sections.</p>
Full article ">Figure 9
<p>Location of loop detector and the GPS data.</p>
Full article ">Figure 10
<p>The amount of GPS trucks whose data are available at the system in each hour of 14th February 2019.</p>
Full article ">Figure 11
<p>Comparison of hourly mean speed from station and GPS data.</p>
Full article ">
17 pages, 10280 KiB  
Article
Is Crowdsourcing a Reliable Method for Mass Data Acquisition? The Case of COVID-19 Spread in Greece During Spring 2020
by Varvara Antoniou, Emmanuel Vassilakis and Maria Hatzaki
ISPRS Int. J. Geo-Inf. 2020, 9(10), 605; https://doi.org/10.3390/ijgi9100605 - 14 Oct 2020
Cited by 17 | Viewed by 4284
Abstract
We present a GIS-based crowdsourcing application that was launched soon after the first COVID-19 cases had been recorded in Greece, motivated by the need for fast, location-wise data acquisition regarding COVID-19 disease spread during spring 2020, due to limited testing. A single question [...] Read more.
We present a GIS-based crowdsourcing application that was launched soon after the first COVID-19 cases had been recorded in Greece, motivated by the need for fast, location-wise data acquisition regarding COVID-19 disease spread during spring 2020, due to limited testing. A single question was posted through a web App, to which the anonymous participants subjectively answered whether or not they had experienced any COVID-19 disease symptoms. Our main goal was to locate geographical areas with increased number of people feeling the symptoms and to determine any temporal changes in the statistics of the survey entries. It was found that the application was rapidly disseminated to the entire Greek territory via social media, having, thus, a great public reception. The higher percentages of participants experiencing symptoms coincided geographically with the highly populated urban areas, having also increased numbers of confirmed cases, while temporal variations were detected that accorded with the restrictions of activities. This application demonstrates that health systems can use crowdsourcing applications that assure anonymity, as an alternative to tracing apps, to identify possible hot spots and to reach and warn the public within a short time interval, increasing at the same time their situational awareness. However, a continuous reminder for participation should be scheduled. Full article
(This article belongs to the Collection Spatial Components of COVID-19 Pandemic)
Show Figures

Figure 1

Figure 1
<p>Interface of the survey in Greek and English (<a href="https://arcg.is/1yiXG4" target="_blank">https://arcg.is/1yiXG4</a>).</p>
Full article ">Figure 2
<p>Dashboard showing the survey responders location, in real time (map on the left) and aggregation analysis for “Not Sure” (map in the middle) and “Yes” (map on the right) per Regional Unit performed every 3 hours (<a href="https://arcg.is/00Kanf" target="_blank">https://arcg.is/00Kanf</a>). Indicators below the maps show the total number of responders, Not Sure and Yes, respectively. Screenshot was taken on 18/03/2020 12:00 (EET).</p>
Full article ">Figure 3
<p>Spatial distribution of survey entries, the majority of which are located in the Regions of Attica and Central Macedonia. Black lines indicate geographical boundaries of the Regions. Green points indicate No, light green Not Sure and red points Yes answers; black circles indicate bigger cities.</p>
Full article ">Figure 4
<p>Temporal evolution for the Greek territory of (<b>a</b>) the number of participants per day with the number of Yes and Not Sure annotated; (<b>b</b>) ratios of the three different answers per day (left axis), compared to the confirmed cases (number/day; right axis); (<b>c</b>) the confirmed cases and deaths per day since the recording of the first confirmed case in Greece. Critical restriction of activities is also indicated with vertical solid or dashed lines.</p>
Full article ">Figure 5
<p>Spatial distribution of survey entries for (<b>a</b>) the Region of Attica and (<b>b</b>) the Athens metropolitan area. Temporal evolution for the Region of Attica of (<b>c</b>) the number of participants per day with the number of Yes and Not Sure annotated, (<b>d</b>) ratios of the three different answers per day (left axis), compared to the confirmed cases (number/day; right axis), (<b>e</b>) the confirmed cases and deaths per day in Attica since the recording of the first confirmed case in Greece. Critical restriction of activities is also indicated with vertical solid or dashed lines.</p>
Full article ">Figure 6
<p>Spatial distribution of survey entries for (<b>a</b>) the Region of Central Macedonia and (<b>b</b>) the greater area of Thessaloniki. Temporal evolution for the Region of Central Macedonia of (<b>c</b>) the number of participants per day with the number of Yes and Not Sure annotated; (<b>d</b>) ratios of the three different answers per day (left axis), compared to the confirmed cases (number/day; right axis); (<b>e</b>) the confirmed cases and deaths per day in Central Macedonia since the recording of the first confirmed case in Greece. Critical restriction of activities is also indicated with vertical solid or dashed lines.</p>
Full article ">Figure 7
<p>Hot spot analysis for Athens metropolitan area and surrounding towns on daily basis for the first eight days of the survey.</p>
Full article ">
25 pages, 15167 KiB  
Article
Reconstruction of Lost Cultural Heritage Sites and Landscapes: Context of Ancient Objects in Time and Space
by Lukáš Brůha, Josef Laštovička, Tomáš Palatý, Eva Štefanová and Přemysl Štych
ISPRS Int. J. Geo-Inf. 2020, 9(10), 604; https://doi.org/10.3390/ijgi9100604 - 14 Oct 2020
Cited by 16 | Viewed by 4464
Abstract
Diachronic studies play a key role in the research and documentation of cultural heritage and its changes, ranging from architectural fragments to landscape. Regarding the reconstructions of lost cultural heritage sites, the determination of landscape conditions in the reconstructed era goes frequently unheeded. [...] Read more.
Diachronic studies play a key role in the research and documentation of cultural heritage and its changes, ranging from architectural fragments to landscape. Regarding the reconstructions of lost cultural heritage sites, the determination of landscape conditions in the reconstructed era goes frequently unheeded. Often, only ruins and detached archeological artefacts remain of the built heritage. Placing them correctly within the reconstructed building complex is of similar importance as placing the lost monument in the context of the landscape at that time. The proposed method harmonizes highly heterogeneous sources to provide such a context. The solution includes the fusion of referential terrain models of different levels of detail (LODs) as well as the fusion of diverse 3D data sources for the reconstruction of the built heritage. Although the combined modeling of large landscapes and small 3D objects of a high detail results in very large datasets, we present a feasible solution, whose data structure is suitable for Geographic Information Systems (GIS) analyses of landscapes and also provides a smooth and clear 3D visualization and inspection of detailed features. The results are demonstrated in the case study of the island monastery, the vanished medieval town of Sekanka, and the surrounding landscape, which is located in Czechia and was the subject of intensive changes over time. Full article
Show Figures

Figure 1

Figure 1
<p>The contemporary state from an aerial view (archeologickyatlas.cz) with marked locations of the lost monastery (red circle), medieval town (green circle), and Church of St. Kilián (blue circle). Source: The Archeology atlas of Czechia (2020; <a href="http://www.archeologickyatlas.cz/cs/lokace/davle_pz_mestecko_sekanka" target="_blank">http://www.archeologickyatlas.cz/cs/lokace/davle_pz_mestecko_sekanka</a>) and Open Street Map (2020; geo: 49.8812, 14.3908).</p>
Full article ">Figure 2
<p>The solution’s workflow summary starting with the conceptual analysis to the employment of technology for distribution of variable thematic data layers.</p>
Full article ">Figure 3
<p>Various historical sources employed: (<b>a</b>) Presentation of artefacts in the National Museum in Prague; (<b>b</b>) The Romanesque pillar of the National Museum collection; (<b>c</b>) The hexagonal tiles of so-called Vyšehrad type; (<b>d</b>) The tombstone of abbot Heřman. Photos: Petr Kříž and Tomáš Palatý (2017).</p>
Full article ">Figure 4
<p>Examples of early maps and land-survey plans originating from archeological research: (<b>a</b>) The Second military mapping; (<b>b</b>) The cadaster mapping (Franziszeische Kataster); (<b>c</b>) The plans differentiating the structures of monastery’s Romanesque and Gothic eras; (<b>d</b>) The plan of the lost town Sekanka (Romanesque era); archaeological drawing of the monastery in Romanesque era (<b>e</b>) and in Gothic era (<b>f</b>). Source: 2nd Military Survey, Section No. 9/2, Austrian State Archive/Military Archive, Vienna (<b>a</b>); State Administration of Land Surveying and Cadaster (<b>b</b>); Miroslav Richter (1982; (<b>c</b>,<b>d</b>)); František Stehlík (1947; original) and Jan Heřman (2009; redrawn with new notes by archeologists; (<b>e</b>,<b>f</b>)).</p>
Full article ">Figure 4 Cont.
<p>Examples of early maps and land-survey plans originating from archeological research: (<b>a</b>) The Second military mapping; (<b>b</b>) The cadaster mapping (Franziszeische Kataster); (<b>c</b>) The plans differentiating the structures of monastery’s Romanesque and Gothic eras; (<b>d</b>) The plan of the lost town Sekanka (Romanesque era); archaeological drawing of the monastery in Romanesque era (<b>e</b>) and in Gothic era (<b>f</b>). Source: 2nd Military Survey, Section No. 9/2, Austrian State Archive/Military Archive, Vienna (<b>a</b>); State Administration of Land Surveying and Cadaster (<b>b</b>); Miroslav Richter (1982; (<b>c</b>,<b>d</b>)); František Stehlík (1947; original) and Jan Heřman (2009; redrawn with new notes by archeologists; (<b>e</b>,<b>f</b>)).</p>
Full article ">Figure 5
<p>(<b>a</b>) The depiction of the triangulation of two datasets with highly different resolution resulting in a sharp divide and unnatural turn. (<b>b</b>) A 3D view on the boundary line and triangles with too small interior angles. (<b>c</b>) An illustration of the different DTMs geometries in the transition zone—low detail, high detail, and the transition zone illustrated with different colors. (<b>d</b>) The figure on the bottom right depicts a 3D view on the resulting geometry of the terrain. The resulting triangulation consists of triangles with far greater interior angles (more equilateral), which follows the objective of the Delaunay triangulation and therefore adheres more closely to the actual terrain course.</p>
Full article ">Figure 6
<p>Photographs of the former island and its merge with mainland from the beginning of the 20th century from northeast (<b>a</b>) and north (<b>b</b>). Source: Josef Dvořák (<a href="http://www.dvorak-davle.cz" target="_blank">http://www.dvorak-davle.cz</a>).</p>
Full article ">Figure 7
<p>The dense cloud representations of the tile processed in Agisoft Metashape at different levels of detail (LODs): (<b>a</b>) 12,782,806 points (<b>b</b>) 2,848,116 points, and (<b>c</b>) 319,603 points.</p>
Full article ">Figure 8
<p>Visualization of the reconstructed middle age town Sekanka.</p>
Full article ">Figure 9
<p>Illustration of the incursion of the Brandenburgers.</p>
Full article ">Figure 10
<p>Visualization of the monastery in the Romanesque era in the summer season.</p>
Full article ">Figure 11
<p>Visualization of Church of St. Kilián on the left bank in the Romanesque era (<b>a</b>) and in the Gothic era (<b>b</b>).</p>
Full article ">Figure 12
<p>Visualization of the Gothic era in the winter season from southwest (<b>a</b>) and southeast (<b>b</b>).</p>
Full article ">Figure 13
<p>Gothic basilica interiors including paintings of the period and tombs. Own work (2020).</p>
Full article ">Figure 14
<p>Interiors of the Gothic quadrature nearby the cloister including paintings of the era and the floor formed from tiles. Own work (2020).</p>
Full article ">Figure 15
<p>The gothic tiled floor created using adapted textures only (<b>a</b>) and using the processed terrestrial laser scanning (TLS) model (<b>b</b>) for comparison.</p>
Full article ">Figure 16
<p>The virtual museum application with hotspots (<b>a</b>) that can be studied in detail in its own pop-up viewer (<b>b</b>).</p>
Full article ">Figure 17
<p>The navigation between historical periods (<b>a</b>) the Romanesque era view and (<b>b</b>) the Gothic era view from the same spot and with the same view frustrum.</p>
Full article ">Figure A1
<p>Illustration of the incursion of the Brandenburgers (the Romanesque era).</p>
Full article ">Figure A2
<p>Visualization of Church of St. Kilián and the monastery in the Romanesque era.</p>
Full article ">Figure A3
<p>Visualization of Church of St. Kilián and the monastery in the Gothic era.</p>
Full article ">Figure A4
<p>Visualization of the monastery in the Gothic era.</p>
Full article ">Figure A5
<p>Visualization of Church of St. Kilián and the monastery in the Gothic era.</p>
Full article ">
32 pages, 7257 KiB  
Article
Detailed Streetspace Modelling for Multiple Applications: Discussions on the Proposed CityGML 3.0 Transportation Model
by Christof Beil, Roland Ruhdorfer, Theresa Coduro and Thomas H. Kolbe
ISPRS Int. J. Geo-Inf. 2020, 9(10), 603; https://doi.org/10.3390/ijgi9100603 - 13 Oct 2020
Cited by 30 | Viewed by 6385
Abstract
In the context of smart cities and digital twins, three-dimensional semantic city models are increasingly used for the analyses of large urban areas. While the representation of buildings, terrain, and vegetation has become standard for most city models, detailed spatio-semantic representations of streetspace [...] Read more.
In the context of smart cities and digital twins, three-dimensional semantic city models are increasingly used for the analyses of large urban areas. While the representation of buildings, terrain, and vegetation has become standard for most city models, detailed spatio-semantic representations of streetspace have played a minor role so far. This is now changing (1) because of data availability, and (2) because recent and emerging applications require having detailed data about the streetspace. The upcoming version 3.0 of the international standard CityGML provides a substantially updated data model regarding the transportation infrastructure, including the representation of the streetspace. However, there already exist a number of other standards and data formats dealing with the representation and exchange of streetspace data. Thus, based on an extensive literature review of potential applications as well as discussions and collaborations with relevant stakeholders, seven key modelling aspects of detailed streetspace models are identified. This allows a structured discussion of representational capabilities of the proposed CityGML3.0 Transportation Model with respect to these aspects and in comparison to the other standards. Subsequently, it is shown that CityGML3.0 meets most of these aspects and that streetspace models can be derived from various data sources and for different cities. Models generated compliant to the CityGML standard are immediately usable for a number of applications. This is demonstrated for some applications, such as land use management, solar potential analyses, and traffic and pedestrian simulations. Full article
(This article belongs to the Special Issue The Applications of 3D-City Models in Urban Studies)
Show Figures

Figure 1

Figure 1
<p>Comparison of different representation types for streetspace modelling [<a href="#B3-ijgi-09-00603" class="html-bibr">3</a>].</p>
Full article ">Figure 2
<p>CityGML 3.0 Transportation Model as presented by the OGC CityGML SWG [<a href="#B69-ijgi-09-00603" class="html-bibr">69</a>]. New classes compared to CityGML2.0 are marked with orange borders.</p>
Full article ">Figure 3
<p>(<b>a</b>) Street network segmentation into Sections (orange) and Intersections (blue); (<b>b</b>) Different possibilities to define Intersection areas.</p>
Full article ">Figure 4
<p>Street Section in different levels of granularity (areal and linear representation).</p>
Full article ">Figure 5
<p>Space concept for Transportation objects in CityGML 3.0 [<a href="#B5-ijgi-09-00603" class="html-bibr">5</a>].</p>
Full article ">Figure 6
<p>Potential inconsistency between linear and areal street representations (<b>left</b>) and proposed solution (<b>right</b>).</p>
Full article ">Figure 7
<p>(<b>a</b>) Comparison of linear representations in granularity = area and granularity = lane; (<b>b</b>) Linear representations of predecessor/successor relations.</p>
Full article ">Figure 8
<p>Predecessor/Successor relations.</p>
Full article ">Figure 9
<p>Affiliations between CityGML 2.0 sub- and top-level features.</p>
Full article ">Figure 10
<p>(<b>a</b>) Object breaklines on DEM, (<b>b</b>) Triangulated surface, (<b>c</b>) Result.</p>
Full article ">Figure 11
<p>Detailed street representation.</p>
Full article ">Figure 12
<p>Streetspace model of Melbourne (near intersection of Flemington Rd and Elizabeth St).</p>
Full article ">Figure 13
<p>Streetspace model including city furniture objects generated from mobile mapping data and traffic simulation.</p>
Full article ">Figure 14
<p>Complex intersection including multi-functional TrafficAreas and Markings.</p>
Full article ">Figure 15
<p>Demo of CityGML 3.0 concepts for an area around TU Munich.</p>
Full article ">Figure 16
<p>Traffic simulation visualization [<a href="#B22-ijgi-09-00603" class="html-bibr">22</a>].</p>
Full article ">Figure 17
<p>Visualization of a global irradiation estimation (kWh/a) for buildings and street objects Maximum Sky View Factor: (1): 0.532, (2): 0.958, (3): 0.254.</p>
Full article ">
19 pages, 4739 KiB  
Article
Concept and Evaluation of Heating Demand Prediction Based on 3D City Models and the CityGML Energy ADE—Case Study Helsinki
by Maxim Rossknecht and Enni Airaksinen
ISPRS Int. J. Geo-Inf. 2020, 9(10), 602; https://doi.org/10.3390/ijgi9100602 - 12 Oct 2020
Cited by 35 | Viewed by 4984
Abstract
This work presents a concept for heating demand and resulting CO2 emissions prediction based on a 3D city model in CityGML format in various scenarios under the consideration of a changing climate. In the case study of Helsinki, the Helsinki Energy and [...] Read more.
This work presents a concept for heating demand and resulting CO2 emissions prediction based on a 3D city model in CityGML format in various scenarios under the consideration of a changing climate. In the case study of Helsinki, the Helsinki Energy and Climate Atlas, that provides detailed information for individual buildings conducting the heating demand, is integrated into the 3D city model using the CityGML Energy Application Domain Extension (Energy ADE) to provide energy-relevant information based on a standardized data model stored in a CityGML database, called 3DCityDB. The simulation environment SimStadt is extended to retrieve the information stored within the Energy ADE schema, use it during simulations, and write simulation results back to the 3DCityDB. Due to climate change, a heating demand reduction of 4% per decade is predicted. By 2035, a reduction of 0.7 TWh is calculated in the normal and of 1.5 TWh in the advanced refurbishment scenario. Including the proposed improvements of the district heating network, heating CO2 emissions are predicted to be reduced by up to 82% by 2035 compared to 1990. The City of Helsinki’s assumed heating demand reduction through the modernization of 2.0 TWh/a by 2035 is not achieved with a 3% refurbishment rate. Furthermore, the reduction of CO2 emissions is mainly achieved through lower CO2 emission factors of the district heating network in Helsinki. Full article
(This article belongs to the Special Issue The Applications of 3D-City Models in Urban Studies)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Building energy predictions approach classification [<a href="#B3-ijgi-09-00602" class="html-bibr">3</a>].</p>
Full article ">Figure 2
<p>Overview of the proposed concept.</p>
Full article ">Figure 3
<p>Principles of the Energy Application Domain Extension (ADE) integration approaches.</p>
Full article ">Figure 4
<p>Simplified CityGML building with used Energy ADE information.</p>
Full article ">Figure 5
<p>SQL statement to insert the floor area to the ng_floorarea table in a 3DCityDB with Energy ADE extension.</p>
Full article ">Figure 6
<p>Expected development of the Helen district heating network CO<math display="inline"><semantics> <msub> <mrow/> <mn>2</mn> </msub> </semantics></math> emissions [<a href="#B25-ijgi-09-00602" class="html-bibr">25</a>].</p>
Full article ">Figure 7
<p>Screenshot of the 3D web visualization, showing the buildings colorized by simulated space heating demand.</p>
Full article ">Figure 8
<p>Distribution of deviations using the usage zone and heated area information for the simulation.</p>
Full article ">Figure 9
<p>Possible incorrectly mapped heat demand of a building.</p>
Full article ">Figure 10
<p>Statistical heating demand (gray) and predicted heating demand (orange) in Helsinki, assuming a 1% refurbishment rate in the BAU scenario.</p>
Full article ">Figure 11
<p>Statistical (gray) and predicted (orange) CO<math display="inline"><semantics> <msub> <mrow/> <mn>2</mn> </msub> </semantics></math> emissions in Helsinki with a 1% refurbishment rate, using the 2019 CO<math display="inline"><semantics> <msub> <mrow/> <mn>2</mn> </msub> </semantics></math> emission factor for district heating.</p>
Full article ">Figure 12
<p>Statistical heating demand (gray) and predicted heating demand (orange) in Helsinki, assuming a 3% refurbishment rate for the actual building stock.</p>
Full article ">Figure 13
<p>Statistical (gray) and predicted (orange) CO<math display="inline"><semantics> <msub> <mrow/> <mn>2</mn> </msub> </semantics></math> emissions in Helsinki with a 3% refurbishment rate and CO<math display="inline"><semantics> <msub> <mrow/> <mn>2</mn> </msub> </semantics></math> efficient district heating network.</p>
Full article ">
15 pages, 2358 KiB  
Article
Semantic Segmentation of Remote-Sensing Imagery Using Heterogeneous Big Data: International Society for Photogrammetry and Remote Sensing Potsdam and Cityscape Datasets
by Ahram Song and Yongil Kim
ISPRS Int. J. Geo-Inf. 2020, 9(10), 601; https://doi.org/10.3390/ijgi9100601 - 12 Oct 2020
Cited by 19 | Viewed by 5232
Abstract
Although semantic segmentation of remote-sensing (RS) images using deep-learning networks has demonstrated its effectiveness recently, compared with natural-image datasets, obtaining RS images under the same conditions to construct data labels is difficult. Indeed, small datasets limit the effective learning of deep-learning networks. To [...] Read more.
Although semantic segmentation of remote-sensing (RS) images using deep-learning networks has demonstrated its effectiveness recently, compared with natural-image datasets, obtaining RS images under the same conditions to construct data labels is difficult. Indeed, small datasets limit the effective learning of deep-learning networks. To address this problem, we propose a combined U-net model that is trained using a combined weighted loss function and can handle heterogeneous datasets. The network consists of encoder and decoder blocks. The convolutional layers that form the encoder blocks are shared with the heterogeneous datasets, and the decoder blocks are assigned separate training weights. Herein, the International Society for Photogrammetry and Remote Sensing (ISPRS) Potsdam and Cityscape datasets are used as the RS and natural-image datasets, respectively. When the layers are shared, only visible bands of the ISPRS Potsdam data are used. Experimental results show that when same-sized heterogeneous datasets are used, the semantic segmentation accuracy of the Potsdam data obtained using our proposed method is lower than that obtained using only the Potsdam data (four bands) with other methods, such as SegNet, DeepLab-V3+, and the simplified version of U-net. However, the segmentation accuracy of the Potsdam images is improved when the larger Cityscape dataset is used. The combined U-net model can effectively train heterogeneous datasets and overcome the insufficient training data problem in the context of RS-image datasets. Furthermore, it is expected that the proposed method can not only be applied to segmentation tasks of aerial images but also to tasks with various purposes of using big heterogeneous datasets. Full article
(This article belongs to the Special Issue Geospatial Big Data and Machine Learning Opportunities and Prospects)
Show Figures

Figure 1

Figure 1
<p>Architecture of the combined U-net model. “Conv2D” denotes the two-dimensional (2D) convolutional layers, and “Conv2DTranspose” denotes a transposed 2D convolutional layer. “Concatenate” denotes a concatenated layer.</p>
Full article ">Figure 2
<p>Framework of the proposed method. The combined U-net model shares encoding blocks and is trained using the combined weighted loss.</p>
Full article ">Figure 3
<p>Example of the International Society for Photogrammetry and Remote Sensing (ISPRS) dataset; the patch number is 2–10: (<b>a</b>) RGB-image, (<b>b</b>) labeling image, (<b>c</b>) enlarged RGB-image, and (<b>d</b>) data labels.</p>
Full article ">Figure 4
<p>Examples from the Cityscape dataset: (<b>a</b>,<b>d</b>) RGB-images, (<b>b</b>,<b>e</b>) original label images, and (<b>c</b>,<b>f</b>) redefined label images.</p>
Full article ">Figure 5
<p>Final label classes of the Potsdam and Cityscape datasets.</p>
Full article ">Figure 6
<p>Learning graph of the overall accuracy (OA) for each epoch. (<b>a</b>) SegNet with original bands of Potsdam dataset, (<b>b</b>) DeepLab-V3+ with original bands of Potsdam dataset, (<b>c</b>) simplified U-net with RGB bands of Potsdam dataset, (<b>d</b>) simplified U-net with original bands of Potsdam dataset, (<b>e</b>) Case 1 wherein training proceeded using both the Potsdam and Cityscape datasets (same sizes) by the combined U-net method. (<b>f</b>) Case 2 wherein training proceeded using the Potsdam and Cityscape datasets by the combined U-net method; however, the Cityscape dataset was about twice as large as the Potsdam dataset.</p>
Full article ">Figure 7
<p>Example of input Potsdam RGB, label, and the resulting semantic segmentation images for the three cases. (<b>a</b>,<b>f</b>,<b>k</b>,<b>q</b>) are the input Potsdam images, (<b>b</b>,<b>g</b>,<b>l</b>,<b>r</b>) are the label images, (<b>c</b>,<b>h</b>,<b>m</b>,<b>s</b>) are the resulting images generated in simplified U-net, (<b>d</b>,<b>i</b>,<b>n</b>,<b>t</b>) are the resulting images generated in Case 1, and (<b>e</b>,<b>j</b>,<b>o</b>,<b>u</b>) are the resulting images generated in Case 2.</p>
Full article ">Figure 8
<p>OA and loss of the test set when (<b>a</b>) the number of Cityscape datasets was varied and when (<b>b</b>) the weight values of the Potsdam and Cityscape were changed</p>
Full article ">
19 pages, 5103 KiB  
Article
Assessing Quality of Life Inequalities. A Geographical Approach
by Antigoni Faka
ISPRS Int. J. Geo-Inf. 2020, 9(10), 600; https://doi.org/10.3390/ijgi9100600 - 12 Oct 2020
Cited by 23 | Viewed by 5086
Abstract
This study proposes an integrated methodology for evaluating and mapping quality of life (QoL) and the quality of a place as residence area, at local level. The QoL assessment was based on the development of composite criteria, using geographical variables that evaluate QoL, [...] Read more.
This study proposes an integrated methodology for evaluating and mapping quality of life (QoL) and the quality of a place as residence area, at local level. The QoL assessment was based on the development of composite criteria, using geographical variables that evaluate QoL, and geographic information systems. The composite criteria are related to the natural and the socioeconomic environment, the housing conditions, the infrastructure and services, and the cultural and recreational facilities. Each criterion was evaluated by a set of variables and each variable was weighted based on the residents’ preferences and the analytical hierarchy process. The criteria were also weighted and combined to assess overall QoL. The methodology was implemented in the Municipality of Katerini, Greece, and QoL mapping led to the zoning of the study area and the identification of areas with low and high QoL. The results revealed the highest level of overall QoL in three out of twenty-nine communities, which provide better housing conditions and access to public services and infrastructures, combining also qualitative natural environment, whereas five mountainous and remote communities scored the lowest level. Mapping QoL may support decision making strategies that target to improve human well-being, increase QoL levels and upgrade living conditions. Full article
(This article belongs to the Special Issue GIS-Based Analysis for Quality of Life and Environmental Monitoring)
Show Figures

Figure 1

Figure 1
<p>Municipality of Katerini and local administrative units (LAUs). LAUs are described by ID in <a href="#ijgi-09-00600-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 2
<p>Flowchart of the proposed methodology.</p>
Full article ">Figure 3
<p>Spatial layers of factors per criterion. (<b>a</b>) spatial layers of public services and infrastructures, (<b>b</b>) spatial layers of cultural and recreational facilities, (<b>c</b>) spatial layers of natural environment.</p>
Full article ">Figure 4
<p>Variables of socioeconomic environment and housing conditions (%), LAUs of the Katerini Municipality. (<b>a</b>) percentage of unemployment, (<b>b</b>) percentage of higher educated population, (<b>c</b>) percentage of employment in upper occupational groups, (<b>d</b>) percentage of employment in lower occupational groups, (<b>e</b>) percentage of households without basic facilities, (<b>f</b>) percentage of households in detached houses, (<b>g</b>) percentage of households in newly built units, (<b>h</b>) housing space per person in m<sup>2</sup>.</p>
Full article ">Figure 5
<p>Variables of public services and infrastructures, cultural and recreational facilities, and natural environment, LAUs of Katerini Municipality. (<b>a</b>) accessibility to medical services in minutes, (<b>b</b>) accessibility to educational facilities in minutes, (<b>c</b>) accessibility to sport facilities in minutes, (<b>d</b>) accessibility to urban amenities in minutes, (<b>e</b>) internet network density (Km/sq. Km), (<b>f</b>) accessibility to leisure facilities in minutes, (<b>g</b>) accessibility to cultural facilities in minutes, (<b>h</b>) accessibility to coastline in minutes, (<b>i</b>) percentage of urban areas, (<b>j</b>) percentage of forestlands.</p>
Full article ">Figure 6
<p>Maps of the QoL criteria, LAUs of Katerini Municipality. (<b>a</b>) socioeconomic environment, (<b>b</b>) housing conditions, (<b>c</b>) public services and infrastructures, (<b>d</b>) cultural and recreational facilities, (<b>e</b>) natural environment.</p>
Full article ">Figure 7
<p>Overall QoL levels, LAUs of the Katerini Municipality.</p>
Full article ">
14 pages, 2222 KiB  
Article
Spatial Patterns of Childhood Obesity Prevalence in Relation to Socioeconomic Factors across England
by Yeran Sun, Xuke Hu, Ying Huang and Ting On Chan
ISPRS Int. J. Geo-Inf. 2020, 9(10), 599; https://doi.org/10.3390/ijgi9100599 - 11 Oct 2020
Cited by 6 | Viewed by 5254
Abstract
To examine to what extent spatial inequalities in childhood obesity are attributable to spatial inequalities in socioeconomic characteristics across a country, we aimed to investigate the spatial associations of socioeconomic characteristics and childhood obesity. We first explored spatial patterns of childhood obesity prevalence, [...] Read more.
To examine to what extent spatial inequalities in childhood obesity are attributable to spatial inequalities in socioeconomic characteristics across a country, we aimed to investigate the spatial associations of socioeconomic characteristics and childhood obesity. We first explored spatial patterns of childhood obesity prevalence, and subsequently investigated the spatial associations of socioeconomic factors and childhood obesity prevalence across England by selecting and estimating appropriate spatial regression models. As the data used are geospatial data, we used two newly developed specifications of spatial regression models to investigate the spatial association of socioeconomic factors and childhood obesity prevalence. As a result, among the two newly developed specifications of spatial regression models, the fast random effects specification of eigenvector spatial filtering (FRES-ESF) model appears to outperform the matrix exponential spatial specification of spatial autoregressive (MESS-SAR) model. Empirical results indicate that positive spatial dependence is found to exist in childhood obesity prevalence across England; and that socioeconomic factors are significantly associated with childhood obesity prevalence across England. In England, children living in areas with lower socioeconomic status are at higher risk of obesity. This study suggests effectively reducing spatial inequalities in socioeconomic status will plays a vital role in mitigating spatial inequalities in childhood obesity prevalence. Full article
Show Figures

Figure 1

Figure 1
<p>Middle Super Output Area (MSOA)-level percentage of obese children in Year 6 (age 10–11 years) across England, 2013/14 to 2015/16.</p>
Full article ">Figure 2
<p>MSOA-level spatial distribution of explanatory variables across England.</p>
Full article ">Figure 2 Cont.
<p>MSOA-level spatial distribution of explanatory variables across England.</p>
Full article ">Figure 3
<p>Global Moran scatterplot of MSOA-level “percentage of obese children”.</p>
Full article ">Figure 4
<p>Clusters and outliers of MSOA-level “percentage of obese children” across the regions of England.</p>
Full article ">
18 pages, 4004 KiB  
Article
Simulating Large-Scale 3D Cadastral Dataset Using Procedural Modelling
by Jernej Tekavec, Anka Lisec and Eugénio Rodrigues
ISPRS Int. J. Geo-Inf. 2020, 9(10), 598; https://doi.org/10.3390/ijgi9100598 - 11 Oct 2020
Cited by 8 | Viewed by 2984
Abstract
Geospatial data and information within contemporary land administration systems are fundamental to manage the territory adequately. 3D land administration systems, often addressed as 3D cadastre, promise several benefits, particularly in managing today’s complex built environment, but these are currently still non-existent in their [...] Read more.
Geospatial data and information within contemporary land administration systems are fundamental to manage the territory adequately. 3D land administration systems, often addressed as 3D cadastre, promise several benefits, particularly in managing today’s complex built environment, but these are currently still non-existent in their full capacity. The development of any complex information and administration system, such as a land administration system, is time-consuming and costly, particularly during the phase of evaluation and testing. In this regard, the process of implementing such systems may benefit from using synthetic data. In this study, the method for simulating the 3D cadastral dataset is presented and discussed. The dataset is generated using a procedural modelling method, referenced to real cadastral data for the Slovenian territory and stored in a spatial database management system (DBMS) that supports storage of 3D spatial data. Spatial queries, related to 3D cadastral data management, are used to evaluate the database performance and storage characteristics, and 3D visualisation options. The results of the study show that the method is feasible for the simulation of large-scale 3D cadastral datasets. Using the developed spatial queries and their performance analysis, we demonstrate the importance of the simulated dataset for developing efficient 3D cadastral data management processes. Full article
Show Figures

Figure 1

Figure 1
<p>Study concept framework.</p>
Full article ">Figure 2
<p>Evolutionary Program for the Space Allocation Problem (EPSAP) algorithm workflow (adapted from ref. [<a href="#B54-ijgi-09-00598" class="html-bibr">54</a>]).</p>
Full article ">Figure 3
<p>Simulated two-storey building model comprising of wall, door and window surfaces.</p>
Full article ">Figure 4
<p>Linking Land Administration Domain Model (LADM) classes (green) to the data model of simulated buildings (white).</p>
Full article ">Figure 5
<p>Solid geometries created from the predefined surfaces.</p>
Full article ">Figure 6
<p>(<b>a</b>) Original placement of generated 3D building models; (<b>b</b>) the same models after their placement to the local coordinate system origin.</p>
Full article ">Figure 7
<p>3D room geometries (random colours), 2D building outline (yellow) and extruded 3D building exterior geometry (transparent blue).</p>
Full article ">Figure 8
<p>Database management system (DBMS) storage of the transformed procedurally generated 3D building models.</p>
Full article ">Figure 9
<p>DBMS storage for the simulated 3D cadastral dataset.</p>
Full article ">Figure 10
<p>Examples of the generated buildings ranging from the single-storey family house to multi-storey mixed-use apartment buildings.</p>
Full article ">Figure 11
<p>3D visualisation of 3D visualisation of georeferenced simulated residential 3D building in Google Earth.</p>
Full article ">
25 pages, 16093 KiB  
Article
Spatial and Temporal Patterns in Volunteer Data Contribution Activities: A Case Study of eBird
by Guiming Zhang
ISPRS Int. J. Geo-Inf. 2020, 9(10), 597; https://doi.org/10.3390/ijgi9100597 - 11 Oct 2020
Cited by 32 | Viewed by 5355
Abstract
Volunteered geographic information (VGI) has great potential to reveal spatial and temporal dynamics of geographic phenomena. However, a variety of potential biases in VGI are recognized, many of which root from volunteer data contribution activities. Examining patterns in volunteer data contribution activities helps [...] Read more.
Volunteered geographic information (VGI) has great potential to reveal spatial and temporal dynamics of geographic phenomena. However, a variety of potential biases in VGI are recognized, many of which root from volunteer data contribution activities. Examining patterns in volunteer data contribution activities helps understand the biases. Using eBird as a case study, this study investigates spatial and temporal patterns in data contribution activities of eBird contributors. eBird sampling efforts are biased in space and time. Most sampling efforts are concentrated in areas of denser populations and/or better accessibility, with the most intensively sampled areas being in proximity to big cities in developed regions of the world. Reported bird species are also spatially biased towards areas where more sampling efforts occur. Temporally, eBird sampling efforts and reported bird species are increasing over the years, with significant monthly fluctuations and notably more data reported on weekends. Such trends are driven by the expansion of eBird and characteristics of bird species and observers. The fitness of use of VGI should be assessed in the context of applications by examining spatial, temporal and other biases. Action may need to be taken to account for the biases so that robust inferences can be made from VGI observations. Full article
(This article belongs to the Special Issue Citizen Science and Geospatial Capacity Building)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of eBird sampling locations in 2019.</p>
Full article ">Figure 2
<p>Covariates used for modeling the spatial pattern in sampling efforts of eBird contributors.</p>
Full article ">Figure 3
<p>Frequency distribution of sampling locations and of the world on the covariates.</p>
Full article ">Figure 4
<p>The number of cumulative eBird sampling events as of 31 December 2019 mapped over 0.25° latitude × 0.25° longitude grid cells. Intervals were determined loosely following quartile classification method.</p>
Full article ">Figure 5
<p>Number of sampling events in each year (2002–2019).</p>
Full article ">Figure 6
<p>Percentage of sampling events in each month relative to the yearly total number of events.</p>
Full article ">Figure 7
<p>Average number of sampling events on each day of the week over the years.</p>
Full article ">Figure 8
<p>Frequency distribution of sampling events by the number of reported species.</p>
Full article ">Figure 9
<p>The cumulative number of observers as of 31 December 2019 mapped over 0.25° latitude × 0.25° longitude grid cells. Intervals were determined loosely following quartile classification method.</p>
Full article ">Figure 10
<p>Number of active observers in each year (2002–2019).</p>
Full article ">Figure 11
<p>Percentage of observers in each month relative to the yearly total number of observers.</p>
Full article ">Figure 12
<p>Average number of active observers on each day of the week over the years.</p>
Full article ">Figure 13
<p>Frequency distribution of observers by year of first observation (<b>upper left</b>), year of last observation (<b>upper right</b>), number of years between the first and last observations (<b>lower left</b>), and number of active dates (<b>lower right</b>).</p>
Full article ">Figure 14
<p>Frequency distribution of observers by number of reported species (<b>left</b>), number of sampling events (<b>center</b>) and number of sampling locations (<b>right</b>).</p>
Full article ">Figure 15
<p>The cumulative number of species reported to eBird as of 31 December 2019 mapped over 0.25° latitude × 0.25° longitude grid cells. Intervals were determined loosely following quartile classification method.</p>
Full article ">Figure 16
<p>Number of species reported in each year (2002–2019).</p>
Full article ">Figure 17
<p>Percentage of species reported in each month relative to the yearly total number of species.</p>
Full article ">Figure 18
<p>Average number of species reported on each day of the week over the years.</p>
Full article ">Figure 19
<p>Distribution of the number of species by the number of observers who reported the same bird species (<b>left</b>) and by the number of sampling events in which the same species was reported (<b>right</b>).</p>
Full article ">Figure 20
<p>Jackknife test of variable importance to the Maxent model based on test AUC.</p>
Full article ">Figure 21
<p>Map of sampling probability of eBird contributors modeled and predicted using Maxent.</p>
Full article ">Figure 22
<p>Number of species (<b>top</b>), incomplete checklists (<b>center</b>), and unapproved observations (<b>bottom</b>) per observer across the days of the week.</p>
Full article ">
16 pages, 42937 KiB  
Article
A Smooth Transition Algorithm for Adjacent Panoramic Viewpoints Using Matched Delaunay Triangular Patches
by Pengcheng Zhao, Qingwu Hu, Zhixiong Tang and Mingyao Ai
ISPRS Int. J. Geo-Inf. 2020, 9(10), 596; https://doi.org/10.3390/ijgi9100596 - 10 Oct 2020
Cited by 5 | Viewed by 3956
Abstract
The unnatural panoramic image transition between two adjacent viewpoints reduces the immersion and interactive experiences of 360° panoramic walkthrough systems. In this paper, a dynamic panoramic image rendering and smooth transition algorithm for adjacent viewpoints is proposed. First, the feature points of adjacent [...] Read more.
The unnatural panoramic image transition between two adjacent viewpoints reduces the immersion and interactive experiences of 360° panoramic walkthrough systems. In this paper, a dynamic panoramic image rendering and smooth transition algorithm for adjacent viewpoints is proposed. First, the feature points of adjacent view images are extracted, a robust matching algorithm is used to establish adjacent point pairs, and the matching triangles are formed by using the homonymous points. Then, a dynamic transition model is formed by the simultaneous linear transitions of shape and texture for each control triangle. Finally, the smooth transition between adjacent viewpoints is implemented by overlaying the dynamic transition model with the 360° panoramic walkthrough scene. Experimental results show that this method has obvious advantages in visual representation with distinct visual movement. It can realize the smooth transition between two indoor panoramic stations with arbitrary station spacing, and its execution efficiency is up to 50 frames per second. It effectively enhances the interactivity and immersion of 360° panoramic walkthrough systems. Full article
(This article belongs to the Special Issue Measuring, Mapping, Modeling, and Visualization of Cities)
Show Figures

Figure 1

Figure 1
<p>Principle schematic diagram of panorama transition guided by matched triangular patches.</p>
Full article ">Figure 2
<p>Proposed technique workflow for smooth transition between adjacent viewpoints.</p>
Full article ">Figure 3
<p>Generation flow chart of triangular patches for controlling transition. The green arrow represents the direction of execution of the local projection, and the red arrow represents the direction of execution of the back projection. The picture at the bottom right shows the extracted control triangular patches.</p>
Full article ">Figure 4
<p>Demonstration of the transition from one tri angle to another.</p>
Full article ">Figure 5
<p>Middle textures generated by the smooth transition algorithm. The panorama scene transitions from the current site to the next with the process ranging from 0.0 to 1.0. The control patches transition from a small piece on the sphere to a hemi-spherical surface as the viewpoint moves forward.</p>
Full article ">Figure 6
<p>The dynamic panoramic ball display strategy. The outer layer is a panoramic sphere for roaming, and the inner layer is the morph layer for smooth site transitions.</p>
Full article ">Figure 7
<p>Visual performance of the 360° panoramic walkthrough system with smooth transitions between adjacent viewpoints. The panorama scene transitions from the current site to the next with the process ranging from 0.0 to 1.0.</p>
Full article ">
26 pages, 13177 KiB  
Article
Hierarchical Instance Recognition of Individual Roadside Trees in Environmentally Complex Urban Areas from UAV Laser Scanning Point Clouds
by Yongjun Wang, Tengping Jiang, Jing Liu, Xiaorui Li and Chong Liang
ISPRS Int. J. Geo-Inf. 2020, 9(10), 595; https://doi.org/10.3390/ijgi9100595 - 10 Oct 2020
Cited by 22 | Viewed by 3477
Abstract
Individual tree segmentation is essential for many applications in city management and urban ecology. Light Detection and Ranging (LiDAR) system acquires accurate point clouds in a fast and environmentally-friendly manner, which enables single tree detection. However, the large number of object categories and [...] Read more.
Individual tree segmentation is essential for many applications in city management and urban ecology. Light Detection and Ranging (LiDAR) system acquires accurate point clouds in a fast and environmentally-friendly manner, which enables single tree detection. However, the large number of object categories and occlusion from nearby objects in complex environment pose great challenges in urban tree inventory, resulting in omission or commission errors. Therefore, this paper addresses these challenges and increases the accuracy of individual tree segmentation by proposing an automated method for instance recognition urban roadside trees. The proposed algorithm was implemented of unmanned aerial vehicles laser scanning (UAV-LS) data. First, an improved filtering algorithm was developed to identify ground and non-ground points. Second, we extracted tree-like objects via labeling on non-ground points using a deep learning model with a few smaller modifications. Unlike only concentrating on the global features in previous method, the proposed method revises a pointwise semantic learning network to capture both the global and local information at multiple scales, significantly avoiding the information loss in local neighborhoods and reducing useless convolutional computations. Afterwards, the semantic representation is fed into a graph-structured optimization model, which obtains globally optimal classification results by constructing a weighted indirect graph and solving the optimization problem with graph-cuts. The segmented tree points were extracted and consolidated through a series of operations, and they were finally recognized by combining graph embedding learning with a structure-aware loss function and a supervoxel-based normalized cut segmentation method. Experimental results on two public datasets demonstrated that our framework achieved better performance in terms of classification accuracy and recognition ratio of tree. Full article
(This article belongs to the Special Issue Advanced Research Based on Multi-Dimensional Point Cloud Analysis)
Show Figures

Figure 1

Figure 1
<p>Pipeline of the proposed method.</p>
Full article ">Figure 2
<p>The example of the filtering result. (<b>a</b>) Raw point cloud; (<b>b</b>) Ground points; (<b>c</b>) Non-ground points after ground removal.</p>
Full article ">Figure 3
<p>The proposed point-wise semantic learning network architecture.</p>
Full article ">Figure 4
<p>Illustration of the structure of the point feature extraction module.</p>
Full article ">Figure 5
<p>Illustration of the structure of the local feature extraction module.</p>
Full article ">Figure 6
<p>Illustration of the whole network architecture for individual tree segmentation. N is the number of points. F is the dimension of the backbone (submanifold convolutional network) output. E is the dimension of the instance embedding. The over-segmentation algorithm is used to cluster the instance embeddings during the inference.</p>
Full article ">Figure 7
<p>Illustration of the aggregator using attention-based KNN.</p>
Full article ">Figure 8
<p>Illustration of the updator part. The skip connection is used to concatenate the output of the aggregator and the input embedding together. Finally, a fully connected (FC) layer follows to update and get the refined output embedding.</p>
Full article ">Figure 9
<p>Classification results of the DFC 3D dataset. (<b>a</b>) the classification result with PointNet; (<b>b</b>) the classification result with the proposed model; (<b>c</b>) the ground truth.</p>
Full article ">Figure 10
<p>Detailed classification results of a certain selected area. (<b>a</b>) the initial classification result; (<b>b</b>) the smoothed classification result; (<b>c</b>) the ground truth.</p>
Full article ">Figure 11
<p>Detailed results from selected area 1 of DFC 3D dataset. (<b>a</b>) the original point clouds; (<b>b</b>) the classification result; (<b>c</b>) the roadside trees extraction; (<b>d</b>) the roadside trees segmentation.</p>
Full article ">Figure 12
<p>Detailed results from selected area 2 of DFC 3D dataset. (<b>a</b>) the original point clouds; (<b>b</b>) the classification result; (<b>c</b>) the roadside trees extraction; (<b>d</b>) the roadside trees segmentation.</p>
Full article ">Figure 13
<p>Detailed results from selected area 1 of DALES dataset. (<b>a</b>) the original point clouds; (<b>b</b>) the classification result; (<b>c</b>) the roadside trees extraction; (<b>d</b>) the roadside trees segmentation.</p>
Full article ">Figure 14
<p>Detailed results from selected area 2 of DALES dataset. (<b>a</b>) the original point clouds; (<b>b</b>) the classification result; (<b>c</b>) the roadside trees extraction; (<b>d</b>) the roadside trees segmentation.</p>
Full article ">Figure 15
<p>Details of the roadside trees recognition outcomes: (<b>a</b>) small trees; (<b>b</b>) incomplete tree.</p>
Full article ">
15 pages, 4010 KiB  
Technical Note
PolySimp: A Tool for Polygon Simplification Based on the Underlying Scaling Hierarchy
by Ding Ma, Zhigang Zhao, Ye Zheng, Renzhong Guo and Wei Zhu
ISPRS Int. J. Geo-Inf. 2020, 9(10), 594; https://doi.org/10.3390/ijgi9100594 - 10 Oct 2020
Cited by 6 | Viewed by 3556
Abstract
Map generalization is a process of reducing the contents of a map or data to properly show a geographic feature(s) at a smaller extent. Over the past few years, the fractal way of thinking has emerged as a new paradigm for map generalization. [...] Read more.
Map generalization is a process of reducing the contents of a map or data to properly show a geographic feature(s) at a smaller extent. Over the past few years, the fractal way of thinking has emerged as a new paradigm for map generalization. A geographic feature can be deemed as a fractal given the perspective of scaling, as its rough, irregular, and unsmooth shape inherently holds a striking scaling hierarchy of far more small elements than large ones. The pattern of far more small things than large ones is a de facto heavy tailed distribution. In this paper, we apply the scaling hierarchy for map generalization to polygonal features. To do this, we firstly revisit the scaling hierarchy of a classic fractal: the Koch Snowflake. We then review previous work that used the Douglas–Peuker algorithm, which identifies characteristic points on a line to derive three types of measures that are long-tailed distributed: the baseline length (d), the perpendicular distance to the baseline (x), and the area formed by x and d (area). More importantly, we extend the usage of the three measures to other most popular cartographical generalization methods; i.e., the bend simplify method, Visvalingam–Whyatt method, and hierarchical decomposition method, each of which decomposes any polygon into a set of bends, triangles, or convex hulls as basic geometric units for simplification. The different levels of details of the polygon can then be derived by recursively selecting the head part of geometric units and omitting the tail part using head/tail breaks, which is a new classification scheme for data with a heavy-tailed distribution. Since there are currently few tools with which to readily conduct the polygon simplification from such a fractal perspective, we have developed PolySimp, a tool that integrates the mentioned four algorithms for polygon simplification based on its underlying scaling hierarchy. The British coastline was selected to demonstrate the tool’s usefulness. The developed tool can be expected to showcase the applicability of fractal way of thinking and contribute to the development of map generalization. Full article
(This article belongs to the Special Issue Geographic Complexity: Concepts, Theories, and Practices)
Show Figures

Figure 1

Figure 1
<p>(Color online) The simplification process of the Koch snowflake guided by head/tail breaks. (Note: The blue polygons in each panel denote the head parts, whereas the red triangles represent the tail part, which needed to be eliminated progressively for generalization purpose).</p>
Full article ">Figure 2
<p>(Color online) Illustration of geometric measures for the (<b>a</b>) point-based simplifying unit (e.g., Douglas–Peucker (DP) algorithm) and (<b>b</b>) areal-based unit (e.g., bend simplify (BS) and hierarchical decomposition of a polygon (HD) algorithm).</p>
Full article ">Figure 3
<p>(Color online) Illustration of the HD algorithm using the Koch snowflake as a working example. Note that the size of derived convex hulls holds a striking scaling hierarchy of far more small ones than large ones.</p>
Full article ">Figure 4
<p>The tool for easily conducting cartographical simplification of polygonal features based on its inherent scaling hierarchies.</p>
Full article ">Figure 5
<p>Flow chart for conducting polygon simplification based on head/tail breaks for each algorithm.</p>
Full article ">Figure 6
<p>(Color online) Basic geometric units of a part of the British coastline (referring to the red box in the right panel) for each polygon simplification algorithm derived by PolySimp. (Note: Panels on the left show a clear scaling hierarchy of far more small ones than large ones, represented by either dot size or patch color).</p>
Full article ">Figure 7
<p>(Color online) The simplified results of the British coastline at different levels of details through four algorithms in terms of the polygonal shape based on x (<b>a</b>), d (<b>b</b>), and area (<b>c</b>); and the corresponding number of points (<b>d</b>) respectively.</p>
Full article ">Figure 8
<p>(Color Online) The shape variation of simplified results of the British coastline at five levels of details, indicated by the polygon’s area (Panels <b>a</b>–<b>c</b>), perimeter (Panels <b>e</b>–<b>g</b>), and shape factor (Panels <b>h</b>–<b>j</b>).</p>
Full article ">Figure 9
<p>(Color online) The first iteration of the simplification process of the Koch snowflake guided by area (note: Red triangles are the tail part that needs to be eliminated at the first level).</p>
Full article ">Figure 10
<p>(Color online) A multiscale representation with the scaling ratio of 1/2 of the simplified result from level 1 to 5 (Panels <b>a</b>–<b>e</b>) of the British coastline using the HD algorithm.</p>
Full article ">Figure A1
<p>(Color online) The illustration of simplifying geometric units among DP, BS, and VW algorithms using Koch snowflake and the generalization result at the first level. (Note: The size of each type of simplifying units, such as points, triangles, and bends, holds the underlying scaling hierarchy that includes a great many smalls (in blue), a few larges (in red), and some in between (in green), all of which can be used for the guidance of cartographical generalization).</p>
Full article ">
25 pages, 63671 KiB  
Article
Spatio-Temporal Relationship between Land Cover and Land Surface Temperature in Urban Areas: A Case Study in Geneva and Paris
by Xu Ge, Dasaraden Mauree, Roberto Castello and Jean-Louis Scartezzini
ISPRS Int. J. Geo-Inf. 2020, 9(10), 593; https://doi.org/10.3390/ijgi9100593 - 10 Oct 2020
Cited by 15 | Viewed by 4427
Abstract
Currently, more than half of the world’s population lives in cities, which leads to major changes in land use and land surface temperature (LST). The associated urban heat island (UHI) effects have multiple impacts on energy consumption and human health. A better understanding [...] Read more.
Currently, more than half of the world’s population lives in cities, which leads to major changes in land use and land surface temperature (LST). The associated urban heat island (UHI) effects have multiple impacts on energy consumption and human health. A better understanding of how different land covers affect LST is necessary for mitigating adverse impacts, and supporting urban planning and public health management. This study explores a distance-based, a grid-based and a point-based analysis to investigate the influence of impervious surfaces, green area and waterbodies on LST, from large (distance and grid based analysis with 400 m grids) to smaller (point based analysis with 30 m grids) scale in the two mid-latitude cities of Paris and Geneva. The results at large scale confirm that the highest LST was observed in the city centers. A significantly positive correlation was observed between LST and impervious surface density. An anticorrelation between LST and green area density was observed in Paris. The spatial lag model was used to explore the spatial correlation among LST, NDBI, NDVI and MNDWI on a smaller scale. Inverse correlations between LST and NDVI and MNDWI, respectively, were observed. We conclude that waterbodies display the greatest mitigation on LST and UHI effects both on the large and smaller scale. Green areas play an important role in cooling effects on the smaller scale. An increase of evenly distributed green area and waterbodies in urban areas is suggested to lower LST and mitigate UHI effects. Full article
(This article belongs to the Special Issue Geodata Science and Spatial Analysis in Urban Studies)
Show Figures

Figure 1

Figure 1
<p>Study area.</p>
Full article ">Figure 2
<p>Land cover data.</p>
Full article ">Figure 3
<p>Four transects in Geneva and Paris.</p>
Full article ">Figure 4
<p>400 m × 400 m grids.</p>
Full article ">Figure 5
<p>Mean Land Surface Temperature (°C) of four seasons in Paris.</p>
Full article ">Figure 6
<p>Mean Land Surface Temperature (°C) of four seasons in Geneva.</p>
Full article ">Figure 6 Cont.
<p>Mean Land Surface Temperature (°C) of four seasons in Geneva.</p>
Full article ">Figure 7
<p>Land surface temperature (°C) of three land cover features for each season in Paris and Geneva.</p>
Full article ">Figure 8
<p>Box charts related to land surface temperature of green area and impervious surfaces for four seasons observed in Paris and Geneva.</p>
Full article ">Figure 9
<p>Seasonal land surface temperature (°C) along the north to south transect in Geneva.</p>
Full article ">Figure 9 Cont.
<p>Seasonal land surface temperature (°C) along the north to south transect in Geneva.</p>
Full article ">Figure 10
<p>Seasonal land surface temperature (°C) along the north to south transect in Paris.</p>
Full article ">Figure 11
<p>Density map of impervious surfaces for four seasons in Paris. Different colored areas and sized circles indicate various temperatures (°C) and different densities respectively. Figures (<b>f</b>–<b>i</b>) present the relationships between LST and density of green area with 95% confidence intervals.</p>
Full article ">Figure 12
<p>Density map of impervious surfaces for four seasons in Geneva. Different colored areas and sized circles indicate various temperatures (°C) and different densities respectively. Figures (<b>f</b>–<b>i</b>) present the relationships between LST and density of green area with 95% confidence intervals.</p>
Full article ">Figure 13
<p>Density map of green area for four seasons in Paris. Different colored areas and sized circles indicate various temperatures (°C) and different density groups respectively. Figures (<b>f</b>–<b>i</b>) present the relationships between LST and density of green area with 95% confidence intervals.</p>
Full article ">Figure 14
<p>Density map of green area for four seasons in Geneva. Different colored areas and sized circles indicate various temperatures (°C) and different density groups respectively. Figures (<b>f</b>–<b>i</b>) present the relationships between LST and density of green area with 95% confidence intervals.</p>
Full article ">Figure 15
<p>Density of green area between 40 and 60 ha/km<sup>2</sup> in winter in Geneva.</p>
Full article ">
20 pages, 5248 KiB  
Technical Note
Combining UAV Imagery, Volunteered Geographic Information, and Field Survey Data to Improve Characterization of Rural Water Points in Malawi
by Marc van den Homberg, Arjen Crince, Jurg Wilbrink, Daniël Kersbergen, Gumbi Gumbi, Simon Tembo and Rob Lemmens
ISPRS Int. J. Geo-Inf. 2020, 9(10), 592; https://doi.org/10.3390/ijgi9100592 - 9 Oct 2020
Cited by 4 | Viewed by 3803
Abstract
As the world is digitizing fast, the increase in Big and Small Data offers opportunities to enrich official statistics for reporting on Sustainable Development Goals (SDG). However, survey data coming from an increased number of organizations (Small Data) and Big Data offer challenges [...] Read more.
As the world is digitizing fast, the increase in Big and Small Data offers opportunities to enrich official statistics for reporting on Sustainable Development Goals (SDG). However, survey data coming from an increased number of organizations (Small Data) and Big Data offer challenges in terms of data heterogeneity. This paper describes a methodology for combining various data sources to create a more comprehensive dataset on SDG 6.1.1. (proportion of population using safely managed drinking water services). We enabled digital volunteers to trace buildings on satellite imagery and used the traces on OpenStreetMap to facilitate visual detection of water points on Unmanned Aerial Vehicle (UAV) imagery and estimate the number of people served per water point. Combining data on water points identified on our UAV imagery with data on water points from field surveys improves the overall quality in terms of removal of inconsistencies and enrichment of attribute information. Satellite imagery enables scaling more easily than UAV imagery but is too costly to acquire at sufficiently high resolution. For small areas, our workflow is cost-effective in creating an up-to-date and consistent water point dataset by combining UAV imagery, Volunteered Geographic Information, and field survey data. Full article
Show Figures

Figure 1

Figure 1
<p>The research framework. Ground Sample Distance (GSD), Unmanned Area Vehicles (UAV). The dark green circles represent organizations that collected data via field surveys on water points: Climate Justice Fund (CJF), Water Point Data Exchange (WPDx), Department of Irrigation and Water Development (DoIWD), Department of Surveys (Dept Surveys).</p>
Full article ">Figure 2
<p>Overview of the case study area, part of the Traditional Authority Makhwira in the district Chikwawa.</p>
Full article ">Figure 3
<p>Images of the UAV mission.</p>
Full article ">Figure 4
<p>Overview of the different imagery for 175 × 125 m<sup>2</sup> for two different locations (top row and bottom row) in the case study area. The red cross refers to water points as identified on the UAV imagery. Green dots represent water points surveyed by Climate Justice Fund, red dots Water Point Data Exchange, yellow dots Department of Surveys, brown dots Madzi Alipo. In some cases, the brown and yellow dots coincide. No water points from the Department of Irrigation and Water Development survey were present in these segments.</p>
Full article ">Figure 5
<p>Protected water point.</p>
Full article ">Figure 6
<p>Non-functional water point.</p>
Full article ">Figure 7
<p>Functional water point.</p>
Full article ">Figure 8
<p>Water point closely located to latrine.</p>
Full article ">Figure 9
<p>Overview of water points (red circle) visually detected on the UAV imagery (red perimeter). Additionally, the water points from Madzi Alipo, Department of Irrigation and Water Development, Climate Justice Fund, Water Point Data Exchange, and Department of Surveys are depicted.</p>
Full article ">Figure 10
<p>Overview of the buffers (50 m, 100 m, 200 m, 500 m) created around the center point of the selected water points.</p>
Full article ">Figure 11
<p>Diagram showing how the number of OpenStreetMap buildings in different buffer zones around a water point can be calculated.</p>
Full article ">
10 pages, 851 KiB  
Review
Usability of IoT and Open Data Repositories for Analyzing Water Pollution. A Case Study in the Czech Republic
by Jan Pavlík, Markéta Hrnčírová, Michal Stočes, Jan Masner and Jiří Vaněk
ISPRS Int. J. Geo-Inf. 2020, 9(10), 591; https://doi.org/10.3390/ijgi9100591 - 8 Oct 2020
Cited by 2 | Viewed by 2574
Abstract
Recently, the process of data opening has intensified, especially thanks to the involvement of many institutions that have not yet shared their data. Some entities provided data to the public long before the trend of open data was pushed to a wider level, [...] Read more.
Recently, the process of data opening has intensified, especially thanks to the involvement of many institutions that have not yet shared their data. Some entities provided data to the public long before the trend of open data was pushed to a wider level, but many institutions have only engaged in this process recently thanks to a systemic state-level effort to make data repositories available to the public. Therefore, there are many new potential sources of data available for research, including the area of water management. This article analyses the current state of available data in the Czech Republic—their content, structure, format, availability, costs and other indicators that affect the usability of these data for independent researchers in the area of water management. The case study was conducted to ascertain the levels of accessibility and usability of data in open data repositories and the possibilities of obtaining data from IoT (Internet of Things) devices such as networked sensors where required data is either not available from existing sources, too costly, or otherwise unsuitable for the research. The goal of the underlying research was to assess the impact/ratio of various watershed factors based on monitored indicators of water pollution in a model watershed. Such information would help propose measures for reducing the volume of pollution resulting in increased security in terms of available drinking water for the capital city Prague. Full article
(This article belongs to the Special Issue Integrating GIS and Internet of Things (IoT) in Sustainable Cities)
Show Figures

Figure 1

Figure 1
<p>Boundaries of BP2 and BP5 profiles made in BNHelp geographical information system (GIS) application.</p>
Full article ">
34 pages, 9315 KiB  
Article
A Comparative Study of Several Metaheuristic Algorithms to Optimize Monetary Incentive in Ridesharing Systems
by Fu-Shiung Hsieh
ISPRS Int. J. Geo-Inf. 2020, 9(10), 590; https://doi.org/10.3390/ijgi9100590 - 8 Oct 2020
Cited by 27 | Viewed by 3154
Abstract
The strong demand on human mobility leads to excessive numbers of cars and raises the problems of serious traffic congestion, large amounts of greenhouse gas emissions, air pollution and insufficient parking space in cities. Although ridesharing is a potential transport mode to solve [...] Read more.
The strong demand on human mobility leads to excessive numbers of cars and raises the problems of serious traffic congestion, large amounts of greenhouse gas emissions, air pollution and insufficient parking space in cities. Although ridesharing is a potential transport mode to solve the above problems through car-sharing, it is still not widely adopted. Most studies consider non-monetary incentive performance indices such as travel distance and successful matches in ridesharing systems. These performance indices fail to provide a strong incentive for ridesharing. The goal of this paper is to address this issue by proposing a monetary incentive performance indicator to improve the incentives for ridesharing. The objectives are to improve the incentive for ridesharing through a monetary incentive optimization problem formulation, development of a solution methodology and comparison of different solution algorithms. A non-linear integer programming optimization problem is formulated to optimize monetary incentive in ridesharing systems. Several discrete metaheuristic algorithms are developed to cope with computational complexity for solving the above problem. These include several discrete variants of particle swarm optimization algorithms, differential evolution algorithms and the firefly algorithm. The effectiveness of applying the above algorithms to solve the monetary incentive optimization problem is compared based on experimental results. Full article
(This article belongs to the Special Issue GIS in Sustainable Transportation)
Show Figures

Figure 1

Figure 1
<p>A flowchart for the discrete particle swarm optimization (PSO) algorithm.</p>
Full article ">Figure 2
<p>A flowchart for the discrete CLPSO algorithm.</p>
Full article ">Figure 3
<p>A flowchart for the discrete CCPSO algorithm.</p>
Full article ">Figure 4
<p>A flowchart for the discrete firefly algorithm.</p>
Full article ">Figure 5
<p>A flowchart for the discrete DE algorithm.</p>
Full article ">Figure 6
<p>The results (obtained with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10 for Test Case 1) displayed on Google Maps.</p>
Full article ">Figure 7
<p>The bar chart for the average fitness function values of discrete PSO, CCPSO, CLPSO, ALPSO and FA algorithms (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10) created based on <a href="#ijgi-09-00590-t012" class="html-table">Table 12</a>.</p>
Full article ">Figure 8
<p>The bar chart for the average number of generations of discrete PSO, CCPSO, CLPSO, ALPSO and FA algorithms (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10) created based on <a href="#ijgi-09-00590-t012" class="html-table">Table 12</a>.</p>
Full article ">Figure 9
<p>The bar chart for the fitness function values of discrete DE algorithms with strategy 1, strategy 2, strategy 3, strategy 4, strategy 5 and strategy 6 (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10) created based on <a href="#ijgi-09-00590-t013" class="html-table">Table 13</a> and <a href="#ijgi-09-00590-t014" class="html-table">Table 14</a>.</p>
Full article ">Figure 10
<p>The bar chart for the average number of generations of discrete DE algorithms with strategy 1, strategy 2, strategy 3, strategy 4, strategy 5 and strategy 6 (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10) created based on <a href="#ijgi-09-00590-t013" class="html-table">Table 13</a> and <a href="#ijgi-09-00590-t014" class="html-table">Table 14</a>.</p>
Full article ">Figure 11
<p>The bar chart for the average fitness function values of discrete PSO, CCPSO, CLPSO, ALPSO and FA algorithms (population size <span class="html-italic">NP</span> = 30) created based on <a href="#ijgi-09-00590-t015" class="html-table">Table 15</a>.</p>
Full article ">Figure 12
<p>The bar chart for the average number of generations of discrete PSO, CCPSO, CLPSO, ALPSO and FA algorithms (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30) created based on <a href="#ijgi-09-00590-t007" class="html-table">Table 7</a>.</p>
Full article ">Figure 13
<p>The bar chart for the fitness function values of discrete DE algorithms with strategy 1, strategy 2, strategy 3, strategy 4, strategy 5 and strategy 6 (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30) created based on <a href="#ijgi-09-00590-t016" class="html-table">Table 16</a> and <a href="#ijgi-09-00590-t017" class="html-table">Table 17</a>.</p>
Full article ">Figure 14
<p>The bar chart for the average number of generations of discrete DE algorithms with strategy 1, strategy 2, strategy 3, strategy 4, strategy 5 and strategy 6 (population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30) created based on <a href="#ijgi-09-00590-t016" class="html-table">Table 16</a> and <a href="#ijgi-09-00590-t017" class="html-table">Table 17</a>.</p>
Full article ">Figure 15
<p>Simulation runs of Test Case 2 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <span class="html-italic">NP</span> = 10.</p>
Full article ">Figure 16
<p>Simulation runs of Test Case 5 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10.</p>
Full article ">Figure 17
<p>Simulation runs of Test Case 7 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10.</p>
Full article ">Figure 18
<p>Simulation runs of Test Case 8 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 10.</p>
Full article ">Figure 19
<p>Simulation runs of Test Case 2 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30.</p>
Full article ">Figure 20
<p>Simulation runs of Test Case 5 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30.</p>
Full article ">Figure 21
<p>Simulation runs of Test Case 7 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30.</p>
Full article ">Figure 22
<p>Simulation runs of Test Case 8 obtained by discrete PSO, CCPSO, CLPSO, ALPSO, FA and DE algorithms with population size <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>P</mi> </mrow> </semantics></math> = 30.</p>
Full article ">
15 pages, 2567 KiB  
Article
Evaluation of the Space Syntax Measures Affecting Pedestrian Density through Ordinal Logistic Regression Analysis
by Özge Öztürk Hacar, Fatih Gülgen and Serdar Bilgi
ISPRS Int. J. Geo-Inf. 2020, 9(10), 589; https://doi.org/10.3390/ijgi9100589 - 7 Oct 2020
Cited by 8 | Viewed by 4154
Abstract
This paper examines the relationship between pedestrian density and space syntax measures in a university campus using ordinal logistic regression analysis. The pedestrian density assumed as the dependent variable of regression analysis was categorised in low, medium, and high classes by using Jenks [...] Read more.
This paper examines the relationship between pedestrian density and space syntax measures in a university campus using ordinal logistic regression analysis. The pedestrian density assumed as the dependent variable of regression analysis was categorised in low, medium, and high classes by using Jenks natural break classification. The data elements of groups were derived from pedestrian counts performed in 22 gates 132 times. The counting period grouped in nominal categories was assumed as an independent variable. Another independent was one of the 15 derived measures of axial analysis and visual graphic analysis. The statistically significant model results indicated that the integration of axial analysis was the most reasonable measure that explained the pedestrian density. Then, the changes in integration values of current and master plan datasets were analysed using paired sample t-test. The calculated p-value of t-test proved that the master plan would change the campus morphology for pedestrians. Full article
(This article belongs to the Special Issue Measuring, Mapping, Modeling, and Visualization of Cities)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Buildings and roads of Davutpasa Campus digitised from orthophotos.</p>
Full article ">Figure 2
<p>Measures obtained from the space syntax (SS) analyses.</p>
Full article ">Figure 3
<p>Box plot graphic for gate-count data.</p>
Full article ">Figure 4
<p>AA Integration values for (<b>a</b>) current and (<b>b</b>) master plan dataset.</p>
Full article ">
30 pages, 20417 KiB  
Article
A Built Heritage Information System Based on Point Cloud Data: HIS-PC
by Florent Poux, Roland Billen, Jean-Paul Kasprzyk, Pierre-Henri Lefebvre and Pierre Hallot
ISPRS Int. J. Geo-Inf. 2020, 9(10), 588; https://doi.org/10.3390/ijgi9100588 - 7 Oct 2020
Cited by 17 | Viewed by 4568
Abstract
The digital management of an archaeological site requires to store, organise, access and represent all the information that is collected on the field. Heritage building information modelling, archaeological or heritage information systems now tend to propose a common framework where all the materials [...] Read more.
The digital management of an archaeological site requires to store, organise, access and represent all the information that is collected on the field. Heritage building information modelling, archaeological or heritage information systems now tend to propose a common framework where all the materials are managed from a central database and visualised through a 3D representation. In this research, we offer the development of a built heritage information system prototype based on a high-resolution 3D point cloud data set. The particularity of the approach is to consider a user-centred development methodology while avoiding meshing/down-sampling operations. The proposed system is initiated by a close collaboration between multi-modal users (managers, visitors, curators) and a development team (designers, developers, architects). The developed heritage information system permits the management of spatial and temporal information, including a wide range of semantics using relational along with NoSQL databases. The semantics used to describe the artifacts are subject to conceptual modelling. Finally, the system proposes a bi-directional communication with a 3D interface able to stream massive point clouds, which is a big step forward to provide a comprehensive site representation for stakeholders while minimising modelling costs. Full article
(This article belongs to the Special Issue BIM for Cultural Heritage (HBIM))
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Computer architecture of the BIMLegacy project (redrawn from [<a href="#B45-ijgi-09-00588" class="html-bibr">45</a>]).</p>
Full article ">Figure 2
<p>UML class diagram representing the data model of the database (from [<a href="#B47-ijgi-09-00588" class="html-bibr">47</a>]).</p>
Full article ">Figure 3
<p>A representation of the palace which dominated the city of Brussels, the object of this research.</p>
Full article ">Figure 4
<p>The actual ruins of the palace mentioned above, standing as an Underground Museum.</p>
Full article ">Figure 5
<p>Methodological approach of user-centred and quality information system design.</p>
Full article ">Figure 6
<p>The four steps of the user-centred approach and their main deliverables: 1.A—user stories, 1.B—specifications, 2.A—conceptual data model, 2.B—IS interfaces and architecture mock-up, 2.C.—data acquisition methodology, 3.A—IS implementation, 3.B—data acquisition, 4.A—training sessions, 4.B—usability report. Black human icons refer to users, and grey human icons refer to developers.</p>
Full article ">Figure 7
<p>Metamodel’s legend.</p>
Full article ">Figure 8
<p>The simplified version of the metamodel adopted for the Coudenberg project.</p>
Full article ">Figure 9
<p>Data Model of the Coudenberg Palace IS.</p>
Full article ">Figure 10
<p>Architecture of the SIA (conceptual model).</p>
Full article ">Figure 11
<p>Use cases diagram of the HIS.</p>
Full article ">Figure 12
<p>Surveyed areas and polygonal of the Coudenberg Museum archaeological site.</p>
Full article ">Figure 13
<p>Point Cloud data after registration.</p>
Full article ">Figure 14
<p>3D point cloud colourised with the internal camera.</p>
Full article ">Figure 15
<p>Laser scanning stations horizontal sections.</p>
Full article ">Figure 16
<p>Geo-referenced point cloud of the Aula Magna.</p>
Full article ">Figure 17
<p>Registration of different photos with laser-scanning data. In white, the different positions of the pictures.</p>
Full article ">Figure 18
<p>Mesh from the Aula Magna—with (<b>right</b>) and without texture (<b>left</b>).</p>
Full article ">Figure 19
<p>Implemented architecture of the Coudenberg Palace IS.</p>
Full article ">Figure 20
<p>The semantic interface of the HIS.</p>
Full article ">Figure 21
<p>The usage of both the measurement and labelling tools on a section of the point cloud within the point cloud Interface.</p>
Full article ">Figure 22
<p>The Geoverse interface.</p>
Full article ">Figure 23
<p>The spatial API scheme connecting Geoverse with PostGIS.</p>
Full article ">Figure 24
<p>Extract of the user manual and video tutorial for the end users.</p>
Full article ">Figure 25
<p>On-site acquisition with an adapted system to reduce problems in low-light conditions.</p>
Full article ">
18 pages, 4627 KiB  
Article
A Machine Learning-Based Approach for Spatial Estimation Using the Spatial Features of Coordinate Information
by Seongin Ahn, Dong-Woo Ryu and Sangho Lee
ISPRS Int. J. Geo-Inf. 2020, 9(10), 587; https://doi.org/10.3390/ijgi9100587 - 6 Oct 2020
Cited by 10 | Viewed by 3558
Abstract
With the development of machine learning technology, research cases for spatial estimation through machine learning approach (MLA) in addition to the traditional geostatistical techniques are increasing. MLA has the advantage that spatial estimation is possible without stationary hypotheses of data, but it is [...] Read more.
With the development of machine learning technology, research cases for spatial estimation through machine learning approach (MLA) in addition to the traditional geostatistical techniques are increasing. MLA has the advantage that spatial estimation is possible without stationary hypotheses of data, but it is possible for the prediction results to ignore spatial autocorrelation. In recent studies, it was considered by using a distance matrix instead of raw coordinates. Although, the performance of spatial estimation could be improved through this approach, the computational complexity of MLA increased rapidly as the number of sample points increased. In this study, we developed a method to reduce the computational complexity of MLA while considering spatial autocorrelation. Principal component analysis is applied to it for extracting spatial features and reducing dimension of inputs. To verify the proposed approach, indicator Kriging was used as a benchmark model, and each performance of MLA was compared when using raw coordinates, distance vector, and spatial features extracted from distance vector as inputs. The proposed approach improved the performance compared to previous MLA and showed similar performance compared with Kriging. We confirmed that extracted features have characteristics of rigid classification in spatial estimation; on this basis, we conclude that the model could improve performance. Full article
Show Figures

Figure 1

Figure 1
<p>Comparison of spatial estimation approaches. Schematic difference between (<b>a</b>) Kriging and (<b>b</b>) the machine-learning approach (MLA) for spatial estimation.</p>
Full article ">Figure 2
<p>Information about Meuse dataset: (<b>a</b>) a spatial distribution of zinc concentrations and (<b>b</b>) a histogram with basic statistics.</p>
Full article ">Figure 3
<p>Information about the Seoul borehole dataset: (<b>a</b>) a spatial distribution of the deposit soil thickness and (<b>b</b>) a histogram with basic statistics.</p>
Full article ">Figure 4
<p>Results of dividing sample data into five folds considering unbiased spatial distribution for (<b>a</b>) Meuse and (<b>b</b>) Seoul borehole datasets.</p>
Full article ">Figure 5
<p>Experimental (dot) and theoretical (line) variograms calculated for thresholds (<b>a</b>) one to (<b>i</b>) nine of the Meuse dataset.</p>
Full article ">Figure 6
<p>Experimental (dot) and theoretical (line) variograms calculated for thresholds (<b>a</b>) one to (<b>i</b>) nine of the Seoul borehole dataset.</p>
Full article ">Figure 7
<p>Results of the performances of spatial estimation according to the number of principal components (PCs): (<b>a</b>) <span class="html-italic">R</span>-squared and (<b>b</b>) RMSE for the Meuse dataset; (<b>c</b>) <span class="html-italic">R</span>-squared and (<b>d</b>) RMSE for the Seoul dataset.</p>
Full article ">Figure 8
<p>Prediction performance results: (<b>a</b>) Meuse dataset and (<b>b</b>) Seoul dataset.</p>
Full article ">Figure 9
<p>Comparison of predictions based on each methodology for the Meuse dataset: predicted zinc concentrations ((<b>a</b>) random forest with coordinate input (RF-Coord); (<b>b</b>) RF with distance input (RF-Dist); (<b>c</b>) RF with principal component analysis (RF-PCA); (<b>d</b>) indicator kriging (IK)) and standard deviation of prediction error ((<b>e</b>) RF-Coord; (<b>f</b>) RF-Dist; (<b>g</b>) RF-PCA; (<b>h</b>) IK).</p>
Full article ">Figure 10
<p>Comparison of predictions based on each methodology for Seoul dataset: predicted deposit soil thickness ((<b>a</b>) RF-Coord; (<b>b</b>) RF-Dist; (<b>c</b>) RF-PCA; (<b>d</b>) IK) and standard deviation of prediction error ((<b>e</b>) RF-Coord; (<b>f</b>) RF-Dist; (<b>g</b>) RF-PCA; (<b>h</b>) IK).</p>
Full article ">Figure 11
<p>Spatial prediction from the RF-PCA using fifteen PCs for the Meuse dataset: (<b>a</b>) Including the third PC and (<b>b</b>) excluding the third PC.</p>
Full article ">Figure 12
<p>Results of five-fold cross-validation of RF-PCA according to the number of PCs increasing for the (<b>a</b>) Meuse dataset and (<b>b</b>) Seoul dataset.</p>
Full article ">
14 pages, 5131 KiB  
Article
Urban Green Accessibility Index: A Measure of Pedestrian-Centered Accessibility to Every Green Point in an Urban Area
by Kee Moon Jang, Jaeman Kim, Hye-Yeong Lee, Hyemin Cho and Youngchul Kim
ISPRS Int. J. Geo-Inf. 2020, 9(10), 586; https://doi.org/10.3390/ijgi9100586 - 6 Oct 2020
Cited by 17 | Viewed by 6217
Abstract
Advancements in remote sensing techniques and urban data analysis tools have enabled the successful monitoring and detection of green spaces in a city. This study aims to develop an index called the urban green accessibility (UGA) index, which measures people’s accessibility to green [...] Read more.
Advancements in remote sensing techniques and urban data analysis tools have enabled the successful monitoring and detection of green spaces in a city. This study aims to develop an index called the urban green accessibility (UGA) index, which measures people’s accessibility to green space and represents the citywide or local characteristics of the distribution pattern of green space. The index is defined as the sum of pedestrians’ accessibility to all vegetation points, which consists of the normalized difference vegetation index (NDVI) with integration and choice values from angular segment analysis. In this study, the proposed index is tested with cases of New York, NY, and San Francisco, CA, in the US. The results reveal differences based on the significance of streets. When analysis ranges are on a neighborhood scale, a few hotspots appear in well-known green areas on commonly accessible streets and in local neighborhood parks on residential blocks. The appearance of high-accessibility points in low-NDVI areas implies the potential of the efficient and proper distribution of green spaces for pedestrians. The proposed measure is expected to help in planning and managing green areas in cities, taking people’s accessibility and spatial relationships into consideration. Full article
(This article belongs to the Special Issue Geodata Science and Spatial Analysis in Urban Studies)
Show Figures

Figure 1

Figure 1
<p>Access situations to green points near a street in a city.</p>
Full article ">Figure 2
<p>Analyses of collected data from San Francisco. (<b>a</b>) NDVI analysis process (aerial image–band manipulation–vectorization–close-up view); (<b>b</b>) Angular segment analysis result (integration).</p>
Full article ">Figure 3
<p>Urban green accessibility distribution with choice index applied as the topological value. (<b>a</b>) New York County (<b>b</b>) San Francisco County.</p>
Full article ">Figure 4
<p>Urban green accessibility distribution with integration index applied as the topological value. (<b>a</b>) New York County (<b>b</b>) San Francisco County.</p>
Full article ">Figure 4 Cont.
<p>Urban green accessibility distribution with integration index applied as the topological value. (<b>a</b>) New York County (<b>b</b>) San Francisco County.</p>
Full article ">Figure 5
<p>Comparison of significant areas in New York. (<b>a</b>) NDVI analysis result; (<b>b</b>) UGA heat map with choice 1600 m setting; (<b>c</b>) UGA heat map with integration 1600 m setting.</p>
Full article ">Figure 6
<p>Comparison of significant areas in San Francisco. (<b>a</b>) NDVI analysis result; (<b>b</b>) UGA heat map with choice 1600 m setting; (<b>c</b>) UGA heat map with integration 1600 m setting.</p>
Full article ">Figure 7
<p>Locations on the maps of (<b>a</b>) New York and (<b>b</b>) San Francisco at the areas of A, B and C.</p>
Full article ">Figure 7 Cont.
<p>Locations on the maps of (<b>a</b>) New York and (<b>b</b>) San Francisco at the areas of A, B and C.</p>
Full article ">
17 pages, 6506 KiB  
Article
Measuring Spatial Accessibility of Urban Fire Services Using Historical Fire Incidents in Nanjing, China
by Kainan Mao, Yuehong Chen, Guohao Wu, Junwang Huang, Wanying Yang and Zelong Xia
ISPRS Int. J. Geo-Inf. 2020, 9(10), 585; https://doi.org/10.3390/ijgi9100585 - 6 Oct 2020
Cited by 17 | Viewed by 3449
Abstract
The measurement of spatial accessibility of fire services is a key task in enhancing fire response efficiency and minimizing property losses and deaths. Recently, the two-step floating catchment area method and its modified versions have been widely applied. However, the circle catchment areas [...] Read more.
The measurement of spatial accessibility of fire services is a key task in enhancing fire response efficiency and minimizing property losses and deaths. Recently, the two-step floating catchment area method and its modified versions have been widely applied. However, the circle catchment areas used in these methods are not suitable for measuring the accessibility of fire services because each fire station is often responsible for the fire incidents within its coverage. Meanwhile, most existing methods take the demographic data and their centroids of residential areas as the demands and locations, respectively, which makes it difficult to reflect the actual demands and locations of fire services. Thus, this paper proposes a fixed-coverage-based two-step floating catchment area (FC2SFCA) method that takes the fixed service coverage of fire stations as the catchment area and the locations and dispatched fire engines of historical fire incidents as the demand location and size, respectively, to measure the spatial accessibility of fire services. Using a case study area in Nanjing, China, the proposed FC2SFCA and enhanced two-step floating catchment area (E2SFCA) are employed to measure and compare the spatial accessibility of fire incidents and fire stations. The results show that (1) the spatial accessibility across Nanjing, China is unbalanced, with relatively high spatial accessibility in the areas around fire stations and the southwest and northeast at the city center area and relatively low spatial accessibility in the periphery and boundary of the service coverage areas and the core of the city center; (2) compared with E2SFCA, FC2SFCA is less influenced by other fire stations and provides greater actual fire service accessibility; (3) the spatial accessibility of fire services is more strongly affected by the number of fire incidents than firefighting capabilities, the area of service coverage, or the average number of crossroads (per kilometer). Suggestions are then made to improve the overall spatial access to fire services. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The spatial distribution of (<b>a</b>) the service coverage of fire stations and (<b>b</b>) historical fire incidents in Nanjing, China.</p>
Full article ">Figure 2
<p>Histogram of the firefighters and fire engines in each fire station.</p>
Full article ">Figure 3
<p>Flowchart of the proposed fixed-coverage-based two-step floating catchment area (FC2SFCA) method.</p>
Full article ">Figure 4
<p>Comparison of E2SFCA and FC2SFCA for generating the catchment area.</p>
Full article ">Figure 5
<p>The probability density distribution histogram of the spatial accessibility score of fire incidents.</p>
Full article ">Figure 6
<p>The spatial accessibility of fire incidents using (<b>a</b>) E2SFCA and (<b>b</b>) FC2SFCA. (<b>c</b>) city center of (<b>a</b>), (<b>d</b>) city center of (<b>b</b>).</p>
Full article ">Figure 7
<p>Optimized Hot Spot Analysis of spatial accessibility of fire incidents using (<b>a</b>) E2SFCA, (<b>b</b>) FC2SFCA.</p>
Full article ">Figure 8
<p>The spatial accessibility of fire stations.</p>
Full article ">Figure 9
<p>The spatial accessibility levels of fire stations using (<b>a</b>) E2SFCA, and (<b>b</b>) FC2SFCA.</p>
Full article ">Figure 10
<p>The correlation between the spatial accessibility score (logarithm) and the firefighting capability, the fire incidents’ count, the area, the average crossroads.</p>
Full article ">
25 pages, 11279 KiB  
Article
ADAtools: Automatic Detection and Classification of Active Deformation Areas from PSI Displacement Maps
by J. A. Navarro, R. Tomás, A. Barra, J. I. Pagán, C. Reyes-Carmona, L. Solari, J. L. Vinielles, S. Falco and M. Crosetto
ISPRS Int. J. Geo-Inf. 2020, 9(10), 584; https://doi.org/10.3390/ijgi9100584 - 6 Oct 2020
Cited by 30 | Viewed by 4123
Abstract
This work describes the set of tools developed, tested, and put into production in the context of the H2020 project Multi-scale Observation and Monitoring of Railway Infrastructure Threats (MOMIT). This project, which ended in 2019, aimed to show how the use of various [...] Read more.
This work describes the set of tools developed, tested, and put into production in the context of the H2020 project Multi-scale Observation and Monitoring of Railway Infrastructure Threats (MOMIT). This project, which ended in 2019, aimed to show how the use of various remote sensing techniques could help to improve the monitoring of railway infrastructures, such as tracks or bridges, and thus, consequently, improve the detection of ground instabilities and facilitate their management. Several lines of work were opened by MOMIT, but the authors of this work concentrated their efforts in the design of tools to help the detection and identification of ground movements using synthetic aperture radar interferometry (InSAR) data. The main output of this activity was a set of tools able to detect the areas labelled active deformation areas (ADA), with the highest deformation rates and to connect them to a geological or anthropogenic process. ADAtools is the name given to the aforementioned set of tools. The description of these tools includes the definition of their targets, inputs, and outputs, as well as details on how the correctness of the applications was checked and on the benchmarks showing their performance. The ADAtools include the following applications: ADAfinder, los2hv, ADAclassifier, and THEXfinder. The toolset is targeted at the analysis and interpretation of InSAR results. Ancillary information supports the semi-automatic interpretation and classification process. Two real use-cases illustrating this statement are included at the end of this paper to show the kind of results that may be obtained with the ADAtools. Full article
Show Figures

Figure 1

Figure 1
<p>The graphical user interface (GUI) version of the ADAfinder tool.</p>
Full article ">Figure 2
<p>An example of the output of ADAfinder in the coastal area of Granada (S Spain). Colors are used to show the quality of the time series (TS) information according to the quality index attribute in the output shapefile. Red: “very reliable”, orange: “reliable”, lime-green: “not so reliable” and purple: “not reliable”.</p>
Full article ">Figure 3
<p>Sketch of the algorithms of ADAclassifier, modified from [<a href="#B4-ijgi-09-00584" class="html-bibr">4</a>]. Note that the Th1–Th11 labels in the diagram correspond to some thresholds described in detail in the ADAclassifier user guide [<a href="#B34-ijgi-09-00584" class="html-bibr">34</a>].</p>
Full article ">Figure 4
<p>The ADAclassifier GUI (options tab).</p>
Full article ">Figure 5
<p>The ADAclassifier GUI (files tab).</p>
Full article ">Figure 6
<p>The THEXfinder GUI (options tab).</p>
Full article ">Figure 7
<p>The THEXfinder GUI (files tab).</p>
Full article ">Figure 8
<p>los2hv GUI.</p>
Full article ">Figure 9
<p>los2v: tesserae and ascending and descending persistent scatterers (PS). The white points represent the PS in ascending orbit; the red points the PS in descending orbit. Source: [<a href="#B25-ijgi-09-00584" class="html-bibr">25</a>].</p>
Full article ">Figure 10
<p>A real example of a los2hv options file.</p>
Full article ">Figure 11
<p>Example of an input shapefile read-map plain text file.</p>
Full article ">Figure 12
<p>ADAclassifier: synthetic dataset samples. (<b>a</b>) represents the checkerboard pattern for ADA, while (<b>b</b>) shows the stripes with values for the horizontal components of the movement. The green rectangle in (<b>c</b>) depicts the area where positive results for the landslides algorithm should be expected since there the conditions set by the algorithms are satisfied. Source: [<a href="#B25-ijgi-09-00584" class="html-bibr">25</a>].</p>
Full article ">Figure 13
<p>Synthetic dataset for los2hv. The green tiles are the only ones where both ascending and descending PS exist, and therefore, the unique areas where the horizontal component of the movement may be computed.</p>
Full article ">Figure 14
<p>(<b>a</b>) Location of the real test cases in the province of Granada (Spain) and Vibo Valentia (Calabria, Italy); (<b>b</b>) Detail of the test area located in the Granada coast. (<b>c</b>) Test site of Vibo Valentia; and (<b>d</b>) Detail of the test area of Vibo Valentia.</p>
Full article ">Figure 15
<p>Classification of ADA in the study area located in the South of Spain (province of Granada, Andalucía).</p>
Full article ">Figure 16
<p>Classification of active deformation areas (ADA) in the study area of Tropea-Zaccanopoli as: (<b>a</b>) landslides; (<b>b</b>) subsidence; (<b>c</b>) consolidation settlements; and (<b>d</b>) sinkholes.</p>
Full article ">
16 pages, 3746 KiB  
Article
Multitemporal Analysis of Deforestation in Response to the Construction of the Tucuruí Dam
by Andres Velastegui-Montoya, Aline de Lima and Marcos Adami
ISPRS Int. J. Geo-Inf. 2020, 9(10), 583; https://doi.org/10.3390/ijgi9100583 - 3 Oct 2020
Cited by 26 | Viewed by 5735
Abstract
The expansion of hydroelectric dams that is planned, and under construction, in the Amazon basin is a proposal to generate “clean” energy, with the purposes of meeting the regional energy demand, and the insertion of Brazil into the international economic market. However, this [...] Read more.
The expansion of hydroelectric dams that is planned, and under construction, in the Amazon basin is a proposal to generate “clean” energy, with the purposes of meeting the regional energy demand, and the insertion of Brazil into the international economic market. However, this type of megaproject can change the dynamics of natural ecosystems. In the present article, the spatiotemporal patterns of deforestation according to distance from the reservoir in the vicinity of the lake of Tucuruí, and within a radius of 30 km from it, are analyzed. A linear spectral mixture model of segmented Landsat-thematic mapper (TM), enhanced thematic mapper plus (ETM+), and operational land imager (OLI) images, and proximity analysis were used for the mapping of the land-cover classes in the vicinity of the artificial lake of Tucuruí. Likewise, landscape metrics were determined with the purpose of quantifying the reduction of primary forest, as a mechanism of loss of ecosystem services in the region. These methods were also used for the evaluation of the influence of the distance from the reservoir on the expansion of anthropogenic activities. This methodology was used for the scenarios of pre-inauguration, completion of phase I, beginning of construction phase II, full completion of the Tucuruí hydroelectric project, and the current scenario of the region. The results showed that the highest deforestation rate occurred in the first period of the analysis, due to the areas submerged by the reservoir and due to the anthropogenic disturbances, such as timber extraction, road construction, and the conversion of forests into large areas of agribusiness. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area: the lake of Tucuruí and the radius of 30 km around the reservoir.</p>
Full article ">Figure 2
<p>Thematic maps of land-cover within the 30 km radius of Lake Tucuruí in the five scenarios from 1984 to 2017.</p>
Full article ">Figure 3
<p>The percentage of forest per year in the six buffer rings surrounding the flooded area. The data belong to the scenarios of 1984, 1988, 1999, 2010, and 2017, for the entire study area, and part of the municipalities of Tucuruí, Breu Branco, Goianésia do Pará, Jacundá, Nova Ipixuna, Novo Repartimento, and Itupiranga. B<sub>0–5</sub>—Buffer from 0 km to 5 km; B<sub>5–10</sub>—Buffer from 5 km to 10 km; B<sub>10–15</sub>—Buffer from 10 km to 15 km; B<sub>15–20</sub>—Buffer from 15 km to 20 km; B<sub>20–25</sub>—Buffer from 20 km to 25 km; and B<sub>25–30</sub>—Buffer from 25 km to 30 km.</p>
Full article ">Figure 4
<p>Landscape change as measured by landscape-level metrics from 1984 to 2017.</p>
Full article ">Figure 5
<p>Habitats loss and fragmentation as measured by class-level metrics from 1984 to 2017.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop