Changeset 1655 for branches/alma/python


Ignore:
Timestamp:
11/05/09 21:40:21 (15 years ago)
Author:
WataruKawasaki
Message:

New Development: Yes

JIRA Issue: Yes (CAS-1433)

Ready to Release: Yes

Interface Changes: Yes

What Interface Changed: A new functionality flag_row() added

Test Programs:

Put in Release Notes: No

Module(s): sdflag

Description: Added a new functionality flag_row() to

enable row_based flagging/unflagging.


File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/alma/python/scantable.py

    r1653 r1655  
    825825        self._add_history("flag", varlist)
    826826
     827    def flag_row(self, rows=[], unflag=False):
     828        """
     829        Flag the selected data in row-based manner.
     830        Parameters:
     831            rows:   list of row numbers to be flagged. Default is no row (must be explicitly specified to execute row-based flagging).
     832            unflag: if True, unflag the data.
     833        """
     834        varlist = vars()
     835        try:
     836            self._flag_row(rows, unflag)
     837        except RuntimeError, msg:
     838            if rcParams['verbose']:
     839                print_log()
     840                asaplog.push(msg.message)
     841                print_log('ERROR')
     842                return
     843            else: raise
     844        self._add_history("flag_row", varlist)
     845       
     846       
    827847    def lag_flag(self, frequency, width=0.0, unit="GHz", insitu=None):
    828848        """
     
    15781598                f.set_function(poly=order)
    15791599
    1580             rows = range(self.nrow())
     1600            rows = range(workscan.nrow())
    15811601            if len(rows) > 0:
    15821602                self.blpars = []
     
    15841604            for r in rows:
    15851605                # take into account flagtra info (CAS-1434)
    1586                 flagtra = self._getmask(r)
     1606                flagtra = workscan._getmask(r)
    15871607                actualmask = mask[:]
    15881608                if len(actualmask) == 0:
     
    15941614                        for i in range(0, len(actualmask)):
    15951615                            actualmask[i] = actualmask[i] and flagtra[i]
    1596                 f.set_scan(self, actualmask)
    1597                 f.x = self._getabcissa(r)
    1598                 f.y = self._getspectrum(r)
     1616                f.set_scan(workscan, actualmask)
     1617                f.x = workscan._getabcissa(r)
     1618                f.y = workscan._getspectrum(r)
    15991619                f.data = None
    16001620                f.fit()
    1601                 fpar = f.get_parameters()
    16021621                if plot:
    16031622                    f.plot(residual=True)
     
    16071626                        continue
    16081627                workscan._setspectrum(f.fitter.getresidual(), r)
    1609                 self.blpars.append(fpar)
     1628                self.blpars.append(f.get_parameters())
     1629
    16101630            if plot:
    16111631                f._p.unmap()
    16121632                f._p = None
    1613             #f.set_scan(self, mask)
    1614             #s = f.auto_fit(insitu, plot=plot)
    1615             ## Save parameters of baseline fits as a class attribute.
    1616             ## NOTICE: It does not reflect changes in scantable!
    1617             #self.blpars = f.blpars
    16181633            workscan._add_history("poly_baseline", varlist)
    16191634            print_log()
     
    17591774            asaplog.push(msg, False)
    17601775
    1761             fpar = f.get_parameters()
    17621776            if plot:
    17631777                f.plot(residual=True)
     
    17671781                    self.masklists.append(None)
    17681782                    continue
     1783
    17691784            workscan._setspectrum(f.fitter.getresidual(), r)
    1770             self.blpars.append(fpar)
     1785            self.blpars.append(f.get_parameters())
    17711786            self.masklists.append(masklist)
    17721787        if plot:
Note: See TracChangeset for help on using the changeset viewer.