... | ... | @@ -36,11 +36,18 @@ for pt in pts : |
|
|
print(new_pt)
|
|
|
```
|
|
|
# Cell fraction to area
|
|
|
## 30 minutes resolution
|
|
|
```
|
|
|
gridhectares = np.zeros((len(latitudes),len(longitudes)))
|
|
|
for lat,latitude in enumerate(latitudes) :
|
|
|
gridhectares[lat,:] = 110 * 110 * math.cos(math.pi*latitude/180) * 1/2 * 1/2 * 100
|
|
|
```
|
|
|
## 5 minutes resolution
|
|
|
```
|
|
|
surfaces = np.zeros_like(cftfracs)
|
|
|
for lat,latitude in enumerate(latitudes) :
|
|
|
surfaces[:,:,lat,:] = 110 * 110 * math.cos(math.pi*latitude/180) * 0.083 * 0.083 * 100
|
|
|
```
|
|
|
# Moving average
|
|
|
```
|
|
|
ma_harvest = []
|
... | ... | |