[GRASS-dev] r.tileset in GRASS 7 defect

Hi

Lately I tried calculating tiles for my region using r.tileset in GRASS 7.

On Ubuntu I get no tiles listed:

On my Windows box I get this error-message:

GRASS 7.0.svn> r.tileset sourceproj=“+proj=utm +no_defs +zone=33 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1” maxcols=25000 maxrows=25000 overlap=10

sh: r.tileset: command not found

On my Ubunut box no tiles are printed at all:

GRASS 7.0.svn (Norge_33N_WGS84):~ > r.tileset sourceproj=“+proj=utm +no_defs +zone=33 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1” maxcols=25000 maxrows=25000 overlap=10

Computing tiling…

GRASS 7.0.svn (Norge_33N_WGS84):~ >

Maybe related to:

http://trac.osgeo.org/grass/ticket/1825

Cheers

Stefan

Hi again

Some update on that issue:

Looks like a declaration of floating type variables I necessary, when the tile_bbox w,s,e,n coordinates are calculated.

Furthermore, the bboxesIntersect() method does not work properly.

I could sent a patch for the named issues (based on my modifications, see below).

But there will be still a rounding oddity, because my regions n,s,e,w have no digits, while the tiles get some:

GRASS 7.0.svn (Norge_33N_WGS84):~ > g.region -up

projection: 1 (UTM)

zone: 33

datum: wgs84

ellipsoid: wgs84

north: 7972675

south: 6415595

west: -97865

east: 1165535

nsres: 10

ewres: 10

rows: 155708

cols: 126340

cells: 19672148720

GRASS 7.0.svn (Norge_33N_WGS84):~ > r.tileset1.py sourceproj=“+proj=utm +no_defs +zone=33 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1” maxcols=25000 maxrows=25000 overlap=10 --v

Projecting destination region into source…

Projecting source bounding box into destination…

Computing length of sides of source bounding box…

Computing tiling…

There will be 6 by 7 tiles each 21057 by 22244 cells

-97864.9638333|6415594.9703|112801.679274|6638134.85547|21067|22254

-97864.9638333|6638034.85552|112801.679274|6860674.74064|21067|22254

-97864.9638333|6860474.74074|112801.679274|7083214.62581|21067|22254

-97864.9638333|7082914.62596|112801.679274|7305754.51097|21067|22254

Cheers

Stefan

···

From: grass-dev-bounces@lists.osgeo.org [mailto:grass-dev-bounces@lists.osgeo.org] On Behalf Of Blumentrath, Stefan
Sent: 2. juni 2013 13:31
To: grass-dev@lists.osgeo.org
Subject: [GRASS-dev] r.tileset in GRASS 7 defect

Hi

Lately I tried calculating tiles for my region using r.tileset in GRASS 7.

On Ubuntu I get no tiles listed:

On my Windows box I get this error-message:

GRASS 7.0.svn> r.tileset sourceproj=“+proj=utm +no_defs +zone=33 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1” maxcols=25000 maxrows=25000 overlap=10

sh: r.tileset: command not found

On my Ubunut box no tiles are printed at all:

GRASS 7.0.svn (Norge_33N_WGS84):~ > r.tileset sourceproj=“+proj=utm +no_defs +zone=33 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1” maxcols=25000 maxrows=25000 overlap=10

Computing tiling…

GRASS 7.0.svn (Norge_33N_WGS84):~ >

Maybe related to:

http://trac.osgeo.org/grass/ticket/1825

Cheers

Stefan