# Axes # Graphics app = vv.use() #vv.clf() # # Superficie como gráfica de una función # f = lambda x, y: x*y xs = np.linspace(-2.0,2.0,21) ys = np.linspace(-2.0,2.0,21) zs = [ [ f(x,y) for x in xs ] for y in ys ] graph_surf = vv.grid(xs, ys, zs, axesAdjust=True) #graph_surf = vv.surf(xs, ys, zs, axesAdjust=True) graph_surf.colormap = vv.CM_JET # # Superficie paramétrica # def paramSurf(u,v): x = 2.0 * cos(u) * sin(v) y = 1.5 * sin(u) * sin(v) z = 2.2 * cos(v) return np.array( (x, y, z) ) NU = 21 NV = 11
#Plot in 3D with matplotlib from mpl_toolkits.mplot3d.axes3d import * import matplotlib.pyplot as plt from matplotlib import cm fig = plt.figure() ax = Axes3D(fig) ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=cm.jet,linewidth=1, antialiased=True) plt.show() #Plot in 3D with visvis import visvis f = visvis.gca() m = visvis.grid(xi,yi,Z) f.daspect = 1,1,10 # z x 10 # draped colors m = visvis.surf(xi,yi,Z) m.colormap = visvis.CM_JET #Save to GeoTiff file nrows,ncols = np.shape(Z) xres = (xmax-xmin)/ncols yres = (ymax-ymin)/nrows geotransform=(xmin,xres,0,ymin,0, yres) driver = gdal.GetDriverByName("GTiff") ds = driver.Create('output.tif', nrows,ncols, 1, gdal.GDT_Float32) ds.SetGeoTransform(geotransform) #Establish its coordinate
#'inverse': 1.0/sqrt((r/self.epsilon)**2 + 1) #x2 = np.linspace(min(x2), max(x2)) #y2 = np.linspace(min(y2), max(y2)) #X2, Y2 = np.meshgrid(x2, y2) # interpolación Z = spline(X, Y) Z1 = spline1(X, Y) #Visualización con visvis f = visvis.gca() m = visvis.plot(x, y, z, lc='k', ls='', mc='g', mw=10, lw=10, ms='.') f.daspect = 1, 1, 10 # z x 10 #m = visvis.surf(xi,yi,Z) m = visvis.surf(xi, yi, Z) m = visvis.grid(xi, yi, Z1) m.colormap = visvis.CM_JET f.axis.visible = True #m = visvis.surf(x2,y2,Z2) #m.colormap = visvis.CM_JET #Volumen en cada interpolación con respecto a los puntos: #volume = Convexhull(xyz).volume # Necesario para que la vista no se cierre app = visvis.use() #Movimiento cámara en 3 dimensiones (FPS) a = visvis.gca()
Keyword arguments ----------------- axesAdjust : bool If True, this function will call axes.SetLimits(), and set the camera type to 3D. If daspectAuto has not been set yet, it is set to False. axes : Axes instance Display the bars in the given axes, or the current axes if not given. Notes ----- * This function should not be confused with the axis grid, see the Axis.showGrid property. * This function is know in Matlab as mesh(), but to avoid confusion with the vv.Mesh class, it is called grid() in visvis. Also see surf() and the solid*() methods. """ m = vv.surf(*args, **kwargs) m.faceShading = None m.edgeShading = 'smooth' m.edgeColor = 'w' return m if __name__ == '__main__': vv.figure() m = vv.grid(vv.peaks()) m.colormap = vv.CM_HOT
----------------- axesAdjust : bool If True, this function will call axes.SetLimits(), and set the camera type to 3D. If daspectAuto has not been set yet, it is set to False. axes : Axes instance Display the bars in the given axes, or the current axes if not given. Notes ----- * This function should not be confused with the axis grid, see the Axis.showGrid property. * This function is know in Matlab as mesh(), but to avoid confusion with the vv.Mesh class, it is called grid() in visvis. Also see surf() and the solid*() methods. """ m = vv.surf(*args, **kwargs) m.faceShading = None m.edgeShading = 'smooth' m.edgeColor = 'w' return m if __name__ == '__main__': vv.figure() m = vv.grid(vv.peaks()) m.colormap = vv.CM_HOT