Changeset 42392 for trunk/extsrc/gpcsw/gpcsrc/fits/burntool/burntool.c

Timestamp:

Feb 14, 2023, 12:48:54 PM (3 years ago)

Author:

tdeboer

Message:

new version of the burntool programme

File:

• trunk/extsrc/gpcsw/gpcsrc/fits/burntool/burntool.c (modified) (7 diffs)

trunk/extsrc/gpcsw/gpcsrc/fits/burntool/burntool.c

@@ -89 +89 @@
    PSF_THRESH   =  5000;        /* Threshold for a star to be a PSF */
 
-   BMASK_GROW =  5.0;           /* Growth of burned boxes in size */
-   RMASK_GROW =  3.0;           /* Growth of burn/star in diameter */
+   BMASK_GROW =  3.0;           /* Growth of burned boxes in pixels */
+   RMASK_GROW =  2.0;           /* Growth of burn/star in diameter */
 
    MIN_PSF_SIZE  =    3;        /* Min box size for stamp selection */
@@ -97 +97 @@
    MAX_PSF_PER_CELL = 20;       /* Max number of PSF stars accepted per cell */
 
-   NEGLIGIBLE_TRAIL = 0.5;      /* Don't sweat less than this * sigma */
+   NEGLIGIBLE_TRAIL = 0.4;      /* Don't sweat less than this * sigma */
    EXPIRE_TRAIL_TIME = 2000;    /* Expire a persist after this [sec] */
 
@@ -564 +564 @@
 
 /* Does this cell look kosher? */
-         if(OTA[cell].rms*OTA[cell].rms > 
-            OTA[cell].sky/MIN_EADU + i*i*MAX_READ_NOISE*MAX_READ_NOISE) {
-            if(VERBOSE > 0) {
-               fprintf(stderr, "logonly: cell %d is unreasonably noisy, skipping...\n", cell);
-               fprintf(stderr, "logonly: bias = %d sky = %d rms = %d satfrac = %.2f\n", 
-                       OTA[cell].bias, OTA[cell].sky, OTA[cell].rms, OTA[cell].satfrac);
-            }
-            continue;
-         }
+//       if(OTA[cell].rms*OTA[cell].rms > 
+//          OTA[cell].sky/MIN_EADU + i*i*MAX_READ_NOISE*MAX_READ_NOISE) {
+//          if(VERBOSE > 0) {
+//             fprintf(stderr, "logonly: cell %d is unreasonably noisy, skipping...\n", cell);
+//             fprintf(stderr, "logonly: bias = %d sky = %d rms = %d satfrac = %.2f\n", 
+//                     OTA[cell].bias, OTA[cell].sky, OTA[cell].rms, OTA[cell].satfrac);
+//          }
+//          continue;
+//       }
 
          if(apply) {
@@ -601 +601 @@
             if(VERBOSE & VERB_NORM) burn_blab(OTA+cell);
 
-/* Fix up the streaks (don't bother if table only) */
-// Why the heck do I get to skip this if table only?  Makes no sense!
-//          if(!tableonly) {
-               persist_fix(naxis1-ovrscan1, naxis2-ovrscan2, naxis1, buf, 
+            persist_fix(naxis1-ovrscan1, naxis2-ovrscan2, naxis1, buf, 
                            OTA+cell);
 /* Tell us about it? */
-               if(VERBOSE & VERB_NORM) persist_blab(OTA+cell);
-//          }
+            if(VERBOSE & VERB_NORM) persist_blab(OTA+cell);
+
          }
 
@@ -667 +664 @@
 {
 /* Make some space for medians */
-   if(2*SKY_MARG*ny > nmedian_buf) {
+/*   if(2*SKY_MARG*ny > nmedian_buf) {
       if(median_buf != NULL) free(median_buf);
       if( (median_buf = (int *)calloc(2*SKY_MARG*ny, sizeof(int))) == NULL) {
@@ -675 +672 @@
 
       nmedian_buf = 2*SKY_MARG*ny;
+   }
+*/
+
+/* Make some space for cell copy */
+   if(nx*ny > nmedian_buf) {
+      if(median_buf != NULL) free(median_buf);
+      if( (median_buf = (int *)calloc(nx*ny, sizeof(DTYPE))) == NULL) {
+       fprintf(stderr, "\rerror: failed to alloc cell copy\n");
+       exit(-668);
+      }
+      nmedian_buf = nx*ny;
    }
 
@@ -728 +736 @@
 /* Get sky stats */
    for(k=n=nsat=0; k<nx*ny; k+=((617*nx)/1000)) {
+//   for(k=n=nsat=0; k<nx*ny; k+=((100*nx)/1000)) {
+//   for(k=n=nsat=0; k<nx*ny; k++) {
       i = k % nx;
       j = k / nx;