IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Apr 2, 2008, 3:57:59 PM (18 years ago)
Author:
Paul Price
Message:

Removing MODE_TARGET mode --- this is now done by generating a fake image (with pmReadoutFakeFromSources) and matching to that with MODE_1 or MODE_2 as appropriate. Changed pmSubtractionMatch and pmSubtractionConvolve so that two output 'convolved' images are always produced. This makes it easier to determine what to subtract from what.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/psModules/src/imcombine/pmSubtractionStamps.c

    r15789 r17297  
    8686    switch (mode) {
    8787      case PM_SUBTRACTION_MODE_1:
    88       case PM_SUBTRACTION_MODE_TARGET:
    8988        maskVal |= PM_SUBTRACTION_MASK_FOOTPRINT_1;
    9089        break;
     
    585584    // Let pmSubtractionStampsSet take care of the rest of the assertions
    586585
    587     const char *format = (mode == PM_SUBTRACTION_MODE_TARGET ? "%f %f %f" : "%f %f"); // Format of file
    588     psArray *data = psVectorsReadFromFile(filename, format);
     586    psArray *data = psVectorsReadFromFile(filename, "%f %f");
    589587    if (!data) {
    590588        psError(PS_ERR_IO, false, "Unable to read stamps file %s", filename);
     
    592590    }
    593591    psVector *x = data->data[0], *y = data->data[1]; // Stamp positions
    594     psVector *flux = (mode == PM_SUBTRACTION_MODE_TARGET ? data->data[2] : NULL); // Stamp fluxes
    595592
    596593    // Correct for IRAF/FITS (unit-offset) positions to C (zero-offset) positions
     
    598595    psBinaryOp(y, y, "-", psScalarAlloc(1.0, PS_TYPE_F32));
    599596
    600     pmSubtractionStampList *stamps = pmSubtractionStampsSet(x, y, flux, image, subMask, region, footprint,
     597    pmSubtractionStampList *stamps = pmSubtractionStampsSet(x, y, NULL, image, subMask, region, footprint,
    601598                                                            spacing, mode);
    602599    psFree(data);
Note: See TracChangeset for help on using the changeset viewer.