Hi again Markus, below is a small example with an extent of 10^12 cells. It does not get to v.rast.stats, as the r.random.surface command fails before that (without error message). I will do more experiments in the following days.
Regards. # Set region for zones g.region n=10000000 s=0 w=0 e=10000000 nsres=2000 ewres=2000 # Create zones map r.mapcalc expr="zones=row()*10000+col()" r.to.vect -b input=zones output=zones column=zone type=area # Set computation region g.region n=10000000 s=0 w=0 e=10000000 nsres=100 ewres=100 # Generate random raster : this command fails without error message r.random.surface output=random seed=420 distance=1000 --overwrite # Stats v.rast.stats map=zones raster=random method=number column_prefix=stats -- Luís ‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐ On Wednesday, May 19, 2021 5:09 PM, Markus Neteler <nete...@osgeo.org> wrote: > > GRASS GIS 7 supports the off_t type, hence it can address an enormous > amount of raster data (like 1e+18 pixels or more), see > https://grasswiki.osgeo.org/wiki/GRASS_GIS_Performance#Large_raster_data_processing > > Bugs and limitations are always possible but there is interest in > fixing them wherever possible. > > > I will try to have a look again at an agnostic replication procedure. > > Alright. > > Best > Markus _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user