Operator Instructions For Ronchical

When an operator is required to process a set of ronchi images they will usually be given a tape by one of the operations people and possibly a note, or verbal information regarding what day the ronchi was recorded and the approximate position of the ronchi data on the tape. The basic steps are then:
  1. Read the tape.
  2. unrice the images.
  3. Use rfits to rename the images and convert them to reals:
        fits_file = "*.fits"        FITS data source
        file_list = ""              File/extensions list
        iraf_file = ""              IRAF filename
      (make_image = yes)            Create an IRAF image?
     (long_header = no)             Print FITS header cards?
    (short_header = yes)            Print short header?
        (datatype = "real")         IRAF data type
           (blank = 0.)             Blank value
           (scale = yes)            Scale the data?
     (oldirafname = yes)            Use old IRAF name in place of iraf_file?
          (offset = 0)              Tape file offset
    This can be accomplished using the following command:
    rfits fits*fits "" "" make_image+ data=real scale+ old+
  4. Identify the calibration and ronchi images.
    These can be found by checking the header keyword INSTMODE in each of the images.
    To process ronchi images you will need to identify the following images:

    Image Type INSTMODE
    east cal 11
    west cal 12
    diffuser 21
    dark 4
    ronchi image 25
    ronchi guide 27

  5. Use imrename to rename the east cal image to east
  6. Use imrename to rename the west cal image to west
  7. Use imrename to rename the diffuser image to diffuser
  8. Use imcombine to combine the east and west calibration images and average the image planes:
    imcombine east,west avgcal proj+
  9. Use imcombine to average the planes of the dark image:
    imcombine dark avgdrk proj+
  10. Identify the ronchi images via:
    hselect *.fits $I,instmode
    This will generate a list containing the image name followed by the value of INSTMODE. Ronchi data sets will have INSTMODE values of 25 or 27. The type 27 images can be deleted (they are images recorded without the ronchi ruling when guiding was re-enabled in order to recenter the solar image.)
  11. Make a list of the input ronchi images to process.
  12. Set the task parameters correctly:
           inlist = "@inlist"       List of input images
             dark = "avgdrk"        Dark image name
              cal = "avgcal         Calibration image name =(east+west)/2.
         diffuser = "diff"          Diffuser image name
          (lc_bn1 = 196)            Last column of band 1
          (lc_bn2 = 452)            Last column of band 2
          (lc_bn3 = 708)            Last column of band 3
           (clmp1 = 11)             First row for clamp correction
           (clmp2 = 40)             Last row for clamp correction
          (vscale = 4.8000000000000E-5) Scale factor for vertical correction
          (hscale = 2.2000000000000E-5) Scale factor for horizontal correction)
        (function = "chebyshev")    Function to fit (legendre|chebyshev)
          (xorder = 5)              Order of function in x for spatial fit
          (yorder = 5)              Order of function in y for spatial fit
     (cross_terms = yes)            Include cross-terms?

  13. Run ronchical: rronchical @inlist
    where inlist is an ASCII file containing the names of the raw ronchi images.
  14. ronchical will produce images named according to: ssrciyymmdd, where ss is the 2-character site abbreviation and yy, mm, dd are the 2-digit year, month, and day of when the data were recorded.

Last updated by Cliff Toner on Tuesday, October 01, 2005, 14:30