from scipy.spatial import ConvexHull
import matplotlib.pyplot as plt
import numpy as np

pts = np.random.rand(30, 2)
hull = ConvexHull(pts)
plt.plot(pts[:,0], pts[:,1], 'o')
for i in hull.simplices:
    plt.plot(pts[i, 0], pts[i, 1], 'k-')

plt.show()

其中simplex為索引點的序號，繪圖之后效果如下

ConvexHull有兩個可選參數(shù)，其中，incremental為布爾型參數(shù)，當(dāng)其為True時，允許添加新的點。

qhull_options的具體參數(shù)可以查看qhull，下面只演示一下QG。

QG

QGn表示將第n個點視為觀察點，在對點集進行凸包劃分后，如果把頂點連接起來，當(dāng)作一個圍墻，那么觀察點可以看得到的點，則標(biāo)記為good，其效果如下所示

pts = np.random.rand(1000, 2)
# 添加一個觀察點
pts = np.vstack([pts, np.array([[2,0.5]])])
hull = ConvexHull(pts, qhull_options='QG1000')
plt.plot(pts[:,0], pts[:,1], '.')
for i in hull.simplices:
    plt.plot(pts[i, 0], pts[i, 1], 'k-')

for i in hull.simplices[hull.good]:
    plt.plot(pts[i, 0],pts[i, 1], lw=5)

plt.show()

效果如圖所示

三維情況

二維情況下的凸包，很明顯是由線構(gòu)成的一個封閉圖形，而三維情況下的凸包，自然應(yīng)該是一個三維幾何體。拓展到任意維度，凸包構(gòu)成的實際上是一個單形，ConvexHull中的simplices便是構(gòu)成單形的點，在原點集中的索引。示例如下

pts = np.random.rand(30, 3)
hull = ConvexHull(pts)
ax = plt.subplot(projection='3d')
ax.scatter(pts[:,0], pts[:,1], pts[:,2])
for i in hull.simplices:
    ax.plot_trisurf(pts[i, 0], pts[i, 1], pts[i,2], alpha=0.5)

???????plt.show()

其中alpha參數(shù)用于調(diào)整三角面的透明度，從而可以透過凸包，看到凸包內(nèi)部的點。

效果如下

ConvexHull屬性

前面已經(jīng)引入了單形的概念，即凸包構(gòu)成的圖形便是單形。作為二維情況下的凸包，是由線段圍成；三維情況下的凸包，則是由平面圍成；推廣到任意維度，可以表述為構(gòu)成凸包的單形，由超曲面圍成。由于超曲面這個概念并沒有邊界，所以具有頂點、邊緣的凸包表面，下文中通稱為單形超表面。

ConvexHull類中常用的屬性如下

points 凸包包圍的點集
vertices 單形頂點在點集中的索引
simplices 單形超表面頂點
neighbors 超表面相鄰超表面的索引
equations 超曲面方程的參數(shù)

三維情況下的超曲面方程示例如下，即每個超曲面有4個參數(shù)

>>> hull.equations
array([[-0.5509472 ,  0.72386104, -0.41530999, -0.36369123],
       [-0.26155355,  0.16210178, -0.95147925,  0.02022163],
       [-0.99132368, -0.0460725 ,  0.12310441,  0.045523  ],
       [-0.98526526, -0.07170442,  0.15527666,  0.04749854],
       [-0.15900968, -0.98529789, -0.06248198,  0.13294496],
   # .......

到此這篇關(guān)于Python實現(xiàn)繪制凸包的示例代碼的文章就介紹到這了,更多相關(guān)Python繪制凸包內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！

您可能感興趣的文章: