Tuesday, April 12, 2011

Detecting branched points in a skeleton

 
The image of a segmented chromosome (top left) is thresholded (top right) then used to calculate a skeleton with mahotas.thin() (lower left). The skeleton branched points  are detected by a hit-or-miss operator (lower right).
The python script depends on mahotas, pymorph, scipy.