[101] | 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)) |
---|