<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN"> <!--Converted with LaTeX2HTML 98.1p1 release (March 2nd, 1998) originally by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds * revised and updated by: Marcus Hennecke, Ross Moore, Herb Swan * with significant contributions from: Jens Lippmann, Marek Rouchal, Martin Wilck and others --> <HTML> <HEAD> <TITLE>Flat Field correction</TITLE> <META NAME="description" CONTENT="Flat Field correction"> <META NAME="keywords" CONTENT="vol2"> <META NAME="resource-type" CONTENT="document"> <META NAME="distribution" CONTENT="global"> <META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1"> <LINK REL="STYLESHEET" HREF="vol2.css"> <LINK REL="previous" HREF="node473.html"> <LINK REL="up" HREF="node469.html"> <LINK REL="next" HREF="node475.html"> </HEAD> <BODY > <!--Navigation Panel--> <A NAME="tex2html7361" HREF="node475.html"> <IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="icons.gif/next_motif.gif"></A> <A NAME="tex2html7358" HREF="node469.html"> <IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="icons.gif/up_motif.gif"></A> <A NAME="tex2html7354" HREF="node473.html"> <IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="icons.gif/previous_motif.gif"></A> <A NAME="tex2html7360" HREF="node1.html"> <IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents" SRC="icons.gif/contents_motif.gif"></A> <BR> <B> Next:</B> <A NAME="tex2html7362" HREF="node475.html">Reduction without Standard Star</A> <B> Up:</B> <A NAME="tex2html7359" HREF="node469.html">Reduction using Standard Stars</A> <B> Previous:</B> <A NAME="tex2html7355" HREF="node473.html">Background Correction</A> <BR> <BR> <!--End of Navigation Panel--> <H3><A NAME="SECTION0027101500000000000000"> Flat Field correction</A> </H3> <P> After successful completion of the command <TT>CALIBRATE/ECHELLE</TT>, the next step in the reduction is the preparation of the <TT>FLAT</TT>s, i.e. the subtraction of the background level and the preparation of the file which is needed to normalise each <TT>FLAT</TT>. <P> The flat field correction is an optional step, it may not be needed in the case of objects exposed in the short wavelength range of the echelle (images centered at <!-- MATH: $\lambda<5000$ --> <IMG WIDTH="91" HEIGHT="41" ALIGN="MIDDLE" BORDER="0" SRC="img956.gif" ALT="$\lambda<5000$">Å). If you do not want to perform a flat field correction, set the echelle keyword <TT>FFOPT</TT> to <TT>NO</TT> and skip this section. Otherwise, type: <P> <BLOCKQUOTE> <TT>SET/ECHELLE FFOPT=YES FLAT=... </TT></BLOCKQUOTE> <P> It is advised to replace the default names ffcorr and blaze for the output files by some readily identifiable ones via the command: <P> <BLOCKQUOTE> <TT>SET/ECHELLE CORRECT=... BLAZE=...</TT></BLOCKQUOTE> <P> All paramters corresponding to the Background, Extraction and Rebin (sections 2, 4, 7) in the <TT>SHOW/ECHELLE</TT> must be set. See <TT>MIDAS</TT> User Manual, Vol. B, Chapter 8 (Echelle Spectra) for more details on the different background and extraction methods. <P> The command to reduce the <TT>FLAT</TT>s is simply: <P> <BLOCKQUOTE> <TT>FLAT/ECHELLE [flat] [correct] [blazecorr]</TT></BLOCKQUOTE> <P> If you are following the tutorial, a plot will appear after a while in order to monitor the quality of the process. You will see a vertical trace of the raw flat field on top of the fitted background. After a short while a frame will appear on the image display showing the residuals of the fit at the grid points which have been used to determine the background. As explained in Section <B>Background definition</B> the background is fitted to points located in the interorder space.This residual map is stored in the intermediate frame <TT>&Z</TT> <P> The sequence described above should be applied both to the standard stars and objects with their corresponding flatfield and wavelength calibration images. Use the command <TT>SAVE/ECHELLE name</TT> to save the relevant files containing the result of the calibrations. The next step is the computation of the instrument response. It will involve the observed standard star(s). If you do not have one available, you could still use this procedure assuming that your object is itself a standard star with flux unity. If you are following the tutorial, the table <TT>LTT1020</TT> with absolute fluxes of this standard will be copied from the area <TT>MID_CASPEC</TT> into your work space by the command <TT>SET/ECHELLE FLUX=LTT1020</TT>. Other tables are also available or you can prepare your own table or use the table <TT>UNITY</TT>. <P> The following parameters have to be supplied: <TT>SLIT</TT>, <TT>OFFSET</TT>, <TT>SAMPLE</TT> and <TT>FLUXTAB</TT>, the identification of the standard star table which is to be used. To determine <TT>SLIT</TT> and <TT>OFFSET</TT> use the same procedure that was used on the <TT>FLAT</TT> to determine <TT>WIDTH1</TT>. To determine the <TT>OFFSET</TT> of the location of the orders for the object with respect to the flatfield the same column should of course be specified. <TT>SAMPLE</TT> determines the step with which the data in the final image will be sampled. When all the parameters are set, give the command: <BLOCKQUOTE> <TT>RESPONSE/ECHELLE</TT></BLOCKQUOTE>This command performs the following tasks. It determines the background in the standard star and subtracts it. Then it divides by the <TT>FLAT</TT>, extracts the orders and rebins them. The counts are reduced to an exposure time of one second and the orders are binned into the same wavelength intervals as the corresponding standard star table (at present this means 12 Å intervals). The table values are divided by the observed counts to give the conversion to absolute flux units, the response. The response frame is lastly interpolated to the resolution required by the parameter <TT>SAMPLE</TT> using a low order polynomial. In order to cancel effects introduced by differences in the <TT>FLAT</TT> for the standard star and the object, the flatfield normalisation is applied. To check the accuracy of the calibration, the standard star is reduced as if it was an object. The response correction is applied to the <TT>STD</TT> and the individual orders are merged to form a one-dimensional spectrum. A description of these last steps is given in the part which describes the command <TT>REDUCE/ECHELLE</TT>. The instrument response is called <TT>RESPONSE</TT> by default, but you could assign a different name, with the command<BLOCKQUOTE> <TT>SET/ECHELLE RESPONSE=yourname</TT></BLOCKQUOTE>before executing <TT>RESPONSE/ECHELLE</TT>. If you do not want to correct for the instrument response, assign the value <TT>NO</TT> to the <TT>RESPOPT</TT> parameter as<BLOCKQUOTE> <TT>SET/ECHELLE RESPOPT=NO</TT></BLOCKQUOTE>and skip the <TT>RESPONSE/ECHELLE</TT> command. <P> <TT>REPONSE/ECHELLE</TT> does not require any user interaction. It is advised, however, to monitor the intermediate results which appear at various stages on the screen and to inspect some of the intermediate files. Two steps are particularly delicate. The first one is the fitting of the background. As already mentioned before one should not specify too high an order for the 2D-polynomial. If <TT>BIAS</TT> and <TT>DARK</TT> frames have been subtracted a mere off-set (degree 0,0) might be sufficient, of course depending on the exposure level. A careful inspection of the residuals of the background fitting as produced by <TT>BACKGROUND/ECHELLE</TT> will help in deciding the optimal choice. Small errors in the background fit will, at the edges of the orders be amplified due to the correction for the blaze. This may critically influence the results on the final, merged spectrum. <P> The second delicate step is the interpolation of the response frame to the required wavelength step. The sampling step of the table is 12 Å. This means about 10 calibration points per order. However discrepant pixels in the response frame which occur either at the end or beginning of an order can influence the fit. The effect is therefore again serious because the correction at the edges of the orders will be larger. <P> <TT>RESPONSE/ECHELLE</TT> will display on the monitor the resulting frame showing the response of the instrument at the resolution of the standard star spectrum specified by <TT>SAMPLE</TT>. Any discrepant pixels can be readily identified. To remove these, use the command<BLOCKQUOTE> <TT>REPEAT/ECHELLE</TT></BLOCKQUOTE>(It is not a <I>nice</I> name but you are repeating some of the steps <IMG WIDTH="30" HEIGHT="21" ALIGN="BOTTOM" BORDER="0" SRC="img957.gif" ALT="$\ldots$">) This command will execute <TT>MODIFY/AREA</TT>. Identify, using the cursor, all pixels which should be excluded from the polynomial fit. When this is done, the command will execute the same set of steps as <TT>RESPONSE/ECHELLE</TT> and ultimately produce a new version of the <TT>RESPONSE</TT> file. <P> Even after this editing, the resulting response might still lead to unsatisfactory results. The reason is quite obvious. The sampling step of the standard star table is too large and therefore the number of reference pixels per order is too small, especially if one moves towards the blue end of the spectrum or if one has used the binned read-out mode. A possible way out of this problem is to use standard star tables sampled at about 3 Å instead of 12 Å . In the current version, standard star tables sampled at lower steps are obtained by interpolating the 12 Å tables. Note that the intrinsic resolution of the tables remains low and narrow (absorption) lines which are present in the standard star observation will not be present in the table. Test runs have shown that the response computed using a 3 Å sampling is superior to the one using the tables sampled at 12 Å . Absorption lines which are not present in the table give of course a wrong response correction but these regions can be easily edited using <TT>REPEAT/ECHELLE</TT>. Real high resolution tables will become available in future versions. <P> At this stage, everything is ready to reduce your observed spectrum. To reduce the object give the command <TT>INIT/ECHELLE</TT> to initialise the tables belonging to the object. Instead of specifying the wavelength step in Å one should now specify a reference file which determines the sampling of the final spectrum. This is normally the response file. Do not forget to specify also the parameter <TT>RESPONSE</TT>. <P> All parameters in sections Background, Extraction, Rebin, Flat-Field Correction, Response Correction, Merging as displayed by <TT>SHOW/ECHELLE</TT> must be checked. <P> Now type: <BLOCKQUOTE> <TT>REDUCE/ECHELLE input output</TT></BLOCKQUOTE>where <TT>input</TT> and <TT>output</TT> are the input object frame and the one-dimensional output spectrum respectively. <TT>REDUCE/ECHELLE</TT> follows essentially the same steps as <TT>RESPONSE/ECHELLE</TT>. The background is determined and subtracted, the image is divided by the flatfield and the echelle orders are extracted and rebinned; the counts are normalised to an exposure time of one second. The orders are not yet corrected for the differential extinction between <TT>STD</TT> and <TT>OBJ</TT>. They are multiplied by the response and by the flatfield normalisation. As in <TT>RESPONSE/ECHELLE</TT> the determination of the background is a delicate step which needs careful monitoring. <P> Individual orders can be extracted using the level 1 command<BLOCKQUOTE> <TT>MERGE/ECHELLE & name order1,order2 NOAPPEND</TT></BLOCKQUOTE>where <TT>&</TT> is a temporary file name (WARNING: This file will be deleted by the system when you initialise MIDAS again) created by <TT>REDUCE/ECHELLE</TT>. The parameter <TT>name</TT> is the root for the output file names to which the order number will be appended, parameters <TT>order1,order2</TT> define the order range and <TT>NOAPPEND</TT> defines this ''method of merging`` orders. You will find all this in the tutorial. The individual orders can be plotted to check the accuracy of the method in the region where the orders overlap. <P> These steps, only four high level commands plus the corresponding <TT>SET/ECHELLE</TT> definitions, are the standard reduction. Depending on the type of observations, the operation can be optimised by running the commands in batch mode. Use the interactive command <TT>CALIBRATE/ECHELLE</TT> at a MIDAS work station and save the session status with the command <BLOCKQUOTE> <TT>SAVE/ECHELLE name</TT></BLOCKQUOTE>where <TT>name </TT> is the session name. The rest of the reduction can be done completely in batch mode at some other time, probably at night, creating a procedure like <BLOCKQUOTE> <TT>INIT/ECHELLE name1 <BR> SHOW/ECHELLE <BR> FLAT/ECHELLE <BR> RESPONSE/ECHELLE <BR> REDUCE/ECHELLE input1 output1 <BR> REDUCE/ECHELLE ... <BR> INIT/ECHELLE name2 <BR> ...</TT></BLOCKQUOTE>where <TT>name1 ...</TT> are names of previously saved sessions, and <TT>input1 ...</TT>, <TT>output1 ...</TT> are the raw image and reduced spectrum, respectively. The command <TT>RESPONSE/ECHELLE</TT> can be excluded from the batch procedure given the interactive editing required. <P> <HR> <!--Navigation Panel--> <A NAME="tex2html7361" HREF="node475.html"> <IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="icons.gif/next_motif.gif"></A> <A NAME="tex2html7358" HREF="node469.html"> <IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="icons.gif/up_motif.gif"></A> <A NAME="tex2html7354" HREF="node473.html"> <IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="icons.gif/previous_motif.gif"></A> <A NAME="tex2html7360" HREF="node1.html"> <IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents" SRC="icons.gif/contents_motif.gif"></A> <BR> <B> Next:</B> <A NAME="tex2html7362" HREF="node475.html">Reduction without Standard Star</A> <B> Up:</B> <A NAME="tex2html7359" HREF="node469.html">Reduction using Standard Stars</A> <B> Previous:</B> <A NAME="tex2html7355" HREF="node473.html">Background Correction</A> <!--End of Navigation Panel--> <ADDRESS> <I>Petra Nass</I> <BR><I>1999-06-15</I> </ADDRESS> </BODY> </HTML>
