title: Python and 1D spectrum - How to access data within fits file
body: <p>I'm very new to Python and have been looking around for an answer I want but haven't found one. I have a 1D spectrum file (.fits) and would like to access its wavelength and flux information. For example, if I plug in a wavelength value I'm trying to write a program that would give me flux at the input wavelength. Could anyone give me a quick guide on how to write such a program, or what modules do I need to begin with? Will numpy and pyfits be enough to do the task? Thank you so much.</p>