[wxWidgets.git] / src / tiff / tools / rgb2ycbcr.c

/* $Id$ */

/*
 * Copyright (c) 1991-1997 Sam Leffler
 * Copyright (c) 1991-1997 Silicon Graphics, Inc.
 *
 * Permission to use, copy, modify, distribute, and sell this software and 
 * its documentation for any purpose is hereby granted without fee, provided
 * that (i) the above copyright notices and this permission notice appear in
 * all copies of the software and related documentation, and (ii) the names of
 * Sam Leffler and Silicon Graphics may not be used in any advertising or
 * publicity relating to the software without the specific, prior written
 * permission of Sam Leffler and Silicon Graphics.
 * 
 * THE SOFTWARE IS PROVIDED "AS-IS" AND WITHOUT WARRANTY OF ANY KIND, 
 * EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY 
 * WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE.  
 * 
 * IN NO EVENT SHALL SAM LEFFLER OR SILICON GRAPHICS BE LIABLE FOR
 * ANY SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY KIND,
 * OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS,
 * WHETHER OR NOT ADVISED OF THE POSSIBILITY OF DAMAGE, AND ON ANY THEORY OF 
 * LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE 
 * OF THIS SOFTWARE.
 */

#include "tif_config.h"

#include <stdio.h>
#include <string.h>
#include <stdlib.h>

#ifdef HAVE_UNISTD_H
# include <unistd.h>
#endif

#include "tiffio.h"

#define	streq(a,b)	(strcmp(a,b) == 0)
#define	CopyField(tag, v) \
    if (TIFFGetField(in, tag, &v)) TIFFSetField(out, tag, v)

#ifndef howmany
#define	howmany(x, y)	(((x)+((y)-1))/(y))
#endif
#define	roundup(x, y)	(howmany(x,y)*((uint32)(y)))

#define	LumaRed		ycbcrCoeffs[0]
#define	LumaGreen	ycbcrCoeffs[1]
#define	LumaBlue	ycbcrCoeffs[2]

uint16	compression = COMPRESSION_PACKBITS;
uint32	rowsperstrip = (uint32) -1;

uint16	horizSubSampling = 2;		/* YCbCr horizontal subsampling */
uint16	vertSubSampling = 2;		/* YCbCr vertical subsampling */
float	ycbcrCoeffs[3] = { .299F, .587F, .114F };
/* default coding range is CCIR Rec 601-1 with no headroom/footroom */
float	refBlackWhite[6] = { 0.F, 255.F, 128.F, 255.F, 128.F, 255.F };

static	int tiffcvt(TIFF* in, TIFF* out);
static	void usage(int code);
static	void setupLumaTables(void);

int
main(int argc, char* argv[])
{
	TIFF *in, *out;
	int c;
	extern int optind;
	extern char *optarg;

	while ((c = getopt(argc, argv, "c:h:r:v:z")) != -1)
		switch (c) {
		case 'c':
			if (streq(optarg, "none"))
			    compression = COMPRESSION_NONE;
			else if (streq(optarg, "packbits"))
			    compression = COMPRESSION_PACKBITS;
			else if (streq(optarg, "lzw"))
			    compression = COMPRESSION_LZW;
			else if (streq(optarg, "jpeg"))
			    compression = COMPRESSION_JPEG;
			else if (streq(optarg, "zip"))
			    compression = COMPRESSION_ADOBE_DEFLATE;
			else
			    usage(-1);
			break;
		case 'h':
			horizSubSampling = atoi(optarg);
			break;
		case 'v':
			vertSubSampling = atoi(optarg);
			break;
		case 'r':
			rowsperstrip = atoi(optarg);
			break;
		case 'z':	/* CCIR Rec 601-1 w/ headroom/footroom */
			refBlackWhite[0] = 16.;
			refBlackWhite[1] = 235.;
			refBlackWhite[2] = 128.;
			refBlackWhite[3] = 240.;
			refBlackWhite[4] = 128.;
			refBlackWhite[5] = 240.;
			break;
		case '?':
			usage(0);
			/*NOTREACHED*/
		}
	if (argc - optind < 2)
		usage(-1);
	out = TIFFOpen(argv[argc-1], "w");
	if (out == NULL)
		return (-2);
	setupLumaTables();
	for (; optind < argc-1; optind++) {
		in = TIFFOpen(argv[optind], "r");
		if (in != NULL) {
			do {
				if (!tiffcvt(in, out) ||
				    !TIFFWriteDirectory(out)) {
					(void) TIFFClose(out);
					return (1);
				}
			} while (TIFFReadDirectory(in));
			(void) TIFFClose(in);
		}
	}
	(void) TIFFClose(out);
	return (0);
}

float	*lumaRed;
float	*lumaGreen;
float	*lumaBlue;
float	D1, D2;
int	Yzero;

static float*
setupLuma(float c)
{
	float *v = (float *)_TIFFmalloc(256 * sizeof (float));
	int i;
	for (i = 0; i < 256; i++)
		v[i] = c * i;
	return (v);
}

static unsigned
V2Code(float f, float RB, float RW, int CR)
{
	unsigned int c = (unsigned int)((((f)*(RW-RB)/CR)+RB)+.5);
	return (c > 255 ? 255 : c);
}

static void
setupLumaTables(void)
{
	lumaRed = setupLuma(LumaRed);
	lumaGreen = setupLuma(LumaGreen);
	lumaBlue = setupLuma(LumaBlue);
	D1 = 1.F/(2.F - 2.F*LumaBlue);
	D2 = 1.F/(2.F - 2.F*LumaRed);
	Yzero = V2Code(0, refBlackWhite[0], refBlackWhite[1], 255);
}

static void
cvtClump(unsigned char* op, uint32* raster, uint32 ch, uint32 cw, uint32 w)
{
	float Y, Cb = 0, Cr = 0;
	uint32 j, k;
	/*
	 * Convert ch-by-cw block of RGB
	 * to YCbCr and sample accordingly.
	 */
	for (k = 0; k < ch; k++) {
		for (j = 0; j < cw; j++) {
			uint32 RGB = (raster - k*w)[j];
			Y = lumaRed[TIFFGetR(RGB)] +
			    lumaGreen[TIFFGetG(RGB)] +
			    lumaBlue[TIFFGetB(RGB)];
			/* accumulate chrominance */
			Cb += (TIFFGetB(RGB) - Y) * D1;
			Cr += (TIFFGetR(RGB) - Y) * D2;
			/* emit luminence */
			*op++ = V2Code(Y,
			    refBlackWhite[0], refBlackWhite[1], 255);
		}
		for (; j < horizSubSampling; j++)
			*op++ = Yzero;
	}
	for (; k < vertSubSampling; k++) {
		for (j = 0; j < horizSubSampling; j++)
			*op++ = Yzero;
	}
	/* emit sampled chrominance values */
	*op++ = V2Code(Cb / (ch*cw), refBlackWhite[2], refBlackWhite[3], 127);
	*op++ = V2Code(Cr / (ch*cw), refBlackWhite[4], refBlackWhite[5], 127);
}
#undef LumaRed
#undef LumaGreen
#undef LumaBlue
#undef V2Code

/*
 * Convert a strip of RGB data to YCbCr and
 * sample to generate the output data.
 */
static void
cvtStrip(unsigned char* op, uint32* raster, uint32 nrows, uint32 width)
{
	uint32 x;
	int clumpSize = vertSubSampling * horizSubSampling + 2;
	uint32 *tp;

	for (; nrows >= vertSubSampling; nrows -= vertSubSampling) {
		tp = raster;
		for (x = width; x >= horizSubSampling; x -= horizSubSampling) {
			cvtClump(op, tp,
			    vertSubSampling, horizSubSampling, width);
			op += clumpSize;
			tp += horizSubSampling;
		}
		if (x > 0) {
			cvtClump(op, tp, vertSubSampling, x, width);
			op += clumpSize;
		}
		raster -= vertSubSampling*width;
	}
	if (nrows > 0) {
		tp = raster;
		for (x = width; x >= horizSubSampling; x -= horizSubSampling) {
			cvtClump(op, tp, nrows, horizSubSampling, width);
			op += clumpSize;
			tp += horizSubSampling;
		}
		if (x > 0)
			cvtClump(op, tp, nrows, x, width);
	}
}

static int
cvtRaster(TIFF* tif, uint32* raster, uint32 width, uint32 height)
{
	uint32 y;
	tstrip_t strip = 0;
	tsize_t cc, acc;
	unsigned char* buf;
	uint32 rwidth = roundup(width, horizSubSampling);
	uint32 rheight = roundup(height, vertSubSampling);
	uint32 nrows = (rowsperstrip > rheight ? rheight : rowsperstrip);
        uint32 rnrows = roundup(nrows,vertSubSampling);

	cc = rnrows*rwidth +
	    2*((rnrows*rwidth) / (horizSubSampling*vertSubSampling));
	buf = (unsigned char*)_TIFFmalloc(cc);
	for (y = height; (int32) y > 0; y -= nrows) {
		uint32 nr = (y > nrows ? nrows : y);
		cvtStrip(buf, raster + (y-1)*width, nr, width);
		nr = roundup(nr, vertSubSampling);
		acc = nr*rwidth +
			2*((nr*rwidth)/(horizSubSampling*vertSubSampling));
		if (!TIFFWriteEncodedStrip(tif, strip++, buf, acc)) {
			_TIFFfree(buf);
			return (0);
		}
	}
	_TIFFfree(buf);
	return (1);
}

static int
tiffcvt(TIFF* in, TIFF* out)
{
	uint32 width, height;		/* image width & height */
	uint32* raster;			/* retrieve RGBA image */
	uint16 shortv;
	float floatv;
	char *stringv;
	uint32 longv;

	TIFFGetField(in, TIFFTAG_IMAGEWIDTH, &width);
	TIFFGetField(in, TIFFTAG_IMAGELENGTH, &height);
	raster = (uint32*)_TIFFmalloc(width * height * sizeof (uint32));
	if (raster == 0) {
		TIFFError(TIFFFileName(in), "No space for raster buffer");
		return (0);
	}
	if (!TIFFReadRGBAImage(in, width, height, raster, 0)) {
		_TIFFfree(raster);
		return (0);
	}

	CopyField(TIFFTAG_SUBFILETYPE, longv);
	TIFFSetField(out, TIFFTAG_IMAGEWIDTH, width);
	TIFFSetField(out, TIFFTAG_IMAGELENGTH, height);
	TIFFSetField(out, TIFFTAG_BITSPERSAMPLE, 8);
	TIFFSetField(out, TIFFTAG_COMPRESSION, compression);
	TIFFSetField(out, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_YCBCR);
	if (compression == COMPRESSION_JPEG)
		TIFFSetField(out, TIFFTAG_JPEGCOLORMODE, JPEGCOLORMODE_RAW);
	CopyField(TIFFTAG_FILLORDER, shortv);
	TIFFSetField(out, TIFFTAG_ORIENTATION, ORIENTATION_TOPLEFT);
	TIFFSetField(out, TIFFTAG_SAMPLESPERPIXEL, 3);
	CopyField(TIFFTAG_XRESOLUTION, floatv);
	CopyField(TIFFTAG_YRESOLUTION, floatv);
	CopyField(TIFFTAG_RESOLUTIONUNIT, shortv);
	TIFFSetField(out, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG);
	{ char buf[2048];
	  char *cp = strrchr(TIFFFileName(in), '/');
	  sprintf(buf, "YCbCr conversion of %s", cp ? cp+1 : TIFFFileName(in));
	  TIFFSetField(out, TIFFTAG_IMAGEDESCRIPTION, buf);
	}
	TIFFSetField(out, TIFFTAG_SOFTWARE, TIFFGetVersion());
	CopyField(TIFFTAG_DOCUMENTNAME, stringv);

	TIFFSetField(out, TIFFTAG_REFERENCEBLACKWHITE, refBlackWhite);
	TIFFSetField(out, TIFFTAG_YCBCRSUBSAMPLING,
	    horizSubSampling, vertSubSampling);
	TIFFSetField(out, TIFFTAG_YCBCRPOSITIONING, YCBCRPOSITION_CENTERED);
	TIFFSetField(out, TIFFTAG_YCBCRCOEFFICIENTS, ycbcrCoeffs);
	rowsperstrip = TIFFDefaultStripSize(out, rowsperstrip);
	TIFFSetField(out, TIFFTAG_ROWSPERSTRIP, rowsperstrip);

	return (cvtRaster(out, raster, width, height));
}

char* stuff[] = {
    "usage: rgb2ycbcr [-c comp] [-r rows] [-h N] [-v N] input... output\n",
    "where comp is one of the following compression algorithms:\n",
    " jpeg\t\tJPEG encoding\n",
    " lzw\t\tLempel-Ziv & Welch encoding\n",
    " zip\t\tdeflate encoding\n",
    " packbits\tPackBits encoding (default)\n",
    " none\t\tno compression\n",
    "and the other options are:\n",
    " -r\trows/strip\n",
    " -h\thorizontal sampling factor (1,2,4)\n",
    " -v\tvertical sampling factor (1,2,4)\n",
    NULL
};

static void
usage(int code)
{
	char buf[BUFSIZ];
	int i;

	setbuf(stderr, buf);
       
 fprintf(stderr, "%s\n\n", TIFFGetVersion());
	for (i = 0; stuff[i] != NULL; i++)
		fprintf(stderr, "%s\n", stuff[i]);
	exit(code);
}

/* vim: set ts=8 sts=8 sw=8 noet: */
Commit	Line	Data
8414a40c VZ	1	/* $Id$ */
	2
	3	/*
	4	* Copyright (c) 1991-1997 Sam Leffler
	5	* Copyright (c) 1991-1997 Silicon Graphics, Inc.
	6	*
	7	* Permission to use, copy, modify, distribute, and sell this software and
	8	* its documentation for any purpose is hereby granted without fee, provided
	9	* that (i) the above copyright notices and this permission notice appear in
	10	* all copies of the software and related documentation, and (ii) the names of
	11	* Sam Leffler and Silicon Graphics may not be used in any advertising or
	12	* publicity relating to the software without the specific, prior written
	13	* permission of Sam Leffler and Silicon Graphics.
	14	*
	15	* THE SOFTWARE IS PROVIDED "AS-IS" AND WITHOUT WARRANTY OF ANY KIND,
	16	* EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY
	17	* WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE.
	18	*
	19	* IN NO EVENT SHALL SAM LEFFLER OR SILICON GRAPHICS BE LIABLE FOR
	20	* ANY SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY KIND,
	21	* OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS,
	22	* WHETHER OR NOT ADVISED OF THE POSSIBILITY OF DAMAGE, AND ON ANY THEORY OF
	23	* LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE
	24	* OF THIS SOFTWARE.
	25	*/
	26
	27	#include "tif_config.h"
	28
	29	#include <stdio.h>
	30	#include <string.h>
	31	#include <stdlib.h>
	32
	33	#ifdef HAVE_UNISTD_H
	34	# include <unistd.h>
	35	#endif
	36
	37	#include "tiffio.h"
	38
	39	#define streq(a,b) (strcmp(a,b) == 0)
	40	#define CopyField(tag, v) \
	41	if (TIFFGetField(in, tag, &v)) TIFFSetField(out, tag, v)
	42
	43	#ifndef howmany
	44	#define howmany(x, y) (((x)+((y)-1))/(y))
	45	#endif
	46	#define roundup(x, y) (howmany(x,y)*((uint32)(y)))
	47
	48	#define LumaRed ycbcrCoeffs[0]
	49	#define LumaGreen ycbcrCoeffs[1]
	50	#define LumaBlue ycbcrCoeffs[2]
	51
	52	uint16 compression = COMPRESSION_PACKBITS;
	53	uint32 rowsperstrip = (uint32) -1;
	54
	55	uint16 horizSubSampling = 2; /* YCbCr horizontal subsampling */
	56	uint16 vertSubSampling = 2; /* YCbCr vertical subsampling */
	57	float ycbcrCoeffs[3] = { .299F, .587F, .114F };
	58	/* default coding range is CCIR Rec 601-1 with no headroom/footroom */
	59	float refBlackWhite[6] = { 0.F, 255.F, 128.F, 255.F, 128.F, 255.F };
	60
	61	static int tiffcvt(TIFF* in, TIFF* out);
	62	static void usage(int code);
	63	static void setupLumaTables(void);
	64
65	int
66	main(int argc, char* argv[])
67	{
68	TIFF in, out;
69	int c;
70	extern int optind;
71	extern char *optarg;
72
73	while ((c = getopt(argc, argv, "c:h:r:v:z")) != -1)
74	switch (c) {
75	case 'c':
76	if (streq(optarg, "none"))
77	compression = COMPRESSION_NONE;
78	else if (streq(optarg, "packbits"))
79	compression = COMPRESSION_PACKBITS;
80	else if (streq(optarg, "lzw"))
81	compression = COMPRESSION_LZW;
82	else if (streq(optarg, "jpeg"))
83	compression = COMPRESSION_JPEG;
84	else if (streq(optarg, "zip"))
85	compression = COMPRESSION_ADOBE_DEFLATE;
86	else
87	usage(-1);
88	break;
89	case 'h':
90	horizSubSampling = atoi(optarg);
91	break;
92	case 'v':
93	vertSubSampling = atoi(optarg);
94	break;
95	case 'r':
96	rowsperstrip = atoi(optarg);
97	break;
98	case 'z': /* CCIR Rec 601-1 w/ headroom/footroom */
99	refBlackWhite[0] = 16.;
100	refBlackWhite[1] = 235.;
101	refBlackWhite[2] = 128.;
102	refBlackWhite[3] = 240.;
103	refBlackWhite[4] = 128.;
104	refBlackWhite[5] = 240.;
105	break;
106	case '?':
107	usage(0);
108	/NOTREACHED/
109	}
110	if (argc - optind < 2)
111	usage(-1);
112	out = TIFFOpen(argv[argc-1], "w");
113	if (out == NULL)
114	return (-2);
115	setupLumaTables();
116	for (; optind < argc-1; optind++) {
117	in = TIFFOpen(argv[optind], "r");
118	if (in != NULL) {
119	do {
120	if (!tiffcvt(in, out) \|\|
121	!TIFFWriteDirectory(out)) {
122	(void) TIFFClose(out);
123	return (1);
124	}
125	} while (TIFFReadDirectory(in));
126	(void) TIFFClose(in);
127	}
128	}
129	(void) TIFFClose(out);
130	return (0);
131	}
132
133	float *lumaRed;
134	float *lumaGreen;
135	float *lumaBlue;
136	float D1, D2;
137	int Yzero;
138
139	static float*
140	setupLuma(float c)
141	{
142	float v = (float )_TIFFmalloc(256 * sizeof (float));
143	int i;
144	for (i = 0; i < 256; i++)
145	v[i] = c * i;
146	return (v);
147	}
148
149	static unsigned
150	V2Code(float f, float RB, float RW, int CR)
151	{
152	unsigned int c = (unsigned int)((((f)*(RW-RB)/CR)+RB)+.5);
153	return (c > 255 ? 255 : c);
154	}
155
156	static void
157	setupLumaTables(void)
158	{
159	lumaRed = setupLuma(LumaRed);
160	lumaGreen = setupLuma(LumaGreen);
161	lumaBlue = setupLuma(LumaBlue);
162	D1 = 1.F/(2.F - 2.F*LumaBlue);
163	D2 = 1.F/(2.F - 2.F*LumaRed);
164	Yzero = V2Code(0, refBlackWhite[0], refBlackWhite[1], 255);
165	}
166
167	static void
168	cvtClump(unsigned char* op, uint32* raster, uint32 ch, uint32 cw, uint32 w)
169	{
170	float Y, Cb = 0, Cr = 0;
171	uint32 j, k;
172	/*
173	* Convert ch-by-cw block of RGB
174	* to YCbCr and sample accordingly.
175	*/
176	for (k = 0; k < ch; k++) {
177	for (j = 0; j < cw; j++) {
178	uint32 RGB = (raster - k*w)[j];
179	Y = lumaRed[TIFFGetR(RGB)] +
180	lumaGreen[TIFFGetG(RGB)] +
181	lumaBlue[TIFFGetB(RGB)];
182	/* accumulate chrominance */
183	Cb += (TIFFGetB(RGB) - Y) * D1;
184	Cr += (TIFFGetR(RGB) - Y) * D2;
185	/* emit luminence */
186	*op++ = V2Code(Y,
187	refBlackWhite[0], refBlackWhite[1], 255);
188	}
189	for (; j < horizSubSampling; j++)
190	*op++ = Yzero;
191	}
192	for (; k < vertSubSampling; k++) {
193	for (j = 0; j < horizSubSampling; j++)
194	*op++ = Yzero;
195	}
196	/* emit sampled chrominance values */
197	op++ = V2Code(Cb / (chcw), refBlackWhite[2], refBlackWhite[3], 127);
198	op++ = V2Code(Cr / (chcw), refBlackWhite[4], refBlackWhite[5], 127);
199	}
200	#undef LumaRed
201	#undef LumaGreen
202	#undef LumaBlue
203	#undef V2Code
204
205	/*
206	* Convert a strip of RGB data to YCbCr and
207	* sample to generate the output data.
208	*/
209	static void
210	cvtStrip(unsigned char* op, uint32* raster, uint32 nrows, uint32 width)
211	{
212	uint32 x;
213	int clumpSize = vertSubSampling * horizSubSampling + 2;
214	uint32 *tp;
215
216	for (; nrows >= vertSubSampling; nrows -= vertSubSampling) {
217	tp = raster;
218	for (x = width; x >= horizSubSampling; x -= horizSubSampling) {
219	cvtClump(op, tp,
220	vertSubSampling, horizSubSampling, width);
221	op += clumpSize;
222	tp += horizSubSampling;
223	}
224	if (x > 0) {
225	cvtClump(op, tp, vertSubSampling, x, width);
226	op += clumpSize;
227	}
228	raster -= vertSubSampling*width;
229	}
230	if (nrows > 0) {
231	tp = raster;
232	for (x = width; x >= horizSubSampling; x -= horizSubSampling) {
233	cvtClump(op, tp, nrows, horizSubSampling, width);
234	op += clumpSize;
235	tp += horizSubSampling;
236	}
237	if (x > 0)
238	cvtClump(op, tp, nrows, x, width);
239	}
240	}
241
242	static int
243	cvtRaster(TIFF* tif, uint32* raster, uint32 width, uint32 height)
244	{
245	uint32 y;
246	tstrip_t strip = 0;
247	tsize_t cc, acc;
248	unsigned char* buf;
249	uint32 rwidth = roundup(width, horizSubSampling);
250	uint32 rheight = roundup(height, vertSubSampling);
251	uint32 nrows = (rowsperstrip > rheight ? rheight : rowsperstrip);
252	uint32 rnrows = roundup(nrows,vertSubSampling);
253
254	cc = rnrows*rwidth +
255	2((rnrowsrwidth) / (horizSubSampling*vertSubSampling));
256	buf = (unsigned char*)_TIFFmalloc(cc);
257	for (y = height; (int32) y > 0; y -= nrows) {
258	uint32 nr = (y > nrows ? nrows : y);
259	cvtStrip(buf, raster + (y-1)*width, nr, width);
260	nr = roundup(nr, vertSubSampling);
261	acc = nr*rwidth +
262	2((nrrwidth)/(horizSubSampling*vertSubSampling));
263	if (!TIFFWriteEncodedStrip(tif, strip++, buf, acc)) {
264	_TIFFfree(buf);
265	return (0);
266	}
267	}
268	_TIFFfree(buf);
269	return (1);
270	}
271
272	static int
273	tiffcvt(TIFF* in, TIFF* out)
274	{
275	uint32 width, height; /* image width & height */
276	uint32* raster; /* retrieve RGBA image */
277	uint16 shortv;
278	float floatv;
279	char *stringv;
280	uint32 longv;
281
282	TIFFGetField(in, TIFFTAG_IMAGEWIDTH, &width);
283	TIFFGetField(in, TIFFTAG_IMAGELENGTH, &height);
284	raster = (uint32)_TIFFmalloc(width height * sizeof (uint32));
285	if (raster == 0) {
286	TIFFError(TIFFFileName(in), "No space for raster buffer");
287	return (0);
288	}
289	if (!TIFFReadRGBAImage(in, width, height, raster, 0)) {
290	_TIFFfree(raster);
291	return (0);
292	}
293
294	CopyField(TIFFTAG_SUBFILETYPE, longv);
295	TIFFSetField(out, TIFFTAG_IMAGEWIDTH, width);
296	TIFFSetField(out, TIFFTAG_IMAGELENGTH, height);
297	TIFFSetField(out, TIFFTAG_BITSPERSAMPLE, 8);
298	TIFFSetField(out, TIFFTAG_COMPRESSION, compression);
299	TIFFSetField(out, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_YCBCR);
300	if (compression == COMPRESSION_JPEG)
301	TIFFSetField(out, TIFFTAG_JPEGCOLORMODE, JPEGCOLORMODE_RAW);
302	CopyField(TIFFTAG_FILLORDER, shortv);
303	TIFFSetField(out, TIFFTAG_ORIENTATION, ORIENTATION_TOPLEFT);
304	TIFFSetField(out, TIFFTAG_SAMPLESPERPIXEL, 3);
305	CopyField(TIFFTAG_XRESOLUTION, floatv);
306	CopyField(TIFFTAG_YRESOLUTION, floatv);
307	CopyField(TIFFTAG_RESOLUTIONUNIT, shortv);
308	TIFFSetField(out, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG);
309	{ char buf[2048];
310	char *cp = strrchr(TIFFFileName(in), '/');
311	sprintf(buf, "YCbCr conversion of %s", cp ? cp+1 : TIFFFileName(in));
312	TIFFSetField(out, TIFFTAG_IMAGEDESCRIPTION, buf);
313	}
314	TIFFSetField(out, TIFFTAG_SOFTWARE, TIFFGetVersion());
315	CopyField(TIFFTAG_DOCUMENTNAME, stringv);
316
317	TIFFSetField(out, TIFFTAG_REFERENCEBLACKWHITE, refBlackWhite);
318	TIFFSetField(out, TIFFTAG_YCBCRSUBSAMPLING,
319	horizSubSampling, vertSubSampling);
320	TIFFSetField(out, TIFFTAG_YCBCRPOSITIONING, YCBCRPOSITION_CENTERED);
321	TIFFSetField(out, TIFFTAG_YCBCRCOEFFICIENTS, ycbcrCoeffs);
322	rowsperstrip = TIFFDefaultStripSize(out, rowsperstrip);
323	TIFFSetField(out, TIFFTAG_ROWSPERSTRIP, rowsperstrip);
324
325	return (cvtRaster(out, raster, width, height));
326	}
327
328	char* stuff[] = {
329	"usage: rgb2ycbcr [-c comp] [-r rows] [-h N] [-v N] input... output\n",
330	"where comp is one of the following compression algorithms:\n",
331	" jpeg\t\tJPEG encoding\n",
332	" lzw\t\tLempel-Ziv & Welch encoding\n",
333	" zip\t\tdeflate encoding\n",
334	" packbits\tPackBits encoding (default)\n",
335	" none\t\tno compression\n",
336	"and the other options are:\n",
337	" -r\trows/strip\n",
338	" -h\thorizontal sampling factor (1,2,4)\n",
339	" -v\tvertical sampling factor (1,2,4)\n",
340	NULL
341	};
342
343	static void
344	usage(int code)
345	{
346	char buf[BUFSIZ];
347	int i;
348
349	setbuf(stderr, buf);
350
351	fprintf(stderr, "%s\n\n", TIFFGetVersion());
352	for (i = 0; stuff[i] != NULL; i++)
353	fprintf(stderr, "%s\n", stuff[i]);
354	exit(code);
355	}
356
357	/* vim: set ts=8 sts=8 sw=8 noet: */