[gdal-dev] gdallocationinfo + WCS question...

Mon Oct 25 09:43:14 PDT 2021

Hi,

Thanks for testing and sending the results back for the community.

-Jukka Rahkonen-

Lähettäjä: gdal-dev <gdal-dev-bounces at lists.osgeo.org> Puolesta Carl Godkin
Lähetetty: maanantai 25. lokakuuta 2021 18.37
Vastaanottaja: gdal-dev at 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 <cgodkin at gmail.com<mailto:cgodkin at 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.rouault at spatialys.com<mailto:even.rouault at 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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20211025/ad5918d8/attachment.html>