Changeset 40483 for branches/czw_branch/20170908/ppTranslate
- Timestamp:
- Jun 27, 2018, 3:31:55 PM (8 years ago)
- Location:
- branches/czw_branch/20170908
- Files:
-
- 5 edited
-
. (modified) (1 prop)
-
ppTranslate/src/ppMops.h (modified) (2 diffs)
-
ppTranslate/src/ppMopsGetSkyChipPsfVersion.c (modified) (1 diff)
-
ppTranslate/src/ppMopsRead.c (modified) (4 diffs)
-
ppTranslate/src/ppMopsWrite.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/czw_branch/20170908
- Property svn:mergeinfo changed
-
branches/czw_branch/20170908/ppTranslate/src/ppMops.h
r40049 r40483 11 11 PM_SOURCE_MODE_CR_LIMIT | PM_SOURCE_MODE_SKY_FAILURE | PM_SOURCE_MODE_DEFECT) // Flags to exclude 12 12 //#define SOURCE_MASK2 (PM_SOURCE_MODE2_DIFF_WITH_DOUBLE | PM_SOURCE_MODE2_ON_SPIKE | PM_SOURCE_MODE2_ON_STARCORE | PM_SOURCE_MODE2_ON_BURNTOOL | PM_SOURCE_MODE2_ON_CONVPOOR) // Flags2 to exclude 13 14 #define SOURCE_MASK2 (PM_SOURCE_MODE2_DIFF_WITH_DOUBLE | PM_SOURCE_MODE2_ON_BURNTOOL | PM_SOURCE_MODE2_ON_STARCORE | \ 15 PM_SOURCE_MODE2_ON_SPIKE) 13 #define SOURCE_MASK2 (PM_SOURCE_MODE2_ON_STARCORE) //Flags2 to exclude 16 14 17 15 #define PS1_DV_FORMAT "PS1_DV%d" … … 68 66 psString fpashutoutc; // FPA shutoutc 69 67 psString fpashutcutc; // FPA shutcutc 68 psString refcat; // Reference catalog used 70 69 } ppMopsDetections; 71 70 -
branches/czw_branch/20170908/ppTranslate/src/ppMopsGetSkyChipPsfVersion.c
r34472 r40483 17 17 psFree(headerSkyChip); 18 18 return 3; 19 } else if (strcmp(version, "PS1_DV4") == 0) { 20 psFree(headerSkyChip); 21 return 4; 22 } else if (strcmp(version, "PS1_DV5") == 0) { 23 psFree(headerSkyChip); 24 return 5; 19 25 } 20 26 psWarning("Unsupported EXTTYPE in SkyChip.psf table: [%s]", version); -
branches/czw_branch/20170908/ppTranslate/src/ppMopsRead.c
r40058 r40483 107 107 det->fpashutoutc = psMemIncrRefCounter(psMetadataLookupStr(NULL, header, "FPA.SHUTOUTC")); 108 108 det->fpashutcutc = psMemIncrRefCounter(psMetadataLookupStr(NULL, header, "FPA.SHUTCUTC")); 109 det->refcat = psMemIncrRefCounter(psMetadataLookupStr(NULL, header, "PSREFCAT")); 109 110 psFree(header); 110 111 … … 242 243 243 244 //MEH -- diff params, may need skyChipPsfVersion 244 psVector *npos = psMetadataLookupVector(NULL, table, "DIFF_NPOS"); 245 // -- unclear for npos, remove for now 246 //psVector *npos = psMetadataLookupVector(NULL, table, "DIFF_NPOS"); 245 247 psVector *fpos = psMetadataLookupVector(NULL, table, "DIFF_FRATIO"); 246 248 psVector *rbad = psMetadataLookupVector(NULL, table, "DIFF_NRATIO_BAD"); … … 270 272 } 271 273 272 // add diff params to discard -- needs to be in config.. -- WSdiff only?274 // MEH -- add diff params to discard -- needs to be in config.. -- WSdiff only? 273 275 // -- needs to be finite to test but unclear if should reject (probably) 274 276 // -- using Napoli 2013 to start with … … 278 280 // DIFF_NRATIO_MASK>0.4 279 281 // DIFF_NRATIO_ALL>0.3 280 // -- unclear about npos still -- 282 // -- unclear about npos still -- so remove from check 281 283 //if (npos->data.S32[row]<=3 || 282 284 if (!isfinite(fpos->data.F32[row]) || !isfinite(rbad->data.F32[row]) || !isfinite(rmask->data.F32[row]) || !isfinite(rall->data.F32[row])) { 283 psTrace("ppMops.read", 10, "Discarding row %ld from input %d because of NAN diff params: 3, 0.6, 0.4, 0.4, 0.3: %d %g %f %f %f ", row, i, npos->data.S32[row],fpos->data.F32[row],rbad->data.F32[row],rmask->data.F32[row],rall->data.F32[row]); 285 //psTrace("ppMops.read", 10, "Discarding row %ld from input %d because of NAN diff params: 3, 0.6, 0.4, 0.4, 0.3: %d %g %f %f %f ", row, i, npos->data.S32[row],fpos->data.F32[row],rbad->data.F32[row],rmask->data.F32[row],rall->data.F32[row]); 286 psTrace("ppMops.read", 10, "Discarding row %ld from input %d because of NAN diff params: 3, 0.6, 0.4, 0.4, 0.3: %g %f %f %f ", row, i, fpos->data.F32[row],rbad->data.F32[row],rmask->data.F32[row],rall->data.F32[row]); 284 287 det->mask->data.U8[row] = 0xFF; 285 288 continue; 286 289 } 287 290 if (fpos->data.F32[row]<=0.6 || rbad->data.F32[row]<=0.4 || rmask->data.F32[row]<=0.4 || rall->data.F32[row]<=0.3) { 288 psTrace("ppMops.read", 10, "Discarding row %ld from input %d because of diff params: 3, 0.6, 0.4, 0.4, 0.3: %d %g %f %f %f ",289 row, i, npos->data.S32[row], fpos->data.F32[row], rbad->data.F32[row], rmask->data.F32[row], rall->data.F32[row]);291 //psTrace("ppMops.read", 10, "Discarding row %ld from input %d because of diff params: 3, 0.6, 0.4, 0.4, 0.3: %d %g %f %f %f ", row, i, npos->data.S32[row], fpos->data.F32[row], rbad->data.F32[row], rmask->data.F32[row], rall->data.F32[row]); 292 psTrace("ppMops.read", 10, "Discarding row %ld from input %d because of diff params: 3, 0.6, 0.4, 0.4, 0.3: %g %f %f %f ", row, i, fpos->data.F32[row], rbad->data.F32[row], rmask->data.F32[row], rall->data.F32[row]); 290 293 det->mask->data.U8[row] = 0xFF; 291 294 continue; -
branches/czw_branch/20170908/ppTranslate/src/ppMopsWrite.c
r38729 r40483 104 104 psMetadataAddStr(header, PS_LIST_TAIL, "FPA.SHUTOUTC", 0, "Time of exposure open (CHIP mod)", det->fpashutoutc); 105 105 psMetadataAddStr(header, PS_LIST_TAIL, "FPA.SHUTCUTC", 0, "Time of exposure close (RAW)", det->fpashutcutc); 106 psMetadataAddStr(header, PS_LIST_TAIL, "PSREFCAT", 0, "Reference catalog used", det->refcat); 106 107 107 108 //field in header that tells about the CMF version
Note:
See TracChangeset
for help on using the changeset viewer.
