Changeset 34046 for branches/meh_branches/ppstack_test/ppStack/src/ppStackPrepare.c

Timestamp:

Jun 20, 2012, 11:53:21 AM (14 years ago)

Author:

mhuber

Message:

adding parts for simple model mean target PSF

File:

• branches/meh_branches/ppstack_test/ppStack/src/ppStackPrepare.c (modified) (3 diffs)

branches/meh_branches/ppstack_test/ppStack/src/ppStackPrepare.c

@@ -306 +306 @@
     psLogMsg("ppStack", PS_LOG_INFO, "Input FWHMs : %f +/- %f", fwhmStats->clippedMean, fwhmStats->clippedStdev);
 
+    //MEH -- add clippedMean/clippedStdev to options, may want to be part of clipFWHMnSig?
+    options->clippedMean =  fwhmStats->clippedMean;
+    options->clippedStdev = fwhmStats->clippedStdev;
+     
+
     if (isfinite(clipFWHMnSig)) {
         float fwhmLimit = fwhmStats->clippedMean + clipFWHMnSig * fwhmStats->clippedStdev;
@@ -318 +323 @@
     }
 
-    //MEH -- by option could reset maxFWHM and reject any > before going into target ppStackPSF however will be masked and not included in conv/unconv stack in end rather than rej from deconv limit
-
-    //MEH -- could use fwhmStats->clippedMean+X to internally set PSF.OUTPUT.FWHM but called in ppStackPSF and would need to pass it
-
 
     //MEH - try runnning every thing that can reject images before target psf and conv - 
@@ -334 +335 @@
     // Generate target PSF
     if (options->convolve) {
-        options->psf = ppStackPSF(config, numCols, numRows, psfs, options->inputMask);
+      //        options->psf = ppStackPSF(config, numCols, numRows, psfs, options->inputMask);
+        options->psf = ppStackPSF(config, numCols, numRows, psfs, options);
         psFree(psfs);
         if (!options->psf) {