Raster size out of range bij downloaden geotiff van AHN3

Beste Lezer,

Nu de TMS service van PDOK is uitgefaseerd probeer ik de raster afbeeldingen via de WCS service te downloaden.

Echter als ik deze service aanspreek via:

http://geodata.nationaalgeoregister.nl/ahn3/wcs?request=GetCoverage&service=WCS&version=2.0.1&coverage=ahn3_05m_dsm&crs=EPSG:4326&response_crs=EPSG:4326&bbox=4.75,52.5,4.76,52.49&width=10&height=10&FORMAT=image/tiff&coverageid=ahn3_05m_dsm

of via:

http://geodata.nationaalgeoregister.nl/ahn3/wcs?request=GetCoverage&service=WCS&version=2.0.1&coverage=ahn3_05m_dsm&crs=EPSG:28992&response_crs=EPSG:28992&bbox=120000,480010,120010,480000&width=20&height=20&FORMAT=image/tiff&coverageid=ahn3_05m_dsm

(verschil alleen in crs)

zie ik een “Raster size out of range” error.

Ik baseer mijn query op de link zoals hier gegeven: https://www.pdok.nl/wcs

Enig idee wat ik hier fout doe? De width en height zijn erg klein en de bounding box ook.

Mijn uiteindelijke doel is om GeoTiffs te maken voor verdere verwerking met Rasterio. Mocht ik dan beter 1 van de andere services gebruiken, dan hoor ik dat natuurlijk ook graag!

Groet
Jorrit

Welke kaart viewer?
De pdok viewer hebben we niet uitgefaseerd…

Hi Wouter, excuus, ik bedoelde de TMS service.

Gelukkig ik schrok al :smiley:

Maar i.p.v. de TMS is er (nog steeds) de WMTS zoals we gebruiken in de pdok viewer.

M.b.t. je voorbeeld requesten lijkt wel iets anders (mogelijk) scheef te zitten, ik zal d’r ff naar kijken…

@jorrit het probleem is dat je WCS aanroept met parameters die bij de WMS horen. Een WMS request komt in het kort neer op “geef me een plaatje van dit gebied met deze afmetingen”. Bij een WCS liggen die afmetingen echter vast, als gevolg van de pixelgrootte in het coverage. Zie bijvoorbeeld dit DescribeCoverage request van ahn3_05m_dsm in de AHN3 WCS. De afmetingen van een pixel in dit coverage zijn een 0.5m x 0.5m:

<gml:offsetVector srsName="http://www.opengis.net/def/crs/EPSG/0/28992">0.500000 0</gml:offsetVector>
<gml:offsetVector srsName="http://www.opengis.net/def/crs/EPSG/0/28992">0 -0.500000</gml:offsetVector>

Dus bij een WCS request van “geef met dit raster bestand van dit gebied” krijg je altijd een raster bestand met dezelfde afmetingen (zolang het gebied en het het bevraagde coverage hetzelfde blijven).

Dit gezegd hebbende is het volgens mij ook mogelijk om via de WCS te resamplen, en daarmee dus een raster op te vragen met een andere pixelgrootte dan de bron coverage in de WCS. Ik heb hier helaase geen voorbeeld van paraat.

Hier is wel een werkend voorbeeld van jouw WCS request:

https://geodata.nationaalgeoregister.nl/ahn3/wcs?service=WCS&Request=GetCoverage&version=2.0.1&CoverageId=ahn3_05m_dsm&format=image/tiff&subset=x(120000,120010)&subset=y(480000,480010)

Dit request resulteert overigens in een multipart response met:

  • out.tif (het geotiff bestand)
  • out.tif.aux.xml (metadata bestand die NODATA specificeert)

Om het bestand te kunnen openen zul je eerst de multipart response moeten opknippen. Ik heb een bash script wat je hiervoor kan gebruiken.

3 likes

Nog een kleine toevoeging; resamplen kan onder andere dmv de ScaleFactor parameter:

https://geodata.nationaalgeoregister.nl/ahn3/wcs?service=WCS&Request=GetCoverage&version=2.0.1&CoverageId=ahn3_05m_dsm&format=image/tiff&subset=x(120000,120010)&subset=y(480000,480010)&ScaleFactor=0.5

Dit request levert dus een GeoTiff bestand op met een lagere resolutie (grotere pixels).

Zie ook de MapServer WCS documentatie voor een overzicht van de parameters van een WCS request.

Hi Anton,

dank voor je hulp!

De url die ik gebruik komt uit dit voorbeeld op de PDOK site en ik gebruik dezelfde parameters, alleen probeer ik ah3 te bevragen.

Ik heb het nu werkend met de volgende parameters:

http://geodata.nationaalgeoregister.nl/ahn3/wcs?request=GetCoverage&service=WCS&version=1.0.0&coverage=ahn3_05m_dsm&crs=EPSG:28992&response_crs=EPSG:28992&bbox=120000,487500,121000,486500&width=2000&height=2000&FORMAT=image/tiff

Dit levet geotiffs op met een dimensie van 2000x2000 pixels bij een coverage van 1000 bij 1000 meter dus een resolutie van 0.5m (welke ik nodig heb).

Het opvragen van een grotere tiff formaat geeft de eerder genoemde foutmelding.

en ik kwam er via via achter dat ik hier nog wel gewoon de kaartbladen kan downloaden:

https://www.arcgis.com/apps/Embed/index.html?appid=a3dfa5a818174aa787392e461c80f781

Wellicht handig voor een volgende google-aar