<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">In the past I've done what Kelle has suggested (and I hope AstroPy will eventually be able to parse all of this, even if I end up writing it myself), but just for reference:<div><br></div><div><a href="http://iraf.net/irafdocs/specwcs.php">http://iraf.net/irafdocs/specwcs.php</a></div><div><br></div><div>Provides a good overview on the different dispersion representations that you'll typically find in FITS headers.</div><div><br></div><div>Cheers,</div><div>Andy</div><div><br></div><div><div>On 14/08/2012, at 11:47 AM, Ian Crossfield wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">

  
    <meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type">
  
  <div bgcolor="#ffffff" text="#000000">
    It's not quite what you're looking for, but I have Python code to
    generate a wavelength array from the "ec" echelle database IRAF
    generates:<br>
    <a class="moz-txt-link-freetext" href="http://www.astro.ucla.edu/~ianc/python/nsdata.html#nsdata.dispeval">http://www.astro.ucla.edu/~ianc/python/nsdata.html#nsdata.dispeval</a><br>
    <br>
    Kelle's suggestion is probably the easiest, though likely not the
    fastest.<br>
    <br>
    -Ian<br>
    <br>
    <br>
    On 8/14/12 5:10 PM, Kevin Gullikson wrote:
    <blockquote cite="mid:CAEUX0Xc+cKcUFLuOPnOTk75x-b=BX010EvaWxgY4rSWyosf+nw@mail.gmail.com" type="cite">Hi,<br>
      <br>
      I have some spectra that were reduced in iraf, and I would like to
      read them in to python to do some fitting that I already have code
      for. However, I am having trouble reading in the wavelength
      calibration. The spectra are echelle, and the wavelengths are
      defined with this thing in the header that looks like this:<br>
      <br>
      "<br>
      WAT2_001= 'wtype=multispec spec1 = "1 32 2 10672.221264362
      0.086046403007761 20'<br>
      WAT2_002= '46 0. 6.65 26.40 1. 0. 1 5 1. 2046. 10761.5346117191
      87.998952480064'<br>
      WAT2_003= ' -1.32409633589482 -0.0165054046288388
      -0.00680394594704411" spec2 ='<br>
      WAT2_004= ' "2 33 2 10348.8582697 0.083445037904806 2046 0. 34.60
      58.07 1. 0. 1'<br>
      WAT2_005= ' 5 1. 2046. 10435.4699172677 85.3379030819273
      -1.2833101652519 -0.01'<br>
      WAT2_006= '53518242619121 -0.0057861452751027" spec3 = "3 34 2
      10043.881872825 '<br>
      ...<br>
      "<br>
      <br>
      The numbers give the conversion from pixel coordinates to
      wavelength coordinates.<br>
      I was wondering if there was a python function that could
      correctly read all that in, parse it, and give me a spectrum in
      flux vs wavelength? I have tried playing with pywcs but that seems
      only to work for images?<br>
      <br>
      Cheers,<br clear="all">
      <br>
      Kevin Gullikson
      <pre wrap=""><fieldset class="mimeAttachmentHeader"></fieldset>
_______________________________________________
AstroPy mailing list
<a class="moz-txt-link-abbreviated" href="mailto:AstroPy@scipy.org">AstroPy@scipy.org</a>
<a class="moz-txt-link-freetext" href="http://mail.scipy.org/mailman/listinfo/astropy">http://mail.scipy.org/mailman/listinfo/astropy</a>
</pre>
    </blockquote>
    <br>
    <br>
    <pre class="moz-signature" cols="72">-- 
Ian Crossfield
Max-Planck-Institüt für Astronomie
</pre>
  </div>

_______________________________________________<br>AstroPy mailing list<br><a href="mailto:AstroPy@scipy.org">AstroPy@scipy.org</a><br>http://mail.scipy.org/mailman/listinfo/astropy<br></blockquote></div><br></body></html>