SunPy Code Tutorial for 12th RHESSI Workshop (Nanjing)

# given on 17-Oct-2011

# plotting a fits file image
import sunpy
aia = sunpy.Map(sunpy.AIA_171_IMAGE)
aia.show()

from matplotlib import colors
aia.show(norm=colors.Normalize(1, 2048))

from sunpy.cm import cm as cm
cm.show_colormaps
aia.show(cmap=cm.sdoaia94)

from matplotlib import cm

type(aia)
aia.header
aia.max()
aia.mean()

sub_aia = aia.submap([-500,500], [-500,500])
sub_aia.show()

from sunpy.net import vso

# create a new VSOClient instance
client = vso.VSOClient()

result = client.query( vso.attrs.Time((2006, 9, 20), (2006, 9, 21,1)), vso.attrs.Instrument('eit') )
result.no_records
result = client.query( vso.attrs.Time((2006, 9, 20), (2006, 9, 21,1)), vso.attrs.Instrument('eit') | vso.attrs.Instrument('trace'))

#now download the data
res = client.get(result[:1]).wait()

eit = sunpy.Map('/Users/schriste/Desktop/efz20100920.011938.fits')

from sunpy.sun import constants as con

con.volume
con.value('volume')
con.unit('volume')
con.print_all()

solar_constants = con.physical_constants
solar_constants.keys()
solar_constants.get('volume')

from sunpy.util import util as util
util.anytim('2003/04/01 01:03:00')

import pidly
idl = pidly.IDL('myidl')
idl('t = anytim_now(/yoh)')
idl.t

idl.x = 2 ** 2
idl.x

idl.interact()
print, x
x = 12
#now type ^D to escape the IDL shell
idl.x

#if time permits
from sunpy.instr import rhessi
image = rhessi.backprojection(sunpy.RHESSI_EVENT_LIST)
from pylab import *
imshow(image)