Changeset 36816 for branches/eam_branches/ipp-20140423/psModules/src/objects

Timestamp:

Jun 5, 2014, 4:02:53 PM (12 years ago)

Author:

eugene

Message:

add some visualization and ability to dump data to pmVisual; clip extreme outliers from psf model

Location:

branches/eam_branches/ipp-20140423/psModules/src/objects

Files:

• pmPSFtryMakePSF.c (modified) (1 diff)
• pmSourceVisual.c (modified) (1 diff)

branches/eam_branches/ipp-20140423/psModules/src/objects/pmPSFtryMakePSF.c

@@ -220 +220 @@
     }
 
+    // weed out extreme e0 outliers here: find the median and exclude points not in the
+    // range MEDIAN / 5 < e0 < 5 * MEDIAN
+    { 
+      psStats *e0stats = psStatsAlloc (PS_STAT_SAMPLE_MEDIAN);
+      if (psVectorStats (e0stats, e0, NULL, srcMask, 0xff)) {
+        float e0med = e0stats->sampleMedian;
+    
+        for (int i = 0; i < sources->n; i++) {
+          // skip any masked sources (failed to fit one of the model steps or get a magnitude)
+          if (srcMask->data.PS_TYPE_VECTOR_MASK_DATA[i]) continue;
+
+          if (e0->data.F32[i] < 0.2*e0med) {
+            srcMask->data.PS_TYPE_VECTOR_MASK_DATA[i] = PSFTRY_MASK_OUTLIER;
+          }
+          if (e0->data.F32[i] > 5.0*e0med) {
+            srcMask->data.PS_TYPE_VECTOR_MASK_DATA[i] = PSFTRY_MASK_OUTLIER;
+          }
+        }
+      }
+      psFree (e0stats);
+    }
+
     // we run 'clipIter' cycles clipping in each of x and y, with only one iteration each.
     // This way, the parameters masked by one of the fits will be applied to the others

branches/eam_branches/ipp-20140423/psModules/src/objects/pmSourceVisual.c

@@ -545 +545 @@
     psFree (model);
 
+    bool dumpData = false;
+
     // pause and wait for user input:
     // continue, save (provide name), ??
-    pmVisualAskUser(&plotPSF);
+retry:
+    pmVisualAskUserOrDump(&plotPSF, &dumpData);
+    if (dumpData) {
+      char name[128];
+      fprintf (stderr, "filename: ");
+      int status = fscanf (stdin, "%127s", name);
+      if (status != 1) {
+        fprintf (stderr, "odd response\n");
+        goto retry;
+      }
+
+      FILE *f = fopen (name, "w");
+      if (!f) {
+        fprintf (stderr, "cannot open %s for output\n", name);
+        goto retry;
+      }
+      for (int i = 0; i < x->n; i++) {
+        float vModel = pmTrend2DEval (trend, x->data.F32[i], y->data.F32[i]);
+        fprintf (f, "%f %f %f %f %d\n", x->data.F32[i], y->data.F32[i], param->data.F32[i], vModel, mask->data.PS_TYPE_VECTOR_MASK_DATA[i]);
+      }
+      fclose (f);
+      goto retry;
+    }
 
     return true;