vascobnunes, вот как я решил эту проблему, используя скрипт Python для последовательного объединения нескольких инструкций ogr2ogr. Вы можете легко преобразовать его в пакетный скрипт, в основном я просто объединяю инструкции ogr2ogr ( cmd
), а затем выполняю их, вызывая os.system(cmd)
, передав команду ogr2ogr, которую я соединял вместе.
Секретным оружием является ( как продемонстрировал capooti ) применение OGR_SQL для наложения имени файла в качестве постоянного значения набора исходных данных, который вы добавляете в свой результат слияния.
В моем примере -sql
флаг обрабатывает это, в коде это выглядит так:
-sql "SELECT \'' + filename + '\' AS filename, * FROM ' + filenameNoExt + '"'
Но это сбивает с толку, потому что мне нужно применять одинарные кавычки и двойные кавычки в результате объединения. Для этого мне нужно избегать одинарных кавычек (то есть \ '), чтобы использовать их «по-настоящему». Так что для удобства чтения, это помогает увидеть его без переменных и escape-последовательностей. Если вы делаете вид, что имя файла было «roads1» для конкретной итерации, результирующая конкатенация будет выглядеть следующим образом в предложении ogr2ogr:
-sql "SELECT 'roads1.shp' AS filename, * FROM roads1"
Этот скрипт .py представляет собой объединение трех трюков, которые я украл из matt wilkie (пустой клон шейп- файла ), j03lar50n (добавление столбца в шейп-файл с использованием ogrinfo и ogr_sql) и capooti (использование ogr_sql для наложения фиксированного значения столбца на всех записях в шейп-файле). Итак, вот полный сценарий:
# merge_shps.py
import os
path = "D:/GIS/01_tutorials/ND_Roads/extracted" # path to your folder of .shp files
merge = "merge_filename" # this will be the name of your merged result
directory = os.listdir(path)
count = 0
for filename in directory:
if ".SHP" in filename.upper() and not ".XML" in filename.upper():
# On the first pass, create a clone and add the filename column.
if count == 0:
# Make a clone (matt wilkie)..
cmd = 'ogr2ogr ' + path + '/' + merge + '.shp ' + path + '/' + filename + ' -where "FID < 0"'
os.system(cmd)
# Add the field (j03lar50n)..
cmd = 'ogrinfo ' + path + '/' + merge + '.shp -sql "ALTER TABLE ' + merge + ' ADD COLUMN filename character(50)"'
os.system(cmd)
# Now populate the data (capooti)..
print "Merging: " + str(filename)
# You'll need the filename without the .shp extension for the OGR_SQL..
filenameNoExt = filename.replace(".shp","")
cmd = 'ogr2ogr -f "esri shapefile" -update -append ' + \
path + '/' + merge + '.shp ' + \
path + '/' + filename + \
' -sql "SELECT \'' + filename + '\' AS filename, * FROM ' + filenameNoExt + '"'
# Uncomment this line to spit the ogr2ogr sentence to the terminal..
#print "\n" + cmd + "\n"
os.system(cmd)
count += 1