Changeset 2269 for trunk/python


Ignore:
Timestamp:
08/11/11 12:19:57 (13 years ago)
Author:
Malte Marquarding
Message:

DO NOT change ORDER of arguments in functions. This WILL BREAK old code

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/python/scantable.py

    r2189 r2269  
    21732173
    21742174    @asaplog_post_dec
    2175     def smooth(self, kernel="hanning", width=5.0, order=2, plot=False, insitu=None):
     2175    def smooth(self, kernel="hanning", width=5.0, order=2, plot=False,
     2176               insitu=None):
    21762177        """\
    21772178        Smooth the spectrum by the specified kernel (conserving flux).
     
    22912292
    22922293    @asaplog_post_dec
    2293     def sinusoid_baseline(self, insitu=None, mask=None, applyfft=None, fftmethod=None, fftthresh=None,
    2294                           addwn=None, rejwn=None, clipthresh=None, clipniter=None, plot=None,
    2295                           getresidual=None, showprogress=None, minnrow=None, outlog=None, blfile=None):
     2294    def sinusoid_baseline(self, insitu=None, mask=None, applyfft=None,
     2295                          fftmethod=None, fftthresh=None,
     2296                          addwn=None, rejwn=None, clipthresh=None,
     2297                          clipniter=None, plot=None,
     2298                          getresidual=None, showprogress=None,
     2299                          minnrow=None, outlog=None, blfile=None):
    22962300        """\
    22972301        Return a scan which has been baselined (all rows) with sinusoidal functions.
     
    26822686
    26832687            #CURRENTLY, PLOT=true UNAVAILABLE UNTIL cubic spline fitting is implemented as a fitter method.
    2684             workscan._auto_cspline_baseline(mask, npiece, clipthresh, clipniter, normalise_edge_param(edge), threshold, chan_avg_limit, getresidual, pack_progress_params(showprogress, minnrow), outlog, blfile)
     2688            workscan._auto_cspline_baseline(mask, npiece, clipthresh,
     2689                                            clipniter,
     2690                                            normalise_edge_param(edge),
     2691                                            threshold,
     2692                                            chan_avg_limit, getresidual,
     2693                                            pack_progress_params(showprogress,
     2694                                                                 minnrow),
     2695                                            outlog, blfile)
    26852696            workscan._add_history("auto_cspline_baseline", varlist)
    26862697           
     
    26942705
    26952706    @asaplog_post_dec
    2696     def poly_baseline(self, insitu=None, mask=None, order=None, plot=None, getresidual=None,
    2697                       showprogress=None, minnrow=None, outlog=None, blfile=None):
     2707    def poly_baseline(self, mask=None, order=None, insitu=None, plot=None,
     2708                      getresidual=None, showprogress=None, minnrow=None,
     2709                      outlog=None, blfile=None):
    26982710        """\
    26992711        Return a scan which has been baselined (all rows) by a polynomial.
     
    27282740            varlist = vars()
    27292741       
    2730             if insitu is None: insitu = rcParams["insitu"]
     2742            if insitu is None:
     2743                insitu = rcParams["insitu"]
    27312744            if insitu:
    27322745                workscan = self
     
    27342747                workscan = self.copy()
    27352748
    2736             if mask         is None: mask         = [True for i in xrange(workscan.nchan())]
     2749            if mask         is None: mask         = [True for i in \
     2750                                                       xrange(workscan.nchan())]
    27372751            if order        is None: order        = 0
    27382752            if plot         is None: plot         = False
     
    27442758
    27452759            if plot:
    2746                 outblfile = (blfile != "") and os.path.exists(os.path.expanduser(os.path.expandvars(blfile)))
    2747                 if outblfile: blf = open(blfile, "a")
     2760                outblfile = (blfile != "") and \
     2761                    os.path.exists(os.path.expanduser(os.path.expandvars(blfile)))
     2762                if outblfile:
     2763                    blf = open(blfile, "a")
    27482764               
    27492765                f = fitter()
     
    27692785                    masklist = workscan.get_masklist(f.mask, row=r, silent=True)
    27702786                    #workscan._append_blinfo(blpars, masklist, f.mask)
    2771                     workscan._setspectrum((f.fitter.getresidual() if getresidual else f.fitter.getfit()), r)
     2787                    workscan._setspectrum((f.fitter.getresidual()
     2788                                           if getresidual else f.fitter.getfit()), r)
    27722789                   
    27732790                    if outblfile:
    27742791                        rms = workscan.get_rms(f.mask, r)
    2775                         dataout = workscan.format_blparams_row(blpars["params"], blpars["fixed"], rms, str(masklist), r, True)
     2792                        dataout = \
     2793                            workscan.format_blparams_row(blpars["params"],
     2794                                                         blpars["fixed"],
     2795                                                         rms, str(masklist),
     2796                                                         r, True)
    27762797                        blf.write(dataout)
    27772798
     
    27812802                if outblfile: blf.close()
    27822803            else:
    2783                 workscan._poly_baseline(mask, order, getresidual, pack_progress_params(showprogress, minnrow), outlog, blfile)
     2804                workscan._poly_baseline(mask, order, getresidual,
     2805                                        pack_progress_params(showprogress,
     2806                                                             minnrow),
     2807                                        outlog, blfile)
    27842808           
    27852809            workscan._add_history("poly_baseline", varlist)
     
    27942818
    27952819    @asaplog_post_dec
    2796     def auto_poly_baseline(self, insitu=None, mask=None, order=None, edge=None, threshold=None, chan_avg_limit=None,
    2797                            plot=None, getresidual=None, showprogress=None, minnrow=None, outlog=None, blfile=None):
     2820    def auto_poly_baseline(self, mask=None, order=None, edge=None,
     2821                           threshold=None, chan_avg_limit=None,
     2822                           plot=None, insitu=None,
     2823                           getresidual=None, showprogress=None,
     2824                           minnrow=None, outlog=None, blfile=None):
    27982825        """\
    27992826        Return a scan which has been baselined (all rows) by a polynomial.
     
    28492876            varlist = vars()
    28502877
    2851             if insitu is None: insitu = rcParams['insitu']
     2878            if insitu is None:
     2879                insitu = rcParams['insitu']
    28522880            if insitu:
    28532881                workscan = self
     
    28702898
    28712899            if plot:
    2872                 outblfile = (blfile != "") and os.path.exists(os.path.expanduser(os.path.expandvars(blfile)))
     2900                outblfile = (blfile != "") and \
     2901                    os.path.exists(os.path.expanduser(os.path.expandvars(blfile)))
    28732902                if outblfile: blf = open(blfile, "a")
    28742903               
    28752904                from asap.asaplinefind import linefinder
    28762905                fl = linefinder()
    2877                 fl.set_options(threshold=threshold,avg_limit=chan_avg_limit)
     2906                fl.set_options(threshold=threshold, avg_limit=chan_avg_limit)
    28782907                fl.set_scan(workscan)
    28792908               
     
    28862915                for r in rows:
    28872916                    idx = 2*workscan.getif(r)
    2888                     fl.find_lines(r, mask_and(mask, workscan._getmask(r)), edge[idx:idx+2])  # (CAS-1434)
     2917                    fl.find_lines(r, mask_and(mask, workscan._getmask(r)),
     2918                                  edge[idx:idx+2])  # (CAS-1434)
    28892919
    28902920                    f.x = workscan._getabcissa(r)
     
    29032933                    masklist = workscan.get_masklist(f.mask, row=r, silent=True)
    29042934                    #workscan._append_blinfo(blpars, masklist, f.mask)
    2905                     workscan._setspectrum((f.fitter.getresidual() if getresidual else f.fitter.getfit()), r)
     2935                    workscan._setspectrum((f.fitter.getresidual() \
     2936                                               if getresidual else f.fitter.getfit()), r)
    29062937
    29072938                    if outblfile:
    29082939                        rms = workscan.get_rms(f.mask, r)
    2909                         dataout = workscan.format_blparams_row(blpars["params"], blpars["fixed"], rms, str(masklist), r, True)
     2940                        dataout = \
     2941                            workscan.format_blparams_row(blpars["params"],
     2942                                                         blpars["fixed"],
     2943                                                         rms, str(masklist),
     2944                                                         r, True)
    29102945                        blf.write(dataout)
    29112946                   
     
    29152950                if outblfile: blf.close()
    29162951            else:
    2917                 workscan._auto_poly_baseline(mask, order, edge, threshold, chan_avg_limit, getresidual, pack_progress_params(showprogress, minnrow), outlog, blfile)
     2952                workscan._auto_poly_baseline(mask, order, edge, threshold,
     2953                                             chan_avg_limit, getresidual,
     2954                                             pack_progress_params(showprogress,
     2955                                                                  minnrow),
     2956                                             outlog, blfile)
    29182957
    29192958            workscan._add_history("auto_poly_baseline", varlist)
Note: See TracChangeset for help on using the changeset viewer.