... | ... | @@ -31,4 +31,11 @@ for pt in pts : |
|
|
p = round(p)+0.25
|
|
|
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 |