--- ray/src/common/tmapcolrs.c 1998/10/22 13:53:19 3.6 +++ ray/src/common/tmapcolrs.c 2003/02/22 02:07:22 3.10 @@ -1,22 +1,80 @@ -/* Copyright (c) 1998 Silicon Graphics, Inc. */ - #ifndef lint -static char SCCSid[] = "$SunId$ SGI"; +static const char RCSid[] = "$Id: tmapcolrs.c,v 3.10 2003/02/22 02:07:22 greg Exp $"; #endif - /* * Routines for tone mapping on Radiance RGBE and XYZE pictures. - * See tonemap.h for detailed function descriptions. + * + * Externals declared in tonemap.h */ +/* ==================================================================== + * The Radiance Software License, Version 1.0 + * + * Copyright (c) 1990 - 2002 The Regents of the University of California, + * through Lawrence Berkeley National Laboratory. All rights reserved. + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in + * the documentation and/or other materials provided with the + * distribution. + * + * 3. The end-user documentation included with the redistribution, + * if any, must include the following acknowledgment: + * "This product includes Radiance software + * (http://radsite.lbl.gov/) + * developed by the Lawrence Berkeley National Laboratory + * (http://www.lbl.gov/)." + * Alternately, this acknowledgment may appear in the software itself, + * if and wherever such third-party acknowledgments normally appear. + * + * 4. The names "Radiance," "Lawrence Berkeley National Laboratory" + * and "The Regents of the University of California" must + * not be used to endorse or promote products derived from this + * software without prior written permission. For written + * permission, please contact radiance@radsite.lbl.gov. + * + * 5. Products derived from this software may not be called "Radiance", + * nor may "Radiance" appear in their name, without prior written + * permission of Lawrence Berkeley National Laboratory. + * + * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED + * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES + * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE + * DISCLAIMED. IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR + * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, + * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT + * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF + * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND + * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, + * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT + * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF + * SUCH DAMAGE. + * ==================================================================== + * + * This software consists of voluntary contributions made by many + * individuals on behalf of Lawrence Berkeley National Laboratory. For more + * information on Lawrence Berkeley National Laboratory, please see + * . + */ + #include +#include #include +#include #include "tmprivat.h" #include "resolu.h" +#ifndef TM_PIC_CTRANS +#define TM_PIC_CTRANS 1 /* transform colors? (expensive) */ +#endif -extern char *tempbuffer(); - #define GAMTSZ 1024 typedef struct { @@ -25,9 +83,13 @@ typedef struct { TMbright inpsfb; /* encoded tm->inpsf */ } COLRDATA; +#ifdef NOPROTO static MEM_PTR colrInit(); static void colrNewSpace(); -extern void free(); +#else +static MEM_PTR colrInit(struct tmStruct *); +static void colrNewSpace(struct tmStruct *); +#endif static struct tmPackage colrPkg = { /* our package functions */ colrInit, colrNewSpace, free }; @@ -35,7 +97,6 @@ static int colrReg = -1; /* our package registration #define LOGISZ 260 static TMbright logi[LOGISZ]; -static BYTE photofact[BMESUPPER-BMESLOWER]; int @@ -54,7 +115,11 @@ int len; returnErr(TM_E_TMINVAL); if (ls == NULL | scan == NULL | len < 0) returnErr(TM_E_ILLEGAL); +#if TM_PIC_CTRANS if (tmNeedMatrix(tmTop)) { /* need floating point */ +#else + if (tmTop->inppri == TM_XYZPRIM) { /* no way around this */ +#endif register COLOR *newscan; newscan = (COLOR *)tempbuffer(len*sizeof(COLOR)); if (newscan == NULL) @@ -70,11 +135,8 @@ int len; for (i = 256; i--; ) logi[i] = TM_BRTSCALE*log((i+.5)/256.) - .5; for (i = 256; i < LOGISZ; i++) - logi[i] = logi[255]; - for (i = BMESLOWER; i < BMESUPPER; i++) - photofact[i-BMESLOWER] = 256. * - (tmLuminance(i) - LMESLOWER) / - (LMESUPPER - LMESLOWER); + logi[i] = 0; + tmMkMesofact(); } if ((cd = (COLRDATA *)tmPkgData(tmTop,colrReg)) == NULL) returnErr(TM_E_NOMEM); @@ -84,29 +146,31 @@ int len; li = ( cd->clfb[RED]*cmon[RED] + cd->clfb[GRN]*cmon[GRN] + cd->clfb[BLU]*cmon[BLU] ) >> 8; - bi = BRT2SCALE*(cmon[EXP]-COLXS) + + bi = BRT2SCALE(cmon[EXP]-COLXS) + logi[li] + cd->inpsfb; - if (bi < MINBRT) { - bi = MINBRT-1; /* bogus value */ - li++; /* avoid li==0 */ + if (li <= 0) { + bi = TM_NOBRT; /* bogus value */ + li = 1; /* avoid li==0 */ } ls[i] = bi; if (cs == TM_NOCHROM) /* no color? */ continue; - if (tmTop->flags & TM_F_BW) - cmon[RED] = cmon[GRN] = cmon[BLU] = li; /* mesopic adj. */ if (tmTop->flags & TM_F_MESOPIC && bi < BMESUPPER) { register int pf, sli = normscot(cmon); if (bi < BMESLOWER) cmon[RED] = cmon[GRN] = cmon[BLU] = sli; else { - pf = photofact[bi-BMESLOWER]; + if (tmTop->flags & TM_F_BW) + cmon[RED] = cmon[GRN] = cmon[BLU] = li; + pf = tmMesofact[bi-BMESLOWER]; sli *= 256 - pf; cmon[RED] = ( sli + pf*cmon[RED] ) >> 8; cmon[GRN] = ( sli + pf*cmon[GRN] ) >> 8; cmon[BLU] = ( sli + pf*cmon[BLU] ) >> 8; } + } else if (tmTop->flags & TM_F_BW) { + cmon[RED] = cmon[GRN] = cmon[BLU] = li; } bi = ( (int4)GAMTSZ*cd->clfb[RED]*cmon[RED]/li ) >> 8; cs[3*i ] = bi>=GAMTSZ ? 255 : cd->gamb[bi]; @@ -146,17 +210,18 @@ register struct radhead *rh; rh->format = FMTCIE; else rh->format = FMTBAD; - return; + return(0); } if (isexpos(s)) { rh->expos *= exposval(s); - return; + return(0); } if (isprims(s)) { primsval(rh->mypri, s); rh->primp = rh->mypri; - return; + return(0); } + return(0); } @@ -186,7 +251,7 @@ FILE *fp; *lpp = NULL; if (cpp != TM_NOCHROMP) *cpp = NULL; info = rhdefault; /* get our header */ - getheader(inpf, headline, (MEM_PTR)&info); + getheader(inpf, headline, (char *)&info); if (info.format == FMTBAD | info.expos <= 0. || fgetresolu(xp, yp, inpf) < 0) { err = TM_E_BADFILE; goto done; @@ -239,6 +304,7 @@ done: /* clean up */ } +#ifdef PCOND int /* run pcond to map picture */ dopcond(psp, xp, yp, flags, monpri, gamval, Lddyn, Ldmax, fname) BYTE **psp; @@ -249,7 +315,7 @@ double gamval, Lddyn, Ldmax; char *fname; { char *funcName = fname; - char cmdbuf[512]; + char cmdbuf[1024]; FILE *infp; register COLR *scan; register BYTE *rp; @@ -259,26 +325,26 @@ char *fname; if (setcolrcor(pow, 1./gamval) < 0) returnErr(TM_E_NOMEM); /* create command */ - strcpy(cmdbuf, "pcond "); + strcpy(cmdbuf, PCOND); if (flags & TM_F_HCONTR) - strcat(cmdbuf, "-s "); + strcat(cmdbuf, " -s"); if (flags & TM_F_MESOPIC) - strcat(cmdbuf, "-c "); + strcat(cmdbuf, " -c"); if (flags & TM_F_LINEAR) - strcat(cmdbuf, "-l "); + strcat(cmdbuf, " -l"); if (flags & TM_F_ACUITY) - strcat(cmdbuf, "-a "); + strcat(cmdbuf, " -a"); if (flags & TM_F_VEIL) - strcat(cmdbuf, "-v "); + strcat(cmdbuf, " -v"); if (flags & TM_F_CWEIGHT) - strcat(cmdbuf, "-w "); - sprintf(cmdbuf+strlen(cmdbuf), - "-p %f %f %f %f %f %f %f %f -d %f -u %f %s", - monpri[RED][CIEX], monpri[RED][CIEY], - monpri[GRN][CIEX], monpri[GRN][CIEY], - monpri[BLU][CIEX], monpri[BLU][CIEY], - monpri[WHT][CIEX], monpri[WHT][CIEY], - Lddyn, Ldmax, fname); + strcat(cmdbuf, " -w"); + if (monpri != stdprims) + sprintf(cmdbuf+strlen(cmdbuf), " -p %f %f %f %f %f %f %f %f", + monpri[RED][CIEX], monpri[RED][CIEY], + monpri[GRN][CIEX], monpri[GRN][CIEY], + monpri[BLU][CIEX], monpri[BLU][CIEY], + monpri[WHT][CIEX], monpri[WHT][CIEY]); + sprintf(cmdbuf+strlen(cmdbuf), " -d %f -u %f %s", Lddyn, Ldmax, fname); /* start pcond */ if ((infp = popen(cmdbuf, "r")) == NULL) returnErr(TM_E_BADFILE); @@ -322,6 +388,7 @@ char *fname; pclose(infp); returnOK; } +#endif int /* map a Radiance picture */ @@ -347,10 +414,12 @@ FILE *fp; if (Lddyn < MINLDDYN) Lddyn = DEFLDDYN; if (Ldmax < MINLDMAX) Ldmax = DEFLDMAX; if (flags & TM_F_BW) monpri = stdprims; +#ifdef PCOND /* check for pcond run */ if (fp == TM_GETFILE && flags & TM_F_UNIMPL) return( dopcond(psp, xp, yp, flags, monpri, gamval, Lddyn, Ldmax, fname) ); +#endif /* initialize tone mapping */ if (tmInit(flags, monpri, gamval) == NULL) returnErr(TM_E_NOMEM);