llamando a gdal_merge.py a otro script de Python que ejecuta procesos GDAL

Quiero usar gdal_merge.pypara fusionar una serie de archivos .tif antes de recortar el archivo fusionado en un límite de shapefile, pero quiero hacerlo como parte de otro script más grande que ejecutará una serie de otros procesos. Soy un novato de Python / GDAL y no estoy seguro de cómo llamar gdal_merge.pya otro script.

No puedo simplemente ejecutarlo gdal_merge.pysolo como es, sino un paso en un script que, con suerte, ejecutará varios procesos. ¿Alguna idea sobre la mejor manera de hacer esto?

La forma más fácil de hacer esto es importando la ruta donde gdal_merge.pyse encuentra, en mi caso, /usr/bin/sustituir la ruta a gdal_mergeen su sistema, que, obviamente, también podría ser una ruta de Windows.

import sys
sys.path.append('/usr/bin/')
import gdal_merge as gm

Ahora tendrá que construir una matriz para sys.argv, como si estuviera llamando gdal_mergedirectamente, por ejemplo,

sys.argv = ['-o','outputfile.tiff','inputfile1.png', 'inputfile2.png', ....'inputfile10.png']
gm.main()

Hay más información sobre esta publicación de desbordamiento de pila

También existe el __init__.pymecanismo , pero esto requiere que el archivo que está importando esté en un subdirectorio desde donde esté ejecutando su archivo de Python.

Editar:

Acabo de leer que esta opción ya se sugirió en los comentarios. De todos modos, por completo. ¿Quizás alguien puede fusionar las respuestas?

Por supuesto, también es posible usar Python subprocess, por ejemplo

import subprocess
def merge(first, second, out_file):
    """
    This utility will automatically mosaic a set of images.
    All the images must be in the same coordinate system and
    have a matching number of bands, but they may be overlapping,
    and at different resolutions. In areas of overlap,
    the last image will be copied over earlier ones.

    :param first:
    :param second:
    :param out_file:
    :return:
    """
    ps = subprocess.Popen(
        ['gdal_merge.py', '-o', out_file,
         '-of', unipath.Path(out_file).ext,
         '-n', NO_DATA_VALUE,
         first, second],
        stdout=subprocess.PIPE
    )
    output = ps.communicate()[0]
    for line in output.splitlines():
        logger.debug("[*] {0}".format(line))