|
Last change
on this file since 1702 was 1313, checked in by mar637, 18 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.
