Si he entendido bien tu pregunta (y lo que se ha publicado en la lista de correo de GRASS), supongo que sí:
- Un vector puntual de asentamientos (asentamientos múltiples)
- Un segundo vector puntual del yacimiento arqueológico (punto único)
- Un mapa rasterizado en el que cada píxel representa el tiempo de paso por esa celda
En este caso creo que r.cost ( manual ) y una simple expresión r.mapcalc hará lo que quieres.
En primer lugar, ejecute r.cost con el ráster de tiempo de recorrido como entrada, utilice el vector de asentamientos como punto_de_inicio y el yacimiento arqueológico como punto_de_fin, y llame a la salida "accum_cost". A continuación, suponiendo que la trama de tiempo de caminata estaba en segundos, ejecute:
r.mapcalc "within_1hr = if(accum_cost<=3600, 1, null())"
Esto debería crear un nuevo raster tal que todas aquellas áreas dentro de 1 h a pie tendrán el valor 1, y todas las demás áreas son nulas.
El valor añadido de r.walk, frente al más sencillo r.cost, es la opción de combinar tanto una capa de elevación y una capa de fricción. Así que si en tu modelo, caminar desde los asentamientos hasta el lugar implica tanto subir colinas como atravesar una jungla, entonces podrías combinar los costes de ambos con r.caminar. Incluso puede establecer los parámetros de velocidad (probablemente ya haya visto la función r.walk manual para esto).
