Hi everyone, I'm wondering whether st_extract() from stars package preserves the order of the input sf dataframe. See the example below; is it a safe approach? Is there perhaps a more "elegant" way to do the same extraction?
library(sf) library(stars) library(geodata) library(magrittr) # Create sf dataframe of cities cities <- data.frame(x=c(2.3522, 5.3698, 4.8357), y=c(48.8566, 43.2965, 45.7640), names=c('Paris', 'Marseille', 'Lyon')) sfdf <- st_as_sf(cities, coords=c('x','y'), crs=st_crs(4326), agr='identity') # Get raster of monthly average tempature tavg <- worldclim_tile(var = 'tavg', lon = 15, lat = 45, path = tempdir()) %>% st_as_stars() # Extract monthly avg temperature for cities tavg_df <- st_extract(tavg, at = sfdf) %>% st_as_sf() %>% st_drop_geometry() # Recombine with original sf dataframe full_df <- cbind(sfdf, tavg_df) Kind regards, Loïc _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo