/*
 * skelaxis.c		Jim Piper		16-05-96
 *
 * compute chromosome medial axis from skeleton
 */

#include <stdio.h>
#include "../h/wstruct.h"

struct object *skelaxis(obj)
struct object *obj;
{
	register struct object *mpol, *skobj;
	struct object *skeleton();
	struct object *apolyfromskel();
	struct valuetable *newvaluetb();
	struct polygondomain *pdom;
	int numbranches;

	/*
	 * skeletonisation
	 */
	skobj = skeleton(obj,6);
	skobj->vdom = newvaluetb(skobj,1,0);
	setval(skobj,1);
	numbranches=afindtipandbranch(skobj);
	if (numbranches == 0) {
		mpol = apolyfromskel(skobj);
		pdom = (struct polygondomain *)mpol->idom;
		if (pdom->nvertices <= 3) {
			freeobj(mpol);
			mpol = NULL;
		} 
	} else
		mpol = NULL;
	freeobj(skobj);
	return(mpol);
}