Я пытаюсь сравнить установленный растр растительного покрова с тем, который я пытаюсь создать из полигональных данных в SpatialPolygonsDataFrame. Я создал пустой растр с тем же размером сетки и т. Д. И хотел бы итеративно извлекать значения из полигонов, лежащих над каждой растровой ячейкой. Для некоторых ячеек я хотел бы извлечь данные из многоугольников, непосредственно перекрывающих эту ячейку, а для некоторых ячеек я бы хотел извлечь данные из любого многоугольника, лежащего над этой ячейкой и его 8 или 24 соседними соседями. Я могу предоставить немного псевдокода, чтобы проиллюстрировать, что я пытаюсь сделать. В пакете R Raster есть много функций, которые используют границы полигонов для извлечения растровых значений, но я пытаюсь сделать наоборот.
for cell in Raster {
layer1 <- extract(cell, SPDF_data$year_built)
layer2 <- extract(cell, SPDF_data$building_type)
layer3 <- extract(cell_8_adjacent_neighbors, number_of_buildings)
}