У меня есть список пространственных буферов (30000 буферов), которые я построил с помощью функции lapply: buff.pts <- lapply(1:nrow(pts.prj), FUN=function(l){ buff <- gBuffer(pts.prj[l,], width=1000) ## 1km return(buff) })) > head(buff.pts) [[1]] class : SpatialPolygons features : 1 extent...
Я перенес код из пакета sp в новый пакет sf. В моем предыдущем коде у меня был полигон SpatialDataFrame (censimentoMap) и SpatialPointDataFrame (indirizzi.sp), и я получил идентификатор ячейки полигона («Cell110») для каждой точки, лежащей внутри, с помощью инструкции ниже: points.data <-...
Я работаю над проектом по экологической эпидемиологии, где у меня есть точечные воздействия (~ 2000 операций по свиноводству - МГО). Эти МГО распыляют на близлежащие поля, но капли воды и запах фекалий могут преодолевать мили. Таким образом, эти точечные экспозиции получают 3-миллиметровые буферы,...
У меня есть два класса, использующих один и тот же CRS (Latitute и Longitude): bolognaQuartieriMap: SpatialPolygonDataFrameсодержащие данные городских районов. crashPoints: SpatialPointsDataFrameсодержащие данные о несчастных случаях. Они хорошо построены с использованием: plot(bolognaQuartieriMap)...
Я хотел бы выполнить пространственное наложение, чтобы идентифицировать полигон, в который попадает набор точек. Однако сначала я хочу буферизовать и растворить многоугольники, чтобы точки, падающие где-либо внутри объединенных многоугольников (но не внутри каких-либо отверстий), были помечены...