-
Notifications
You must be signed in to change notification settings - Fork 0
/
fuel_moisture
51 lines (42 loc) · 1.25 KB
/
fuel_moisture
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
39
40
41
42
43
44
45
46
47
48
49
50
51
var bbox = ee.Geometry.Rectangle({coords:[-180,-60,180,84],geodesic:false});
var start = ee.Date('2015-01-01');
var end = ee.Date('2018-12-31');
var MOD16A2 = ee.ImageCollection("MODIS/006/MOD16A2").filterDate(start,end);
var calc_alpha = function(image) {
return image.expression(
'(b1)/b2', {
'b1': image.select('ET'),
'b2': image.select('PET')
}).unmask(-9999)};
var alpha = MOD16A2.map(calc_alpha);
// set null values as -999
print(alpha,'alpha');
Map.addLayer(alpha);
var shapes = ee.FeatureCollection("users/dsandovalhdh/fuelmoisture_extract");
Map.addLayer(shapes);
var values = alpha.map(function(image) {
return image.reduceRegions({
collection: shapes,
reducer: ee.Reducer.mean(),
crs:'EPSG:4326',
tileScale: 16,
scale: 500
}).map(function(f){
return f.set('date', image.get('system:index'));
});
}).flatten();
print(values.limit(20),'values');
//; error, multiple points same coordinates
Export.table.toDrive({
collection:values,
folder: "fuelmoisture_db",
selectors:(["code","date","mean",".geo"]),
});
/*
Export.table.toDrive(values,
"fuelmoisture_db",
"fuelmoisture_db",
"fuelmoisture_db");
//Map.addLayer(PML,
// {min: [0, 0, 0], max: [10, 10, 20], bands: ['Es' , 'Ec','ET_water']}, 'PML');
*/