Changeset 5565 for trunk/psastro/src/psastroIO.c

Timestamp:

Nov 21, 2005, 5:33:37 PM (21 years ago)

Author:

eugene

Message:

reorg

File:

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

trunk/psastro/src/psastroIO.c

@@ -124 +124 @@
 }
 
+// elixir-style pseudo FITS table (header + ascii list)
+bool testWriteCMP (psMetadata *header, char *filename, psArray *sources) {
+
+    int i;
+    psMetadataItem *mdi;
+
+    psMetadataAdd (header, PS_LIST_TAIL, "NSTARS", PS_DATA_S32    | PS_META_REPLACE, "foo", sources->n);
+
+    // create file, write-out header
+    unlink (filename);
+    psFits *fits = psFitsAlloc (filename);
+
+    // set NAXIS to 0 : CFITSIO requires isolated header to have NAXIS = 0
+    mdi = psMetadataLookup (header, "NAXIS");
+    mdi->data.S32 = 0;
+    mdi->type = PS_DATA_S32;
+
+    psFitsWriteHeader (header, fits);
+    psFree (fits);
+
+    // re-open, add data to end of file
+    FILE *f = fopen (filename, "a+");
+    if (f == NULL) {
+        psLogMsg ("WriteSourceOBJ", 3, "can't open output file for output %s\n", filename);
+        return false;
+    }
+    fseek (f, 0, SEEK_END);
+
+    for (i = 0; i < sources->n; i++) {
+        
+        pmAstromObj *star = sources->data[i];
+
+        fprintf (f, "%6.1f %6.1f %6.3f %03d %2d %3.1f %6.3f %6.3f %6.2f %6.2f %5.1f\n", 
+                 star->pix.x, star->pix.y, star->Mag, (int)(star->dMag*1000), 1, 2.5, 11.0, 12.0, 2.0, 0.5, 15.0);
+    }
+    fclose (f);
+    return true;
+}
+
+// write out raw objects
+bool testWriteRaw (char *filename, psArray *sources) {
+
+    int i;
+
+    // re-open, add data to end of file
+    FILE *f = fopen (filename, "w");
+    if (f == NULL) {
+        psLogMsg ("WriteSourceOBJ", 3, "can't open output file for output %s\n", filename);
+        return false;
+    }
+
+    for (i = 0; i < sources->n; i++) {
+        
+        pmAstromObj *star = sources->data[i];
+
+        fprintf (f, "%8.2f %8.2f   %8.2f %8.2f   %8.2f %8.2f   %8.2f %8.2f   %8.2f %8.2f   %10.6f %10.6f   %8.2f %8.2f\n", 
+                 star->pix.x, star->pix.y, 
+                 star->cell.x, star->cell.y, 
+                 star->chip.x, star->chip.y, 
+                 star->FP.x, star->FP.y, 
+                 star->TP.x, star->TP.y, 
+                 star->sky.r*DEG_RAD, star->sky.d*DEG_RAD, 
+                 star->Mag, star->dMag);
+    }
+    fclose (f);
+    return true;
+}
+
+// elixir-style pseudo FITS table (header + ascii list)
+bool testWriteRef (char *filename, psArray *sources) {
+
+    int i;
+
+    // re-open, add data to end of file
+    FILE *f = fopen (filename, "w");
+    if (f == NULL) {
+        psLogMsg ("WriteSourceOBJ", 3, "can't open output file for output %s\n", filename);
+        return false;
+    }
+
+    for (i = 0; i < sources->n; i++) {
+        
+        pmAstromObj *star = sources->data[i];
+
+        fprintf (f, "%10.7f %10.7f %6.3f %6.3f\n", star->sky.r, star->sky.d+10*RAD_DEG/3600, 10.0, 0.05);
+    }
+    fclose (f);
+    return true;
+}
+
+void psastroDumpStars (psArray *sources) {
+
+    for (int i = 0; i < sources->n; i++) {
+        pmAstromObj *star = sources->data[i];
+
+        fprintf (stderr, "%8.2f %8.2f   %8.2f %8.2f   %8.2f %8.2f   %8.2f %8.2f   %8.2f %8.2f   %10.6f %10.6f   %8.2f %8.2f\n", 
+                 star->pix.x, star->pix.y, 
+                 star->cell.x, star->cell.y, 
+                 star->chip.x, star->chip.y, 
+                 star->FP.x, star->FP.y, 
+                 star->TP.x, star->TP.y, 
+                 star->sky.r*DEG_RAD, star->sky.d*DEG_RAD, 
+                 star->Mag, star->dMag);
+    }
+}