Changeset 3890


Ignore:
Timestamp:
Apr 12, 2019 9:30:04 AM (3 years ago)
Author:
vondreele
Message:

contour maps finished

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIphsGUI.py

    r3888 r3890  
    62406240                drawingData['unitCellBox'] = unitCellBox.GetValue()
    62416241                G2plt.PlotStructure(G2frame,data)
     6242                wx.CallAfter(UpdateDrawOptions)
    62426243   
    62436244            def OnShowHyd(event):
     
    63446345            showSizer.Add(line2Sizer)
    63456346           
    6346             line3Sizer = wx.BoxSizer(wx.HORIZONTAL)
    6347            
    6348             showCS = wx.CheckBox(drawOptions,-1,label=' Show contour slice?')
    6349             showCS.Bind(wx.EVT_CHECKBOX, OnShowSlice)
    6350             showCS.SetValue(drawingData['showSlice'])           
    6351             line3Sizer.Add(showCS,0,WACV)
    6352            
    6353             showSizer.Add(line3Sizer)
     6347            if generalData['Map']['rhoMax'] and drawingData['unitCellBox']:
     6348                line3Sizer = wx.BoxSizer(wx.HORIZONTAL)
     6349           
     6350                showCS = wx.CheckBox(drawOptions,-1,label=' Show contour slice? NB: can be very slow for large structures!')
     6351                showCS.Bind(wx.EVT_CHECKBOX, OnShowSlice)
     6352                showCS.SetValue(drawingData['showSlice'])           
     6353                line3Sizer.Add(showCS,0,WACV)
     6354               
     6355                showSizer.Add(line3Sizer)
    63546356           
    63556357            return showSizer
  • trunk/GSASIIplot.py

    r3889 r3890  
    75587558    wxOrange = wx.Colour(255,128,0)
    75597559    uBox = np.array([[0,0,0],[1,0,0],[1,1,0],[0,1,0],[0,0,1],[1,0,1],[1,1,1],[0,1,1]])
    7560     eBox = np.array([[.125,.875],[.125,.125],[.9,.125],[.9,.875],])
     7560    eBox = np.array([[0,1],[0,0],[1,0],[1,1],])
    75617561    eplane = np.array([[-1,-1,0],[-1,1,0],[1,1,0],[1,-1,0]])
    75627562    uEdges = np.array([
     
    82328232        GL.glBindTexture(GL.GL_TEXTURE_2D, txID)
    82338233        GL.glPixelStorei(GL.GL_UNPACK_ALIGNMENT,1)
     8234        GL.glEnable(GL.GL_BLEND)
    82348235        GL.glBlendFunc(GL.GL_SRC_ALPHA,GL.GL_ONE_MINUS_SRC_ALPHA)
    8235         GL.glEnable(GL.GL_BLEND)
    82368236        GL.glEnable(GL.GL_TEXTURE_2D)
    82378237        GL.glPushMatrix()
     
    82418241        GL.glFrontFace(GL.GL_CW)
    82428242        GL.glTexEnvf(GL.GL_TEXTURE_ENV, GL.GL_TEXTURE_ENV_MODE, GL.GL_REPLACE)
    8243 #        GL.glTexEnvf(GL.GL_TEXTURE_ENV, GL.GL_ALPHA_SCALE, 1.0)
     8243        GL.glTexEnvf(GL.GL_TEXTURE_ENV, GL.GL_ALPHA_SCALE, 1.0)
    82448244        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_BASE_LEVEL, 0)
    82458245        GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_MAX_LEVEL, 0)
     
    86808680                    RenderPlane(plane,color)
    86818681            if drawingData['showSlice']:
    8682                 rho = generalData['Map']['rho']
    8683                 if not len(rho):
     8682                if len(D4mapData.get('rho',[])):        #preferentially select 4D map if there
     8683                    rho = D4mapData['rho'][:,:,:,int(G2frame.tau*10)]   #pick current tau 3D slice
     8684                elif len(mapData['rho']):               #ordinary 3D map
     8685                    rho = mapData['rho']
     8686                else:
    86848687                    return
    86858688                from matplotlib.backends.backend_agg import FigureCanvasAgg
     
    86978700                plt.contour(Z,colors='k',linewidths=1)
    86988701                plt.axis("off")
     8702                plt.subplots_adjust(bottom=0.,top=1.,left=0.,right=1.,wspace=0.,hspace=0.)
    86998703                canvas = plt.get_current_fig_manager().canvas
    87008704                agg = canvas.switch_backends(FigureCanvasAgg)
     
    87028706                img, (width, height) = agg.print_to_buffer()
    87038707                Zimg = np.frombuffer(img, np.uint8).reshape((height, width, 4))
     8708                Zimg[:,:,3] = 64            #sets alpha to 25%
    87048709                RenderViewPlane(msize*eplane,Zimg,width,height)
    87058710               
Note: See TracChangeset for help on using the changeset viewer.