1 | It appears that asap has some way of marking scans as bad not
|
---|
2 | and including them in the summary. However, the scans are in the
|
---|
3 | scantable. From a 50 row SDFITS bintable::
|
---|
4 |
|
---|
5 | In [1]: len(s_table.get_sourcename())
|
---|
6 | Out[1]: 50
|
---|
7 |
|
---|
8 | In [2]: np.array(s_table.get_sourcename())
|
---|
9 | Out[2]:
|
---|
10 | array(['n159 D', 'n159 D', 'G309.92', 'N159 C', 'N159 C', 'N159 C',
|
---|
11 | 'N159 C', 'N159 C', 'N159 C', 'N159 A', '0637-752', '0637-752',
|
---|
12 | '0637-752', '0637-752', 'G309.92', 'G309.92', 'G309.92', 'G309.92',
|
---|
13 | '1226+023', 'G339.88', 'G339.88', 'G339.88', 'G339.88', 'G188.95',
|
---|
14 | 'G188.95', 'G188.95', 'G323.74', 'G323.74', 'SGR B2', 'SGR B2',
|
---|
15 | 'SGR B2', 'SGR B2', 'G1.6', 'G1.6', 'G1.6', 'G1.6', 'SGR AB',
|
---|
16 | 'SGR AB', 'SGR A 2E0N', 'Sgr A meth pk1', 'Sgr A meth pk1', '30dor',
|
---|
17 | 'lmc-co-1', 'lmc-co-1', 'lmc-co-1', 'lmc-co-1', 'lmc-co-1',
|
---|
18 | 'lmc-co-1', 'lmc-co-1', 'lmc-co-1'],
|
---|
19 | dtype='|S14')
|
---|
20 |
|
---|
21 | In [3]: s_table.data.field('OBJECT')
|
---|
22 | Out[3]:
|
---|
23 | chararray(['n159 D', 'n159 D', 'G309.92', 'N159 C', 'N159 C', 'N159 C',
|
---|
24 | 'N159 C', 'N159 C', 'N159 C', 'N159 A', '0637-752', '0637-752',
|
---|
25 | '0637-752', '0637-752', 'G309.92', 'G309.92', 'G309.92', 'G309.92',
|
---|
26 | '1226+023', 'G339.88', 'G339.88', 'G339.88', 'G339.88', 'G188.95',
|
---|
27 | 'G188.95', 'G188.95', 'G323.74', 'G323.74', 'SGR B2', 'SGR B2',
|
---|
28 | 'SGR B2', 'SGR B2', 'G1.6', 'G1.6', 'G1.6', 'G1.6', 'SGR AB',
|
---|
29 | 'SGR AB', 'SGR A 2E0N', 'Sgr A meth pk1', 'Sgr A meth pk1', '30dor',
|
---|
30 | 'lmc-co-1', 'lmc-co-1', 'lmc-co-1', 'lmc-co-1', 'lmc-co-1',
|
---|
31 | 'lmc-co-1', 'lmc-co-1', 'lmc-co-1'],
|
---|
32 | dtype='|S16')
|
---|
33 |
|
---|
34 | So the scantable.get_sourcename() agrees with the SDFITS extension.
|
---|
35 | Likewise, the other gets give the right number of values::
|
---|
36 |
|
---|
37 | In [7]: len(s_table.get_azimuth())
|
---|
38 | Out[7]: 50
|
---|
39 |
|
---|
40 | In [8]: len(s_table.get_direction())
|
---|
41 | Out[8]: 50
|
---|
42 |
|
---|
43 | In [9]: len(s_table.get_inttime())
|
---|
44 | Out[9]: 50
|
---|
45 |
|
---|
46 | In [10]: len(s_table.get_time())
|
---|
47 | Out[10]: 50
|
---|
48 |
|
---|
49 | In [11]: len(s_table.get_tsys())
|
---|
50 | Out[11]: 50
|
---|
51 |
|
---|
52 | In [12]: len(s_table.getscannos())
|
---|
53 | Out[12]: 45
|
---|
54 |
|
---|
55 | In [16]: np.array(s_table.getscannos())
|
---|
56 | Out[16]:
|
---|
57 | array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
|
---|
58 | 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
|
---|
59 | 20, 21, 22, 23, 24, 25, 26, 27, 28, 29,
|
---|
60 | 30, 31, 32, 33, 34, 39,
|
---|
61 | 43, 44, 45, 46, 47, 48, 49,
|
---|
62 | 50, 51])
|
---|
63 |
|
---|
64 | In [18]: s35 = s_table.get_scan([35])
|
---|
65 | Couldn't find any match.
|
---|
66 |
|
---|
67 | In [19]: s51 = s_table.get_scan([51])
|
---|
68 |
|
---|
69 | In [20]:
|
---|
70 |
|
---|
71 | These absent scans cannot be extracted and are missing from
|
---|
72 | s_table.summary(). Also, there are fanthom scan numbers in
|
---|
73 | the sense that 50 scans counted from 0 should end at 49.
|
---|
74 | So we can't use the above scan numbers to index into the arrays
|
---|
75 | of the scantable.
|
---|
76 |
|
---|
77 | Here's the clue::
|
---|
78 |
|
---|
79 | In [22]: for i in s_table.getscannos():
|
---|
80 | ....: s = s_table.get_scan([i])
|
---|
81 | ....: print i, s.get_sourcename()
|
---|
82 | ....:
|
---|
83 | 0 ['n159 D']
|
---|
84 | 1 ['n159 D']
|
---|
85 | 2 ['G309.92']
|
---|
86 | 3 ['N159 C']
|
---|
87 | 4 ['N159 C']
|
---|
88 | 5 ['N159 C']
|
---|
89 | 6 ['N159 C']
|
---|
90 | 7 ['N159 C']
|
---|
91 | 8 ['N159 C']
|
---|
92 | 9 ['N159 A']
|
---|
93 | 10 ['0637-752', '0637-752', 'G309.92']
|
---|
94 | 11 ['0637-752', 'G309.92', 'G309.92']
|
---|
95 | 12 ['0637-752', 'G309.92']
|
---|
96 | 13 ['1226+023']
|
---|
97 | 14 ['G339.88']
|
---|
98 | 15 ['G339.88']
|
---|
99 | 16 ['G339.88']
|
---|
100 | 17 ['G339.88']
|
---|
101 | 18 ['G188.95']
|
---|
102 | 19 ['G188.95']
|
---|
103 | 20 ['G188.95']
|
---|
104 | 21 ['G323.74']
|
---|
105 | 22 ['G323.74']
|
---|
106 | 23 ['SGR B2']
|
---|
107 | 24 ['SGR B2']
|
---|
108 | 25 ['SGR B2']
|
---|
109 | 26 ['SGR B2']
|
---|
110 | 27 ['G1.6']
|
---|
111 | 28 ['G1.6']
|
---|
112 | 29 ['G1.6']
|
---|
113 | 30 ['G1.6']
|
---|
114 | 31 ['SGR AB']
|
---|
115 | 32 ['SGR AB']
|
---|
116 | 33 ['SGR A 2E0N']
|
---|
117 | 34 ['Sgr A meth pk1']
|
---|
118 | 39 ['Sgr A meth pk1']
|
---|
119 | 43 ['30dor']
|
---|
120 | 44 ['lmc-co-1']
|
---|
121 | 45 ['lmc-co-1']
|
---|
122 | 46 ['lmc-co-1']
|
---|
123 | 47 ['lmc-co-1']
|
---|
124 | 48 ['lmc-co-1']
|
---|
125 | 49 ['lmc-co-1']
|
---|
126 | 50 ['lmc-co-1']
|
---|
127 | 51 ['lmc-co-1']
|
---|
128 |
|
---|
129 | Aha! Fetching scan by number for 10, 11 and 12 returns scantables
|
---|
130 | with more than one row. Eight scans in three scantables reduces
|
---|
131 | the number of rows by five. Why should some rows be combined into
|
---|
132 | mini-scantables?
|
---|
133 |
|
---|
134 | In [23]: s10 = s_table.get_scan([10])
|
---|
135 |
|
---|
136 | In [25]: s10.get_direction()
|
---|
137 | Out[25]: ['06:37:23.5 -75.13.36.8', '00:00:00.0 +00.00.00.0', '06:37:23.5 -61.20.17.2']
|
---|
138 |
|
---|
139 | In [26]: s11 = s_table.get_scan([11])
|
---|
140 |
|
---|
141 | In [27]: s11.get_direction()
|
---|
142 | Out[27]: ['06:37:23.5 -75.13.36.8', '13:47:13.6 -61.20.17.2', '06:37:23.5 -61.20.17.2']
|
---|
143 |
|
---|
144 | In [28]: s12 = s_table.get_scan([12])
|
---|
145 |
|
---|
146 | In [29]: s12.get_direction()
|
---|
147 | Out[29]: ['06:37:23.5 -75.13.36.8', '06:37:23.5 -61.20.17.2']
|
---|