aboutsummaryrefslogtreecommitdiff
path: root/noao/onedspec/ecidentify/ecpeak.x
blob: f797fbac880088d38c694165cc576d3602225709 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
include	"ecidentify.h"

# EC_PEAK -- Find the peak value above continuum.

double procedure ec_peak (ec, pix)

pointer	ec			# ID pointer
double	pix			# Pixel position
double	peak			# Peak value

int	c, l, u

begin
	if (IS_INDEFD(pix))
	    return (INDEFD)

	c = nint (pix)
	l = max (1, nint (pix - EC_FWIDTH(ec)))
	u = min (EC_NPTS(ec), nint (pix + EC_FWIDTH(ec)))
	peak = IMDATA(ec,c) - (IMDATA(ec,l) +
	    IMDATA(ec,u)) / 2.

	return (peak)
end