Index: /branches/eam_branches/ipp-20220316/fpcamera/src/fpcameraMetadataStats.c
===================================================================
--- /branches/eam_branches/ipp-20220316/fpcamera/src/fpcameraMetadataStats.c	(revision 42220)
+++ /branches/eam_branches/ipp-20220316/fpcamera/src/fpcameraMetadataStats.c	(revision 42220)
@@ -0,0 +1,96 @@
+/** @file fpcameraMetaDataStats.c
+ *
+ *  @brief
+ *
+ *  @ingroup fpcamera
+ *
+ *  @author IfA
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2009-02-07 02:03:34 $
+ *  Copyright 2009 Institute for Astronomy, University of Hawaii
+ */
+
+# include "fpcameraInternal.h"
+
+bool fpcameraMetadataStats (pmConfig *config, psMetadata *stats) {
+
+    bool status;
+
+    char *filename = psMetadataLookupStr (&status, config->arguments, "STATS");
+    if (!filename) return true; // stats not requested
+
+    // Extract statistics from the last output fpa
+    pmFPAfile *output = psMetadataLookupPtr(&status, config->files, "FPCAMERA.OUTPUT");
+
+    if (!output) {
+        psError(PS_ERR_UNEXPECTED_NULL, true, "Unable to find output file (FPCAMERA.OUTPUT).");
+        return false;
+    }
+
+    // extract stats for the complete fpa
+    pmFPAview *view = pmFPAviewAlloc(0);
+
+    // Entries in the FPCAMERA.HEADER in the FPA don't get seen by ppStats
+    psMetadata *header = psMetadataLookupMetadata(&status, output->fpa->analysis, "FPCAMERA.HEADER");
+    if (status && header) {
+        psMetadataItemSupplement(&status, stats, header, "ZPT_OBS");
+        psMetadataItemSupplement(&status, stats, header, "ZPT_ERR");
+        psMetadataItemSupplement(&status, stats, header, "NASTRO");
+    }
+
+    if (!ppStatsMetadata(stats, output->fpa, view, 0, config)) {
+        psError(PS_ERR_UNEXPECTED_NULL, false, "Unable to generate stats for image.");
+        psFree(view);
+        return false;
+    }
+
+    psFree(view);
+
+    // if we did not request any specific stats, the structure is empty
+    if (stats && stats->list->n == 0) {
+        psWarning ("stats output specified, but no requested stats entries in headers");
+        return true;
+    }
+
+    // convert the stats MDC to a string
+    char *statsMDC = psMetadataConfigFormat(stats);
+    if (!statsMDC || strlen(statsMDC) == 0) {
+        psError(PS_ERR_IO, false, "Unable to serialize stats metadata.\n");
+        return false;
+    }
+
+    // convert to a real UNIX filename
+    psString resolved = pmConfigConvertFilename(filename, config, true, false); // Resolved filename
+    if (!resolved) {
+        psError(psErrorCodeLast(), false, "Unable to resolve statistics filename: %s", filename);
+        psFree(statsMDC);
+        return false;
+    }
+    FILE *statsFile = fopen (resolved, "w");
+    if (!statsFile) {
+        psError(PS_ERR_IO, true, "Unable to open statistics file %s for writing.\n", resolved);
+        psFree(statsMDC);
+        psFree(resolved);
+        return false;
+    }
+
+    // write the stats MDC to a file
+    if (fprintf(statsFile, "%s", statsMDC) != strlen(statsMDC)) {
+        psError(PS_ERR_IO, false, "Unable to write statistics file %s", resolved);
+        psFree(statsMDC);
+        psFree(resolved);
+        return false;
+    }
+    psFree(statsMDC);
+    if (fclose(statsFile) == EOF) {
+        psError(PS_ERR_IO, false, "Unable to write statistics file %s", resolved);
+        psFree(resolved);
+        return false;
+    }
+
+    psFree(resolved);
+
+    pmConfigRunFilenameAddWrite(config, "STATS", filename);
+
+    return true;
+}