Mapview problems with CRS, WKT,rgdal with projections

[fabeit]

I am not even sure where to start but mapview is not working anymore for me, I cannot plot anything anymore and I wonder if it's related to the recent changes to PROJ6.

When I try to use it with a raster I get

y <- raster(ncol=180, nrow=75, xmn=-180, xmx=180, ymn=-60, ymx=90)
y=setValues(y,1)
class      : RasterLayer 
dimensions : 75, 180, 13500  (nrow, ncol, ncell)
resolution : 2, 2  (x, y)
extent     : -180, 180, -60, 90  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : layer 
values     : 0.0104, 7.2  (min, max)

Error in rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[,  : 
  target crs creation failed: generic error of unknown origin
In addition: Warning messages:
1: In showSRID(uprojargs, format = "PROJ", multiline = "NO") :
  Discarded ellps WGS 84 in CRS definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs
2: In showSRID(uprojargs, format = "PROJ", multiline = "NO") :
  Discarded datum WGS_1984 in CRS definition
3: In rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[,  :
  48 projected point(s) not finite
4: In rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[,  :
  48 projected point(s) not finite
5: In showSRID(uprojargs, format = "PROJ", multiline = "NO") :
  Discarded ellps WGS 84 in CRS definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs
6: In showSRID(uprojargs, format = "PROJ", multiline = "NO") :
  Discarded datum WGS_1984 in CRS definition
7: In CPL_crs_from_input(x) :
  GDAL Error 4: Unable to open EPSG support file gcs.csv.  Try setting the GDAL_DATA environment variable to point to the directory containing EPSG csv files.
8: In wkt(pto) : CRS object has no comment
9: In rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[,  :
  Using PROJ not WKT2 strings 

When I try with a polygon

t=rasterToPolygons(y)
class       : SpatialPolygonsDataFrame 
features    : 3845 
extent      : -180, 180, -58, 90  (xmin, xmax, ymin, ymax)
crs         : +proj=longlat +datum=WGS84 +no_defs 
variables   : 1
names       :              layer 
min values  : 0.0104000000283122 
max values  :   7.19999980926514 
> mapview(t)

Error in st_crs.character(comment(x)) : invalid crs: GEOGCRS["unknown",
    DATUM["World Geodetic System 1984",
        ELLIPSOID["WGS 84",6378137,298.257223563,
            LENGTHUNIT["metre",1]],
        ID["EPSG",6326]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433],
        ID["EPSG",8901]],
    CS[ellipsoidal,2],
        AXIS["longitude",east,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]],
        AXIS["latitude",north,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]]]

I have just reinstalled all packages

> packageVersion('mapview')
[1] ‘2.9.0’
> packageVersion('raster')
[1] ‘3.3.13’
> packageVersion('sp')
[1] ‘1.4.2’
> packageVersion('rgdal')
[1] ‘1.5.16’

What's happening all of a sudden? I never had these problems.

[tim-salabim]

Are you on windows? What's your sessionInfo()?
I have a feeling the errors you see are related to an outdated GDAL or PROJ version.

[tim-salabim]

Here's what I see (with the commit above):

library(mapview)
#> Warning: replacing previous import 'vctrs::data_frame' by 'tibble::data_frame'
#> when loading 'dplyr'
#> GDAL version >= 3.1.0 | setting mapviewOptions(fgb = TRUE)
library(raster)
#> Loading required package: sp

mapviewOptions(fgb = FALSE)

y <- raster(ncol=180, nrow=75, xmn=-180, xmx=180, ymn=-60, ymx=90)
y = setValues(y, 1)

y
#> class      : RasterLayer 
#> dimensions : 75, 180, 13500  (nrow, ncol, ncell)
#> resolution : 2, 2  (x, y)
#> extent     : -180, 180, -60, 90  (xmin, xmax, ymin, ymax)
#> crs        : +proj=longlat +datum=WGS84 +no_defs 
#> source     : memory
#> names      : layer 
#> values     : 1, 1  (min, max)

mapview(y)
#> Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
#> ellps WGS 84 in CRS definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0
#> +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs
#> Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
#> datum WGS_1984 in CRS definition
#> Warning in rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, : 48
#> projected point(s) not finite

#> Warning in rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, : 48
#> projected point(s) not finite
#> Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
#> ellps WGS 84 in CRS definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0
#> +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs
#> Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO"): Discarded
#> datum WGS_1984 in CRS definition

t = rasterToPolygons(y)

mapview(t)

^{Created on 2020-10-07 by the reprex package (v0.3.0)}

My sessionInfo

R version 4.0.2 (2020-06-22)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Pop!_OS 18.04 LTS

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.7.1
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.7.1

locale:
 [1] LC_CTYPE=en_GB.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_GB.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_GB.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

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

other attached packages:
[1] leaflet_2.0.3 raster_3.3-13 sp_1.4-2      mapview_2.9.4

loaded via a namespace (and not attached):
 [1] xfun_0.18               tidyselect_1.1.0        purrr_0.3.4            
 [4] sf_0.9-6                lattice_0.20-41         colorspace_1.4-1       
 [7] vctrs_0.3.4             generics_0.0.2          leaflet.providers_1.9.0
[10] htmltools_0.5.0         stats4_4.0.2            yaml_2.2.1             
[13] base64enc_0.1-3         rlang_0.4.7             leafpop_0.0.6          
[16] e1071_1.7-3             pillar_1.4.6            glue_1.4.2             
[19] DBI_1.1.0               gdtools_0.2.2           RColorBrewer_1.1-2     
[22] uuid_0.1-4              lifecycle_0.2.0         munsell_0.5.0          
[25] htmlwidgets_1.5.1       evaluate_0.14           codetools_0.2-16       
[28] knitr_1.29              callr_3.4.4             ps_1.3.4               
[31] crosstalk_1.1.0.1       class_7.3-17            leafem_0.1.3           
[34] Rcpp_1.0.5              KernSmooth_2.23-17      clipr_0.7.0            
[37] scales_1.1.1            classInt_0.4-3          satellite_1.0.2        
[40] webshot_0.5.2           jsonlite_1.7.1          farver_2.0.3           
[43] systemfonts_0.2.3       fs_1.5.0                brew_1.0-6             
[46] png_0.1-7               digest_0.6.25           processx_3.4.4         
[49] dplyr_1.0.0             grid_4.0.2              rgdal_1.5-12           
[52] tools_4.0.2             magrittr_1.5            tibble_3.0.3           
[55] whisker_0.4             crayon_1.3.4            pkgconfig_2.0.3        
[58] ellipsis_0.3.1          reprex_0.3.0            rmarkdown_2.3          
[61] svglite_1.2.3.2         rstudioapi_0.11         R6_2.4.1               
[64] units_0.6-7             compiler_4.0.2  

[fabeit]

Thanks Tim, I will reupdate everything again and get back with you.

[fabeit]

After updating R itself and a fresh install of all package mapview works again. Thanks!