MODIStsp failed in to reproject from hdf to WGS84

[cedlfc44]

HI:
Im working into a researh, so i need to change crs from MODIS images during downloads.
So, i working with MODIStsp function, but after finish to create hdf files in my cpu, a fail is occuring. BUT, when, i put #output_proj, download works fine.

thanks for your help

Script
 rm(list=ls(all=TRUE))
 library(MODIStsp)  
 library(bench)
 library(sp)
 library(rgdal)   
 
 MTSP_fechas <- read.csv(file = "./fechas_imagenes_MODIS_tsp.txt")
 class(MTSP_fechas[1,1])
 
 for (ll in 1:24) {
 MODIStsp_get_prodlayers("M*D13Q1")
 spatial_file <- system.file("testdata/lakeshapes/cercado.shp", package = "MODIStsp")
 
 MTSP_descarga<-MODIStsp(gui             = FALSE,
         spatmeth        = "file",
         spafile         = spatial_file,  
         out_folder      = "$tempdir",
         selprod         = "Vegetation Indexes_16Days_250m (M*D13Q1)",
         bandsel         = c("EVI"),
         output_proj = st_as_text(st_crs(4326)),
         # output_proj     = "EPSG:4326",
         # output_proj     = "+proj=longlat +datum=WGS84 +no_defs",
         user            = "********" ,
         password        = "**********",
         start_date      = MTSP_fechas[ll,1],  
         end_date        = MTSP_fechas[ll,1],
         verbose         = FALSE)

 out_fold <- file.path(tempdir(), "MODIStsp/cercado/VI_16Days_250m_v6/") 

FAIL
Error in gdal_utils("buildvrt", gdalfile, output.vrt, opts) :
gdal_utils buildvrt: an error occured
Además: Warning message:
In CPL_gdalbuildvrt(source, destination, options, oo, quiet) :
GDAL Error 5: Invalid user extent

[ranghetti]

Hi, can you attach the shapefile cercado.shp (among with all the accessoriy files) and paste the content of file fechas_imagenes_MODIS_tsp.txt?
Moreover, provide the output of:

sessionInfo()
packageVersion("sen2r")

[cedlfc44]

Hi:
this is fechas_imagenes_MODIS_tsp.txt

Fechas_MODIS
2015.01.01
2015.01.17
2015.02.02
2015.02.18
2015.03.06
2015.03.22
2015.04.07
2015.04.23
2015.05.09
2015.05.25
2015.06.10
2015.06.26
2015.07.12

sessionInfo()
R version 4.0.3 (2020-10-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 18363)

Matrix products: default

locale:
[1] LC_COLLATE=Spanish_Mexico.1252 LC_CTYPE=Spanish_Mexico.1252
[3] LC_MONETARY=Spanish_Mexico.1252 LC_NUMERIC=C
[5] LC_TIME=Spanish_Mexico.1252

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] sen2r_1.3.9 sf_0.9-6 rgdal_1.5-18 sp_1.4-4 bench_1.1.1 MODIStsp_2.0.4

loaded via a namespace (and not attached):
[1] gdalUtilities_1.1.0 zoo_1.8-8 tidyselect_1.1.0 purrr_0.3.4
[5] lattice_0.20-41 V8_3.2.0 vctrs_0.3.4 generics_0.0.2
[9] stars_0.4-4 XML_3.99-0.5 rlang_0.4.8 e1071_1.7-4
[13] pillar_1.4.6 foreign_0.8-80 glue_1.4.2 httpcode_0.3.0
[17] DBI_1.1.0 foreach_1.5.0 lifecycle_0.2.0 stringr_1.4.0.9000
[21] geojsonio_0.9.2 rgeos_0.5-3 raster_3.3-13 RcppTOML_0.1.7
[25] codetools_0.2-16 doParallel_1.0.15 maptools_1.0-1 parallel_4.0.3
[29] curl_4.3 class_7.3-17 xts_0.12-0 Rcpp_1.0.5
[33] KernSmooth_2.23-17 classInt_0.4-3 lwgeom_0.2-5 jsonlite_1.7.1
[37] abind_1.4-5 stringi_1.5.3 dplyr_1.0.0 grid_4.0.3
[41] jqr_1.1.0 tools_4.0.3 bitops_1.0-6 magrittr_1.5
[45] lazyeval_0.2.2 geojson_0.3.4 tibble_3.0.4 crul_1.0.0
[49] crayon_1.3.4 pkgconfig_2.0.3 ellipsis_0.3.1 data.table_1.12.8
[53] xml2_1.3.2 assertthat_0.2.1 httr_1.4.2 rstudioapi_0.11
[57] iterators_1.0.12 R6_2.5.0 units_0.6-7 compiler_4.0.3

packageVersion("sen2r")
[1] ‘1.3.9’

[cedlfc44]

Cercado shape
cercado.zip

[ranghetti]

Sorry for the long waiting. It should be fixed, but I cannon test it since NASA server is currently down. Please install the develop package version with the command below, retry your code and let me know if the bug was fixed.

[ranghetti]

Sorry, I forgot to write the command:

install.packages("remotes")
remotes::install_github("ropensci/MODIStsp", ref = "devel")

[cedlfc44]

ok, thanks im going to try right now

[cedlfc44]

hi:
MODIStsp command is working perfectly. bug is fixed.
thanks for your help