Changeset 17368 for trunk/ppSub/src/ppSubReadout.c

Timestamp:

Apr 7, 2008, 12:35:31 PM (18 years ago)

Author:

Paul Price

Message:

Add renormalisation of weight map.

File:

• trunk/ppSub/src/ppSubReadout.c (modified) (2 diffs)

trunk/ppSub/src/ppSubReadout.c

@@ -65 +65 @@
     int optOrder = psMetadataLookupS32(&mdok, config->arguments, "OPTIMUM.ORDER"); // Order for search
     bool dual = psMetadataLookupBool(&mdok, config->arguments, "DUAL"); // Dual convolution?
+    bool renorm = psMetadataLookupBool(&mdok, config->arguments, "RENORM"); // Renormalise weights?
+    int renormWidth = psMetadataLookupS32(&mdok, config->arguments, "RENORM.WIDTH"); // Width for renormalise
 
     pmSubtractionMode mode = dual ? PM_SUBTRACTION_MODE_DUAL : PM_SUBTRACTION_MODE_UNSURE; // Subtraction mode
@@ -188 +190 @@
     }
 
+    if (renorm) {
+        psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS, 0); // Random number generator
+        if (!pmReadoutWeightRenorm(outRO, maskBlank, PS_STAT_ROBUST_MEDIAN, PS_STAT_ROBUST_STDEV,
+                                   renormWidth, rng)) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to renormalise weights for photometry.");
+            psFree(outRO);
+            return false;
+        }
+    }
+
     if (psMetadataLookupBool(NULL, config->arguments, "PHOTOMETRY")) {
         pmFPAfile *photFile = psMetadataLookupPtr(NULL, config->files, "PSPHOT.INPUT");