Changeset 19511

Timestamp:

Sep 11, 2008, 3:08:04 PM (18 years ago)

Author:

eugene

Message:

adding up and down coords to corner dump

File:

• trunk/psastro/src/psastroDemoDump.c (modified) (2 diffs)

trunk/psastro/src/psastroDemoDump.c

@@ -195 +195 @@
 }
 
-bool psastroDumpCorners (char *filename, pmFPA *fpa) {
+bool psastroDumpCorners (char *filenameU, char *filenameD, pmFPA *fpa) {
 
   // XXX test output of chip corners based on model
-  FILE *f = fopen (filename, "w");
+  FILE *fu = fopen (filenameU, "w");
+  FILE *fd = fopen (filenameD, "w");
 
   pmFPAview *view = pmFPAviewAlloc (0);
+
+  float fpaAngle = PM_DEG_RAD * atan2 (fpa->toTPA->y->coeff[1][0], fpa->toTPA->x->coeff[1][0]);
+
+  fprintf (fu, "# boresite: %f, %f @ %f\n", fpa->toSky->R*PS_DEG_RAD, fpa->toSky->D*PS_DEG_RAD, fpaAngle);
+  fprintf (fd, "# boresite: %f, %f @ %f\n", fpa->toSky->R*PS_DEG_RAD, fpa->toSky->D*PS_DEG_RAD, fpaAngle);
 
   pmChip *chip = NULL;
@@ -211 +217 @@
         psSphere ptSky;
 
+        // UP 0,0
         ptCP.x = region->x0; ptCP.y = region->y0;
         psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
         psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
         psDeproject (&ptSky, &ptTP, fpa->toSky);
-        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);
-
+        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);
+
+        // DOWN 0,0
+        psProject (&ptTP, &ptSky, fpa->toSky);
+        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
+        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
+        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);
+
+        // UP 1,0
         ptCP.x = region->x1; ptCP.y = region->y0;
         psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
         psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
         psDeproject (&ptSky, &ptTP, fpa->toSky);
-        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);
-
+        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);
+        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);
+
+        // DOWN 1,0
+        psProject (&ptTP, &ptSky, fpa->toSky);
+        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
+        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
+        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);
+        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);
+
+        // UP 1,1
         ptCP.x = region->x1; ptCP.y = region->y1;
         psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
         psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
         psDeproject (&ptSky, &ptTP, fpa->toSky);
-        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);
-
+        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);
+        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);
+
+        // DOWN 1,1
+        psProject (&ptTP, &ptSky, fpa->toSky);
+        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
+        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
+        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);
+        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);
+
+        // UP 0,1
         ptCP.x = region->x0; ptCP.y = region->y1;
         psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
         psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
         psDeproject (&ptSky, &ptTP, fpa->toSky);
-        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);
-
+        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);
+        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);
+
+        // DOWN 0,1
+        psProject (&ptTP, &ptSky, fpa->toSky);
+        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
+        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
+        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);
+        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);
+
+        // UP 0,0
         ptCP.x = region->x0; ptCP.y = region->y0;
         psPlaneTransformApply (&ptFP, chip->toFPA, &ptCP);
         psPlaneTransformApply (&ptTP, fpa->toTPA, &ptFP);
         psDeproject (&ptSky, &ptTP, fpa->toSky);
-        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);
+        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);
+
+        // DOWN 0,0
+        psProject (&ptTP, &ptSky, fpa->toSky);
+        psPlaneTransformApply (&ptFP, fpa->fromTPA, &ptTP);
+        psPlaneTransformApply (&ptCP, chip->fromFPA, &ptFP);
+        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);
 
         psFree (region);
   }
 
-  fclose (f);
+  fclose (fu);
+  fclose (fd);
   psFree (view);
   return true;