我正在合并一些shapefile，但在QGIS中这样做时遇到了一些问题，因此我直接使用了ogr2ogr。我正在这样做（批量）：
ogr2ogr -overwrite %destination% %n1%
ogr2ogr -update -append %destination% %n2% -nln all_new
ogr2ogr -update -append %destination% %n3% -nln all_new
ogr2ogr -update -append %destination% %n4% -nln all_new它工作正常，但现在我需要在生成的shapefile中有一个字段，其中包含我合并的原始shapefile的名称。听起来并不困难，但我没有设法做到这一点。
有人可以帮忙吗？谢谢！
Answers:
使用小的脚本，这将是可行的。使用类似以下的方法，您应该能够将列添加到文件夹中所有shapefile中的shapefile中，并将它们合并到merged.shp文件中
for %f in (*.shp) do (
ogrinfo %f -sql "ALTER TABLE %f ADD COLUMN filename character(15)"
ogrinfo %f -sql "UPDATE TABLE %f filename = '%f'"
ogr2ogr -update -append merged.shp %f -f “esri shapefile” -nln merge
)编辑：与Bash脚本相同，但进行一些更改以使其起作用：
for f in *.shp
do
base=${f%.shp}
ogrinfo $f -sql "ALTER TABLE $base ADD COLUMN filename character(15)"
ogrinfo $f -dialect SQLite -sql "UPDATE $base SET filename = '$base'"
ogr2ogr -update -append merged.shp $f
done我将使用-sql选项，并通过以下方式导入shapefile：
ogr2ogr -update -append %destination% %n2% -sql 'SELECT "%n2%" as SHAPE_ORIG, field1, field2, ... FROM %n2%'有一些合并shapefile的方法。
mkdir merged
for %f in (*.shp) do (
if not exist merged\merged.shp (
ogr2ogr -f “esri shapefile” merged\merged.shp %f) else (
ogr2ogr -f “esri shapefile” -update -append merged\merged.shp %f -nln Merged )
)并将属性添加到shapefile @dango Directon是好的。您可以使用layer.CreateField（field_name）创建一个新列，该列从
import os
shapeFileName = os.path.splitext("your_shape_file_path")[0]我希望它可以帮助您...
vascobnunes，这是我通过使用Python脚本将多个ogr2ogr指令菊花链式链接在一起的方法来解决此问题。您可以轻松地将其转换为批处理脚本，基本上我只是将ogr2ogr指令（cmd）连接在一起，然后通过调用os.system(cmd)，传入我串联在一起的ogr2ogr命令来执行它们。
秘密武器是（如capooti所示）应用OGR_SQL将文件名强加为要追加到合并结果中的源数据集的常量值。
在我的示例中，-sql标志在代码中处理如下：
-sql "SELECT \'' + filename + '\' AS filename, * FROM ' + filenameNoExt + '"'但这令人困惑，因为我需要在结果串联中应用单引号和双引号。为此，我必须转义单引号（即\'）以“实际”使用它们。因此，为了提高可读性，有助于查看没有变量和转义序列的代码。如果您将文件名假装为特定迭代的“ roads1”，则在ogr2ogr语句中产生的连接看起来像这样：
-sql "SELECT 'roads1.shp' AS filename, * FROM roads1"这个.py脚本是我从Matt wilkie（一个空的shapefile的克隆），j03lar50n（使用ogrinfo和ogr_sql将一个列添加到shapefile中）和capooti（使用ogr_sql施加一个固定的列值）中窃取的三个技巧的组合在shapefile中的所有记录上）。所以这是完整的脚本：
# 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从shapefile的文件夹中添加带有源文件名的列。需要GDAL 1.10dev，我删除.shp扩展名的尝试无效，但总的来说可以。-我想可以将其添加到与OGR合并的行中。
for f in *.shp;
do
name=${f%.shp}
/Users/you/gdal_src/bin/ogrinfo $f -sql "ALTER TABLE $name ADD COLUMN filename character(21)"
/Users/you/gdal_src/bin/ogrinfo $f -dialect SQLite -sql "UPDATE $name SET filename = '$f'"
done;在QGIS内部，您可以添加Merge Shapefile插件。有一个选项“添加带有文件名的列”
JaaKL答案的略微修改版本。请注意，-append foo.shp和-nln foo需要匹配。另外，请注意使用SQLite方言（GDAL显然不接受关键字“ Update”，因此必须立即使用SQLite方言），并在单词“ UPDATE”后缺少关键字“ TABLE”（不需要）或被SQLite接受）。
for %%f in (*.shp) do (
if not "%%f" == "merge.shp" (
ogrinfo %%f -sql "ALTER TABLE %%~nf ADD COLUMN fname character(15)"
ogrinfo %%f -dialect SQLite -sql "UPDATE %%~nf SET fname = '%%~nf'"
ogr2ogr -update -append merge.shp %%f -f "ESRI SHAPEFILE" -nln merge
)
)