<div dir="ltr">Hi again,<div><br></div><div>Like Jukka, I get a 404 from the National Map server. (I kept removing the cache directory just in case that mattered.)</div><div><br></div><div>I just asked a colleague in Texas to try the National Map WCS command and it's not working for him either so Ari and Even are apparently lucky!</div><div><br></div><div>Thanks,</div><div>carl</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Oct 21, 2021 at 1:32 PM Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi">jukka.rahkonen@maanmittauslaitos.fi</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="FI" style="overflow-wrap: break-word;">
<div class="gmail-m_-5380455332683267470WordSection1">
<p class="MsoNormal"><span>Hi Even,<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">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.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">The <a href="http://nationalmap.gov" target="_blank">nationalmap.gov</a> 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<br>
SUBSET=x(-20037507.0672000013,-13590403.131594259)&SUBSET=y(5303978.1067894027,5304490.1067842878)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">and gdallocationinfo error is this:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Report:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> Location: (6446955P,13504941L)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> Band 1:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 1: HTTP error code : 404<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 1: Failed to get coverage data:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 1: C:\Users\JRAHKONEN\.gdal\wcs_cache\QpYSX.xml, band 1: IReadBlock failed at X offset 6295, Y offset 26376: Failed to get coverage data:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka-<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<div>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> Even Rouault <<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>> <br>
<b>Lähetetty:</b> torstai 21. lokakuuta 2021 22.50<br>
<b>Vastaanottaja:</b> Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>>; Carl Godkin <<a href="mailto:cgodkin@gmail.com" target="_blank">cgodkin@gmail.com</a>>; <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<b>Aihe:</b> Re: [gdal-dev] gdallocationinfo + WCS question...<u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<p>Jukka,<u></u><u></u></p>
<p>I don't understand how you get a non-error on the below request on your WCS service:<u></u><u></u></p>
<p><a href="https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?SERVICE=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=korkeusmalli__korkeusmalli&FORMAT=text/xml" target="_blank">https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?SERVICE=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=korkeusmalli__korkeusmalli&FORMAT=text/xml</a>
advertizes 3 bands<u></u><u></u></p>
<p>but the GeoTIFF returned by<u></u><u></u></p>
<p><a href="https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?SERVICE=WCS&REQUEST=GetCoverage&VERSION=2.0.1&COVERAGEID=korkeusmalli__korkeusmalli&SUBSET=E%28498656,500704%29&SUBSET=N%287542384,7543408%29&Format=image/tiff" target="_blank">https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?SERVICE=WCS&REQUEST=GetCoverage&VERSION=2.0.1&COVERAGEID=korkeusmalli__korkeusmalli&SUBSET=E%28498656,500704%29&SUBSET=N%287542384,7543408%29&Format=image/tiff</a><u></u><u></u></p>
<p>has just one band.<u></u><u></u></p>
<p>Or maybe you have in your ~/gdal/wcs_cache a previous version of the DescribeCoverage response that advertizes just one band ?<u></u><u></u></p>
<p><u></u> <u></u></p>
<p>And for what is worth,<u></u><u></u></p>
<p><span lang="EN-US" style="font-family:"Courier New",serif">gdallocationinfo "WCS:</span><span style="font-family:"Courier New",serif"><a href="https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation" target="_blank"><span lang="EN-US">https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation</span></a></span><span lang="EN-US" style="font-family:"Courier New",serif">"
-wgs84 -122.086 42.948 --debug ON</span><span lang="EN-US"><u></u><u></u></span></p>
<p><span style="font-family:"Courier New",serif">returns 1882.00549316406 for me</span><u></u><u></u></p>
<p>The multipart answer doesn't seem to be an issue.<u></u><u></u></p>
<p><u></u> <u></u></p>
<div>
<p class="MsoNormal">Le 21/10/2021 à 20:16, Rahkonen Jukka (MML) a écrit :<u></u><u></u></p>
</div>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<p class="MsoNormal"><span>Hi,</span><u></u><u></u></p>
<p class="MsoNormal"><span> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US">At least it seems to work with our WCS service<br>
<br>
gdallocationinfo WCS:<a href="https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?version=2.0.1&coverage=korkeusmalli__korkeusmalli" target="_blank">https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?version=2.0.1&coverage=korkeusmalli__korkeusmalli</a>
-wgs84 27 68</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US">Report:</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> Location: (228000P,119567L)</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> Band 1:</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> Value: 303.463989257812</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka Rahkonen-</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> gdal-dev <a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">
<gdal-dev-bounces@lists.osgeo.org></a> <b>Puolesta </b>Carl Godkin<br>
<b>Lähetetty:</b> torstai 21. lokakuuta 2021 18.47<br>
<b>Vastaanottaja:</b> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<b>Aihe:</b> [gdal-dev] gdallocationinfo + WCS question...<u></u><u></u></p>
</div>
<p class="MsoNormal"> <u></u><u></u></p>
<div>
<p class="MsoNormal">Hi,<u></u><u></u></p>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">I am trying to use gdallocationinfo to query point locations from a WCS server and not having any luck. I don't see anything in the docs that says this is unsupported but I can't make it work.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">First of all, is there a better way to get an elevation from a given lat,lon using WCS?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Here's an example query using the USGS National Map server:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New",serif">gdallocationinfo "WCS:<a href="https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation" target="_blank">https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation</a>"
-wgs84 -122.086 42.948 --debug ON</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">This is a point in the middle of Crater Lake which I picked since it's flat and I know what to expect.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Anyway, looking at the debug output it appears that the program is querying for a GeoTIFF of a VERY large area and that's failing.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">I actually want to do this with my own WCS server but thought it might be mis-configured somehow so I tried the USGS server. With my own server, gdallocationinfo seems to be requesting a strange shape from the server (as I can see through
the debug output).<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks for any help,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">carl<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"><br>
<br>
<u></u><u></u></p>
<pre>_______________________________________________<u></u><u></u></pre>
<pre>gdal-dev mailing list<u></u><u></u></pre>
<pre><a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><u></u><u></u></pre>
<pre><a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><u></u><u></u></pre>
</blockquote>
<pre>-- <u></u><u></u></pre>
<pre><a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><u></u><u></u></pre>
<pre>My software is free, but my time generally not.<u></u><u></u></pre>
</div>
</div>
</blockquote></div>