Форумы GIS-Lab.info

Код: Выделить всё

	TIFF *tif=XTIFFOpen(fname_1, "r");
	TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, &width);
	TIFFGetField(tif, TIFFTAG_IMAGELENGTH, &height);

	uint32 npixels=width*height;
	uint32 *raster=(uint32 *) malloc(npixels *sizeof(uint32));

	TIFFReadRGBAImage(tif, width, height, raster, 0);
	TIFFReadRGBAImageOriented(tif, width, height, raster, ORIENTATION_TOPLEFT, 0);
	TIFFClose(tif);

	unsigned char *data = (unsigned char *)raster; 

	TIFF *tif_=XTIFFOpen(fname_2, "w");
	GTIF *gtif=GTIFNew(tif_);

	TIFFSetField(tif_,TIFFTAG_IMAGEWIDTH, width);
	TIFFSetField(tif_,TIFFTAG_IMAGELENGTH, height);

	int sampleperpixel = 1;
	TIFFSetField(tif_, TIFFTAG_SAMPLESPERPIXEL, sampleperpixel);

	TIFFSetField(tif_,TIFFTAG_BITSPERSAMPLE, 8);
	TIFFSetField(tif_, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG);
	TIFFSetField(tif_,TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_PALETTE);

	tsize_t linebytes = sampleperpixel * width;

	unsigned char *buf = NULL;
	if (TIFFScanlineSize(tif_))
		buf =(unsigned char *)_TIFFmalloc(linebytes);
	else
		buf = (unsigned char *)_TIFFmalloc(TIFFScanlineSize(tif_));

	// We set the strip size of the file to be size of one row of pixels
	TIFFSetField(tif_, TIFFTAG_ROWSPERSTRIP, TIFFDefaultStripSize(tif_, linebytes));

	for (uint32 row = 0; row < height; row++)
	{
		memcpy(buf, &data[row*linebytes], linebytes);
		if (TIFFWriteScanline(tif_, buf, row, 0) < 0)
			break;
	}

	GTIFWriteKeys(gtif);
	GTIFFree(gtif);
	XTIFFClose(tif_);

Код: Выделить всё

int sampleperpixel = 1;
TIFFSetField(tif_,TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_PALETTE);

Код: Выделить всё

int sampleperpixel = 3;
TIFFSetField(tif_,TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_RGB);