src/px/slide.c

#ifndef lint
static const char       RCSid[] = "$Id$";
#endif
/*
 *  slide.c - driver for dicomed film recorder using GRAFPAC under VMS.
 *
 *     2/19/86
 *
 *      Version for color slide ouput.
 *      Link with grafpac driver "DS":
 *              link slide.obj,color.obj,sys_grafpac:ds/lib
 */

#include  <stdio.h>

#include  "color.h"

#include  "random.h"


#define  XSIZE          4096            /* device width */
#define  YSIZE          4096            /* device length */

#define  NCOLS          4018            /* maximum # columns for output */
#define  NROWS          2712            /* maximum # rows for output */

#define  XCENTER        2048            /* center for x */
#define  YCENTER        2048            /* center for y */

int  nrows;                             /* number of rows for output */
int  ncols;                             /* number of columns for output */

int  xoffset;                           /* offset for x */
int  yoffset;                           /* offset for y */

#define  RSIZ           17              /* size of response map */

float  rmap[2*RSIZ] = {                         /* response map */
        .0,             .0,
        .0004,          .059,
        .0008,          .122,
        .0018,          .184,
        .0033,          .247,
        .0061,          .310,
        .0119,          .373,
        .0202,          .435,
        .0373,          .498,
        .065,           .561,
        .113,           .624,
        .199,           .686,
        .300,           .749,
        .432,           .812,
        .596,           .875,
        .767,           .937,
        1.000,          1.000,
};


main(argc, argv)
int  argc;
char  **argv;
{
        FILE  *fin;
        COLOR  scanline[NCOLS];
        int  mult;
        char  sbuf[128];
        int  i, j;

/*
        for (i = 1; i < argc; i++)
                if (argv[i][0] == '-')
                        switch (argv[i][1]) {
                        default:
                                fprintf(stderr, "unknown option: %s\n", argv[i]);
                                exit(1);
                                break;
                        }
                else
                        break;
                        
*/
        argc = 2;
        argv[1] = "pixfile";
        i = 1;

        if (argc - i < 1)
                fin = 0;
        else if ((fin = fopen(argv[i], "r")) == NULL) {
                        fprintf(stderr, "%s: can't open file \"%s\"\n",
                                                argv[0], argv[i]);
                        exit(1);
                }
        
                                /* discard header */
        while (fgets(sbuf, sizeof(sbuf), fin) != NULL && sbuf[0] != '\n')
                ;
                                /* get picture size */
        if (fgets(sbuf, sizeof(sbuf), fin) == NULL ||
                        sscanf(sbuf, "-Y %d +X %d\n", &nrows, &ncols) != 2) {
                fprintf(stderr, "%s: bad picture size\n", argv[0]);
                exit(1);
        }

        if (ncols > NCOLS || nrows > NROWS) {
                fprintf(stderr, "%s: resolution mismatch\n", argv[0]);
                exit(1);
        }

        mult = NCOLS/ncols;
        if (NROWS/nrows < mult)
                mult = NROWS/nrows;

        xoffset = XCENTER/mult - ncols/2;
        yoffset = YCENTER/mult - nrows/2;

        tvinit();
        i = XSIZE/mult; j = YSIZE/mult; tvvrs(&i, &j);
        i = 12; tvset("ifile", &i);
        tvset("nrco", "rgb");

        for (i = nrows-1; i >= 0; i--) {
                if (freadscan(scanline, ncols, fin) < 0) {
                        fprintf(stderr, "%s: read error in row %d\n",
                                        argv[0], i);
                        exit(1);
                }
                scanout(scanline, i);
        }

        tvsend();
        tvend();

        return(0);
}


scanout(scan, y)                        /* output scan line */
register COLOR  *scan;
int  y;
{
        static int  zero;
        static float  rgb[3][NCOLS];
        int  i;
        register int  j;
        
        zero = xoffset;
        y += yoffset;

        for (j = 0; j < ncols; j++) {
                mapcolor(scan[j], scan[j]);
                for (i = 0; i < 3; i++)
                        rgb[i][j] = scan[j][i] + (.5-frandom())/255.0;
        }
        tvrstc(&zero, &y, rgb[0], rgb[1], rgb[2], &ncols);
}


mapcolor(c1, c2)                /* map c1 into c2 */
register COLOR  c1;
COLOR  c2;
{
        static float  *mp[3] = {rmap, rmap, rmap};
        double  x;
        register int  i;

        for (i = 0; i < 3; i++) {
                if (colval(c1,i) >= 1.0) {
                        colval(c2,i) = 1.0;
                        continue;
                }
                while (colval(c1,i) >= mp[i][2])
                        mp[i] += 2;
                while (colval(c1,i) < mp[i][0])
                        mp[i] -= 2;
                x = (colval(c1,i) - mp[i][0]) / (mp[i][2] - mp[i][0]);
                colval(c2,i) = (1.0-x)*mp[i][1] + x*mp[i][3];
        }
}
Revision:	1.3
Committed:	Sat Feb 22 02:07:28 2003 UTC (21 years, 2 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R6P1, rad3R5, rad3R6
Changes since 1.2:	+1 -5 lines
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	Content
1	#ifndef lint
2	static const char RCSid[] = "$Id$";
3	#endif
4	/*
5	* slide.c - driver for dicomed film recorder using GRAFPAC under VMS.
6	*
7	* 2/19/86
8	*
9	* Version for color slide ouput.
10	* Link with grafpac driver "DS":
11	* link slide.obj,color.obj,sys_grafpac:ds/lib
12	*/
13
14	#include <stdio.h>
15
16	#include "color.h"
17
18	#include "random.h"
19
20
21	#define XSIZE 4096 /* device width */
22	#define YSIZE 4096 /* device length */
23
24	#define NCOLS 4018 /* maximum # columns for output */
25	#define NROWS 2712 /* maximum # rows for output */
26
27	#define XCENTER 2048 /* center for x */
28	#define YCENTER 2048 /* center for y */
29
30	int nrows; /* number of rows for output */
31	int ncols; /* number of columns for output */
32
33	int xoffset; /* offset for x */
34	int yoffset; /* offset for y */
35
36	#define RSIZ 17 /* size of response map */
37
38	float rmap[2RSIZ] = { / response map */
39	.0, .0,
40	.0004, .059,
41	.0008, .122,
42	.0018, .184,
43	.0033, .247,
44	.0061, .310,
45	.0119, .373,
46	.0202, .435,
47	.0373, .498,
48	.065, .561,
49	.113, .624,
50	.199, .686,
51	.300, .749,
52	.432, .812,
53	.596, .875,
54	.767, .937,
55	1.000, 1.000,
56	};
57
58
59	main(argc, argv)
60	int argc;
61	char **argv;
62	{
63	FILE *fin;
64	COLOR scanline[NCOLS];
65	int mult;
66	char sbuf[128];
67	int i, j;
68
69	/*
70	for (i = 1; i < argc; i++)
71	if (argv[i][0] == '-')
72	switch (argv[i][1]) {
73	default:
74	fprintf(stderr, "unknown option: %s\n", argv[i]);
75	exit(1);
76	break;
77	}
78	else
79	break;
80
81	*/
82	argc = 2;
83	argv[1] = "pixfile";
84	i = 1;
85
86	if (argc - i < 1)
87	fin = 0;
88	else if ((fin = fopen(argv[i], "r")) == NULL) {
89	fprintf(stderr, "%s: can't open file \"%s\"\n",
90	argv[0], argv[i]);
91	exit(1);
92	}
93
94	/* discard header */
95	while (fgets(sbuf, sizeof(sbuf), fin) != NULL && sbuf[0] != '\n')
96	;
97	/* get picture size */
98	if (fgets(sbuf, sizeof(sbuf), fin) == NULL \|\|
99	sscanf(sbuf, "-Y %d +X %d\n", &nrows, &ncols) != 2) {
100	fprintf(stderr, "%s: bad picture size\n", argv[0]);
101	exit(1);
102	}
103
104	if (ncols > NCOLS \|\| nrows > NROWS) {
105	fprintf(stderr, "%s: resolution mismatch\n", argv[0]);
106	exit(1);
107	}
108
109	mult = NCOLS/ncols;
110	if (NROWS/nrows < mult)
111	mult = NROWS/nrows;
112
113	xoffset = XCENTER/mult - ncols/2;
114	yoffset = YCENTER/mult - nrows/2;
115
116	tvinit();
117	i = XSIZE/mult; j = YSIZE/mult; tvvrs(&i, &j);
118	i = 12; tvset("ifile", &i);
119	tvset("nrco", "rgb");
120
121	for (i = nrows-1; i >= 0; i--) {
122	if (freadscan(scanline, ncols, fin) < 0) {
123	fprintf(stderr, "%s: read error in row %d\n",
124	argv[0], i);
125	exit(1);
126	}
127	scanout(scanline, i);
128	}
129
130	tvsend();
131	tvend();
132
133	return(0);
134	}
135
136
137	scanout(scan, y) /* output scan line */
138	register COLOR *scan;
139	int y;
140	{
141	static int zero;
142	static float rgb[3][NCOLS];
143	int i;
144	register int j;
145
146	zero = xoffset;
147	y += yoffset;
148
149	for (j = 0; j < ncols; j++) {
150	mapcolor(scan[j], scan[j]);
151	for (i = 0; i < 3; i++)
152	rgb[i][j] = scan[j][i] + (.5-frandom())/255.0;
153	}
154	tvrstc(&zero, &y, rgb[0], rgb[1], rgb[2], &ncols);
155	}
156
157
158	mapcolor(c1, c2) /* map c1 into c2 */
159	register COLOR c1;
160	COLOR c2;
161	{
162	static float *mp[3] = {rmap, rmap, rmap};
163	double x;
164	register int i;
165
166	for (i = 0; i < 3; i++) {
167	if (colval(c1,i) >= 1.0) {
168	colval(c2,i) = 1.0;
169	continue;
170	}
171	while (colval(c1,i) >= mp[i][2])
172	mp[i] += 2;
173	while (colval(c1,i) < mp[i][0])
174	mp[i] -= 2;
175	x = (colval(c1,i) - mp[i][0]) / (mp[i][2] - mp[i][0]);
176	colval(c2,i) = (1.0-x)mp[i][1] + xmp[i][3];
177	}
178	}