GEOTIFF

Abrindo e visualizando uma Fmask no QGIS usando um código simples em Python no PyQGIS

Hoje, eu vou apresentar para vocês um gist (fragmento de código) que desenvolvi, para separar as informações de máscaras contidas dentro de uma Fmask em diferentes bandas, com o objetivo de visualizá-las separadamente. Nós vamos visualizá-las no QGIS, usando o Python Console do QGIS, também conhecido como PyQGIS.

Quando se faz um mosaico de dois rasters, o tamanho do raster resultante é quase sempre maior do que a soma dos tamanhos dos rasters originais. Por quê?

Resposta curta: seu software de SIG precisa preencher os espaços em branco no espaço em branco com pixels de valor NoData. Resposta longa: vou demonstrar o que acontece, por meio de um exemplo no QGIS 3.18, e mostrar porque o raster final tende a ser maior do que os outros dois somados.