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.
ERROR 1: Attempt to write non-multipoint (POINT) geometry to multipoint shapefile. ERROR 1: Unable to write feature 0 from layer points_layer. ERROR 1: Terminating translation prematurely after failed translation of layer points_layer (use -skipfailures to skip errors)
Neste post, vou mostrar como realizar este processo usando duas ferramentas diferentes, uma delas é nativa do QGIS e a outra é um plugin disponível no repositório de plugins. Em seguida, vou comentar sobre como amostrar um raster usando interpolação do tipo cúbica, spline cúbica ou bilinear.
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.
[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.
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.
Hoje, vou escrever sobre como abrir, editar e salvar arquivos raster usando Python, se você já estiver familiarizado com Python (principalmente Anaconda). Talvez você precise executar uma função Python em cada pixel de um raster, ou talvez você deseje fazer operações personalizadas que não seriam possíveis em softwares com interfaces gráficas como o QGIS. O Python também pode ser usado para fazer processamento em lotes. Para o processamento de hoje, usaremos o GDAL no Python, instalado através do Conda.
Uma das maneiras mais diretas, e por outro lado, uma das mais complicadas, de fazer processamento em lotes é usando scripts Bash. Bash é um tipo de Shell, o que significa que é executado diretamente em uma interface de usuário que acessa os serviços do sistema operacional.
Na semana passada, eu falei um pouco sobre como usar o SAGA na linha de comando do Prompt. Hoje, vou discutir soluções para o problema “Error: Catchment Area” no SAGA GIS. Ou “Error: seu raster de entrada”, ou “Input layers do not have the same grid extent”.