IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 15208


Ignore:
Timestamp:
Oct 3, 2007, 4:55:24 PM (19 years ago)
Author:
eugene
Message:

updates to the db interations

Location:
trunk/Ohana/src/opihi
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/opihi/dvo/dbExtractMeasures.c

    r15038 r15208  
    136136      value = measure[0].airmass;
    137137      break;
     138    case MEAS_AZ: /* OK */
     139      value = measure[0].az;
     140      break;
    138141    case MEAS_EXPTIME: /* OK */
    139142      value = pow (10.0, measure[0].dt * 0.4);
    … …  
    157160      value = measure[0].theta;
    158161      break;
    159 //    case MEAS_DOPHOT: /* OK */
    160 //      value = measure[0].dophot;
    161 //      break;
    162     case MEAS_FLAGS: /* ? */
     162    case MEAS_DOPHOT: /* OK */
     163      value = measure[0].dophot;
     164      break;
     165    case MEAS_DB_FLAGS: /* ? */
    163166      value = measure[0].dbFlags;
     167      break;
     168    case MEAS_PHOT_FLAGS: /* ? */
     169      value = measure[0].photFlags;
    164170      break;
    165171    case MEAS_XCCD: /* OK */
    … …  
    205211      value = y;
    206212      break;
     213
     214    case MEAS_SKY: /* OK */
     215      value = measure[0].Sky;
     216      break;
     217    case MEAS_dSKY: /* OK */
     218      value = measure[0].dSky;
     219      break;
     220    case MEAS_DET_ID: /* OK */
     221      value = measure[0].detID;
     222      break;
     223    case MEAS_OBJ_ID: /* OK */
     224      value = measure[0].averef;
     225      break;
     226    case MEAS_IMAGE_ID: /* OK */
     227      value = measure[0].imageID;
     228      break;
     229    case MEAS_PSF_QF: /* OK */
     230      value = measure[0].qPSF;
     231      break;
     232    case MEAS_PSF_PROB: /* OK */
     233      value = measure[0].psfProb;
     234      break;
     235    case MEAS_CR_NSIGMA: /* OK */
     236      value = measure[0].crNsigma;
     237      break;
     238    case MEAS_EXT_NSIGMA: /* OK */
     239      value = measure[0].extNsigma;
     240      break;
     241    case MEAS_STARGAL: /* OK */
     242      value = measure[0].stargal;
     243      break;
    207244  }
    208245  return (value);

  • trunk/Ohana/src/opihi/dvo/dbFields.c

    r14590 r15208  
    103103  field->name  = strcreate (fieldName);
    104104
    105   if (!strcasecmp (fieldName, "RA"))       ESCAPE (MEAS_RA,             MAG_NONE);
    106   if (!strcasecmp (fieldName, "DEC"))      ESCAPE (MEAS_DEC,            MAG_NONE);
    107   if (!strcasecmp (fieldName, "RA:AVE"))   ESCAPE (MEAS_RA_AVE,         MAG_NONE);
    108   if (!strcasecmp (fieldName, "DEC:AVE"))  ESCAPE (MEAS_DEC_AVE,        MAG_NONE);
    109   if (!strcasecmp (fieldName, "RA:ERR"))   ESCAPE (MEAS_RA_AVE_ERR,     MAG_NONE);
    110   if (!strcasecmp (fieldName, "DEC:ERR"))  ESCAPE (MEAS_DEC_AVE_ERR,    MAG_NONE);
    111   if (!strcasecmp (fieldName, "uRA"))      ESCAPE (MEAS_U_RA,           MAG_NONE);
    112   if (!strcasecmp (fieldName, "uDEC"))     ESCAPE (MEAS_U_DEC,          MAG_NONE);
    113   if (!strcasecmp (fieldName, "duRA"))     ESCAPE (MEAS_U_RA_ERR,       MAG_NONE);
    114   if (!strcasecmp (fieldName, "duDEC"))    ESCAPE (MEAS_U_DEC_ERR,      MAG_NONE);
    115   if (!strcasecmp (fieldName, "PAR"))      ESCAPE (MEAS_PAR,            MAG_NONE);
    116   if (!strcasecmp (fieldName, "dPAR"))     ESCAPE (MEAS_PAR_ERR,        MAG_NONE);
    117   if (!strcasecmp (fieldName, "dR"))       ESCAPE (MEAS_RA_OFFSET,      MAG_NONE);
    118   if (!strcasecmp (fieldName, "dD"))       ESCAPE (MEAS_DEC_OFFSET,     MAG_NONE);
    119   if (!strcasecmp (fieldName, "dR:ERR"))   ESCAPE (MEAS_RA_OFFSET_ERR,  MAG_NONE);
    120   if (!strcasecmp (fieldName, "dD:ERR"))   ESCAPE (MEAS_DEC_OFFSET_ERR, MAG_NONE);
    121   if (!strcasecmp (fieldName, "xp"))       ESCAPE (MEAS_XP,             MAG_NONE);
    122   if (!strcasecmp (fieldName, "nmeas"))    ESCAPE (MEAS_NMEAS,          MAG_NONE);
    123   if (!strcasecmp (fieldName, "nmiss"))    ESCAPE (MEAS_NMISS,          MAG_NONE);
    124   if (!strcasecmp (fieldName, "objflags")) ESCAPE (MEAS_OBJFLAGS,       MAG_NONE);
    125   if (!strcasecmp (fieldName, "AIRMASS"))  ESCAPE (MEAS_AIRMASS,        MAG_NONE);
    126   if (!strcasecmp (fieldName, "ALT"))      ESCAPE (MEAS_ALT,            MAG_NONE);
    127   if (!strcasecmp (fieldName, "AZ"))       ESCAPE (MEAS_AZ,             MAG_NONE);
    128   if (!strcasecmp (fieldName, "EXPTIME"))  ESCAPE (MEAS_EXPTIME,        MAG_NONE);
    129   if (!strcasecmp (fieldName, "PHOTCODE")) ESCAPE (MEAS_PHOTCODE,       MAG_NONE);
    130   if (!strcasecmp (fieldName, "TIME"))     ESCAPE (MEAS_TIME,           MAG_NONE);
    131   if (!strcasecmp (fieldName, "FWHM"))     ESCAPE (MEAS_FWHM,           MAG_NONE);
    132   if (!strcasecmp (fieldName, "FWHM_MAJ")) ESCAPE (MEAS_FWHM_MAJ,       MAG_NONE);
    133   if (!strcasecmp (fieldName, "FWHM_MIN")) ESCAPE (MEAS_FWHM_MIN,       MAG_NONE);
    134   if (!strcasecmp (fieldName, "THETA"))    ESCAPE (MEAS_THETA,          MAG_NONE);
    135   if (!strcasecmp (fieldName, "FLAGS"))    ESCAPE (MEAS_FLAGS,          MAG_NONE);
    136   if (!strcasecmp (fieldName, "XCCD"))     ESCAPE (MEAS_XCCD,           MAG_NONE);
    137   if (!strcasecmp (fieldName, "YCCD"))     ESCAPE (MEAS_YCCD,           MAG_NONE);
    138   if (!strcasecmp (fieldName, "XMOSAIC"))  ESCAPE (MEAS_XMOSAIC,        MAG_NONE);
    139   if (!strcasecmp (fieldName, "YMOSAIC"))  ESCAPE (MEAS_YMOSAIC,        MAG_NONE);
    140   if (!strcasecmp (fieldName, "XCHIP"))    ESCAPE (MEAS_XCCD,           MAG_NONE);
    141   if (!strcasecmp (fieldName, "YCHIP"))    ESCAPE (MEAS_YCCD,           MAG_NONE);
    142   if (!strcasecmp (fieldName, "XFPA"))     ESCAPE (MEAS_XMOSAIC,        MAG_NONE);
    143   if (!strcasecmp (fieldName, "YFPA"))     ESCAPE (MEAS_YMOSAIC,        MAG_NONE);
    144   if (!strcasecmp (fieldName, "DETID"))    ESCAPE (MEAS_DETID,          MAG_NONE);
    145   if (!strcasecmp (fieldName, "OBJID"))    ESCAPE (MEAS_OBJID,          MAG_NONE);
    146   if (!strcasecmp (fieldName, "IMAGEID"))  ESCAPE (MEAS_IMAGEID,        MAG_NONE);
    147   if (!strcasecmp (fieldName, "stargal"))  ESCAPE (MEAS_IMAGEID,        MAG_NONE);
    148   if (!strcasecmp (fieldName, "psfqf"))    ESCAPE (MEAS_IMAGEID,        MAG_NONE);
    149   if (!strcasecmp (fieldName, "sky"))      ESCAPE (MEAS_IMAGEID,        MAG_NONE);
    150   if (!strcasecmp (fieldName, "sky_err"))  ESCAPE (MEAS_IMAGEID,        MAG_NONE);
     105  if (!strcasecmp (fieldName, "RA"))         ESCAPE (MEAS_RA,           MAG_NONE);
     106  if (!strcasecmp (fieldName, "DEC"))        ESCAPE (MEAS_DEC,          MAG_NONE);
     107  if (!strcasecmp (fieldName, "RA:AVE"))     ESCAPE (MEAS_RA_AVE,       MAG_NONE);
     108  if (!strcasecmp (fieldName, "DEC:AVE"))    ESCAPE (MEAS_DEC_AVE,      MAG_NONE);
     109  if (!strcasecmp (fieldName, "RA:ERR"))     ESCAPE (MEAS_RA_AVE_ERR,   MAG_NONE);
     110  if (!strcasecmp (fieldName, "DEC:ERR"))    ESCAPE (MEAS_DEC_AVE_ERR,  MAG_NONE);
     111  if (!strcasecmp (fieldName, "uRA"))        ESCAPE (MEAS_U_RA,         MAG_NONE);
     112  if (!strcasecmp (fieldName, "uDEC"))       ESCAPE (MEAS_U_DEC,        MAG_NONE);
     113  if (!strcasecmp (fieldName, "duRA"))       ESCAPE (MEAS_U_RA_ERR,     MAG_NONE);
     114  if (!strcasecmp (fieldName, "duDEC"))      ESCAPE (MEAS_U_DEC_ERR,    MAG_NONE);
     115  if (!strcasecmp (fieldName, "PAR"))        ESCAPE (MEAS_PAR,          MAG_NONE);
     116  if (!strcasecmp (fieldName, "dPAR"))       ESCAPE (MEAS_PAR_ERR,      MAG_NONE);
     117  if (!strcasecmp (fieldName, "dR"))         ESCAPE (MEAS_RA_OFFSET,    MAG_NONE);
     118  if (!strcasecmp (fieldName, "dD"))         ESCAPE (MEAS_DEC_OFFSET,   MAG_NONE);
     119  if (!strcasecmp (fieldName, "dR:ERR"))     ESCAPE (MEAS_RA_OFFSET_ERR,  MAG_NONE);
     120  if (!strcasecmp (fieldName, "dD:ERR"))     ESCAPE (MEAS_DEC_OFFSET_ERR, MAG_NONE);
     121  if (!strcasecmp (fieldName, "xp"))         ESCAPE (MEAS_XP,           MAG_NONE);
     122  if (!strcasecmp (fieldName, "nmeas"))      ESCAPE (MEAS_NMEAS,        MAG_NONE);
     123  if (!strcasecmp (fieldName, "nmiss"))      ESCAPE (MEAS_NMISS,        MAG_NONE);
     124  if (!strcasecmp (fieldName, "objflags"))   ESCAPE (MEAS_OBJFLAGS,     MAG_NONE);
     125  if (!strcasecmp (fieldName, "AIRMASS"))    ESCAPE (MEAS_AIRMASS,      MAG_NONE);
     126  if (!strcasecmp (fieldName, "ALT"))        ESCAPE (MEAS_ALT,          MAG_NONE);
     127  if (!strcasecmp (fieldName, "AZ"))         ESCAPE (MEAS_AZ,           MAG_NONE);
     128  if (!strcasecmp (fieldName, "EXPTIME"))    ESCAPE (MEAS_EXPTIME,      MAG_NONE);
     129  if (!strcasecmp (fieldName, "PHOTCODE"))   ESCAPE (MEAS_PHOTCODE,     MAG_NONE);
     130  if (!strcasecmp (fieldName, "TIME"))       ESCAPE (MEAS_TIME,         MAG_NONE);
     131  if (!strcasecmp (fieldName, "FWHM"))       ESCAPE (MEAS_FWHM,         MAG_NONE);
     132  if (!strcasecmp (fieldName, "FWHM_MAJ"))   ESCAPE (MEAS_FWHM_MAJ,     MAG_NONE);
     133  if (!strcasecmp (fieldName, "FWHM_MIN"))   ESCAPE (MEAS_FWHM_MIN,     MAG_NONE);
     134  if (!strcasecmp (fieldName, "THETA"))      ESCAPE (MEAS_THETA,        MAG_NONE);
     135  if (!strcasecmp (fieldName, "DOPHOT"))     ESCAPE (MEAS_DOPHOT,       MAG_NONE);
     136  if (!strcasecmp (fieldName, "DB_FLAGS"))   ESCAPE (MEAS_DB_FLAGS,     MAG_NONE);
     137  if (!strcasecmp (fieldName, "PHOT_FLAGS")) ESCAPE (MEAS_PHOT_FLAGS,   MAG_NONE);
     138  if (!strcasecmp (fieldName, "XCCD"))       ESCAPE (MEAS_XCCD,         MAG_NONE);
     139  if (!strcasecmp (fieldName, "YCCD"))       ESCAPE (MEAS_YCCD,         MAG_NONE);
     140  if (!strcasecmp (fieldName, "XMOSAIC"))    ESCAPE (MEAS_XMOSAIC,      MAG_NONE);
     141  if (!strcasecmp (fieldName, "YMOSAIC"))    ESCAPE (MEAS_YMOSAIC,      MAG_NONE);
     142  if (!strcasecmp (fieldName, "XCHIP"))      ESCAPE (MEAS_XCCD,         MAG_NONE);
     143  if (!strcasecmp (fieldName, "YCHIP"))      ESCAPE (MEAS_YCCD,         MAG_NONE);
     144  if (!strcasecmp (fieldName, "XFPA"))       ESCAPE (MEAS_XMOSAIC,      MAG_NONE);
     145  if (!strcasecmp (fieldName, "YFPA"))       ESCAPE (MEAS_YMOSAIC,      MAG_NONE);
     146  if (!strcasecmp (fieldName, "DETID"))      ESCAPE (MEAS_DET_ID,       MAG_NONE);
     147  if (!strcasecmp (fieldName, "OBJID"))      ESCAPE (MEAS_OBJ_ID,       MAG_NONE);
     148  if (!strcasecmp (fieldName, "IMAGEID"))    ESCAPE (MEAS_IMAGE_ID,     MAG_NONE);
     149  if (!strcasecmp (fieldName, "PSF_QF"))     ESCAPE (MEAS_PSF_QF,       MAG_NONE);
     150  if (!strcasecmp (fieldName, "PSF_PROB"))   ESCAPE (MEAS_PSF_PROB,     MAG_NONE);
     151  if (!strcasecmp (fieldName, "CR_NSIGMA"))  ESCAPE (MEAS_CR_NSIGMA,    MAG_NONE);
     152  if (!strcasecmp (fieldName, "EXT_NSIGMA")) ESCAPE (MEAS_EXT_NSIGMA,   MAG_NONE);
     153  if (!strcasecmp (fieldName, "SKY"))        ESCAPE (MEAS_SKY,          MAG_NONE);
     154  if (!strcasecmp (fieldName, "SKY_ERR"))    ESCAPE (MEAS_dSKY,         MAG_NONE);
     155  if (!strcasecmp (fieldName, "STARGAL"))    ESCAPE (MEAS_STARGAL,      MAG_NONE);
    151156
    152157  // for words that don't parse, try a photcode

  • trunk/Ohana/src/opihi/dvo/gstar.c

    r15038 r15208  
    1111  double Mcat, Mrel;
    1212  float *RA, *DEC;
    13   int i, j, k, m, N, *N1, Nsecfilt, NPTS, QUIET;
     13  int i, j, k, m, N, *N1, Nsecfilt, NPTS, QUIET, FULL_OUTPUT;
    1414  int Nstars, found, GetMeasures, Nlo, Nhi;
    1515  int SaveVectors;
    … …  
    4242  }
    4343
     44  FULL_OUTPUT = FALSE;
     45  if ((N = get_argument (argc, argv, "-full"))) {
     46    FULL_OUTPUT = TRUE;
     47    remove_argument (N, &argc, argv);
     48  }
     49
    4450  GetMeasures = FALSE;
    4551  if ((N = get_argument (argc, argv, "-m"))) {
     52    GetMeasures = TRUE;
     53    remove_argument (N, &argc, argv);
     54  }
     55  if ((N = get_argument (argc, argv, "-meas"))) {
    4656    GetMeasures = TRUE;
    4757    remove_argument (N, &argc, argv);
    … …  
    139149      if (!QUIET) {
    140150        gprint (GP_LOG, "star: %d\n", k);
    141         gprint (GP_LOG, "%11.7f %11.7f %3d of %3d  %4.1f %5d\n", catalog.average[k].R, catalog.average[k].D,
    142                  catalog.average[k].Nm, catalog.average[k].Nn + catalog.average[k].Nm,
    143                  0.01*catalog.average[k].Xp, catalog.average[k].code);
     151        gprint (GP_LOG, "%11.7f ", catalog.average[k].R);
     152        gprint (GP_LOG, "%11.7f ", catalog.average[k].D);
     153        gprint (GP_LOG, "%3d   ",  catalog.average[k].Nm);
     154        gprint (GP_LOG, "%4.1f ",  0.01*catalog.average[k].Xp);
     155        gprint (GP_LOG, "%5d",     catalog.average[k].code);
     156
     157        if (FULL_OUTPUT) {
     158            gprint (GP_LOG, "%f",     catalog.average[k].dR);
     159            gprint (GP_LOG, "%f",     catalog.average[k].dD);
     160            gprint (GP_LOG, "%f",     catalog.average[k].uR);
     161            gprint (GP_LOG, "%f",     catalog.average[k].uD);
     162            gprint (GP_LOG, "%f",     catalog.average[k].duR);
     163            gprint (GP_LOG, "%f",     catalog.average[k].duD);
     164            gprint (GP_LOG, "%f",     catalog.average[k].P);
     165            gprint (GP_LOG, "%f",     catalog.average[k].dP);
     166            gprint (GP_LOG, "%x",     catalog.average[k].objID_hi);
     167            gprint (GP_LOG, "%x",     catalog.average[k].objID_lo);
     168        }
     169
     170        gprint (GP_LOG, "\n");
    144171     
    145172        /* filter names */
    … …  
    169196          if (GetMeasures && !QUIET) {
    170197            date = ohana_sec_to_date (catalog.measure[m].t);
    171             gprint (GP_LOG, "%6.3f %6.3f %5.3f  %20s  %5.2f %5.2f %2d %3x %3d %-20s\n",
    172                      Mcat, Mrel, catalog.measure[m].dM,
    173                      date, catalog.measure[m].dR, catalog.measure[m].dD,
    174                      catalog.measure[m].dophot, catalog.measure[m].dbFlags,
    175                      catalog.measure[m].photcode, GetPhotcodeNamebyCode (catalog.measure[m].photcode));
     198            gprint (GP_LOG, "%6.3f ",  Mcat);
     199            gprint (GP_LOG, "%6.3f ",  Mrel);
     200            gprint (GP_LOG, "%5.3f  ", catalog.measure[m].dM);
     201            gprint (GP_LOG, "%20s  ",  date);
     202            gprint (GP_LOG, "%7.4f ",  catalog.measure[m].dR);
     203            gprint (GP_LOG, "%7.4f",   catalog.measure[m].dD);
     204            gprint (GP_LOG, "%2d ",    catalog.measure[m].dophot);
     205            gprint (GP_LOG, "%3x ",    catalog.measure[m].dbFlags);
     206            gprint (GP_LOG, "%3d ",    catalog.measure[m].photcode);
     207            gprint (GP_LOG, "%-20s  ", GetPhotcodeNamebyCode (catalog.measure[m].photcode));
     208
     209            if (FULL_OUTPUT) {
     210                gprint (GP_LOG, "%f", catalog.measure[m].Mcal);
     211                gprint (GP_LOG, "%f", catalog.measure[m].Mgal);
     212                gprint (GP_LOG, "%f", pow(10.0, 0.4*catalog.measure[m].dt));
     213                gprint (GP_LOG, "%f", 1.0 + catalog.measure[m].airmass);
     214                gprint (GP_LOG, "%f", catalog.measure[m].az);
     215                gprint (GP_LOG, "%f", catalog.measure[m].Xccd);
     216                gprint (GP_LOG, "%f", catalog.measure[m].Yccd);
     217                gprint (GP_LOG, "%d", catalog.measure[m].dXccd);
     218                gprint (GP_LOG, "%d", catalog.measure[m].dYccd);
     219                gprint (GP_LOG, "%f", catalog.measure[m].Sky);
     220                gprint (GP_LOG, "%f", catalog.measure[m].dSky);
     221                gprint (GP_LOG, "%d", catalog.measure[m].averef);
     222                gprint (GP_LOG, "%d", catalog.measure[m].detID);
     223                gprint (GP_LOG, "%d", catalog.measure[m].imageID);
     224                gprint (GP_LOG, "%f", catalog.measure[m].qPSF);
     225                gprint (GP_LOG, "%f", catalog.measure[m].psfProb);
     226                gprint (GP_LOG, "%f", catalog.measure[m].crNsigma);
     227                gprint (GP_LOG, "%f", catalog.measure[m].extNsigma);
     228                gprint (GP_LOG, "%f", 0.01*catalog.measure[m].FWx);
     229                gprint (GP_LOG, "%f", 0.01*catalog.measure[m].FWy);
     230                gprint (GP_LOG, "%f", (360.0/(float)0xffff)*catalog.measure[m].theta);
     231
     232                gprint (GP_LOG, "%x", catalog.measure[m].photFlags);
     233                gprint (GP_LOG, "%d", catalog.measure[m].stargal);
     234            }
     235            gprint (GP_LOG, "\n");
     236
    176237            free (date);
    177238          }

  • trunk/Ohana/src/opihi/dvo/photometry.c

    r15038 r15208  
    191191  if (!strcasecmp (parname, "fwhm"))     param = MEAS_FWHM;
    192192  // if (!strcasecmp (parname, "dophot"))        param = MEAS_DOPHOT;
    193   if (!strcasecmp (parname, "FLAGS"))    param = MEAS_FLAGS;
     193  if (!strcasecmp (parname, "FLAGS"))    param = MEAS_DB_FLAGS;
    194194  if (!strcasecmp (parname, "XCCD"))     param = MEAS_XCCD;
    195195  if (!strcasecmp (parname, "YCCD"))     param = MEAS_YCCD;
    … …  
    11681168      value = 0.01*measure[0].FWx;
    11691169      break;
    1170     case MEAS_FLAGS: /* ? */
     1170    case MEAS_DB_FLAGS: /* ? */
    11711171      value = measure[0].dbFlags;
    11721172      break;

  • trunk/Ohana/src/opihi/include/dvoshell.h

    r14590 r15208  
    5959      MEAS_FWHM_MIN,
    6060      MEAS_THETA,
    61       MEAS_FLAGS,
     61      MEAS_DOPHOT,
     62      MEAS_DB_FLAGS,
     63      MEAS_PHOT_FLAGS,
    6264      MEAS_XCCD,
    6365      MEAS_YCCD,
    6466      MEAS_XMOSAIC,
    6567      MEAS_YMOSAIC,
    66       MEAS_DETID,
    67       MEAS_OBJID,
    68       MEAS_IMAGEID,
     68      MEAS_SKY,
     69      MEAS_dSKY,
     70      MEAS_DET_ID,
     71      MEAS_OBJ_ID,
     72      MEAS_IMAGE_ID,
     73      MEAS_PSF_QF,
     74      MEAS_PSF_PROB,
     75      MEAS_CR_NSIGMA,
     76      MEAS_EXT_NSIGMA,
     77      MEAS_STARGAL,
    6978};
    7079
Note: See TracChangeset for help on using the changeset viewer.