In 3-voud (QGis, FME en Postgis): hoeveel punten in elke polygoon?

Voor een boek waar Travelingo mee bezig is moesten we het aantal POIS (points of interest) per Nederlandse provincie weten. Meerdere tools zijn dan mogelijk. Er bestaan er ongetwijfeld veel meer, maar de makkelijkste zijn voor mij QGIS, FME en Postgis.

“In 3-voud (QGis, FME en Postgis): hoeveel punten in elke polygoon?” verder lezen

Labels workflow in QGis and Inkscape

However matured the labeling engine in Qgis is, there is still a moment when you think beyond QGis in order to position the labels and icons as you want. I managed to setup a workflow with QGis and Inkscape with satisfactory results. Both are opensource. But if anybody has better options I would like to hear from you!

“Labels workflow in QGis and Inkscape” verder lezen

Clustering van bushaltes in QGis met postgis

Bushaltes bevinden zich vaak in paren, aan beide zijden van de weg. Op een kaart is het vaak beter om daar één symbool met één label van te maken. Dat kan met postgis en de DBManager in Qgis.

Soms ook gaat het om meer dan 2 bushaltes met dezelfde namen. Kijk hieronder naar het plein voor NS-treinstation Zwolle. Het wemelt van de bushaltes.

bushalte noncentroids

Beter zou zijn om haltes met dezelfde naam te clusteren. Dat kan door een laag aan te maken in de DBManager van Qgis. Niet in de Postgresql  database zelf! Het is beter een laag met de sql-code te koppelen, omdat de st_makeenveloppe functie 4 waarden heeft, corresponderend met de kaartextent, die per project verschillen. Het gaat om de volgende SQL-code:

select name, st_centroid(st_union(way)) as waygeom
from planet_osm_point
where way && ST_MakeEnvelope(202248, 501898, 205180, 504072, 28992) and highway = 'bus_stop'
group by name;

De st_union wordt gedaan met behulp van de group by functie, waarbij records worden samengevoegd tot multipoints, gebaseerd op dezelfde naam.

P.S. De ST_MakeEnveloppe functie beperkt het samenvoegen op naam tot de extent van de kaart in Qgis. Immers, een bushalte met ‘Markt’ als naam kan overal in Nederland voorkomen en tot vreemde centroids leiden.

Daarna doet de st_centroid functie de multipoints weer tot één punt smelten. Het resultaat, zie hieronder. Een stuk mooier.

bushalte centroids

Ook een gewoon bushalten’paar’ oogt beter als één symbool. Zie het verschil hieronder.

bushalte noncentroids bushalte centroids

Neutraal Moresnet

Een interessant boek over Neutraal Moresnet. het vroegere ministaatje ten zuiden van Vaals, leidde ertoe dat ik de grenzen van Moresnet eens precies ‘in kaart’ wilde zien. Een shapefile is niet zo gemakkelijk te vinden. Wat wel te achterhalen is, zijn de grensstenen, die nog steeds in het veld staan en die op een kaart weergegeven, de grens prima aangeven. “Neutraal Moresnet” verder lezen

Hillshading, Ahn2 en FME

AHN2 en hillshading

Hillshading in FME blijkt heel makkelijk. De RasterHillShader transformer is voldoende. Gebruik als Reader een raster-TIF van de PDOK AHN2 maaiveld 5 meter hoogtebestanden. De default parameters zijn  voldoende. Importeer de TIF Writer in de bovenste laag van QGis met blending mode op ‘multiply’.

Het is maar een klein deel van Zuid-Limburg, zoals je ziet. De kaart zelf is aangemaakt met behulp van OSM en Top10NL. Het is een ongoing project, maar ik vind het er al aardig uit zien…. Het oogt wel eigenaardig dat de Nederlandse Hoogtebestanden bij de grens ook onmiddellijk ophouden. De heuvels aldaar niet …

Fietsknooppuntennetwerk van OSM en LFP vergeleken

Om kaarten te maken met de fietsknooppunten in Nederland en de fietsknooppuntlijnen daartussen, zijn er twee bronnen die je als appmaker of cartograaf kunt gebruiken.

1. Openstreetmap (OSM)
2. Het Landelijk Fietsplatform (LFP, verantwoordelijk voor het beheer van o.a. de fietsknooppuntnetwerken)

“Fietsknooppuntennetwerk van OSM en LFP vergeleken” verder lezen

Schema naam aanpassingen van QGis projectfiles via XSLT

Dankzij de dames en heren van NLExtract is het sinds vorig jaar mogelijk een Postgresql backup te krijgen van o.a. de BAG en de Top10nl kaarten. Via een restore heb je dan weer de meest actuele cartografische gegevens. Klein ongemak is dat bij de restore de databases in een apart schema worden gezet (‘bagactueel’ voor de BAG en ‘ttnl´ voor de top10nl), maar dat de QGis projectfiles die ik de afgelopen anderhalf jaar heb in elkaar geknutseld, een verwijzing bevatten naar het standaard ‘public’ schema.

“Schema naam aanpassingen van QGis projectfiles via XSLT” verder lezen

Generate a polygons QGIS layer to use as a Atlas coverage layer’

Proposition: B&Bs, Hotels and other potential clients are looking for simple A4 printed maps of their region, including walking- and bicycling routes, tourist attractions, pubs, but also things like medical provisions, shops etc.. The clients location can be used as a focal point and there should be 4 maps for each block surrounding the focal point (NW, SW, NE and SE), but also a 5th central map for general use.

Solution: If I can make 5 polygons in one QGIS layer, I could use this as a coverage layer for generating an QGis Atlas. See script below…

“Generate a polygons QGIS layer to use as a Atlas coverage layer’” verder lezen

Antipode van Nederland dicht bij Nieuw Zeeland

In 1997 vlogen wij voor het eerst naar Nieuw Zeeland. Voor een lange fietstocht van drie maanden. Toen pas besefte ik hoe ver dat is. De halve aardbol rond gaat je als passagier niet in je koude kleren zitten, wat een lange zit! Het zal in die jaren zijn geweest dat ik benieuwd was waar je zou uitkomen als je echt door het middelpunt van de aarde zou gaan. In de geografie is dit de antipode van een punt: het tegenoverliggende punt op de aardbol. De afgelopen dagen heb ik daar een QGis-kaart van gemaakt. Dat ging niet zonder slag of stoot…

“Antipode van Nederland dicht bij Nieuw Zeeland” verder lezen