[gdal-dev] Retrieving 3DEP data via GDAL's WCS driver...

[Carl Godkin]

Hi GDAL list,

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).

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.

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?

Does anyone use this service regularly with GDAL tools (or just via ArcGIS
or some other tools)?

Thank you very much for your help,
carl


On Sat, Jun 13, 2020 at 5:27 PM Carl Godkin <cgodkin at gmail.com> wrote:

> Hi,
>
> This page https://viewer.nationalmap.gov/services/ includes a link under
> "WCS Services" that lists a number of subdatasets.  This command
>
> gdalnfo "
> https://elevation.nationalmap.gov/arcgis/services/3DEPElevation/ImageServer/WCSServer?request=GetCapabilities&service=WCS
> "
>
> works fine.  I have been retrieving elevation data from the first
> subdataset:
>
>   SUBDATASET_1_NAME=WCS:
> https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation
>   SUBDATASET_1_DESC=DEP3Elevation
>
> and most of the time it works but I think I must be doing something wrong
> since sometimes there are large missing areas.
>
> For instance this area (Boston) works just fine:
>
> gdal_translate "WCS:
> https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation"
> test.tif -projwin -71.07 42.36 -71.03 42.33 -projwin_srs EPSG:4326 -tr 100
> 100
>
> but this command with only a different "projwin" (in Colorado) is about
> 75% missing:
>
> gdal_translate "WCS:
> https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation"
> test.tif -projwin -107.03 37.28 -107 37.25 -projwin_srs EPSG:4326 -tr 100
> 100
>
> Can anyone help me understand what I'm doing wrong here and how to fix
> it?
>
> I am using a very recent OSGeo4W featuring GDAL 3.0.4, released
> 2020/01/28.
>
> Thanks very much for your help,
>
> carl
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20201101/1a317300/attachment.html>