Changeset 26259 for trunk/psastro/src/psastroConvert.c

Timestamp:

Nov 22, 2009, 2:57:41 PM (16 years ago)

Author:

eugene

Message:

various fixes to psastro:

1) added bootstrap resampling to zero point error analysis
2) added iterative clump removal from refstars and rawstars
3) added unique reference match option
4) some improved visualizations
5) improved mosaic iterations

File:

• trunk/psastro/src/psastroConvert.c (modified) (2 diffs)

trunk/psastro/src/psastroConvert.c

@@ -127 +127 @@
 
     psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.RAWSTARS", PS_DATA_ARRAY, "astrometry objects", rawStars);
+
     psLogMsg ("psastro", 4, "loaded %ld sources, using %ld of %ld good sources (inst mag: %f to %f)\n", sources->n, rawStars->n, inStars->n, mMin, mMax);
     psLogMsg ("psastro", 4, "skip reasons: mode: %d, faint: %d, bright: %d, inf: %d\n", nModeSkip, nFaintSkip, nBrightSkip, nInfSkip);
@@ -153 +154 @@
 
         psF32 *PAR = model->params->data.F32;
+        psF32 *dPAR = model->dparams->data.F32;
 
         pmAstromObj *obj = pmAstromObjAlloc ();
 
         // is the source magnitude calibrated in any sense?
-        obj->pix->x = PAR[PM_PAR_XPOS];
-        obj->pix->y = PAR[PM_PAR_YPOS];
+        obj->pix->x    = PAR[PM_PAR_XPOS];
+        obj->pix->y    = PAR[PM_PAR_YPOS];
+        obj->pix->xErr = dPAR[PM_PAR_XPOS];
+        obj->pix->yErr = dPAR[PM_PAR_YPOS];
         obj->Mag = source->psfMag;
+        obj->dMag = source->errMag;
 
         // XXX do we have the information giving the readout and cell offset?