Changeset 12445 for trunk/psModules/src/objects
- Timestamp:
- Mar 15, 2007, 10:44:25 AM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/objects/pmSourceIO_PS1_DEV_0.c
r12402 r12445 3 3 * @author EAM, IfA 4 4 * 5 * @version $Revision: 1. 1$ $Name: not supported by cvs2svn $6 * @date $Date: 2007-03-1 1 18:56:38$5 * @version $Revision: 1.2 $ $Name: not supported by cvs2svn $ 6 * @date $Date: 2007-03-15 20:44:25 $ 7 7 * 8 8 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 45 45 psArray *table; 46 46 psMetadata *row; 47 int i , type;47 int i; 48 48 psF32 *PAR, *dPAR; 49 float lsky;50 bool status;51 49 psEllipseShape shape; 52 50 psEllipseAxes axes; 53 54 // find config information for output header55 float ZERO_POINT = psMetadataLookupF32 (&status, header, "ZERO_PT");56 51 57 52 table = psArrayAllocEmpty (sources->n); … … 66 61 dPAR = model->dparams->data.F32; 67 62 68 type = pmSourceDophotType (source);69 lsky = (source->sky < 1.0) ? 0.0 : log10(source->sky);70 71 63 shape.sx = PAR[PM_PAR_SXX]; 72 64 shape.sy = PAR[PM_PAR_SYY]; … … 74 66 axes = psEllipseShapeToAxes (shape, 20.0); 75 67 76 peakMag = (source->peak->flux > 0) ? -2.5*log(source->peak->flux) : NAN;77 nImageOverlap = 1;78 ID = 0; // XXX need to figure out how to generate this68 float peakMag = (source->peak->flux > 0) ? -2.5*log(source->peak->flux) : NAN; 69 psS16 nImageOverlap = 1; 70 psS32 ID = 0; // XXX need to figure out how to generate this 79 71 80 72 row = psMetadataAlloc (); 81 psMetadataAdd (row, PS_LIST_TAIL, "IPP_IDET", PS_DATA_ S64, "IPP detection identifier index", ID);73 psMetadataAdd (row, PS_LIST_TAIL, "IPP_IDET", PS_DATA_U32, "IPP detection identifier index", ID); 82 74 psMetadataAdd (row, PS_LIST_TAIL, "X_PSF", PS_DATA_F32, "PSF x coordinate", PAR[PM_PAR_XPOS]); 83 75 psMetadataAdd (row, PS_LIST_TAIL, "Y_PSF", PS_DATA_F32, "PSF y coordinate", PAR[PM_PAR_YPOS]); … … 99 91 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight); 100 92 // XXX not sure how to get this : need to load Nimages with weight 101 psMetadataAdd (row, PS_LIST_TAIL, "N_FRAMES", PS_DATA_S32, "Number of frames overlapping source center", nImageOverlap); 93 psMetadataAdd (row, PS_LIST_TAIL, "N_FRAMES", PS_DATA_U16, "Number of frames overlapping source center", nImageOverlap); 94 psMetadataAdd (row, PS_LIST_TAIL, "DUMMY", PS_DATA_U16, "padding", 0); 102 95 103 96 // XXX these calibrated values are not supplied by psphot analysis … … 138 131 psEllipseShape shape; 139 132 psEllipseAxes axes; 140 float lsky;141 133 142 134 // define PSF model type … … 147 139 modelType = pmModelSetType (PSF_NAME); 148 140 } 149 150 // find config information for output header151 float ZERO_POINT = psMetadataLookupF32 (&status, header, "ZERO_PT");152 if (!status)153 ZERO_POINT = 25.0;154 141 155 142 psArray *table = psFitsReadTable (fits);
Note:
See TracChangeset
for help on using the changeset viewer.
