Changeset 42820 for branches/eam_branches/ipp-20230313/extsrc/gpcsw/gpcsrc/fits/burntool/stardetect.c

Timestamp:

May 8, 2025, 4:10:23 PM (15 months ago)

Author:

eugene

Message:

merge changes from trunk

Location:

branches/eam_branches/ipp-20230313

Files:

• . (modified) (1 prop)
• extsrc/gpcsw/gpcsrc/fits/burntool/stardetect.c (modified) (3 diffs)

branches/eam_branches/ipp-20230313/extsrc/gpcsw/gpcsrc/fits/burntool/stardetect.c

@@ -34 +34 @@
    starcut = STAR_THRESH/2 + cell->sky + cell->bias;
 
-   /*try a new value for thresh_lo which is 1 sigmas above the background, but with a safety limit of 1000 counts*/
+   /*try a new value for thresh_lo which is 10 sigmas above the background, but with a safety limit of 1000 counts*/
         trailthresh = cell->sky + (10*cell->rms);
         if(trailthresh <= 1000) trailthresh = 1000;
         trailthresh = trailthresh + cell->bias;
+        thresh_lo = trailthresh;
+
+   /*try a new value for burnthresh which is twice the bias above the sky background, but with a safety limit of 30000 counts*/
+        burnthresh = cell->sky + (2*cell->bias);
+        if(burnthresh <= BURN_THRESH) burnthresh = BURN_THRESH;
+        thresh_hi = burnthresh;
 
 /* Look at all the pixels which pass the burn threshold */
@@ -81 +87 @@
                 boxbuf[nbox].ey = k-1;
                 grow_box(nx, ny, NX, data, thresh_lo, boxbuf+nbox);
+                
            /* Fill in max and center info */
                 boxbuf[nbox].max = 0;
@@ -100 +107 @@
                              boxbuf[nbox].ex, boxbuf[nbox].ey,
                              boxbuf[nbox].max);
+                        fflush(stdout);
+                  }
+/* Mark the veto mask so as not to trigger on this one again */
+                  for(l=boxbuf[nbox].sy; l<=boxbuf[nbox].ey; l++) {
+                         for(k=boxbuf[nbox].sx; k<=boxbuf[nbox].ex; k++) {
+                          veto[k+l*NX] = boxbuf[nbox].max;
+                         }
+                  }
+                  xon = -1;
+                  continue;
+                }
+
+/* A box that is composed of a single pixel should not be considered a burn */
+                if((boxbuf[nbox].ex - boxbuf[nbox].sx) <= 2 && (boxbuf[nbox].ey - boxbuf[nbox].sy) <= 2) {
+                  if(VERBOSE & VERB_DETECT) {
+                        printf("Ditching box %d %d %d %d  for being too small\n",
+                             boxbuf[nbox].sx, boxbuf[nbox].sy,
+                             boxbuf[nbox].ex, boxbuf[nbox].ey);
                         fflush(stdout);
                   }