Geolibre (en AHN)

Ik probeer een workflow waarin ik Geolibre wil combineren met een online AHN‑bron. Het idee is om hoogte‑informatie direct te streamen in plaats van lokaal te downloaden, zodat iets makkeler te visualiseren is met online 3D tiles. Zonder de nodige Arcgis licenties :stuck_out_tongue:

Gebruikt iemand hier dit platform? Als ik een tif download en wil gebruiken krijg ik de melding terug dat een niet ondersteunde EPSG is Daarmee zou

Zijn er hier mensen die al ervaring hebben met het streamen van AHN‑data in combinatie met Geolibre? Tips over geschikte endpoints, performance‑ervaringen of valkuilen zijn welkom.

1 like

Ik ken GeoLibre niet, maar weet dat deze intern MapLibre gebruikt. Om AHN DEM data te gebruiken zal deze m.i. geconverteerd moeten worden. Naar GeoTIFF COG of “Terrarium Tiles”. In EPSG:3857 of EPSG:4326. EPSG:28992 lijkt mij (nog) niet.

Ik ben met het laatste bezig, work in progress voor https://map5topo.nl : via het project Mapterhorn genereer ik XYZ PMTiles in RGB Terrarium en serveer via Martin in Maplibre. Voor NL (AHN5+AHN6 50cm) draag ik PRs bij, o.a. Updates for Dutch DEM AHN data sources 5m and 0.5m - see issue #247 by justb4 · Pull Request #287 · mapterhorn/mapterhorn · GitHub .

Zie screenshots:

Als je BGT met de OGC api methode benaderd komt niet alles over heb ik het idee. Toch blijft het een leuk initiatief, zo’n webbased viewer.

(edit) doooohh. nu zie ik dat er een max features staat op : 1000

1 like

Topic neemt even een afslag hier…

Om mijn eerdere antwoord boven aan te vullen: Terrarium Tiles is dus een techniek/encoding om hoogte informatie naar clients zoals ook CesiumJS, te streamen. Een andere encoding (ipv RGB) is ‘quantized-mesh’ (preferred door CesiumJS).

Om topic-afslag te beantwoorden: een OGC API (Features) service zal altijd een geconfigureerd maximum aan Features teruggeven. Echter kent deze API ook twee andere faciliteiten: “Paging” en “BBOX Query”. Een “slimme client” zoals bijvoorbeeld QGIS combineert deze twee technieken om toch snel Features parallel binnen te halen en te renderen. Ik weet niet hoe/of MapLibre dit doet.

Geolibre heeft nog geen methode om the vullen op gebied of een boundary. Een verzoek neergelegd bij de maker

Wel kun je landxml als 3D formaat visualieren maar dit heeft zo z’n beperkingen. Mooiste zou zijn als je inderdaad AHN (ondersteuning nationale EPSG code) kunt gebruiken zodat je niet afhankelijk bent van onnauwkeurige wereldwijde diensten.

Fout melding AHN als WCS dienst

PDOK WCS https://service.pdok.nl/rws/ahn/wcs/v1_0 doet het nu (GeoLibre main branch), voornaamste probleem was dat DescribeCoverages (https://service.pdok.nl/rws/ahn/wcs/v1_0?SERVICE=WCS&VERSION=1.0.0&REQUEST=DescribeCoverage&COVERAGE=dtm_05m) alleen EPSG:28922 aangeeft bij SupportedCRSs, maar EPSG:3857/EPSG:4326 werken ook.

Nice! Nu is de kunst om dit als DTM weer te geven. Ik bedoel 3D.
Dat heb ik even gemeld.

De AHN via PDOK is versie 4. Het wordt spannender als je de WCS van AHN 5 en 6 uit de dataroom van AHN probeert. Deze serveert namelijk een andere indeling/formaat, dus ik ben wel benieuwd of dat lukt.

He Anton, ik heb even getest met https://service.pdok.nl/rws/actueel-hoogtebestand-nederland/wcs/v1_0 , werkt goed.

Ja, dat is dus het verwarrende, ook de actueel is gewoon AHN4 en wordt niet meer geupdate op PDOK. AHN heeft besloten om de data alleen nog via hun eigen dataroom te verspreiden.

GeoLibre Mapterhorn plugin (zonder WCS) - eerste versie

geolibre_mapterhorn.gif

Hi there, sorry to reply in english.

In case it is useful, I just wanted to post a link to an app specifically made for terrain, here is a terrain-viewer link to load the AHN raster DTM and DSM from WCS official endpoint, and compute the nDSM = DSM - DTM on-the-fly. The visualization modes docs is a useful place to better understand what modes exist, how to toggle them in 3D, 2D, load your own datasets, extract snapshots or COGs, comapre historical imagery etc. Hope this will be useful to the discussion!

1 like

This is truly an amazing initiative. Earlier this year someone pointed me to “Godzilla” from TU Delft. An inittiative to get fast DSM and DTM data for NL. Since it went offlineI tried to get in touch to learn the story behind it and why it ended, but i could not retrieve the right personsto get a clear answer. This feels very much like the sequel :slight_smile: Great! Can it do geotiff download of some area as well?

Godzilla AHN download (landsdekkend) - YouTube

Is it possible to supply this as a feature request? Geojson ánd Geotiff

Dooooohh!.. Nevermind. I Found the right button.

For the Dutch readers, I put together this short demo to highlight some of the viewer’s core functionality. It’s a seriously useful tool and well worth exploring.

@jo-chemla, could you explain how additional services can be added?

I tried loading a service from the Province of Wallonia, but it doesn’t appear in the viewer. I assume the correct service URL, including the required BBOX and coverage parameters, needs to be used, but I haven’t been able to figure out the correct configuration.

Relief de la Wallonie - Modèle Numérique de Terrain (MNT) 2021-2022 | Géoportail de la Wallonie


Any guidance or example would be greatly appreciated. Thanks in advance!

Thanks for the kind words!
Indeed DEM Export is accessible from the (usually folded) section, Download and Snapshot - hte export buttons in the drawings section are only for the drawn feature and layers. You can choose your DEM Geotiff export resolution, your screen resolution/viewport will be very fast - can even export all the other visualization modes layers like slope from terrain-analysis, relief visualization SVF etc, as raw values or color-ramped. Or you can choose a higher resolution, it will load more tiles, and can go from my tests up to 8k.

When you activate (pretty hidden, under source-info, last section, coverage overlays, mapterhorn), you see maptehorn combined source already includes the Relief de la Wallonie MNT, see screen
0.5 m · Relief de la Wallonie - Modèle Numérique de Terrain (MNT) 2021-2022 (Service public de Wallonie (SPW)) · “bewallonie”

It seems the terrain dataset is not served as raw elevation values, so it can only be added as an overlay.

Hope this helps! Note you can now even choose, in split-mode: overlay or side, to desync the visualization modes, for example to get a 3x2 map views, each with its own visualizatino mode - eg top row mapterhorn, bottom row another terrain, applying basemap on first column, hillshade on second, and openness on third column for example.