Changeset 6727 for trunk/psphot/src/psphotModelTest.c
- Timestamp:
- Mar 29, 2006, 11:08:05 AM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotModelTest.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotModelTest.c
r6715 r6727 1 1 # include "psphot.h" 2 # include "psEllipse.h"3 2 static char DEFAULT_MODE[] = "EXT"; 4 3 5 bool psphotModelTest (pmReadout *readout, recipe) {4 bool psphotModelTest (pmReadout *readout, psMetadata *recipe) { 6 5 7 6 bool status; … … 97 96 fprintf (stderr, "sum: %f @ (%f, %f)\n", source->moments->Sum, source->moments->x, source->moments->y); 98 97 99 EllipseMoments moments;98 psEllipseMoments moments; 100 99 moments.x2 = source->moments->Sx; 101 100 moments.y2 = source->moments->Sy; 102 101 moments.xy = source->moments->Sxy; 103 EllipseAxes axes =EllipseMomentsToAxes (moments);102 psEllipseAxes axes = psEllipseMomentsToAxes (moments); 104 103 105 104 fprintf (stderr, "axes: %f @ (%f, %f)\n", axes.theta*180/M_PI, axes.major, axes.minor); … … 137 136 138 137 // list model input shape 139 EllipseShape shape;138 psEllipseShape shape; 140 139 shape.sx = 1.4 / model->params->data.F32[4]; 141 140 shape.sy = 1.4 / model->params->data.F32[5]; 142 141 shape.sxy = model->params->data.F32[6]; 143 axes = EllipseShapeToAxes (shape);142 axes = psEllipseShapeToAxes (shape); 144 143 145 144 fprintf (stderr, "guess: %f @ (%f, %f)\n", axes.theta*180/M_PI, axes.major, axes.minor); … … 152 151 153 152 // define the pixels used for the fit 154 psImageKeepCircle (source->mask, xObj, yObj, RADIUS, "OR", P SPHOT_MASK_MARKED);153 psImageKeepCircle (source->mask, xObj, yObj, RADIUS, "OR", PM_SOURCE_MASK_MARKED); 155 154 156 155 char *fitset = psMetadataLookupStr (&status, recipe, "TEST_FIT_SET"); … … 163 162 164 163 // measure the source mags 165 pmSourcePhotometry (&fitMag, &obsMag, model, source->pixels, source->mask); 164 pmSourcePhotometryModel (&fitMag, model); 165 pmSourcePhotometryAper (&obsMag, model, source->pixels, source->mask); 166 166 fprintf (stderr, "ap: %f, fit: %f, apmifit: %f\n", obsMag, fitMag, obsMag - fitMag); 167 167
Note:
See TracChangeset
for help on using the changeset viewer.
