> 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
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo