Changeset 19311

Timestamp:

Sep 2, 2008, 9:09:51 AM (18 years ago)

Author:

eugene

Message:

measure the gradients but do not apply: this function is now called as part of the distortion model fit

File:

• trunk/psastro/src/psastroMosaicGradients.c (modified) (4 diffs)

trunk/psastro/src/psastroMosaicGradients.c

@@ -2 +2 @@
 static int nPass = 0;
 
-bool psastroMosaicGradients (pmFPA *fpa, psMetadata *recipe) {
+bool psastroMosaicDistortionFromGradients (pmFPA *fpa, psMetadata *recipe) {
 
     bool status;
@@ -9 +9 @@
     pmReadout *readout = NULL;
     psArray *gradients = NULL;
+
+    // Measure the gradient as a function of position.  This must be performed between the
+    // corrected ref->TP and the observed raw->FP, for which the distortion is a perturbation.
 
     pmFPAview *view = pmFPAviewAlloc (0);
@@ -59 +62 @@
     }
 
+    // Fit the gradient field and convert to the distortion terms.
+
     // allocate mosaic-level polynomial transformation and set masks needed by DVO
     int order = psMetadataLookupF32 (&status, recipe, "PSASTRO.MOSAIC.ORDER");
@@ -95 +100 @@
     }
         
-    if (!psastroMosaicSetAstrom (fpa)) {
-        psError(PSASTRO_ERR_UNKNOWN, false, "failed to apply mosaic distortion terms\n");
-        psFree (gradients);
-        psFree (view);
-        return false;
-    }
-
     psFree (gradients);
     psFree (view);
     return true;
 }
+