SunPy: Python for Solar Physics
Software in Solar Physics Data Analysis
Many IDL
IDL + "SolarSoft", other people starting to use MATLAB or PDL and writing their own code to get the job done.
The code base is becoming more fragmented, and even harder to share code.
SunPy: A Gateway into Scientific Python
SunPy aims to make it easy for solar physicists to utilise the power of scientific python. (Or at least that is the idea!)
SunPy's Tools
Data Acquisition
In [6]:
from sunpy.net import vso
In [7]:
client = vso.VSOClient()
In [8]:
resp = client.query(vso.attrs.Time('2011/9/12T20:30:00', '2011/9/12T21:00:00'),
vso.attrs.Wave(304,304), vso.attrs.Instrument('AIA'),
vso.attrs.Sample(30))
print len(resp)
In [4]:
results = client.get(resp, methods=('URL-FILE_Rice'))
In [5]:
results.wait(progress=True)
Out[5]:
SunPy's Database Module
SunPy has a database moudle, which allows easy management of local data.
In [18]:
import sunpy.database
from sunpy.database.tables import display_entries
In [22]:
db = sunpy.database.Database('sqlite:///')
In [23]:
db.fetch(vso.attrs.Time('2011/9/12T20:30:00', '2011/9/12T21:00:00'),
vso.attrs.Wave(171,171), vso.attrs.Instrument('AIA'),
vso.attrs.Sample(300))
db.commit()
In [24]:
print display_entries(db,
['id', 'observation_time_start', 'observation_time_end', 'instrument', 'wavemin', 'wavemax'])
Query the database:
In [27]:
rs = db.query(vso.attrs.Time('2011/9/12T20:45:00', '2011/9/12T21:00:00'))
In [28]:
print display_entries(rs,
['id', 'observation_time_start', 'observation_time_end', 'instrument', 'wavemin', 'wavemax'])
Data Types
Map
In [20]:
import sunpy.map
In [3]:
aia_map = sunpy.map.Map('/home/stuart/sunpy/data/aia_lev1_171a_2012_07_23t00_00_11_34z_image_lev1_fits')
aia_map = aia_map.rotate()
In [6]:
fig = plt.figure(figsize=(10,10))
im = aia_map.plot()
li = aia_map.draw_limb()
gr = aia_map.draw_grid()
Map Cubes
In [10]:
mapc = sunpy.map.Map('/home/stuart/sunpy/data/aia_lev1_304a*', cube=True)
# Crop the images
mapc = sunpy.map.Map([amap.submap([-1200, -500], [-200, 700]) for amap in mapc.maps], cube=True)
In [18]:
fig = plt.figure(figsize=(10,12))
ani = mapc.plot()
ani.save("test.mp4", bitrate=5000)
plt.close()
LightCurve
LightCurve is about to be re-written! This is using SunPy 0.5!
In [21]:
import sunpy.lightcurve
In [22]:
goeslc = sunpy.lightcurve.GOESLightCurve.create('2014/6/11', '2014/6/12')
In [23]:
fig = plt.figure(figsize=(15,8))
ax = goeslc.plot(figure=fig)
Spectra / Spectrogram
In [21]:
from sunpy.spectra.sources.callisto import CallistoSpectrogram
In [23]:
spec = CallistoSpectrogram.read(sunpy.CALLISTO_IMAGE)
In [32]:
plt.figure(figsize=(12,8))
_ = spec.plot()