通过GDAL OSGeo4W中的shapefile循环

Question

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

首先，我定义坐标系（因为它们没有为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中运行的批处理脚本，或者在MSYS shell中运行的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')

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