IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Mar 27, 2009, 1:26:55 PM (17 years ago)
Author:
Paul Price
Message:

Can now combine images without convolving, but it doesn't apply any rejection yet.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ppStack/src/ppStackCombineInitial.c

    r23573 r23576  
    2323    psTimerStart("PPSTACK_INITIAL");
    2424
    25     psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, PPSTACK_RECIPE); // ppStack recipe
    26     psAssert(recipe, "We've thrown an error on this before.");
    27 
    2825    psMetadata *ppsub = psMetadataLookupMetadata(NULL, config->recipes, "PPSUB"); // PPSUB recipe
    2926    int overlap = 2 * psMetadataLookupS32(NULL, ppsub,
    3027                                          "KERNEL.SIZE"); // Overlap by kernel size between consecutive scans
    3128
    32     pmFPAview *view = NULL;             // View to readout
    33 
    34     int row0, col0;                 // Offset for readout
    35     int numCols, numRows;           // Size of readout
    36     ppStackThread *thread = stack->threads->data[0]; // Representative thread
    37     if (!pmReadoutStackSetOutputSize(&col0, &row0, &numCols, &numRows, thread->readouts)) {
    38         psError(PS_ERR_UNKNOWN, false, "problem setting output readout size.");
    39         return false;
    40     }
    41 
    42     pmFPAfileActivate(config->files, false, NULL);
    43     ppStackFileActivation(config, PPSTACK_FILES_COMBINE, true);
    44     view = ppStackFilesIterateDown(config);
    45     if (!view) {
    46         return false;
    47     }
    48 
    49     pmCell *outCell = pmFPAfileThisCell(config->files, view, "PPSTACK.OUTPUT"); // Output cell
    50     options->outRO = pmReadoutAlloc(outCell); // Output readout
    51     psFree(view);
    52 
    53     psString maskBadStr = psMetadataLookupStr(NULL, recipe, "MASK.BAD"); // Name of bits to mask for bad
    54     psImageMaskType maskBad = pmConfigMaskGet(maskBadStr, config); // Bits to mask for bad pixels
    55     if (!pmReadoutStackDefineOutput(options->outRO, col0, row0, numCols, numRows, true, true, maskBad)) {
    56         psError(PS_ERR_UNKNOWN, false, "Unable to prepare output.");
    57         return false;
    58     }
    59 
    60     psFree(options->cells); options->cells = NULL;
    6129
    6230    bool status;                    // Status of read
Note: See TracChangeset for help on using the changeset viewer.