Message-ID: <VI1P194MB0558FAAE37D58CE5FB3961F7B4AD0@VI1P194MB0558.EURP194.PROD.OUTLOOK.COM>
Date: 2020-04-29T22:33:45Z
From: Jesús Rojo
Subject: Help: assign crs to RasterBrick
In-Reply-To: <VI1P194MB0558E097A213BB0E1680DA94B4AE0@VI1P194MB0558.EURP194.PROD.OUTLOOK.COM>
> Hello everyone
> I have a flipped RasterBrick without crs assigned (crs = NA):
> plot(a[[100]]
> <PastedGraphic-2.tiff>
> 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")
> <PastedGraphic-1.tiff>
> 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