Waterschapsgrenzen download FOUT, WMS service goed

Datasets
#1

Ik wil graag de waterschapsgrenzen gebruiken.

De dataset is te vinden in het nationaal georegister: http://www.nationaalgeoregister.nl/geonetwork/srv/dut/catalog.search#/metadata/2d4ec7e6-50ae-4c71-bd20-5b9211e2135f

Als ik de service gebruik (QGIS) zie ik: image
te controleren met: http://geodata.nationaalgeoregister.nl/wsaeenhedeninspire/wms?language=dut&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=10000,300000,280000,625000&CRS=EPSG:28992&WIDTH=227&HEIGHT=273&LAYERS=waterschapsgrens&STYLES=&FORMAT=image/png&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE

Maar als ik de download gebruik (http://geodata.nationaalgeoregister.nl/wsaeenhedeninspire/extract/Waterschapsgrenzen_4258.zip) dan krijg ik iets heel anders als ik dit inlees in QGIS:
image

Er lijkt hier iets mis te gaan met de gekoppelde download, ook gezien de omschrijving op de pagina: image
Hier staat dat deze data in RD is, terwijl het in EPSG:4258 is.

Wie kan mij helpen aan de echte waterschapsgrenzen als data (WFS mag, enig geografisch data bestand ook, liever niet de data zoals door ESRI gemaakt ivm aanpassingen die zij uitvoeren).

#2

Beste Marco,

Ik heb de zelfde data gedownload van PDOK en in QGIS 3.10.2 gesleept, met een netjes resultaat.

image
image1295×796 156 KB

Wellicht dat er iets anders aan de hand is.

groet, Erkan

#3

Dank voor je hulp.

Hetzelfde nog een keer gedaan, en weer die vierkanten.

Ook ik gebruik 3.10.2 … Er is iets geks aan de hand.

#4

Als ik met ogr2ogr omzet naar een shapefile, dan is het resultaat wel goed. Dit is echt heel gek …

#5

Het GML-bestand bevat per feature ook een bounding box. (<gml:boundedBy><gml:Envelope …). Klaarblijkelijk wordt die door je QGIS versie gebruikt i.p.v. de daadwerkelijke geometrie. Wat zou kunnen is dat jij het XSD-bestand niet in dezelfde directory hebt staan, maar Erkan dat wel heeft. In dat geval wordt een GFS-bestand aangemaakt en gaat GDAL zelf bepalen wat de geometriekolom is. Wanneer het XSD-bestand wel aanwezig is, gebeurt dat niet. Tenminste niet bij mij, maar ik zit nog steeds op QGIS 3.8.0. Overigens zie ik wel de juiste polygonen als de XSD niet aanwezig is.

Zou jij de inhoud van het GFS-bestand willen posten, indien die aangemaakt wordt? Welke versie van GDAL wordt bij jouw QGIS-versie meegeleverd? Is dat GDAL 2 of 3? En welke versie gebruik jij met ogr2ogr? Mogelijk is dat de oorzaak van het verschil. En is bij ogr2ogr de XSD wel aanwezig, of wordt daar ook een GFS-bestand aangemaakt?

1 like
#6

Ha Frank,

MET xsd wordt er inderdaad geen .gfs bestand aangemaakt. De bounding boxen worden dan weergegeven zoals op eerder geposte screenshot.

ZONDER xsd word er wel een .gfs bestand aangemaakt. Helaas worden dan ook de bounding boxen weergegeven zoals op eerder screenshot.

Bij ogr2ogr gebruikte ik versie 1.11. Er werd dan geen .gfs bestand aangemaakt (voor zover ik mij kan herinneren; inmiddels gdal geupgrade, dus ik kan dat niet meer reproduceren).

De volledige QGIS versie info zoals ik die gebruik:


QGIS version
3.10.2-A Coruña
QGIS code revision
616ad4531b
Compiled against Qt
5.9.5
Running against Qt
5.9.5
Compiled against GDAL/OGR
3.0.2
Running against GDAL/OGR
3.0.2
Compiled against GEOS
3.8.0-CAPI-1.13.1
Running against GEOS
3.8.0-CAPI-1.13.1 
Compiled against SQLite
3.22.0
Running against SQLite
3.22.0
PostgreSQL Client Version
10.10 (Ubuntu 10.10-0ubuntu0.18.04.1)
SpatiaLite Version
4.3.0a
QWT Version
6.1.3
QScintilla2 Version
2.10.2
Compiled against PROJ
6.2.1
Running against PROJ
Rel. 6.2.1, November 1st, 2019
OS Version
Ubuntu 18.04.3 LTS
Active python plugins
maplibrary; 
qgisnetworklogger; 
QGIS_QuickPrint; 
plugin_reloader; 
RuimtelijkePlannen; 
pdokservicesplugin; 
pdok_locatieserver_locator_filter; 
BGTImport; 
QuickWKT; 
QGIS-NL-base; 
MetaSearch; 
processing; 
db_manager

De .gfs zoals gegenereerd wanneer de xsd niet wordt gebruikt:

<GMLFeatureClassList>
  <GMLFeatureClass>
    <Name>wsh.waterschapsgrens</Name>
    <ElementPath>wsh.waterschapsgrens</ElementPath>
    <!--MULTIPOLYGON-->
    <GeometryType>6</GeometryType>
    <SRSName>GEOGCS["ETRS89",DATUM["European_Terrestrial_Reference_System_1989",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4258"]]</SRSName>
    <DatasetSpecificInfo>
      <FeatureCount>21</FeatureCount>
      <ExtentXMin>50.75037</ExtentXMin>
      <ExtentXMax>53.49892</ExtentXMax>
      <ExtentYMin>3.35838</ExtentYMin>
      <ExtentYMax>7.22750</ExtentYMax>
    </DatasetSpecificInfo>
    <PropertyDefn>
      <Name>gag_id</Name>
      <ElementPath>gag_id</ElementPath>
      <Type>String</Type>
      <Width>16</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>admin_code</Name>
      <ElementPath>admin_code</ElementPath>
      <Type>String</Type>
      <Width>4</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>land_code</Name>
      <ElementPath>land_code</ElementPath>
      <Type>String</Type>
      <Width>2</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>waterschapsnaam</Name>
      <ElementPath>waterschapsnaam</ElementPath>
      <Type>String</Type>
      <Width>36</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>laatste_wijziging</Name>
      <ElementPath>laatste_wijziging</ElementPath>
      <Type>String</Type>
      <Width>10</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>hierarchieniveau</Name>
      <ElementPath>hierarchieniveau</ElementPath>
      <Type>Integer</Type>
    </PropertyDefn>
    <PropertyDefn>
      <Name>hierarchienaam</Name>
      <ElementPath>hierarchienaam</ElementPath>
      <Type>String</Type>
      <Width>25</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>inspire_id_p</Name>
      <ElementPath>inspire_id_p</ElementPath>
      <Type>String</Type>
      <Width>47</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>inspire_id_l</Name>
      <ElementPath>inspire_id_l</ElementPath>
      <Type>String</Type>
      <Width>47</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>inspire_id_v</Name>
      <ElementPath>inspire_id_v</ElementPath>
      <Type>String</Type>
      <Width>47</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>publicerende_waterbeheerder</Name>
      <ElementPath>publicerende_waterbeheerder</ElementPath>
      <Type>Integer</Type>
    </PropertyDefn>
    <PropertyDefn>
      <Name>wbh_code_objectidentificatie</Name>
      <ElementPath>wbh_code_objectidentificatie</ElementPath>
      <Type>Integer</Type>
    </PropertyDefn>
    <PropertyDefn>
      <Name>einde_levensduur</Name>
      <ElementPath>einde_levensduur</ElementPath>
      <Type>String</Type>
      <Width>10</Width>
    </PropertyDefn>
    <PropertyDefn>
      <Name>sde_id</Name>
      <ElementPath>sde_id</ElementPath>
      <Type>Integer</Type>
    </PropertyDefn>
  </GMLFeatureClass>
</GMLFeatureClassList>