rspatial / terra

R package for spatial data handling https://rspatial.github.io/terra/reference/terra-package.html
GNU General Public License v3.0
536 stars 89 forks source link

Let plet() work with "local" CRS #1570

Closed aloboa closed 2 months ago

aloboa commented 2 months ago

Given

> s <- rast(system.file("ex/logo.tif", package="terra"))
> crs(s)
[1] "ENGCRS[\"Cartesian (Meter)\",\n    EDATUM[\"Unknown engineering datum\"],\n    CS[Cartesian,2],\n        AXIS[\"(E)\",east,\n            ORDER[1],\n            LENGTHUNIT[\"metre\",1,\n                ID[\"EPSG\",9001]]],\n        AXIS[\"(N)\",north,\n            ORDER[2],\n            LENGTHUNIT[\"metre\",1,\n                ID[\"EPSG\",9001]]]]"

I get:

> plet(s[[1]])
Error in h(simpleError(msg, call)) : 
  error in evaluating the argument 'x' in selecting a method for function 'ext': error in evaluating the argument 'x' in selecting a method for function 'project': [project] Cannot do this transformation
In addition: Warning message:
Cannot find coordinate operations from `ENGCRS["Cartesian (Meter)",EDATUM["Unknown engineering datum"],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]' to `+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 +type=crs' (GDAL error 6) 

I would greatly appreciate the zooming capabilities of plet() for objects with arbitrary ("local") coordinate systems.

rhijmans commented 2 months ago

Done.

aloboa commented 2 months ago

I confirm it works, thianks