Changeset 34825

Timestamp:

Dec 15, 2012, 6:49:04 AM (14 years ago)

Author:

eugene

Message:

adding notes on projection cell / skycell mapping

Location:

branches/eam_branches/ipp-20121130/Ohana/src

Files:

• addstar/src/findskycell.c (modified) (4 diffs)
• relphot/src/BoundaryTreeOps.c (modified) (1 diff)

branches/eam_branches/ipp-20121130/Ohana/src/addstar/src/findskycell.c

@@ -15 +15 @@
 
 void usage (void) {
-  fprintf (stderr, "USAGE: findcell -mktree (tree) (catdir)\n");
+  fprintf (stderr, "USAGE: findcell -mktree (tree) (catdir) [-nx Nx] [-ny Ny]\n");
   fprintf (stderr, "USAGE: findcell -tree (tree) (datafile)\n");
   fprintf (stderr, "   (datafile) should contain a list of RA,DEC pairs\n");
@@ -26 +26 @@
 int main (int argc, char **argv) {
 
-  int N;
+  int N, NX_SUB, NY_SUB;
   char *treefile = NULL;
 
@@ -37 +37 @@
   if (get_argument (argc, argv, "-help")) usage ();
   if (get_argument (argc, argv, "-h")) usage ();
+
+  NX_SUB = NY_SUB = 1;
+  if ((N = get_argument (argc, argv, "-nx"))) {
+    remove_argument (N, &argc, argv);
+    NX_SUB = atof (argv[N]);
+    remove_argument (N, &argc, argv);
+  }
+  if ((N = get_argument (argc, argv, "-ny"))) {
+    remove_argument (N, &argc, argv);
+    NY_SUB = atof (argv[N]);
+    remove_argument (N, &argc, argv);
+  }
 
   /* extra error messages */
@@ -326 +338 @@
 
     fprintf (stdout, "%10.6f %10.6f  %3d %3d  %s\n", ra, dec, zone, band, tree->name[zone][band]);
+
+    // I know the projection cell (band,zone), but I need to find the skycell within that projection cell.
+    // the proj cell is divided into Nx, Ny bits.
+    // (ra,dec) for (Ro,Do) -> (x,y).  given (Xo,Yo),(dX,dY) I can find ix,iy
+    // I currently track Ro,Do (tree->RA_origin[zone], tree->...)
+
+
   }
 

branches/eam_branches/ipp-20121130/Ohana/src/relphot/src/BoundaryTreeOps.c

@@ -22 +22 @@
   
   snprintf (primaryCellName, DVO_MAX_PATH, "RINGS.V3.%s", tree->name[zone][band]);
+
+  // I have ra, dec, and the primary projection cell.  In order to choose the primary skycell,
+  // I just need to project to ra,dec to X,Y based on the center of the cell and then get the subdivision right.
+
   return TRUE;
 }