Changeset 34862

Timestamp:

Dec 19, 2012, 2:09:54 PM (14 years ago)

Author:

bills

Message:

Add recipe value to turn off stacked background construction

Location:

tags/ipp-20121218

Files:

• ippconfig/recipes/ppStack.config (modified) (1 diff, 1 prop)
• ppStack (modified) (1 prop)
• ppStack/src (modified) (1 prop)
• ppStack/src/ppStackCombinePrepare.c (modified) (2 diffs)
• ppStack/src/ppStackLoop.c (modified) (3 diffs)

tags/ipp-20121218/ippconfig/recipes/ppStack.config

@@ -32 +32 @@
 TARGET.FRAC     F32     0.1             # Minimum flux fraction for target PSF
 TARGET.MINMAG   F32     3.0             # Minimum magnitude difference to tolerate in stamp
+
+BACKGROUND.MODEL    BOOL   FALSE        # Construct a stacked version of the warp stage background
 
 RENORM          BOOL    TRUE            # Renormalise weight maps?

tags/ipp-20121218/ppStack/src

@@ -6 +6 @@
 *.la
 ppStack
+ppStackMedian
 config.h
 config.h.in
@@ -12 +13 @@
 ppStackErrorCodes.c
 ppStackErrorCodes.h
+

@@ -6 +6 @@
 *.la
 ppStack
+ppStackMedian
 config.h
 config.h.in
@@ -12 +13 @@
 ppStackErrorCodes.c
 ppStackErrorCodes.h
+

tags/ipp-20121218/ppStack/src/ppStackCombinePrepare.c

@@ -34 +34 @@
 /*       options->expRO = NULL; */
 /*     } */
-    pmCell *bkgCell;
-    int bkg_r0 = 0, bkg_c0 = 0;
-    int bkg_nC = 0, bkg_nR = 0;
+    int bkg_r0 = 0,bkg_c0 = 0;
+    int bkg_nC = 0,bkg_nR = 0;
     if (bkgName) {
       ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, true, 0);
       pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PPSTACK.INPUT.BKGMODEL", 0);
+      if (!file)  {
+        psError(PPSTACK_ERR_ARGUMENTS, false, "Output background model selected but no inputs provided.");
+        return false;
+      }
       pmFPAview *view = ppStackFilesIterateDown(config);
       pmReadout *ro = pmFPAviewThisReadout(view,file->fpa);
@@ -47 +50 @@
       bkg_nC = ro->image->numCols;
       bkg_nR = ro->image->numRows;
-      bkgCell = pmFPAfileThisCell(config->files, view, bkgName); // Bkg cell
+      pmCell *bkgCell = pmFPAfileThisCell(config->files, view, bkgName); // Bkg cell
       
       options->bkgRO = pmReadoutAlloc(bkgCell); // BKG readout

tags/ipp-20121218/ppStack/src/ppStackLoop.c

@@ -41 +41 @@
     if (options->quality) return true; // Can't do anything else
 
+    psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, PPSTACK_RECIPE); // ppStack recipe
+
     // Ensure sufficient inputs
     {
         int numGood = stackSummary(options, "initial combination");
-        psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, PPSTACK_RECIPE); // ppStack recipe
         bool safe = psMetadataLookupBool(NULL, recipe, "SAFE"); // Be safe when combining
         if (safe && numGood <= 1) {
@@ -62 +63 @@
     }
 
+    bool doBackground = psMetadataLookupBool(NULL, recipe, "BACKGROUND.MODEL");
+
     // Prepare for combination
-    if (!ppStackCombinePrepare("PPSTACK.OUTPUT", "PPSTACK.OUTPUT.EXP", "PPSTACK.OUTPUT.BKGMODEL", PPSTACK_FILES_STACK, stack, options, config)) {
+    if (!ppStackCombinePrepare("PPSTACK.OUTPUT", "PPSTACK.OUTPUT.EXP", doBackground ? "PPSTACK.OUTPUT.BKGMODEL" : NULL, PPSTACK_FILES_STACK, stack, options, config)) {
         psError(psErrorCodeLast(), false, "Unable to prepare for combination.");
         psFree(stack);
@@ -143 +146 @@
 
     // Generate median background stack here.
-    if (!ppStackCombineBackground(stack, options, config)) {
+    if (doBackground && !ppStackCombineBackground(stack, options, config)) {
       psError(psErrorCodeLast(), false, "Unable to generate median of background images.");
       psFree(stack);

@@ -6 +6 @@
 *.la
 ppStack
+ppStackMedian
 config.h
 config.h.in
@@ -12 +13 @@
 ppStackErrorCodes.c
 ppStackErrorCodes.h
+