<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="auto"><div>Hello Michael,<div dir="auto"><br></div><div dir="auto">Yes, both the grid data and the shapefile after the reprojection are longitude latitude on wgs84.</div><div dir="auto"><br></div><div dir="auto">Thanks for your reply. I'll try to reproject the shape by the qgis, extract and compare the results from one layer.</div><div dir="auto"><br></div><div dir="auto">Cheers.</div><br><div class="gmail_extra"><br><div class="gmail_quote">Em 5 de nov de 2020 23:30, Michael Sumner <mdsumner@gmail.com> escreveu:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">It probably doesn't matter at all. Your data is longitude,latitude on datum WGS84? Those strings are equivalent (some softwares will flip the order of lon,lat as that's what the standard dictates but easy to spot that usually). <div><br></div><div>raster::projection(), raster::crs(), sp::proj4string() all are fairly equivalent functions, they just get the projection metadata - they are just overlapping and inconsistent because of historical vagarie, they differ in ideological stance and how much checking and churn they do to the metadata and in efforts to interpret your intentions. But on the face of it I don't think you need to worry about your code here. </div><div><br></div><div>Just keep an eye on whether your data looks right, before and after reprojection. It takes care, but the details are yours alone for checking unless you can pin down a reproducible example that raises questions. </div><div><br></div><div>Good luck. </div><div><br></div><div>fwiw you seem to not be doing anything with ncdf4 (or ncdf4::nc_open) here. the raster package *uses* ncdf4 on your behalf via raster()/stack()/brick() functions. </div><div><br></div><div>Cheers, Mike. </div><div><br></div><div><br></div></div><br><div class="elided-text"><div dir="ltr">On Fri, Nov 6, 2020 at 12:11 PM Pietro Andre Telatin Paschoalino <<a href="mailto:Pietro_telato@hotmail.com">Pietro_telato@hotmail.com</a>> wrote:<br></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb( 204 , 204 , 204 );padding-left:1ex">Dear,<br>
<br>
I have a question that is possibly extremely easy for you.<br>
<br>
I'm extracting data from layers of grid to a shape file.<br>
<br>
To do this, first, I know that I need to reproject the coordinates reference system of my shape to the same CRS as my grid. The code for this is:<br>
<br>
library(ncdf4)<br>
<br>
fn<-file.path(path of my data)<br>
<br>
ncfile <- nc_open(fn)<br>
<br>
rasbrick <- stack(fn)<br>
<br>
rasbrick <- brick[[70:81]] (the layers that I want)<br>
<br>
shape2<- spTransform(shape, crs(rasbrick)) in lowercase, right?  It is my first question.<br>
<br>
The rest of the extraction code works without any problems. If I open this file in qgis the CRS is:<br>
<br>
OGC:CRS84 - WGS 84 (CRS84) - Geographic<br>
<br>
The problem is that if I write a raster of a layer using:<br>
<br>
rasbrick <- brick[[70]]<br>
<br>
writeRaster(rasbrick, "test.tif", format="GTiff")<br>
<br>
and open the raster in qgis the CRS is:<br>
<br>
EPSG:4326 - WGS 84 - Geographic<br>
<br>
So, my second question is why on the raster it looks like EPSG and on the shape it looks like OGC? And could this cause any problems in extracting data from layers for the shape? If yes, what I need to change to reproject the shape for the correct CRS?<br>
<br>
Thanks.<br>
<br>
Pietro Andre Telatin Paschoalino<br>
Doutorando em Ciências Econômicas da Universidade Estadual de Maringá - PCE.<br>
<br>
        [[alternative HTML version deleted]]<br>
<br>
_______________________________________________<br>
R-sig-Geo mailing list<br>
<a href="mailto:R-sig-Geo@r-project.org">R-sig-Geo@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr">Michael Sumner<br>Software and Database Engineer<br>Australian Antarctic Division<br>Hobart, Australia<br>e-mail: <a href="mailto:mdsumner@gmail.com">mdsumner@gmail.com</a></div>
</blockquote></div><br></div></div></div>