Changeset 27989 for trunk/pstamp/src/ppstampMakeStamp.c
- Timestamp:
- May 17, 2010, 5:49:51 PM (16 years ago)
- File:
-
- 1 edited
-
trunk/pstamp/src/ppstampMakeStamp.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/pstamp/src/ppstampMakeStamp.c
r27461 r27989 36 36 double d_before = check->sky->d; 37 37 psLogMsg("ppstampMakeStamp", 2, "Before fit Center Pixel RA: %f DEC: %f (degrees)\n", 38 r_before * PS_DEG_RAD, d_before * PS_DEG_RAD);38 r_before * PS_DEG_RAD, d_before * PS_DEG_RAD); 39 39 40 40 if (outFPA->toSky->type != PS_PROJ_TAN) { 41 // we need to make astrometry terms for the output which are centered on the output chip center,42 // but we keep the original plate scale43 psFree(outFPA->toSky);44 outFPA->toSky = psProjectionAlloc (r_before, d_before, inFPA->toSky->Xs, inFPA->toSky->Ys, PS_PROJ_TAN);45 46 if (!pmAstromLinearizeToSky(inFPA, inChip, outFPA, outChip, roi)) {41 // we need to make astrometry terms for the output which are centered on the output chip center, 42 // but we keep the original plate scale 43 psFree(outFPA->toSky); 44 outFPA->toSky = psProjectionAlloc (r_before, d_before, inFPA->toSky->Xs, inFPA->toSky->Ys, PS_PROJ_TAN); 45 46 if (!pmAstromLinearizeToSky(inFPA, inChip, outFPA, outChip, roi)) { 47 47 psFree(check); 48 48 psError(PS_ERR_UNKNOWN, false, "Failed to linearize astrometry\n"); … … 56 56 57 57 psLogMsg("ppstampMakeStamp", 2, "After fit: Center Pixel RA: %f DEC: %f (degrees)\n", 58 r_after * PS_DEG_RAD, d_after * PS_DEG_RAD);59 psLogMsg("ppstampMakeStamp", 2, "Error in fit to astrometry %.2f %.2f (arcseconds)\n", 60 (r_after - r_before) * PS_DEG_RAD * 3600, (d_after - d_before) * PS_DEG_RAD * 3600);58 r_after * PS_DEG_RAD, d_after * PS_DEG_RAD); 59 psLogMsg("ppstampMakeStamp", 2, "Error in fit to astrometry %.2f %.2f (arcseconds)\n", 60 (r_after - r_before) * PS_DEG_RAD * 3600, (d_after - d_before) * PS_DEG_RAD * 3600); 61 61 62 62 // XXX: should we fail if the fit is bad ?? … … 238 238 return PS_EXIT_DATA_ERROR; 239 239 } 240 char *fpaName = psMetadataLookupStr(NULL, input->fpa->concepts, "FPA.OBS"); // Name of FPA241 240 pmFPAview *outview = pmFPAviewAlloc(0); 242 pmFPAAddSourceFromView(output->fpa, fpaName,outview, output->format);241 pmFPAAddSourceFromView(output->fpa, outview, output->format); 243 242 244 243 outview->chip = 0;
Note:
See TracChangeset
for help on using the changeset viewer.
