Mosaic_Landcover.py

# mosaic

import rasterio
from rasterio.merge import merge
from pathlib import Path
import glob

# Caminho para a pasta onde estão os GeoTIFFs baixados
input_dir = Path("/home/adriano/Doutorado/landcover/EarthEngine_Exports/")  # ajuste o nome da pasta se necessário

# Lista todos os arquivos .tif
tif_files = sorted(glob.glob(str(input_dir / "*.tif")))

# Abre todos os arquivos como objetos rasterio
src_files_to_mosaic = [rasterio.open(fp) for fp in tif_files]

# Realiza a fusão
mosaic, out_transform = merge(src_files_to_mosaic)

# Usa metadados do primeiro arquivo como base
out_meta = src_files_to_mosaic[0].meta.copy()
out_meta.update({
    "driver": "GTiff",
    "height": mosaic.shape[1],
    "width": mosaic.shape[2],
    "transform": out_transform,
    "compress": "lzw"
})

# Caminho para salvar o mosaico final
output_file = "landcover_peru_merged.tif"

# Salva o arquivo final
with rasterio.open(output_file, "w", **out_meta) as dest:
    dest.write(mosaic)

print(f"Mosaico salvo como {output_file}")