Создание векторизованного массива из списка начальных/конечных индексов

У меня есть двухколоночная матрица M, которая содержит начальные/конечные индексы связки интервалов:

startInd   EndInd
1          3
6          10
12         12
15         16

Как я могу сгенерировать вектор всех индексов интервалов:

v = [1 2 3 6 7 8 9 10 12 15 16];

Я делаю выше, используя циклы, но мне интересно, есть ли более элегантное векторное решение?

v = [];
for i=1:size(M,1)
    v = [v M(i,1):M(i,2)];
end

Ответ 1

Здесь векторизованное решение, которое я хотел бы использовать для этой конкретной задачи, используя функцию cumsum:

v = zeros(1, max(endInd)+1);  % An array of zeroes
v(startInd) = 1;              % Place 1 at the starts of the intervals
v(endInd+1) = v(endInd+1)-1;  % Add -1 one index after the ends of the intervals
v = find(cumsum(v));          % Perform a cumulative sum and find the nonzero entries

Ответ 2

cell2mat(arrayfun(@colon,M(:,1)',M(:,2)','UniformOutput',false))

У меня нет IMFILL, но на моей машине этот метод быстрее, чем другие предложения, и я думаю, будет бить метод IMFILL из-за использования find.

Это можно сделать еще быстрее, если M настроен транспонированным (и мы корректируем третий и четвертый аргументы arrayfun).

Ответ 3

Вероятно, даже лучшее решение, которое я как-то не вижу, но здесь версия с IMFILL

startInd = [1,6,12,15];
endInd = [3,10,12,16];

%# create a logical vector with starts and ends set to true to prepare for imfill
tf = false(endInd(end),1);
tf([startInd,endInd]) = true;

%# fill at startInd+1 wherever startInd is not equal endInd
tf = imfill(tf,startInd(startInd~=endInd)'+1); %' SO formatting

%# use find to get the indices
v = find(tf)'  %' SO formatting

v =
     1     2     3     6     7     8     9    10    12    15    16

Ответ 4

Очень странное решение IMHO, создающее временные строки и использующее EVAL. Может быть также однострочным.

tmp = cellstr(strcat(num2str(M(:,1)),{':'},num2str(M(:,2)),{' '}));
v = eval(['[' cell2mat(tmp') ']']);

Я знаю, что он, вероятно, не будет работать на большой матрице. Просто для удовольствия.