윤곽선 경로 가끔 비움

Question

나는 윤곽선의 주위와 영역을 얻으려고 노력해왔다. (솔리드/하이브리드 로켓 모터의 단점을 모델링하기 위해 레벨 설정 방법을 사용하고있다.) 여기

matplotlib - extracting values from contour lines

오류 재현 최소한의 예입니다 : 내가 여기 시스템 권장 접근 방식을 촬영했습니다

import numpy as np 
import pylab as pl 
import skfmm 

n_pts = 1000 
r_f = 98 #mm (final radius) 
fin_l = 20 #mm (fin length) 
fin_w = 10 #mm (fin length) 

#set up grid 
X, Y = np.meshgrid(np.linspace(-r_f,r_f,n_pts), np.linspace(-r_f,r_f,n_pts)) 
initial_grid = -1 * np.ones_like(X) 

#set up initial shape 
initial_grid[np.logical_and(np.abs(Y) < fin_l, np.abs(X) < fin_w/2.)] = 1 
initial_grid[np.logical_and(np.abs(Y) < fin_w/2., np.abs(X) < fin_l)] = 1 

#obtain 3D Array 
contour_grid = skfmm.distance(initial_grid, dx=1e-2) 

perims = np.array([]) 
webs = np.linspace(0, r_f, n_pts/10) 
for web in webs: 
    web_ctr = pl.contour(X, Y, contour_grid, [web]) 
    no_error_here = p = web_ctr.collections[0].get_paths() 
    out_of_range_error = web_ctr.collections[0].get_paths()[0] 

내가 오류 "범위를 벗어 목록 지수"는 얻을, 그리고 그것은 밝혀 그 no_error_here 빈 목록입니다. 나는 조금 둘러 보았고, matplotlib가 이것을하는 이유를 발견 할 수 없다. 팁?

Answer 1

기본적으로 contour은 데이터 범위를 벗어나는 값으로 윤곽선을 그리면 빈 컬렉션을 반환합니다.

윤곽선 지정 값이 데이터 범위 내에 있는지 확인하거나 빈 콜렉션이 있는지 확인하십시오.

또한 현재는 지정된 값에 윤곽선이 하나만 있다고 가정합니다. 아마 나쁜 가정 일 겁니다. 단일 값을 컨투어링하기 때문에 collections[0]을 사용하려고하지만 그 윤곽에는 둘 이상의 경로가있을 수 있습니다.

---

결과를 플롯하지 않으려면 윤곽 값을 직접 추출하지 않는 이유는 무엇입니까? skimage.measure.find_contours은 자신이하는 일에 더 잘 어울립니다. 예를 들어

는 :

import numpy as np 
from skimage.measure import find_contours 

data = np.random.random((10,10)) 
contours = find_contours(data, 0.5) 
for xy in contours: 
    x, y = xy.T 
    perim = np.hypot(np.diff(x), np.diff(y)).sum()