Changeset 37477 for branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/include/fakeastro.h

Timestamp:

Oct 9, 2014, 3:44:07 PM (12 years ago)

Author:

eugene

Message:

adding code to fit the images to the reference positions

File:

• branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/include/fakeastro.h (modified) (6 diffs)

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/include/fakeastro.h

@@ -9 +9 @@
 
 typedef struct {
+  int   NfakeImage;
+  int   NFAKEIMAGE;
+  Image *fakeImage;
+
+  int   NtrueImage;
+  int   NTRUEIMAGE;
+  Image *trueImage;
+} ImageInfo;
+
+typedef struct {
   double R, D;
   StarPar starpar;
@@ -16 +26 @@
 
 typedef struct {
+  double Rref; // "reference" coordinate of the fake stars
+  double Dref;
   Average  average;
   Measure  measure;
@@ -21 +33 @@
   int found;
 } Stars;
+
+// structure to hold coordinate fitting terms
+typedef struct {
+    int Npts;
+    int Nterms;
+    int Norder;
+    int Nsums;
+    int Nelems;
+    double **sum;
+    double **xsum;
+    double **ysum;
+    double **xfit;
+    double **yfit;
+} CoordFit;
 
 /* used in find_matches, find_matches_refstars */
@@ -33 +59 @@
 # define DVO_MAX_PATH 1024
 
-char   ImageCat[DVO_MAX_PATH];
 char   GSCFILE[DVO_MAX_PATH];
 char   CATDIR[DVO_MAX_PATH];
@@ -116 +141 @@
 Image *make_fake_images (Image *image, int *nfakeImage);
 
-Image *fakeastro_images_region (Image *image, int *nimage, int *NIMAGE, Image *refImage, int NrefImage, SkyTable *skyTableInput, SkyTable *skyTableOutput, SkyRegion *region);
+int fakeastro_images_region (ImageInfo *imageInfo, Image *refImage, int NrefImage, SkyTable *skyTableInput, SkyTable *skyTableOutput, SkyRegion *innerRegion);
 
 SkyRegion *get_image_patch (Image *image);
@@ -135 +160 @@
 float airmass (float secz_image, double ra, double dec, double st, double latitude);
 float azimuth (double ha, double dec, double latitude);
+
+int fit_fake_stars (Stars *stars, int Nstars, Image *image);
+
+double **array_init (int Nx, int Ny);
+void array_free (double **array, int Nx);
+CoordFit *fit_init (int order);
+void fit_free (CoordFit *fit);
+void fit_add (CoordFit *fit, double x1, double y1, double x2, double y2, double wt);
+int fit_eval (CoordFit *fit);
+void fit_apply (CoordFit *fit, double *x2, double *y2, double x1, double y1);
+double **poly2d_dx (double **poly, int Nx, int Ny);
+double **poly2d_dy (double **poly, int Nx, int Ny);
+double **poly2d_copy (double **poly, int Nx, int Ny);
+double poly2d_eval (double **poly, int Nx, int Ny, double x, double y);
+int fit_apply_coords (CoordFit *fit, Coords *coords, int keepRef);
+
+int save_astrom_table ();
+AstromOffsetTable *get_astrom_table ();