[GRASSLIST:2734] r.neighbors with size>25

GRASS GRASS user
1

I want to assign to each grid cell the number of grid cells with a certain characteristics in a neighborhood of 30x30 cells. The problem with r.neighbors is that it allows only neighborhood sizes up to 25.
Any ideas how to do it (without typing in 30x30 coordinates of the neighborhood cells in r.mapcalc...)?

Christof

2

Christof Bigler wrote:

I want to assign to each grid cell the number of grid cells with a
certain characteristics in a neighborhood of 30x30 cells. The problem
with r.neighbors is that it allows only neighborhood sizes up to 25.
Any ideas how to do it (without typing in 30x30 coordinates of the
neighborhood cells in r.mapcalc...)?

If you have the source code, remove the line:

  parm.size->options = "1,3,5,7,9,11,13,15,17,19,21,23,25" ;

from src/raster/r.neighbors/cmd/main.c and re-compile.

--
Glynn Clements <glynn.clements@virgin.net>