Просто додайте кілька додаткових думок, які можуть допомогти іншим із нерегулярними проблемами типу домену. У ситуації, коли у користувача є три вектори / списки, x, y, z, що представляють собою 2D-рішення, де z має бути нанесено на прямокутну сітку як поверхню, застосовуються коментарі ArtifixR 'plot_trisurf ()'. Аналогічний приклад, але з не прямокутним доменом:
import matplotlib.pyplot as plt
from matplotlib import cm
from mpl_toolkits.mplot3d import Axes3D
# problem parameters
nu = 50; nv = 50
u = np.linspace(0, 2*np.pi, nu,)
v = np.linspace(0, np.pi, nv,)
xx = np.zeros((nu,nv),dtype='d')
yy = np.zeros((nu,nv),dtype='d')
zz = np.zeros((nu,nv),dtype='d')
# populate x,y,z arrays
for i in range(nu):
for j in range(nv):
xx[i,j] = np.sin(v[j])*np.cos(u[i])
yy[i,j] = np.sin(v[j])*np.sin(u[i])
zz[i,j] = np.exp(-4*(xx[i,j]**2 + yy[i,j]**2)) # bell curve
# convert arrays to vectors
x = xx.flatten()
y = yy.flatten()
z = zz.flatten()
# Plot solution surface
fig = plt.figure(figsize=(6,6))
ax = Axes3D(fig)
ax.plot_trisurf(x, y, z, cmap=cm.jet, linewidth=0,
antialiased=False)
ax.set_title(r'trisurf example',fontsize=16, color='k')
ax.view_init(60, 35)
fig.tight_layout()
plt.show()
Вищевказаний код виробляє:
Однак це може не вирішити всіх проблем, зокрема, коли проблема визначена в неправильній області. Крім того, у випадку, коли домен має одну або більше увігнутих областей, триангуляція делауна може призвести до генерування помилкових трикутників, зовнішніх до домену. У таких випадках ці шахрайські трикутники повинні бути вилучені з триангуляції, щоб досягти правильного подання поверхні. У цих ситуаціях користувачеві, можливо, доведеться явно включити розрахунок триангуляції делауної, щоб ці трикутники можна було програмно видалити. За цих обставин наступний код може замінити попередній код сюжету:
import matplotlib.tri as mtri
import scipy.spatial
# plot final solution
pts = np.vstack([x, y]).T
tess = scipy.spatial.Delaunay(pts) # tessilation
# Create the matplotlib Triangulation object
xx = tess.points[:, 0]
yy = tess.points[:, 1]
tri = tess.vertices # or tess.simplices depending on scipy version
#############################################################
# NOTE: If 2D domain has concave properties one has to
# remove delaunay triangles that are exterior to the domain.
# This operation is problem specific!
# For simple situations create a polygon of the
# domain from boundary nodes and identify triangles
# in 'tri' outside the polygon. Then delete them from
# 'tri'.
# <ADD THE CODE HERE>
#############################################################
triDat = mtri.Triangulation(x=pts[:, 0], y=pts[:, 1], triangles=tri)
# Plot solution surface
fig = plt.figure(figsize=(6,6))
ax = fig.gca(projection='3d')
ax.plot_trisurf(triDat, z, linewidth=0, edgecolor='none',
antialiased=False, cmap=cm.jet)
ax.set_title(r'trisurf with delaunay triangulation',
fontsize=16, color='k')
plt.show()
Нижче наведено приклади графіків, що ілюструють рішення 1) із помилковими трикутниками та 2) де вони були видалені:
Я сподіваюсь, що вищесказане може допомогти людям із ситуацією увігнутості в даних рішення.
