1 | from __future__ import absolute_import
|
---|
2 |
|
---|
3 | import collections
|
---|
4 | import numpy as np
|
---|
5 |
|
---|
6 |
|
---|
7 | class FTS(object):
|
---|
8 | """Class to describe the Fourier Transform Spectrometer.
|
---|
9 | """
|
---|
10 |
|
---|
11 | def __init__(self, parameters):
|
---|
12 | self.parameters = parameters
|
---|
13 |
|
---|
14 | print parameters['substages']['FTSpectrograph'].keys()
|
---|
15 |
|
---|
16 | # read params from FTSpectrograph sheet
|
---|
17 | ftsparm = parameters['substages']['FTSpectrograph']
|
---|
18 |
|
---|
19 | self.result = collections.OrderedDict()
|
---|
20 |
|
---|
21 | # slightly awkward getting values as they are keyed by
|
---|
22 | # the row number in the parent spreadsheet
|
---|
23 | row = ftsparm['wn min [cm-1]'].keys()[0]
|
---|
24 | self.result['wnmin'] = wnmin = ftsparm['wn min [cm-1]'][row]
|
---|
25 | self.result['wnmax'] = wnmax = ftsparm['wn max [cm-1]'][row]
|
---|
26 | self.result['specres'] = specres = ftsparm['Spectral Res'][row]
|
---|
27 |
|
---|
28 | # read params from FTSMechanical sheet
|
---|
29 | ftsparm = parameters['substages']['FTSMechanical']
|
---|
30 | row = ftsparm['V drive [cm/s]'].keys()[0]
|
---|
31 | self.result['vdrive'] = ftsparm['V drive [cm/s]'][row]
|
---|
32 |
|
---|
33 | # derived parameters
|
---|
34 | self.result['delta_wn'] = delta_wn = wnmin / specres
|
---|
35 | self.result['delta_opd'] = 1.0 / (2.0 * wnmax * 100.0)
|
---|
36 | self.result['ftsnsample'] = nsample = np.ceil(wnmax / delta_wn)
|
---|
37 | self.result['opd_max'] = 1.0 / (2.0 * delta_wn)
|
---|
38 |
|
---|
39 | # spectral points sampled by FTS (wn=wavenumber[cm-1]
|
---|
40 | # lambda=wavelength[m])
|
---|
41 | i = np.arange(int(nsample))
|
---|
42 | fts_wn = wnmax * i / float(nsample-1)
|
---|
43 | # crop the spectrum below wnmin. This does not seem to be done
|
---|
44 | # in the matlab version.
|
---|
45 | fts_wn = fts_wn[fts_wn > wnmin]
|
---|
46 | fts_lambda = 1.0 / (fts_wn * 100.0)
|
---|
47 |
|
---|
48 | self.result['fts_wn'] = fts_wn
|
---|
49 | self.result['fts_lambda'] = fts_lambda
|
---|
50 |
|
---|
51 | print 'more to add in FTSDrive'
|
---|
52 |
|
---|
53 | def run(self):
|
---|
54 | print 'FTS.run'
|
---|
55 | return self.result
|
---|
56 |
|
---|
57 | def __repr__(self):
|
---|
58 | return '''
|
---|
59 | FTS:
|
---|
60 | wnmin : {wnmin}
|
---|
61 | wnmax : {wnmax}
|
---|
62 | specres : {specres}
|
---|
63 | delta wn : {delta_wn}
|
---|
64 | delta opd: {delta_opd}
|
---|
65 | nsample : {nsample}
|
---|
66 | max opd : {opd_max}
|
---|
67 | '''.format(
|
---|
68 | wnmin=self.result['wnmin'],
|
---|
69 | wnmax=self.result['wnmax'],
|
---|
70 | specres=self.result['specres'],
|
---|
71 | delta_wn=self.result['delta_wn'],
|
---|
72 | delta_opd=self.result['delta_opd'],
|
---|
73 | nsample=self.result['nsample'],
|
---|
74 | opd_max=self.result['opd_max'])
|
---|
75 |
|
---|