<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>

    <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
  </head>
  <body bgcolor="#ffffff" text="#000000">
    Hello,<br>
    <br>
    This error comes from a problem with how the coordinate lists are
    created within 'tweakshifts' and only shows up when there is only 1
    identified source in a chip given the object finding parameters
    provided when that task was run.  This task will be replaced shortly
    (end of this summer?) with a new Python, WCSLIB-based task, so no
    effort will be going into fixing any bugs found in using this task. 
    There are a couple of parameters in 'tweakshift' that provide some
    control over the object finding algorithm in 'daofind' (and
    'sextractor', depending on which you use); namely, fwhmpsf, sigma,
    datamin, datamax, threshold and nsigma.  Varying these can usually
    result in more than a single object being detected, assuming more
    than 1 object exists in the field of view at all.   This will then
    allow 'tweakshifts' to at least get past this problem.  <br>
    <br>
    As a side note, there are many other ways of specifying the names of
    the input images other than building an association (ASN) table,
    most of which are easier to use and manage.  Please review the help
    provided for this task, something which is made easiest by using the
    EPAR GUI and clicking on the 'tweakshift help' button at the top. 
    For example, you can run 'tweakshifts' with '*flt.fits' as the
    input. The ASN table actually gets stored as a binary FITS table for
    cases where you want to keep track of the shifts computed in the
    same file, rather relying on than the more convenient ASCII shift
    file that gets written by 'tweakshifts'.  <br>
    <br>
    If you have any further questions or problems with this task, please
    feel free to email them '<a class="moz-txt-link-abbreviated" href="mailto:help@stsci.edu">help@stsci.edu</a>' and I will be glad to
    continue working with you to help you use this task as successfully
    as possible.<br>
    <br>
    Warren Hack<br>
    Science Software Branch<br>
    Space Telescope Science Institute<br>
    <br>
    <blockquote type="cite" style="color: rgb(0, 0, 0);">From: Steve
      White <a class="moz-txt-link-rfc2396E"
        href="mailto:stevan.white@googlemail.com"><stevan.white@googlemail.com></a>
      <br>
      Date: June 20, 2011 5:44:16 AM EDT
      <br>
      To: <a class="moz-txt-link-rfc2396E"
        href="mailto:astropy@scipy.org"><astropy@scipy.org></a>
      <br>
      Subject: [AstroPy] tweakshifts in writeCoordFile()--IndexError:
      invalid index to scalar variable
      <br>
      <br>
      Hi!
      <br>
      <br>
      I'm trying to create a shifts file along the lines of the
      <br>
          Multidrizzle Handbook 5.5.3.6 Tweakshifts
      <br>
          <a class="moz-txt-link-freetext"
href="http://www.stsci.edu/hst/HST_overview/documents/multidrizzle/ch56.html#251052">http://www.stsci.edu/hst/HST_overview/documents/multidrizzle/ch56.html#251052</a>
      <br>
      <br>
      Tweakshifts gets stuck in a strange place deep in the python code.
      <br>
      I can't tell if I'm missing some step or if there is a bug in the
      code.
      <br>
      <br>
      Iraf 2.14.1 is installed (but also tried 2.15.a1) with stsdas and
      <br>
      tables v 3.13 and pytools SVN revision 11593.
      <br>
      Here, for example, the input files are from HST proposal 6431.
      <br>
      <br>
      First I determined which files correspond to a single pointing,
      and feed that
      <br>
      list into
      <br>
          asnt = asnutil.ASNTable( files, output = out_prefix ),
      <br>
             asnt.create()
      <br>
             asnt.write()
      <br>
      to create asn table files.
      <br>
      <br>
      When I run tweakshifts on the ASN files thus produced,
      <br>
          --> tweakshifts u34l390250_asn.fits undistort=yes mode=h
      <br>
      it does some stuff then dies, as listed below.
      <br>
      <br>
      I also attach the ASN file (which seems to consist of mixed text
      and binary).
      <br>
      <br>
      Thanks!
      <br>
      ----------------------------------------------------------
      <br>
      --> tweakshifts u34l390250_asn.fits undistort=yes mode=h
      <br>
      xyxin:  INDEF    xyyin:  INDEF
      <br>
      fitgeometry:  rscale
      <br>
      Tweakshifts Version  0.7.3 (19 Aug 2010)
      <br>
      <br>
      + MAKEWCS Version 1.1.7 (6 Jul 2010)
      <br>
      -Updating image  u34l3901m_c0m.fits[sci,1]
      <br>
      - IDCTAB: Distortion model from row 25 for chip 1 : F555W and
      CLEAR
      <br>
      - OFFTAB: Offset interpolated from rows 9 and 13
      <br>
      - IDCTAB: Distortion model from row 165 for chip 3 : F555W and
      CLEAR
      <br>
      - OFFTAB: Offset interpolated from rows 11 and 15
      <br>
      ...
      <br>
      Setting up output name:  u34l390250_drz.fits
      <br>
      Starting PyDrizzle Version  6.3.7 (3-Jan-2011)  at  10:02:52
      (20/06/2011)
      <br>
      - IDCTAB: Distortion model from row 25 for chip 1 : F555W and
      CLEAR
      <br>
      - OFFTAB: Offset interpolated from rows 9 and 13
      <br>
      - IDCTAB: Distortion model from row 95 for chip 2 : F555W and
      CLEAR
      <br>
      ...
      <br>
      Drizzle parameters have been calculated. Ready to .run()...
      <br>
      Finished calculating parameters at  10:02:55 (20/06/2011)
      <br>
      <br>
      <blockquote type="cite" style="color: rgb(0, 0, 0);">WARNING:
        Pixel stack overflow at position 493,431
        <br>
      </blockquote>
      <br>
         Computing undistorted target positions for  322 objects.
      <br>
      ...
      <br>
         Computing undistorted target positions for  1 objects.
      <br>
      Traceback (innermost last):
      <br>
       File "<CL script CL1>", line 1, in <module>
      <br>
       File "<i class="moz-txt-slash"><span class="moz-txt-tag">/</span>opt/iraf/extern/stsdas/pkg/analysis/dither<span
          class="moz-txt-tag">/</span></i>tweakshifts_iraf.py",
      <br>
      line 217, in tweak_iraf
      <br>
         fitbox        = _fitbox)
      <br>
       File "/opt/iraf/extern/stsdas/python/tweak.py", line 625, in run
      <br>
         shift_dict = runCatalog(flist,refname,match_pars)
      <br>
       File "/opt/iraf/extern/stsdas/python/tweak.py", line 861, in
      runCatalog
      <br>
         writeCoordFile(cxylist, fluxes, coordname)
      <br>
       File "/opt/iraf/extern/stsdas/python/tweak.py", line 252, in
      writeCoordFile
      <br>
         cstr = str(coord[0])+'    '+str(coord[1])+'    '+str(flux)
      +'\n'
      <br>
      IndexError: invalid index to scalar variable.
      <br>
    </blockquote>
    <blockquote type="cite" style="color: rgb(0, 0, 0);">_______________________________________________
      <br>
      AstroPy mailing list
      <br>
      <a class="moz-txt-link-abbreviated"
        href="mailto:AstroPy@scipy.org">AstroPy@scipy.org</a>
      <br>
      <a class="moz-txt-link-freetext"
        href="http://mail.scipy.org/mailman/listinfo/astropy">http://mail.scipy.org/mailman/listinfo/astropy</a>
      <br>
    </blockquote>
  </body>
</html>