Добавьте значения к полигонам в SpatialPolygonsDataFrame для универсального кригинга
У меня есть кусок кода, который определяет пустую сетку для обычного прогноза кригинга:
library(automap)
library(rgdal)
library(raster)
library(rgeos)
library(dismo)
CP <- as(extent(649500,651000,5591000,5592000), "SpatialPolygons") # Cropping area
proj4string(CP) <- CRS(proj4string(LosRios)) # Reproject
out <- gIntersection(LosRios, CP, byid=TRUE) # Intersect with Los Rios Region Shapefile
grid <- raster(extent(out)) # Create the grid with a raster
res(grid) <- 25 # Grid resolution in meters
proj4string(grid)<-proj4string(out) # Reproject the grid
gridpolygon <- rasterToPolygons(grid) # Turn to polygon
dry.grid <- intersect(out, gridpolygon) # Intersect with crop
Теперь я понимаю, что dry.grid это SpatialPolygonsDataFrame, но как я могу назначить значения (в нашем случае, STREET а также LENGTH) каждому полигону?
Прошу прощения за то, что я не публикую данные, но для воспроизведения примера вам понадобится шейп-файл. Может кто-нибудь попытаться ответить на собственном опыте?
Заранее благодарим и, пожалуйста, попросите любую дополнительную информацию, которая может оказаться необходимой для ответа на вопрос.