Changeset 7067 for trunk/ppMerge/src/ppMergeCombine.c
- Timestamp:
- May 4, 2006, 3:53:10 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/ppMerge/src/ppMergeCombine.c (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppMerge/src/ppMergeCombine.c
r7061 r7067 4 4 #include <psmodules.h> 5 5 6 #include "ppMerge.h" 6 7 #include "ppMergeData.h" 7 8 #include "ppMergeCombine.h" … … 21 22 // Sanity checks 22 23 assert(!options->scale || scales); 23 assert(!scales || !*scales || ((*scales)->type.type == PS_TYPE_F64 &&24 (*scales)->numCols == data->numCells &&25 (*scales)->numRows == filenames->n));24 assert(!scales || (scales->type.type == PS_TYPE_F64 && 25 scales->numCols == data->numCells && 26 scales->numRows == filenames->n)); 26 27 assert(!options->zero || zeros); 27 assert(!zeros || !*zeros || ((*zeros)->type.type == PS_TYPE_F64 &&28 (*zeros)->numCols == data->numCells &&29 (*zeros)->numRows == filenames->n));28 assert(!zeros || (zeros->type.type == PS_TYPE_F64 && 29 zeros->numCols == data->numCells && 30 zeros->numRows == filenames->n)); 30 31 31 32 // Iterate over the FPA … … 38 39 continue; 39 40 } 40 pmChipWrite(chip, outFile, config->database, false);41 pmChipWrite(chip, data->outFile, config->database, false); 41 42 psArray *cells = chip->cells; // Array of output cells 42 43 for (int j = 0; j < cells->n; j++) { … … 46 47 } 47 48 cellNum++; 48 pmCellWrite(cell, outFile, config->database, false);49 pmCellWrite(cell, data->outFile, config->database, false); 49 50 pmReadout *readout = pmReadoutAlloc(cell); // Output readout of interest 50 51 psArray *stack = psArrayAlloc(filenames->n); // Stack of readouts to combine … … 53 54 cellScales = psImageCol(NULL, scales, cellNum); 54 55 } 55 psVector cellZeros = NULL; // Zeros for this cell56 psVector *cellZeros = NULL; // Zeros for this cell 56 57 if (zeros) { 57 58 cellZeros = psImageCol(NULL, zeros, cellNum); … … 81 82 82 83 // Only reading and writing the first readout in each cell (plane 0) 83 stillReading |= pmReadoutReadNext(stack->data[k], fits, 0, options->rows);84 stillReadingRows |= pmReadoutReadNext(stack->data[k], fits, 0, options->rows); 84 85 psFitsClose(fits); 85 86 } 86 87 87 pmReadoutCombine(readout, stack, cellZeros, cellScales, combineParams);88 pmReadoutWriteNext(readout, outFile, 0);88 pmReadoutCombine(readout, stack, cellZeros, cellScales, options->combine); 89 pmReadoutWriteNext(readout, data->outFile, 0); 89 90 90 91 } while (stillReadingRows); 91 92 92 93 // Write the pixels 93 pmCellWrite(cell, outFile, config->database, true);94 pmCellWrite(cell, data->outFile, config->database, true); 94 95 95 96 // Blow away the cell data … … 104 105 105 106 // Write the pixels 106 pmChipWrite(chip, outFile, config->database, true);107 pmChipWrite(chip, data->outFile, config->database, true); 107 108 108 109 // Blow away the chip data … … 117 118 118 119 // Write the pixels 119 pmFPAWrite( fpa,outFile, config->database, true);120 pmFPAWrite(data->out, data->outFile, config->database, true); 120 121 121 122 // Blow away the FPA data … … 124 125 pmFPAFreeData(fpaIn); 125 126 } 126 pmFPAFreeData( fpa);127 pmFPAFreeData(data->out); 127 128 128 129 return true;
Note:
See TracChangeset
for help on using the changeset viewer.
