IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jun 7, 2014, 6:33:11 AM (12 years ago)
Author:
eugene
Message:

merge changes from eam_branches/ipp-20140423: add dcr-color option, ref_color, bad detection clipping to relastro; fix RA,DEC for mosaics in imextract; allow photdbc to restrict images to output region; change PHOTCODE_SKIP_LIST to PHOTCODE_KEEP_LIST in photdbc; add -reset-to-photcode option to relastro; ensure relastro child processes get the complete set of options; allow relastro region hosts to continue after load-object clients are done, but before the data is loaded; invalidate extreme proper motion; limit relastro catalog density on a grid; require more detections for higher order fits; option to delete duplicate images based on exp time/date and photcode; have kapa use XSync instead of XFlush to avoid flashing; handle some cases to avoid crashing kapa; extend kapa connection timeout; allow kapa to handle > 2GB images; ensure dvomerge works even if incoming db has no images; add dvo fuctions coordmosaic & psastro_model

Location:
trunk/Ohana/src/libdvo
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/libdvo/include/dvo.h

    r36630 r36833  
    171171  ID_STAR_USE_PAR   = 0x00200000, // parallax used (not AVE or PM)
    172172  ID_STAR_NO_ASTROM = 0x00400000, // mean astrometry could not be measured
     173  ID_STAR_BAD_PM    = 0x00800000, // mean astrometry could not be measured
    173174  ID_OBJ_EXT        = 0x01000000, // extended in our data (eg, PS)
    174175  ID_OBJ_EXT_ALT    = 0x02000000, // extended in external data (eg, 2MASS)
    … …  
    481482  uint32_t         measureOffset;        // offset to first psf measurement
    482483  uint32_t         missingOffset;        // offset to first missing obs
    483   uint32_t         extendOffset;         // offset to first extended measurement
     484  float            refColor;            // offset to first extended measurement
    484485  uint32_t         flags;                // average object flags (star; ghost; etc)
    485486  uint32_t         photFlagsUpper;       // upper bit of 2 bit summary of per-measure photflags

  • trunk/Ohana/src/libdvo/include/dvodb.h

    r36630 r36833  
    222222      AVE_EXTID_HI,
    223223      AVE_EXTID_LO,
     224      AVE_REF_COLOR,
    224225};
    225226
    … …  
    289290      IMAGE_NFIT_ASTROM,
    290291      IMAGE_NLINK_PHOTOM,
    291       IMAGE_NLINK_ASTROM
     292      IMAGE_NLINK_ASTROM,
     293      IMAGE_REF_COLOR
    292294};
    293295

  • trunk/Ohana/src/libdvo/src/dbExtractAverages.c

    r36630 r36833  
    161161    case AVE_STARGAL:
    162162      value.Flt = average[0].stargal;
     163      break;
     164    case AVE_REF_COLOR:
     165      value.Flt = average[0].refColor;
    163166      break;
    164167

  • trunk/Ohana/src/libdvo/src/dbExtractImages.c

    r35102 r36833  
    6969      if (!haveCelestial) {
    7070        if (!FindMosaicForImage (image, Nimage, N)) return value;
    71         x = 0.5*image[N].NX;
    72         y = 0.5*image[N].NY;
     71        if (!strcmp(&image[N].coords.ctype[4], "-DIS")) {
     72          x = 0.0;
     73          y = 0.0;
     74        } else {
     75          x = 0.5*image[N].NX;
     76          y = 0.5*image[N].NY;
     77        }
    7378        XY_to_RD (&RAo, &DECo, x, y, &image[N].coords);
    7479        haveCelestial = TRUE;
    … …  
    7984      if (!haveCelestial) {
    8085        if (!FindMosaicForImage (image, Nimage, N)) return value;
    81         x = 0.5*image[N].NX;
    82         y = 0.5*image[N].NY;
     86        if (!strcmp(&image[N].coords.ctype[4], "-DIS")) {
     87          x = 0.0;
     88          y = 0.0;
     89        } else {
     90          x = 0.5*image[N].NX;
     91          y = 0.5*image[N].NY;
     92        }
    8393        XY_to_RD (&RAo, &DECo, x, y, &image[N].coords);
    8494        haveCelestial = TRUE;
    … …  
    90100        if (!haveCelestial) {
    91101          if (!FindMosaicForImage (image, Nimage, N)) return value;
    92           x = 0.5*image[N].NX;
    93           y = 0.5*image[N].NY;
     102          if (!strcmp(&image[N].coords.ctype[4], "-DIS")) {
     103            x = 0.0;
     104            y = 0.0;
     105          } else {
     106            x = 0.5*image[N].NX;
     107            y = 0.5*image[N].NY;
     108          }
    94109          XY_to_RD (&RAo, &DECo, x, y, &image[N].coords);
    95110          haveCelestial = TRUE;
    … …  
    104119        if (!haveCelestial) {
    105120          if (!FindMosaicForImage (image, Nimage, N)) return value;
    106           x = 0.5*image[N].NX;
    107           y = 0.5*image[N].NY;
     121          if (!strcmp(&image[N].coords.ctype[4], "-DIS")) {
     122            x = 0.0;
     123            y = 0.0;
     124          } else {
     125            x = 0.5*image[N].NX;
     126            y = 0.5*image[N].NY;
     127          }
    108128          XY_to_RD (&RAo, &DECo, x, y, &image[N].coords);
    109129          haveCelestial = TRUE;
    … …  
    118138        if (!haveCelestial) {
    119139          if (!FindMosaicForImage (image, Nimage, N)) return value;
    120           x = 0.5*image[N].NX;
    121           y = 0.5*image[N].NY;
     140          if (!strcmp(&image[N].coords.ctype[4], "-DIS")) {
     141            x = 0.0;
     142            y = 0.0;
     143          } else {
     144            x = 0.5*image[N].NX;
     145            y = 0.5*image[N].NY;
     146          }
    122147          XY_to_RD (&RAo, &DECo, x, y, &image[N].coords);
    123148          haveCelestial = TRUE;
    … …  
    132157        if (!haveCelestial) {
    133158          if (!FindMosaicForImage (image, Nimage, N)) return value;
    134           x = 0.5*image[N].NX;
    135           y = 0.5*image[N].NY;
     159          if (!strcmp(&image[N].coords.ctype[4], "-DIS")) {
     160            x = 0.0;
     161            y = 0.0;
     162          } else {
     163            x = 0.5*image[N].NX;
     164            y = 0.5*image[N].NY;
     165          }
    136166          XY_to_RD (&RAo, &DECo, x, y, &image[N].coords);
    137167          haveCelestial = TRUE;
    … …  
    344374      value.Int = image[N].nLinkAstrom;
    345375      break;
     376
     377    case IMAGE_REF_COLOR:
     378      value.Flt = image[N].refColor;
     379      break;
    346380  }
    347381  return (value);

  • trunk/Ohana/src/libdvo/src/dbFields.c

    r36630 r36833  
    381381  if (!strcasecmp (fieldName, "EXTID_HI"))    ESCAPE (AVE_EXTID_HI,    MAG_NONE, OPIHI_INT);
    382382  if (!strcasecmp (fieldName, "EXTID_LO"))    ESCAPE (AVE_EXTID_LO,    MAG_NONE, OPIHI_INT);
     383  if (!strcasecmp (fieldName, "REF_COLOR"))   ESCAPE (AVE_REF_COLOR,   MAG_NONE, OPIHI_FLT);
    383384
    384385  // check for code:mode in photcode name
    … …  
    522523  if (!strcasecmp (fieldName, "UCDIST")) ESCAPE (IMAGE_UBERCAL_DIST,MAG_NONE, OPIHI_INT);
    523524
    524   if (!strcasecmp (fieldName, "NFIT_PHOTOM")) ESCAPE (IMAGE_NFIT_PHOTOM,MAG_NONE, OPIHI_INT);
    525   if (!strcasecmp (fieldName, "NFIT_ASTROM")) ESCAPE (IMAGE_NFIT_ASTROM,MAG_NONE, OPIHI_INT);
    526   if (!strcasecmp (fieldName, "NLINK_PHOTOM")) ESCAPE (IMAGE_NLINK_PHOTOM,MAG_NONE, OPIHI_INT);
    527   if (!strcasecmp (fieldName, "NLINK_ASTROM")) ESCAPE (IMAGE_NLINK_ASTROM,MAG_NONE, OPIHI_INT);
     525  if (!strcasecmp (fieldName, "NFIT_PHOTOM"))  ESCAPE (IMAGE_NFIT_PHOTOM,  MAG_NONE, OPIHI_INT);
     526  if (!strcasecmp (fieldName, "NFIT_ASTROM"))  ESCAPE (IMAGE_NFIT_ASTROM,  MAG_NONE, OPIHI_INT);
     527  if (!strcasecmp (fieldName, "NLINK_PHOTOM")) ESCAPE (IMAGE_NLINK_PHOTOM, MAG_NONE, OPIHI_INT);
     528  if (!strcasecmp (fieldName, "NLINK_ASTROM")) ESCAPE (IMAGE_NLINK_ASTROM, MAG_NONE, OPIHI_INT);
     529  if (!strcasecmp (fieldName, "REF_COLOR"))    ESCAPE (IMAGE_REF_COLOR,    MAG_NONE, OPIHI_FLT);
    528530
    529531  // for words that don't parse, try a photcode

  • trunk/Ohana/src/libdvo/src/dvo_catalog.c

    r36630 r36833  
    100100  average->measureOffset   = -1;
    101101  average->missingOffset   = -1;
    102   average->extendOffset    = -1;
     102  average->refColor        = NAN;
    103103
    104104  average->flags           = 0;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_DEV_2.c

    r35162 r36833  
    119119    out[i].measureOffset = in[i].measureOffset;
    120120    out[i].missingOffset = in[i].missingOffset;
    121     out[i].extendOffset  = in[i].extendOffset;
     121    out[i].refColor  = in[i].refColor;
    122122    out[i].objID         = in[i].objID;
    123123    out[i].catID         = in[i].catID;
    … …  
    153153    out[i].measureOffset = in[i].measureOffset;
    154154    out[i].missingOffset = in[i].missingOffset;
    155     out[i].extendOffset  = in[i].extendOffset;
     155    out[i].refColor  = in[i].refColor;
    156156    out[i].objID         = in[i].objID;
    157157    out[i].catID         = in[i].catID;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V1.c

    r35162 r36833  
    140140    out[i].measureOffset = in[i].measureOffset;
    141141    out[i].missingOffset = in[i].missingOffset;
    142     out[i].extendOffset  = in[i].extendOffset;
     142    out[i].refColor  = in[i].refColor;
    143143    out[i].flags         = in[i].flags;   
    144144    out[i].objID         = in[i].objID;
    … …  
    175175    out[i].measureOffset = in[i].measureOffset;
    176176    out[i].missingOffset = in[i].missingOffset;
    177     out[i].extendOffset  = in[i].extendOffset;
     177    out[i].refColor  = in[i].refColor;
    178178    out[i].flags         = in[i].flags;   
    179179    out[i].objID         = in[i].objID;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V2.c

    r35162 r36833  
    145145    out[i].measureOffset = in[i].measureOffset;
    146146    out[i].missingOffset = in[i].missingOffset;
    147     out[i].extendOffset  = in[i].extendOffset;
     147    out[i].refColor  = in[i].refColor;
    148148    out[i].flags         = in[i].flags;   
    149149    out[i].objID         = in[i].objID;
    … …  
    184184    out[i].measureOffset = in[i].measureOffset;
    185185    out[i].missingOffset = in[i].missingOffset;
    186     out[i].extendOffset  = in[i].extendOffset;
     186    out[i].refColor  = in[i].refColor;
    187187    out[i].flags         = in[i].flags;   
    188188    out[i].objID         = in[i].objID;
    … …  
    269269    out[i].DECo             = in[i].DECo;
    270270    out[i].Radius           = in[i].Radius;
     271    out[i].refColor         = in[i].refColor;
    271272
    272273    out[i].detection_limit  = in[i].detection_limit;
    … …  
    331332    out[i].DECo             = in[i].DECo;
    332333    out[i].Radius           = in[i].Radius;
     334    out[i].refColor         = in[i].refColor;
    333335
    334336    out[i].detection_limit  = in[i].detection_limit;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V3.c

    r35162 r36833  
    145145    out[i].measureOffset = in[i].measureOffset;
    146146    out[i].missingOffset = in[i].missingOffset;
    147     out[i].extendOffset  = in[i].extendOffset;
     147    out[i].refColor  = in[i].refColor;
    148148    out[i].flags         = in[i].flags;   
    149149    out[i].photFlagsUpper = in[i].photFlagsUpper;   
    … …  
    186186    out[i].measureOffset = in[i].measureOffset;
    187187    out[i].missingOffset = in[i].missingOffset;
    188     out[i].extendOffset  = in[i].extendOffset;
     188    out[i].refColor  = in[i].refColor;
    189189    out[i].flags         = in[i].flags;   
    190190    out[i].photFlagsUpper = in[i].photFlagsUpper;   
    … …  
    279279    out[i].DECo             = in[i].DECo;
    280280    out[i].Radius           = in[i].Radius;
     281    out[i].refColor         = in[i].refColor;
    281282
    282283    out[i].detection_limit  = in[i].detection_limit;
    … …  
    341342    out[i].DECo             = in[i].DECo;
    342343    out[i].Radius           = in[i].Radius;
     344    out[i].refColor         = in[i].refColor;
    343345
    344346    out[i].detection_limit  = in[i].detection_limit;

  • trunk/Ohana/src/libdvo/src/dvo_convert_PS1_V4.c

    r36630 r36833  
    166166    out[i].measureOffset = in[i].measureOffset;
    167167    out[i].missingOffset = in[i].missingOffset;
    168     out[i].extendOffset  = in[i].extendOffset;
     168    out[i].refColor      = in[i].refColor;
    169169    out[i].flags         = in[i].flags;   
    170170    out[i].photFlagsUpper = in[i].photFlagsUpper;   
    … …  
    210210    out[i].measureOffset = in[i].measureOffset;
    211211    out[i].missingOffset = in[i].missingOffset;
    212     out[i].extendOffset  = in[i].extendOffset;
     212    out[i].refColor      = in[i].refColor;
    213213    out[i].flags         = in[i].flags;   
    214214    out[i].photFlagsUpper = in[i].photFlagsUpper;   
    … …  
    317317    out[i].DECo             = in[i].DECo;
    318318    out[i].Radius           = in[i].Radius;
     319    out[i].refColor         = in[i].refColor;
    319320
    320321    out[i].detection_limit  = in[i].detection_limit;
    … …  
    379380    out[i].DECo             = in[i].DECo;
    380381    out[i].Radius           = in[i].Radius;
     382    out[i].refColor         = in[i].refColor;
    381383
    382384    out[i].detection_limit  = in[i].detection_limit;
    … …  
    680682    out[i].measureOffset = in[i].measureOffset;
    681683    out[i].missingOffset = in[i].missingOffset;
    682     out[i].extendOffset  = in[i].extendOffset;
     684    out[i].refColor      = in[i].refColor;
    683685    out[i].flags         = in[i].flags;   
    684686    out[i].photFlagsUpper = in[i].photFlagsUpper;   

  • trunk/Ohana/src/libdvo/src/skyregion_io.c

    r33649 r36833  
    7878  gfits_free_matrix (&matrix);
    7979  gfits_free_header (&theader);
     80  fclose (f);
    8081
    8182  return (skytable);
Note: See TracChangeset for help on using the changeset viewer.