Changeset 20307 for trunk/psModules
- Timestamp:
- Oct 21, 2008, 4:11:36 PM (18 years ago)
- Location:
- trunk/psModules/src
- Files:
-
- 3 edited
-
camera/pmFPAMaskWeight.c (modified) (1 diff)
-
objects/pmModel.c (modified) (3 diffs)
-
objects/pmSourcePhotometry.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/camera/pmFPAMaskWeight.c
r20095 r20307 453 453 psImage *weight = readout->weight; // Weight map 454 454 455 psImageInterpolat eOptions *interp = psImageInterpolateOptionsAlloc(mode, image, weight, mask, maskVal,456 NAN, NAN, maskBad, maskPoor, poorFrac);455 psImageInterpolation *interp = psImageInterpolationAlloc(mode, image, weight, mask, maskVal, 456 NAN, NAN, maskBad, maskPoor, poorFrac, 0); 457 457 interp->shifting = false; // Turn off "exact shifts" so we get proper interpolation 458 458 -
trunk/psModules/src/objects/pmModel.c
r15977 r20307 6 6 * @author EAM, IfA 7 7 * 8 * @version $Revision: 1.2 2$ $Name: not supported by cvs2svn $9 * @date $Date: 2008- 01-02 20:36:36$8 * @version $Revision: 1.23 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2008-10-22 02:11:08 $ 10 10 * 11 11 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 222 222 psImage *myRx = NULL; 223 223 psImage *myRy = NULL; 224 psImageInterpolat eOptions*Ro = NULL;225 psImageInterpolat eOptions*Rx = NULL;226 psImageInterpolat eOptions*Ry = NULL;224 psImageInterpolation *Ro = NULL; 225 psImageInterpolation *Rx = NULL; 226 psImageInterpolation *Ry = NULL; 227 227 if (model->residuals && (mode & (PM_MODEL_OP_RES0 | PM_MODEL_OP_RES1))) { 228 228 // if the residual image and object image don't match, … … 274 274 } 275 275 276 Ro = psImageInterpolat eOptionsAlloc(PS_INTERPOLATE_BILINEAR, myRo, NULL, mask, 0, 0.0, 0.0, 1, 0, 0.0);277 Rx = psImageInterpolat eOptionsAlloc(PS_INTERPOLATE_BILINEAR, myRx, NULL, NULL, 0, 0.0, 0.0, 1, 0, 0.0);278 Ry = psImageInterpolat eOptionsAlloc(PS_INTERPOLATE_BILINEAR, myRy, NULL, NULL, 0, 0.0, 0.0, 1, 0, 0.0);276 Ro = psImageInterpolationAlloc(PS_INTERPOLATE_BILINEAR, myRo, NULL, mask, 0, 0.0, 0.0, 1, 0, 0.0, 0); 277 Rx = psImageInterpolationAlloc(PS_INTERPOLATE_BILINEAR, myRx, NULL, NULL, 0, 0.0, 0.0, 1, 0, 0.0, 0); 278 Ry = psImageInterpolationAlloc(PS_INTERPOLATE_BILINEAR, myRy, NULL, NULL, 0, 0.0, 0.0, 1, 0, 0.0, 0); 279 279 280 280 } -
trunk/psModules/src/objects/pmSourcePhotometry.c
r20077 r20307 3 3 * @author EAM, IfA; GLG, MHPCC 4 4 * 5 * @version $Revision: 1.4 5$ $Name: not supported by cvs2svn $6 * @date $Date: 2008-10- 13 01:59:05$5 * @version $Revision: 1.46 $ $Name: not supported by cvs2svn $ 6 * @date $Date: 2008-10-22 02:11:08 $ 7 7 * 8 8 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 122 122 if (source->modelFits) { 123 123 for (int i = 0; i < source->modelFits->n; i++) { 124 pmModel *model = source->modelFits->data[i];125 status = pmSourcePhotometryModel (&model->mag, model);124 pmModel *model = source->modelFits->data[i]; 125 status = pmSourcePhotometryModel (&model->mag, model); 126 126 } 127 127 if (source->modelEXT) { 128 source->extMag = source->modelEXT->mag;128 source->extMag = source->modelEXT->mag; 129 129 } 130 130 } else { 131 131 if (source->modelEXT) { 132 status = pmSourcePhotometryModel (&source->extMag, source->modelEXT);132 status = pmSourcePhotometryModel (&source->extMag, source->modelEXT); 133 133 } 134 134 } … … 177 177 178 178 if (!psImageShiftMask(&flux, &mask, source->pixels, source->maskObj, maskVal, dx, dy, 179 NAN, 0xff, PS_INTERPOLATE_BI CUBE)) {179 NAN, 0xff, PS_INTERPOLATE_BIQUADRATIC)) { 180 180 // Not much we can do about it 181 181 psErrorClear(); … … 256 256 } 257 257 if (mode & PM_SOURCE_PHOT_APCORR) { 258 rflux = pow (10.0, 0.4*source->psfMag);259 source->apMag -= PS_SQR(model->radiusFit)*rflux * psf->skyBias + psf->skySat / rflux;258 rflux = pow (10.0, 0.4*source->psfMag); 259 source->apMag -= PS_SQR(model->radiusFit)*rflux * psf->skyBias + psf->skySat / rflux; 260 260 } 261 261 } … … 333 333 } 334 334 if (apSum <= 0) { 335 *apMag = NAN;335 *apMag = NAN; 336 336 return true; 337 337 }
Note:
See TracChangeset
for help on using the changeset viewer.
