【问题标题】:Get a contour (x,y) coordinates of a shape in an image with python使用python获取图像中形状的轮廓(x,y)坐标
【发布时间】:2016-02-09 03:45:49
【问题描述】:

我需要用python获取下图轮廓坐标(x,y)的矩阵。

我尝试使用 opencv canny 检测器并找到轮廓,但我得到了很多轮廓,但我不知道如何获得我想要的。

import numpy as np
from matplotlib import pyplot as plt
import cv2
#from skimage import measure, feature, io
#from skimage import img_as_ubyte

x1 = 330
xf = 690
y1 = 0
yf = 400

img = cv2.imread('test.tif')
img = img[y1:yf, x1:xf]
edge = cv2.Canny(img, 100, 200)

image, contours, hierarchy = cv2.findContours(edge, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_NONE)

我只需要一个包含轮廓 (x,y) 坐标的数组。我认为它在cv2.findContours() 的轮廓输出中,但我没有找到我想要的轮廓……

我也尝试了matplotlib.pyplot.contour 函数:

import cv2
import numpy as np
from matplotlib import pyplot as plt

img = cv2.imread('test.tif', 0) # read image
img = cv2.threshold(img, 0, 255, cv2.THRESH_OTSU)[1] # threshold image
img = cv2.medianBlur(img, 15)  # remove noise

# skeletonize 
size = np.size(img)  # get number of pixels
skel = np.zeros(img.shape, np.uint8) # create an array of zeros with the same shape as the image and 256 gray levels

element = cv2.getStructuringElement(cv2.MORPH_ELLIPSE, (3, 3)) # create a structurant element (cross)
done = False

while(not done):
    eroded = cv2.erode(img, element)
    temp = cv2.dilate(eroded, element)
    temp = cv2.subtract(img, temp)
    skel = cv2.bitwise_or(skel, temp)
    img = eroded.copy()
    zeros = size - cv2.countNonZero(img)
    if zeros == size:
        done = True

cs = plt.contour(skel, 1)
p = cs.collections[0].get_paths()[0]
v = p.vertices
x = v[:, 0]
y = v[:, 1]

但我只有闭合轮廓,而不是从图像左侧到右侧的开放轮廓。

非常感谢您的回答。

【问题讨论】:

  • “轮廓”是什么意思?你的预期输出是什么?具体一点。
  • 感谢您的评论@barny。我所说的等高线是指一条从左到右并穿过强度梯度最大的点的线。我需要矩阵中这条线的(x,y)坐标。

标签: python opencv numpy matplotlib


【解决方案1】:

您几乎找到了问题的答案。首先,边缘检测轮廓检测是有区别的。从根本上说,边缘检测会导致您(不正确地)所谓的“开放轮廓”(即边缘),而轮廓检测会导致您所谓的“闭合轮廓”(即轮廓)。

Canny 边缘检测 是一种流行的边缘检测算法。由于您想以数组的形式检测边缘,其中 (x,y) 坐标从图像的左侧到右侧,因此 Canny 边缘检测是一个好主意。

答案是edge,它不是想要的格式。

import numpy as np
import matplotlib.pyplot as plt
import cv2

img = cv2.imread('test.tif')
edge = cv2.Canny(img, 100, 200)

ans = []
for y in range(0, edge.shape[0]):
    for x in range(0, edge.shape[1]):
        if edge[y, x] != 0:
            ans = ans + [[x, y]]
ans = np.array(ans)

print(ans.shape)
print(ans[0:10, :])

数组ans(形状等于(n, 2))存储构成检测边缘的n 像素的(x,y)坐标。这就是您要寻找的结果。

这是我用白色绘制的n 像素的图像:

希望对你有帮助。

【讨论】:

  • 非常感谢巴蒂斯特的详细解释,这正是我想要的!
猜你喜欢
  • 2016-06-08
  • 1970-01-01
  • 1970-01-01
  • 2021-07-10
  • 1970-01-01
  • 2020-09-10
  • 1970-01-01
  • 2021-01-27
  • 1970-01-01
相关资源
最近更新 更多