Changeset 27605

Timestamp:

Apr 5, 2010, 2:31:06 PM (16 years ago)

Author:

Paul Price

Message:

Output ds9 region files.

Location:

trunk/ppViz/src/ppCoord

Files:

• ppCoord.h (modified) (2 diffs)
• ppCoordArguments.c (modified) (2 diffs)
• ppCoordData.c (modified) (2 diffs)
• ppCoordLoop.c (modified) (1 diff)

trunk/ppViz/src/ppCoord/ppCoord.h

@@ -2 +2 @@
 #define PP_VIZPSF_H
 
+#include <stdio.h>
 #include <pslib.h>
 #include <psmodules.h>
+
 
 #define PPVIZPSF_RECIPE "PPVIZPSF"      // Recipe name
@@ -17 +19 @@
     bool radians;                       // RA,Dec are in radians?
     bool all;                           // Only all coordinates?
+    psString ds9name;                   // Name of ds9 region file
+    FILE *ds9;                          // ds9 output file handle
+    float ds9radius;                    // Radius of ds9 regions
+    psString ds9color;                  // Color of ds9 regions
 } ppCoordData;
 

trunk/ppViz/src/ppCoord/ppCoordArguments.c

@@ -53 +53 @@
     psMetadataAddBool(arguments, PS_LIST_TAIL, "-radians", 0, "RA,Dec in radians?", NULL);
     psMetadataAddBool(arguments, PS_LIST_TAIL, "-all", 0, "Output all coordinates?", NULL);
+    psMetadataAddStr(arguments, PS_LIST_TAIL, "-ds9", 0, "Output ds9 region file", NULL);
+    psMetadataAddF32(arguments, PS_LIST_TAIL, "-ds9-radius", 0, "ds9 region radii", 5.0);
+    psMetadataAddStr(arguments, PS_LIST_TAIL, "-ds9-color", 0, "ds9 region color", "red");
     if (!psArgumentParse(arguments, &argc, argv) || argc != 1) {
         usage(argv[0], arguments, data);
@@ -65 +68 @@
     data->all = psMetadataLookupBool(NULL, arguments, "-all");
 //    psMetadataAddStr(data->config->arguments, PS_LIST_TAIL, "OUTPUT", 0, "Output root name", data->outRoot);
+
+    data->ds9name = psMemIncrRefCounter(psMetadataLookupStr(NULL, arguments, "-ds9"));
+    if (data->ds9name) {
+        data->ds9 = fopen(data->ds9name, "w");
+        if (!data->ds9) {
+            psError(PS_ERR_IO, true, "Unable to open region file %s", data->ds9name);
+            return false;
+        }
+    }
+    data->ds9radius = psMetadataLookupF32(NULL, arguments, "-ds9-radius");
+    data->ds9color = psMetadataLookupStr(NULL, arguments, "-ds9-color");
 
     psTrace("ppCoord", 1, "Done reading command-line arguments\n");

trunk/ppViz/src/ppCoord/ppCoordData.c

@@ -19 +19 @@
     psFree(data->radecName);
     psFree(data->config);
+    psFree(data->ds9name);
+    fclose(data->ds9);
+    psFree(data->ds9color);
     return;
 }
@@ -34 +37 @@
     data->radecName = NULL;
     data->config = NULL;
+    data->radians = false;
+    data->all = false;
+    data->ds9name = NULL;
+    data->ds9 = NULL;
+    data->ds9radius = NAN;
+    data->ds9color = NULL;
 
     return data;

trunk/ppViz/src/ppCoord/ppCoordLoop.c

@@ -314 +314 @@
                     rawFile ? " " : "",
                     rawFile && cellName ? cellName : (rawFile ? "UNKNOWN" : ""));
+            if (!rawFile && data->ds9) {
+                // Region file is only appropriate if we're not mapping all the way back to cell coordinates
+                fprintf(data->ds9, "image;circle(%f,%f,%f) # color=%s",
+                        xPix->data.F32[i], yPix->data.F32[i], data->ds9radius, data->ds9color);
+            }
         }
     }