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
