src/px/pcwarp.c

#ifndef lint
static const char       RCSid[] = "$Id: pcwarp.c,v 3.7 2019/07/19 17:37:56 greg Exp $";
#endif
/*
 * Warp colors in Radiance picture to correct for input/output changes.
 */

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

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;

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