Pseudo-code for a medial axis transform algorithm based on Voronoi diagram.

gistfile1.txt

//Get point cloud S from polygon P
//by running a densify operation
var S = P.Densify ();

//Create Voronoi diagram V for
//point cloud S
var V = Voronoi (S);

//Get intersection C of Voronoi diagram V
//with polygon P
var C = V.Intersection (P);

//Retrieve medial axis by selecting edges
//in C which do not touch the perimeter
//of polygon P
var skeleton = from edge in C
  where (not (edge.Touches (P.Boundary))
  and edge.ContainedBy (F))
  select edge;