Changeset 20718 for trunk/psModules/src/imcombine/pmSubtractionMatch.c
- Timestamp:
- Nov 12, 2008, 10:01:19 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/imcombine/pmSubtractionMatch.c
r20717 r20718 210 210 pmSubtractionStampList *stamps = NULL; // Stamps for matching PSF 211 211 pmSubtractionKernels *kernels = NULL; // Kernel basis functions 212 psMetadata *analysis = psMetadataAlloc(); // QA data 212 213 213 214 int numCols = ro1->image->numCols, numRows = ro1->image->numRows; // Image dimensions 214 215 215 216 psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS, 0); // Random number generator 216 217 {218 psStats *bg = psStatsAlloc(PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV); // Statistics for backgroun219 if (!psImageBackground(bg, NULL, ro1->image, ro1->mask, maskVal, rng)) {220 psError(PS_ERR_UNKNOWN, false, "Unable to measure background statistics.");221 psFree(bg);222 psFree(rng);223 goto MATCH_ERROR;224 }225 threshold = bg->robustMedian + threshold * bg->robustStdev;226 psFree(bg);227 }228 217 229 218 memCheck("start"); … … 233 222 if (!subMask) { 234 223 psError(PS_ERR_UNKNOWN, false, "Unable to generate subtraction mask."); 235 psFree(weight); 236 return false; 224 goto MATCH_ERROR; 237 225 } 238 226 … … 250 238 } 251 239 252 psMetadata *analysis = psMetadataAlloc(); // QA data 240 { 241 psStats *bg = psStatsAlloc(PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV); // Statistics for backgroun 242 if (!psImageBackground(bg, NULL, ro1->image, ro1->mask, maskVal, rng)) { 243 psError(PS_ERR_UNKNOWN, false, "Unable to measure background statistics."); 244 psFree(bg); 245 psFree(rng); 246 goto MATCH_ERROR; 247 } 248 threshold = bg->robustMedian + threshold * bg->robustStdev; 249 psFree(bg); 250 } 253 251 254 252 // Iterate over iso-kernel regions
Note:
See TracChangeset
for help on using the changeset viewer.
