绘制多边形

我试图绘制多边形是这样的：绘制多边形

In [1]: canvas = numpy.zeros((12, 12), dtype=int) 

In [2]: mahotas.polygon.fill_polygon(
    ...: [(1, 1), (1, 10), (10, 10), (10, 1)], 
    ...: canvas) 

In [3]: canvas 
Out[3]: 
array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
     [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]])

我希望下面的输出，但：

array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0], 
     [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]])

为什么[(10,2):(10:10)]仍为零？有没有另一种方法来绘制一个填充多边形到数组？

来源

2011-04-07 Attila O.

这是一个奇怪的结果。我发现如果你颠倒点的顺序，它会画出完整的图。换句话说：

# this is broken 
pts = [(1, 1), (1, 10), (10, 10), (10, 1)] 
# this works 
pts = [(1, 1), (10, 1), (10, 10), (1, 10)]

这是一个测试程序：

import numpy 
import mahotas.polygon 

def run(n, reverse=0): 
    canvas = numpy.zeros((n, n), dtype=int) 
    lim = n-2 
    print '\n%d x %d, lim=%d reverse=%d' % (n, n, lim, reverse) 
    pts = [(1, 1), (1, lim), (lim, lim), (lim, 1), (1, 1)] 
    if reverse: 
     pts.reverse() 
    mahotas.polygon.fill_polygon(pts, canvas) 
    return canvas 

for rev in (0, 1): 
    for n in range(3, 14): 
     print run(n, rev)

例子：

6 x 6, lim=4 reverse=0 
[[0 0 0 0 0 0] 
[0 1 0 0 1 0] 
[0 1 1 1 1 0] 
[0 1 1 1 1 0] 
[0 1 0 0 0 0] 
[0 0 0 0 0 0]] 

6 x 6, lim=4 reverse=1 
[[0 0 0 0 0 0] 
[0 1 1 1 1 0] 
[0 1 1 1 1 0] 
[0 1 1 1 1 0] 
[0 1 1 1 1 0] 
[0 0 0 0 0 0]]

来源

2011-04-07 22:36:32 samplebias

令人印象深刻！唯一的是，我不知道提前订购。我可能做的是以两种方式绘制多边形，然后将其中包含更多“1”的版本作为正确结果。感谢你的回答！ – 2011-04-08 14:26:57

如果多边形总是矩形，那么我们只需要两点：

import numpy 
canvas = numpy.zeros((12, 12), dtype=int) 
points = [(1, 1), (1, 10), (10, 10), (10, 1)] 
start_pt, end_pt = min(points), max(points) 
canvas[start_pt[1]:end_pt[1]+1, start_pt[0]:end_pt[0]+1] = 1

来源

2011-04-08 23:13:37 Benjamin

这只是一个例子，我的多边形是来自OpenStreetMap数据库的建筑物（通常不是矩形）。但谢谢你的答案。 – 2011-04-09 19:18:05

我认为逆转是因为numpy引用数组坐标为y，x（行，列），其中大多数其他程序将坐标指定为x，y。

来源

2013-06-03 08:27:05

回答

相关问题