... | ... | @@ -49,3 +49,17 @@ for y,year in enumerate(years[window:]): |
|
|
ma_harvest.append(np.average(harvests_regions[id_region][y-window:y]))
|
|
|
ma_harvest = np.asarray(ma_harvest)
|
|
|
```
|
|
|
# Regional selection (biome in brazil)
|
|
|
|
|
|
```
|
|
|
def region_selection (latitudes,longitudes,region):
|
|
|
nc = Dataset("/home/marie/ownCloud/Documents/spitfire grassland paper/visualization/python_toolbox/brazil_biomes/brazilian_biomes_newversion.nc", mode='r')
|
|
|
lat = nc.variables['latitude'][:]
|
|
|
lon = nc.variables['longitude'][:]
|
|
|
biomes = nc.variables['biome'][:]
|
|
|
nc.close()
|
|
|
|
|
|
new_latitudes,new_longitudes = np.where(region==biomes)
|
|
|
|
|
|
return (new_latitudes,new_longitudes)
|
|
|
``` |
|
|
\ No newline at end of file |