Nothing Special   »   [go: up one dir, main page]

Next Issue
Volume 13, January
Previous Issue
Volume 12, November
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 12, Issue 12 (December 2023) – 38 articles

Cover Story (view full-size image): Three-dimensional indoor models are a crucial component for simulating pedestrian evacuations. However, existing 3D indoor models cannot fully represent indoor environments to simulate 3D pedestrian motions in evacuations since spaces above/below some physical components (e.g., desks, chairs) are largely overlooked. This paper presents a conceptual space model that fully captures six 3D pedestrian motions (e.g., low crawling, bent-over walking). We first present the definition and parameterisation of the motions. Next, the definition and specifications of three types of space components are articulated based on the motions. Finally, the concepts are implemented using a voxel-based approach and demonstrated using an illustrative example. This work advances 3D indoor modelling towards realistically simulating 3D pedestrian motions. View this paper
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Select all
Export citation of selected articles as:
19 pages, 5755 KiB  
Article
Hyperspectral Image Classification Network Based on 3D Octave Convolution and Multiscale Depthwise Separable Convolution
by Qingqing Hong, Xinyi Zhong, Weitong Chen, Zhenghua Zhang and Bin Li
ISPRS Int. J. Geo-Inf. 2023, 12(12), 505; https://doi.org/10.3390/ijgi12120505 - 17 Dec 2023
Cited by 3 | Viewed by 2328
Abstract
Hyperspectral images (HSIs) are pivotal in various fields due to their rich spectral–spatial information. While convolutional neural networks (CNNs) have notably enhanced HSI classification, they often generate redundant spatial features. To address this, we introduce a novel HSI classification method, OMDSC, employing 3D [...] Read more.
Hyperspectral images (HSIs) are pivotal in various fields due to their rich spectral–spatial information. While convolutional neural networks (CNNs) have notably enhanced HSI classification, they often generate redundant spatial features. To address this, we introduce a novel HSI classification method, OMDSC, employing 3D Octave convolution combined with multiscale depthwise separable convolutional networks. This method initially utilizes 3D Octave convolution for efficient spectral–spatial feature extraction from HSIs, thereby reducing spatial redundancy. Subsequently, multiscale depthwise separable convolution is used to further improve the extraction of spatial features. Finally, the HSI classification results are output by softmax classifier. This work compares the method with other methods on three publicly available datasets in order to confirm its efficacy. The outcomes show that the method performs better in terms of classification. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the 3D Octave convolution.</p>
Full article ">Figure 2
<p>(<b>a</b>) Diagram of the depthwise separable convolution process. (<b>b</b>) Diagram of the ordinary convolution process.</p>
Full article ">Figure 3
<p>Overall schematic diagram of the proposed method.</p>
Full article ">Figure 4
<p>Multiscale depthwise separable convolution-specific parameter map.</p>
Full article ">Figure 5
<p>Technical flowchart of the method.</p>
Full article ">Figure 6
<p>Classification maps generated by all of the competing methods on the Indian Pines data with 10% training samples. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">Figure 6 Cont.
<p>Classification maps generated by all of the competing methods on the Indian Pines data with 10% training samples. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">Figure 7
<p>Classification maps generated by all of the competing methods on the University of Pavia data with 5% training samples. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">Figure 8
<p>Classification maps generated by all of the competing methods on the WHU-Hi-LongKou dataset with 1% training samples. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">Figure 9
<p>Confusion matrix of different methods for the Indian Pines dataset. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">Figure 10
<p>Confusion matrix of different methods for the University of Pavia dataset. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">Figure 11
<p>Confusion matrix of different methods for the Salinas Scene dataset. (<b>a</b>) 2DCNN, (<b>b</b>) 3DCNN, (<b>c</b>) M3D-DCNN, (<b>d</b>) HybridSN, (<b>e</b>) Vit, (<b>f</b>) SATNet, (<b>g</b>) SSFTT, (<b>h</b>) OMDSC.</p>
Full article ">
23 pages, 13895 KiB  
Article
A Web-Based Geodesign Tool for Evaluating the Integration of Transport Infrastructure, Public Spaces, and Human Activities
by Liu Yang
ISPRS Int. J. Geo-Inf. 2023, 12(12), 504; https://doi.org/10.3390/ijgi12120504 - 17 Dec 2023
Viewed by 2354
Abstract
The need for addressing the adverse impacts of transport infrastructure on public spaces and human activities (TSH) emphasizes the importance of designing integrated TSH system, thereby necessitating tailored planning support systems (PSS). This study begins by assessing the demand for PSS using surveys [...] Read more.
The need for addressing the adverse impacts of transport infrastructure on public spaces and human activities (TSH) emphasizes the importance of designing integrated TSH system, thereby necessitating tailored planning support systems (PSS). This study begins by assessing the demand for PSS using surveys and interviews to uncover the need for robust analysis and evaluation support, particularly through the use of geographical information systems (GIS). On this basis, a prototype GIS platform is proposed for analyzing and evaluating the integration of the TSH system at the block scale. This user-friendly geodesign tool encompasses a customizable evaluation index (includes seven KPAs and KPIs), allowing for combined quantitative and qualitative assessments. Notably, it introduces a buffer effect index to quantify transport–space interaction. The proposed tool serves as a dedicated platform for evaluating TSH systems, offering 2D/3D visualization capabilities and two analysis units and facilitating cross-platform collaboration. Applied to a case study in Nanjing, China, it effectively assessed the interdependence among different TSH system components and block integration around expressways, railways, and main roads. This tool holds promise in offering invaluable insights into urban planning and (re)development, thereby enhancing the integration of transport infrastructure and public spaces. Full article
Show Figures

Figure 1

Figure 1
<p>Tool demands of different urban planning stages.</p>
Full article ">Figure 2
<p>(<b>left</b>) Needs of quantitative methods in different urban design stages; (<b>right</b>) The most used tools (word size reflects the number of respondents who mentioned it).</p>
Full article ">Figure 3
<p>The indicators that should be quantified in an urban plan evaluation.</p>
Full article ">Figure 4
<p>The architecture of the tool.</p>
Full article ">Figure 5
<p>The operation process of the tool.</p>
Full article ">Figure 6
<p>The indicator set for evaluating the TSH system, adapted from Yang et al. [<a href="#B33-ijgi-12-00504" class="html-bibr">33</a>]. Note: Bold red items are KPAs selected for this study.</p>
Full article ">Figure 7
<p>The interface and use process of the Web-based GIS tool.</p>
Full article ">Figure 8
<p>Location of the case study blocks in Nanjing, China (coordinate system: WGS 84/UTM zone 50N).</p>
Full article ">Figure 9
<p>An example of the evaluation index/3D model (<b>left</b>) and evaluation result/2D data (<b>right</b>).</p>
Full article ">Figure 10
<p>Integration scores of all blocks in the three study areas and position of 24 selected blocks (coordinate system: WGS 84/UTM zone 50N).</p>
Full article ">Figure 11
<p>Score and KPI values of each area (<b>top</b>) and the 24 selected blocks (<b>bottom</b>).</p>
Full article ">Figure 12
<p>The correlation among the KPIs. Note: Significance level: *** <span class="html-italic">p</span> &lt; 0.01, ** <span class="html-italic">p</span> &lt; 0.05, and * <span class="html-italic">p</span> &lt; 0.1. The orange border highlights strong positive correlations among KPIs.</p>
Full article ">
36 pages, 15617 KiB  
Article
Machine-Learning-Based Forest Classification and Regression (FCR) for Spatial Prediction of Liver Fluke Opisthorchis viverrini (OV) Infection in Small Sub-Watersheds
by Benjamabhorn Pumhirunroj, Patiwat Littidej, Thidarut Boonmars, Kanokwan Bootyothee, Atchara Artchayasawat, Phusit Khamphilung and Donald Slack
ISPRS Int. J. Geo-Inf. 2023, 12(12), 503; https://doi.org/10.3390/ijgi12120503 - 14 Dec 2023
Cited by 4 | Viewed by 2719
Abstract
Infection of liver flukes (Opisthorchis viverrini) is partly due to their suitability for habitats in sub-basin areas, which causes the intermediate host to remain in the watershed system in all seasons. The spatial monitoring of fluke at the small basin scale [...] Read more.
Infection of liver flukes (Opisthorchis viverrini) is partly due to their suitability for habitats in sub-basin areas, which causes the intermediate host to remain in the watershed system in all seasons. The spatial monitoring of fluke at the small basin scale is important because this can enable analysis at the level of the factors involved that influence infections. A spatial mathematical model was weighted by the nine spatial factors X1 (index of land-use types), X2 (index of soil drainage properties), X3 (distance index from the road network, X4 (distance index from surface water resources), X5 (distance index from the flow accumulation lines), X6 (index of average surface temperature), X7 (average surface moisture index), X8 (average normalized difference vegetation index), and X9 (average soil-adjusted vegetation index) by dividing the analysis into two steps: (1) the sub-basin boundary level was analyzed with an ordinary least square (OLS) model used to select the spatial criteria of liver flukes aimed at analyzing the factors related to human liver fluke infection according to sub-watersheds, and (2) we used the infection risk positional analysis level through machine-learning-based forest classification and regression (FCR) to display the predictive results of infection risk locations along stream lines. The analysis results show four prototype models that import different independent variable factors. The results show that Model 1 and Model 2 gave the most AUC (0.964), and the variables that influenced infection risk the most were the distance to stream lines and the distance to water bodies; the NDMI and NDVI factors rarely affected the accuracy. This FCR machine-learning application approach can be applied to the analysis of infection risk areas at the sub-basin level, but independent variables must be screened with a preliminary mathematical model weighted to the spatial units in order to obtain the most accurate predictions. Full article
Show Figures

Figure 1

Figure 1
<p>The boundaries of the study area show the proximity of freshwater bodies that are fish habitats to the Mekong River.</p>
Full article ">Figure 2
<p>Percentage of people infected with liver flukes during 2019–2021 according to the 8th Regional Health Province (R8) near the Mekong River (adapted from R8, [<a href="#B80-ijgi-12-00503" class="html-bibr">80</a>]. <a href="https://r8way.moph.go.th/r8-primary/" target="_blank">https://r8way.moph.go.th/r8-primary/</a> (accessed on 20 June 2022).</p>
Full article ">Figure 3
<p>(<b>a</b>) The number of infected populations of each village with population density (persons/sq·km) (<b>b</b>) Infected percentage of each sub-basin and the Nongharn Lake boundary in the rainy season.</p>
Full article ">Figure 4
<p>The framework of OLS modeling finds the relationship between liver fluke occurrence and spatial factors at the sub-basin level.</p>
Full article ">Figure 5
<p>Pareto plots for screening sub-basin area units to determine the location of machine learning for importing forest-based classification and regression models: (<b>a</b>) Y (% of OV), (<b>b</b>) <span class="html-italic">X</span><sub>1</sub> (Land use), (<b>c</b>) <span class="html-italic">X</span><sub>2</sub> (Soil), (<b>d</b>) <span class="html-italic">X</span><sub>3</sub> (Road), (<b>e</b>) <span class="html-italic">X</span><sub>4</sub> (Water), (<b>f</b>) <span class="html-italic">X</span><sub>5</sub> (Stream), (<b>g</b>) <span class="html-italic">X</span><sub>6</sub> (Surface Temperature), (<b>h</b>) <span class="html-italic">X</span><sub>7</sub> (NDMI), (<b>i</b>) <span class="html-italic">X</span><sub>8</sub> (NDVI), and (<b>j</b>) <span class="html-italic">X</span><sub>9</sub> (SAVI).</p>
Full article ">Figure 5 Cont.
<p>Pareto plots for screening sub-basin area units to determine the location of machine learning for importing forest-based classification and regression models: (<b>a</b>) Y (% of OV), (<b>b</b>) <span class="html-italic">X</span><sub>1</sub> (Land use), (<b>c</b>) <span class="html-italic">X</span><sub>2</sub> (Soil), (<b>d</b>) <span class="html-italic">X</span><sub>3</sub> (Road), (<b>e</b>) <span class="html-italic">X</span><sub>4</sub> (Water), (<b>f</b>) <span class="html-italic">X</span><sub>5</sub> (Stream), (<b>g</b>) <span class="html-italic">X</span><sub>6</sub> (Surface Temperature), (<b>h</b>) <span class="html-italic">X</span><sub>7</sub> (NDMI), (<b>i</b>) <span class="html-italic">X</span><sub>8</sub> (NDVI), and (<b>j</b>) <span class="html-italic">X</span><sub>9</sub> (SAVI).</p>
Full article ">Figure 6
<p>The forest-based classification and regression (FCR) for liver fluke (<span class="html-italic">Opisthorchis viverrini</span>)-infection prediction.</p>
Full article ">Figure 7
<p>The whole dataset for training and testing.</p>
Full article ">Figure 8
<p>Comparison of the standard residuals of OLS alternative models: (<b>a</b>) OLS Model 1, (<b>b</b>) OLS Model 2, (<b>c</b>) OLS Model 3, and (<b>d</b>) OLS Model 4.</p>
Full article ">Figure 9
<p>Sub-basin boundary map obtained from the analysis of the DEM data.</p>
Full article ">Figure 10
<p>Raster mapping radius of <span class="html-italic">OV</span>-infected points using a heatmap: (<b>a</b>) radius: 2 km, (<b>b</b>) radius: 4 km.</p>
Full article ">Figure 11
<p>Maps of independent variable indexes, where (<b>a</b>) is <span class="html-italic">X</span><sub>1</sub> (index of land-use types), (<b>b</b>) is <span class="html-italic">X</span><sub>2</sub> (index of soil drainage properties), (<b>c</b>) is <span class="html-italic">X</span><sub>3</sub> (the distance index from the road network, (<b>d</b>) is <span class="html-italic">X</span><sub>4</sub> (distance index from surface water sources), (<b>e</b>) is <span class="html-italic">X</span><sub>5</sub> (distance index from the flow accumulation lines), (<b>f</b>) is <span class="html-italic">X</span><sub>6</sub> (index of average surface temperature), (<b>g</b>) is <span class="html-italic">X</span><sub>7</sub> (average surface moisture index), (<b>h</b>) is <span class="html-italic">X</span><sub>8</sub> (average normalized difference vegetation index), and (<b>i</b>) is <span class="html-italic">X</span><sub>9</sub> (average soil-adjusted vegetation index).</p>
Full article ">Figure 11 Cont.
<p>Maps of independent variable indexes, where (<b>a</b>) is <span class="html-italic">X</span><sub>1</sub> (index of land-use types), (<b>b</b>) is <span class="html-italic">X</span><sub>2</sub> (index of soil drainage properties), (<b>c</b>) is <span class="html-italic">X</span><sub>3</sub> (the distance index from the road network, (<b>d</b>) is <span class="html-italic">X</span><sub>4</sub> (distance index from surface water sources), (<b>e</b>) is <span class="html-italic">X</span><sub>5</sub> (distance index from the flow accumulation lines), (<b>f</b>) is <span class="html-italic">X</span><sub>6</sub> (index of average surface temperature), (<b>g</b>) is <span class="html-italic">X</span><sub>7</sub> (average surface moisture index), (<b>h</b>) is <span class="html-italic">X</span><sub>8</sub> (average normalized difference vegetation index), and (<b>i</b>) is <span class="html-italic">X</span><sub>9</sub> (average soil-adjusted vegetation index).</p>
Full article ">Figure 11 Cont.
<p>Maps of independent variable indexes, where (<b>a</b>) is <span class="html-italic">X</span><sub>1</sub> (index of land-use types), (<b>b</b>) is <span class="html-italic">X</span><sub>2</sub> (index of soil drainage properties), (<b>c</b>) is <span class="html-italic">X</span><sub>3</sub> (the distance index from the road network, (<b>d</b>) is <span class="html-italic">X</span><sub>4</sub> (distance index from surface water sources), (<b>e</b>) is <span class="html-italic">X</span><sub>5</sub> (distance index from the flow accumulation lines), (<b>f</b>) is <span class="html-italic">X</span><sub>6</sub> (index of average surface temperature), (<b>g</b>) is <span class="html-italic">X</span><sub>7</sub> (average surface moisture index), (<b>h</b>) is <span class="html-italic">X</span><sub>8</sub> (average normalized difference vegetation index), and (<b>i</b>) is <span class="html-italic">X</span><sub>9</sub> (average soil-adjusted vegetation index).</p>
Full article ">Figure 12
<p>Display receiver operating characteristic (ROC) curve and area under the ROC curve (AUC), comparing four models: (<b>a</b>) Model 1, (<b>b</b>) Model 2, (<b>c</b>) Model 3, (<b>d</b>) Model 4, and (<b>e</b>) all models.</p>
Full article ">Figure 13
<p>The results of the predictions for FCR for Models 1–4: (<b>a</b>) Model 1, (<b>b</b>) Model 2, (<b>c</b>) Model 3, and (<b>d</b>) Model 4.</p>
Full article ">Figure 14
<p>Residual plot and fit-plot graphs of variable correlations: <span class="html-italic">X</span><sub>5</sub> (<b>a1</b>,<b>a2</b>); <span class="html-italic">X</span><sub>6</sub> (<b>b1</b>,<b>b2</b>); <span class="html-italic">X</span><sub>7</sub> (<b>c1</b>,<b>c2</b>); <span class="html-italic">X</span><sub>8</sub> (<b>d1</b>,<b>d2</b>); and <span class="html-italic">X</span><sub>9</sub> (<b>e1</b>,<b>e2</b>), selected from correlation analysis.</p>
Full article ">Figure 14 Cont.
<p>Residual plot and fit-plot graphs of variable correlations: <span class="html-italic">X</span><sub>5</sub> (<b>a1</b>,<b>a2</b>); <span class="html-italic">X</span><sub>6</sub> (<b>b1</b>,<b>b2</b>); <span class="html-italic">X</span><sub>7</sub> (<b>c1</b>,<b>c2</b>); <span class="html-italic">X</span><sub>8</sub> (<b>d1</b>,<b>d2</b>); and <span class="html-italic">X</span><sub>9</sub> (<b>e1</b>,<b>e2</b>), selected from correlation analysis.</p>
Full article ">
22 pages, 8890 KiB  
Article
Construction of a Real-Time Ship Trajectory Prediction Model Based on Ship Automatic Identification System Data
by Daping Xi, Yuhao Feng, Wenping Jiang, Nai Yang, Xini Hu and Chuyuan Wang
ISPRS Int. J. Geo-Inf. 2023, 12(12), 502; https://doi.org/10.3390/ijgi12120502 - 13 Dec 2023
Cited by 1 | Viewed by 2356
Abstract
The extraction of ship behavior patterns from Automatic Identification System (AIS) data and the subsequent prediction of travel routes play crucial roles in mitigating the risk of ship accidents. This study focuses on the Wuhan section of the dendritic river system in the [...] Read more.
The extraction of ship behavior patterns from Automatic Identification System (AIS) data and the subsequent prediction of travel routes play crucial roles in mitigating the risk of ship accidents. This study focuses on the Wuhan section of the dendritic river system in the middle reaches of the Yangtze River and the partial reticulated river system in the northern part of the Zhejiang Province as its primary investigation areas. Considering the structure and attributes of AIS data, we introduce a novel algorithm known as the Combination of DBSCAN and DTW (CDDTW) to identify regional navigation characteristics of ships. Subsequently, we develop a real-time ship trajectory prediction model (RSTPM) to facilitate real-time ship trajectory predictions. Experimental tests on two distinct types of river sections are conducted to assess the model’s reliability. The results indicate that the RSTPM exhibits superior prediction accuracy when compared to conventional trajectory prediction models, achieving an approximate 20 m prediction accuracy for ship trajectories on inland waterways. This showcases the advancements made by this model. Full article
Show Figures

Figure 1

Figure 1
<p>Technical roadmap for the RSTPM.</p>
Full article ">Figure 2
<p>Segmentation results of Ball-Tree with different segmentation thresholds (<span class="html-italic">N</span> = 3, <b>left</b>) (<span class="html-italic">N</span> = 2, <b>right</b>).</p>
Full article ">Figure 3
<p>The influence of two clustering parameters on the number of clusters and the number of noise points.</p>
Full article ">Figure 4
<p>Improvement of the algorithm.</p>
Full article ">Figure 5
<p>Two different ship sailing tracks.</p>
Full article ">Figure 6
<p>Gray similarity distance matrix and matching path of two trajectories before and after the improvement of the DTW algorithm.</p>
Full article ">Figure 7
<p>Structure of the RSTPM.</p>
Full article ">Figure 8
<p>Nonforked section of the dendritic river system in the Yangtze River’s middle reaches in Wuhan.</p>
Full article ">Figure 9
<p>Some multiforked river sections of reticulated river systems in northern Zhejiang Province.</p>
Full article ">Figure 10
<p>Study area and numbering.</p>
Full article ">Figure 11
<p>Mean squared error for different batch sizes in the training set and validation set.</p>
Full article ">Figure 12
<p>Variation in mean squared error and accuracy for the model with the number of neurons.</p>
Full article ">Figure 13
<p>Variation of mean squared error and accuracy of the model with step size.</p>
Full article ">Figure 14
<p>Channel names and ship tracks in the study area.</p>
Full article ">Figure 15
<p>Vessel trajectory clustering results.</p>
Full article ">Figure 16
<p>Normal trajectory prediction results and prediction error.</p>
Full article ">Figure 17
<p>Prediction of trajectory within 30 min periods in different positions.</p>
Full article ">Figure 18
<p>The average prediction error and the maximum error for the 30 min period across the entire trajectory.</p>
Full article ">Figure 19
<p>Location of multifork river in the study area.</p>
Full article ">Figure 20
<p>Ship history track.</p>
Full article ">Figure 21
<p>Predicted trajectory on multiple segments when the ship reaches the bifurcation.</p>
Full article ">Figure 22
<p>Predicted trajectory of the ship after passing the bifurcation.</p>
Full article ">
23 pages, 9627 KiB  
Article
The Spatial Effects of Regional Poverty: Spatial Dependence, Spatial Heterogeneity and Scale Effects
by Mengxiao Liu, Yong Ge, Shan Hu and Haiguang Hao
ISPRS Int. J. Geo-Inf. 2023, 12(12), 501; https://doi.org/10.3390/ijgi12120501 - 13 Dec 2023
Cited by 2 | Viewed by 2439
Abstract
Recognizing the spatial effects of regional poverty is essential for achieving sustainable poverty alleviation. This study investigates these spatial effects and their determinants across three distinct administrative levels within Hubei Province, China. To analyze the spatial patterns and heterogeneity of multi-scale regional poverty, [...] Read more.
Recognizing the spatial effects of regional poverty is essential for achieving sustainable poverty alleviation. This study investigates these spatial effects and their determinants across three distinct administrative levels within Hubei Province, China. To analyze the spatial patterns and heterogeneity of multi-scale regional poverty, we employed various spatial analysis techniques, including the global and local Moran’s I statistics, the Lineman, Merenda, and Gold (LMG) method, as well as Multiscale Geographically Weighted Regression (MGWR). We found that: (1) Regional poverty exhibits significant spatial dependence across various scales, with a higher level of spatial dependence observed at higher administrative levels. (2) The spatial distribution of poverty is primarily influenced by geographical factors, encompassing first-, second-, and third-nature geographical elements. Notably, first-nature geographical factors make substantial contributions, accounting for 36.99%, 42.23%, and 23.79% at the county, township, and village levels, respectively. (3) The influence of geographical factors varies with scale. Global effects of various factors may transcend scales or remain confined to specific scales, while the local impacts of different factors also exhibit variations across scales. These results underscore the necessity for collaborative efforts among government entities at different levels with the anti-poverty measures tailored to local contexts. Full article
Show Figures

Figure 1

Figure 1
<p>The location of Hubei Province in China; the county and township levels poverty incidence in Hubei province; and the village–level poverty incidence in Yunyang County in 2013.</p>
Full article ">Figure 2
<p>Conceptual framework of the formation of the spatial effects of regional poverty.</p>
Full article ">Figure 3
<p>Variation in LISA of regional poverty at different levels.</p>
Full article ">Figure 4
<p>Relative importance of geographical variables explaining poverty at different levels (metrics first normalized to add up to 100% and then multiplied by the adjusted R<sup>2</sup> to represent their real contribution at each scale).</p>
Full article ">Figure 5
<p>Bivariate relations between poverty incidence and individual significant geographical factors.</p>
Full article ">Figure 5 Cont.
<p>Bivariate relations between poverty incidence and individual significant geographical factors.</p>
Full article ">Figure 6
<p>Ranges of the regression coefficient of different geographical factors in MGWR models at different scales.</p>
Full article ">Figure 7
<p>Local estimates for (<b>a</b>) average elevation, (<b>b</b>) distance to the nearest administrative center, (<b>c</b>) proportion of the population having access to the internet at different levels.</p>
Full article ">
22 pages, 12065 KiB  
Article
Geo-Referencing and Analysis of Entities Extracted from Old Drawings and Photos Using Computer Vision and Deep Learning Algorithms
by Liat David, Motti Zohar and Ilan Shimshoni
ISPRS Int. J. Geo-Inf. 2023, 12(12), 500; https://doi.org/10.3390/ijgi12120500 - 13 Dec 2023
Cited by 2 | Viewed by 1925
Abstract
This study offers a quantitative solution that automates the creation of a historical timeline starting with old drawings from the beginning of the 18th century and ending with present-day photographs of the Old City of Jerusalem. This is performed using GIScience approaches, computer [...] Read more.
This study offers a quantitative solution that automates the creation of a historical timeline starting with old drawings from the beginning of the 18th century and ending with present-day photographs of the Old City of Jerusalem. This is performed using GIScience approaches, computer vision, and deep learning. The motivation to select the Old City of Jerusalem is the substantial availability of old archival drawings and photographs, owing to the area’s significance throughout the years. This task is challenging, as drawings, old photographs, and new photographs exhibit distinct characteristics. Our method encompasses several key components for the analyses: a 2D location recommendation engine, which detects an approximate location in the image of 3D landmarks; 2D landmarks to 3D conversion; and 2D polygonal areas to 3D GIS polylines conversion. This is applied to the segmentation of built areas. To achieve more accurate results, Meta’s Segment Anything model was utilized, which eliminates the need for extensive data preparation, training, and validation, thus optimizing the process. Using such techniques enabled us to examine the landscape development throughout the last three centuries and gain deeper insights concerning the evolution of prominent landmarks and features such as built area over time. Full article
Show Figures

Figure 1

Figure 1
<p>A diagram describing the various components of the algorithms presented in the paper.</p>
Full article ">Figure 2
<p>For every 2D pixel <span class="html-italic">p</span> (landmark) in a drawing/photo, <math display="inline"><semantics> <mover accent="true"> <mi>P</mi> <mo>˜</mo> </mover> </semantics></math> is calculated from the smallest group of <math display="inline"><semantics> <mi>α</mi> </semantics></math> with the shortest <span class="html-italic">d</span>.</p>
Full article ">Figure 3
<p>Segmentation using the Segment Anything model and standard semantic segmentation. (<b>a</b>,<b>c</b>) shows squares that mark all the buildings in the photo/drawing. (<b>b</b>,<b>d</b>) Segmentation results after applying the Segment Anything model to the photo/drawing. (<b>e</b>) Example of results in <a href="http://apeer.com" target="_blank">apeer.com</a> (accessed on 27 September 2023), U-Net, FPN, and LinkNet semantic segmentation networks on the drawings dataset that was created.</p>
Full article ">Figure 4
<p>Step-by-step process for creating 3D GIS areas from photo/drawing buildings. The figure visually presents the sequential steps involved in transforming 2D representation into a 3D model.</p>
Full article ">Figure 5
<p>Evolution of the Church of the Redeemer (Z) and Custodia Terrae Sanctae (Q). (<b>a</b>) Vue du Mont des Oliviers, centre, Panorama (1870–1871) does not include the landmarks Z and Q. (<b>b</b>) Panorama photo American Colony (1899–1904), landmarks Z and Q can be observed.</p>
Full article ">Figure 6
<p>Evolution of Dormition Abbey (A2). (<b>a</b>) Panorama photo—American Colony taken between 1899–1904 does not include landmark A2. (<b>b</b>) American Colony photo taken in 1920, we can observe the landmark A2.</p>
Full article ">Figure 7
<p>Two examples of good results from the recommendation engine. Al-Hanka minaret (M1) marked in (<b>a</b>) DSC_0518 (2020) and suggested in (<b>b</b>) American Colony (1899–1904) and landmark Al-Omar minaret (I) marked in (<b>c</b>) DSC_0518 (2020) and suggested in (<b>d</b>) Gotfryd, Bernard, photographer (1971).</p>
Full article ">Figure 8
<p>Part (<b>b</b>) is a 2D drawing with marked landmarks (in yellow) and calculated landmarks (in black), and part (<b>a</b>) is these 2D landmarks converted to 3D GIS (yellow and black) and real landmarks location in white.</p>
Full article ">Figure 9
<p>Illustrates the built areas (without the wall), measured in square meters, within the Old City of Jerusalem. These calculations are based on drawings and photos that capture the entire view of the built areas without big voids.</p>
Full article ">Figure 10
<p>Calculated GIS 3D built areas and landmarks in the Jerusalem (cartographic material) drawing (1782). The annotated landmarks are in yellow, the recommended landmarks are in black, the ground truth landmarks are in white, and the drawing built areas are in blue.</p>
Full article ">Figure A1
<p>The 40 noted features in and around the Old City of Jerusalem. True geotagged locations from which modern photographs were taken (squares); and calculated viewpoints of old artworks and modern photographs (circles) [<a href="#B23-ijgi-12-00500" class="html-bibr">23</a>].</p>
Full article ">Figure A2
<p>The Jerusalem cartographic drawing from 1782 [<a href="#B26-ijgi-12-00500" class="html-bibr">26</a>].</p>
Full article ">Figure A3
<p>Visual comparison of built areas in different drawings: (<b>a</b>) Trion’s drawing (1732). (<b>b</b>) Jerusalem (Cartographic Material) drawing (1782). (<b>c</b>) Henniker drawing (1823).</p>
Full article ">Figure A4
<p>Comprehensive comparison between the 2D drawings and 3D renderings of Trion and De Bruyn drawings. It highlights that the Trion drawing (1732) appears to be based on the De Bruyn drawing (1700), with noticeable similarities observed in their 3D representations. Part (<b>a</b>) shows De Bruyn’s drawing and part (<b>b</b>) shows Tirion’s drawing. Part (<b>c</b>) shows the building area in these two drawings.</p>
Full article ">Figure A5
<p>Three types of drawings/photos: cut, full-view, and panoramic.“Real” landmarks are in yellow, calculated landmarks are in black, and the distance between “real” and calculated landmarks is in green. (<b>a</b>) De Bruyn’s. (<b>b</b>) Trion’s. (<b>c</b>) Vue Générale de Jérusalem. (<b>d</b>) Vue du Mont des Oliviers. (<b>e</b>) Gotfryd, Bernard.</p>
Full article ">
21 pages, 7098 KiB  
Article
What Drives the Spatial Heterogeneity of Urban Leisure Activity Participation? A Multisource Big Data-Based Metrics in Nanjing, China
by Shaojun Liu, Xiawei Chen, Fengji Zhang, Yiyan Liu and Junlian Ge
ISPRS Int. J. Geo-Inf. 2023, 12(12), 499; https://doi.org/10.3390/ijgi12120499 - 12 Dec 2023
Cited by 2 | Viewed by 2328
Abstract
With the rapid pace of urbanization, enhancing the quality of life has become an urgent demand for the general public in both developed and developing countries. This study addresses the pressing need to understand the spatial distribution and underlying mechanisms of urban leisure [...] Read more.
With the rapid pace of urbanization, enhancing the quality of life has become an urgent demand for the general public in both developed and developing countries. This study addresses the pressing need to understand the spatial distribution and underlying mechanisms of urban leisure activity participation. To achieve this, we propose a novel methodological framework that integrates diverse big data sources, including mobile phone signaling data, urban geospatial data, and web-crawled data. By applying this framework to the urban area of Nanjing, our study reveals both the temporal and spatial patterns of urban leisure activity participation in the city. Notably, leisure activity participation is significantly higher on weekends, with distinctive daily peaks. Moreover, we identify spatial heterogeneity in leisure activity participation across the study area. Leveraging the OLS regression model, we design and quantify a comprehensive set of 12 internal and external indicators to explore the formation mechanisms of leisure participation for different leisure activity types. Our findings offer valuable guidance for urban planners and policymakers to optimize the allocation of resources, enhance urban street environments, and develop leisure resources in a rational and inclusive manner. Ultimately, this study contributes to the ongoing efforts to improve the quality of urban life and foster vibrant and sustainable cities. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) is the case study area of Nanjing, China, and (<b>b</b>) shows the distribution of ULRs in our study area. The study area includes one main urban area and three sub-districts of Jiangbei, Xianlin, and Dongshan.</p>
Full article ">Figure 2
<p>All collected points of BMSV images in the left figure and street map view image examples of four viewing angles.</p>
Full article ">Figure 3
<p>Framework for analysis of urban leisure participation patterns and their driving mechanisms.</p>
Full article ">Figure 4
<p>The example for the spatial association of ULRs and station buffers. (<b>a</b>) is the example stop regions and ULRs. (<b>b</b>) represents a stop region with only one ULR associated with it. (<b>c</b>,<b>d</b>) shows a stop region associated with multiple ULRs.</p>
Full article ">Figure 5
<p>Examples of the semantic segmentation results for BMSV images.</p>
Full article ">Figure 6
<p>The temporal pattern of LAP during April in the study area of Nanjing. (<b>a</b>) shows the hourly changes of LAP on weekdays and weekends; (<b>b</b>,<b>c</b>) shows the frequency distribution of leisure duration on weekdays and weekends, the curves are the cumulative percentage change.</p>
Full article ">Figure 7
<p>The spatial distribution of LAP during April in the study area of Nanjing.</p>
Full article ">Figure 8
<p>Histogram of standardized residuals.</p>
Full article ">Figure 9
<p>Streetscape images of two attractive cultural venues in Nanjing (front and rear view). (<b>a</b>) is the street scene near Nanjing Culture and Art Center and (<b>b</b>) is the street scene near Nanjing Museum.</p>
Full article ">
19 pages, 34010 KiB  
Article
Visibility-Based R-Tree Spatial Index for Consistent Visualization in Indoor and Outdoor Scenes
by Chengpeng Li, Xi Kuai, Biao He, Zhigang Zhao, Haojia Lin, Wei Zhu, Yu Liu and Renzhong Guo
ISPRS Int. J. Geo-Inf. 2023, 12(12), 498; https://doi.org/10.3390/ijgi12120498 - 12 Dec 2023
Cited by 1 | Viewed by 2124
Abstract
(1) Background: The smart city management system, with GIS technology as its core, is based on realistic visualization of multiple types of 3D model data syntheses. However, the efficiency barriers to achieving smooth and continuous visualization from outdoor scenes to small indoor scenes [...] Read more.
(1) Background: The smart city management system, with GIS technology as its core, is based on realistic visualization of multiple types of 3D model data syntheses. However, the efficiency barriers to achieving smooth and continuous visualization from outdoor scenes to small indoor scenes remain a challenge. (2) Methods: This paper uses the visibility prediction method to obtain potential visual sets at three levels—outdoor, indoor and outdoor connection, and indoor—and constructs an R-tree spatial index structure for organizing potential visual sets. By integrating these potential visible sets with spatial indexes, scene visualization can be carried out effectively. (3) Results: A near-reality indoor and outdoor scene was used for experimentation, resulting in stable 10% fluctuation visual frame rates around 90 FPS. (4) Conclusions: Spatial indexing methods that combine potential visible sets can effectively solve the continuity and stability problem of indoor and outdoor scene visualization in smart city management systems. Full article
Show Figures

Figure 1

Figure 1
<p>Viewpoint space and the visual world: (<b>a</b>) the spaces and entities; (<b>b</b>) the visual results of a single direction; and (<b>c</b>) the entire 3D visual world.</p>
Full article ">Figure 2
<p>The relationships among each component of VESI.</p>
Full article ">Figure 3
<p>The process of restructuring an indoor entity.</p>
Full article ">Figure 4
<p>PVS sampling strategies: (<b>a</b>) calculating visible building shells, (<b>b</b>) calculating visible wall elements for VS voxels and the other side, and (<b>c</b>) calculating visible objects in dense spaces with an omnidirectional view for one viewpoint.</p>
Full article ">Figure 5
<p>Multi-frame recovery to obtain visible objects.</p>
Full article ">Figure 6
<p>The organization of VESI: (<b>a</b>) the viewpoint space voxels that record the visible objects and (<b>b</b>) the R-tree for voxels and the associated visible objects.</p>
Full article ">Figure 7
<p>The pseudo-code for “Data Scheduling Using VESI”.</p>
Full article ">Figure 8
<p>The pseudo-code of “Find <span class="html-italic">Nextvoxel</span>”.</p>
Full article ">Figure 9
<p>A scheduling sample: (<b>a</b>) two VS voxels that record the visible objects with intersect parts and (<b>b</b>) incremental objects when the camera moves from Voxel5 to Voxel6.</p>
Full article ">Figure 10
<p>Continuous indoor and outdoor scenes: (<b>a</b>) a physical scene as a photo and (<b>b</b>) a digital scene as 3D models.</p>
Full article ">Figure 11
<p>The structure of the VESI file in Json format (show in Notepad++). (<b>a</b>) The list of visible sets in VESI. (<b>b</b>) The tree structure of VESI.</p>
Full article ">Figure 12
<p>The preprocessing speed of VESI: (<b>a</b>) the cost limited by the threshold of the number of triangles (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>H</mi> </mrow> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </semantics></math>) in a single object and (<b>b</b>) the R-tree for voxels and the associated visible objects.</p>
Full article ">Figure 13
<p>The visual effects test: (<b>a</b>) based on VESI, (<b>b</b>) based on the spatial index, and (<b>c</b>) based on PVS.</p>
Full article ">Figure 13 Cont.
<p>The visual effects test: (<b>a</b>) based on VESI, (<b>b</b>) based on the spatial index, and (<b>c</b>) based on PVS.</p>
Full article ">Figure 14
<p>Visual stability statistics: (<b>a</b>) based on VESI, (<b>b</b>) based on the spatial index, and (<b>c</b>) based on PVS.</p>
Full article ">
17 pages, 4062 KiB  
Article
Location Planning of Emergency Medical Facilities Using the p-Dispersed-Median Modeling Approach
by Changwha Oh, Yongwan Chun and Hyun Kim
ISPRS Int. J. Geo-Inf. 2023, 12(12), 497; https://doi.org/10.3390/ijgi12120497 - 12 Dec 2023
Viewed by 1926
Abstract
This research employs a spatial optimization approach customized for addressing equitable emergency medical facility location problems through the p-dispersed-median problem (p-DIME). The p-DIME integrates two conflicting classes of spatial optimization problems, dispersion and median problems, aiming to identify the [...] Read more.
This research employs a spatial optimization approach customized for addressing equitable emergency medical facility location problems through the p-dispersed-median problem (p-DIME). The p-DIME integrates two conflicting classes of spatial optimization problems, dispersion and median problems, aiming to identify the optimal locations for emergency medical facilities to achieve an equitable spatial distribution of emergency medical services (EMS) while effectively serving demand. To demonstrate the utility of the p-DIME model, we selected Gyeongsangbuk-do in South Korea, recognized as one of the most challenging areas for providing EMS to the elderly population (aged 65 and over). This challenge arises from the significant spatial disparity in the distribution of emergency medical facilities. The results of the model assessment gauge the spatial disparity of EMS, provide significantly enhanced solutions for a more equitable EMS distribution in terms of service coverage, and offer policy implications for future EMS location planning. In addition, to address the computational challenges posed by p-DIME’s inherent complexity, involving mixed-integer programming, this study introduces a solution technique through constraint formulations aimed at tightening the lower bounds of the problem’s solution space. The computational results confirm the effectiveness of this approach in ensuring reliable computational performance, with significant reductions in solution times, while still producing optimal solutions. Full article
Show Figures

Figure 1

Figure 1
<p>Research extent and the distribution of the candidate medical facilities for LEMIs. Note: green dots represent the location of current LEMIs.</p>
Full article ">Figure 2
<p>Population distribution of the research area. Note: the natural breaks method was used as a classification method to visualize the distribution of the registered population.</p>
Full article ">Figure 3
<p>Comparison of geographic coverage by the current 31 LEMIs (<b>a</b>) and the LEMIs by the <span class="html-italic">p</span>-DIME model when <span class="html-italic">p</span> = 31 (<b>b</b>).</p>
Full article ">Figure 4
<p>Change in two objective values of the <span class="html-italic">p</span>-DIME model with <span class="html-italic">p</span>. (<b>a</b>) Change in the objective function values on the maxisum dispersion terms (locations). (<b>b</b>) Change in the objective function values on the median terms (allocations).</p>
Full article ">Figure 5
<p>Comparison of solution times: (a) <span class="html-italic">p</span>-DIME model without APRIL (a) vs. (b) the <span class="html-italic">p</span>-DIME with APRIL constraint.</p>
Full article ">Figure 6
<p>Trade-off between dispersion and median of LEMIs in the <span class="html-italic">p</span>-DIME model (<span class="html-italic">p</span> = 31). Note: for a better visualization of Pareto curve, the population-weighted demand values are rescaled by 1/10,000 [<a href="#B45-ijgi-12-00497" class="html-bibr">45</a>].</p>
Full article ">
15 pages, 3162 KiB  
Article
Research on Traffic Accident Risk Prediction Method Based on Spatial and Visual Semantics
by Wenjing Li and Zihao Luo
ISPRS Int. J. Geo-Inf. 2023, 12(12), 496; https://doi.org/10.3390/ijgi12120496 - 11 Dec 2023
Cited by 2 | Viewed by 2657
Abstract
Predicting traffic accidents involves analyzing historical data, determining the relevant factors affecting the occurrence of traffic accidents, and predicting the likelihood of future traffic accidents. Most of the previous studies used statistical methods or single deep learning network model prediction methods while ignoring [...] Read more.
Predicting traffic accidents involves analyzing historical data, determining the relevant factors affecting the occurrence of traffic accidents, and predicting the likelihood of future traffic accidents. Most of the previous studies used statistical methods or single deep learning network model prediction methods while ignoring the visual effects of the city landscape on the drivers and the zero-inflation problem, resulting in poor prediction performance. Therefore, this paper constructs a city traffic accident risk prediction model that incorporates spatial and visual effects on drivers. The improved STGCN model is used in the model, a CNN and GRU replace the origin space–time convolution layer, two layers of a GCN are added to extract the city landscape similarity of different regions, and a BN layer is added to solve the gradient explosion problem. Finally, the features extracted from the time–space correlation module, the city landscape similarity module and the spatial correlation module are fused. The model is trained with the self-made Chicago dataset and compared with the existing network model. The comparison experiment proves that the prediction effect of the model in both the full time period and the high-frequency time period is better than that of the existing model. The ablation experiment proves that the city landscape similarity module added in this paper performs well in the high-frequency area. Full article
Show Figures

Figure 1

Figure 1
<p>The architecture of SVSNet.</p>
Full article ">Figure 2
<p>The road of Chicago.</p>
Full article ">Figure 3
<p>An example of semantics segmentation.</p>
Full article ">Figure 4
<p>An example of the same pixel type and different distribution. The number of pixels with the same color in (<b>a</b>,<b>b</b>) is identical. The difference between these two images lies in the distribution of colors. The color distribution in (<b>a</b>) is noticeably more complex than in (<b>b</b>). Therefore, even though the colors in the two images are the same, the perceived complexity by the human eye is different due to the distinct color distributions.</p>
Full article ">Figure 5
<p>Traffic scene complexity distribution map.</p>
Full article ">Figure 6
<p>An example of geographical correlations. The red area indicates the region where a traffic accident occurred at a certain moment, while the yellow area represents the possibility that vehicles, due to the occurrence of a traffic accident in the red area, may choose to detour through the yellow area. This diversion results in increased traffic flow in the yellow area, potentially impacting the probability of traffic accidents occurring in the yellow region.</p>
Full article ">Figure 7
<p>Traffic accident risk chart for four consecutive weeks.</p>
Full article ">Figure 8
<p>The prediction results.</p>
Full article ">
21 pages, 10820 KiB  
Article
An Aggregated Shape Similarity Index: A Case Study of Comparing the Footprints of OpenStreetMap and INSPIRE Buildings
by Renata Ďuračiová
ISPRS Int. J. Geo-Inf. 2023, 12(12), 495; https://doi.org/10.3390/ijgi12120495 - 9 Dec 2023
Viewed by 2098
Abstract
The mutual identification of spatial objects is a fundamental issue when updating geographic data with other data sets. Representations of spatial objects in different sources may not have the same identifiers, which would unambiguously assign them to each other. Intersections of spatial objects [...] Read more.
The mutual identification of spatial objects is a fundamental issue when updating geographic data with other data sets. Representations of spatial objects in different sources may not have the same identifiers, which would unambiguously assign them to each other. Intersections of spatial objects can be used for this purpose, but this does not allow for the detection of possible changes and their quantification. The aim of this paper is to propose a simple, applicable procedure for calculating the shape similarity measure, which should be able to efficiently identify different representations of spatial objects in two data sources, even though they may be changed or generalised. The main result is the aggregated index of shape similarity and instructions for its calculation and implementation. The shape similarity index is based on the calculation of the set similarity, the distance of the boundaries, and the differences in the area, perimeter, and number of the vertices of areal spatial objects. In the case study, the footprints of the building complexes in Dúbravka (part of the city of Bratislava, the capital of Slovakia) are compared using data from OpenStreetMap and INSPIRE (Infrastructure for Spatial Information in Europe) Buildings. A contribution to the quality check of the OpenStreetMap data is then a secondary result. The proposed method can be effectively used in the semi-automatic integration of heterogeneous data sources, updating the data source with other spatial data, or in their quality control. Full article
Show Figures

Figure 1

Figure 1
<p>An example of the representation of buildings in two data sources: (<b>a</b>) objects with a different level of detail (different number of vertices), (<b>b</b>) objects with a different number of polygons, (<b>c</b>) moved position of objects, (<b>d</b>–<b>f</b>) objects whose identity or possible change needs to be assessed individually.</p>
Full article ">Figure 2
<p>Case study area—Bratislava–Dúbravka (Slovakia, Europe), data sources: [<a href="#B35-ijgi-12-00495" class="html-bibr">35</a>,<a href="#B46-ijgi-12-00495" class="html-bibr">46</a>].</p>
Full article ">Figure 3
<p>A sample of the OSM and INSPIRE Buildings footprints with their similarity indices, data sources: [<a href="#B35-ijgi-12-00495" class="html-bibr">35</a>,<a href="#B46-ijgi-12-00495" class="html-bibr">46</a>].</p>
Full article ">Figure 4
<p>Example of visualisation of aggregated shape similarity indices of OSM and INSPIRE Buildings in QGIS.</p>
Full article ">Figure 5
<p>Example of the classification of building footprints into the categories of object similarity, data sources: [<a href="#B35-ijgi-12-00495" class="html-bibr">35</a>,<a href="#B46-ijgi-12-00495" class="html-bibr">46</a>].</p>
Full article ">Figure 6
<p>A sample of identical building footprints represented by different numbers of areal objects. (<b>a</b>) buildings that are very close to each other but not connected; (<b>b</b>) multiple buildings being connected into one, for example, by a walkway.</p>
Full article ">
21 pages, 5886 KiB  
Article
A Motion-Based Conceptual Space Model to Support 3D Evacuation Simulation in Indoor Environments
by Ruihang Xie, Sisi Zlatanova, Jinwoo (Brian) Lee and Mitko Aleksandrov
ISPRS Int. J. Geo-Inf. 2023, 12(12), 494; https://doi.org/10.3390/ijgi12120494 - 8 Dec 2023
Cited by 2 | Viewed by 2184
Abstract
Three-dimensional (3D) indoor models are a crucial component to simulate pedestrian evacuations realistically in indoor environments. However, existing 3D indoor models cannot fully represent realistic indoor environments to enable the simulation of 3D pedestrian motions in evacuations because spaces above/below some physical components [...] Read more.
Three-dimensional (3D) indoor models are a crucial component to simulate pedestrian evacuations realistically in indoor environments. However, existing 3D indoor models cannot fully represent realistic indoor environments to enable the simulation of 3D pedestrian motions in evacuations because spaces above/below some physical components (e.g., desks, chairs) have been largely overlooked. Thus, this paper introduces a conceptual space model to advance a space identification and classification scheme that can fully capture 3D pedestrian motions. This paper first proposes the definition and parameterisation of different 3D pedestrian motions. Then, the definition and specifications of three categories of space components are elaborated on based on the motions. Finally, a voxel-based approach is introduced to identify and classify the space components, which are demonstrated by an illustrative example. This work contributes to advancing 3D indoor modelling to enable a more realistic simulation of 3D pedestrian motions. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of 3D pedestrian motions: (1) 3D movements: (<b>a</b>) low crawling, (<b>b</b>) knee and hand crawling, (<b>c</b>) bent-over walking and (<b>d</b>) walking upright; and (2) and 3D actions: (<b>e</b>) stepping up/down, (<b>f</b>) jumping up/down and (<b>g</b>) climbing up/down.</p>
Full article ">Figure 2
<p>Three-dimensional cuboids for the four types of three-dimensional movements: (<b>a</b>) low crawling, (<b>b</b>) knee and hand crawling, (<b>c</b>) bent-over walking, (<b>d</b>) walking upright.</p>
Full article ">Figure 3
<p>(<b>a</b>) A populated room with static, movable and dynamic objects. (<b>b</b>) The room with space components: P-space (green), C-space (yellow) and N-space (red).</p>
Full article ">Figure 4
<p>Spaces are categorised into three types according to the definitions of P-, C- and N-spaces.</p>
Full article ">Figure 5
<p>Four types of C-spaces. (<b>a</b>) C<sub>l</sub>-spaces. (<b>b</b>) C<sub>k</sub>-spaces. (<b>c</b>) C<sub>b</sub>-spaces. (<b>d</b>) C<sub>u</sub>-spaces.</p>
Full article ">Figure 6
<p>A C<sub>l</sub>-space under a table. (<b>a</b>) <span class="html-italic">h<sub>l</sub></span> and <span class="html-italic">l<sub>l</sub></span> of the 3D cuboid for low crawling fit within the space. (<b>b</b>) Adjacent P-spaces accommodate a portion of the <span class="html-italic">w<sub>l</sub></span>.</p>
Full article ">Figure 7
<p>Schematic diagram of a parameter for the clearance <span class="html-italic">h<sub>n</sub></span>.</p>
Full article ">Figure 8
<p>Flowchart of the space identification and classification based on the conceptual space model.</p>
Full article ">Figure 9
<p>Location of the four 3D cuboids above an occupied voxel in the X, Y and Z dimensions.</p>
Full article ">Figure 10
<p>BIM IFC model of a furnished room.</p>
Full article ">Figure 11
<p>Identification and classification process of space components.</p>
Full article ">Figure 12
<p>Three-dimensional voxel model of the furnished room with identified and classified space components: (<b>a</b>) P-spaces (green voxels), (<b>b</b>) C<sub>u</sub>-spaces (blue voxels), C<sub>k</sub>-spaces (yellow voxels), C<sub>l</sub>-spaces (cyan voxels) and N-spaces (red voxels).</p>
Full article ">Figure 13
<p>Diagram of evacuation paths in the furniture room. (<b>a</b>) A 2D path. (<b>b</b>) A 3D path using P-spaces and C<sub>u</sub>-spaces. (<b>c</b>) A 3D path using P-spaces, C<sub>k</sub>-spaces and C<sub>l</sub>-spaces. N-spaces (red voxels), C-spaces (yellow and blue voxels) and P-spaces (green voxels).</p>
Full article ">
23 pages, 14593 KiB  
Article
Risk Assessment of Landslide Collapse Disasters along National Highways Based on Information Quantity and Random Forest Coupling Methods: A Case Study of the G331 National Highway
by Zuoquan Nie, Qiuling Lang, Yichen Zhang, Jiquan Zhang, Yanan Chen and Zengkai Pan
ISPRS Int. J. Geo-Inf. 2023, 12(12), 493; https://doi.org/10.3390/ijgi12120493 - 6 Dec 2023
Cited by 2 | Viewed by 2212
Abstract
Based on the data from two field surveys in 2015 and 2022, this paper calculates the weight of values using the entropy weight method and the variation coefficient method, and evaluates risk using the information quantity method. The information quantities of four levels [...] Read more.
Based on the data from two field surveys in 2015 and 2022, this paper calculates the weight of values using the entropy weight method and the variation coefficient method, and evaluates risk using the information quantity method. The information quantities of four levels of criteria (hazards, exposure, vulnerability, emergency responses, and capability of recovery) were extracted and inputted into a random forest model. After optimizing the hyperparameters of the random forest using GridSearchCV, the risk assessment was performed again. Finally, the accuracy of the two evaluation results was verified using an ROC curve, and the model with the higher AUC value was selected to create a risk map. Compared with previous studies, this paper considers the factors of emergency responses and recovery capability, which makes the risk assessment more comprehensive. Our findings show that the evaluation results based on the coupling model are more accurate than the evaluation results of the information method, as the coupling model had an AUC value of 0.9329. After considering the indices of emergency responses and capability of recovery, the risk level of the highest-risk area in the study area decreased. Full article
(This article belongs to the Topic Geotechnics for Hazard Mitigation)
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area.</p>
Full article ">Figure 2
<p>Workflow of the study.</p>
Full article ">Figure 3
<p>Spatial distribution of the hazard indicators. ((<b>a</b>) Slope, (<b>b</b>) aspect, (<b>c</b>) curvature, (<b>d</b>) NDVI, (<b>e</b>) mean annual precipitation, (<b>f</b>) distance from fault, (<b>g</b>) lithology, (<b>h</b>) distance from road).</p>
Full article ">Figure 3 Cont.
<p>Spatial distribution of the hazard indicators. ((<b>a</b>) Slope, (<b>b</b>) aspect, (<b>c</b>) curvature, (<b>d</b>) NDVI, (<b>e</b>) mean annual precipitation, (<b>f</b>) distance from fault, (<b>g</b>) lithology, (<b>h</b>) distance from road).</p>
Full article ">Figure 4
<p>Spatial distribution of the exposure indices. ((<b>a</b>) Population density, (<b>b</b>) road density, (<b>c</b>) building density).</p>
Full article ">Figure 5
<p>Spatial distribution of the vulnerability indices. ((<b>a</b>) Age of the structure, (<b>b</b>) road classification, (<b>c</b>) building types).</p>
Full article ">Figure 6
<p>Spatial distribution of the emergency responses and recovery capability indices. ((<b>a</b>) GDP, (<b>b</b>) educational status, (<b>c</b>) number of medical staff).</p>
Full article ">Figure 7
<p>Drawing of each criterion layer. ((<b>a</b>) Hazard map, (<b>b</b>) exposure map, (<b>c</b>) vulnerability map, (<b>d</b>) emergency responses and recovery capability map).</p>
Full article ">Figure 7 Cont.
<p>Drawing of each criterion layer. ((<b>a</b>) Hazard map, (<b>b</b>) exposure map, (<b>c</b>) vulnerability map, (<b>d</b>) emergency responses and recovery capability map).</p>
Full article ">Figure 8
<p>ROC curve for the information quantity method.</p>
Full article ">Figure 9
<p>ROC curve for the random forest model. ((<b>a</b>) Initial model, (<b>b</b>) optimizing n_estimators, (<b>c</b>) final model).</p>
Full article ">Figure 10
<p>Risk map.</p>
Full article ">Figure 11
<p>Risk map excluding emergency responses and recovery capability.</p>
Full article ">
24 pages, 6275 KiB  
Article
Exploring the Spatial Relationship between Street Crime Events and the Distribution of Urban Greenspace: The Case of Porto, Portugal
by Miguel Saraiva and Bárbara Teixeira
ISPRS Int. J. Geo-Inf. 2023, 12(12), 492; https://doi.org/10.3390/ijgi12120492 - 6 Dec 2023
Cited by 2 | Viewed by 3445
Abstract
In post-pandemic, climate-changing societies, the presence of urban greenspace assumes paramount functions, at the same time that socio-economic crises and shocks augment vulnerabilities and insecurities. The recent literature on environmental criminology argues that the geography of crime is not random, and that the [...] Read more.
In post-pandemic, climate-changing societies, the presence of urban greenspace assumes paramount functions, at the same time that socio-economic crises and shocks augment vulnerabilities and insecurities. The recent literature on environmental criminology argues that the geography of crime is not random, and that the presence of greenery, due to its impact on well-being and the environment, can have positive associations with feeling safe; although the opposite effect can occur if spaces are not properly designed or maintained. In this paper, the case study of Porto, Portugal, is presented; one of the municipalities with higher crime rates, that also pledged to double the available greenspace in the near future. As a way to support decision-making, the aim of this study was to present an overall exploratory diagnosis of how street crime patterns, of different typologies, spatially co-exist with greenspaces. Using a 10-year street crime dataset at the segment level, descriptive quantitative methods with the support of GIS have been applied to plot crime’s spatial distribution over time, as well as the walking accessibility to greenspaces. The results confirm crime’s geographical non-randomness, with distinct categories occupying specific locations, even though there was a consistently proportional distribution in the different distance bands. On the contrary, the cumulative effect of the proximity to greenspaces was variable. Almost half of the city’s street crimes (46%) were within a 5 min walking distance of greenspaces, but they were much closer to smaller inner-city urban gardens, with higher densities of street crimes (hot spots), than to larger municipal parks, where lower densities (cold spots) were seen. Full article
Show Figures

Figure 1

Figure 1
<p>Porto Municipality, Portugal.</p>
Full article ">Figure 2
<p>The city of Porto and the locations of the analysed green spaces, by typology.</p>
Full article ">Figure 3
<p>Walking distance (5, 10, 15 min) to the closest greenspaces of Porto.</p>
Full article ">Figure 4
<p>Kernel Density Estimation, considering a 50 m cell size and crime at street segment, for all reported street crimes between 2009 and 2018 in Porto (source: own; based on untreated raw data from the Public Safety Police).</p>
Full article ">Figure 5
<p>Kernel Density Estimation, with 50 m cell size and crime at the street segment level, for three of the thirteen categories analysed: (<b>a</b>) pickpocketing, (<b>b</b>) thefts in motorized vehicles, and (<b>c</b>) thefts in the street, between 2009 and 2018 in Porto (source: own; based on untreated raw data from the Public Safety Police).</p>
Full article ">Figure 6
<p>Emerging hot spot analysis for reported street crimes between 2009 and 2018 with a 1-year time step (source: own; based on untreated raw data from the Public Safety Police).</p>
Full article ">Figure 7
<p>Distribution (in percentages) of the total amount of each typology of registered crime, by walking distance to overall greenspaces, gardens, and parks.</p>
Full article ">Figure 8
<p>Distribution (in percentages) of how total street crimes, by walking distance, are divided into crime typologies, for overall greenspaces, gardens, and parks.</p>
Full article ">
19 pages, 68820 KiB  
Article
Multiscale Analysis of Spatial Accessibility to Acute Hospitals in Carinthia, Austria
by Changzhen Wang, Michael Leitner and Gernot Paulus
ISPRS Int. J. Geo-Inf. 2023, 12(12), 491; https://doi.org/10.3390/ijgi12120491 - 6 Dec 2023
Cited by 2 | Viewed by 2032
Abstract
Health care accessibility studies are well established in the US but lacking in Austria, even though both experience high costs and have hospital care as the largest contributor to health care spending. This study aims to examine multiscale spatial accessibility to acute hospitals [...] Read more.
Health care accessibility studies are well established in the US but lacking in Austria, even though both experience high costs and have hospital care as the largest contributor to health care spending. This study aims to examine multiscale spatial accessibility to acute hospitals in Carinthia, Austria. Using the most recent data at census block and 250 meter grid levels, we refine proximity and generalized two-step floating catchment area (G2SFCA) methods while accounting for the modifiable areal unit problem (MAUP) and edge effects. For census blocks and 250 meter grids, the mean travel times to the nearest acute hospitals are 16 and 21 min, respectively, covering 58.8% and 76.2% of the population, which, however, increases to 25 and 31 min to the three nearest hospitals with similar populations. People bypassing the nearest hospital to seek hospitals at a longer distance, termed “bypass behavior”, is more influential, as 20% more of the population living in mountainous or rural areas need to travel 30 min longer. The G2SFCA method with a more pronounced distance decay results in a more decentralized polycentric structure of accessibility and identifies poorer access areas. While urban advantage is most evident in Klagenfurt and Villach, not all areas near hospitals enjoy the highest accessibility. A combination of the proximity and G2SFCA methods identifies less accessible areas. The MAUP overestimates accessibility at a coarse level and in less populous areas. Edge effects occur at the border when using proximity only, but they are more sensitive when considering bypass behavior or a weak distance decay effect. This study contributes to our understanding of acute hospitals’ accessibility in Carinthia and highlights the need to improve low-accessible areas in addition to universal health coverage. Cautions need to be exercised when using different geographic units or considering edge effects for health care planning and management. Full article
Show Figures

Figure 1

Figure 1
<p>The study area: The Austrian province of Carinthia located in southern Austria.</p>
Full article ">Figure 2
<p>(<b>a</b>) Population density of census blocks and acute hospital beds in Carinthia; (<b>b</b>) 250 meter grid population and acute hospital beds in Carinthia.</p>
Full article ">Figure 2 Cont.
<p>(<b>a</b>) Population density of census blocks and acute hospital beds in Carinthia; (<b>b</b>) 250 meter grid population and acute hospital beds in Carinthia.</p>
Full article ">Figure 3
<p>Boxplots of travel times to acute hospitals across census blocks and 250 meter grids. The horizontal red dash line represents the mean value of travel times in each category.</p>
Full article ">Figure 4
<p>Travel times to the nearest acute hospital across (<b>a</b>) census blocks and (<b>b</b>) 250 meter grids; and average travel times to the three nearest acute hospitals across (<b>c</b>) census blocks and (<b>d</b>) 250 meter grids in Carinthia.</p>
Full article ">Figure 5
<p>(<b>a</b>) Differences in travel time to the nearest acute hospital; and (<b>b</b>) differences in average travel time to the three nearest acute hospitals between grids and blocks.</p>
Full article ">Figure 6
<p>Scatter plots and maps of block- and grid-based accessibilities using (<b>a</b>,<b>b</b>) β = 1; (<b>c</b>,<b>d</b>) β = 1.1; (<b>e</b>,<b>f</b>) β = 1.2; (<b>g</b>,<b>h</b>) β = 1.3; (<b>i</b>,<b>j</b>) β = 1.4; (<b>k</b>,<b>l</b>) β = 1.5; (<b>m</b>,<b>n</b>) β = 1.6 in the G2SFCA method. Note: all district names from Carinthia can be found in <a href="#ijgi-12-00491-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 6 Cont.
<p>Scatter plots and maps of block- and grid-based accessibilities using (<b>a</b>,<b>b</b>) β = 1; (<b>c</b>,<b>d</b>) β = 1.1; (<b>e</b>,<b>f</b>) β = 1.2; (<b>g</b>,<b>h</b>) β = 1.3; (<b>i</b>,<b>j</b>) β = 1.4; (<b>k</b>,<b>l</b>) β = 1.5; (<b>m</b>,<b>n</b>) β = 1.6 in the G2SFCA method. Note: all district names from Carinthia can be found in <a href="#ijgi-12-00491-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 6 Cont.
<p>Scatter plots and maps of block- and grid-based accessibilities using (<b>a</b>,<b>b</b>) β = 1; (<b>c</b>,<b>d</b>) β = 1.1; (<b>e</b>,<b>f</b>) β = 1.2; (<b>g</b>,<b>h</b>) β = 1.3; (<b>i</b>,<b>j</b>) β = 1.4; (<b>k</b>,<b>l</b>) β = 1.5; (<b>m</b>,<b>n</b>) β = 1.6 in the G2SFCA method. Note: all district names from Carinthia can be found in <a href="#ijgi-12-00491-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 7
<p>(<b>a</b>,<b>b</b>) Edge effect of travel time across census blocks and grids in Carinthia; (<b>c</b>,<b>d</b>) edge effect of census block-based and block-based accessibility by G2SFCA using β = 1; and (<b>e</b>,<b>f</b>) edge effect of census block-based and block-based accessibility by G2SFCA using β = 1.6 in Carinthia.</p>
Full article ">
18 pages, 7077 KiB  
Article
On the Definition of Standard Parallels in Map Projections
by Miljenko Lapaine
ISPRS Int. J. Geo-Inf. 2023, 12(12), 490; https://doi.org/10.3390/ijgi12120490 - 6 Dec 2023
Cited by 3 | Viewed by 2367
Abstract
The article belongs to the field of theoretical research on map projections. It is observed that there is no unique and generally accepted definition of standard parallels in the cartographic literature. For some authors, a standard line is a line along which there [...] Read more.
The article belongs to the field of theoretical research on map projections. It is observed that there is no unique and generally accepted definition of standard parallels in the cartographic literature. For some authors, a standard line is a line along which there is no distortion, and for others, it is a line along which there is no distortion of length. At the same time, it is forgotten that the length distortions at any point generally change and depend on the direction. The main goal of this article is very simple: the sentence “linear deformation is zero in all directions” is expressed using a mathematical formula. Besides that, the paper introduces equidistance in a broader sense. This is a novelty in the theory of map projections. Equidistance is defined at a point, along a line and in an area, especially in the direction of the parallels and especially in the direction of the meridian. This enables an unambiguous definition of standard parallels. Theoretical considerations are illustrated with examples of cylindrical projections. The practical value of the proposed approach is manifested in the possibility of a better understanding of the distribution of distortions in any map projection used. Full article
Show Figures

Figure 1

Figure 1
<p>Tissot’s indicatrix when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> at a point. At that point, the linear distortion is zero in the direction of the meridian. The direction of the meridian image is drawn with a part of a vertical line.</p>
Full article ">Figure 2
<p>Tissot’s indicatrix when <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> at a point. At that point, the linear distortion is zero in the direction of the parallel. The direction of the parallel image is drawn with a part of a horizontal line.</p>
Full article ">Figure 3
<p>Tissot’s indicatrix is a unit circle when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> at a point. We say that if a point was mapped equidistantly in all directions, the point is <span class="html-italic">standard</span>. At that point, the linear distortion is equal to zero in all directions.</p>
Full article ">Figure 4
<p><math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <mfenced separators="|"> <mrow> <mi>π</mi> <mi>φ</mi> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> </mrow> </mfenced> </mrow> <mrow> <mn>2</mn> <mi>π</mi> </mrow> </mfrac> </mrow> </semantics></math>. Horizontal axis in degrees.</p>
Full article ">Figure 5
<p><math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>y</mi> </mrow> <mrow> <mi>d</mi> <mi>φ</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <mfenced separators="|"> <mrow> <mi>π</mi> <mo>+</mo> <mn>2</mn> <mi>φ</mi> </mrow> </mfenced> </mrow> <mrow> <mn>2</mn> <mi>π</mi> </mrow> </mfrac> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mrow> <mrow> <mn>2</mn> <mi mathvariant="normal">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>φ</mi> </mrow> </mrow> </mrow> </mfrac> </mrow> </semantics></math> (red). Horizontal axis in degrees.</p>
Full article ">Figure 6
<p>World map in cylindrical projection <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mrow> <mrow> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mn>2</mn> </mrow> </mfrac> </mrow> <mo>⁡</mo> <mrow> <mo>·</mo> </mrow> </mrow> <mi>λ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <mfenced separators="|"> <mrow> <mi>π</mi> <mi>φ</mi> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> </mrow> </mfenced> </mrow> <mrow> <mn>2</mn> <mi>π</mi> </mrow> </mfrac> </mrow> </semantics></math>. Equidistantly mapped parallels −30° and 30°. Standard parallel −30°.</p>
Full article ">Figure 7
<p><math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <mfenced separators="|"> <mrow> <mi>π</mi> <mi>φ</mi> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> </mrow> </mfenced> </mrow> <mrow> <mn>5</mn> <mi>π</mi> </mrow> </mfrac> </mrow> </semantics></math>. Horizontal axis in degrees.</p>
Full article ">Figure 8
<p><math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>y</mi> </mrow> <mrow> <mi>d</mi> <mi>φ</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <mfenced separators="|"> <mrow> <mi>π</mi> <mo>+</mo> <mn>2</mn> <mi>φ</mi> </mrow> </mfenced> </mrow> <mrow> <mn>5</mn> <mi>π</mi> </mrow> </mfrac> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mrow> <mrow> <mn>2</mn> <mi mathvariant="normal">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>φ</mi> </mrow> </mrow> </mrow> </mfrac> </mrow> </semantics></math> (red). Horizontal axis in degrees.</p>
Full article ">Figure 9
<p>World map in cylindrical projection <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mrow> <mrow> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mn>2</mn> </mrow> </mfrac> </mrow> <mo>⁡</mo> <mrow> <mo>·</mo> </mrow> </mrow> <mi>λ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mn>3</mn> <mfenced separators="|"> <mrow> <mi>π</mi> <mi>φ</mi> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> </mrow> </mfenced> </mrow> <mrow> <mn>5</mn> <mi>π</mi> </mrow> </mfrac> </mrow> </semantics></math>. Equidistantly mapped parallels −30° and 30°. No standard parallels.</p>
Full article ">Figure 10
<p><math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> </mrow> <mrow> <mn>5</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msup> <mo>+</mo> <mfrac> <mrow> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msup> <mo>+</mo> <mi>e</mi> <mi>φ</mi> </mrow> </semantics></math>; <span class="html-italic">a</span>, <span class="html-italic">c</span> and <span class="html-italic">e</span> from (35).</p>
Full article ">Figure 11
<p><math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>y</mi> </mrow> <mrow> <mi>d</mi> <mi>φ</mi> </mrow> </mfrac> <mo>=</mo> <mi>a</mi> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msup> <mo>+</mo> <mi>c</mi> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mo>+</mo> <mi>e</mi> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mrow> <mrow> <mn>2</mn> <mi mathvariant="normal">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>φ</mi> </mrow> </mrow> </mrow> </mfrac> </mrow> </semantics></math> (red); <span class="html-italic">a</span>, <span class="html-italic">c</span> and <span class="html-italic">e</span> from (35).</p>
Full article ">Figure 12
<p>Map of the world in cylindrical projection <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mrow> <mrow> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mn>2</mn> </mrow> </mfrac> </mrow> <mo>⁡</mo> <mrow> <mo>·</mo> </mrow> </mrow> <mi>λ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> </mrow> <mrow> <mn>5</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msup> <mo>+</mo> <mfrac> <mrow> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msup> <mo>+</mo> <mi>e</mi> <mi>φ</mi> </mrow> </semantics></math>; <span class="html-italic">a</span>, <span class="html-italic">c</span> and <span class="html-italic">e</span> from (35). Two standard parallels <math display="inline"><semantics> <mrow> <mfenced separators="|"> <mrow> <mo>±</mo> <mn>30</mn> <mo>°</mo> </mrow> </mfenced> </mrow> </semantics></math>.</p>
Full article ">Figure 13
<p><math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> </mrow> <mrow> <mn>5</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msup> <mo>+</mo> <mfrac> <mrow> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msup> <mo>+</mo> <mi>e</mi> <mi>φ</mi> </mrow> </semantics></math>; a, c and e from (36).</p>
Full article ">Figure 14
<p><math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>y</mi> </mrow> <mrow> <mi>d</mi> <mi>φ</mi> </mrow> </mfrac> <mo>=</mo> <mi>a</mi> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msup> <mo>+</mo> <mi>c</mi> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mo>+</mo> <mi>e</mi> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mrow> <mrow> <mn>2</mn> <mi mathvariant="normal">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>φ</mi> </mrow> </mrow> </mrow> </mfrac> </mrow> </semantics></math> (red); <span class="html-italic">a</span>, <span class="html-italic">c</span> and <span class="html-italic">e</span> from (36).</p>
Full article ">Figure 15
<p>The world map in cylindrical projection <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mrow> <mrow> <mfrac> <mrow> <msqrt> <mn>3</mn> </msqrt> </mrow> <mrow> <mn>2</mn> </mrow> </mfrac> </mrow> <mo>⁡</mo> <mrow> <mo>·</mo> </mrow> </mrow> <mi>λ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> </mrow> <mrow> <mn>5</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msup> <mo>+</mo> <mfrac> <mrow> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </mfrac> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msup> <mo>+</mo> <mi>e</mi> <mi>φ</mi> </mrow> </semantics></math>; <span class="html-italic">a</span>, <span class="html-italic">c</span> and <span class="html-italic">e</span> from (36). Two parallels are mapped equidistantly in the direction of parallels and the other two equidistantly in the direction of meridians. But there are no standard parallels.</p>
Full article ">Figure 16
<p><math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mi>φ</mi> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msub> </mrow> </mfenced> </mrow> </mfenced> </mrow> </mfenced> </mrow> </semantics></math> for the Patterson projection.</p>
Full article ">Figure 17
<p><math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>y</mi> </mrow> <mrow> <mi>d</mi> <mi>φ</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mfenced separators="|"> <mrow> <msub> <mrow> <mn>5</mn> <mi>c</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mfenced separators="|"> <mrow> <msub> <mrow> <mn>7</mn> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mn>9</mn> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msub> </mrow> </mfenced> </mrow> </mfenced> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> </mrow> <mrow> <mrow> <mrow> <mi mathvariant="normal">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>φ</mi> </mrow> </mrow> </mrow> </mfrac> </mrow> </semantics></math> (red) for the Patterson projection. There are no standard parallels because <math display="inline"><semantics> <mrow> <mi>h</mi> <mfenced separators="|"> <mrow> <mn>0</mn> </mrow> </mfenced> <mo>=</mo> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>&gt;</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 18
<p>The map of the world in the Patterson projection: <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mi>λ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mi>φ</mi> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <msub> <mrow> <mi>c</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msub> </mrow> </mfenced> </mrow> </mfenced> </mrow> </mfenced> </mrow> </semantics></math>. There are no standard parallels.</p>
Full article ">Figure 19
<p><math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mi>φ</mi> </mrow> </semantics></math> for the equidistant cylindrical projection.</p>
Full article ">Figure 20
<p><math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>y</mi> </mrow> <mrow> <mi>d</mi> <mi>φ</mi> </mrow> </mfrac> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> </mrow> <mrow> <mrow> <mrow> <mi mathvariant="normal">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>φ</mi> </mrow> </mrow> </mrow> </mfrac> </mrow> </semantics></math> (red) for the equidistant cylindrical projection. The equator is the standard parallel.</p>
Full article ">Figure 21
<p>World map in the equidistant cylindrical projection <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mi>λ</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>=</mo> <mi>φ</mi> </mrow> </semantics></math>. The standard parallel is the equator.</p>
Full article ">
19 pages, 7724 KiB  
Article
High-Temporal-Resolution Prediction of Malaria Transmission Risk in South Kivu, Democratic Republic of the Congo, Based on Multi-Criteria Evaluation Using Geospatial Data
by Ryunosuke Komura and Masayuki Matsuoka
ISPRS Int. J. Geo-Inf. 2023, 12(12), 489; https://doi.org/10.3390/ijgi12120489 - 6 Dec 2023
Cited by 1 | Viewed by 2020
Abstract
Malaria is a major public health concern, and accurate mapping of malaria risk is essential to effectively managing the disease. However, current models are unable to predict malaria risk with high temporal and spatial resolution. This study describes a climate-based model that can [...] Read more.
Malaria is a major public health concern, and accurate mapping of malaria risk is essential to effectively managing the disease. However, current models are unable to predict malaria risk with high temporal and spatial resolution. This study describes a climate-based model that can predict malaria risk in South Kivu, Democratic Republic of the Congo, daily at a resolution of 2 km using meteorological (relative humidity, precipitation, wind speed, and temperature) and elevation data. We used the multi-criteria evaluation technique to develop the model. For the weighting of factors, we employed the analytical hierarchy process and linear regression techniques to compare expert knowledge-driven and mathematical methods. Using climate data from the prior 2 weeks, the model successfully mapped regions with high malaria case numbers, enabling accurate prediction of high-risk regions. This research may contribute to the development of a sustainable malaria risk forecasting system, which has been a longstanding challenge. Overall, this study provides insights into model development to predict malaria risk with high temporal and spatial resolution, supporting malaria control and management efforts. Full article
(This article belongs to the Topic Spatial Epidemiology and GeoInformatics)
Show Figures

Figure 1

Figure 1
<p>Map of Democratic Republic of the Congo (<b>a</b>), and South Kivu (<b>b</b>). Reprinted with permission from Ref. [<a href="#B21-ijgi-12-00489" class="html-bibr">21</a>]. Copyright 2023 Google.</p>
Full article ">Figure 2
<p>Research workflow: This figure illustrates the sequential steps of the research process. The orange segment represents the preprocessing phase, encompassing data preparation. The red segmented phase involves the computation of malaria case numbers. The blue segment represents that data from the period of 1 January 2018, to 31 December 2020, was used as the training period, while the period from 1 January 2021, to 31 December 2021, was selected as the validation period. Subsequently, the green segment calculates malaria risks. The purple phase denotes the application of the LR method, while the yellow phase corresponds to the AHP method, responsible for factor weighting. Finally, the white phase involves a comparative analysis of the two methods.</p>
Full article ">Figure 3
<p>Relationships of 2-week average values of each factor with malaria risk over the period of 1 January 2018–31 December 2020. (<b>a</b>) Elevation (<b>b</b>) Humidity (<b>c</b>) Precipitaion (<b>d</b>) Temperature and (<b>e</b>) Wind speed.</p>
Full article ">Figure 4
<p>Relationship of malaria case number with the normalized factor value for the training period (left axis) and representation with the f function (right axis). (<b>a</b>) Normalized elevation (<b>b</b>) Normalized humidity (<b>c</b>) Normalized precipitaion (<b>d</b>) Normalized temperature (<b>e</b>) Normalized wind speed.</p>
Full article ">Figure 5
<p>Map of population (colored dots indicate locations with more than 0 people within a 20 m × 20 m grid) and elevation.</p>
Full article ">Figure 6
<p>Relationship of malaria risk (<span class="html-italic">S</span>) with malaria cases for AHP method for the training period (<b>a</b>), LR method for the training period (<b>b</b>), AHP method for the validation period (<b>c</b>), and LR method for the validation period (<b>d</b>). The training period was from 1 January 2018 to 31 December 2020, and the validation period was from 1 January 2021 to 31 December 2021.</p>
Full article ">Figure 7
<p>Relation of every 0.1 increment of malaria risk (<span class="html-italic">S</span>) and malaria cases of top 1–5% from 1 January 2018–31 December 2020, for AHP (<b>a</b>) and LR (<b>b</b>).</p>
Full article ">Figure 8
<p>Relationships of the top 1%, 5%, 10%, and 20% of malaria cases with malaria risk, shown for the training period (1 January 2018–31 December 2020) in the top row and validation period (1 January 2021–31 December 2021) in the bottom row with the AHP method in the left column and LR method in the right column. (<b>a</b>) AHP method for the training period (2018–2020), (<b>b</b>) LR method for the training period (2018–2020) (<b>c</b>) AHP method for the validation period (2021) (<b>d</b>) LR method for the validation period (2021).</p>
Full article ">Figure 9
<p>Comparison of malaria risk map and malaria case map. Observation dates are 1 April 2021 (rainy season) for the top panels and 21 June 2021 (dry season) for the bottom panels. To improve readability of the malaria maps, the maximum point is set to the 98th percentile of the 21 June 2021 (dry season) data. (<b>a</b>) Risk map obtained using AHP on 1 April 2021 (rainy season), (<b>b</b>) Risk map obtained using LR on 1 April 2021 (rainy season), (<b>c</b>) Malaria case map on 1 April 2021 (rainy season), (<b>d</b>) Risk map obtained using AHP on 21 June 2021 (dry season), (<b>e</b>) Risk map obtained using LR on 21 June 2021 (dry season), and (<b>f</b>) Malaria case map on 21 June 2021 (dry season).</p>
Full article ">Figure 10
<p>Zoomed map of malaria risk predicted and malaria cases observed at populated area on 1 April 2021 (rainy season) and 21 June 2021 (dry season). For the malaria map due to its visibility, the maximum point is set to the 98th percentile of the 21 June 2021 (dry season) data. (<b>a</b>) Risk map obtained using AHP on 1 April 2021 (rainy season), (<b>b</b>) Risk map obtained using LR on 1 April 2021 (rainy season) (<b>c</b>) Malaria case map on 1 April 2021 (rainy season) (<b>d</b>) Risk map obtained using AHP on 21 June 2021 (dry season) (<b>e</b>) Risk map obtained using LR on 21 June 2021 (dry season) (<b>f</b>) Malaria case map on 21 June 2021 (dry season).</p>
Full article ">
25 pages, 5028 KiB  
Article
Applying Dynamic Human Activity to Disentangle Property Crime Patterns in London during the Pandemic: An Empirical Analysis Using Geo-Tagged Big Data
by Tongxin Chen, Kate Bowers and Tao Cheng
ISPRS Int. J. Geo-Inf. 2023, 12(12), 488; https://doi.org/10.3390/ijgi12120488 - 6 Dec 2023
Viewed by 2390
Abstract
This study aimed to evaluate the relationships between different groups of explanatory variables (i.e., dynamic human activity variables, static variables of social disorganisation and crime generators, and combinations of both sets of variables) and property crime patterns across neighbourhood areas of London during [...] Read more.
This study aimed to evaluate the relationships between different groups of explanatory variables (i.e., dynamic human activity variables, static variables of social disorganisation and crime generators, and combinations of both sets of variables) and property crime patterns across neighbourhood areas of London during the pandemic (from 2020 to 2021). Using the dynamic human activity variables sensed from mobile phone GPS big data sets, three types of ‘Least Absolute Shrinkage and Selection Operator’ (LASSO) regression models (i.e., static, dynamic, and static and dynamic) differentiated into explanatory variable groups were developed for seven types of property crime. Then, the geographically weighted regression (GWR) model was used to reveal the spatial associations between distinct explanatory variables and the specific type of crime. The findings demonstrated that human activity dynamics impose a substantially stronger influence on specific types of property crimes than other static variables. In terms of crime type, theft obtained particularly high relationships with dynamic human activity compared to other property crimes. Further analysis revealed important nuances in the spatial associations between property crimes and human activity across different contexts during the pandemic. The result provides support for crime risk prediction that considers the impact of dynamic human activity variables and their varying influences in distinct situations. Full article
Show Figures

Figure 1

Figure 1
<p>Footfall-generation process. The stay detection retrieves stays from different users’ raw mobile phone GPS trajectories. Next, the footfalls in geospatial units are aggregated by the detected stays.</p>
Full article ">Figure 2
<p>All property crime and footfall monthly change (per LSOA) from 2020 to 2021.</p>
Full article ">Figure 3
<p>Standardised local coefficient values of explanatory variables in property crime (S + D) models. Each model’s training set is 19 months of data from January 2020 to July 2021.</p>
Full article ">Figure 4
<p>Standardised local coefficient values of explanatory variables in ‘Short-term theft from a person (S + D) models’. Each model’s training set is one month from January 2020, April 2020, August 2020, and November 2020, respectively.</p>
Full article ">Figure 5
<p>Maps of the spatial association (local coefficient values) between crime rates of theft from a person and the selected human activity variables in fit GWRs for February 2020, April 2020, August 2020, and November 2020. The bandwidths of the four GWRs are 14, 80, 23, and 34, respectively. The global <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> of the four GWRs are 0.92, 0.36, 0.82, and 0.74, respectively.</p>
Full article ">Figure A1
<p>The map of the residential population in London LSOAs.</p>
Full article ">Figure A2
<p>The relationship (<math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math>) between the ten types of human activity variables and resident population in LSOAs. The R-squared values were calculated from the fit OLS models (without constant) between the resident population and the dynamic human activity variables of the LSOAs during 24 months. All the <span class="html-italic">p</span>-values of the coefficients (human activity variables) are statistically significant (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure A3
<p>Standardisedlocal coefficient values of explanatory variables in ‘Short-term shoplifting (S + D) models’. Each model’s training set is one month from January 2020, April 2020, August 2020, and November 2020, respectively.</p>
Full article ">Figure A4
<p>Maps of the spatial association (local coefficient values) between crime rates of shoplifting and the selected human activity variables in fit GWRs for February 2020, April 2020, August 2020, and November 2020. The bandwidths of four GWRs are 41, 71, 31, and 30, respectively. The global <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> of four GWRs are 0.69, 0.34, 0.71, and 0.61, respectively.</p>
Full article ">
22 pages, 4022 KiB  
Article
Vertical vs. Horizontal Fractal Dimensions of Roads in Relation to Relief Characteristics
by Klemen Prah and Ashton M. Shortridge
ISPRS Int. J. Geo-Inf. 2023, 12(12), 487; https://doi.org/10.3390/ijgi12120487 - 30 Nov 2023
Cited by 1 | Viewed by 1835
Abstract
This paper investigated the surface length of roads from both horizontal and vertical perspectives using the theory of fractal dimension of surfaces and curves. Three progressive experiments were conducted. The first demonstrated the magnitude of the differences between the planar road length and [...] Read more.
This paper investigated the surface length of roads from both horizontal and vertical perspectives using the theory of fractal dimension of surfaces and curves. Three progressive experiments were conducted. The first demonstrated the magnitude of the differences between the planar road length and the DTM-derived surface road length and assessed its correlation with the DTM-calculated road slope. The second investigated the road distance complexity through the fractal dimension in both planar and vertical dimensions. The third related the vertical with the horizontal fractal dimension of roads across a range of distinct physiographic regions. The study contributed theoretically by linking the planimetric complexity to vertical complexity, with clear applications for advanced transportation studies and network analyses. The core methodology used geographic information systems (GIS) to integrate a high resolution (1 × 1 m) digital terrain model (DTM) with a road network layer. A novel concept, the vertical fractal dimension of roads was introduced. Both the vertical and horizontal fractal dimensions of the roads were calculated using the box-counting methodology. We conducted an investigation into the relationship between the two fractal dimensions using fourteen study areas within four distinct physiographic regions across Slovenia. We found that the average slope of a three-dimensional (3D) road was directly related to the length difference between 3D and two-dimensional (2D) roads. The calculated values for the vertical fractal dimension in the study areas were only slightly above 1, while the maximum horizontal fractal dimension of 1.1837 reflected the more sinuous properties of the road in plan. Variations in the vertical and horizontal fractal dimensions of the roads varied between the different physiographic regions. Full article
Show Figures

Figure 1

Figure 1
<p>The research design was organized into three general, progressive stages (orange ellipses) and associated analysis details (blue rectangles).</p>
Full article ">Figure 2
<p>Vertical profile of a five-meter road section overlain on a one-meter cell size raster DTM processed using ArcGIS, according to the bilinear interpolation method, to obtain the Z, slope, and surface length information.</p>
Full article ">Figure 3
<p>Relationship between the length difference and average slope of the five-meter road sections in the municipality of Gornji Grad in Slovenia.</p>
Full article ">Figure 4
<p>Vertical profile of the sample road section. A one-hundred-meter-long detail of the road illustrated its vertical roughness.</p>
Full article ">Figure 5
<p>Procedure of box counting to obtain D<sub>2</sub>, the horizontal fractal dimension of a road in the Gornji Grad region. Ten different square side lengths were calculated, of which four are presented to illustrate the approach: (<b>a</b>) 1024 m, (<b>b</b>) 256 m, (<b>c</b>) 128 m, and (<b>d</b>) 16 m. The last example (<b>d</b>) presents in detail a 128-by-128 m area. In all the examples, the squares that intersect with the road are colored light blue. An elevation map provides context.</p>
Full article ">Figure 6
<p>Vertical and horizontal fractal dimensions for the Gornji Grad road section. The red square features represent the vertical profile, while the round dark blue features represent the horizontal profile. The dots indicate the paired log(ƞ):log(N) values that were empirically derived at different mesh side lengths. The regression model lines are depicted for each set of points, and each regression equation with parameter estimates and R<sup>2</sup> are reported. The vertical and horizontal fractal dimensions are given by the absolute value of the line slope in bold.</p>
Full article ">Figure 7
<p>Fourteen study areas, represented by numbered dark blue squares, spanning all four Slovenian macro-regions, including plains, low hills, high hills, and mountains. Source: Perko 1998, ARSO, GURS.</p>
Full article ">Figure 8
<p>Vertical and horizontal fractal dimensions of the road sections for all fourteen study sites (site numbers). The linear regression line and expression quantified the strength of the relationship between these dimensions.</p>
Full article ">Figure 9
<p>6-by-6 km excerpts from selected study sites with roads (red lines): (<b>a</b>) plain (site 1); (<b>b</b>) low hills (site 3); (<b>c</b>) hills (site 6); and (<b>d</b>) mountains (site 10). The shaded relief was generated from the DTM used in the study. The streams and rivers (blue lines) provide context. Source: ARSO, GURS, NAVTEQ.</p>
Full article ">Figure 10
<p>Vertical profiles of the road sections presented in <a href="#ijgi-12-00487-f009" class="html-fig">Figure 9</a>: (<b>a</b>) plains of study area 1; (<b>b</b>) low hills of study area 3; (<b>c</b>) hills of study area 6; and (<b>d</b>) mountains of study area 10. In each case, the lowest grid line represents the lowest elevation of the road section (also labeled on the axis). For better representation, the ratio between the horizontal and vertical axes was 1:2.</p>
Full article ">
21 pages, 12463 KiB  
Article
Measuring the Multiple Functions and Tradeoffs among Streets: A New Framework Using the Deep Learning Method
by Shihang Fu, Ying Fang, Nannan Wang, Zhaomin Tong and Yaolin Liu
ISPRS Int. J. Geo-Inf. 2023, 12(12), 486; https://doi.org/10.3390/ijgi12120486 - 29 Nov 2023
Viewed by 2080
Abstract
With the sustainable and coordinated development of cities, the formulation of urban street policies requires multiangle analysis. In regard to the existing street research, a large number of studies have focused on specific landscapes or accessibility of streets, and there is a lack [...] Read more.
With the sustainable and coordinated development of cities, the formulation of urban street policies requires multiangle analysis. In regard to the existing street research, a large number of studies have focused on specific landscapes or accessibility of streets, and there is a lack of research on the multiple functions of streets. Recent advances in sensor technology and digitization have produced a wealth of data and methods. Thus, we may comprehensively understand streets in a less labor-intensive way, not just single street functions. This paper defines an index system of the multiple functions of urban streets and proposes a framework for multifunctional street measurement. Via the application of deep learning to Baidu Street View (BSV) imagery, we generate three functions, namely, landscape, traffic, and economic functions. The results indicate that street facilities and features are suitably identified. According to the multifunctional perspective, this paper further classifies urban streets into multifunctional categories and provides targeted policy recommendations for urban street planning. There exist correlations among the various street functions, and the correlation between the street landscape and economic functions is highly significant. This framework can be widely applied in other countries and cities to better understand street differences in various cities. Full article
Show Figures

Figure 1

Figure 1
<p>Study area and sampled streets.</p>
Full article ">Figure 2
<p>The workflow of the analysis of multiple street functions.</p>
Full article ">Figure 3
<p>Attenuation coefficient with the POI distance.</p>
Full article ">Figure 4
<p>Semantic segmentation results of street view images.</p>
Full article ">Figure 5
<p>Distribution of the green vegetation factor (<b>a</b>), sky view factor (<b>b</b>), sidewalk factor (<b>c</b>), and landscape function (<b>d</b>) of streets.</p>
Full article ">Figure 6
<p>Distribution of street traffic safety (<b>a</b>) and traffic function (<b>b</b>).</p>
Full article ">Figure 7
<p>Spatial characteristics (<b>a</b>,<b>b</b>) and distribution (<b>c</b>,<b>d</b>) of TPBt and NQPD in the study area.</p>
Full article ">Figure 8
<p>Distribution of the street economic function.</p>
Full article ">Figure 9
<p>Coupling degree (<b>a</b>) and coordination degree (<b>b</b>) of the multiple functions of streets.</p>
Full article ">Figure 10
<p>Multifunctional street classification.</p>
Full article ">
19 pages, 29056 KiB  
Article
A Head/Tail Breaks-Based Approach to Characterizing Space-Time Risks of COVID-19 Epidemic in China’s Cities
by Tingting Wu, Bisong Hu, Jin Luo and Shuhua Qi
ISPRS Int. J. Geo-Inf. 2023, 12(12), 485; https://doi.org/10.3390/ijgi12120485 - 29 Nov 2023
Viewed by 2117
Abstract
The novel coronavirus pneumonia (COVID-19) pandemic has caused enormous impacts around the world. Characterizing the risk dynamics for urgent epidemics such as COVID-19 is of great benefit to epidemic control and emergency management. This article presents a novel approach to characterizing the space-time [...] Read more.
The novel coronavirus pneumonia (COVID-19) pandemic has caused enormous impacts around the world. Characterizing the risk dynamics for urgent epidemics such as COVID-19 is of great benefit to epidemic control and emergency management. This article presents a novel approach to characterizing the space-time risks of the COVID-19 epidemic. We analyzed the heavy-tailed distribution and spatial hierarchy of confirmed COVID-19 cases in 367 cities from 20 January to 12 April 2020, and population density data for 2019, and modelled two parameters, COVID-19 confirmed cases and population density, to measure the risk value of each city and assess the epidemic from the perspective of spatial and temporal changes. The evolution pattern of high-risk areas was assessed from a spatial and temporal perspective. The number of high-risk cities decreased from 57 in week 1 to 6 in week 12. The results show that the risk measurement model based on the head/tail breaks approach can describe the spatial and temporal evolution characteristics of the risk of COVID-19, and can better predict the risk trend of future epidemics in each city and identify the risk of future epidemics even during low incidence periods. Compared with the traditional risk assessment method model, it pays more attention to the differences in the spatial level of each city and provides a new perspective for the assessment of the risk level of epidemic transmission. It has generality and flexibility and provides a certain reference for the prevention of infectious diseases as well as a theoretical basis for government implementation strategies. Full article
(This article belongs to the Collection Spatial Components of COVID-19 Pandemic)
Show Figures

Figure 1

Figure 1
<p>Development of the COVID-19 epidemic in China ((<b>a</b>) shows the time-series variation in new cases of COVID-19 in China; (<b>b</b>) shows the spatial distribution of the cumulative number of confirmed cases of COVID-19 as of 12 April 2020).</p>
Full article ">Figure 2
<p>Rank-size plot (note: the x-axis is the sort and the y-axis is the size distribution of the corresponding sorted data values. The mean is m and the corresponding sort is R(m), with the head (e.g., 13%) and the tail (e.g., 87%)).</p>
Full article ">Figure 3
<p>The recursive flow of head/tail break (the 1001 numbers [1, 1/2, …, 1/1001] are divided into 5 levels, [1/134, 1/135, …, 1/1001], [1/25, 1/26, …, 1/133], [1/7, 1/8 …, 1/24], [1/3, 1/4, …, 1/6], [1, 1/2], with 5 inherent levels: L1, L2, L3, L4, L5).</p>
Full article ">Figure 4
<p>Population density and the heavy-tailed distribution of COVID-19 confirmed cases ((<b>a</b>,<b>b</b>) show the first three levels of rank-size plots for population density and first-week COVID-19 confirmed case head/tail break results, respectively; (<b>c</b>) shows the power-law variation for 12 weeks of COVID-19 confirmed cases; (<b>d</b>,<b>e</b>) show the trend of power-law index α and goodness-of-fit p, respectively).</p>
Full article ">Figure 5
<p>Time-series variation in population density and ht index of COVID-19 confirmed cases.</p>
Full article ">Figure 6
<p>Geographical distributions of the spatial hierarchy of population density and COVID-19 confirmed cases ((<b>a</b>) shows the hierarchical structure of the spatial hierarchy of population density; (<b>b</b>–<b>d</b>) show the change in the hierarchical structure of the spatial hierarchy in the 1st, 6th, and 12th weeks, respectively).</p>
Full article ">Figure 7
<p>Geographical distributions of high-risk cities (note: (<b>a</b>–<b>l</b>) represent the spatial distributions of high-risk cities from week 1 to week 12, respectively; the darker the colour, the higher the risk).</p>
Full article ">Figure 8
<p>Numbers of high-risk cities and statistics of the risk assessment indicator during 12 weeks (note: (<b>a</b>–<b>d</b>) show the number of cities at risk, mean, standard deviation, and maximum risk value changes per week, respectively).</p>
Full article ">Figure 9
<p>Time-series variation in risk, incidence, and the ht index rank of population migration intensity in six cities ((<b>a</b>–<b>f</b>) represent: Beijing, Guangzhou, Wuhan, Shanghai, Wenzhou, Chongqing; blue indicates risk values <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">R</mi> </mrow> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> </mrow> </semantics></math>, red indicates prevalence, black indicates population migration intensity ht index rank).</p>
Full article ">
27 pages, 870 KiB  
Article
Internet in the Middle of Nowhere: Performance of Geoportals in Rural Areas According to Core Web Vitals
by Karol Król and Wojciech Sroka
ISPRS Int. J. Geo-Inf. 2023, 12(12), 484; https://doi.org/10.3390/ijgi12120484 - 29 Nov 2023
Cited by 2 | Viewed by 2121
Abstract
The spatial planning system in Poland is undergoing a fundamental reform. It emphasises the digital representation of spatial data. Low performance of geoportals, no Internet access, or poor connectivity can contribute to the exclusion from the spatial planning process, and consequently to the [...] Read more.
The spatial planning system in Poland is undergoing a fundamental reform. It emphasises the digital representation of spatial data. Low performance of geoportals, no Internet access, or poor connectivity can contribute to the exclusion from the spatial planning process, and consequently to the exclusion from a specific part of public life. Considering these developments, the present study seems relevant by pointing out the issue with geoportal performance and availability of quality Internet in rural areas. The primary contribution of the article is (1) results of performance measurements for selected geoportals; (2) presentation of measuring tools and performance indices combined with methods for ad-hoc performance measuring; and (3) presentation of potential actions to improve geoportal performance on the device with which it is used. The article offers case studies where the performance of selected geoportals was tested in rural mountainous areas with limited Internet access. Five geoportals were tested with PageSpeed Insights (PSI), WebPageTest, GTmetrix, Pingdom, and GiftOfSpeed. Core Web Vitals indices were analysed: Largest Contentful Paint (LCP), First Input Delay (FID), Cumulative Layout Shift (CLS), and First Contentful Paint (FCP). The author verified values of the Speed Index and Fully Loaded Time along with other performance indices, like GTmetrix Structure. The study failed to provide unambiguous evidence that radio link users in rural areas could experience problems with geoportal performance, although the results seem to suggest it indirectly. PSI Lab Data and Field Data tests revealed a relatively low performance of the geoportals. The Performance index remained below 50 in most cases, which is ‘Poor’ according to the PSI scale. The fully loaded time exceeded 10 s for all the geoportals and 20 s in some cases (Lab Data). It means that the perceived performance of the tested geoportals on a radio link in rural areas is most probably even lower. The case studies demonstrated further that the user has limited possibilities to speed up map applications. It is possible to slightly improve the geoportal experience through the optimisation of the device locally, but the responsibility to ensure geoportal performance is mainly the publisher’s. Full article
Show Figures

Figure 1

Figure 1
<p>Research method conceptual diagram.</p>
Full article ">Figure 2
<p>Illustrative figure with the measurement location and Internet signal strength. Source: original work based on: RFBenchmark [<a href="#B64-ijgi-12-00484" class="html-bibr">64</a>].</p>
Full article ">
30 pages, 11203 KiB  
Article
Mixed-Methods Approach to Land Use Renewal Strategies in and around Abandoned Airports: The Case of Beijing Nanyuan Airport
by Haoxian Cai and Wei Duan
ISPRS Int. J. Geo-Inf. 2023, 12(12), 483; https://doi.org/10.3390/ijgi12120483 - 28 Nov 2023
Viewed by 2599
Abstract
Urban airports are typically large infrastructures with significant cultural, economic, and ecological impacts; meanwhile, abandoned airports are common worldwide. However, there is limited knowledge regarding transformation strategies for the renewal of abandoned airports and their surrounding regions in historically and culturally rich areas. [...] Read more.
Urban airports are typically large infrastructures with significant cultural, economic, and ecological impacts; meanwhile, abandoned airports are common worldwide. However, there is limited knowledge regarding transformation strategies for the renewal of abandoned airports and their surrounding regions in historically and culturally rich areas. We use Beijing’s Nanyuan Airport as a case study, combining the historic urban landscape approach, land use and land cover change, and counterfactual simulations of land use patterns to construct a comprehensive analytical framework. Our framework was used to analyze the long-term land use patterns of the study area, determine its value, and improve perception from a macro- and multi-perspective. We discovered that the traditional knowledge and planning systems in the study area have largely disappeared, but Nanyuan Airport’s impact on the surrounding land use patterns is unique and significant. By considering the characteristics and mechanisms of land use in the study area, we aimed to find a balance point between the historical context and future potential. As such, we propose optimized recommendations with the theme of connection and development engines. Our findings supplement the planning knowledge of relevant areas and provide a springboard for interdisciplinary research in landscape planning. Full article
Show Figures

Figure 1

Figure 1
<p>The present scenario in the vicinity of Nanyuan Airport (left photo shows the ongoing demolition and clearance work in Nanyuan, Beijing, to make way for future construction; right photo shows the protective forest area around the airport, taken by the authors in 2023).</p>
Full article ">Figure 2
<p>Location and scope of the Nanyuan research area and Nanyuan Airport.</p>
Full article ">Figure 3
<p>Flowchart of the cellular automata (CA)-Markov simulation workflow.</p>
Full article ">Figure 4
<p>Changes in the airport’s influence on the surrounding land use pattern were evaluated using conventional comparison (a) and counterfactual simulation (b).</p>
Full article ">Figure 5
<p>Research framework.</p>
Full article ">Figure 6
<p>A diagram of the evolution of the extent of Nanyuan throughout history [<a href="#B55-ijgi-12-00483" class="html-bibr">55</a>].</p>
Full article ">Figure 7
<p>The Qianlong Emperor Hunting Hare by Giuseppe Castiglione (1755).</p>
Full article ">Figure 8
<p>Zoning map of the Nanyuan site during the Kangxi period of the Qing dynasty (adapted from the First Historical Archives of China, “The Complete Map of Nanyuan in 1699”).</p>
Full article ">Figure 9
<p>Map of the royal hunting grounds in Nanyuan, Beijing, painted on wood at the end of the Qing dynasty.</p>
Full article ">Figure 10
<p>The main functions of Nanyuan during the Qing dynasty (self-drawn by the authors).</p>
Full article ">Figure 11
<p>The map of Nanyuan (39 × 53 cm). Scale: 1:25,000 (By G.Bouillard Ingènieur, Beijing, 1923 in Carte Des Environs De Peking).</p>
Full article ">Figure 12
<p>The 2019 land simulation and realistic land use classification results.</p>
Full article ">Figure 13
<p>(<b>a</b>) Keyhole satellite images and land use and land cover maps based on (<b>b</b>) 1984 TM images; (<b>c</b>) 1993 TM images; (<b>d</b>) 2001 ETM images; (<b>e</b>) 2011 ETM images; (<b>f</b>) 2019 OLI_TIRS images; (<b>g</b>) 2022 OLI_TIRS images in Nanyuan Airport study area.</p>
Full article ">Figure 14
<p>Real and simulated changes of various land use and land cover classifications for 1984, 1993, 2001, 2011, 2019, and 2022.</p>
Full article ">Figure 15
<p>Dependency wheels for land use flow changes among different years in the study area.</p>
Full article ">Figure 16
<p>The 2022 land simulation and realistic land use classification results.</p>
Full article ">Figure 17
<p>Land use change analysis under different methods; the traditional method compares real land use patterns in 2019 and 2022, and the counterfactual method compares real 2019 and simulated 2022 land use patterns. Note: Data shown in the figure are simulated or real 2022 LULC land use data minus real 2019 LULC land use data in hectares.</p>
Full article ">
27 pages, 8450 KiB  
Article
Impacts of Urban Morphology on Seasonal Land Surface Temperatures: Comparing Grid- and Block-Based Approaches
by Gyuwon Jeon, Yujin Park and Jean-Michel Guldmann
ISPRS Int. J. Geo-Inf. 2023, 12(12), 482; https://doi.org/10.3390/ijgi12120482 - 28 Nov 2023
Cited by 4 | Viewed by 2404
Abstract
Climate change is expected to result in increased occurrences of extreme weather events such as heat waves and cold spells. Urban planning responses are crucial for improving the capacity of cities and communities to deal with significant temperature variations across seasons. This study [...] Read more.
Climate change is expected to result in increased occurrences of extreme weather events such as heat waves and cold spells. Urban planning responses are crucial for improving the capacity of cities and communities to deal with significant temperature variations across seasons. This study aims to investigate the relationship between urban temperature fluctuations and urban morphology throughout the four seasons. Through quadrant and statistical analyses, built-environment factors are identified that moderate or exacerbate seasonal land surface temperatures (LSTs). The focus is on Seoul, South Korea, as a case study, and seasonal LST values are calculated at both the grid (100 m × 100 m) and street block levels, incorporating factors such as vegetation density, land use patterns, albedo, two- and three-dimensional building forms, and gravity indices for large forests and water bodies. The quadrant analysis reveals a spatial segregation between areas demonstrating high LST adaptability (cooler summers and warmer winters) and those displaying LST vulnerability (hotter summers and colder winters), with significant differences in vegetation and building forms. Spatial regression analyses demonstrate that higher vegetation density and proximity to water bodies play key roles in moderating LSTs, leading to cooler summers and warmer winters. Building characteristics have a constant impact on LSTs across all seasons: horizontal expansion increases the LST, while vertical expansion reduces the LST. These findings are consistent for both grid- and block-level analyses. This study emphasizes the flexible role of the natural environment in moderating temperatures. Full article
Show Figures

Figure 1

Figure 1
<p>Geographic location of the study area.</p>
Full article ">Figure 2
<p>Comparison of grid-based and street block-based units: (<b>a</b>) grid unit, (<b>b</b>) block unit.</p>
Full article ">Figure 3
<p>Seasonal distributions of land surface temperature (LST) in 2017, Seoul, South Korea: (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 4
<p>Partition of four quadrants as the intersection of summer and winter LSTs.</p>
Full article ">Figure 5
<p>Gravity index for urban forests (GIUF): (<b>a</b>) grid level, (<b>b</b>) block level.</p>
Full article ">Figure 6
<p>Gravity index for water bodies (GIWB): (<b>a</b>) grid level, (<b>b</b>) block level.</p>
Full article ">Figure 7
<p>Spatial distribution of the summer–winter LST quadrants in Seoul: (<b>a</b>) grid cell level; (<b>b</b>) street block level results.</p>
Full article ">
24 pages, 7987 KiB  
Article
Mapping Gross Domestic Product Distribution at 1 km Resolution across Thailand Using the Random Forest Area-to-Area Regression Kriging Model
by Yan Jin, Yong Ge, Haoyu Fan, Zeshuo Li, Yaojie Liu and Yan Jia
ISPRS Int. J. Geo-Inf. 2023, 12(12), 481; https://doi.org/10.3390/ijgi12120481 - 27 Nov 2023
Cited by 2 | Viewed by 2538
Abstract
Accurate spatial distribution of gridded gross domestic product (GDP) data is crucial for revealing regional disparities within administrative units, thus facilitating a deeper understanding of regional economic dynamics, industrial distribution, and urbanization trends. The existing GDP spatial models often rely on prediction residuals [...] Read more.
Accurate spatial distribution of gridded gross domestic product (GDP) data is crucial for revealing regional disparities within administrative units, thus facilitating a deeper understanding of regional economic dynamics, industrial distribution, and urbanization trends. The existing GDP spatial models often rely on prediction residuals for model evaluation or utilize residual distribution to improve the final accuracy, frequently overlooking the modifiable areal unit problem within residual distribution. This paper introduces a hybrid downscaling model that combines random forest and area-to-area kriging to map gridded GDP. Employing Thailand as a case study, GDP distribution maps were generated at a 1 km spatial resolution for the year 2015 and compared with five alternative downscaling methods and an existing GDP product. The results demonstrate that the proposed approach yields higher accuracy and greater precision in detailing GDP distribution, as evidenced by the smallest mean absolute error and root mean squared error values, which stand at USD 256.458 and 699.348 ten million, respectively. Among the four different sets of auxiliary variables considered, one consistently exhibited a higher prediction accuracy. This particular set of auxiliary variables integrated classification-based variables, illustrating the advantages of incorporating such integrated variables into modeling while accounting for classification characteristics. Full article
Show Figures

Figure 1

Figure 1
<p>Elevation and regional divisions (<b>a</b>), LULC (<b>b</b>), agriculture statistical GDP (<b>c</b>), and non-agriculture statistical GDP (<b>d</b>) of Thailand at the province level. (GDP in constant 2011 international USD).</p>
Full article ">Figure 2
<p>The result of KDE bandwidth determination.</p>
Full article ">Figure 3
<p>Tuning curves of random forest modeling for agricultural and non-agricultural GDP with four different scenarios. (The horizontal coordinate represents mtree).</p>
Full article ">Figure 4
<p>Flowchart of the GDP spatialization.</p>
Full article ">Figure 5
<p>The spatial distribution of downscaled GDP using RF, SVR, and MLR methods.</p>
Full article ">Figure 6
<p>The spatial distribution of downscaled GDP by using RFATARK, SVATARK, and MLATARK methods.</p>
Full article ">Figure 7
<p>The linear relationship between statistical GDP and spatialized GDP at the provincial scale. GDP derived by (<b>a</b>) RFATARK, (<b>b</b>) SVATARK, (<b>c</b>) MLATARK using MC group of auxiliary variables, and (<b>d</b>) G_GDP product. The black dashed line represents the 1:1 line, and the orange line represents the line fitted through the scatter points.</p>
Full article ">Figure 8
<p>The spatial distribution of downscaled GDP by using (<b>a</b>) RFIDW, (<b>b</b>) RF_MC, and (<b>c</b>) XGBoost methods and MC group auxiliary variables.</p>
Full article ">
21 pages, 3568 KiB  
Article
Research on Approximate Spatial Keyword Group Queries Based on Differential Privacy and Exclusion Preferences in Road Networks
by Liping Zhang, Jing Li and Song Li
ISPRS Int. J. Geo-Inf. 2023, 12(12), 480; https://doi.org/10.3390/ijgi12120480 - 26 Nov 2023
Cited by 1 | Viewed by 1757
Abstract
A new spatial keyword group query method is proposed in this paper to address the existing issue of user privacy leakage and exclusion of preferences in road networks. The proposed query method is based on the IGgram-tree index and minimum hash set. To [...] Read more.
A new spatial keyword group query method is proposed in this paper to address the existing issue of user privacy leakage and exclusion of preferences in road networks. The proposed query method is based on the IGgram-tree index and minimum hash set. To deal with this problem effectively, this paper proposes a query method based on the IGgram-tree index and minimum hash set. The IGgram-tree index is proposed for the first time to deal with the approximate keyword query problem in the road network. This index significantly improves the efficiency of calculating the road network distance and querying approximate keywords. Considering that spatial keyword group queries are caused by NP-hard problems with high time complexity, this paper proposes a data structure that uses the minimum hash set, which can efficiently search for the result set. To address the problem that the traditional spatial keyword group query does not consider user privacy leakage and the limitations of existing privacy protection techniques, this method proposes a differential privacy-based allocation method to better protect the privacy of data. The theoretical study and experimental analysis show that the proposed method can better handle the approximate spatial keyword group query problem based on its use of differential privacy and exclusion preferences in road networks. Full article
Show Figures

Figure 1

Figure 1
<p>Algorithm relationship and data-processing flow.</p>
Full article ">Figure 2
<p>Illustration of road network environment.</p>
Full article ">Figure 3
<p>Illustration of road network environmental partition.</p>
Full article ">Figure 4
<p>The IGgram-tree index.</p>
Full article ">Figure 5
<p>Illustration of MH’s structure<b>.</b></p>
Full article ">Figure 6
<p>Effect of the number of query positive keywords on the efficiency of the algorithm and the number of expansion nodes.</p>
Full article ">Figure 7
<p>Effect of the number of query rejection keywords on the efficiency of the algorithm and the number of expansion nodes<b>.</b></p>
Full article ">Figure 8
<p>Algorithm accuracy<b>.</b></p>
Full article ">Figure 9
<p>Availability comparison.</p>
Full article ">Figure 10
<p>Impact of the number of keywords under different privacy-preserving methods on accuracy.</p>
Full article ">
20 pages, 4054 KiB  
Article
Understanding Map Misinterpretation: Factors Influencing Correct Map Reading and Common Errors
by Csaba Szigeti-Pap, Dávid Kis and Gáspár Albert
ISPRS Int. J. Geo-Inf. 2023, 12(12), 479; https://doi.org/10.3390/ijgi12120479 - 26 Nov 2023
Cited by 2 | Viewed by 2646
Abstract
Misinterpreting maps can have serious consequences, especially in situations requiring quick decisions like using car navigation systems. Studies indicate that a map reader’s experience is crucial for understanding maps, but factors such as age, education, and gender can also influence interpretation. However, understanding [...] Read more.
Misinterpreting maps can have serious consequences, especially in situations requiring quick decisions like using car navigation systems. Studies indicate that a map reader’s experience is crucial for understanding maps, but factors such as age, education, and gender can also influence interpretation. However, understanding only the proportion of correctly interpreted information is not enough. It is essential to investigate the types of mistakes made and their causes. To address this, we conducted a study available in six languages with 511 participants who completed an online questionnaire testing their map reading skills. The questions focused on scale usage, mental rotation, and recognizing map categories (relief, line and point symbols, and geographic names). Gender had significant relation with one skill, qualification with two and age with three. Experience was associated to the highest number of skills, a total of four, confirming previous findings. When making mistakes, participants tended to overestimate distances and struggled with conceptual similarities in symbol recognition. Experienced readers often misplaced reference locations of geographic names. The results of the research could be used in the design of large-scale maps (e.g., car navigation), as they allow to reduce typical map reading errors by careful selection of symbol types and placements. Full article
Show Figures

Figure 1

Figure 1
<p>The relief map for the test used for questions Q1–3.</p>
Full article ">Figure 2
<p>The relief and hydrography map for the test used for question Q4.</p>
Full article ">Figure 3
<p>The simplified topographic map for the test used for question Q5, and Q6.</p>
Full article ">Figure 4
<p>The last map for the test showing linear and polygonal map symbols and geographic names was used for question Q7, and Q8. Note: the map legend was not included in the test.</p>
Full article ">Figure 5
<p>The visual depiction of the study procedure.</p>
Full article ">Figure 6
<p>The relative proportion of votes given to incorrect answer options in the case of four-choice questions (the skills tested by the questions are shown on the left). Capital letters indicate the answer variations shown in <a href="#sec2dot1-ijgi-12-00479" class="html-sec">Section 2.1</a>. The answer “I don’t know” was not counted as one of the wrong answers and was therefore not included in the analysis.</p>
Full article ">Figure 7
<p>The relative frequency distribution of mistakes grouped by the different levels of erroneous performance.</p>
Full article ">
14 pages, 3856 KiB  
Article
Spatial Accessibility of Public Electric Vehicle Charging Services in China
by Yu Chen, Yuehong Chen and Yuqi Lu
ISPRS Int. J. Geo-Inf. 2023, 12(12), 478; https://doi.org/10.3390/ijgi12120478 - 25 Nov 2023
Cited by 1 | Viewed by 2888
Abstract
Decarbonizing the transport sector using electric vehicles (EVs) is a vital pathway for China to achieve the carbon peak and carbon neutrality goals. Despite the unprecedented growth of EV diffusion in China, little information is available for the spatial accessibility of public electric [...] Read more.
Decarbonizing the transport sector using electric vehicles (EVs) is a vital pathway for China to achieve the carbon peak and carbon neutrality goals. Despite the unprecedented growth of EV diffusion in China, little information is available for the spatial accessibility of public electric vehicle charging services (EVCSs). This study developed an applicable accessibility measurement framework to examine the city-level accessibility of EVCSs in China using the Gaussian two-step floating catchment area (G2SFCA) method. G2SFCA takes the EV charging stations with charging piles as supply and the EV ownership data as demand. The results indicate that (1) the eastern region of China has the highest density of EV charging stations (69.1%), followed by the central region, while the western region has the lowest density; (2) the spatial accessibility of EVCSs has a different pattern, where the central region has the highest accessibility, followed by the eastern and western regions; (3) the spatial mismatch between EVCSs and EV diffusion in the eastern region is larger than that of the other two regions, which may be attributed to the suboptimal layout of EV charging stations and the inconsistent pace between EV penetration and EV charging station construction; and (4) there is a significant spatial inequity in the accessibility of EVCSs across both all three regions and the entirety of China, with the western region exhibiting the highest inequity, followed by the central and eastern regions. Based on these findings, policy implications are drawn for different regions in China, which may aid policymakers in crafting strategic policies and subsidy programs to foster the advancement of EVCSs. Full article
Show Figures

Figure 1

Figure 1
<p>Development of EVCSs in China (2006–2021).</p>
Full article ">Figure 2
<p>Experimental data. (<b>a</b>) Spatial distribution of public EV charging stations and (<b>b</b>) city-level ownership data of EVs.</p>
Full article ">Figure 3
<p>Research framework.</p>
Full article ">Figure 4
<p>EV charging stations and piles in each province.</p>
Full article ">Figure 5
<p>EV charging stations per 1000 EV cars (<b>a</b>) and their LISA (<b>b</b>).</p>
Full article ">Figure 6
<p>City-level spatial accessibility of EVCSs (<b>a</b>) and its LISA in China (<b>b</b>).</p>
Full article ">
14 pages, 4193 KiB  
Article
Analysis of Spatial and Temporal Distribution Patterns of Traditional Opera Culture along the Beijing–Hangzhou Grand Canal
by Jiayi Yang, Di Hu, Zihan Chen, Yicheng Xu, Zewei Zou and Ying Zhu
ISPRS Int. J. Geo-Inf. 2023, 12(12), 477; https://doi.org/10.3390/ijgi12120477 - 25 Nov 2023
Cited by 4 | Viewed by 1934
Abstract
As an exquisite asset of Chinese traditional culture, traditional opera occupies a place of high esteem within the world’s cultural and artistic treasury. The impact of emerging cultures has threatened the future of traditional opera culture, necessitating a thorough examination of the historical [...] Read more.
As an exquisite asset of Chinese traditional culture, traditional opera occupies a place of high esteem within the world’s cultural and artistic treasury. The impact of emerging cultures has threatened the future of traditional opera culture, necessitating a thorough examination of the historical context of the Grand Canal and traditional opera. There is insufficient research on the spatial evolution of the traditional opera culture along the Grand Canal; thus, this study takes ancient opera stages, a representative cultural relic of traditional opera, as an entry point and employs methods such as kernel density analysis and standard deviation ellipse analysis to analyze the spatial and temporal distribution patterns of the traditional opera culture along the Grand Canal. The results showed that: (i) Nationwide, opera stages in the areas along the Grand Canal exhibit a significant clustering characteristic. (ii) The changes in the number and locations of opera stages in the areas along the Grand Canal are closely related to the rise and fall of the Canal. The opera stages emerged along the Canal, gradually prospered with the development of the Canal, and finally clustered in a band-like cluster along the Grand Canal. (iii) From the Ming Dynasty to the founding of the People’s Republic of China, the opera stages in the areas along the Grand Canal spread in the “southeast–northwest” direction, which was consistent with the main direction of the Grand Canal, indicating its driving influence. (iv) On the centennial scale, from the 14th century to the 20th century, the evolution characteristics of the distribution centroid of opera stages in the areas along the Grand Canal were closely related to the key time nodes of Grand Canal construction and basin expansion. This study reveals the relationship between the Grand Canal and the spatial pattern evolution of traditional opera culture, aiming to promote the construction of the Grand Canal cultural belt. Full article
Show Figures

Figure 1

Figure 1
<p>Distribution of traditional opera stages in China.</p>
Full article ">Figure 2
<p>Spatial distribution of opera stages in China and in the eight provinces and cities.</p>
Full article ">Figure 3
<p>Standard deviation ellipses for different regions during different periods.</p>
Full article ">Figure 3 Cont.
<p>Standard deviation ellipses for different regions during different periods.</p>
Full article ">Figure 4
<p>The kernel density distribution of opera stages in the eight provinces and cities along the Grand Canal in different periods.</p>
Full article ">Figure 5
<p>Standard deviation ellipses of important time nodes from before the 14th century to the 20th century.</p>
Full article ">
19 pages, 4063 KiB  
Article
Effects of Spatial Reference Frames, Map Dimensionality, and Navigation Modes on Spatial Orientation Efficiency
by Hongyun Guo, Nai Yang, Zhong Wang and Hao Fang
ISPRS Int. J. Geo-Inf. 2023, 12(12), 476; https://doi.org/10.3390/ijgi12120476 - 23 Nov 2023
Viewed by 1893
Abstract
How can the interactive mode of a map be optimized to facilitate efficient positioning and improve cognitive efficiency? This paper addresses this crucial aspect of map design. It explores the impact of spatial reference frames, map dimensionality, and navigation modes on spatial orientation [...] Read more.
How can the interactive mode of a map be optimized to facilitate efficient positioning and improve cognitive efficiency? This paper addresses this crucial aspect of map design. It explores the impact of spatial reference frames, map dimensionality, and navigation modes on spatial orientation efficiency, as well as their interactions, through empirical eye-movement experiments. The results demonstrate the following: (1) When using a 2D fixed map in an allocentric reference frame, participants exhibit a high correct rate, a low cognitive load, and a short reaction time. In contrast, when operating within an egocentric reference frame using a 2D rotating map, participants demonstrate a higher correct rate, a reduced cognitive load, and a quicker reaction time. (2) The simplicity of 2D maps, despite their reduced authenticity compared to 3D maps, diminishes users’ cognitive load and enhances positioning efficiency. (3) The fixed map aligns more closely with the cognitive habits of participants in the allocentric reference frame, while the rotating map corresponds better to the cognitive habits of participants in the egocentric reference frame, thereby improving their cognitive efficiency. This study offers insights that can inform the optimization design of spatial orientation efficiency. Full article
Show Figures

Figure 1

Figure 1
<p>Demonstration diagram of animal testing.</p>
Full article ">Figure 2
<p>Animal test results of two different participants.</p>
Full article ">Figure 3
<p>Four experimental materials.</p>
Full article ">Figure 4
<p>Interaction plots of spatial reference frame and map dimensionality on accuracy rate.</p>
Full article ">Figure 5
<p>Interaction plots of spatial reference frame and navigation method on accuracy rate.</p>
Full article ">Figure 6
<p>Interaction plots of navigation mode and map dimensionality on reaction time.</p>
Full article ">Figure 7
<p>Interaction plots of spatial reference frame and navigation mode on reaction time.</p>
Full article ">Figure 8
<p>Interaction plots of navigation mode and map dimensionality on fixation duration.</p>
Full article ">Figure 9
<p>Interaction plots of spatial reference frame and map dimensionality on fixation duration.</p>
Full article ">Figure 10
<p>Interaction plots of spatial reference frame and navigation mode on fixation duration.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop