[R-sig-Geo] Help: assign crs to RasterBrick

Jesús Rojo jru85yebe @end|ng |rom hotm@||@com
Sun Apr 26 17:21:23 CEST 2020


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 +x_0=4000000 +y_0=2800000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs”

I try to solve with the next commands for the raster library:

In a first step I have flipped the map, and it has been right with:

a <- t(flip(a, direction='y' ))

But I have problems with the assignments of the suitable crs.

crs(a) <- "+init=EPSG:3034 +proj=lcc +lat_1=35 +lat_2=65 +lat_0=50 +lon_0=10 +x_0=4000000 +y_0=2800000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs”

Because although the map assigns crs = “+init = EPSG:3034”, the map doesn’t seem georreferenced, because the starting point of the map is (0,0). My final objective is  to transform the crs to "+proj=longlat +datum=WGS84”

A <- projectRaster(a, crs = "+proj=longlat +datum=WGS84")

I will be grateful some kind of help. Do I need more additional information? Any step is wrong?
Thank you so much
Jesús Rojo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20200426/58ab109b/attachment.html>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-2.tiff
Type: image/tiff
Size: 15216 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20200426/58ab109b/attachment.tiff>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-1.tiff
Type: image/tiff
Size: 15528 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20200426/58ab109b/attachment-0001.tiff>


More information about the R-sig-Geo mailing list