Hi again, In case anyone else is interested, I followed one of the links that Michael provided and found that I should be able to use WCS directly:
gdalinfo "WCS: https://elevation.nationalmap.gov:443/arcgis/services/3DEPElevation/ImageServer/WCSServer?version=2.0.1&coverage=DEP3Elevation " This returns one band of Float32 elevation data. This is exactly what I want! Thanks to both of you for the suggestions, carl On Thu, Jun 4, 2020 at 7:55 AM Carl Godkin <cgod...@gmail.com> wrote: > Hi Michael, > > That seems perfect. I'll experiment with that. > > Thank you very much! > > carl > > On Thu, Jun 4, 2020 at 7:45 AM Michael Smith <michael.smith.e...@gmail.com> > wrote: > >> Looks like they don’t provide the data from that service just the >> extents. To get the data, they have the product view >> https://viewer.nationalmap.gov/basic/#productGroupSearch >> >> that links to files for download >> >> >> >> Or you can access them directly via gdal with something like >> >> >> >> gdalinfo /vsizip//vsicurl/ >> https://prd-tnm.s3.amazonaws.com/StagedProducts/Elevation/1/GridFloat/USGS_NED_1_n35w117_GridFloat.zip/usgs_ned_1_n35w117_gridfloat.flt >> >> >> >> Mike >> >> >> >> >> >> -- >> >> Michael Smith >> >> US Army Corps of Engineers >> >> Remote Sensing/GIS Center >> >> >> >> >> >> *From: *Carl Godkin <cgod...@gmail.com> >> *Date: *Thursday, June 4, 2020 at 10:28 AM >> *To: *<michael.smith.e...@gmail.com> >> *Cc: *<gdal-dev@lists.osgeo.org> >> *Subject: *Re: [gdal-dev] Using GDAL WMS with USGS 3DEPElevationIndex >> >> >> >> Hi Michael, >> >> >> >> Thanks a lot for your reply. >> >> >> >> I managed to come up with something sort of like that by using >> gdal_translate to a WMS service description file and then using that. >> >> >> >> But I don't see how to get elevations from there. Have you figured that >> out? I know how to get imagery, but I'm interested in elevation values. >> >> >> >> Thanks a lot, >> >> carl >> >> >> >> >> >> On Wed, Jun 3, 2020 at 5:42 PM <michael.smith.e...@gmail.com> wrote: >> >> Carl, >> >> >> >> Try something like >> >> >> >> gdal_translate -f gtiff "<GDAL_WMS><Service name=\"AGS\"><ServerUrl> >> https://index.nationalmap.gov/arcgis/rest/services/3DEPElevationIndex/MapServer</ServerUrl><BBoxOrder>xyXY</BBoxOrder><SRS>3857</SRS></Service><DataWindow><UpperLeftX>-20037508.34</UpperLeftX><UpperLeftY>20037508.34</UpperLeftY><LowerRightX>20037508.34</LowerRightX><LowerRightY>-20037508.34</LowerRightY><SizeX>2500</SizeX><SizeY>2500</SizeY></DataWindow></GDAL_WMS>" >> test1.tif -projwin ulx uly lrx lry -projwin_srs <srs of your window> >> -t_srs <srs you want for the output> >> >> >> >> If you want a different size than 2500x2500 just change that size >> >> >> >> Michael Smith >> >> US Army Corps >> >> >> >> >> >> On Jun 3, 2020, at 8:27 PM, Carl Godkin <cgod...@gmail.com> wrote: >> >> >> >> Hi, >> >> >> >> I'm unable to figure out how to retrieve elevation data (grids) from the >> 3DEP program's endpoints. >> >> >> >> Starting here: >> https://index.nationalmap.gov/arcgis/rest/services/3DEPElevationIndex/MapServer >> >> >> >> >> it would seem to be possible but I just cannot get anything back other >> than solid white images. Most of the help from that page leads off to >> ArcGIS-land which I don't have or want to use. And Googling for 3DEP leads >> only to more government sites that aren't helpful. >> >> >> >> If you have some experience with this (or a related public elevation data >> source), would you please send an example "gdal_translate" command line? >> >> >> >> I'd really like to figure out how to download terrain (as perhaps a >> GeoTIFF with floating point elevations) for an arbitrary CONUS lat/lon >> range. >> >> >> >> Thanks in advance for any help! >> >> >> >> carl >> >> _______________________________________________ >> gdal-dev mailing list >> gdal-dev@lists.osgeo.org >> https://lists.osgeo.org/mailman/listinfo/gdal-dev >> >>
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
