--- ray/src/px/ra_tiff.c 2003/07/16 01:32:53 2.24 +++ ray/src/px/ra_tiff.c 2019/07/19 17:37:56 2.36 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: ra_tiff.c,v 2.24 2003/07/16 01:32:53 greg Exp $"; +static const char RCSid[] = "$Id: ra_tiff.c,v 2.36 2019/07/19 17:37:56 greg Exp $"; #endif /* * Program to convert between RADIANCE and TIFF files. @@ -7,10 +7,11 @@ static const char RCSid[] = "$Id: ra_tiff.c,v 2.24 200 * Added white-balance adjustment 10/01 (GW). */ -#include #include #include -#include + +#include "rtio.h" +#include "platform.h" #include "tiffio.h" #include "color.h" #include "resolu.h" @@ -28,6 +29,8 @@ static const char RCSid[] = "$Id: ra_tiff.c,v 2.24 200 #define C_TWRD 0x80 /* TIFF data is 16-bit */ #define C_PRIM 0x100 /* has assigned primaries */ +typedef void colcvf_t(uint32 y); + struct { uint16 flags; /* conversion flags (defined above) */ char capdate[20]; /* capture date/time */ @@ -56,7 +59,7 @@ struct { float *fp; /* float pointer */ char *p; /* generic pointer */ } t; /* TIFF scanline */ - void (*tf)(); /* translation procedure */ + colcvf_t *tf; /* translation procedure */ } cvts = { /* conversion structure */ 0, "", "", COMPRESSION_NONE, PHOTOMETRIC_RGB, PLANARCONFIG_CONTIG, GAMCOR, 0, 1, 1., 1., @@ -67,10 +70,20 @@ struct { #define CLR(f) (cvts.flags &= ~(f)) #define TGL(f) (cvts.flags ^= (f)) -void Luv2Color(), L2Color(), RGB2Colr(), Gry2Colr(); -void Color2Luv(), Color2L(), Colr2RGB(), Colr2Gry(); -void RRGGBB2Color(), GGry2Color(), Color2RRGGBB(), Color2GGry(); +static colcvf_t Luv2Color, L2Color, RGB2Colr, Gry2Colr; +static colcvf_t Color2Luv, Color2L, Colr2RGB, Colr2Gry; +static colcvf_t RRGGBB2Color, GGry2Color, Color2RRGGBB, Color2GGry; +static gethfunc headline; +static void quiterr(char *err); +static void allocbufs(void); +static void initfromtif(void); +static void tiff2ra(int ac, char *av[]); +static void initfromrad(void); +static void ra2tiff(int ac, char *av[]); + + + #define RfGfBf2Color Luv2Color #define Gryf2Color L2Color #define Color2Gryf Color2L @@ -95,9 +108,8 @@ char OWNSTR[] = "OWNER="; char *progname; -main(argc, argv) -int argc; -char *argv[]; +int +main(int argc, char *argv[]) { int reverse = 0; int i; @@ -161,13 +173,14 @@ doneopts: if (i != argc-2) goto userr; /* consistency checks */ - if (CHK(C_GRY)) + if (CHK(C_GRY)) { if (cvts.phot == PHOTOMETRIC_RGB) cvts.phot = PHOTOMETRIC_MINISBLACK; else { cvts.phot = PHOTOMETRIC_LOGL; cvts.comp = COMPRESSION_SGILOG; } + } if (CHK(C_TWRD|C_TFLT) == (C_TWRD|C_TFLT)) goto userr; @@ -177,17 +190,19 @@ doneopts: exit(0); userr: fprintf(stderr, - "Usage: %s [-z|-L|-l|-f|-w][-b][-e +/-stops][-g gamma] {in.pic|-} out.tif\n", + "Usage: %s [-z|-L|-l|-f|-w][-b][-e +/-stops][-g gamma] {in.hdr|-} out.tif\n", progname); fprintf(stderr, - " Or: %s -r [-x][-e +/-stops][-g gamma] in.tif [out.pic|-]\n", + " Or: %s -r [-x][-e +/-stops][-g gamma] in.tif [out.hdr|-]\n", progname); exit(1); } -quiterr(err) /* print message and exit */ -char *err; +static void +quiterr( /* print message and exit */ + char *err +) { if (err != NULL) { fprintf(stderr, "%s: %s\n", progname, err); @@ -197,7 +212,8 @@ char *err; } -allocbufs() /* allocate scanline buffers */ +static void +allocbufs(void) /* allocate scanline buffers */ { int rsiz, tsiz; @@ -207,12 +223,13 @@ allocbufs() /* allocate scanline buffers */ (CHK(C_GRY) ? 1 : 3); cvts.r.p = (char *)malloc(rsiz*cvts.xmax); cvts.t.p = (char *)malloc(tsiz*cvts.xmax); - if (cvts.r.p == NULL | cvts.t.p == NULL) + if ((cvts.r.p == NULL) | (cvts.t.p == NULL)) quiterr("no memory to allocate scanline buffers"); } -initfromtif() /* initialize conversion from TIFF input */ +static void +initfromtif(void) /* initialize conversion from TIFF input */ { uint16 hi; char *cp; @@ -257,7 +274,7 @@ initfromtif() /* initialize conversion from TIFF inpu cpcolormat(cvts.cmat, xyz2rgbmat); SET(C_CXFM|C_GAMUT); } else if (cvts.comp == COMPRESSION_SGILOG) - SET(C_GAMUT); + SET(C_GAMUT); /* may be outside XYZ gamut */ if (cvts.pconf != PLANARCONFIG_CONTIG) quiterr("cannot handle separate Luv planes"); TIFFSetField(cvts.tif, TIFFTAG_SGILOGDATAFMT, @@ -344,7 +361,7 @@ initfromtif() /* initialize conversion from TIFF inpu quiterr("unknown input image resolution"); if (!TIFFGetField(cvts.tif, TIFFTAG_STONITS, &cvts.stonits)) - cvts.stonits = 1.; + cvts.stonits = -1.; if (!TIFFGetField(cvts.tif, TIFFTAG_DATETIME, &cp)) cvts.capdate[0] = '\0'; @@ -362,13 +379,15 @@ initfromtif() /* initialize conversion from TIFF inpu if (cvts.pixrat < .99 || cvts.pixrat > 1.01) fputaspect(cvts.pixrat, cvts.rfp); if (CHK(C_XYZE)) { - fputexpos(pow(2.,(double)cvts.bradj)/cvts.stonits, cvts.rfp); + if (cvts.stonits > .0) + fputexpos(pow(2.,(double)cvts.bradj)/cvts.stonits, cvts.rfp); fputformat(CIEFMT, cvts.rfp); } else { if (CHK(C_PRIM)) fputprims(cvts.prims, cvts.rfp); - fputexpos(WHTEFFICACY*pow(2.,(double)cvts.bradj)/cvts.stonits, - cvts.rfp); + if (cvts.stonits > .0) + fputexpos(WHTEFFICACY*pow(2.,(double)cvts.bradj)/cvts.stonits, + cvts.rfp); fputformat(COLRFMT, cvts.rfp); } if (cvts.capdate[0]) @@ -380,9 +399,11 @@ initfromtif() /* initialize conversion from TIFF inpu } -tiff2ra(ac, av) /* convert TIFF image to Radiance picture */ -int ac; -char *av[]; +static void +tiff2ra( /* convert TIFF image to Radiance picture */ + int ac, + char *av[] +) { int32 y; /* open TIFF input */ @@ -393,6 +414,7 @@ char *av[]; cvts.rfp = stdout; else if ((cvts.rfp = fopen(av[ac+1], "w")) == NULL) quiterr("cannot open Radiance output picture"); + SET_FILE_BINARY(cvts.rfp); /* start output header */ newheader("RADIANCE", cvts.rfp); printargs(ac, av, cvts.rfp); @@ -410,13 +432,15 @@ char *av[]; } -int -headline(s) /* process Radiance input header line */ -char *s; +static int +headline( /* process Radiance input header line */ + char *s, + void *p +) { static int tmstrlen = 0; static int ownstrlen = 0; - char fmt[32]; + char fmt[MAXFMTLEN]; if (!tmstrlen) tmstrlen = strlen(TMSTR); @@ -469,7 +493,8 @@ char *s; } -initfromrad() /* initialize input from a Radiance picture */ +static void +initfromrad(void) /* initialize input from a Radiance picture */ { int i1, i2, po; /* read Radiance header */ @@ -546,6 +571,7 @@ initfromrad() /* initialize input from a Radiance pi SAMPLEFORMAT_IEEEFP); cvts.tf = Color2RfGfBf; SET(C_RFLT); + CLR(C_GAMUT); } else cvts.tf = Colr2RGB; break; @@ -596,9 +622,11 @@ initfromrad() /* initialize input from a Radiance pi } -ra2tiff(ac, av) /* convert Radiance picture to TIFF image */ -int ac; -char *av[]; +static void +ra2tiff( /* convert Radiance picture to TIFF image */ + int ac, + char *av[] +) { uint32 y; /* open Radiance file */ @@ -606,6 +634,7 @@ char *av[]; cvts.rfp = stdin; else if ((cvts.rfp = fopen(av[ac], "r")) == NULL) quiterr("cannot open Radiance input picture"); + SET_FILE_BINARY(cvts.rfp); /* open TIFF file */ if ((cvts.tif = TIFFOpen(av[ac+1], "w")) == NULL) quiterr("cannot open TIFF output"); @@ -620,15 +649,19 @@ char *av[]; } -void -Luv2Color(y) /* read/convert/write Luv->COLOR scanline */ -uint32 y; +static void +Luv2Color( /* read/convert/write Luv->COLOR scanline */ + uint32 y +) { register int x; if (CHK(C_RFLT|C_TWRD|C_TFLT|C_GRY) != (C_RFLT|C_TFLT)) quiterr("internal error 1 in Luv2Color"); + if (cvts.pconf != PLANARCONFIG_CONTIG) + quiterr("cannot handle separate 32-bit color planes"); + if (TIFFReadScanline(cvts.tif, cvts.t.p, y, 0) < 0) quiterr("error reading TIFF input"); /* also works for float RGB */ @@ -655,17 +688,21 @@ uint32 y; } -void -RRGGBB2Color(y) /* read/convert/write RGB16->COLOR scanline */ -uint32 y; +static void +RRGGBB2Color( /* read/convert/write RGB16->COLOR scanline */ + uint32 y +) { - int dogamma = cvts.gamcor < 0.99 | cvts.gamcor > 1.01; + int dogamma = (cvts.gamcor < 0.99) | (cvts.gamcor > 1.01); register double d; register int x; if (CHK(C_RFLT|C_TWRD|C_TFLT|C_GRY) != (C_TWRD|C_RFLT)) quiterr("internal error 1 in RRGGBB2Color"); + if (cvts.pconf != PLANARCONFIG_CONTIG) + quiterr("cannot handle separate 16-bit color planes"); + if (TIFFReadScanline(cvts.tif, cvts.t.p, y, 0) < 0) quiterr("error reading TIFF input"); @@ -682,9 +719,6 @@ uint32 y; if (CHK(C_CXFM)) colortrans(cvts.r.colors[x], cvts.cmat, cvts.r.colors[x]); - if (CHK(C_GAMUT)) - clipgamut(cvts.r.colors[x], cvts.t.fp[3*x + 1], - CGAMUT_LOWER, cblack, cwhite); } if (cvts.bradj) { d = pow(2.,(double)cvts.bradj); @@ -697,9 +731,10 @@ uint32 y; } -void -L2Color(y) /* read/convert/write Lfloat->COLOR scanline */ -uint32 y; +static void +L2Color( /* read/convert/write Lfloat->COLOR scanline */ + uint32 y +) { float m = pow(2., (double)cvts.bradj); register int x; @@ -720,9 +755,10 @@ uint32 y; } -void -RGB2Colr(y) /* read/convert/write RGB->COLR scanline */ -uint32 y; +static void +RGB2Colr( /* read/convert/write RGB->COLR scanline */ + uint32 y +) { COLOR ctmp; register int x; @@ -776,9 +812,10 @@ readerr: } -void -Gry2Colr(y) /* read/convert/write G8->COLR scanline */ -uint32 y; +static void +Gry2Colr( /* read/convert/write G8->COLR scanline */ + uint32 y +) { register int x; @@ -802,11 +839,12 @@ uint32 y; } -void -GGry2Color(y) /* read/convert/write G16->COLOR scanline */ -uint32 y; +static void +GGry2Color( /* read/convert/write G16->COLOR scanline */ + uint32 y +) { - int dogamma = cvts.gamcor < 0.99 | cvts.gamcor > 1.01; + int dogamma = (cvts.gamcor < 0.99) | (cvts.gamcor > 1.01); double m; register double d; register int x; @@ -832,11 +870,12 @@ uint32 y; } -void -Color2GGry(y) /* read/convert/write COLOR->G16 scanline */ -uint32 y; +static void +Color2GGry( /* read/convert/write COLOR->G16 scanline */ + uint32 y +) { - int dogamma = cvts.gamcor < 0.99 | cvts.gamcor > 1.01; + int dogamma = (cvts.gamcor < 0.99) | (cvts.gamcor > 1.01); float m = pow(2.,(double)cvts.bradj); register int x; @@ -866,9 +905,10 @@ uint32 y; } -void -Color2L(y) /* read/convert/write COLOR->Lfloat scanline */ -uint32 y; +static void +Color2L( /* read/convert/write COLOR->Lfloat scanline */ + uint32 y +) { float m = pow(2.,(double)cvts.bradj); register int x; @@ -888,9 +928,10 @@ uint32 y; } -void -Color2Luv(y) /* read/convert/write COLOR->Luv scanline */ -uint32 y; +static void +Color2Luv( /* read/convert/write COLOR->Luv scanline */ + uint32 y +) { register int x; @@ -921,11 +962,12 @@ uint32 y; } -void -Color2RRGGBB(y) /* read/convert/write COLOR->RGB16 scanline */ -uint32 y; +static void +Color2RRGGBB( /* read/convert/write COLOR->RGB16 scanline */ + uint32 y +) { - int dogamma = cvts.gamcor < 0.99 | cvts.gamcor > 1.01; + int dogamma = (cvts.gamcor < 0.99) | (cvts.gamcor > 1.01); float m = pow(2.,(double)cvts.bradj); register int x, i; @@ -935,7 +977,14 @@ uint32 y; if (freadscan(cvts.r.colors, cvts.xmax, cvts.rfp) < 0) quiterr("error reading Radiance picture"); - for (x = cvts.xmax; x--; ) + for (x = cvts.xmax; x--; ) { + if (CHK(C_CXFM)) { + colortrans(cvts.r.colors[x], cvts.cmat, + cvts.r.colors[x]); + if (CHK(C_GAMUT)) + clipgamut(cvts.r.colors[x], bright(cvts.r.colors[x]), + CGAMUT_LOWER, cblack, cwhite); + } for (i = 3; i--; ) { register float f = m*colval(cvts.r.colors[x],i); if (f <= 0) @@ -948,15 +997,17 @@ uint32 y; else cvts.t.wp[3*x + i] = (int)((float)(1L<<16)*f); } + } if (TIFFWriteScanline(cvts.tif, cvts.t.p, y, 0) < 0) quiterr("error writing TIFF output"); } -void -Colr2Gry(y) /* read/convert/write COLR->RGB scanline */ -uint32 y; +static void +Colr2Gry( /* read/convert/write COLR->RGB scanline */ + uint32 y +) { register int x; @@ -980,9 +1031,10 @@ uint32 y; } -void -Colr2RGB(y) /* read/convert/write COLR->RGB scanline */ -uint32 y; +static void +Colr2RGB( /* read/convert/write COLR->RGB scanline */ + uint32 y +) { COLOR ctmp; register int x;