src/px/pcwarp.c

#ifndef lint
static const char       RCSid[] = "$Id: pcwarp.c,v 3.6 2018/08/02 18:33:44 greg Exp $";
#endif
/*
 * Warp colors in Radiance picture to correct for input/output changes.
 */

#include "rtio.h"
#include "resolu.h"
#include "color.h"
#include "warp3d.h"

char    *progname;                      /* global argv[0] */

FILE    *infp = NULL;                   /* input stream */
int     xres, yres;                     /* input picture resolution */

WARP3D  *cwarp;                         /* our warp map */
int     iclip = CGAMUT_UPPER;           /* input value gamut clipping */
int     oclip = CGAMUT_LOWER;           /* output value gamut clipping */

static void syserror(char *s);
static void picwarp(void);


int
main(
        int     argc,
        char    *argv[]
)
{
        static char     picfmt[MAXFMTLEN] = PICFMT;
        int     cwflags = 0;
        int     rval;
        int     i;

        progname = argv[0];
        infp = stdin;
                                        /* get options */
        for (i = 1; i < argc && argv[i][0] == '-'; i++)
                switch (argv[i][1]) {
                case 'e':
                        cwflags = W3EXACT;
                        break;
                case 'f':
                        cwflags = W3FAST;
                        break;
                case 'i':
                        iclip = 0;
                        break;
                case 'o':
                        oclip = CGAMUT;
                        break;
                default:
                        goto userr;
                }
                                        /* load warp map */
        if (i >= argc)
                goto userr;
        if ((cwarp = load3dw(argv[i++], NULL)) == NULL)
                syserror("load3dw");
        set3dwfl(cwarp, cwflags);
                                        /* open input and output pictures */
        if (i < argc && (infp = fopen(argv[i], "r")) == NULL)
                syserror(argv[i]);
        if (i < argc-1 && freopen(argv[i+1], "w", stdout) == NULL)
                syserror(argv[i+1]);
                                        /* transfer header */
        if ((rval = checkheader(infp, picfmt, stdout)) < 0) {
                fprintf(stderr, "%s: input not a Radiance picture\n",
                                progname);
                exit(1);
        }
        if (rval)
                fputformat(picfmt, stdout);
                                        /* add new header info. */
        printargs(i, argv, stdout);
        putchar('\n');
                                        /* get picture size */
        if ((rval = fgetresolu(&xres, &yres, infp)) < 0) {
                fprintf(stderr, "%s: bad picture size\n", progname);
                exit(1);
        }
                                        /* new picture size the same */
        fputresolu(rval, xres, yres, stdout);
                                        /* warp those colors! */
        picwarp();
        exit(0);
userr:
        fprintf(stderr,
                "Usage: %s [-i][-o][-e|-f] map.cwp [input.hdr [output.hdr]]\n",
                        progname);
        exit(1);
}


static void
syserror(                       /* print system error and exit */
        char    *s
)
{
        fprintf(stderr, "%s: ", progname);
        perror(s);
        exit(2);
}


static void
picwarp(void)                   /* warp our picture scanlines */
{
        register COLOR  *scan;
        long    ngamut = 0;
        int     rval;
        int     y;
        register int    x;

        scan = (COLOR *)malloc(xres*sizeof(COLOR));
        if (scan == NULL)
                syserror("picwarp");
        for (y = 0; y < yres; y++) {
                if (freadscan(scan, xres, infp) < 0) {
                        fprintf(stderr, "%s: error reading input picture\n",
                                        progname);
                        exit(1);
                }
                for (x = 0; x < xres; x++) {
                        if (iclip)
                                clipgamut(scan[x], bright(scan[x]), iclip,
                                                cblack, cwhite);
                        rval = warp3d(scan[x], scan[x], cwarp);
                        if (rval & W3ERROR)
                                syserror("warp3d");
                        if (rval & W3BADMAP) {
                                fprintf(stderr, "%s: singular color mapping\n",
                                                progname);
                                exit(1);
                        }
                        if (rval & W3GAMUT)
                                ngamut++;
                        if (oclip)
                                clipgamut(scan[x], bright(scan[x]), oclip,
                                                cblack, cwhite);
                }
                if (fwritescan(scan, xres, stdout) < 0) {
                        fprintf(stderr, "%s: error writing output picture\n",
                                        progname);
                        exit(1);
                }
        }
        if (ngamut >= (long)xres*yres/100)
                fprintf(stderr, "%s: warning - %ld%% of pixels out of gamut\n",
                                progname, 100*ngamut/((long)xres*yres));
        free((void *)scan);
}
Revision:	3.7
Committed:	Fri Jul 19 17:37:56 2019 UTC (4 years, 9 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad5R4, rad5R3, HEAD
Changes since 3.6:	+2 -3 lines
Log Message:	Moved declarations and definitions for header.c from resolu.h to rtio.h
#	User	Rev	Content
1	greg	3.1	#ifndef lint
2	greg	3.7	static const char RCSid[] = "$Id: pcwarp.c,v 3.6 2018/08/02 18:33:44 greg Exp $";
3	greg	3.1	#endif
4			/*
5			* Warp colors in Radiance picture to correct for input/output changes.
6			*/
7
8	greg	3.7	#include "rtio.h"
9	schorsch	3.4	#include "resolu.h"
10	greg	3.1	#include "color.h"
11			#include "warp3d.h"
12
13			char progname; / global argv[0] */
14
15	gwlarson	3.2	FILE infp = NULL; / input stream */
16	greg	3.1	int xres, yres; /* input picture resolution */
17
18			WARP3D cwarp; / our warp map */
19			int iclip = CGAMUT_UPPER; /* input value gamut clipping */
20			int oclip = CGAMUT_LOWER; /* output value gamut clipping */
21
22	schorsch	3.4	static void syserror(char *s);
23			static void picwarp(void);
24
25	greg	3.1
26	schorsch	3.4	int
27			main(
28			int argc,
29			char *argv[]
30			)
31	greg	3.1	{
32	greg	3.6	static char picfmt[MAXFMTLEN] = PICFMT;
33	greg	3.1	int cwflags = 0;
34			int rval;
35			int i;
36
37			progname = argv[0];
38	gwlarson	3.2	infp = stdin;
39	greg	3.1	/* get options */
40			for (i = 1; i < argc && argv[i][0] == '-'; i++)
41			switch (argv[i][1]) {
42			case 'e':
43			cwflags = W3EXACT;
44			break;
45			case 'f':
46			cwflags = W3FAST;
47			break;
48			case 'i':
49			iclip = 0;
50			break;
51			case 'o':
52			oclip = CGAMUT;
53			break;
54			default:
55			goto userr;
56			}
57			/* load warp map */
58			if (i >= argc)
59			goto userr;
60			if ((cwarp = load3dw(argv[i++], NULL)) == NULL)
61			syserror("load3dw");
62			set3dwfl(cwarp, cwflags);
63			/* open input and output pictures */
64			if (i < argc && (infp = fopen(argv[i], "r")) == NULL)
65			syserror(argv[i]);
66			if (i < argc-1 && freopen(argv[i+1], "w", stdout) == NULL)
67			syserror(argv[i+1]);
68			/* transfer header */
69			if ((rval = checkheader(infp, picfmt, stdout)) < 0) {
70			fprintf(stderr, "%s: input not a Radiance picture\n",
71			progname);
72			exit(1);
73			}
74			if (rval)
75			fputformat(picfmt, stdout);
76			/* add new header info. */
77			printargs(i, argv, stdout);
78			putchar('\n');
79			/* get picture size */
80			if ((rval = fgetresolu(&xres, &yres, infp)) < 0) {
81			fprintf(stderr, "%s: bad picture size\n", progname);
82			exit(1);
83			}
84			/* new picture size the same */
85			fputresolu(rval, xres, yres, stdout);
86			/* warp those colors! */
87			picwarp();
88			exit(0);
89			userr:
90			fprintf(stderr,
91	greg	3.5	"Usage: %s [-i][-o][-e\|-f] map.cwp [input.hdr [output.hdr]]\n",
92	greg	3.1	progname);
93			exit(1);
94			}
95
96
97	schorsch	3.4	static void
98			syserror( /* print system error and exit */
99			char *s
100			)
101	greg	3.1	{
102			fprintf(stderr, "%s: ", progname);
103			perror(s);
104			exit(2);
105			}
106
107
108	schorsch	3.4	static void
109			picwarp(void) /* warp our picture scanlines */
110	greg	3.1	{
111			register COLOR *scan;
112			long ngamut = 0;
113			int rval;
114			int y;
115			register int x;
116
117			scan = (COLOR )malloc(xressizeof(COLOR));
118			if (scan == NULL)
119			syserror("picwarp");
120			for (y = 0; y < yres; y++) {
121			if (freadscan(scan, xres, infp) < 0) {
122			fprintf(stderr, "%s: error reading input picture\n",
123			progname);
124			exit(1);
125			}
126			for (x = 0; x < xres; x++) {
127			if (iclip)
128			clipgamut(scan[x], bright(scan[x]), iclip,
129			cblack, cwhite);
130			rval = warp3d(scan[x], scan[x], cwarp);
131			if (rval & W3ERROR)
132			syserror("warp3d");
133			if (rval & W3BADMAP) {
134			fprintf(stderr, "%s: singular color mapping\n",
135			progname);
136			exit(1);
137			}
138			if (rval & W3GAMUT)
139			ngamut++;
140			if (oclip)
141			clipgamut(scan[x], bright(scan[x]), oclip,
142			cblack, cwhite);
143			}
144			if (fwritescan(scan, xres, stdout) < 0) {
145			fprintf(stderr, "%s: error writing output picture\n",
146			progname);
147			exit(1);
148			}
149			}
150			if (ngamut >= (long)xres*yres/100)
151	schorsch	3.4	fprintf(stderr, "%s: warning - %ld%% of pixels out of gamut\n",
152	greg	3.1	progname, 100ngamut/((long)xresyres));
153	greg	3.3	free((void *)scan);
154	greg	3.1	}