/*
 * chromaxis3.c		26-05-86		Jim Piper
 */

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

#define TIP 3
#define STEP 3
#define TIPSTRAIGHT 3

#define XPAR 6
#define YPAR 40
#define XEXCPAR	30


/*
 * Find chromosome axis of vertical chromosome.
 * First find "poor man's skeleton".  If this is both long and bent,
 * find axis from skeleton if skeleton is unbranched.
 */
struct object *chromaxis3(obj)
struct object *obj;
{
	register struct object *mpol, *skpol;
	struct object *pmsaxis(), *skelaxis();
	struct intervaldomain *intrp8(),*uspaxis8();
	struct polygondomain *pdom;
	register int tip, yr;

	/*
	 * find mid points polygon by new poor man's skeleton routine
	 */
	tip = (obj->idom->lastln - obj->idom->line1)/4;
	if (tip < 2)
		tip = 2;
	if (tip > TIP)
		tip = TIP;
	mpol = pmsaxis(obj,tip,STEP,TIPSTRAIGHT);
	/*
	 * features of this raw polygon for deciding whether to use skeleton:
	 * (1) height
	 * (2) x range
	 * (3) whether any points except end points lie outside x bounds
	 * (4) ratio of total x excursion to height
	 */
	yr = obj->idom->lastln - obj->idom->line1;
	if ((yr > YPAR && xrange(mpol) > XPAR) ||
			100*xexcursion(mpol)/yr > XEXCPAR ||
			axoutside(mpol,obj) ) {
		/*
		 * skeletonisation
		 */
		skpol = skelaxis(obj);
		if (skpol != NULL) {
			freeobj(mpol);
			mpol = skpol;
		}
	}
	pdom = (struct polygondomain *)mpol->idom;
	axissmooth(pdom,2);
	mpol->idom = uspaxis8(pdom);
	free(pdom);
#ifndef MIZAR
	pdom = (struct polygondomain *)mpol->idom;
	mpol->idom = intrp8(pdom);
	free(pdom);
#endif MIZAR
	return(mpol);
}