<div dir="ltr">Hi Tanya,<div><br></div><div>Thanks for the suggestion.  I wondered about that too but there's only one endpoint for "3DEP" in that WCS link so I don't know how to switch<br>products with WCS.  It appears from this view</div><div><br></div><div><a href="https://viewer.nationalmap.gov/basic/?basemap=b1&category=ned&q=&zoom=9&bbox=-107.61108398,36.81368259,-105.38909912,38.17775097&preview=&avail=DEM">https://viewer.nationalmap.gov/basic/?basemap=b1&category=ned&q=&zoom=9&bbox=-107.61108398,36.81368259,-105.38909912,38.17775097&preview=&avail=DEM</a> 1/9 arc-second&refpoly=<br></div><div><br></div><div>that I might be near the edge of the 1/9 arc second availability but I would have thought that the service would fill in from the next higher resolution but I can find little in the way of documentation of this assumption.  I have also written to the The National Map help email and I'm hoping that someone either here or there can set me straight.</div><div><br></div><div>Thank you so much for the response,</div><div>carl</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 13, 2020 at 6:20 PM TC Haddad <<a href="mailto:tchaddad@gmail.com" target="_blank">tchaddad@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"><div><br></div><div>Hi Carl,</div><div><br></div><div>Is it possible that you are just encountering the edges of the subdataset? <br></div><div><br></div><div>That URL you provided in your first email had the list of indexes for the collection of data sets, maybe you can cross reference the shapes it was showing with the edges you are getting, to see if you need to switch to a different product when you are in certain areas:<br></div><div><br></div><div><a href="http://www.arcgis.com/home/webmap/viewer.html?url=https%3A%2F%2Findex.nationalmap.gov%2Farcgis%2Frest%2Fservices%2F3DEPElevationIndex%2FMapServer&source=sd" target="_blank">http://www.arcgis.com/home/webmap/viewer.html?url=https%3A%2F%2Findex.nationalmap.gov%2Farcgis%2Frest%2Fservices%2F3DEPElevationIndex%2FMapServer&source=sd</a></div><div><br></div><div>Tanya<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 13, 2020 at 5:28 PM Carl Godkin <<a href="mailto:cgodkin@gmail.com" target="_blank">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>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div>
</blockquote></div>