Changeset 226 for trunk/python/scantable.py

Timestamp:

01/19/05 18:24:16 (19 years ago)

Author:

mar637

Message:

added rcParams to support rc style default parameters, read from .asaprc

File:

• trunk/python/scantable.py (modified) (23 diffs)

trunk/python/scantable.py

@@ -1 +1 @@
 from asap._asap import sdtable
+from asap import rcParams
 from numarray import ones,zeros
 import sys
@@ -21 +22 @@
                          scantable
         """
-        self._vb = True
+        self._vb = rcParams['verbose']
         self._p = None
         from os import stat as st
@@ -28 +29 @@
             sdtable.__init__(self, filename)            
         else:
-            mode = st(filename)[stat.ST_MODE]
+            try:
+                mode = st(filename)[stat.ST_MODE]
+            except OSError:
+                print "File not found"
+                return
             if stat.S_ISDIR(mode):
                 # crude check if asap table
@@ -35 +40 @@
                 else:
                     print 'The given file is not a valid asap table'
-            else:            
+                    return
+            else:
+                autoav = rcParams['scantable.autoaverage']
+
                 from asap._asap import sdreader
                 r = sdreader(filename)
@@ -41 +49 @@
                 r.read([-1])
                 tbl = r.getdata()
-                
-                from asap._asap import average
-                tmp = tuple([tbl])
-                print 'Auto averging integrations...'
-                tbl2 = average(tmp,(),True,'none')
-                sdtable.__init__(self,tbl2)
-                del r,tbl
-
-    def save(self, name="", format='ASAP'):
+                if autoav:
+                    from asap._asap import average
+                    tmp = tuple([tbl])
+                    print 'Auto averaging integrations...'
+                    tbl2 = average(tmp,(),True,'none')
+                    sdtable.__init__(self,tbl2)
+                    del r,tbl
+                else:
+                    sdtable.__init__(self,tbl)
+
+    def save(self, name="", format=None):
         """
         Store the scantable on disk. This can be a asap file or SDFITS/MS2.
@@ -60 +70 @@
                                        'FITS' (saves each row as a FITS Image) 
                                        'ASCII' (saves as ascii text file)
-                                       'MS2' (saves as an aips++ MeasurementSet V2)
+                                       'MS2' (saves as an aips++
+                                              MeasurementSet V2)
         Example:
             scan.save('myscan.asap')
             scan.save('myscan.sdfits','SDFITS')
         """
+        if format is None: format = rcParams['scantable.save']
         if format == 'ASAP':
             self._save(name)
@@ -72 +84 @@
             w.write(self, name)
         return
-
-    def _verbose(self, *args):
-        """
-        Set the verbose level True or False, to indicate if output
-        should be printed as well as returned.
-        """
-        if len(args) == 0:
-            return self._vb
-        elif type(args[0]) is bool:
-            self._vb = args[0]
-            return
 
     def copy(self):
@@ -166 +167 @@
             print " Cursor selection"
             print "--------------------------------------------------"
-            out = 'Beam=%d IF=%d Pol=%d '% (i,j,k)
+            out = 'Beam=%d IF=%d Pol=%d ' % (i,j,k)
             print out
         return i,j,k
 
-    def stats(self, stat='stddev', mask=None, all=True):
+    def stats(self, stat='stddev', mask=None, all=None):
         """
         Determine the specified statistic of the current beam/if/pol
@@ -180 +181 @@
             mask:    an optional mask specifying where the statistic
                      should be determined.
-            all:     optional flag to show all (default) or a selected
-                     spectrum of Beam/IF/Pol
+            all:     optional flag to show all (default or .asaprc) or a
+                     cursor selected spectrum of Beam/IF/Pol
 
         Example:
@@ -188 +189 @@
             scan.stats(stat='mean', mask=m)
         """
+        if all is None: all = rcParams['scantable.allaxes']
         from asap._asap import stats as _stats
         if mask == None:
@@ -210 +212 @@
                             out += '= %3.3f\n' % (statval[l])
                 out += "--------------------------------------------------\n"
-
             if self._vb:
                 print "--------------------------------------------------"
@@ -232 +233 @@
                 out += '= %3.3f\n' % (statval[l])
                 out +=  "--------------------------------------------------\n"
+
             if self._vb:
                 print "--------------------------------------------------"
@@ -239 +241 @@
             return statval
 
-    def stddev(self,mask=None, all=True):
+    def stddev(self,mask=None, all=None):
         """
         Determine the standard deviation of the current beam/if/pol
@@ -247 +249 @@
             mask:    an optional mask specifying where the standard
                      deviation should be determined.
-            all:     optional flag to show all or a selected
-                     spectrum of Beam/IF/Pol
+            all:     optional flag to show all or a cursor selected
+                     spectrum of Beam/IF/Pol. Default is all or taken
+                     from .asaprc
 
         Example:
@@ -255 +258 @@
             scan.stddev(mask=m)
         """
+        if all is None: all = rcParams['scantable.allaxes']
         return self.stats(stat='stddev',mask=mask, all=all);
 
-    def get_tsys(self, all=True):
+    def get_tsys(self, all=None):
         """
         Return the System temperatures.
@@ -268 +272 @@
             a list of Tsys values.
         """
+        if all is None: all = rcParams['scantable.allaxes']
         if all:
             out = ''
@@ -287 +292 @@
                             out += '= %3.3f\n' % (ts[l])
                 out+= "--------------------------------------------------\n"
-
             if self._vb:
                 print "--------------------------------------------------"
@@ -308 +312 @@
                 out += '= %3.3f\n' % (ts[l])
                 out += "--------------------------------------------------\n"
-
             if self._vb:
                 print "--------------------------------------------------"
@@ -345 +348 @@
         if self._p: self.plot()
 
-    def set_freqframe(self, frame='LSRK'):
+    def set_freqframe(self, frame=None):
         """
         Set the frame type of the Spectral Axis.
@@ -353 +356 @@
             scan.set_freqframe('BARY')
         """
+        if not frame: frame = rcParams['scantable.freqframe']
         valid = ['REST','TOPO','LSRD','LSRK','BARY', \
                    'GEO','GALACTO','LGROUP','CMB']
@@ -379 +383 @@
         selected Beam/IF/Pol
         Parameters:
-            none
+            rowno:    an optional row number in the scantable. Default is the
+                      first row, i.e. rowno=0
         Returns:
             The abcissa values and it's format string.
@@ -501 +506 @@
             self._p.set_panels(rows=npan)
         xlab,ylab,tlab = None,None,None
-        vb = self._verbose()
-        self._verbose(False)
-        sel = self.get_selection()
+        self._vb = False
+        sel = self.get_selection()        
         for i in range(npan):
-            self._p.subplot(i)
+            if npan > 1:
+                self._p.subplot(i)
             for j in range(validcol[col]):
                 x = None
@@ -540 +545 @@
         self._p.release()
         self.set_selection(sel[0],sel[1],sel[2])
-        self._verbose(vb)
+        self._vb = rcParams['verbose']
         return