Mam listę buforów przestrzennych (30000 buforów), które zbudowałem z funkcją 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 : 307941.8,...
Migruję kod z pakietu sp do nowszego pakietu sf. Mój poprzedni kod miałem wielokąta SpatialDataFrame (censimentoMap) i SpatialPointDataFrame (indirizzi.sp) i otrzymałem identyfikator komórki wielokąta („Cell110”) dla każdego punktu leżącego w obrębie poniższej instrukcji: points.data <-...
Pracuję nad projektem epidemiologii środowiskowej, w którym mam narażenia punktowe (~ 2000 operacji wieprzowych - IHO). Te IHO rozpylają na pobliskie pola, ale krople wody i zapach kału mogą pokonywać kilometry. Te ekspozycje punktowe otrzymują 3 miliony buforów i chcę poznać liczbę ekspozycji IHO...
Mam dwie klasy dzielące ten sam CRS (Latitit i Longitude): bolognaQuartieriMap: SpatialPolygonDataFramezawierające dane gmin miejskich. crashPoints: SpatialPointsDataFramezawierający dane o wypadkach. Są dobrze wykreślone przy użyciu: plot(bolognaQuartieriMap) title("Crash per...
Chciałbym wykonać nakładkę przestrzenną, aby zidentyfikować wielokąt, w którym mieści się zbiór punktów. Najpierw jednak chcę buforować i rozpuszczać wielokąty, tak aby punkty mieszczące się w dowolnym miejscu w połączonych wielokątach (ale nie w żadnych dziurach) zostaną podobnie oznaczone za...