Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [spatial-statistics]

The field of study concerning statistical methods that use space and spatial relationships (such as distance, area, volume, length, height, orientation, centrality and/or other spatial characteristics of data) directly in their mathematical computations.

1
vote
0answers
13 views

Raster correlation: Error (arg6) in raster.modified.ttest (R)

I am performing pairwise spatial correlation of 6 rasters using the raster.modified.ttest() function in R, as per this very useful post. I have consistently had the following error with my data; ...
0
votes
0answers
41 views

Seeking ideas for geostatistics applied assignment? [closed]

Any ideas for Geo statistics applied assignment, including comparison between two methods?
0
votes
0answers
9 views

Extract random effects from CARBayes output

We conducted a study at the U.S. census tract level looking at the relationship between environmental factors and social mobility in the ten largest U.S. cities. We used CARBayes::S.CARmultilevel in R ...
0
votes
0answers
46 views

Kriging in a large geographical area in QGIS

I am aware that already there are some answers about Kriging, but those answers couldn't solve my problem. I am trying to use Kriging to interpolate my data in a large area; in fact, the area is 400 ...
2
votes
1answer
30 views

How do I plot spatialpolygon dataframe layer on a leaflet map

I am not sure why this piece of code works library(leaflet) library(raster) adm <- getData('GADM', country='SGP', level=1) popup <- paste0("<strong>Name: </strong>", ...
1
vote
0answers
24 views

BIOMOD2: error message: Error in .local(sp, env, …) : object 'mask.tmp' not found

I am trying to model a plant species (composed by 298 occurrences) with BIOMOD2 using the consensus landcover as explanatory variables. However, R shows me an error message during the formatting of ...
0
votes
1answer
52 views

Calculating Moran's I for buffers of forest cover?

I'm trying to figure out if buffers I have created in ArcMap (measuring percent forest cover in a 1 km distance from a land use map within the buffer) around points at which I will be collecting ...
1
vote
1answer
29 views

Extracting vector/raster data for CSV data points [closed]

I am using QGIS and I have a CSV file of lat and long coordinates which represents presence points of a species in the USA. I also have vector (the state extents/names) as well as raster information (...
1
vote
1answer
28 views

Seeking ArcGIS Pro tool to calculate network distance from point to nearest line?

I am trying to calculate the network distance from a point to the nearest poly line feature as pointed in the image below. Is there any ArcGIS Pro tool that can calculate network distance?
0
votes
0answers
21 views

Understanding Local Moran's I Maps?

I want to better understand the local Moran's I chloropleth maps which show areas (zones) as either in high-high, high-low or low-low clusters. As I assume the "high"/"low" is just relative to the ...
2
votes
0answers
73 views

Rank-based Spatial Clustering in Python

I was wondering if there was any python (or R) packages out there that can do spatial clustering on ranked data? For each US state I spatially varying value for education rankings. For example, for ...
1
vote
2answers
27 views

Specifying which point buffers fall within a particualar land type (polygon)

I have a point buffer layer which represents tick occurrence and a polygon layer which represents landcover types (derived from Corine land cover). I want to define which point buffers (1km) fall ...
1
vote
0answers
22 views

How to create a surface layer from gridded point data of air temperature and Precipitation?

I want to aggregate monthly averages of temperature and precipitation over a 10-year period into weighted seasonal averages (DJF) for use in interpolation. I have downloaded Shapefiles for each month (...
0
votes
1answer
38 views

Measure distance for spatialLine in R returns bad value

I have this spatialLine in R: library(sp) library(maptools) library(SDraw) library(leaflet) x <- c(153.9, 152.67072,147.9879,144.25364,135.45557,116.03998,76.92708,45.10622,37.02522,36.76888,36....
0
votes
0answers
21 views

Using spatial clustering or statistical tools for attribute data

I have a data-set of well test results with attributes like iron and manganese concentrations. What are the best tools to try and determine if there are any spatial trends for each attribute? I would ...
1
vote
1answer
77 views

Plot trajectory in R

Does anyone knows a library/technique in R where you can plot a trajectory using multiple timesteps (points) where each point does have a lat/lon coordinate? I have this data x <- read.table("...
0
votes
0answers
8 views

Adding output value to attribute table of a feature class in model builder (ArcGIS 10.4)

I was doing spatial autocorrelation analysis in ArcGIS 10.4. So is there any way by which I can input the obtained Moran's index value in attribute table of other feature class if the objectID is ...
1
vote
1answer
48 views

Seeking ArcGIS TREND tool equivalent in R?

Arc features the TREND function with linear interpolation up to 12th order. Is there an equivalent approach in R (preferably) or python? I'm trying to de-trend a raster surface using polynomial ...
1
vote
0answers
53 views

Geographic scale of model and covariate significance using ArcGIS Spatial Analyst and R?

A few months back I created a spatial logistic regression model of terrestrial skunk rabies in Arizona/Colorado/New Mexico/Utah. A link to the project can be found here (https://docs.wixstatic.com/ugd/...
1
vote
2answers
41 views

Understanding Neighborhood Majority Statistic in ArcGIS Desktop?

Here both 12 and 13 are equal so would the neighborhood majority be null for this 3 by 3 neighborhood?
0
votes
1answer
66 views

Using raster stack function time in R?

I successfully ran a function over a smaller test raster stack (48 layers) in about 30 minutes but cannot run it over my large stack (1462 layers). The last one I tried over the large stack ran from ...
3
votes
2answers
94 views

Heuristic relationship between number of satellites and GPS precision?

I'm a statistician conducting an analysis on animal movement based on data collected using ATS Wildlink 500 collars. As part of this analysis, I'm estimating GPS error assuming that the error in the ...
1
vote
0answers
60 views

Reducing “noise” in spatially regressed points

I'm working on spatial data on a grid. I have two variables that I think to be correlated and I would like to build a model (see Fig attached). I did some exploratory analyses following the ...
0
votes
1answer
212 views

How kriging variance is calculated in R gstat?

How does exactly the krige function (a wrapper to gstat and predict functions) from package gstat calculate kriging variance (minimized estimation variance) in Ordinary Kriging? I wanted to reproduce ...
0
votes
0answers
251 views

Estimate spatial kernel density from raster and output kernel values in a grid

I would like to use a raster map (find below) to estimate or create a kernel density distribution. The objective is to extrapolate the kernel density function and re-use it in order to distribute a ...
0
votes
0answers
27 views

Method to spatially analyse and compare sets of points with eachother

I have 3 sets of points represented by the image below (3 different shapefiles). Each point has a coordinate. All points represent the same type of real world feature (veterinary surgeries) Using ...
3
votes
1answer
173 views

Standard deviation for percentiles interval in Google Earth Engine

I'm working on a GEE script aiming to extract DEM-based elevation statistics for regions defined by polygons. Specifically, I want to obtain the mean, median, standard deviation, and min/max values of ...
1
vote
1answer
33 views

How to create two Windows in spatstat for the same dataset?

Spatial patterns: I have a dataset of coordinates and marks (diameter) attributed to the points. I am trying to create a window (owin()) for each plot (I have 2 plots), but I don't know how to set the ...
1
vote
0answers
106 views

Error in calculating covariance matrix and fitting variogram model

When constructing variogram model, usually the spatial points are given, and we build a variogram model and conduct parameter estimation. But here in my question the points are not given while ...
0
votes
1answer
44 views

Mismatch in area calculation using R and ArcMap?

I am new to spatial statistics using R. In this scenario, I am trying to calculate the polygons area with in a shapefile by usingrgeos::gArea and compared with ArcMap calculate geometry value. both ...
0
votes
0answers
55 views

Statistical analysis of polygons QGIS : find areas of statistical difference

I have polygons for the country divided into E.D. areas. These polygons illustrate the average annual rainfall in their ED area. Is there any way to statistically tell if one polygon area has a ...
5
votes
1answer
114 views

Cokriging: zero distance cross-semivariance in undersampled cases [gstat error]

Trying to cokrige two variables that are not perfectly colocated (one has sparser measurements than the other), I faced an issue that I'll illustrate with the following MRE. With the meuse dataset, ...
0
votes
1answer
37 views

Install GRASS GIS extension r.pi including r.pi.prox and implement the Addon in QGIS

I would like to install the GRASS GIS addon r.pi within QGIS, I only have access to basic GRASS commands I installed QGIS using the osgeo4w-setup-x86_64.exe advanced installation using the http://...
2
votes
1answer
69 views

How can / should sf spatial objects be translated (i.e. moved within same projection)?

Big fan of the sf (and previously, sp) spatial packages in R. I see spTransform and st_transform as methods to reproject a spatial object into a new space. However, I have a case where I have a ...
2
votes
0answers
68 views

Any advantage to use elongated focal window for raster Pearson's Correlation?

I am performing Pearson's Correlation with 2 rasters (same res, size, extent), and will use a focal window, but is there any advantage or robustness to using a focal window that better reflects the ...
2
votes
1answer
165 views

What is a relevant geospatial modelling approach when no spatial autocorrelation is found?

Although I have found on StackOverflow all answers I ever had on GIS matters using R, ArcGIS and QGIS in the past, I am currently faced with one question I can't find an answer to online. For my ...
1
vote
1answer
136 views

Moran's I: Selecting the best option for conceptualizing the spatial relationship between # births & county locations

I'm using Moran's I to see whether there are clusters of births rates present among the 83 Michigan counties, and, if so, whether the clustering is statistically significant. The inputs are the ...
0
votes
0answers
24 views

Subsampling data to reduce autocorrelation/pseudo-replication

I am interested in the effect subsampling data has on reducing autocorrelation/pseudo-replication. My application is biologically-inspired: I am examining trends in genetic variation within species ...
1
vote
1answer
76 views

LecoS batch processing - Is there an easier way?

I'm using the LecoS software: https://github.com/Martin-Jung/LecoS I've clipped my raster to reduce processing time to focus in on the areas I'm interested in. The entire raster has 11 classes, but ...
1
vote
0answers
59 views

How to calculate sum of values in a point shapefile within a moving window

I have a point shapefile that has a 'rank' attribute with numbers between 0-100. I need to calculate the sum of this value for points within a 10km circular moving window anywhere within my area of ...
1
vote
0answers
109 views

Calculating Moran's I on large point dataframe in R

I have a dataset of about 77,000 points of which I want to test autocorrelation. I have tried to compute the test statistic using the testSpatialAutocorrelation tool in the DHARMa package, which uses "...
0
votes
0answers
70 views

Points to a 1km x 1km cell grid

I have a large database of point data throughout all of Eastern Australia. Each point represents species presence/absence data and has attached about 18 predictor variables. To account for sampling ...
0
votes
1answer
82 views

correlation matrix and p-value in grass

I'd like to know if there is a way to obtain statistical significance and p-value for correlation coefficients calculated for different raster in grass with the command r.covar All that I obtain is ...
2
votes
1answer
167 views

Is there a statistical test to calculate difference/similarity between two sets of lat/lng data points?

I have two sets of data points where each data point is a latitude-longitude combination. The two sets represent two different time frames, let's say 2010-2014 and 2015-2018. The number of data points ...
0
votes
0answers
25 views

Visualizing two different datasets amicably when they are measuring two different things

I am working on a story map in AGOL that is comparing the percentage of a population that lives in a house that has been deemed transit dependent versus the concentration of a population that lives in ...
0
votes
0answers
13 views

test.W ade4 package FUNCTION NOT FOUND

I am trying to perform a MEM analyses, Morans eigenvectors spatial analyses, test.W function Selects among the MEM variables constructed on Delaunay Triangulation I am trying to use the function test....
0
votes
1answer
91 views

How do you find the appropriate spatial weight to use for a spatial analysis? [closed]

I am currently about to run a spatial regression (OLS, Lag, Error etc) for my dataset. My dataset contains 1000 randomly generated points throughout the utilisation distribution for my animals (so, ...
1
vote
1answer
260 views

How to extract specific values with point coordinates from Kriging interpolations made in R?

By using R version 3.4.2 and the library "geoR", I made kriging interpolations for different variables (bellow I give an example of my process). I also made a matrix with the coordinates for 305 trees ...
0
votes
1answer
73 views

Is there a tool in ArcToolbox that can compare types of shapefiles?

I have a some point shapefiles that represent categories of folklore story, witch, mermaid etc. I also have some point shapefiles of landmarks like portal dolmens, ancient woods etc. Is there a tool ...
0
votes
0answers
47 views

How to assess significant land cover changes using Chi-square test?

I have a series of land use maps at different times as a result of a classification process on Landsat images. I am mapping 3 different thematic classes : bare soils, vegetation and croplands. I ...