def plot_tunnels(tunnels,horiz,fignum=1,keep_lim=True,clear_lines=True,**kwargs): '''plot current tunnels in specified figure (axes[0]) re-zoom to limits at time of call after plotting if keep_lim clear previous plotted lines in axes[0] if clear_lines ''' ax = pylab.figure(fignum).axes[0] prev_ylim = ax.get_ybound() prev_xlim = ax.get_xbound() ax.lines = [] x,y = Util.dezip(horiz) pylab.plot(x,y,'k--',lw=2) for c,k in iplot.subspec_enum(sorted(tunnels.keys())): x,y = Util.dezip(tunnels[k]) if not '.' in k: pylab.plot(x,y,'w:',lw=2) else: pylab.plot(x,y,c=c,lw=2) pylab.xlim(*prev_xlim) pylab.ylim(*prev_ylim[::-1])
reload(measure_burrows_from_image) cm_pts,horiz = measure_burrows_from_image.get_calibration(2) x,y = Util.dezip(horiz) plot(x,y,'k--') tunnels2 = {} reload(measure_burrows_from_image) reload(measure_burrows_from_image) tunnel_oris = {} measure_burrows_from_image.add_straight_tunnel(tunnels2,tunnel_oris,'MRhill') measure_burrows_from_image.add_straight_tunnel(tunnels2,tunnel_oris,'MR1.1','MRhill') measure_burrows_from_image.add_curved_tunnel(tunnels2,tunnel_oris,'MR1.2','MR1.1') measure_burrows_from_image.add_curved_tunnel(tunnels2,tunnel_oris,'MR1.2','MR1.1') measure_burrows_from_image.add_straight_tunnel(tunnels2,tunnel_oris,'MR1.3','MR1.2') tunnel_oris tunnels2 for c,k in iplot.subspec_enum(tunnels2): x,y = Util.dezip(tunnels2[k]) if not '.' in k: plot(x,y,'w--',lw=2) else: plot(x,y,c=c,lw=2) figure(2).axes[0].lines figure(2).axes[0].lines = [] for c,k in iplot.subspec_enum(tunnels2): x,y = Util.dezip(tunnels2[k]) if not '.' in k: plot(x,y,'w--',lw=2) else: plot(x,y,c=c,lw=2)