Come posso convertire i punti lon-lat in funzioni semplici (sfg) e poi inserirli in una raccolta di funzioni semplici (sfc)?
Ecco un MWE che non funziona ma è il più vicino a cui sono arrivato.
library(data.table)
library(sf)
# The DT data.table is the data I have (but 10,000s of rows, each row is a point)
DT <- data.table(
place=c("Finland", "Canada", "Tanzania", "Bolivia", "France"),
longitude=c(27.472918, -90.476303, 34.679950, -65.691146, 4.533465),
latitude=c(63.293001, 54.239631, -2.855123, -13.795272, 48.603949),
crs="+proj=longlat +datum=WGS84")
DT[, rowid:=1:.N]
# The following two rows do not work
DT[, place.sfg:=st_point(x=c(longitude, latitude), dim="XY"), by=rowid]
places.sfc <- st_sfc(DT[, place.sfg], crs=DT[, crs])
# This should result in five points, which it doesn't
plot(places.sfc)
Sto cercando di apprendere le funzionalità semplici (motivo per cui non voglio usare la libreria sp), e in seguito ho bisogno di eseguire st_buffer su sfc.
Forse è meglio creare lo sfc direttamente, senza sfg per punto?
Uso data.table per motivi di velocità (10.000 di migliaia di punti che vengono analizzati anche senza aspetti geografici).
Penso di aver bisogno di uno sfc di sfg-points e non di un MULTIPOINT-sfg.