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

[Michael Smith]

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/109a927a/attachment-0001.html>