99 wetlands

Retrieving wetland boundaries from the National Wetlands Inventory (NWI)

The National Wetlands Inventory (NWI) is a comprehensive geospatial inventory of U.S. wetlands. It is a publicly available resource that provides detailed information on the abundance, characteristics, and distribution of wetlands in the United States. The NWI dataset is maintained by the U.S. Fish and Wildlife Service (USFWS) and is used by a wide range of stakeholders, including federal, state, and local agencies, researchers, and conservation organizations.

The notebook demonstrates how to retrieve wetland boundaries from the NWI using the leafmap Python package. The NWI dataset is available as a web service, which allows users to access the data programmatically. The leafmap package provides a simple interface for querying the NWI web service and visualizing the wetland boundaries on an interactive map.

Uncomment the following line to install leafmap if needed.

In [1]:

# %pip install -U leafmap

# %pip install -U leafmap

In [2]:

import leafmap

import leafmap

In [3]:

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14)
m.add_basemap("HYBRID")
m

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14) m.add_basemap("HYBRID") m

Out[3]:

Using point geometry¶

In [4]:

point_geometry = {"x": -99.907986, "y": 47.216359}

point_geometry = {"x": -99.907986, "y": 47.216359}

In [5]:

gdf = leafmap.get_nwi(point_geometry)

gdf = leafmap.get_nwi(point_geometry)

In [6]:

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14)
m.add_basemap("HYBRID")
m.add_nwi(gdf, add_legend=True)
m

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14) m.add_basemap("HYBRID") m.add_nwi(gdf, add_legend=True) m

Out[6]:

Using bounding box¶

In [7]:

bbox_geometry = {"xmin": -99.9023, "ymin": 47.211, "xmax": -99.8556, "ymax": 47.2325}

bbox_geometry = {"xmin": -99.9023, "ymin": 47.211, "xmax": -99.8556, "ymax": 47.2325}

In [8]:

gdf = leafmap.get_nwi(bbox_geometry)

gdf = leafmap.get_nwi(bbox_geometry)

In [9]:

m = leafmap.Map()
m.add_basemap("HYBRID")
m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True)
m

m = leafmap.Map() m.add_basemap("HYBRID") m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True) m

Out[9]:

Using GeoDataFrame¶

In [10]:

bbox = [-99.8653, 47.3952, -99.7498, 47.4441]
bbox_geometry = leafmap.bbox_to_gdf(bbox)
bbox_geometry

bbox = [-99.8653, 47.3952, -99.7498, 47.4441] bbox_geometry = leafmap.bbox_to_gdf(bbox) bbox_geometry

Out[10]:

	geometry
0	POLYGON ((-99.8653 47.3952, -99.8653 47.4441, ...

In [11]:

gdf = leafmap.get_nwi(bbox_geometry)

gdf = leafmap.get_nwi(bbox_geometry)

In [12]:

m = leafmap.Map()
m.add_basemap("HYBRID")
m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True)
m

m = leafmap.Map() m.add_basemap("HYBRID") m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True) m

Out[12]:

Using HUC8 watershed boundary¶

In [13]:

huc8 = "11120104"

huc8 = "11120104"

In [14]:

gdf = leafmap.get_nwi_by_huc8(huc8, quiet=False)

gdf = leafmap.get_nwi_by_huc8(huc8, quiet=False)

Downloading...
From: https://documentst.ecosphere.fws.gov/wetlands/downloads/watershed/HU8_11120104_Watershed.zip
To: /tmp/HU8_11120104_Watershed.zip

  0%|          | 0.00/4.36M [00:00<?, ?B/s]

100%|██████████| 4.36M/4.36M [00:00<00:00, 50.9MB/s]

Extracting files...

In [15]:

m = leafmap.Map()
m.add_basemap("HYBRID")
m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True)
m

m = leafmap.Map() m.add_basemap("HYBRID") m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True) m

Out[15]: