rmbl-sdp / rsdp Goto Github PK

## bring in landcover raster
lc <- sdp_get_raster(catalog_id="R3D018") 
## bring in 8 polygons as sf
field_planning = st_read("C:/Users/AmandaHenderson/SynologyDrive/UER_vegmap/output/sampling_areas_draft.gpkg")
## extract landcover data for polygons of interest
fp_lc = sdp_extract_data(raster = lc, locations = field_planning)
#[1] "Re-projecting locations to coordinate system of the raster."
#[1] "Extracting data at 8 locations for 1 raster layers."
#Error: [cbind] nrow does not match

location_df <- data.frame(SiteName=c("Roaring Judy","Gothic","Galena Lake"),
                          Lat=c(38.716995,38.958446,39.021644),
                          Lon=c(-106.853186,-106.988934,-107.072569))
location_sv <- terra::vect(location_df,geom=c("Lon","Lat"),crs="EPSG:4327")
days <- seq(as.Date("2007-10-01"),as.Date("2007-10-31"),by="day")

cl <- parallel::makeCluster(2L)
doParallel::registerDoParallel(cl)

extr_list4 <- foreach::foreach(i=1:length(days),.packages=c("terra","devtools","rSDP")) %dopar% {
  #devtools::load_all() ## During package development.
  tmax4 <- rSDP::sdp_get_raster("R4D004",date_start=days[i],date_end=days[i],verbose=FALSE)
  extr_dat <- rSDP::sdp_extract_data(tmax4,locations_proj,verbose=FALSE)[,3]
  (extr_dat)
}
parallel::stopCluster(cl)

: NULL value as symbol address.