search

USEPA / elevatr

An R package for accessing elevation data
Other
203 stars 26 forks source link

Error with get_elev_raster #76

Closed ACalleros closed 1 year ago

ACalleros commented 2 years ago

I was trying to replicate this repository and when I try to download the location data it gives me the following error:

Accessing raster elevation [=========================] 100% Mosaicing & Projecting Error in sf::gdal_utils(util = "warp", source = files, destination = destfile, : gdal_utils warp: an error occured In addition: Warning messages: 1: In CPL_gdalwarp(source, destination, options, oo, doo, quiet, "-overwrite" %in% : GDAL Error 1: Error: Failed to build program executable! Build Log: Error in processing command line: Don't understand command line argument "inff"! 2: In CPL_gdalwarp(source, destination, options, oo, doo, quiet, "-overwrite" %in% : GDAL Error 1: Error at file /build/gdal/src/gdal-3.5.1/alg/gdalwarpkernel_opencl.cpp line 2479: CL_BUILD_PROGRAM_FAILURE 3: In CPL_gdalwarp(source, destination, options, oo, doo, quiet, "-overwrite" %in% : GDAL Error 1: OpenCL routines reported failure (-11) on line 4695.

What can I do?

jhollist commented 2 years ago

@ACalleros Hard to say without more details. It might be something to do with your local version of GDAL?

Can you provide the details of your sessionInfo() and sf::sf_extSoftVersion()?

ACalleros commented 2 years ago

Sorry, I forgot to include my session info

> sessionInfo()
R version 4.2.1 (2022-06-23)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Manjaro Linux

Matrix products: default
BLAS:   /usr/lib/libopenblasp-r0.3.21.so
LAPACK: /usr/lib/liblapack.so.3.10.1

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=es_MX.UTF-8       
 [4] LC_COLLATE=en_US.UTF-8     LC_MONETARY=es_MX.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                  LC_ADDRESS=C              
[10] LC_TELEPHONE=C             LC_MEASUREMENT=es_MX.UTF-8 LC_IDENTIFICATION=C       

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

other attached packages:
 [1] forcats_0.5.2      stringr_1.4.0      dplyr_1.0.10       purrr_0.3.5        readr_2.1.3       
 [6] tidyr_1.2.1        tibble_3.1.8       ggplot2_3.4.0      tidyverse_1.3.2    terra_1.6-17      
[11] showtext_0.9-5     showtextdb_3.0     sysfonts_0.8.8     metR_0.13.0        MexBrewer_0.0.1   
[16] elevatr_0.4.3.9999

loaded via a namespace (and not attached):
 [1] fs_1.5.2            sf_1.0-9            usethis_2.1.6       lubridate_1.9.0     devtools_2.4.4     
 [6] httr_1.4.4          tools_4.2.1         profvis_0.3.7       backports_1.4.1     utf8_1.2.2         
[11] R6_2.5.1            KernSmooth_2.23-20  DBI_1.1.3           colorspace_2.0-3    withr_2.5.0        
[16] urlchecker_1.0.1    sp_1.5-1            tidyselect_1.2.0    prettyunits_1.1.1   processx_3.8.0     
[21] compiler_4.2.1      progressr_0.11.0    rvest_1.0.3         cli_3.4.1           xml2_1.3.3         
[26] scales_1.2.1        checkmate_2.1.0     classInt_0.4-8      callr_3.7.3         proxy_0.4-27       
[31] digest_0.6.30       pkgconfig_2.0.3     htmltools_0.5.3     sessioninfo_1.2.2   dbplyr_2.2.1       
[36] fastmap_1.1.0       htmlwidgets_1.5.4   rlang_1.0.6         readxl_1.4.1        rstudioapi_0.14    
[41] shiny_1.7.3         generics_0.1.3      jsonlite_1.8.3      googlesheets4_1.0.1 magrittr_2.0.3     
[46] Rcpp_1.0.9          munsell_0.5.0       fansi_1.0.3         lifecycle_1.0.3     stringi_1.7.8      
[51] pkgbuild_1.3.1      grid_4.2.1          promises_1.2.0.1    crayon_1.5.2        miniUI_0.1.1.1     
[56] lattice_0.20-45     haven_2.5.1         hms_1.1.2           ps_1.7.2            pillar_1.8.1       
[61] codetools_0.2-18    pkgload_1.3.1       reprex_2.0.2        glue_1.6.2          data.table_1.14.4  
[66] remotes_2.4.2       modelr_0.1.9        vctrs_0.5.0         tzdb_0.3.0          httpuv_1.6.6       
[71] cellranger_1.1.0    gtable_0.3.1        assertthat_0.2.1    cachem_1.0.6        mime_0.12          
[76] xtable_1.8-4        broom_1.0.1         e1071_1.7-12        later_1.3.0         class_7.3-20       
[81] googledrive_2.0.0   gargle_1.2.1        memoise_2.0.1       units_0.8-0         timechange_0.1.1   
[86] ellipsis_0.3.2     
> sf::sf_extSoftVersion()
          GEOS           GDAL         proj.4 GDAL_with_GEOS     USE_PROJ_H           PROJ 
      "3.11.0"        "3.5.1"        "9.0.1"         "true"         "true"        "9.0.1"
jhollist commented 2 years ago

@ACalleros I am unable to recreate the error, I am on slightly different versions of GDAL and PROJ. One quick question, do you possibly have other versions of GDAL installed on your machines, perhaps via pip or from source? I am getting out of my depth a bit here but I do vaguely recall others having issues with multiple or failed GDAL installs.

alexkrohn commented 1 year ago

Connecting this issue back to #62 as the errors are the same. I'm having this same problem in Ubuntu 22.04, whether I use sf 1.0.9 (current) or 1.0.6 as suggested in #62.

My sf-associated packages are:

> sf_extSoftVersion()
          GEOS           GDAL         proj.4 GDAL_with_GEOS     USE_PROJ_H           PROJ 
      "3.10.2"        "3.4.1"        "8.2.1"         "true"         "true"        "8.2.1"

Might this be a problem with GDAL versions rather than sf? (Downgrading to sf 1.0.6 via source does not seem to also downgrade the dependencies.)

jhollist commented 1 year ago

I am trying again on my PC with GDAL 3.5.2 and have it running right now on Debian with GDAL 3.2.2 too.

Can you run this with just a single state? Different Zoom level?

alexkrohn commented 1 year ago

Wow! Very odd! When I run:

states <- us_states(states = c("AL", "FL"))

se.elevations <- get_elev_raster(locations = states, z = 10, clip = "locations")

It works without a problem. In fact, the whole plot works if I decrease z to 7 instead of 10. I guess the larger plot failing was a memory problem rather than a GDAL problem. Interesting that R didn't crash, but instead gave the gdalwarp error.

Thanks for helping me figure out the issue.

jhollist commented 1 year ago

Yep! On my Ubuntu machine it failed on the larger ones. And I got different errors depending on how many states I included. On a mid size one it reported back that it needed 20+ GB to store something. Windows and Debian (actually a chromebook) the full thing worked fine. May be a GDAL version thing since that was different across all three.

At a minimum, this at least seems to be a work around that'll at least get you 1/2 way there!

On Tue, Dec 13, 2022 at 11:55 AM alexkrohn @.***> wrote:

Wow! Very odd! When I run:

states <- us_states(states = c("AL", "FL"))

se.elevations <- get_elev_raster(locations = states, z = 10, clip = "locations")

It works without a problem. In fact, the whole plot works if I decrease z to 7 instead of 10. I guess the larger plot failing was a memory problem rather than a GDAL problem. Interesting that R didn't crash, but instead gave the gdalwarp error.

Thanks for helping me figure out the issue.

— Reply to this email directly, view it on GitHub https://github.com/jhollist/elevatr/issues/76#issuecomment-1349057953, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABJPYSZJBFE4AWHSWJJTA73WNCTAZANCNFSM6AAAAAAR45CLUU . You are receiving this because you commented.Message ID: @.***>

-- Jeffrey W. Hollister email: @.*** cell: 401 556 4087 https://jwhollister.com

alexkrohn commented 1 year ago

Interesting! I got it to work with all the states by taking the zoom level down to 8. Thanks again for your help.

jhollist commented 1 year ago

Definitely a memory thing... Feels like we made progress. You OK if I close the issue?

On Tue, Dec 13, 2022 at 4:07 PM alexkrohn @.***> wrote:

Interesting! I got it to work with all the states by taking the zoom level down to 8. Thanks again for your help.

— Reply to this email directly, view it on GitHub https://github.com/jhollist/elevatr/issues/76#issuecomment-1349689614, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABJPYS6UCJOZC2NBLGC2DCTWNDQQHANCNFSM6AAAAAAR45CLUU . You are receiving this because you commented.Message ID: @.***>

-- Jeffrey W. Hollister email: @.*** cell: 401 556 4087 https://jwhollister.com

alexkrohn commented 1 year ago

My issue is resolved. Feel free to close. Thanks!

On Tue, Dec 13, 2022 at 4:43 PM Jeffrey W Hollister < @.***> wrote:

Definitely a memory thing... Feels like we made progress. You OK if I close the issue?

On Tue, Dec 13, 2022 at 4:07 PM alexkrohn @.***> wrote:

Interesting! I got it to work with all the states by taking the zoom level down to 8. Thanks again for your help.

— Reply to this email directly, view it on GitHub https://github.com/jhollist/elevatr/issues/76#issuecomment-1349689614, or unsubscribe < https://github.com/notifications/unsubscribe-auth/ABJPYS6UCJOZC2NBLGC2DCTWNDQQHANCNFSM6AAAAAAR45CLUU

. You are receiving this because you commented.Message ID: @.***>

-- Jeffrey W. Hollister email: @.*** cell: 401 556 4087 https://jwhollister.com

— Reply to this email directly, view it on GitHub https://github.com/jhollist/elevatr/issues/76#issuecomment-1349784062, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADAEIHH4QMYE4P3JTCFOBOLWNDUWTANCNFSM6AAAAAAR45CLUU . You are receiving this because you commented.Message ID: @.***>