Problema con GRASS i.tasscap

Question

Estoy tratando de llevar a cabo las transformaciones de Tasseled Cap en GRASS GIS. Tengo una imagen Landsat 8 con las bandas 2-7 que son necesarias para el cálculo. He importado las bandas a GRASS y he intentado realizar el cálculo, pero cuando lo hago, sale un tamaño de imagen de 1x1. He comprobado todas las bandas de entrada y todas son de 8202 x 8113 (filas x columnas).

La siguiente captura de pantalla es la salida de r.info en el archivo de salida. Creo que está haciendo el cálculo en todo el mapa en lugar de píxel por píxel, pero no tengo idea de por qué. ¿Alguien sabe cómo arreglar esto?

Estoy usando Fedora 27 y GRASS GIS 7.2.2. He comprobado en Internet pero nadie más parece tener el mismo problema.

r.in.gdal input=LC08_L1TP_208023_20170408_20170414_01_T1_B2_2157.TIF output=B2

r.in.gdal input=LC08_L1TP_208023_20170408_20170414_01_T1_B3_2157.TIF output=B3

r.in.gdal input=LC08_L1TP_208023_20170408_20170414_01_T1_B4_2157.TIF output=B4

r.in.gdal input=LC08_L1TP_208023_20170408_20170414_01_T1_B5_2157.TIF output=B5

r.in.gdal input=LC08_L1TP_208023_20170408_20170414_01_T1_B6_2157.TIF output=B6

r.in.gdal input=LC08_L1TP_208023_20170408_20170414_01_T1_B7_2157.TIF output=B7

i.tasscap sensor=landsat8_oli input=B2,B3,B4,B5,B6,B7 output=tct_test

Answer 1

Parece que no has puesto el GRASS región correctamente. Después de importar sus datos, debe establecer la región a la extensión y resolución de una de sus bandas, por ejemplo:

g.region rast=B2

y luego tratar de ejecutar el proceso de nuevo.

Puede comprobar que la configuración regional es correcta:

g.region -p 

y deberías encontrar que son consistentes con la salida de:

gdalinfo LC08_L1TP_208023_20170408_20170414_01_T1_B2_2157.TIF