Tạo lưới Delaunay thông thường với scipy

Question

Có phương pháp nào để có được một tam giác trong 2D được đặt hàng nhiều hơn như Matlab Delaunay tạo ra không? Đây là một ví dụ về tam giác 2D Delaunay của Matlab.

Sử dụng mã này:

xPoints = np.arange(0,11,1) 
yPoints = np.arange(0,11,1) 
gridPoints = np.array([[x,y] for y in yPoints for x in xPoints]) 
tri = Delaunay(gridPoints) 
plt.triplot(gridPoints[:,0],gridPoints[:,1],tri.simplices.copy()) 
plt.plot(gridPoints[:,0],gridPoints[:,1],'bo') 
plt.title("Triangulation Visualization") 

tôi nhận được tam giác dưới đây:

Chú ý cách cung đường chéo trong Matlab kết quả tất cả đều có độ dốc tương tự; nhưng những người trong kết quả scipy là khác nhau. Kể từ Matlab và Scipy cả hai sử dụng QHull nội bộ, tôi đoán có một số phương pháp để bắt chước kết quả Matlab.

Answer 1

Bạn có thể thử tam giác thay vì Delaunay:

import matplotlib.tri as tri 

xlen = 10 
ylen = 16 
xPoints = np.arange(0,xlen+1,1) 
yPoints = np.arange(0,ylen+1,1) 

gridPoints = np.array([[[x,y] for y in yPoints] for x in xPoints]) 
a = [[i+j*(ylen+1),(i+1)+j*(ylen+1),i+(j+1)*(ylen+1)] for i in range(ylen) for j in range(xlen)] 
triang = tri.Triangulation(gridPoints[:,:,0].flatten(), gridPoints[:,:,1].flatten(),a) 

plt.triplot(triang) 
plt.plot(gridPoints[:,:,0],gridPoints[:,:,1],'bo') 
plt.title("Triangulation Visualization")