Changeset 1192 for trunk/python

Timestamp:

08/28/06 16:33:58 (18 years ago)

Author:

mar637

Message:

added lag_flag - flagging of frequencies in fft space

Location:

trunk/python

Files:

• __init__.py (modified) (1 diff)
• scantable.py (modified) (1 diff)

trunk/python/__init__.py

@@ -409 +409 @@
             set_restfreqs   - set a list of rest frequencies
             flag            - flag selected channels in the data
+            lag_flag        - flag specified frequency in the data
             save            - save the scantable to disk as either 'ASAP',
                               'SDFITS' or 'ASCII'

trunk/python/scantable.py

@@ -639 +639 @@
             else: raise
         self._add_history("flag", varlist)
+
+    def lag_flag(self, frequency, unit="GHz", width=0, insitu=None):
+        """
+        Flag the data in 'lag' space by providing a frequency to remove.
+        Flagged data in the scantable gets set to 0.0 before the fft.
+        No taper is applied.
+        Parameters:
+            frequency:    the frequency to remove
+            unit:         the frequency unit ()default "GHz")
+            width:        the number of 'lag' channels to extent the frequency
+                          by, i.e. the interval
+                          [bandwidth/frequency-width, bandwidth/frequency+width]
+        """
+        if insitu is None: insitu = rcParams['insitu']
+        self._math._setinsitu(insitu)
+        varlist = vars()
+        base = { "GHz": 100000000., "MHz": 1000000., "kHz": 1000., "Hz": 1. }
+        if not base.has_key(unit):
+            raise ValueError("%s is not a valid unit." % unit)
+        try:
+            s = scantable(self._math._lag_flag(self, frequency*base[unit], width))
+        except RuntimeError, msg:
+            if rcParams['verbose']:
+                print msg
+                return
+            else: raise
+        s._add_history("lag_flag", varlist)
+        print_log()
+        if insitu:
+            self._assign(s)
+        else:
+            return s