Changeset 27989 for trunk/pstamp
- Timestamp:
- May 17, 2010, 5:49:51 PM (16 years ago)
- Location:
- trunk/pstamp/src
- Files:
-
- 2 edited
-
ppstampMakeStamp.c (modified) (3 diffs)
-
ppstampMosaic.c (modified) (1 diff)
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; -
trunk/pstamp/src/ppstampMosaic.c
r18712 r27989 24 24 pmChip *inChip = pmFPAviewThisChip(view, input->fpa); 25 25 if (!mChip->hdu && !mChip->parent->hdu) { 26 const char *name = psMetadataLookupStr(&status, input->fpa->concepts, "FPA.OBS"); // Name of FPA 27 pmFPAAddSourceFromView(mosaic->fpa, name, mosaicView, mosaic->format); 26 pmFPAAddSourceFromView(mosaic->fpa, mosaicView, mosaic->format); 28 27 } 29 28 psFree(mosaicView);
Note:
See TracChangeset
for help on using the changeset viewer.
