[gdal-dev] Converting WFS to shapefile

Paul Meems bontepaarden at gmail.com
Fri Jun 18 02:07:03 PDT 2021


Hi List,

I'm having trouble converting the result of a WFS request to a shapefile.
I'm using two servers, one is working the other not.

My command for the working server:
ogr2ogr -f "ESRI Shapefile" -s_srs "EPSG:4326" -t_srs "EPSG:28992"
"Perceel-bbox.shp" "WFS:
https://geodata.nationaalgeoregister.nl/kadastralekaart/wfs/v4_0?service=WFS&request=GetFeature&version=2.0.0&typename=kadastralekaartv4:perceel&srsname=EPSG:4326&outputFormat=application/json&bbox=163049,467594,163256,468020,EPSG:28992
"

And the not working server:
ogr2ogr -f "ESRI Shapefile" --debug ON  -s_srs "EPSG:4326" -t_srs
"EPSG:28992" "Perceel-aan.shp" "WFS:
https://geodata.nationaalgeoregister.nl/aan/wfs?service=WFS&request=GetFeature&version=2.0.0&typename=aan&srsname=EPSG:4326&outputFormat=application/json&bbox=163049,467594,163256,468020,EPSG:28992
"

I don't get an error, the file is just not created.
The debug messages:
WFS:
https://geodata.nationaalgeoregister.nl/aan/wfs?SERVICE=WFS&REQUEST=GetCapabilities&version=2.0.0&srsname=EPSG:4326&bbox=163049,467594,163256,468020,EPSG:28992
HTTP: Fetch(
https://geodata.nationaalgeoregister.nl/aan/wfs?SERVICE=WFS&REQUEST=GetCapabilities&version=2.0.0&srsname=EPSG:4326&bbox=163049,467594,163256,468020,EPSG:28992
)
HTTP: libcurl/7.74.0-DEV OpenSSL/1.1.1i zlib/1.2.11
WFS: Paging support with page size 15000
WFS: No transaction support
GDAL: GDALOpen(WFS:
https://geodata.nationaalgeoregister.nl/aan/wfs?service=WFS&request=GetFeature&version=2.0.0&typename=aan&srsname=EPSG:4326&outputFormat=application/json&bbox=163049,467594,163256,468020,EPSG:28992,
this=0000018A6EA691F0) succeeds as WFS.
GDAL: GDALDriver::Create(ESRI
Shapefile,Perceel-aan.shp,0,0,0,Unknown,0000000000000000)
GDAL: GDALClose(WFS:
https://geodata.nationaalgeoregister.nl/aan/wfs?service=WFS&request=GetFeature&version=2.0.0&typename=aan&srsname=EPSG:4326&outputFormat=application/json&bbox=163049,467594,163256,468020,EPSG:28992,
this=0000018A6EA691F0)
GDAL: GDALClose(Perceel-aan.shp, this=0000018A70565430)
GDAL: In GDALDestroy - unloading GDAL shared library.

When I open the URL in the browser, I do get back JSON data. When I copy
this and paste it in QGis I get the expected polygons.

I'm not sure how I need to change my command to get this working for the
second server as well.
Please advice.

With regards,

Paul Meems
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20210618/fadb3e07/attachment.html>


More information about the gdal-dev mailing list