# !pip install geemap

import ee
import geemap.foliumap as geemap

Map = geemap.Map(center=[40, -100], zoom=4)

# Select the seven NLCD epoches after 2000.
years = ['2001', '2004', '2006', '2008', '2011', '2013', '2016']

# Get an NLCD image by year.
def getNLCD(year):
    # Import the NLCD collection.
    dataset = ee.ImageCollection('USGS/NLCD_RELEASES/2016_REL')

    # Filter the collection by year.
    nlcd = dataset.filter(ee.Filter.eq('system:index', year)).first()

    # Select the land cover band.
    landcover = nlcd.select('landcover')
    return landcover

for year in years:
    data = getNLCD(year)
    Map.addLayer(data, {}, 'NLCD ' + year)

# Add the NLCD legend to the map.
Map.add_legend(legend_title='NLCD Land Cover Classification', builtin_legend='NLCD')

Map

Map.publish(name='National Land Cover Database (NLCD)')