Constructing the medial axis of a discrete level-set curve

In (mathematical) physics one can encounter problems with a fase change, the so called Stefan problems. Numerical methods to solve these problems often use a level-set function which indicates the distance to the fase boundary. Problems can arise in the medial axis points of the zero set. Therefore it can be usefull to calculate this medial axis...