113 image area

Uncomment the following line to install geemap if needed.

In [ ]:

# !pip install geemap

# !pip install geemap

In [ ]:

import ee
import geemap

import ee import geemap

Add ESA Land Cover data.

In [ ]:

Map = geemap.Map()
dataset = ee.ImageCollection("ESA/WorldCover/v100").first()
Map.addLayer(dataset, {'bands': ['Map']}, 'ESA Land Cover')
Map.add_legend(builtin_legend='ESA_WorldCover')
Map

Map = geemap.Map() dataset = ee.ImageCollection("ESA/WorldCover/v100").first() Map.addLayer(dataset, {'bands': ['Map']}, 'ESA Land Cover') Map.add_legend(builtin_legend='ESA_WorldCover') Map

Calculate the area of each land cover type.

In [ ]:

df = geemap.image_area_by_group(
    dataset, scale=1000, denominator=1e6, decimal_places=4, verbose=True
)
df

df = geemap.image_area_by_group( dataset, scale=1000, denominator=1e6, decimal_places=4, verbose=True ) df

Save the results to a CSV.

In [ ]:

df.to_csv('esa_area.csv')

df.to_csv('esa_area.csv')

Add NLCD land cover data.

In [ ]:

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

nlcd = ee.Image('USGS/NLCD_RELEASES/2019_REL/NLCD/2019')
landcover = nlcd.select('landcover')

Map.addLayer(landcover, {}, 'NLCD Land Cover 2019')
Map.add_legend(
    title="NLCD Land Cover Classification", builtin_legend='NLCD', height='465px'
)
Map

Map = geemap.Map(center=[40, -100], zoom=4) Map.add_basemap('HYBRID') nlcd = ee.Image('USGS/NLCD_RELEASES/2019_REL/NLCD/2019') landcover = nlcd.select('landcover') Map.addLayer(landcover, {}, 'NLCD Land Cover 2019') Map.add_legend( title="NLCD Land Cover Classification", builtin_legend='NLCD', height='465px' ) Map

Calculate the area of each land cover type.

In [ ]:

df = geemap.image_area_by_group(
    landcover, scale=1000, denominator=1e6, decimal_places=4, verbose=True
)
df

df = geemap.image_area_by_group( landcover, scale=1000, denominator=1e6, decimal_places=4, verbose=True ) df

Save the results to a CSV.

In [ ]:

df.to_csv('nlcd_area.csv')

df.to_csv('nlcd_area.csv')