繪製多邊形

我試圖繪製多邊形是這樣的：繪製多邊形

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

回答

相關問題