Data-treatment.md

@@ -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