Kadastrale percelen genereren vanuit BRK

In deze blog beschrijf ik hoe je met QGIS “echte” vlakgerichte Kadastrale Percelen kunt genereren wanneer de BRK gegevens (de Kadastrale Kaart) alleen lijngerichte kadastrale grenzen en perceel punten blijkt te bevatten.

polygonized

Sinds 1 januari 2016 is ook de BRK (Basis Registratie Kadaster) beschikbaar via PDOK. Je kunt deze in GML formaat per provincie downloaden en direkt in QGIS openen.

Deze BRK dataset bevat gegevens over percelen. Elk perceel heeft een unieke ID en bij het kadaster ligt ook vast wie de eigenaar is en de historie. Hier heeft het Kadaster zijn bestaan en belangrijke rol in de samenleving aan te danken! De gegevens in de BRK zijn een momentopname van de perceelsgrenzen en bevat geen eigenaarsgegevens of historie.

Ik had zelf verwacht dat de BRK een polygoonbestand met perceelsgegevens zou bevatten. De ZIP-file bevat echter de bestanden Kadastralegrens.gml (lijnen), Perceel.gml (punten met perceelsgegevens), Annotatie.gml (punten) en Bebouwing.gml (lijnen). Hiervan zijn Annotatie.gml en Bebouwing.gml minder interessant (BAG en BGT bevatten immers al de omtrek van gebouwen en die labels kunnen we prima zelf binnen het vlak genereren) .

Het proces van het omzetten naar perceelvlakken met perceelgegevens bestaat globaal uit de volgende stappen:

1) inlezen van GML bestanden en opslaan als Shape bestanden
2) toevoegen van spatiale indexen aan Shape bestanden
3) omzetten van kadastrale grenzen naar kadastrale vlakken
4) toevoegen van de attributen van perceelpunten aan kadastrale vlakken

1) Download eerst de BRK van Flevoland (of andere provincie) van de downloadpagina:

https://www.pdok.nl/nl/producten/pdok-downloads/basis-registratie-kadaster/kadastrale-kaart

2) in dit voorbeeld pak ik de zip-file van Flevoland uit in de folder C:\data\brk\input\ en maak ik de folder C:\data\brk\output aan waarin ik de resultaatbestanden zal aanmaken.

3) Activeer de Processing Toolbox via het menu Processing >> Toolbox. Zet rechtsonderin QGIS de Render functie (Her)teken uit en zet de interface processing om van Simplified interface naar Advanced interface.
Render of and advanced interface

4) laad de gml-files Kadastralegrens.gml en Perceel.gml
inlezen_gml

5) Sla allebei de bijgeladen GML bestanden op in het shape-formaat, dit kan door deze in het paneel “Lagen” te selecteren en via het menu Kaartlagen >> Opslaan Als… op te slaan als shape-bestand. Dit kan ook via Processing. Zoek in de Processing Toolbox met het zoekveld naar het GDAL/OGR algoritme “Convert format” met convert. Je moet die van OGR hebben (voor omzetten van vectordata). Dubbelklik op het gevonden algoritme en de dialoog opent waarmee je kaartlagen kunt opslaan in verschillende bestandsformaten.

GDAL_convert

6) Zorg ervoor dat de aangemaakte shapebestanden zijn geopend in QGIS. Het toevoegen van een index, doe je via de Laag Eigenschappen. Dubbelklik op de geselecteerde kaartlaag KadastraleGrens en Druk op de knop [Ruimtelijke index aanmaken] (Create spatial index).
Create Spatial Index

7) Het aanmaken van perceelsvlakken uit perceelsgrenzen doen we vervolgens met de Processing Tool, Polygonize. Zoek in de Processing Toolbox met het zoekveld naar “Polygonize”. Dubbelklik op het gevonden algoritme en de  volgende dialoog opent. Wanneer je de optie Create geometry columns aanvinkt , dan worden de attribuutvelden oppervlakte en omtrek aangemaakt in de output layer KadastraleVlak.shp.
grens_naar_vlak

8) Tijdens de uitvoering verschijnt het logbestand van het omzetten van perceelsgrenzen naar vlakken, bij 40% staat deze erg lang stil en lijkt QGIS nergens op te reageren. (sterker nog het scherm van het processing algoritme kan tijdelijk zijn verdwenen). Dit is een proces wat echter veel tijd nodig heeft! Gewoon even geduld hebben, laat QGIS zijn werk doen.
log_polygonize

9) Gebruik het QGIS Processing Algoritme Join attribute by location, Voor het toevoegen van de attributen van Perceel.shp aan KadastraalVlak.shp om deze vervolgens op te slaan in KadastralePercelen.shp. Op de volgende dialoog door hierop te dubbelklikken en vul deze als volgt in. Dit algoritme is een stuk sneller en wanneer je hierna alleen de KadastralePercelen.shp laad en bekijkt, zie je Flevoland zoals getoond aan het begin van deze blog (wellicht in andere kleur). Vergeet ook niet de renderfunctie (Her)teken rechtsonderin de applicatie weer aan te zetten.
join_attributes_to_polygon

Wanneer je een kadastraal perceel selecteert en bekijkt heeft deze de volgende velden.
Eigenschappen Kadastraal Perceel

Door: Diethard Jansen (GIS-hulp)

17 thoughts on “Kadastrale percelen genereren vanuit BRK”

  1. Hallo Diethard Jansen, veel dank voor werkwijze om BRK.GML om te zetten naar shape. Echter bij het omzetten van line naar polygon krijg ik de foutmelding “No shapely geometry can be created from null value See log for more details”
    Misschien doe ik iets verkeerds via QGIS. Ik heb het probleem zelf opgelost door via ARCGIS “line naar polygon” te gebruiken en dat werkte wel.

    1. Dag Willem,

      Ik denk dat ik die fout heb gezien op het moment dat ik zelf probeerde het GML bestand te gebruiken als input voor de Polygonize tool. Wellicht goed om eerst alleen het aangemaakte shapebestand te openen in QGIS en dan de Polygonize tool te gebruiken. Laat mij anders weten welke versie van QGIS en welke provincie u heeft gebruikt.

  2. Dag Diethard,

    Ik krijg bij het polygoniseren van de shapefile van de kadastrale grens de foutmelding float division by zero.
    De shapefile van Perceel kan ik geeneens polygoniseren omdat deze is te kiezen als invoerlaag.

    Groeten Marc

  3. Hoi Diethard,

    Je reactie op Willem werkt. Ik had hetzelfde probleem en dat werkt nu. Waar ik zelf nog wel mee worstel is hoe de lagen te exporteren naar een autocad bestand. Als ik nu opslaan als, in kaartlagen, uitvoer en dan naar Autocad DXF, verdwijnt er een hoop informatie. Iemand een idee hoe dit te voorkomen?

    Groetjes,

    John Haan

    1. Alhoewel er in de laatste versie flink gesleuteld is aan de DXF-export, (zie: QGIS changelog 2.10 ) verlies je uiteraard veel data wanneer je van een echt GIS formaat (geometrie + gegevens) gegevens exporteert naar een CAD formaat. De Geometrie + symbologie blijft behouden maar de gegevens raak je helaas kwijt!

  4. Dag Diethard,

    Ik had zelf verwacht dat de BRK een polygoonbestand met perceelsgegevens zou bevatten

    Laat dat toch inderdaad het geval zijn !
    Perceel.gml bevat per perceel een vlak en een label (beiden!).
    Bij het inlezen van een GML bestand wordt er vooraf een analyse gedaan (perceel.gfs), waarbij er maar één geometrietype wordt geaccepteerd (POINT in dit geval).
    Vervang in het bestand perceel.gfs
    1
    door:
    3
    en lees het GML bestand opnieuw in.
    De vlakken worden nu waarschijnlijk in één keer verwerkt.
    Anders heb ik nog wel een iets andere perceel.gfs voor je.

    Zoals je zelf opmerkte kun je labels prima door QGIS laten toevoegen.

    1. Dag Fred,

      Ik heb in perceel.gfs het GeometryType gewijzigd van 1 naar 3. Als ik perceel.gml vervolgens opnieuw inlaad blijft het scherm leeg. Wat doe ik verkeerd?

      groet
      Sander

      1. Dag Sander,

        Ik had dezelfde vraag gesteld aan Fred en hij laat op zijn blog-site de inhoud van een werkende Perceel.gfs zien. Ik zal hier vanavond op deze site een nieuwe blog online zetten!

        1. Dag Fred en Diethard,

          Met het perceel.gfs bestand van de blog-site van Fred werkt het nu inderdaad wel. Super oplossing!

          groet
          Sander

  5. Ik ben een beginner op gebied van PDOK, QGIS en AutoCad. Ik heb AutoCad2012 tot mijn beschikking en kan het kadastrale GML bestand zodoende niet inlezen. Via QGIS wil ik de GML bestanden uit BRK omzetten naar DXF om in AutoCad verder te werken.
    Bij inlezen van Annotatie.gml in QGIS zie ik alleen punten en bij DXF-Export zijn ze in AutoCad niet aanwezig / niet zichtbaar.
    Voor mij volstaat de polygoon informatie uit BRK, echter ook heb ik de leesbare Annotatie gegevens in AutoCad nodig.
    Hoe krijg ik dat voor elkaar?

    1. Binnen QGIS heb ik het voor elkaar: annotatie / waarde van label “tekst” is zichtbaar in gewenste weergave en de (invoeg)punten zijn uitgeschakeld.
      Probleem blijft dat de annotatie / label “tekst” in AutoCad niet zichtbaar is na DXF export uit QGIS.
      Hoe zorg ik ervoor dat ze wel zichtbaar zijn in ACAD resp met de DXF-export meekomen?

  6. Kan iemand de omgezette bestanden ergens wegzetten? Ik heb met name de kadastrale polygonen van Gelderland en Noord-Brabant nodig.

    Ik heb zelf meerdere pogingen gedaan, maar ook na een nacht ‘polygonizeren’ blijft het proces op 40% staan. Bovendien loopt mijn laptop vast, waardoor het proces telkens wordt verstoord.

    Ik hoor graag.

    1. Dag Jildert,

      Excuses voor het late antwoord.
      Een interessant idee om de bestanden om te zetten, wellicht dat iemand dit oppakt. Wanneer het genereren van Kadastrale percelen niet lukt, lees dan vooral het vervolgartikel “Kadastrale percelen en het gfs bestand” wat zeker werkt! Sterker nog ik heb even zitten twijfelen of ik deze blog niet moest verwijderen, want het geeft wel goed weer hoe je met processing vlakken kunt genereren uit lijnen, maar het aanpassen van het gfs bestand werkt veel beter!

Laat een reactie achter bij Diethard Jansen Reactie annuleren

Je e-mailadres wordt niet gepubliceerd. Verplichte velden zijn gemarkeerd met *

De volgende HTML tags en attributen zijn toegestaan: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>