IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 42018


Ignore:
Timestamp:
Feb 7, 2022, 10:56:09 AM (4 years ago)
Author:
eugene
Message:

fixing UKIDSS transformations; adding measure->extID accessor

Location:
branches/eam_branches/ipp-20211108/Ohana/src
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • branches/eam_branches/ipp-20211108/Ohana/src/addstar/include/ukirt_uhs.h

    r41341 r42018  
    6464char *iparse_csv_rpt (int *X,    int Nwant, int Nlast, char *line, int *status);
    6565char *jparse_csv_rpt (uint64_t *X, int Nwant, int Nlast, char *line, int *status);
     66
     67float psfQFfromXClass (int xClass);

  • branches/eam_branches/ipp-20211108/Ohana/src/addstar/src/loadukirt_uhs_readstars.c

    r41341 r42018  
    4545}
    4646
     47float psfQFfromXClass (int xClass) {
     48
     49  float value = NAN;
     50  switch (xClass) {
     51    case -1: value = 0.99; break; // star
     52    case -2: value = 0.90; break; // probable star
     53    case -3: value = 0.40; break; // probable galaxy
     54    case +1: value = 0.20; break; // galaxy
     55    case  0: value = 0.10; break; // noise
     56    default: break;
     57  }
     58
     59  return value;
     60}
     61

  • branches/eam_branches/ipp-20211108/Ohana/src/addstar/src/loadukirt_uhs_readstars_ugcs.c

    r41341 r42018  
    605605    stars[Nstars].measure[0].Mkron      = yAperMag4[i];
    606606    stars[Nstars].measure[0].dMkron     = yAperMag4Err[i];
    607     stars[Nstars].measure[0].posangle   = yPA[i]; // XXX units
    608     stars[Nstars].measure[0].psfQF      = yClass[i];
     607    stars[Nstars].measure[0].posangle   = ToShortDegrees(yPA[i]); // XXX units
     608    stars[Nstars].measure[0].psfQF      = psfQFfromXClass(yClass[i]);
    609609    stars[Nstars].measure[0].photFlags  = yppErrBits[i];
    610610    stars[Nstars].measure[0].photFlags2 = photFlags2;
    … …  
    625625    stars[Nstars].measure[1].Mkron      = zAperMag4[i];
    626626    stars[Nstars].measure[1].dMkron     = zAperMag4Err[i];
    627     stars[Nstars].measure[1].posangle   = zPA[i]; // XXX units
    628     stars[Nstars].measure[1].psfQF      = zClass[i];
     627    stars[Nstars].measure[1].posangle   = ToShortDegrees(zPA[i]); // XXX units
     628    stars[Nstars].measure[1].psfQF      = psfQFfromXClass(zClass[i]);
    629629    stars[Nstars].measure[1].photFlags  = zppErrBits[i];
    630630    stars[Nstars].measure[1].photFlags2 = photFlags2;
    … …  
    645645    stars[Nstars].measure[2].Mkron      = jAperMag4[i];
    646646    stars[Nstars].measure[2].dMkron     = jAperMag4Err[i];
    647     stars[Nstars].measure[2].posangle   = jPA[i]; // XXX units
    648     stars[Nstars].measure[2].psfQF      = jClass[i];
     647    stars[Nstars].measure[2].posangle   = ToShortDegrees(jPA[i]); // XXX units
     648    stars[Nstars].measure[2].psfQF      = psfQFfromXClass(jClass[i]);
    649649    stars[Nstars].measure[2].photFlags  = jppErrBits[i];
    650650    stars[Nstars].measure[2].photFlags2 = photFlags2;
    … …  
    665665    stars[Nstars].measure[3].Mkron      = hAperMag4[i];
    666666    stars[Nstars].measure[3].dMkron     = hAperMag4Err[i];
    667     stars[Nstars].measure[3].posangle   = hPA[i]; // XXX units
    668     stars[Nstars].measure[3].psfQF      = hClass[i];
     667    stars[Nstars].measure[3].posangle   = ToShortDegrees(hPA[i]); // XXX units
     668    stars[Nstars].measure[3].psfQF      = psfQFfromXClass(hClass[i]);
    669669    stars[Nstars].measure[3].photFlags  = hppErrBits[i];
    670670    stars[Nstars].measure[3].photFlags2 = photFlags2;
    … …  
    685685    stars[Nstars].measure[4].Mkron      = k1AperMag4[i];
    686686    stars[Nstars].measure[4].dMkron     = k1AperMag4Err[i];
    687     stars[Nstars].measure[4].posangle   = k1PA[i]; // XXX units
    688     stars[Nstars].measure[4].psfQF      = k1Class[i];
     687    stars[Nstars].measure[4].posangle   = ToShortDegrees(k1PA[i]); // XXX units
     688    stars[Nstars].measure[4].psfQF      = psfQFfromXClass(k1Class[i]);
    689689    stars[Nstars].measure[4].photFlags  = k1ppErrBits[i];
    690690    stars[Nstars].measure[4].photFlags2 = photFlags2;
    … …  
    705705    stars[Nstars].measure[5].Mkron      = k2AperMag4[i];
    706706    stars[Nstars].measure[5].dMkron     = k2AperMag4Err[i];
    707     stars[Nstars].measure[5].posangle   = k2PA[i]; // XXX units
    708     stars[Nstars].measure[5].psfQF      = k2Class[i];
     707    stars[Nstars].measure[5].posangle   = ToShortDegrees(k2PA[i]); // XXX units
     708    stars[Nstars].measure[5].psfQF      = psfQFfromXClass(k2Class[i]);
    709709    stars[Nstars].measure[5].photFlags  = k2ppErrBits[i];
    710710    stars[Nstars].measure[5].photFlags2 = photFlags2;

  • branches/eam_branches/ipp-20211108/Ohana/src/addstar/src/loadukirt_uhs_readstars_ugps.c

    r41341 r42018  
    516516    stars[Nstars].measure[0].Mkron      = jAperMag4[i];
    517517    stars[Nstars].measure[0].dMkron     = jAperMag4Err[i];
    518     stars[Nstars].measure[0].posangle   = jPA[i]; // XXX units
    519     stars[Nstars].measure[0].psfQF      = jClass[i];
     518    stars[Nstars].measure[0].posangle   = ToShortDegrees(jPA[i]); // XXX units
     519    stars[Nstars].measure[0].psfQF      = psfQFfromXClass(jClass[i]);
    520520    stars[Nstars].measure[0].photFlags  = jppErrBits[i];
    521521    stars[Nstars].measure[0].photFlags2 = photFlags2;
    … …  
    536536    stars[Nstars].measure[1].Mkron      = hAperMag4[i];
    537537    stars[Nstars].measure[1].dMkron     = hAperMag4Err[i];
    538     stars[Nstars].measure[1].posangle   = hPA[i]; // XXX units
    539     stars[Nstars].measure[1].psfQF      = hClass[i];
     538    stars[Nstars].measure[1].posangle   = ToShortDegrees(hPA[i]); // XXX units
     539    stars[Nstars].measure[1].psfQF      = psfQFfromXClass(hClass[i]);
    540540    stars[Nstars].measure[1].photFlags  = hppErrBits[i];
    541541    stars[Nstars].measure[1].photFlags2 = photFlags2;
    … …  
    556556    stars[Nstars].measure[2].Mkron      = k1AperMag4[i];
    557557    stars[Nstars].measure[2].dMkron     = k1AperMag4Err[i];
    558     stars[Nstars].measure[2].posangle   = k1PA[i]; // XXX units
    559     stars[Nstars].measure[2].psfQF      = k1Class[i];
     558    stars[Nstars].measure[2].posangle   = ToShortDegrees(k1PA[i]); // XXX units
     559    stars[Nstars].measure[2].psfQF      = psfQFfromXClass(k1Class[i]);
    560560    stars[Nstars].measure[2].photFlags  = k1ppErrBits[i];
    561561    stars[Nstars].measure[2].photFlags2 = photFlags2;
    … …  
    576576    stars[Nstars].measure[3].Mkron      = k2AperMag4[i];
    577577    stars[Nstars].measure[3].dMkron     = k2AperMag4Err[i];
    578     stars[Nstars].measure[3].posangle   = k2PA[i]; // XXX units
    579     stars[Nstars].measure[3].psfQF      = k2Class[i];
     578    stars[Nstars].measure[3].posangle   = ToShortDegrees(k2PA[i]); // XXX units
     579    stars[Nstars].measure[3].psfQF      = psfQFfromXClass(k2Class[i]);
    580580    stars[Nstars].measure[3].photFlags  = k2ppErrBits[i];
    581581    stars[Nstars].measure[3].photFlags2 = photFlags2;
    … …  
    596596    stars[Nstars].measure[4].Mkron      = h2AperMag4[i];
    597597    stars[Nstars].measure[4].dMkron     = h2AperMag4Err[i];
    598     stars[Nstars].measure[4].posangle   = h2PA[i]; // XXX units
    599     stars[Nstars].measure[4].psfQF      = h2Class[i];
     598    stars[Nstars].measure[4].posangle   = ToShortDegrees(h2PA[i]); // XXX units
     599    stars[Nstars].measure[4].psfQF      = psfQFfromXClass(h2Class[i]);
    600600    stars[Nstars].measure[4].photFlags  = h2ppErrBits[i];
    601601    stars[Nstars].measure[4].photFlags2 = photFlags2;

  • branches/eam_branches/ipp-20211108/Ohana/src/addstar/src/loadukirt_uhs_readstars_uhs.c

    r41341 r42018  
    297297    stars[Nstars].measure[0].Mkron      = jAperMag4[i];
    298298    stars[Nstars].measure[0].dMkron     = jAperMag4Err[i];
    299     stars[Nstars].measure[0].posangle   = jPA[i]; // XXX units
    300     stars[Nstars].measure[0].psfQF      = jClass[i];
     299    stars[Nstars].measure[0].posangle   = ToShortDegrees(jPA[i]); // jPA is in degrees, posangle is short
     300    stars[Nstars].measure[0].psfQF      = psfQFfromXClass(jClass[i]);
    301301    stars[Nstars].measure[0].photFlags  = jppErrBits[i];
    302302    stars[Nstars].measure[0].photFlags2 = photFlags2;

  • branches/eam_branches/ipp-20211108/Ohana/src/addstar/src/loadukirt_uhs_readstars_ulas.c

    r41341 r42018  
    543543    stars[Nstars].measure[0].Mkron      = yAperMag4[i];
    544544    stars[Nstars].measure[0].dMkron     = yAperMag4Err[i];
    545     stars[Nstars].measure[0].posangle   = yPA[i]; // XXX units
    546     stars[Nstars].measure[0].psfQF      = yClass[i];
     545    stars[Nstars].measure[0].posangle   = ToShortDegrees(yPA[i]); // XXX units
     546    stars[Nstars].measure[0].psfQF      = psfQFfromXClass(yClass[i]);
    547547    stars[Nstars].measure[0].photFlags  = yppErrBits[i];
    548548    stars[Nstars].measure[0].photFlags2 = photFlags2;
    … …  
    563563    stars[Nstars].measure[1].Mkron      = j1AperMag4[i];
    564564    stars[Nstars].measure[1].dMkron     = j1AperMag4Err[i];
    565     stars[Nstars].measure[1].posangle   = j1PA[i]; // XXX units
    566     stars[Nstars].measure[1].psfQF      = j1Class[i];
     565    stars[Nstars].measure[1].posangle   = ToShortDegrees(j1PA[i]); // XXX units
     566    stars[Nstars].measure[1].psfQF      = psfQFfromXClass(j1Class[i]);
    567567    stars[Nstars].measure[1].photFlags  = j1ppErrBits[i];
    568568    stars[Nstars].measure[1].photFlags2 = photFlags2;
    … …  
    583583    stars[Nstars].measure[2].Mkron      = j2AperMag4[i];
    584584    stars[Nstars].measure[2].dMkron     = j2AperMag4Err[i];
    585     stars[Nstars].measure[2].posangle   = j2PA[i]; // XXX units
    586     stars[Nstars].measure[2].psfQF      = j2Class[i];
     585    stars[Nstars].measure[2].posangle   = ToShortDegrees(j2PA[i]); // XXX units
     586    stars[Nstars].measure[2].psfQF      = psfQFfromXClass(j2Class[i]);
    587587    stars[Nstars].measure[2].photFlags  = j2ppErrBits[i];
    588588    stars[Nstars].measure[2].photFlags2 = photFlags2;
    … …  
    603603    stars[Nstars].measure[3].Mkron      = hAperMag4[i];
    604604    stars[Nstars].measure[3].dMkron     = hAperMag4Err[i];
    605     stars[Nstars].measure[3].posangle   = hPA[i]; // XXX units
    606     stars[Nstars].measure[3].psfQF      = hClass[i];
     605    stars[Nstars].measure[3].posangle   = ToShortDegrees(hPA[i]); // XXX units
     606    stars[Nstars].measure[3].psfQF      = psfQFfromXClass(hClass[i]);
    607607    stars[Nstars].measure[3].photFlags  = hppErrBits[i];
    608608    stars[Nstars].measure[3].photFlags2 = photFlags2;
    … …  
    623623    stars[Nstars].measure[4].Mkron      = kAperMag4[i];
    624624    stars[Nstars].measure[4].dMkron     = kAperMag4Err[i];
    625     stars[Nstars].measure[4].posangle   = kPA[i]; // XXX units
    626     stars[Nstars].measure[4].psfQF      = kClass[i];
     625    stars[Nstars].measure[4].posangle   = ToShortDegrees(kPA[i]); // XXX units
     626    stars[Nstars].measure[4].psfQF      = psfQFfromXClass(kClass[i]);
    627627    stars[Nstars].measure[4].photFlags  = kppErrBits[i];
    628628    stars[Nstars].measure[4].photFlags2 = photFlags2;

  • branches/eam_branches/ipp-20211108/Ohana/src/libdvo/include/dvodb.h

    r41341 r42018  
    240240      MEAS_OBJ_ID,
    241241      MEAS_CAT_ID,
     242      MEAS_EXT_ID,
    242243      MEAS_IMAGE_ID,
    243244      MEAS_PSF_QF,
    … …  
    248249      MEAS_CR_NSIGMA,
    249250      MEAS_EXT_NSIGMA,
    250       MEAS_EXTERN_ID,
     251      MEAS_IMAGE_EXTERN_ID, // return image.externID
    251252      MEAS_EXPNAME_AS_INT,
    252253      MEAS_MCAL_OFFSET_PSF, // make this a dvoMagOption?

  • branches/eam_branches/ipp-20211108/Ohana/src/libdvo/src/dbExtractMeasures.c

    r40463 r42018  
    696696      value.Int = measure[0].imageID;
    697697      break;
     698    case MEAS_EXT_ID: /* OK */
     699      value.Int = measure[0].extID;
     700      break;
    698701    case MEAS_PSF_QF: /* OK */
    699702      value.Flt = measure[0].psfQF;
    … …  
    714717      value.Flt = measure[0].extNsigma;
    715718      break;
    716     case MEAS_EXTERN_ID: /* OK */
     719    case MEAS_IMAGE_EXTERN_ID: /* OK */
    717720      if (REMOTE_CLIENT) {
    718721        ImageMetadata *image = MatchImageMetadataDVO (measure[0].imageID);

  • branches/eam_branches/ipp-20211108/Ohana/src/libdvo/src/dbFields.c

    r40463 r42018  
    603603  if (!strcasecmp (fieldName, "CATID"))          ESCAPE (MEAS_CAT_ID,         OPIHI_INT);
    604604  if (!strcasecmp (fieldName, "IMAGEID"))        ESCAPE (MEAS_IMAGE_ID,       OPIHI_INT);
    605   if (!strcasecmp (fieldName, "EXTERNID"))       ESCAPE (MEAS_EXTERN_ID,      OPIHI_INT);
     605  if (!strcasecmp (fieldName, "IMAGE_EXTERN_ID")) ESCAPE (MEAS_IMAGE_EXTERN_ID, OPIHI_INT);
     606  if (!strcasecmp (fieldName, "IMAGE_EXT_ID"))    ESCAPE (MEAS_IMAGE_EXTERN_ID, OPIHI_INT);
     607  if (!strcasecmp (fieldName, "EXTERNID"))       ESCAPE (MEAS_EXT_ID,         OPIHI_INT);
     608  if (!strcasecmp (fieldName, "EXT_ID"))         ESCAPE (MEAS_EXT_ID,         OPIHI_INT);
    606609  if (!strcasecmp (fieldName, "EXPNAME"))        ESCAPE (MEAS_EXPNAME_AS_INT, OPIHI_INT);
    607610  if (!strcasecmp (fieldName, "PSF_QF"))         ESCAPE (MEAS_PSF_QF,         OPIHI_FLT);
Note: See TracChangeset for help on using the changeset viewer.