- Timestamp:
- Sep 20, 2012, 4:25:06 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/czw_branch/20120906/ppSkycell/src/ppSkycellLoop.c
r34454 r34456 133 133 psArray *imageRegions = psArrayAlloc(data->numInputs); // Region for image 134 134 psArray *regionHDUs = psArrayAlloc(data->numInputs); 135 // Determine which projection cells we have to deal with. 135 136 for (int i = 0; i < data->numInputs; i++) { 136 137 pmFPAfile *file = pmFPAfileSelectSingle(data->config->files, "PPSKYCELL.IMAGE", i); // File to examine … … 184 185 185 186 pmFPAfileActivate(data->config->files, false, NULL); 186 187 // Loop over projections 187 188 for (int i = 0; i < numProj; i++) { 188 189 psRegion *projRegion = projRegions->data[i]; // Region for skycell projection … … 211 212 int modify_wcs1 = 1; 212 213 int modify_wcs2 = 1; 214 // Loop over inputs to this projection. 213 215 for (int j = 0; j < data->numInputs; j++) { 214 216 if (target->data.S32[j] != i) { … … 249 251 // Offsets for image on skycell 250 252 int xOffset1 = (-imageRegion->x0 + projRegion->x0) / (float)data->bin1; 251 int yOffset1 = ( imageRegion->y0 - projRegion->y0) / (float)data->bin1;253 int yOffset1 = (-imageRegion->y1 + projRegion->y1) / (float)data->bin1; 252 254 int xOffset2 = xOffset1 / (float)data->bin2, yOffset2 = yOffset1 / (float)data->bin2; 253 255 fprintf(stderr,"OFF: %d %d %d %g %g %g %g %g %g %g %g\n",j,xOffset1,yOffset1, 256 imageRegion->x0,imageRegion->y0, 257 imageRegion->x1,imageRegion->y1, 258 projRegion->x0,projRegion->y0, 259 projRegion->x1,projRegion->y1); 254 260 // XXX Completely neglecting rotations 255 261 // The skycells are divided up neatly with them all having the same orientation
Note:
See TracChangeset
for help on using the changeset viewer.
