In a method of image segmentation, a tree is derived from input data. A no-fill ordering of a matrix is produced, the matrix representing the tree. A system of linear equations is solved with respect to the matrix, obtaining a solution. The solution is used to define a family of segmentations. A segmentation is selected from the family of segmentations based on a metric of segment quality.