循环遍历GDAL OSGeo4W中的shapefile

Question

我有一个要转换为geoJSON的形状文件列表。我正在使用GDAL OSGeo4W外壳进行转换。

首先，我定义了坐标系(因为它们不是为shapefile定义的)，然后我转换它们。

我可以很容易地为单个文件做这件事，但是因为我有超过200个shapefile，所以我想遍历它们。

我的个人命令是：

ogr2ogr -s_srs http://spatialreference.org/ref/epsg/27700/ srs_def x.shp

ogr2ogr -f GeoJSON -t_srs EPSG:27700 x.geojson x.shp

Answer 1

创建一个批处理脚本以在主OSGeo4w shell中运行，或创建一个shell脚本以在MSYS shell中运行。要了解它们的外观，请参阅https://gis.stackexchange.com/a/7806/1872

Answer 2

批量处理gdal或ogr命令的一个简单方法是使用python (您可以通过启动python直接在OSGeo4w中使用它)：

import os
os.system('ogr2ogr -s_srs http://spatialreference.org/ref/epsg/27700/ srs_def x.shp')
os.system('ogr2ogr -f GeoJSON -t_srs EPSG:27700 x.geojson x.shp')

然后，你可以只做一个简单的循环。