Maaiveldlijn uit AHN voor dwarsprofielen

Beste Topic lezer,
In de AHN viewer zoeken wij regelmatig hoogtes van enkele punten op. Dit werkt niet echt handig, omdat je met de muis naar het gewenst XY coördinaat moet zoeken en het aanklikken. Bovendien willen we graag een lijn uitzetten en van alle punten op die lijn de Z coördinaten ontvangen in .csv of .dxf.
Sinds kort hebben we QGIS geinstalleerd en de PDOK plugin, waarmee we de AHN gegevens zichtbaar kunnen maken. Maar het lukt nog niet de gewenste info eruit te halen. De tutorials die we hierover gevonden hebben, zijn erop gericht maaivelden in 3D dxf bestanden om te zetten.

Weet iemand een manier om specifiek de Z coördinaten van een lijn te exporteren naar .csv of .dxf?
De lijn is meestal een rechte lijn, maar liever nog zouden we ook van een lijn opgebouwd uit rechte en kromme stukken de Z-coördinaten binnen kunnen halen.

We gebruiken deze dan om dwarsprofielen van de bodem te maken, met daarin de maaiveldhoogtes. Dit is voor ons vooral van belang wanneer zich in het dwarsprofiel sloten of wegen bevinden.

In de AHN-viewer kan je het resulterende hoogteprofiel exporteren als csv, waarbij je dan een lijst krijgt met coördinaten (X,Y,Z). Wellicht voldoet dit? zie afbeelding:

Als je op zoekt bent naar mogelijkheden om eenvoudiger een lijn in te tekenen op basis van vooraf bekende coördinaten, dan kan je beter gaan kijken naar professionele GIS-software (commercieel of open source). Of een alternatief kan nog zijn om direct een hoogteprofiel op te vragen door de API te gebruiken die in de AHN-viewer wordt gebruikt. Als je de coordinaten van de lijn waarvoor je het profiel wilt berekenen al weet, kan je die direct sturen naar de API, zie voorbeeldje:

https://ahn.arcgisonline.nl/arcgis/rest/services/Geoprocessing/Profile_AHN3/GPServer/Profile/execute?InputLineFeatures=%7B%22fields%22%3A%5B%7B%22name%22%3A%22OID%22%2C%22type%22%3A%22esriFieldTypeObjectID%22%2C%22alias%22%3A%22OID%22%7D%5D%2C%22geometryType%22%3A%22esriGeometryPolyline%22%2C%22features%22%3A%5B%7B%22geometry%22%3A%7B%22paths%22%3A%5B%5B%5B123487%2C454666%5D%2C%5B123520%2C454588%5D%5D%5D%2C%22spatialReference%22%3A%7B%22wkid%22%3A28992%7D%7D%2C%22attributes%22%3A%7B%22OID%22%3A1%7D%7D%5D%2C%22sr%22%3A%7B%22wkid%22%3A28992%7D%7D&ProfileIDField=OID&DEMResolution=FINEST&MaximumSampleDistance=&MaximumSampleDistanceUnits=Meters&env%3AoutSR=28992&env%3AprocessSR=&returnZ=true&returnM=true&returnTrueCurves=false&returnFeatureCollection=false&f=pjson

Beste JBak, dit komt inderdaad aardig in de buurt, maar:

  1. Het handmatig met de muis selecteren van de XY coördinaten is te tijdrovend. Vandaar dat we al de stap gemaakt hebben naar QGIS. Echter hebben we de juiste tutorial nog niet kunnen vinden om hier goed mee te kunnen werken
  2. de maaivelden die we op deze manier vanuit de viewer gedownload hebben, zijn niet gecorrigeerd. Dus ongeacht welke layer (ruw/gefilterd, dynamisch/statisch, AHN2/AHN3) geselecteerd zijn, blijft het hoogteprofiel dezelfde resultaten leveren. Vaak zien we dat bij bomen rare pieken optreden. De verwachting is door vanuit QGIS te werken, wel de juist gefilterde AHN bestanden te kunnen raadplegen.
  3. Op deze manier vanuit de viewer zijn geen krommen in te tekenen.

We zijn dus nog op zoek en staan open voor verdere suggestie.

Beste JBak, heb jouw bericht nog eens goed gelezen, nu vooral het derde alternatief. Echter lijkt de link niet te werken, of snap ik het gewoon niet. Ik krijg via de link een stuk code te zien, maar zie niet hoe ik daarmee tot een .csv of eventueel een .dxf kom.
Kun je het toelichten?

Een ontwikkelaar kan via de API een request sturen en zo een json response krijgen wat de ontwikkelaar vervolgens bijv. zelf weer verwerkt / opslaat / mee verder werkt in een maatwerk applicatie / tool. Volgens mij is dit niet iets om verder naar te kijken gezien je vraag.

Terug naar je behoefte om Z-waarden te krijgen voor een lijn (in uiteindelijk CSV of DXF formaat). Hoe dat met QGIS gaat kan ik je niet vertellen maar wellicht iemand anders op dit forum. Wel kan ik je meegeven hoe je zoiets binnen ons product ArcGIS zou kunnen doen met de Stack Profile functie: Stack Profile (3D Analyst)—ArcGIS Pro | Documentation

De stappen in het kort zijn:

  • Voeg de AHN2 of AHN3 kaart vanuit ArcGIS Online toe aan ArcGIS Pro (je hoeft dus geen data lokaal te hebben)
  • Zorg dat er geen raster processing template aan staat (zet die voor de laag dus op “None” zodat je met de hoogtewaardes kunt gaan rekenen)
  • Gebruik de Stack Profile tool
    ** Teken een nieuwe lijn of geef een bestaande laag met lijnen op.
    ** Geef 1 of meerdere AHN2/AHN3 rasters op als profile target(s) - zo kan je dus bijvoorbeeld eenvoudig verschillende opname momenten vergelijken of bijvoorbeeld het ruwe en maaiveld bestand
    ** De output is een tabel. Hierin staat de afstand vanaf het begin van de lijn en de Z-waarde. Deze tabel kan je gebruiken om gegevens in een grafiek weer te geven (zie voorbeeld hieronder), maar de tabel kan je natuurlijk dan ook weer exporteren naar alle formaten die ArcGIS Pro je toestaat te exporteren (waaronder ook csv).

Mocht je interesse hebben om het uit te proberen en geen beschikking hebben over ArcGIS kan je altijd een trial aanvragen: Hoe ga ik om met ECW-bestanden in ArcGIS for Desktop en ArcGIS for Server?.

Bison (https://www.food4rhino.com/app/bison) heeft een pluginn ontwikkeld voor Rhino/Grasshopper. Hiermee kunnen geotiff(float32) bestanden worden geïmporteerd in rhino/grasshopper. Op basis van deze input kunnen tal van oppervlakteanalyses worden uitgevoerd en gevisualiseerd. (https://www.bison.la/releases/Bisonv09Reference.pdf)
De geïmporteerde informatie kan vervolgens in vele formaten worden omgezet. (Rhino - Supported File Formats)

image

image

image