<div dir="ltr"><div>Hi GDAL list,</div><div><br></div><div>I sent the following question in June about using gdal_translate with the USGS' WCS Services.  One kind person tried to help but their suggestion did not lead anywhere.  Since then I've also written a few times to the USGS help addresses and all they can tell me is that the areas look okay to them (through their ESRI interfaces, I think).  </div><div><br></div><div>I just tried it again with (the latest OSGeo4W) GDAL 3.1.4 and found the same problem.  Since the USGS thinks the data sets are fine, it's either a usage mistake on my part (still the leading candidate!) or possibly a bug in the WCS driver.</div><div><br></div><div>Can anyone suggest how I might investigate this since I'm not very familiar with the inner workings of gdal_translate and the WCS driver?   </div><div><br></div><div>Does anyone use this service regularly with GDAL tools (or just via ArcGIS or some other tools)?</div><div><br></div><div>Thank you very much for your help,</div><div>carl</div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 13, 2020 at 5:27 PM Carl Godkin <<a href="mailto:cgodkin@gmail.com">cgodkin@gmail.com</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 dir="ltr">Hi,<div><br></div><div>This page <a href="https://viewer.nationalmap.gov/services/" target="_blank">https://viewer.nationalmap.gov/services/</a> includes a link under "WCS Services" that lists a number of subdatasets.  This command </div><div><br></div><div><font face="monospace">gdalnfo "<a href="https://elevation.nationalmap.gov/arcgis/services/3DEPElevation/ImageServer/WCSServer?request=GetCapabilities&service=WCS" target="_blank">https://elevation.nationalmap.gov/arcgis/services/3DEPElevation/ImageServer/WCSServer?request=GetCapabilities&service=WCS</a>"</font></div><div><br></div><div>works fine.  I have been retrieving elevation data from the first subdataset:</div><div><br></div><div><font face="monospace">  SUBDATASET_1_NAME=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><br>  SUBDATASET_1_DESC=DEP3Elevation</font><br></div><div><br></div><div>and most of the time it works but I think I must be doing something wrong since sometimes there are large missing areas.</div><div><br></div><div>For instance this area (Boston) works just fine:</div><div><br></div><div><font face="monospace">gdal_translate "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>" test.tif -projwin -71.07 42.36 -71.03 42.33 -projwin_srs EPSG:4326 -tr 100 100</font><br></div><div><font face="monospace"><br></font></div><div><font face="arial, sans-serif">but this command with only a different "projwin" (in Colorado) is about 75% missing:</font></div><div><br></div><div><font face="monospace">gdal_translate "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>" test.tif -projwin -107.03 37.28 -107 37.25 -projwin_srs EPSG:4326 -tr 100 100</font><br></div><div><font face="monospace"><br></font></div><div><font face="arial, sans-serif">Can anyone help me understand what I'm doing wrong here and how to fix it?  </font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I am using a very recent OSGeo4W featuring </font>GDAL 3.0.4, released 2020/01/28.</div><div><br></div><div>Thanks very much for your help,</div><div><br></div><div>carl</div></div>
</blockquote></div></div>