source: trunk/python/asaplotbase.py @ 1553

Last change on this file since 1553 was 1553, checked in by Malte Marquarding, 15 years ago

Fix for Ticket #157; numpy api changed for mask/data access in MaskedArray?

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 27.9 KB
RevLine 
[705]1"""
2ASAP plotting class based on matplotlib.
3"""
4
5import sys
6from re import match
7
8import matplotlib
9
10from matplotlib.figure import Figure, Text
[1147]11from matplotlib.font_manager import FontProperties as FP
[705]12from matplotlib.numerix import sqrt
13from matplotlib import rc, rcParams
[710]14from asap import rcParams as asaprcParams
[1259]15from matplotlib.ticker import OldScalarFormatter
[1019]16from matplotlib.ticker import NullLocator
[1147]17
[1425]18# API change in mpl >= 0.98
19try:
20    from matplotlib.transforms import blended_transform_factory
21except ImportError:
22    from matplotlib.transforms import blend_xy_sep_transform  as blended_transform_factory
23
[1095]24if int(matplotlib.__version__.split(".")[1]) < 87:
25    print "Warning: matplotlib version < 0.87. This might cause errors. Please upgrade."
[1019]26
[1259]27#class MyFormatter(OldScalarFormatter):
28#    def __call__(self, x, pos=None):
29#        last = len(self.locs)-2
30#        if pos==0:
31#            return ''
32#        else: return OldScalarFormatter.__call__(self, x, pos)
[1019]33
[705]34class asaplotbase:
35    """
36    ASAP plotting base class based on matplotlib.
37    """
38
39    def __init__(self, rows=1, cols=0, title='', size=(8,6), buffering=False):
[1019]40        """
41        Create a new instance of the ASAPlot plotting class.
[705]42
[1019]43        If rows < 1 then a separate call to set_panels() is required to define
44        the panel layout; refer to the doctext for set_panels().
45        """
[705]46        self.is_dead = False
[1019]47        self.figure = Figure(figsize=size, facecolor='#ddddee')
[705]48        self.canvas = None
49
[1019]50        self.set_title(title)
51        self.subplots = []
52        if rows > 0:
53            self.set_panels(rows, cols)
[705]54
[710]55        # Set matplotlib default colour sequence.
56        self.colormap = "green red black cyan magenta orange blue purple yellow pink".split()
[1019]57
[710]58        c = asaprcParams['plotter.colours']
59        if isinstance(c,str) and len(c) > 0:
60            self.colormap = c.split()
61
62        self.lsalias = {"line":  [1,0],
63                        "dashdot": [4,2,1,2],
64                        "dashed" : [4,2,4,2],
65                        "dotted" : [1,2],
66                        "dashdotdot": [4,2,1,2,1,2],
67                        "dashdashdot": [4,2,4,2,1,2]
68                        }
69
70        styles = "line dashed dotted dashdot".split()
71        c = asaprcParams['plotter.linestyles']
72        if isinstance(c,str) and len(c) > 0:
73            styles = c.split()
74        s = []
75        for ls in styles:
76            if self.lsalias.has_key(ls):
77                s.append(self.lsalias.get(ls))
78            else:
79                s.append('-')
80        self.linestyles = s
81
[705]82        self.color = 0;
[710]83        self.linestyle = 0;
[1019]84        self.attributes = {}
85        self.loc = 0
[705]86
[1019]87        self.buffering = buffering
[705]88
89    def clear(self):
[1019]90        """
[1147]91        Delete all lines from the plot.  Line numbering will restart from 0.
[1019]92        """
[705]93
[1019]94        for i in range(len(self.lines)):
95           self.delete(i)
96        self.axes.clear()
97        self.color = 0
98        self.lines = []
[705]99
[710]100    def palette(self, color, colormap=None, linestyle=0, linestyles=None):
[705]101        if colormap:
[710]102            if isinstance(colormap,list):
103                self.colormap = colormap
104            elif isinstance(colormap,str):
105                self.colormap = colormap.split()
[705]106        if 0 <= color < len(self.colormap):
107            self.color = color
[710]108        if linestyles:
109            self.linestyles = []
110            if isinstance(linestyles,list):
111                styles = linestyles
112            elif isinstance(linestyles,str):
113                styles = linestyles.split()
114            for ls in styles:
115                if self.lsalias.has_key(ls):
116                    self.linestyles.append(self.lsalias.get(ls))
117                else:
118                    self.linestyles.append(self.lsalias.get('line'))
119        if 0 <= linestyle < len(self.linestyles):
120            self.linestyle = linestyle
[705]121
122    def delete(self, numbers=None):
[1019]123        """
124        Delete the 0-relative line number, default is to delete the last.
125        The remaining lines are NOT renumbered.
126        """
[705]127
[1019]128        if numbers is None: numbers = [len(self.lines)-1]
[705]129
[1019]130        if not hasattr(numbers, '__iter__'):
131            numbers = [numbers]
[705]132
[1019]133        for number in numbers:
134            if 0 <= number < len(self.lines):
135                if self.lines[number] is not None:
136                    for line in self.lines[number]:
137                        line.set_linestyle('None')
138                        self.lines[number] = None
139        self.show()
[705]140
141    def get_line(self):
[1019]142        """
143        Get the current default line attributes.
144        """
145        return self.attributes
[705]146
147
[1086]148    def hist(self, x=None, y=None, fmt=None, add=None):
[1019]149        """
150        Plot a histogram.  N.B. the x values refer to the start of the
151        histogram bin.
[705]152
[1019]153        fmt is the line style as in plot().
154        """
[1086]155        from matplotlib.numerix import array
156        from matplotlib.numerix.ma import MaskedArray
[1553]157       
[1019]158        if x is None:
159            if y is None: return
[1023]160            x = range(len(y))
[705]161
[1019]162        if len(x) != len(y):
163            return
164        l2 = 2*len(x)
[1023]165        x2 = range(l2)
[1086]166        y2 = range(12)
[1023]167        y2 = range(l2)
168        m2 = range(l2)
[1553]169        ymsk = None
170        ydat = None
171        if hasattr(y, "raw_mask"):
172            # numpy < 1.1
173            ymsk = y.raw_mask()
174            ydat = y.raw_data()
175        else:
176            ymsk = y.mask
177            ydat = y.data
178
[1023]179        for i in range(l2):
[1019]180            x2[i] = x[i/2]
[1086]181            m2[i] = ymsk[i/2]
[705]182
[1023]183        y2[0] = 0.0
[1019]184        for i in range(1,l2):
[1086]185            y2[i] = ydat[(i-1)/2]
[705]186
[1086]187        self.plot(x2, MaskedArray(y2,mask=m2,copy=0), fmt, add)
[705]188
189
190    def hold(self, hold=True):
[1019]191        """
192        Buffer graphics until subsequently released.
193        """
194        self.buffering = hold
[705]195
196
197    def legend(self, loc=None):
[1019]198        """
199        Add a legend to the plot.
[705]200
[1019]201        Any other value for loc else disables the legend:
202             1: upper right
203             2: upper left
204             3: lower left
205             4: lower right
206             5: right
207             6: center left
208             7: center right
209             8: lower center
210             9: upper center
211            10: center
[705]212
[1019]213        """
[1095]214        if isinstance(loc, int):
[1098]215            self.loc = None
216            if 0 <= loc <= 10: self.loc = loc
[1095]217        else:
218            self.loc = None
219        #self.show()
[705]220
221
[1086]222    def plot(self, x=None, y=None, fmt=None, add=None):
[1019]223        """
224        Plot the next line in the current frame using the current line
225        attributes.  The ASAPlot graphics window will be mapped and raised.
[705]226
[1019]227        The argument list works a bit like the matlab plot() function.
228        """
229        if x is None:
230            if y is None: return
231            x = range(len(y))
[705]232
[1019]233        elif y is None:
234            y = x
235            x = range(len(y))
[1086]236        if fmt is None:
237            line = self.axes.plot(x, y)
[1019]238        else:
[1086]239            line = self.axes.plot(x, y, fmt)
[705]240
[1019]241        # Add to an existing line?
[1086]242        i = None
[1019]243        if add is None or len(self.lines) < add < 0:
244            # Don't add.
245            self.lines.append(line)
246            i = len(self.lines) - 1
247        else:
248            if add == 0: add = len(self.lines)
249            i = add - 1
250            self.lines[i].extend(line)
[705]251
[1019]252        # Set/reset attributes for the line.
253        gotcolour = False
254        for k, v in self.attributes.iteritems():
255            if k == 'color': gotcolour = True
256            for segment in self.lines[i]:
257                getattr(segment, "set_%s"%k)(v)
[705]258
[1019]259        if not gotcolour and len(self.colormap):
260            for segment in self.lines[i]:
261                getattr(segment, "set_color")(self.colormap[self.color])
[710]262                if len(self.colormap)  == 1:
263                    getattr(segment, "set_dashes")(self.linestyles[self.linestyle])
[1086]264
[1019]265            self.color += 1
266            if self.color >= len(self.colormap):
267                self.color = 0
[705]268
[710]269            if len(self.colormap) == 1:
270                self.linestyle += 1
[1019]271            if self.linestyle >= len(self.linestyles):
272                self.linestyle = 0
[710]273
[1019]274        self.show()
[705]275
276
277    def position(self):
[1019]278        """
279        Use the mouse to get a position from a graph.
280        """
[705]281
[1019]282        def position_disable(event):
283            self.register('button_press', None)
284            print '%.4f, %.4f' % (event.xdata, event.ydata)
[705]285
[1019]286        print 'Press any mouse button...'
287        self.register('button_press', position_disable)
[705]288
289
[1546]290    def get_region(self):
291        pos = []
292        print "Please select the bottom/left point"
293        pos.append(self.figure.ginput(n=1, show_clicks=False)[0])
294        print "Please select the top/right point"
295        pos.append(self.figure.ginput(n=1, show_clicks=False)[0])
296        return pos
297
298    def get_point(self):
299        print "Please select the point"
[1553]300        pt = self.figure.ginput(n=1, show_clicks=False)
301        if pt:
302            return pt[0]
303        else:
304            return None
[1546]305
[705]306    def region(self):
[1019]307        """
308        Use the mouse to get a rectangular region from a plot.
[705]309
[1019]310        The return value is [x0, y0, x1, y1] in world coordinates.
311        """
[705]312
[1019]313        def region_start(event):
314            height = self.canvas.figure.bbox.height()
315            self.rect = {'fig': None, 'height': height,
316                         'x': event.x, 'y': height - event.y,
317                         'world': [event.xdata, event.ydata,
318                                   event.xdata, event.ydata]}
319            self.register('button_press', None)
320            self.register('motion_notify', region_draw)
321            self.register('button_release', region_disable)
[705]322
[1019]323        def region_draw(event):
324            self.canvas._tkcanvas.delete(self.rect['fig'])
325            self.rect['fig'] = self.canvas._tkcanvas.create_rectangle(
326                                self.rect['x'], self.rect['y'],
327                                event.x, self.rect['height'] - event.y)
[705]328
[1019]329        def region_disable(event):
330            self.register('motion_notify', None)
331            self.register('button_release', None)
[705]332
[1019]333            self.canvas._tkcanvas.delete(self.rect['fig'])
[705]334
[1019]335            self.rect['world'][2:4] = [event.xdata, event.ydata]
336            print '(%.2f, %.2f)  (%.2f, %.2f)' % (self.rect['world'][0],
337                self.rect['world'][1], self.rect['world'][2],
338                self.rect['world'][3])
[705]339
[1019]340        self.register('button_press', region_start)
[705]341
[1019]342        # This has to be modified to block and return the result (currently
343        # printed by region_disable) when that becomes possible in matplotlib.
[705]344
[1019]345        return [0.0, 0.0, 0.0, 0.0]
[705]346
347
348    def register(self, type=None, func=None):
[1019]349        """
350        Register, reregister, or deregister events of type 'button_press',
351        'button_release', or 'motion_notify'.
[705]352
[1019]353        The specified callback function should have the following signature:
[705]354
[1019]355            def func(event)
[705]356
[1019]357        where event is an MplEvent instance containing the following data:
[705]358
[1019]359            name                # Event name.
360            canvas              # FigureCanvas instance generating the event.
361            x      = None       # x position - pixels from left of canvas.
362            y      = None       # y position - pixels from bottom of canvas.
363            button = None       # Button pressed: None, 1, 2, 3.
364            key    = None       # Key pressed: None, chr(range(255)), shift,
365                                  win, or control
366            inaxes = None       # Axes instance if cursor within axes.
367            xdata  = None       # x world coordinate.
368            ydata  = None       # y world coordinate.
[705]369
[1019]370        For example:
[705]371
[1019]372            def mouse_move(event):
373                print event.xdata, event.ydata
[705]374
[1019]375            a = asaplot()
376            a.register('motion_notify', mouse_move)
[705]377
[1019]378        If func is None, the event is deregistered.
[705]379
[1019]380        Note that in TkAgg keyboard button presses don't generate an event.
381        """
[705]382
[1019]383        if not self.events.has_key(type): return
[705]384
[1019]385        if func is None:
386            if self.events[type] is not None:
387                # It's not clear that this does anything.
388                self.canvas.mpl_disconnect(self.events[type])
389                self.events[type] = None
[705]390
[1019]391                # It seems to be necessary to return events to the toolbar.
392                if type == 'motion_notify':
393                    self.canvas.mpl_connect(type + '_event',
394                        self.figmgr.toolbar.mouse_move)
395                elif type == 'button_press':
396                    self.canvas.mpl_connect(type + '_event',
397                        self.figmgr.toolbar.press)
398                elif type == 'button_release':
399                    self.canvas.mpl_connect(type + '_event',
400                        self.figmgr.toolbar.release)
[705]401
[1019]402        else:
403            self.events[type] = self.canvas.mpl_connect(type + '_event', func)
[705]404
405
406    def release(self):
[1019]407        """
408        Release buffered graphics.
409        """
410        self.buffering = False
411        self.show()
[705]412
413
[1095]414    def save(self, fname=None, orientation=None, dpi=None, papertype=None):
[1019]415        """
416        Save the plot to a file.
[705]417
[1019]418        fname is the name of the output file.  The image format is determined
419        from the file suffix; 'png', 'ps', and 'eps' are recognized.  If no
420        file name is specified 'yyyymmdd_hhmmss.png' is created in the current
421        directory.
422        """
[1095]423        from asap import rcParams
424        if papertype is None:
425            papertype = rcParams['plotter.papertype']
[1019]426        if fname is None:
427            from datetime import datetime
428            dstr = datetime.now().strftime('%Y%m%d_%H%M%S')
429            fname = 'asap'+dstr+'.png'
[705]430
[1019]431        d = ['png','.ps','eps']
[705]432
[1019]433        from os.path import expandvars
434        fname = expandvars(fname)
[705]435
[1019]436        if fname[-3:].lower() in d:
437            try:
[705]438                if fname[-3:].lower() == ".ps":
[1020]439                    from matplotlib import __version__ as mv
[1479]440                    w = self.figure.get_figwidth()
441                    h = self.figure.get_figheight()
[1019]442
[705]443                    if orientation is None:
[1147]444                        # oriented
[705]445                        if w > h:
446                            orientation = 'landscape'
447                        else:
448                            orientation = 'portrait'
[1095]449                    from matplotlib.backends.backend_ps import papersize
450                    pw,ph = papersize[papertype.lower()]
[1025]451                    ds = None
452                    if orientation == 'landscape':
[1095]453                        ds = min(ph/w, pw/h)
[1025]454                    else:
[1095]455                        ds = min(pw/w, ph/h)
[1025]456                    ow = ds * w
457                    oh = ds * h
[1479]458                    self.figure.set_size_inches((ow, oh))
[1095]459                    self.figure.savefig(fname, orientation=orientation,
460                                        papertype=papertype.lower())
[1479]461                    self.figure.set_size_inches((w, h))
[705]462                    print 'Written file %s' % (fname)
[1019]463                else:
[705]464                    if dpi is None:
465                        dpi =150
[1025]466                    self.figure.savefig(fname,dpi=dpi)
[705]467                    print 'Written file %s' % (fname)
[1019]468            except IOError, msg:
469                print 'Failed to save %s: Error msg was\n\n%s' % (fname, err)
470                return
471        else:
472            print "Invalid image type. Valid types are:"
473            print "'ps', 'eps', 'png'"
[705]474
475
476    def set_axes(self, what=None, *args, **kwargs):
[1019]477        """
478        Set attributes for the axes by calling the relevant Axes.set_*()
479        method.  Colour translation is done as described in the doctext
480        for palette().
481        """
[705]482
[1019]483        if what is None: return
484        if what[-6:] == 'colour': what = what[:-6] + 'color'
[705]485
[1153]486        key = "colour"
487        if kwargs.has_key(key):
488            val = kwargs.pop(key)
489            kwargs["color"] = val
[705]490
[1153]491        getattr(self.axes, "set_%s"%what)(*args, **kwargs)
[705]492
[1153]493        self.show(hardrefresh=False)
[705]494
[1019]495
[705]496    def set_figure(self, what=None, *args, **kwargs):
[1019]497        """
498        Set attributes for the figure by calling the relevant Figure.set_*()
499        method.  Colour translation is done as described in the doctext
500        for palette().
501        """
[705]502
[1019]503        if what is None: return
504        if what[-6:] == 'colour': what = what[:-6] + 'color'
505        #if what[-5:] == 'color' and len(args):
506        #    args = (get_colour(args[0]),)
[705]507
[1019]508        newargs = {}
509        for k, v in kwargs.iteritems():
510            k = k.lower()
511            if k == 'colour': k = 'color'
512            newargs[k] = v
[705]513
[1019]514        getattr(self.figure, "set_%s"%what)(*args, **newargs)
[1153]515        self.show(hardrefresh=False)
[705]516
517
518    def set_limits(self, xlim=None, ylim=None):
[1019]519        """
520        Set x-, and y-limits for each subplot.
[705]521
[1019]522        xlim = [xmin, xmax] as in axes.set_xlim().
523        ylim = [ymin, ymax] as in axes.set_ylim().
524        """
525        for s in self.subplots:
526            self.axes  = s['axes']
527            self.lines = s['lines']
[705]528            oldxlim =  list(self.axes.get_xlim())
529            oldylim =  list(self.axes.get_ylim())
530            if xlim is not None:
531                for i in range(len(xlim)):
532                    if xlim[i] is not None:
533                        oldxlim[i] = xlim[i]
[1019]534            if ylim is not None:
[705]535                for i in range(len(ylim)):
536                    if ylim[i] is not None:
537                        oldylim[i] = ylim[i]
538            self.axes.set_xlim(oldxlim)
539            self.axes.set_ylim(oldylim)
540        return
541
542
543    def set_line(self, number=None, **kwargs):
[1019]544        """
545        Set attributes for the specified line, or else the next line(s)
546        to be plotted.
[705]547
[1019]548        number is the 0-relative number of a line that has already been
549        plotted.  If no such line exists, attributes are recorded and used
550        for the next line(s) to be plotted.
[705]551
[1019]552        Keyword arguments specify Line2D attributes, e.g. color='r'.  Do
[705]553
[1019]554            import matplotlib
555            help(matplotlib.lines)
[705]556
[1019]557        The set_* methods of class Line2D define the attribute names and
558        values.  For non-US usage, "colour" is recognized as synonymous with
559        "color".
[705]560
[1019]561        Set the value to None to delete an attribute.
[705]562
[1019]563        Colour translation is done as described in the doctext for palette().
564        """
[705]565
[1019]566        redraw = False
567        for k, v in kwargs.iteritems():
568            k = k.lower()
569            if k == 'colour': k = 'color'
[705]570
[1019]571            if 0 <= number < len(self.lines):
572                if self.lines[number] is not None:
573                    for line in self.lines[number]:
574                        getattr(line, "set_%s"%k)(v)
575                    redraw = True
576            else:
577                if v is None:
578                    del self.attributes[k]
579                else:
580                    self.attributes[k] = v
[705]581
[1153]582        if redraw: self.show(hardrefresh=False)
[705]583
584
585    def set_panels(self, rows=1, cols=0, n=-1, nplots=-1, ganged=True):
[1019]586        """
587        Set the panel layout.
[705]588
[1019]589        rows and cols, if cols != 0, specify the number of rows and columns in
590        a regular layout.   (Indexing of these panels in matplotlib is row-
591        major, i.e. column varies fastest.)
[705]592
[1019]593        cols == 0 is interpreted as a retangular layout that accomodates
594        'rows' panels, e.g. rows == 6, cols == 0 is equivalent to
595        rows == 2, cols == 3.
[705]596
[1019]597        0 <= n < rows*cols is interpreted as the 0-relative panel number in
598        the configuration specified by rows and cols to be added to the
599        current figure as its next 0-relative panel number (i).  This allows
600        non-regular panel layouts to be constructed via multiple calls.  Any
601        other value of n clears the plot and produces a rectangular array of
602        empty panels.  The number of these may be limited by nplots.
603        """
604        if n < 0 and len(self.subplots):
605            self.figure.clear()
606            self.set_title()
[705]607
[1019]608        if rows < 1: rows = 1
[705]609
[1019]610        if cols <= 0:
611            i = int(sqrt(rows))
612            if i*i < rows: i += 1
613            cols = i
[705]614
[1019]615            if i*(i-1) >= rows: i -= 1
616            rows = i
[705]617
[1019]618        if 0 <= n < rows*cols:
619            i = len(self.subplots)
620            self.subplots.append({})
[705]621
[1019]622            self.subplots[i]['axes']  = self.figure.add_subplot(rows,
623                                            cols, n+1)
624            self.subplots[i]['lines'] = []
[705]625
[1019]626            if i == 0: self.subplot(0)
[705]627
[1019]628            self.rows = 0
629            self.cols = 0
[705]630
[1019]631        else:
632            self.subplots = []
[705]633
[1019]634            if nplots < 1 or rows*cols < nplots:
635                nplots = rows*cols
[1025]636            if ganged:
637                hsp,wsp = None,None
638                if rows > 1: hsp = 0.0001
639                if cols > 1: wsp = 0.0001
640                self.figure.subplots_adjust(wspace=wsp,hspace=hsp)
[1019]641            for i in range(nplots):
642                self.subplots.append({})
[1153]643                self.subplots[i]['lines'] = []
644                if not ganged:
645                    self.subplots[i]['axes'] = self.figure.add_subplot(rows,
[1019]646                                                cols, i+1)
[1513]647                    if asaprcParams['plotter.xaxisformatting'] == 'mpl':
648                        self.subplots[i]['axes'].xaxis.set_major_formatter(OldScalarFormatter())
[1153]649                else:
650                    if i == 0:
651                        self.subplots[i]['axes'] = self.figure.add_subplot(rows,
652                                                cols, i+1)
[1513]653                        if asaprcParams['plotter.xaxisformatting'] != 'mpl':
654                           
655                            self.subplots[i]['axes'].xaxis.set_major_formatter(OldScalarFormatter())
[1153]656                    else:
657                        self.subplots[i]['axes'] = self.figure.add_subplot(rows,
658                                                cols, i+1,
659                                                sharex=self.subplots[0]['axes'],
660                                                sharey=self.subplots[0]['axes'])
[1259]661
[705]662                    # Suppress tick labelling for interior subplots.
663                    if i <= (rows-1)*cols - 1:
664                        if i+cols < nplots:
665                            # Suppress x-labels for frames width
666                            # adjacent frames
[1153]667                            for tick in self.subplots[i]['axes'].xaxis.majorTicks:
668                                tick.label1On = False
[1019]669                            self.subplots[i]['axes'].xaxis.label.set_visible(False)
[705]670                    if i%cols:
671                        # Suppress y-labels for frames not in the left column.
672                        for tick in self.subplots[i]['axes'].yaxis.majorTicks:
673                            tick.label1On = False
674                        self.subplots[i]['axes'].yaxis.label.set_visible(False)
[1025]675                    # disable the first tick of [1:ncol-1] of the last row
[1153]676                    #if i+1 < nplots:
677                    #    self.subplots[i]['axes'].xaxis.majorTicks[0].label1On = False
[1019]678                self.rows = rows
679                self.cols = cols
680            self.subplot(0)
[705]681
[1153]682    def tidy(self):
683        # this needs to be exceuted after the first "refresh"
684        nplots = len(self.subplots)
685        if nplots == 1: return
686        for i in xrange(nplots):
687            ax = self.subplots[i]['axes']
688            if i%self.cols:
689                ax.xaxis.majorTicks[0].label1On = False
690            else:
691                if i != 0:
692                    ax.yaxis.majorTicks[-1].label1On = False
693
694
[705]695    def set_title(self, title=None):
[1019]696        """
697        Set the title of the plot window.  Use the previous title if title is
698        omitted.
699        """
700        if title is not None:
701            self.title = title
[705]702
[1019]703        self.figure.text(0.5, 0.95, self.title, horizontalalignment='center')
[705]704
705
[1153]706    def show(self, hardrefresh=True):
[1019]707        """
708        Show graphics dependent on the current buffering state.
709        """
[1153]710        if not hardrefresh: return
[1019]711        if not self.buffering:
712            if self.loc is not None:
[1086]713                for sp in self.subplots:
[1019]714                    lines  = []
715                    labels = []
716                    i = 0
[1086]717                    for line in sp['lines']:
[1019]718                        i += 1
719                        if line is not None:
720                            lines.append(line[0])
721                            lbl = line[0].get_label()
722                            if lbl == '':
723                                lbl = str(i)
724                            labels.append(lbl)
[705]725
[1019]726                    if len(lines):
[1147]727                        fp = FP(size=rcParams['legend.fontsize'])
728                        fsz = fp.get_size_in_points() - len(lines)
729                        fp.set_size(max(fsz,6))
[1086]730                        sp['axes'].legend(tuple(lines), tuple(labels),
[1147]731                                          self.loc, prop=fp)
[1019]732                    else:
[1086]733                        sp['axes'].legend((' '))
[705]734
[1086]735            from matplotlib.artist import setp
[1147]736            fp = FP(size=rcParams['xtick.labelsize'])
737            xts = fp.get_size_in_points()- (self.cols)/2
738            fp = FP(size=rcParams['ytick.labelsize'])
739            yts = fp.get_size_in_points() - (self.rows)/2
[1086]740            for sp in self.subplots:
741                ax = sp['axes']
742                s = ax.title.get_size()
743                tsize = s-(self.cols+self.rows)
744                ax.title.set_size(tsize)
[1147]745                fp = FP(size=rcParams['axes.labelsize'])
[1086]746                setp(ax.get_xticklabels(), fontsize=xts)
747                setp(ax.get_yticklabels(), fontsize=yts)
[1147]748                origx =  fp.get_size_in_points()
749                origy = origx
[1086]750                off = 0
751                if self.cols > 1: off = self.cols
752                xfsize = origx-off
753                ax.xaxis.label.set_size(xfsize)
754                off = 0
755                if self.rows > 1: off = self.rows
756                yfsize = origy-off
757                ax.yaxis.label.set_size(yfsize)
[705]758
759    def subplot(self, i=None, inc=None):
[1019]760        """
761        Set the subplot to the 0-relative panel number as defined by one or
762        more invokations of set_panels().
763        """
764        l = len(self.subplots)
765        if l:
766            if i is not None:
767                self.i = i
[705]768
[1019]769            if inc is not None:
770                self.i += inc
[705]771
[1019]772            self.i %= l
773            self.axes  = self.subplots[self.i]['axes']
774            self.lines = self.subplots[self.i]['lines']
[705]775
776    def text(self, *args, **kwargs):
[1019]777        """
778        Add text to the figure.
779        """
780        self.figure.text(*args, **kwargs)
781        self.show()
[1147]782
783    def vline_with_label(self, x, y, label,
784                         location='bottom', rotate=0.0, **kwargs):
785        """
786        Plot a vertical line with label.
787        It takes "world" values fo x and y.
788        """
789        ax = self.axes
790        # need this to suppress autoscaling during this function
791        self.axes.set_autoscale_on(False)
792        ymin = 0.0
793        ymax = 1.0
794        valign = 'center'
795        if location.lower() == 'top':
796            y = max(0.0, y)
797        elif location.lower() == 'bottom':
798            y = min(0.0, y)
799        lbloffset = 0.06
800        # a rough estimate for the bb of the text
801        if rotate > 0.0: lbloffset = 0.03*len(label)
802        peakoffset = 0.01
[1535]803        xy = None
804        xy0 = None
805        # matplotlib api change 0.98 is using transform now
806        if hasattr(ax.transData, "inverse_xy_tup"):
807            # get relative coords
808            xy0 = ax.transData.xy_tup((x,y))
809            xy = ax.transAxes.inverse_xy_tup(xy0)
810        else:
811            xy0 = ax.transData.transform((x,y))
812            # get relative coords
813            xy = ax.transAxes.inverted().transform(xy0)
[1147]814        if location.lower() == 'top':
815            ymax = 1.0-lbloffset
816            ymin = xy[1]+peakoffset
817            valign = 'bottom'
818            ylbl = ymax+0.01
819        elif location.lower() == 'bottom':
820            ymin = lbloffset
821            ymax = xy[1]-peakoffset
822            valign = 'top'
823            ylbl = ymin-0.01
[1425]824        trans = blended_transform_factory(ax.transData, ax.transAxes)
[1147]825        l = ax.axvline(x, ymin, ymax, color='black', **kwargs)
826        t = ax.text(x, ylbl ,label, verticalalignment=valign,
827                                    horizontalalignment='center',
828                    rotation=rotate,transform = trans)
829        self.axes.set_autoscale_on(True)
Note: See TracBrowser for help on using the repository browser.