Changeset 41610

Timestamp:

May 25, 2021, 2:07:05 PM (5 years ago)

Author:

eugene

Message:

added code to output a basic dump of the grid corrections

Location:

branches/eam_branches/relphot.20210521/src

Files:

• GridIO.c (added)
• GridOps.c (modified) (2 diffs)
• relphot_images.c (modified) (1 diff)

branches/eam_branches/relphot.20210521/src/GridOps.c

@@ -16 +16 @@
  */
 
-typedef struct {
-  unsigned short photcode;
-  float **Mgrid; // grid of average corrections
-  float **dMgrid; // grid of stdev of corrections
-    int **nMgrid; // grid of number of stars for corrections
-  int Nx; // bin = int(Xchip * (Nx / NxChip))
-  int Ny; // bin = int(Ychip * (Ny / NyChip))
-  float dX; // bin = int(Xchip * dX), dX = Nx / NxChip
-  float dY; // bin = int(Ychip * dY), dY = Ny / NyChip
-  // NxChip, NyChip = 4900,4900 for now
-} GridCorrectionType;
-
 static GridCorrectionType **GridCorr = NULL;
 static int                 NGridCorr = 0;
@@ -95 +83 @@
 
   return;
+}
+
+GridCorrectionType *getGridCorr (int *Nlast) {
+
+  if (GridCorr == NULL) return NULL;
+
+  if (*Nlast >= NGridCorr) return NULL;
+
+  if (*Nlast < 0) *Nlast = -1;
+
+  GridCorrectionType *result = NULL;
+  for (int i = *Nlast + 1; i < NGridCorr; i++) {
+    if (GridCorr[i] == NULL) continue;
+    result = GridCorr[i];
+    *Nlast = i;
+  }
+  return result;
 }
 

branches/eam_branches/relphot.20210521/src/relphot_images.c

@@ -143 +143 @@
     freeMosaicBins (Ncatalog, TRUE);
     freeTGroupBins (Ncatalog);
+
+    GridCorrectionSave ("test.grid.fits");
     // end of if (NLOOP > 0) block : this loop determines the offsets per chip 
   } else {