Hoi @Diede,
Er zit een bug in MapServer (wat PDOK gebruikt als server backend voor de WCS) die ervoor zorgt dat het niet mogelijk is een GetCoverage request te doen met de dekking van 1 pixel. Zie het draadje op de Mapserver mailing list ([mapserver-users] WCS GetCoverage request one by one pixel).
Een workaround kan zijn om een coverage van 2x2 pixels op te vragen en vervolgens de waarde uit te lezen van de pixel die overlapt met je POI.
Dit request werkt namelijk wel:
Merk overigens op dat de WCS een multipart response terug geeft. Ik heb hier een bash script voor het splitten van multipart responses (Bash script for splitting multipart response files in seperate files #wcs #bash #multipart · GitHub).
Zie het onderstaande bash script voor het uitlezen van de hoogte van een punt locatie (in RD coordinaten).
#!/usr/bin/env bash
x=196396.84
y=450032.55
x_lower_bound=$(python3 -c "print(($x//0.25)*0.25)")
# retrieving coverage 2x2, pixel size is 0.5
x_upper_bound=$(python3 -c "print($x_lower_bound + 1.0)")
y_lower_bound=$(python3 -c "print(($y//0.25)*0.25)")
y_upper_bound=$(python3 -c "print($y_lower_bound + 1.0)")
curl "https://geodata.nationaalgeoregister.nl/ahn3/wcs?service=WCS&Request=GetCoverage&version=2.0.1&CoverageId=ahn3_05m_dsm&format=image/tiff&subset=x($x_lower_bound,$x_upper_bound)&subset=y($y_lower_bound,$y_upper_bound)" -o wcs_response
split-mulitpart wcs_response wcs_response_split
gdallocationinfo wcs_response_split/coverage/out.tif -geoloc $x $y
Zie overigens ook deze vraag hier: point - How to query elevation data from "Digitales Geländemodell (DGM) NRW", maybe using WCS? - Geographic Information Systems Stack Exchange
Vermoedelijk gebruikt die WCS service ook MapServer als service backend.