Fixed the handling of 16bit TIFF files for cinema compression. Modified "convert.c".

Fixed the parameters used for cinema compression (9-7 transform used instead of 5-3). Modified "image_to_j2k.c"

ChangeLog

@@ -5,6 +5,21 @@ What's New for OpenJPEG
 ! : changed
 + : added
 
+May 31, 2007
+* [FOD] Fixed the handling of 16bit TIFF files for cinema compression. Modified "convert.c"
+* [FOD] Fixed the parameters used for cinema compression (9-7 transform used instead of 5-3). Modified "image_to_j2k.c"
+
+May 24, 2007
+* [FOD] Bug fixed by Sylvain Munaut. Change in the reading of the POC marker. Since COD/COC can be anywhere in the header, the decoder cannot always know while decoding the POC marker the value of numlayers and numresolution.
+
+May 23, 2007
+! [FOD] Patch suggested by Callum Lerwick <seg@haxxed.com>: "This makes the t1 data arrays dynamic, which greatly reduces cache thrashing. Also, some minor cleanup to prevent unnecessary casts"
+
+May 22, 2007
+! [FOD] Patch suggested by Callum Lerwick <seg@haxxed.com>: "Some formatting cleanups, so that the long function definitions and calls fit on screen. Use of prefix increment which is theoretically faster, in practice any sane compiler can optimize a postfix increment but its best not to count on such things. Consolidation of some redundant calculations in the inner loops, which becomes very useful in the future autovectorize patch."
+! [FOD] Patch suggested by Callum Lerwick <seg@haxxed.com>: "This changes the flag code in t1 to use a flag_t type, which can then be changed to reduce RAM usage. It is now typedef to a short."
+! [FOD] Patch suggested by Callum Lerwick <seg@haxxed.com>: "This patch makes the t1 LUTs static. I actually intend this as a prelude to possibly eliminating some or all of the LUTs entirely."
+
 May 11, 2007
 - [FOD]	JP3D library Binaries removed from SVN.
 ! [FOD] MJ2 codec updated from LibOpenJPEG version 0.97 to LibOpenJPEG version 1. Hence, the MJ2 codec will now take advantage of all the improvements and optimizations done on the LibOpenJPEG library.

codec/convert.c

@@ -1082,7 +1082,7 @@ typedef struct tiff_infoheader{
 }tiff_infoheader_t;
 
 int imagetotif(opj_image_t * image, const char *outfile) {
-	int width, height, imgsize;
+	int width, height;
 	int bps,index;
 	TIFF *tif;
 	tdata_t buf;
@@ -1108,7 +1108,6 @@ int imagetotif(opj_image_t * image, const char *outfile) {
 
 			width	= image->comps[0].w;
 			height= image->comps[0].h;
-			imgsize = image->comps[0].w * image->comps[0].h ;
 			bps		= image->comps[0].prec;
 			/* Set tags */
 			TIFFSetField(tif, TIFFTAG_IMAGEWIDTH, width);
@@ -1130,74 +1129,34 @@ int imagetotif(opj_image_t * image, const char *outfile) {
 				int i;
 				dat8 = buf;
 				if (image->comps[0].prec == 8){
-					for (i=0; i < strip_size; i+=3) {	// 8 bits per pixel 
+					for (i=0; i<TIFFStripSize(tif); i+=3) {	// 8 bits per pixel 
-						int r = 0,g = 0,b = 0;
+						dat8[i+0] = image->comps[0].data[index] ;	// R 
-						if(index < imgsize){
+						dat8[i+1] = image->comps[1].data[index] ;	// G 
-							r = image->comps[0].data[index];
+						dat8[i+2] = image->comps[2].data[index] ;	// B 
-							r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-							g = image->comps[1].data[index];
-							g += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-							b = image->comps[2].data[index];
-							b += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-
-							dat8[i+0] = r ;	// R 
-							dat8[i+1] = g ;	// G 
-							dat8[i+2] = b ;	// B 
 						index++;
-						}else
-							break;
 					}
 				}else if (image->comps[0].prec == 12){
-					for (i=0; i<strip_size; i+=9) {	// 12 bits per pixel 
+					for (i=0; i<TIFFStripSize(tif); i+=9) {	// 12 bits per pixel 
-						int r = 0,g = 0,b = 0;
+						dat8[i+0] = (image->comps[0].data[index]>>8)<<4 | (image->comps[0].data[index]>>4);
-						int r1 = 0,g1 = 0,b1 = 0;
+						dat8[i+1] = (image->comps[0].data[index]<<4)|((image->comps[1].data[index]>>8)& 0x0f);
-						if(index < imgsize){
+						dat8[i+2] = (image->comps[1].data[index]);
-							r = image->comps[0].data[index];
+						dat8[i+3] = (image->comps[2].data[index]>>8)<<4 | (image->comps[2].data[index]>>4);
-							r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
+						dat8[i+4] = (image->comps[2].data[index]<<4)|((image->comps[1].data[index+1]>>8)& 0x0f);
-							g = image->comps[1].data[index];
+						dat8[i+5] = (image->comps[0].data[index+1]);
-							g += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
+						dat8[i+6] = (image->comps[1].data[index+1]>>8)<<4 | (image->comps[1].data[index+1]>>4);
-							b = image->comps[2].data[index];
+						dat8[i+7] = (image->comps[1].data[index+1]<<4)|((image->comps[2].data[index+1]>>8)& 0x0f);
-							b += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
+						dat8[i+8] = (image->comps[2].data[index+1]);
-							r1 = image->comps[0].data[index+1];
-							r1 += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-							g1 = image->comps[1].data[index+1];
-							g1 += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-							b1 = image->comps[2].data[index+1];
-							b1 += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-
-							dat8[i+0] = (r >> 4);
-							dat8[i+1] = ((r & 0x0f) << 4 )|((g >> 8)& 0x0f);
-							dat8[i+2] = g ;
-							dat8[i+3] = (b >> 4);
-							dat8[i+4] = ((b & 0x0f) << 4 )|((r1 >> 8)& 0x0f);
-							dat8[i+5] = (r1);
-							dat8[i+6] = (g1 >> 4);
-							dat8[i+7] = ((g1 & 0x0f)<< 4 )|((b1 >> 8)& 0x0f);
-							dat8[i+8] = (b1);
 						index+=2;
-						}else
-							break;
 					}
 				}else if (image->comps[0].prec == 16){
-					for (i=0; i<strip_size; i+=6) {	// 16 bits per pixel
+					for (i=0; i<TIFFStripSize(tif); i+=6) {	// 16 bits per pixel 
-						int r = 0,g = 0,b = 0;
+						dat8[i+0] =  image->comps[0].data[index];//LSB
-						if(index < imgsize){
+						dat8[i+1] = (image->comps[0].data[index]>> 8);//MSB	 
-							r = image->comps[0].data[index];
+						dat8[i+2] =  image->comps[1].data[index]; 
-							r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
+						dat8[i+3] = (image->comps[1].data[index]>> 8);	
-							g = image->comps[1].data[index];
+						dat8[i+4] =  image->comps[2].data[index];	 
-							g += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
+						dat8[i+5] = (image->comps[2].data[index]>> 8); 
-							b = image->comps[2].data[index];
-							b += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
-
-							dat8[i+0] =  r;//LSB
-							dat8[i+1] = (r >> 8);//MSB	 
-							dat8[i+2] =  g; 
-							dat8[i+3] = (g >> 8);	
-							dat8[i+4] =  b;	 
-							dat8[i+5] = (b >> 8); 
 						index++;
-						}else
-							break;
 					}
 				}else{
 					fprintf(stderr,"Bits=%d, Only 8,12,16 bits implemented\n",image->comps[0].prec);
@@ -1290,7 +1249,6 @@ opj_image_t* tiftoimage(char *filename, opj_cparameters_t *parameters)
 	OPJ_COLOR_SPACE color_space;
 	opj_image_cmptparm_t cmptparm[3];
 	opj_image_t * image = NULL;
-	int imgsize;
 
 	tif = TIFFOpen(filename, "r");
 
@@ -1320,8 +1278,13 @@ opj_image_t* tiftoimage(char *filename, opj_cparameters_t *parameters)
 		/* initialize image components*/ 
 		memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));
 		for(j = 0; j < numcomps; j++) {
+			if ((parameters->cp_cinema) & (Info.tiBps== 16)){
+				cmptparm[j].prec = 12;
+				cmptparm[j].bpp = 12;
+			}else{
 				cmptparm[j].prec = Info.tiBps;
 				cmptparm[j].bpp = Info.tiBps;
+			}
 			cmptparm[j].sgnd = 0;
 			cmptparm[j].dx = subsampling_dx;
 			cmptparm[j].dy = subsampling_dy;
@@ -1345,7 +1308,6 @@ opj_image_t* tiftoimage(char *filename, opj_cparameters_t *parameters)
 		strip_size=0;
 		strip_size=TIFFStripSize(tif);
 		index = 0;
-		imgsize = image->comps[0].w * image->comps[0].h ;
 		/* Read the Image components*/
 		for (strip = 0; strip < TIFFNumberOfStrips(tif); strip++) {
 			unsigned char *dat8;
@@ -1355,7 +1317,6 @@ opj_image_t* tiftoimage(char *filename, opj_cparameters_t *parameters)
 
 			if (Info.tiBps==12){
 				for (i=0; i<ssize; i+=9) {	/*12 bits per pixel*/
-					if(index < imgsize){
 					image->comps[0].data[index]   = ( dat8[i+0]<<4 )		|(dat8[i+1]>>4);
 					image->comps[1].data[index]   = ((dat8[i+1]& 0x0f)<< 8)	| dat8[i+2];
 					image->comps[2].data[index]   = ( dat8[i+3]<<4)			|(dat8[i+4]>>4);
@@ -1363,30 +1324,27 @@ opj_image_t* tiftoimage(char *filename, opj_cparameters_t *parameters)
 					image->comps[1].data[index+1] = ( dat8[i+6] <<4)		|(dat8[i+7]>>4);
 					image->comps[2].data[index+1] = ((dat8[i+7]& 0x0f)<< 8)	| dat8[i+8];
 					index+=2;
-					}else
-						break;
 				}
 			}
 			else if( Info.tiBps==16){
 				for (i=0; i<ssize; i+=6) {	/* 16 bits per pixel */
-					if(index < imgsize){
 					image->comps[0].data[index] = ( dat8[i+1] << 8 ) | dat8[i+0];	// R 
 					image->comps[1].data[index] = ( dat8[i+3] << 8 ) | dat8[i+2];	// G 
 					image->comps[2].data[index] = ( dat8[i+5] << 8 ) | dat8[i+4];	// B 
+					if(parameters->cp_cinema){/* Rounding to 12 bits*/
+						image->comps[0].data[index] = (image->comps[0].data[index] + 0x08) >> 4 ;
+						image->comps[1].data[index] = (image->comps[1].data[index] + 0x08) >> 4 ;
+						image->comps[2].data[index] = (image->comps[2].data[index] + 0x08) >> 4 ;
+                    }
 					index++;
-					}else
-						break;
 				}
 			}
 			else if ( Info.tiBps==8){
 				for (i=0; i<ssize; i+=3) {	/* 8 bits per pixel */
-					if(index < imgsize){
 					image->comps[0].data[index] = dat8[i+0];	// R 
 					image->comps[1].data[index] = dat8[i+1];	// G 
 					image->comps[2].data[index] = dat8[i+2];	// B 
 					index++;
-					}else
-						break;
 				}
 			}
 			else{

codec/image_to_j2k.c

@@ -465,6 +465,10 @@ void cinema_parameters(opj_cparameters_t *parameters){
 	parameters->roi_compno = -1;
 
 	parameters->subsampling_dx = 1;		parameters->subsampling_dy = 1;
+
+	/* 9-7 transform */
+	parameters->irreversible = 1;
+
 }
 
 void cinema_setup_encoder(opj_cparameters_t *parameters,opj_image_t *image){