[GRASS-user] some questions on r.tile

Margherita Di Leo diregola at gmail.com
Thu Aug 11 01:09:13 PDT 2016


Ciao Markus! Thanks for your answer,

On Wed, Aug 10, 2016 at 8:00 PM, Markus Neteler <neteler at osgeo.org> wrote:

> Hi
>
> On Aug 10, 2016 4:34 PM, "Margherita Di Leo" <diregola at gmail.com> wrote:
> >
> > ERROR: Illegal latitude for North
>
> Which projection are you using?
> I suspect latlong?
>
yes, correct.

> > How should I calculate my overlap and my tiles in order to avoid this
> problem?
>
> Please let us know more.. Perhaps also the command with parameters.
>
The data source I'm using is GMTED 2010 at 30 arc second spatial
resolution, global coverage --> created VRT --> linked via r.external -->
set g.region rast=<the whole mosaic>

g.region rast=bln300_mosaic -ap
projection: 3 (Latitude-Longitude)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      89:59:59.5N
south:      70:00:00.5S
west:       179:59:59.5E
east:       179:59:59.5E
nsres:      0:00:30
ewres:      0:00:30
rows:       19200
cols:       43200
cells:      829440000

I wanted to try create 6 large tiles as a first tentative, to see how
r.tile works.
Then i calculate them like this:
rows: 19200 / 2 = 9600
cols: 43200 / 3 = 14400

r.tile input=bln300_mosaic output=tiled width=14400 height=9600 overlap=4
--v
ERROR: Illegal latitude for North

hm..



-- 
Margherita Di Leo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160811/fcf55108/attachment.html>


More information about the grass-user mailing list