België en Duitsland in RD met QGIS

De Basisregistratie Topografie (BRT/Top10NL) houdt ongeveer een kilometer na de landsgrens al op. Velen, bijvoorbeeld brandweermensen, hebben behoefte aan meer kaartbeeld over de landsgrens. Die kaartdata is er, en nog kosteloos ook. De Open Street Map (OSM) van de Europese landen is als Shapefiles te downloaden, o.a. via de Duitse GeoFabrik, http://download.geofabrik.de/europe.html. Ze wordt aangeboden in WGS84 (latitude, longitude). QGIS 1.8 biedt een eenvoudige mogelijkheid om die te converteren naar het Rijksdriehoeksstelsel (EPSG:28992): maak een kopie van de Shapefile onder een andere naam (eventueel alleen het relevante geselecteerde gebied langs de grens “selectie opslaan als…”), en kies daarbij voor het geselecteerde CRS “Amersfoort New/RD”. QGIS maakt een nieuwe shapefile aan met het gewenste RD-coördinaatstelsel. Die Shapefile kun je toevoegen als kaartlaag in je QGIS RD-projectfile. Die Shapefile sluit aan de landsgrens mooi aan op de BRT kaartlagen.
Voorbeeldje van de stad Emden, bij Delfzijl:

OSM beeld van Emden, geconverteerd naar RD stelsel, en opgemaakt in het QGIS JW kleurenschema.

Voor de kleurenopmaak heb ik gebruikt gemaakt van de kleurcodes die ik eerder op dit forum heb gedeeld. De OSM versies (QML formaat, voor in QGIS) vind je hier. Per OSM sublaag (natural, landuse, buildings, roads, etc) vind je een stijlbestand. Let wel, de wegen heb je minimaal 2x nodig: één keer voor de zwarte omranding, en één keer voor de kleuren daar bovenop. En omdat de OSM roads ook een “tunnel” indicator hebben, kun je die ook als aparte kaartlaag toevoegen met een eigen QML opmaak.
Natuurlijk is OSM vooral gericht op wegen, en minder op terrein. Mij valt op dat er in de OSM terreinvlakken nogal wat ‘gaten’ zitten, afhankelijk per landsgebied, provincie of zelfs gemeente. Zo heeft Emden z’n zaakjes goed voor elkaar, maar de buurgemeenten hebben veel minder terreinvlakken beschikbaar (zie links bovenin bij het voorbeeld).
 

QGIS en WFS caching

Via de PDOK-servicesplugin is het mogelijk om de verschillende PDOK kaartservices te gebruiken.
WFS is op zich een interessante service, echter wanneer je een laag als ‘wfs nwb wegvakken’ aanroept in QGIS via de PDOKservicesPlugin krijg je het volgende resultaat:

Het (in eerste instantie) vreemde resultaat verkrijg je ook wanneer je ingezoomed bent op een gemeente. Je zou verwachten dat je dan de features krijgt binnen de huidige Extent/BoundingBox, maar je blijkt toch weer diezelfde balk te krijgen….
Maar wij kunnen / QGIS kan dit fixen.
Continue reading “QGIS en WFS caching”

BAG gebruiksfuncties op Top10NL in QGIS

De pandcontouren van de BAG vormen een elegant alternatief voor de gebouwvlakken van de Top10NL, het verschil in actualiteit van beide bronnen daargelaten. Met de BAG is het mogelijk om in een stad snel de locaties van winkelcentra, uitgaanscentra (“bijeenkomst”), kantoorparken en industriegebieden zichtbaar te maken op de Top10NL. Maar de gebruiksfuncties zitten niet in een geometriekolom van de BAG (PostGIS formaat), omdat een verblijfsobject meerdere gebruikersfuncties kan hebben (1:n relatie). Via de gratis BAG converter van Geon is het echter mogelijk om een shapefile van verblijfsobjecten te maken waarin per verblijfsobject twee gebruiksfuncties als attribuut worden meegenomen.
In QGIS is deze punt-shapefile (vaak de voordeur van een pand, maar niet altijd) te plotten bovenop de BAG pandcontouren die op de Top10NL zijn geplakt. Hieronder een voorbeeld van het centrum van Almere, met de volgende legenda van gebruiksfuncties (exclusief ‘woonfunctie’, anders staat de kaart vol met stipjes):

  • Groen: winkelfunctie
  • Paars: bijeenkomstfunctie
  • Geel: onderwijsfunctie
  • Roze: gezondheidsfunctie
  • Blauw: kantoorfunctie
  • Bruin: industriefuntie
  • Groenblauw (cyaan): sportfunctie
  • Rood: logiesfunctie
Almere Top10NL + BAG pandcontouren + BAG VBO verblijfsfuncties
QGIS screenshot: Almere Top10NL + BAG pandcontouren + BAG VBO verblijfsfuncties

De volledige plaat met legenda is hier te zien.
De Topografische kaart geeft zo snel inzicht in waar de winkels zijn, waar het uitgaansgebied is, waar de kantoorparken zijn, waar de scholen zijn etc. Een volgende stap zou zijn om in QGIS een heel pandvlak te kleuren in plaats van alleen een stipje. Dat vergt echter wat complexere attribuut-joining berekeningen. Wizards met een slimme oplossing hiervoor zijn van harte uitgenodigd om dat hier te melden!
 

Nieuwe versie 0.4 PDOK-servicesplugin

Ik heb net een nieuwe versie gemaakt van de PDOK-servicesplugin. Deze plugin is beschikbaar via de standaard plugin-repository in QGIS. Laad de plugin door via het menu the kiezen voor ‘Plugins/Python plugins ophalen’ en dan te zoeken met de term ‘pdok’. Zie verder ook het eerdere artikel Eenvoudige PDOK-services plugin voor QGIS. Sinds de laatste update van de PDOK-services werkte die namelijk niet meer.
Deze nieuwe update van de PDOK-services voegt een paar leuke nieuwe kaartlagen toe, bijvoorbeeld de fietsknooppuntroutes in Nederland. Daarnaast zijn er echter wat technische veranderingen doorgevoerd, waardoor de manier van laagnamen in de Capabilities niet meer strookte met het systeem wat hiervoor in de lucht was. Voor de liefhebbers: in de Capabilities stonden voorheen laagnamen als “natura2000”, terwijl de laag moest worden opgehaald met “natura2000:natura2000”. Dit is nu gefixed.
Er zijn natuurlijk altijd een paar mogelijke verbeterpunten toe noemen. Ik ventileer ze hier maar:

  • Fietsknooppunten: Persoonlijk zou ik iets meer liefde aan de visualisatie toevoegen 🙂 De fietsknooppunten worden nu in een soort neon-groen getoond, met de labels neongroen met een witte halo. Het zou mooier (en nog niet eens zo heel veel werk) zijn als dat gewoon het groen zou zijn van de knooppuntbordjes zelf, en de labels groene tekst in een groenomrand wit cirkeltje (net als de bordjes).
  • Alle CBS services: waarom hebben alle lagen van het cbs, zoals de gemeenten, provincies, wijken etc eigenlijk een vulling? Die zullen toch nooit als ondergrond worden gebruikt, maar veel meer als deel van een meer uitgebreider kaart?

  • TopXRaster-services: elke top-kaart is in een aparte (schaalafhankelijke) service opgenomen. Met als gevolg dat iemand die op 1:50.000 niveau de Top25Raster service laadt NIETS ziet. Waarom niet een van de volgende mogelijkheden:
    1. Een (extra?) service die je TopRaster noemt, en waarin alle TopRaster-kaarten (schaalafhankelijk) zitten. Zodat je tenminste in het gebruik op alle zoomniveau’s een kaart hebt.
    2. Als minder alternatief: laat desnoods een lichte contour van Nederland of de provincies zien op die niveau’s dat je nu niets ziet.

En voor de plugin zelf: op dit moment gebruikt de plugin een kopie van een JSON-bestand die alleen opgehaald kan worden door mensen die in een whitelist staan (gemeenten etc). Het zou natuurlijk veel mooier zijn als zo’n bestand ergens openbaar zou staan, zodat de plugin live aanpassingen meteen zou tonen.
Gelukkig: er is nog genoeg werk dus 😉

Proportionele symbolen (2)

In het artikel Gebieden inkleuren werd de waarde van een attribuut met QGIS gepresenteerd door gebieden in te kleuren. Voor relatieve gegevens, zoals percentages of verhoudingsgetallen, is dat een goede optie. Voor absolute gegevens, zoals inwonertal, zijn symbolen die in grootte variëren geschikter. Beide kaartvormen kunnen overigens ook goed met elkaar worden gecombineerd. Dit is een herziene versie van een ouder artikel.

Methoden

QGIS 1.8 bevat vier mogelijkheden om proportionele symbolen weer te geven. Die zijn allemaal te vinden in het dialoogvenster ‘Layer properties’, op verschillende tabbladen. De volgende tabel geeft de voor- en nadelen weer. ‘Old symbology’ wordt daarbij buiten beschouwing gelaten. Continue reading “Proportionele symbolen (2)”

Video: hoe maak je razendsnel kaarten met de PDOK Plugin en je eigen data in Qgis

Op Eenvoudige PDOK services plugin voor QGIS beschreef Richard al hoe je met de PDOK plugin de services van PDOK snel in Qgis kan gebruiken.

Deze video http://www.youtube.com/watch?v=NjLG4KtpdYk demonstreert hoe je met deze plugin in enkele minuten eigen data (GML in dit geval) kan combineren met PDOK services.

Eenvoudige PDOK services plugin voor QGIS

Om de PDOK-services te kunnen bekijken is er een eenvoudige plugin beschikbaar via de standaard plugin-repositories in QGIS. Laad de plugin door via het menu the kiezen voor ‘Plugins/Python plugins ophalen’ en dan te zoeken met de term ‘pdok’.
Met de plugin krijgt u een nieuw scherm met een lijst met de beschikbare services die PDOK aanbiedt. Met een klik op een regel laadt u die service in QGIS. Het is ook mogelijk om in de lijst te zoeken door het filter te gebruiken.

Continue reading “Eenvoudige PDOK services plugin voor QGIS”

Python one-liners for adding wms, wmts, wfs layers

In the Netherlands, currently there is a movement to make more and more data available for the general public. This is done via different services: WMS, WMTS and WFS services via a service called PDOK (Publieke Data Op de Kaart).
While creating a plugin to list and load all those layers in QGIS I had to dive into the details of the different api-possibilities to add raster and vector layers to the map via the Python inteface.
Because I think this is generally helpfull, here is some code to try.

Continue reading “Python one-liners for adding wms, wmts, wfs layers”

Iets als lyr stijl bestand voor QGIS

Afgelopen weekend was er een posting op de gebruikers mailinglist met de vraag of er in QGIS iets dergelijks was als het .lyr bestand in de ESRI wereld.
Hoewel ik wel wist dat je met QGIS stylen kunt opslaan en hergebruiken (ze worden dan opgeslagen als .qml bestanden). Was ik me niet bewust van het feit dat ze zich OOK gedragen als de .lyr files in de zin dat als er een .lyr file staat naast een shape file met exact dezelfde naam, dat QGIS dan automatisch die stijl laadt en toewijst aan de laag.
Altijd mooi zo’n ‘leermomentje’ als begin van de week, bedankt Morten.

BRT

Het Kadaster levert sinds kort de Basisregistratie Topografie (BRT) in gml-formaat uit. In QGIS kun je deze bestanden direct inlezen en gebruiken. Hieronder een beschrijving.
Klik op Vectorlaag Toevoegen (ctrl-shift-v) en blader naar het gml-bestand. Zorg dat rechtsonder het bestandstype op [GML] of op “Alle bestanden” staat. Selecteer alle lagen die QGIS in het gml-bestand vindt en open deze. De lagen worden nu met een standaard (random) symbool toegevoegd aan je kaart.

QGIS-scherm direct na openen gml-bestand

Continue reading “BRT”