Changeset 3946
 Timestamp:
 May 2, 2019 2:51:39 PM (2 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/GSASIIplot.py
r3945 r3946 117 117 import numpy.ma as ma 118 118 import numpy.linalg as nl 119 from scipy.ndimage.interpolation import map_coordinates 119 120 # Don't depend on wx/matplotlib for scriptable 120 121 try: … … 7712 7713 neqv = 0 7713 7714 for xyz,atyp in zip(XYZs,Types): 7714 equiv = G2spc.GenAtom(xyz,SGData,All=True,Move=False)7715 equiv = list(G2spc.GenAtom(xyz,SGData,All=True,Move=False)) 7715 7716 neqv = max(neqv,len(equiv)) 7716 7717 for item in equiv: … … 8893 8894 return 8894 8895 Rmax = np.max(rho)*drawingData['contourMax'] 8895 from scipy.ndimage.interpolation import map_coordinates8896 8896 from matplotlib.backends.backend_agg import FigureCanvasAgg 8897 8897 import matplotlib.pyplot as plt … … 8907 8907 if FourD: 8908 8908 SXYZT = np.vstack((SXYZ.T,np.inner(SXYZ,modQ)+G2frame.tau)).T 8909 Z = np.reshape(map_coordinates(rho4D,(SXYZT *rho4D.shape).T,order=1,mode='wrap'),(npts,npts))8909 Z = np.reshape(map_coordinates(rho4D,(SXYZT%1.*rho4D.shape).T,order=1,mode='wrap'),(npts,npts)) 8910 8910 else: 8911 Z = np.reshape(map_coordinates(rho,(SXYZ *rho.shape).T,order=1,mode='wrap'),(npts,npts))8911 Z = np.reshape(map_coordinates(rho,(SXYZ%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts)) 8912 8912 Z = np.where(Z<=Rmax,Z,Rmax) 8913 8913 plt.rcParams['figure.facecolor'] = (1.,1.,1.,.5)
Note: See TracChangeset
for help on using the changeset viewer.