convert_tif_to_tiles.py

import os
import subprocess
import tempfile

def convert_tif_to_tiles(input_tif, output_directory):

    # Ensure the output directory exists
    os.makedirs(output_directory, exist_ok=True)

    # Create a temporary file for the `gdal_translate` output
    with tempfile.NamedTemporaryFile(suffix=".tif", delete=False) as temp_tif:
        temp_tif_path = temp_tif.name

    try:
        # Step 1: Run gdal_translate to set the 4th band as the transparency mask
        translate_command = ['gdal_translate', '-mask', '4', input_tif, temp_tif_path]
        subprocess.run(translate_command, check=True)

        # Step 2: Run gdal2tiles on the result of gdal_translate
        #todo: -x necessary to avoid producing empty tiles?
        #todo: best zoom level ?
        #todo: avoid creating html with ---webviewer none

        tiles_command = ['gdal2tiles.py', '-z', '5-10', '-v', '-r', 'bilinear', temp_tif_path, output_directory]
        subprocess.run(tiles_command, check=True)

    except subprocess.CalledProcessError as e:
        print(f"An error occurred: {e}")

    finally:
        # Clean up the temporary file
        os.remove(temp_tif_path)

    # todo: Generate a Folder List File: Create a JSON file (or another format)
    #  containing all folder names, such as tile_folders.json, to be put in the amazon S3 bucket.

if __name__ == "__main__":
    input_tif = "/media/maffe/nvme/iceboost_global_deploy/RGI70G/rgi11/mapbox/RGI2000-v7.0-G-11-02596_Mapbox.tif"
    glacier_name = input_tif.split('/')[-1].replace('.tif', '').replace('_Mapbox', '')
    output_dir = f"/home/maffe/PycharmProjects/iceboost_webapp/tiles/RGI70G/{glacier_name}/"
    convert_tif_to_tiles(input_tif, output_dir)