Changeset 42338 for trunk/psModules/src/detrend/pmOverscan.c

Timestamp:

Jan 30, 2023, 9:53:40 AM (3 years ago)

Author:

eugene

Message:

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

Location:

trunk/psModules

Files:

• . (modified) (1 prop)
• src/detrend/pmOverscan.c (modified) (3 diffs)

trunk/psModules/src/detrend/pmOverscan.c

@@ -136 +136 @@
         break;
       case PM_FIT_SPLINE:
-        // XXX I don't think psSpline1D is up to scratch yet --- it has no mask, and requires an
-        // input spline
-        overscanOpts->spline = psVectorFitSpline1D(reduced, ordinate);
-        fitted = psSpline1DEvalVector(overscanOpts->spline, ordinate);
+
+        // XXX I don't think psSpline1D is up to scratch yet --- it has no mask, and it assumes
+        // a knot for every input point.  it needs an argument like 'number of knots' for the
+        // output spline.  EAM: still true 2023.01.22
+
+        // overscanOpts->spline = psVectorFitSpline1D(reduced, ordinate);
+        // fitted = psSpline1DEvalVector(overscanOpts->spline, ordinate);
+        psError(PS_ERR_UNKNOWN, true, "Spline overscan fitting is broken\n");
         break;
       default:
@@ -185 +189 @@
       }
       case PM_FIT_SPLINE: {
+        /*
           psSpline1D *spline = overscanOpts->spline; // The spline
           for (int i = 0; i < spline->n; i++) {
@@ -197 +202 @@
               comment = NULL;
           }
+        */
           // write metadata header value
           psMetadataAddF32(hdu->header, PS_LIST_TAIL, "OVER_VAL", PS_META_REPLACE,