MODIStsp export RasterStack without CRS

[zhihaojin]

I updated MODIStsp package today and tried to process MOD13A3 NDVI products with the default output projection. After the processing was done, the CRS of the rasterstack was NA. The warning was: "In .newCRS(value) : +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +R=6371007.181 +units=m +no_defs is not a valid PROJ.4 CRS string In .gd_SetProject(object, ...) : NOT UPDATED FOR PROJ >= 6""

I guess the point is about package "PORJ", so I downloaded this package and it didn't help.

There are also some other problems like I cannot change the folder paths with the button "browse". I can provide more information if needed.

Thank you for your time Here is the session info:

R version 3.6.3 (2020-02-29) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows >= 8 x64 (build 9200)

Matrix products: default

locale: [1] LC_COLLATE=Chinese (Simplified)_China.936 LC_CTYPE=Chinese (Simplified)_China.936
[3] LC_MONETARY=Chinese (Simplified)_China.936 LC_NUMERIC=C
[5] LC_TIME=Chinese (Simplified)_China.936

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

other attached packages: [1] gWidgetsRGtk2_0.0-86 cairoDevice_2.28 gWidgets_0.0-54.1 RGtk2_2.20.36
[5] MODIStsp_1.4.0 lubridate_1.7.8 tmaptools_3.0 tmap_3.0
[9] rgdal_1.5-8 rasterVis_0.47 latticeExtra_0.6-29 lattice_0.20-38
[13] raster_3.1-5 sp_1.3-1

loaded via a namespace (and not attached): [1] Rcpp_1.0.1 png_0.1-7 class_7.3-15 zoo_1.8-5 digest_0.6.20
[6] mime_0.7 R6_2.4.0 stats4_3.6.3 e1071_1.7-1 httr_1.4.1
[11] rlang_0.4.4 data.table_1.12.8 hexbin_1.28.1 webshot_0.5.2 stringr_1.4.0
[16] htmlwidgets_1.5.1 munsell_0.5.0 shiny_1.4.0 compiler_3.6.3 httpuv_1.5.2
[21] base64enc_0.1-3 htmltools_0.4.0 gdalUtilities_1.1.0 mapedit_0.6.0 codetools_0.2-16
[26] mapview_2.7.0 XML_3.99-0.3 viridisLite_0.3.0 later_1.0.0 sf_0.9-3
[31] bitops_1.0-6 grid_3.6.3 jsonlite_1.6 lwgeom_0.2-1 satellite_1.0.2
[36] xtable_1.8-4 DBI_1.1.0 magrittr_1.5 units_0.6-5 scales_1.0.0
[41] KernSmooth_2.23-16 stringi_1.4.3 promises_1.1.0 leafsync_0.1.0 leaflet_2.0.3
[46] xml2_1.2.2 xts_0.12-0 generics_0.0.2 RColorBrewer_1.1-2 tools_3.6.3
[51] dichromat_2.0-0 leafem_0.1.1 crosstalk_1.0.0 jpeg_0.1-8.1 abind_1.4-5
[56] parallel_3.6.3 fastmap_1.0.1 yaml_2.2.0 colorspace_1.4-1 stars_0.4-1
[61] classInt_0.4-3

[lbusett]

@zhihaojin

thanks for reporting this. I'll have a look as soon as possible.

Lorenzo

[lbusett]

@zhihaojin

unfortunately at the moment I was not able to replicate the problem.

Could you please save your processing options to JSON from the GUI and cut and paste the content of the file here?

Lorenzo

[lbusett]

Hi. Are you still having problems? Otherwise I'll close this.

[zhihaojin]

Sorry for replying late. I found another way out. Thank you very much. You can close this.