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

[Blumentrath, 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 at lists.osgeo.org [mailto:grass-dev-bounces at lists.osgeo.org] On Behalf Of Blumentrath, Stefan
Sent: 2. juni 2013 13:31
To: grass-dev at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20130603/d1994087/attachment-0001.html>