1 | from scantable import scantable |
---|
2 | def average_scan(scan): |
---|
3 | """ |
---|
4 | Return a (time) averaged a scan, i.e. all correlator cycles |
---|
5 | are averaged into one "scan". |
---|
6 | """ |
---|
7 | from asap._asap import average as _av |
---|
8 | return scantable(_av(scan)) |
---|
9 | |
---|
10 | |
---|
11 | def average_scans(*args, **kwargs): |
---|
12 | """ |
---|
13 | Return the (time) average of a list of scans. [in channels only] |
---|
14 | Parameters: |
---|
15 | a comma separated list of scans |
---|
16 | mask: an optional mask |
---|
17 | Example: |
---|
18 | scanav = average_scans(scana,scanb) |
---|
19 | # return a time averaged scan from scan and scanb |
---|
20 | # without using a mask |
---|
21 | """ |
---|
22 | |
---|
23 | if len(args) < 2: |
---|
24 | print "Please give at least two scantables" |
---|
25 | return |
---|
26 | |
---|
27 | from asap._asap import averages as _av |
---|
28 | d = [args[0].nbeam(),args[0].nif(),args[0].npol(),args[0].nchan()] |
---|
29 | for s in args: |
---|
30 | if not isinstance(s,scantable): |
---|
31 | print "Please give a list of scantables" |
---|
32 | return |
---|
33 | dim = [s.nbeam(),s.nif(),s.npol(),s.nchan()] |
---|
34 | if (dim != d): |
---|
35 | print "All scans have to have the same numer of Beams/IFs/Pols/Chans" |
---|
36 | return |
---|
37 | if kwargs.has_key('mask'): |
---|
38 | return scantable(_av(args, kwargs.get('mask'))) |
---|
39 | else: |
---|
40 | from numarray import ones |
---|
41 | mask = list(ones(d[3])) |
---|
42 | return scantable(_av((args), mask)) |
---|
43 | |
---|
44 | def quotient(source, reference): |
---|
45 | """ |
---|
46 | Return the quotient of a 'source' scan and a 'reference' scan |
---|
47 | Parameters: |
---|
48 | source: the 'on' scan |
---|
49 | reference: the 'off' scan |
---|
50 | """ |
---|
51 | from asap._asap import quotient as _quot |
---|
52 | return scantable(_quot(source, reference)) |
---|
53 | |
---|
54 | def scale(scan, factor): |
---|
55 | """ |
---|
56 | Return a scan where all spectra are scaled by the give 'factor' |
---|
57 | Parameters: |
---|
58 | scan: a scantable |
---|
59 | factor: the sclaing factor |
---|
60 | Note: |
---|
61 | This currently applies the all beams/IFs/pols |
---|
62 | """ |
---|
63 | from asap._asap import scale as _scale |
---|
64 | return scantable(_scale(scan, factor)) |
---|
65 | |
---|
66 | |
---|
67 | def bin(scan, binwidth=5): |
---|
68 | """ |
---|
69 | """ |
---|
70 | from asap._asap import bin as _bin |
---|
71 | return scantable(_bin(scan, binwidth)) |
---|