IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 24862


Ignore:
Timestamp:
Jul 20, 2009, 2:39:58 PM (17 years ago)
Author:
Paul Price
Message:

Mask low pixels in the convolved images to avoid holes in subtrahend from creating positive sources in difference image.

Location:
trunk/ppSub/src
Files:
1 added
1 deleted
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/ppSub/src/Makefile.am

    r24672 r24862  
    4545        ppSubReadoutStats.c             \
    4646        ppSubReadoutSubtract.c          \
    47         ppSubSetMasks.c
     47        ppSubSetMasks.c                 \
     48        ppSubThreshold.c
    4849
    4950ppSubKernel_CPPFLAGS = $(PSLIB_CFLAGS) $(PSMODULE_CFLAGS) $(PPSUB_CFLAGS)

  • trunk/ppSub/src/ppSub.h

    r24672 r24862  
    7676/// looking up the parameters in the recipe and supplying them to the function pmSubtractionMatch()
    7777bool ppSubMatchPSFs(ppSubData *data     ///< Processing data
     78    );
     79
     80/// Threshold low pixels in image
     81bool ppSubLowThreshold(ppSubData *data  ///< Processing data
    7882    );
    7983

  • trunk/ppSub/src/ppSubLoop.c

    r24672 r24862  
    6262    if (!ppSubFilesIterateUp(config, PPSUB_FILES_INPUT)) {
    6363        psError(PPSUB_ERR_IO, false, "Unable to close input files.");
     64        return false;
     65    }
     66
     67    if (!ppSubLowThreshold(data)) {
     68        psError(PS_ERR_UNKNOWN, false, "Unable to threshold images.");
    6469        return false;
    6570    }

  • trunk/ppSub/src/ppSubMatchPSFs.c

    r24833 r24862  
    1818#include <pslib.h>
    1919#include <psmodules.h>
    20 #include <psphot.h>
    2120
    2221#include "ppSub.h"

  • trunk/ppSub/src/ppSubReadoutSubtract.c

    r24155 r24862  
    4444        subtrahend = pmFPAfileThisReadout(config->files, view, "PPSUB.REF.CONV");
    4545    }
    46 
    47 // XXX this should be replaced by SAVE options to the inConv and refConf pmFPAfiles
    48 #ifdef TESTING
    49     {
    50         pmReadoutMaskApply(minuend, maskVal);
    51         psFits *fits = psFitsOpen("minuend.fits", "w");
    52         psFitsWriteImage(fits, NULL, minuend->image, 0, NULL);
    53         psFitsClose(fits);
    54     }
    55     {
    56         pmReadoutMaskApply(subtrahend, maskVal);
    57         psFits *fits = psFitsOpen("subtrahend.fits", "w");
    58         psFitsWriteImage(fits, NULL, subtrahend->image, 0, NULL);
    59             psFitsClose(fits);
    60     }
    61 #endif
    6246
    6347    // Do the actual subtraction
Note: See TracChangeset for help on using the changeset viewer.