-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathget_geometry.py
113 lines (102 loc) · 2.61 KB
/
get_geometry.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
import ee
def get_bounds(region_name):
if region_name == 'afghan':
geom = ee.Geometry.Polygon(
[
[
70.27809406204916,
34.043798365898766
],
[
70.79994464798666,
34.043798365898766
],
[
70.79994464798666,
34.33233747415387
],
[
70.27809406204916,
34.33233747415387
],
[
70.27809406204916,
34.043798365898766
]
], 'EPSG:4326', False
)
elif region_name == 'sudan':
geom = ee.Geometry.Polygon(
[
[
[
33.09961652873555,
14.09063032622901
],
[
33.6022410404543,
14.09063032622901
],
[
33.6022410404543,
14.646707097666118
],
[
33.09961652873555,
14.646707097666118
],
[
33.09961652873555,
14.09063032622901
]
]
], 'EPSG:4326', False
)
elif region_name == 'iran':
geom = ee.Geometry.Polygon(
[
[
[
48.556574946261954,
32.54634289750261
],
[
48.018244868136954,
32.54634289750261
],
[
48.03060448727758,
32.04485995394506
],
[
48.46731103024633,
32.04485995394506
],
[
48.559321528293204,
32.08501075292506
],
[
48.58043587765844,
32.10827850580867
],
[
48.59914696774633,
32.13212180060732
],
[
48.59708703122289,
32.424130475331516
],
[
48.57236779294164,
32.47280364604964
],
[
48.556574946261954,
32.54634289750261
]
]
], 'EPSG:4326', False
)
return geom