VTKPlots.py

import animate_helper
import warnings
import vtk
import vcs
import vcs2vtk
import numpy
from vtk.util import numpy_support as VN
import meshfill,boxfill,isofill,isoline
import os, traceback, sys
import cdms2
import DV3D
import MV2
import cdtime

def smooth(x,beta,window_len=11):
   """ kaiser window smoothing """
   # extending the data at beginning and at the end
   # to apply the window at the borders
   s = numpy.r_[x[window_len-1:0:-1],x,x[-1:-window_len:-1]]
   w = numpy.kaiser(window_len,beta)
   y = numpy.convolve(w/w.sum(),s,mode='valid')
   return y[(window_len/2):-(window_len/2)]

class VCSInteractorStyle(vtk.vtkInteractorStyleUser):

  def __init__(self,parent):
      self.AddObserver("LeftButtonPressEvent", parent.leftButtonPressEvent )
      self.AddObserver("LeftButtonReleaseEvent", parent.leftButtonReleaseEvent )
      self.AddObserver( "ModifiedEvent", parent.configureEvent )
      self.AddObserver( "ConfigureEvent", parent.configureEvent )
      if sys.platform == "darwin":
          self.AddObserver( "RenderEvent", parent.renderEvent )

class VTKVCSBackend(object):
  def __init__(self,canvas,renWin=None, debug=False,bg=None):
    self._lastSize = None
    self.canvas = canvas
    self.renWin = renWin
    self.debug = debug
    self.bg = bg
    self.type = "vtk"
    self.plotApps = {}
    self.plotRenderers = set()
    self.renderer = None
    self._renderers = {}
    self._plot_keywords = ['renderer','vtk_backend_grid','vtk_backend_geo']
    self.numberOfPlotCalls = 0
    self.renderWindowSize=None
    if renWin is not None:
      self.renWin = renWin
      if renWin.GetInteractor() is None and self.bg is False:
        self.createDefaultInteractor()
    self.createLogo()

    if sys.platform == "darwin":
        self.reRender = False

  def setAnimationStepper( self, stepper ):
      for plot in self.plotApps.values():
        plot.setAnimationStepper( stepper )

  def interact(self,*args,**kargs):
      if self.renWin is None:
          warnings.warn("Cannot interact if you did not open the canvas yet")
          return
      interactor = self.renWin.GetInteractor()
      ## Mac seems to handle events a bit differently
      ## Need to add observers on renWin
      ## Linux is fine w/o it so no need to do it
      if sys.platform == "darwin":
          self.renWin.AddObserver( "RenderEvent", self.renderEvent )
          self.renWin.AddObserver("LeftButtonPressEvent", self.leftButtonPressEvent )
          self.renWin.AddObserver("LeftButtonReleaseEvent", self.leftButtonReleaseEvent )
          self.renWin.AddObserver( "ModifiedEvent", self.configureEvent )
          self.renWin.AddObserver( "ConfigureEvent", self.configureEvent )
          self.renWin.AddObserver( "EndEvent",self.endEvent)
      if interactor is None:
          warnings.warn("Cannot start interaction. Blank plot?")
          return
      warnings.warn("Press 'Q' to exit interactive mode and continue script execution")
      interactor.Start()

  def endEvent(self,obj,event):
    if self.renWin is not None:
      if self.reRender:
        self.reRender = False
        self._lastSize = None
        self.renWin.Render()

  def renderEvent(self,caller,evt):
    renwin = self.renWin if (caller == None) else caller
    window_size = renwin.GetSize()
    if ( window_size <> self.renderWindowSize ):
      self.configureEvent(caller,evt)
      self.renderWindowSize = window_size

  def leftButtonPressEvent(self,obj,event):
    xy = self.renWin.GetInteractor().GetEventPosition()
    sz = self.renWin.GetSize()
    x = float(xy[0])/sz[0]
    y = float(xy[1])/sz[1]
    st = ""
    for dnm in self.canvas.display_names:
      d=vcs.elements["display"][dnm]
      if d.array[0] is None:
        continue
      t=vcs.elements["template"][d.template]
      gm = vcs.elements[d.g_type][d.g_name]
      if t.data.x1<=x<=t.data.x2 and t.data.y1<=y<=t.data.y2:
        ## Ok we clicked within template
        if numpy.allclose(gm.datawc_x1,1.e20):
          x1 = d.array[0].getAxis(-1)[0]
        else:
          x1 = gm.datawc_x1
        if numpy.allclose(gm.datawc_x2,1.e20):
          x2 = d.array[0].getAxis(-1)[-1]
        else:
          x2 = gm.datawc_x2
        if numpy.allclose(gm.datawc_y1,1.e20):
          y1 = d.array[0].getAxis(-2)[0]
        else:
          y1 = gm.datawc_y1
        if numpy.allclose(gm.datawc_y2,1.e20):
          y2 = d.array[0].getAxis(-2)[-1]
        else:
          y2 = gm.datawc_y2

        X = (x-t.data.x1)/(t.data.x2-t.data.x1)*(x2-x1)+x1
        Y = (y-t.data.y1)/(t.data.y2-t.data.y1)*(y2-y1)+y1
        # Ok we now have the X/Y values we need to figure out the indices
        try:
            I = d.array[0].getAxis(-1).mapInterval((X,X,'cob'))[0]
            try:
                J = d.array[0].getAxis(-2).mapInterval((Y,Y,'cob'))[0]
                # Values at that point
                V = d.array[0][...,J,I]
            except:
                V = d.array[0][...,I]
            if isinstance(V,numpy.ndarray):
              V=V.flat[0]
            try:
                st+="Var: %s\nX[%i] = %g\nY[%i] = %g\nValue: %g" % (d.array[0].id,I,X,J,Y,V)
            except:
                st+="Var: %s\nX = %g\nY[%i] = %g\nValue: %g" % (d.array[0].id,X,I,Y,V)
        except:
            st+="Var: %s\nX=%g\nY=%g\nValue = N/A" % (d.array[0].id,X,Y)
    ren = vtk.vtkRenderer()
    ren.SetBackground(.96,.96,.86)
    ren.SetViewport(x,y,min(x+.2,1.),min(y+.2,1))
    ren.SetLayer(self.renWin.GetNumberOfLayers()-1)
    a=vtk.vtkTextActor()
    a.SetInput(st)
    p=a.GetProperty()
    p.SetColor(0,0,0)
    bb = [0,0,0,0]
    a.GetBoundingBox(ren,bb)
    ps=vtk.vtkPlaneSource()
    ps.SetCenter(bb[0],bb[2],0.)
    ps.SetPoint1(bb[1],bb[2],0.)
    ps.SetPoint2(bb[0],bb[3],0.)
    ps.Update()
    m2d=vtk.vtkPolyDataMapper2D()
    m2d.SetInputConnection(ps.GetOutputPort())
    a2d=vtk.vtkActor2D()
    a2d.SetMapper(m2d)
    a2d.GetProperty().SetColor(.93,.91,.67)
    ren.AddActor(a2d)
    ren.AddActor(a)
    ren.ResetCamera()
    self.clickRenderer= ren
    self.renWin.AddRenderer(ren)
    self.renWin.Render()

  def leftButtonReleaseEvent(self,obj,event):
    self.clickRenderer.RemoveAllViewProps()
    self.renWin.RemoveRenderer(self.clickRenderer)
    self.renWin.Render()

  def configureEvent(self,obj,ev):
    sz = self.renWin.GetSize()
    if self._lastSize == sz:
      # We really only care about resize event
      # this is mainly to avoid segfault vwith Vistraisl which does
      # not catch configure Events but only modifiedEvents....
      if self.renWin is not None and sys.platform == "darwin":
        self.renWin.Render()
      return
    interactor = self.renWin.GetInteractor()
    self._lastSize = sz
    plots_args = []
    key_args =[]
    for dnm in self.canvas.display_names:
      d=vcs.elements["display"][dnm]
      parg = []
      for a in d.array:
        if a is not None:
          parg.append(a)
      parg.append(d._template_origin)
      parg.append(d.g_type)
      parg.append(d.g_name)
      plots_args.append(parg)
      if d.ratio is not None:
          key_args.append({"ratio":d.ratio})
      else:
          key_args.append({})
    self.canvas.clear()
    for i, pargs in enumerate(plots_args):
      self.canvas.plot(*pargs,**key_args[i])
    if self.renWin.GetSize()!=(0,0):
      self.scaleLogo()
    if self.renWin is not None and sys.platform == "darwin":
      self.renWin.Render()
    if sys.platform == "darwin":
        ## ON mac somehow we need to issue an extra Render after resize
        self.reRender = True

  def clear(self):
    if self.renWin is None: #Nothing to clear
          return
    renderers = self.renWin.GetRenderers()
    renderers.InitTraversal()
    ren = renderers.GetNextItem()
    hasValidRenderer = True if ren is not None else False
    while ren is not None:
        if not ren in self.plotRenderers:
            ren.RemoveAllViewProps()
            if not ren.GetLayer()==0:
              self.renWin.RemoveRenderer(ren)
        ren = renderers.GetNextItem()
    if hasValidRenderer and self.renWin.IsDrawable():
        self.renWin.Render()
    self.numberOfPlotCalls = 0
    self.createLogo()
    self._renderers = {}

  def createDefaultInteractor( self, ren=None ):
    defaultInteractor = self.renWin.GetInteractor()
    if defaultInteractor is None:
      defaultInteractor = vtk.vtkRenderWindowInteractor()
    self.vcsInteractorStyle = VCSInteractorStyle(self)
    if ren:
      self.vcsInteractorStyle.SetCurrentRenderer( ren )
    defaultInteractor.SetInteractorStyle( self.vcsInteractorStyle )
    defaultInteractor.SetRenderWindow(self.renWin)
    self.vcsInteractorStyle.On()

  def createRenWin(self,*args,**kargs):
    if self.renWin is None:
      # Create the usual rendering stuff.
      self.renWin = vtk.vtkRenderWindow()
      self.renWin.SetWindowName("VCS Canvas %i" % self.canvas._canvas_id)
      self.renWin.SetAlphaBitPlanes(1)
      ## turning off antialiasing by default
      ## mostly so that pngs are same accross platforms
      self.renWin.SetMultiSamples(0)
      self.initialSize()
      
    if self.renderer == None:
      self.renderer = self.createRenderer()
      if self.bg is False:
          self.createDefaultInteractor(self.renderer)
      self.renWin.AddRenderer(self.renderer)

  def createRenderer(self, *args, **kargs):
      # For now always use the canvas background
      ren = vtk.vtkRenderer()
      r,g,b = self.canvas.backgroundcolor
      ren.SetBackground(r/255., g/255., b/255.)
      return ren

  def update(self, *args, **kargs):
    self._lastSize=-1
    if self.renWin:
      self.configureEvent(None,None)

  def canvasinfo(self):
    if self.renWin is None:
      mapstate = False
      height = self.canvas.bgY
      width = self.canvas.bgX
      depth = None
      x=0
      y=0
    else:
      try: #mac but not linux
        mapstate = self.renWin.GetWindowCreated()
      except:
        mapstate = True
      width, height = self.renWin.GetSize()
      depth=self.renWin.GetDepthBufferSize()
      try: #mac not linux
        x,y = self.renWin.GetPosition()
      except:
        x,y = 0,0
    info = {
        "mapstate":mapstate,
        "height":height,
        "width":width,
        "depth":depth,
        "x":x,
        "y":y,
        }
    return info

  def orientation(self,*args,**kargs):
    if self.renWin is None:
      return "landscape"
    w,h = self.renWin.GetSize()
    if w>h:
      return "landscape"
    else:
      return "portrait"

  def portrait(self,W,H,x,y,clear):
      if clear:
          self.clear()
      if self.renWin is None:
          if W!=-99:
              self.canvas.bgX = W
              self.canvas.bgY = H
          else:
              W = self.canvas.bgX
              self.canvas.bgX = self.canvas.bgY
              self.canvas.bgY = W
      else:
          self.renWin.SetSize(W,H)

  def initialSize(self):
      #screenSize = self.renWin.GetScreenSize()
      self.renWin.SetSize(self.canvas.bgX,self.canvas.bgY)

  def open(self):
    self.createRenWin( open=True )

  def close(self):
    if self.renWin is None:
      return
    self.clear()
    self.renWin.Finalize()
    self.renWin = None

  def geometry(self,x,y,*args):
      self.renWin.SetSize(x,y)

  def flush(self):
      if self.renWin is not None:
          self.renWin.Render()

  def plot(self,data1,data2,template,gtype,gname,bg,*args,**kargs):
    self.numberOfPlotCalls+=1
    ## these are keyargs that can be reused later by the backend.
    returned = {}
    if self.bg is None:
      if bg:
        self.bg= True
      else:
        self.bg= False
    self.createRenWin(**kargs)
    if self.bg:
        self.renWin.SetOffScreenRendering(True)
        self.renWin.SetSize(self.canvas.bgX,self.canvas.bgY)
    self.cell_coordinates=kargs.get( 'cell_coordinates', None )
    if gtype == "text":
      tt,to = gname.split(":::")
      tt = vcs.elements["texttable"][tt]
      to = vcs.elements["textorientation"][to]
      gm=tt
    else:
      gm = vcs.elements[gtype][gname]
    tpl = vcs.elements["template"][template]

    if kargs.get("renderer",None) is None:
        if ( gtype in ["3d_scalar", "3d_dual_scalar", "3d_vector"] ) and (self.renderer <> None):
            ren = self.renderer
        else:
            #print "Calling 1"
            #ren = self.createRenderer()
            #if not (vcs.issecondaryobject(gm) and gm.priority==0):
            #    self.setLayer(ren,tpl.data.priority)
            #    self.renderer = ren
            #    self.renWin.AddRenderer(ren)
            pass
        #ren.SetPreserveDepthBuffer(True)
    else:
      ren = kargs["renderer"]

    vtk_backend_grid = kargs.get("vtk_backend_grid",None)
    vtk_backend_geo = kargs.get("vtk_backend_geo",None)
    if gtype in ["boxfill","meshfill","isofill","isoline"]:
      returned.update(self.plot2D(data1,data2,tpl,gm,vtk_backend_grid=vtk_backend_grid,vtk_backend_geo=vtk_backend_geo))
    elif gtype in ["3d_scalar", "3d_dual_scalar", "3d_vector"]:
      cdms_file = kargs.get( 'cdmsfile', None )
      cdms_var = kargs.get( 'cdmsvar', None )
      if not cdms_var is None:
          raise Exception()
      if not cdms_file is None:
          gm.addPlotAttribute( 'file', cdms_file )
          gm.addPlotAttribute( 'filename', cdms_file )
          gm.addPlotAttribute( 'url', cdms_file )
      returned.update(self.plot3D(data1,data2,tpl,gm,ren,**kargs))
    elif gtype in ["text"]:
      if tt.priority!=0:
        if not (None,None,None) in self._renderers.keys():
            ren = self.createRenderer()
            self.renWin.AddRenderer(ren)
            self.setLayer(ren,1)
            self._renderers[(None,None,None)]=ren
        else:
            ren = self._renderers[(None,None,None)]
        vcs2vtk.genTextActor(ren,to=to,tt=tt)
        self.setLayer(ren,tt.priority)
    elif gtype=="line":
      if gm.priority!=0:
        actors = vcs2vtk.prepLine(self.renWin,gm)
        for act,geo in actors:
            ren = self.fitToViewport(act,gm.viewport,wc=gm.worldcoordinate,geo=geo)
            ren.AddActor(act)
            self.setLayer(ren,gm.priority)
    elif gtype=="marker":
      if gm.priority!=0:
        actors = vcs2vtk.prepMarker(self.renWin,gm)
        for g,gs,pd,act,geo in actors:
            ren = self.fitToViewport(act,gm.viewport,wc=gm.worldcoordinate,geo=geo)
            ren.AddActor(act)
            self.setLayer(ren,gm.priority)
            # Add a transform to correct the glyph's aspect ratio:
            if pd is None and act.GetUserTransform():
              # Invert the scale of the actor's transform.
              glyphTransform = vtk.vtkTransform()
              scale = act.GetUserTransform().GetScale()
              xComp = scale[0]
              scale = [xComp / float(val) for val in scale]
              glyphTransform.Scale(scale)

              glyphFixer = vtk.vtkTransformPolyDataFilter()
              glyphFixer.SetTransform(glyphTransform)

              if pd is None:
                glyphFixer.SetInputConnection(gs.GetOutputPort())
              else:
                glyphFixer.SetInputData(pd)
                g.SetSourceData(None)

              g.SetSourceConnection(glyphFixer.GetOutputPort())
    elif gtype=="fillarea":
      if gm.priority!=0:
        actors = vcs2vtk.prepFillarea(self.renWin,gm)
        for act,geo in actors:
            ren = self.fitToViewport(act,gm.viewport,wc=gm.worldcoordinate,geo=geo)
            ren.AddActor(act)
            self.setLayer(ren,gm.priority)
    elif gtype=="1d":
      #self.renWin.AddRenderer(ren)
      returned.update(self.plot1D(data1,data2,tpl,gm))
    elif gtype=="vector":
      returned.update(self.plotVector(data1,data2,tpl,gm,vtk_backend_grid=vtk_backend_grid,vtk_backend_geo=vtk_backend_geo))
    else:
      raise Exception,"Graphic type: '%s' not re-implemented yet" % gtype
    self.scaleLogo()
    if not kargs.get("donotstoredisplay",False):
      self.renWin.Render()
    return returned


  def plot1D(self,data1,data2,tmpl,gm):
    Y = self.trimData1D(data1)
    if data2 is None:
      X=Y.getAxis(0)
    else:
      X=Y
      data1._yname = data2.id
      Y=self.trimData1D(data2)

    if gm.flip:
      tmp = Y
      Y = X
      X = tmp

    if gm.smooth is not None:
        Y = smooth(Y,gm.smooth)
    l = self.canvas.createline()
    Xs = X[:].tolist()
    Ys = Y[:].tolist()
    xs = []
    ys = []
    prev = None
    for i,v in enumerate(Ys):
        if v is not None and Xs[i] is not None: # Valid data
            if prev is None:
                prev=[]
                prev2 = []
            prev.append(Xs[i])
            prev2.append(v)
        else:
            if prev is not None:
                xs.append(prev)
                ys.append(prev2)
                prev = None
    if prev is not None:
        xs.append(prev)
        ys.append(prev2)
    l._x = xs
    l._y = ys
    l.color=gm.linecolor
    if gm.linewidth>0:
        l.width = gm.linewidth
    else:
        l.priority=0
    l.type = gm.line
    l._viewport = [tmpl.data.x1,tmpl.data.x2,tmpl.data.y1,tmpl.data.y2]
    # Also need to make sure it fills the whole space
    x1,x2,y1,y2 = vcs.utils.getworldcoordinates(gm,X,Y)
    if numpy.allclose(y1,y2):
        y1-=.0001
        y2+=.0001
    if numpy.allclose(x1,x2):
        x1-=.0001
        x2+=.0001
    l._worldcoordinate = [x1,x2,y1,y2]
    m=self.canvas.createmarker()
    m.type = gm.marker
    m.color = gm.markercolor
    if gm.markersize>0:
        m.size = gm.markersize
    else:
        m.priority=0
    m._x = l.x
    m._y = l.y
    m._viewport=l.viewport
    m._worldcoordinate = l.worldcoordinate

    if not (Y[:].min()>max(y1,y2) or Y[:].max()<min(y1,y2) \
            or X[:].min()>max(x1,x2) or X[:].max()<min(x1,x2)):
    	if l.priority>0:
            self.canvas.plot(l,donotstoredisplay=True)
        if m.priority>0:
            self.canvas.plot(m,donotstoredisplay=True)
    ren2 = self.createRenderer()
    self.renWin.AddRenderer(ren2)
    tmpl.plot(self.canvas,data1,gm,bg=self.bg,renderer=ren2,X=X,Y=Y)
    if hasattr(data1,"_yname"):
      del(data1._yname)
    del(vcs.elements["line"][l.name])
    del(vcs.elements["marker"][m.name])

    if tmpl.legend.priority>0:
        legd = self.canvas.createline()
        legd.x = [tmpl.legend.x1, tmpl.legend.x2]
        legd.y = [tmpl.legend.y1, tmpl.legend.y1]
        legd.color = l.color
        legd.width = l.width
        legd.type  = l.type
        t=self.canvas.createtext(To_source=tmpl.legend.textorientation,Tt_source=tmpl.legend.texttable)
        t.x=tmpl.legend.x2
        t.y=tmpl.legend.y2
        t.string=data1.id
        self.canvas.plot(t,donotstoredisplay=True)
        sp = t.name.split(":::")
        del(vcs.elements["texttable"][sp[0]])
        del(vcs.elements["textorientation"][sp[1]])
        del(vcs.elements["textcombined"][t.name])
        self.canvas.plot(legd,donotstoredisplay=True)
        del(vcs.elements["line"][legd.name])
    return {}

  def setLayer(self,renderer,priority):
    n = self.numberOfPlotCalls + (priority-1)*10000+1
    nMax = max(self.renWin.GetNumberOfLayers(),n+1)
    self.renWin.SetNumberOfLayers(nMax)
    renderer.SetLayer(n)
    pass

  def plot3D(self,data1,data2,tmpl,gm,ren,**kargs):
      from DV3D.Application import DV3DApp
      requiresFileVariable = True
      self.canvas.drawLogo = False
      if ( data1 is None ) or ( requiresFileVariable and not ( isinstance(data1, cdms2.fvariable.FileVariable ) or isinstance(data1, cdms2.tvariable.TransientVariable ) ) ):
          traceback.print_stack()
          raise Exception, "Error, must pass a cdms2 variable object as the first input to the dv3d gm ( found '%s')" % ( data1.__class__.__name__ )
      g = self.plotApps.get( gm, None )
      if g == None:
          g = DV3DApp( self.canvas, self.cell_coordinates )
          n_overview_points = 500000
          grid_coords = ( None, None, None, None )
          var_proc_op = None
          interface = None
          roi = None # ( 0, 0, 50, 50 )
          g.gminit( data1, data2, roi=roi, axes=gm.axes, n_overview_points=n_overview_points, n_cores=gm.NumCores, renwin=ren.GetRenderWindow(), plot_attributes=gm.getPlotAttributes(), gmname=gm.g_name, cm=gm.cfgManager, **kargs  ) #, plot_type = PlotType.List  )
          self.plotApps[ gm ] = g
          self.plotRenderers.add( g.plot.renderer )
      else:
          g.update( tmpl )
      return {}

  def onClosing(self):
      for plot in self.plotApps.values():
          if hasattr( plot, 'onClosing' ):
              plot.onClosing()

  def plotVector(self,data1,data2,tmpl,gm,vtk_backend_grid=None,vtk_backend_geo=None):
    #Preserve time and z axis for plotting these inof in rendertemplate
    returned = {}
    taxis = data1.getTime()
    if data1.ndim>2:
        zaxis = data1.getAxis(-3)
    else:
        zaxis = None
    data1 = self.trimData2D(data1) # Ok get3 only the last 2 dims
    data2 = self.trimData2D(data2)
    gridGenDict = vcs2vtk.genGridOnPoints(data1,data2,gm,deep=False,grid=vtk_backend_grid,geo=vtk_backend_geo)
    for k in ['vtk_backend_grid','xm','xM','ym','yM','continents','wrap','geo']:
        exec("%s = gridGenDict['%s']" % (k,k))
    returned["vtk_backend_grid"]=vtk_backend_grid
    returned["vtk_backend_geo"]=geo
    missingMapper = vcs2vtk.putMaskOnVTKGrid(data1,vtk_backend_grid,None,False,deep=False)
    returned["vtk_backend_missing_mapper"]=missingMapper

    u=numpy.ma.ravel(data1)
    v=numpy.ma.ravel(data2)
    sh = list(u.shape)
    sh.append(1)
    u = numpy.reshape(u,sh)
    v = numpy.reshape(v,sh)
    z = numpy.zeros(u.shape)
    w = numpy.concatenate((u,v),axis=1)
    w = numpy.concatenate((w,z),axis=1)

    # HACK The grid returned by vtk2vcs.genGrid is not the same size as the
    # data array. I'm not sure where the issue is...for now let's just zero-pad
    # data array so that we can at least test rendering until Charles gets
    # back from vacation:
    wLen = len(w)
    numPts = vtk_backend_grid.GetNumberOfPoints()
    if wLen != numPts:
        warnings.warn("!!! Warning during vector plotting: Number of points does not "\
              "match the number of vectors to be glyphed (%s points vs %s "\
              "vectors). The vectors will be padded/truncated to match for "\
              "rendering purposes, but the resulting image should not be "\
              "trusted."%(numPts, wLen))
        newShape = (numPts,) + w.shape[1:]
        w = numpy.ma.resize(w, newShape)

    w = vcs2vtk.numpy_to_vtk_wrapper(w,deep=False)
    w.SetName("vectors")
    vtk_backend_grid.GetPointData().AddArray(w)

    ## Vector attempt
    l = gm.line
    if l is None:
        l = "default"
    try:
      l = vcs.getline(l)
      lwidth = l.width[0]
      lcolor = l.color[0]
      lstyle = l.type[0]
    except:
      lstyle = "solid"
      lwidth = 1.
      lcolor = 0
    if gm.linewidth is not None:
        lwidth = gm.linewidth
    if gm.linecolor is not None:
        lcolor = gm.linecolor

    # Strip out masked points.
    if vtk_backend_grid.IsA("vtkStructuredGrid"):
        if vtk_backend_grid.GetCellBlanking():
            visArray = vtk_backend_grid.GetCellVisibilityArray()
            visArray.SetName("BlankingArray")
            vtk_backend_grid.GetCellData().AddArray(visArray)
            thresh = vtk.vtkThreshold()
            thresh.SetInputData(vtk_backend_grid)
            thresh.ThresholdByUpper(0.5)
            thresh.SetInputArrayToProcess(0, 0, 0,
                                          "vtkDataObject::FIELD_ASSOCIATION_CELLS",
                                          "BlankingArray")
            thresh.Update()
            vtk_backend_grid = thresh.GetOutput()
        elif vtk_backend_grid.GetPointBlanking():
            visArray = vtk_backend_grid.GetPointVisibilityArray()
            visArray.SetName("BlankingArray")
            vtk_backend_grid.GetPointData().AddArray(visArray)
            thresh = vtk.vtkThreshold()
            thresh.SetInputData(vtk_backend_grid)
            thresh.SetUpperThreshold(0.5)
            thresh.SetInputArrayToProcess(0, 0, 0,
                                          "vtkDataObject::FIELD_ASSOCIATION_POINTS",
                                          "BlankingArray")
            thresh.Update()
            vtk_backend_grid = thresh.GetOutput()

    arrow = vtk.vtkGlyphSource2D()
    arrow.SetGlyphTypeToArrow()
    arrow.FilledOff()

    glyphFilter = vtk.vtkGlyph2D()
    glyphFilter.SetSourceConnection(arrow.GetOutputPort())
    glyphFilter.SetVectorModeToUseVector()

    # Rotate arrows to match vector data:
    glyphFilter.OrientOn()

    # Scale to vector magnitude:
    glyphFilter.SetScaleModeToScaleByVector()

    # These are some unfortunately named methods. It does *not* clamp the scale
    # range to [min, max], but rather remaps the range [min, max]-->[0,1]. Bump
    # up min so that near-zero vectors will not be rendered, as these tend to
    # come out randomly oriented.
    glyphFilter.ClampingOn()
    glyphFilter.SetRange(0.01, 1.0)

    glyphFilter.SetInputArrayToProcess(1,0,0,0,"vectors")
    glyphFilter.SetScaleFactor(2.*gm.scale)

    glyphFilter.SetInputData(vtk_backend_grid)

    mapper = vtk.vtkPolyDataMapper()
    mapper.SetInputConnection(glyphFilter.GetOutputPort())
    act = vtk.vtkActor()
    act.SetMapper(mapper)
    try:
      cmap = vcs.elements["colormap"][cmap]
    except:
      cmap = vcs.elements["colormap"][self.canvas.getcolormapname()]
    r,g,b = cmap.index[lcolor]
    act.GetProperty().SetColor(r/100.,g/100.,b/100.)
    x1,x2,y1,y2 = vcs.utils.getworldcoordinates(gm,data1.getAxis(-1),data1.getAxis(-2))
    act = vcs2vtk.doWrap(act,[x1,x2,y1,y2],wrap)
    ren = self.fitToViewport(act,[tmpl.data.x1,tmpl.data.x2,tmpl.data.y1,tmpl.data.y2],[x1,x2,y1,y2])
    if tmpl.data.priority!=0:
        ren.AddActor(act)
        self.setLayer(ren,tmpl.data.priority)
    ren.AddActor(act)
    returned.update(self.renderTemplate(tmpl,data1,gm,taxis,zaxis))
    if self.canvas._continents is None:
      continents = False
    if continents:
        projection = vcs.elements["projection"][gm.projection]
        self.plotContinents(x1,x2,y1,y2,projection,wrap,tmpl)
    return returned

  def plot2D(self,data1,data2,tmpl,gm,vtk_backend_grid=None,vtk_backend_geo=None):
    #Preserve time and z axis for plotting these inof in rendertemplate
    returned = {}
    t = data1.getTime()
    if data1.ndim>2:
        z = data1.getAxis(-3)
    else:
        z = None
    data1 = self.trimData2D(data1) # Ok get3 only the last 2 dims
    if gm.g_name!="Gfm":
      data2 = self.trimData2D(data2)
    gridGenDict = vcs2vtk.genGrid(data1,data2,gm,deep=False,grid=vtk_backend_grid,geo=vtk_backend_geo)
    for k in ['vtk_backend_grid','xm','xM','ym','yM','continents','wrap','geo','cellData']:
        exec("%s = gridGenDict['%s']" % (k,k))
    returned["vtk_backend_grid"]=vtk_backend_grid
    returned["vtk_backend_geo"]=geo
    #Now applies the actual data on each cell
    if isinstance(gm,boxfill.Gfb) and gm.boxfill_type=="log10":
        data1=numpy.ma.log10(data1)
    data = vcs2vtk.numpy_to_vtk_wrapper(data1.filled(0.).flat, deep=False)
    if cellData:
        vtk_backend_grid.GetCellData().SetScalars(data)
    else:
        vtk_backend_grid.GetPointData().SetScalars(data)

    try:
      cmap = vcs.elements["colormap"][cmap]
    except:
      cmap = vcs.elements["colormap"][self.canvas.getcolormapname()]

    color = getattr(gm,"missing",None)
    if color is not None:
        color = cmap.index[color]
    missingMapper = vcs2vtk.putMaskOnVTKGrid(data1,vtk_backend_grid,color,cellData,deep=False)
    returned["vtk_backend_missing_mapper"]=missingMapper
    lut = vtk.vtkLookupTable()
    mn,mx=vcs.minmax(data1)
    #Ok now we have grid and data let's use the mapper
    mapper = vtk.vtkPolyDataMapper()
    legend = None
    if isinstance(gm,(meshfill.Gfm,boxfill.Gfb)):
      geoFilter = vtk.vtkDataSetSurfaceFilter()
      if cellData:
          p2c = vtk.vtkPointDataToCellData()
          p2c.SetInputData(vtk_backend_grid)
          geoFilter.SetInputConnection(p2c.GetOutputPort())
      else:
        geoFilter.SetInputData(vtk_backend_grid)
      geoFilter.Update()

    if isinstance(gm,(isofill.Gfi,isoline.Gi,meshfill.Gfm)) or \
        (isinstance(gm,boxfill.Gfb) and gm.boxfill_type=="custom"):

      #Now this filter seems to create the good polydata
      sFilter = vtk.vtkDataSetSurfaceFilter()
      if cellData:
          # Sets data to point instead of just cells
          c2p = vtk.vtkCellDataToPointData()
          c2p.SetInputData(vtk_backend_grid)
          c2p.Update()
          if self.debug:
            vcs2vtk.dump2VTK(c2p)
          #For contouring duplicate points seem to confuse it
          if vtk_backend_grid.IsA("vtkUntructuredGrid"):
              cln = vtk.vtkCleanUnstructuredGrid()
              cln.SetInputConnection(c2p.GetOutputPort())
              if self.debug:
                vcs2vtk.dump2VTK(cln)
              sFilter.SetInputConnection(cln.GetOutputPort())
          else:
              sFilter.SetInputConnection(c2p.GetOutputPort())
      else:
          sFilter.SetInputData(vtk_backend_grid)
      sFilter.Update()
      if self.debug:
        vcs2vtk.dump2VTK(sFilter)
      if isinstance(gm,isoline.Gi):
        cot = vtk.vtkContourFilter()
        if cellData:
          cot.SetInputData(sFilter.GetOutput())
        else:
          cot.SetInputData(vtk_backend_grid)


      levs = gm.levels
      if (isinstance(gm,isoline.Gi) and numpy.allclose( levs[0],[0.,1.e20])) or numpy.allclose(levs,1.e20):
        levs = vcs.mkscale(mn,mx)
        if len(levs)==1: # constant value ?
          levs = [levs[0],levs[0]+.00001]
        Ncolors = len(levs)
        if isinstance(gm,(isofill.Gfi,meshfill.Gfm)):
          levs2 = vcs.mkscale(mn,mx)
          if len(levs2)==1: # constant value ?
            levs2 = [levs2[0],levs2[0]+.00001]
          levs=[]
          for i in range(len(levs2)-1):
            levs.append([levs2[i],levs2[i+1]])
      else:
        if isinstance(gm.levels[0],(list,tuple)):
          if isinstance(gm,isoline.Gi):
            levs = [x[0] for x in gm.levels]
          else:
            levs = gm.levels
        else:
          levs = []
          levs2=gm.levels
          if numpy.allclose(levs2[0],1.e20):
            levs2[0]=-1.e20
          for i in range(len(levs2)-1):
            levs.append([levs2[i],levs2[i+1]])
          if isinstance(gm,isoline.Gi):
            levs = levs2
      Nlevs=len(levs)
      Ncolors = Nlevs
      ## Figure out colors
      if isinstance(gm,boxfill.Gfb):
        cols = gm.fillareacolors
        if cols is None:
          cols = vcs.getcolors(levs2,split=0)
      elif isinstance(gm,(isofill.Gfi,meshfill.Gfm)):
        cols = gm.fillareacolors
        if cols==[1,]:
          cols = vcs.getcolors(levs2,split=0)
          if isinstance(cols,(int,float)):
              cols=[cols,]
      elif isinstance(gm,isoline.Gi):
        cols = gm.linecolors

      if isinstance(gm,isoline.Gi):
        cot.SetNumberOfContours(Nlevs)
        if levs[0]==1.e20:
          levs[0]=-1.e20
        for i in range(Nlevs):
          cot.SetValue(i,levs[i])
        cot.SetValue(Nlevs,levs[-1])
        cot.Update()
        mapper.SetInputConnection(cot.GetOutputPort())
        mappers = []
      else:
        mappers = []
        LEVS = []
        INDX = []
        COLS = []
        indices = gm.fillareaindices
        if indices is None:
            indices=[1,]
        while len(indices)<len(cols):
            indices.append(indices[-1])
        for i,l in enumerate(levs):
            if i==0:
                C = [cols[i],]
                if numpy.allclose(levs[0][0],-1.e20):
                    ## ok it's an extension arrow
                    L=[mn-1.,levs[0][1]]
                else:
                    L = list(levs[i])
                I = [indices[i],]
            else:
                if l[0] == L[-1] and I[-1]==indices[i]:
                    # Ok same type lets keep going
                    if numpy.allclose(l[1],1.e20):
                        L.append(mx+1.)
                    else:
                        L.append(l[1])
                    C.append(cols[i])
                else: # ok we need new contouring
                    LEVS.append(L)
                    COLS.append(C)
                    INDX.append(I)
                    C = [cols[i],]
                    L = levs[i]
                    I = [indices[i],]
        LEVS.append(L)
        COLS.append(C)
        INDX.append(I)


        for i,l in enumerate(LEVS):
          # Ok here we are trying to group together levels can be, a join will happen if:
          # next set of levels contnues where one left off AND pattern is identical

          if isinstance(gm,isofill.Gfi):
              mapper = vtk.vtkPolyDataMapper()
              lut = vtk.vtkLookupTable()
              cot = vtk.vtkBandedPolyDataContourFilter()
              cot.ClippingOn()
              cot.SetInputData(sFilter.GetOutput())
              cot.SetNumberOfContours(len(l))
              cot.SetClipTolerance(0.)
              for j,v in enumerate(l):
                cot.SetValue(j,v)
              #cot.SetScalarModeToIndex()
              cot.Update()
              mapper.SetInputConnection(cot.GetOutputPort())
              lut.SetNumberOfTableValues(len(COLS[i]))
              for j,color in enumerate(COLS[i]):
                  r,g,b = cmap.index[color]
                  lut.SetTableValue(j,r/100.,g/100.,b/100.)
              mapper.SetLookupTable(lut)
              mapper.SetScalarRange(0,len(l)-1)
              mapper.SetScalarModeToUseCellData()
          else:
              for j,color in enumerate(COLS[i]):
                  mapper = vtk.vtkPolyDataMapper()
                  lut = vtk.vtkLookupTable()
                  th = vtk.vtkThreshold()
                  th.ThresholdBetween(l[j],l[j+1])
                  th.SetInputConnection(geoFilter.GetOutputPort())
                  geoFilter2 = vtk.vtkDataSetSurfaceFilter()
                  geoFilter2.SetInputConnection(th.GetOutputPort())
                  mapper.SetInputConnection(geoFilter2.GetOutputPort())
                  lut.SetNumberOfTableValues(1)
                  r,g,b = cmap.index[color]
                  lut.SetTableValue(0,r/100.,g/100.,b/100.)
                  mapper.SetLookupTable(lut)
                  mapper.SetScalarRange(l[j],l[j+1])
                  mappers.append([mapper,])

          #png = vtk.vtkPNGReader()
          #png.SetFileName("/git/uvcdat/Packages/vcs/Share/uvcdat_texture.png")
          #T=vtk.vtkTexture()
          #T.SetInputConnection(png.GetOutputPort())
          if isinstance(gm,isofill.Gfi):
              mappers.append([mapper,])

    else: #Boxfill (non custom)/Meshfill
      if isinstance(gm,boxfill.Gfb):
        if numpy.allclose(gm.level_1,1.e20) or numpy.allclose(gm.level_2,1.e20):
          levs = vcs.mkscale(mn,mx)
          if len(levs)==1: # constant value ?
              levs = [levs[0],levs[0]+.00001]
          legend = vcs.mklabels(levs)
          dx = (levs[-1]-levs[0])/(gm.color_2-gm.color_1+1)
          levs = numpy.arange(levs[0],levs[-1]+dx,dx)
        else:
          if gm.boxfill_type=="log10":
              levslbls = vcs.mkscale(numpy.ma.log10(gm.level_1),numpy.ma.log10(gm.level_2))
              levs = vcs.mkevenlevels(numpy.ma.log10(gm.level_1),
                      numpy.ma.log10(gm.level_2),
                      nlev=(gm.color_2-gm.color_1)+1)
          else:
              levslbls = vcs.mkscale(gm.level_1,gm.level_2)
              levs = vcs.mkevenlevels(gm.level_1,gm.level_2,nlev=(gm.color_2-gm.color_1)+1)
          if len(levs)>25:
              ## Too many colors/levels need to prettyfy this for legend
              legend = vcs.mklabels(levslbls)
              ## Make sure extremes are in
              legd2=vcs.mklabels([levs[0],levs[-1]])
              legend.update(legd2)
          else:
              legend = vcs.mklabels(levs)
          if gm.boxfill_type=="log10":
              for k in legend.keys():
                  legend[float(numpy.ma.log10(legend[k]))] = legend[k]
                  del(legend[k])

        cols = range(gm.color_1,gm.color_2+1)
      else:
        if numpy.allclose(gm.levels,1.e20):
          levs = vcs.mkscale(mn,mx)
        else:
          levs = gm.levels
          if numpy.allclose(levs[0],1.e20):
            levs[0]=-1.e20
        cols = gm.fillareacolors
        if cols==[1,]:
          cols = vcs.getcolors(levs)
      Nlevs = len(levs)
      Ncolors = Nlevs-1
      #Prep mapper
      mappers=[]
      mapper = vtk.vtkPolyDataMapper()
      thr = vtk.vtkThreshold()
      thr.SetInputConnection(geoFilter.GetOutputPort())
      if not gm.ext_1 in ["y",1,True]  and not gm.ext_2 in ["y",1,True] :
          thr.ThresholdBetween(levs[0],levs[-1])
      elif gm.ext_1 in ["y",1,True]  and not gm.ext_2 in ["y",1,True] :
          thr.ThresholdByLower(levs[-1])
      elif not gm.ext_1 in ["y",1,True]  and gm.ext_2 in ["y",1,True] :
          thr.ThresholdByUpper(levs[0])
      thr.Update()
      geoFilter2 = vtk.vtkDataSetSurfaceFilter()
      geoFilter2.SetInputConnection(thr.GetOutputPort())
      if gm.ext_1 in ["y",1,True]  and gm.ext_2 in ["y",1,True] :
          mapper.SetInputConnection(geoFilter.GetOutputPort())
      else:
          mapper.SetInputConnection(geoFilter2.GetOutputPort())

    if mappers == []: # ok didn't need to have special banded contours
      mappers=[mapper,]
      ## Colortable bit
      # make sure length match
      while len(cols)<Ncolors:
        cols.append(cols[-1])

      lut.SetNumberOfTableValues(Ncolors)
      for i in range(Ncolors):
        r,g,b = cmap.index[cols[i]]
        lut.SetTableValue(i,r/100.,g/100.,b/100.)

      mapper.SetLookupTable(lut)
      if numpy.allclose(levs[0],-1.e20):
        lmn = mn-1.
      else:
        lmn= levs[0]
      if numpy.allclose(levs[-1],1.e20):
        lmx = mx+1.
      else:
        lmx= levs[-1]
      mapper.SetScalarRange(lmn,lmx)

    if missingMapper is not None:
      if isinstance(gm,meshfill.Gfm):
        mappers.append(missingMapper)
      else:
        mappers.insert(0,missingMapper)

    if isinstance(gm,meshfill.Gfm):
        x1,x2,y1,y2 = vcs2vtk.getRange(gm,xm,xM,ym,yM)
    else:
        x1,x2,y1,y2 = vcs.utils.getworldcoordinates(gm,data1.getAxis(-1),data1.getAxis(-2))

    # Add a second mapper for wireframe meshfill:
    if isinstance(gm, meshfill.Gfm) and gm.mesh:
      lineMappers = []
      wireLUT = vtk.vtkLookupTable()
      wireLUT.SetNumberOfTableValues(1)
      wireLUT.SetTableValue(0,0,0,0)
      for polyMapper in mappers:
        if isinstance(polyMapper, list):
            polyMapper = polyMapper[0]
        lineMapper = vtk.vtkPolyDataMapper()
        lineMapper.SetInputConnection(polyMapper.GetInputConnection(0, 0))
        #lineMapper.SetInputConnection(geoFilter.GetOutputPort())
        lineMapper._useWireFrame = True

        # Setup depth resolution so lines stay above points:
        polyMapper.SetResolveCoincidentTopologyPolygonOffsetParameters(0, 1)
        polyMapper.SetResolveCoincidentTopologyToPolygonOffset()
        lineMapper.SetResolveCoincidentTopologyPolygonOffsetParameters(1, 1)
        lineMapper.SetResolveCoincidentTopologyToPolygonOffset()
        lineMapper.SetLookupTable(wireLUT)

        lineMappers.append(lineMapper)
      mappers.extend(lineMappers)

    if tmpl.data.priority != 0:
      # And now we need actors to actually render this thing
      for mapper in mappers:
        act = vtk.vtkActor()
        if isinstance(mapper,list):
          act.SetMapper(mapper[0])
        else:
          mapper.Update()
          act.SetMapper(mapper)
          if hasattr(mapper, "_useWireFrame"):
            prop = act.GetProperty()
            # Makes wireframed
            prop.SetRepresentationToWireframe()
        if geo is None:
          #If using geofilter on wireframed does not get wrppaed not sure why so sticking to many mappers
          act = vcs2vtk.doWrap(act,[x1,x2,y1,y2],wrap)
        if isinstance(mapper,list):
          ## This is the sport to add patterns
          #act.GetMapper().ScalarVisibilityOff()
          #act.SetTexture(mapper[1])
          pass
        # create a new renderer for this mapper
        # (we need one for each mapper because of cmaera flips)
        ren = self.fitToViewport(act,[tmpl.data.x1,tmpl.data.x2,tmpl.data.y1,tmpl.data.y2],wc=[x1,x2,y1,y2],geo=geo)
        if tmpl.data.priority>0:
            ren.AddActor(act)
            self.setLayer(ren,tmpl.data.priority)

    if isinstance(gm,meshfill.Gfm):
      tmpl.plot(self.canvas,data1,gm,
                bg=self.bg,
                X=numpy.arange(xm,xM*1.1,(xM-xm)/10.),
                Y=numpy.arange(ym,yM*1.1,(yM-ym)/10.))
    else:
      returned.update(self.renderTemplate(tmpl,data1,gm,t,z))
    if isinstance(gm,(isofill.Gfi,meshfill.Gfm,boxfill.Gfb)):
      if getattr(gm,"legend",None) is not None:
        legend = gm.legend
      if gm.ext_1 in ["y",1,True] and not numpy.allclose(levs[0],-1.e20):
          if isinstance(levs,numpy.ndarray):
              levs=levs.tolist()
          if not (isinstance(levs[0],list) and numpy.less_equal(levs[0][0],-1.e20)):
            levs.insert(0,-1.e20)
      if gm.ext_2 in ["y",1,True] and not numpy.allclose(levs[-1],1.e20):
          if isinstance(levs,numpy.ndarray):
              levs=levs.tolist()
          if not (isinstance(levs[-1],list) and numpy.greater_equal(levs[-1][-1],1.e20)):
            levs.append(1.e20)

      returned.update(self.renderColorBar(tmpl,levs,cols,legend,cmap))
    if self.canvas._continents is None:
      continents = False
    if continents:
        projection = vcs.elements["projection"][gm.projection]
        self.plotContinents(x1,x2,y1,y2,projection,wrap,tmpl)
    return returned

  def plotContinents(self,x1,x2,y1,y2,projection,wrap,tmpl):
      contData = vcs2vtk.prepContinents(self.canvas._continents)
      contMapper = vtk.vtkPolyDataMapper()
      contMapper.SetInputData(contData)
      contActor = vtk.vtkActor()
      contActor.SetMapper(contMapper)
      contActor.GetProperty().SetColor(0.,0.,0.)
      contActor = vcs2vtk.doWrap(contActor,[x1,x2,y1,y2],wrap,fastClip=False)
      if projection.type!="linear":
          contData=contActor.GetMapper().GetInput()
          cpts = contData.GetPoints()
          geo, gcpts = vcs2vtk.project(cpts,projection,[x1,x2,y1,y2])
          contData.SetPoints(gcpts)
          contMapper = vtk.vtkPolyDataMapper()
          contMapper.SetInputData(contData)
          contActor = vtk.vtkActor()
          contActor.SetMapper(contMapper)
          contActor.GetProperty().SetColor(0.,0.,0.)
      else:
          geo=None
      contMapper.SetResolveCoincidentTopologyPolygonOffsetParameters(1, 1)
      contMapper.SetResolveCoincidentTopologyToPolygonOffset()

      ren = self.fitToViewport(contActor,[tmpl.data.x1,tmpl.data.x2,tmpl.data.y1,tmpl.data.y2],wc=[x1,x2,y1,y2],geo=geo)
      if tmpl.data.priority!=0:
        self.setLayer(ren,tmpl.data.priority)
        ren.AddActor(contActor)
        col = ren.GetActors()
        col.InitTraversal()
        n=col.GetNextActor()
      return {}

  def renderTemplate(self,tmpl,data,gm,taxis,zaxis):
    tmpl.plot(self.canvas,data,gm,bg=self.bg)
    if taxis is not None:
        tstr = str(cdtime.reltime(taxis[0],taxis.units).tocomp(taxis.getCalendar()))
        #ok we have a time axis let's display the time
        crdate = vcs2vtk.applyAttributesFromVCStmpl(tmpl,"crdate")
        crdate.string = tstr.split()[0].replace("-","/")
        crtime = vcs2vtk.applyAttributesFromVCStmpl(tmpl,"crtime")
        crtime.string = tstr.split()[1]
        if not (None,None,None) in self._renderers.keys():
            ren = self.createRenderer()
            self.renWin.AddRenderer(ren)
            self.setLayer(ren,1)
            self._renderers[(None,None,None)]=ren
        else:
            ren = self._renderers[(None,None,None)]
        tt,to = crdate.name.split(":::")
        tt = vcs.elements["texttable"][tt]
        to = vcs.elements["textorientation"][to]
        if crdate.priority>0:
            vcs2vtk.genTextActor(ren,to=to,tt=tt)
        del(vcs.elements["texttable"][tt.name])
        del(vcs.elements["textorientation"][to.name])
        del(vcs.elements["textcombined"][crdate.name])
        tt,to = crtime.name.split(":::")
        tt = vcs.elements["texttable"][tt]
        to = vcs.elements["textorientation"][to]
        if crtime.priority>0:
            vcs2vtk.genTextActor(ren,to=to,tt=tt)
        del(vcs.elements["texttable"][tt.name])
        del(vcs.elements["textorientation"][to.name])
        del(vcs.elements["textcombined"][crtime.name])
    if zaxis is not None:
        # ok we have a zaxis to draw
        zname = vcs2vtk.applyAttributesFromVCStmpl(tmpl,"zname")
        zname.string=zaxis.id
        zvalue = vcs2vtk.applyAttributesFromVCStmpl(tmpl,"zvalue")
        if zaxis.isTime():
            zvalue.string = str(zaxis.asComponentTime()[0])
        else:
            zvalue.string= "%g" % zaxis[0]
        if not (None,None,None) in self._renderers.keys():
            ren = self.createRenderer()
            self.renWin.AddRenderer(ren)
            self.setLayer(ren,1)
            self._renderers[(None,None,None)]=ren
        else:
            ren = self._renderers[(None,None,None)]
        tt,to = zname.name.split(":::")
        tt = vcs.elements["texttable"][tt]
        to = vcs.elements["textorientation"][to]
        if zname.priority>0:
            vcs2vtk.genTextActor(ren,to=to,tt=tt)
        del(vcs.elements["texttable"][tt.name])
        del(vcs.elements["textorientation"][to.name])
        del(vcs.elements["textcombined"][zname.name])
        if hasattr(zaxis,"units"):
            zunits = vcs2vtk.applyAttributesFromVCStmpl(tmpl,"zunits")
            zunits.string=zaxis.units
            if zunits.priority>0:
                tt,to = zunits.name.split(":::")
                tt = vcs.elements["texttable"][tt]
                to = vcs.elements["textorientation"][to]
                vcs2vtk.genTextActor(ren,to=to,tt=tt)
                del(vcs.elements["texttable"][tt.name])
                del(vcs.elements["textorientation"][to.name])
                del(vcs.elements["textcombined"][zunits.name])
        tt,to = zvalue.name.split(":::")
        tt = vcs.elements["texttable"][tt]
        to = vcs.elements["textorientation"][to]
        if zvalue.priority>0:
            vcs2vtk.genTextActor(ren,to=to,tt=tt)
        del(vcs.elements["texttable"][tt.name])
        del(vcs.elements["textorientation"][to.name])
        del(vcs.elements["textcombined"][zvalue.name])
    return {}


  def renderColorBar(self,tmpl,levels,colors,legend,cmap):
    if tmpl.legend.priority>0:
      tmpl.drawColorBar(colors,levels,x=self.canvas,legend=legend,cmap=cmap)
    return {}

  def cleanupData(self,data):
      data[:] = numpy.ma.masked_invalid(data,numpy.nan)
      return data

  def trimData1D(self,data):
    if data is None:
      return None
    while len(data.shape)>1:
      data = data[0]
    return self.cleanupData(data)

  #ok now trying to figure the actual data to plot
  def trimData2D(self,data):
    if data is None:
      return None
    try:
      g=data.getGrid()
      gaxes=list(g.getAxisList())
      daxes=list(data.getAxisList())
      if daxes[len(daxes)-len(gaxes):] == gaxes:
        # Ok it is gridded and the grid axes are last
        return self.cleanupData(data(*(slice(0,1),)*(len(daxes)-len(gaxes)),squeeze=1))
      else:
        # Ok just return the last two dims
        return self.cleanupData(data(*(slice(0,1),)*(len(daxes)-2),squeeze=1))
    except Exception,err: # ok no grid info
      daxes=list(data.getAxisList())
      if cdms2.isVariable(data):
        return self.cleanupData( data(*(slice(0,1),)*(len(daxes)-2)))
      else: #numpy arrays are not callable
        op = ()
        for i in range(numpy.rank(data)-2):
          op.append(slice(0,1))
        return self.cleanupData(data[op])

  def put_png_on_canvas(self,filename,zoom=1,xOffset=0,yOffset=0,*args,**kargs):
      return self.put_img_on_canvas(filename,zoom,xOffset,yOffset,*args,**kargs)

  def put_img_on_canvas(self,filename,zoom=1,xOffset=0,yOffset=0,*args,**kargs):
    readerFactory = vtk.vtkImageReader2Factory()
    reader = readerFactory.CreateImageReader2(filename)
    reader.SetFileName(filename)
    reader.Update()
    imageData = reader.GetOutput()
    a = vtk.vtkImageActor()
    a.GetMapper().SetInputConnection(reader.GetOutputPort())
    origin = imageData.GetOrigin()
    spc = imageData.GetSpacing()
    ext = imageData.GetExtent()
    ren = self.createRenderer()
    cam = ren.GetActiveCamera()
    cam.ParallelProjectionOn()
    width = (ext[1]-ext[0])*spc[0]
    height = (ext[3]-ext[2])*spc[1]
    xoff = width*xOffset/zoom/200.
    yoff = height*yOffset/zoom/200.
    xc  = origin[0] + .5*(ext[0]+ext[1])*spc[0]
    yc  = origin[1] + .5*(ext[2]+ext[3])*spc[1]
    yd = (ext[3]-ext[2])*spc[1]
    d = cam.GetDistance()
    cam.SetParallelScale(.5*yd/zoom)
    cam.SetFocalPoint(xc+xoff,yc+yoff,0.)
    cam.SetPosition(xc+xoff,yc+yoff,d)
    ren.AddActor(a)
    self.renWin.AddRenderer(ren)
    self.renWin.Render()
    return

  def vectorGraphics(self, output_type, file, width=None, height=None, units=None):
    if self.renWin is None:
      raise Exception("Nothing on Canvas to dump to file")

    gl  = vtk.vtkGL2PSExporter()
    gl.SetInput(self.renWin)
    gl.SetCompress(0) # Do not compress
    gl.SetFilePrefix(".".join(file.split(".")[:-1]))
    if output_type=="svg":
        gl.SetFileFormatToSVG()
    elif output_type == "ps":
        gl.SetFileFormatToPS()
    elif output_type=="pdf":
        gl.SetFileFormatToPDF()
    else:
        raise Exception("Unknown format: %s" % output_type)
    gl.Write()

  def postscript(self, file, width=None, height=None, units=None,left=None,right=None,top=None,bottom=None):
      if right is not None:
          warnings.warn("the right_margin keyword for postscript has been deprecated in 2.0 and is being ignored")
      if left is not None:
          warnings.warn("the left_margin keyword for postscript has been deprecated in 2.0 and is being ignored")
      if top is not None:
          warnings.warn("the top_margin keyword for postscript has been deprecated in 2.0 and is being ignored")
      if bottom is not None:
          warnings.warn("the bottom_margin keyword for postscript has been deprecated in 2.0 and is being ignored")

      return self.vectorGraphics("ps", file, width, height, units)

  def pdf(self, file, width=None, height=None, units=None):
      return self.vectorGraphics("pdf", file, width, height, units)

  def svg(self, file, width=None, height=None, units=None):
      return self.vectorGraphics("svg", file, width, height, units)

  def gif(self,filename='noname.gif', merge='r', orientation=None, geometry='1600x1200'):
    raise RuntimeError("gif method not implemented in VTK backend yet")

  def png(self, file, width=None,height=None,units=None,draw_white_background = True, **args ):

        if self.renWin is None:
          raise Exception,"Nothing to dump aborting"

        if not file.split('.')[-1].lower() in ['png']:
            file+='.png'

        try:
          os.remove(file)
        except:
          pass

        #if width is not None and height is not None:
        #  self.renWin.SetSize(width,height)
          #self.renWin.Render()
        imgfiltr = vtk.vtkWindowToImageFilter()
        imgfiltr.SetInput(self.renWin)
#        imgfiltr.SetMagnification(3)
        ignore_alpha = args.get( 'ignore_alpha', False )
        if ignore_alpha or draw_white_background:
          imgfiltr.SetInputBufferTypeToRGB()
        else:
          imgfiltr.SetInputBufferTypeToRGBA()
        imgfiltr.Update()
        writer = vtk.vtkPNGWriter()
        writer.SetInputConnection(imgfiltr.GetOutputPort())
        writer.SetFileName(file)
        writer.Write()

  def cgm(self,file):
        if self.renWin is None:
          raise Exception,"Nothing to dump aborting"

        if not file.split('.')[-1].lower() in ['cgm']:
            file+='.cgm'

        try:
          os.remove(file)
        except:
          pass

        writer = vtk.vtkIOCGM.vtkCGMWriter()
        writer.SetFileName(file)
        R = self.renWin.GetRenderers()
        r=R.GetFirstRenderer()
        A = r.GetActors()
        A.InitTraversal()
        a = A.GetNextActor()
        while a is not None:
          m = a.GetMapper()
          m.Update()
          writer.SetInputData(m.GetInput())
          writer.Write()
          a=A.GetNextActor()
  def Animate(self,*args,**kargs):
    return VTKAnimate(*args,**kargs)

  def gettextextent(self,textorientation,texttable):
      warnings.warn("Please implement gettextextent for VTK Backend")

  def getantialiasing(self):
    if self.renWin is None:
      return 0
    else:
      return self.renWin.GetMultiSamples()

  def setantialiasing(self,antialiasing):
    if self.renWin is None:
      warnings.warn("no RenderWindow ready, skipping setantialiasing call, please reissue at a later time")
    else:
      self.renWin.SetMultiSamples(antialiasing)

  def createLogo(self):
    if self.canvas.drawLogo:
        defaultLogoFile = os.path.join(sys.prefix,"share","vcs","uvcdat.png")
        reader = vtk.vtkPNGReader()
        reader.SetFileName( defaultLogoFile )
        reader.Update()
        logo_input = reader.GetOutput()
        self.logoRepresentation = vtk.vtkLogoRepresentation()
        self.logoRepresentation.SetImage(logo_input)
        self.logoRepresentation.ProportionalResizeOn ()
        self.logoRepresentation.SetPosition( 0.882, 0.0 )
        self.logoRepresentation.SetPosition2( 0.10, 0.05 )
        self.logoRepresentation.GetImageProperty().SetOpacity( .8 )
        self.logoRepresentation.GetImageProperty().SetDisplayLocationToBackground()
        self.logoRenderer = vtk.vtkRenderer()
        self.logoRenderer.AddViewProp(self.logoRepresentation)
        self.logoRepresentation.SetRenderer(self.logoRenderer)

  def scaleLogo(self):
    if self.canvas.drawLogo:
        if self.renWin is not None:
            self.setLayer(self.logoRenderer,1)
            self.renWin.AddRenderer(self.logoRenderer)

  def fitToViewport(self,Actor,vp,wc=None,geo=None):
      ## Data range in World Coordinates
      vp=tuple(vp)
      if wc is None:
        Xrg = list(Actor.GetXRange())
        Yrg = list(Actor.GetYRange())
      else:
        Xrg=[float(wc[0]),float(wc[1])]
        Yrg=[float(wc[2]),float(wc[3])]
      if Yrg[0]>Yrg[1]:
        #Yrg=[Yrg[1],Yrg[0]]
        #T.RotateY(180)
        Yrg=[Yrg[1],Yrg[0]]
        flipY = True
      else:
        flipY = False
      if Xrg[0]>Xrg[1]:
        Xrg=[Xrg[1],Xrg[0]]
        flipX=True
      else:
        flipX=False

      if geo is not None:
       pt = vtk.vtkPoints()
       pt.SetNumberOfPoints(1)
       Xrg2 = [1.e20,-1.e20]
       Yrg2 = [1.e20,-1.e20]
       Npts=50.
       for x in numpy.arange(Xrg[0],Xrg[1],(Xrg[1]-Xrg[0])/Npts):
         for y in numpy.arange(Yrg[0],Yrg[1],(Yrg[1]-Yrg[0])/Npts):
           pt.SetPoint(0,x,y,0)
           pts = vtk.vtkPoints()
           geo.TransformPoints(pt,pts)
           b = pts.GetBounds()
           xm,xM,ym,yM=b[:4]
           if xm!=-numpy.inf:
             Xrg2[0]=min(Xrg2[0],xm)
           if xM!=numpy.inf:
             Xrg2[1]=max(Xrg2[1],xM)
           if ym!=-numpy.inf:
             Yrg2[0]=min(Yrg2[0],ym)
           if yM!=numpy.inf:
             Yrg2[1]=max(Yrg2[1],yM)
       Xrg=Xrg2
       Yrg=Yrg2
      sc = self.renWin.GetSize()
      wRatio = float(sc[0])/float(sc[1])
      dRatio = (Xrg[1]-Xrg[0])/(Yrg[1]-Yrg[0])
      vRatio = float(vp[1]-vp[0])/float(vp[3]-vp[2])


      if wRatio>1.: #landscape orientated window
          yScale = 1.
          xScale = vRatio*wRatio/dRatio
      else:
          xScale = 1.
          yScale = dRatio/(vRatio*wRatio)
      ## Ok now we know scaling and vp, let's see if we did this already.
      if (vp,xScale,yScale) in self._renderers.keys():
        #yep already have one, we will use this Renderer
        Renderer = self._renderers[(vp,xScale,yScale)]
        didRenderer = True
      else:
        Renderer = self.createRenderer()
        self.renWin.AddRenderer(Renderer)
        self._renderers[(vp,xScale,yScale)]=Renderer
        Renderer.SetViewport(vp[0],vp[2],vp[1],vp[3])
        didRenderer = False

      T = vtk.vtkTransform()
      T.Scale(xScale,yScale,1.)

      Actor.SetUserTransform(T)

      mapper = Actor.GetMapper()
      planeCollection = mapper.GetClippingPlanes()

      # We have to transform the hardware clip planes as well
      if (planeCollection is not None):
          planeCollection.InitTraversal()
          plane = planeCollection.GetNextItem()
          while (plane):
              origin = plane.GetOrigin()
              inOrigin = [origin[0], origin[1], origin[2], 1.0]
              outOrigin = [origin[0], origin[1], origin[2], 1.0]

              normal = plane.GetNormal()
              inNormal = [normal[0], normal[1], normal[2], 0.0]
              outNormal = [normal[0], normal[1], normal[2], 0.0]

              T.MultiplyPoint(inOrigin, outOrigin)
              if (outOrigin[3] != 0.0):
                  outOrigin[0] /= outOrigin[3]
                  outOrigin[1] /= outOrigin[3]
                  outOrigin[2] /= outOrigin[3]
              plane.SetOrigin(outOrigin[0], outOrigin[1], outOrigin[2])

              # For normal matrix, compute the transpose of inverse
              normalTransform = vtk.vtkTransform()
              normalTransform.DeepCopy(T)
              mat = vtk.vtkMatrix4x4()
              normalTransform.GetTranspose(mat)
              normalTransform.GetInverse(mat)
              normalTransform.SetMatrix(mat)
              normalTransform.MultiplyPoint(inNormal, outNormal)
              if (outNormal[3] != 0.0):
                  outNormal[0] /= outNormal[3]
                  outNormal[1] /= outNormal[3]
                  outNormal[2] /= outNormal[3]
              plane.SetNormal(outNormal[0], outNormal[1], outNormal[2])
              plane = planeCollection.GetNextItem()

      if not didRenderer:
        xc = xScale*float(Xrg[1]+Xrg[0])/2.
        yc = yScale*float(Yrg[1]+Yrg[0])/2.
        xd = xScale*float(Xrg[1]-Xrg[0])/2.
        yd = yScale*float(Yrg[1]-Yrg[0])/2.
        cam = Renderer.GetActiveCamera()
        cam.ParallelProjectionOn()
        cam.SetParallelScale(yd)
        cd = cam.GetDistance()
        cam.SetPosition(xc,yc,cd)
        cam.SetFocalPoint(xc,yc,0.)
        if geo is None:
          if flipY:
            cam.Elevation(180.)
            cam.Roll(180.)
            pass
          if flipX:
            cam.Azimuth(180.)
      return Renderer

class VTKAnimate(animate_helper.AnimationController):
   pass

# class VTKAnimate(animate_helper.animate_obj):
#   def __init__(self,*args,**kargs):
#     animate_helper.animate_obj.__init__(self,*args,**kargs)
#     self._initial_blink_done = False
#   def draw2(self,frame):
#     if self.create_flg == 1:
#         self.current_frame = frame
#         kargs = {}
#         if self._initial_blink_done:
#           kargs["noblink"]=True
#         else:
#           self._initial_blink_done = True
#         self.vcs_self.backend.clear()
#         self.vcs_self.put_png_on_canvas(self.animation_files[frame],
#                 self.zoom_factor, self.vertical_factor, self.horizontal_factor,**kargs)
#         if animate_helper.hasPyQt:
#           self.signals.drew.emit()