IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Oct 9, 2008, 2:39:01 PM (18 years ago)
Author:
beaumont
Message:

merged mainline into branch. resolved conflicts. added plots.

Location:
branches/cnb_branch_20080830/psastro/src
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • branches/cnb_branch_20080830/psastro/src

    • Property svn:ignore

      •  

        old new  
        1515psastroModel
        1616gpcModel
         17psastroModelFit

  • branches/cnb_branch_20080830/psastro/src/psastroDemoDump.c

    r16073 r20033  
    128128        if (!chip->process || !chip->file_exists) continue;
    129129       
     130        char *chipName = psMetadataLookupStr(NULL, chip->concepts, "CHIP.NAME");
     131
    130132        while ((cell = pmFPAviewNextCell (view, fpa, 1)) != NULL) {
    131133            psTrace ("psastro", 4, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
    … …  
    153155
    154156                    pmAstromObj *raw = rawstars->data[match->raw];
    155                     fprintf (f, "%f %f  %f %f  %f %f  %f %f  %f   |   ", 
    156                              DEG_RAD*raw->sky->r, DEG_RAD*raw->sky->d,
     157                    fprintf (f, "%s  %f %f  %f %f  %f %f  %f %f  %f   |   ", 
     158                             chipName, DEG_RAD*raw->sky->r, DEG_RAD*raw->sky->d,
    157159                             raw->TP->x, raw->TP->y,
    158160                             raw->FP->x, raw->FP->y,
    … …  
    193195}
    194196
    195 bool psastroDumpCorners (char *filename, pmFPA *fpa) {
     197bool psastroDumpCorners (char *filenameU, char *filenameD, pmFPA *fpa) {
    196198
    197199  // XXX test output of chip corners based on model
    198   FILE *f = fopen (filename, "w");
     200  FILE *fu = fopen (filenameU, "w");
     201  FILE *fd = fopen (filenameD, "w");
    199202
    200203  pmFPAview *view = pmFPAviewAlloc (0);
     204
     205  float fpaAngle = PM_DEG_RAD * atan2 (fpa->toTPA->y->coeff[1][0], fpa->toTPA->x->coeff[1][0]);
     206
     207  fprintf (fu, "# boresite: %f, %f @ %f\n", fpa->toSky->R*PS_DEG_RAD, fpa->toSky->D*PS_DEG_RAD, fpaAngle);
     208  fprintf (fd, "# boresite: %f, %f @ %f\n", fpa->toSky->R*PS_DEG_RAD, fpa->toSky->D*PS_DEG_RAD, fpaAngle);
    201209
    202210  pmChip *chip = NULL;
    … …  
    209217        psSphere ptSky;
    210218
     219        // UP 0,0
    211220        ptCP.x = region->x0; ptCP.y = region->y0;
    212221        psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
    213222        psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
    214223        psDeproject (&ptSky, &ptTP, fpa->toSky);
    215         fprintf (f, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
    216 
     224        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     225
     226        // DOWN 0,0
     227        psProject (&ptTP, &ptSky, fpa->toSky);
     228        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
     229        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
     230        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     231
     232        // UP 1,0
    217233        ptCP.x = region->x1; ptCP.y = region->y0;
    218234        psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
    219235        psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
    220236        psDeproject (&ptSky, &ptTP, fpa->toSky);
    221         fprintf (f, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
    222 
     237        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     238        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     239
     240        // DOWN 1,0
     241        psProject (&ptTP, &ptSky, fpa->toSky);
     242        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
     243        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
     244        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     245        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     246
     247        // UP 1,1
    223248        ptCP.x = region->x1; ptCP.y = region->y1;
    224249        psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
    225250        psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
    226251        psDeproject (&ptSky, &ptTP, fpa->toSky);
    227         fprintf (f, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
    228 
     252        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     253        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     254
     255        // DOWN 1,1
     256        psProject (&ptTP, &ptSky, fpa->toSky);
     257        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
     258        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
     259        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     260        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     261
     262        // UP 0,1
    229263        ptCP.x = region->x0; ptCP.y = region->y1;
    230264        psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
    231265        psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
    232266        psDeproject (&ptSky, &ptTP, fpa->toSky);
    233         fprintf (f, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
    234 
     267        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     268        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     269
     270        // DOWN 0,1
     271        psProject (&ptTP, &ptSky, fpa->toSky);
     272        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
     273        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
     274        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     275        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     276
     277        // UP 0,0
    235278        ptCP.x = region->x0; ptCP.y = region->y0;
    236279        psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
    237280        psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
    238281        psDeproject (&ptSky, &ptTP, fpa->toSky);
    239         fprintf (f, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     282        fprintf (fu, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
     283
     284        // DOWN 0,0
     285        psProject (&ptTP, &ptSky, fpa->toSky);
     286        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
     287        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
     288        fprintf (fd, "%10.6f %10.6f  %8.1f %8.1f  %8.1f %8.1f  %8.1f %8.1f\n", ptSky.r, ptSky.d, ptTP.x, ptTP.y, ptFP.x, ptFP.y, ptCP.x, ptCP.y);
    240289
    241290        psFree (region);
    242291  }
    243292
    244   fclose (f);
     293  fclose (fu);
     294  fclose (fd);
    245295  psFree (view);
    246296  return true;
Note: See TracChangeset for help on using the changeset viewer.