This page (revision-13) was last changed on 07-Dec-2016 14:14 by Peter Young

This page was created on 24-Feb-2010 18:15 by PeterYoung

Only authorized users are allowed to rename pages.

Only authorized users are allowed to delete pages.

Page revision history

Version Date Modified Size Author Changes ... Change note
13 07-Dec-2016 14:14 2 KB Peter Young to previous
12 24-Aug-2012 22:18 2 KB Peter Young to previous | to last
11 24-Aug-2012 22:15 2 KB Peter Young to previous | to last
10 28-Jun-2011 19:17 2 KB Peter Young to previous | to last
9 04-Aug-2010 21:40 1 KB PeterYoung to previous | to last
8 07-Jul-2010 19:34 1 KB PeterYoung to previous | to last
7 17-Mar-2010 21:33 1 KB PeterYoung to previous | to last
6 24-Feb-2010 22:22 1 KB PeterYoung to previous | to last
5 24-Feb-2010 22:19 1 KB PeterYoung to previous | to last
4 24-Feb-2010 22:09 1 KB PeterYoung to previous | to last
3 24-Feb-2010 18:32 1 KB PeterYoung to previous | to last
2 24-Feb-2010 18:29 1 KB PeterYoung to previous | to last
1 24-Feb-2010 18:15 1 KB PeterYoung to last

Page References

Incoming links Outgoing links

Version management

Difference between version and

At line 1 changed one line
!!!Worksheet 2 - Reading and browsing data
!!!Worksheet 2 - Calibrating EIS data (eis_prep)
At line 3 changed one line
The EIS FITS file to be used in all the worksheets is from 9-Dec-2007, 11:30 UT:
Generally it is recommended that you calibrate your data file before you begin taking a look at it. This is because EIS data are affected by warm pixels on the CCD which make it difficult to study the images.
At line 5 changed one line
filename=’eis_l0_20061209_113031.fits’
To calibrate your data, do the standard call to eis_prep is
At line 7 changed one line
(NOTE: you may need to add the path name, depending on where the file is located.)
IDL> eis_prep,filename,/default,/save,/retain,/quiet
At line 9 changed one line
To read an EIS FITS file from the command line, type:
To see what these keywords do (and see other options), go to the wiki section on eis_prep keywords.
At line 11 changed one line
IDL> data=obj_new(‘eis_data’,filename)
The keyword /default makes the routine use the default processing options, and the keyword /save makes eis_prep create two files in your current working directory. These are the level-1 FITS file, and the error FITS file.
At line 13 changed one line
data is an IDL object. You can extract information from the object using various methods. E.g.,
TIP: it is a good idea to compress these files with gzip and move them to a directory called /level1.
At line 15 changed 3 lines
IDL> exp=data->getexp()\\
IDL> xpos=data->getxpos()\\
IDL> wd0=data->getvar(0)\\
NOTE: “bad” data (referred to as missing) are flagged with a value of -100 in the error FITS file arrays (not the level-1 arrays). The missing data in the level-1 arrays are replaced with the median value of the surrounding pixels. This is useful for display purposes, but for scientific analysis you should not use any pixels flagged as missing.
With the /default option, eis_prep does the following things:
At line 19 changed one line
wd0 is a 3D array containing the window data for the specified window (0 in this case). You can get a list of all methods by doing:
* Flag any saturated data as missing (saturated data have a DN value of 16383)
* Remove the CCD pedestal and dark current
* Flag hot pixels as missing
* Flag warm pixels as missing (not implemented at present)
* Flag cosmic rays as missing
* Calibrate data, and create error arrays
At line 21 removed 15 lines
IDL> data->display_methods
To extract the header into a string array, do:
IDL> hdr_obj=data->gethdr()\\
IDL> text=hdr_obj->getbte_hdr()
A number of useful header tags can be sent to a structure by doing:
IDL> eis_get_hdr_struc,data,hdrstr
To find which wavelength windows are contained in data, do:
IDL> wininfo=eis_get_wininfo(filename,/list)
At line 26 added one line
#Using commands from Worksheet 2, how many pixels in the Fe XII 195 data array have been flagged as missing by eis_prep?
At line 38 removed 8 lines
#Using eis_get_wininfo and the getvar() object method, extract an image in Fe XII 195 by summing all wavelength pixels in the window. Plot the image. What do you have to do to the image to orient it correctly? (You may need to refer to XRT images from Day 1 of the workshop.)
#Using the following commands:
IDL> aux_data=data->getaux_data()\\
IDL> exptimes=data->sec_from_obs_start(*(aux_data.ti_1))\\
IDL> xpos=data->getxpos()\\
work out in which direction EIS rasters (east-west, or west-east).