Changeset 17254
- Timestamp:
- Mar 31, 2008, 2:19:39 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/opihi/dvo/dbExtractMeasures.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/opihi/dvo/dbExtractMeasures.c
r17199 r17254 78 78 Image *image; 79 79 Coords *mosaic; 80 PhotCode *equiv; 80 81 81 82 value = NAN; … … 83 84 switch (field->ID) { 84 85 case MEAS_MAG: /* magnitudes are already determined above */ 86 equiv = GetPhotcodeEquivbyCode (measure[0].photcode); 87 85 88 // we return the magnitude for this measure if: 86 89 if (field->photcode->type == PHOT_MAG) goto valid_photcode; 87 90 if ((field->photcode->type == PHOT_REF) && (measure[0].photcode == field->photcode->code)) goto valid_photcode; 88 91 if ((field->photcode->type == PHOT_DEP) && (measure[0].photcode == field->photcode->code)) goto valid_photcode; 89 if ((field->photcode->type == PHOT_SEC) && (GetPhotcodeEquivCodebyCode (measure[0].photcode) == field->photcode->code)) goto valid_photcode; 92 93 if ((equiv != NULL) && (field->photcode->type == PHOT_SEC) && (equiv[0].code == field->photcode->code)) goto valid_photcode; 90 94 break; 91 95 … … 105 109 break; 106 110 case MAG_CAL: 107 value = PhotCal (measure, average, secfilt, measure, GetPhotcodeEquivbyCode (measure[0].photcode));111 value = PhotCal (measure, average, secfilt, measure, equiv); 108 112 break; 109 113 case MAG_AVE: 110 value = PhotAve ( GetPhotcodeEquivbyCode (measure[0].photcode), average, secfilt);114 value = PhotAve (equiv, average, secfilt); 111 115 break; 112 116 case MAG_REF: 113 value = PhotRef ( GetPhotcodeEquivbyCode (measure[0].photcode), average, secfilt, measure);117 value = PhotRef (equiv, average, secfilt, measure); 114 118 break; 115 119 case MAG_ERR: 116 120 if (field->photcode->type == PHOT_SEC) { 117 value = PhotdM ( GetPhotcodeEquivbyCode (measure[0].photcode), average, secfilt);121 value = PhotdM (equiv, average, secfilt); 118 122 } else { 119 123 value = measure[0].dM; … … 121 125 break; 122 126 case MAG_CHISQ: 123 value = PhotXm ( GetPhotcodeEquivbyCode (measure[0].photcode), average, secfilt);127 value = PhotXm (equiv, average, secfilt); 124 128 break; 125 129 case MAG_NCODE: 126 Nsec = GetPhotcodeNsec (GetPhotcodeEquivCodebyCode (measure[0].photcode)); 130 if (equiv == NULL) { 131 value = NAN; 132 break; 133 } 134 Nsec = GetPhotcodeNsec (equiv); 127 135 if (Nsec == -1) break; 128 136 value = secfilt[Nsec].Ncode; 129 137 break; 130 138 case MAG_NPHOT: 131 Nsec = GetPhotcodeNsec (GetPhotcodeEquivCodebyCode (measure[0].photcode)); 139 if (equiv == NULL) { 140 value = NAN; 141 break; 142 } 143 Nsec = GetPhotcodeNsec (equiv); 132 144 if (Nsec == -1) break; 133 145 value = secfilt[Nsec].Nused;
Note:
See TracChangeset
for help on using the changeset viewer.
