Changeset 26662

Timestamp:

Jan 22, 2010, 2:02:53 AM (16 years ago)

Author:

bills

Message:

add option to save the expanded background model file

Location:

trunk/ppbgrestore/src

Files:

• ppbgrestoreArguments.c (modified) (2 diffs)
• ppbgrestoreLoop.c (modified) (5 diffs)
• ppbgrestoreParseCamera.c (modified) (1 diff)

trunk/ppbgrestore/src/ppbgrestoreArguments.c

@@ -8 +8 @@
 pmConfig *ppbgrestoreArguments(int argc, char *argv[])
 {
+    int argnum;
     if (argc == 1) {
         usage();
@@ -32 +33 @@
         usage();
     }
+
+    if ((argnum = psArgumentGet(argc, argv, "-savebg"))) {
+        psArgumentRemove(argnum, &argc, argv);
+        psMetadataAddBool(config->arguments, PS_LIST_TAIL, "SAVE_BACKGROUND", 0, "save background innmage", true);
+    }
+
     if (argc != 2) usage ();
 

trunk/ppbgrestore/src/ppbgrestoreLoop.c

@@ -30 +30 @@
         exit(PS_EXIT_PROG_ERROR);
     }
+    pmFPAfile *background= psMetadataLookupPtr(&status, config->files, "PPBGRESTORE.BACKGROUND");
+    if (!status) {
+        psErrorStackPrint(stderr, "Can't find output data!\n");
+        ppbgrestoreCleanup(config, options);
+        exit(PS_EXIT_PROG_ERROR);
+    }
+    bool save_background = psMetadataLookupBool(&status, config->arguments, "SAVE_BACKGROUND");
 
     pmFPAview *view = pmFPAviewAlloc(0);// View for level of interest
@@ -90 +97 @@
         }
         output->save = true;
+        background->save = save_background;
 
         psImageBinning *binning = psphotBackgroundBinning(image, config);
@@ -96 +104 @@
             return false;
         }
-        pmReadout *background = pmFPAfileDefineInternal (config->files, "PSPHOT.BACKGND", image->numCols, image->numRows, PS_TYPE_F32);
+        pmReadout *backgroundRO = pmFPAviewThisReadout(&roView, background->fpa);
+//        pmReadout *backgroundRO = pmFPAfileDefineInternal (config->files, "PSPHOT.BACKGND", image->numCols, image->numRows, PS_TYPE_F32);
+        backgroundRO->image = psImageAlloc(image->numCols, image->numRows, PS_TYPE_F32);
 
         // linear interpolation to full-scale
-        if (!psImageUnbin (background->image, modelRO->image, binning)) {
+        if (!psImageUnbin (backgroundRO->image, modelRO->image, binning)) {
             psError (PSPHOT_ERR_PROG, true, "inconsistent sizes for unbinning");
             psFree(binning);
@@ -105 +115 @@
         }
         psFree(binning);
-        psF32 **backData = background->image->data.F32;
+        psF32 **backData = backgroundRO->image->data.F32;
 
-        // Do the background subtraction
+        // Undo the background subtraction
         int numCols = image->numCols, numRows = image->numRows; // Size of image
         long nancount = 0;
@@ -137 +147 @@
             psWarning("%ld infinite pixels found in background %ld pixels were finite\n", nancount, finitecount);
         }
-        pmFPAfileDropInternal(config->files, "PSPHOT.BACKGND");
+//        pmFPAfileDropInternal(config->files, "PSPHOT.BACKGND");
 
         // Close chip

trunk/ppbgrestore/src/ppbgrestoreParseCamera.c

@@ -23 +23 @@
     pmFPAfile *outImage = pmFPAfileDefineOutput(config, input->fpa, "PPBGRESTORE.OUTPUT");
     if (!outImage) {
-        psError(PS_ERR_IO, false, _("Unable to generate output file from PPIMAGE.OUTPUT"));
+        psError(PS_ERR_IO, false, _("Unable to generate output file from PPBGRESTORE.OUTPUT"));
+        psFree(options);
+        return NULL;
+    }
+    pmFPAfile *outBG = pmFPAfileDefineOutput(config, input->fpa, "PPBGRESTORE.BACKGROUND");
+    if (!outBG) {
+        psError(PS_ERR_IO, false, _("Unable to generate bacgkround file from PPBGRESTORE.BACKGROUND"));
         psFree(options);
         return NULL;