Dan wordt de resolutie van de data verlaagd naar 2 x 2 meter. Hoe zorg ik ervoor dat ik toch 0.5 x 0.5 binnenhaal? Ik heb wat zitten zoeken naar/in de documentatie, maar dacht daaruit af te leiden dat het met width en height opgelost zou moeten zijn.
(ik heb ook een workbench voor het downloaden via kaartbladen, maar bij stukken die over meerdere kaartbladen vallen is de wcs vaak veel sneller)
Uit de documentatie: “Specifies the pixel size to be returned, in the X and Y dimensions. For gridCRS UTM values this is returned in meters. For gridCRS Lat/Long values this is returned in degrees/pixel. Minimum values for G-EGD are 0.5m and .0000045 degrees; these equate to full resolution data for WV01”.
Volgens mij komt dit omdat de width en height niet goed zijn berekend. De resolutie die je o.b.v. de bounding box en width/height krijgt, kun je als volgt berekenen:
x res: (max x - min x) / width => (223094.446 - 200246.970) / 11424 = 2.000 (afgerond)
y res: (max y - min y) / height => (507561.821 - 501229.934) / 3166 = 2.000 (afgerond)
Ofwel, je moet een width van (223094.446 - 200246.970) / 0,5 = 45695 en een height van (507561.821 - 501229.934) / 0,5 = 12664 gebruiken. Het is maar de vraag of dat in één keer gaat lukken, omdat dat een erg grote hoeveelheid data betreft.
Dank je, dat heb ik inderdaad fout gedaan. Ik heb het aangepast, maar het maakt geen verschil voor het resultaat. Ook niet met gridofsets of RESX etc… aan de link toegvoegd
Ik zie al iets wat foput gaat. Blijkbaar wordt het automatisch teruggezet naar 800x600. Ik krijg deze melding:
The retrieved value for ‘WIDTH’ is unnaturally large or 0, width and height will be set to 800x600. This can be overrided by specifying your own values
Het rare is dat ik natuurlijk al mijn eigen width en height heb aangegeven.
Ik los het nu op door het invoerbestand waar ik de hoogtegegvens van wil hebben in FME te tilen met max 600x800 per tile en dan de losse tiles down te loaden en te mosaicen.
aanvulling op wat @fsteggink zei. Ik heb hier ook lang mee zitten stoeien. Uiteindelijk werkte het wel goed door hele coördinaten te gebruiken i.p.v. werken met decimalen.
Ik krijg bij dit request juist een exception report waarin wordt gezegd dat de max width en height 4000 is:
<?xml version='1.0' encoding="UTF-8" ?>
msWCSGetCoverage(): WCS server error. Raster size out of range, width and height of resulting coverage must be no more than MAXSIZE=4000.
V.w.b. het downloaden van meerdere tiles zou ik dan ook deze waarde aanhouden. 1 tile van 4000x4000 is vergelijkbaar met ca. 33 tiles van 800x600. Je vraagt daarmee ineens gebieden van 2x2 km op.
En afronden zoals @Jellest is ook aan te raden, zodat je beter kunt garanderen dat de tiles naadloos op elkaar zullen aansluiten. Je moet dan afronden op de resolutie die je wilt opvragen, dus als je doelresolutie 0,5 meter is, kun je afronden op 0,5 meter.
Maar als je dan toch data nodig hebt van een gebied van 23 x 6 km, kun je beter de AHN-bladen downloaden en die door FME aan elkaar laten mosaicken. Voor kleinere gebieden zal de WCS sneller zijn, maar voor grote gebieden zal het downloaden en aan elkaar plakken van losse bladen waarschijnlijk sneller zijn.
