У меня есть набор точек вдоль овала. Как создать заполненную двоичную маску

Question

Я пытаюсь получить заполненную двоичную маску контура этого изображения.

Я взглянул на этот вопрос SciPy Create 2D Polygon Mask; однако, похоже, мне не нравится мой набор данных.

import numpy as np 
from matplotlib.nxutils import points_inside_poly 

nx, ny = 10, 10 
poly_verts = [(1,1), (5,1), (5,9),(3,2),(1,1)] 

# Create vertex coordinates for each grid cell... 
# (<0,0> is at the top left of the grid in this system) 
x, y = np.meshgrid(np.arange(nx), np.arange(ny)) 
x, y = x.flatten(), y.flatten() 

points = np.vstack((x,y)).T 

grid = points_inside_poly(points, poly_verts) 
grid = grid.reshape((ny,nx)) 

print grid 

Интересно, есть ли другой способ, что я могу попытаться вернуть двоичную маску или кто-то, чтобы объяснить ограничения points_inside_poly

, потому что, кажется, в конечном итоге что-то вроде этого

Answer 1

I не уверен, что вы черчения в конце, но ваш пример работает для меня:

import numpy as np 
import matplotlib.pyplot as plt 
from matplotlib.nxutils import points_inside_poly 
from itertools import product, compress 

pv = [(1,1),(5,1),(5,9),(3,2),(1,1)] 

x, y = np.meshgrid(np.arange(10),np.arange(10)) 
x, y = x.flatten(), y.flatten() 

xy = np.vstack((x,y)).T 

grid = points_inside_poly(xy,pv) 

xv, yv = zip(*pv) 
xp, yp = zip(*compress(xy,grid)) 

plt.plot(xp,yp,'o',color='red',label='points') 
plt.plot(xv,yv,'o',color='blue',label='vertices') 
plt.xlim((0,10)) 
plt.ylim((0,10)) 
plt.legend() 
plt.show()