Changeset 20999 for trunk/psModules/src/imcombine/pmPSFEnvelope.c

Timestamp:

Dec 15, 2008, 3:58:30 PM (17 years ago)

Author:

Paul Price

Message:

Attempting to get correct scaling for stacked images.

File:

• trunk/psModules/src/imcombine/pmPSFEnvelope.c (modified) (4 diffs)

trunk/psModules/src/imcombine/pmPSFEnvelope.c

@@ -118 +118 @@
         psf->residuals = NULL;
         if (!pmReadoutFakeFromSources(fakeRO, fakeSize, fakeSize, fakes, xOffset, yOffset, psf,
-                                      NAN, radius, true)) {
+                                      NAN, radius, true, true)) {
             psError(PS_ERR_UNKNOWN, false, "Unable to generate fake readout.");
             psFree(envelope);
@@ -148 +148 @@
             // Get the radius
             pmModel *model = pmModelFromPSFforXY(psf, x, y, PEAK_FLUX); // Model for source
-            psAssert (model, "failed to generate model: should this be an error or not?");
+            psAssert (model, "failed to generate model: should this be an error or not?");
             float srcRadius = model->modelRadius(model->params, PS_SQR(WEIGHT_VAL)); // Radius for source
             if (srcRadius > maxRadius) {
@@ -299 +299 @@
         pmModel *fakeModel = pmModelFromPSFforXY(psf, (float)numCols / 2.0, (float)numRows / 2.0,
                                                  1.0); // Fake model, with central intensity of 1.0
-        psAssert (fakeModel, "failed to generate model: should this be an error or not?");
+        psAssert (fakeModel, "failed to generate model: should this be an error or not?");
         float flux0 = fakeModel->modelFlux(fakeModel->params); // Flux for central intensity of 1.0
         for (int i = 0; i < numFakes; i++) {
@@ -308 +308 @@
         pmReadout *generated = pmReadoutAlloc(NULL); // Generated image
         pmReadoutFakeFromSources(generated, numCols, numRows, fakes, NULL, NULL, psf, NAN, radius,
-                                 false);
+                                 false, true);
         {
             psFits *fits = psFitsOpen("psf_field_model.fits", "w");