Извлечение значений из растровых полигонов

Question

Я пытаюсь найти среднесуточную температуру для округов в Южной Дакоте по сеткам растров (файлы 'bil'), найденным в http://prism.oregonstate.edu/. Я получаю границы округов из пакета 'maps'.

library(maps)
library(raster)
sd_counties <- map('county','south dakota')
sd_raster <- raster('file_path')

Как извлечь ячейки сетки в каждом округе?Я думаю, что для этого нужно превратить каждый округ в отдельный полигон, но как?Тогда я должен быть в состоянии сделать что-то вроде следующего.Любая помощь будет принята с благодарностью.

values <- extract(raster, list of polygons)  
polygon_means <- unlist(lapply(values, FUN=mean))

Answer 1

Я не знаком с пакетом maps или функцией map, но похоже, что он предназначен исключительно для визуализации, а не для геопространственных операций.

Хотя может быть способ конвертироватьmap объект к фактическим многоугольникам, вот простой способ спеть raster * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 100 * * * * * * * * * * * * 1.Вы можете извлечь пиксели для каждого округа, используя extract(r,sd_counties), где r - желаемый растр.

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