Hoe WCS van AHN 6 aan te roepen in CRS 28992

[Anton]

De WCS van AHN 6 werkt anders dan de WCS van AHN 4 via PDOK. Ik wil deze AHN 6 WCS aanroepen vanuit een eigen applicatie zoals ik ook AHN 4 kan ophalen. Maar ik stuit hierbij op issues die ik niet eenvoudig kan oplossen, dus hopelijk is er iemand die me op weg kan helpen.

Als ik via PDOK een stukje ophaal uit AHN 4 dan gebruik ik deze url:

https://service.pdok.nl/rws/ahn/wcs/v1_0?service=wcs&version=2.0.1&request=GetCoverage&CoverageId=dtm_05m&Format=image/tiff&SubSet=X(246779.06,248603.06)&SubSet=Y(561427.98,562695.98)&ScaleSize=X(1824)&ScaleSize=Y(1268)

Vraag ik exact hetzelfde stukje op van AHN 6 dan zou ik deze url moeten aanroepen:

https://api.ellipsis-drive.com/v3/ogc/wcs/22f5740b-1276-42e0-9d2d-dd2a9745e928?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&FORMAT=GeoTIFF&COVERAGE=c936131a-2c4f-44ce-b321-a085c6c0ce3b&BBOX=246779.06,248603.06,561427.98,562695.98&CRS=EPSG:28992&RESPONSE_CRS=EPSG:28992&WIDTH=1824&HEIGHT=1268

Via de GetCapabilities kwam ik er niet uit, niet alle parameters lijken beschreven, dus via QGIS gekeken hoe de aanvraag tot stand komt.

QGIS gebruikt EPSG:3857 WGS84/Pseudo-Mercator, ondanks dat het project EPSG:28992 hanteert. Als ik deze CRS instel bij de source, dan wordt er niets opgehaald uit de WCS. Laat ik de CRS op 2867, dan worden Mercator coördinaten gehanteerd in de Bounding Box en de CRS en Response_CRS verwijzen dan naar 3857.

Als ik deze gegenereerde url aanroep waarbij ik 3857 aan heb gepast naar 28992 en RD coördinaten instel, dan krijg ik wel een GeoTiff van pakweg 160KB maar die ziet er zo uit:

Haal ik hetzelfde gebied uit AHN 4, dan krijg ik een GeoTiff van 6.5MB die er zo uitziet:

Het is niet dat het gebied in AHN 6 nog niet beschikbaar is, daar heb ik op gelet.

Volgens de Capabilities van AHN 6 wordt 28992 ondersteund:

Maar wat doe ik hier fout? Iemand een idee?

[sbjager]

Ik heb al eens eerder gesignaleerd dat de WCS-services vanuit ellipsis drive niet lekker werken in QGis: Ellipsis WCS-services werken niet lekker in QGis
Helaas is daar nooit een reactie op gekomen, maar het lijkt er op dat WCS nog steeds niet echt goed ondersteund word helaas.

Edit: Ik heb dus geen idee wat er fout gaat

[Anton]

Haha Toch bedankt, dacht al dat ik onnozel was, dat lijkt dan weer mee te vallen

[Anton]

Iemand die hier nog iets over kan zeggen? @AHN misschien?

[AHN]

In de BBOX lijken de coördinaten omgedraaid en de URL was van het oppervalktemodel (DSM), niet het maaiveldmodel (DTM). Onderstaande URL geeft wel het gevraagde beeld:
https://api.ellipsis-drive.com/v3/ogc/wcs/231a4c66-9de9-4f41-91af-e4e991df9160?requestedEpsg=28992&SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&FORMAT=GeoTIFF&COVERAGE=05df775a-ae39-4bb0-a52e-d3044c993d1e&BBOX=246779.06,561427.98,248603.06,562695.98&CRS=EPSG:28992&RESPONSE_CRS=EPSG:28992&WIDTH=1824&HEIGHT=1268

[Anton]

Ah, dat was inderdaad een foutje, waarschijnlijk bij het testen de volgorde van de SubSet van AHN 4 aangehouden.

Hoewel het gebied nu wel klopt, kan ik helaas niet veel met het resultaat. Het lijkt wel of de GeoTiff niet correct wordt gegenereerd, of in een ander tiff formaat dan vroeger wordt aangeleverd. Ik gebruik de AHN om Surfaces te maken in Civil 3D. Dit is het resultaat van de AHN 6 WCS:

Terwijl dit het resultaat van AHN 4 WCS is:

In QGIS kan ik de gegenereerde GeoTiff wel inlezen. Het lijkt erop dat QGIS wat meer formaten ondersteunt dan Civil 3D. Wel blijft het raar (wat Stefan ook al schreef in zijn topic) dat de naamgeving en de dubbele band beter kan:

Bovenstaand is No Description (tenzij gedownload, dan is het een Guid). Onderste is AHN 4 en geeft keurig de bron en de hoogten.

Verzoek aan AHN. Is het mogelijk om:

• Het GeoTiff formaat gelijk te houden aan hoe de AHN 4 wordt geleverd? Er zijn meer applicaties dan alleen QGIS of ArcGis die iets met hoogten doen.
• Maar 1 band te leveren in plaats van 2?
• Omschrijvingen van de bron en de band in duidelijke taal?

Dus eigenlijk om de levering van AHN via de dataroom op dezelfde wijze in te richten als AHN via PDOK.