Changeset 10801
- Timestamp:
- Dec 17, 2006, 12:49:47 AM (20 years ago)
- Location:
- trunk/psphot/src
- Files:
-
- 13 edited
-
psphotApResid.c (modified) (2 diffs)
-
psphotBasicDeblend.c (modified) (1 diff)
-
psphotBlendFit.c (modified) (1 diff)
-
psphotChoosePSF.c (modified) (5 diffs)
-
psphotEnsemblePSF.c (modified) (3 diffs)
-
psphotFindPeaks.c (modified) (6 diffs)
-
psphotGuessModels.c (modified) (2 diffs)
-
psphotImageMedian.c (modified) (4 diffs)
-
psphotMagnitudes.c (modified) (1 diff)
-
psphotReplaceUnfit.c (modified) (1 diff)
-
psphotRoughClass.c (modified) (2 diffs)
-
psphotSkyReplace.c (modified) (2 diffs)
-
psphotSourceStats.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotApResid.c
r10555 r10801 115 115 Npsf ++; 116 116 } 117 psLogMsg ("psphot.apresid", 4, "measure aperture residuals : %f secfor %d objects (%d skipped, %d failed, %ld invalid)\n",117 psLogMsg ("psphot.apresid", PS_LOG_DETAIL, "measure aperture residuals for %d objects (%d skipped, %d failed, %ld invalid)\n", 118 118 psTimerMark ("psphot"), Npsf, Nskip, Nfail, sources->n - Npsf - Nskip - Nfail); 119 119 … … 333 333 psMetadataAdd (recipe, PS_LIST_TAIL, "NAPMIFIT", PS_DATA_S32 | PS_META_REPLACE, "ap residual scatter", psf->nApResid); 334 334 335 psLogMsg ("psphot.apresid", 3, "measure full-frame aperture residual: %f sec\n", psTimerMark ("psphot"));336 psLogMsg ("psphot.apresid", 4, "aperture residual: %f +/- %f : %f bias, %f skysat (%d of %d used)\n",335 psLogMsg ("psphot.apresid", PS_LOG_INFO, "measure full-frame aperture residual: %f sec\n", psTimerMark ("psphot")); 336 psLogMsg ("psphot.apresid", PS_LOG_DETAIL, "aperture residual: %f +/- %f : %f bias, %f skysat (%d of %d used)\n", 337 337 psf->ApResid, psf->dApResid, psf->skyBias, psf->skySat, Nkeep, Npsf); 338 psLogMsg ("psphot.apresid", 4, "apresid trends: %f %f %f %f %f\n",338 psLogMsg ("psphot.apresid", PS_LOG_MINUTIA, "apresid trends: %f %f %f %f %f\n", 339 339 1e3*psf->ApTrend->coeff[1][0][0][0], 340 340 1e6*psf->ApTrend->coeff[2][0][0][0], -
trunk/psphot/src/psphotBasicDeblend.c
r9734 r10801 120 120 psFree (contour); 121 121 } 122 psLogMsg ("psphot.deblend", 3, "identified %d blended objects (%f sec)\n", Nblend, psTimerMark ("psphot"));122 psLogMsg ("psphot.deblend", PS_LOG_INFO, "identified %d blended objects (%f sec)\n", Nblend, psTimerMark ("psphot")); 123 123 124 124 psFree (SN); -
trunk/psphot/src/psphotBlendFit.c
r9888 r10801 86 86 } 87 87 88 psLogMsg ("psphot.psphotBlendFit", 3, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %ld skipped)\n", psTimerMark ("psphot"), Nfit, Npsf, Next, Nfail, sources->n - Nfit);88 psLogMsg ("psphot.psphotBlendFit", PS_LOG_INFO, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %ld skipped)\n", psTimerMark ("psphot"), Nfit, Npsf, Next, Nfail, sources->n - Nfit); 89 89 return (true); 90 90 } -
trunk/psphot/src/psphotChoosePSF.c
r10267 r10801 75 75 if (source->mode & PM_SOURCE_MODE_PSFSTAR) psArrayAdd (stars, 200, source); 76 76 } 77 psLogMsg ("psphot.pspsf", 4, "selected candidate %ld PSF objects\n", stars->n);77 psLogMsg ("psphot.pspsf", PS_LOG_DETAIL, "selected candidate %ld PSF objects\n", stars->n); 78 78 79 79 if (stars->n == 0) { … … 173 173 psVectorStats (stats, dSx, NULL, NULL, 0); 174 174 float dSxo = stats->sampleMean; 175 psLogMsg ("psphot.choosePSF", 3, "dSx: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSx->n);175 psLogMsg ("psphot.choosePSF", PS_LOG_INFO, "dSx: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSx->n); 176 176 177 177 psVectorStats (stats, dSy, NULL, NULL, 0); 178 178 float dSyo = stats->sampleMean; 179 psLogMsg ("psphot.choosePSF", 3, "dSy: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSy->n);179 psLogMsg ("psphot.choosePSF", PS_LOG_INFO, "dSy: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSy->n); 180 180 181 181 for (int i = 0; i < dSN->n; i++) { … … 186 186 psVectorStats (stats, dSx, NULL, NULL, 0); 187 187 float nSx = stats->sampleStdev; 188 psLogMsg ("psphot.choosePSF", 3, "dSx: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSx->n);188 psLogMsg ("psphot.choosePSF", PS_LOG_INFO, "dSx: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSx->n); 189 189 psVectorStats (stats, dSy, NULL, NULL, 0); 190 190 float nSy = stats->sampleStdev; 191 psLogMsg ("psphot.choosePSF", 3, "dSy: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSy->n);191 psLogMsg ("psphot.choosePSF", PS_LOG_INFO, "dSy: mean: %f, median: %f, stdev: %f, npts: %ld", stats->sampleMean, stats->sampleMedian, stats->sampleStdev, dSy->n); 192 192 193 193 psFree (Sx); … … 234 234 psMetadataConfigWrite (psfData, "psfmodel.dat"); 235 235 psFree (psfData); 236 psLogMsg ("psphot.choosePSF", 3, "wrote out psf-subtracted image, psf data, exiting\n");236 psLogMsg ("psphot.choosePSF", PS_LOG_INFO, "wrote out psf-subtracted image, psf data, exiting\n"); 237 237 exit (0); 238 238 } … … 247 247 248 248 modelName = pmModelGetType (psf->type); 249 psLogMsg ("psphot.pspsf", 3, "select psf model: %f sec\n", psTimerMark ("psphot"));250 psLogMsg ("psphot.pspsf", 3, "selectedpsf model %s, ApResid: %f +/- %f\n", modelName, psf->ApResid, psf->dApResid);249 psLogMsg ("psphot.pspsf", PS_LOG_INFO, "select psf model: %f sec\n", psTimerMark ("psphot")); 250 psLogMsg ("psphot.pspsf", PS_LOG_INFO, "psf model %s, ApResid: %f +/- %f\n", modelName, psf->ApResid, psf->dApResid); 251 251 252 252 return (psf); -
trunk/psphot/src/psphotEnsemblePSF.c
r10171 r10801 105 105 psFree (fitSource); 106 106 } 107 psLogMsg ( __func__, 4, "built fitSources: %f (%ld objects)\n", psTimerMark ("psphot"), refSources->n);107 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built fitSources: %f (%ld objects)\n", psTimerMark ("psphot"), refSources->n); 108 108 109 109 // vectors to store stats for each object … … 178 178 } 179 179 psSparseResort (sparse); 180 psLogMsg ( __func__, 4, "built matrix: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem);180 psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built matrix: %f (%d elements)\n", psTimerMark ("psphot"), sparse->Nelem); 181 181 182 182 // set the sky, sky_x, sky_y components of border matrix … … 242 242 psFree (border); 243 243 244 psLogMsg ( __func__, 3, "measure ensemble of PSFs: %f sec\n", psTimerMark ("psphot"));244 psLogMsg ("psphot.ensemble", PS_LOG_INFO, "measure ensemble of PSFs: %f sec\n", psTimerMark ("psphot")); 245 245 return true; 246 246 } -
trunk/psphot/src/psphotFindPeaks.c
r10295 r10801 20 20 psImage *smooth_im = psImageCopy (NULL, readout->image, PS_TYPE_F32); 21 21 psImageSmooth (smooth_im, SIGMA, NSIGMA); 22 psLogMsg ("psphot", 4, "smooth image: %f sec\n", psTimerMark ("psphot"));22 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth image: %f sec\n", psTimerMark ("psphot")); 23 23 24 24 psImage *smooth_wt = psImageCopy (NULL, readout->weight, PS_TYPE_F32); … … 26 26 smooth_wt = (psImage*)psBinaryOp(smooth_wt, smooth_wt, "*", 27 27 psScalarAlloc(1/(4*M_PI*PS_SQR(SIGMA)), PS_TYPE_F32)); 28 psLogMsg ("psphot", 4, "smooth weight: %f sec\n", psTimerMark ("psphot"));28 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth weight: %f sec\n", psTimerMark ("psphot")); 29 29 30 30 psImage *mask = readout->mask; … … 36 36 } 37 37 38 psTimerStart ("psphot");39 38 for (int j = 0; j < smooth_im->numRows; j++) { 40 39 for (int i = 0; i < smooth_im->numCols; i++) { … … 47 46 } 48 47 } 49 psLogMsg ("psphot", 4, "built S/N map: %f sec\n", psTimerMark ("psphot")); 50 51 psTimerStart ("psphot"); 48 psLogMsg ("psphot", PS_LOG_MINUTIA, "built S/N map: %f sec\n", psTimerMark ("psphot")); 52 49 53 50 // set peak threshold … … 55 52 // threshold = NSIGMA*sky->sampleStdev + sky->sampleMean; 56 53 threshold = PS_SQR(NSIGMA); 57 psLogMsg ("psphot", 4, "threshold: %f DN\n", threshold);54 psLogMsg ("psphot", PS_LOG_DETAIL, "threshold: %f DN\n", threshold); 58 55 59 56 // find the peaks in the smoothed image … … 68 65 pmPeaksWriteText (peaks, output); 69 66 } 70 psLogMsg ("psphot", 3, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("psphot"));67 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("psphot")); 71 68 72 69 return (peaks); -
trunk/psphot/src/psphotGuessModels.c
r10052 r10801 3 3 // construct an initial PSF model for each object 4 4 bool psphotGuessModels (pmReadout *readout, psArray *sources, psMetadata *recipe, pmPSF *psf) { 5 6 psTimerStart ("psphot"); 5 7 6 8 // setup the PSF fit radius details … … 35 37 source->modelPSF = modelPSF; 36 38 } 37 psLogMsg ("psphot. emsemble", 4, "built models: %f (%ld objects)\n", psTimerMark ("psphot"), sources->n);39 psLogMsg ("psphot.models", 4, "built models: %f (%ld objects)\n", psTimerMark ("psphot"), sources->n); 38 40 return true; 39 41 } -
trunk/psphot/src/psphotImageMedian.c
r10553 r10801 199 199 } 200 200 } 201 psLogMsg ("psphot", 3, "build median image: %f sec\n", psTimerMark ("psphot"));201 psLogMsg ("psphot", PS_LOG_MINUTIA, "build median image: %f sec\n", psTimerMark ("psphot")); 202 202 203 203 // measure background stats and save for later output … … 210 210 psMetadataAdd (recipe, PS_LIST_TAIL, "SKY_NX", PS_DATA_S32 | PS_META_REPLACE, "sky model size (x)", nx); 211 211 psMetadataAdd (recipe, PS_LIST_TAIL, "SKY_NY", PS_DATA_S32 | PS_META_REPLACE, "sky model size (y)", ny); 212 psLogMsg ("psphot", 3, "background sky : min %f mean %f max %f stdev %f", statsBck->min, statsBck->sampleMean, statsBck->max, statsBck->sampleStdev); 212 psLogMsg ("psphot", PS_LOG_INFO, "background sky : min %f mean %f max %f stdev %f", 213 statsBck->min, statsBck->sampleMean, statsBck->max, statsBck->sampleStdev); 213 214 psFree (statsBck); 214 215 … … 232 233 // linear interpolation to full-scale 233 234 psImageUnbin (background->image, model->image, DX, DY, dx, dy); 234 psLogMsg ("psphot", 3, "build resampled image: %f sec\n", psTimerMark ("psphot"));235 psLogMsg ("psphot", PS_LOG_MINUTIA, "build resampled image: %f sec\n", psTimerMark ("psphot")); 235 236 236 237 // back-sub image pixels, from output background file (don't create if not requested) … … 267 268 # endif 268 269 269 psLogMsg ("psphot", 3, "subtracted background model: %f sec\n", psTimerMark ("psphot"));270 psLogMsg ("psphot", PS_LOG_INFO, "subtracted background model: %f sec\n", psTimerMark ("psphot")); 270 271 psFree(stats); 271 272 psFree(rng); -
trunk/psphot/src/psphotMagnitudes.c
r10185 r10801 48 48 } 49 49 50 psLogMsg ("psphot.magnitudes", 4, "measure magnitudes : %f sec for %ld objects (%d with apertures)\n", psTimerMark ("psphot"), sources->n, Nap);50 psLogMsg ("psphot.magnitudes", PS_LOG_DETAIL, "measure magnitudes : %f sec for %ld objects (%d with apertures)\n", psTimerMark ("psphot"), sources->n, Nap); 51 51 return true; 52 52 } -
trunk/psphot/src/psphotReplaceUnfit.c
r10201 r10801 50 50 source->mode &= ~PM_SOURCE_MODE_SUBTRACTED; 51 51 } 52 psLogMsg ("psphot.replace", 3, "replace models: %f sec (%ld objects)\n", psTimerMark ("psphot"), sources->n);52 psLogMsg ("psphot.replace", PS_LOG_INFO, "replace models: %f sec (%ld objects)\n", psTimerMark ("psphot"), sources->n); 53 53 return true; 54 54 } -
trunk/psphot/src/psphotRoughClass.c
r10053 r10801 5 5 6 6 pmPSFClump psfClump; 7 8 psTimerStart ("psphot"); 7 9 8 10 psfClump = pmSourcePSFClump (sources, recipe); … … 22 24 psphotDumpMoments (recipe, sources); 23 25 26 psLogMsg ("psphot.roughclass", PS_LOG_INFO, "rough classification (%f sec)\n", psTimerMark ("psphot")); 27 24 28 return true; 25 29 } -
trunk/psphot/src/psphotSkyReplace.c
r9270 r10801 4 4 // in order to successfully replace the sky, we must define a corresponding file... 5 5 bool psphotSkyReplace (pmConfig *config, pmFPAview *view) { 6 7 psTimerStart ("psphot"); 6 8 7 9 // find the currently selected readout … … 26 28 } 27 29 } 30 psLogMsg ("psphot.sky", PS_LOG_DETAIL, "replace background flux : %f sec\n", psTimerMark ("psphot")); 28 31 return true; 29 32 } -
trunk/psphot/src/psphotSourceStats.c
r10053 r10801 74 74 } 75 75 76 psLogMsg ("psphot", 3, "%ld moments: %f sec\n", sources->n, psTimerMark ("psphot"));76 psLogMsg ("psphot", PS_LOG_INFO, "%ld moments: %f sec\n", sources->n, psTimerMark ("psphot")); 77 77 78 78 return (sources);
Note:
See TracChangeset
for help on using the changeset viewer.
