找到中心像素周围半径内所有邻居的线性索引？ [复制]

Question

我有一个线性索引数组，对于每个线性索引，我想找到半径为5-pixels 的相邻像素的线性索引。我找到了以下代码，它为8-connected neighborhood 完成了这项工作。但是，如何实现它以找到120 neighbors 的线性索引，用于5-pixel 邻域。

%# target_array: array where pixels are marked
%# idx: linear index of a marked pixel
[M,N] = size(target_array)

neighbor_offsets=[-M-1 -M -M+1 1 M+1 M M-1 -1];

neighbors = bsxfun(@plus, idx, neighbor_offsets);

Answer 1

您提到的代码会找到像素周围的线性索引，只要该像素不太靠近目标数组的边界即可。

在你的情况下，我建议你逐个循环遍历像素，find 邻域：

[M,N] = size(target_array);

SE = strel('disk',5,inf);

%# the linear indices are stored in idxList
nIndices = length(idxList);

neighbors = cell(nIndices);

for ii = 1:nIndices
    idx = idxList(ii);
    bw = false(M,N);
    bw(idx) = true;
    bw = imdilate(bw,SE);

    %# mask the center
    bw(idx) = false;

    neighbors{ii} = find(bw);
end

如果你知道没有一个邻域重叠或接触，你可以简化上面的：

bw = false(M,N);
bw(idxList = true;
bw = imdilate(bw,SE);
bw(idxList) = false;
cc = bwconncomp(bw);
neighbors = cc.PixelIdxList;

Answer 2

您可以使用 meshgrid 和 sub2ind 来实现

假设一个名为T的目标数组和索引点（m，n）其中[m, n] = ind2sub(size(T), ind);

[X, Y] = meshgrid(m-5:m+5, n-5:n+5);
I = sub2ind(size(T), X(:), Y(:));

但如果你担心边缘（你应该是），那么使用 min 和 max 如下：

[M, N] = size(T);
[X, Y] = meshgrid(max(1,m-5):min(M,m+5), max(1,n-5):min(N,n+5));
I = sub2ind(size(T), X(:), Y(:));

另请注意，这将包括中心点，但很容易通过使用sub2ind(size(T), m, n);找到它的线性索引来删除它