[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Issue
Volume 7, December
Previous Issue
Volume 7, October
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 

ISPRS Int. J. Geo-Inf., Volume 7, Issue 11 (November 2018) – 39 articles

Cover Story (view full-size image): The increase in deprived living conditions in many cities of the Global South contradicts efforts to make them inclusive, safe, resilient, and sustainable places. Global statistics on the “proportion of urban population living in slums, informal settlements, or inadequate housing” (SDG indicator 11.1.1) are highly uncertain due to measurement and reporting problems. In addition, the data cannot be localized to effectively support upgrading. Using examples from across the globe, we show the scope of earth observation-based mapping of deprived living conditions and related population estimations in support of providing consistent information for the SDG indicator 11.1.1, and address uncertainties and geo-ethics. 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:
26 pages, 6900 KiB  
Article
Place versus Space: From Points, Lines and Polygons in GIS to Place-Based Representations Reflecting Language and Culture
by Thomas Blaschke, Helena Merschdorf, Pablo Cabrera-Barona, Song Gao, Emmanuel Papadakis and Anna Kovacs-Györi
ISPRS Int. J. Geo-Inf. 2018, 7(11), 452; https://doi.org/10.3390/ijgi7110452 - 19 Nov 2018
Cited by 34 | Viewed by 13023
Abstract
Around the globe, Geographic Information Systems (GISs) are well established in the daily workflow of authorities, businesses and non-profit organisations. GIS can effectively handle spatial entities and offer sophisticated analysis and modelling functions to deal with space. Only a small fraction of the [...] Read more.
Around the globe, Geographic Information Systems (GISs) are well established in the daily workflow of authorities, businesses and non-profit organisations. GIS can effectively handle spatial entities and offer sophisticated analysis and modelling functions to deal with space. Only a small fraction of the literature in Geographic Information Science—or GIScience in short—has advanced the development of place, addressing entities with an ambiguous boundary and relying more on the human or social attributes of a location rather than on crisp geographic boundaries. While the GIScience developments support the establishment of the digital humanities, GISs were never designed to handle subjective or vague data. We, an international group of authors, juxtapose place and space in English language and in several other languages and discuss potential consequences for Geoinformatics and GIScience. In particular, we address the question of whether linguistic and cultural settings play a role in the perception of place. We report on some facts revealed by this multi-language and multi-cultural dialogue, and what particular aspects of place we were able to discern regarding the few languages addressed. Full article
(This article belongs to the Special Issue Place-Based Research in GIScience and Geoinformatics)
Show Figures

Figure 1

Figure 1
<p>The perception of place.</p>
Full article ">Figure 2
<p>“Place, which place?”.</p>
Full article ">Figure 3
<p>Crowdsourced data from Twitter and Flickr to define the concept of “Downtown Santa Barbara.” Legend: <span class="html-italic">sb</span> stands for Santa Barbara, SDE for <span class="html-italic">Standard Deviational Ellipse</span>, an analysis method to identify the significant points, which is robust to outliers and summarises the central tendency and directional trend of point distributions.</p>
Full article ">Figure 4
<p>Placemaking: A concept that calls for a place-based GIS representation.</p>
Full article ">
14 pages, 14802 KiB  
Article
Direct Impacts of Climate Change and Indirect Impacts of Non-Climate Change on Land Surface Phenology Variation across Northern China
by Zhaohui Luo, Qingmei Song, Tao Wang, Huanmu Zeng, Tao He, Hengjun Zhang and Wenchen Wu
ISPRS Int. J. Geo-Inf. 2018, 7(11), 451; https://doi.org/10.3390/ijgi7110451 - 19 Nov 2018
Cited by 8 | Viewed by 3881
Abstract
Land surface phenology (LSP) is a sensitive indicator of climate change. Understanding the variation in LSP under various impacts can improve our knowledge on ecosystem dynamics and biosphere-atmosphere interactions. Over recent decades, LSP derived from remote sensing data and climate change-related variation of [...] Read more.
Land surface phenology (LSP) is a sensitive indicator of climate change. Understanding the variation in LSP under various impacts can improve our knowledge on ecosystem dynamics and biosphere-atmosphere interactions. Over recent decades, LSP derived from remote sensing data and climate change-related variation of LSP have been widely reported at the regional and global scales. However, the smoothing methods of the vegetation index (i.e., NDVI) are diverse, and discrepancies among methods may result in different results. Additionally, LSP is affected by climate change and non-climate change simultaneously. However, few studies have focused on the isolated impacts of climate change and the impacts of non-climate change on LSP variation. In this study, four methods were applied to reconstruct the MODIS enhanced vegetation index (EVI) dataset to choose the best smoothing result to estimate LSP. Subsequently, the variation in the start of season (SOS) and end of season (EOS) under isolated impacts of climate change were analyzed. Furthermore, the indirect effects of isolated impacts of non-climate change were conducted based on the differences between the combined impact (the impacts of both climate change and non-climate change) and isolated impacts of climate change. Our results indicated that the Savitzky-Golay method is the best method of the four for smoothing EVI in Northern China. Additionally, SOS displayed an advanced trend under the impacts of both climate change and non-climate change (hereafter called the combined impact), isolated impacts of climate change, and isolated impacts of non-climate change, with mean values of −0.26, −0.07, and −0.17 days per year, respectively. Moreover, the trend of SOS continued after 2000, but the magnitudes of changes in SOS after 2000 were lower than those that were estimated over the last two decades of the twentieth century (previous studies). EOS showed a delayed trend under the combined impact and isolated impacts of non-climate change, with mean values of 0.41 and 0.43 days per year, respectively. However, EOS advanced with a mean value of −0.16 days per year under the isolated impacts of climate change. Furthermore, the absolute mean values of SOS and EOS trends under the isolated impacts of non-climate change were larger than that of the isolated impacts of climate change, indicating that the effect of non-climate change on LSP variation was larger than that of climate change. With regard to the relative contribution of climatic factors to the variation in SOS and EOS, the proportion of solar radiation was the largest for both SOS and EOS, followed by precipitation and temperature. Full article
Show Figures

Figure 1

Figure 1
<p>Trends of the start of season (SOS) (<b>a</b>) and end of season (EOS) (<b>c</b>) under the impacts of climate change and the impacts of non-climate change. (<b>b</b>,<b>d</b>) show the corresponding frequency distributions of SOS and EOS, respectively.</p>
Full article ">Figure 2
<p>Isolated impacts of climatic factors (<b>a</b>) and non-climate change (<b>c</b>) on the variation in the SOS. (<b>b</b>,<b>d</b>) are the corresponding frequency distributions for climatic factors (<b>a</b>) and non-climate change (<b>c</b>), respectively.</p>
Full article ">Figure 3
<p>Isolated impacts of climatic factors (<b>a</b>) and non-climate change (<b>c</b>) on the variation in EOS. (<b>b</b>,<b>d</b>) are the corresponding frequency distributions for climatic factors (<b>a</b>) and non-climate change (<b>c</b>), respectively.</p>
Full article ">Figure 4
<p>Relative contributions of climatic factors to the variation in the SOS (<b>a</b>) and EOS (<b>b</b>), respectively. Tem, temperature; Pre, precipitation; Rad, solar radiation.</p>
Full article ">
23 pages, 10554 KiB  
Article
Real-Time Efficient Exploration in Unknown Dynamic Environments Using MAVs
by Haytham Mohamed, Adel Moussa, Mohamed Elhabiby and Naser El-Sheimy
ISPRS Int. J. Geo-Inf. 2018, 7(11), 450; https://doi.org/10.3390/ijgi7110450 - 18 Nov 2018
Viewed by 3045
Abstract
Micro aerial vehicles (MAVs) have been acknowledged as an influential technology for indoor search and rescue operations. The time constraint is a crucial factor in most search and rescue operations. The employed MAVs in indoor environments are characterized by short endurance flight time [...] Read more.
Micro aerial vehicles (MAVs) have been acknowledged as an influential technology for indoor search and rescue operations. The time constraint is a crucial factor in most search and rescue operations. The employed MAVs in indoor environments are characterized by short endurance flight time and limited payload weights. Hence, adding more batteries to extend the flight time is practically not feasible. Typically, most of the indoor missions’ environments might not be accessed and remain unknown. Working in such environments requires effective exploration and information gathering to save time and maximize the coverage area. Furthermore, due to the dynamism of such environments, choosing the least risky trajectory is an important task. This paper proposes a real-time active exploration technique which is capable of efficiently generating paths that minimize the vehicle’s risk and maximize the coverage area. Furthermore, it accomplishes real-time monitoring of sudden changes in the estimated map, due to the dynamic objects, by reevaluating at real-time the destination and trajectory to minimize the risk on the chosen path and simultaneously preserving the maximization of the coverage area. Ultimately, recording the implemented trajectory of the vehicle also assists in time-saving as the vehicle depends on this trajectory during the exit process. The performance of the technique is studied under static and dynamic environments and is also compared with different algorithms. Full article
(This article belongs to the Special Issue Applications and Potential of UAV Photogrammetric Survey)
Show Figures

Figure 1

Figure 1
<p>Overview of the system structure.</p>
Full article ">Figure 2
<p>Distance transform (DT) method using Euclidean distance.</p>
Full article ">Figure 3
<p>Binary image for indoor environment.</p>
Full article ">Figure 4
<p>The result of the DT method using Euclidean distance.</p>
Full article ">Figure 5
<p>Part of the constructed map I: (<b>a</b>) binary representation; and (<b>b</b>) the corresponding result using the DT method.</p>
Full article ">Figure 6
<p>Part of the constructed map II: (<b>a</b>) binary representation; and (<b>b</b>) the corresponding result using the DT method.</p>
Full article ">Figure 7
<p>Representation of the dynamic movement using successive Gaussian kernels: (<b>a</b>) a long period of the history information; and (<b>b</b>) a short period of the history information.</p>
Full article ">Figure 8
<p>A part of the vehicle’s trajectory representation using Gaussian kernel.</p>
Full article ">Figure 9
<p>A* algorithm.</p>
Full article ">Figure 10
<p>Monitoring of the dynamic objects.</p>
Full article ">Figure 11
<p>Generated path, for Dataset I, using: (<b>a</b>) A* algorithm and APF method; (<b>b</b>) first iteration of the RRT algorithm; (<b>c</b>) second iteration of the RRT algorithm; (<b>d</b>) first iteration of the bidirectional RRT algorithm; and (<b>e</b>) second iteration of the bidirectional RRT algorithm.</p>
Full article ">Figure 12
<p>The cost map of the static representation for Dataset I.</p>
Full article ">Figure 13
<p>The generated path, for Dataset I, within the extension of the static obstacles.</p>
Full article ">Figure 14
<p>2D representation of the generated path, for Dataset I, using the adjusted A* algorithm in a static environment.</p>
Full article ">Figure 15
<p>Representation of the generated and smoothed paths.</p>
Full article ">Figure 16
<p>3D representation of the generated path, for Dataset I, using the adjusted A* algorithm in a static environment.</p>
Full article ">Figure 17
<p>Generated path, for Dataset II, using: (<b>a</b>) A* algorithm; (<b>b</b>) first iteration of APF; (<b>c</b>) second iteration of APF; (<b>d</b>) RRT; and (<b>e</b>) bidirectional RRT.</p>
Full article ">Figure 18
<p>Adjusted A* algorithm in a static environment: (<b>a</b>) cost map; (<b>b</b>) 2D representation; and (<b>c</b>) 3D representation.</p>
Full article ">Figure 19
<p>Generated path, for Dataset III, using: (<b>a</b>) A* algorithm and APF method; (<b>b</b>) RRT algorithm; and (<b>c</b>) bidirectional RRT algorithm.</p>
Full article ">Figure 20
<p>Adjusted A* algorithm in a dynamic environment: (<b>a</b>) cost map; (<b>b</b>) 2D representation; and (<b>c</b>) 3D representation.</p>
Full article ">Figure 21
<p>Generated path, for Dataset IV, using: (<b>a</b>) A* algorithm and APF method; (<b>b</b>) second iteration of APF method; (<b>c</b>) RRT algorithm; and (<b>d</b>) bidirectional RRT algorithm.</p>
Full article ">Figure 22
<p>Adjusted A* algorithm in a static environment: (<b>a</b>) cost map; (<b>b</b>) 2D representation; and (<b>c</b>) 3D representation.</p>
Full article ">Figure 23
<p>Generated path, for Dataset V, using: (<b>a</b>) A* algorithm and APF method; (<b>b</b>) RRT algorithm; and (<b>c</b>) bidirectional RRT algorithm.</p>
Full article ">
13 pages, 6705 KiB  
Article
Estimation of Hourly Link Population and Flow Directions from Mobile CDR
by Ko Ko Lwin, Yoshihide Sekimoto and Wataru Takeuchi
ISPRS Int. J. Geo-Inf. 2018, 7(11), 449; https://doi.org/10.3390/ijgi7110449 - 17 Nov 2018
Cited by 14 | Viewed by 5071
Abstract
The rise in big data applications in urban planning and transport management is now widening and becoming a part of local government decision-making processes. Understanding people flow inside the city helps urban and transport planners build a healthy and lively city. Many flow [...] Read more.
The rise in big data applications in urban planning and transport management is now widening and becoming a part of local government decision-making processes. Understanding people flow inside the city helps urban and transport planners build a healthy and lively city. Many flow maps are based on origin-and-destination points with crossing lines, which reduce the map’s readability and overall appearance. Today, with the emergence of geolocation-enabled handheld devices with wireless communication and networking capabilities, human mobility and the resulting events can be captured and stored as text-based geospatial big data. In this paper, we used one-week mobile-call-detail records (CDR) and a GIS road network model to estimate hourly link population and flow directions, based on mobile-call activities of origin–destination pairs with a shortest-path analysis for the whole city. Moreover, to gain the actual population size from the number of mobile-call users, we introduced a home-based magnification factor (h-MF) by integrating with the national census. Therefore, the final output link data have both magnitude (actual population) and flow direction at one-hour intervals between 06:00 and 21:00. The hourly link population and flow direction dataset are intended to optimize bus routes, solve traffic congestion problems, and enhance disaster and emergency preparedness. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Yangon city road network patterns (left) and census population in 2014 (right).</p>
Full article ">Figure 2
<p>Data processing and research flow. MPT—Myanma Posts and Telecommunications; CDR—call-detail records; OD—origin-and-destination.</p>
Full article ">Figure 3
<p>Formatted MPT CDR for both voice and data with base transceiver station (BTS) locations.</p>
Full article ">Figure 4
<p>Graphical illustration of home-based magnification factor, link counts, population, and flow direction computational steps. (<b>a</b>) Home user extraction; (<b>b</b>) Finding PID with maximum call frequencies; (<b>c</b>) Counting total home users by Cell-ID; (<b>d</b>) Disaggregation of Cell-ID population with its home users; (<b>e</b>) Illustration of link counts, population and flow directions computation.</p>
Full article ">Figure 5
<p>Moving cell tower location to nearest road node to synchronize the directions and shortest path analysis. (<b>a</b>) Before moving Cell-ID; (<b>b</b>) After moving Cell-ID; (<b>c</b>) find the shortest path between two successive calls.</p>
Full article ">Figure 6
<p>Link population (magnitude) between 06:00 and 07:00.</p>
Full article ">Figure 7
<p>Link population (magnitude) between 17:00 and 18:00.</p>
Full article ">Figure 8
<p>Hourly link flow direction between 06:00 and 07:00.</p>
Full article ">Figure 9
<p>Hourly link flow direction between 17:00 and 18:00.</p>
Full article ">Figure 10
<p>(<b>a</b>) Hourly link flow direction with population values for both directions between 06:00 and 07:00. (<b>b</b>) Comparison of hourly flow magnitude for both southbound and northbound at Hledan Junction. (<b>c</b>) Ground-truth data collection for result validation.</p>
Full article ">
16 pages, 2593 KiB  
Article
Toward Model-Generated Household Listing in Low- and Middle-Income Countries Using Deep Learning
by Robert Chew, Kasey Jones, Jennifer Unangst, James Cajka, Justine Allpress, Safaa Amer and Karol Krotki
ISPRS Int. J. Geo-Inf. 2018, 7(11), 448; https://doi.org/10.3390/ijgi7110448 - 16 Nov 2018
Cited by 8 | Viewed by 4530
Abstract
While governments, researchers, and NGOs are exploring ways to leverage big data sources for sustainable development, household surveys are still a critical source of information for dozens of the 232 indicators for the Sustainable Development Goals (SDGs) in low- and middle-income countries (LMICs). [...] Read more.
While governments, researchers, and NGOs are exploring ways to leverage big data sources for sustainable development, household surveys are still a critical source of information for dozens of the 232 indicators for the Sustainable Development Goals (SDGs) in low- and middle-income countries (LMICs). Though some countries’ statistical agencies maintain databases of persons or households for sampling, conducting household surveys in LMICs is complicated due to incomplete, outdated, or inaccurate sampling frames. As a means to develop or update household listings in LMICs, this paper explores the use of machine learning models to detect and enumerate building structures directly from satellite imagery in the Kaduna state of Nigeria. Specifically, an object detection model was used to identify and locate buildings in satellite images. In the test set, the model attained a mean average precision (mAP) of 0.48 for detecting structures, with relatively higher values in areas with lower building density (mAP = 0.65). Furthermore, when model predictions were compared against recent household listings from fieldwork in Nigeria, the predictions showed high correlation with household coverage (Pearson = 0.70; Spearman = 0.81). With the need to produce comparable, scalable SDG indicators, this case study explores the feasibility and challenges of using object detection models to help develop timely enumerated household lists in LMICs. Full article
(This article belongs to the Special Issue Geo-Information and the Sustainable Development Goals (SDGs))
Show Figures

Figure 1

Figure 1
<p>Example of an 100 × 100 m image annotated with human-labeled bounding boxes.</p>
Full article ">Figure 2
<p>Predicted boxes: low building density.</p>
Full article ">Figure 3
<p>Predicted boxes: high building density.</p>
Full article ">Figure 4
<p>Scatterplot of household counts vs. predicted building counts with a 1:1 diagonal line.</p>
Full article ">Figure 5
<p>Difference between Alive and Thrive (A&amp;T) survey household counts and predicted buildings per SGU vs. (<b>a</b>) household counts per SGU and (<b>b</b>) predicted building counts per SGU.</p>
Full article ">Figure 6
<p>(<b>a</b>) SGC grid area with 40 predicted buildings and 6 households; (<b>b</b>) SGC grid area with 48 predicted buildings and 91 households.</p>
Full article ">
22 pages, 7952 KiB  
Article
Modelling, Validation and Quantification of Climate and Other Sensitivities of Building Energy Model on 3D City Models
by Syed Monjur Murshed, Solène Picard and Andreas Koch
ISPRS Int. J. Geo-Inf. 2018, 7(11), 447; https://doi.org/10.3390/ijgi7110447 - 15 Nov 2018
Cited by 13 | Viewed by 5306
Abstract
New planning tools are required to depict the complete building stock in a city and investigate detailed measures on reaching local and global targets to improve energy efficiency and reduce greenhouse gas emissions. To pursue this objective, ISO (the International Organization for Standardization) [...] Read more.
New planning tools are required to depict the complete building stock in a city and investigate detailed measures on reaching local and global targets to improve energy efficiency and reduce greenhouse gas emissions. To pursue this objective, ISO (the International Organization for Standardization) 13790:2008 monthly heating and cooling energy calculation method is implemented using geometric information from 3D city models (e.g., CityGML format) in an open source software architecture. A model is developed and applied in several urban districts with different number of 3D buildings in various cities. The model is validated with the simulation software TRNSYS. We also perform a sensitivity analysis to quantify the impact of climate change and other physical and behavioral factors on modelling results. The proposed approach can help to perform city or district-wide analysis of the building energy needs and prepare different renovation plans to support decision-making, which finally will enhance the livability of a city and the quality of life of the citizens. Full article
Show Figures

Figure 1

Figure 1
<p>Overall research methodology of this study.</p>
Full article ">Figure 2
<p>Overview of different software, tools and datasets required to implement CityBEM model (modified after Reference [<a href="#B30-ijgi-07-00447" class="html-bibr">30</a>]).</p>
Full article ">Figure 3
<p>Implementation steps of the CityBEM model: Different input data, related software packages and the tables created in the databases.</p>
Full article ">Figure 4
<p>Illustration of annual heating energy need (kWh/m<sup>2</sup>) in a district in the city of Karlsruhe in Germany.</p>
Full article ">Figure 5
<p>Annual specific building heating and cooling energy needs in 33 building typologies and corresponding number of buildings in each typology (SFH = Single Family House, DFH = Double Family House, MFH = Multi Family House, AB = Apartment Block, HRAB = High Rise Apartment Block, OB = Office Building, WB = Workshop Building).</p>
Full article ">Figure 6
<p>Monthly specific (<b>a</b>) heating and (<b>b</b>) cooling energy needs for seven building types and (<b>c</b>) heating energy needs for five building age classes.</p>
Full article ">Figure 6 Cont.
<p>Monthly specific (<b>a</b>) heating and (<b>b</b>) cooling energy needs for seven building types and (<b>c</b>) heating energy needs for five building age classes.</p>
Full article ">Figure 7
<p>The CityGML building extracted for simulation in TRNSYS (<b>left</b>) and corresponding geometric properties (<b>right</b>).</p>
Full article ">Figure 8
<p>Comparison of energy simulation results obtained with TRNSYS and ISO-based CityBEM for an office building (1975).</p>
Full article ">Figure 9
<p>Comparison of energy simulation results obtained with TRNSYS and ISO-based CityBEM for a residential building (1985).</p>
Full article ">Figure 10
<p>Sensitivity analyses approach for the CityBEM monthly model.</p>
Full article ">Figure 11
<p>Sensitivity of U-values on monthly energy needs for heating (<b>left</b>) and cooling (<b>right</b>).</p>
Full article ">Figure 12
<p>Sensitivity of set point temperature on monthly energy needs for heating (<b>left</b>) and cooling (<b>right</b>).</p>
Full article ">Figure 13
<p>Sensitivity indices for monthly heating and cooling energy needs considering the local method.</p>
Full article ">Figure 14
<p>Decrease of annual heating energy needs in 2050 (<b>left</b>) and in 2100 (<b>right</b>), compared to the reference scenarios.</p>
Full article ">Figure 15
<p>Increase of annual cooling energy needs in 2050 (<b>left</b>) and in 2100 (<b>right</b>), compared to the reference scenarios.</p>
Full article ">Figure 16
<p>Impact of future climate change on specific cooling energy needs of different building typologies in an urban district in Karlsruhe.</p>
Full article ">Figure 17
<p>Impact of future climate change on specific heating energy needs of different building typologies in an urban district in Karlsruhe.</p>
Full article ">
18 pages, 2070 KiB  
Article
Building a Framework of Usability Patterns for Web Applications in Spatial Data Infrastructures
by Christin Henzen
ISPRS Int. J. Geo-Inf. 2018, 7(11), 446; https://doi.org/10.3390/ijgi7110446 - 15 Nov 2018
Cited by 5 | Viewed by 4122
Abstract
Web applications in spatial data infrastructures (SDIs) should provide robust and user-friendly user interfaces for geoinformation (GI) discovery, analysis, and usage. Poor usability, e.g., caused by unsuitable information presentation or inappropriate (non) availability of functions, can result in inefficient or faulty usage and [...] Read more.
Web applications in spatial data infrastructures (SDIs) should provide robust and user-friendly user interfaces for geoinformation (GI) discovery, analysis, and usage. Poor usability, e.g., caused by unsuitable information presentation or inappropriate (non) availability of functions, can result in inefficient or faulty usage and can increase the acceptance of the application and provided geoinformation. Until now, a number of usability problems in GI web applications were identified; however, methods to summarize these problems, to provide (software-independent) solutions for them, and to find pairs of problems and related solutions hardly exist. We propose an adapted usability pattern concept for web applications in SDIs to map and categorize usability problems and best practice solutions and we enable a GI context-specific creation and discovery of these problems and solutions. The concept includes developed pattern types, relationships, and rules on how to use the relationships for the different pattern types. Full article
Show Figures

Figure 1

Figure 1
<p>Aspects of usability in geoinformation (GI) web applications.</p>
Full article ">Figure 2
<p>Sub-attributes of the pattern context.</p>
Full article ">Figure 3
<p>Eight usability patterns for GI web applications and their relationships.</p>
Full article ">
16 pages, 4334 KiB  
Article
Species-Level Vegetation Mapping in a Himalayan Treeline Ecotone Using Unmanned Aerial System (UAS) Imagery
by Niti B. Mishra, Kumar P. Mainali, Bharat B. Shrestha, Jackson Radenz and Debendra Karki
ISPRS Int. J. Geo-Inf. 2018, 7(11), 445; https://doi.org/10.3390/ijgi7110445 - 14 Nov 2018
Cited by 29 | Viewed by 5890
Abstract
Understanding ecological patterns and response to climate change requires unbiased data on species distribution. This can be challenging, especially in biodiverse but extreme environments like the Himalaya. This study presents the results of the first ever application of Unmanned Aerial Systems (UAS) imagery [...] Read more.
Understanding ecological patterns and response to climate change requires unbiased data on species distribution. This can be challenging, especially in biodiverse but extreme environments like the Himalaya. This study presents the results of the first ever application of Unmanned Aerial Systems (UAS) imagery for species-level mapping of vegetation in the Himalaya following a hierarchical Geographic Object Based Image Analysis (GEOBIA) method. The first level of classification separated green vegetated objects from the rest with overall accuracy of 95%. At the second level, seven cover types were identified (including four woody vegetation species). For this, the suitability of various spectral, shape and textural features were tested for classifying them using an ensemble decision tree algorithm. Spectral features alone yielded ~70% accuracy (kappa 0.66) whereas adding textural and shape features marginally improved the accuracy (73%) but at the cost of a substantial increase in processing time. Contrast in plant morphological traits was the key to distinguishing nearby stands as different species. Hence, broad-leaved versus fine needle leaved vegetation were mapped more accurately than structurally similar classes such as Rhododendron anthopogon versus non-photosynthetic vegetation. Results highlight the potential and limitations of the suggested UAS-GEOBIA approach for detailed mapping of plant communities and suggests future research directions. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Study area location in Nepal Himalaya showing elevation range; (<b>b</b>) location of study area within Langtang National Park, Nepal; (<b>c</b>) field vegetation sampling locations overlaid on multispectral orthomosaic for the study area created from UAS imagery; and (<b>d</b>) ground control point.</p>
Full article ">Figure 2
<p>(<b>a</b>) Calibration panel used for estimating surface reflectance; (<b>b</b>) Sequoia camera that acquires images at fixed wavelengths including red-edge and near-infrared; (<b>c</b>) the entire UAS setup with remote and tablet used for mission planning; and (<b>d</b>) setting up ground control points (GCPs) prior to UAS mission in the treeline ecotone in Langtang national park (image date June, 2016). Also labeled are the woody species mapped in this study.</p>
Full article ">Figure 3
<p>Workflow of data collection, pre-processing and analysis to produce a species classification map.</p>
Full article ">Figure 4
<p>Results of the segmentation method used in this study (background image band combination: Blue: green band, Green: red band and Red: red edge band).</p>
Full article ">Figure 5
<p>GEOBIA classification outputs for level 1 and level 2 classification.</p>
Full article ">Figure 6
<p>Variable importance of various feature used for random forest classification. Only the top 15 variables are shown.</p>
Full article ">
22 pages, 16957 KiB  
Article
Multi-Criteria Decision Making (MCDM) Model for Seismic Vulnerability Assessment (SVA) of Urban Residential Buildings
by Mohsen Alizadeh, Mazlan Hashim, Esmaeil Alizadeh, Himan Shahabi, Mohammad Reza Karami, Amin Beiranvand Pour, Biswajeet Pradhan and Hassan Zabihi
ISPRS Int. J. Geo-Inf. 2018, 7(11), 444; https://doi.org/10.3390/ijgi7110444 - 14 Nov 2018
Cited by 53 | Viewed by 9088
Abstract
Earthquakes are among the most catastrophic natural geo-hazards worldwide and endanger numerous lives annually. Therefore, it is vital to evaluate seismic vulnerability beforehand to decrease future fatalities. The aim of this research is to assess the seismic vulnerability of residential houses in an [...] Read more.
Earthquakes are among the most catastrophic natural geo-hazards worldwide and endanger numerous lives annually. Therefore, it is vital to evaluate seismic vulnerability beforehand to decrease future fatalities. The aim of this research is to assess the seismic vulnerability of residential houses in an urban region on the basis of the Multi-Criteria Decision Making (MCDM) model, including the analytic hierarchy process (AHP) and geographical information system (GIS). Tabriz city located adjacent to the North Tabriz Fault (NTF) in North-West Iran was selected as a case study. The NTF is one of the major seismogenic faults in the north-western part of Iran. First, several parameters such as distance to fault, percent of slope, and geology layers were used to develop a geotechnical map. In addition, the structural construction materials, building materials, size of building blocks, quality of buildings and buildings-floors were used as key factors impacting on the building’s structural vulnerability in residential areas. Subsequently, the AHP technique was adopted to measure the priority ranking, criteria weight (layers), and alternatives (classes) of every criterion through pair-wise comparison at all levels. Lastly, the layers of geotechnical and spatial structures were superimposed to design the seismic vulnerability map of buildings in the residential area of Tabriz city. The results showed that South and Southeast areas of Tabriz city exhibit low to moderate vulnerability, while some regions of the north-eastern area are under severe vulnerability conditions. In conclusion, the suggested approach offers a practical and effective evaluation of Seismic Vulnerability Assessment (SVA) and provides valuable information that could assist urban planners during mitigation and preparatory phases of less examined areas in many other regions around the world. Full article
(This article belongs to the Special Issue Natural Hazards and Geospatial Information)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Geographical location of the East Azarbaijan province and Tabriz city in relation to the political provincial and national border.</p>
Full article ">Figure 2
<p>Flowchart of the methodology illustrating the various stages of the analysis for the preparation of the seismic vulnerability map of the study area.</p>
Full article ">Figure 3
<p>Distribution of (<b>a</b>) Fault systems (<b>b</b>) Slope gradients and (<b>c</b>) Geology map of Tabriz city. Mmg2 = Interlayer of greenish grey marl associated with an interlayer of gypsum- bring sandy marl. Msc5 = Interbedded red conglomerate with sandstone and red marl. Msm4 = Sandstone and red marl. Pldt = Diatomic and fish interbedded with fine particles sediment. Plqc = Interlayer of semi-hard conglomerate associated with sandstone and pumice. Plqc = Interlayer of semi-hard conglomerate associated with sandstone and pumice. Qal = Quaternary alluvium. Qt2 = Young terrace and alluvium deposits.</p>
Full article ">Figure 4
<p>Distribution of (<b>a</b>) buildings materials (<b>b</b>) size of buildings density (<b>c</b>) quality of buildings and (<b>d</b>) buildings floor in Tabriz city.</p>
Full article ">Figure 4 Cont.
<p>Distribution of (<b>a</b>) buildings materials (<b>b</b>) size of buildings density (<b>c</b>) quality of buildings and (<b>d</b>) buildings floor in Tabriz city.</p>
Full article ">Figure 5
<p>Residential buildings vulnerability distribution showing geotechnical vulnerability in Tabriz municipality.</p>
Full article ">Figure 6
<p>Residential buildings vulnerability distribution showing structural vulnerability in Tabriz municipality.</p>
Full article ">Figure 7
<p>Residential buildings vulnerability distribution showing structural and geotechnical vulnerability in Tabriz municipality.</p>
Full article ">
15 pages, 12105 KiB  
Article
Effective Identification of Terrain Positions from Gridded DEM Data Using Multimodal Classification Integration
by Ling Jiang, Dequan Ling, Mingwei Zhao, Chun Wang, Qiuhua Liang and Kai Liu
ISPRS Int. J. Geo-Inf. 2018, 7(11), 443; https://doi.org/10.3390/ijgi7110443 - 14 Nov 2018
Cited by 16 | Viewed by 5616
Abstract
Terrain positions are widely used to describe the Earth’s topographic features and play an important role in the studies of landform evolution, soil erosion and hydrological modeling. This work develops a new multimodal classification system with enhanced classification performance by integrating different approaches [...] Read more.
Terrain positions are widely used to describe the Earth’s topographic features and play an important role in the studies of landform evolution, soil erosion and hydrological modeling. This work develops a new multimodal classification system with enhanced classification performance by integrating different approaches for terrain position identification. The adopted classification approaches include local terrain attribute (LA)-based and regional terrain attribute (RA)-based, rule-based and supervised, and pixel-based and object-oriented methods. Firstly, a double-level definition scheme is presented for terrain positions. Then, utilizing a hierarchical framework, a multimodal approach is developed by integrating different classification techniques. Finally, an assessment method is established to evaluate the new classification system from different aspects. The experimental results, obtained at a Loess Plateau region in northern China on a 5 m digital elevation model (DEM), show reasonably positional relationship, and larger inter-class and smaller intra-class variances. This indicates that identified terrain positions are consistent with the actual topography from both overall and local perspectives, and have relatively good integrity and rationality. This study demonstrates that the current multimodal classification system, developed by taking advantage of various classification methods, can reflect the geographic meanings and topographic features of terrain positions from different levels. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location and digital elevation model (DEM) of the study area.</p>
Full article ">Figure 2
<p>Flowchart of the proposed classification method.</p>
Full article ">Figure 3
<p>Illustration of terrain positions of hills (Modified from Ruhe [<a href="#B38-ijgi-07-00443" class="html-bibr">38</a>]).</p>
Full article ">Figure 4
<p>Distribution of RPI.</p>
Full article ">Figure 5
<p>Relationship between the average slope and RPI: (<b>a</b>) lower-slope; (<b>b</b>) upper-slope.</p>
Full article ">Figure 6
<p>Terrain positions at the first level.</p>
Full article ">Figure 7
<p>Patch skeleton lines (local area).</p>
Full article ">Figure 8
<p>Classification of hillocks and ridges.</p>
Full article ">Figure 9
<p>Terrain positions of the proposed classification method.</p>
Full article ">Figure 10
<p>Terrain positions classified by the pixel-based classification method using the rules at the second level.</p>
Full article ">
28 pages, 27783 KiB  
Article
Complying with Privacy Legislation: From Legal Text to Implementation of Privacy-Aware Location-Based Services
by Mehrnaz Ataei, Auriol Degbelo, Christian Kray and Vitor Santos
ISPRS Int. J. Geo-Inf. 2018, 7(11), 442; https://doi.org/10.3390/ijgi7110442 - 13 Nov 2018
Cited by 9 | Viewed by 4803
Abstract
An individual’s location data is very sensitive geoinformation. While its disclosure is necessary, e.g., to provide location-based services (LBS), it also facilitates deep insights into the lives of LBS users as well as various attacks on these users. Location privacy threats can be [...] Read more.
An individual’s location data is very sensitive geoinformation. While its disclosure is necessary, e.g., to provide location-based services (LBS), it also facilitates deep insights into the lives of LBS users as well as various attacks on these users. Location privacy threats can be mitigated through privacy regulations such as the General Data Protection Regulation (GDPR), which was introduced recently and harmonises data privacy laws across Europe. While the GDPR is meant to protect users’ privacy, the main problem is that it does not provide explicit guidelines for designers and developers about how to build systems that comply with it. In order to bridge this gap, we systematically analysed the legal text, carried out expert interviews, and ran a nine-week-long take-home study with four developers. We particularly focused on user-facing issues, as these have received little attention compared to technical issues. Our main contributions are a list of aspects from the legal text of the GDPR that can be tackled at the user interface level and a set of guidelines on how to realise this. Our results can help service providers, designers and developers of applications dealing with location information from human users to comply with the GDPR. Full article
(This article belongs to the Special Issue Recent Trends in Location Based Services and Science)
Show Figures

Figure 1

Figure 1
<p>Words mentioned at least 200 times in the GDPR reference document [<a href="#B5-ijgi-07-00442" class="html-bibr">5</a>] (stopwords removed). It helps see most prominent concepts of the regulation: personal data, the processing of it, and the interaction between controllers and (data) subjects. Word frequencies were obtained by using Wordart [<a href="#B32-ijgi-07-00442" class="html-bibr">32</a>].</p>
Full article ">Figure 2
<p>Recommended guidelines based on legal requirements and expert suggestions.</p>
Full article ">Figure 3
<p>GeoFreebie (left to right): (<b>a</b>) approximate location of the users shown in a blue circle; (<b>b</b>) marked location info as private; (<b>c</b>) setting options to adjust location tracking and location sharing; (<b>d</b>) notification pop-up if the location sharing and tracking are enabled.</p>
Full article ">Figure 4
<p>TourChamp (left to right): (<b>a</b>) notice; (<b>b</b>) visual indications of friends’ presence when location sharing is enabled; (<b>c</b>) setting options to adjust location sharing (GPS) for public or only friends; (<b>d</b>) layered setting adjustment to enable/disable location sharing.</p>
Full article ">Figure 5
<p>Recommended guidelines reflected in the two LBS developed during the take-home study: “G” and “T” indicate factors that were realised in the GeoFreebie and TourChamp LBS, respectively.</p>
Full article ">Figure 6
<p>The correspondence between the guidelines and the original GDPR document.</p>
Full article ">
27 pages, 22984 KiB  
Article
Change Detection in Coral Reef Environment Using High-Resolution Images: Comparison of Object-Based and Pixel-Based Paradigms
by Zhenjin Zhou, Lei Ma, Tengyu Fu, Ge Zhang, Mengru Yao and Manchun Li
ISPRS Int. J. Geo-Inf. 2018, 7(11), 441; https://doi.org/10.3390/ijgi7110441 - 12 Nov 2018
Cited by 31 | Viewed by 6163
Abstract
Despite increases in the spatial resolution of satellite imagery prompting interest in object-based image analysis, few studies have used object-based methods for monitoring changes in coral reefs. This study proposes a high accuracy object-based change detection (OBCD) method intended for coral reef environment, [...] Read more.
Despite increases in the spatial resolution of satellite imagery prompting interest in object-based image analysis, few studies have used object-based methods for monitoring changes in coral reefs. This study proposes a high accuracy object-based change detection (OBCD) method intended for coral reef environment, which uses QuickBird and WorldView-2 images. The proposed methodological framework includes image fusion, multi-temporal image segmentation, image differencing, random forests models, and object-area-based accuracy assessment. For validation, we applied the method to images of four coral reef study sites in the South China Sea. We compared the proposed OBCD method with a conventional pixel-based change detection (PBCD) method by implementing both methods under the same conditions. The average overall accuracy of OBCD exceeded 90%, which was approximately 20% higher than PBCD. The OBCD method was free from salt-and-pepper effects and was less prone to images misregistration in terms of change detection accuracy and mapping results. The object-area-based accuracy assessment reached a higher overall accuracy and per-class accuracy than the object-number-based and pixel-number-based accuracy assessment. Full article
(This article belongs to the Special Issue GEOBIA in a Changing World)
Show Figures

Figure 1

Figure 1
<p>Four coral reef study sites: (<b>a</b>,<b>b</b>) QuickBird satellite images of Taiping Island and Zhongye Island (Bands 3, 2, and 1 in RGB); (<b>c</b>,<b>d</b>) WorldView-2 satellite images of Barque Canada Reef (Bands 5, 3, and 2 in RGB), where red rectangles indicate the two study sites in Barque Canada reef.</p>
Full article ">Figure 2
<p>Schematic of the procedures implemented in coral reef change detection using the pixel-based change detection (PBCD) and object-based change detection (OBCD) methods and the comparison of their results.</p>
Full article ">Figure 3
<p>Rate of change of local variance (ROC-LV) curve of multiresolution segmentation of Zhongye Island images.</p>
Full article ">Figure 4
<p>Change detection maps and reference maps of Zhongye Island and Taiping Island.</p>
Full article ">Figure 5
<p>Change detection maps and reference maps of the two study sites on Barque Canada Reef.</p>
Full article ">Figure 6
<p>Overall accuracy and Kappa coefficients of the PBCD and OBCD results of the four study sites.</p>
Full article ">Figure 7
<p>The impact of registration error on PBCD and OBCD change detection methods. (Pixels between neighboring segments were marked by red circle in PBCD change detection maps).</p>
Full article ">
15 pages, 7419 KiB  
Article
Low-Power LoRa Signal-Based Outdoor Positioning Using Fingerprint Algorithm
by Wongeun Choi, Yoon-Seop Chang, Yeonuk Jung and Junkeun Song
ISPRS Int. J. Geo-Inf. 2018, 7(11), 440; https://doi.org/10.3390/ijgi7110440 - 9 Nov 2018
Cited by 64 | Viewed by 8558
Abstract
Positioning is an essential element in most Internet of Things (IoT) applications. Global Positioning System (GPS) chips have high cost and power consumption, making it unsuitable for long-range (LoRa) and low-power IoT devices. Alternatively, low-power wide-area (LPWA) signals can be used for simultaneous [...] Read more.
Positioning is an essential element in most Internet of Things (IoT) applications. Global Positioning System (GPS) chips have high cost and power consumption, making it unsuitable for long-range (LoRa) and low-power IoT devices. Alternatively, low-power wide-area (LPWA) signals can be used for simultaneous positioning and communication. We summarize previous studies related to LoRa signal-based positioning systems, including those addressing proximity, a path loss model, time difference of arrival (TDoA), and fingerprint positioning methods. We propose a LoRa signal-based positioning method that uses a fingerprint algorithm instead of a received signal strength indicator (RSSI) proximity or TDoA method. The main objective of this study was to evaluate the accuracy and usability of the fingerprint algorithm for large areas in the real world. We estimated the locations using probabilistic means based on three different algorithms that use interpolated fingerprint RSSI maps. The average accuracy of the three proposed algorithms in our experiments was 28.8 m. Our method also reduced the battery consumption significantly compared with that of existing GPS-based positioning methods. Full article
(This article belongs to the Special Issue Geospatial Applications of the Internet of Things (IoT))
Show Figures

Figure 1

Figure 1
<p>Telecommunication methods [<a href="#B2-ijgi-07-00440" class="html-bibr">2</a>].</p>
Full article ">Figure 2
<p>Comparison of positioning technologies [<a href="#B6-ijgi-07-00440" class="html-bibr">6</a>].</p>
Full article ">Figure 3
<p>Link budget calculation to measure RSSI values [<a href="#B2-ijgi-07-00440" class="html-bibr">2</a>].</p>
Full article ">Figure 4
<p>Semtech Corporation’s TDoA-based geolocation architecture [<a href="#B6-ijgi-07-00440" class="html-bibr">6</a>].</p>
Full article ">Figure 5
<p>Comparisons of packet flow between Wi-Fi and LoRa positioning methods: (<b>a</b>) downlink packet from AP to smartphone and (<b>b</b>) uplink packet from end-device to LoRa gateways.</p>
Full article ">Figure 6
<p>Positioning using LoRa fingerprint.</p>
Full article ">Figure 7
<p>LoRa RSSI data acquisition and fingerprint map generation.</p>
Full article ">Figure 8
<p>Interpolation of f(x,y) using Gaussian RBF with 25 sample points: (<b>a</b>) Random data points. (<b>b</b>) the RBF centered at the collocation points (<b>c</b>) shows the interpolated surface [<a href="#B21-ijgi-07-00440" class="html-bibr">21</a>].</p>
Full article ">Figure 9
<p>Comparison of interpolated fingerprint maps: (<b>a</b>) cubic, (<b>b</b>) Gaussian, (<b>c</b>) quintic, (<b>d</b>) linear, and (<b>e</b>) thin plate interpolation methods.</p>
Full article ">Figure 10
<p>Probability map generation.</p>
Full article ">Figure 11
<p>Comparison of candidate regions. We filter candidate regions with 3 dBm. Dark areas represent candidates region with (<b>a</b>) strong and (<b>b</b>) weak RSSI signals.</p>
Full article ">Figure 12
<p>Signal-strength distribution at a fixed location. [<a href="#B22-ijgi-07-00440" class="html-bibr">22</a>].</p>
Full article ">Figure 13
<p>Vector space RSSI data collection.</p>
Full article ">Figure 14
<p>Experimental area (pink: gateway; green: training dataset; red: test dataset).</p>
Full article ">Figure 15
<p>Cumulative probability distribution of errors.</p>
Full article ">Figure 16
<p>Candidate regions extracted with buffer.</p>
Full article ">Figure 17
<p>Probability maps and errors for three algorithms.</p>
Full article ">
15 pages, 8183 KiB  
Article
Automatic Parametrization of Urban Areas Using ALS Data: The Case Study of Santiago de Compostela
by Mario Soilán, Belén Riveiro, Patricia Liñares and Andrea Pérez-Rivas
ISPRS Int. J. Geo-Inf. 2018, 7(11), 439; https://doi.org/10.3390/ijgi7110439 - 9 Nov 2018
Cited by 5 | Viewed by 3659
Abstract
Nowadays, gathering accurate and meaningful information about the urban environment with the maximum efficiency in terms of cost and time has become more relevant for city administrations, as this information is essential if the sustainability or the resilience of the urban structure has [...] Read more.
Nowadays, gathering accurate and meaningful information about the urban environment with the maximum efficiency in terms of cost and time has become more relevant for city administrations, as this information is essential if the sustainability or the resilience of the urban structure has to be improved. This work presents a methodology for the automatic parametrization and characterization of different urban typologies, for the specific case study of Santiago de Compostela (Spain), using data from Aerial Laser Scanners (ALS). This methodology consists of a number of sequential processes of point cloud data, using exclusively their geometric coordinates. Three of the main elements of the urban structure are assessed in this work: intersections, building blocks, and streets. Different geometric and contextual metrics are automatically extracted for each of the elements, defining the urban typology of the studied area. The accuracy of the measurements is validated against a manual reference, obtaining average errors of less than 3%, proving that the input data is valid for this assessment. Full article
Show Figures

Figure 1

Figure 1
<p>Methodology workflow.</p>
Full article ">Figure 2
<p>Point cloud classification. (<b>a</b>) Point cloud coloured by elevation. (<b>b</b>) Elevation-based filtering. Points that do not belong to roofs or to the ground are removed from the point cloud. (<b>c</b>) After performing the triangulation <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mrow> <mi>P</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math>, two groups of points are labelled as either ground or roof. (<b>d</b>) A region growing algorithm classifies the complete point cloud.</p>
Full article ">Figure 3
<p>Building block boundaries. (<b>a</b>) Raster image of the point cloud <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="bold-italic">P</mi> <mi>r</mi> </msub> </mrow> </semantics></math> (<math display="inline"><semantics> <mrow> <msub> <mi mathvariant="bold-italic">I</mi> <mi>b</mi> </msub> <mo stretchy="false">)</mo> </mrow> </semantics></math>. (<b>b</b>) Image <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <msub> <mi mathvariant="bold-italic">I</mi> <mi>b</mi> </msub> </mrow> <mo stretchy="true">¯</mo> </mover> </mrow> </semantics></math>, defined to remove noisy elements. (<b>c</b>) External and internal building boundaries, coloured in red and green, respectively.</p>
Full article ">Figure 4
<p>Triangulation <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="bold-italic">T</mi> <mi mathvariant="bold-italic">B</mi> </msub> </mrow> </semantics></math>, where each triangle is coloured according to the number of building blocks that are in contact with its vertices. Red triangles are in contact with 3 blocks (intersections), green triangles are in contact with two blocks (streets) and blue triangles are in contact with only one block (building blocks).</p>
Full article ">Figure 5
<p>Urban parametrization. (<b>a</b>) Façade orientation labels. (<b>b</b>) The street width (<math display="inline"><semantics> <mrow> <msub> <mi>W</mi> <mi>s</mi> </msub> </mrow> </semantics></math>) is obtained computing the distance between each point in <math display="inline"><semantics> <mrow> <msub> <mi>P</mi> <mrow> <mi>B</mi> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> and the closest point in <math display="inline"><semantics> <mrow> <msub> <mi>P</mi> <mrow> <mi>B</mi> <mi>j</mi> </mrow> </msub> </mrow> </semantics></math>. The street axis is obtained averaging each pair of those points. (<b>c</b>) Street orientation labels.</p>
Full article ">Figure 6
<p>Three different classes with their respective properties are defined for storing information of the intersections, building blocks, and streets.</p>
Full article ">Figure 7
<p>Case study data. Three well-differentiated urban areas from Santiago de Compostela (highlighted in red, green and blue) are chosen for the validation of the proposed methodology.</p>
Full article ">Figure 8
<p>Comparison of intersection surfaces for the three case study areas.</p>
Full article ">Figure 9
<p>(<b>a</b>) Comparison of building block surface and built surface for the three case study areas. (<b>b</b>) Although individual buildings are smaller in the city centre, in many cases they are built without physical separation, hence they are considered the same building block, having an impact on the interpretation of the results.</p>
Full article ">Figure 10
<p>(<b>a</b>) Comparison of building block heights for the three case study areas, (<b>b</b>) comparison of façade orientation, having in account the four orientation labels that were defined for each point.</p>
Full article ">Figure 11
<p>(<b>a</b>) Comparison of street areas. (<b>b</b>) Comparison of street widths.</p>
Full article ">Figure 12
<p>Manual reference. (<b>a</b>) Area of a building block. (<b>b</b>) Area of an intersection. (<b>c</b>) Width of a street.</p>
Full article ">Figure 13
<p>Visualization of the results as a GIS shape layer which overlaps an orthophoto of the case study area. (<b>a</b>) Building block parameters can be visualized either point-wise (azimuth and orientation label) or block wise (area, volume, etc.). (<b>b</b>) Street parameters can be similarly visualized.</p>
Full article ">
29 pages, 16129 KiB  
Article
Landslide Susceptibility Mapping Using Logistic Regression Analysis along the Jinsha River and Its Tributaries Close to Derong and Deqin County, Southwestern China
by Xiaohui Sun, Jianping Chen, Yiding Bao, Xudong Han, Jiewei Zhan and Wei Peng
ISPRS Int. J. Geo-Inf. 2018, 7(11), 438; https://doi.org/10.3390/ijgi7110438 - 8 Nov 2018
Cited by 85 | Viewed by 8278
Abstract
The objective of this study was to identify the areas that are most susceptible to landslide occurrence, and to find the key factors associated with landslides along Jinsha River and its tributaries close to Derong and Deqin County. Thirteen influencing factors, including (a) [...] Read more.
The objective of this study was to identify the areas that are most susceptible to landslide occurrence, and to find the key factors associated with landslides along Jinsha River and its tributaries close to Derong and Deqin County. Thirteen influencing factors, including (a) lithology, (b) slope angle, (c) slope aspect, (d) TWI, (e) curvature, (f) SPI, (g) STI, (h) topographic relief, (i) rainfall, (j) vegetation, (k) NDVI, (l) distance-to-river, (m) and distance-to-fault, were selected as the landslide conditioning factors in landslide susceptibility mapping. These factors were mainly obtained from the field survey, digital elevation model (DEM), and Landsat 4–5 imagery using ArcGIS software. A total of 40 landslides were identified in the study area from field survey and aerial photos’ interpretation. First, the frequency ratio (FR) method was used to clarify the relationship between the landslide occurrence and the influencing factors. Then, the principal component analysis (PCA) was used to eliminate multiple collinearities between the 13 influencing factors and to reduce the dimension of the influencing factors. Subsequently, the factors that were reselected using the PCA were introduced into the logistic regression analysis to produce the landslide susceptibility map. Finally, the receiver operating characteristic (ROC) curve was used to evaluate the accuracy of the logistic regression analysis model. The landslide susceptibility map was divided into the following five classes: very low, low, moderate, high, and very high. The results showed that the ratios of the areas of the five susceptibility classes were 23.14%, 22.49%, 18.00%, 19.08%, and 17.28%, respectively. And the prediction accuracy of the model was 83.4%. The results were also compared with the FR method (79.9%) and the AHP method (76.9%), which meant that the susceptibility model was reasonable. Finally, the key factors of the landslide occurrence were determined based on the above results. Consequently, this study could serve as an effective guide for further land use planning and for the implementation of development. Full article
(This article belongs to the Special Issue Natural Hazards and Geospatial Information)
Show Figures

Figure 1

Figure 1
<p>Geographical position and landslide inventory of the study area.</p>
Full article ">Figure 2
<p>Flow chart of this study.</p>
Full article ">Figure 3
<p>The relationship between elevation and average annual precipitation based on the nine precipitation stations.</p>
Full article ">Figure 4
<p>Influencing factors maps of the study area: (<b>a</b>) lithology; (<b>b</b>) slope angle; (<b>c</b>) slope aspect; and (<b>d</b>) topographic wetness index (TWI).</p>
Full article ">Figure 5
<p>Influencing factor maps of the study area: (<b>a</b>) curvature; (<b>b</b>) steam power index (SPI); (<b>c</b>) sediment transport index (STI); and (<b>d</b>) topographic relief.</p>
Full article ">Figure 6
<p>Influencing factor maps of the study area: (<b>a</b>) rainfall; (<b>b</b>) vegetation; (<b>c</b>) normalized difference vegetation index (NDVI); and (<b>d</b>) distance-to-river.</p>
Full article ">Figure 7
<p>Influencing factors maps of the study area: (<b>a</b>) distance-to-fault and (<b>b</b>) elevation.</p>
Full article ">Figure 8
<p>Influencing factor maps selected using principal component analysis (PCA): (<b>a</b>) Factor 1; (2) Factor 2; (3) Factor 3; and (<b>b</b>) Factor 4.</p>
Full article ">Figure 9
<p>Influencing factors maps selected by PCA: (<b>a</b>) Factor 5 and (<b>b</b>) Factor 6.</p>
Full article ">Figure 10
<p>Landslide susceptibility map: (<b>a</b>) PCA-logistic regression (LR) method; (<b>b</b>) FR method; and (<b>c</b>) analytic hierarchy process (AHP) method.</p>
Full article ">Figure 11
<p>Receiver operating characteristic (ROC) curve of the model.</p>
Full article ">
21 pages, 2296 KiB  
Article
The Elephant in the Room: Informality in Tanzania’s Rural Waterscape
by Jesper Katomero and Yola Georgiadou
ISPRS Int. J. Geo-Inf. 2018, 7(11), 437; https://doi.org/10.3390/ijgi7110437 - 8 Nov 2018
Cited by 4 | Viewed by 6651
Abstract
Informality is pervasive in Tanzania’s rural waterscape, but not acknowledged by development partners (donors and beneficiaries), despite persistent warnings by development scholars. Informality is thus the proverbial elephant in the room. In this paper, we examine a case of superior rural water access [...] Read more.
Informality is pervasive in Tanzania’s rural waterscape, but not acknowledged by development partners (donors and beneficiaries), despite persistent warnings by development scholars. Informality is thus the proverbial elephant in the room. In this paper, we examine a case of superior rural water access in two geographical locales—Hai and Siha districts—in Tanzania, where actors not only acknowledge, but actively harness informality to provide access to water to rural populations. We employ concepts from organization and institutional theory to show that when informal programs and related informal sanctions/rewards complement their formal counterparts, chances for achieving the Sustainable Development Goals (SDG) target 6.1 ‘By 2030, achieve universal and equitable access to safe and affordable drinking water for all’ are significantly increased. Full article
(This article belongs to the Special Issue Geo-Information and the Sustainable Development Goals (SDGs))
Show Figures

Figure 1

Figure 1
<p>An “either-or” view of (in)formality and development, reproduced from [<a href="#B14-ijgi-07-00437" class="html-bibr">14</a>].</p>
Full article ">Figure 2
<p>Organisational levels in the rural water supply in Hai and Siha (Source: Fieldwork data, 2017).</p>
Full article ">Figure 3
<p>Part of O&amp;M infrastructure in Hai &amp; Siha (Source: Field work data).</p>
Full article ">Figure 4
<p>Appraisal form used to evaluate the WTs [<a href="#B37-ijgi-07-00437" class="html-bibr">37</a>].</p>
Full article ">Figure 5
<p>The Chief engineer of Hai and Siha pictured left, in the middle and right engaging Islamic and evangelical church leaders on rural water supply issues [Fieldwork data, 2017] [<a href="#B37-ijgi-07-00437" class="html-bibr">37</a>].</p>
Full article ">
19 pages, 4986 KiB  
Article
Crisis Maps—Observed Shortcomings and Recommendations for Improvement
by Ana Kuveždić Divjak and Miljenko Lapaine
ISPRS Int. J. Geo-Inf. 2018, 7(11), 436; https://doi.org/10.3390/ijgi7110436 - 7 Nov 2018
Cited by 15 | Viewed by 4531
Abstract
Cartographic communication through crisis maps takes place in a unique environment characterised by the immediate risks of considerable loss and stress. Many such maps are designed by practitioners with limited resources, pressured for time, and who often fail to pay the necessary attention [...] Read more.
Cartographic communication through crisis maps takes place in a unique environment characterised by the immediate risks of considerable loss and stress. Many such maps are designed by practitioners with limited resources, pressured for time, and who often fail to pay the necessary attention to map graphics. This can reduce map clarity and make orientation to and understanding of essential crisis information difficult. To identify the most frequent shortcomings that may compromise the interpretation of depicted objects, phenomena presented, and actions required, we assessed the map graphics of 106 maps specifically designed for communication and action in crises. The results showed that they were often visually overloaded. Crisis data were not conveyed by appropriate cartographic representations, and due to the inappropriate use of visual variables, the associative and selective properties of cartographic symbols were overlooked, and their ordered and quantitative features ignored. The use of colour was often not adapted to conventional visual language, and colour symbolism was not always taken into account. The cartographic symbols used were often incomprehensible, illegible, ambiguous, and unclassified, and they lacked symbolism and hierarchical organisation. The article aims to address these problems by proposing guidelines which do not require much time or expertise, but which would ensure that cartographically correct crisis maps are well designed. Objects, phenomena or actions specific to crisis management would be indicated using appropriate map graphics and their importance highlighted, so as to make interpretation easier for all participants in a crisis event, and so facilitate crisis communication and response. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Maps grouped by crisis type. The percentage of maps analysed for each group is given above the column.</p>
Full article ">Figure 2
<p>Maps divided according to scale into five groups (very large-scale maps, large-scale maps, medium-scale maps, small-scale maps and very small-scale maps).</p>
Full article ">Figure 3
<p>The relationship between crisis type and scale on the maps collected. Crisis communication maps are designed on a large or very large scale, regardless of the type of the crisis depicted.</p>
Full article ">Figure 4
<p>Examples of different background maps for the representation of cartographic symbols in crisis communication: (1) aerial photograph (Map No. 19 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (2) topographic map (Map No. 94 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (3) relief map (Map No. 98 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (4) city map (Map No. 78 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>) and (5) informative map (Map No. 14 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>).</p>
Full article ">Figure 5
<p>When optimal contrast between the background map and thematic content about a crisis event is achieved, optimal layering will also be achieved. (1) Too little light-dark and colour contrast. (2) Sufficient contrast between the background map and thematic content about the crisis event. Grey on a background map forms a neutral ground for other colours.</p>
Full article ">Figure 6
<p>Geometric and graphic elements analysed divided into three groups (<span class="html-italic">points</span>, <span class="html-italic">lines</span>, <span class="html-italic">areas</span>). The percentage of elements used for each group is given above the column.</p>
Full article ">Figure 7
<p>Examples of area symbols used to represent: (1) a wildfire-burned area (Map No. 42 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (2) a flooded area (Map No. 10 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), and (3) an area affected by a hazardous gas leak (Map No. 31 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>).</p>
Full article ">Figure 8
<p>Examples of line symbols used to represent: (1) fire front progression (Map No. 4 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (2) oil pollution (Map No. 85 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), and (3) evacuation routes during a crisis (Map No. 12 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>).</p>
Full article ">Figure 9
<p>Examples of point symbols used to represent: (1) objects and actions during a fire (Map No. 45 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (2) calls for emergency assistance after a hurricane (Map No. 84 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), and (3) the number of refugees due to a humanitarian crisis (Map No. 3 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>).</p>
Full article ">Figure 10
<p>Good examples of map graphics on the crisis maps analysed (1) (Map No. 35 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (2) (Map No. 1 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (3) (Map No. 102 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), (4) (Map No. 103 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>), and (5) (Map No. 38 in <a href="#app1-ijgi-07-00436" class="html-app">Figure S1</a>). Correctly used map graphics contribute not only to clarity, legibility and layout, but also facilitate crisis communication and response.</p>
Full article ">
12 pages, 8017 KiB  
Article
Spatial Distribution Estimates of the Urban Population Using DSM and DEM Data in China
by Junlin Zhang, Wei Xu, Lianjie Qin and Yugang Tian
ISPRS Int. J. Geo-Inf. 2018, 7(11), 435; https://doi.org/10.3390/ijgi7110435 - 7 Nov 2018
Cited by 10 | Viewed by 4293
Abstract
Spatial distribution and population density are important parameters in studies on urban development, resource allocation, emergency management, and risk analysis. High-resolution height data can be used to estimate the total or spatial pattern of the urban population for small study areas, e.g., the [...] Read more.
Spatial distribution and population density are important parameters in studies on urban development, resource allocation, emergency management, and risk analysis. High-resolution height data can be used to estimate the total or spatial pattern of the urban population for small study areas, e.g., the downtown area of a city or a community. However, there has been no case of population estimation for large areas. This paper tries to estimate the urban population of prefectural cities in China using building height data. Building height in urban population settlement (Mdiffs) was first extracted using the digital surface model (DSM), digital elevation model (DEM), and land use data. Then, the relationships between the census-based urban population density (CPD) and the Mdiffs density (MDD) for different regions were regressed. Using these results, the urban population for prefectural cities of China was finally estimated. The results showed that a good linear correlation was found between Mdiffs and the census data in each type of region, as all the adjusted R2 values were above 0.9 and all the models passed the significance test (95% confidence level). The ratio of the estimated population to the census population (PER) was between 0.7 and 1.3 for 76% of the cities in China. This is the first attempt to estimate the urban population using building height data for prefectural cities in China. This method produced reasonable results and can be effectively used for spatial distribution estimates of the urban population in large scale areas. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Technical flow of the proposed method.</p>
Full article ">Figure 2
<p>Distribution of <span class="html-italic">k.</span></p>
Full article ">Figure 3
<p>Distribution of the urban estimation and census population at the prefectural city level. (<b>a</b>) Estimated urban population size. (<b>b</b>) Census urban population size. (<b>c</b>) Estimated urban population density. (<b>d</b>) Census urban population density.</p>
Full article ">Figure 4
<p>Distribution of <span class="html-italic">PER</span> (EP/CP).</p>
Full article ">Figure 5
<p>The histogram of <span class="html-italic">PER</span> (EP/CP).</p>
Full article ">
18 pages, 2458 KiB  
Article
An Automatic Recognition and Positioning Method for Point Source Targets on Satellite Images
by Kai Li, Yongsheng Zhang, Zhenchao Zhang and Ying Yu
ISPRS Int. J. Geo-Inf. 2018, 7(11), 434; https://doi.org/10.3390/ijgi7110434 - 7 Nov 2018
Cited by 3 | Viewed by 3138
Abstract
Currently, the geometric and radiometric calibration of on-board satellite sensors utilizes different ground targets using some form of manual intervention. Point source targets provide high precision geometric and radiometric information and have the potential to become a new tool for joint geometric and [...] Read more.
Currently, the geometric and radiometric calibration of on-board satellite sensors utilizes different ground targets using some form of manual intervention. Point source targets provide high precision geometric and radiometric information and have the potential to become a new tool for joint geometric and radiometric calibration. In this paper, an automatic recognition and positioning method for point source target images is proposed. First, the template matching method was used to effectively reduce nonpoint source target image pixels in the satellite imagery. The point source target images were then identified using particular feature parameters. Using the template matching method, the weighted centroid method, and the Gaussian fitting method, the positions of the centroid of the point source target images were calculated. The maximum position detection error obtained using the three methods was 0.07 pixels, which is comparably better than artificial targets currently in use. The experimental results show point source targets provide high precision geometric information, which can become a suitable alternative for automatic joint geometric and radiometric calibration of spaceborne optical sensors. Full article
Show Figures

Figure 1

Figure 1
<p>Reflective point source.</p>
Full article ">Figure 2
<p>Degraded images of point sources when phases of point sources are (0, 0) (<b>a</b>) and (0.6, 0.8) (<b>b</b>).</p>
Full article ">Figure 3
<p>Different distribution of pixel grayscale caused by sampling phase.</p>
Full article ">Figure 4
<p>Experimental data: (<b>a</b>) Satellite imagery; (<b>b</b>) point source target image; (<b>c</b>) large area grayscale target image; and (<b>d</b>) image patch containing a point source target.</p>
Full article ">Figure 5
<p>(<b>a</b>) Edge Spread Function (ESF) raw sample values, average values within each window, and smoothed curves in the horizontal direction and (<b>b</b>) ESF raw sample values, average values within each window, and smoothed curves in the vertical direction.</p>
Full article ">Figure 6
<p>(<b>a</b>) Edge Line Spread Function (LSF) curve and fitted Gaussian curve in the horizontal direction and (<b>b</b>) Edge LSF curve and fitted Gaussian curve in the vertical direction.</p>
Full article ">Figure 7
<p>Number of point source image candidates retained in the 100 image patches.</p>
Full article ">Figure 8
<p>(<b>a</b>) Point source target image, (<b>b</b>) nonpoint source target image, and (<b>c</b>) common matching template of (<b>a</b>,<b>b</b>).</p>
Full article ">Figure 9
<p>Curve of feature parameters as a function of candidate point source image number. (<b>a</b>–<b>d</b>) correspond to <math display="inline"><semantics> <mover accent="true"> <mi>σ</mi> <mo stretchy="false">^</mo> </mover> </semantics></math> and <math display="inline"><semantics> <mover accent="true"> <mi>ξ</mi> <mo stretchy="false">^</mo> </mover> </semantics></math>, <math display="inline"><semantics> <mover accent="true"> <mi>b</mi> <mo stretchy="false">^</mo> </mover> </semantics></math>, <math display="inline"><semantics> <mrow> <mrow> <mrow> <mo stretchy="false">(</mo> <mover accent="true"> <mi>K</mi> <mo stretchy="false">^</mo> </mover> <mo>+</mo> <mover accent="true"> <mi>b</mi> <mo stretchy="false">^</mo> </mover> <mo stretchy="false">)</mo> </mrow> <mo>/</mo> <mover accent="true"> <mi>b</mi> <mo stretchy="false">^</mo> </mover> </mrow> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>ε</mi> <mrow> <mi>min</mi> </mrow> </msub> </mrow> </semantics></math>, respectively.</p>
Full article ">
35 pages, 6541 KiB  
Article
Impaired Water Hazard Zones: Mapping Intersecting Environmental Health Vulnerabilities and Polluter Disproportionality
by Raoul S. Liévanos
ISPRS Int. J. Geo-Inf. 2018, 7(11), 433; https://doi.org/10.3390/ijgi7110433 - 6 Nov 2018
Cited by 10 | Viewed by 7208
Abstract
This study advanced a rigorous spatial analysis of surface water-related environmental health vulnerabilities in the California Bay-Delta region, USA, from 2000 to 2006. It constructed a novel hazard indicator—“impaired water hazard zones’’—from regulatory estimates of extensive non-point-source (NPS) and point-source surface water pollution, [...] Read more.
This study advanced a rigorous spatial analysis of surface water-related environmental health vulnerabilities in the California Bay-Delta region, USA, from 2000 to 2006. It constructed a novel hazard indicator—“impaired water hazard zones’’—from regulatory estimates of extensive non-point-source (NPS) and point-source surface water pollution, per section 303(d) of the U.S. Clean Water Act. Bivariate and global logistic regression (GLR) analyses examined how established predictors of surface water health-hazard exposure vulnerability explain census block groups’ proximity to impaired water hazard zones in the Bay-Delta. GLR results indicate the spatial concentration of Black disadvantage, isolated Latinx disadvantage, low median housing values, proximate industrial water pollution levels, and proximity to the Chevron oil refinery—a disproportionate, “super emitter”, in the Bay-Delta—significantly predicted block group proximity to impaired water hazard zones. A geographically weighted logistic regression (GWLR) specification improved model fit and uncovered spatial heterogeneity in the predictors of block group proximity to impaired water hazard zones. The modal GWLR results in Oakland, California, show how major polluters beyond the Chevron refinery impair the local environment, and how isolated Latinx disadvantage was the lone positively significant population vulnerability factor. The article concludes with a discussion of its scholarly and practical implications. Full article
(This article belongs to the Special Issue Geoprocessing in Public and Environmental Health)
Show Figures

Figure 1

Figure 1
<p>Map (<b>A</b>) California and Map (<b>B</b>) the six-county Bay-Delta study area.</p>
Full article ">Figure 2
<p>The critical physical geography “methodological four-square” [<a href="#B51-ijgi-07-00433" class="html-bibr">51</a>] adapted to the present study’s examination of impaired water hazard zones and the sociospatial dimensions of water injustice in the Bay-Delta. CAWMHS, cumulative areal-weighted modeled hazard score.</p>
Full article ">Figure 3
<p>Illustration of Map (<b>A</b>) the impaired minor water bodies and their buffers, Map (<b>B</b>) the impaired major water bodies and their buffers, and Map (<b>C</b>) the impaired water hazard zones that result from intersecting the buffers for impaired minor and major water bodies in the Bay-Delta, 2006.</p>
Full article ">Figure 4
<p>The spatial distribution of Black disadvantage in the Bay-Delta, 2000 (N = 3064).</p>
Full article ">Figure 5
<p>The spatial distribution of isolated Latinx disadvantage in the Bay-Delta, 2000 (N = 3061).</p>
Full article ">Figure 6
<p>The spatial distribution of relative median housing values in the Bay-Delta, 2000 (N = 3041).</p>
Full article ">Figure 7
<p>The spatial distribution of percentile ranges for surface water CAWMHS associated with surface water toxic releases from 2000 to 2006 in Map (<b>A</b>) the Bay-Delta region (N = 3073) and Map (<b>B</b>) in the Richmond-Berkeley area.</p>
Full article ">Figure 8
<p>GWLR model fit statistic: Percent local deviance explained by block group.</p>
Full article ">Figure 9
<p>Residuals from Map (<b>A</b>) GLR analysis and Map (<b>B</b>) GWLR analysis.</p>
Full article ">Figure 10
<p>Results from spatial autocorrelation analysis of residuals from the GLR analysis and the GWLR analysis by nearest neighbor spatial weights matrix. <span class="html-italic">Note</span>: Square markers on trend lines indicate significant Moran’s I value (<span class="html-italic">p</span> &lt; 0.001; 9999 permutations).</p>
Full article ">Figure 11
<p>GWLR <span class="html-italic">t</span>-values and odds ratios for Map (<b>A</b>) Black disadvantage and Map (<b>B</b>) isolated Latinx disadvantage.</p>
Full article ">Figure 12
<p>GWLR <span class="html-italic">t</span>-values and odds ratios for Map (<b>A</b>) surface water CAWMHS, 2000−2006 (100,000s) and Map (<b>B</b>) kilometers to the Richmond Chevron refinery.</p>
Full article ">Figure 13
<p>Modal GWLR results in the Map (<b>A</b>) Bay-Delta and Map (<b>B</b>) the Oakland-Alameda-Piedmont Area.</p>
Full article ">
26 pages, 12449 KiB  
Article
Winter Is Coming: A Socio-Environmental Monitoring and Spatiotemporal Modelling Approach for Better Understanding a Respiratory Disease
by Lukas Marek, Malcolm Campbell, Michael Epton, Simon Kingham and Malina Storer
ISPRS Int. J. Geo-Inf. 2018, 7(11), 432; https://doi.org/10.3390/ijgi7110432 - 6 Nov 2018
Cited by 6 | Viewed by 4566
Abstract
Chronic Obstructive Pulmonary Disease is a progressive lung disease affecting the respiratory function of every sixth New Zealander and over 300 million people worldwide. In this paper, we explored how the combination of social, demographical and environmental conditions (represented by increased winter air [...] Read more.
Chronic Obstructive Pulmonary Disease is a progressive lung disease affecting the respiratory function of every sixth New Zealander and over 300 million people worldwide. In this paper, we explored how the combination of social, demographical and environmental conditions (represented by increased winter air pollution) affected hospital admissions due to COPD in an urban area of Christchurch (NZ). We juxtaposed the hospitalisation data with dynamic air pollution data and census data to investigate the spatiotemporal patterns of hospital admissions. Spatial analysis identified high-risk health hot spots both overall and season specific, exhibiting higher rates in winter months not solely due to air pollution, but rather as a result of its combination with other factors that initiate deterioration of breathing, increasing impairments and lead to the hospitalisation of COPD patients. From this we found that socioeconomic deprivation and air pollution, followed by the age and ethnicity structure contribute the most to the increased winter hospital admissions. This research shows the continued importance of including both individual (composition) and area level (composition) factors when examining and analysing disease patterns. Full article
(This article belongs to the Special Issue Geoprocessing in Public and Environmental Health)
Show Figures

Figure 1

Figure 1
<p>Age and gender structure of people hospitalised due to COPD in Christchurch, August 2014–October 2017 (<b>a</b>), and absolute frequency of hospital admissions due to COPD (<b>b</b>).</p>
Full article ">Figure 2
<p>Comparison of average PM<sub>2.5</sub> concentrations in Christchurch—annually (<b>a</b>), in winter months (<b>b</b>), and during non-winter months (<b>c</b>).</p>
Full article ">Figure 3
<p>Incidence of hospitalisations due to COPD exacerbations annually (<b>a</b>), in winter (<b>b</b>) and non-winter (<b>c</b>) months.</p>
Full article ">Figure 4
<p>SIR of hospitalisation due to Chronic Obstructive Pulmonary Disease (COPD) exacerbations annually (<b>a</b>), in winter (<b>b</b>) and non-winter (<b>c</b>) months.</p>
Full article ">Figure 5
<p>Local <span class="html-italic">G<sub>i</sub>*</span> of average and season-specific standardised incidence ratios (SIRs) of COPD hospitalisation (clusters of high SIRs are in red, clusters of low SIRs are in blue)—annually (<b>a</b>), in winter (<b>b</b>) and non-winter (<b>c</b>) months.</p>
Full article ">Figure 6
<p>Spatiotemporal clusters of SIRs of hospitalisations due to COPD in Census Area Units (CAU).</p>
Full article ">Figure 7
<p>Spearman’s correlation among season specific SIRs of hospitalisations due to COPD and likely triggering factors. Statistically significant correlations (<span class="html-italic">p</span>-value &lt; 0.05) are emphasised by the square coloured and sized accordingly to the strength of correlation.</p>
Full article ">Figure 8
<p>Local correlation among average and season-specific SIRs and characteristics of CAUs.</p>
Full article ">Figure 9
<p>Boxplot of characteristics of CAU with lines characterising clusters. (Group colours match the group colours in <a href="#ijgi-07-00432-f010" class="html-fig">Figure 10</a>).</p>
Full article ">Figure 10
<p>The location of clusters of CAUs in the study area. (Group colours match the group colours in <a href="#ijgi-07-00432-f009" class="html-fig">Figure 9</a> and <a href="#ijgi-07-00432-t0A1" class="html-table">Table A1</a>).</p>
Full article ">Figure 11
<p>The estimates of GWPR regression parameters (log-odds) in CAUs (classified in sextiles during) the winter season. The coefficients estimates are log-odds. Positive values of the log-odds indicate positive relationships between the explanatory variable and COPD hospitalisations, while negative values of the log-odds indicate negative relationships. Dots indicate significant associations (<span class="html-italic">p</span> &lt; 0.05) based on pseudo <span class="html-italic">t</span>-value.</p>
Full article ">Figure 12
<p>The estimates of GWPR regression parameters (log-odds) in CAUs (classified in sextiles) during the non-winter season. The coefficients estimates are log-odds. Positive values of the log-odds indicate positive relationships between the explanatory variable and COPD hospitalisations, while negative values of the log-odds indicate negative relationships. Dots indicate significant associations (<span class="html-italic">p</span> &lt; 0.05) based on pseudo <span class="html-italic">t</span>-value.</p>
Full article ">Figure A1
<p>Census area units and location of air pollution monitoring sites in Christchurch.</p>
Full article ">Figure A2
<p>Geographic distribution of COPD hospital admissions in Christchurch—annually (<b>a</b>), in winter (<b>b</b>) and non-winter (<b>c</b>) months.</p>
Full article ">
21 pages, 8188 KiB  
Article
Intact Planar Abstraction of Buildings via Global Normal Refinement from Noisy Oblique Photogrammetric Point Clouds
by Qing Zhu, Feng Wang, Han Hu, Yulin Ding, Jiali Xie, Weixi Wang and Ruofei Zhong
ISPRS Int. J. Geo-Inf. 2018, 7(11), 431; https://doi.org/10.3390/ijgi7110431 - 6 Nov 2018
Cited by 6 | Viewed by 3666
Abstract
Oblique photogrammetric point clouds are currently one of the major data sources for the three-dimensional level-of-detail reconstruction of buildings. However, they are severely noise-laden and pose serious problems for the effective and automatic surface extraction of buildings. In addition, conventional methods generally use [...] Read more.
Oblique photogrammetric point clouds are currently one of the major data sources for the three-dimensional level-of-detail reconstruction of buildings. However, they are severely noise-laden and pose serious problems for the effective and automatic surface extraction of buildings. In addition, conventional methods generally use normal vectors estimated in a local neighborhood, which are liable to be affected by noise, leading to inferior results in successive building reconstruction. In this paper, we propose an intact planar abstraction method for buildings, which explicitly handles noise by integrating information in a larger context through global optimization. The information propagates hierarchically from a local to global scale through the following steps: first, based on voxel cloud connectivity segmentation, single points are clustered into supervoxels that are enforced to not cross the surface boundary; second, each supervoxel is expanded to nearby supervoxels through the maximal support region, which strictly enforces planarity; third, the relationships established by the maximal support regions are injected into a global optimization, which reorients the local normal vectors to be more consistent in a larger context; finally, the intact planar surfaces are obtained by region growing using robust normal and point connectivity in the established spatial relations. Experiments on the photogrammetric point clouds obtained from oblique images showed that the proposed method is effective in reducing the influence of noise and retrieving almost all of the major planar structures of the examined buildings. Full article
Show Figures

Figure 1

Figure 1
<p>Different region selection and geometric kernel for normal estimation. Principal component analysis (PCA) (<b>a</b>), jet fitting (<b>b</b>) and MLS Sphere fitting (<b>c</b>) generally use local regions, and select planes, quadratic surfaces, and spheres as the geometric kernel, respectively. In addition, the geometric kernel can also be learned from exemplar datasets (<b>d</b>). However, for noisy photogrammetric point clouds in urban environments, a large context should be selected to keep sharp features, such as the sharp edges of a building.</p>
Full article ">Figure 2
<p>Processing pipeline of the proposed method.</p>
Full article ">Figure 3
<p>Directional expansion of the maximal support regions. (<b>a</b>) Photogrammetric point cloud, (<b>b</b>) supervoxels, (<b>c</b>) growing region of a supervoxel located in a plane, and (<b>d</b>) growing region of a supervoxel located in a sharp feature.</p>
Full article ">Figure 4
<p>Constraints used during region growing of the support region. (<b>a</b>) Planarity constraint determined from the spectral features. (<b>b</b>) Angle deviation in the supervoxels.</p>
Full article ">Figure 5
<p>Mutual validation of the connectivity test. The dots represent supervoxels and the rounded rectangles represent support regions. (<b>a</b>) Only one side of the region is contained, <span class="html-italic">c<sub>j</sub></span> ∈ <span class="html-italic">S<sub>i</sub></span> but <span class="html-italic">c<sub>i</sub></span> <math display="inline"><semantics> <mo>∉</mo> </semantics></math> <span class="html-italic">S<sub>j</sub></span>. (<b>b</b>) The shaded dots are mutually contained, <span class="html-italic">c<sub>i</sub></span> ∈ <span class="html-italic">S<sub>j</sub></span> and <span class="html-italic">c<sub>j</sub></span> ∈ <span class="html-italic">S<sub>i</sub></span>.</p>
Full article ">Figure 6
<p>Visualization of normal vectors. (<b>a</b>) shows a cube model with 240,000 points. (<b>c</b>) shows the same model but with 0.5% noise. (<b>b</b>,<b>d</b>) show the results from the proposed method, the RMSE are 0.0022 and 0.0231 respectively. (<b>c</b>,<b>f</b>) show the results estimated by local information only, the RMSE are 0.0707 and 0.433 respectively. The proposed method produces more consistent normal vectors.</p>
Full article ">Figure 7
<p>Normal vectors visualization. The left column represents photogrammetric point clouds, the middle column represents normal vectors optimized by IPA, and the right column represents initial normal vectors estimated by the PCA.</p>
Full article ">Figure 8
<p>Comparison of planar surface extraction from the point cloud of tile 1.</p>
Full article ">Figure 9
<p>Comparison of planar surface extraction from the point cloud of tile 2.</p>
Full article ">Figure 10
<p>Comparison of planar surface extraction from the point cloud of tile 3.</p>
Full article ">Figure 11
<p>Comparison of the planar abstraction quality of building 1.</p>
Full article ">Figure 12
<p>Comparison of the planar abstraction quality of building 2. (<b>a</b>) Building point cloud, (<b>b</b>) planar abstractions, (<b>c</b>) extracted planes approximated by a set of planar polygons, and (<b>d</b>) magnified images of the circle marked regions in (<b>c</b>).</p>
Full article ">Figure 13
<p>Comparison of the planar abstraction quality of building 3. (<b>a</b>) Building point cloud, (<b>b</b>) planar abstractions, and (<b>c</b>) extracted planes approximated by a set of planar polygons.</p>
Full article ">Figure 14
<p>Efficiency analysis of algorithms.</p>
Full article ">
13 pages, 9800 KiB  
Article
Use of a Multilayer Perceptron to Automate Terrain Assessment for the Needs of the Armed Forces
by Krzysztof Pokonieczny
ISPRS Int. J. Geo-Inf. 2018, 7(11), 430; https://doi.org/10.3390/ijgi7110430 - 6 Nov 2018
Cited by 9 | Viewed by 4445
Abstract
The classification of terrain in terms of passability plays a significant role in the process of military terrain assessment. It involves classifying selected terrain to specific classes (GO, SLOW-GO, NO-GO). In this article, the problem of terrain classification to the respective category of [...] Read more.
The classification of terrain in terms of passability plays a significant role in the process of military terrain assessment. It involves classifying selected terrain to specific classes (GO, SLOW-GO, NO-GO). In this article, the problem of terrain classification to the respective category of passability was solved by applying artificial neural networks (multilayer perceptron) to generate a continuous Index of Passability (IOP). The neural networks defined this factor for primary fields in two sizes (1000 × 1000 m and 100 × 100 m) based on the land cover elements obtained from Vector Smart Map (VMap) Level 2 and Shuttle Radar Topography Mission (SRTM). The work used a feedforward neural network consisting of three layers. The paper presents a comprehensive analysis of the reliability of the neural network parameters, taking into account the number of neurons, learning algorithm, activation functions and input data configuration. The studies and tests carried out have shown that a well-trained neural network can automate the process of terrain classification in terms of passability conditions. Full article
Show Figures

Figure 1

Figure 1
<p>The study area.</p>
Full article ">Figure 2
<p>Thematic categories and feature classes in VMap Level 2 (A: Area, L: Line, P: Point object) [<a href="#B23-ijgi-07-00430" class="html-bibr">23</a>].</p>
Full article ">Figure 3
<p>Example of gathering land cover element data for a 1000 × 1000 m grid.</p>
Full article ">Figure 4
<p>Example of a vector model conversion to a grid for the built-up area feature class.</p>
Full article ">Figure 5
<p>McCulloch and Pitts artificial neuron model [<a href="#B27-ijgi-07-00430" class="html-bibr">27</a>]<sup>.</sup></p>
Full article ">Figure 6
<p>Structure of used neural network.</p>
Full article ">Figure 7
<p>Configuration of parameters taken into account when generating artificial neural networks.</p>
Full article ">Figure 8
<p>Average quality of the validation sample based on the number of neurons in the hidden layer for two sizes of primary field.</p>
Full article ">Figure 9
<p>Average quality of the validation sample based on the number of learning iterations for two sizes of primary field.</p>
Full article ">Figure 10
<p>Part of a passability map generated using a primary field of 1000 × 1000 m with avenues of approach and the most significant elements of land cover.</p>
Full article ">Figure 11
<p>Part of a passability map generated using a primary field of 100 × 100 m with avenues of approach and the most significant elements of land cover.</p>
Full article ">Figure 12
<p>Comparison of map generation time.</p>
Full article ">
19 pages, 3115 KiB  
Article
Importance of Remotely-Sensed Vegetation Variables for Predicting the Spatial Distribution of African Citrus Triozid (Trioza erytreae) in Kenya
by Kyalo Richard, Elfatih M. Abdel-Rahman, Samira A. Mohamed, Sunday Ekesi, Christian Borgemeister and Tobias Landmann
ISPRS Int. J. Geo-Inf. 2018, 7(11), 429; https://doi.org/10.3390/ijgi7110429 - 3 Nov 2018
Cited by 27 | Viewed by 5713
Abstract
Citrus is considered one of the most important fruit crops globally due to its contribution to food and nutritional security. However, the production of citrus has recently been in decline due to many biological, environmental, and socio-economic constraints. Amongst the biological ones, pests [...] Read more.
Citrus is considered one of the most important fruit crops globally due to its contribution to food and nutritional security. However, the production of citrus has recently been in decline due to many biological, environmental, and socio-economic constraints. Amongst the biological ones, pests and diseases play a major role in threatening citrus quantity and quality. The most damaging disease in Kenya, is the African citrus greening disease (ACGD) or Huanglongbing (HLB) which is transmitted by the African citrus triozid (ACT), Trioza erytreae. HLB in Kenya is reported to have had the greatest impact on citrus production in the highlands, causing yield losses of 25% to 100%. This study aimed at predicting the occurrence of ACT using an ecological habitat suitability modeling approach. Specifically, we tested the contribution of vegetation phenological variables derived from remotely-sensed (RS) data combined with bio-climatic and topographical variables (BCL) to accurately predict the distribution of ACT in citrus-growing areas in Kenya. A MaxEnt (maximum entropy) suitability modeling approach was used on ACT presence-only data. Forty-seven (47) ACT observations were collected while 23 BCL and 12 RS covariates were used as predictor variables in the MaxEnt modeling. The BCL variables were extracted from the WorldClim data set, while the RS variables were predicted from vegetation phenological time-series data (spanning the years 2014–2016) and annually-summed land surface temperature (LST) metrics (2014–2016). We developed two MaxEnt models; one including both the BCL and the RS variables (BCL-RS) and another with only the BCL variables. Further, we tested the relationship between ACT habitat suitability and the surrounding land use/land cover (LULC) proportions using a random forest regression model. The results showed that the combined BCL-RS model predicted the distribution and habitat suitability for ACT better than the BCL-only model. The overall accuracy for the BCL-RS model result was 92% (true skills statistic: TSS = 0.83), whereas the BCL-only model had an accuracy of 85% (TSS = 0.57). Also, the results revealed that the proportion of shrub cover surrounding citrus orchards positively influenced the suitability probability of the ACT. These results provide a resourceful tool for precise, timely, and site-specific implementation of ACGD control strategies. Full article
Show Figures

Figure 1

Figure 1
<p>Study area (major citrus-growing regions) in Kenya where the African citrus triozid (<span class="html-italic">Trioza erytreae</span>) presence data were collected.</p>
Full article ">Figure 2
<p>Collinearity matrix for predictor variables. Darker shades of blue and red colors indicate high variable collinearity while light shades indicate low collinearity between variables.</p>
Full article ">Figure 3
<p>A 20-m spatial resolution land use/land cover map for the study area generated by the Climate Change Initiative (CCI) Land Cover (LC) team. Using yearly Sentinel-2 observations. (<b>a</b>–<b>c</b>) represent zoomed buffers of 1500 m radius each around certain representative African citrus triozid (ACT) occurrence points.</p>
Full article ">Figure 4
<p>Jackknife variable importance test of regulated gains for the BCL model. The dark blue shades show the regularized training gain for the specific variable, light blue shows the relevance when the variable is omitted, while red shows the regularized training gain with all the variables combined.</p>
Full article ">Figure 5
<p>Jackknife variable importance test of regulated gains for the BCL-RS model. The dark blue shades show the regularized training gain for the specific variable, light blue illustrates gains without the variable, while red shows the regularized training gain with all the variables combined.</p>
Full article ">Figure 6
<p>Predicted distribution suitability map for African citrus triozid (<span class="html-italic">Trioza erytreae</span>) using environmental (BCL model) variables (<b>a</b>), and environmental and remotely-sensed (BCL-RS model) variables (<b>b</b>). Blue indicates low distribution suitability, while red represents high distribution suitability.</p>
Full article ">Figure 7
<p>The relevance of the four major land use/land cover classes to the habitat suitability of African citrus triozid (<span class="html-italic">Trioza erytreae</span>) using a random forest variable importance rank.</p>
Full article ">
28 pages, 15767 KiB  
Review
The Scope of Earth-Observation to Improve the Consistency of the SDG Slum Indicator
by Monika Kuffer, Jiong Wang, Michael Nagenborg, Karin Pfeffer, Divyani Kohli, Richard Sliuzas and Claudio Persello
ISPRS Int. J. Geo-Inf. 2018, 7(11), 428; https://doi.org/10.3390/ijgi7110428 - 1 Nov 2018
Cited by 65 | Viewed by 10397
Abstract
The continuous increase in deprived living conditions in many cities of the Global South contradicts efforts to make cities inclusive, safe, resilient, and sustainable places. Using examples of Asian, African, and Latin American cities, this study shows the scope and limits of earth [...] Read more.
The continuous increase in deprived living conditions in many cities of the Global South contradicts efforts to make cities inclusive, safe, resilient, and sustainable places. Using examples of Asian, African, and Latin American cities, this study shows the scope and limits of earth observation (EO)-based mapping of deprived living conditions in support of providing consistent global information for the SDG indicator 11.1.1 “proportion of urban population living in slums, informal settlements or inadequate housing”. At the technical level, we compare several EO-based methods and imagery for mapping deprived living conditions, discussing their ability to map such areas including differences in terms of accuracy and performance at the city scale. At the operational level, we compare available municipal maps showing identified deprived areas with the spatial extent of morphological mapped areas of deprived living conditions (using EO) at the city scale, discussing the reasons for inconsistencies between municipal and EO-based maps. We provide an outlook on how EO-based mapping of deprived living conditions could contribute to a global spatial information base to support targeting of deprived living conditions in support of the SDG Goal 11.1.1 indicator, when uncertainties and ethical considerations on data provision are well addressed. Full article
(This article belongs to the Special Issue Geo-Information and the Sustainable Development Goals (SDGs))
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Romani settlements in Europe: Novi Sad, Serbia (<b>left</b>) and deprived settlement in Paris (yellow outlines) using the abandoned Petite Ceinture railway (<b>right</b>) (Image source: Google Earth).</p>
Full article ">Figure 2
<p>Overview of methodological steps and aspects to utilize remote sensing for the SDG indicator 11.1.1.</p>
Full article ">Figure 3
<p>The ability of different sensor systems to capture a small slum pocket in Bangalore, India (image dimension: 400 × 200 m), with slum pocket in yellow (area: 2832 m<sup>2</sup>) (Image sources: DigitalGlobe).</p>
Full article ">Figure 4
<p>Comparing VHR satellite imagery with UAV images for the case of Dar es Salaam: (<b>1</b>) and (<b>3</b>) show an informal settlement surrounded by formal areas ((<b>1</b>): VHR and (<b>3</b>): UAV image); (<b>2</b>) shows the level of detail available in VHR imagery; (<b>4</b>,<b>5</b>) show the level of detail available in UAV images (Image source: Google Earth and Ramani Huria).</p>
Full article ">Figure 5
<p>Citywide-level mapping of deprived areas in Bangalore, India: (<b>a</b>) 8 km × 8 km coverage of downtown Bangalore with numbered sample areas for illustration; (<b>b</b>) reference data with sample areas highlighted by boxes; (<b>c</b>) predicted deprived areas with sample areas highlighted by boxes; (<b>d</b>) zoomed-in sample area containing the training samples; and (<b>e</b>–<b>h</b>) zoomed-in sample areas of predictions with reference data superimposed (Image source: WorldView-2 images provided by DigitalGlobe).</p>
Full article ">Figure 6
<p>Deprived area mapping in Dar es Salaam, Tanzania: (<b>a</b>) reference data of deprived areas showing the center and its adjacent western suburban areas; (<b>b</b>) deprived area mapping using training area 1 labelled in (<b>a</b>); (<b>c</b>) deprived area mapping using training area 2 labelled in (<b>a</b>); and (<b>d</b>) gradient of deprived area configurations from suburban to urban areas labelled in (<b>a</b>) by blue boxes from left to right (Image source: Pleiades image provided via the ESA Third Party data grant).</p>
Full article ">Figure 7
<p>The complexity of drawing boundaries. (<b>Left</b>) the example of Rio de Janeiro, Brazil (red line municipal favela boundary in 2016); (<b>Right</b>) visual delineations of boundaries of deprived areas in Ahmedabad, India by several geospatial experts (different colors refer to different interpreters) (Image source: Google Earth).</p>
Full article ">Figure 8
<p>Spatiotemporal dynamics of a deprived area in Bangalore (Image source: Google Earth)—the red boundaries represent the slum boundaries delineated in a city-wide mapping campaign by a local team in 2017 for the Dynaslum project.</p>
Full article ">Figure 9
<p>Main groups of accuracy assessment methods, which are relevant to assess mapping results of deprived areas.</p>
Full article ">Figure 10
<p>Difference between non-deprived (<b>1</b>) and deprived kampungs (<b>2</b>) (Image source: Google Earth image and ground photos J. Pratomo—adapted from [<a href="#B72-ijgi-07-00428" class="html-bibr">72</a>]).</p>
Full article ">Figure 11
<p>Difference between slum maps produced by the local mapping team of the Indonesian Government (<b>left</b>) and a combination of ground survey and image interpretation (<b>right</b>) [<a href="#B71-ijgi-07-00428" class="html-bibr">71</a>] (Image source: Pleiades provided via the ESA Third Party data grant).</p>
Full article ">Figure 12
<p>Modelling the urban population, Kalyan-Dombivali, India: (<b>a</b>) building footprints and slum boundaries, (<b>b</b>) DSM extracted from Cartosat-1 stereo-images (Image source: ISRO) combined with a Google Earth image, (<b>c</b>) estimated building height, and (<b>d</b>) estimated population density surface.</p>
Full article ">Figure 13
<p>Disappearance and growth of deprived areas driven by infrastructure and formal housing developments (red dot: example of an evicted area causing the incremental growth of deprived areas in the surrounding—purple dot) (Image source: Google Earth).</p>
Full article ">Figure 14
<p>Different spatial information levels of deprived areas: (<b>a</b>) deprived areas mapped together with other land cover/use; (<b>b</b>) detailed urban deprivation classes; (<b>c</b>) point locations of deprived areas and size; and (<b>d</b>) percentage of deprived areas within 200 × 200 m grid cells (Image source: WorldView-2 provided by DigitalGlobe).</p>
Full article ">
15 pages, 3106 KiB  
Article
Effect of Size, Shape and Map Background in Cartographic Visualization: Experimental Study on Czech and Chinese Populations
by Zdeněk Stachoň, Čeněk Šašinka, Jiří Čeněk, Stephan Angsüsser, Petr Kubíček, Zbyněk Štěrba and Martina Bilíková
ISPRS Int. J. Geo-Inf. 2018, 7(11), 427; https://doi.org/10.3390/ijgi7110427 - 1 Nov 2018
Cited by 19 | Viewed by 5662
Abstract
This paper deals with the issue of the perceptual aspects of selected graphic variables (specifically shape and size) and map background in cartographic visualization. The continued experimental study is based on previous findings and the presupposed cross-cultural universality of shape and size as [...] Read more.
This paper deals with the issue of the perceptual aspects of selected graphic variables (specifically shape and size) and map background in cartographic visualization. The continued experimental study is based on previous findings and the presupposed cross-cultural universality of shape and size as a graphic variable. The results bring a new perspective on the usage of shape, size and presence/absence of background as graphic variables, as well as a comparison to previous studies. The results suggest that all examined variables influence the speed of processing. Respondents (Czech and Chinese, N = 69) identified target stimuli faster without a map background, with larger stimuli, and with triangular and circular shapes. Czech respondents were universally faster than Chinese respondents. The implications of our research were discussed, and further directions were outlined. Full article
Show Figures

Figure 1

Figure 1
<p>Bertin’s six basic variables ([<a href="#B3-ijgi-07-00427" class="html-bibr">3</a>]; depiction adapted from [<a href="#B20-ijgi-07-00427" class="html-bibr">20</a>]).</p>
Full article ">Figure 2
<p>Minimum size (mm) at a reading distance of 40 cm (adapted from [<a href="#B34-ijgi-07-00427" class="html-bibr">34</a>]).</p>
Full article ">Figure 3
<p>Map reading speed (adapted from [<a href="#B34-ijgi-07-00427" class="html-bibr">34</a>]).</p>
Full article ">Figure 4
<p>Differences in the stimuli used with and without a map background (task instruction: find the shape shown at the left side of the screen).</p>
Full article ">Figure 5
<p>Schema of the procedure used for testing.</p>
Full article ">Figure 6
<p>Comparison of mean detection times with background and without background (in seconds, N = 69).</p>
Full article ">Figure 7
<p>Mean detection times regarding point symbol shape (in seconds, N = 69).</p>
Full article ">
19 pages, 26247 KiB  
Article
Application of Open-Source Software in Community Heritage Resources Management
by Jihn-Fa Jan
ISPRS Int. J. Geo-Inf. 2018, 7(11), 426; https://doi.org/10.3390/ijgi7110426 - 31 Oct 2018
Cited by 7 | Viewed by 6664
Abstract
In this paper, we present a case study of community heritage resources investigation and management, which was a collaborative project conducted by researchers and participants from rural communities. Geotagged photos were obtained using smart phones, and 360-degree panoramas were acquired using a robotic [...] Read more.
In this paper, we present a case study of community heritage resources investigation and management, which was a collaborative project conducted by researchers and participants from rural communities. Geotagged photos were obtained using smart phones, and 360-degree panoramas were acquired using a robotic camera system. These images were then uploaded to a web-based GIS (WebGIS) developed using Arches-Heritage Inventory Package (HIP), an open-source geospatial software system for cultural heritage inventory and management. By providing various tools for resources annotation, data exploration, mapping, geovisualization, and spatial analysis, the WebGIS not only serves as a platform for heritage resources database management, but also empowers the community residents to acquire, share, interpret, and analyze the data. The results show that this type of collaborative working model between researcher and community can promote public awareness of the importance of heritage conservation and achieve the research goal more effectively and efficiently. Full article
(This article belongs to the Special Issue Data Acquisition and Processing in Cultural Heritage)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The MAP VIEW of Nouli Community Arches platform showing locations of heritage sites on Bing Maps.</p>
Full article ">Figure 2
<p>The study site of this research. The overview image of Taiwan is from Google Earth, and the base map is from a WMS server created by the National Land Surveying and Mapping Center, Ministry of the Interior, Taiwan.</p>
Full article ">Figure 3
<p>Panoramic image of the of Hakka Culture Exhibition Center, which was obtained by using LizardQ camera system. The image is on Google Cloud so that viewers can observe the scene from every viewing angle. <a href="https://bit.ly/2PYxE7U" target="_blank">https://bit.ly/2PYxE7U</a>.</p>
Full article ">Figure 4
<p>A screen capture of a 360<math display="inline"><semantics> <msup> <mrow/> <mo>∘</mo> </msup> </semantics></math> video derived from panorama of the Hakka Culture Exhibition Center viewed in YouTube app. The video is on YouTube so that viewers can use VR headset to experience immersive virtual reality without the need to visit the site. <a href="https://youtu.be/x4PnMW1FG9U" target="_blank">https://youtu.be/x4PnMW1FG9U</a>.</p>
Full article ">Figure 5
<p>Definition of a placemark in a KML file.</p>
Full article ">Figure 6
<p>Microsoft Excel spreadsheets were used to enter heritage resource data, including site IDs, names, geometries, town IDs, and description.</p>
Full article ">Figure 7
<p>An example of “.arches” file, which contains attributes and geometry of the “Toyota Culture and History Exhibition Hall”.</p>
Full article ">Figure 8
<p>An example of “.relations” file, which defines relationships among heritage resources.</p>
Full article ">Figure 9
<p>The overall architecture of the Arches heritage management system.</p>
Full article ">Figure 10
<p>The Resource Data Manager of Arches-HIP.</p>
Full article ">Figure 11
<p>Heritage resources overlaid on historic map produced in 1897.</p>
Full article ">Figure 12
<p>The detail information of a heritage resource and the relation graph.</p>
Full article ">Figure 13
<p>The website uses scalable vector graphics (SVG) to show dynamic statistical chart for open data downloaded from government agencies. This figure depicts the immigration of Hakka people to the Nouli Community spanning a period of more than a century. <a href="http://arches.nccu.edu.tw/hakka" target="_blank">http://arches.nccu.edu.tw/hakka</a>.</p>
Full article ">
15 pages, 5922 KiB  
Article
Street Centralities and Land Use Intensities Based on Points of Interest (POI) in Shenzhen, China
by Shuai Wang, Gang Xu and Qingsheng Guo
ISPRS Int. J. Geo-Inf. 2018, 7(11), 425; https://doi.org/10.3390/ijgi7110425 - 31 Oct 2018
Cited by 29 | Viewed by 7045
Abstract
Urban land use and transportation are closely associated. Previous studies have investigated the spatial interrelationship between street centralities and land use intensities using land cover data, thus neglecting the social functions of urban land. Taking the city of Shenzhen, China, as a case [...] Read more.
Urban land use and transportation are closely associated. Previous studies have investigated the spatial interrelationship between street centralities and land use intensities using land cover data, thus neglecting the social functions of urban land. Taking the city of Shenzhen, China, as a case study, we used reclassified points of interest (POI) data to represent commercial, public service, and residential land, and then investigated the varying interrelationships between the street centralities and different types of urban land use intensities. We calculated three global centralities (“closeness”, “betweenness”, and “straightness”) as well as local centralities (1-km, 2-km, 3-km, and 5-km searching radiuses), which were transformed into raster frameworks using kernel density estimation (KDE) for correlation analysis. Global closeness and straightness are high in the urban core area, and roads with high global betweenness outline the skeleton of the street network. The spatial patterns of the local centralities are distinguished from the global centralities, reflecting local location advantages. High intensities of commercial and public service land are concentrated in the urban core, while residential land is relatively scattered. The bivariate correlation analysis implies that commercial and public service land are more dependent on centralities than residential land. Closeness and straightness have stronger abilities in measuring the location advantages than betweenness. The centralities and intensities are more positively correlated on a larger scale (census block). These findings of the spatial patterns and interrelationships of the centralities and intensities have major implications for urban land use and transportation planning. Full article
Show Figures

Figure 1

Figure 1
<p>Location of Shenzhen City in China and its spatial extent.</p>
Full article ">Figure 2
<p>Spatial distribution of the road net and road lengths in Shenzhen City (<b>a</b>) and density of nodes in Shenzhen City (<b>b</b>). The nodes density used the kernel density estimation (KDE) in a default searching radius.</p>
Full article ">Figure 3
<p>Numbers and proportions of 15 primary categories of points of interest (POIs) in Shenzhen City. These POIs are divided into commercial and business facilities (61.6%), residential sites (6.9%), administration and public services (8.1%), and others (23.4%).</p>
Full article ">Figure 4
<p>Flowchart of analysis of interrelationship between street centralities and POI-based urban land use intensities.</p>
Full article ">Figure 5
<p>Spatial distributions of global centrality indicators of street networks in Shenzhen. The centrality indicators of each edge are averaged from two nodes linking to this edge, and they are divided into five levels using natural breaks in ArcGIS10.2 (<b>a</b>,<b>c</b>,<b>e</b>); kernel density estimation of centrality indicators at nodes (<b>b</b>,<b>d</b>,<b>f</b>).</p>
Full article ">Figure 6
<p>Frequency distribution of global centrality indicators at each node, global closeness (<b>a</b>), global betweenness (<b>b</b>), and global straightness (<b>c</b>).</p>
Full article ">Figure 7
<p>Spatial distributions of intensities of commercial and business facilities (<b>a</b>), public and administrative services (<b>b</b>), and residential sites (<b>c</b>). These represent commercial land use intensities, public service land use intensities, and residential land use intensities, respectively.</p>
Full article ">Figure 8
<p>Frequency distributions of land use intensities using kernel density estimation (KDE) in a 1 km grid. (<b>a</b>) Commercial land use intensity, (<b>b</b>) public service land use intensity, (<b>c</b>) residential land use intensity.</p>
Full article ">Figure 9
<p>Scatter plots of street centralities and land use intensities. The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively. The same horizontal row shares the same land use intensity: commercial land, public service land, and residential land, from top to bottom, respectively. The straight lines are linear regression lines.</p>
Full article ">Figure 10
<p>Spatial patterns of KDE of the local street centralities using searching radiuses of 1 km, 2 km, 3 km, and 5 km. The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively. The same horizontal row shares the same searching radius: 1 km, 2 km, 3 km, and 5 km, from top to bottom, respectively.</p>
Full article ">Figure 11
<p>Spatial distributions of street centralities and land use intensities at the subdistrict level (census block) in Shenzhen City. (<b>a</b>) Closeness at 1 km searching radius, (<b>b</b>) Commercial land use intensity, (<b>c</b>) Global Betweenness, (<b>d</b>) Public service land use intensity, (<b>e</b>) Straightness at 2 km searching radius, (<b>f</b>) Residential land use intensity.</p>
Full article ">
26 pages, 15439 KiB  
Article
Assessment of Segmentation Parameters for Object-Based Land Cover Classification Using Color-Infrared Imagery
by Ozgun Akcay, Emin Ozgur Avsar, Melis Inalpulat, Levent Genc and Ahmet Cam
ISPRS Int. J. Geo-Inf. 2018, 7(11), 424; https://doi.org/10.3390/ijgi7110424 - 31 Oct 2018
Cited by 15 | Viewed by 5821
Abstract
Using object-based image analysis (OBIA) techniques for land use-land cover classification (LULC) has become an area of interest due to the availability of high-resolution data and segmentation methods. Multi-resolution segmentation in particular, statistically seen as the most used algorithm, is able to produce [...] Read more.
Using object-based image analysis (OBIA) techniques for land use-land cover classification (LULC) has become an area of interest due to the availability of high-resolution data and segmentation methods. Multi-resolution segmentation in particular, statistically seen as the most used algorithm, is able to produce non-identical segmentations depending on the required parameters. The total effect of segmentation parameters on the classification accuracy of high-resolution imagery is still an open question, though some studies were implemented to define the optimum segmentation parameters. However, recent studies have not properly considered the parameters and their consequences on LULC accuracy. The main objective of this study is to assess OBIA segmentation and classification accuracy according to the segmentation parameters using different overlap ratios during image object sampling for a predetermined scale. With this aim, we analyzed and compared (a) high-resolution color-infrared aerial images of a newly-developed urban area including different land use types; (b) combinations of multi-resolution segmentation with different shape, color, compactness, bands, and band-weights; and (c) accuracies of classifications based on varied segmentations. The results of various parameters in the study showed an explicit correlation between segmentation accuracies and classification accuracies. The effect of changes in segmentation parameters using different sample selection methods for five main LULC types was studied. Specifically, moderate shape and compactness values provided more consistency than lower and higher values; also, band weighting demonstrated substantial results due to the chosen bands. Differences in the variable importance of the classifications and changes in LULC maps were also explained. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study area: (<b>a</b>–<b>c</b>) Google Earth imagery, (<b>d</b>) ortho-mosaic image.</p>
Full article ">Figure 2
<p>Sample distribution of LULC classes.</p>
Full article ">Figure 3
<p>Workflow diagram.</p>
Full article ">Figure 4
<p>Number of segments according to segmentation parameters.</p>
Full article ">Figure 5
<p>Comparison of segmentation attempts: (<b>a</b>) southeast of the study area, (<b>b</b>) SN 1, (<b>c</b>) SN 73, (<b>d</b>) SN 97, (<b>e</b>) SN 98, and (<b>f</b>) SN 101.</p>
Full article ">Figure 6
<p>Defining training segments using different criteria.</p>
Full article ">Figure 7
<p>Class, sample and segment relations.</p>
Full article ">Figure 8
<p>(<b>a</b>) Over segmentation, (<b>b</b>) area fit index, (<b>c</b>) 1-Quality rate, and (<b>d</b>) root mean square error of segmentations.</p>
Full article ">Figure 8 Cont.
<p>(<b>a</b>) Over segmentation, (<b>b</b>) area fit index, (<b>c</b>) 1-Quality rate, and (<b>d</b>) root mean square error of segmentations.</p>
Full article ">Figure 8 Cont.
<p>(<b>a</b>) Over segmentation, (<b>b</b>) area fit index, (<b>c</b>) 1-Quality rate, and (<b>d</b>) root mean square error of segmentations.</p>
Full article ">Figure 9
<p>Number of sample-segments for each class in classifications based on SSM 3.</p>
Full article ">Figure 10
<p>User accuracies for each class in classifications based on SSM 3.</p>
Full article ">Figure 11
<p>Producer accuracies for each class in classifications based on SSM 3.</p>
Full article ">Figure 12
<p>Kappa values according to segmentation numbers.</p>
Full article ">Figure 13
<p>Comparison of classification accuracies of equal and unequal number of sample-segment classifications based on SSM 3.</p>
Full article ">Figure 14
<p>Mean decrease accuracy for variables.</p>
Full article ">Figure 15
<p>Classification examples.</p>
Full article ">Figure 16
<p>Percentage of change in class assignment in classifications based on SSM 3.</p>
Full article ">Figure 17
<p>Classification kappa values and segmentation RMSE values according to segment selection method, shape and compactness criteria and Euclidean distances (<b>a</b>) SSM 1, (<b>b</b>) SSM 2, (<b>c</b>) SSM 3, and (<b>d</b>) legend.</p>
Full article ">Figure 17 Cont.
<p>Classification kappa values and segmentation RMSE values according to segment selection method, shape and compactness criteria and Euclidean distances (<b>a</b>) SSM 1, (<b>b</b>) SSM 2, (<b>c</b>) SSM 3, and (<b>d</b>) legend.</p>
Full article ">
18 pages, 6298 KiB  
Article
A Task-Oriented Knowledge Base for Geospatial Problem-Solving
by Can Zhuang, Zhong Xie, Kai Ma, Mingqiang Guo and Liang Wu
ISPRS Int. J. Geo-Inf. 2018, 7(11), 423; https://doi.org/10.3390/ijgi7110423 - 31 Oct 2018
Cited by 6 | Viewed by 5112
Abstract
In recent years, the rapid development of cloud computing and web technologies has led to a significant advancement to chain geospatial information services (GI services) in order to solve complex geospatial problems. However, the construction of a problem-solving workflow requires considerable expertise for [...] Read more.
In recent years, the rapid development of cloud computing and web technologies has led to a significant advancement to chain geospatial information services (GI services) in order to solve complex geospatial problems. However, the construction of a problem-solving workflow requires considerable expertise for end-users. Currently, few studies design a knowledge base to capture and share geospatial problem-solving knowledge. This paper abstracts a geospatial problem as a task that can be further decomposed into multiple subtasks. The task distinguishes three distinct granularities: Geooperator, Atomic Task, and Composite Task. A task model is presented to define the outline of problem solution at a conceptual level that closely reflects the processes for problem-solving. A task-oriented knowledge base that leverages an ontology-based approach is built to capture and share task knowledge. This knowledge base provides the potential for reusing task knowledge when faced with a similar problem. Conclusively, the details of implementation are described through using a meteorological early-warning analysis as an example. Full article
(This article belongs to the Special Issue Big Data Computing for Geospatial Applications)
Show Figures

Figure 1

Figure 1
<p>Sequences of the meteorological early-warning (MEW) process.</p>
Full article ">Figure 2
<p>Relationships between Task, AtomicTask, CompositeTask, and Geooperator.</p>
Full article ">Figure 3
<p>(<b>a</b>) Different perspectives on Geooperator [<a href="#B41-ijgi-07-00423" class="html-bibr">41</a>] (<b>b</b>) Description elements of Geooperator [<a href="#B17-ijgi-07-00423" class="html-bibr">17</a>].</p>
Full article ">Figure 4
<p>An example of task decomposition.</p>
Full article ">Figure 5
<p>The relationships of ontologies in the knowledge base.</p>
Full article ">Figure 6
<p>Data type specifications.</p>
Full article ">Figure 7
<p>Interface for annotating Task and Geooperator, and Data Type for specifying Interface.</p>
Full article ">Figure 8
<p>An excerpt of ontologies where (<b>a</b>) depicts the classes of ontologies, and (<b>b</b>) illustrates the object properties between classes.</p>
Full article ">Figure 9
<p>Snippets of owl notation using a universal restriction.</p>
Full article ">Figure 10
<p>The task instance of an atomic task (EffectiveRainfallCalTask).</p>
Full article ">Figure 11
<p>The task instance of a composite task (EarlyWarningAnalysisTask).</p>
Full article ">Figure 12
<p>The workflow model of EWATask. Data access services are represented in elliptical shapes, and geoprocessing services are represented in rectangular shapes.</p>
Full article ">Figure 13
<p>The graphical user interface of the prototype system.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop