A fast and robust segmentation model for piecewise smooth images is provided. Local statistics in an energy formulation are provided as a functional. The shape gradient of this new functional gives a contour evolution controlled by local averaging of image intensities inside and outside the contour. Fast computation is realized by expressing terms as the result of convolutions implemented via recursive filters. Results are similar to the general Mumford-Shah model but realized faster without having to solve a Poisson partial differential equation at each iteration. Examples are provided. A system to implement segmentation methods is also provided.