Hoogtestatistieken op specifieke locatie met Python

Je kan de processing tool ook al handmatig installeren in QGIS; zie de instructies hier in de README.

Dat klopt, daarbij moet je wel rekening houden met de volgende twee zaken (zie ook deze reactie):

  • de WCS service kan geen raster serveren van 1x1 pixel (bekende bug)
  • als je een WCS 2.0.1 request doet krijg je een multipart response terug waar een geotiff in zit, wil je direct een geotiff response gebruik dan WCS versie 1.0.0 (voorbeeld 1.0.0 GetCoverage request). Het capabilities document geeft aan welke WCS versie’s ondersteund worden:
    image

Dat gezegd hebbende, als je de hoogte van een gebouw wilt bepalen met de AHN3, zal je nog wel de terrein hoogte van de gevonden waarde af moeten trekken op die positie. Aangezien de AHN de absolute hoogte bevat t.ov. NAP. Zie ook de uitleg over het verschil tussen de DTM en DSM coverages in de AHN3 op de download pagina.

Persoonlijk zou ik de PDOK 3D hoogtestatistieken gebruiken om gebouw hoogte op te halen. Dan hoef je niet eerst te bepalen hoe hoog het maaiveld is rond het gebouw om de hoogte van het gebouw te bepalen.