Cómo crear una grilla de puntos espaciales

grid geospatial

Question

library(reshape2)
library(data.table)
library(dplyr)
library(magrittr)
library(ggplot2)
library(scales)
library(gstat)
library(DescTools)
library(sp)



#I want a colorado grid#
data("colorado.grid")

#making cordinates into spatial points dataframe#
coordinates(Gold_tracer_kri) <- ~ long_orig + lat_orig
#attempt at kriging but no grid#
lzn.kriged <- krige(Au ~ 1, Gold_tracer_kri, colorado.grid, model=lzn.fit)

lzn.kriged %>% as.data.frame %>%
ggplot(aes(long_orig=long_orig, lat_orig=lat_orig)) +     geom_tile(aes(fill=var1.pred)) + coord_equal() +
scale_fill_gradient(low = "yellow", high="red") +
scale_x_continuous(labels=comma) + scale_y_continuous(labels=comma) +
theme_bw()

cargar el dominio espacial para interpolar sobre

data("meuse.grid")

Intento usar métodos kriging en R pero estoy atascado debido a que no puedo encontrar una grilla para mis datos. Mis datos abarcan todo el estado estadounidense de Colorado y me gustaría poder obtener una cuadrícula para relacionar mis datos. Similar al meuse.grid que se usa en el ejemplo que estoy siguiendo.

Cualquier ayuda sería apreciada

Answer 1

puedes crear una grilla usando sp::makegrid

library(sp)
library(rgdal)
library(raster)

# load some spatial data. Administrative Boundary
us <- getData(''GADM'', country = ''US'', level = 1)
us$NAME_1
colorado <- us[us$NAME_1 == "Colorado",]


# check the CRS to know which map units are used
proj4string(colorado)
# "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

# Create a grid of points within the bbox of the SpatialPolygonsDataFrame 
# colorado with decimal degrees as map units
grid <- makegrid(colorado, cellsize = 0.1)

# grid is a data.frame. To change it to a spatial data set we have to
grid <- SpatialPoints(grid, proj4string = CRS(proj4string(colorado)))

plot(colorado)
plot(grid, pch = ".", add = T)

Y otro ejemplo con Austria (código GADM ''AUT'' ).

Puntos de rejilla solo dentro del polígono:

Para extraer solo los puntos dentro del polígono, use `[` para subdividir los puntos según la ubicación de esta manera:

grid <- grid[colorado, ]