|
Last change
on this file was
1313,
checked in by mar637, 17 years ago
|
|
scantable plugin to flag negative values in the bandpass (MOPS)
|
|
File size:
1.0 KB
|
| Line | |
|---|
| 1 | def _vanvleck_flag(scan, level=5.0): |
|---|
| 2 | """ |
|---|
| 3 | Find level crossings and flag all data before (and after it). |
|---|
| 4 | This rectifies the problem caused by the MOPS bandpass being negative. |
|---|
| 5 | Use levels greater > 0.0 to get a better result. |
|---|
| 6 | """ |
|---|
| 7 | sel = selector() |
|---|
| 8 | saved = scan.get_selection() |
|---|
| 9 | for i in xrange(scan.nrow()): |
|---|
| 10 | sel.reset() |
|---|
| 11 | sel.set_scans(scan.getscan(i)) |
|---|
| 12 | sel.set_beams(scan.getbeam(i)) |
|---|
| 13 | sel.set_ifs(scan.getif(i)) |
|---|
| 14 | sel.set_polarisations(scan.getpol(i)) |
|---|
| 15 | sel.set_cycles(scan.getcycle(i)) |
|---|
| 16 | scan.set_selection(sel) |
|---|
| 17 | spec = scan._getspectrum() |
|---|
| 18 | msk = mask_not(scan._getmask()) |
|---|
| 19 | n = len(msk) |
|---|
| 20 | modified = False |
|---|
| 21 | for j in xrange(n): |
|---|
| 22 | if j < n and spec[j] <= level: |
|---|
| 23 | msk[j] = 1 |
|---|
| 24 | modified = True |
|---|
| 25 | if modified: |
|---|
| 26 | scan.flag(msk) |
|---|
| 27 | scan.set_selection(saved) |
|---|
| 28 | |
|---|
| 29 | if not "vanvleck_flag" in dir(scantable): |
|---|
| 30 | scantable.vanvleck_flag = _vanvleck_flag |
|---|
| 31 | del _vanvleck_flag |
|---|
Note: See
TracBrowser
for help on using the repository browser.
