-
Notifications
You must be signed in to change notification settings - Fork 1
/
davitpy_sd_fov_plot.py
46 lines (41 loc) · 2.3 KB
/
davitpy_sd_fov_plot.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
import matplotlib
matplotlib.use("Agg")
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
import datetime as dt
from davitpy.utils import plotUtils
from davitpy.utils.plotUtils import textHighlighted
from davitpy.pydarn.plotting.mapOverlay import overlayRadar, overlayFov
dateTime = dt.datetime(2015, 1, 1)
fig = plt.figure(figsize=(10,5), dpi=300)
ax = fig.add_subplot(121)
m = plotUtils.mapObj(projection='npstere', lon_0=270, boundinglat=30,
coords="geo", dateTime=dateTime, ax=ax, fill_alpha=.7, )
codes = ['bks', 'fhe','fhw','cve','cvw',"wal", "ade", "adw", "hok", "hkw"]
overlayFov(m, codes=codes, dateTime=dateTime, fovColor="darkorange", fovAlpha=0.4, maxGate=75, lineWidth=0.5)
codes = ['kap','sas','pgr',"gbr", "ksr","kod","sto","pyk","han"]
overlayFov(m, codes=codes, dateTime=dateTime, fovColor="aqua", fovAlpha=0.4, maxGate=75, lineWidth=0.5)
codes = ["inv","cly","rkn","lyr"]
overlayFov(m, codes=codes, dateTime=dateTime, fovColor="lime", fovAlpha=0.4, maxGate=75, lineWidth=0.5)
codes = ['bks', 'fhe','fhw','cve','cvw',"wal", "ade", "adw", "hok", "hkw"]
overlayRadar(m, fontSize=30, codes=codes, dateTime=dateTime)
codes = ['kap','sas','pgr',"gbr", "ksr","kod","sto","pyk","han"]
overlayRadar(m, fontSize=30, codes=codes, dateTime=dateTime)
codes = ["inv","cly","rkn","lyr"]
overlayRadar(m, fontSize=30, codes=codes, dateTime=dateTime)
ax.set_title("Northern Hemisphere", fontdict={"size":20,})
ax = fig.add_subplot(122)
m = plotUtils.mapObj(projection='spstere', lon_0=270, boundinglat=-35,
coords="geo", dateTime=dateTime, ax=ax, fill_alpha=.7, )
codes = ["bpk", "tig", "unw", "fir"]
overlayRadar(m, fontSize=30, codes=codes, dateTime=dateTime)
overlayFov(m, codes=codes, dateTime=dateTime, fovColor="darkorange", fovAlpha=0.4, maxGate=75, lineWidth=0.5)
codes = ["ker","zho","sye","sys","sps","san","dce"]
overlayRadar(m, fontSize=30, codes=codes, dateTime=dateTime)
overlayFov(m, codes=codes, dateTime=dateTime, fovColor="aqua", fovAlpha=0.4, maxGate=75, lineWidth=0.5)
codes = ["mcm"]
overlayRadar(m, fontSize=30, codes=codes, dateTime=dateTime)
overlayFov(m, codes=codes, dateTime=dateTime, fovColor="lime", fovAlpha=0.4, maxGate=75, lineWidth=0.5)
ax.set_title("Southern Hemisphere", fontdict={"size":20,})
fig.savefig("data/overview.png", bbox_inches="tight")