मैं एसपी पैकेज से नए एसएफ पैकेज के लिए कोड माइग्रेट कर रहा हूं। मेरा पिछला कोड मेरे पास एक बहुभुज SpatialDataFrame (censimentoMap) और एक SpatialPointDataFrame (indirizzi.sp) था और मुझे नीचे दिए गए निर्देश के साथ प्रत्येक बिंदु पर बहुभुज सेल आईडी ("सेल110") मिला:
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 की तलाश कर रहा हूँ ... यह हो:
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"]