python-microscopy 25.11.20__cp313-cp313-win_amd64.whl

PYME/LMVis/glcanvas_core.py

@@ -0,0 +1,1157 @@
+#!/usr/bin/python
+
+##################
+# gl_render3D.py
+#
+# Copyright David Baddeley, Michael Graff
+# graff@hm.edu
+#
+# This program is free software: you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation, either version 3 of the License, or
+# (at your option) any later version.
+#
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with this program.  If not, see <http://www.gnu.org/licenses/>.
+#
+##################
+from math import floor
+import time
+
+import numpy
+import numpy as np
+import wx
+import wx.glcanvas
+import PYME.LMVis.mv_math as mm
+from OpenGL import GL, GLU
+#from OpenGL import GLU
+
+from PYME.LMVis.layers import AxesOverlayLayer, LUTOverlayLayer, ScaleBarOverlayLayer, SelectionOverlayLayer
+    #QuadTreeRenderLayer, VertexRenderLayer, , ShadedPointRenderLayer, \
+    #TetrahedraRenderLayer
+
+from PYME.LMVis.gl_offScreenHandler import OffScreenHandler
+from wx.glcanvas import GLCanvas
+
+from PYME.LMVis.shader_programs.ShaderProgramFactory import ShaderProgramFactory
+from PYME.LMVis.shader_programs.DefaultShaderProgram import DefaultShaderProgram
+
+from PYME import config
+from PYME.LMVis.views import View
+
+try:
+    from PYME.Analysis.points.gen3DTriangs import gen3DTriangs, gen3DBlobs
+except:
+    pass
+
+try:
+    # location in Python 2.7 and 3.1
+    from weakref import WeakSet
+except ImportError:
+    # separately installed py 2.6 compatibility
+    from weakrefset import WeakSet
+
+# import time
+
+from warnings import warn
+import logging
+logger = logging.getLogger(__name__)
+
+import sys
+
+# if sys.platform == 'darwin':
+#     # osx gives us LOTS of scroll events
+#     # ajust the mag in smaller increments
+#     ZOOM_FACTOR = 1.1
+# else:
+#     ZOOM_FACTOR = 2.0
+ZOOM_FACTOR = config.get('pymevis-zoom-factor', 1.1)
+
+# import statusLog
+
+name = 'ball_glut'
+
+from . import views
+from PYME.contrib import dispatch
+from PYME.ui import selection
+
+class SelectionSettings(selection.Selection):
+    def __init__(self):
+        selection.Selection.__init__(self)
+        self.colour = [1, 1, 0]
+        self.show = False
+
+class LMGLShaderCanvas(GLCanvas):
+    LUTOverlayLayer = None
+    AxesOverlayLayer = None
+    ScaleBarOverlayLayer = None
+    ScaleBoxOverlayLayer = None
+    _is_initialized = False
+
+    _stencil_shader = None
+
+    core_profile=True
+
+    def __init__(self, parent, show_lut=True, display_mode='2D', view=None):
+        print("New Canvas")
+        attribute_list = [wx.glcanvas.WX_GL_RGBA, wx.glcanvas.WX_GL_STENCIL_SIZE, 8, wx.glcanvas.WX_GL_DOUBLEBUFFER, wx.glcanvas.WX_GL_CORE_PROFILE]
+        self._num_antialias_samples = int(config.get('VisGUI-antialias_samples', 4))
+        if self._num_antialias_samples > 0:
+            attribute_list.extend([wx.glcanvas.WX_GL_SAMPLE_BUFFERS, 1, wx.glcanvas.WX_GL_SAMPLES, self._num_antialias_samples])
+        
+        GLCanvas.__init__(self, parent, -1, attribList=attribute_list)
+        self.Bind(wx.EVT_PAINT, self.OnPaint)
+        self.Bind(wx.EVT_SIZE, self.OnSize)
+        self.Bind(wx.EVT_MOUSEWHEEL, self.OnWheel)
+        self.Bind(wx.EVT_LEFT_DOWN, self.OnLeftDown)
+        self.Bind(wx.EVT_LEFT_UP, self.OnLeftUp)
+        self.Bind(wx.EVT_MIDDLE_DOWN, self.OnMiddleDown)
+        self.Bind(wx.EVT_MIDDLE_UP, self.OnMiddleUp)
+        self.Bind(wx.EVT_RIGHT_DOWN, self.OnMiddleDown)
+        self.Bind(wx.EVT_RIGHT_UP, self.OnMiddleUp)
+        self.Bind(wx.EVT_MOTION, self.OnMouseMove)
+        self.Bind(wx.EVT_KEY_DOWN, self.OnKeyPress)
+        # wx.EVT_MOVE(self, self.OnMove)
+        try:
+            self.gl_context = wx.glcanvas.GLContext(self)
+        except:
+            logger.exception('Error creating OpenGL context, try modifying the number of anti-aliasing samples using the "VisGUI-antialias_samples" config setting (0 to disable antialiasing)')
+            raise
+
+        self.displayMode = display_mode  # 3DPersp' #one of 3DPersp, 3DOrtho, 2D
+
+        self.parent = parent
+
+        self._scaleBarLength = 1000
+        self.clear_colour = [0,0,0,1.0]
+
+        self.LUTDraw = show_lut
+
+        self.zmin = -10
+        self.zmax = 10
+
+        self.angup = 0
+        self.angright = 0
+        
+        if view is None:
+            self.view = views.View()
+        else:
+            self.view = view
+
+        self.sx = 100
+        self.sy = 100
+
+        self.zc_o = 0
+
+        self.stereo = False
+
+        self.eye_dist = .01
+
+        self.dragging = False
+        self.panning = False
+
+        self.selectionSettings = SelectionSettings()
+        self.selectionDragging = False
+
+        self.layers = []
+        self.overlays = []
+        self.underlays = [] #draw these before data (assuming depth-testing is disabled)
+
+        self.wantViewChangeNotification = WeakSet()
+        self.pointSelectionCallbacks = []
+
+        self.on_screen = True
+        self.view_port_size = (self.Size[0], self.Size[1])
+        
+        self._old_bbox = None
+
+        self.layer_added = dispatch.Signal()
+
+        self._stencil_vao = None
+        self._stencil_vbo = None
+
+        self._stencil_size = (None,None)
+    
+    @property
+    def xc(self):
+        warn('Use view.translation[0] instead', DeprecationWarning, stacklevel=2)
+        return self.view.translation[0]
+    
+    @property
+    def yc(self):
+        warn('Use view.translation[1] instead', DeprecationWarning, stacklevel=2)
+        return self.view.translation[1]
+
+    @property
+    def zc(self):
+        warn('Use view.translation[2] instead', DeprecationWarning, stacklevel=2)
+        return self.view.translation[2]
+
+    @property
+    def scale(self):
+        warn('Use view.scale instead', DeprecationWarning, stacklevel=2)
+        return self.view.scale
+    
+    @property
+    def bounds(self):
+        return self.view.clipping
+
+    @property
+    def scaleBarLength(self):
+        return self._scaleBarLength
+
+    @scaleBarLength.setter
+    def scaleBarLength(self, value):
+        self._scaleBarLength = value
+        self.ScaleBarOverlayLayer.set_scale_bar_length(value)
+
+    def OnPaint(self, events):
+        if not self.IsShown():
+            print('ns')
+            return
+        
+        #wx.PaintDC(self)
+        #self.gl_context.SetCurrent(self)
+        self.SetCurrent(self.gl_context)
+
+        if not self._is_initialized:
+            self.initialize()
+        else:
+            self.OnDraw()
+        return
+    
+    @property
+    def bbox(self):
+        """ Bounding box in format [x0,y0,z0, x1, y1, z1]
+        
+        Calculated as the spanning box of all visible layers
+        """
+        
+        bb = np.zeros(6)
+        bb[:3] = 1e12
+        bb[3:] = -1e12
+        
+        nLayers = 0
+        
+        for l in self.layers:
+            # if getattr(l, 'visible', True):
+            if True:
+                bbl = l.bbox
+                
+                if not bbl is None:
+                    bb[:3] = np.minimum(bb[:3], bbl[:3])
+                    bb[3:] = np.maximum(bb[3:], bbl[3:])
+                    nLayers += 1
+                
+        if nLayers > 0:
+            #print('bbox: %s' % bb)
+            if not self._old_bbox is None and not np.allclose(self._old_bbox, bb):
+                #bbox has changed - update our cliping region
+                self.bounds['x'] = [bb[0]-1, bb[3]+1]
+                self.bounds['y'] = [bb[1]-1, bb[4]+1]
+                self.bounds['z'] = [bb[2]-1, bb[5]+1]
+            self._old_bbox = bb
+            return bb
+        else:
+            return None
+        
+
+    def initialize(self):
+        from .layers.ScaleBoxOverlayLayer import ScaleBoxOverlayLayer
+        self.InitGL()
+        self.ScaleBarOverlayLayer = ScaleBarOverlayLayer()
+        self.ScaleBoxOverlayLayer = ScaleBoxOverlayLayer()
+
+        self.LUTOverlayLayer = LUTOverlayLayer()
+        self.AxesOverlayLayer = AxesOverlayLayer()
+        
+        self.overlays.append(SelectionOverlayLayer(self.selectionSettings,))
+        self.underlays.append(self.ScaleBoxOverlayLayer)
+
+        self._is_initialized = True
+
+    def OnSize(self, event):
+        self.view_port_size = (self.Size[0], self.Size[1])
+        if self._is_initialized:
+            self.OnDraw()
+        self.Refresh()
+
+        #self.interlace_stencil()
+
+    def OnMove(self, event):
+        self.Refresh()
+        
+    def add_layer(self, layer, connect=True):
+        self.layers.append(layer)
+        if connect:
+            layer.on_update.connect(self.refresh)
+            
+        self.layer_added.send(self)
+
+    def get_session_info(self):
+        return {'layers': [l.serialise() for l in self.layers],
+                'view': self.view.to_json()}
+
+    def setOverlayMessage(self, message=''):
+        # self.messageOverlay.set_message(message)
+        # if self._is_initialized:
+        #     self.Refresh()
+        pass
+
+    def interlace_stencil(self):
+        window_width = self.view_port_size[0]
+        window_height = self.view_port_size[1]
+
+        #high dpi screens
+        sc = self.content_scale_factor
+
+        #print('scale factor:', sc)
+
+        # setting screen-corresponding geometry
+        GL.glViewport(0, 0, int(sc*window_width), int(sc*window_height))
+
+        with self.stencil_shader as sp:
+        
+            #mv = np.eye(4)
+            proj = mm.ortho(0, window_width -1, 0, window_height - 1, -1000, 1000)
+
+            sp.set_modelviewprojectionmatrix(proj)
+            
+
+            # clearing and configuring stencil drawing
+            if self.on_screen:
+                GL.glDrawBuffer(GL.GL_BACK)
+
+            GL.glEnable(GL.GL_STENCIL_TEST)
+            GL.glClearStencil(0)
+            GL.glClear(GL.GL_STENCIL_BUFFER_BIT)
+            GL.glStencilOp(GL.GL_REPLACE, GL.GL_REPLACE, GL.GL_REPLACE)  # colorbuffer is copied to stencil
+            GL.glDisable(GL.GL_DEPTH_TEST)
+            GL.glStencilFunc(GL.GL_ALWAYS, 1, 1)  # to avoid interaction with stencil content
+
+            # drawing stencil pattern
+            #GL.glColor4f(1, 1, 1, 0)  # alfa is 0 not to interfere with alpha tests
+
+            if self._stencil_size != (window_width, window_height):
+                self._stencil_size = (window_width, window_height)
+                start = self.ScreenPosition[1] % 2
+
+                # TODO - cache for a given screen size???
+                verts = np.hstack([[0, y, 0, window_width, y, 0] for y in range(start, int(window_height), 2)]).astype('f')
+                # alfa is 0 not to interfere with alpha tests
+                cols = np.hstack([[1, 1, 1, 0] for y in range(start, int(window_height), 2)]).astype('f')
+                
+                if self._stencil_vao is None:
+                    self._stencil_vao = GL.glGenVertexArrays(1)
+                    self._stencil_vbo = GL.glGenBuffers(3)
+
+                GL.glBindVertexArray(self._stencil_vao)
+                GL.glBindBuffer(GL.GL_ARRAY_BUFFER, self._stencil_vbo[0])
+                GL.glBufferData(GL.GL_ARRAY_BUFFER, verts, GL.GL_STATIC_DRAW)
+                GL.glVertexAttribPointer(0, 3, GL.GL_FLOAT, GL.GL_FALSE, 0, None)
+                GL.glEnableVertexAttribArray(0)
+                GL.glBindBuffer(GL.GL_ARRAY_BUFFER, self._stencil_vbo[1])
+                GL.glBufferData(GL.GL_ARRAY_BUFFER, 0*verts, GL.GL_STATIC_DRAW)
+                GL.glVertexAttribPointer(1, 3, GL.GL_FLOAT, GL.GL_FALSE, 0, None)
+                GL.glEnableVertexAttribArray(1)
+                GL.glBindBuffer(GL.GL_ARRAY_BUFFER, self._stencil_vbo[2])
+                GL.glBufferData(GL.GL_ARRAY_BUFFER, cols, GL.GL_STATIC_DRAW)
+                GL.glVertexAttribPointer(2, 4, GL.GL_FLOAT, GL.GL_FALSE, 0, None)
+                GL.glEnableVertexAttribArray(2)
+
+                self._stencil_count = len(verts)//3
+            else:
+                GL.glBindVertexArray(self._stencil_vao)
+
+            GL.glDrawArrays(GL.GL_LINES, 0, self._stencil_count)
+
+            GL.glStencilOp(GL.GL_KEEP, GL.GL_KEEP, GL.GL_KEEP)  # // disabling changes in stencil buffer
+        GL.glFlush()
+
+        # print 'is'
+
+    @property
+    def _stereo_views(self):
+        if self.displayMode == '2D':
+            return ['2D']
+        elif self.stereo:
+            return ['left', 'right']
+        else:
+            return ['centre']
+    
+    def OnDraw(self):
+        t0 = time.time()
+        self.SetCurrent(self.gl_context)
+        self.interlace_stencil()
+        GL.glEnable(GL.GL_DEPTH_TEST)
+        
+        GL.glClearColor(*self.clear_colour)
+        GL.glClear(GL.GL_COLOR_BUFFER_BIT)
+
+        #aspect ratio of window
+        ys = float(self.view_port_size[1]) / float(self.view_port_size[0])
+
+        for stereo_view in self._stereo_views:
+            if stereo_view == 'left':
+                eye = -self.eye_dist
+                GL.glStencilFunc(GL.GL_NOTEQUAL, 1, 1)
+            elif stereo_view == 'right':
+                eye = +self.eye_dist
+                GL.glStencilFunc(GL.GL_EQUAL, 1, 1)
+            else:
+                eye = 0
+
+            GL.glClear(GL.GL_DEPTH_BUFFER_BIT)
+
+            
+            #proj = glm.mat4(1.0)
+            # GL.glMatrixMode(GL.GL_MODELVIEW)
+            # GL.glLoadIdentity()
+            # GL.glMatrixMode(GL.GL_PROJECTION)
+            # GL.glLoadIdentity()
+
+            if self.displayMode == '3DPersp':
+                # our object will be be scaled to fit a 2x2x2 box at z=10 - see translate and scale calls below
+                #GL.glFrustum(-1 + eye, 1 + eye, ys, -ys, 8.5, 11.5)
+                proj = mm.frustum(-1 + eye, 1 + eye, ys, -ys, 8.5, 11.5)
+            else:
+                proj = mm.ortho(-1, 1, ys, -ys, -1000, 1000)
+                #GL.glOrtho(-1, 1, ys, -ys, -1000, 1000)
+
+            self.proj = proj
+
+            #GL.glMatrixMode(GL.GL_MODELVIEW)
+
+            self.mv = np.eye(4)#glm.mat4(1.0)
+
+            #stereo offset
+            self.mv = mm.translate(self.mv, eye, 0, 0) #glm.translate(self.mv, glm.vec3(eye, 0.0, 0.0))
+            #GL.glTranslatef(eye, 0.0, 0.0)
+
+            # move our object to be centred at -10
+            if self.displayMode == '3DPersp':
+                #GL.glTranslatef(0, 0, -10)
+                self.mv = mm.translate(self.mv, 0, 0, -10)
+
+            if not self.displayMode == '2D':
+                self.AxesOverlayLayer.render(self)
+
+            # scale object to fit a 2x2x2 box
+            #GL.glScalef(self.view.scale, self.view.scale, self.view.scale)
+            self.mv = mm.scale(self.mv, self.view.scale, self.view.scale, self.view.scale)
+
+            #GL.glPushMatrix()
+            _mv = self.mv
+
+            try:
+                # rotate object
+                #GL.glMultMatrixf(self.object_rotation_matrix)
+                self.mv = np.dot(self.mv,  self.object_rotation_matrix.T)
+                #GL.glTranslatef(-self.view.translation[0], -self.view.translation[1], -self.view.translation[2])
+                self.mv = mm.translate(self.mv, *(-self.view.translation))
+
+                self.mvp = np.dot(proj, self.mv)
+                
+                for l in self.underlays:
+                    l.render(self)
+    
+                
+                oit_layers = [] #layers which support order independent transparency - render these all together
+                for l in self.layers:
+                    if hasattr(l, 'engine') and (getattr(l.engine, 'use_oit', lambda l: False)(l) or getattr(l.engine, 'oit', False)):
+                        #defer rendering
+                        oit_layers.append(l)
+                    else:
+                        l.render(self)
+                        
+                if len(oit_layers) > 0:
+                    self.render_oit_layers(oit_layers)
+                    
+    
+                for o in self.overlays:
+                    o.render(self)
+            finally:
+                self.mv = _mv
+                #GL.glPopMatrix()
+
+            self.mvp = np.dot(proj, self.mv)
+            #scale bar gets drawn without the rotation
+            self.ScaleBarOverlayLayer.render(self)
+
+            if self.LUTDraw:
+                # set us up to draw in pixel coordinates
+                #GL.glMatrixMode(GL.GL_PROJECTION)
+                #GL.glLoadIdentity()
+                #GL.glOrtho(0, self.Size[0], self.Size[1], 0, -1000, 1000)
+                proj = mm.ortho(0, self.Size[0], self.Size[1], 0, -1000, 1000)
+                self.mvp = proj
+
+                #GL.glMatrixMode(GL.GL_MODELVIEW)
+                #GL.glLoadIdentity()
+                
+                self.LUTOverlayLayer.render(self)
+
+        GL.glFlush()
+
+        self.SwapBuffers()
+
+        #logger.debug('Frame rendered in: %3.4fs' % (time.time() - t0))
+
+    @property
+    def normal_matrix(self):
+        return np.linalg.inv(np.array(self.mv)[:3,:3].T)
+
+    def init_oit(self):
+        self._fb = GL.glGenFramebuffers(1)
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, self._fb)
+    
+        self._accumT = GL.glGenTextures(1)
+        GL.glBindTexture(GL.GL_TEXTURE_2D, self._accumT)
+    
+        self._revealT = GL.glGenTextures(1)
+        GL.glBindTexture(GL.GL_TEXTURE_2D, self._revealT)
+    
+        self._db = GL.glGenRenderbuffers(1)
+        GL.glBindRenderbuffer(GL.GL_RENDERBUFFER, self._db)
+    
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, 0)
+
+    def cleanup_oit(self):
+        GL.glDeleteFramebuffers(1, self._fb)
+        GL.glDeleteTextures(1, self._accumT)
+        GL.glDeleteTextures(1, self._revealT)
+        GL.glDeleteRenderbuffers(1, self._db)
+    
+    def resize_gl_oit(self, w, h):
+        #print('resize_gl')
+        self._oit_w = w
+        self._oit_h = h
+
+        # fix scaling error
+        # TODO - should this check be on wx, or OpenGL (or potentially python)
+        # TODO - do we need to do anything differnt on win?
+        sc = self.content_scale_factor
+
+        #sc = 1
+
+        w = int(sc*w)
+        h = int(sc*h)
+             
+        GL.glViewport(0, 0, w, h)
+        #self._accumdata = np.zeros([4,w,h], 'f')
+    
+        #print('bind fb')
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, self._fb)
+    
+        #print('bind accum texture')
+        GL.glBindTexture(GL.GL_TEXTURE_2D, self._accumT)
+        GL.glTexImage2D(GL.GL_TEXTURE_2D, 0, GL.GL_RGBA32F, w, h, 0, GL.GL_RGBA, GL.GL_FLOAT, None)
+    
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_MAG_FILTER, GL.GL_LINEAR)
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_MIN_FILTER, GL.GL_LINEAR)
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_WRAP_S, GL.GL_CLAMP_TO_EDGE)
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_WRAP_T, GL.GL_CLAMP_TO_EDGE)
+    
+        #print('framebuffer texture')
+        GL.glFramebufferTexture2D(GL.GL_FRAMEBUFFER, GL.GL_COLOR_ATTACHMENT0, GL.GL_TEXTURE_2D, self._accumT, 0)
+    
+        #print('bind reveal texture')
+        GL.glBindTexture(GL.GL_TEXTURE_2D, self._revealT)
+        GL.glTexImage2D(GL.GL_TEXTURE_2D, 0, GL.GL_R32F, w, h, 0, GL.GL_RED, GL.GL_FLOAT, None)
+    
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_MAG_FILTER, GL.GL_LINEAR)
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_MIN_FILTER, GL.GL_LINEAR)
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_WRAP_S, GL.GL_CLAMP_TO_EDGE)
+        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_WRAP_T, GL.GL_CLAMP_TO_EDGE)
+    
+        GL.glFramebufferTexture2D(GL.GL_FRAMEBUFFER, GL.GL_COLOR_ATTACHMENT1, GL.GL_TEXTURE_2D, self._revealT, 0)
+    
+        GL.glBindTexture(GL.GL_TEXTURE_2D, 0)
+    
+        #print('bind render buffer')
+        GL.glBindRenderbuffer(GL.GL_RENDERBUFFER, self._db)
+        GL.glRenderbufferStorage(GL.GL_RENDERBUFFER, GL.GL_DEPTH_COMPONENT, w, h)
+        GL.glFramebufferRenderbuffer(GL.GL_FRAMEBUFFER, GL.GL_DEPTH_ATTACHMENT, GL.GL_RENDERBUFFER, self._db)
+        GL.glBindRenderbuffer(GL.GL_RENDERBUFFER, 0)
+    
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, 0)
+    
+    def render_oit_layers(self, layers):
+        # initialise OIT if not already done
+        if not hasattr(self, '_fb'):
+            self.init_oit()
+            
+        # reallocate OIT textures if window has changed size
+        w, h = self.view_port_size
+        if not ((getattr(self, '_oit_w', None) == w) and (getattr(self,'_oit_h', None) == h)):
+            self.resize_gl_oit(w, h)
+    
+    
+        # draw to an offscreen framebuffer
+        current_fb = GL.glGetIntegerv(GL.GL_DRAW_FRAMEBUFFER_BINDING) #keep a reference to the current framebuffer so we can restore after we're done
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, self._fb) #bind our offscreen framebuffer
+    
+        #clear the offscreen framebuffer
+    
+        #glClearBufferfv(GL.GL_COLOR, 0, (0., 0., 0., 1.))
+        #glClearBufferfv(GL.GL_COLOR, 1, (1., 0., 0., 0.))
+        GL.glClearColor(0.0, 0.0, 0.0, 1.0)
+        GL.glClearDepth(1.0)
+        GL.glClear(GL.GL_COLOR_BUFFER_BIT | GL.GL_DEPTH_BUFFER_BIT)
+    
+        #Copy depth buffer from opaque stuff rendered previously
+        GL.glBindFramebuffer(GL.GL_DRAW_FRAMEBUFFER, self._fb)
+        GL.glBindFramebuffer(GL.GL_READ_FRAMEBUFFER, 0)
+        GL.glBlitFramebuffer(0, 0, self._oit_w, self._oit_h,
+                          0, 0, self._oit_w, self._oit_h,
+                          GL.GL_DEPTH_BUFFER_BIT, GL.GL_NEAREST)
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, self._fb)
+    
+        #glDrawBuffer(GL.GL_COLOR_ATTACHMENT0)
+        #glClear(GL.GL_COLOR_BUFFER_BIT) #| GL.GL_DEPTH_BUFFER_BIT)
+    
+        #clear the second colour attachment buffer
+        GL.glDrawBuffers(1, [GL.GL_COLOR_ATTACHMENT1, ])
+        GL.glClearColor(1.0, 0.0, 0.0, 1.0)
+        GL.glClear(GL.GL_COLOR_BUFFER_BIT)
+    
+        GL.glDrawBuffers(2, [GL.GL_COLOR_ATTACHMENT0, GL.GL_COLOR_ATTACHMENT1])
+        
+        #draw the OIT layers
+        for l in layers:
+            l.render(self)
+
+        # set the framebuffer back
+        GL.glBindFramebuffer(GL.GL_FRAMEBUFFER, current_fb)
+
+        #now do the compositing pass
+        with self.composite_shader as c:
+            # bind our pre-rendered textures
+
+            GL.glActiveTexture(GL.GL_TEXTURE0)
+            #GL.glEnable(GL.GL_TEXTURE_2D)
+            GL.glBindTexture(GL.GL_TEXTURE_2D, self._accumT)
+            self._acc_buf = GL.glGetTexImage(GL.GL_TEXTURE_2D, 0, GL.GL_RGBA, GL.GL_FLOAT)
+            GL.glUniform1i(c.get_uniform_location('accum_t'), 0)
+
+            GL.glActiveTexture(GL.GL_TEXTURE1)
+            #GL.glEnable(GL.GL_TEXTURE_2D)
+            GL.glBindTexture(GL.GL_TEXTURE_2D, self._revealT)
+            self._reveal_buf = GL.glGetTexImage(GL.GL_TEXTURE_2D, 0, GL.GL_RED, GL.GL_FLOAT)
+            GL.glUniform1i(c.get_uniform_location('reveal_t'), 1)
+
+            GL.glEnable(GL.GL_BLEND)
+            GL.glBlendFunc(GL.GL_ONE_MINUS_SRC_ALPHA, GL.GL_SRC_ALPHA)
+
+            # Draw triangles to display texture on
+            #GL.glColor4f(1., 1., 1., 1.)
+
+            vao = GL.glGenVertexArrays(1)
+            vbo = GL.glGenBuffers(1)
+            GL.glBindVertexArray(vao)
+            GL.glBindBuffer(GL.GL_ARRAY_BUFFER, vbo)
+            # GL.glBufferData(GL.GL_ARRAY_BUFFER, numpy.array([-1, -1, 0,  
+            #                                                  1, -1, 0, 
+            #                                                  1, 1, 0, 
+            #                                                  -1, 1, 0], 'f'), GL.GL_STATIC_DRAW)
+            GL.glBufferData(GL.GL_ARRAY_BUFFER, numpy.array([-1, -1, 0,  
+                                                             1, -1, 0, 
+                                                             1, 1, 0,
+                                                             1, 1, 0, 
+                                                             -1, 1, 0,
+                                                             -1, -1, 0], 'f'), GL.GL_STATIC_DRAW)
+            GL.glVertexAttribPointer(0, 3, GL.GL_FLOAT, GL.GL_FALSE, 0, None)
+            GL.glEnableVertexAttribArray(0)
+
+            #GL.glDrawArrays(GL.GL_TRIANGLE_STRIP, 0, 4)
+            GL.glDrawArrays(GL.GL_TRIANGLES, 0, 6)
+
+            GL.glBindBuffer(GL.GL_ARRAY_BUFFER, 0)
+            GL.glBindVertexArray(0)
+            GL.glDeleteBuffers(1, [vbo, ])
+            GL.glDeleteVertexArrays(1, [vao, ])
+
+            # GL.glBegin(GL.GL_QUADS)
+            # GL.glTexCoord2f(0., 0.) # lower left corner of image */
+            # GL.glVertex3f(-1, -1, 0.0)
+            # GL.glTexCoord2f(1., 0.) # lower right corner of image */
+            # GL.glVertex3f(1, -1, 0.0)
+            # GL.glTexCoord2f(1.0, 1.0) # upper right corner of image */
+            # GL.glVertex3f(1, 1, 0.0)
+            # GL.glTexCoord2f(0.0, 1.0) # upper left corner of image */
+            # GL.glVertex3f(-1, 1, 0.0)
+            # GL.glEnd()
+
+            #unbind our textures
+            GL.glActiveTexture(GL.GL_TEXTURE0)
+            GL.glBindTexture(GL.GL_TEXTURE_2D, 0)
+            #GL.glDisable(GL.GL_TEXTURE_2D)
+
+            GL.glActiveTexture(GL.GL_TEXTURE1)
+            GL.glBindTexture(GL.GL_TEXTURE_2D, 0)
+            #GL.glDisable(GL.GL_TEXTURE_2D)
+                
+    @property
+    def composite_shader(self):
+        from PYME.LMVis.shader_programs.ShaderProgramFactory import ShaderProgramFactory
+        from PYME.LMVis.shader_programs.oit_compositor import OITCompositorProgram
+        if not hasattr(self, '_composite_shader'):
+            self._composite_shader = ShaderProgramFactory.get_program(OITCompositorProgram, self.gl_context, self)
+            
+        return self._composite_shader
+
+    @property
+    def object_rotation_matrix(self):
+        """
+        The transformation matrix used to map coordinates in real space to our 3D view space. Currently implements
+        rotation, defined by 3 vectors (up, right, and back). Does not include scaling or projection.
+        
+        Returns
+        -------
+        a 4x4 matrix describing the rotation of the points within our 3D world
+        
+        """
+        if not self.displayMode == '2D':
+            return numpy.array(
+                [numpy.hstack((self.view.vec_right, 0)), numpy.hstack((self.view.vec_up, 0)), numpy.hstack((self.view.vec_back, 0)),
+                 [0, 0, 0, 1]], 'f')
+        else:
+            return numpy.eye(4, dtype='f')
+
+    def InitGL(self):
+        print('OpenGL - Version: {}'.format(GL.glGetString(GL.GL_VERSION)))
+        self.glsl_version = GL.glGetString(GL.GL_SHADING_LANGUAGE_VERSION).decode().replace('.', '')
+        print('Shader - Version: {}'.format(self.glsl_version))
+
+        #supported_glsl_versions = [GL.glGetStringi(GL.GL_SHADING_LANGUAGE_VERSION, i) for i in range(GL.glGetIntegerv(GL.GL_NUM_SHADING_LANGUAGE_VERSIONS))]
+        #print('Supported GLSL versions: {}'.format(supported_glsl_versions))
+        
+        max_samples = GL.glGetInteger(GL.GL_MAX_SAMPLES)
+        n_samples = GL.glGetInteger(GL.GL_SAMPLES)
+        print('GL.GL_MAX_SAMPLES: %d, GL.GL_SAMPLES: %d' % (max_samples, n_samples))
+        
+        if (max_samples >= 4) and (n_samples < 4):
+            logger.info('Your machine supports OpenGL antialiasing, but antialiasing disabled - enable by setting the "VisGUI-antialias_samples" PYME config setting to 4 or higher')
+            
+        GL.glEnable(GL.GL_DEPTH_TEST)
+        #GL.glEnable(GL.GL_NORMALIZE)
+        GL.glEnable(GL.GL_MULTISAMPLE)
+
+        #GL.glLoadIdentity()
+        #GL.glOrtho(self.xmin, self.xmax, self.ymin, self.ymax, self.zmin, self.zmax)
+        GL.glClearColor(0.0, 0.0, 0.0, 1.0)
+        GL.glClearDepth(1.0)
+
+        #GL.glEnableClientState(GL.GL_VERTEX_ARRAY)
+        #GL.glEnableClientState(GL.GL_COLOR_ARRAY)
+        #GL.glEnableClientState(GL.GL_NORMAL_ARRAY)
+        GL.glEnable(GL.GL_BLEND)
+        GL.glBlendFunc(GL.GL_SRC_ALPHA, GL.GL_ONE_MINUS_SRC_ALPHA)
+        #GL.glEnable(GL.GL_POINT_SMOOTH)
+
+        self.ResetView()
+
+    @property
+    def stencil_shader(self):
+        if self._stencil_shader is None:
+            self._stencil_shader =ShaderProgramFactory.get_program(DefaultShaderProgram, self.gl_context, self)
+        return self._stencil_shader
+
+    def ResetView(self):
+        self.view.vec_up = numpy.array([0, 1, 0])
+        self.view.vec_right = numpy.array([1, 0, 0])
+        self.view.vec_back = numpy.array([0, 0, 1])
+
+        self.Refresh()
+
+    @property
+    def xmin(self):
+        return self.view.translation[0] - 0.5 * self.pixelsize * self.view_port_size[0]
+
+    @property
+    def xmax(self):
+        return self.view.translation[0] + 0.5 * self.pixelsize * self.view_port_size[0]
+
+    @property
+    def ymin(self):
+        return self.view.translation[1] - 0.5 * self.pixelsize * self.view_port_size[1]
+
+    @property
+    def ymax(self):
+        return self.view.translation[1] + 0.5 * self.pixelsize * self.view_port_size[1]
+
+    @property
+    def pixelsize(self):
+        return 2. / (self.view.scale * self.view_port_size[0])
+    
+    @property
+    def content_scale_factor(self):
+        # version-safe way to get the content scale factor for high DPI screens
+        vmajor,vminor = wx.VERSION[:2]
+        if vmajor >= 4 and vminor >=1: # 4.1 or later 
+            sc = self.GetContentScaleFactor()
+        else:
+            sc = 1
+        return sc
+    
+    def _view_changed(self):
+        """Notify anyone who is synced to our view"""
+        for callback in self.wantViewChangeNotification:
+            if callback:
+                callback.Refresh()
+
+    def setView(self, xmin, xmax, ymin, ymax):
+
+        self.view.translation[0] = (xmin + xmax) / 2.0
+        self.view.translation[1] = (ymin + ymax) / 2.0
+        self.view.translation[2] = self.zc_o  # 0#z.mean()
+
+        self.view.scale = 2. / (xmax - xmin)
+
+        self.Refresh()
+        if 'OnGLViewChanged' in dir(self.parent):
+            self.parent.OnGLViewChanged()
+
+        self._view_changed()
+
+    def moveView(self, dx, dy, dz=0):
+        return self.pan(dx, dy, dz)
+
+    def pan(self, dx, dy, dz=0):
+        self.view.translation[0] += dx
+        self.view.translation[1] += dy
+        self.view.translation[2] += dz
+
+        self.Refresh()
+        self._view_changed()
+
+    def OnWheel(self, event):
+        rot = event.GetWheelRotation()
+        xp, yp = self._ScreenCoordinatesToNm(event.GetX(), event.GetY())
+
+        dx, dy = (xp - self.view.translation[0]), (yp - self.view.translation[1])
+        dx_, dy_, dz_, c_ = numpy.dot(self.object_rotation_matrix, [dx, dy, 0, 0])
+        xp_, yp_, zp_ = (self.view.translation[0] + dx_), (self.view.translation[1] + dy_), (self.view.translation[2] + dz_)
+
+        if event.MiddleIsDown():
+            self.WheelFocus(rot, xp_, yp_, zp_)
+        else:
+            self.WheelZoom(rot, xp_, yp_, zp_)
+
+    def WheelZoom(self, rot, xp, yp, zp=0):
+        dx = xp - self.view.translation[0]
+        dy = yp - self.view.translation[1]
+        dz = zp - self.view.translation[2]
+
+        if rot > 0:
+            # zoom out
+            self.view.scale *= ZOOM_FACTOR
+
+            self.view.translation[0] += dx * (1. - 1. / ZOOM_FACTOR)
+            self.view.translation[1] += dy * (1. - 1. / ZOOM_FACTOR)
+            self.view.translation[2] += dz * (1. - 1. / ZOOM_FACTOR)
+
+        if rot < 0:
+            # zoom in
+            self.view.scale /= ZOOM_FACTOR
+
+            self.view.translation[0] += dx * (1. - ZOOM_FACTOR)
+            self.view.translation[1] += dy * (1. - ZOOM_FACTOR)
+            self.view.translation[2] += dz * (1. - ZOOM_FACTOR)
+
+        self.Refresh()
+        self._view_changed()
+
+    def WheelFocus(self, rot, xp, yp, zp=0):
+        if rot > 0:
+            self.view.translation[2] -= 1.
+
+        if rot < 0:
+            self.view.translation[2] += 1.
+
+        self.Refresh()
+        self._view_changed()
+
+    def OnLeftDown(self, event):
+        if not self.displayMode == '2D':
+            # dragging the mouse rotates the object
+            self.xDragStart = event.GetX()
+            self.yDragStart = event.GetY()
+
+            self.angyst = self.angup
+            self.angxst = self.angright
+
+            self.dragging = True
+        else:  # 2D
+            # dragging the mouse sets an ROI
+            xp, yp = self._ScreenCoordinatesToNm(event.GetX(), event.GetY())
+
+            self.selectionDragging = True
+            self.selectionSettings.show = True
+
+            self.selectionSettings.start = (xp, yp, None)
+            self.selectionSettings.finish = (xp, yp, None)
+            self.selectionSettings.trace.clear()
+            self.selectionSettings.trace.append((xp, yp))
+
+        event.Skip()
+
+    def OnLeftUp(self, event):
+        self.dragging = False
+
+        if self.selectionDragging:
+            xp, yp = self._ScreenCoordinatesToNm(event.GetX(), event.GetY())
+
+            self.selectionSettings.finish = (xp, yp, None)
+            self.selectionSettings.trace.append((xp, yp))
+            self.selectionDragging = False
+
+            self.Refresh()
+            self.Update()
+
+        event.Skip()
+
+    def OnMiddleDown(self, event):
+        self.xDragStart = event.GetX()
+        self.yDragStart = event.GetY()
+
+        self.panning = True
+        event.Skip()
+
+    def OnMiddleUp(self, event):
+        self.panning = False
+        event.Skip()
+
+    def _ScreenCoordinatesToNm(self, x, y):
+        # FIXME!!!
+        sc = self.content_scale_factor
+        x_ = self.pixelsize * (x - 0.5 * float(self.view_port_size[0])) + self.view.translation[0]
+        y_ = self.pixelsize * (y - 0.5 * float(self.view_port_size[1])) + self.view.translation[1]
+        # print x_, y_
+        return x_, y_
+
+    def OnMouseMove(self, event):
+        x = event.GetX()
+        y = event.GetY()
+
+        if self.selectionDragging:
+            sx, sy = self._ScreenCoordinatesToNm(x, y)
+            self.selectionSettings.finish = (sx, sy, None)
+            self.selectionSettings.trace.append((sx, sy))
+            self.Refresh()
+            event.Skip()
+
+        elif self.dragging:
+            angx = numpy.pi * (x - self.xDragStart) / 180
+            angy = numpy.pi * (y - self.yDragStart) / 180
+
+            r_mat1 = numpy.matrix(
+                [[numpy.cos(angx), 0, numpy.sin(angx)], [0, 1, 0], [-numpy.sin(angx), 0, numpy.cos(angx)]])
+            r_mat = r_mat1 * numpy.matrix(
+                [[1, 0, 0], [0, numpy.cos(angy), numpy.sin(angy)], [0, -numpy.sin(angy), numpy.cos(angy)]])
+
+            vec_right_n = numpy.array(r_mat * numpy.matrix(self.view.vec_right).T).squeeze()
+            vec_up_n = numpy.array(r_mat * numpy.matrix(self.view.vec_up).T).squeeze()
+            vec_back_n = numpy.array(r_mat * numpy.matrix(self.view.vec_back).T).squeeze()
+
+            self.view.vec_right = vec_right_n
+
+            self.view.vec_up = vec_up_n
+            self.view.vec_back = vec_back_n
+
+            self.xDragStart = x
+            self.yDragStart = y
+
+            self.Refresh()
+            event.Skip()
+
+        elif self.panning:
+            dx = self.pixelsize * (x - self.xDragStart)
+            dy = self.pixelsize * (y - self.yDragStart)
+
+            dx_, dy_, dz_, c_ = numpy.dot(self.object_rotation_matrix, [dx, dy, 0, 0])
+
+            self.xDragStart = x
+            self.yDragStart = y
+
+            self.pan(-dx_, -dy_, -dz_)
+
+            event.Skip()
+
+    def OnKeyPress(self, event):
+        # print event.GetKeyCode()
+        if event.GetKeyCode() == 83:  # S - toggle stereo
+            self.stereo = not self.stereo
+            self.Refresh()
+        elif event.GetKeyCode() == 67:  # C - centre
+            self.recenter_bbox()
+            self.Refresh()
+
+        elif event.GetKeyCode() == 91:  # [ decrease eye separation
+            self.eye_dist /= 1.5
+            self.Refresh()
+
+        elif event.GetKeyCode() == 93:  # ] increase eye separation
+            self.eye_dist *= 1.5
+            self.Refresh()
+
+        elif event.GetKeyCode() == 82:  # R reset view
+            self.ResetView()
+            self.Refresh()
+
+        elif event.GetKeyCode() == 314:  # left
+            pos = numpy.array([self.view.translation[0], self.view.translation[1], self.view.translation[2]], 'f')
+            pos -= 300 * self.view.vec_right
+            self.view.translation[0], self.view.translation[1], self.view.translation[2] = pos
+            # print 'l'
+            self.Refresh()
+
+        elif event.GetKeyCode() == 315:  # up
+            pos = numpy.array([self.view.translation[0], self.view.translation[1], self.view.translation[2]])
+            pos -= 300 * self.view.vec_back
+            self.view.translation[0], self.view.translation[1], self.view.translation[2] = pos
+            self.Refresh()
+
+        elif event.GetKeyCode() == 316:  # right
+            pos = numpy.array([self.view.translation[0], self.view.translation[1], self.view.translation[2]])
+            pos += 300 * self.view.vec_right
+            self.view.translation[0], self.view.translation[1], self.view.translation[2] = pos
+            self.Refresh()
+
+        elif event.GetKeyCode() == 317:  # down
+            pos = numpy.array([self.view.translation[0], self.view.translation[1], self.view.translation[2]])
+            pos += 300 * self.view.vec_back
+            self.view.translation[0], self.view.translation[1], self.view.translation[2] = pos
+            self.Refresh()
+
+        else:
+            event.Skip()
+
+    def getSnapshot(self, mode=GL.GL_RGB):
+        # GL.glBindFramebuffer(GL.GL_READ_FRAMEBUFFER, 0)
+        # width, height = self.view_port_size[0], self.view_port_size[1]
+        # snap = GL.glReadPixelsf(0, 0, width, height, mode)
+        # snap = snap.ravel().reshape(width, height, -1, order='F')
+        #
+        # if mode == GL.GL_LUMINANCE:
+        #     # snap.strides = (4, 4 * snap.shape[0])
+        #     pass
+        # elif mode == GL.GL_RGB:
+        #     snap.strides = (12, 12 * snap.shape[0], 4)
+        # else:
+        #     raise RuntimeError('{} is not a supported format.'.format(mode))
+        # img.show()
+        self.on_screen = False
+        sc = self.content_scale_factor
+        view_port_size = (int(self.view_port_size[0]*sc), int(self.view_port_size[1]*sc))
+        off_screen_handler = OffScreenHandler(view_port_size, mode, self._num_antialias_samples)
+        with off_screen_handler:
+            self.OnDraw()
+        snap = off_screen_handler.get_snap()
+        self.on_screen = True
+        return  (255*np.asarray(snap)).astype('uint8')
+
+    def getIm(self, pixel_size=None, mode=GL.GL_RGB, image_bounds=None):
+        if ((pixel_size is None) or (abs(1 - pixel_size) < 0.001) and image_bounds is None):  # use current pixel size
+            return self.getSnapshot(mode=mode)
+        else:
+            # set size before moving the view, since self.setView includes a self.Refresh() call
+            old_view = self.get_view()
+            old_size = self.view_port_size
+
+            try:
+                if image_bounds is None:
+                    # for snapshots, just scale the viewport by the ratio of desired and current pixel sizes
+                    # as this works well with 3D rotated (and perspective) views
+                    # NOTE: rounding to the nearest pixel causes the pixel size to be inexact, and images produced this way should NOT
+                    # be used for quantification.
+                    sc = pixel_size/self.pixelsize
+                    self.view_port_size = (int(old_size[0]/sc), int(old_size[1]/sc))
+
+                else:
+                    # legacy fallback for CurrentRenderer
+                    # NOTE - does not work for 3D views - strictly 2D only
+                    x0, y0, x1, y1, z0, z1 = image_bounds.bounds 
+
+            
+                    self.view_port_size = (int((x1 - x0) / pixel_size),
+                                    int((y1 - y0) / pixel_size))
+                    logging.debug('viewport size %s' % (self.view_port_size,))
+                
+                    self.setView(x0, x0 + self.view_port_size[0]*pixel_size,
+                            y0, y0 + self.view_port_size[1]*pixel_size)
+                    
+                    #enforce top-down view (2D)
+                    self.view.top()
+
+                    print(pixel_size, self.pixelsize)
+                    assert(self.pixelsize == pixel_size)
+                
+                snap = self.getSnapshot(mode=mode)
+
+            finally:
+                self.view_port_size = old_size
+                self.set_view(old_view)
+            return snap
+
+    def recenter(self, x, y):
+        self.view.translation[0] = x.mean()
+        self.view.translation[1] = y.mean()
+        self.view.translation[2] = 0  # z.mean()
+
+        self.sx = x.max() - x.min()
+        self.sy = y.max() - y.min()
+        self.sz = 0  # z.max() - z.min()
+
+        self.view.scale = 2. / (max(self.sx, self.sy))
+        
+    def recenter_bbox(self):
+        bb = self.bbox
+        if bb is None:
+            return
+        
+        #print('recenter, bbox:', bb)
+        
+        centre = 0.5*(bb[:3] + bb[3:])
+        
+        self.view.translation[0], self.view.translation[1], self.view.translation[2] = centre
+        
+    def fit_bbox(self):
+        bb = self.bbox
+        if bb is None:
+            return
+    
+        sz = (bb[3:] - bb[:3])
+    
+        self.view.scale = 1.5/max(np.abs(sz))
+        self.recenter_bbox()
+
+    def set_view(self, view):
+        self.view=view
+        self.Refresh()
+
+    def get_view(self, view_id='id'):
+        view = View.copy(self.view)
+        view.view_id = view_id
+        return view
+    
+    def refresh(self, *args, **kwargs):
+        bb = self.bbox #force an update of our bounding box
+        self.Refresh()
+
+    
+
+
+def showGLFrame():
+    f = wx.Frame(None, size=(800, 800))
+    #sizer = wx.BoxSizer(wx.VERTICAL)
+    c = LMGLShaderCanvas(f)
+    #sizer.Add(c, 1, wx.EXPAND, 0)
+    #f.SetSizer(sizer)
+    f.Show()
+    return c