Changeset 5350 for trunk/psphot/src/psphotOutput.c
- Timestamp:
- Oct 17, 2005, 8:07:54 AM (21 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotOutput.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotOutput.c
r5128 r5350 186 186 187 187 // write necessary information to output header 188 psMetadataAdd (imdata->header, PS_LIST_TAIL, "NSTARS", PS_ META_S32 | PS_META_REPLACE, "NUMBER OF STARS", sources->n);189 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ZERO_PT", PS_ META_F32 | PS_META_REPLACE, "zero point", ZERO_POINT);190 psMetadataAdd (imdata->header, PS_LIST_TAIL, "APMIFIT", PS_ META_F32 | PS_META_REPLACE, "aperture residual", psf->ApResid);191 psMetadataAdd (imdata->header, PS_LIST_TAIL, "dAPMIFIT", PS_ META_F32 | PS_META_REPLACE, "ap residual scatter", psf->dApResid);192 psMetadataAdd (imdata->header, PS_LIST_TAIL, "SKYBIAS", PS_ META_F32 | PS_META_REPLACE, "aperture sky bias", psf->skyBias);193 psMetadataAdd (imdata->header, PS_LIST_TAIL, "PHOTCODE", PS_ META_STR| PS_META_REPLACE, "photometry code", PHOTCODE);194 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_X", PS_ META_F32 | PS_META_REPLACE, "PSF FWHM X", 0.0);195 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_Y", PS_ META_F32 | PS_META_REPLACE, "PSF FWHM Y", 0.0);196 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ANGLE", PS_ META_F32 | PS_META_REPLACE, "PSF ANGLE", 0.0);197 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FSATUR", PS_ META_F32 | PS_META_REPLACE, "SATURATION MAG", 0.0);198 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FLIMIT", PS_ META_F32 | PS_META_REPLACE, "COMPLETENESS MAG", 0.0);188 psMetadataAdd (imdata->header, PS_LIST_TAIL, "NSTARS", PS_DATA_S32 | PS_META_REPLACE, "NUMBER OF STARS", sources->n); 189 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ZERO_PT", PS_DATA_F32 | PS_META_REPLACE, "zero point", ZERO_POINT); 190 psMetadataAdd (imdata->header, PS_LIST_TAIL, "APMIFIT", PS_DATA_F32 | PS_META_REPLACE, "aperture residual", psf->ApResid); 191 psMetadataAdd (imdata->header, PS_LIST_TAIL, "dAPMIFIT", PS_DATA_F32 | PS_META_REPLACE, "ap residual scatter", psf->dApResid); 192 psMetadataAdd (imdata->header, PS_LIST_TAIL, "SKYBIAS", PS_DATA_F32 | PS_META_REPLACE, "aperture sky bias", psf->skyBias); 193 psMetadataAdd (imdata->header, PS_LIST_TAIL, "PHOTCODE", PS_DATA_STRING | PS_META_REPLACE, "photometry code", PHOTCODE); 194 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_X", PS_DATA_F32 | PS_META_REPLACE, "PSF FWHM X", 0.0); 195 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_Y", PS_DATA_F32 | PS_META_REPLACE, "PSF FWHM Y", 0.0); 196 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ANGLE", PS_DATA_F32 | PS_META_REPLACE, "PSF ANGLE", 0.0); 197 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FSATUR", PS_DATA_F32 | PS_META_REPLACE, "SATURATION MAG", 0.0); 198 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FLIMIT", PS_DATA_F32 | PS_META_REPLACE, "COMPLETENESS MAG", 0.0); 199 199 200 200 // create file, write-out header … … 205 205 mdi = psMetadataLookup (imdata->header, "NAXIS"); 206 206 mdi->data.S32 = 0; 207 mdi->type = PS_ META_S32;208 209 // psMetadataAdd (imdata->header, PS_LIST_HEAD, "NAXIS", PS_ META_S32 | PS_META_REPLACE, "", 0);207 mdi->type = PS_DATA_S32; 208 209 // psMetadataAdd (imdata->header, PS_LIST_HEAD, "NAXIS", PS_DATA_S32 | PS_META_REPLACE, "", 0); 210 210 psFitsWriteHeader (imdata->header, fits); 211 211 psFree (fits); … … 275 275 276 276 // write necessary information to output header 277 psMetadataAdd (imdata->header, PS_LIST_TAIL, "NSTARS", PS_ META_S32 | PS_META_REPLACE, "NUMBER OF STARS", sources->n);278 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ZERO_PT", PS_ META_F32 | PS_META_REPLACE, "zero point", ZERO_POINT);279 psMetadataAdd (imdata->header, PS_LIST_TAIL, "APMIFIT", PS_ META_F32 | PS_META_REPLACE, "aperture residual", psf->ApResid);280 psMetadataAdd (imdata->header, PS_LIST_TAIL, "dAPMIFIT", PS_ META_F32 | PS_META_REPLACE, "ap residual scatter", psf->dApResid);281 psMetadataAdd (imdata->header, PS_LIST_TAIL, "SKYBIAS", PS_ META_F32 | PS_META_REPLACE, "aperture sky bias", psf->skyBias);282 psMetadataAdd (imdata->header, PS_LIST_TAIL, "PHOTCODE", PS_ META_STR| PS_META_REPLACE, "photometry code", PHOTCODE);283 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_X", PS_ META_F32 | PS_META_REPLACE, "PSF FWHM X", 0.0);284 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_Y", PS_ META_F32 | PS_META_REPLACE, "PSF FWHM Y", 0.0);285 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ANGLE", PS_ META_F32 | PS_META_REPLACE, "PSF ANGLE", 0.0);286 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FSATUR", PS_ META_F32 | PS_META_REPLACE, "SATURATION MAG", 0.0);287 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FLIMIT", PS_ META_F32 | PS_META_REPLACE, "COMPLETENESS MAG", 0.0);277 psMetadataAdd (imdata->header, PS_LIST_TAIL, "NSTARS", PS_DATA_S32 | PS_META_REPLACE, "NUMBER OF STARS", sources->n); 278 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ZERO_PT", PS_DATA_F32 | PS_META_REPLACE, "zero point", ZERO_POINT); 279 psMetadataAdd (imdata->header, PS_LIST_TAIL, "APMIFIT", PS_DATA_F32 | PS_META_REPLACE, "aperture residual", psf->ApResid); 280 psMetadataAdd (imdata->header, PS_LIST_TAIL, "dAPMIFIT", PS_DATA_F32 | PS_META_REPLACE, "ap residual scatter", psf->dApResid); 281 psMetadataAdd (imdata->header, PS_LIST_TAIL, "SKYBIAS", PS_DATA_F32 | PS_META_REPLACE, "aperture sky bias", psf->skyBias); 282 psMetadataAdd (imdata->header, PS_LIST_TAIL, "PHOTCODE", PS_DATA_STRING | PS_META_REPLACE, "photometry code", PHOTCODE); 283 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_X", PS_DATA_F32 | PS_META_REPLACE, "PSF FWHM X", 0.0); 284 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FWHM_Y", PS_DATA_F32 | PS_META_REPLACE, "PSF FWHM Y", 0.0); 285 psMetadataAdd (imdata->header, PS_LIST_TAIL, "ANGLE", PS_DATA_F32 | PS_META_REPLACE, "PSF ANGLE", 0.0); 286 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FSATUR", PS_DATA_F32 | PS_META_REPLACE, "SATURATION MAG", 0.0); 287 psMetadataAdd (imdata->header, PS_LIST_TAIL, "FLIMIT", PS_DATA_F32 | PS_META_REPLACE, "COMPLETENESS MAG", 0.0); 288 288 289 289 // set NAXIS to 0 … … 291 291 mdi = psMetadataLookup (imdata->header, "NAXIS"); 292 292 mdi->data.S32 = 0; 293 mdi->type = PS_ META_S32;293 mdi->type = PS_DATA_S32; 294 294 295 295 table = psArrayAlloc (sources->n); … … 310 310 311 311 row = psMetadataAlloc (); 312 psMetadataAdd (row, PS_LIST_HEAD, "X_PIX", PS_ META_F32, "", PAR[2]);313 psMetadataAdd (row, PS_LIST_HEAD, "Y_PIX", PS_ META_F32, "", PAR[3]);314 psMetadataAdd (row, PS_LIST_HEAD, "MAG_RAW", PS_ META_F32, "", source->fitMag + ZERO_POINT);315 psMetadataAdd (row, PS_LIST_HEAD, "MAG_ERR", PS_ META_F32, "", (int)(1000*dmag));316 psMetadataAdd (row, PS_LIST_HEAD, "MAG_GAL", PS_ META_F32, "", type);317 psMetadataAdd (row, PS_LIST_HEAD, "MAG_AP", PS_ META_F32, "", lsky);318 psMetadataAdd (row, PS_LIST_HEAD, "LOG_SKY", PS_ META_F32, "", 32.0); // should be 'Mgal'319 psMetadataAdd (row, PS_LIST_HEAD, "FWHM_X", PS_ META_F32, "", source->apMag + ZERO_POINT);320 psMetadataAdd (row, PS_LIST_HEAD, "FWHM_Y", PS_ META_F32, "", PAR[4]); // should be 'FHWM x'321 psMetadataAdd (row, PS_LIST_HEAD, "THETA", PS_ META_F32, "", PAR[5]); // should be 'FHWM y'312 psMetadataAdd (row, PS_LIST_HEAD, "X_PIX", PS_DATA_F32, "", PAR[2]); 313 psMetadataAdd (row, PS_LIST_HEAD, "Y_PIX", PS_DATA_F32, "", PAR[3]); 314 psMetadataAdd (row, PS_LIST_HEAD, "MAG_RAW", PS_DATA_F32, "", source->fitMag + ZERO_POINT); 315 psMetadataAdd (row, PS_LIST_HEAD, "MAG_ERR", PS_DATA_F32, "", (int)(1000*dmag)); 316 psMetadataAdd (row, PS_LIST_HEAD, "MAG_GAL", PS_DATA_F32, "", type); 317 psMetadataAdd (row, PS_LIST_HEAD, "MAG_AP", PS_DATA_F32, "", lsky); 318 psMetadataAdd (row, PS_LIST_HEAD, "LOG_SKY", PS_DATA_F32, "", 32.0); // should be 'Mgal' 319 psMetadataAdd (row, PS_LIST_HEAD, "FWHM_X", PS_DATA_F32, "", source->apMag + ZERO_POINT); 320 psMetadataAdd (row, PS_LIST_HEAD, "FWHM_Y", PS_DATA_F32, "", PAR[4]); // should be 'FHWM x' 321 psMetadataAdd (row, PS_LIST_HEAD, "THETA", PS_DATA_F32, "", PAR[5]); // should be 'FHWM y' 322 322 psMetadataAdd (row, PS_LIST_HEAD, "DOPHOT", PS_TYPE_S8, "", 0); // should be theta 323 // psMetadataAdd (header, PS_LIST_HEAD, "DUMMY", PS_ META_STR, "", NULL);323 // psMetadataAdd (header, PS_LIST_HEAD, "DUMMY", PS_DATA_STRING, "", NULL); 324 324 325 325 psArrayAdd (table, 100, row); … … 334 334 335 335 theader = psMetadataAlloc (); 336 psMetadataAdd (theader, PS_LIST_HEAD, "EXTNAME", PS_ META_STR, "extension name", "SMPFILE");336 psMetadataAdd (theader, PS_LIST_HEAD, "EXTNAME", PS_DATA_STRING, "extension name", "SMPFILE"); 337 337 338 338 unlink (filename);
Note:
See TracChangeset
for help on using the changeset viewer.
