include	<smw.h>
include	"identify.h"

# ID_CENTER -- Locate the center of a feature.

double procedure id_center (id, x, width, type)

pointer	id				# ID pointer
double	x				# Initial guess
real	width				# Feature width
int	type				# Feature type

int	np1
real	value

real	center1dx()
double	smw_c1trand()

begin
	np1 = NP1(ID_SH(id)) - 1
	value = smw_c1trand (ID_PL(id), x) - np1
 	value = center1dx (value, IMDATA(id,1), ID_NPTS(id),
	    width, type, ID_CRADIUS(id), ID_THRESHOLD(id))
	if (IS_INDEF(value))
	    return (INDEFD)
	else
	    return (smw_c1trand (ID_LP(id), double(value+np1)))
end