Cómo cortar una imagen ráster en QGIS-GRASS usando g.region sin preservar la resolución original

Si se trabaja en el entorno de QGIS existen varias maneras de cortar una imagen ráster. Sin embargo, si se importan los ráster al entorno de GRASS no existe ningún módulo que incluya la opción “clip”, “cut” o “crop” que nmemotécnicamente nos señale una acción de corte. Sin embargo, esa tarea no es difícil de ejecutar en GRASS si se recuerda que el alcance de ejecución de los módulos viene circunscrita a la llamada región de trabajo de GRASS. Esta puede ser modificada manualmente por arrastre del ratón haciendo click en el icono de configuración de la región de GRASS (último en la barra de herramientas de GRASS) o mediante cónsola de GRASS usando el comando de propósitos generales g.región. Esto último es lo que se va a utilizar porque permite tener un control preciso sobre la extensión de la región de corte.

Calcular media aritmética por columnas en un ráster usando QGIS y lenguaje R

Para calcular medias aritméticas, por columnas, en un ráster es necesario pasar primero el ráster a formato ascii. Para ello, una vez cargado el ráster (yo empleo generalmente el formato *.tif), se apela a Ráster -> Conversión -> Traducir (convertir formato) y en la ventana de diálogo que aparece escoger el nombre adecuado del ráster (con la terminación *.asc), seleccionando en “Ficheros de tipo” la opción [GDAL] Arc/Info ASCII Grid (*.asc, *.ASC) y luego Guardar -> Aceptar (se asume que la proyección es la misma que la de la imagen fuente).

Script de bash para copiar columnas de tablas externas con db.execute (GRASS-QGIS)

En el artículo pasado se consideró la elaboración de un script de bash para determinar raíces cuadradas (extensible a otras funciones especiales como logaritmos o trigonométricas) en registros de tablas atributivas con db.execute (GRASS-QGIS). La razón de esta búsqueda fue que no se pueden realizar este tipo de operaciones con el previsible v.db.update desde el entorno de GRASS-QGISS (el driver dbf no lo soporta; ni siquiera para los comandos UPDATE de SQL). Las raíces cuadradas aparecen en la fórmula de la distancia y los logaritmos en las funciones “p” para manejar concentraciones muy pequeñas; entre otras. Aunque la calculadora de campos de QGIS funciona bien en las determinaciones con funciones especiales lo ideal sería hacer todas las operaciones dentro de GRASS y no externamente para no tener que, posteriormente, importar los vectoriales al formato de GRASS con la eventual pérdida de tiempo y evitar la posibilidad de automatizar procesos repetitivos mediante scripts en ambiente de GRASS.

Script de bash para determinar raíz cuadrada en registros con db.execute

Hasta ahora, no he encontrado la manera de calcular raíces cuadradas (y otras funciones especiales como logaritmos, funciones trigonométricas, etc) usando la opción v.db.update de GRASS-QGIS para actualizar simultáneamente todos los registros de un campo. Sin embargo, mediante un script de bash, los cálculos con raíces cuadradas (y otras funciones especiales) pueden realizarse con la librería de funciones (-l) de la calculadora de precisión arbitraria bc y las actualizaciones simultáneas, en cada uno de los registros, mediante la sentencia UPDATE de SQL empleando el módulo db.execute.

Script de bash para determinar promedio y desviación estándar en columnas con v.db.select (GRASS-QGIS)

Como una extensión y generalización del script para la suma de atributos en columnas se propone éste para determinar promedio y desviación estándar con la posibilidad de escoger la database y la columna sobre la que se desean tales operaciones. Para el despliegue de las databases (incluye la de los vectoriales) se usa db.tables y para las columnas db.columns. El script se encuentra a continuación:

Script de bash para determinar valor máximo en columnas con v.db.select (GRASS-QGIS)

En artículos pasados se propuso la determinación del índice de alargamiento de una cuenca hidrográfica cuyo valor se estimaba como la relación del largo del río entre la anchura máxima de la cuenca perpendicular al río. Para determinar de manera expedita el largo del río se creo un script de bash. Para la anchura máxima se va a hacer lo mismo pero antes se va exponer un recuento del procedimiento para justificar la estructura del script.

Script de bash para sumar atributos en columnas con v.db.select (GRASS-QGIS)

En un artículo reciente tuve la necesidad de obtener la suma de todos los valores de un campo en una tabla atributiva para estimar un índice. Por premura, usé la salida del comando: