mirror of
https://github.com/google/earthengine-api.git
synced 2025-12-08 19:26:12 +00:00
50 lines
1.3 KiB
Python
50 lines
1.3 KiB
Python
#!/usr/bin/env python
|
|
"""Map an expression.
|
|
|
|
Computes the mean NDVI and SAVI by mapping an expression over a collection
|
|
and taking the mean. This intentionally exercises both variants of
|
|
Image.expression.
|
|
"""
|
|
|
|
import datetime
|
|
import ee
|
|
import ee.mapclient
|
|
|
|
ee.Initialize()
|
|
|
|
# Filter the L7 collection to a single month.
|
|
collection = (ee.ImageCollection('LANDSAT/LE07/C01/T1_TOA')
|
|
.filterDate(datetime.datetime(2002, 11, 1),
|
|
datetime.datetime(2002, 12, 1)))
|
|
|
|
|
|
def NDVI(image):
|
|
"""A function to compute NDVI."""
|
|
return image.expression('float(b("B4") - b("B3")) / (b("B4") + b("B3"))')
|
|
|
|
|
|
def SAVI(image):
|
|
"""A function to compute Soil Adjusted Vegetation Index."""
|
|
return ee.Image(0).expression(
|
|
'(1 + L) * float(nir - red)/ (nir + red + L)',
|
|
{
|
|
'nir': image.select('B4'),
|
|
'red': image.select('B3'),
|
|
'L': 0.2
|
|
})
|
|
|
|
vis = {
|
|
'min': 0,
|
|
'max': 1,
|
|
'palette': [
|
|
'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163',
|
|
'99B718', '74A901', '66A000', '529400', '3E8601',
|
|
'207401', '056201', '004C00', '023B01', '012E01',
|
|
'011D01', '011301'
|
|
]}
|
|
|
|
ee.mapclient.centerMap(-93.7848, 30.3252, 11)
|
|
ee.mapclient.addToMap(collection.map(NDVI).mean(), vis)
|
|
ee.mapclient.addToMap(collection.map(SAVI).mean(), vis)
|
|
|