Changeset 16836
- Timestamp:
- Mar 5, 2008, 12:54:48 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/ppStack/src/ppStackLoop.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppStack/src/ppStackLoop.c
r16717 r16836 162 162 static bool writeImage(const char *name, // Name of image 163 163 psMetadata *header, // Header 164 const psImage *image // Image 164 const psImage *image, // Image 165 pmConfig *config // Configuration 165 166 ) 166 167 { … … 168 169 assert(image); 169 170 170 psFits *fits = psFitsOpen(name, "w"); 171 psString resolved = pmConfigConvertFilename(name, config, true); // Resolved file name 172 psFits *fits = psFitsOpen(resolved, "w"); 171 173 if (!fits) { 172 psError(PS_ERR_IO, false, "Unable to open FITS file to write image."); 174 psError(PS_ERR_IO, false, "Unable to open FITS file %s to write image.", resolved); 175 psFree(resolved); 173 176 return false; 174 177 } 175 178 if (!psFitsWriteImage(fits, header, image, 0, NULL)) { 176 psError(PS_ERR_IO, false, "Unable to write FITS image."); 179 psError(PS_ERR_IO, false, "Unable to write FITS image %s.", resolved); 180 psFitsClose(fits); 181 psFree(resolved); 177 182 return false; 178 183 } 179 184 psFitsClose(fits); 185 psFree(resolved); 180 186 return true; 181 187 } … … 344 350 pmHDU *hdu = readout->parent->parent->parent->hdu; // HDU for convolved image 345 351 assert(hdu); 346 writeImage(imageNames->data[i], hdu->header, readout->image );347 writeImage(maskNames->data[i], hdu->header, readout->mask );348 writeImage(weightNames->data[i], hdu->header, readout->weight );352 writeImage(imageNames->data[i], hdu->header, readout->image, config); 353 writeImage(maskNames->data[i], hdu->header, readout->mask, config); 354 writeImage(weightNames->data[i], hdu->header, readout->weight, config); 349 355 #endif 350 356 … … 389 395 psArray *weightFits = psArrayAlloc(num); 390 396 for (int i = 0; i < num; i++) { 391 imageFits->data[i] = psFitsOpen(imageNames->data[i], "r"); 392 maskFits->data[i] = psFitsOpen(maskNames->data[i], "r"); 393 weightFits->data[i] = psFitsOpen(weightNames->data[i], "r"); 397 // Resolved names 398 psString imageResolved = pmConfigConvertFilename(imageNames->data[i], config, false); 399 psString maskResolved = pmConfigConvertFilename(maskNames->data[i], config, false); 400 psString weightResolved = pmConfigConvertFilename(weightNames->data[i], config, false); 401 imageFits->data[i] = psFitsOpen(imageResolved, "r"); 402 maskFits->data[i] = psFitsOpen(maskResolved, "r"); 403 weightFits->data[i] = psFitsOpen(weightResolved, "r"); 404 psFree(imageResolved); 405 psFree(maskResolved); 406 psFree(weightResolved); 407 if (!imageFits->data[i] || !maskFits->data[i] || !weightFits->data[i]) { 408 psError(PS_ERR_UNKNOWN, false, "Unable to open convolved files %s, %s, %s", 409 (char*)imageNames->data[i], (char*)maskNames->data[i], (char*)weightNames->data[i]); 410 psFree(cells); 411 psFree(subKernels); 412 psFree(subRegions); 413 return false; 414 } 394 415 } 395 416
Note:
See TracChangeset
for help on using the changeset viewer.
