11

M一連の間隔の開始/終了インデックスを含む2列の行列があります。

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
4

4 に答える 4

9

これは、関数を使用して、この特定の問題に使用するのが好きなベクトル化されたソリューションです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
于 2010-05-11T04:04:22.617 に答える
6 
cell2mat(arrayfun(@colon,M(:,1)',M(:,2)','UniformOutput',false))

私はIMFILLを持っていませんが、私のマシンでは、このメソッドは他の提案よりも高速であり、findを使用するためIMFILLメソッドよりも優れていると思います。

Mを転置して設定すると、さらに高速になります（arrayfunの3番目と4番目の引数を調整します）。

于 2011-01-31T20:51:35.217 に答える
1

私がどういうわけか見ていませんが、おそらくもっと良い解決策がありますが、これは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
于 2010-05-11T01:11:59.347 に答える
0

非常に奇妙な解決策IMHOが一時的な文字列を作成し、EVALを使用します。ワンライナーにすることもできます。

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

私はそれがおそらく大きなマトリックスでは機能しないことを知っています。ただ楽しみのために。

于 2010-05-11T01:19:56.393 に答える