我有一张二进制图片。我想找到像素值= 1并将其标记为当前像素。然后,我想将其8个邻居像素值相加。如果当前像素的8个相邻像素值的总和= 1,则用标记标记该当前像素。二进制图像的某些部分如下:
0 0 0 0 0
0 1 0 0 0
0 0 1 1 0
0 0 0 0 1
0 0 0 0 0
我尝试了以下matlab代码,但它有一些错误(在此行 - > Sums = sum(currentPix, nOffsets);)
。我该如何解决?
Sums = 0;
S = size(BW,1);
nOffsets = [S, S+1, 1, -S+1, -S, -S-1, -1, S-1]'; %8-neighbors offsets
BW_Out = BW;
for row=1:S
for col=1:S
if BW(row,col),
break;
end
end
idx = sub2ind(size(BW),row,col);
neighbors = bsxfun(@plus, idx, nOffsets);
currentPix = find(BW==1); %if found 1, define it as current pixel
while ~isempty(currentPix)
% new current pixel list is set of neighbors of current list.
currentPix = bsxfun(@plus, currentPix, nOffsets);
currentPix = currentPix(:);
Sums = sum(currentPix, nOffsets); %error at this line
if (Sums==1) %if the sum of 8-neighbor values = 1, mark ROI
plot(currentPix,'r*','LineWidth',1);
end
% Remove from the current pixel list pixels that are already
currentPix(BW_Out(currentPix)) = [];
% Remove duplicates from the list.
currentPix = unique(currentPix);
end
end
答案 0 :(得分:3)
我认为你可以在一行中实现这一点(在定义内核之后)
I = [0 0 0 0 0
0 1 0 0 0
0 0 1 1 0
0 0 0 0 1
0 0 0 0 0];
K = [1 1 1;
1 0 1;
1 1 1;];
(conv2(I,K,'same')==1) & I
ans =
0 0 0 0 0
0 1 0 0 0
0 0 0 0 0
0 0 0 0 1
0 0 0 0 0
打破这个局面:
M = conv2(I,K, 'same'); %// convolving with this specific kernel sums up the 8 neighbours excluding the central element (i.e. the 0 in the middle)
(M==1) & I %// Only take results where there was a 1 in the original image.