Identifying time for prolong changes in a sequence

binarySequence = [0 1 0 0 0 1 1 0 0 0 0 0 1 1 1 1 1 0 0 0 0 1 0 0 0 1 1 0 0 1];
props = regionprops(logical(binarySequence), 'Area');
allLengths = [props.Area]
% Count number of regions longer than some number:
maxLength = max(allLengths)
for k = 1 : maxLength
	longerThanN = sum(allLengths == k);
	fprintf('There are %d sequences of %d.\n', longerThanN, k);
end

You get:

There are 3 sequences of 1.

There are 2 sequences of 2.

There are 0 sequences of 3.

There are 0 sequences of 4.

There are 1 sequences of 5.

10 个评论
显示 8更早的评论隐藏 8更早的评论

Image Analyst 2021-7-11

在 MATLAB Online 中打开

Actually you can just ask regionprops() for the locations:

binarySequence = [0 1 0 0 0 1 1 0 0 0 0 0 1 1 1 1 1 0 0 0 0 1 0 0 0 1 1 0 0 1];
% Label the data
props = regionprops(logical(binarySequence), 'Area', 'PixelIdxList');
allLengths = [props.Area]
% Count number of regions longer than some number:
maxLength = max(allLengths)
for k = 1 : maxLength
	% Get the labels for this length
	labels = allLengths == k;
	longerThanN = sum(labels);
	fprintf('There are %d sequences of %d.\n', longerThanN, k);
	% Get 1's wherever the run has a length of k.
	locations = zeros(size(binarySequence));
	for k2 = 1 : length(labels)
		if labels(k2)
			locations(props(k2).PixelIdxList) = 1;
		end
	end
	% Report to the command window:
	locations
end

If you want locations in a 2-D array, you can do this:

binarySequence = [0 1 0 0 0 1 1 0 0 0 0 0 1 1 1 1 1 0 0 0 0 1 0 0 0 1 1 0 0 1];
% Label the data
props = regionprops(logical(binarySequence), 'Area', 'PixelIdxList');
allLengths = [props.Area]
% Count number of regions longer than some number:
maxLength = max(allLengths)
locations = zeros(maxLength, length(binarySequence));
for k = 1 : maxLength
	% Get the labels for this length
	labels = allLengths == k;
	longerThanN = sum(labels);
	fprintf('There are %d sequences of %d.\n', longerThanN, k);
	% Get 1's wherever the run has a length of k.
	for k2 = 1 : length(labels)
		if labels(k2)
			locations(k, props(k2).PixelIdxList) = 1;
		end
	end
end
% Report to the command window:
locations

and you get

locations =

Columns 1 through 23

0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0

0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0

Columns 24 through 30

0 0 0 0 0 0 1

0 0 1 1 0 0 0

0 0 0 0 0 0 0

Note the length is the row number. So row 2 has 1's wherever those locations are in a run of 2 digits long, namely 1's in columns 6 & 7, and 26 & 27.

Tino 2021-7-12

Hi Image analyst.

If I wanted only onces to show (and there group occurences in the sequence) in the result instead of both 1 and zero

for instance if I have 1111100001111000000111111

I will like the results to be

There are 1 sequences of 4 = second occurences

There are 2 sequences of 5 = first and third occurences

I will appreciate it if you answer the question for me as you ve done in the past

Best wishes

Image Analyst 2021-7-12

在 MATLAB Online 中打开

Try this:

binarySequence = [0 1 0 0 0 1 1 0 0 0 0 0 1 1 1 1 1 0 0 0 0 1 0 0 0 1 1 0 0 1]
% Measure the data
props = regionprops(logical(binarySequence), 'Area', 'PixelIdxList');
allLengths = [props.Area]
% Count number of regions longer than some number:
maxLength = max(allLengths)
for k = 1 : maxLength
	% Get the labels for this length
	labels = allLengths == k;
	longerThanN = sum(labels);
	fprintf('There are %d sequences of %d.\n', longerThanN, k);
	indexes{k} = find(labels);
	for k2 = 1 : length(indexes{k})
		fprintf('    Sequence %d\n', indexes{k}(k2));
	end
	% Get 1's wherever the run has a length of k.
	locations = zeros(size(binarySequence));
	for k2 = 1 : length(labels)
		if labels(k2)
			locations(props(k2).PixelIdxList) = 1;
		end
	end
	% Report to the command window:
	%locations
end

There are 3 sequences of 1.

Sequence 1

Sequence 4

Sequence 6

There are 2 sequences of 2.

Sequence 2

Sequence 5

There are 0 sequences of 3.

There are 0 sequences of 4.

There are 1 sequences of 5.

Sequence 3

请先登录，再进行评论。

Identifying time for prolong changes in a sequence

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

10 个评论
显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Identifying time for prolong changes in a sequence

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

10 个评论 显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

10 个评论
显示 8更早的评论隐藏 8更早的评论