[Anuga-user] How to set up topography from ADF file downloaded from USGS

[Fei L. <Fe...@NJ...>]

Hi,

1) I have an ADF file downloaded from USGS website:

https://viewer.nationalmap.gov/basic/#productSearch

https://prd-tnm.s3.amazonaws.com/StagedProducts/Elevation/1/ArcGrid/USGS_NED_1_n43w084_ArcGrid.zip

2) After unzipping this file, the topography file is this 55MB adf data 
file

file grdn43w084_1/w001001.adf
grdn43w084_1/w001001.adf: ESRI Shapefile version 0 length 27774971 type 
Null Shape

ls -la grdn43w084_1/w001001.adf
-rw-r--r-- 1 feiliu feiliu 55549942 Apr 12 13:52 grdn43w084_1/w001001.adf

3) In your example you can convert

import zipfile as zf
if project.verbose: print 'Reading ASC from cairns.zip'
zf.ZipFile(project.name_stem+'.zip').extract(project.name_stem+'.asc')

# Create DEM from asc data
anuga.asc2dem(project.name_stem+'.asc', use_cache=project.cache,
               verbose=project.verbose)

# Create pts file for onshore DEM
anuga.dem2pts(project.name_stem+'.dem', use_cache=project.cache,
               verbose=project.verbose)

4) My initial reading of the user manual only found reference about 
arcgis file on page 45.

.import_ungenerate_file(ofile, tag=None, region tag=None) Module: 
pmesh.mesh This method is used to import a polygon file in the 
ungenerate format, which is used by arcGIS. The polygons from the file 
are converted to vertices and segments. ofile is the name of the polygon 
file. tag is the tag given to all the polygon’s segments. If tag is not 
supplied then the segment will not effect the water flow, it will only 
effect the mesh generation. region_tag is the tag given to all the 
polygon’s segments. If it is a string the tag will be assigned to all 
regions. If it is a list the first value in the list will be applied to 
the first polygon etc. This function can be used to import building 
footprints.

5) I am not sure how I could convert the elevation (topography) ADF file 
from USGS to DEM file. Any help would be appreciated. Thanks,

- Fei