Changeset 34470 for branches/czw_branch/20120906/pswarp/src/pswarpMapGrid.c

Timestamp:

Sep 25, 2012, 5:33:33 PM (14 years ago)

Author:

watersc1

Message:

Hardcoded version that gets the right solution. I wanted to save this before attempting to clean up the code.

File:

• branches/czw_branch/20120906/pswarp/src/pswarpMapGrid.c (modified) (1 diff)

branches/czw_branch/20120906/pswarp/src/pswarpMapGrid.c

@@ -52 +52 @@
     return grid;
 }
+/**
+ * pswarpMapGridFromImageForBackground builds a set (a grid) of locally-linear maps which convert the source
+ * coordinates (src) to destination coordinates (dest).  we construct a grid with superpixel
+ * spacing of nXpix, nYpix.  The transformation for each grid cell is valid for the superpixel.
+ * The grid over-fills the source image so every source image pixel is guaranteed to have a map.
+ */
+pswarpMapGrid *pswarpMapGridFromImageForBackground (pmReadout *dest, pmReadout *src, int nXpix, int nYpix) {
+
+    int i, ni;
+    int j, nj;
+
+    // start counting from the center of the superpixels
+    double xMin = 0.5*nXpix;
+    double yMin = 0.5*nYpix;
+
+    // the map is defined for coordinates in the image parent frame.
+    int Nx = src->image->numCols + src->image->col0;
+    int Ny = src->image->numRows + src->image->row0;
+
+    // always allocate an extra superpixel to handle spillover
+    int nXpts = (int)(Nx / nXpix) + 1;
+    int nYpts = (int)(Ny / nYpix) + 1;
+
+    // create the grid of maps
+    pswarpMapGrid *grid = pswarpMapGridAlloc (nXpts, nYpts);
+
+    // measure the map for the center of each superpixel
+    for (ni = 0, i = xMin; ni < nXpts; i += nXpix, ni++) {
+        for (nj = 0, j = yMin; nj < nYpts; j += nYpix, nj++) {
+            pswarpMapSetLocalModel (grid->maps[ni][nj], dest, src, i, j);
+/*          grid->maps[ni][nj]->Xo += (13.0 * 400.0 - 4846.0) / 2.0; */
+/*          grid->maps[ni][nj]->Yo += (13.0 * 400.0 - 4868.0) / 2.0; */
+        }
+    }
+
+    grid->nXpix = nXpix;
+    grid->nYpix = nYpix;
+    grid->xMin = xMin;
+    grid->yMin = yMin;
+    return grid;
+}
 
 /**