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