Download GPKG from WFS example

dl_from_wfs.py

import sys

try:
    from osgeo import ogr, osr, gdal
except:
    sys.exit('ERROR: cannot find GDAL/OGR modules')

# Set the driver (optional)
wfs_drv = ogr.GetDriverByName('WFS')

# Speeds up querying WFS capabilities for services with alot of layers
gdal.SetConfigOption('OGR_WFS_LOAD_MULTIPLE_LAYER_DEFN', 'NO')

# Set config for paging. Works on WFS 2.0 services and WFS 1.0 and 1.1 with some other services.
gdal.SetConfigOption('OGR_WFS_PAGING_ALLOWED', 'YES')
gdal.SetConfigOption('OGR_WFS_PAGE_SIZE', '1000')

# Open the webservice
url = 'https://wxs.ign.fr/ortho/geoportail/wfs?SERVICE=WFS'
wfs_ds = wfs_drv.Open('WFS:' + url)
if not wfs_ds:
    sys.exit('ERROR: cannot open WFS datasource')
else:
    pass

# Get a specific layer
layer = wfs_ds.GetLayerByName("ORTHOIMAGERY.ORTHOPHOTOS.GRAPHE-MOSAIQUAGE:graphe_bdortho")
dr = ogr.GetDriverByName( 'GPKG' )
ds = dr.CreateDataSource( 'graphe_bdortho.gpkg' )
ds.CopyLayer(layer, 'local_copy')