找出哪个 blob 具有像素位置 [x,y]？

Question

我有一个使用 bwlabel 标记的 blob 图像，我想找到一个具有像素位置 [x,y] 的 blob 并通过删除其余 blob 来显示它。 这是我写的代码，但它没有给出正确的答案，请修复这个

[y, x] = ginput(1);
x = round(x);
y = round(y); % here x and y is a location of blob i want to keep

    BW = bwlabel(newImgg,4) ; % labelled image contains several blobs
%     figure, imshow(BW, [])
    props = regionprops(logical(BW),'all');
    while(1)
        for k = 2:length(props)
            if ismember([x,y],props(k,1).PixelList) == [1, 1];
                keeperIndex = k;
                break
            end
        end
        break
    end

    keeperBlobsImage = ismember(BW, keeperIndex);
    keeperBlobsImage = imfill(keeperBlobsImage,'holes');
    figure, imshow(keeperBlobsImage,[])

谢谢，

戈皮

Answer 1

我目前没有 MATLAB 许可证，所以我无法在我的机器上测试它，我也有一段时间没有使用 MATLAB 语法了。这是一个想法：

From MATLAB's documentation、PixelList 是一个数组，其中每一行的格式为 [x,y,...]，具体取决于您的尺寸。

使用您的图片我假设PixelList 的格式为[x,y]

循环通过PixelList，跟踪要丢弃的索引。如果你测量了n 像素：

discardList = []
for i = 1:n
    if (PixelList(i) != [target_x,target_y]
        discardList=[discardList,i]
    end
end
newPixelList = PixelList
newPixelList(discardList) = []

再次重申，我已经有很长时间没有使用 MATLAB，因此对于语法（括号、循环和条件）中的任何问题，我深表歉意

编辑/更新：

根据 MA​​TLAB 的文档，它显示 bwlabel 仅用于 BW 图像。所以请确保你这样做，我猜。

此外，在regionprops 的输出中，您应该有WeightedCentroid。

从您的ginput 中，找到质心最近的区域。

我的建议是使用vision.BlobAnalysis System Object

[y,x] = ginput(1)
bA = vision.BlobAnalysis;
centroids = step(bA,BWImage);

使用文档确保关闭系统对象的所有“输出端口”，并保持质心输出端口打开。

d = 1e10;
d2 = 0;
dArr = [x,y;0,0]
cIndex=0;
for i = 1:length(centroids)
    dArr(2,:) = centroids(i,:);
    d2 = pdist(dArr);
    if (d2<d)
        d = d2;
        cIndex = i;
    end
end

变量cIndex 将包含您需要的blob 的索引。您可以运行 blob 分析并将其与其他部分隔离