Optimale manier om hoogte data op te halen?

Beste Halam,
Ik las jouw bericht aan EkkoNap: ‘ik heb een andere manier om een dtm te maken met ahn voor in een autodesk produkt. Als je wilt kan ik je het recept wel doormailen’
Zou je het mij ook kunnen mailen?
Groet,
Martien.

Beste R Duivenvoorde,
Ik las jouw bericht in post ‘Optimale manier om hoogte data op te halen?’ op Geoforum.
Ik zou heel graag via WCS het Z coordinaat van een reeks XY coordinaten uit het AHN3 ophalen en indien mogelijk in .csv formaat downloaden.
Ben al aan het orienteren hoe Postman werkt of dat dit via QGIS kan, maar het is allemaal nieuw voor mij, maar kan wel wat hulp gebruiken.
Heb je misschien tips, of kun jij misschien dit misschien heel makkelijk voor een reeks XY doen? Als ik één keer de code heb, zou ik wellicht heel gemakkelijk die code voor mijn volgende reeksen XY coordinaten kunnen gebruiken…
Bij voorbaat dank,
Martien.

Ik heb het link met het documentje hieronder nog eens geplakt. Een trukje om punten te filteren is om met notepad++ een scriptje te maken die 99 regels delete uit de cvs en er een overslaat. Zo hou je 1 procent over, een behapbaar cvs bestand dat toch een redelijk goed beeld geeft van het oppervlakte. Omgezet naar DWG kun je ze overal wel in gebruiken.

AHN in Revit verwerken.pdf (1,0 MB)

Beste Jelle

Deze scripts zijn volgens mij erg handig!
Ik heb enkele pogingen gedaan om je scripts werkend te krijgen maar daar heb ik nog wat moeite mee.
Daarom heb ik op een ander forum wat hulp gevraagd om met R studio dit draaiend te krijgen

Zou jij hier naar willen kijken om potentiële gebruikers klein beetje op weg kunnen helpen?

BVD

Dit 'demo’tje laat zien hoe je xyz direct goed uit Qgis kunt exporteren.
XYZ kunt je in de meeste DWG programma’s of in Revit importeren als topo surface,

Hoi

Je bent me net voor want ik was bezig om een Nederlandstalige post op deze forum hierover te plaatsen. Deze kun je nu hier vinden.

in deze post staat ook uitgelegd hoe de package preceis werkt en hoe je hem moet installeren. Als je een foutmelding krijgt dat hij de functie ahn_point() niet kan vinden, dan heb je waarschijnlijk de package nog niet gedownload EN geladen.

Super!
Ik denk dat het handig is de discussie hierover hier voort te zetten

Installatie

library(raster)
library(devtools)
devtools::install_github(“Jellest/rAHNextract”)
library(rAHNextract)

Dit resultaat


package ‘R.methodsS3’ successfully unpacked and MD5 sums checked
package ‘shiny’ successfully unpacked and MD5 sums checked
package ‘manipulateWidget’ successfully unpacked and MD5 sums checked
package ‘miniUI’ successfully unpacked and MD5 sums checked
package ‘base64enc’ successfully unpacked and MD5 sums checked
package ‘webshot’ successfully unpacked and MD5 sums checked
package ‘httpuv’ successfully unpacked and MD5 sums checked
package ‘xtable’ successfully unpacked and MD5 sums checked
package ‘sourcetools’ successfully unpacked and MD5 sums checked
package ‘fastmap’ successfully unpacked and MD5 sums checked
package ‘glue’ successfully unpacked and MD5 sums checked
Error: Failed to install ‘rAHNextract’ from GitHub:
** (converted from warning) cannot remove prior installation of package ‘glue’**

Dus de install verloopt niet vlekkeloos
Functie wordt niet herkend

ahn_point(name = “Utrecht point”, X = 136550, Y = 456060)
Error in ahn_point(name = “Utrecht point”, X = 136550, Y = 456060) :
could not find function "ahn_point"

Dit probleem kom ik ook wel eens tegen met het installeren van andere packages. Snap dat zelf ook nog niet zo goed. Bij mij lukte het toen om de package waar het mis ging eerst te verwijderen en deze handmatig te installeren.

In dit geval:

remove.packages(“glue”)
install.packages(“glue”)
en dan weer
devtools::install_github(“Jellest/rAHNextract”)

kun je dat eens proberen?

> install.packages("glue")
WARNING: Rtools is required to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding:

https://cran.rstudio.com/bin/windows/Rtools/
Installing package into ‘C:/Users/Admin/Documents/R/win-library/4.0’
(as ‘lib’ is unspecified)

  There is a binary version available but the source version is
  later:
     binary source needs_compilation
glue  1.4.0  1.4.1              TRUE

  Binaries will be installed
trying URL 'https://cran.rstudio.com/bin/windows/contrib/4.0/glue_1.4.0.zip'
Content type 'application/zip' length 158875 bytes (155 KB)
downloaded 155 KB

package ‘glue’ successfully unpacked and MD5 sums checked

The downloaded binary packages are in
	C:\Users\Admin\AppData\Local\Temp\RtmpshAVvb\downloaded_packages



> devtools::install_github(“Jellest/rAHNextract”)
Error: unexpected input in "devtools::install_github(“"

Het is niet gelukt.
Heb andere varianten van de tekst geprobeerd maar het werkt niet.
Is Glue nu geinstalleerd of niet. Nogal onduidelijk om alleen te vertellen dat het gedownload is…

Vervolgens lijkt Rstudio me te zeggen dat het Jelle/rAHNextract misschien niet kan bereiken. Lastige materie.

Beste Jelle, ik ben een klein stapje verder. Blijkbaar mag je in Rstudio geef copy paste in de command line uitvoeren. ik heb de tekst eens uitgetikt en dat werkt de install wel. Heel maf.

Nu de volgende …

> ahn_point(name = “Utrecht point”, X = 136550, Y = 456060)
Error: unexpected input in “ahn_point(name = “”

Hij snapt de naamgeving niet. Is dit een verwijzing naar een ID in een database om weg te schrijven…
BVD Hans

Dat is bijzonder. Goed om te lezen dat het nu wel werkt.

Ik denk dat het nu ook fout gaat met het kopiëren en plakken. Probeer het eens zelf in te voeren… de " wordt namelijk verkeerd gekopieerd. Het is overigen niet verplicht om de ‘name’ parameter te gebruiken. Dit is slechts een naam die hij gebruikt als je bestanden wilt wegschrijven.

Ik ben overigens nog een klein foutje tegen gekomen als je de ‘resolution’ parameter niet zelf invult. Ik heb het aangepast. Je kunt het beste de package opnieuw installeren via devtools:install_github("Jellest/rAHNextract")

Als je nog problemen of gekke dingen tegenkomt, dan hoor ik het graag van je. Deze package in namelijk nog in ontwikkeling (pre-release).

1 like

wordt er een beetje moedeloos van :frowning:

devtools::install_github(“Jellest/rAHNextract”)
Error in loadNamespace(j ← i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) :
there is no package called ‘glue’

Glue weer opnieuw installeren denk ik

remove.packages(“glue”)
Removing package from ‘C:/Users/Admin/Documents/R/win-library/4.0’
(as ‘lib’ is unspecified)
Error in remove.packages : there is no package called ‘glue’

package ‘glue’ successfully unpacked and MD5 sums checked
Error: Failed to install ‘rAHNextract’ from GitHub:
(converted from warning) cannot remove prior installation of package ‘glue’

ja ja
Ik denk dat we beter de discussie kunnen voeren op het Rstudio forum.
Even volhouden… :stuck_out_tongue:

Beste Jelle!
Eureka. Rstudio heb ik nu zovaak getweaked en gespeeld dat het volgens mij AHN bespeelt
Wat ik me alleen afvraag is hoe deze tot waarde +58 boven NAP komt. +4.60 lijkt me eerder een reeele waarde. Kun ik de verwijzing laten lopen naar data die wél gefilterd is?

ahn_point(X = 136550, Y = 456073)
AHN3 0.5 resolution DSM (ruw) selected.
trying URL ‘https://geodata.nationaalgeoregister.nl/ahn3/wcs?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&COVERAGE=ahn3_05m_DSM&BBOX=136549.5,456072.5,136551,456074&CRS=EPSG:28992&RESPONSE_CRS=EPSG:28992&FORMAT=GEOTIFF_FLOAT32&WIDTH=3&HEIGHT=3
Content type ‘image/tiff’ length 562 bytes
downloaded 562 bytes

[1] “Download raster image succeeded.”
[1] “Intersecting raster. Getting elevation…”
[1] "Elevation of AHNelevation: 58.98 m."
[1] 58.98

ahn_point(X = 136550, Y = 456063)
AHN3 0.5 resolution DSM (ruw) selected.
trying URL ‘https://geodata.nationaalgeoregister.nl/ahn3/wcs?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&COVERAGE=ahn3_05m_DSM&BBOX=136549.5,456062.5,136551,456064&CRS=EPSG:28992&RESPONSE_CRS=EPSG:28992&FORMAT=GEOTIFF_FLOAT32&WIDTH=3&HEIGHT=3
Content type ‘image/tiff’ length 562 bytes
downloaded 562 bytes

[1] “Download raster image succeeded.”
[1] “Intersecting raster. Getting elevation…”
[1] “Elevation of AHNelevation: 4.61 m.”
[1] 4.61

De 58.98 m is de hoogte van het gebouw dat op Neude in Utrecht te vinden is. Zoals je ziet heeft hij de DSM geselecteerd, en neemt hij dus de hoogte van gebouwen mee. Als je het maaiveld hoogte wilt hebben, moet je `dem=“DTM” gebruiken als parameter. Bij deze specifieke coördinaten, heeft de DTM hier geen hoogte opgenomen.

ahn_point(X = 136550, Y = 456073, dem = "DTM")
AHN3 0.5 resolution DTM (maaiveld) selected.
trying URL 'https://geodata.nationaalgeoregister.nl/ahn3/wcs?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&COVERAGE=ahn3_05m_DTM&BBOX=136549.5,456072.5,136551,456074&CRS=EPSG:28992&RESPONSE_CRS=EPSG:28992&FORMAT=GEOTIFF_FLOAT32&WIDTH=3&HEIGHT=3'
Content type 'image/tiff' length 562 bytes
downloaded 562 bytes

[1] "Download raster image succeeded."
[1] "Intersecting raster. Getting elevation..."
[1] "Elevation of AHNelevation: NA m."
[1] NA
1 like

Jelle, ik begrijp dat er een manier is om een csv bestandje voorzien van X,Y coordinaen via Rstudio en rAHNextract te laten ‘vullen’ met Z waarden van het maailveld. De documentatie is erg vrij cryptisch en mager en ik kom er niet achter hoe je dit werkend krijgt. Kun jij daar iets over toelichten? Het is namelijk wel ontzettend intressant voor legio toepassingen !

Jelle, ik kan me voorstellen dat je het erg druk hebt en dit niet direct top priority is.
Zou fjin zijn als je een beetje uit de doeken kunt doen hoe bijvoorbeeld de csv file eruit ziet die in 262 wordt aangehaald.

C:/Users/jelle/Documents/coordinaten_Zwaanshoek.csv"

Nog mooier zou zijn als er een soort ‘specify file’ > browse functie zou worden gebruikt in Rstudio.
Zijn er andere mensen bekend met Rstudio die dit werkend kunnen krijgen?Zou leuk zijn :slight_smile:

Schermopname van deze WCS method die ik graag met wat hulp meer generiek wil krijgen voor ‘POINT SERIES’

=> INPUT XY RANGE
=> OUTPUT XYZ RANGE

Hoi ,

Excuses dat ik nu pas reageer. Heb dit even gemist in combinatie met een paar dagen er tussenuit. Nogmaals mijn excuses hiervoor.

R packages zijn vaak zo gemaakt dat ze werken met andere packages en dat ze om algemene R kennis vragen. Het inlezen van een .csv bestand en een functie in een loop draaien zijn voorbeelden daarvan. Mijn documentatie legt slechts uit hoe je de functies van mijn package moet gebruiken.

het ahn_hoogte.R bestand dat je hebt gevonden op de github pagina is inderdaad een voorbeeld hoe het moet. Dit was echter een eerste versie van het script dat iemand anders toen heeft gebruikt. Let wel, gebruik niet de rest van de code want dat is verouderd. Ik raad om het via de package te doen. Wel kun je de techniek gebruiken om een csv in te lezen (regel 262 in het bestand ahn_hoogte.R).

Wat bedoel je precies met ‘specify file’ > browse functie’? In R zul je alles moeten doen met inlezen van bestanden vai allerlei commands/functies.

2 likes

Hey Jelle, Ik probeer momenteel voor een project hoogtes op te halen van punten binnen amsterdam. Mijn vraag is eigenlijk: hoe weet jij de correspodentie tussen punten op het raster van de TIFF en de long/lats?

Alvast bedankt,

Pim

Hoi Pim,

De TIFF bestanden van de AHN zijn GeoTIFF. Deze tiff bestanden bevat een attribuut met de extent (xmin,ymin,xmax,ymax) waar het is. Als je deze tiff bestanden als raster in je script/programma laadt, dan weet hij meteen waar het is via de extent attribuut. Beantwoord dit een beetje je vraag?