- Timestamp:
- Oct 18, 2009, 10:23:28 AM (17 years ago)
- Location:
- branches/eam_branches/20090820
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
psModules/src/objects/pmPSF_IO.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/20090820
-
branches/eam_branches/20090820/psModules/src/objects/pmPSF_IO.c
r21183 r25870 321 321 pmPSFClump psfClump; 322 322 323 324 int nRegions = psMetadataLookupS32 (&status, recipe, "PSF.CLUMP.NREGIONS");325 326 327 328 329 psMetadata *regionMD = psMetadataLookupPtr (&status, recipe, regionName);330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 323 // we now save clump parameters for each region : need to save all of those 324 int nRegions = psMetadataLookupS32 (&status, analysis, "PSF.CLUMP.NREGIONS"); 325 psMetadataAddS32 (header, PS_LIST_TAIL, "PSF_CLN", PS_META_REPLACE, "number of psf clump regions", nRegions); 326 for (int i = 0; i < nRegions; i++) { 327 char regionName[64]; 328 snprintf (regionName, 64, "PSF.CLUMP.REGION.%03d", i); 329 psMetadata *regionMD = psMetadataLookupPtr (&status, analysis, regionName); 330 331 psfClump.X = psMetadataLookupF32 (&status, regionMD, "PSF.CLUMP.X"); assert (status); 332 psfClump.Y = psMetadataLookupF32 (&status, regionMD, "PSF.CLUMP.Y"); assert (status); 333 psfClump.dX = psMetadataLookupF32 (&status, regionMD, "PSF.CLUMP.DX"); assert (status); 334 psfClump.dY = psMetadataLookupF32 (&status, regionMD, "PSF.CLUMP.DY"); assert (status); 335 336 char key[16]; 337 snprintf (key, 9, "CLX_%03d", i); 338 psMetadataAddF32 (header, PS_LIST_TAIL, key, PS_META_REPLACE, "psf clump center", psfClump.X); 339 snprintf (key, 9, "CLY_%03d", i); 340 psMetadataAddF32 (header, PS_LIST_TAIL, key, PS_META_REPLACE, "psf clump center", psfClump.Y); 341 snprintf (key, 9, "CLDX_%03d", i); 342 psMetadataAddF32 (header, PS_LIST_TAIL, key, PS_META_REPLACE, "psf clump size", psfClump.dX); 343 snprintf (key, 9, "CLDY_%03d", i); 344 psMetadataAddF32 (header, PS_LIST_TAIL, key, PS_META_REPLACE, "psf clump size", psfClump.dY); 345 } 346 346 347 347 // save the dimensions of each parameter … … 667 667 int nRegions = psMetadataLookupS32 (&status, header, "PSF_CLN"); 668 668 if (!status) { 669 670 671 672 673 psMetadataAddS32 (recipe, PS_LIST_TAIL, "PSF.CLUMP.NREGIONS", PS_META_REPLACE, "psf clump regions", 1);674 675 psMetadata *regionMD = psMetadataLookupPtr (&status, recipe, regionName);676 677 678 psMetadataAddMetadata (recipe, PS_LIST_TAIL, regionName, PS_META_REPLACE, "psf clump region", regionMD);679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 669 // read old-style psf clump data 670 671 char regionName[64]; 672 snprintf (regionName, 64, "PSF.CLUMP.REGION.000"); 673 psMetadataAddS32 (analysis, PS_LIST_TAIL, "PSF.CLUMP.NREGIONS", PS_META_REPLACE, "psf clump regions", 1); 674 675 psMetadata *regionMD = psMetadataLookupPtr (&status, analysis, regionName); 676 if (!regionMD) { 677 regionMD = psMetadataAlloc(); 678 psMetadataAddMetadata (analysis, PS_LIST_TAIL, regionName, PS_META_REPLACE, "psf clump region", regionMD); 679 psFree (regionMD); 680 } 681 682 // psf clump data 683 pmPSFClump psfClump; 684 685 psfClump.X = psMetadataLookupF32 (&status, header, "PSF_CLX" ); assert(status); 686 psfClump.Y = psMetadataLookupF32 (&status, header, "PSF_CLY" ); assert(status); 687 psfClump.dX = psMetadataLookupF32 (&status, header, "PSF_CLDX"); assert(status); 688 psfClump.dY = psMetadataLookupF32 (&status, header, "PSF_CLDY"); assert(status); 689 690 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.X", PS_META_REPLACE, "psf clump center", psfClump.X); 691 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.Y", PS_META_REPLACE, "psf clump center", psfClump.Y); 692 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.DX", PS_META_REPLACE, "psf clump center", psfClump.dX); 693 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.DY", PS_META_REPLACE, "psf clump center", psfClump.dY); 694 694 } else { 695 696 697 698 699 700 701 702 psMetadata *regionMD = psMetadataLookupPtr (&status, recipe, regionName);703 704 705 psMetadataAddMetadata (recipe, PS_LIST_TAIL, regionName, PS_META_REPLACE, "psf clump region", regionMD);706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 695 psMetadataAddS32 (recipe, PS_LIST_TAIL, "PSF.CLUMP.NREGIONS", PS_META_REPLACE, "psf clump regions", nRegions); 696 697 for (int i = 0; i < nRegions; i++) { 698 char key[10]; 699 char regionName[64]; 700 snprintf (regionName, 64, "PSF.CLUMP.REGION.%03d", i); 701 702 psMetadata *regionMD = psMetadataLookupPtr (&status, analysis, regionName); 703 if (!regionMD) { 704 regionMD = psMetadataAlloc(); 705 psMetadataAddMetadata (analysis, PS_LIST_TAIL, regionName, PS_META_REPLACE, "psf clump region", regionMD); 706 psFree (regionMD); 707 } 708 709 // psf clump data 710 pmPSFClump psfClump; 711 712 snprintf (key, 9, "CLX_%03d", i); 713 psfClump.X = psMetadataLookupF32 (&status, header, key); assert(status); 714 snprintf (key, 9, "CLY_%03d", i); 715 psfClump.Y = psMetadataLookupF32 (&status, header, key); assert(status); 716 snprintf (key, 9, "CLDX_%03d", i); 717 psfClump.dX = psMetadataLookupF32 (&status, header, key); assert(status); 718 snprintf (key, 9, "CLDY_%03d", i); 719 psfClump.dY = psMetadataLookupF32 (&status, header, key); assert(status); 720 721 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.X", PS_META_REPLACE, "psf clump center", psfClump.X); 722 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.Y", PS_META_REPLACE, "psf clump center", psfClump.Y); 723 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.DX", PS_META_REPLACE, "psf clump center", psfClump.dX); 724 psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.DY", PS_META_REPLACE, "psf clump center", psfClump.dY); 725 } 726 726 } 727 727 … … 871 871 psFitsReadImageBuffer(psf->residuals->Ry, file->fits, fullImage, 2); // Desired pixels 872 872 } 873 873 // XXX notice that we are not saving the resid->mask 874 874 } 875 875
Note:
See TracChangeset
for help on using the changeset viewer.
