A quantitative framework for analysing long term spatial clustering and vegetation fragmentation in an urban landscape using multi-temporal landsat data

Rapid urbanization threatens urban green spaces and vegetation, demonstrated by a decrease in connectivity and higher levels of fragmentation. Understanding historic spatial and temporal patterns of such fragmentation is important for habitat and biological conservation, ecosystem management and urban planning. Despite their potential value, Local Indicators of Spatial Autocorrelation (LISA) measures have not been sufficiently exploited in monitoring the spatial and temporal variability in clustering and fragmentation of vegetation patterns in urban areas. LISA statistics are an important structural measure that indicates the presence of outliers, zones of similarity (hot spots) and of dissimilarity (cold spots) at proximate locations, hence they could be used to explicitly capture spatial patterns that are clustered, dispersed or random. In this study, we applied landscape metrics, LISA indices to analyse the temporal variability in clustering and fragmentation patterns of vegetation patches in Harare metropolitan city, Zimbabwe using Landsat series data for 1994, 2001 and 2017. Analysis of landscape metrics showed an increase in the fragmentation of vegetation patches between 1994–2017 as shown by the decrease in mean patch size, an increase in number of patches, edge density and shape complexity of vegetation patches. The study further demonstrates the utility of LISA indices in identifying key hot spot and cold spots. Comparatively, the highly vegetated northern parts of the city were characterised by significantly high positive spatial autocorrelation (p < 0.05) of vegetation patches. Conversely, more dispersed vegetation patches were found in the highly and densely urbanized western, eastern and southern parts of the city. This suggest that with increasing vegetation fragmentation, small and isolated vegetation patches do not spatially cluster but are dispersed geographically. The findings of the study underline the potential of LISA measures as a valuable spatially explicit method for the assessment of spatial clustering and fragmentation of urban vegetation patterns.     

Urban event detection with big data of taxi OD trips: A time series decomposition approach

Big urban mobility data, such as taxi trips, cell phone records, and geo‐social media check‐ins, offer great opportunities for analyzing the dynamics, events, and spatiotemporal trends of the urban social landscape. In this article, we present a new approach to the detection of urban events based on location‐specific time series decomposition and outlier detection. The approach first extracts long‐term temporal trends and seasonal periodicity patterns. Events are defined as anomalies that deviate significantly from the prediction with the discovered temporal patterns, i.e., trend and periodicity. Specifically, we adopt the STL approach, i.e., seasonal and trend decomposition using LOESS (locally weighted scatterplot smoothing), to decompose the time series for each location into three components: long‐term trend, seasonal periodicity, and the remainder. Events are extracted from the remainder component for each location with an outlier detection method. We analyze over a billion taxi trips for over seven years in Manhattan (New York City) to detect and map urban events at different temporal resolutions. Results show that the approach is effective and robust in detecting events and revealing urban dynamics with both holistic understandings and location‐specific interpretations.     

Sequence-based mapping approach to spatio-temporal snow patterns from MODIS time-series applied to Scotland

Snow cover and its monitoring are important because of the impact on important environmental variables, hydrological circulation and ecosystem services. For regional snow cover mapping and monitoring, the MODIS satellite sensors are particularly appealing. However cloud presence is an important limiting factor. This study addressed the problem of cloud cover for time-series in a boreal-Atlantic region where melting and re-covering of snow often do not follow the usual alpine-like patterns. A key requirement in this context was to apply improved methods to deal with the high cloud cover and the irregular spatio-temporal snow occurrence, through exploitation of space-time correlation of pixel values. The information contained in snow presence sequences was then used to derive summary indices to describe the time series patterns. Finally it was tested whether the derived indices can be considered an accurate summary of the snow presence data by establishing and evaluating their statistical relations with morphology and the landscape. The proposed cloud filling method had a good agreement (between 80 and 99%) with validation data even with a large number of pixels missing. The sequence analysis algorithm proposed takes into account the position of the states to fully consider the temporal dimension, i.e. the order in which a certain state appears in an image sequence compared to its neighbourhoods. The indices that were derived from the sequence of snow presence proved useful for describing the general spatio-temporal patterns of snow in Scotland as they were well related (more than 60% of explained deviance) with environmental information such as morphology supporting their use as a summary of snow patterns over time. The use of the derived indices is an advantage because of data reduction, easier interpretability and capture of sequence position-wise information (e.g. importance of short term fall/melt cycles). The derived seven clusters took into account the temporal patterns of the snow presence and they were well separated both spatially and according to the snow patterns and the environmental information. In conclusion, the use of sequences proved useful for analysing different spatio-temporal patterns of snow that could be related to other environmental information to characterize snow regimes regions in Scotland and to be integrated with ground measures for further hydrological and climatological analysis as baseline data for climate change models.     

Analytical methods for error propagation in planar space–time prisms

The space–time prism demarcates all locations in space–time that a mobile object or person can occupy during an episode of potential or unobserved movement. The prism is central to time geography as a measure of potential mobility and to mobile object databases as a measure of location possibilities given sampling error. This paper develops an analytical approach to assessing error propagation in space–time prisms and prism–prism intersections. We analyze the geometry of the prisms to derive a core set of geometric problems involving the intersection of circles and ellipses. Analytical error propagation techniques such as the Taylor linearization method based on the first-order partial derivatives are not available since explicit functions describing the intersections and their derivatives are unwieldy. However, since we have implicit functions describing prism geometry, we modify this approach using an implicit function theorem that provides the required first-order partials without the explicit expressions. We describe the general method as well as details for the two spatial dimensions case and provide example calculations.     

Detection of clusters in traffic networks based on spatio‐temporal flow modeling

Spatio‐temporal clustering is a highly active research topic and a challenging issue in spatio‐temporal data mining. Many spatio‐temporal clustering methods have been designed for geo‐referenced time series. Under some special circumstances, such as monitoring traffic flow on roads, existing methods cannot handle the temporally dynamic and spatially heterogeneous correlations among road segments when detecting clusters. Therefore, this article develops a spatio‐temporal flow‐based approach to detect clusters in traffic networks. First, a spatio‐temporal flow process is modeled by combining network topology relations with real‐time traffic status. On this basis, spatio‐temporal neighborhoods are captured by considering traffic time‐series similarity in spatio‐temporal flows. Spatio‐temporal clusters are further formed by successive connection of spatio‐temporal neighbors. Experiments on traffic time series of central London's road network on both weekdays and weekends are performed to demonstrate the effectiveness and practicality of the proposed method.     

Dynamics of vegetation indices in tropical and subtropical savannas defined by ecoregions and Moderate Resolution Imaging Spectroradiometer (MODIS) land cover

In this study, we explored the capacity of vegetation indices derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance products to characterize global savannas in Australia, Africa and South America. The savannas were spatially defined and subdivided using the World Wildlife Fund (WWF) global ecoregions and MODIS land cover classes. Average annual profiles of Normalized Difference Vegetation Index, shortwave infrared ratio (SWIR32), White Sky Albedo (WSA) and the Structural Scattering Index (SSI) were created. Metrics derived from average annual profiles of vegetation indices were used to classify savanna ecoregions. The response spaces between vegetation indices were used to examine the potential to derive structural and fractional cover measures. The ecoregions showed distinct temporal profiles and formed groups with similar structural properties, including higher levels of woody vegetation, similar forest–savanna mixtures and similar grassland predominance. The potential benefits from the use of combinations of indices to characterize savannas are discussed.     

Multi‐scale approach to mining significant spatial co‐location patterns

Spatial co‐location pattern mining aims to discover a collection of Boolean spatial features, which are frequently located in close geographic proximity to each other. Existing methods for identifying spatial co‐location patterns usually require users to specify two thresholds, i.e. the prevalence threshold for measuring the prevalence of candidate co‐location patterns and distance threshold to search the spatial co‐location patterns. However, these two thresholds are difficult to determine in practice, and improper thresholds may lead to the misidentification of useful patterns and the incorrect reporting of meaningless patterns. The multi‐scale approach proposed in this study overcomes this limitation. Initially, the prevalence of candidate co‐location patterns is measured statistically by using a significance test, and a non‐parametric model is developed to construct the null distribution of features with the consideration of spatial auto‐correlation. Next, the spatial co‐location patterns are explored at multi‐scales instead of single scale (or distance threshold) discovery. The validity of the co‐location patterns is evaluated based on the concept of lifetime. Experiments on both synthetic and ecological datasets show that spatial co‐location patterns are discovered correctly and completely by using the proposed method; on the other hand, the subjectivity in discovery of spatial co‐location patterns is reduced significantly.     

A Back‐Propagation Neural Network‐Based Approach for Multi‐Represented Feature Matching in Update Propagation

Spatial data infrastructures, which are characterized by multi‐represented datasets, are prevalent throughout the world. The multi‐represented datasets contain different representations for identical real‐world entities. Therefore, update propagation is useful and required for maintaining multi‐represented datasets. The key to update propagation is the detection of identical features in different datasets that represent corresponding real‐world entities and the detection of changes in updated datasets. Using polygon features of settlements as examples, this article addresses these key problems and proposes an approach for multi‐represented feature matching based on spatial similarity and a back‐propagation neural network (BPNN). Although this approach only utilizes the measures of distance, area, direction and length, it dynamically and objectively determines the weight of each measure through intelligent learning; in contrast, traditional approaches determine weight using expertise. Therefore, the weight may be variable in different data contexts but not for different levels of expertise. This approach can be applied not only to one‐to‐one matching but also to one‐to‐many and many‐to‐many matching. Experiments are designed using two different approaches and four datasets that encompass an area in China. The goals are to demonstrate the weight differences in different data contexts and to measure the performance of the BPNN‐based feature matching approach.     

Comparative study on two methods for calculating the gravitational potential of a prism

The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy’s approach provided a closed expression, one does not know whether it is valid, due to the fact that this approach has not been confirmed in literatures. In this paper, a study on the comparison of Nagy’s approach with discrete element method is presented. The results show that Nagy’s formulas for determining the gravitational potential of a prism are valid in the domain both inside and outside the prism.     

Comparative study on two methods for calculating the gravitational potential of a prism

The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy’s approach provi...     

Measuring differential access to facilities between population groups using spatial Lorenz curves and related indices

Access to certain types of facilities can promote health and well‐being. When population and facilities are not uniformly distributed across the landscape, inequities in accessibility may occur. Current research into GIS‐based accessibility measures has focused primarily on spatial inequities between different geographic locations but not directly on differences in accessibility between subgroups of the population. The research presented here develops a new method for measuring differential accessibility to facilities between various segments of the population. The method extends concepts and techniques in spatial point pattern analysis that account for the spatial structure of demand and its relationship to supply. In this approach, the traditional Lorenz curve and its associated indices, the Gini coefficient and the dissimilarity index, which are used to measure inequality, are recast in spatial terms for measuring differences in accessibility between population subgroups. An analysis of spatial accessibility to grocery stores in Akron, OH illustrates the value of the spatial Lorenz curve and its associated indices compared to other methods.     

Development of a GIS-based Methodology for Quantifying Predation Risk in a Spatial Context

We used GIS‐based techniques and inductive‐spatial modelling to analyse spatial interactions between predator and prey populations, a unique example of an extension to an existing ecological approach. The approach identifies spatial predictors of predation and is demonstrated for braided riverbed habitat in the Mackenzie Basin, New Zealand. This habitat is highly modified and is the focus of intensive management because of the presence of introduced mammalian predators and their native prey (primarily nests of ground‐nesting native birds that are extremely vulnerable to predators). Native prey are currently in decline due to predation; GIS techniques therefore have the potential to augment conventional ecological techniques and management practices. Predation risk to nests of ground‐nesting birds was interpolated across an area using inverse distance weighting. A linear regression model for spatio‐temporal predictors of nest predation was deter‐mined using a predation risk surface and three separate surfaces representing predator abundance, prey abundance, and habitat structure. These measures can enhance the non‐spatial results from conventional ecological inquiry and help solve pest management problems affecting the natural environment.     

STAMP: spatial–temporal analysis of moving polygons

Research questions regarding temporal change in spatial patterns are increasingly common in geographical analysis. In this research, we explore and extend an approach to the spatial–temporal analysis of polygons that are spatially distinct and experience discrete changes though time. We present five new movement events for describing spatial processes: displacement, convergence, divergence, fragmentation and concentration. Spatial–temporal measures of events for size and direction are presented for two time periods, and multiple time periods. Size change metrics are based on area overlaps and a modified cone-based model is used for calculating polygon directional relationships. Quantitative directional measures are used to develop application specific metrics, such as an estimation of the concentration parameter for a von Mises distribution, and the directional rate of spread. The utility of the STAMP methods are demonstrated by a case study on the spread of a wildfire in northwestern Montana.      

Discovering Spatial Patterns in Origin‐Destination Mobility Data

Mobility and spatial interaction data have become increasingly available due to the wide adoption of location‐aware technologies. Examples of mobility data include human daily activities, vehicle trajectories, and animal movements, among others. In this article we focus on a special type of mobility data, i.e. origin‐destination pairs, and present a new approach to the discovery and understanding of spatio‐temporal patterns in the movements. Specifically, to extract information from complex connections among a large number of point locations, the approach involves two steps: (1) spatial clustering of massive GPS points to recognize potentially meaningful places; and (2) extraction and mapping of the flow measures of clusters to understand the spatial distribution and temporal trends of movements. We present a case study with a large dataset of taxi trajectories in Shenzhen, China to demonstrate and evaluate the methodology. The contribution of the research is two‐fold. First, it presents a new methodology for detecting location patterns and spatial structures embedded in origin‐destination movements. Second, the approach is scalable to large data sets and can summarize massive data to facilitate pattern extraction and understanding.     

Spatio‐Temporal Building Population Estimation for Highly Urbanized Areas Using GIS

Detailed population information is crucial for the micro‐scale modeling and analysis of human behavior in urban areas. Since it is not available on the basis of individual persons, it has become necessary to derive data from aggregated census data. A variety of approaches have been published in the past, yet they are not entirely suitable for use in the micro‐scale context of highly urbanized areas, due mainly to their broad spatial scale and missing temporal scale. Here we introduce an enhanced approach for the spatio‐temporal estimation of building populations in highly urbanized areas. It builds upon other estimation methodologies, but extends them by introducing multiple usage categories and the temporal dimension. This allows for a more realistic representation of human activities in highly urbanized areas and the fact that populations change over time as a result of these activities. The model makes use of a variety of micro‐scale data sets to operationalize the activities and their spatio‐temporal representations. The outcome of the model provides estimated population figures for all buildings at each time step and thereby reveals spatio‐temporal behavior patterns. It can be used in a variety of applications concerning the implications of human behavior in urban areas.     

Multidimensional Similarity Measuring for Semantic Trajectories

Most existing approaches aiming at measuring trajectory similarity are focused on two‐dimensional sequences of points, called raw trajectories. However, recent proposals have used background geographic information and social media data to enrich these trajectories with a semantic dimension, giving rise to the concept of semantic trajectories. Only a few works have proposed similarity measures for semantic trajectories or multidimensional sequences, having limitations such as predefined weight of the dimensions, sensitivity to noise, tolerance for gaps with different sizes, and the prevalence of the worst dimension similarity. In this article we propose MSM, a novel similarity measure for multidimensional sequences that overcomes the aforementioned limitations by considering and weighting the similarity in all dimensions. MSM is evaluated through an extensive experimental study that, based on a seed trajectory, creates sets of semantic trajectories with controlled transformations to introduce different kinds and levels of dissimilarity. For each set, we compute the similarity between the seed and the transformed trajectories, using different measures. The results showed that MSM was more robust and efficient than related approaches in the domain of semantic trajectories.     

A Space‐Time Raster GIS Data Model for Spatiotemporal Analysis of Vegetation Responses to a Freeze Event

Many past space‐time GIS data models viewed the world mainly from a spatial perspective. They attached a time stamp to each state of an entity or the entire area of study. This approach is less efficient for certain spatio‐temporal analyses that focus on how locations change over time, which require researchers to view each location from a temporal perspective. In this article, we present a data model to organize multi‐temporal remote sensing datasets and track their changes at the individual pixel level. This data model can also integrate raster datasets from heterogeneous sources under a unified framework. The proposed data model consists of several object classes under a hierarchical structure. Each object class is associated with specific properties and behaviors to facilitate efficient spatio‐temporal analyses. We apply this data model to a case study of analyzing the impact of the 2007 freeze in Knoxville, Tennessee. The characteristics of different vegetation clusters before, during, and after the 2007 freeze event are compared. Our findings indicate that the majority of the study area is impacted by this freeze event, and different vegetation types show different response patterns to this freeze.     

Measures of spatio-temporal accuracy for time series land cover data

Remote sensing is a useful tool for monitoring changes in land cover over time. The accuracy of such time-series analyses has hitherto only been assessed using confusion matrices. The matrix allows global measures of user, producer and overall accuracies to be generated, but lacks consideration of any spatial aspects of accuracy. It is well known that land cover errors are typically spatially auto-correlated and can have a distinct spatial distribution. As yet little work has considered the temporal dimension and investigated the persistence or errors in both geographic and temporal dimensions. Spatio-temporal errors can have a profound impact on both change detection and on environmental monitoring and modelling activities using land cover data. This study investigated methods for describing the spatio-temporal characteristics of classification accuracy. Annual thematic maps were created using a random forest classification of MODIS data over the Jakarta metropolitan areas for the period of 2001–2013. A logistic geographically weighted model was used to estimate annual spatial measures of user, producer and overall accuracies. A principal component analysis was then used to extract summaries of the multi-temporal accuracy. The results showed how the spatial distribution of user and producer accuracy varied over space and time, and overall spatial variance was confirmed by the principal component analysis. The results indicated that areas of homogeneous land cover were mapped with relatively high accuracy and low variability, and areas of mixed land cover with the opposite characteristics. A multi-temporal spatial approach to accuracy is shown to provide more informative measures of accuracy, allowing map producers and users to evaluate time series thematic maps more comprehensively than a standard confusion matrix approach. The need to identify suitable properties for a temporal kernel are discussed.     

A Space-Time GIS Approach to Exploring Large Individual-based Spatiotemporal Datasets

The increasing number of large individual-based spatiotemporal datasets in various research fields has challenged the GIS community to develop analysis tools that can efficiently help researchers explore the datasets in order to uncover useful information. Rooted in Hägerstrand's time geography, this study presents a generalized space-time path (GSTP) approach to facilitating visualization and exploration of spatiotemporal changes among individuals in a large dataset. The fundamental idea of this approach is to derive a small number of representative space-time paths (i.e. GSTPs) from the raw dataset by identifying spatial cluster centers of observed individuals at different time periods and connecting them according to their temporal sequence. A space-time GIS environment is developed to implement the GSTP concept. Different methods of handling temporal data aggregation and the creation of GSTPs are discussed in this article. Using a large individual-based migration history dataset, this study successfully develops an operational space-time GIS prototype in ESRI's ArcScene and ArcMap to provide a proof-of-concept study of this approach. This space-time GIS system demonstrates that the proposed GSTP approach can provide a useful exploratory analysis and geovisualization environment to help researchers effectively search for hidden patterns and trends in such datasets.