[gdal-dev] Using GDAL WMS with USGS 3DEPElevationIndex

[Carl Godkin]

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.erdc at 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 <cgodkin at gmail.com>
> *Date: *Thursday, June 4, 2020 at 10:28 AM
> *To: *<michael.smith.erdc at gmail.com>
> *Cc: *<gdal-dev at 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.erdc at 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 <cgodkin at 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 at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200604/592ef694/attachment.html>