From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- noao/rv/rvidlines/idpeak.x | 23 +++++++++++++++++++++++ 1 file changed, 23 insertions(+) create mode 100644 noao/rv/rvidlines/idpeak.x (limited to 'noao/rv/rvidlines/idpeak.x') diff --git a/noao/rv/rvidlines/idpeak.x b/noao/rv/rvidlines/idpeak.x new file mode 100644 index 00000000..9cba49c4 --- /dev/null +++ b/noao/rv/rvidlines/idpeak.x @@ -0,0 +1,23 @@ +include "identify.h" + +# ID_PEAK -- Find the peak value above continuum. + +double procedure id_peak (id, pix) + +pointer id # 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 - ID_FWIDTH(id))) + u = min (ID_NPTS(id), nint (pix + ID_FWIDTH(id))) + peak = IMDATA(id,c) - (IMDATA(id,l) + IMDATA(id,u)) / 2. + + return (peak) +end -- cgit