[GRASS-user] v.rast.stats error: "Unable to seek"

[Luí­s Moreira de Sousa]

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 <neteler at 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