Changeset 1995


Ignore:
Timestamp:
02/16/11 21:06:10 (13 years ago)
Author:
Kana Sugimoto
Message:

New Development: No

JIRA Issue: Yes (CAS-2775)

Ready for Test: Yes

Interface Changes: No

What Interface Changed:

Test Programs: plot spectra with FLAGROW=1

Put in Release Notes: No

Module(s): asapplotter, sdplot

Description:

Proper handling of FLAGROW values in sdplot/asapplotter.
Also, proper handling of plotrange for totally flagged spectra.


File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/python/asapplotter.py

    r1989 r1995  
    8585            #from asap.casatoolbar import CustomFlagToolbarTkAgg
    8686            #return CustomFlagToolbarTkAgg(self)
    87         else: return None
     87        return None
    8888
    8989    @asaplog_post_dec
     
    949949                else:
    950950                    y = scan._getspectrum(r)
    951                 m = scan._getmask(r)
    952                 from numpy import logical_not, logical_and
    953                 if self._maskselection and len(self._usermask) == len(m):
    954                     if d[self._stacking](r) in self._maskselection[self._stacking]:
    955                         m = logical_and(m, self._usermask)
     951                # flag application
     952                mr = scan._getflagrow(r)
    956953                from numpy import ma, array
     954                if mr:
     955                    y = ma.masked_array(y,mask=mr)
     956                else:
     957                    m = scan._getmask(r)
     958                    from numpy import logical_not, logical_and
     959                    if self._maskselection and len(self._usermask) == len(m):
     960                        if d[self._stacking](r) in self._maskselection[self._stacking]:
     961                            m = logical_and(m, self._usermask)
     962                    y = ma.masked_array(y,mask=logical_not(array(m,copy=False)))
     963
    957964                x = array(scan._getabcissa(r))
    958965                if self._offset:
    959966                    x += self._offset
    960                 y = ma.masked_array(y,mask=logical_not(array(m,copy=False)))
    961967                if self._minmaxx is not None:
    962968                    s,e = self._slice_indeces(x)
     
    978984                plotit = self._plotter.plot
    979985                if self._hist: plotit = self._plotter.hist
    980                 if len(x) > 0:
     986                if len(x) > 0 and not mr:
    981987                    plotit(x,y)
    982988                    xlim= self._minmaxx or [min(x),max(x)]
Note: See TracChangeset for help on using the changeset viewer.