(f) Plot EUVI Subimage Time Sequence + GOES Light Curve
For an overview of the evolution of a selected region on the Sun
(i.e., an active region) it is convenient to plot a time sequence
of EUVI subimages along with the GOES soft X-ray light curve.
We can do this taks with the IDL procedure EUVI_EVOL.PRO.
IDL>
datadir='$secchi/lz/L0/'
;disk or data archive
date ='20080325'
;'YYYYMMDD' date of observation
tstart ='184000'
;'HHMMSS' start time interval for image selection
tend ='192000'
;'HHMMSS' end time interval for image selection
sc ='a'
;spacecraft "a" or 'b'
wave ='171'
;wavelength '171','195','284', or '304' (Angstroem)
color=0
;0=automated scaling of dynamic range, 1=interactive
movie=0
;0=automated movie play, 1=interactive stepping of images
fov=[100,500,600,1100]
;[i1,j1,i2,j2] pixels of field-of-view
mismax =0
;maximum offset (in pixels) for misalignment correction by cross-correlation
iseq=[4,5,6,7]
;indices of images selected for display (iseq=[-1] for equidistant times)
nmin =3
;3 or 5 min smoothing of GOES light curve
id =-1
;ID number of event in EUVI event catalog (id=-1 if not used)
nsme =0
;(currently unused parameter)
area =[0.5,0.4,0.75,0.9]
;fractional field-of-view area of EUV light curve [x1,y1,x2,y2]
filename=date+'_'+tstart+'_'+wave+'_'+sc+'_'
;filename of plot
dir ='~/work/'
;work directory
euvi_files,datadir,date,tstart,tend,sc,wave,datafiles
euvi_display,datafiles,color,movie,images,index
euvi_fov,images,color,fov,images_fov,index,mismax
euvi_evol,date,tstart,tend,index,images_fov,iseq,sc,wave,nmin,id,nsme,area,pos,dir,filename
In this time interval we obtain a sequence 31 images.
From these 31 images (with IDL index 0,1,...,30) we choose the subset of
ISEQ=[4,5,6,7] for our plot. Let us suppose that you clicked the bottom
left corner (BLC) and the top right corner (TRC) of the FOV that encompasses
the active region at the south-east limb, say at FOV=[0,500,500,1000].
The resulting plot looks then like this:
The upper panel shows the GOES soft X-ray light curves (thin solid curves),
the time derivative, which is a proxy for the hard X-ray light curve according
to the Neupert effect (thick solid curve). The times of the 31 images are
marked with vertical dotted lines, the displayed 4 images with vertical
solid lines, and the EUV light curve (averaged over
the fractional image area (0.5-0.75 of the x-range and 0.4-0.9 of the y-range)
is shown with a dashed curve and diamond symbols. Note the EUV dimming
by 25.3% over the selected area with a size of 146 Mm in x-direction and
291 Mm in y-direction (information is given in title of plot).
The second panel shows the 4 images (of the chosen sequence
numbers ISEQ), and the third panel shows difference plots, which indicate
the dynamic changes in the flaring region (white is flux increase, while black
denotes flux decreases). Note that the difference images nicely enhance
the expanding CME bubble, which is not visible in the original images.
You should also find a postscript file of this plot
in your directory with the filename: 20080325_184000_171_a_evol_col.ps