43 extract values to points

Uncomment the following line to install geemap if needed.

In [1]:

# !pip install geemap

# !pip install geemap

In [2]:

import os
import ee
import geemap

import os import ee import geemap

Upgrade geemap to the latest version

In [3]:

# geemap.update_package()

# geemap.update_package()

In [4]:

Map = geemap.Map()
Map

Map = geemap.Map() Map

Out[4]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [5]:

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

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

# Add Earth Eninge 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('LE7_TOA_5YEAR/1999_2003') # Set visualization parameters. vis_params = { 'min': 0, 'max': 4000, 'palette': ['006633', 'E5FFCC', '662A00', 'D8D8D8', 'F5F5F5'], } # Add Earth Eninge 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 [6]:

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)

Downloading https://github.com/giswqs/data/raw/main/us/us_cities.zip ...
Unzipping us_cities.zip ...
Data downloaded to: /home/runner/Downloads/us_cities

In [7]:

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 [8]:

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)

Generating URL ...
Downloading data from https://earthengine.googleapis.com/v1alpha/projects/earthengine-legacy/tables/e7cc35bb2456e91e4e610421708c47d8-34aecdcdaa5f3b4422e6c54778b9db5f:getFeatures
Please wait ...
Data downloaded to /home/runner/Downloads/dem.shp

Export pixel values as a csv

In [9]:

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)

Generating URL ...
Downloading data from https://earthengine.googleapis.com/v1alpha/projects/earthengine-legacy/tables/8269818bfe0b6b2e8fcc62f574f89136-acfc86e26b483bd7f826de4abb37d9be:getFeatures
Please wait ...
Data downloaded to /home/runner/Downloads/landsat.csv

---

Last update: 2022-03-14