Я хотел бы объединить полигональные и растровые данные в один фрейм данных, чтобы затем использовать пакет randomForests в R.
Это включает в себя сначала извлечение среднего значения растра на полигон.

Пока у меня есть следующее:

#load libraries
library(raster)
library(rgdal)
library(sp)
library(maptools)

#import raster data 
r <- raster("myRasterdata.tif")

#import polygon data 
p <- readShapePoly("myPolydata.shp")

#extract mean raster value for each polygon
ExtractMyData <- extract(r, p, small=TRUE, fun=mean, na.rm=TRUE, df=FALSE,  nl=1, sp=TRUE)
# note I have also tried this with df=TRUE and sp=FALSE

Выходные данные представляют собой матрицу, которую я могу записать в фрейм данных. Но у него нет пространственных координат или исходных идентификаторов полигонов, поэтому я не знаю, как объединить вывод в одну и ту же базу данных. Я думал, что аргумент sp = TRUE сделает это, но, похоже, он не работает.

Обратите внимание, что на самом деле мне нужно будет преобразовать полигоны в точки (используя метод центроида?) Для целей RandomForests, чтобы я мог догадаться, что я действительно хочу, это соединить средние значения растра, соединенные с точками, а не с полигонами.

Любые предложения будут ценны. Спасибо!!