Changeset 6812

Timestamp:

Apr 7, 2006, 8:41:21 AM (20 years ago)

Author:

eugene

Message:

updating getstar dvo interactions

Location:

trunk/Ohana/src/getstar

Files:

• include/getstar.h (modified) (2 diffs)
• src/Shutdown.c (added)
• src/getstar.c (modified) (4 diffs)
• src/mkcatalog.c (added)
• src/select_by_image.c (added)
• src/select_by_region.c (added)
• src/wcatalog.c (added)

trunk/Ohana/src/getstar/include/getstar.h

@@ -12 +12 @@
 };
 
+# if (0)
 typedef struct {
   double X;
@@ -24 +25 @@
   int found;
 } Stars;
+# endif
 
 /* global variables set in parameter file */

trunk/Ohana/src/getstar/src/getstar.c

@@ -12 +12 @@
   set_db (&db);
 
+  sky = SkyTableLoadOptimal (CATDIR, SKY_TABLE, GSCFILE, SKY_DEPTH, VERBOSE);
+  SkyTableSetFilenames (sky, CATDIR, "cpt");
+
+  /* choose elements to load based on options */
+  catalog.catflags = LOAD_AVES | LOAD_MEAS | LOAD_MISS | LOAD_SECF;
+
   Nstars = 0;
   stars = NULL;
@@ -20 +26 @@
     case BY_RADIUS:
 
-      /* load corresponding sky regions, associated images */
-      status = dvo_image_lock (&db, ImageCat, 60.0, LCK_SOFT);
-      if (!status) Shutdown ("ERROR: failure to lock image catalog %s", db.filename);
+      /* load corresponding sky regions */
+      skylist = SkyListByPatch (sky, -1, &REGION);
+      for (i = 0; i < skylist[0].Nregions; i++) {
+        catalog.filename = skylist[0].filename[i];
+        load_pt_catalog (&catalog, skylist[0].regions[i]);
+        unlock_catalog (&catalog);
 
-      skylist = SkyListByPatch (sky, -1, &REGION);
-      overlap = gimages (&db, &image, MOSAIC, &Noverlap);
-      dvo_image_unlock (&db);
+        /* skip empty catalogs */
+        if (catalog.Nave_disk == 0) continue;
+        select_by_region (&output, &catalog, &REGION, 0, 0);
+      }
+      break;
 
+    case BY_IMLIST:
+      /* load image list */
+    case BY_IMAGE:
+
+      # if (0)
+      /* load corresponding sky regions */
+      skylist = SkyListByImage (sky, -1, &image);
       for (i = 0; i < skylist[0].Nregions; i++) {
-          
           catalog.filename = skylist[0].filename[i];
           load_pt_catalog (&catalog, skylist[0].regions[i]);
@@ -36 +53 @@
           /* skip empty catalogs */
           if (catalog.Nave_disk == 0) continue;
-
-          stars = find_matches (&catalog, &image, 0, 0, stars, &Nstars);
-
+          stars = select_by_image (&catalog, &image, 0, 0, stars, &Nstars);
+      }
+      # endif
       break;
 
@@ -45 +62 @@
       break;
 
-    case BY_IMAGE:
-
-      break;
-
-    case BY_IMLIST:
-
+      /* not implemented */
       break;
   }
 
-  image = gimages (argv[1]);
-
-  region = gregions (&image[0], &Nregions);
-
-  stars = (Stars *) NULL;
-  Nstars = 0;
-
-  for (i = 0; i < Nregions; i++) {
-    if (VERBOSE) fprintf (stderr, "extracting from %s\n", region[i].filename);
-    catalog.filename = region[i].filename;  /* don't free region before catalog! */
-    switch (lock_catalog (&catalog, LCK_SOFT)) {
-    case 0:
-      fprintf (stderr, "ERROR: can't lock file\n");
-      exit (1);
-    case 1:
-      gcatalog (&catalog);
-      break;
-    case 2:
-      fprintf (stderr, "no data in %s, skipping\n", catalog.filename);
-      catalog.Naverage = 0;
-      catalog.Nmeasure = 0;
-      break;
-    default:
-      fprintf (stderr, "weird lock_catalog exit state\n");
-      exit (1);
-    }
-    unlock_catalog (&catalog);
-    stars = find_matches (&catalog, &image[0], 0, 0, stars, &Nstars);
-  }
-
   /* write out the selected stars */
-  clean_stars (stars, &Nstars, image);
-  wstars (argv[2], stars, Nstars, image);
+  wcatalog (argv[2], output);
   fprintf (stderr, "SUCCESS\n");
   exit (0);