Nikos wrote:

>> Import using r.in.gdal, _without_ any of `-l` or `-a` and then I get the
>> closest to the reported spatial resolution. Else, with `-a`, for
>> example, the spatial resolution is not as close to the "original" one.
>> Makes sense?

Markus M:

>what would be the output resolution and extends with r.in.gdal -a?
>Generally r.in.gdal -a provides the best results.

Importing without `-a`:

import
360 degree EW extent is exceeded by 0.00019226 cells
north=80.0223214291667
south=-79.9776823855556
east=179.977687153889
west=-180.022321429167
nsres=0.0446428582072328
ewres=0.0446428582072241
rows=3584
cols=8064
cells=28901376
datatype=CELL
ncats=0

Importing with `-a`:

import_with_a
north=80.0223333333333
south=-79.9634444444444

Markus M:

something went wrong here:
east=179.945666666667
east should be 180 - ew_res / 2

just like
west=-180.022333333333
i.e west = -180 - ew_res / 2

I will investigate

Reported here https://trac.osgeo.org/grass/ticket/3672#ticket.

nsres=0.0446388888888889
ewres=0.0446388888888889
rows=3584
cols=8064
cells=28901376
datatype=CELL
ncats=0

Nikos

Attachment: signature.asc
Description: PGP signature

_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to