<div dir="ltr"><br><br>On Tue, Sep 25, 2018 at 4:51 PM Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net">nik@nikosalexandris.net</a>> wrote:<br>><br>> * Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>> [2018-09-25 08:24:44 -0300]:<br>><br>> >AFAIU, plate carree (lat long grid, no meters) is deprecated, and you<br>> >should use latlong instead. I had a similar issue with modis ocean color<br>> >products.<br>> ><br>> >Do you at least get the same number of row and columns that is described by<br>> >gdalinfo when you import?<br>><br>> No, it is 512^2 (see below) against<br>> ```<br>> g.regio -p<br>> ..<br>> rows:       3584<br>> cols:       8064<br>> ..<br>> ```<br>><br>> Markus' hint, not *directly* on the netCDF file, rather on one of the<br>> layers:<br>><br>> gdalinfo NETCDF:"<a href="http://g2_BIOPAR_LST_201606220100_GLOBE_GEO_V1.2.nc">g2_BIOPAR_LST_201606220100_GLOBE_GEO_V1.2.nc</a>":LST<br>> Driver: netCDF/Network Common Data Format<br>> Files: <a href="http://g2_BIOPAR_LST_201606220100_GLOBE_GEO_V1.2.nc">g2_BIOPAR_LST_201606220100_GLOBE_GEO_V1.2.nc</a><br>> Size is 8064, 3584<br>> Coordinate System is:<br>> GEOGCS["unknown",<br>>     DATUM["unknown",<br>>         SPHEROID["Spheroid",6378137,298.2572326660156]],<br>>     PRIMEM["Greenwich",0],<br>>     UNIT["degree",0.0174532925199433,<br>>         AUTHORITY["EPSG","9122"]]]<br>> Origin = (-180.022321429103613,80.022321429103613)<br>> Pixel Size = (0.044642858207226,-0.044642858207226)<br><div>[...]<br></div><div>> Corner Coordinates:</div>> Upper Left  (-180.0223214,  80.0223214) (180d 1'20.36"W, 80d 1'20.36"N)<br>> Lower Left  (-180.0223214, -79.9776824) (180d 1'20.36"W, 79d58'39.66"S)<br>> Upper Right ( 179.9776872,  80.0223214) (179d58'39.67"E, 80d 1'20.36"N)<br>> Lower Right ( 179.9776872, -79.9776824) (179d58'39.67"E, 79d58'39.66"S)<br>> Center      (  -0.0223171,   0.0223195) (  0d 1'20.34"W,  0d 1'20.35"N)<br><div><br></div><div>apparently the center of the left-most column is -180, therefore forcing the corner to -180 would shift the grid by half a cell.<br></div><div><br></div>><br>> Import using r.in.gdal, _without_ any of `-l` or `-a` and then I get the<br>> closest to the reported spatial resolution. Else, with `-a`, for<br>> example, the spatial resolution is not as close to the "original" one.<br>> Makes sense?<br><div><br></div><div>what would be the output resolution and extends with r.in.gdal -a? Generally r.in.gdal -a provides the best results.<br></div>><br>> Better to cut off the west side (?):<br>> ```<br>> g.region raster=g2_BIOPAR_LST_201606220100_GLOBE_GEO_V1.2.nc_LST -pag w=-180 e=180<br><div><br></div><div>this shifts the grid by half a cell to the east. The -a flag does not make sense because 1) you want to force new extends, 2) there is no resolution given for use with the -a flag.</div><div><br></div><div>Using g.region -pg defaults to g.region -g. Use either -p or -g.<br></div><div><br></div><div>Markus M<br></div></div>