buildParticleListAutoDP.m

function [indexVec,meanSegLength] = buildParticleListAutoDP(tol,dt,xPos,yPos,lengthMin)


temp = size(xPos);
nParticles = temp(2);
indexVec = [];

for ii = 1:nParticles
    t = 0:dt:(length(xPos{ii})-1)*dt;
    %%%%%%%%%%%% FIX SCALING FACTOR ABOVE!!!!!!
    [ps,ix] = dpsimplify([xPos{ii},yPos{ii},(t')/1E4],tol);
    nSegments = length(ix)-1;
    lengthT = 0;
    for jj = 1:nSegments
        lengthT = lengthT+sqrt((ps(jj+1,1)-ps(jj,1))^2+(ps(jj+1,2)-ps(jj,2))^2);
    end
    if lengthT>1E-6
	indexVec = [indexVec, ii];
    end
    meanSegLength(ii) = mean(diff(ix));
end

indexVec2 = [];
for ii = indexVec
    if length(xPos{ii})<200
    else
        indexVec2 = [indexVec2,ii];
    end
end
indexVec = indexVec2;