Quelqu'un a-t-il un moyen éloquent d'empiler plusieurs fichiers .tif dans une pile à plusieurs bandes en utilisant Rasterio et / ou GDAL?
Je cherche un moyen d'éviter d'utiliser une commande de sous-processus comme gdal_merge.py et plutôt de l'avoir dans le cadre de mon script python.
Je sais que Rasterio et GDAL lisent les fichiers .tif sous forme de tableaux, mais comment puis-je empiler ces tableaux et écrire le résultat sous forme de bandes distinctes?