Generate parcels areas from parcel boundaries

Hi! In this blog I describe how you can create proper parcels with polygon geometry in from polylines (parcel boundaries) and points (Parcel point with parcel attributes placed inside parcel boundaries).

polygonized

Since the 1st of januari 2016 a dataset named, BRK (Basis Registratie Kadaster) is available from PDOK. You can download these in GML format for each province and open directly in QGIS.

The BRK holds information on parcels, land areas. Every parcel is registrated with an unique ID at Kadaster, the Dutch Land Registration authority. They also register who owns the parcel of land and the historie of previous owners. De dataset BRK provides a momentary view of the parcels at a recent point in time.

I expected to receive polygons but the BRK included only lines and points. Fortunately, QGIS provides the tools to provide what I want.

The process of converting parcel boundaries with parcel points to parcel polygons includes globally following steps:

1) Read parcels provided in GML and save as Shape
2) Add spatial indexes to shape files
3) Convert parcel boundaries (lines) to areas (polygones)
4) Add attributes from parcel points to areas

1) First Download the dataset BRK of the Flevoland Province from the download webpage:

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

2) In this case, I unpack the downloaded zip-file in folder C:\data\brk\input\ and create the folder C:\data\brk\output for (intermediate) results.

3) Activate the Processing Toolbox from the menu Processing >> Toolbox. Turn rendering of in QGIS (right below in QGIS application window) and switch the interface of processing from Simplified interface to Advanced interface.
Render of and advanced interface

4) load both gml-files Kadastralegrens.gml and Perceel.gml
inlezen_gml

5) Save both GML files loaded to shape-format, in het layer panel you can select a layer and select in het menu Layers >> Save As…. Using the Processing Toolbox you can use Convert format (Translate format) [OGR] from the GDAL/OGR group of algorithms. When you double click on the algorithm the dialog opens that can be used to convert to many vector formats.
GDAL_convert

6) Make sure the created shapefiles are opened in QGIS first. To add a spatial index double click on a selected layer which opens the Layer Properties dialog. Use the button [Create spatial index] .
Create Spatial Index

7) To create parcel areas from boundaries, use the Processing Algorithm Polygonize. Search in the Processing Toolbox with “Polygonize”.
grens_naar_vlak

8) This algorithm takes a lot of time! At 40% it seems to stop, but please be patient. Let QGIS do its job.
log_polygonize

9) Now add the attributes of the point layer Perceel.shp to created area layer KadastraalVlak.shp and I saved it to KadastralePercelen.shp.
join_attributes_to_polygon

When you select a parcel using the identify tool you should see following fields.
Eigenschappen Kadastraal Perceel

That’s it.

Regards, Diethard Jansen (GIS-hulp)

15 thoughts on “Generate parcels areas from parcel boundaries”

  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.

Leave a Reply

Your email address will not be published. Required fields are marked *

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>