Thanks Roger for your feedback and clarification.
Best regards.
El lun., 27 abr. 2020 a las 5:04, Roger Bivand ()
escribió:
> On Sat, 25 Apr 2020, Jose Ramon Martinez Batlle wrote:
>
> > Dear Anaïs.
> >
> > I am sure more experienced members will give you a better answer, but
> until
> > that
You need to also setExtent() of the raster. That's the coordinates that the
crs refers to xmin,xmax,ymin,ymax values in metres. Sometimes visual cues
can be used to derive these from an image, sometimes the grid is a commonly
used convention (a certain dimension, I.e. extent +resolution in a
> Hello everyone
> I have a flipped RasterBrick without crs assigned (crs = NA):
> plot(a[[100]]
>
> I would like to fill the map and assign the correct coordinate systems, that
> it has been said to me to be specifically:
>
> "+init=EPSG:3034 +proj=lcc +lat_1=35 +lat_2=65 +lat_0=50 +lon_0=10
Similar to sp and rgdal, upcoming sf CRAN windows binaries will use GDAL
3.0.4, PROJ 6.3.1, and GEOS 3.8.0. You can try a beta-version of sf
0.9-3 now by installing it with:
install.packages("sf", repos = "https://edzer.github.io/drat;, type =
"win.binary")
Comments are welcome here, or as
Roger and Mathias,
I tried the latest version (rev 962) and it works! Perfect!
Thanks a lot!
Cheers,
Thorsten
29.04.2020 11:27:18 Roger Bivand :
> Thanks, as I thought. Committed to R-forge as SVN rev. 962. Really we'd need
> a test rig for reading, manipulation in sp/rgdal and raster, and
Thanks, as I thought. Committed to R-forge as SVN rev. 962. Really we'd need a
test rig for reading, manipulation in sp/rgdal and raster, and writing, to be
sure that all int assumptions that need conversion to R_xlen_t.
Roger
--
Roger Bivand
Norwegian School of Economics
Helleveien 30, 5045
