jhollist
commented
3 years ago I am not on my machine right now so can't look too closely at this, but I noticed you are on an old version of elevatr. Current version on CRAN is 0.3.4. can you update this and confirm that you are still having this issue?
On Fri, Feb 26, 2021, 7:05 AM James P. Duffy notifications@github.com wrote:
Hello,
I'm afraid this is another issue likely related to the recent large changes to PROJ & GDAL. I get the following error when running:
elevation <- elevatr::get_elev_raster(lake,z = 9) Error in h(simpleError(msg, call)) : error in evaluating the argument 'CRSobj' in selecting a method for function 'spTransform': NA
I've narrowed it down to this line of code:
https://github.com/jhollist/elevatr/blob/4c45d548231400284f72b92b42dff6d9fd9a0928/R/internal.R#L196
...and found that changing it to:
locations <- sp::spTransform(locations, sp::CRS("+proj=longlat +datum=WGS84 +no_defs"))
...resolves the issue.
I think this shouldn't be a breaking change and wondered whether you might consider implementing it?
Here is my sessioninfo / rgdal loading text if useful:
Platform: x86_64-w64-mingw32/x64 (64-bit)Running under: Windows 10 x64 (build 18363) Matrix products: default locale: [1] LC_COLLATE=English_United Kingdom.1252 [2] LC_CTYPE=English_United Kingdom.1252 [3] LC_MONETARY=English_United Kingdom.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United Kingdom.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] rgdal_1.5-23 elevatr_0.3.1 raster_3.4-5 sp_1.4-5 sf_0.9-7 loaded via a namespace (and not attached): [1] Rcpp_1.0.5 magrittr_2.0.1 hms_1.0.0 progress_1.2.2 [5] units_0.6-7 tidyselect_1.1.0 lattice_0.20-41 R6_2.5.0 [9] rlang_0.4.10 httr_1.4.2 dplyr_1.0.4 tools_4.0.3 [13] grid_4.0.3 KernSmooth_2.23-18 e1071_1.7-4 DBI_1.1.1 [17] ellipsis_0.3.1 class_7.3-18 tibble_3.0.6 lifecycle_0.2.0 [21] crayon_1.4.1 purrr_0.3.4 codetools_0.2-18 vctrs_0.3.6 [25] curl_4.3 glue_1.4.2 compiler_4.0.3 pillar_1.4.7 [29] prettyunits_1.1.1 generics_0.1.0 classInt_0.4-3 pkgconfig_2.0.3 ```
when loading rgdal
Loading required package: sprgdal: version: 1.5-23, (SVN revision 1121)Geospatial Data Abstraction Library extensions to R successfully loadedLoaded GDAL runtime: GDAL 3.2.1, released 2020/12/29Path to GDAL shared files: C:/R/R403/library/rgdal/gdalGDAL binary built with GEOS: TRUE Loaded PROJ runtime: Rel. 7.2.1, January 1st, 2021, [PJ_VERSION: 721]Path to PROJ shared files: C:/R/R403/library/rgdal/projPROJ CDN enabled: FALSELinking to sp version:1.4-5To mute warnings of possible GDAL/OSR exportToProj4() degradation,use options("rgdal_show_exportToProj4_warnings"="none") before loading rgdal.Overwritten PROJ_LIB was C:/R/R403/library/rgdal/proj
P.S I think in the much longer term, switching to sf/stars/terra would help solve such issues, but understand that would be a massive undertaking. Oh and thank you for making a very useful package!
— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/jhollist/elevatr/issues/48, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABJPYSZ5GIJT3NFTYDO55JTTA6FCFANCNFSM4YIKODTA .
everydayduffy
arielfuentes
xtimbeau
fipoucat
Hello,
I'm afraid this is another issue likely related to the recent large changes to PROJ & GDAL. I get the following error when running:
I've narrowed it down to this line of code:
https://github.com/jhollist/elevatr/blob/4c45d548231400284f72b92b42dff6d9fd9a0928/R/internal.R#L196
...and found that changing it to:
locations <- sp::spTransform(locations, sp::CRS("+proj=longlat +datum=WGS84 +no_defs"))...resolves the issue.
I think this shouldn't be a breaking change and wondered whether you might consider implementing it?
Here is my sessioninfo / rgdal loading text if useful:
P.S I think in the much longer term, switching to sf/stars/terra would help solve such issues, but understand that would be a massive undertaking. Oh and thank you for making a very useful package!