[GRASS-dev] Retiling functionality for 2D raster?

GRASS GIS GRASS GIS developer
1

Hi,

for voxels, we have the r3.retile module to subset a volume:
http://grass.osgeo.org/grass70/manuals/r3.retile.html

I am lacking a "r.retile" module to easily split a big 2D raster
map into tiles, e.g. for parallelized analysis on a cluster.
To go the v.mkgrid way might be a bit tedious.

Any ideas? (Ab)use the segment lib to achieve this?

thanks,
Markus

2

Markus Neteler wrote:

I am lacking a "r.retile" module to easily split a big 2D raster
map into tiles, e.g. for parallelized analysis on a cluster.
To go the v.mkgrid way might be a bit tedious.

r.tile added in r55511.

--
Glynn Clements <glynn@gclements.plus.com>

3

On Mon, Mar 25, 2013 at 3:55 AM, Glynn Clements
<glynn@gclements.plus.com> wrote:

Markus Neteler wrote:

I am lacking a "r.retile" module to easily split a big 2D raster
map into tiles, e.g. for parallelized analysis on a cluster.
To go the v.mkgrid way might be a bit tedious.

r.tile added in r55511.

Wonderful, thanks.

How about renaming r3.retile to r3.tile (to gain name consistency)?

Markus

4

Markus Neteler wrote:

>> I am lacking a "r.retile" module to easily split a big 2D raster
>> map into tiles, e.g. for parallelized analysis on a cluster.
>> To go the v.mkgrid way might be a bit tedious.
>
> r.tile added in r55511.

Wonderful, thanks.

How about renaming r3.retile to r3.tile (to gain name consistency)?

They do different things. r3.retile generates a single output map, but
with the specified tiling (3d rasters are tile-based). r.tile
generates a separate raster for each tile (equivalent to running
g.region+r.resample in a double loop).

--
Glynn Clements <glynn@gclements.plus.com>

5

On Mon, Mar 25, 2013 at 3:55 AM, Glynn Clements
<glynn@gclements.plus.com> wrote:

Markus Neteler wrote:

I am lacking a "r.retile" module to easily split a big 2D raster
map into tiles, e.g. for parallelized analysis on a cluster.
To go the v.mkgrid way might be a bit tedious.

r.tile added in r55511.

I found a small issue:

g.region rast=elevation -p
r.tile input=elevation output=elev_tile width=135 height=150
WARNING: category information for [elev_tile] in [user1] missing or invalid
WARNING: category information for [elev_tile] in [user1] missing or invalid
...

Could that be fixed/suppressed?

thanks
Markus