しきい値を超えてデータを抽出する

Question

したがって、行1がタイムスタンプである2000 x 4行列がある場合、行2-4はさまざまな変数のデータです。

1つの変数（列2など）のデータを抽出するのは特定のしきい値を超えていますが、値が元に戻ってから前の100個のサンプルと次の100個のサンプルを抽出します再び閾値。

overThresh=sampleMatrix(:,columnOfInterest)>threshold; 

これが閾値よりも大きいいる「真」を示す要素を使用してマトリックスと同じ長さの論理的な列ベクトルを作成します値が閾値を超えている場合

Answer 1

を検出することを含みます。

今100個のサンプルでこれらの領域を拡大するために、私は画像処理のトリックを使用します：画像膨張：

SE = strel('arbitrary', ones(100,1)) 
extractVector=imdilate(overThresh,SE) 

今extractVectorは（それらの隣にしきい値+ 100を介して）したいサンプルのために真であります。

interestingSamples=sampleMatrix(extractVector,:) 

これは、あなたのすべての列の値を持つ行列を与えるが、興味のある行のみ

編集：。 あなたが使用して、しきい値違反の各イベントにラベル（整数）を割り当てることができますが：

overthreshLabeled=bwlabel(overthresh); 

ここで、各領域は「真」の代わりに数字で構成されます。各番号は隣接するイベントの特定の領域に対応し、閾値を超えている。 次に、リージョン番号をループして値を抽出し、セルにすべて入れることができます（長さが異なるため、セルにする必要があります）。

ROIs{1,max(labeledThreshold)}=[]; % preallocation 

for i=1:max(labeledThreshold) 
    region=labeledThreshold==i; 
    extractVector=imdilate(overThresh,SE) % Note: set SE before loop 
    ROIs{1,i}=sampleMatrix(extractVector,:) 
end