Je migre le code du package sp vers le package sf plus récent. Mon code précédent, j'avais un polygone SpatialDataFrame (censimentoMap) et un SpatialPointDataFrame (indirizzi.sp) et j'ai obtenu l'ID de cellule de polygone ("Cell110") pour chaque point se trouvant à l'intérieur avec les instructions ci-dessous:
points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"])
En fait, j'ai créé deux objets sf:
shape_sf <- st_read(dsn = shape_dsn)
shape_sf <- st_transform(x=shape_sf, crs=crs_string)
et
indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string)
Et je cherche l'équivalent sf de l'instruction ci-dessus ... Migth ce soit:
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"]