Changeset 29515 for branches/eam_branches/ipp-20100823/ippToPsps/src/ippToPspsBatchStack.c

Timestamp:

Oct 21, 2010, 2:45:13 PM (16 years ago)

Author:

eugene

Message:

merge changes from trunk

File:

• branches/eam_branches/ipp-20100823/ippToPsps/src/ippToPspsBatchStack.c (modified) (2 diffs)

branches/eam_branches/ipp-20100823/ippToPsps/src/ippToPspsBatchStack.c

@@ -9 +9 @@
 #include "ippToPsps.h"
 #include "ippToPspsConfig.h"
+#include "ippToPspsStackEnums.h"
+
 
 /**
@@ -15 +17 @@
 int ippToPsps_batchStack(IppToPsps *this) {
 
-//    char extensionName[20];
-    uint32_t skycell = 0; // TODO
+    int status = 0;
+    fitsfile *fitsIn;
 
-    for (uint16_t i = 0; i<this->numOfInputFiles; i++) {
+    if (fits_open_file(&fitsIn, this->inputFiles[0], READONLY, &status)) {
 
-        skycell = i; // TODO
+        fits_report_error(stderr, status);
+        return PS_EXIT_SYS_ERROR;
+    }
+
+    long removeList[MAXDETECT];
+    float instMag[MAXDETECT];
+    float floatnull = -999.0;
+    int anynull = 0;
+
+    // get primary header and pull stuff out for later
+    int nKeys = 0;
+    fits_get_hdrspace(fitsIn, &nKeys, NULL, &status);
+    float exposureTime; status=0; fits_read_key(fitsIn, TFLOAT, "EXPTIME", &exposureTime, NULL, &status);
+    char filterType[20]; status=0; fits_read_key(fitsIn, TSTRING, "FPA.FILTERID", filterType, NULL, &status);
 
 
+    char** assocDate = (char**)calloc(MAXDETECT, sizeof(char**));
+    for (uint32_t i=0; i<MAXDETECT;i++) assocDate[i] = (char*)calloc(20,sizeof(char));
+    int8_t filterIDs[MAXDETECT], surveyIDs[MAXDETECT];
+
+
+    // write StackMeta
+    ippToPspsConfig_writeTable(this->config, fitsIn, this->fitsOut, 1, "StackMeta", true);
+    fits_write_col(this->fitsOut, TLONG, STACKMETA_SKYCELLID, 1, 1, 1, &this->id, &status);
+
+    int8_t filterID = -1;
+    if (!ippToPspsConfig_getFilterId(this->config, filterType, &filterID)) {return PS_EXIT_DATA_ERROR;}
+    fits_write_col(this->fitsOut, TBYTE, STACKMETA_FILTERID, 1, 1, 1, &filterID, &status);
+
+    fits_write_col(this->fitsOut, TBYTE, STACKMETA_SURVEYID, 1, 1, 1, &this->surveyID, &status);
+
+
+    // psf detections
+    char extensionName[15];
+    sprintf(extensionName, "SkyChip.psf");
+    if (fits_movnam_hdu(fitsIn, BINARY_TBL, extensionName, 0, &status)) {
+        psError(PS_ERR_IO, false, "Can't move to extension: %s\n", extensionName);
+        
     }
+    else {
+   
+        // some stuff is the same for all detections so we can populate here
+        for (long s = 0; s<MAXDETECT; s++) {
+
+            filterIDs[s] = filterID;
+            surveyIDs[s] = this->surveyID;
+            strcpy(assocDate[s], this->todaysDate);
+        }
+
+
+        long nDet = 0;
+        if (fits_get_num_rows(fitsIn, &nDet, &status)) {
+            fits_report_error(stderr, status);
+        }
+
+        int instMagNum;
+        status=0;fits_get_colnum(fitsIn, CASESEN, "PSF_INST_MAG", &instMagNum, &status);
+        if (status) psError(PS_ERR_IO, false, "Unable to read col num for PSF_INST_MAG");
+        fits_read_col(fitsIn, TFLOAT, instMagNum, 1, 1, nDet, &floatnull, instMag, &anynull, &status);
+
+
+        printf("Looping through %ld psf detections\n", nDet);
+        float mag;
+        long unmatched = 0, totalDetections = 0, numOfDuplicates = 0, numInvalidFlux = 0, numDetectionsOut = 0;
+
+        for (long s = 0; s<nDet; s++) {
+
+            // TODO implement this match in DVO
+            if (1) {
+
+                mag = instMag[s];
+                if (!isfinite(mag) || mag < -998.0) {
+
+                    removeList[numOfDuplicates+numInvalidFlux] = s+1;
+                    numInvalidFlux++;
+                }
+
+                totalDetections++;
+            }
+            else {
+
+                unmatched++;
+                continue;
+            }
+        }
+
+        numDetectionsOut = totalDetections - numInvalidFlux;
+
+        if (numDetectionsOut > 0) {
+
+            ippToPspsConfig_writeTable(this->config, fitsIn, this->fitsOut, nDet, "StackDetection", false);
+            fits_write_col(this->fitsOut, TLONG, STACKDETECTION_SKYCELLID, 1, 1, 1, &this->id, &status);
+            fits_write_col(this->fitsOut, TBYTE, STACKDETECTION_FILTERID, 1, 1, nDet, filterIDs, &status);
+            fits_write_col(this->fitsOut, TBYTE, STACKDETECTION_SURVEYID, 1, 1, nDet, surveyIDs, &status);
+            fits_write_col(this->fitsOut, TSTRING, STACKDETECTION_ASSOCDATE, 1, 1, nDet, assocDate, &status);
+
+            if (numInvalidFlux) fits_delete_rowlist(this->fitsOut, removeList, numInvalidFlux, &status);
+
+        }
+        psLogMsg("ippToPsps", PS_LOG_INFO,
+                "+---------------+---------+----------+------------------+---------------+--------------+\n"
+                "|   Extension   | Rows in | Rows out | Missing from DVO | Duplicate IDs | Invalid Flux |\n"
+                "|  %12s |  %5ld  |   %5ld  |      %5ld       |    %5ld      |    %5ld     |\n",
+                extensionName, nDet, numDetectionsOut, unmatched, numOfDuplicates, numInvalidFlux);
+
+    }
+
+
+
+    // extended source
+    sprintf(extensionName, "SkyChip.xsrc");
+    if (fits_movnam_hdu(fitsIn, BINARY_TBL, extensionName, 0, &status)) {
+        psError(PS_ERR_IO, false, "Can't move to extension: %s\n", extensionName);
+
+    }
+    else {
+
+        long nDet = 0;
+        if (fits_get_num_rows(fitsIn, &nDet, &status)) {
+            fits_report_error(stderr, status);
+        }
+
+        printf("Looping through %ld extended source detections\n", nDet);
+        for (long s = 0; s<nDet; s++) {
+
+
+
+        }
+
+        ippToPspsConfig_writeTable(this->config, fitsIn, this->fitsOut, nDet, "StackApFlx", false);
+        fits_write_col(this->fitsOut, TBYTE, STACKAPFLX_FILTERID, 1, 1, nDet, filterIDs, &status);
+        fits_write_col(this->fitsOut, TBYTE, STACKAPFLX_SURVEYID, 1, 1, nDet, surveyIDs, &status);
+    }
+
+    status=0;
+    if (fits_close_file(fitsIn, &status)) fits_report_error(stderr, status);
+
 
     return PS_EXIT_SUCCESS;