IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Feb 6, 2013, 3:26:27 PM (13 years ago)
Author:
eugene
Message:

support to use the Xfix,Yfix values instead of Xraw,Yraw for calculating astrometry

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/relastro/src/BrightCatalog.c

    r34405 r35105  
    7676    GET_COLUMN(Xccd,      "X_CCD",      float);
    7777    GET_COLUMN(Yccd,      "Y_CCD",      float);
     78    GET_COLUMN(Xfix,      "X_FIX",      float);
     79    GET_COLUMN(Yfix,      "Y_FIX",      float);
    7880    GET_COLUMN(dt,        "EXPTIME",    float);
    7981    GET_COLUMN(t,         "TIME",       int);
    … …  
    98100      measure[i].Xccd      = Xccd[i];
    99101      measure[i].Yccd      = Yccd[i];
     102      measure[i].Xfix      = Xfix[i];
     103      measure[i].Yfix      = Yfix[i];
    100104      measure[i].dt        = dt[i];
    101105      measure[i].t         = t[i];
    … …  
    117121    free (Xccd    );
    118122    free (Yccd    );
     123    free (Xfix    );
     124    free (Yfix    );
    119125    free (dt      );
    120126    free (t       );
    … …  
    340346    gfits_define_bintable_column (&theader, "E", "MAG_ERR",  "magnitude (err)",             NULL,    1.0, 0.0);
    341347    gfits_define_bintable_column (&theader, "E", "AIRMASS",  "airmass",                     NULL,    1.0, 0.0);
    342     gfits_define_bintable_column (&theader, "E", "X_CCD",    "ccd x coord",                "pix",    1.0, 0.0);
    343     gfits_define_bintable_column (&theader, "E", "Y_CCD",    "ccd y coord",                "pix",    1.0, 0.0);
     348    gfits_define_bintable_column (&theader, "E", "X_CCD",    "ccd x raw coord",            "pix",    1.0, 0.0);
     349    gfits_define_bintable_column (&theader, "E", "Y_CCD",    "ccd y raw coord",            "pix",    1.0, 0.0);
     350    gfits_define_bintable_column (&theader, "E", "X_FIX",    "ccd x fixed coord",          "pix",    1.0, 0.0);
     351    gfits_define_bintable_column (&theader, "E", "Y_FIX",    "ccd y fiex coord",           "pix",    1.0, 0.0);
    344352    gfits_define_bintable_column (&theader, "E", "EXPTIME",  "-2.5 * log (exposure time)", "sec",    1.0, 0.0);
    345353    gfits_define_bintable_column (&theader, "J", "TIME",     "time of exp",                "sec",    1.0, 1.0*0x8000);
    … …  
    363371    float *Xccd      ; ALLOCATE (Xccd     ,  float, catalog->Nmeasure);
    364372    float *Yccd      ; ALLOCATE (Yccd     ,  float, catalog->Nmeasure);
     373    float *Xfix      ; ALLOCATE (Xfix     ,  float, catalog->Nmeasure);
     374    float *Yfix      ; ALLOCATE (Yfix     ,  float, catalog->Nmeasure);
    365375    float *dt        ; ALLOCATE (dt       ,  float, catalog->Nmeasure);
    366376    int   *t         ; ALLOCATE (t        ,  int  , catalog->Nmeasure);
    … …  
    383393      Xccd[i]     = measure[i].Xccd     ;
    384394      Yccd[i]     = measure[i].Yccd     ;
     395      Xfix[i]     = measure[i].Xfix     ;
     396      Yfix[i]     = measure[i].Yfix     ;
    385397      dt[i]       = measure[i].dt       ;
    386398      t[i]        = measure[i].t        ;
    … …  
    402414    gfits_set_bintable_column (&theader, &ftable, "X_CCD",      Xccd,      catalog->Nmeasure);
    403415    gfits_set_bintable_column (&theader, &ftable, "Y_CCD",      Yccd,      catalog->Nmeasure);
     416    gfits_set_bintable_column (&theader, &ftable, "X_FIX",      Xfix,      catalog->Nmeasure);
     417    gfits_set_bintable_column (&theader, &ftable, "Y_FIX",      Yfix,      catalog->Nmeasure);
    404418    gfits_set_bintable_column (&theader, &ftable, "EXPTIME",    dt,        catalog->Nmeasure);
    405419    gfits_set_bintable_column (&theader, &ftable, "TIME",       t,         catalog->Nmeasure);
    … …  
    419433    free (Xccd    );
    420434    free (Yccd    );
     435    free (Xfix    );
     436    free (Yfix    );
    421437    free (dt      );
    422438    free (t       );
    … …  
    453469    gfits_define_bintable_column (&theader, "J", "MEAN_EPOCH",  "mean epoch (PM-PAR ref)",                                "", 1.0, 1.0*0x8000);
    454470    gfits_define_bintable_column (&theader, "J", "TIME_RANGE",  "mean epoch (PM-PAR ref)",                                "", 1.0, 1.0*0x8000);
    455     gfits_define_bintable_column (&theader, "E", "SIGMA_POS",   "position scatter",                                       "", 1.0, 0.0);   
     471    gfits_define_bintable_column (&theader, "E", "STARGAL_SEP", "star/galaxy separator",                                  "", 1.0, 0.0);   
    456472    gfits_define_bintable_column (&theader, "I", "NUMBER_POS",  "number of detections used for astrometry",               "", 1.0, 1.0*0x80);
    457473    gfits_define_bintable_column (&theader, "I", "NMEASURE",    "number of psf measurements",                             "", 1.0, 1.0*0x80);
Note: See TracChangeset for help on using the changeset viewer.