mass download vector data from WFS service. Generates a sh script with ogr2ogr calls.

wfs_download.py

from owslib.fes import *
from owslib.etree import etree
from owslib.wfs import WebFeatureService
wfs11 = WebFeatureService(url='http://example.com:8080/geoserver/wfs/', version='1.1.0')

f = open('download.sh', 'wb')
f.write('#!/bin/bash'+'\n')
f.write(''+'\n')


for layername in wfs11.contents:
    if 'namepart'.lower() in layername.lower():     
        cmd='ogr2ogr -progress -f "GeoJSON" {}.geojson WFS:"http://example.com:8080/geoserver/wfs" {}'.format(layername, layername)
        print cmd
        f.write(cmd+'\n')


f.close