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/twodspec/apextract/apcenter.x | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) create mode 100644 noao/twodspec/apextract/apcenter.x (limited to 'noao/twodspec/apextract/apcenter.x') diff --git a/noao/twodspec/apextract/apcenter.x b/noao/twodspec/apextract/apcenter.x new file mode 100644 index 00000000..88f089d1 --- /dev/null +++ b/noao/twodspec/apextract/apcenter.x @@ -0,0 +1,26 @@ +include + +# AP_CENTER -- Locate the center of an emission profile. This is done +# using the CENTER1D algorithm. The procedure gets the centering +# parameters using CL queries. If the center is not found because of the +# RADIUS or THRESHOLD centering criteria then INDEF is returned. + +real procedure ap_center (x, data, npts) + +real x # Initial guess +real data[npts] # Data +int npts # Number of data points + +real width # Centering width +real radius # Centering radius +real threshold # Detection threshold + +real apgetr(), center1d() + +begin + width = apgetr ("width") + radius = apgetr ("radius") + threshold = apgetr ("threshold") + + return (center1d (x, data, npts, width, EMISSION, radius, threshold)) +end -- cgit