-
Notifications
You must be signed in to change notification settings - Fork 50
/
reduce_resolution.py
38 lines (30 loc) · 1.08 KB
/
reduce_resolution.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
import ee
import geemap
# Create a map centered at (lat, lon).
Map = geemap.Map(center=[40, -100], zoom=4)
# Load a MODIS EVI image.
modis = ee.Image(ee.ImageCollection('MODIS/006/MOD13A1').first()) \
.select('EVI')
# Display the EVI image near La Honda, California.
Map.setCenter(-122.3616, 37.5331, 12)
Map.addLayer(modis, {'min': 2000, 'max': 5000}, 'MODIS EVI')
# Get information about the MODIS projection.
modisProjection = modis.projection()
print('MODIS projection:', modisProjection.getInfo())
# Load and display forest cover data at 30 meters resolution.
forest = ee.Image('UMD/hansen/global_forest_change_2015') \
.select('treecover2000')
Map.addLayer(forest, {'max': 80}, 'forest cover 30 m')
# Get the forest cover data at MODIS scale and projection.
forestMean = forest \
.reduceResolution(**{
'reducer': ee.Reducer.mean(),
'maxPixels': 1024
}) \
.reproject(**{
'crs': modisProjection
})
# Display the aggregated, reprojected forest cover data.
Map.addLayer(forestMean, {'max': 80}, 'forest cover at MODIS scale')
# Display the map.
Map