Hi, Thanks for testing and sending the results back for the community.
-Jukka Rahkonen- Lähettäjä: gdal-dev <gdal-dev-boun...@lists.osgeo.org> Puolesta Carl Godkin Lähetetty: maanantai 25. lokakuuta 2021 18.37 Vastaanottaja: gdal-dev@lists.osgeo.org Aihe: Re: [gdal-dev] gdallocationinfo + WCS question... Hi, I just built GDAL 3.3.3 RC1 here on Windows with VS2019 and PROJ 8.1.1 and it built easily. I retried the gdallocationinfo queries that I was having trouble with beforehand (with 3.3.1) and both the USGS National Map and my local WCS server queries now return the expected values. I know others had tested the fix Even mentioned but I wanted to try it myself. Thanks again for the explanation and the fix, carl On Thu, Oct 21, 2021 at 2:00 PM Carl Godkin <cgod...@gmail.com<mailto:cgod...@gmail.com>> wrote: Thanks, Even. That makes sense. (Finally!) I can wait a little bit to come back to this so I'll either build 3.3 HEAD on my own or await the next release. It's great that this isn't a random server load thing or something! The weird result on my own WCS server also seemed to be caused by a bad "xmin" value which is how I started down this road. Thanks again, carl On Thu, Oct 21, 2021 at 1:49 PM Even Rouault <even.roua...@spatialys.com<mailto:even.roua...@spatialys.com>> wrote: Le 21/10/2021 à 22:32, Rahkonen Jukka (MML) a écrit : Hi Even, You are right, re-running the request -oo clear_cache gives me the same error about 3 bands. Obviously our server was configured right when I used it with GDAL last time. Now gdallocationinfo does not work even for me because I wiped the working configuration from my cache. We will fix the service metadata when we go from beta to official soon. The nationalmap.gov<http://nationalmap.gov> service does not still work for me. I cleared the cache but I still get a 404 error from the server. Subsets in the reques are these SUBSET=x(-20037507.0672000013,-13590403.131594259)&SUBSET=y(5303978.1067894027,5304490.1067842878) The xmin value in the request you get is obviously wrong. The subset I get with GDAL master and head of 3.3 branch is: SUBSET=x%(-13591427.131584032,-13590403.131594259)&SUBSET=y(5303978.1067894027,5304490.1067842878) I suspect you might use a GDAL release version, that hasn't yet received the following recent fix that is likely related to that issue and will be in 3.3.3 and 3.4.0: https://github.com/OSGeo/gdal/commit/8af67b64371c275a3770e75da1bd0d4bb7329295 (I get the same error as you with the 3.2 branch, so I strongly suspect the above is indeed the fix for that) Even -- http://www.spatialys.com My software is free, but my time generally not.
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev