Changeset 26764 for branches/eam_branches/20091201/extsrc/gpcsw/gpcsrc/fits/burntool/trailfit.c

Timestamp:

Feb 2, 2010, 7:06:18 PM (17 years ago)

Author:

watersc1

Message:

new burntool code

File:

• branches/eam_branches/20091201/extsrc/gpcsw/gpcsrc/fits/burntool/trailfit.c (modified) (3 diffs)

branches/eam_branches/20091201/extsrc/gpcsw/gpcsrc/fits/burntool/trailfit.c

@@ -137 +137 @@
    box->xfit = (int *)calloc(xe-xs+1, sizeof(int));
    box->yfit = (int *)calloc(xe-xs+1, sizeof(int));
+   if(box->zero == NULL || box->xfit == NULL || box->yfit == NULL) {
+      fprintf(stderr, "\rerror: failed to alloc box memory\n");
+      exit(-679);
+   }
 /* Some defaults */
    box->slope = 0.0;
@@ -226 +230 @@
 /* FIXME: sanity check fits */
    if(slope >= FIT_MAX_SLOPE || slope < FIT_MIN_SLOPE) {
+/* Check whether it's a significant trail (but with wrong slope) or just 
+ * noise.  In the former case negate box->nfit, 
+ * in the latter case set it to zero. */
+      linearrms(nfit, ybuf+y1, zbuf+y1, slope, zero, &trial);
       box->slope = slope;
-      box->nfit = 0;
+      if(trial > 2*rms) box->nfit = -box->nfit;
+      else box->nfit = 0;
       box->fiterr = FIT_SLOPE_ERROR;
       return(-1);
@@ -268 +277 @@
 /* FIXME: what's a really good criterion for negligible fit? */
 /* 100 pixels up fit is zsum or ~zsum/e */
-      if(zsum > NEGLIGIBLE_TRAIL*rms) {
+//      if(zsum > NEGLIGIBLE_TRAIL*rms) {
+      if(zsum > NEGLIGIBLE_TRAIL*rms ||
+/* 100113: but also evaluate at the end of the fit for stubby trails */
+         zbuf[y2]*zsum > NEGLIGIBLE_TRAIL*rms) {
 
 /* Ascertain the starting point of where the fit is good */