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)

---

Last update: 2023-04-06
Created: 2020-03-31