Optimale manier om hoogte data op te halen?

halam 2020-05-12 20:34:41 UTC #21

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

https://community.rstudio.com/t/assistance-with-running-r-code-ahnextract/65686

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

BVD

halam 2020-05-12 20:38:38 UTC #22

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,

Jellest 2020-05-13 15:13:27 UTC #23

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.

halam 2020-05-13 19:52:18 UTC #24

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”

Jellest 2020-05-14 16:18:58 UTC #25

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?

halam 2020-05-14 19:08:37 UTC #26

> 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.

halam 2020-05-15 14:17:34 UTC #27

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

Jellest 2020-05-15 17:41:15 UTC #28

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).

halam 2020-05-15 18:39:06 UTC #29

wordt er een beetje moedeloos van

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…

halam 2020-05-18 21:14:54 UTC #30

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

Jellest 2020-05-19 11:08:30 UTC #31

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

halam 2020-05-23 21:41:48 UTC #32

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 !

https://github.com/Martien1973/rAHNextract

Mogelijk gemaakt door Discourse, best bekeken met JavaScript ingeschakeld