def test_trace_edge(self):
(header, data1, targs, ssl, msl, asl) = \
IO.readfits_all("/users/npk/desktop/c9/m110326_3242.fits")
data = data1
ssl = ssl[ssl["Slit_Number"] != ' ']
numslits = np.round(
np.array(ssl["Slit_length"], dtype=np.float) / 7.02)
for i in range(len(ssl)):
print ssl[i]["Target_Name"], numslits[i]
def test_trace_edge(self):
(header, data1, targs, ssl, msl, asl) = \
IO.readfits_all("/users/npk/desktop/c9/m110326_3242.fits")
data = data1
ssl = ssl[ssl["Slit_Number"] != ' ']
numslits = np.round(np.array(ssl["Slit_length"],
dtype=np.float) / 7.02)
for i in range(len(ssl)):
print ssl[i]["Target_Name"], numslits[i]
npk April 14th 2011
'''
import MOSFIRE
import time
from MOSFIRE import Fit, IO
import numpy as np, pylab as pl
reload(Fit)
reload(IO)
if __name__ == "__main__":
(header, data1, targs, ssl, msl, asl) = IO.readfits_all("/users/npk/desktop/c9/m110326_3242.fits")
data = data1
ssl = ssl[ssl["Slit_Number"] != ' ']
numslits = np.round(np.array(ssl["Slit_length"], dtype=np.float) / 7.02)
for i in range(len(ssl)):
print ssl[i]["Target_Name"], numslits[i]
Outputs:
xposs []: Array of x positions along the slit edge [pix]
yposs []: The fitted y positions of the "top" edge of the slit [pix]
import pylab as pl
import scipy as sp
from MOSFIRE import IO, Fit, Bspline
reload(IO)
reload(Fit)
reload(Bspline)
np.set_printoptions(precision=3)
# use the following file and id'd spectrum to guess
(header, data, targs, ssl, msl, asl) = IO.readfits_all(
"/users/npk/desktop/c9/m110319_1949.fits")
band = 'H'
pl.ion()
DRAW = True
def fit_spec(data, pos, alpha, sinbeta, gamma, delta, band):
global DRAW
ar_h_lines = np.array([1.465435, 1.474317, 1.505062, 1.517684, 1.530607,
1.533334, 1.540685, 1.590403, 1.599386, 1.618444, 1.644107,
1.674465, 1.744967, 1.791961])
Ar_K_lines = np.array([1.982291, 1.997118, 2.032256, 2.057443,
2.062186, 2.099184, 2.133871, 2.15409, 2.204558, 2.208321,
2.313952, 2.385154])
