我想要三角化的三维点阵列为openFoam创建stls。我想要创造的是一个带有波纹的板材,它的长度通常在板式换热器中找到。对于stl的创建,我使用numpy stl,而对于三角剖分,我使用了来自枕的delauny三角测量。三角剖分本身创造了一个很好的标准,但它不考虑边缘,因此我不能创建波纹。
import numpy as np
from stl import mesh
import stl
from mpl_toolkits import mplot3d
from matplotlib import pyplot
import math as math
from scipy.spatial
我试图使用spatstat函数' Delaunay‘在R中生成delaunay三角剖分。但是,我查看了文档,似乎没有设置最大长度的参数。我注意到了这篇文章:How to set maximum length of triangle side in Delaunay triangulation in R?这似乎和我想做的事情是一样的,但因为我的点模式很大,所以我更喜欢一个简单而快速的解决方案。谢谢! 下面是我的代码: pts <- data.frame(readMat(paste('./TMA - Coordinates/HE_Rescaled_Coords/',
我使用Delaunay三角剖分来插值在一个规则的四维网格上的一组参数上计算的函数值。有时,当一个参数值发生少量的变化,从而使它变成一个新的单纯形时,单形中的多个点会发生变化。我希望,由于我不断地改变一个参数,我会从单纯形移动到单纯形,一次只改变单形中的一个点(通常在我的代码中也是这样)。相反,请考虑以下脚本:
import numpy as np
from scipy.spatial import Delaunay
# hideous construction to get the desired 4d grid of points
# with points at [-1, -0.5, 0
我只想计算两个点的Voronoi图: from scipy.spatial import Voronoi
vor = Voronoi([[0, 0], [0, 1]]) 但会得到以下错误: vor = Voronoi(points)
File "qhull.pyx", line 2518, in scipy.spatial.qhull.Voronoi.__init__
File "qhull.pyx", line 354, in scipy.spatial.qhull._Qhull.__init__
scipy.spatial.qhull.Qh
我在工作中遇到了一个任务,我需要使用Python在几个不同的多边形中创建三角形网格。虽然多边形是凸壳的分段,但每个分段本身都是凹形的。当我尝试使用scipy.spatial.Delaunay()时,该算法似乎假设多边形是凸的,并最终在凹壳区域之外生成三角形。代码和结果图如下所示。
for i in range(len(mp)):
dtri = np.array(mp[i])
plt.scatter(dtri[:,0],dtri[:,1])
plt.plot(np.array(fdc_seg[i])[:,0],np.array(fdc_seg[i])[:,1])