--- ray/src/common/tonemap.c 2011/04/22 14:35:54 3.36 +++ ray/src/common/tonemap.c 2021/02/24 20:35:22 3.45 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: tonemap.c,v 3.36 2011/04/22 14:35:54 greg Exp $"; +static const char RCSid[] = "$Id: tonemap.c,v 3.45 2021/02/24 20:35:22 greg Exp $"; #endif /* * Tone mapping functions. @@ -14,6 +14,7 @@ static const char RCSid[] = "$Id: tonemap.c,v 3.36 201 #include #include #include +#include #include "tmprivat.h" #include "tmerrmsg.h" @@ -66,7 +67,7 @@ double gamval tmnew->mongam = gamval; /* set color divisors */ for (i = 0; i < 3; i++) - tmnew->cdiv[i] = 256.*pow(tmnew->clf[i], 1./tmnew->mongam); + tmnew->cdiv[i] = TM_BRES*pow(tmnew->clf[i], 1./tmnew->mongam); /* set input transform */ tmnew->inppri = tmnew->monpri; @@ -142,7 +143,7 @@ MEM_PTR dat tms->cmat[i][j] *= tms->inpsf; /* set color divisors */ for (i = 0; i < 3; i++) - tms->cdiv[i] = 256.*pow(tms->clf[i] < .001 ? .001 : + tms->cdiv[i] = TM_BRES*pow(tms->clf[i] < .001 ? .001 : tms->clf[i], 1./tms->mongam); /* notify packages */ for (i = tmNumPkgs; i--; ) @@ -252,13 +253,13 @@ int tmCvColors( /* convert float colors */ TMstruct *tms, TMbright *ls, -BYTE *cs, +uby8 *cs, COLOR *scan, int len ) { static const char funcName[] = "tmCvColors"; - static BYTE gamtab[1024]; + static uby8 gamtab[1024]; static double curgam = .0; COLOR cmon; float lum, slum, d; @@ -374,12 +375,12 @@ int wt horig = HISTI(tms->hbrmin); hlen = HISTI(tms->hbrmax) + 1 - horig; if (hlen > oldlen) { /* (re)allocate histogram */ - int *newhist = (int *)calloc(hlen, sizeof(int)); + HIST_TYP *newhist = (HIST_TYP *)calloc(hlen, sizeof(HIST_TYP)); if (newhist == NULL) returnErr(TM_E_NOMEM); if (oldlen) { /* copy and free old */ - for (i = oldlen, j = i+oldorig-horig; i; ) - newhist[--j] = tms->histo[--i]; + memcpy((MEM_PTR)(newhist+(oldorig-horig)), + (MEM_PTR)tms->histo, oldlen*sizeof(HIST_TYP)); free((MEM_PTR)tms->histo); } tms->histo = newhist; @@ -435,10 +436,12 @@ double gamval gamval = tms->mongam; d = log(expmult/tms->inpsf); for (i = tms->mbrmax-tms->mbrmin+1; i--; ) { - double val = 256. * exp( + double val = TM_BRES * exp( ( d + (tms->mbrmin+i)*(1./TM_BRTSCALE) ) / gamval); - tms->lumap[i] = val >= (double)0xffff ? 0xffff : (int)val; + tms->lumap[i] = val; + if (sizeof(TMAP_TYP) == 2 && val >= 0xffff) + tms->lumap[i] = 0xffff; } returnOK; } @@ -453,11 +456,11 @@ double Ldmax ) { static const char funcName[] = "tmComputeMapping"; - int *histo; + HIST_TYP *histo; float *cumf; - int brt0, histlen, threshold, ceiling, trimmings; - double logLddyn, Ldmin, Ldavg, Lwavg, Tr, Lw, Ld; - int32 histot; + int brt0, histlen; + HIST_TYP threshold, ceiling, trimmings, histot; + double logLddyn, Ldmin, Lwavg, Tr, Lw, Ld; double sum; double d; int i, j; @@ -471,7 +474,6 @@ double Ldmax /* compute handy values */ Ldmin = Ldmax/Lddyn; logLddyn = log(Lddyn); - Ldavg = sqrt(Ldmax*Ldmin); i = HISTI(tms->hbrmin); brt0 = HISTV(i); histlen = HISTI(tms->hbrmax) + 1 - i; @@ -482,7 +484,7 @@ double Ldmax histot += tms->histo[i]; sum += (double)(j -= HISTEP) * tms->histo[i]; } - threshold = histot*0.005 + .5; + threshold = histot*0.002 + .5; if (!histot) returnErr(TM_E_TMFAIL); Lwavg = tmLuminance( (double)sum / histot ); @@ -491,22 +493,22 @@ double Ldmax tms->hbrmax - tms->hbrmin < TM_BRTSCALE*logLddyn) goto linearmap; /* clamp histogram */ - histo = (int *)malloc(histlen*sizeof(int)); + histo = (HIST_TYP *)malloc(histlen*sizeof(HIST_TYP)); cumf = (float *)malloc((histlen+2)*sizeof(float)); if ((histo == NULL) | (cumf == NULL)) returnErr(TM_E_NOMEM); cumf[histlen+1] = 1.; /* guard for assignment code */ - for (i = histlen; i--; ) /* make malleable copy */ - histo[i] = tms->histo[i]; + /* make malleable copy */ + memcpy((MEM_PTR)histo, (MEM_PTR)tms->histo, histlen*sizeof(HIST_TYP)); do { /* iterate to solution */ sum = 0; /* cumulative probability */ for (i = 0; i < histlen; i++) { - cumf[i] = (double)sum/histot; - sum += histo[i]; + cumf[i] = sum/histot; + sum += (double)histo[i]; } cumf[histlen] = 1.; Tr = histot * (double)(tms->hbrmax - tms->hbrmin) / - ((double)histlen*TM_BRTSCALE) / logLddyn; + ((double)TM_BRTSCALE*histlen*logLddyn); ceiling = Tr + 1.; trimmings = 0; /* clip to envelope */ for (i = histlen; i--; ) { @@ -538,16 +540,16 @@ double Ldmax d -= (double)j; Ld = Ldmin*exp(logLddyn*((1.-d)*cumf[j]+d*cumf[j+1])); d = (Ld - Ldmin)/(Ldmax - Ldmin); - tms->lumap[i] = 256.*pow(d, 1./gamval); + tms->lumap[i] = TM_BRES*pow(d, 1./gamval); } free((MEM_PTR)histo); /* clean up and return */ free((MEM_PTR)cumf); returnOK; linearmap: /* linear tone-mapping */ if (tms->flags & TM_F_HCONTR) - d = htcontrs(Ldavg) / htcontrs(Lwavg); + d = htcontrs(sqrt(Ldmax*Ldmin)) / htcontrs(Lwavg); else - d = Ldavg / Lwavg; + d = Ldmax / tmLuminance(tms->hbrmax); return(tmFixedMapping(tms, tms->inpsf*d/Ldmax, gamval)); } @@ -555,30 +557,36 @@ linearmap: /* linear tone-mapping */ int tmMapPixels( /* apply tone-mapping to pixel(s) */ TMstruct *tms, -BYTE *ps, +uby8 *ps, TMbright *ls, -BYTE *cs, +uby8 *cs, int len ) { static const char funcName[] = "tmMapPixels"; - int32 li, pv; + TMbright lv; + TMAP_TYP li; + int pv; if (tms == NULL || tms->lumap == NULL) returnErr(TM_E_TMINVAL); if ((ps == NULL) | (ls == NULL) | (len < 0)) returnErr(TM_E_ILLEGAL); while (len--) { - if ((li = *ls++) < tms->mbrmin) { + if ((lv = *ls++) < tms->mbrmin) { li = 0; } else { - if (li > tms->mbrmax) - li = tms->mbrmax; - li = tms->lumap[li - tms->mbrmin]; + if (lv > tms->mbrmax) + lv = tms->mbrmax; + li = tms->lumap[lv - tms->mbrmin]; } - if (cs == TM_NOCHROM) - *ps++ = li>255 ? 255 : li; - else { + if (cs == TM_NOCHROM) { +#if !(TM_BRES & 0xff) + *ps++ = li>=TM_BRES ? 255 : li/(TM_BRES>>8); +#else + *ps++ = li>=TM_BRES ? 255 : (li<<8)/TM_BRES; +#endif + } else { pv = *cs++ * li / tms->cdiv[RED]; *ps++ = pv>255 ? 255 : pv; pv = *cs++ * li / tms->cdiv[GRN]; @@ -592,7 +600,7 @@ int len TMstruct * -tmDup( /* duplicate top tone mapping */ +tmDup( /* duplicate tone mapping */ TMstruct *tms ) { @@ -608,18 +616,17 @@ TMstruct *tms *tmnew = *tms; /* copy everything */ if (tmnew->histo != NULL) { /* duplicate histogram */ len = HISTI(tmnew->hbrmax) + 1 - HISTI(tmnew->hbrmin); - tmnew->histo = (int *)malloc(len*sizeof(int)); + tmnew->histo = (HIST_TYP *)malloc(len*sizeof(HIST_TYP)); if (tmnew->histo != NULL) - for (i = len; i--; ) - tmnew->histo[i] = tms->histo[i]; + memcpy((MEM_PTR)tmnew->histo, (MEM_PTR)tms->histo, + len*sizeof(HIST_TYP)); } if (tmnew->lumap != NULL) { /* duplicate luminance mapping */ len = tmnew->mbrmax-tmnew->mbrmin+1; - tmnew->lumap = (unsigned short *)malloc( - len*sizeof(unsigned short) ); + tmnew->lumap = (TMAP_TYP *)malloc(len*sizeof(TMAP_TYP)); if (tmnew->lumap != NULL) - for (i = len; i--; ) - tmnew->lumap[i] = tms->lumap[i]; + memcpy((MEM_PTR)tmnew->lumap, (MEM_PTR)tms->lumap, + len*sizeof(TMAP_TYP)); } /* clear package data */ for (i = tmNumPkgs; i--; ) @@ -651,7 +658,7 @@ TMstruct *tms; /******************** Shared but Private library routines *********************/ -BYTE tmMesofact[BMESUPPER-BMESLOWER]; +uby8 tmMesofact[BMESUPPER-BMESLOWER]; void tmMkMesofact() /* build mesopic lookup factor table */ @@ -683,7 +690,7 @@ TMstruct *tms if (tms->mbrmin > tms->mbrmax) return 0; if (tms->lumap == NULL) - tms->lumap = (unsigned short *)malloc(sizeof(unsigned short)* + tms->lumap = (TMAP_TYP *)malloc(sizeof(TMAP_TYP)* (tms->mbrmax-tms->mbrmin+1)); return(tms->lumap != NULL); }