Changeset 34470 for branches/czw_branch/20120906/pswarp/src/pswarpTransformTile.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/pswarpTransformTile.c (modified) (3 diffs)

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

@@ -86 +86 @@
     for (int y = yMin; y < yMax; y++) {
         for (int x = xMin; x < xMax; x++) {
-/*        if ((x == 931)||(x == 717)||(x == 924)) { */
-/* /\*      if ((y == 1348)||(y == 1249)||(y == 1081)) { *\/ */
-/* /\*        fprintf(stderr,"BEcause: %d %d [%d %d %d %d]\n",x,y,yMin,yMax,xMin,xMax); *\/ */
-/* /\*      } *\/ */
-/*        } */
-            if (((x == 931)&&(y == 1348))||
-                ((x == 717)&&(y == 1249))||
-                ((x == 924)&&(y == 1081))) {
-              psTrace("pswarp",3,"TT: A\n");
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-
-            }
-
             // Only transform those pixels requested
             if (region && region->data.PS_TYPE_IMAGE_MASK_DATA[y][x]) {
                 continue;
             }
-            if (((x == 931)&&(y == 1348))||
-                ((x == 717)&&(y == 1249))||
-                ((x == 924)&&(y == 1081))) {
-              psTrace("pswarp",3,"TT: B\n");
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-
-            }
-
             // pswarpMapApply converts the output coordinate (x,y) to the input coordinate.
             // both are in the parent frames of the input and output images.
             double xIn, yIn;            // Input pixel coordinates
             pswarpMapApply(&xIn, &yIn, map, x + 0.5, y + 0.5);
-            if (((x == 931)&&(y == 1348))||
-                ((x == 717)&&(y == 1249))||
-                ((x == 924)&&(y == 1081))) {
-              psTrace("pswarp",3,"TT: Ct\n");
-              psTrace("pswarp",3,"TT: %d %d %d %d %g %g %d %d\n",xMin,xMax,yMin,yMax,xIn,yIn,inNumCols,inNumRows);
 
-            }
+            if (args->interp->mode == 8) {
+              xIn += 177.0 / 400.0;
+              yIn += 166.0 / 400.0;
 
+              if ((xIn > inNumCols - 177.0/ 400.0)||
+                  (yIn > inNumRows - 166.0/ 400.0)||
+                  (xIn < 177.0/ 400.0)||
+                  (yIn < 166.0/ 400.0)
+                  
+                  ) {
+                fprintf(stderr,"\n");
+                continue;
+              }
+            }           
+            
             if (xIn < 0 || xIn >= inNumCols || yIn < 0 || yIn >= inNumRows) {
                 continue;
             }
-            if (((x == 931)&&(y == 1348))||
-                ((x == 717)&&(y == 1249))||
-                ((x == 924)&&(y == 1081))) {
-              psTrace("pswarp",3,"TT: C\n");
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-
-            }
 
             // psImagePixelInterpolate determines the value at pixel coordinate (x,y) in child coordinates
@@ -144 +124 @@
             int xOut = x - outCol0, yOut = y - outRow0; ///< Position on output image
 
-            if ((xIn >= args->interp->image->numCols - 1)||
-                (yIn >= args->interp->image->numRows - 1)) {
-              imageValue = NAN;
+            if (((xOut % 500) == 0)&&((yOut % 500) == 0)) {
+              fprintf(stderr,"%d %d %g %g %d %d (%d %d %d %d) %g %g %g\n",
+                      x,y,xIn,yIn,xOut,yOut,xMin,xMax,yMin,yMax,imageValue,
+                      (4846 - 400 * (12)) / 2.0,
+                      (4868 - 400 * (12)) / 2.0
+                      
+                      );
             }
+            
+/*          if ((xIn >= args->interp->image->numCols - 1)|| */
+/*              (yIn >= args->interp->image->numRows - 1)) { */
+/*            imageValue = NAN; */
+/*          } */
             if (outImageData) {
                 outImageData[yOut][xOut] = imageValue * jacobian;
@@ -159 +148 @@
 
             goodPixels++;
-
-            if ((x == 931)||(x == 717)||(x == 924)) {
-              if ((y == 1348)||(y == 1249)||(y == 1081)) {
-              psTrace("pswarp",3,"TT: %d %d %g %g %d %d %g %g (%g) %ld %g %g %g %g\n",x,y,xIn,yIn,xOut,yOut,imageValue,jacobian,
-                      args->input->image->data.F32[(int) yIn][(int) xIn],(long) args->interp,
-                      0.0,0.0,0.0,0.0);
-/*                    inImage->data.F32[(int) yIn][(int) xIn], */
-/*                    inImage->data.F32[(int) yIn+1][(int) xIn], */
-/*                    inImage->data.F32[(int) yIn+1][(int) xIn+1], */
-/*                    inImage->data.F32[(int) yIn][(int) xIn+1]); */
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-              }
-            }
-            
-            if (((xOut == 931)&&(yOut == 1348))||
-                ((xOut == 717)&&(yOut == 1249))||
-                ((xOut == 924)&&(yOut == 1081))) {
-              psTrace("pswarp",3,"TT: %d %d %g %g %d %d %g %g (%g) %ld %g %g %g %g\n",x,y,xIn,yIn,xOut,yOut,imageValue,jacobian,
-                      args->input->image->data.F32[(int) yIn][(int) xIn],(long) args->interp,
-                                      0.0,0.0,0.0,0.0);
-/*                    inImage->data.F32[(int) yIn][(int) xIn], */
-/*                    inImage->data.F32[(int) yIn+1][(int) xIn], */
-/*                    inImage->data.F32[(int) yIn+1][(int) xIn+1], */
-/*                    inImage->data.F32[(int) yIn][(int) xIn+1]); */
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-
-            }
-            if ((xOut ==  443)&&(yOut == 659)) {
-              psTrace("pswarp",3,"TT: %d %d %g %g %d %d %g %g (%g) %ld\n",x,y,xIn,yIn,xOut,yOut,imageValue,jacobian,
-                      args->input->image->data.F32[(int) yIn][(int) xIn],(long) args->interp);
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-
-            }
-            if ((xOut ==  524)&&(yOut == 576)) {
-              psTrace("pswarp",3,"TT: %d %d %g %g %d %d %g %g (%g) %ld\n",x,y,xIn,yIn,xOut,yOut,imageValue,jacobian,
-                      args->input->image->data.F32[(int) yIn][(int) xIn],(long) args->interp);
-              psTrace("pswarp",3,"TT: %d %d %d %d\n",xMin,xMax,yMin,yMax);
-              
-            }
         }
     }