Question

我正在使用Python 2.7。我必须在图片上定义一些感兴趣的区域（AoI）。基本上，我试图在图片的特定部分绘制椭圆（或更多），并获得其轮廓的坐标（x; y）。我想将这些坐标保存在文件中，以便稍后使用它们来查看我的数据是否在此区域内。

这是我的代码：

import matplotlib.pyplot as plt
import numpy as np
from matplotlib.patches import Ellipse, Circle
from matplotlib.path import Path

# Get an example image

img = imread('sposa.png')

# Create a figure. Equal aspect so circles look circular
fig,ax = plt.subplots(1)
ax.set_aspect('equal')

# Show the image
ax.imshow(img)
ax.set_xlim(0,1600)
ax.set_ylim(0,1200)
# Now, loop through coord arrays, and create a circle at each x,y pair
ellipse = Ellipse((1000, 400), width=400, height=100, edgecolor='white',facecolor='none',linewidth=2)

ax.add_patch(ellipse)
path = ellipse.get_path()

# Show the image
plt.show()

当我运行代码时，我得到了这个（这正是我想要的）：

然而，当我打印路径以便检查它时，我得到以下输出，（我猜）它与椭圆完全相关。

Path(array([[ 0.        , -1.        ],
   [ 0.2652031 , -1.        ],
   [ 0.51957987, -0.89463369],
   [ 0.70710678, -0.70710678],
   [ 0.89463369, -0.51957987],
   [ 1.        , -0.2652031 ],
   [ 1.        ,  0.        ],
   [ 1.        ,  0.2652031 ],
   [ 0.89463369,  0.51957987],
   [ 0.70710678,  0.70710678],
   [ 0.51957987,  0.89463369],
   [ 0.2652031 ,  1.        ],
   [ 0.        ,  1.        ],
   [-0.2652031 ,  1.        ],
   [-0.51957987,  0.89463369],
   [-0.70710678,  0.70710678],
   [-0.89463369,  0.51957987],
   [-1.        ,  0.2652031 ],
   [-1.        ,  0.        ],
   [-1.        , -0.2652031 ],
   [-0.89463369, -0.51957987],
   [-0.70710678, -0.70710678],
   [-0.51957987, -0.89463369],
   [-0.2652031 , -1.        ],
   [ 0.        , -1.        ],
   [ 0.        , -1.        ]]), array([ 1,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,
    4,  4,  4,  4,  4,  4,  4,  4, 79], dtype=uint8))

但是，我需要一个与图像像素相关的椭圆坐标列表（1600 X 1200）。我可能使用了错误的功能，或者图片和椭圆之间存在不匹配的东西。

我应该得到这样的东西（这是前一个实验的一个例子）：

[ Path(array([[ 1599.        ,   868.86791294],
   [ 1598.        ,   868.87197971],
   [ 1597.        ,   868.8801087 ],
   ..., 
   [ 1597.        ,   675.30378536],
   [ 1598.        ,   675.31373204],
   [ 1599.        ,   675.31870792]]), None)]
665

任何人都可以帮助我吗？先感谢您， [R

Answer 1

路径数组似乎是一个粗略的规范化圆圈 - 我忽略它

你已经有了椭圆信息 ellipse = Ellipse((1000, 400), width=400, height=100, ...)

如果你想明确地绘制一个高分辨率椭圆，我会根据Ellipse((1000, 400), width=400, height=100中的前几个数字做一个sin，cos paramaterized椭圆，它是中心点和轴长度

成员资格考试(x - x_0)^2/a^2 + (y - y_0)^2/b^2 <= 1

可能最适合a，b 是width=400, height=100

的1/2

当然你只需要在边界矩形

中测试像素索引

Answer 2

您应该使用Ellipse.get_path().vertices，但是它不在正确的坐标系中。要对其进行转换，请对其应用ellipse.get_patch_transform().transform。参见下面的工作示例

import matplotlib.pyplot as plt
from matplotlib.patches import Ellipse
from matplotlib.path import Path
from matplotlib.patches import PathPatch

img = plt.imread("image.jpg")

fig, ax = plt.subplots(1)
ax.set_aspect('equal')

ax.imshow(img)

# Create the base ellipse
ellipse = Ellipse((300, 300), width=400, height=100,
                  edgecolor='white', facecolor='none', linewidth=2)

# Get the path
path = ellipse.get_path()
# Get the list of path vertices
vertices = path.vertices.copy()
# Transform the vertices so that they have the correct coordinates
vertices = ellipse.get_patch_transform().transform(vertices)

# You can then save the vertices array to a file: csv, pickle... It's up to you

plt.show()

在图上绘制一个椭圆并获得coordinates_Python

2 个答案: