IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Feb 14, 2023, 12:59:09 PM (3 years ago)
Author:
tdeboer
Message:

burntool running updates

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/ppImage/src/ppImageBurntoolApply.c

    r35685 r42395  
    4949    psF32 satVal = psMetadataLookupF32(NULL, readout->parent->concepts, "CELL.SATURATION");
    5050
     51    /* The new burntool tables and method coming into play in 2022 has a slightly different format (an extra column). Need to differentiate */
     52    psS16 BURNTOOL_STATE_GOOD = psMetadataLookupS16(NULL, config->camera, "BURNTOOL.STATE.GOOD");
     53
    5154    for (int p = 0; p < burntoolData->npersist; p++) {
    5255        OBJBOX *box = &(burntoolData->persist[p]);
    … …  
    6164
    6265        y0 = box->y0;
    63         y1 = (box->up) ? ny-1 : 0;
    64         dy = (y0<y1) ? +1 : -1;
    65         if(box->up) {
    66             ys = box->sy;
     66        if (BURNTOOL_STATE_GOOD >= 15) {
     67          y1 = box->eyfit;
    6768        } else {
    68             ys = box->ey;
     69          y1 = (box->up) ? ny-1 : 0;
     70        }
     71        dy = (box->up) ? +1 : -1;
     72        if (BURNTOOL_STATE_GOOD >= 15) {
     73          ys = (box->up) ? box->sy : box->y0;
     74        } else {
     75          ys = (box->up) ? box->y1m : box->y1p;
    6976        }
    7077
    … …  
    9198        for(k=0; k<box->nfit; k++) {
    9299            i = box->xfit[k];
    93             for(j=box->yfit[k]; dy*j<=dy*y1; j+=dy) {
     100            for(j=ys; dy*j<=dy*y1; j+=dy) {
    94101                psF32 inputPixel = image->data.F32[j][i];
    95102                if (false && inputPixel >= satVal) {
    … …  
    150157            return NULL;
    151158        }
    152         if (persist_read(OTA, burntoolFileName, 1)) {
     159
     160        /* The new burntool tables and method coming into play in 2022 has a slightly different format (an extra column). Need to differentiate */
     161        psS16 BURNTOOL_STATE_GOOD = psMetadataLookupS16(NULL, config->camera, "BURNTOOL.STATE.GOOD");
     162        int oldfile = 1;
     163        if (BURNTOOL_STATE_GOOD >= 15) oldfile = 0;
     164
     165        if (persist_read(OTA, burntoolFileName, 1,oldfile)) {
    153166            psError(PS_ERR_UNKNOWN, "true", "failed to read burntool file");
    154167            return NULL;
    … …  
    174187  /* Set the maskValue from the recipes. */
    175188  psImageMaskType maskValue = options->burntoolMask;
     189
     190  /* The new burntool tables and method coming into play in 2022 has a slightly different format (an extra column). Need to differentiate */
     191  psS16 BURNTOOL_STATE_GOOD = psMetadataLookupS16(NULL, config->camera, "BURNTOOL.STATE.GOOD");
    176192
    177193  burntool_cell = (view->cell % 8) * 8 + (view->cell - (view->cell % 8)) / 8;
    … …  
    206222                  );
    207223#endif
     224
     225
     226        /* do separate masking strategy for old (14 and lower) and new burntool tables*/
     227        if (BURNTOOL_STATE_GOOD >= 15) {
     228         printf("DET: New burntool style %d\n",BURNTOOL_STATE_GOOD);
    208229          for (int i = box->sxfit; i<= box->exfit; i++) {
    209 
     230              if (box->up == 0) {
     231                  for (int j = box->eyfit; j <= box->y0; j++) {
     232                    #if PPIMAGE_BURNTOOL_DEBUG
     233                    psLogMsg("ppImageBurntoolMask", 4, "Noisy!: %d %d %d %d\n",
     234                         i,j,image->data.PS_TYPE_IMAGE_MASK_DATA[j][i],maskValue);
     235                    #endif
     236                      image->data.PS_TYPE_IMAGE_MASK_DATA[j][i] |= maskValue;
     237                  }
     238              } else {
     239                  for (int j = box->sy; j < box->eyfit ; j++) {
     240                    #if PPIMAGE_BURNTOOL_DEBUG
     241                      psLogMsg("ppImageBurntoolMask", 4, "Noisy!: %d %d %d %d\n",
     242                       i,j,image->data.PS_TYPE_IMAGE_MASK_DATA[j][i],maskValue);
     243                    #endif
     244                      image->data.PS_TYPE_IMAGE_MASK_DATA[j][i] |= maskValue;
     245                  }
     246              }
     247          }
     248
     249        } else {
     250         printf("DET: Old burntool style %d\n",BURNTOOL_STATE_GOOD);
     251          for (int i = box->sxfit; i<= box->exfit; i++) {
    210252              if (box->up == 0) {
    211253                  for (int j = 0; j <= box->y1p; j++) {
    212 #if PPIMAGE_BURNTOOL_DEBUG
     254                    #if PPIMAGE_BURNTOOL_DEBUG
    213255                    psLogMsg("ppImageBurntoolMask", 4, "Noisy!: %d %d %d %d\n",
    214256                         i,j,image->data.PS_TYPE_IMAGE_MASK_DATA[j][i],maskValue);
    215 #endif
     257                    #endif
    216258                      image->data.PS_TYPE_IMAGE_MASK_DATA[j][i] |= maskValue;
    217259                  }
    218260              } else {
    219261                  for (int j = box->y1m; j < image->numRows ; j++) {
    220 #if PPIMAGE_BURNTOOL_DEBUG
     262                    #if PPIMAGE_BURNTOOL_DEBUG
    221263                      psLogMsg("ppImageBurntoolMask", 4, "Noisy!: %d %d %d %d\n",
    222264                       i,j,image->data.PS_TYPE_IMAGE_MASK_DATA[j][i],maskValue);
    223 #endif
     265                    #endif
    224266                      image->data.PS_TYPE_IMAGE_MASK_DATA[j][i] |= maskValue;
    225                  }
    226              }
    227          }
     267                  }
     268              }
     269          }
     270        }
     271
     272
    228273      }
    229274  }
    … …  
    235280/****************************************************************/
    236281/* persist_read(): Read all the persistence trails from a file */
    237 STATIC int persist_read(CELL *cell, const char *infile, int apply)
     282STATIC int persist_read(CELL *cell, const char *infile, int apply, int oldfile)
    238283{
    239284   int i, k, nbox=0;
    … …  
    257302   while(fgets(line, 1024, fp) != NULL) {
    258303      if(line[0] == '#') continue;
    259       sscanf(line, "%d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %lf %d %d %d %d",
     304
     305      /* read-in of new burntool tables (including the eyfit column)*/
     306      if(!oldfile) {
     307        sscanf(line, "%d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %lf %d %d %d %d %d",
     308             &boxbuf[nbox].cell, &boxbuf[nbox].time,
     309             &boxbuf[nbox].cx, &boxbuf[nbox].cy,
     310             &boxbuf[nbox].max, &boxbuf[nbox].y0,
     311             &boxbuf[nbox].sx, &boxbuf[nbox].sy,
     312             &boxbuf[nbox].ex, &boxbuf[nbox].ey,
     313             &boxbuf[nbox].y0m, &boxbuf[nbox].y0p,
     314             &boxbuf[nbox].y1m, &boxbuf[nbox].y1p,
     315             &boxbuf[nbox].x0m, &boxbuf[nbox].x0p,
     316             &boxbuf[nbox].x1m, &boxbuf[nbox].x1p,
     317             &boxbuf[nbox].func, &boxbuf[nbox].up,
     318             &boxbuf[nbox].slope, &boxbuf[nbox].nfit,
     319             &boxbuf[nbox].sxfit, &boxbuf[nbox].exfit,
     320             &boxbuf[nbox].fiterr,&boxbuf[nbox].eyfit);
     321
     322      } else {
     323        sscanf(line, "%d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %d %lf %d %d %d %d",
    260324             &boxbuf[nbox].cell, &boxbuf[nbox].time,
    261325             &boxbuf[nbox].cx, &boxbuf[nbox].cy,
    … …  
    271335             &boxbuf[nbox].sxfit, &boxbuf[nbox].exfit,
    272336             &boxbuf[nbox].fiterr);
     337      }
     338
    273339      if(boxbuf[nbox].nfit > 0) {
    274340         boxbuf[nbox].zero = (double *)calloc(boxbuf[nbox].nfit, sizeof(double));
    … …  
    431497               box[j].sxfit = MIN(box[j].sxfit, box[k].sxfit);
    432498               box[j].exfit = MAX(box[j].exfit, box[k].exfit);
     499               box[j].eyfit = MIN(box[j].eyfit, box[k].eyfit);
    433500
    434501               box[k].exfit = box[k].sxfit - 1;
    … …  
    459526/****************************************************************/
    460527/* persist_write(): Write all the persistence data for the next image */
    461 STATIC int persist_write(CELL *cell, const char *outfile)
     528STATIC int persist_write(CELL *cell, const char *outfile, int oldfile)
    462529{
    463530   int i, k, j, err=0;
    … …  
    474541   fprintf(fp,  "# Cell: %d  sky= %d   rms= %d   bias= %d\n",
    475542           cell[0].cell, cell[0].sky, cell[0].rms, cell[0].bias);
    476    fprintf(fp, "#Cell time    cx  cy  max   y0   sx  sy   ex  ey  y0m y0p y1m y1p x0m x0p x1m x1p F up    slope nfit sxf exf\n");
     543
     544   /*write new burntool tables (including the eyfit column)*/
     545   if(!oldfile) {
     546     fprintf(fp, "#Cell time    cx  cy  max   y0   sx  sy   ex  ey  y0m y0p y1m y1p x0m x0p x1m x1p F up    slope nfit sxf exf fiterr eyf\n");
     547   } else {
     548     fprintf(fp, "#Cell time    cx  cy  max   y0   sx  sy   ex  ey  y0m y0p y1m y1p x0m x0p x1m x1p F up    slope nfit sxf exf fiterr\n");
     549   }
    477550
    478551   for(j=0; j<MAXCELL; j++) {
    … …  
    494567            if(cell[j].persist[k].nfit <= 0) continue;
    495568         }
    496          fprintf(fp, "%3d %7d  %3d %3d %5d %3d  %3d %3d  %3d %3d  %3d %3d %3d %3d %3d %3d %3d %3d  %1d %1d %9.6f %3d %3d %3d %d\n",
     569 
     570         /*write new burntool tables (including the eyfit column)*/
     571         if(!oldfile) {
     572           fprintf(fp, "%3d %7d  %3d %3d %5d %3d  %3d %3d  %3d %3d  %3d %3d %3d %3d %3d %3d %3d %3d  %1d %1d %9.6f %3d %3d %3d %d %3d\n",
     573                 j, cell[j].persist[k].time,
     574                 cell[j].persist[k].cx, cell[j].persist[k].cy,
     575                 cell[j].persist[k].max, cell[j].persist[k].y0,
     576                 cell[j].persist[k].sx, cell[j].persist[k].sy,
     577                 cell[j].persist[k].ex, cell[j].persist[k].ey,
     578                 cell[j].persist[k].y0m, cell[j].persist[k].y0p,
     579                 cell[j].persist[k].y1m, cell[j].persist[k].y1p,
     580                 cell[j].persist[k].x0m, cell[j].persist[k].x0p,
     581                 cell[j].persist[k].x1m, cell[j].persist[k].x1p,
     582                 cell[j].persist[k].func, cell[j].persist[k].up,
     583                 cell[j].persist[k].slope, cell[j].persist[k].nfit,
     584                 cell[j].persist[k].sxfit, cell[j].persist[k].exfit,
     585                 cell[j].persist[k].fiterr, cell[j].persist[k].eyfit);
     586         } else {
     587           fprintf(fp, "%3d %7d  %3d %3d %5d %3d  %3d %3d  %3d %3d  %3d %3d %3d %3d %3d %3d %3d %3d  %1d %1d %9.6f %3d %3d %3d %d\n",
    497588                 j, cell[j].persist[k].time,
    498589                 cell[j].persist[k].cx, cell[j].persist[k].cy,
    … …  
    508599                 cell[j].persist[k].sxfit, cell[j].persist[k].exfit,
    509600                 cell[j].persist[k].fiterr);
     601         }
     602
    510603         for(i=0; i<cell[j].persist[k].nfit; i++) {
    511604            fprintf(fp, "%3d %3d %8.4f\n", cell[j].persist[k].xfit[i],
    … …  
    531624
    532625         i = (cell[j].burn[k].ex - cell[j].burn[k].sx + 1) / 2;
    533          fprintf(fp, "%3d %7d  %3d %3d %5d %3d  %3d %3d  %3d %3d  %3d %3d %3d %3d %3d %3d %3d %3d  %1d %1d %9.6f %3d %3d %3d %d\n",
     626
     627         if(!oldfile) {
     628           fprintf(fp, "%3d %7d  %3d %3d %5d %3d  %3d %3d  %3d %3d  %3d %3d %3d %3d %3d %3d %3d %3d  %1d %1d %9.6f %3d %3d %3d %d %3d\n",
     629                 j, cell[j].burn[k].time,
     630                 cell[j].burn[k].cx, cell[j].burn[k].cy,
     631                 cell[j].burn[k].max, cell[j].burn[k].y0,
     632                 cell[j].burn[k].sx, cell[j].burn[k].sy,
     633                 cell[j].burn[k].ex, cell[j].burn[k].ey,
     634                 cell[j].burn[k].y0m, cell[j].burn[k].y0p,
     635                 cell[j].burn[k].y1m, cell[j].burn[k].y1p,
     636                 cell[j].burn[k].x0m, cell[j].burn[k].x0p,
     637                 cell[j].burn[k].x1m, cell[j].burn[k].x1p,
     638                 cell[j].burn[k].func, cell[j].burn[k].up,
     639                 cell[j].burn[k].slope, cell[j].burn[k].nfit,
     640                 cell[j].burn[k].sxfit, cell[j].burn[k].exfit,
     641                 cell[j].burn[k].fiterr, cell[j].burn[k].eyfit);
     642         } else {
     643           fprintf(fp, "%3d %7d  %3d %3d %5d %3d  %3d %3d  %3d %3d  %3d %3d %3d %3d %3d %3d %3d %3d  %1d %1d %9.6f %3d %3d %3d %d\n",
    534644                 j, cell[j].burn[k].time,
    535645                 cell[j].burn[k].cx, cell[j].burn[k].cy,
    … …  
    545655                 cell[j].burn[k].sxfit, cell[j].burn[k].exfit,
    546656                 cell[j].burn[k].fiterr);
     657         }
     658
    547659         for(i=0; i<cell[j].burn[k].nfit; i++) {
    548660            fprintf(fp, "%3d %3d %8.4f\n", cell[j].burn[k].xfit[i],
Note: See TracChangeset for help on using the changeset viewer.