Update Data treatment authored by Marie Brunel's avatar Marie Brunel
......@@ -32,3 +32,10 @@ for pt in pts :
new_pt.append(p)
print(new_pt)
```
# Cell fraction to area
```
grid_surface = np.zeros((len(latitudes),len(longitudes))) # dim [30,33]
for lat,latitude in enumerate(latitudes) :
for lon,longitude in enumerate(longitudes) :
grid_surface[lat][lon] = 110 * 110 * math.cos(math.pi*latitude/180) * 1/2 * 1/2 * 100 #in ha
```
\ No newline at end of file