Changeset 31275
- Timestamp:
- Apr 12, 2011, 3:48:58 PM (15 years ago)
- Location:
- branches/eam_branches/ipp-20110404/Ohana/src/relphot/src
- Files:
-
- 3 edited
-
ImageOps.c (modified) (1 diff)
-
MosaicOps.c (modified) (6 diffs)
-
StarOps.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20110404/Ohana/src/relphot/src/ImageOps.c
r31273 r31275 126 126 } 127 127 } 128 fprintf (stderr, "matched %d detections to images\n", Nmatch);128 // fprintf (stderr, "matched %d detections to images\n", Nmatch); 129 129 } 130 130 -
branches/eam_branches/ipp-20110404/Ohana/src/relphot/src/MosaicOps.c
r31273 r31275 285 285 } 286 286 } 287 fprintf (stderr, "matched %d detections to mosaics\n", Nmatch);287 // fprintf (stderr, "matched %d detections to mosaics\n", Nmatch); 288 288 return (TRUE); 289 289 } … … 432 432 mark = (N < IMAGE_TOOFEW) || (N < IMAGE_GOOD_FRACTION*Nlist[i]); 433 433 if (mark) { 434 if ( 1 ||VERBOSE2) { fprintf (stderr, "marked mosaic %s ("OFF_T_FMT"), ("OFF_T_FMT" < %d) || ("OFF_T_FMT" < %f*"OFF_T_FMT")\n", image[imlist[i][0]].name, i, N, IMAGE_TOOFEW, N, IMAGE_GOOD_FRACTION, Nlist[i]); }434 if (VERBOSE2) { fprintf (stderr, "marked mosaic %s ("OFF_T_FMT"), ("OFF_T_FMT" < %d) || ("OFF_T_FMT" < %f*"OFF_T_FMT")\n", image[imlist[i][0]].name, i, N, IMAGE_TOOFEW, N, IMAGE_GOOD_FRACTION, Nlist[i]); } 435 435 mosaic[i].flags |= ID_IMAGE_PHOTOM_FEW; 436 436 Nfew ++; … … 499 499 imageOffset[i] = 0.0; 500 500 501 if ( fabs(mosaic[i].Mcal) < CLOUD_TOLERANCE) {501 if (VERBOSE2 && (fabs(mosaic[i].Mcal) < CLOUD_TOLERANCE)) { 502 502 fprintf (stderr, "cloud-free: %s : %f\n", image[imlist[i][0]].name, mosaic[i].Mcal); 503 503 } 504 if ( mosaic[i].Mcal < -CLOUD_TOLERANCE) {504 if (VERBOSE2 && (mosaic[i].Mcal < -CLOUD_TOLERANCE)) { 505 505 imageOffset[i] = -mosaic[i].Mcal; 506 506 // NOTE the negative sign: down below, we are going to add in the negative of Mcal … … 508 508 fprintf (stderr, "anti-clouds: %s : %f\n", image[imlist[i][0]].name, mosaic[i].Mcal); 509 509 } 510 if ( mosaic[i].Mcal > CLOUD_TOLERANCE) {510 if (VERBOSE2 && (mosaic[i].Mcal > CLOUD_TOLERANCE)) { 511 511 fprintf (stderr, "cloudy : %s : %f\n", image[imlist[i][0]].name, mosaic[i].Mcal); 512 512 } … … 626 626 } 627 627 imageOffset[i] = dM / NSlist[i]; 628 fprintf (stderr, "adjust image: %s : (%f %d) : %f\n", image[imlist[i][0]].name, dM, NSlist[i], imageOffset[i]);629 } 630 631 for (i = 0; i < Nmosaic; i++) {632 fprintf (stderr, "correction: %s : %f\n", image[imlist[i][0]].name, imageOffset[i]);633 }628 // fprintf (stderr, "adjust image: %s : (%f %d) : %f\n", image[imlist[i][0]].name, dM, NSlist[i], imageOffset[i]); 629 } 630 631 // for (i = 0; i < Nmosaic; i++) { 632 // fprintf (stderr, "correction: %s : %f\n", image[imlist[i][0]].name, imageOffset[i]); 633 // } 634 634 635 635 // apply all offset values to the mosaics … … 747 747 n++; 748 748 } 749 fprintf (stderr, "Nmosaic: "OFF_T_FMT", n: "OFF_T_FMT"\n", Nmosaic, n);749 // fprintf (stderr, "Nmosaic: "OFF_T_FMT", n: "OFF_T_FMT"\n", Nmosaic, n); 750 750 751 751 liststats (list, dlist, n, &stats); -
branches/eam_branches/ipp-20110404/Ohana/src/relphot/src/StarOps.c
r31273 r31275 598 598 float Mcal, Mmos, Mgrid; 599 599 StatType stats; 600 int N1, N2, N3, N4, N0;601 N1 = N2 = N3 = N4 = N0 = 0;600 // int N1, N2, N3, N4, N0; 601 // N1 = N2 = N3 = N4 = N0 = 0; 602 602 603 603 int Nsecfilt = GetPhotcodeNsecfilt (); … … 616 616 617 617 /* calculate the average value for a single star */ 618 if (catalog[i].secfilt[Nsecfilt*j+Nsec].flags & STAR_BAD) { N1++;continue; }618 if (catalog[i].secfilt[Nsecfilt*j+Nsec].flags & STAR_BAD) { continue; } 619 619 m = catalog[i].averageT[j].measureOffset; 620 620 … … 622 622 for (k = 0; k < catalog[i].averageT[j].Nmeasure; k++, m++) { 623 623 int ecode = GetPhotcodeEquivCodebyCode (catalog[i].measureT[m].photcode); 624 if (ecode != seccode) { N0++;continue;}624 if (ecode != seccode) { continue;} 625 625 Mcal = getMcal (m, i); 626 if (isnan(Mcal)) { N2++;continue;}626 if (isnan(Mcal)) { continue;} 627 627 Mmos = getMmos (m, i); 628 if (isnan(Mmos)) { N3++;continue; }628 if (isnan(Mmos)) { continue; } 629 629 Mgrid = getMgrid (m, i); 630 if (isnan(Mgrid)) { N4++;continue;}630 if (isnan(Mgrid)) { continue;} 631 631 N++; 632 632 } … … 638 638 } 639 639 640 fprintf (stderr, "N1: %d, N2: %d, N3: %d, N4: %d, N0: %d\n", N1, N2, N3, N4, N0);640 // fprintf (stderr, "N1: %d, N2: %d, N3: %d, N4: %d, N0: %d\n", N1, N2, N3, N4, N0); 641 641 liststats (list, dlist, n, &stats); 642 642 free (list);
Note:
See TracChangeset
for help on using the changeset viewer.
