This is great code! It's very easy to read and taught me a lot about mean shift.
I found that the runtime was decreased on my machine and data by replacing the lines:
sqDistToAll = sum((repmat(myMean,1,numPts) - dataPts).^2);
and
beenVisitedFlag(myMembers) = 1;
with:
sqDistToAll = sum(bsxfun(@minus,myMean,dataPts).^2);
and
beenVisitedFlag(inInds) = 1;
respectively.
Thanks for the code!
Comment only