Changeset 4585 for trunk/Ohana/src/opihi/dvo/region_list.c

Timestamp:

Jul 20, 2005, 4:27:25 AM (21 years ago)

Author:

eugene

Message:

substantial cleanup of photometry extraction functions

File:

• trunk/Ohana/src/opihi/dvo/region_list.c (modified) (3 diffs)

trunk/Ohana/src/opihi/dvo/region_list.c

@@ -1 +1 @@
 # include "dvo1.h"
+
+static int RegionSelect;
+
+int GetRegionSelection () {
+  return (RegionSelect);
+}
+
+int SetRegionSelection (int *argc, char **argv, char **RegionName, char **RegionList) {
+  
+  int N;
+
+  RegionSelect = FALSE;
+
+  /* check for Region selection */
+  *RegionName = NULL;
+  if (N = get_argument (*argc, argv, "-cpt")) {
+    remove_argument (N, argc, argv);
+    *RegionName = strcreate (argv[N]);
+    remove_argument (N, argc, argv);
+    return (TRUE);
+  }    
+  /* check for Region list */
+  *RegionList = NULL;
+  if (N = get_argument (*argc, argv, "-cptlist")) {
+    remove_argument (N, argc, argv);
+    *RegionList = strcreate (argv[N]);
+    remove_argument (N, argc, argv);
+    return (TRUE);
+  } 
+  if ((*RegionName == NULL) && (*RegionList == NULL)) {
+    RegionSelect = TRUE;
+    return (TRUE);
+  }
+
+  fprintf (stderr, "-cpt and -cptlist are incompatible\n");
+  free (*RegionName);
+  free (*RegionList);
+  *RegionName = NULL;
+  *RegionList = NULL;
+  return (FALSE);
+}
+
+RegionFile *SelectRegions (char *RegionName, char *RegionList, int *nregions) {
+
+  double Radius;
+  int Ngraph;
+  Graphdata graphsky;
+  RegionFile *regions;
+
+  /* determine region-file names */
+  if (RegionName != NULL) {
+    *nregions = 1;
+    ALLOCATE (regions, RegionFile, 1);
+    strcpy (regions[0].name, RegionName);
+    free (RegionName);
+    RegionName = NULL;
+    return (regions);
+  } 
+  if (RegionList != NULL) {
+    regions = region_list (RegionList, nregions);
+    free (RegionList);
+    RegionList = NULL;
+    return (regions);
+  }
+
+  Ngraph = 0;
+  if (!GetGraphData (&graphsky, NULL, &Ngraph)) {
+    fprintf (stderr, "region display not available\n");
+    return (NULL);
+  }
+
+  Radius = MAX (fabs(graphsky.xmax), fabs(graphsky.ymax));
+  regions = find_regions (graphsky.coords.crval1, graphsky.coords.crval2, Radius, nregions);
+  return (regions);
+}
 
 /* returns a list of region files names from file */
@@ -12 +87 @@
     fprintf (stderr, "ERROR: can't find region list file %s\n", filename);
     *Nregions = 0;
-    return ((RegionFile *) NULL);
+    return (NULL);
   }
   
@@ -28 +103 @@
 
   *Nregions = nregion;
-
   return (regions);
-
 }