Changeset 25020
- Timestamp:
- Aug 7, 2009, 12:42:40 PM (17 years ago)
- Location:
- branches/eam_branches/20090715/Ohana/src
- Files:
-
- 19 edited
- 6 copied
-
addstar/Makefile (modified) (2 diffs)
-
addstar/doc/notes.txt (modified) (1 diff)
-
addstar/include/2mass.h (modified) (1 diff)
-
addstar/include/addstar.h (modified) (3 diffs)
-
addstar/src/ConfigInit.c (modified) (1 diff)
-
addstar/src/GetZeroPointExposure.c (copied) (copied from trunk/Ohana/src/addstar/src/GetZeroPointExposure.c )
-
addstar/src/LoadData.c (modified) (1 diff)
-
addstar/src/ReadImageHeader.c (modified) (4 diffs)
-
addstar/src/SkyRegionUtils.c (modified) (2 diffs)
-
addstar/src/args.c (modified) (1 diff)
-
addstar/src/args_load2mass.c (modified) (2 diffs)
-
addstar/src/load2mass.c (modified) (3 diffs)
-
addstar/src/load2mass_as_rawdata.c (modified) (4 diffs)
-
getstar/Makefile (modified) (2 diffs)
-
getstar/include/dvoImagesAtCoords.h (copied) (copied from trunk/Ohana/src/getstar/include/dvoImagesAtCoords.h )
-
getstar/src/ConfigInit_coords.c (copied) (copied from trunk/Ohana/src/getstar/src/ConfigInit_coords.c )
-
getstar/src/MatchCoords.c (copied) (copied from trunk/Ohana/src/getstar/src/MatchCoords.c )
-
getstar/src/args_coords.c (copied) (copied from trunk/Ohana/src/getstar/src/args_coords.c )
-
getstar/src/dvoImagesAtCoords.c (copied) (copied from trunk/Ohana/src/getstar/src/dvoImagesAtCoords.c )
-
libdvo/include/dvo.h (modified) (3 diffs)
-
libdvo/src/dvo_photcode_ops.c (modified) (1 diff)
-
libdvo/src/skyregion_ops.c (modified) (1 diff)
-
opihi/dvo/dbExtractMeasures.c (modified) (1 diff)
-
opihi/dvo/dbFields.c (modified) (1 diff)
-
opihi/include/dvoshell.h (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/20090715/Ohana/src/addstar/Makefile
r21508 r25020 58 58 $(SRC)/greference.$(ARCH).o \ 59 59 $(SRC)/grefstars.$(ARCH).o \ 60 $(SRC)/GetZeroPointExposure.$(ARCH).o \ 60 61 $(SRC)/LoadStars.$(ARCH).o \ 61 62 $(SRC)/LoadHeaders.$(ARCH).o \ … … 199 200 $(SRC)/replace_match.$(ARCH).o \ 200 201 $(SRC)/update_coords.$(ARCH).o \ 202 $(SRC)/SkyRegionUtils.$(ARCH).o \ 201 203 $(SRC)/StarOps.$(ARCH).o \ 202 204 $(SRC)/ConfigInit.$(ARCH).o \ -
branches/eam_branches/20090715/Ohana/src/addstar/doc/notes.txt
r21508 r25020 1 2 2009.08.02 3 4 ZERO POINTS: addstar currently has a couple of inconsistent ways that it 5 measures or uses image zero points: 6 7 * the standard zero point: data in the database is saved as M_inst + 8 25.0. This value is written to the image header and checked when 9 the image is loaded. This is completely archaic, but removing 10 this concept from the DVO system will require some work and some 11 care. This value is NOT reflected in the image->Mcal values 12 13 * on-the-fly calibration: addstar was used by CFHT skyprobe to 14 measure the sky transparency. This calculation was done by 15 measuring the zero point between the Tycho stars in the database 16 and the stars in the image. Star data is saved in 17 AddToCalibration as they are found (in find_matches, or 18 equivalent). At the end of the match processing, the zero point 19 for the image was calculated with FindCalibration and saved to the 20 image table but NOT the measure table. For skyprobe, the stars 21 are normally not saved to the database -- the Tycho stars are kept 22 in the measure & average tables, and used as the reference for the 23 zero point calculation above. 24 25 For Pan-STARRS, the zero point analysis is performed by psastro 26 before we run addstar. However, there are a few ways we could 27 choose to apply the zero point calculation: 28 29 * use the per-chip zero point reported in the individual chip 30 headers (easy, but fairly wrong) 31 32 * calculate the per-exposure zero point from the chip headers and 33 supply. 34 35 * calculate the per-exposure zero point from the MATCHED_REFS table 36 37 We also need to consider the zero point error. We currently have 38 only dM (poisson error) + dMcal (calibration error). We need to 39 divide dMsys into dMcal (photometry calibration error) and dMap 40 (aperture correction error). 41 42 * zero point options: 43 * NOMINAL : set Mcal = 0.0 44 * CHIP_HEADER : each chip has zero point in header 45 * PHU_HEADER : a single exposure-wide value in the header 46 * CHIP_AVERAGE : determine zero point for exposure from chip zero points 47 * MATCHED_REFS : recalculate from table 1 48 2 49 2009.02.11 -
branches/eam_branches/20090715/Ohana/src/addstar/include/2mass.h
r17439 r25020 33 33 e_time get2mass_date (char *ptr, int Nbound, int Nmax); 34 34 35 int load2mass_as_rawdata (Sky Table*skytable, char *filename, AddstarClientOptions options);35 int load2mass_as_rawdata (SkyList *skytable, char *filename, AddstarClientOptions options); 36 36 SkyTable *load2mass_acc (char *path, char *accel); 37 37 int get2mass_3star (Stars *star, char *line, int Nmax); -
branches/eam_branches/20090715/Ohana/src/addstar/include/addstar.h
r24400 r25020 73 73 char TYCHO_DIR[256]; 74 74 char SubpixDatafile[256]; 75 char *USE_NAME; 75 76 char PASSWORD[80]; 76 77 char HOSTNAME[80]; … … 118 119 double FAKE_THETA; // boresite angle for fake images 119 120 121 char ZERO_POINT_OPTION[64]; 122 float ZERO_POINT_OFFSET; 123 120 124 // carries the mosaic into gstars 121 125 … … 187 191 HeaderSet *MatchHeaders PROTO((int **extsize, int *nimage, int mode, Header **headers, int Nheaders)); 188 192 int LoadData PROTO((FILE *f, char *file, Image **images, int *nvalid, Stars **stars, int *Nstars, Header **headers, int *extsize, HeaderSet *headerSets, int NheaderSets)); 193 int GetZeroPointExposure PROTO((Header **headers, HeaderSet *headerSets, int Nimages)); 189 194 190 195 int in_image PROTO((double r, double d, Image *image)); -
branches/eam_branches/20090715/Ohana/src/addstar/src/ConfigInit.c
r21508 r25020 77 77 if (!ScanConfig (config, "EXTNAME-KEYWORD", "%s", 0, ExtnameKeyword)) { 78 78 strcpy (ExtnameKeyword, "EXTNAME"); 79 } 80 81 // if this config variable is not set, do not set any zero point offset 82 if (!ScanConfig (config, "ZERO_POINT_OPTION", "%s", 0, ZERO_POINT_OPTION)) { 83 strcpy (ZERO_POINT_OPTION, "NOMINAL"); 79 84 } 80 85 -
branches/eam_branches/20090715/Ohana/src/addstar/src/LoadData.c
r21508 r25020 23 23 24 24 // find image rootname 25 name = filebasename (file); 25 if (USE_NAME) { 26 name = filebasename (USE_NAME); 27 } else { 28 name = filebasename (file); 29 } 30 31 // if zero points are calculated for the full exposure using more than just the matched chip header, 32 // we need to perform that analysis here 33 GetZeroPointExposure (headers, headerSets, Nimages); 26 34 27 35 // now run through the images, interpret the headers and read the stars -
branches/eam_branches/20090715/Ohana/src/addstar/src/ReadImageHeader.c
r21508 r25020 8 8 double tmp, sec, Cerror, ZeroPt, FWHM_X, FWHM_Y; 9 9 char *c, photname[64], line[80]; 10 PhotCode *photcodeData = NULL; 10 11 11 12 // zero out the entire image structure … … 104 105 return (FALSE); 105 106 } 106 photcode = GetPhotcodeCodebyName (photname);107 if (photcode == 0) {107 photcodeData = GetPhotcodebyName (photname); 108 if (photcodeData == NULL) { 108 109 fprintf (stderr, "photcode %s not found in photcode table\n", photname); 109 110 return (FALSE); 110 111 } 112 photcode = photcodeData[0].code; 111 113 } 112 114 if (photcode == 0) { … … 191 193 } 192 194 195 // XXX this is archaic: we used to set a fixed zero point of 25 to shift data into the range 196 // 0 - 32 so it would fit in an unsigned int. This is also needed because some programs like 197 // sextractor will put in an arbitrary zero point that we need to understand to get back to 198 // instrumental mags. 193 199 gfits_scan (header, "ZERO_PT", "%lf", 1, &ZeroPt); 194 200 if (ZeroPt != GetZeroPoint()) { … … 197 203 } 198 204 205 // if it exists, find the 206 if (!strcasecmp(ZERO_POINT_OPTION, "CHIP_HEADER")) { 207 float ZPT_OBS; 208 if (!photcodeData || !gfits_scan (header, "ZPT_OBS", "%f", 1, &ZPT_OBS)) { 209 fprintf (stderr, "zero point not supplied in header\n"); 210 ZERO_POINT_OFFSET = 0.0; 211 } else { 212 ZERO_POINT_OFFSET = 0.001*photcodeData[0].C - ZPT_OBS; 213 } 214 } 215 199 216 /* secz is in units milli-airmass */ 200 image[0].Mcal = 0.0;217 image[0].Mcal = ZERO_POINT_OFFSET; 201 218 image[0].Xm = NAN_S_SHORT; 202 219 image[0].flags = 0; -
branches/eam_branches/20090715/Ohana/src/addstar/src/SkyRegionUtils.c
r14401 r25020 11 11 ALLOCATE (subset, SkyList, 1); 12 12 ALLOCATE (subset[0].regions, SkyRegion *, NSUBSET); 13 ALLOCATE (subset[0].filename, char *, NSUBSET); 13 14 subset[0].ownElements = FALSE; // free these elements when freeing the list 14 15 … … 20 21 21 22 subset[0].regions[Nsubset] = input[0].regions[i]; 23 subset[0].filename[Nsubset] = input[0].filename[i]; 22 24 Nsubset ++; 23 CHECK_REALLOCATE (subset[0].regions, SkyRegion *, NSUBSET, Nsubset, 100); 25 if (Nsubset >= NSUBSET) { 26 NSUBSET += 100; 27 REALLOCATE (subset[0].regions, SkyRegion *, NSUBSET); 28 REALLOCATE (subset[0].filename, char *, NSUBSET); 29 } 24 30 subset[0].Nregions = Nsubset; 25 31 } -
branches/eam_branches/20090715/Ohana/src/addstar/src/args.c
r21508 r25020 46 46 if ((N = get_argument (argc, argv, "-list"))) { 47 47 options.filelist = TRUE; 48 remove_argument (N, &argc, argv); 49 } 50 51 USE_NAME = NULL; 52 if ((N = get_argument (argc, argv, "-use-name"))) { 53 remove_argument (N, &argc, argv); 54 USE_NAME = strcreate (argv[N]); 48 55 remove_argument (N, &argc, argv); 49 56 } -
branches/eam_branches/20090715/Ohana/src/addstar/src/args_load2mass.c
r18374 r25020 38 38 } 39 39 40 /* only add to existing regions */ 41 options.existing_regions = FALSE; 42 if ((N = get_argument (argc, argv, "-existing-regions"))) { 43 options.existing_regions = TRUE; 44 remove_argument (N, &argc, argv); 45 } 40 46 /* only add to existing objects */ 41 47 options.only_match = FALSE; … … 62 68 options.timeref = 0; 63 69 options.mosaic = FALSE; 64 options.existing_regions = FALSE;65 70 options.skip_missed = FALSE; 66 71 options.closest = FALSE; -
branches/eam_branches/20090715/Ohana/src/addstar/src/load2mass.c
r7780 r25020 7 7 int i; 8 8 SkyTable *sky, *sky2mass; 9 SkyList *skylist = NULL; 10 SkyList *overlap = NULL; 9 11 AddstarClientOptions options; 10 12 … … 13 15 options = args_load2mass (argc, argv, options); 14 16 17 // load the full sky description table: 15 18 sky = SkyTableLoadOptimal (CATDIR, SKY_TABLE, GSCFILE, SKY_DEPTH, VERBOSE); 16 19 SkyTableSetFilenames (sky, CATDIR, "cpt"); 17 20 21 // generate the subset matching the user-selected region 22 skylist = SkyListByPatch (sky, -1, &UserPatch); 23 24 // if we only match to existing (already populated) regions, limit the select to those regions: 25 if (options.existing_regions) { 26 SkyList *tmp; 27 tmp = SkyListExistingSubset (skylist, CATDIR); 28 SkyListFree (skylist); 29 skylist = tmp; 30 } 31 18 32 path = TWO_MASS_DIR_AS; 19 33 … … 23 37 24 38 for (i = 0; i < sky2mass[0].Nregions; i++) { 39 // check if any of the skylist entries overlap this 2mass catalog: 40 overlap = SkyListByBounds_List (skylist, -1, sky2mass[0].regions[i].Rmin, sky2mass[0].regions[i].Rmax, sky2mass[0].regions[i].Dmin, sky2mass[0].regions[i].Dmax); 41 if (overlap[0].Nregions == 0) { 42 SkyListFree (overlap); 43 continue; 44 } 45 25 46 fprintf (stderr, "loading %s\n", sky2mass[0].filename[i]); 26 load2mass_as_rawdata (sky, sky2mass[0].filename[i], options); 47 load2mass_as_rawdata (overlap, sky2mass[0].filename[i], options); 48 SkyListFree (overlap); 27 49 } 28 50 exit (0); -
branches/eam_branches/20090715/Ohana/src/addstar/src/load2mass_as_rawdata.c
r21508 r25020 10 10 # define DEBUG 0 11 11 12 int load2mass_as_rawdata (Sky Table*skytable, char *filename, AddstarClientOptions options) {12 int load2mass_as_rawdata (SkyList *skytable, char *filename, AddstarClientOptions options) { 13 13 14 14 int i, j, verbose; … … 109 109 if (tstars[i].D > UserPatch.Dmax) continue; 110 110 111 // identify the relevant catalog 112 skylist = SkyRegionByPoint_List (skytable, -1, tstars[i].R, tstars[i].D); 113 if (skylist[0].Nregions == 0) { 114 SkyListFree (skylist); 115 continue; 116 } 117 region = skylist[0].regions[0]; 118 if (DEBUG) fprintf (stderr, "writing to %s\n", skylist[0].filename[0]); 119 111 120 // collect array of (Stars *) stars in a new output catalog 112 121 Nstars = 0; … … 114 123 ALLOCATE (stars, Stars *, NSTARS); 115 124 116 // identify the relevant catalog 117 skylist = SkyRegionByPoint (skytable, -1, tstars[i].R, tstars[i].D); 118 region = skylist[0].regions[0]; 119 if (DEBUG) fprintf (stderr, "writing to %s\n", skylist[0].filename[0]); 120 125 // loop over stars in this 2mass region that are also in this output region 121 126 for (j = i; j < Ntstars; j++) { 122 127 if (tstars[j].flag) continue; … … 166 171 } 167 172 168 if (DEBUG) fprintf (stderr, "selected %d stars (%10.6f - %10.6f, %10.6f - %10.6f \n", Nstars,173 if (DEBUG) fprintf (stderr, "selected %d stars (%10.6f - %10.6f, %10.6f - %10.6f)\n", Nstars, 169 174 region[0].Rmin, region[0].Rmax, region[0].Dmin, region[0].Dmax); 170 175 -
branches/eam_branches/20090715/Ohana/src/getstar/Makefile
r20983 r25020 23 23 dvoImageOverlaps.install: $(DESTBIN)/dvoImageOverlaps 24 24 25 dvoImagesAtCoords: $(BIN)/dvoImagesAtCoords.$(ARCH) 26 dvoImagesAtCoords.install: $(DESTBIN)/dvoImagesAtCoords 27 25 28 dvoImageExtract: $(BIN)/dvoImageExtract.$(ARCH) 26 29 dvoImageExtract.install: $(DESTBIN)/dvoImageExtract 27 30 28 all: getstar dvoImageOverlaps dvoImageExtract 29 install: getstar.install dvoImageOverlaps.install dvoImageExtract.install 31 all: getstar dvoImageOverlaps dvoImageExtract dvoImagesAtCoords 32 install: getstar.install dvoImageOverlaps.install dvoImageExtract.install dvoImagesAtCoords.install 30 33 31 34 GETSTAR = \ … … 61 64 $(BIN)/dvoImageOverlaps.$(ARCH): $(OVERLAPS) 62 65 66 COVERLAPS = \ 67 $(SRC)/dvoImagesAtCoords.$(ARCH).o \ 68 $(SRC)/args_coords.$(ARCH).o \ 69 $(SRC)/ConfigInit_coords.$(ARCH).o \ 70 $(SRC)/Shutdown.$(ARCH).o \ 71 $(SRC)/SetSignals.$(ARCH).o \ 72 $(SRC)/ReadImageFiles.$(ARCH).o \ 73 $(SRC)/ReadImageHeader.$(ARCH).o \ 74 $(SRC)/GetFileMode.$(ARCH).o \ 75 $(SRC)/MatchCoords.$(ARCH).o 76 77 $(COVERLAPS): $(INC)/dvoImagesAtCoords.h $(DVO_INCS) 78 $(BIN)/dvoImagesAtCoords.$(ARCH): $(COVERLAPS) 79 63 80 EXTRACT = \ 64 81 $(SRC)/dvoImageExtract.$(ARCH).o \ -
branches/eam_branches/20090715/Ohana/src/libdvo/include/dvo.h
r24748 r25020 311 311 float PhotInst (Measure *measure); 312 312 float PhotCat (Measure *measure); 313 float PhotAper (Measure *measure); 313 314 float PhotSys (Measure *measure, Average *average, SecFilt *secfilt); 314 315 float PhotRel (Measure *measure, Average *average, SecFilt *secfilt); … … 441 442 SkyList *SkyListByBounds PROTO((SkyTable *table, int depth, double Rmin, double Rmax, double Dmin, double Dmax)); 442 443 SkyList *SkyListChildrenByBounds PROTO((SkyTable *table, int No, int depth, double Rmin, double Rmax, double Dmin, double Dmax)); 444 443 445 int SkyListMerge PROTO((SkyList **outlist, SkyList *newlist)); 444 446 int SkyListFree PROTO((SkyList *list)); … … 447 449 int SkyTableSetFilenames PROTO((SkyTable *sky, char *path, char *ext)); 448 450 451 SkyList *SkyRegionByPoint_List PROTO((SkyList *inList, int depth, double ra, double dec)); 452 SkyList *SkyListByBounds_List PROTO((SkyList *table, int depth, double Rmin, double Rmax, double Dmin, double Dmax)); 453 SkyList *SkyListChildrenByBounds_List PROTO((SkyList *table, int depth, double Rmin, double Rmax, double Dmin, double Dmax)); 454 449 455 /* dvo-specific sorting functions */ 450 456 void sortave (Average *ave, int N); -
branches/eam_branches/20090715/Ohana/src/libdvo/src/dvo_photcode_ops.c
r17253 r25020 249 249 code = &photcodes[0].code[Np]; 250 250 Mcat = measure[0].M - ZERO_POINT + code[0].K*(measure[0].airmass - 1.000) + SCALE*code[0].C; 251 252 return (Mcat); 253 } 254 255 float PhotAper (Measure *measure) { 256 257 int Np; 258 float Mcat; 259 PhotCode *code; 260 261 Np = photcodes[0].hashcode[measure[0].photcode]; 262 if (Np == -1) return (NAN); 263 264 if (photcodes[0].code[Np].type == PHOT_REF) { 265 Mcat = measure[0].Map; 266 return (Mcat); 267 } 268 code = &photcodes[0].code[Np]; 269 Mcat = measure[0].Map - ZERO_POINT + code[0].K*(measure[0].airmass - 1.000) + SCALE*code[0].C; 251 270 252 271 return (Mcat); -
branches/eam_branches/20090715/Ohana/src/libdvo/src/skyregion_ops.c
r24747 r25020 300 300 REALLOCATE (list[0].filename, char *, NNEW); 301 301 } 302 } 303 } 304 305 list[0].Nregions = Nnew; 306 return (list); 307 } 308 309 /* find region which overlaps c at given depth (-1 : populated ) */ 310 SkyList *SkyRegionByPoint_List (SkyList *inList, int depth, double ra, double dec) { 311 312 int i; 313 SkyRegion **region; 314 SkyList *list; 315 316 ALLOCATE (list, SkyList, 1); 317 ALLOCATE (list[0].regions, SkyRegion *, 1); 318 ALLOCATE (list[0].filename, char *, 1); 319 list[0].Nregions = 0; 320 list[0].ownElements = FALSE; // this list is only holding a view to the elements 321 322 region = inList[0].regions; 323 324 for (i = 0; i < inList[0].Nregions; i++) { 325 326 // once we pass our desired depth, quit 327 if ((depth > -1) && (depth < region[i][0].depth)) break; 328 329 // skip tables that do not overlap 330 if (ra < region[i][0].Rmin) continue; 331 if (ra > region[i][0].Rmax) continue; 332 if (dec < region[i][0].Dmin) continue; 333 if (dec > region[i][0].Dmax) continue; 334 335 // skip tables that are not at our depth 336 if ((depth > -1) && (depth > region[i][0].depth)) continue; 337 if ((depth == -1) && (region[i][0].table == FALSE)) continue; 338 339 list[0].regions[0] = region[i]; 340 list[0].filename[0] = inList[0].filename[i]; 341 list[0].Nregions = 1; 342 break; 343 } 344 return (list); 345 } 346 347 SkyList *SkyListByBounds_List (SkyList *table, int depth, double Rmin, double Rmax, double Dmin, double Dmax) { 348 349 int i, j, Ns; 350 SkyList *list, *extra; 351 352 Rmin = ohana_normalize_angle (Rmin); 353 Rmax = ohana_normalize_angle (Rmax); 354 355 /* handle 0,360 boundary requests */ 356 /* this is probably wrong: I will get duplicates for all Dec bands... */ 357 if (Rmin > Rmax) { 358 list = SkyListChildrenByBounds_List (table, depth, 0.0, Rmax, Dmin, Dmax); 359 360 extra = SkyListChildrenByBounds_List (table, depth, Rmin, 360.0, Dmin, Dmax); 361 REALLOCATE (list[0].regions, SkyRegion *, list[0].Nregions + extra[0].Nregions); 362 REALLOCATE (list[0].filename, char *, list[0].Nregions + extra[0].Nregions); 363 Ns = list[0].Nregions; 364 for (i = 0; i < extra[0].Nregions; i++) { 365 // search for pre-existing match 366 for (j = 0; j < list[0].Nregions; j++) { 367 if (list[0].regions[j] == extra[0].regions[i]) { 368 goto skip; 369 } 370 } 371 list[0].regions[Ns] = extra[0].regions[i]; 372 list[0].filename[Ns] = extra[0].filename[i]; 373 Ns ++; 374 skip: 375 continue; 376 } 377 list[0].Nregions = Ns; 378 SkyListFree (extra); 379 } else { 380 list = SkyListChildrenByBounds_List (table, depth, Rmin, Rmax, Dmin, Dmax); 381 } 382 383 return (list); 384 } 385 386 SkyList *SkyListChildrenByBounds_List (SkyList *table, int depth, double Rmin, double Rmax, double Dmin, double Dmax) { 387 388 int i, j, Ns, Ne, Nnew, NNEW; 389 int append; 390 SkyList *children; 391 SkyList *list; 392 SkyRegion **region; 393 394 Nnew = 0; 395 NNEW = 50; 396 ALLOCATE (list, SkyList, 1); 397 ALLOCATE (list[0].regions, SkyRegion *, NNEW); 398 ALLOCATE (list[0].filename, char *, NNEW); 399 list[0].ownElements = FALSE; // this list is only holding a view to the elements 400 401 region = table[0].regions; 402 403 // can we assume the skylist is globally sorted by depth? i think so... 404 for (i = 0; i < table[0].Nregions; i++) { 405 406 // once we pass our desired depth, quit 407 if ((depth > -1) && (depth < region[i][0].depth)) break; 408 409 // skip tables that do not overlap 410 if (Rmax <= region[i][0].Rmin) continue; 411 if (Rmin >= region[i][0].Rmax) continue; 412 if (Dmax <= region[i][0].Dmin) continue; 413 if (Dmin >= region[i][0].Dmax) continue; 414 415 // skip tables that are not at our depth 416 if ((depth > -1) && (depth > region[i][0].depth)) continue; 417 if ((depth == -1) && (region[i][0].table == FALSE)) continue; 418 419 list[0].regions[Nnew] = region[i]; 420 list[0].filename[Nnew] = table[0].filename[i]; 421 Nnew ++; 422 if (Nnew >= NNEW) { 423 NNEW += 50; 424 REALLOCATE (list[0].regions, SkyRegion *, NNEW); 425 REALLOCATE (list[0].filename, char *, NNEW); 302 426 } 303 427 } -
branches/eam_branches/20090715/Ohana/src/opihi/dvo/dbExtractMeasures.c
r21508 r25020 117 117 case MAG_REF: 118 118 value.Flt = PhotRef (equiv, average, secfilt, measure); 119 break; 120 case MAG_APER: 121 value.Flt = PhotAper (measure); 119 122 break; 120 123 case MAG_ERR: -
branches/eam_branches/20090715/Ohana/src/opihi/dvo/dbFields.c
r21508 r25020 37 37 if (!strcasecmp (string, "ave")) return (MAG_AVE); 38 38 if (!strcasecmp (string, "ref")) return (MAG_REF); 39 if (!strcasecmp (string, "ap")) return (MAG_APER); 40 if (!strcasecmp (string, "aper")) return (MAG_APER); 39 41 if (!strcasecmp (string, "err")) return (MAG_ERR); 40 42 if (!strcasecmp (string, "photflags")) return (MAG_PHOT_FLAGS); -
branches/eam_branches/20090715/Ohana/src/opihi/include/dvoshell.h
r21508 r25020 27 27 MAG_INST, 28 28 MAG_CAT, 29 MAG_APER, 29 30 MAG_SYS, 30 31 MAG_REL,
Note:
See TracChangeset
for help on using the changeset viewer.
