Saya memigrasi kode dari paket sp ke paket sf yang lebih baru. Kode saya sebelumnya, saya memiliki poligon SpatialDataFrame (censimentoMap) dan SpatialPointDataFrame (indirizzi.sp) dan saya mendapatkan id sel poligon ("Cell110") untuk setiap titik yang diletakkan dengan instruksi di bawah ini:
points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"])
Sebenarnya saya membuat dua objek sf:
shape_sf <- st_read(dsn = shape_dsn)
shape_sf <- st_transform(x=shape_sf, crs=crs_string)
dan
indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string)
Dan saya mencari setara sf dari instruksi di atas ... Migth it be:
ids<-sapply(st_intersects(x=indirizzi_sf,y=shshape_sfpeCrif), function(z) if (length(z)==0) NA_integer_ else z[1])
cell_ids <- shape_sf[ids,"Cell110"]