Changeset 21016
- Timestamp:
- Dec 16, 2008, 2:22:25 PM (17 years ago)
- File:
-
- 1 edited
-
trunk/ppStack/src/ppStackSources.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppStack/src/ppStackSources.c
r20995 r21016 87 87 } 88 88 89 zp->data.F32[i] = airmassTerm * airmass +2.5 * log10(exptime);89 zp->data.F32[i] = airmassTerm * airmass - 2.5 * log10(exptime); 90 90 sumExpTime += exptime; 91 91 } … … 131 131 } 132 132 133 // M = m + c0 + c1 * airmass +2.5log(t) + transparency133 // M = m + c0 + c1 * airmass - 2.5log(t) + transparency 134 134 // Want sources to have m corresponding to airmass = 1 and t = sumExpTime and transparency = 0 135 135 // m_0 + c1 * airmass_0 + 2.5log(t_0) - trans_0 = m_1 + c1 * airmass_1 + 2.5log(t_1) - trans_1 … … 138 138 for (int i = 0; i < num; i++) { 139 139 psArray *sources = sourceLists->data[i]; // Sources of interest 140 float magCorr = airmassTerm +2.5*log10(sumExpTime) - zp->data.F32[i] - trans->data.F32[i];140 float magCorr = airmassTerm - 2.5*log10(sumExpTime) - zp->data.F32[i] - trans->data.F32[i]; 141 141 psLogMsg("ppStack", PS_LOG_INFO, "Applying magnitude correction to image %d: %f\n", i, magCorr); 142 142 … … 150 150 } 151 151 psFree(trans); 152 153 152 154 153 #ifdef TESTING
Note:
See TracChangeset
for help on using the changeset viewer.
