Changeset 4375 for trunk/psphot/src/psphot-utils.c

Timestamp:

Jun 24, 2005, 5:38:37 AM (21 years ago)

Author:

eugene

Message:

adding model tests

File:

• trunk/psphot/src/psphot-utils.c (modified) (3 diffs)

trunk/psphot/src/psphot-utils.c

@@ -77 +77 @@
         params = model->params;
         dparams = model->dparams;
-        fprintf (f, "%7.1f %7.1f  %5.1f %7.1f  %7.4f  ", 
+        fprintf (f, "%7.1f %7.1f  %5.1f %7.3f  %7.4f  ", 
                  params[0].data.F32[2], params[0].data.F32[3], 
-                 params[0].data.F32[0], params[0].data.F32[1], (dparams[0].data.F32[1]/params[0].data.F32[1]));
-        for (j = 0; j < model->params->n - 4; j++) {
-            fprintf (f, "%9.3g ", params[0].data.F32[j+4]);
-        }
-        for (j = 0; j < model->params->n - 4; j++) {
-            fprintf (f, "%9.3g ", dparams[0].data.F32[j+4]);
+                 params[0].data.F32[0], -2.5*log10(params[0].data.F32[1]), (dparams[0].data.F32[1]/params[0].data.F32[1]));
+        for (j = 0; j < model->params->n - 4; j++) {
+            fprintf (f, "%9.6f ", params[0].data.F32[j+4]);
+        }
+        for (j = 0; j < model->params->n - 4; j++) {
+            fprintf (f, "%9.6f ", dparams[0].data.F32[j+4]);
         }
         fprintf (f, ": %2d %7.4f %7.4f %7.2f (%d %d)\n", source[0].type, log10(model[0].chisq), source[0].moments->SN, model[0].radius, model[0].nDOF, model[0].nIter);
+    }
+    fclose (f);
+    return true;
+}
+
+// dump the sources to an output file
+bool DumpModelFLT (psArray *sources, char *filename) 
+{
+
+    int i, j;
+    FILE *f;
+    psVector *params;
+    psVector *dparams;
+    psModel  *model;
+
+    f = fopen (filename, "w");
+    if (f == NULL) {
+        psLogMsg ("DumpObjects", 3, "can't open output file for moments%s\n", filename);
+        return false;
+    }
+
+    // write sources with models first 
+    for (i = 0; i < sources->n; i++) {
+        psSource *source = (psSource *) sources->data[i];
+        model = (psModel  *) source->modelFLT;
+        if (model == NULL) continue;
+        params = model->params;
+        dparams = model->dparams;
+        fprintf (f, "%7.1f %7.1f  %5.1f %7.3f  %7.4f  ", 
+                 params[0].data.F32[2], params[0].data.F32[3], 
+                 params[0].data.F32[0], -2.5*log10(params[0].data.F32[1]), (dparams[0].data.F32[1]/params[0].data.F32[1]));
+        for (j = 0; j < model->params->n - 4; j++) {
+            fprintf (f, "%9.6f ", params[0].data.F32[j+4]);
+        }
+        for (j = 0; j < model->params->n - 4; j++) {
+            fprintf (f, "%9.6f ", dparams[0].data.F32[j+4]);
+        }
+        fprintf (f, ": %2d %7.4f %7.4f %7.2f (%d %d)\n", source[0].type, log10(model[0].chisq), source[0].moments->SN, model[0].radius, model[0].nDOF, model[0].nIter);
+    }
+    fclose (f);
+    return true;
+}
+
+// dump the sources to an output file
+bool DumpModelNULL (psArray *sources, char *filename) 
+{
+
+    int i, j;
+    FILE *f;
+    psVector *params;
+    psVector *dparams;
+    psModel  *model;
+
+    f = fopen (filename, "w");
+    if (f == NULL) {
+        psLogMsg ("DumpObjects", 3, "can't open output file for moments%s\n", filename);
+        return false;
+    }
+
+    // write sources with models first 
+    for (i = 0; i < sources->n; i++) {
+        psSource *source = (psSource *) sources->data[i];
+        if (source->modelFLT != NULL) continue;
+        if (source->modelPSF != NULL) continue;
+        if (source->moments == NULL) continue;
+        fprintf (f, "%5d %5d  %7.1f  %7.1f %7.1f  %6.3f %6.3f  %8.1f %7.1f %7.1f %7.1f  %4d %2d\n", 
+                 source->peak->x, source->peak->y, source->peak->counts, 
+                 source->moments->x, source->moments->y, 
+                 source->moments->Sx, source->moments->Sy, 
+                 source->moments->Sum, source->moments->Peak, 
+                 source->moments->Sky, source->moments->SN, 
+                 source->moments->nPixels, source->type); 
     }
     fclose (f);
@@ -118 +190 @@
                  params[0].data.F32[0], params[0].data.F32[1]);
         for (j = 0; j < model->params->n - 4; j++) {
-            fprintf (f, "%9.3g ", params[0].data.F32[j+4]);
-        }
-        for (j = 0; j < model->params->n - 4; j++) {
-            fprintf (f, "%9.3g ", dparams[0].data.F32[j+4]);
+            fprintf (f, "%9.6f ", params[0].data.F32[j+4]);
+        }
+        for (j = 0; j < model->params->n - 4; j++) {
+            fprintf (f, "%9.6f ", dparams[0].data.F32[j+4]);
         }
         fprintf (f, ":  %7.1f (%d %d)\n", model[0].chisq, model[0].nDOF, model[0].nIter);
@@ -138 +210 @@
 }
 
+bool pmSourcePhotometry (float *fitMag, float *obsMag, psModel *model, psImage *image, psImage *mask) {
+
+    float obsSum = 0;
+    float fitSum = 0;
+    float sky = model->params->data.F32[0];
+
+    psModelFlux modelFluxFunc = psModelFlux_GetFunction (model->type);
+    fitSum = modelFluxFunc (model->params);
+
+    for (int ix = 0; ix < image->numCols; ix++) {
+        for (int iy = 0; iy < image->numRows; iy++) {
+            if (mask->data.U8[iy][ix]) continue;
+            obsSum += image->data.F32[iy][ix] - sky;
+            // fitSum += psModelEval (model, image, ix, iy) - sky;
+        }
+    }
+    if (obsSum <= 0) return false;
+    if (fitSum <= 0) return false;
+
+    *fitMag = -2.5*log10(fitSum);
+    *obsMag = -2.5*log10(obsSum);
+    return (true);
+} 
+
 # if (0)
 void pmSourceMaskSaturated (psSource *source, float saturate) {