source: trunk/python/asaplotbase.py@ 1554

Last change on this file since 1554 was 1553, checked in by Malte Marquarding, 16 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.