Changeset 38533 for branches/eam_branches/ipp-20150616/Ohana/src/addstar/src/LoadData.c

Timestamp:

Jun 24, 2015, 9:21:59 AM (11 years ago)

Author:

eugene

Message:

converting to find_matches

File:

• branches/eam_branches/ipp-20150616/Ohana/src/addstar/src/LoadData.c (modified) (2 diffs)

branches/eam_branches/ipp-20150616/Ohana/src/addstar/src/LoadData.c

@@ -10 +10 @@
 
   Catalog *catalog = NULL;
-  // ALLOCATE (catalog, Catalog, 1);
-  // dvo_catalog_init (catalog, TRUE);
+  ALLOCATE (catalog, Catalog, 1);
+  dvo_catalog_init (catalog, TRUE);
+  ALLOCATE (catalog->average, Average, 1);
+  ALLOCATE (catalog->secfilt, SecFilt, 1);
+  ALLOCATE (catalog->measure, Measure, 1);
+  ALLOCATE (catalog->lensing, Lensing, 1);
 
   if (images[0] == NULL) {
@@ -103 +107 @@
     // replace full input catalog newcat with subset version
     newcat = FilterStars (newcat, &images[0][Nvalid], Nvalid, options);
-    catalog = newcat;
 
-    // XXX I need to use find_matches_closest here, but for testing:
-    // find_matches_closest (region, catalog, newcat, options);
-    // dvo_catalog_free (newcat);
-    // free (newcat);
+    double Rref, Dref;
+    XY_to_RD (&Rref, &Dref, 0.0, 0.0, &images[0][Nvalid].coords);
+    SkyRegion region;
+    region.Rmin = Rref - 180.0;
+    region.Rmax = Rref + 180.0;
+    region.Dmin = Dref -  90.0;
+    region.Dmax = Dref +  90.0;
+    
+    AddstarClientOptions matchOptions = *options;
+    matchOptions.radius = 0.1; // tight radius at this stage
+    matchOptions.calibrate = FALSE;
+    matchOptions.only_match = FALSE;
+    matchOptions.nosort = FALSE;
+
+    find_matches_closest (&region, newcat, catalog, matchOptions);
+    dvo_catalog_free (newcat);
+    free (newcat);
 
     Nvalid++;