[wxWidgets.git] / src / tiff / tif_tile.c

/* $Header$ */

/*
 * 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.
 */

/*
 * TIFF Library.
 *
 * Tiled Image Support Routines.
 */
#include "tiffiop.h"

static uint32
summarize(TIFF* tif, size_t summand1, size_t summand2, const char* where)
{
	uint32	bytes = summand1 + summand2;

	if (bytes - summand1 != summand2) {
		TIFFError(tif->tif_name, "Integer overflow in %s", where);
		bytes = 0;
	}

	return (bytes);
}

static uint32
multiply(TIFF* tif, size_t nmemb, size_t elem_size, const char* where)
{
	uint32	bytes = nmemb * elem_size;

	if (elem_size && bytes / elem_size != nmemb) {
		TIFFError(tif->tif_name, "Integer overflow in %s", where);
		bytes = 0;
	}

	return (bytes);
}

/*
 * Compute which tile an (x,y,z,s) value is in.
 */
ttile_t
TIFFComputeTile(TIFF* tif, uint32 x, uint32 y, uint32 z, tsample_t s)
{
	TIFFDirectory *td = &tif->tif_dir;
	uint32 dx = td->td_tilewidth;
	uint32 dy = td->td_tilelength;
	uint32 dz = td->td_tiledepth;
	ttile_t tile = 1;

	if (td->td_imagedepth == 1)
		z = 0;
	if (dx == (uint32) -1)
		dx = td->td_imagewidth;
	if (dy == (uint32) -1)
		dy = td->td_imagelength;
	if (dz == (uint32) -1)
		dz = td->td_imagedepth;
	if (dx != 0 && dy != 0 && dz != 0) {
		uint32 xpt = TIFFhowmany(td->td_imagewidth, dx); 
		uint32 ypt = TIFFhowmany(td->td_imagelength, dy); 
		uint32 zpt = TIFFhowmany(td->td_imagedepth, dz); 

		if (td->td_planarconfig == PLANARCONFIG_SEPARATE) 
			tile = (xpt*ypt*zpt)*s +
			     (xpt*ypt)*(z/dz) +
			     xpt*(y/dy) +
			     x/dx;
		else
			tile = (xpt*ypt)*(z/dz) + xpt*(y/dy) + x/dx;
	}
	return (tile);
}

/*
 * Check an (x,y,z,s) coordinate
 * against the image bounds.
 */
int
TIFFCheckTile(TIFF* tif, uint32 x, uint32 y, uint32 z, tsample_t s)
{
	TIFFDirectory *td = &tif->tif_dir;

	if (x >= td->td_imagewidth) {
		TIFFError(tif->tif_name, "Col %ld out of range, max %lu",
		    (long) x, (u_long) td->td_imagewidth);
		return (0);
	}
	if (y >= td->td_imagelength) {
		TIFFError(tif->tif_name, "Row %ld out of range, max %lu",
		    (long) y, (u_long) td->td_imagelength);
		return (0);
	}
	if (z >= td->td_imagedepth) {
		TIFFError(tif->tif_name, "Depth %ld out of range, max %lu",
		    (long) z, (u_long) td->td_imagedepth);
		return (0);
	}
	if (td->td_planarconfig == PLANARCONFIG_SEPARATE &&
	    s >= td->td_samplesperpixel) {
		TIFFError(tif->tif_name, "Sample %d out of range, max %u",
		    (int) s, td->td_samplesperpixel);
		return (0);
	}
	return (1);
}

/*
 * Compute how many tiles are in an image.
 */
ttile_t
TIFFNumberOfTiles(TIFF* tif)
{
	TIFFDirectory *td = &tif->tif_dir;
	uint32 dx = td->td_tilewidth;
	uint32 dy = td->td_tilelength;
	uint32 dz = td->td_tiledepth;
	ttile_t ntiles;

	if (dx == (uint32) -1)
		dx = td->td_imagewidth;
	if (dy == (uint32) -1)
		dy = td->td_imagelength;
	if (dz == (uint32) -1)
		dz = td->td_imagedepth;
	ntiles = (dx == 0 || dy == 0 || dz == 0) ? 0 :
	    multiply(tif, multiply(tif, TIFFhowmany(td->td_imagewidth, dx),
				   TIFFhowmany(td->td_imagelength, dy),
				   "TIFFNumberOfTiles"),
		     TIFFhowmany(td->td_imagedepth, dz), "TIFFNumberOfTiles");
	if (td->td_planarconfig == PLANARCONFIG_SEPARATE)
		ntiles = multiply(tif, ntiles, td->td_samplesperpixel,
				  "TIFFNumberOfTiles");
	return (ntiles);
}

/*
 * Compute the # bytes in each row of a tile.
 */
tsize_t
TIFFTileRowSize(TIFF* tif)
{
	TIFFDirectory *td = &tif->tif_dir;
	tsize_t rowsize;
	
	if (td->td_tilelength == 0 || td->td_tilewidth == 0)
		return ((tsize_t) 0);
	rowsize = multiply(tif, td->td_bitspersample, td->td_tilewidth,
			   "TIFFTileRowSize");
	if (td->td_planarconfig == PLANARCONFIG_CONTIG)
		rowsize = multiply(tif, rowsize, td->td_samplesperpixel,
				   "TIFFTileRowSize");
	return ((tsize_t) TIFFhowmany8(rowsize));
}

/*
 * Compute the # bytes in a variable length, row-aligned tile.
 */
tsize_t
TIFFVTileSize(TIFF* tif, uint32 nrows)
{
	TIFFDirectory *td = &tif->tif_dir;
	tsize_t tilesize;

	if (td->td_tilelength == 0 || td->td_tilewidth == 0 ||
	    td->td_tiledepth == 0)
		return ((tsize_t) 0);
#ifdef YCBCR_SUPPORT
	if (td->td_planarconfig == PLANARCONFIG_CONTIG &&
	    td->td_photometric == PHOTOMETRIC_YCBCR &&
	    !isUpSampled(tif)) {
		/*
		 * Packed YCbCr data contain one Cb+Cr for every
		 * HorizontalSampling*VerticalSampling Y values.
		 * Must also roundup width and height when calculating
		 * since images that are not a multiple of the
		 * horizontal/vertical subsampling area include
		 * YCbCr data for the extended image.
		 */
		tsize_t w =
		    TIFFroundup(td->td_tilewidth, td->td_ycbcrsubsampling[0]);
		tsize_t rowsize =
		    TIFFhowmany8(multiply(tif, w, td->td_bitspersample,
					  "TIFFVTileSize"));
		tsize_t samplingarea =
		    td->td_ycbcrsubsampling[0]*td->td_ycbcrsubsampling[1];
		nrows = TIFFroundup(nrows, td->td_ycbcrsubsampling[1]);
		/* NB: don't need TIFFhowmany here 'cuz everything is rounded */
		tilesize = multiply(tif, nrows, rowsize, "TIFFVTileSize");
		tilesize = summarize(tif, tilesize,
				     multiply(tif, 2, tilesize / samplingarea,
					      "TIFFVTileSize"),
				     "TIFFVTileSize");
	} else
#endif
		tilesize = multiply(tif, nrows, TIFFTileRowSize(tif),
				    "TIFFVTileSize");
	return ((tsize_t)
	    multiply(tif, tilesize, td->td_tiledepth, "TIFFVTileSize"));
}

/*
 * Compute the # bytes in a row-aligned tile.
 */
tsize_t
TIFFTileSize(TIFF* tif)
{
	return (TIFFVTileSize(tif, tif->tif_dir.td_tilelength));
}

/*
 * Compute a default tile size based on the image
 * characteristics and a requested value.  If a
 * request is <1 then we choose a size according
 * to certain heuristics.
 */
void
TIFFDefaultTileSize(TIFF* tif, uint32* tw, uint32* th)
{
	(*tif->tif_deftilesize)(tif, tw, th);
}

void
_TIFFDefaultTileSize(TIFF* tif, uint32* tw, uint32* th)
{
	(void) tif;
	if (*(int32*) tw < 1)
		*tw = 256;
	if (*(int32*) th < 1)
		*th = 256;
	/* roundup to a multiple of 16 per the spec */
	if (*tw & 0xf)
		*tw = TIFFroundup(*tw, 16);
	if (*th & 0xf)
		*th = TIFFroundup(*th, 16);
}
Commit	Line	Data
b47c832e RR	1	/* $Header$ */
	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	/*
	28	* TIFF Library.
	29	*
	30	* Tiled Image Support Routines.
	31	*/
	32	#include "tiffiop.h"
	33
00cb87b4 VZ	34	static uint32
	35	summarize(TIFF* tif, size_t summand1, size_t summand2, const char* where)
	36	{
	37	uint32 bytes = summand1 + summand2;
	38
	39	if (bytes - summand1 != summand2) {
	40	TIFFError(tif->tif_name, "Integer overflow in %s", where);
	41	bytes = 0;
	42	}
	43
	44	return (bytes);
	45	}
	46
	47	static uint32
	48	multiply(TIFF* tif, size_t nmemb, size_t elem_size, const char* where)
	49	{
	50	uint32 bytes = nmemb * elem_size;
	51
	52	if (elem_size && bytes / elem_size != nmemb) {
	53	TIFFError(tif->tif_name, "Integer overflow in %s", where);
	54	bytes = 0;
	55	}
	56
	57	return (bytes);
	58	}
	59
b47c832e RR	60	/*
	61	* Compute which tile an (x,y,z,s) value is in.
	62	*/
	63	ttile_t
	64	TIFFComputeTile(TIFF* tif, uint32 x, uint32 y, uint32 z, tsample_t s)
	65	{
	66	TIFFDirectory *td = &tif->tif_dir;
	67	uint32 dx = td->td_tilewidth;
	68	uint32 dy = td->td_tilelength;
	69	uint32 dz = td->td_tiledepth;
	70	ttile_t tile = 1;
	71
	72	if (td->td_imagedepth == 1)
	73	z = 0;
	74	if (dx == (uint32) -1)
	75	dx = td->td_imagewidth;
	76	if (dy == (uint32) -1)
	77	dy = td->td_imagelength;
	78	if (dz == (uint32) -1)
	79	dz = td->td_imagedepth;
	80	if (dx != 0 && dy != 0 && dz != 0) {
	81	uint32 xpt = TIFFhowmany(td->td_imagewidth, dx);
	82	uint32 ypt = TIFFhowmany(td->td_imagelength, dy);
	83	uint32 zpt = TIFFhowmany(td->td_imagedepth, dz);
	84
	85	if (td->td_planarconfig == PLANARCONFIG_SEPARATE)
	86	tile = (xptyptzpt)*s +
	87	(xptypt)(z/dz) +
	88	xpt*(y/dy) +
	89	x/dx;
	90	else
00cb87b4	91	tile = (xptypt)(z/dz) + xpt*(y/dy) + x/dx;
b47c832e RR	92	}
	93	return (tile);
	94	}
	95
	96	/*
	97	* Check an (x,y,z,s) coordinate
	98	* against the image bounds.
	99	*/
	100	int
	101	TIFFCheckTile(TIFF* tif, uint32 x, uint32 y, uint32 z, tsample_t s)
	102	{
	103	TIFFDirectory *td = &tif->tif_dir;
	104
	105	if (x >= td->td_imagewidth) {
	106	TIFFError(tif->tif_name, "Col %ld out of range, max %lu",
	107	(long) x, (u_long) td->td_imagewidth);
	108	return (0);
	109	}
	110	if (y >= td->td_imagelength) {
	111	TIFFError(tif->tif_name, "Row %ld out of range, max %lu",
	112	(long) y, (u_long) td->td_imagelength);
	113	return (0);
	114	}
	115	if (z >= td->td_imagedepth) {
	116	TIFFError(tif->tif_name, "Depth %ld out of range, max %lu",
	117	(long) z, (u_long) td->td_imagedepth);
	118	return (0);
	119	}
	120	if (td->td_planarconfig == PLANARCONFIG_SEPARATE &&
	121	s >= td->td_samplesperpixel) {
	122	TIFFError(tif->tif_name, "Sample %d out of range, max %u",
	123	(int) s, td->td_samplesperpixel);
	124	return (0);
	125	}
	126	return (1);
	127	}
	128
	129	/*
	130	* Compute how many tiles are in an image.
	131	*/
	132	ttile_t
	133	TIFFNumberOfTiles(TIFF* tif)
	134	{
	135	TIFFDirectory *td = &tif->tif_dir;
	136	uint32 dx = td->td_tilewidth;
	137	uint32 dy = td->td_tilelength;
	138	uint32 dz = td->td_tiledepth;
	139	ttile_t ntiles;
	140
	141	if (dx == (uint32) -1)
	142	dx = td->td_imagewidth;
	143	if (dy == (uint32) -1)
	144	dy = td->td_imagelength;
	145	if (dz == (uint32) -1)
	146	dz = td->td_imagedepth;
	147	ntiles = (dx == 0 \|\| dy == 0 \|\| dz == 0) ? 0 :
00cb87b4 VZ	148	multiply(tif, multiply(tif, TIFFhowmany(td->td_imagewidth, dx),
	149	TIFFhowmany(td->td_imagelength, dy),
	150	"TIFFNumberOfTiles"),
	151	TIFFhowmany(td->td_imagedepth, dz), "TIFFNumberOfTiles");
b47c832e	152	if (td->td_planarconfig == PLANARCONFIG_SEPARATE)
00cb87b4 VZ	153	ntiles = multiply(tif, ntiles, td->td_samplesperpixel,
00cb87b4 VZ	154	"TIFFNumberOfTiles");
b47c832e RR	155	return (ntiles);
	156	}
	157
	158	/*
	159	* Compute the # bytes in each row of a tile.
	160	*/
	161	tsize_t
	162	TIFFTileRowSize(TIFF* tif)
	163	{
	164	TIFFDirectory *td = &tif->tif_dir;
	165	tsize_t rowsize;
	166
	167	if (td->td_tilelength == 0 \|\| td->td_tilewidth == 0)
	168	return ((tsize_t) 0);
00cb87b4 VZ	169	rowsize = multiply(tif, td->td_bitspersample, td->td_tilewidth,
00cb87b4 VZ	170	"TIFFTileRowSize");
b47c832e	171	if (td->td_planarconfig == PLANARCONFIG_CONTIG)
00cb87b4 VZ	172	rowsize = multiply(tif, rowsize, td->td_samplesperpixel,
	173	"TIFFTileRowSize");
	174	return ((tsize_t) TIFFhowmany8(rowsize));
b47c832e RR	175	}
	176
	177	/*
	178	* Compute the # bytes in a variable length, row-aligned tile.
	179	*/
	180	tsize_t
	181	TIFFVTileSize(TIFF* tif, uint32 nrows)
	182	{
	183	TIFFDirectory *td = &tif->tif_dir;
	184	tsize_t tilesize;
	185
	186	if (td->td_tilelength == 0 \|\| td->td_tilewidth == 0 \|\|
	187	td->td_tiledepth == 0)
	188	return ((tsize_t) 0);
	189	#ifdef YCBCR_SUPPORT
	190	if (td->td_planarconfig == PLANARCONFIG_CONTIG &&
	191	td->td_photometric == PHOTOMETRIC_YCBCR &&
	192	!isUpSampled(tif)) {
	193	/*
	194	* Packed YCbCr data contain one Cb+Cr for every
	195	* HorizontalSampling*VerticalSampling Y values.
	196	* Must also roundup width and height when calculating
	197	* since images that are not a multiple of the
	198	* horizontal/vertical subsampling area include
	199	* YCbCr data for the extended image.
	200	*/
	201	tsize_t w =
	202	TIFFroundup(td->td_tilewidth, td->td_ycbcrsubsampling[0]);
00cb87b4 VZ	203	tsize_t rowsize =
	204	TIFFhowmany8(multiply(tif, w, td->td_bitspersample,
	205	"TIFFVTileSize"));
b47c832e RR	206	tsize_t samplingarea =
	207	td->td_ycbcrsubsampling[0]*td->td_ycbcrsubsampling[1];
	208	nrows = TIFFroundup(nrows, td->td_ycbcrsubsampling[1]);
	209	/* NB: don't need TIFFhowmany here 'cuz everything is rounded */
00cb87b4 VZ	210	tilesize = multiply(tif, nrows, rowsize, "TIFFVTileSize");
	211	tilesize = summarize(tif, tilesize,
	212	multiply(tif, 2, tilesize / samplingarea,
	213	"TIFFVTileSize"),
	214	"TIFFVTileSize");
b47c832e RR	215	} else
b47c832e RR	216	#endif
00cb87b4 VZ	217	tilesize = multiply(tif, nrows, TIFFTileRowSize(tif),
	218	"TIFFVTileSize");
	219	return ((tsize_t)
	220	multiply(tif, tilesize, td->td_tiledepth, "TIFFVTileSize"));
b47c832e RR	221	}
	222
	223	/*
	224	* Compute the # bytes in a row-aligned tile.
	225	*/
	226	tsize_t
	227	TIFFTileSize(TIFF* tif)
	228	{
	229	return (TIFFVTileSize(tif, tif->tif_dir.td_tilelength));
	230	}
	231
	232	/*
	233	* Compute a default tile size based on the image
	234	* characteristics and a requested value. If a
	235	* request is <1 then we choose a size according
	236	* to certain heuristics.
	237	*/
	238	void
	239	TIFFDefaultTileSize(TIFF* tif, uint32* tw, uint32* th)
	240	{
	241	(*tif->tif_deftilesize)(tif, tw, th);
	242	}
	243
	244	void
	245	_TIFFDefaultTileSize(TIFF* tif, uint32* tw, uint32* th)
	246	{
	247	(void) tif;
	248	if ((int32) tw < 1)
	249	*tw = 256;
	250	if ((int32) th < 1)
	251	*th = 256;
	252	/* roundup to a multiple of 16 per the spec */
	253	if (*tw & 0xf)
	254	tw = TIFFroundup(tw, 16);
	255	if (*th & 0xf)
	256	th = TIFFroundup(th, 16);
	257	}