Changeset 31336 for branches/eam_branches/ipp-20110404/psModules/src/objects/pmPSF_IO.c

Timestamp:

Apr 21, 2011, 12:15:30 PM (15 years ago)

Author:

eugene

Message:

add PSF_APERTURE and PSF_FIT_RADIUS to psf I/O file

File:

• branches/eam_branches/ipp-20110404/psModules/src/objects/pmPSF_IO.c (modified) (2 diffs)

branches/eam_branches/ipp-20110404/psModules/src/objects/pmPSF_IO.c

@@ -428 +428 @@
         psMetadataAddF32 (header, PS_LIST_TAIL, "SKY_BIAS", PS_DATA_F32, "sky bias level", psf->skyBias);
 
+        float PSF_APERTURE =  psMetadataLookupF32(&status, roAnalysis, "PSF_APERTURE");
+        if (status) {
+            psMetadataAddF32 (header, PS_LIST_TAIL, "PSF_APERTURE", PS_DATA_F32, "aperture for psf objects", PSF_APERTURE);
+        }
+        float PSF_FIT_RADIUS =  psMetadataLookupF32(&status, roAnalysis, "PSF_FIT_RADIUS");
+        if (status) {
+            psMetadataAddF32 (header, PS_LIST_TAIL, "PSF_FIT_RADIUS", PS_DATA_F32, "aperture for psf objects", PSF_FIT_RADIUS);
+        }
+
         // build a FITS table of the PSF parameters
         psArray *psfTable = psArrayAllocEmpty (100);
@@ -933 +942 @@
     // XXX can we drop this now?
     psf->skyBias   = psMetadataLookupF32 (&status, header, "SKY_BIAS");
+
+    if (roAnalysis) {
+        float PSF_APERTURE =  psMetadataLookupF32(&status, header, "PSF_APERTURE");
+        if (status) {
+            psMetadataAddF32 (roAnalysis, PS_LIST_TAIL, "PSF_APERTURE", PS_DATA_F32, "aperture for psf objects", PSF_APERTURE);
+        }
+        float PSF_FIT_RADIUS =  psMetadataLookupF32(&status, header, "PSF_FIT_RADIUS");
+        if (status) {
+            psMetadataAddF32 (roAnalysis, PS_LIST_TAIL, "PSF_FIT_RADIUS", PS_DATA_F32, "aperture for psf objects", PSF_FIT_RADIUS);
+        }
+    } else {
+        psWarning ("unable to read PSF_APERTURE or PSF_FIT_RADIUS");
+    }
 
     // read the raw table data