<div dir="auto"><div>You need to also setExtent() of the raster. That's the coordinates that the crs refers to xmin,xmax,ymin,ymax values in metres. Sometimes visual cues can be used to derive these from an image, sometimes the grid is a commonly used convention (a certain dimension, I.e. extent +resolution in a region). It really needs the file source itself or some info from the provider where the data was already a map.</div><div dir="auto"><br></div><div dir="auto">Not much can be done from a picture alone though. </div><div dir="auto"><br></div><div dir="auto">I'll have a look around sometimes wild guesses work because projection-use is rarely bespoke.</div><div dir="auto"><br></div><div dir="auto">Best </div><div dir="auto"><div dir="auto"><br></div><div dir="auto"><br></div><br><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Mon., 27 Apr. 2020, 01:21 Jesús Rojo, <<a href="mailto:jru85yebe@hotmail.com" target="_blank" rel="noreferrer">jru85yebe@hotmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;line-break:after-white-space">Hello everyone<div>I have a flipped RasterBrick without crs assigned (crs = NA):</div><div>plot(a[[100]]</div><div><img id="m_-6738428942250779720m_-6344601436895321837D81E60BC-F239-4012-BC54-8D6FE32FA074" src="cid:6779E725-735B-4AD3-B80F-4EBDC3619C94"></div><div>I would like to fill the map and assign the correct coordinate systems, that it has been said to me to be specifically:</div><div><br></div><div>"+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”</div><div><br></div><div>I try to solve with the next commands for the raster library:</div><div><br></div><div>In a first step I have flipped the map, and it has been right with:</div><div><br></div><div>a <- t(flip(a, direction='y' ))</div><div><br></div><div>But I have problems with the assignments of the suitable crs.</div><div><br></div><div>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”</div><div><br></div><div>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”</div><div><br></div><div>A <- projectRaster(a, crs = "+proj=longlat +datum=WGS84")</div><div><img id="m_-6738428942250779720m_-6344601436895321837F12669EC-B3AA-434A-B2C3-C8EAA86A7FA3" src="cid:361FA390-927A-42D2-8FA8-2826A07D5F65"></div><div>I will be grateful some kind of help. Do I need more additional information? Any step is wrong?</div><div>Thank you so much</div><div>Jesús Rojo</div></div>_______________________________________________<br>
R-sig-Geo mailing list<br>
<a href="mailto:R-sig-Geo@r-project.org" rel="noreferrer noreferrer" target="_blank">R-sig-Geo@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" rel="noreferrer noreferrer noreferrer" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a><br>
</blockquote></div></div></div>