The code is reasonably fast due to (optional) randomization and full code vectorization. However, as the algorithm needs to compute pairwise point distances, it can be quite memory intensive. If you get out of memory errors, either downsample the input image or somehow decrease the number of non-zero points in it. It can deal with big amount of noise but can have severe problem with occlusions (major axis end points need to be visible)