colored wireframe plot in matplotlib
Asked Answered
M

3

22

I am trying to color a wireframe plot according to the z-value. I can't find any code examples on the internet.

Here is an example of a surface plot that has the colors I want and a wireframe plot where I can't manage to get the colors on the lines:

import numpy as np
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
import matplotlib.pyplot as plt

# some numbers for the data
P=12000 #W
Q=1     #kg/s
DT=3    #K
cp=4169.32  #J/kgK

dDT=np.logspace(-2,0,20,endpoint=True)
dQ=Q*np.logspace(-3,-1,20,endpoint=True)

# the plotting data
m1,m2=np.meshgrid(dDT,dQ)
err=cp*np.sqrt((m1*Q)**2+(m2*DT)**2)/P

# the wiremesh plot that i need fixed
fig=plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_wireframe(m1, m2, err, color=err/err.max(),cmap='jet')
ax.set_xlabel('dDT')
ax.set_ylabel('DQ')
ax.set_zlabel('relative error')

# the surface plot that has the colors i want
fig = plt.figure()
ax = fig.gca(projection='3d')

surf = ax.plot_surface(m1, m2, err,rstride=1,  cstride=1, cmap=cm.jet,
    linewidth=0.1, antialiased=False)

fig.colorbar(surf, shrink=0.5, aspect=5)

ax.set_xlabel('dDT')
ax.set_ylabel('DQ')
ax.set_zlabel('relative error')
plt.show()

Thank you for any help!

Male answered 28/2, 2013 at 11:11 Comment(1)
Possible duplicate of #24909756 ?Carpic
A
30

When you use plot_wireframe, each line can only have one color. Instead, you can use plot_surface. To get plot_surface to set the edgecolors, you need to give it facecolors. Then you can set the alpha of facecolors to zero.

from mpl_toolkits.mplot3d import axes3d
import matplotlib.pyplot as plt
from matplotlib import cm

X, Y, Z = axes3d.get_test_data(0.2)

# Normalize to [0,1]
norm = plt.Normalize(Z.min(), Z.max())
colors = cm.viridis(norm(Z))
rcount, ccount, _ = colors.shape

fig = plt.figure()
ax = fig.add_subplot(projection='3d')
surf = ax.plot_surface(X, Y, Z, rcount=rcount, ccount=ccount,
                       facecolors=colors, shade=False)
surf.set_facecolor((0,0,0,0))
plt.show()

color wireframe

Amathiste answered 5/7, 2017 at 23:34 Comment(6)
Thanks..Works well, but pretty slow, about the same as using light source. Pity that one has to normalize Z to get it to work, but I guess one could avoid that by tweaking the color map handling?Olvan
@Olvan I edited the code so that only the colors are normalized.Amathiste
this code does not render wireframe, but you posted wireframe picTerryterrye
@GrzegorzKrug did you edit the code? I get the same figure using the same code.Amathiste
@Amathiste I missed last command, sorry for inconveniance :) it worksTerryterrye
Hm, I get TypeError: FigureBase.gca() got an unexpected keyword argument 'projection' error...Merryman
N
4

I had a similar problem with coloring and sizing circles according to a variable which did not work either. So my workaround was to bin the values of the variables and loop over the bins. I masked the data such that the array mask did only contain the data with values in that bin.

ax.plot_wireframe(mask[i], ..., color="red")
ax.plot_wireframe(mask[i], ..., color="blue") 
etc.

I know it's not very elegant but in my case it did the job ;)

Nonpros answered 31/3, 2013 at 23:23 Comment(0)
N
3

Maybe you need to use plot_surface instead?

import matplotlib.pylab as plt
from matplotlib import cm 
from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure(figsize=(8, 8))
ax = fig.gca(projection='3d')

t = np.linspace(-3, 2, 31)
s = np.linspace(-3, 2, 31)

T, S = np.meshgrid(t, s)

ax.plot_surface(T * T, sqrt2 * T * S, S * S, cmap=cm.jet, rstride=1, cstride=1)

ax.set_xlabel('$t^2$')
ax.set_ylabel('$\sqrt{2} s t$')
ax.set_zlabel('$s^2$')

ax.set_title('line $s = t$ in $\cal F$')

plt.show()

enter image description here

Noshow answered 29/12, 2014 at 20:54 Comment(0)

© 2022 - 2024 — McMap. All rights reserved.