Как объединить множество растров со скриптом Python в QGIS, используя GDAL? - PullRequest
0 голосов
/ 01 июля 2019

У меня есть много растров, которые я пытаюсь объединить, чтобы получить контуры.Они в формате .bil, так что я понимаю, что мне нужно сначала перевести их, используя gdal, но ни одна из моих попыток не сработала.

Я пытался использовать gdal_translate и gdal_merge безуспешно.Ниже приведен пример двух .bils, которые я попытался объединить.

import subprocess
import sys
import gdal

#imagea ='C:\\Users\\jeswa\\Desktop\\New folder\\N34_E\\n34_e024_1arc_v3.bil'
#imageb ='C:\\Users\\jeswa\\Desktop\\New folder\\N34_E\\n34_e025_1arc_v3.bil'
#output = 'C:\\Users\\jeswa\\Desktop\\test'
#subprocess.call(['gdal_merge', '-o', output, imagea, imageb])

Traceback (most recent call last):
  File "C:\OSGEO4~1\apps\Python37\lib\code.py", line 90, in runcode
    exec(code, self.locals)
  File "<input>", line 1, in <module>
  File "C:\OSGEO4~1\apps\Python37\lib\subprocess.py", line 304, in call
    with Popen(*popenargs, **kwargs) as p:
  File "C:\OSGEO4~1\apps\Python37\lib\subprocess.py", line 756, in __init__
    restore_signals, start_new_session)
  File "C:\OSGEO4~1\apps\Python37\lib\subprocess.py", line 1155, in _execute_child
    startupinfo)
FileNotFoundError: [WinError 2] The system cannot find the file specified
...