这篇文章主要介绍“Python德劳内三角怎么实现”的相关知识，小编通过实际案例向大家展示操作过程，操作方法简单快捷，实用性强，希望这篇“Python德劳内三角怎么实现”文章能帮助大家解决问题。
对于熟悉matplotlib三维画图的人来说，最常用的应该是plot_surface，但这个函数的绘图逻辑是，将xy平面映射到z轴，所以没法一次性绘制球，只能把球分成两半，上半球和下半球分别绘制。
如果想一次性绘制封闭图形，则可通过tri_surface，其绘图逻辑便是将图形拆分成一个个三角面，然后在对这些三角面进行绘制。所以，将一个曲面拆分成三角面，便构成了一个非常现实的问题，德劳内三角剖分便是建立在这个问题背景之下的。
scipy.spatial中提供了Delaunay类，下面以二维散点为例，来初步认识一下。
计算
安全
数据库
网络和加速
企业服务
import numpy as np
import matplotlib.pyplot as plt
from scipy.spatial import Delaunay
pts = np.array([[0, 0], [0, 1], [2, 0], [2, 1]])
tri = Delaunay(pts)
plt.triplot(pts[:,0], pts[:,1], tri.simplices)
plt.plot(pts[:,0], pts[:,1], 'o')
plt.show()效果如下
Delaunay的构造函数如下
Delaunay(points, furthest_site=False, incremental=False, qhull_options=None)各参数含义为
points 输入散点
furthest_site 为True时，计算最远点
incremental 为True时，允许增量添加点
qhull_options 为qhull参数，具体可参考qhull
在Delaunay对象中，有下面几个必须知道的常用属性
points 即输入的点集
simplices 三角面顶点在点集中的序号
neighbors 三角面相邻三角面的序号
equations 三角面方程
想要画个球，第一步是要得到一个球
# N为点数
def getBall(N):
pts = []
while len(pts) < N:
while True:
u = np.random.uniform(-1, 1)
v = np.random.uniform(-1, 1)
r2 = u**2 + v**2
if r2 < 1:
break
x = 2*u*np.sqrt(1-r2)
y = 2*v*np.sqrt(1-r2)
z = 1 - 2*r2
pts.append((x,y,z))
return np.vstack(pts)下面测试一下
pts = getBall(200)
ax = plt.subplot(projection='3d')
ax.scatter(pts[:,0], pts[:,1], pts[:,2])
plt.show()
接下来将这些随机点生成三角面，并进行绘图
tri = Delaunay(pts) ax = plt.subplot(projection='3d') for i in tri.simplices: ax.plot_trisurf(pts[i, 0], pts[i, 1], pts[i,2]) plt.show()
效果如下
看上去花花绿绿的这些三角形，便是通过德劳内三角剖分得到的，其equations属性可以查看这些三角面的方程参数
>>> tri.equations
array([[-2.35739179e-16, -1.64155539e-15, -1.54600295e-15,
-1.00000000e+00, 2.41181971e-16],
[-2.35739179e-16, -1.64155539e-15, -1.54600295e-15,
-1.00000000e+00, 2.41181971e-16],
[-2.35739179e-16, -1.64155539e-15, -1.54600295e-15,
-1.00000000e+00, 2.41181971e-16],
...,
[-2.35739179e-16, -1.64155539e-15, -1.54600295e-15,
-1.00000000e+00, 2.41181971e-16],
[-2.35739179e-16, -1.64155539e-15, -1.54600295e-15,
-1.00000000e+00, 2.41181971e-16],
[-2.35739179e-16, -1.64155539e-15, -1.54600295e-15,
-1.00000000e+00, 2.41181971e-16]])
``关于“Python德劳内三角怎么实现”的内容就介绍到这里了，感谢大家的阅读。如果想了解更多行业相关的知识，可以关注亿速云行业资讯频道，小编每天都会为大家更新不同的知识点。
亿速云「云服务器」，即开即用、新一代英特尔至强铂金CPU、三副本存储NVMe SSD云盘，价格低至29元/月。点击查看>>
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
原文链接：https://blog.csdn.net/m0_37816922/article/details/130140602
