Changeset 37457

Timestamp:

Oct 2, 2014, 2:19:23 PM (12 years ago)

Author:

eugene

Message:

adding fake images and stars

Location:

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src

Files:

• make_fake_images.c (modified) (1 diff)
• make_fake_stars.c (added)

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fake_images.c

@@ -1 +1 @@
 # include "fakeastro.h"
 
-Images *make_boresites (Images *images, int Nimages; int *nbores) {
+Image *make_fake_images (Image *image, int *nfakeImage) {
 
+  // given an image (really, the WCS) which defines an exposure, generate the image
+  // parameters for a full set of chips in an exposure
   
+  int NfakeImage = 65;
+  Image *fakeImage = NULL;
 
+  ALLOCATE (fakeImage, Image, NfakeImage);
+
+  // first image is the PHU, copy basic data:
+
+  fakeImage[0] = image[0];
+
+  int N = 1;
+
+  int ix, iy;
+  for (ix = 0; ix < 8; ix++) {
+    for (iy = 0; iy < 8; iy++) {
+      
+      fakeImage[N] = fakeImage[0];
+
+      // coords map relative to image center:
+      fakeImage[N].coords.crval1 = 0.0;
+      fakeImage[N].coords.crval2 = 0.0;
+
+      if (ix < 4) {
+        fakeImage[N].coords.crpix1 = (ix - 3)*4900;
+        fakeImage[N].coords.crpix2 = (iy - 3)*4900;
+        fakeImage[N].coords.pc1_1 = 1.0
+        fakeImage[N].coords.pc2_2 = 1.0
+      } else {
+        fakeImage[N].coords.crpix1 = (3 - ix)*4900;
+        fakeImage[N].coords.crpix2 = (4 - iy)*4900;
+        fakeImage[N].coords.pc1_1 = -1.0
+        fakeImage[N].coords.pc2_2 = -1.0
+      } 
+
+      fakeImage[N].coords.cdelt1 = 1.0;
+      fakeImage[N].coords.cdelt2 = 1.0;
+
+      fakeImage[N].coords.pc1_2 = 0.0;
+      fakeImage[N].coords.pc1_2 = 0.0;
+
+      strcpy (fakeImage[N].coords.ctype, "DEC--WRP");
+      fakeImage[N].coords.Npolyterms = 1;
+
+      fakeImage[N].coords.mosaic = &fakeImage[0].coords;
+
+      fakeImage[N].NX = 4850;
+      fakeImage[N].NY = 4850;
+
+      fakeImage[N].photcode = 10000 + iy*10 + ix;
+      fakeImage[N].exptime = 45.0;
+      N++;
+    }
+  }
+
+  *nfakeImage = NfakeImage;
+  return fakeImage;
 }
+