通过Python脚本调用OSGeo4W Shell执行多轮ogr2ogr请求的问题

没问题，你的思路完全可行——通过subprocess调用OSGeo4W的bat文件，在交互式Shell里批量执行ogr2ogr命令，完美绕开无法安装osgeo库的权限限制。下面给你细化实现方案和代码示例：

核心思路解析 #

• 调用OSGeo4W.bat会启动一个配置好GDAL/OGR环境的Shell，我们可以通过stdin给这个Shell传入命令，stdout捕获输出
• 循环生成WFS请求的ogr2ogr命令，逐个传入Shell执行
• 注意处理Shell的交互模式，避免命令执行阻塞

完整代码示例 #

import subprocess
import time

def run_ogr2ogr_wfs_commands(wfs_urls, output_dir):
    # 启动OSGeo4W Shell进程，替换成你实际的bat文件路径
    osgeo_bat_path = r'C:\Program Files\QGIS 2.18\OSGeo4W.bat'
    proc = subprocess.Popen(
        osgeo_bat_path,
        stdin=subprocess.PIPE,
        stdout=subprocess.PIPE,
        stderr=subprocess.STDOUT,
        text=True,
        bufsize=1,
        universal_newlines=True
    )

    # 先切换到输出目录（可选，方便管理输出文件）
    proc.stdin.write(f'cd /d "{output_dir}"\n')
    proc.stdin.flush()

    # 循环执行每个WFS请求的ogr2ogr命令
    for idx, wfs_url in enumerate(wfs_urls):
        # 生成ogr2ogr命令，这里以导出GeoJSON为例，可按需调整格式和参数
        output_file = f'wfs_layer_{idx}.geojson'
        ogr_cmd = f'ogr2ogr -f GeoJSON "{output_file}" "{wfs_url}"\n'
        
        print(f'正在执行: {ogr_cmd.strip()}')
        proc.stdin.write(ogr_cmd)
        proc.stdin.flush()

        # 根据数据大小调整等待时间，也可通过读取stdout判断执行完成
        time.sleep(5)

    # 执行完所有命令后关闭Shell
    proc.stdin.write('exit\n')
    proc.stdin.flush()

    # 打印所有执行输出
    for line in proc.stdout:
        print(line.strip())

    proc.wait()

# 示例用法
if __name__ == '__main__':
    # 替换成你的WFS服务URL列表
    wfs_urls = [
        'http://your-wfs-server.com/wfs?service=WFS&version=1.0.0&request=GetFeature&typename=your_layer_1',
        'http://your-wfs-server.com/wfs?service=WFS&version=1.0.0&request=GetFeature&typename=your_layer_2'
    ]
    # 替换成你的输出目录
    output_directory = r'C:\temp\wfs_export'
    run_ogr2ogr_wfs_commands(wfs_urls, output_directory)

关键注意事项 #

• 路径校验：确保OSGeo4W.bat的路径正确，QGIS 3.x版本的路径通常是C:\Program Files\QGIS 3.x\OSGeo4W.bat
• 命令合法性：WFS URL要包含service=WFS、request=GetFeature、typename等必要参数，可先在OSGeo4W Shell里单独测试命令是否能运行
• 执行等待：time.sleep()是简单的等待方式，如果需要更精准的判断，可以循环读取proc.stdout直到出现命令执行完成的标识
• 权限与编码：确保输出目录有写入权限，使用text=True避免文本编码乱码问题

这样调整后的代码可以稳定批量执行WFS请求，完全不需要依赖osgeo库，完美适配你的权限受限场景。

内容的提问来源于stack exchange，提问作者user1851594