Changeset 12492 for trunk/psastro/src/psastroDemoDump.c

Timestamp:

Mar 18, 2007, 12:28:16 PM (19 years ago)

Author:

eugene

Message:

a variety of fixes to get psastro working in chip and mosaic mode

• various memory leaks
• when loading mosaic data, the modification of the WCS for the common focal-plane scale was getting the offset wrong (fixed in psModules)
• added some debugging dump; more plots still needed
• error handling in psastroMosaicAstrom
• loading chip and fpa dimensions from concepts
• now reading chip concepts when reading objects (in psModules)
• now supplying region and subdivisions for gradient measurements
• grid search is now optional

-

File:

• trunk/psastro/src/psastroDemoDump.c (modified) (1 diff)

trunk/psastro/src/psastroDemoDump.c

@@ -1 +1 @@
 # include "psastro.h"
+
+// this function is used for test purposes (-trace psastro.dump.psastroAstromGuess 1)
+bool psastroDumpStars (psArray *stars, char *filename) {
+
+    FILE *f = fopen (filename, "w");
+
+    for (int i = 0; i < stars->n; i++) {
+        pmAstromObj *obj = stars->data[i];
+
+        // write out the upward projections
+        fprintf (f, "%d  %f %f  %f  %f %f  %f %f  %f %f\n", i,
+                 obj->sky->r, obj->sky->d, obj->Mag, 
+                 obj->TP->x, obj->TP->y, 
+                 obj->FP->x, obj->FP->y, 
+                 obj->chip->x, obj->chip->y);
+    }
+    fclose (f);
+    return true;
+}
 
 // this function is used for test purposes (-trace psastro.dump.psastroAstromGuess 1)
 bool psastroDumpRawstars (psArray *rawstars, pmFPA *fpa, pmChip *chip) {
 
-    FILE *f1 = fopen ("rawstars.dat", "w");
-    FILE *f2 = fopen ("rawstars.down.dat", "w");
+    char *filename = NULL;
+    char *chipname = psMetadataLookupStr (NULL, chip->concepts, "CHIP.NAME");
+
+    psStringAppend (&filename, "rawstars.up.%s.dat", chipname);
+    FILE *f1 = fopen (filename, "w");
+    psFree (filename);
+    filename = NULL;
+
+    psStringAppend (&filename, "rawstars.dn.%s.dat", chipname);
+    FILE *f2 = fopen (filename, "w");
+    psFree (filename);
+    filename = NULL;
 
     for (int i = 0; i < rawstars->n; i++) {