最大の閉じた3D等高線プロットを見つける

Question

私は散布図と多くの等高線図を持っています。多くの輪郭があり、一部は閉じていない。 [this]から最大の閉じた輪郭の値を特定するにはどうしたらいいですか？ 1

Answer 1

最大の閉じた輪郭を見つけることは次のように行うことができます。ただし、これは「最大」とは最大のポイントツーポイント距離を意味することに注意してください。

他のサイズのメトリックも簡単に置換できます。

extend3d=Trueバージョンを使用している場合は、作成したPoly3DCollectionを少し操作する必要があります。

from mpl_toolkits.mplot3d import axes3d, art3d 
import matplotlib.pyplot as plt 
from matplotlib import cm 
from scipy.spatial import distance_matrix 
import numpy as np  

fig = plt.figure() 
ax = fig.gca(projection='3d') 
X, Y, Z = axes3d.get_test_data(0.05) 
cset = ax.contour(X, Y, Z, cmap=cm.coolwarm) 
ax.clabel(cset, fontsize=9, inline=1) 

maxsize = 0 
# Iterate over all the contour segments and find the largest 
for i, segs in enumerate(cset.allsegs): 
    for j, seg in enumerate(segs): 
     # First make sure it's closed 
     if (seg[0]-seg[-1]).any(): 
      continue 
     # Now get it's size 
     size = distance_matrix(seg, seg).max() 
     if size > maxsize: 
      maxsize = size 
      maxseg = (i, j) 

# Now highlight the "biggest" closed contour 
cset.collections[maxseg[0]].set_color('y') 
cset.collections[maxseg[0]].set_lw(5) 
pts2d = cset.allsegs[maxseg[0]][maxseg[1]] 
z = cset.levels[maxseg[0]] 
coords = np.c_[ pts2d, z*np.ones(pts2d.shape[0]) ] 
print(coords) 

plt.show() 

これは、その結果：あなたは閉じた輪郭がどのように大規模な定義はどうすればよい