fixed "tiffio" header inclusion to use user version on Win32 system only (otherwise...
[openjpeg.git] / codec / convert.c
index d074909fdff10325876e122c4fc34a5478d54f65..a9ffb0da1e08c9554f18743c6cbdace3f523a498 100644 (file)
 #include <stdio.h>
 #include <stdlib.h>
 #include <string.h>
-#include "openjpeg.h"
+#ifdef WIN32
 #include "../libs/libtiff/tiffio.h"
+#else
+#include <tiffio.h>
+#endif /* WIN32 */
+#include "openjpeg.h"
 #include "convert.h"
 
 /*
@@ -99,7 +103,7 @@ typedef struct tga_header
 } tga_header;
 #pragma pack(pop) // Return to normal structure packing alignment.
 
-int tga_readheader(FILE *fp, int *bits_per_pixel, int *width, int *height, int *flip_image)
+int tga_readheader(FILE *fp, uint32 *bits_per_pixel, uint32 *width, uint32 *height, int *flip_image)
 {
        int palette_size;
        tga_header tga ;
@@ -310,15 +314,15 @@ int imagetotga(opj_image_t * image, const char *outfile) {
    }
        }
 
-       width  = int_ceildiv(image->x1-image->x0, image->comps[0].dx);
-       height = int_ceildiv(image->y1-image->y0, image->comps[0].dy);
+       width = image->comps[0].w;
+       height = image->comps[0].h; 
 
        // Mono with alpha, or RGB with alpha.
        write_alpha = (image->numcomps==2) || (image->numcomps==4);   
 
        // Write TGA header 
        bpp = write_alpha ? 32 : 24;
-       if (!tga_writeheader(fdest, bpp, width, height, true))
+       if (!tga_writeheader(fdest, bpp, width , height, true))
                return 1;
 
        alpha_channel = image->numcomps-1; 
@@ -764,7 +768,7 @@ opj_image_t* bmptoimage(const char *filename, opj_cparameters_t *parameters) {
 }
 
 int imagetobmp(opj_image_t * image, const char *outfile) {
-       int w, wr, h, hr;
+       int w, h;
        int i, pad;
        FILE *fdest = NULL;
        int adjustR, adjustG, adjustB;
@@ -786,42 +790,39 @@ int imagetobmp(opj_image_t * image, const char *outfile) {
                        return 1;
                }
            
-               w = image->comps[0].w;
-               wr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);
-           
+               w = image->comps[0].w;      
                h = image->comps[0].h;
-               hr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);
            
                fprintf(fdest, "BM");
            
                /* FILE HEADER */
                /* ------------- */
                fprintf(fdest, "%c%c%c%c",
-                       (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) + 54) & 0xff,
-                       (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) >> 8) & 0xff,
-                       (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) >> 16) & 0xff,
-                       (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) >> 24) & 0xff);
+                       (unsigned char) (h * w * 3 + 3 * h * (w % 2) + 54) & 0xff,
+                       (unsigned char) ((h * w * 3 + 3 * h * (w % 2) + 54)     >> 8) & 0xff,
+                       (unsigned char) ((h * w * 3 + 3 * h * (w % 2) + 54)     >> 16) & 0xff,
+                       (unsigned char) ((h * w * 3 + 3 * h * (w % 2) + 54)     >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,((54) >> 16) & 0xff, ((54) >> 24) & 0xff);
            
                /* INFO HEADER   */
                /* ------------- */
                fprintf(fdest, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,     ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
-               fprintf(fdest, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
-                       (unsigned char) ((wr) >> 8) & 0xff,
-                       (unsigned char) ((wr) >> 16) & 0xff,
-                       (unsigned char) ((wr) >> 24) & 0xff);
-               fprintf(fdest, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
-                       (unsigned char) ((hr) >> 8) & 0xff,
-                       (unsigned char) ((hr) >> 16) & 0xff,
-                       (unsigned char) ((hr) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) ((w) & 0xff),
+                       (unsigned char) ((w) >> 8) & 0xff,
+                       (unsigned char) ((w) >> 16) & 0xff,
+                       (unsigned char) ((w) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) ((h) & 0xff),
+                       (unsigned char) ((h) >> 8) & 0xff,
+                       (unsigned char) ((h) >> 16) & 0xff,
+                       (unsigned char) ((h) >> 24) & 0xff);
                fprintf(fdest, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
                fprintf(fdest, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);
                fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
-               fprintf(fdest, "%c%c%c%c", (unsigned char) (3 * hr * wr + 3 * hr * (wr % 2)) & 0xff,
-                       (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> 8) & 0xff,
-                       (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> 16) & 0xff,
-                       (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) (3 * h * w + 3 * h * (w % 2)) & 0xff,
+                       (unsigned char) ((h * w * 3 + 3 * h * (w % 2)) >> 8) & 0xff,
+                       (unsigned char) ((h * w * 3 + 3 * h * (w % 2)) >> 16) & 0xff,
+                       (unsigned char) ((h * w * 3 + 3 * h * (w % 2)) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
@@ -846,24 +847,24 @@ int imagetobmp(opj_image_t * image, const char *outfile) {
                else 
                        adjustB = 0;
 
-               for (i = 0; i < wr * hr; i++) {
+               for (i = 0; i < w * h; i++) {
                        unsigned char rc, gc, bc;
                        int r, g, b;
                                                        
-                       r = image->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+                       r = image->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
                        r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
                        rc = (unsigned char) ((r >> adjustR)+((r >> (adjustR-1))%2));
-                       g = image->comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+                       g = image->comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
                        g += (image->comps[1].sgnd ? 1 << (image->comps[1].prec - 1) : 0);
                        gc = (unsigned char) ((g >> adjustG)+((g >> (adjustG-1))%2));
-                       b = image->comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+                       b = image->comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
                        b += (image->comps[2].sgnd ? 1 << (image->comps[2].prec - 1) : 0);
                        bc = (unsigned char) ((b >> adjustB)+((b >> (adjustB-1))%2));
 
                        fprintf(fdest, "%c%c%c", bc, gc, rc);
                        
-                       if ((i + 1) % wr == 0) {
-                               for (pad = (3 * wr) % 4 ? 4 - (3 * wr) % 4 : 0; pad > 0; pad--) /* ADD */
+                       if ((i + 1) % w == 0) {
+                               for (pad = (3 * w) % 4 ? 4 - (3 * w) % 4 : 0; pad > 0; pad--)   /* ADD */
                                        fprintf(fdest, "%c", 0);
                        }
                }
@@ -875,20 +876,17 @@ int imagetobmp(opj_image_t * image, const char *outfile) {
                <<-- <<-- <<-- <<-- */
 
                fdest = fopen(outfile, "wb");
-               w = image->comps[0].w;
-               wr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);
-           
+               w = image->comps[0].w;      
                h = image->comps[0].h;
-               hr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);
            
                fprintf(fdest, "BM");
            
                /* FILE HEADER */
                /* ------------- */
-               fprintf(fdest, "%c%c%c%c", (unsigned char) (hr * wr + 54 + 1024 + hr * (wr % 2)) & 0xff,
-                       (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2)) >> 8) & 0xff,
-                       (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2)) >> 16) & 0xff,
-                       (unsigned char) ((hr * wr + 54 + 1024 + wr * (wr % 2)) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) (h * w + 54 + 1024 + h * (w % 2)) & 0xff,
+                       (unsigned char) ((h * w + 54 + 1024 + h * (w % 2)) >> 8) & 0xff,
+                       (unsigned char) ((h * w + 54 + 1024 + h * (w % 2)) >> 16) & 0xff,
+                       (unsigned char) ((h * w + 54 + 1024 + w * (w % 2)) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (54 + 1024) & 0xff, ((54 + 1024) >> 8) & 0xff, 
                        ((54 + 1024) >> 16) & 0xff,
@@ -897,21 +895,21 @@ int imagetobmp(opj_image_t * image, const char *outfile) {
                /* INFO HEADER */
                /* ------------- */
                fprintf(fdest, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,     ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
-               fprintf(fdest, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
-                       (unsigned char) ((wr) >> 8) & 0xff,
-                       (unsigned char) ((wr) >> 16) & 0xff,
-                       (unsigned char) ((wr) >> 24) & 0xff);
-               fprintf(fdest, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
-                       (unsigned char) ((hr) >> 8) & 0xff,
-                       (unsigned char) ((hr) >> 16) & 0xff,
-                       (unsigned char) ((hr) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) ((w) & 0xff),
+                       (unsigned char) ((w) >> 8) & 0xff,
+                       (unsigned char) ((w) >> 16) & 0xff,
+                       (unsigned char) ((w) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) ((h) & 0xff),
+                       (unsigned char) ((h) >> 8) & 0xff,
+                       (unsigned char) ((h) >> 16) & 0xff,
+                       (unsigned char) ((h) >> 24) & 0xff);
                fprintf(fdest, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
                fprintf(fdest, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff);
                fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
-               fprintf(fdest, "%c%c%c%c", (unsigned char) (hr * wr + hr * (wr % 2)) & 0xff,
-                       (unsigned char) ((hr * wr + hr * (wr % 2)) >> 8) &      0xff,
-                       (unsigned char) ((hr * wr + hr * (wr % 2)) >> 16) &     0xff,
-                       (unsigned char) ((hr * wr + hr * (wr % 2)) >> 24) & 0xff);
+               fprintf(fdest, "%c%c%c%c", (unsigned char) (h * w + h * (w % 2)) & 0xff,
+                       (unsigned char) ((h * w + h * (w % 2)) >> 8) &  0xff,
+                       (unsigned char) ((h * w + h * (w % 2)) >> 16) & 0xff,
+                       (unsigned char) ((h * w + h * (w % 2)) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
                fprintf(fdest, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
@@ -920,24 +918,25 @@ int imagetobmp(opj_image_t * image, const char *outfile) {
                if (image->comps[0].prec > 8) {
                        adjustR = image->comps[0].prec - 8;
                        printf("BMP CONVERSION: Truncating component 0 from %d bits to 8 bits\n", image->comps[0].prec);
-               }
+               }else 
+                       adjustR = 0;
 
                for (i = 0; i < 256; i++) {
                        fprintf(fdest, "%c%c%c%c", i, i, i, 0);
                }
 
-               for (i = 0; i < wr * hr; i++) {
+               for (i = 0; i < w * h; i++) {
                        unsigned char rc;
                        int r;
                        
-                       r = image->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+                       r = image->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
                        r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
                        rc = (unsigned char) ((r >> adjustR)+((r >> (adjustR-1))%2));
                        
                        fprintf(fdest, "%c", rc);
 
-                       if ((i + 1) % wr == 0) {
-                               for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--)     /* ADD */
+                       if ((i + 1) % w == 0) {
+                               for (pad = w % 4 ? 4 - w % 4 : 0; pad > 0; pad--)       /* ADD */
                                        fprintf(fdest, "%c", 0);
                        }
                }
@@ -1098,7 +1097,7 @@ opj_image_t* pgxtoimage(const char *filename, opj_cparameters_t *parameters) {
 }
 
 int imagetopgx(opj_image_t * image, const char *outfile) {
-       int w, wr, h, hr;
+       int w, h;
        int i, j, compno;
        FILE *fdest = NULL;
 
@@ -1135,12 +1134,9 @@ int imagetopgx(opj_image_t * image, const char *outfile) {
       }
 
                w = image->comps[compno].w;
-               wr = int_ceildivpow2(image->comps[compno].w, image->comps[compno].factor);
-           
                h = image->comps[compno].h;
-               hr = int_ceildivpow2(image->comps[compno].h, image->comps[compno].factor);
            
-               fprintf(fdest, "PG ML %c %d %d %d\n", comp->sgnd ? '-' : '+', comp->prec, wr, hr);
+               fprintf(fdest, "PG ML %c %d %d %d\n", comp->sgnd ? '-' : '+', comp->prec, w, h);
                if (comp->prec <= 8) {
                        nbytes = 1;
                } else if (comp->prec <= 16) {
@@ -1148,8 +1144,8 @@ int imagetopgx(opj_image_t * image, const char *outfile) {
                } else {
                        nbytes = 4;
                }
-               for (i = 0; i < wr * hr; i++) {
-                       int v = image->comps[compno].data[i / wr * w + i % wr];
+               for (i = 0; i < w * h; i++) {
+                       int v = image->comps[compno].data[i];
                        for (j = nbytes - 1; j >= 0; j--) {
                                char byte = (char) (v >> (j * 8));
                                fwrite(&byte, 1, 1, fdest);
@@ -1267,7 +1263,7 @@ opj_image_t* pnmtoimage(const char *filename, opj_cparameters_t *parameters) {
 }
 
 int imagetopnm(opj_image_t * image, const char *outfile) {
-       int w, wr, wrr, h, hr, hrr, max;
+       int w, wr, h, hr, max;
        int i, compno;
        int adjustR, adjustG, adjustB, adjustX;
        FILE *fdest = NULL;
@@ -1297,18 +1293,16 @@ int imagetopnm(opj_image_t * image, const char *outfile) {
 
                w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx);
                wr = image->comps[0].w;
-               wrr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);
         
                h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy);
                hr = image->comps[0].h;
-               hrr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);
            
                max = image->comps[0].prec > 8 ? 255 : (1 << image->comps[0].prec) - 1;
            
                image->comps[0].x0 = int_ceildivpow2(image->comps[0].x0 - int_ceildiv(image->x0, image->comps[0].dx), image->comps[0].factor);
                image->comps[0].y0 = int_ceildivpow2(image->comps[0].y0 -       int_ceildiv(image->y0, image->comps[0].dy), image->comps[0].factor);
 
-               fprintf(fdest, "P6\n%d %d\n%d\n", wrr, hrr, max);
+               fprintf(fdest, "P6\n%d %d\n%d\n", wr, hr, max);
 
                if (image->comps[0].prec > 8) {
                        adjustR = image->comps[0].prec - 8;
@@ -1330,18 +1324,18 @@ int imagetopnm(opj_image_t * image, const char *outfile) {
                        adjustB = 0;
 
 
-               for (i = 0; i < wrr * hrr; i++) {
+               for (i = 0; i < wr * hr; i++) {
                        int r, g, b;
                        unsigned char rc,gc,bc;
-                       r = image->comps[0].data[i / wrr * wr + i % wrr];
+                       r = image->comps[0].data[i];
                        r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
                        rc = (unsigned char) ((r >> adjustR)+((r >> (adjustR-1))%2));
 
-                       g = image->comps[1].data[i / wrr * wr + i % wrr];
+                       g = image->comps[1].data[i];
                        g += (image->comps[1].sgnd ? 1 << (image->comps[1].prec - 1) : 0);
                        gc = (unsigned char) ((g >> adjustG)+((g >> (adjustG-1))%2));
                        
-                       b = image->comps[2].data[i / wrr * wr + i % wrr];
+                       b = image->comps[2].data[i];
                        b += (image->comps[2].sgnd ? 1 << (image->comps[2].prec - 1) : 0);
                        bc = (unsigned char) ((b >> adjustB)+((b >> (adjustB-1))%2));
                        
@@ -1371,18 +1365,16 @@ int imagetopnm(opj_image_t * image, const char *outfile) {
             
                        w = int_ceildiv(image->x1 - image->x0, image->comps[compno].dx);
                        wr = image->comps[compno].w;
-                       wrr = int_ceildivpow2(image->comps[compno].w, image->comps[compno].factor);
                        
                        h = int_ceildiv(image->y1 - image->y0, image->comps[compno].dy);
                        hr = image->comps[compno].h;
-                       hrr = int_ceildivpow2(image->comps[compno].h, image->comps[compno].factor);
                        
                        max = image->comps[compno].prec > 8 ? 255 : (1 << image->comps[compno].prec) - 1;
                        
                        image->comps[compno].x0 = int_ceildivpow2(image->comps[compno].x0 - int_ceildiv(image->x0, image->comps[compno].dx), image->comps[compno].factor);
                        image->comps[compno].y0 = int_ceildivpow2(image->comps[compno].y0 - int_ceildiv(image->y0, image->comps[compno].dy), image->comps[compno].factor);
                        
-                       fprintf(fdest, "P5\n%d %d\n%d\n", wrr, hrr, max);
+                       fprintf(fdest, "P5\n%d %d\n%d\n", wr, hr, max);
                        
                        if (image->comps[compno].prec > 8) {
                                adjustX = image->comps[0].prec - 8;
@@ -1391,10 +1383,10 @@ int imagetopnm(opj_image_t * image, const char *outfile) {
                        else 
                                adjustX = 0;
                        
-                       for (i = 0; i < wrr * hrr; i++) {
+                       for (i = 0; i < wr * hr; i++) {
                                int l;
                                unsigned char lc;
-                               l = image->comps[compno].data[i / wrr * wr + i % wrr];
+                               l = image->comps[compno].data[i];
                                l += (image->comps[compno].sgnd ? 1 << (image->comps[compno].prec - 1) : 0);
                                lc = (unsigned char) ((l >> adjustX)+((l >> (adjustX-1))%2));
                                fprintf(fdest, "%c", lc);
@@ -1423,8 +1415,9 @@ typedef struct tiff_infoheader{
 }tiff_infoheader_t;
 
 int imagetotif(opj_image_t * image, const char *outfile) {
-       int width, height;
-       int bps,index;
+       int width, height, imgsize;
+       int bps,index,adjust = 0;
+       int last_i=0;
        TIFF *tif;
        tdata_t buf;
        tstrip_t strip;
@@ -1448,7 +1441,8 @@ int imagetotif(opj_image_t * image, const char *outfile) {
                        }
 
                        width   = image->comps[0].w;
-                       height= image->comps[0].h;
+                       height  = image->comps[0].h;
+                       imgsize = width * height ;
                        bps             = image->comps[0].prec;
                        /* Set tags */
                        TIFFSetField(tif, TIFFTAG_IMAGEWIDTH, width);
@@ -1461,43 +1455,166 @@ int imagetotif(opj_image_t * image, const char *outfile) {
                        TIFFSetField(tif, TIFFTAG_ROWSPERSTRIP, 1);
 
                        /* Get a buffer for the data */
-                       buf = _TIFFmalloc(TIFFStripSize(tif));
-                       index=0;
-                       strip_size=0;
                        strip_size=TIFFStripSize(tif);
+                       buf = _TIFFmalloc(strip_size);
+                       index=0;                
+                       adjust = image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0;
                        for (strip = 0; strip < TIFFNumberOfStrips(tif); strip++) {
                                unsigned char *dat8;
-                               int i;
+                               int i, ssize;
+                               ssize = TIFFStripSize(tif);
                                dat8 = buf;
                                if (image->comps[0].prec == 8){
-                                       for (i=0; i<TIFFStripSize(tif); i+=3) { // 8 bits per pixel 
-                                               dat8[i+0] = image->comps[0].data[index] ;       // R 
-                                               dat8[i+1] = image->comps[1].data[index] ;       // G 
-                                               dat8[i+2] = image->comps[2].data[index] ;       // B 
-                                               index++;
+                                       for (i=0; i<ssize-2; i+=3) {    // 8 bits per pixel 
+                                               int r = 0,g = 0,b = 0;
+                                               if(index < imgsize){
+                                                       r = image->comps[0].data[index];
+                                                       g = image->comps[1].data[index];
+                                                       b = image->comps[2].data[index];
+                                                       if (image->comps[0].sgnd){                      
+                                                               r += adjust;
+                                                               g += adjust;
+                                                               b += adjust;
+                                                       }
+                                                       dat8[i+0] = r ; // R 
+                                                       dat8[i+1] = g ; // G 
+                                                       dat8[i+2] = b ; // B 
+                                                       index++;
+                                                       last_i = i+3;
+                                               }else
+                                                       break;
+                                       }
+                                       if(last_i < ssize){
+                                               for (i=last_i; i<ssize; i+=3) { // 8 bits per pixel 
+                                                       int r = 0,g = 0,b = 0;
+                                                       if(index < imgsize){
+                                                               r = image->comps[0].data[index];
+                                                               g = image->comps[1].data[index];
+                                                               b = image->comps[2].data[index];
+                                                               if (image->comps[0].sgnd){                      
+                                                                       r += adjust;
+                                                                       g += adjust;
+                                                                       b += adjust;
+                                                               }
+                                                               dat8[i+0] = r ; // R 
+                                                               if(i+1 <ssize) dat8[i+1] = g ;  else break;// G 
+                                                               if(i+2 <ssize) dat8[i+2] = b ;  else break;// B 
+                                                               index++;
+                                                       }else
+                                                               break;
+                                               }
                                        }
                                }else if (image->comps[0].prec == 12){
-                                       for (i=0; i<TIFFStripSize(tif); i+=9) { // 12 bits per pixel 
-                                               dat8[i+0] = (image->comps[0].data[index]>>8)<<4 | (image->comps[0].data[index]>>4);
-                                               dat8[i+1] = (image->comps[0].data[index]<<4)|((image->comps[1].data[index]>>8)& 0x0f);
-                                               dat8[i+2] = (image->comps[1].data[index]);
-                                               dat8[i+3] = (image->comps[2].data[index]>>8)<<4 | (image->comps[2].data[index]>>4);
-                                               dat8[i+4] = (image->comps[2].data[index]<<4)|((image->comps[0].data[index+1]>>8)& 0x0f);
-                                               dat8[i+5] = (image->comps[0].data[index+1]);
-                                               dat8[i+6] = (image->comps[1].data[index+1]>>8)<<4 | (image->comps[1].data[index+1]>>4);
-                                               dat8[i+7] = (image->comps[1].data[index+1]<<4)|((image->comps[2].data[index+1]>>8)& 0x0f);
-                                               dat8[i+8] = (image->comps[2].data[index+1]);
-                                               index+=2;
+                                       for (i=0; i<ssize-8; i+=9) {    // 12 bits per pixel 
+                                               int r = 0,g = 0,b = 0;
+                                               int r1 = 0,g1 = 0,b1 = 0;
+                                               if((index < imgsize)&(index+1 < imgsize)){
+                                                       r  = image->comps[0].data[index];
+                                                       g  = image->comps[1].data[index];
+                                                       b  = image->comps[2].data[index];
+                                                       r1 = image->comps[0].data[index+1];
+                                                       g1 = image->comps[1].data[index+1];
+                                                       b1 = image->comps[2].data[index+1];
+                                                       if (image->comps[0].sgnd){                                                                                                              
+                                                               r  += adjust;
+                                                               g  += adjust;
+                                                               b  += adjust;
+                                                               r1 += adjust;
+                                                               g1 += adjust;
+                                                               b1 += adjust;
+                                                       }
+                                                       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;
+                                                       last_i = i+9;
+                                               }else
+                                                       break;
+                                       }
+                                       if(last_i < ssize){
+                                               for (i= last_i; i<ssize; i+=9) {        // 12 bits per pixel 
+                                                       int r = 0,g = 0,b = 0;
+                                                       int r1 = 0,g1 = 0,b1 = 0;
+                                                       if((index < imgsize)&(index+1 < imgsize)){
+                                                               r  = image->comps[0].data[index];
+                                                               g  = image->comps[1].data[index];
+                                                               b  = image->comps[2].data[index];
+                                                               r1 = image->comps[0].data[index+1];
+                                                               g1 = image->comps[1].data[index+1];
+                                                               b1 = image->comps[2].data[index+1];
+                                                               if (image->comps[0].sgnd){                                                                                                              
+                                                                       r  += adjust;
+                                                                       g  += adjust;
+                                                                       b  += adjust;
+                                                                       r1 += adjust;
+                                                                       g1 += adjust;
+                                                                       b1 += adjust;
+                                                               }
+                                                               dat8[i+0] = (r >> 4);
+                                                               if(i+1 <ssize) dat8[i+1] = ((r & 0x0f) << 4 )|((g >> 8)& 0x0f); else break;
+                                                               if(i+2 <ssize) dat8[i+2] = g ;                  else break;
+                                                               if(i+3 <ssize) dat8[i+3] = (b >> 4);    else break;
+                                                               if(i+4 <ssize) dat8[i+4] = ((b & 0x0f) << 4 )|((r1 >> 8)& 0x0f);else break;
+                                                               if(i+5 <ssize) dat8[i+5] = r1;                  else break;
+                                                               if(i+6 <ssize) dat8[i+6] = (g1 >> 4);   else break;
+                                                               if(i+7 <ssize) dat8[i+7] = ((g1 & 0x0f)<< 4 )|((b1 >> 8)& 0x0f);else break;
+                                                               if(i+8 <ssize) dat8[i+8] = b1;                  else break;
+                                                               index+=2;
+                                                       }else
+                                                               break;
+                                               }
                                        }
                                }else if (image->comps[0].prec == 16){
-                                       for (i=0; i<TIFFStripSize(tif); i+=6) { // 16 bits per pixel 
-                                               dat8[i+0] =  image->comps[0].data[index];//LSB
-                                               dat8[i+1] = (image->comps[0].data[index]>> 8);//MSB      
-                                               dat8[i+2] =  image->comps[1].data[index]; 
-                                               dat8[i+3] = (image->comps[1].data[index]>> 8);  
-                                               dat8[i+4] =  image->comps[2].data[index];        
-                                               dat8[i+5] = (image->comps[2].data[index]>> 8); 
-                                               index++;
+                                       for (i=0 ; i<ssize-5 ; i+=6) {  // 16 bits per pixel 
+                                               int r = 0,g = 0,b = 0;
+                                               if(index < imgsize){
+                                                       r = image->comps[0].data[index];
+                                                       g = image->comps[1].data[index];
+                                                       b = image->comps[2].data[index];
+                                                       if (image->comps[0].sgnd){
+                                                       r += adjust;
+                                                       g += adjust;
+                                                       b += adjust;
+                                                       }
+                                                       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++;
+                                                       last_i = i+6;
+                                               }else
+                                                       break; 
+                                       }
+                                       if(last_i < ssize){
+                                               for (i=0 ; i<ssize ; i+=6) {    // 16 bits per pixel 
+                                                       int r = 0,g = 0,b = 0;
+                                                       if(index < imgsize){
+                                                               r = image->comps[0].data[index];
+                                                               g = image->comps[1].data[index];
+                                                               b = image->comps[2].data[index];
+                                                               if (image->comps[0].sgnd){
+                                                                       r += adjust;
+                                                                       g += adjust;
+                                                                       b += adjust;
+                                                               }
+                                                               dat8[i+0] =  r;//LSB
+                                                               if(i+1 <ssize) dat8[i+1] = (r >> 8);else break;//MSB     
+                                                               if(i+2 <ssize) dat8[i+2] =  g;          else break;
+                                                               if(i+3 <ssize) dat8[i+3] = (g >> 8);else break;
+                                                               if(i+4 <ssize) dat8[i+4] =  b;          else break;
+                                                               if(i+5 <ssize) dat8[i+5] = (b >> 8);else break;
+                                                               index++;
+                                                       }else
+                                                               break; 
+                                               }                                               
                                        }
                                }else{
                                        fprintf(stderr,"Bits=%d, Only 8,12,16 bits implemented\n",image->comps[0].prec);
@@ -1520,7 +1637,8 @@ int imagetotif(opj_image_t * image, const char *outfile) {
                        }
 
                        width   = image->comps[0].w;
-                       height= image->comps[0].h;
+                       height  = image->comps[0].h;
+                       imgsize = width * height;
                        bps             = image->comps[0].prec;
 
                        /* Set tags */
@@ -1534,34 +1652,59 @@ int imagetotif(opj_image_t * image, const char *outfile) {
                        TIFFSetField(tif, TIFFTAG_ROWSPERSTRIP, 1);
 
                        /* Get a buffer for the data */
-                       buf = _TIFFmalloc(TIFFStripSize(tif));
-                       index = 0;
-                       strip_size = 0;
                        strip_size = TIFFStripSize(tif);
+                       buf = _TIFFmalloc(strip_size);
+                       index = 0;                      
                        for (strip = 0; strip < TIFFNumberOfStrips(tif); strip++) {
                                unsigned char *dat8;
                                int i;
                                dat8 = buf;
                                if (image->comps[0].prec == 8){
                                        for (i=0; i<TIFFStripSize(tif); i+=1) { // 8 bits per pixel 
-                                               dat8[i+0] = image->comps[0].data[index] ;
-                                               index++;
+                                               if(index < imgsize){
+                                                       int r = 0;
+                                                       r = image->comps[0].data[index];
+                                                       if (image->comps[0].sgnd){
+                                                               r  += adjust;
+                                                       }
+                                                       dat8[i+0] = r;
+                                                       index++;
+                                               }else
+                                                       break; 
                                        }
                                }else if (image->comps[0].prec == 12){
                                        for (i = 0; i<TIFFStripSize(tif); i+=3) {       // 12 bits per pixel 
-                                               dat8[i+0] = (image->comps[0].data[index]>>8)<<4 | (image->comps[0].data[index]>>4);
-                                               dat8[i+1] = (image->comps[0].data[index]<<4)|((image->comps[0].data[index+1]>>8)& 0x0f);
-                                               dat8[i+2] = (image->comps[0].data[index+1]);
-                                               index+=2;
+                                               if(index < imgsize){
+                                                       int r = 0, r1 = 0;
+                                                       r  = image->comps[0].data[index];
+                                                       r1 = image->comps[0].data[index+1];
+                                                       if (image->comps[0].sgnd){
+                                                               r  += adjust;
+                                                               r1 += adjust;
+                                                       }
+                                                       dat8[i+0] = (r >> 4);
+                                                       dat8[i+1] = ((r & 0x0f) << 4 )|((r1 >> 8)& 0x0f);
+                                                       dat8[i+2] = r1 ;
+                                                       index+=2;
+                                               }else
+                                                       break; 
                                        }
                                }else if (image->comps[0].prec == 16){
                                        for (i=0; i<TIFFStripSize(tif); i+=2) { // 16 bits per pixel 
-                                               dat8[i+0] =  image->comps[0].data[index];
-                                               dat8[i+1] = (image->comps[0].data[index]>> 8);
-                                               index++;
+                                               if(index < imgsize){
+                                                       int r = 0;
+                                                       r = image->comps[0].data[index];
+                                                       if (image->comps[0].sgnd){
+                                                               r  += adjust;
+                                                       }
+                                                       dat8[i+0] = r;
+                                                       dat8[i+1] = r >> 8;
+                                                       index++;
+                                               }else
+                                                       break; 
                                        }
                                }else{
-                                       fprintf(stderr,"Bits=%d, Only 8,12,16 bits implemented\n",image->comps[0].prec);
+                                       fprintf(stderr,"TIFF file creation. Bits=%d, Only 8,12,16 bits implemented\n",image->comps[0].prec);
                                        fprintf(stderr,"Aborting\n");
                                        return 1;
                                }
@@ -1570,7 +1713,7 @@ int imagetotif(opj_image_t * image, const char *outfile) {
                        _TIFFfree(buf);
                        TIFFClose(tif);
                }else{
-                       fprintf(stderr,"False color format. Only RGB & Grayscale has been implemented\n");
+                       fprintf(stderr,"TIFF file creation. Bad color format. Only RGB & Grayscale has been implemented\n");
                        fprintf(stderr,"Aborting\n");
                        return 1;
                }
@@ -1590,6 +1733,7 @@ opj_image_t* tiftoimage(const char *filename, opj_cparameters_t *parameters)
        OPJ_COLOR_SPACE color_space;
        opj_image_cmptparm_t cmptparm[3];
        opj_image_t * image = NULL;
+       int imgsize = 0;
 
        tif = TIFFOpen(filename, "r");
 
@@ -1649,6 +1793,7 @@ opj_image_t* tiftoimage(const 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;
@@ -1658,43 +1803,52 @@ opj_image_t* tiftoimage(const char *filename, opj_cparameters_t *parameters)
 
                        if (Info.tiBps==12){
                                for (i=0; i<ssize; i+=9) {      /*12 bits per pixel*/
-                                       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);
-                                       image->comps[0].data[index+1] = ((dat8[i+4]& 0x0f)<< 8) | dat8[i+5];
-                                       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;
+                                       if((index < imgsize)&(index+1 < 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);
+                                               image->comps[0].data[index+1] = ((dat8[i+4]& 0x0f)<< 8) | dat8[i+5];
+                                               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 */
-                                       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++;
+                                       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 */
-                                       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 
-                                       if(parameters->cp_cinema){/* Rounding to 12 bits*/
-                                               image->comps[0].data[index] = image->comps[0].data[index] << 4 ;
-                                               image->comps[1].data[index] = image->comps[1].data[index] << 4 ;
-                                               image->comps[2].data[index] = image->comps[2].data[index] << 4 ;
-                                       }
-                                       index++;
+                                       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 
+                                               if(parameters->cp_cinema){/* Rounding to 12 bits*/
+                                                       image->comps[0].data[index] = image->comps[0].data[index] << 4 ;
+                                                       image->comps[1].data[index] = image->comps[1].data[index] << 4 ;
+                                                       image->comps[2].data[index] = image->comps[2].data[index] << 4 ;
+                                               }
+                                               index++;
+                                       }else
+                                               break;
                                }
                        }
                        else{
-                               fprintf(stderr,"Bits=%d, Only 8,12,16 bits implemented\n",Info.tiBps);
+                               fprintf(stderr,"TIFF file creation. Bits=%d, Only 8,12,16 bits implemented\n",Info.tiBps);
                                fprintf(stderr,"Aborting\n");
                                return NULL;
                        }
@@ -1735,6 +1889,7 @@ opj_image_t* tiftoimage(const 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;
@@ -1744,25 +1899,34 @@ opj_image_t* tiftoimage(const char *filename, opj_cparameters_t *parameters)
 
                        if (Info.tiBps==12){
                                for (i=0; i<ssize; i+=3) {      /* 12 bits per pixel*/
-                                       image->comps[0].data[index] = ( dat8[i+0]<<4 )                          |(dat8[i+1]>>4) ;
-                                       image->comps[0].data[index] = ((dat8[i+1]& 0x0f)<< 8)   | dat8[i+2];
-                                       index+=2;
+                                       if(index < imgsize){
+                                               image->comps[0].data[index]   = ( dat8[i+0]<<4 )                |(dat8[i+1]>>4) ;
+                                               image->comps[0].data[index+1] = ((dat8[i+1]& 0x0f)<< 8) | dat8[i+2];
+                                               index+=2;
+                                       }else
+                                               break;
                                }
                        }
                        else if( Info.tiBps==16){
                                for (i=0; i<ssize; i+=2) {      /* 16 bits per pixel */
-                                       image->comps[0].data[index] = ( dat8[i+1] << 8 ) | dat8[i+0];
-                                       index++;
+                                       if(index < imgsize){
+                                               image->comps[0].data[index] = ( dat8[i+1] << 8 ) | dat8[i+0];
+                                               index++;
+                                       }else
+                                               break;
                                }
                        }
                        else if ( Info.tiBps==8){
                                for (i=0; i<ssize; i+=1) {      /* 8 bits per pixel */
-                                       image->comps[0].data[index] = dat8[i+0];
-                                       index++;
+                                       if(index < imgsize){
+                                               image->comps[0].data[index] = dat8[i+0];
+                                               index++;
+                                       }else
+                                               break;
                                }
                        }
                        else{
-                               fprintf(stderr,"Bits=%d, Only 8,12,16 bits implemented\n",Info.tiBps);
+                               fprintf(stderr,"TIFF file creation. Bits=%d, Only 8,12,16 bits implemented\n",Info.tiBps);
                                fprintf(stderr,"Aborting\n");
                                return NULL;
                        }
@@ -1771,7 +1935,7 @@ opj_image_t* tiftoimage(const char *filename, opj_cparameters_t *parameters)
                _TIFFfree(buf);
                TIFFClose(tif);
        }else{
-               fprintf(stderr,"False color format. Only RGB & Grayscale has been implemented\n");
+               fprintf(stderr,"TIFF file creation. Bad color format. Only RGB & Grayscale has been implemented\n");
                fprintf(stderr,"Aborting\n");
                return NULL;
        }
@@ -1795,7 +1959,7 @@ opj_image_t* rawtoimage(const char *filename, opj_cparameters_t *parameters, raw
        opj_image_t * image = NULL;
        unsigned short ch;
        
-       if((raw_cp->rawWidth * raw_cp->rawHeight * raw_cp->rawComp * raw_cp->rawBitDepth) == 0)
+       if((! (raw_cp->rawWidth & raw_cp->rawHeight & raw_cp->rawComp & raw_cp->rawBitDepth)) == 0)
        {
                fprintf(stderr,"\nError: invalid raw image parameters\n");
                fprintf(stderr,"Please use the Format option -F:\n");
@@ -1834,7 +1998,6 @@ opj_image_t* rawtoimage(const char *filename, opj_cparameters_t *parameters, raw
                fclose(f);
                return NULL;
        }
-
        /* set image offset and reference grid */
        image->x0 = parameters->image_offset_x0;
        image->y0 = parameters->image_offset_y0;
@@ -1854,24 +2017,34 @@ opj_image_t* rawtoimage(const char *filename, opj_cparameters_t *parameters, raw
                        }
                }
        }
-       else
+       else if(raw_cp->rawBitDepth <= 16)
        {
-               unsigned short value = 0;
+               unsigned short value;
                for(compno = 0; compno < numcomps; compno++) {
                        for (i = 0; i < w * h; i++) {
-                               if (!fread(&value, 2, 1, f)) {
+                               unsigned char temp;
+                               if (!fread(&temp, 1, 1, f)) {
                                        fprintf(stderr,"Error reading raw file. End of file probably reached.\n");
                                        return NULL;
                                }
+                               value = temp << 8;
+                               if (!fread(&temp, 1, 1, f)) {
+                                       fprintf(stderr,"Error reading raw file. End of file probably reached.\n");
+                                       return NULL;
+                               }
+                               value += temp;
                                image->comps[compno].data[i] = raw_cp->rawSigned?(short)value:value;
                        }
                }
        }
+       else {
+               fprintf(stderr,"OpenJPEG cannot encode raw components with bit depth higher than 16 bits.\n");
+               return NULL;
+       }
 
        if (fread(&ch, 1, 1, f)) {
                fprintf(stderr,"Warning. End of raw file not reached... processing anyway\n");
        }
-
        fclose(f);
 
        return image;
@@ -1880,7 +2053,10 @@ opj_image_t* rawtoimage(const char *filename, opj_cparameters_t *parameters, raw
 int imagetoraw(opj_image_t * image, const char *outfile)
 {
        FILE *rawFile = NULL;
-       int compno, pixelsToWrite, offset, cont;
+       int compno;
+       int w, h;
+       int line, row;
+       int *ptr;
 
        if((image->numcomps * image->x1 * image->y1) == 0)
        {
@@ -1901,8 +2077,8 @@ int imagetoraw(opj_image_t * image, const char *outfile)
                fprintf(stdout,"Component %d characteristics: %dx%dx%d %s\n", compno, image->comps[compno].w,
                        image->comps[compno].h, image->comps[compno].prec, image->comps[compno].sgnd==1 ? "signed": "unsigned");
 
-               pixelsToWrite = image->comps[compno].w * image->comps[compno].h;
-               offset = 0;
+               w = image->comps[compno].w;
+               h = image->comps[compno].h;
 
                if(image->comps[compno].prec <= 8)
                {
@@ -1910,20 +2086,26 @@ int imagetoraw(opj_image_t * image, const char *outfile)
                        {
                                signed char curr;
                                int mask = (1 << image->comps[compno].prec) - 1;
-                               for(cont = 0; cont < pixelsToWrite; cont++)
-                               {                               
-                                       curr = (signed char) (image->comps[compno].data[cont] & mask);
-                                       fwrite(&curr, sizeof(signed char), 1, rawFile);
+                               ptr = image->comps[compno].data;
+                               for (line = 0; line < h; line++) {
+                                       for(row = 0; row < w; row++)    {                               
+                                               curr = (signed char) (*ptr & mask);
+                                               fwrite(&curr, sizeof(signed char), 1, rawFile);
+                                               ptr++;
+                                       }
                                }
                        }
                        else if(image->comps[compno].sgnd == 0)
                        {
                                unsigned char curr;
                                int mask = (1 << image->comps[compno].prec) - 1;
-                               for(cont = 0; cont < pixelsToWrite; cont++)
-                               {                               
-                                       curr = (unsigned char) (image->comps[compno].data[cont] & mask);
-                                       fwrite(&curr, sizeof(unsigned char), 1, rawFile);
+                               ptr = image->comps[compno].data;
+                               for (line = 0; line < h; line++) {
+                                       for(row = 0; row < w; row++)    {       
+                                               curr = (unsigned char) (*ptr & mask);
+                                               fwrite(&curr, sizeof(unsigned char), 1, rawFile);
+                                               ptr++;
+                                       }
                                }
                        }
                }
@@ -1933,31 +2115,45 @@ int imagetoraw(opj_image_t * image, const char *outfile)
                        {
                                signed short int curr;
                                int mask = (1 << image->comps[compno].prec) - 1;
-                               for(cont = 0; cont < pixelsToWrite; cont++)
-                               {                               
-                                       curr = (signed short int) (image->comps[compno].data[cont] & mask);
-                                       fwrite(&curr, sizeof(signed short int), 1, rawFile);
+                               ptr = image->comps[compno].data;
+                               for (line = 0; line < h; line++) {
+                                       for(row = 0; row < w; row++)    {                                       
+                                               unsigned char temp;
+                                               curr = (signed short int) (*ptr & mask);                                                
+                                               temp = curr >> 8;
+                                               fwrite(&temp, 1, 1, rawFile);
+                                               temp = curr;
+                                               fwrite(&temp, 1, 1, rawFile);
+                                               ptr++;
+                                       }
                                }
                        }
                        else if(image->comps[compno].sgnd == 0)
                        {
                                unsigned short int curr;
                                int mask = (1 << image->comps[compno].prec) - 1;
-                               for(cont = 0; cont < pixelsToWrite; cont++)
-                               {                               
-                                       curr = (unsigned short int) (image->comps[compno].data[cont] & mask);
-                                       fwrite(&curr, sizeof(unsigned short int), 1, rawFile);
+                               ptr = image->comps[compno].data;
+                               for (line = 0; line < h; line++) {
+                                       for(row = 0; row < w; row++)    {                               
+                                               unsigned char temp;
+                                               curr = (unsigned short int) (*ptr & mask);                                              
+                                               temp = curr >> 8;
+                                               fwrite(&temp, 1, 1, rawFile);
+                                               temp = curr;
+                                               fwrite(&temp, 1, 1, rawFile);
+                                               ptr++;
+                                       }
                                }
                        }
                }
                else if (image->comps[compno].prec <= 32)
                {
-
-
+                       fprintf(stderr,"More than 16 bits per component no handled yet\n");
+                       return 1;
                }
                else
                {
-                       fprintf(stderr,"\nError: invalid precision\n");
+                       fprintf(stderr,"Error: invalid precision: %d\n", image->comps[compno].prec);
                        return 1;
                }
        }