sp 패키지에서 최신 sf 패키지로 코드를 마이그레이션하고 있습니다. 내 이전 코드에는 다각형 SpatialDataFrame (censimentoMap)과 SpatialPointDataFrame (indirizzi.sp)이 있었고 아래 지침에 따라 각 점에 대한 다각형 셀 ID ( "Cell110")를 얻었습니다.
points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"])
실제로 두 개의 sf 객체를 만들었습니다.
shape_sf <- st_read(dsn = shape_dsn)
shape_sf <- st_transform(x=shape_sf, crs=crs_string)
과
indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string)
그리고 위의 명령과 동등한 sf를 찾고 있습니다 ... Migth it :
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"]