Summary: Two numerical approximation schemes for minimising the Mumford-Shah functional for unit vector fields are proposed, analyzed, and compared. The first uses a projection strategy, the second a penalization strategy to enforce the sphere constraint. Both schemes are then applied to the segmentation of colour images using the Chromaticity and Brightness colour model.