how to count the number of zeros between 2 one's
6 views (last 30 days)
Show older comments
i am using the code
num = [0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0];
num = num(find(num==1, 1, 'first'):find(num==1, 1, 'last'));
% Measure the lengths of the stretches/runs of zeros.
measurements = regionprops(logical(~num), 'Area');
zerospan = [measurements.Area]
a = max(zerospan)
this is giving me answer as 9 which is correct for middle part. but i would also count the end around number of zeros. that is zeros starting from second last position and ending at fourth position. which should give me answer as 6.
thanks monica
0 Comments
Accepted Answer
Andrei Bobrov
on 24 Apr 2016
Edited: Andrei Bobrov
on 24 Apr 2016
num = [0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0]
measurements = regionprops(logical(~num), 'Area')
m = [measurements.Area];
out = [sum(m([1,end])), max(m(2:end-1))]
1 Comment
Image Analyst
on 24 Apr 2016
"which should give me answer as 6." Just be aware that this answer gives [6,9].
More Answers (2)
Azzi Abdelmalek
on 24 Apr 2016
Edited: Azzi Abdelmalek
on 24 Apr 2016
nu=~num;
s=cumsum(num)+~num(1);
ii=strfind(num,[1 0]);
jj=zeros(size(num));
jj(ii)=1;
s=s+cumsum(jj)
f=accumarray(s',(1:numel(s))',[],@(x) sum(nu(x)));
out=[ nonzeros(f(2:end-1))' sum([f(1) f(end)])]
2 Comments
Image Analyst
on 24 Apr 2016
You almost had it. You just needed to ask for the 'PixelIdxList' when you called regionprops. This gives you the index of every element in each grouping of 0's.
num = [0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0];
% Measure the lengths of the stretches/runs of zeros.
measurements = regionprops(logical(~num), 'Area', 'PixelIdxList');
zerospan = [measurements.Area]
% Get index of the largest run in the lengths array:
[longestRun, indexOfLongestRun] = max(zerospan)
% Get the index in the original num array:
indexOfStartOfLongestRun = measurements(indexOfLongestRun).PixelIdxList(1)
You will see that indexOfLongestRun = 9, which is where that run starts at.
Image Analyst
on 24 Apr 2016
Just add the first and last area:
result = zerospan(1) + zerospan(end);
0 Comments
See Also
Categories
Find more on Logical in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!