Changeset 5837 for trunk/psphot/src/psphotTest.c

Timestamp:

Dec 23, 2005, 1:58:45 PM (21 years ago)

Author:

eugene

Message:

fixing output functions

File:

• trunk/psphot/src/psphotTest.c (modified) (1 diff)

trunk/psphot/src/psphotTest.c

@@ -3 +3 @@
 int main (int argc, char **argv) {
 
-    bool status;
+    char line[80];
+    psImage *image;
 
-    psMetadata *config = psphotTestArguments (&argc, argv);
-    eamReadout *imdata = psphotSetup (config);
-
-    int x = psMetadataLookupF32 (&status, config, "TEST_FIT_X");
-    int y = psMetadataLookupF32 (&status, config, "TEST_FIT_Y");
-    float fraction = psMetadataLookupF32 (&status, config, "TEST_FRACTION");
-    float threshold = psMetadataLookupF32 (&status, config, "TEST_THRESHOLD");
-
-    if (threshold == -1) {
-        threshold = fraction * imdata->image->data.F32[y][x];
+    image = psImageAlloc (512, 512, PS_DATA_F32);
+    for (int j = 0; j < 512; j++) {
+        for (int i = 0; i < 512; i++) {
+            image->data.F32[j][i] = i;
+        }
     }
 
-    psArray *contour = pmSourceContour_EAM (imdata->image, x, y, threshold);
-    if (contour == NULL) {
-        fprintf (stderr, "error: invalid contour\n");
-        exit (1);
+    psMetadata *row = NULL;
+
+    psArray *table = psArrayAlloc (100);
+    table->n = 0;
+    
+    for (int i = 0; i < 100; i++) {
+        row = psMetadataAlloc ();
+        sprintf (line, "line %03d", i);
+        psMetadataAdd (row, PS_LIST_HEAD, "FLOAT_DATA", PS_DATA_F32,    "", 0.2*i);
+        psMetadataAdd (row, PS_LIST_HEAD, "INT_DATA",   PS_DATA_S32,    "", i);
+        psMetadataAdd (row, PS_LIST_HEAD, "CHAR_DATA",  PS_DATA_STRING, "", line);
+        psArrayAdd (table, 100, row);
+        // psFree (row);
     }
 
-    psVector *xV = contour->data[0];
-    psVector *yV = contour->data[1];
+    psFits *fits = psFitsOpen ("test.fits", "w");
+    psFitsWriteImage (fits, NULL, image, 0);
+    psFitsWriteTable (fits, NULL, table);
+    psFitsClose (fits);
 
-    for (int i = 0; i < xV->n; i++) {
-        fprintf (stdout, "%d %f %f\n", i, xV->data.F32[i], yV->data.F32[i]);
-    }
-
-    // psSparseMatrixTest ();
     exit (0);
 }