pygeohydro.levee#

Accessing National Flood Hazard Layers (NLD) through web services.

Module Contents#

class pygeohydro.levee.NLD(layer, outfields='*', crs=4326)#

Access National Levee Database (NLD) services.

Notes

For more info visit: https://geospatial.sec.usace.army.mil/server/rest/services/NLD2_PUBLIC/FeatureServer

Parameters:

  • layer (str, optional) – A valid service layer. Valid layers are:

    • boreholes

    • crossings

    • levee_stations

    • piezometers

    • pump_stations

    • relief_wells

    • alignment_lines

    • closure_structures

    • cross_sections

    • embankments

    • floodwalls

    • frm_lines

    • pipe_gates

    • toe_drains

    • leveed_areas

    • system_routes

    • pipes

    • channels

  • outfields (str or list, optional) – Target field name(s), default to “*” i.e., all the fields.

  • crs (str, int, or pyproj.CRS, optional) – Target spatial reference, default to EPSG:4326.

bygeom(geom, geo_crs=4326, sql_clause='', distance=None, return_m=False, return_geom=True)#

Get features within a geometry that can be combined with a SQL where clause.

byids(field, fids, return_m=False, return_geom=True)#

Get features by object IDs.

bysql(sql_clause, return_m=False, return_geom=True)#

Get features using a valid SQL 92 WHERE clause.

Examples

>>> from pygeohydro import NLD
>>> nld = NLD("levee_stations")
>>> levees = nld.bygeom((-105.914551, 37.437388, -105.807434, 37.522392))
>>> levees.shape
(1838, 12)