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

[Carl Godkin]

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/20200613/ded689ab/attachment.html>