Changeset 1901 for trunk/psModules/src/pmSubtractBias.c

Timestamp:

Sep 25, 2004, 11:03:05 AM (22 years ago)

Author:

gusciora

Message:

Comments.

File:

• trunk/psModules/src/pmSubtractBias.c (modified) (14 diffs)

trunk/psModules/src/pmSubtractBias.c

@@ -8 +8 @@
  *  @author George Gusciora, MHPCC
  *
- *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-09-24 23:47:21 $
+ *  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-09-25 21:03:05 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -18 +18 @@
 #endif
 
-#include<stdio.h>
-#include<math.h>
-#include "pslib.h"
 #include "pmSubtractBias.h"
 
@@ -57 +54 @@
             (in->image)->data.F32[i][j]-=
                 (bias->image)->data.F32[i+in->row0-bias->row0][j+in->col0-bias->col0];
+
             ((in->mask)->data.U8[i][j])|=
                 (bias->mask)->data.U8[i+in->row0-bias->row0][j+in->col0-bias->col0];
@@ -146 +144 @@
 /******************************************************************************
 ScaleOverscanVector(): this routine takes as input an arbitrary vector,
-creates anew vector of length n, and fills the new vector with the
+creates a new vector of length n, and fills the new vector with the
 interpolated values of the old vector.  The type of interpolation is:
     PM_FIT_POLYNOMIAL: fit a polynomial to the entire input vector data.
@@ -159 +157 @@
                               pmFit fit)
 {
+    if (n == overscanVector) {
+        return(overscanVector);
+    }
+
     psVector *newVec = psVectorAlloc(n, PS_TYPE_F32);
     psPolynomial1D *myPoly;
@@ -166 +168 @@
 
     if (fit == PM_FIT_POLYNOMIAL) {
+        // Fit a polynomial to the old overscan vector.
         myPoly = (psPolynomial1D *) fitSpec;
         myPoly = psVectorFitPolynomial1D(myPoly, NULL, overscanVector, NULL);
+
+        // For each element of the new vector, convert the x-ordinate to that
+        // of the old vector, use the fitted polynomial to determine the
+        // interpolated value at that point, and set the new vector.
         for (i=0;i<n;i++) {
-            x = (((float) i) / ((float) n)) * ((float) overscanVector->n);
+            x = ((float) i) * ((float) overscanVector->n) / ((float) n);
             newVec->data.F32[i] = psPolynomial1DEval(x, myPoly);
         }
     } else if (fit == PM_FIT_SPLINE) {
+        // Fit a spline to the old overscan vector.
         mySpline = (psSpline1D *) fitSpec;
         mySpline = psVectorFitSpline1D(mySpline, NULL, overscanVector, NULL);
+
+        // For each element of the new vector, convert the x-ordinate to that
+        // of the old vector, use the fitted polynomial to determine the
+        // interpolated value at that point, and set the new vector.
         for (i=0;i<n;i++) {
-            x = (((float) i) / ((float) n)) * ((float) overscanVector->n);
+            x = ((float) i) * ((float) overscanVector->n) / ((float) n);
             newVec->data.F32[i] = psSpline1DEval(mySpline, x);
         }
@@ -254 +266 @@
 
     // XXX: How to use multiple overscans?  Currently, only the first is used.
+
     if ((overScanAxis == PM_OVERSCAN_ROWS) || (overScanAxis == PM_OVERSCAN_COLUMNS)) {
         myOverscanImage = (psImage *) ((overscans->head)->data);
@@ -271 +284 @@
                         tmpRow->data.F32[j] = myOverscanImage->data.F32[i][j];
                     }
-                    stat = psVectorStats(stat, tmpRow, NULL, 0xffffffff);
+                    stat = psVectorStats(stat, tmpRow, NULL, 0);
                     p_psGetStatValue(stat, &statValue);
                     overscanVector->data.F32[i] = statValue;
@@ -295 +308 @@
                         tmpRow->data.F32[j] = myOverscanImage->data.F32[i][j];
                     }
-                    stat = psVectorStats(stat, tmpRow, NULL, 0xffffffff);
+                    stat = psVectorStats(stat, tmpRow, NULL, 0);
                     p_psGetStatValue(stat, &statValue);
                     overscanVector->data.F32[i] = statValue;
                 }
-                //                tmpOverscan = tmpOverscan->next;
+                // tmpOverscan = tmpOverscan->next;
                 tmpOverscan = NULL;
             }
@@ -306 +319 @@
         }
 
+        //
+        // Re-bin the overscan vector (change its length).
+        //
         if ((nBin > 0) && (nBin < overscanVector->n)) {
             numBins = 1+((overscanVector->n)/nBin);
@@ -327 +343 @@
 
             }
+
             // Change the effective size of overscanVector.
             overscanVector->n = numBins;
@@ -338 +355 @@
         }
 
+        //
+        // Fit a polynomial to the overscan vector.
+        //
         if (!((fitSpec == NULL) || (fit == PM_FIT_NONE))) {
             if (fit == PM_FIT_POLYNOMIAL) {
@@ -352 +372 @@
         }
 
+        //
+        // Subtract overscan vector row-wise from the image.
+        //
         if (overScanAxis == PM_OVERSCAN_ROWS) {
             for (i=0;i<(in->image)->numCols;i++) {
@@ -360 +383 @@
         }
 
+        //
+        // Subtract overscan vector column-wise from the image.
+        //
         if (overScanAxis == PM_OVERSCAN_COLUMNS) {
             for (i=0;i<(in->image)->numRows;i++) {