[GRASS-dev] [GRASS GIS] #1825: r.tileset only returns one tile

GRASS GIS trac at osgeo.org
Tue Jun 4 00:52:29 PDT 2013

#1825: r.tileset only returns one tile
 Reporter:  pvanbosgeo   |       Owner:  grass-dev@…              
     Type:  defect       |      Status:  new                      
 Priority:  normal       |   Milestone:  7.0.0                    
Component:  Raster       |     Version:  svn-trunk                
 Keywords:  r.tileset    |    Platform:  Unspecified              
      Cpu:  Unspecified  |  

Comment(by sbl):

 Some variant of the same problem I guess:

 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):~ >

 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.

 The provided patch should fix these issues.

 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

Ticket URL: <http://trac.osgeo.org/grass/ticket/1825#comment:2>
GRASS GIS <http://grass.osgeo.org>

More information about the grass-dev mailing list