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

Jesús Rojo jru85yebe @end|ng |rom hotm@||@com
Thu Apr 30 00:33:45 CEST 2020


> 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



More information about the R-sig-Geo mailing list