Creates path to be used by GeoAxes.
Args:
minlon (float): minimum longitude
maxlon (float): maximum longitude
minlat (float): minimum latitude
maxlat (float): maximum latitude
n (int, optional): Interpolation points for each segment. Defaults to 50.
Returns:
boundary_path (mpath.Path): Boundary Path in flat projection
Source code in jetutils/plots.py
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248 | def make_boundary_path(
minlon: float, maxlon: float, minlat: float, maxlat: float, n: int = 50
) -> mpath.Path:
"""Creates path to be used by GeoAxes.
Args:
minlon (float): minimum longitude
maxlon (float): maximum longitude
minlat (float): minimum latitude
maxlat (float): maximum latitude
n (int, optional): Interpolation points for each segment. Defaults to 50.
Returns:
boundary_path (mpath.Path): Boundary Path in flat projection
"""
boundary_path = []
# North (E->W)
if maxlon > minlon:
linspace = np.linspace(minlon, maxlon, n)
else:
linspace = (180 + np.linspace(minlon % 360, maxlon % 360, n)) % 360 - 180
edge = [linspace, np.full(n, maxlat)]
boundary_path += [[i, j] for i, j in zip(*edge)]
# West (N->S)
edge = [np.full(n, maxlon), np.linspace(maxlat, minlat, n)]
boundary_path += [[i, j] for i, j in zip(*edge)]
# South (W->E)
edge = [linspace[::-1], np.full(n, minlat)]
boundary_path += [[i, j] for i, j in zip(*edge)]
# East (S->N)
edge = [np.full(n, minlon), np.linspace(minlat, maxlat, n)]
boundary_path += [[i, j] for i, j in zip(*edge)]
boundary_path = mpath.Path(boundary_path)
return boundary_path
|