43 extract values to points

Uncomment the following line to install geemap if needed.

In [ ]:

# !pip install geemap

# !pip install geemap

In [ ]:

import os
import ee
import geemap

import os import ee import geemap

Upgrade geemap to the latest version

In [ ]:

# geemap.update_package()

# geemap.update_package()

In [ ]:

Map = geemap.Map()
Map

Map = geemap.Map() Map

In [ ]:

# Add Earth Engine dataset
dem = ee.Image("USGS/SRTMGL1_003")
landsat7 = ee.Image("LANDSAT/LE7_TOA_5YEAR/1999_2003")

# Set visualization parameters.
vis_params = {
    "min": 0,
    "max": 4000,
    "palette": ["006633", "E5FFCC", "662A00", "D8D8D8", "F5F5F5"],
}

# Add Earth Engine layers to Map
Map.addLayer(
    landsat7, {"bands": ["B4", "B3", "B2"], "min": 20, "max": 200}, "Landsat 7"
)
Map.addLayer(dem, vis_params, "SRTM DEM", True, 1)

# Add Earth Engine dataset dem = ee.Image("USGS/SRTMGL1_003") landsat7 = ee.Image("LANDSAT/LE7_TOA_5YEAR/1999_2003") # Set visualization parameters. vis_params = { "min": 0, "max": 4000, "palette": ["006633", "E5FFCC", "662A00", "D8D8D8", "F5F5F5"], } # Add Earth Engine layers to Map Map.addLayer( landsat7, {"bands": ["B4", "B3", "B2"], "min": 20, "max": 200}, "Landsat 7" ) Map.addLayer(dem, vis_params, "SRTM DEM", True, 1)

Download sample data

In [ ]:

work_dir = os.path.expanduser("~/Downloads")
in_shp = os.path.join(work_dir, "us_cities.shp")
if not os.path.exists(in_shp):
    data_url = "https://github.com/giswqs/data/raw/main/us/us_cities.zip"
    geemap.download_from_url(data_url, out_dir=work_dir)

work_dir = os.path.expanduser("~/Downloads") in_shp = os.path.join(work_dir, "us_cities.shp") if not os.path.exists(in_shp): data_url = "https://github.com/giswqs/data/raw/main/us/us_cities.zip" geemap.download_from_url(data_url, out_dir=work_dir)

In [ ]:

in_fc = geemap.shp_to_ee(in_shp)
Map.addLayer(in_fc, {}, "Cities")

in_fc = geemap.shp_to_ee(in_shp) Map.addLayer(in_fc, {}, "Cities")

Export pixel values as a shapefile

In [ ]:

out_shp = os.path.join(work_dir, "dem.shp")
geemap.extract_values_to_points(in_fc, dem, out_shp)

out_shp = os.path.join(work_dir, "dem.shp") geemap.extract_values_to_points(in_fc, dem, out_shp)

Export pixel values as a csv

In [ ]:

out_csv = os.path.join(work_dir, "landsat.csv")
geemap.extract_values_to_points(in_fc, landsat7, out_csv)

out_csv = os.path.join(work_dir, "landsat.csv") geemap.extract_values_to_points(in_fc, landsat7, out_csv)