Estoy usando GRASS 6.4.3 en Linuxmint 16. Tengo un bash script para GRASS que uso para procesar casi 300 archivos raster. El script primero importa algunos archivos raster ( r.in.gdal ), luego hace algunos cálculos y clasificación (usando r.mapcalc y r.reclass ). Y finalmente transforma el raster resultante en vector ( r.to.vect ) y lo guarda ( v.out.ogr , utilizando type=area ).
Digamos que cada archivo vectorial tiene un campo denominado categoría y que sólo toma dos valores 1 y 2. Como sólo me interesa el área total de los polígonos con categoría = 2 y con un área superior a 10.000 m2, abro cada archivo vectorial en QGIS, edito su tabla de atributos, calculo el área de cada polígono y los elimino todos con categoría=1 y $superficie<10.000 .
Pero tengo que hacerlo manualmente,
¿Hay alguna manera de hacerlo en un proceso por lotes dentro de QGIS?
O mejor aún, ¿es posible hacerlo en GRASS?
¿Es posible editar la tabla de atributos dentro de la consola de comandos en GRASS?
¿Podría hacerlo en el mismo bash script que uso para hacer todos los demás procesos en GRASS?
He estado leyendo este pero no estoy seguro de cómo aplicarlo a mi problema. No he podido encontrar cómo calcular el área de los polígonos y añadirla a una nueva columna usando GRASS. Creo que podría utilizar db.select para buscar los polígonos que quiero eliminar, pero no sé cómo hacerlo.
Gracias por su ayuda.
