gọi gdal_merge.py vào một tập lệnh python khác đang chạy các quy trình GDAL

Tôi muốn sử dụng gdal_merge.pyđể hợp nhất một loạt các tệp .tif trước khi cắt tệp đã hợp nhất thành một ranh giới shapefile, nhưng muốn thực hiện như một phần của một tập lệnh lớn hơn sẽ thực thi một số quy trình khác. Tôi là một người mới chơi python / GDAL và không biết làm thế nào để gọi về gdal_merge.pymột kịch bản khác.

Tôi không thể gdal_merge.pytự mình chạy một mình mà chỉ một bước trong một kịch bản hy vọng sẽ thực thi một số quy trình. Bất kỳ suy nghĩ về cách tốt nhất để làm điều này?

Cách dễ nhất để làm điều này là bằng cách nhập đường dẫn gdal_merge.pyđược đặt, trong trường hợp của tôi, /usr/bin/- thay thế bằng đường dẫn đến gdal_mergetrên hệ thống của bạn, rõ ràng cũng có thể là đường dẫn Windows.

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

Bây giờ bạn sẽ phải xây dựng một mảng cho sys.argv, như thể bạn đang gọi gdal_mergetrực tiếp, ví dụ:

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

Có nhiều thông tin hơn về bài viết Stack Overflow này

Ngoài ra còn có __init__.pycơ chế , nhưng điều này đòi hỏi tệp bạn đang nhập phải nằm trong thư mục con của bất cứ nơi nào bạn đang chạy tệp python của mình.

Biên tập:

Chỉ cần đọc rằng tùy chọn này đã được đề xuất trong các ý kiến ​​đã. Dù sao, vì lợi ích hoàn toàn. Có lẽ ai đó có thể hợp nhất các câu trả lời?

Tất nhiên cũng có thể sử dụng python's subprocess, vd

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))