Changeset 3889
- Timestamp:
- Apr 11, 2019 9:54:43 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIplot.py
r3888 r3889 7480 7480 7481 7481 # PlotStructure initialization here 7482 global mcsaXYZ,mcsaTypes,mcsaBonds 7482 global mcsaXYZ,mcsaTypes,mcsaBonds,txID 7483 7483 global cell, Vol, Amat, Bmat, A4mat, B4mat 7484 txID = 0 7484 7485 ForthirdPI = 4.0*math.pi/3.0 7485 7486 generalData = data['General'] … … 8222 8223 GL.glShadeModel(GL.GL_SMOOTH) 8223 8224 8224 def RenderViewPlane(plane, color,Z,width,height):8225 fade = list(color) + [.5,]8225 def RenderViewPlane(plane,Z,width,height): 8226 global txID 8226 8227 GL.glShadeModel(GL.GL_FLAT) 8227 ID = GL.glGenTextures(1) 8228 GL.glBindTexture(GL.GL_TEXTURE_2D, ID) 8228 newTX = False 8229 if not txID: 8230 txID = GL.glGenTextures(1) 8231 newTX = True 8232 GL.glBindTexture(GL.GL_TEXTURE_2D, txID) 8229 8233 GL.glPixelStorei(GL.GL_UNPACK_ALIGNMENT,1) 8230 8234 GL.glBlendFunc(GL.GL_SRC_ALPHA,GL.GL_ONE_MINUS_SRC_ALPHA) … … 8237 8241 GL.glFrontFace(GL.GL_CW) 8238 8242 GL.glTexEnvf(GL.GL_TEXTURE_ENV, GL.GL_TEXTURE_ENV_MODE, GL.GL_REPLACE) 8239 GL.glBindTexture(GL.GL_TEXTURE_2D, ID)8243 # GL.glTexEnvf(GL.GL_TEXTURE_ENV, GL.GL_ALPHA_SCALE, 1.0) 8240 8244 GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_BASE_LEVEL, 0) 8241 8245 GL.glTexParameteri(GL.GL_TEXTURE_2D, GL.GL_TEXTURE_MAX_LEVEL, 0) 8242 GL.glTexImage2D(GL.GL_TEXTURE_2D,0,GL.GL_RGBA,width,height,0,GL.GL_RGBA,GL.GL_UNSIGNED_BYTE,Z) 8246 if newTX: 8247 GL.glTexImage2D(GL.GL_TEXTURE_2D,0,GL.GL_RGBA,width,height,0,GL.GL_RGBA,GL.GL_UNSIGNED_BYTE,Z) 8248 else: 8249 GL.glTexSubImage2D(GL.GL_TEXTURE_2D,0,0,0,width,height,GL.GL_RGBA,GL.GL_UNSIGNED_BYTE,Z) 8243 8250 GL.glBegin(GL.GL_POLYGON) 8244 8251 for vertex,evertex in zip(plane,eBox): … … 8246 8253 GL.glVertex3fv(vertex) 8247 8254 GL.glEnd() 8248 # GL.glDrawPixels(width,height,GL.GL_RGBA,GL.GL_UNSIGNED_BYTE,Z)8249 8255 GL.glPopMatrix() 8250 8256 GL.glDisable(GL.GL_TEXTURE_2D) … … 8688 8694 SXYZ = np.reshape(np.inner(SXYZ,invModel[:3,:3].T)+VP[nxs,nxs,:],(-1,3)) 8689 8695 Z = np.reshape(G2mth.getRhos(SXYZ,rho),(npts,npts)) 8696 plt.cla() 8690 8697 plt.contour(Z,colors='k',linewidths=1) 8691 8698 plt.axis("off") … … 8695 8702 img, (width, height) = agg.print_to_buffer() 8696 8703 Zimg = np.frombuffer(img, np.uint8).reshape((height, width, 4)) 8697 RenderViewPlane(msize*eplane, Wt,Zimg,width,height)8704 RenderViewPlane(msize*eplane,Zimg,width,height) 8698 8705 8699 8706 # print time.time()-time0
Note: See TracChangeset
for help on using the changeset viewer.