Changeset 42340

Timestamp:

Jan 30, 2023, 9:55:36 AM (3 years ago)

Author:

eugene

Message:

merge from branches/eam_branches/ppImage.20230123: add new non-linearity correction using new psSpline

Location:

trunk/ppImage

Files:

• . (modified) (1 prop)
• src/Makefile.am (modified) (1 diff)
• src/ppImage.h (modified) (2 diffs)
• src/ppImageArguments.c (modified) (2 diffs)
• src/ppImageDetrendFree.c (modified) (1 diff)
• src/ppImageDetrendNewNonLinear.c (copied) (copied from branches/eam_branches/ppImage.20230123/src/ppImageDetrendNewNonLinear.c )
• src/ppImageDetrendNonLinear.c (modified) (1 diff)
• src/ppImageDetrendReadout.c (modified) (2 diffs)
• src/ppImageDetrendRecord.c (modified) (1 diff)
• src/ppImageOptions.c (modified) (2 diffs)
• src/ppImageParseCamera.c (modified) (1 diff)

trunk/ppImage/src/Makefile.am

@@ -23 +23 @@
         ppImageDetrendRecord.c \
         ppImageDetrendNonLinear.c \
+        ppImageDetrendNewNonLinear.c \
         ppImageDetrendFringe.c \
         ppImageDetrendFree.c \

trunk/ppImage/src/ppImage.h

@@ -35 +35 @@
     bool doAuxMask;                     // apply auxillary mask
     bool doNonLin;                      // Non-linearity correction
+    bool doNewNonLin;                   // Non-linearity correction
     bool doOverscan;                    // Overscan subtraction
     bool doNoiseMap;                    // Bias subtraction
@@ -175 +176 @@
 bool ppImageDetrendBias(pmReadout *inputReadout, pmReadout *bias, pmReadout *dark, ppImageOptions *options);
 
+bool ppImageDetrendNewNonLinear(pmReadout *input, pmFPAview *linearity, pmConfig *config);
+
 bool ppImageDetrendNonLinear(pmReadout *input, pmFPAview *linearity, pmConfig *config);
 bool ppImageDetrendNonLinearLookup(pmReadout *input, psMetadataItem *dataItem);

trunk/ppImage/src/ppImageArguments.c

@@ -25 +25 @@
     fprintf(stderr, "\t-fringe/-fringelist: Fringe image and data.\n");
     fprintf(stderr, "\t-linearity/-linearlist: linearity correction file.\n");
+    fprintf(stderr, "\t-newnonlin/-newnonlinlist: non-linearity correction file (v 2023.01).\n");
     fprintf(stderr, "\n");
     exit (2);
@@ -116 +117 @@
     pmConfigFileSetsMD (config->arguments, &argc, argv, "FRINGE", "-fringe", "-fringelist");
     pmConfigFileSetsMD (config->arguments, &argc, argv, "LINEARITY", "-linearity", "-linearlist");
+    pmConfigFileSetsMD (config->arguments, &argc, argv, "NEWNONLIN", "-newnonlin", "-newnonlinlist");
     pmConfigFileSetsMD (config->arguments, &argc, argv, "PATTERN.ROW.AMP", "-pattern-row-amplitude", "-not-defined");
 

trunk/ppImage/src/ppImageDetrendFree.c

@@ -14 +14 @@
     "PPIMAGE.SHUTTER",
     "PPIMAGE.LINEARITY",
+    "PPIMAGE.NEWNONLIN",
     NULL
 };

trunk/ppImage/src/ppImageDetrendNonLinear.c

@@ -1 @@
-#ifdef HAVE_CONFIG_H
-#include <config.h>
-#endif
-
 #include "ppImage.h"
 

trunk/ppImage/src/ppImageDetrendReadout.c

@@ -124 +124 @@
       // psLogMsg ("ppImage", 6, "nonlinear correction: %f sec\n", psTimerMark ("detrend.readout"));
     }
+    // New Non-linearity correction (exclusive of the above)
+    if (options->doNewNonLin) {
+      if (!ppImageDetrendNewNonLinear(input, detview, config)) {
+        psError(PS_ERR_UNKNOWN, false, "Unable to correct Non-Linearity with new version (2023)");
+        psFree(detview);
+        return(false);
+      }
+      // psLogMsg ("ppImage", 6, "nonlinear correction: %f sec\n", psTimerMark ("detrend.readout"));
+    }
 
     // set up the dark and bias
@@ -188 +197 @@
             psImageUnbin (noiseImage, noiseMap->image, binning);
             psFree (binning);
-            // Stolen from pmSkySubtract.c
-            // CZW: Unneeded, as psImageUnbin does the bilinear interpolation?  It still looks blocky. 
-/*          psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, RECIPE_NAME); // Recipe */
-/*          psAssert(recipe, "Should be there!"); */
-
-/*          psS32 xBinFactor = psMetadataLookupS32(NULL,recipe,"NOISE.XBIN"); */
-/*          psS32 yBinFactor = psMetadataLookupS32(NULL,recipe,"NOISE.YBIN"); */
-/*          psImageInterpolateOptions *interp = psImageInterpolateOptionsAlloc(PS_INTERPOLATE_BILINEAR, */
-/*                                                                             noiseMap->image, */
-/*                                                                             NULL, NULL, */
-/*                                                                             0, 0.0, 0.0, 0, 0, 0.0); */
-/*          for (psS32 row = 0; row < input->image->numRows; row++) { */
-/*            for (psS32 col = 0; col < input->image->numCols; col++) { */
-/*              // We calculate the F32 value of the pixel coordinates in the */
-/*              // binned image and then use a pixel interpolation routine to */
-/*              // determine the value of the pixel at that location. */
-/*              psF32 binRowF64 = ((psF32) row) / ((psF32) yBinFactor); */
-/*              psF32 binColF64 = ((psF32) col) / ((psF32) xBinFactor); */
-/*              // We add 0.5 to the pixel locations since the pixel */
-/*              // interpolation routine defines the location of pixel */
-/*              // (i, j) as (i+0.5, j+0.5). */
-/*              binRowF64+= 0.5; */
-/*              binColF64+= 0.5; */
-
-/*              double binPixel; */
-/*              if (!psImagePixelInterpolate(&binPixel, NULL, NULL, binColF64, binRowF64, interp)) { */
-/*                psError(PS_ERR_UNKNOWN, false, "Unable to interpolate image."); */
-/*                psFree(interp); */
-/*                psFree(noiseImage); */
-/*                return NULL; */
-/*              } */
-/*              noiseImage->data.F32[row][col] = binPixel; */
-/*            } */
-/*          } */
-/*          psFree(interp); */
         }
 

trunk/ppImage/src/ppImageDetrendRecord.c

@@ -75 +75 @@
     detrendRecord(options->doFringe,   detrend, config, view, "PPIMAGE.FRINGE",   "DETREND.FRINGE",   "Fringe filename");
 
-    detrendRecord(options->doNonLin,   detrend, config, view, "PPIMAGE.LINEARITY","DETREND.NONLIN",   "Non-linearity table filename");
+    detrendRecord(options->doNonLin,    detrend, config, view, "PPIMAGE.LINEARITY","DETREND.NONLIN",   "Non-linearity table filename");
+    detrendRecord(options->doNewNonLin, detrend, config, view, "PPIMAGE.NEWNONLIN","DETREND.NEWNONLIN","Non-linearity table filename (v2023)");
 
     detrendRecord(options->doDark & options->useVideoDark, detrend, config, view, "PPIMAGE.VIDEODARK", "DETREND.VIDEODARK", "VideoDark filename");

trunk/ppImage/src/ppImageOptions.c

@@ -28 +28 @@
     options->doAuxMask       = false;   // apply auxillary mask
     options->doNonLin        = false;   // Non-linearity correction
+    options->doNewNonLin     = false;   // New Non-linearity correction (v2023)
     options->doOverscan      = false;   // Overscan subtraction
     options->doNoiseMap      = false;   // Apply Read Noise Map
@@ -175 +176 @@
     // XXX PAP: The overscan stuff needs to be updated following the reworked API
 
+    // New Non-linearity (v 2023) recipe options
+    // non-linearity corrections are loaded from a file defined in the detrend system or on the command-line
+    if (psMetadataLookupBool(NULL, recipe, "NEWNONLIN")) {
+        options->doNewNonLin = true;
+    }
+
+    // XXX PAP: The overscan stuff needs to be updated following the reworked API
+
     // Overscan recipe options
     // XXX EAM : we should abort on invalid options. default options?

trunk/ppImage/src/ppImageParseCamera.c

@@ -104 +104 @@
 
     if (options->doNonLin) {
-      if (!ppImageDefineFile(config, input->fpa, "PPIMAGE.LINEARITY", "LINEARITY",
-                             PM_FPA_FILE_LINEARITY, PM_DETREND_TYPE_LINEARITY)) {
-        psError(PS_ERR_IO, false, "Can't find a non-linearity correction source");
-        psFree(options);
-        return NULL;
-      }
+        if (!ppImageDefineFile(config, input->fpa, "PPIMAGE.LINEARITY", "LINEARITY",
+                               PM_FPA_FILE_LINEARITY, PM_DETREND_TYPE_LINEARITY)) {
+            psError(PS_ERR_IO, false, "Can't find a non-linearity correction source");
+            psFree(options);
+            return NULL;
+        }
+    }
+    if (options->doNewNonLin) {
+        // if the file has been specified on the command-line (-newnonlin file), then the file
+        // is identified by the NEWNONLIN entry in config->arguments.  otherwise, load from the
+        // detrend system as type NEWNONLIN
+        if (!ppImageDefineFile(config, input->fpa, "PPIMAGE.NEWNONLIN", "NEWNONLIN",
+                               PM_FPA_FILE_NEWNONLIN, PM_DETREND_TYPE_NEWNONLIN)) {
+            psError(PS_ERR_IO, false, "Can't find a new non-linearity correction source");
+            psFree(options);
+            return NULL;
+        }
     }