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.
Soluções usando rasterio funcionam, todos nós sabemos disso. No entanto, instalar o rasterio e o GDAL no mesmo ambiente do conda pode ser um desafio… Então, geralmente, eu tento não instalar os dois juntos. Vou mostrar para vocês a maneira que eu acho que é a mais fácil para salvar dados raster usando o GDAL no Python.
Hoje, estou apresentando a criação de uma camada virtual como uma alternativa ao processo de mosaicar rasters. Bem, a principal razão é que um arquivo de layer virtual ocupa muito menos espaço no disco rígido do que um mosaico raster. Se os meus rasters tiverem cada um 80 KB, isso não é um problema. Mas e se eles tiverem 1 GB cada?
Eu estou usando o pacote pygrib para abrir os arquivos. Em sua documentação, eles descrevem para duas maneiras de instalar o pygrib, usando pip ou conda. No entanto, só consegui fazer funcionar usando conda. Execute isso em seu prompt do Anaconda para instalar o pygrib
Esta semana, estou discutindo um erro que costumava ser um simples aviso (warning) no antigo rasterio - e eu ignorei. Aprendi a lição! Eu estava salvando rasters GeoTIFF, como eu sempre faço, mas em um ambiente Python atualizado (com o rasterio 1.2.6) quando obtive este erro.
A Calculadora Raster do SAGA é útil porque ela permite fazer alguns cálculos que a Calculadora Raster nativa do QGIS não faz. Especialmente, ela permite usar as coordenadas X e Y de um raster nos cálculos. Vamos lá?
[RESOLVIDO] ERROR 3: Free disk space available is 85802630536 bytes, whereas 730874237760027488 are at least necessary. You can disable this check by defining the CHECK_DISK_FREE_SPACE configuration option to FALSE.
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.
Hoje, vou mostrar aqui no blog como georreferenciar uma imagem usando o QGIS 3.18. Talvez você queira georreferenciar uma imagem de um artigo de revista, ou talvez seja uma fotografia aérea. Não importa! Se você conhece as coordenadas de pelo menos dois pontos não colineares, é possível adicionar geolocalização à sua imagem!
Hoje, vamos discutir diferentes maneiras de cortar (recortar) um raster no QGIS 3.X e o que elas fazem. Em especial, estou comparando as ferramentas do QGIS provenientes do GDAL “Clip raster by mask layer” (Recortar raster por camada máscara), “Clip raster by extent” (Recortar raster pela extensão) e “Warp (reproject)”, para que servem, e o que exatamente fazem com os dados raster originais.