This notebook will demonstrate how to do basic SuperDARN data plotting.



In [1]:

    
%pylab inline
import datetime
import os
import matplotlib.pyplot as plt
from davitpy import pydarn









    



Populating the interactive namespace from numpy and matplotlib
Loaded davitpyrc file from DaViTpy installation directory. Path: /usr/local/lib/python2.7/dist-packages/davitpy-0.5-py2.7-linux-x86_64.egg/davitpy



In [2]:

    
sTime = datetime.datetime(2008,2,22)
eTime = datetime.datetime(2008,2,23)
radar = 'bks'
beam  = 7

Remote File RTI Plots



In [3]:

    
#The following command will print the docstring for the plotRti routine:
#pydarn.plotting.rti.plotRti?



In [3]:

    
fig = plt.figure(figsize=(14,12)) #Define a figure with a custom size.
pydarn.plotting.rti.plotRti(sTime, radar, eTime=eTime, bmnum=beam)
plt.show()









    



Warning: This function is being deprecated. Use plot_rti instead.
Calling plot_rti.
Found cached file: /tmp/sd/20080222.000000.20080223.000000.bks.fitex

reached end of data
plotting took: 0:00:12.922522






    



/usr/local/lib/python2.7/dist-packages/matplotlib/figure.py:397: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
  "matplotlib is currently using a non-GUI backend, "






    





<matplotlib.figure.Figure at 0x7f8503e4b590>



In [4]:

    
#Now save as a PNG to your home folder...
home = os.getenv('HOME')
filename = os.path.join(home,'rti.png')
fig.savefig(filename)



In [5]:

    
fig.clear() #Clear the figure from memory.

Local File RTI Plot

You can also plot data stored in a local file. Just change the variables in the cell below.



In [5]:

    
fileName = '/tmp/sd/20080222.000000.20080223.000000.bks.fitex'
fileType = 'fitex'
radar = 'bks'
beam  = 7
sTime = datetime.datetime(2008,2,22)
eTime = datetime.datetime(2008,2,23)

fig = plt.figure(figsize=(14,12)) #Define a figure with a custom size.
pydarn.plotting.rti.plotRti(sTime, radar, eTime=eTime, bmnum=beam, fileName=fileName,fileType=fileType)
plt.show()









    



Warning: This function is being deprecated. Use plot_rti instead.
Calling plot_rti.
problem reading /tmp/sd/20080222.000000.20080223.000000.bks.fitex :file does not exist
error, your pointer does not point to any data
error, no data available for the requested time/radar/filetype combination






    





<matplotlib.figure.Figure at 0x7f00575a2fd0>



In [8]:

    
fig.clear() #Clear the figure from memory.

Fan Plots

Right now we don't have plotFan setup to accept local files. But, we will add that in shortly.

Geopgraphic Coordinates



In [13]:

    
import datetime
import os
import matplotlib.pyplot as plt
from davitpy import pydarn

pydarn.plotting.fan.plotFan(datetime.datetime(2013,3,16,16,30),['fhe','fhw'],param='power',gsct=False)









    



Found cached file: /tmp/sd/20130316.163000.20130316.163100.fhe.fitex
Found cached file: /tmp/sd/20130316.163000.20130316.163100.fhw.fitex
Warning, setting dateTime to datetime
Warning, setting dateTime to datetime
0:00:04.784423

reached end of data

reached end of data





    



/usr/local/lib/python2.7/dist-packages/matplotlib/artist.py:221: MatplotlibDeprecationWarning: This has been deprecated in mpl 1.5, please use the
axes property.  A removal date has not been set.
  warnings.warn(_get_axes_msg, mplDeprecation, stacklevel=1)
/usr/local/lib/python2.7/dist-packages/matplotlib/figure.py:397: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
  "matplotlib is currently using a non-GUI backend, "






    









    





<matplotlib.figure.Figure at 0x7efff893ad50>

Magnetic Coordinates

Magnetic coordinates still need a little work. For instance, high latitude continent lines don't always plot. Also, we are working on getting Simon's new AACGM system in place (http://dx.doi.org/doi/10.1002/2014JA020264). Not there yet...



In [14]:

    
import datetime
import os
import matplotlib.pyplot as plt
from davitpy import pydarn

pydarn.plotting.fan.plotFan(datetime.datetime(2014,9,25,6,2),['cvw','cve','fhe'],param='velocity',gsct=False,coords='mag')









    



Found cached file: /tmp/sd/20140925.060200.20140925.060300.cvw.fitex
Found cached file: /tmp/sd/20140925.060200.20140925.060300.cve.fitex
Found cached file: /tmp/sd/20140925.060200.20140925.060300.fhe.fitex
Warning, setting dateTime to datetime
Warning, setting dateTime to datetime
0:00:04.400967

reached end of data

reached end of data

reached end of data






    





<matplotlib.figure.Figure at 0x7efff7a57b50>

Convection Plotting



In [7]:

    
import datetime
import matplotlib.pyplot as plt
import davitpy.pydarn.plotting.plotMapGrd
from davitpy.utils import *

fig = plt.figure(figsize=(15,15))
ax = fig.add_subplot(111)

sdate = datetime.datetime(2014,9,25,6,2)
mObj = plotUtils.mapObj(boundinglat=50.,gridLabels=True, coords='mag')
mapDatObj = davitpy.pydarn.plotting.plotMapGrd.MapConv(sdate, mObj, ax)
mapDatObj.overlayMapFitVel()
#mapDatObj.overlayCnvCntrs()
mapDatObj.overlayHMB()









    



Warning, datetime/dateTime not specified, using current time.
Found cached file: /tmp/sd/20140925.060200.20140925.060400.north.grdex
Found cached file: /tmp/sd/20140925.060200.20140925.060400.north.mapex



In [ ]: