imageio.c

Go to the documentation of this file.

#include <string.h>
#include <ctype.h>
#include "imageio.h"

#ifdef USE_LIBPNG
#include <zlib.h>
#include <png.h>
#if PNG_LIBPNG_VER < 10400
/* For compatibility with older libpng */
#define png_set_expand_gray_1_2_4_to_8  png_set_gray_1_2_4_to_8
#endif
#endif
#ifdef USE_LIBTIFF
#include <tiffio.h>
#endif
#ifdef USE_LIBJPEG
#include <jpeglib.h>
#include <setjmp.h>
#endif

#define FILE_BUFFER_CAPACITY    (1024*4)

#define ROUNDCLAMPF(x)   ((x < 0.0f) ? 0 : \
    ((x > 1.0f) ? 255 : (uint8_t)(255.0f*(x) + 0.5f)))
#define ROUNDCLAMP(x)   ((x < 0.0) ? 0 : \
    ((x > 1.0) ? 255 : (uint8_t)(255.0*(x) + 0.5)))


static int string_ends_with(const char *string, const char *suffix)
{
    unsigned string_length = strlen(string), suffix_length = strlen(suffix);
    unsigned i;
    
    if(string_length < suffix_length)
        return 0;
    
    string += string_length - suffix_length;
    
    for(i = 0; i < suffix_length; i++)
        if(tolower(string[i]) != tolower(suffix[i]))
            return 0;
    
    return 1;
}


static void fill_image(uint32_t *image, int width, int height, uint32_t color)
{
    int x, y;
    
    if(image)
        for(y = 0; y < height; y++, image += width)
            for(x = 0; x < width; x++)
                image[x] = color;
}


static uint32_t *get_image_palette(int *num_colors, int *use_color,
    int *use_alpha, const uint32_t *image, int width, int height)
{
    const int max_colors = 256;
    uint32_t *palette = NULL;
    uint32_t pixel;
    int x, y, i, red, green, blue, alpha;
    
    
    if(!use_color || !num_colors || !use_alpha)
        return NULL;
    else if(!image
        || !(palette = (uint32_t *)Malloc(sizeof(uint32_t)*max_colors)))
    {
        *num_colors = -1;
        *use_color = *use_alpha = 1;
        return NULL;
    }
    
    *num_colors = *use_color = *use_alpha = 0;
   
    for(y = 0; y < height; y++)
    {
        for(x = 0; x < width; x++)
        {
            pixel = *(image++);
            red = ((uint8_t *)&pixel)[0];
            green = ((uint8_t *)&pixel)[1];
            blue = ((uint8_t *)&pixel)[2];
            alpha = ((uint8_t *)&pixel)[3];
                        
            if(red != green || red != blue)     /* Check color */
                *use_color = 1;
            
            if(alpha != 255)                    /* Check alpha */
                *use_alpha = 1;
            
            /* Check palette colors (if *num_colors != -1) */
            for(i = 0; i < *num_colors; i++)
                if(pixel == palette[i])
                    break;
            
            if(i == *num_colors)
            {
                if(i < max_colors)
                {   /* Add new color to palette */
                    palette[i] = pixel;
                    (*num_colors)++;
                }
                else
                {   /* Maximum size for palette exceeded */
                    Free(palette);
                    palette = NULL;
                    *num_colors = -1;    /* Don't check palette colors */
                }
            }
        }
    }
    
    return palette;
}


static uint16_t read_u16_le(FILE *file)
{
    uint16_t w;
    w = (uint16_t) getc(file);
    w |= ((uint16_t) getc(file) << 8);
    return w;
}


static uint32_t read_u32_le(FILE *file)
{
    uint32_t dw;
    dw = (uint32_t) getc(file);
    dw |= ((uint32_t) getc(file) << 8);
    dw |= ((uint32_t) getc(file) << 16);
    dw |= ((uint32_t) getc(file) << 24);
    return dw;
}


static void write_u16_le(uint16_t w, FILE *file)
{
    putc(w & 0xFF, file);
    putc((w & 0xFF00) >> 8, file);
}


static void write_u32_le(uint32_t dw, FILE *file)
{
    putc(dw & 0xFF, file);
    putc((dw & 0xFF00) >> 8, file);
    putc((dw & 0xFF0000) >> 16, file);
    putc((dw & 0xFF000000) >> 24, file);
}


static int read_bmp_1bpp(uint32_t *image, int width, int height,
    FILE *file, const uint32_t *palette)
{
    int row_padding = (-(width+7)/8)&3;
    int x, y, bit;
    unsigned code;
    
    image += ((long int)width)*((long int)height - 1);
    
    for(y = height; y; y--, image -= width)
    {
        if(feof(file))
            return 0;
        
        for(x = 0; x < width;)
        {
            code = getc(file);
            
            for(bit = 7; bit >= 0 && x < width; bit--, code <<= 1)
                image[x++] = palette[(code & 0x80) ? 1:0];
        }
        
        for(x = row_padding; x; x--)
            getc(file); /* Skip padding bytes at the end of the row */
    }
    
    return 1;
}


static int read_bmp_4bpp(uint32_t *image, int width, int height,
    FILE *file, const uint32_t *palette)
{
    int row_padding = (-(width+1)/2)&3;
    int x, y;
    unsigned code;

    image += ((long int)width)*((long int)height - 1);
    
    for(y = height; y; y--, image -= width)
    {
        if(feof(file))
            return 0;
        
        for(x = 0; x < width;)
        {
            code = getc(file);
            image[x++] = palette[(code & 0xF0) >> 4];
            
            if(x < width)
                image[x++] = palette[code & 0x0F];
        }
        
        for(x = row_padding; x; x--)
            getc(file); /* Skip padding bytes at the end of the row */
    }
    
    return 1;
}


static int read_bmp_4bpp_rle(uint32_t *image, int width, int height,
    FILE *file, const uint32_t *palette)
{
    int x, y, dy, k;
    unsigned count, value;
    uint32_t color_high, color_low;
    
    fill_image(image, width, height, palette[0]);
    image += ((long int)width)*((long int)height - 1);
    
    for(x = 0, y = height; y;)
    {
        if(feof(file))
            return 0;
        
        count = getc(file);
        value = getc(file);
        
        if(!count)
        {   /* count = 0 is the escape code */
            switch(value)
            {
            case 0:     /* End of line */
                image -= width;
                x = 0;
                y--;
                break;
            case 1:     /* End of bitmap */
                return 1;
            case 2:     /* Delta */
                x += getc(file);
                dy = getc(file);
                y -= dy;
                image -= dy*width;
                
                if(x >= width || y < 0)
                    return 0;
                break;
            default:
                /* Read a run of uncompressed data (value = length of run) */
                count = k = value;
                
                if(x >= width)
                    return 0;

                do
                {
                    value = getc(file);
                    image[x++] = palette[(value & 0xF0) >> 4];
                    
                    if(x >= width)
                        break;
                        
                    if(--k)
                    {
                        image[x++] = palette[value & 0x0F];
                        k--;
                        
                        if(x >= width)
                            break;
                    }
                }while(k);
                
                if(((count + 1)/2) & 1)
                    getc(file); /* Padding for word align */
            }
        }
        else
        {   /* Run of pixels (count = length of run) */
            color_high = palette[(value & 0xF0) >> 4];
            color_low = palette[value & 0xF];
            
            if(x >= width)
                return 0;
            
            do
            {
                image[x++] = color_high;
                count--;
                
                if(x >= width)
                    break;
                
                if(count)
                {
                    image[x++] = color_low;
                    count--;
                    
                    if(x >= width)
                        break;
                }
            }while(count);
        }
    }
    
    return 1;
}


static int read_bmp_8bpp(uint32_t *image, int width, int height,
    FILE *file, const uint32_t *palette)
{
    int row_padding = (-width)&3;
    int x, y;
    
    image += ((long int)width)*((long int)height - 1);
    
    for(y = height; y; y--, image -= width)
    {
        if(feof(file))
            return 0;
        
        for(x = 0; x < width; x++)
            image[x] = palette[getc(file) & 0xFF];
        
        for(x = row_padding; x; x--)
            getc(file); /* Skip padding bytes at the end of the row */
    }
    
    return 1;
}


static int read_bmp_8bpp_rle(uint32_t *image, int width, int height,
    FILE *file, const uint32_t *palette)
{
    int x, y, dy, k;
    unsigned count, value;
    uint32_t color;
    
    fill_image(image, width, height, palette[0]);
    image += ((long int)width)*((long int)height - 1);
    
    for(x = 0, y = height; y;)
    {
        if(feof(file))
            return 0;
            
        count = getc(file);
        value = getc(file);
        
        if(!count)
        {   /* count = 0 is the escape code */
            switch(value)
            {
            case 0:     /* End of line */
                image -= width;
                x = 0;
                y--;
                break;
            case 1:     /* End of bitmap */
                return 1;
            case 2:     /* Delta */
                x += getc(file);
                dy = getc(file);
                y -= dy;
                image -= dy*width;
                
                if(x >= width || y < 0)
                    return 0;
                break;
            default:
                /* Read a run of uncompressed data (value = length of run) */
                count = k = value;
                
                do
                {
                    if(x >= width)
                        break;
                    
                    image[x++] = palette[getc(file) & 0xFF];
                }while(--k);
                
                if(count&1)
                    getc(file); /* Padding for word align */
            }
        }
        else
        {   /* Run of pixels equal to value (count = length of run) */
            color = palette[value & 0xFF];
            
            do
            {
                if(x >= width)
                    break;
                
                image[x++] = color;
            }while(--count);
        }
    }
    
    return 1;
}


static int read_bmp_24bpp(uint32_t *image, int width, int height, FILE *file)
{
    uint8_t *image_ptr = (uint8_t *)image;
    int row_padding = (-3*width)&3;
    int x, y;
    
    width <<= 2;
    image_ptr += ((long int)width)*((long int)height - 1);
    
    for(y = height; y; y--, image_ptr -= width)
    {
        if(feof(file))
            return 0;
        
        for(x = 0; x < width; x += 4)
        {
            image_ptr[x+3] = 255;        /* Set alpha            */
            image_ptr[x+2] = getc(file); /* Read blue component  */
            image_ptr[x+1] = getc(file); /* Read green component */
            image_ptr[x+0] = getc(file); /* Read red component   */
        }
        
        for(x = row_padding; x; x--)
            getc(file); /* Skip padding bytes at the end of the row */
    }
    
    return 1;
}

static void get_mask_shifts(uint32_t mask, int *left_shift, int *right_shift)
{
    int shift = 0, bitcount = 0;
    
    if(!mask)
    {
        *left_shift = 0;
        *right_shift = 0;
        return;
    }
    
    while(!(mask & 1))  /* Find the first true bit */
    {
        mask >>= 1;
        ++shift;
    }
    
    /* Adjust the result for scaling to 8-bit quantities */
    while(mask & 1)     /* count the number of true bits */
    {
        mask >>= 1;
        ++bitcount;
    }
    
    /* Compute a signed shift (right is positive) */
    shift += bitcount - 8;
    
    if(shift >= 0)
    {
        *left_shift = 0;
        *right_shift = shift;
    }
    else
    {
        *left_shift = -shift;
        *right_shift = 0;
    }
}

static int read_bmp_16bpp(uint32_t *image, int width, int height, FILE *file,
    uint32_t redmask, uint32_t greenmask,
    uint32_t bluemask, uint32_t alphamask)
{
    uint8_t *image_ptr = (uint8_t *)image;
    uint32_t code;
    int row_padding = (-2*width)&3;
    int redleft_shift, greenleft_shift, blueleft_shift, alphaleft_shift;
    int redright_shift, greenright_shift, blueright_shift, alpharight_shift;
    int x, y;
    
    get_mask_shifts(redmask, &redleft_shift, &redright_shift);
    get_mask_shifts(greenmask, &greenleft_shift, &greenright_shift);
    get_mask_shifts(bluemask, &blueleft_shift, &blueright_shift);
    get_mask_shifts(alphamask, &alphaleft_shift, &alpharight_shift);
    width <<= 2;
    image_ptr += ((long int)width)*((long int)height - 1);
    
    for(y = height; y; y--, image_ptr -= width)
    {
        if(feof(file))
            return 0;
        
        for(x = 0; x < width; x += 4)
        {
            code = read_u16_le(file);
            /* By the Windows 4.x BMP specification, masks must be contiguous
               <http://www.fileformat.info/format/bmp/egff.htm>.  So we can
               decode bitfields by bitshifting and bitwise AND. */
            image_ptr[x + 3] = ((code & alphamask) >> alpharight_shift)
                << alphaleft_shift;
            image_ptr[x + 2] = ((code & bluemask ) >> blueright_shift )
                << blueleft_shift;
            image_ptr[x + 1] = ((code & greenmask) >> greenright_shift)
                << greenleft_shift;
            image_ptr[x + 0] = ((code & redmask  ) >> redright_shift  )
                << redleft_shift;
        }
        
        for(x = row_padding; x; x--)
            getc(file); /* Skip padding bytes at the end of the row */
    }

    return 1;
}


static int read_bmp_32bpp(uint32_t *image, int width, int height, FILE *file,
    uint32_t redmask, uint32_t greenmask,
    uint32_t bluemask, uint32_t alphamask)
{
    uint8_t *image_ptr;
    uint32_t code;
    int redleft_shift, greenleft_shift, blueleft_shift, alphaleft_shift;
    int redright_shift, greenright_shift, blueright_shift, alpharight_shift;
    int x, y;
    
    get_mask_shifts(redmask, &redleft_shift, &redright_shift);
    get_mask_shifts(greenmask, &greenleft_shift, &greenright_shift);
    get_mask_shifts(bluemask, &blueleft_shift, &blueright_shift);
    get_mask_shifts(alphamask, &alphaleft_shift, &alpharight_shift);
    width <<= 2;
    image_ptr = (uint8_t *)image + ((long int)width)*((long int)height - 1);
    
    for(y = height; y; y--, image_ptr -= width)
    {
        if(feof(file))
            return 0;
        
        for(x = 0; x < width; x += 4)
        {
            code = read_u32_le(file);
            image_ptr[x + 3] = ((code & alphamask) >> alpharight_shift)
                << alphaleft_shift;
            image_ptr[x + 2] = ((code & bluemask ) >> blueright_shift )
                << blueleft_shift;
            image_ptr[x + 1] = ((code & greenmask) >> greenright_shift)
                << greenleft_shift;
            image_ptr[x + 0] = ((code & redmask  ) >> redright_shift  )
                << redleft_shift;
        }
    }

    return 1;
}

static int read_bmp(uint32_t **image, int *width, int *height, FILE *file)
{
    uint32_t *palette = NULL;
    uint8_t *palette_ptr;
    long int image_data_offset, info_size;
    unsigned i, num_planes, bits_per_pixel, compression, num_colors;
    uint32_t redmask, greenmask, bluemask, alphamask;
    int success = 0, os2bmp;
    uint8_t magic[2];
    
    *image = NULL;
    *width = *height = 0;
    fseek(file, 0, SEEK_SET);

    magic[0] = getc(file);
    magic[1] = getc(file);

    if(!(magic[0] == 0x42 && magic[1] == 0x4D) /* Verify the magic numbers */
        || fseek(file, 8, SEEK_CUR))         /* Skip the reserved fields */
    {
        ErrorMessage("Invalid BMP header.\n");
        goto fail;
    }
    
    image_data_offset = read_u32_le(file);
    info_size = read_u32_le(file);
    
    /* Read the info header */
    if(info_size < 12)
    {
        ErrorMessage("Invalid BMP info header.\n");
        goto fail;
    }
    
    if((os2bmp = (info_size == 12)))  /* This is an OS/2 V1 infoheader */
    {
        *width = (int)read_u16_le(file);
        *height = (int)read_u16_le(file);
        num_planes = (unsigned)read_u16_le(file);
        bits_per_pixel = (unsigned)read_u16_le(file);
        compression = 0;
        num_colors = 0;
        redmask = 0x00FF0000;
        greenmask = 0x0000FF00;
        bluemask = 0x000000FF;
        alphamask = 0xFF000000;
    }
    else
    {
        *width = abs((int)read_u32_le(file));
        *height = abs((int)read_u32_le(file));
        num_planes = (unsigned)read_u16_le(file);
        bits_per_pixel = (unsigned)read_u16_le(file);
        compression = (unsigned)read_u32_le(file);
        fseek(file, 12, SEEK_CUR);
        num_colors = (unsigned)read_u32_le(file);
        fseek(file, 4, SEEK_CUR);
        redmask = read_u32_le(file);
        greenmask = read_u32_le(file);
        bluemask = read_u32_le(file);
        alphamask = read_u32_le(file);
    }
    
    /* Check for problems or unsupported compression modes */
    if(*width > MAX_IMAGE_SIZE || *height > MAX_IMAGE_SIZE)
    {
        ErrorMessage("image dimensions exceed MAX_IMAGE_SIZE.\n");
        goto fail;
    }
    
    if(feof(file) || num_planes != 1 || compression > 3)
        goto fail;
    
    /* Allocate the image data */
    if(!(*image = (uint32_t *)Malloc(
        sizeof(uint32_t)*((long int)*width)*((long int)*height))))
        goto fail;
    
    /* Read palette */
    if(bits_per_pixel <= 8)
    {
        fseek(file, 14 + info_size, SEEK_SET);
        
        if(!num_colors)
            num_colors = 1 << bits_per_pixel;
        
        if(!(palette = (uint32_t *)Malloc(sizeof(uint32_t)*256)))
            goto fail;
        
        for(i = 0, palette_ptr = (uint8_t *)palette; i < num_colors; i++)
        {
            palette_ptr[3] = 255;          /* Set alpha            */
            palette_ptr[2] = getc(file);   /* Read blue component  */
            palette_ptr[1] = getc(file);   /* Read green component */
            palette_ptr[0] = getc(file);   /* Read red component   */
            palette_ptr += 4;
            
            if(!os2bmp)
                getc(file); /* Skip extra byte (for non-OS/2 bitmaps) */
        }
        
        /* Fill the rest of the palette with the first color */
        for(; i < 256; i++)
            palette[i] = palette[0];
    }
    
    if(fseek(file, image_data_offset, SEEK_SET) || feof(file))
    {
        ErrorMessage("file error.\n");
        goto fail;
    }
    
    /*** Read the bitmap image data ***/
    switch(compression)
    {
        case 0: /* Uncompressed data */
            switch(bits_per_pixel)
            {
            case 1: /* Read 1-bit uncompressed indexed data */
                success = read_bmp_1bpp(
                    *image, *width, *height, file, palette);
                break;
            case 4: /* Read 4-bit uncompressed indexed data */
                success = read_bmp_4bpp(
                    *image, *width, *height, file, palette);
                break;
            case 8: /* Read 8-bit uncompressed indexed data */
                success = read_bmp_8bpp(
                    *image, *width, *height, file, palette);
                break;
            case 24: /* Read 24-bit BGR image data */
                success = read_bmp_24bpp(*image, *width, *height, file);
                break;
            case 16: /* Read 16-bit data */
                success = read_bmp_16bpp(*image, *width, *height, file,
                    0x001F << 10, 0x001F << 5, 0x0001F, 0);
                break;
            case 32: /* Read 32-bit BGRA image data */
                success = read_bmp_32bpp(*image, *width, *height, file,
                    0x00FF0000, 0x0000FF00, 0x000000FF, 0xFF000000);
                break;
            }
            break;
        case 1: /* 8-bit RLE */
            if(bits_per_pixel == 8)
                success = read_bmp_8bpp_rle(
                    *image, *width, *height, file, palette);
            break;
        case 2: /* 4-bit RLE */
            if(bits_per_pixel == 4)
                success = read_bmp_4bpp_rle(
                    *image, *width, *height, file, palette);
            break;
        case 3: /* bitfields data */
            switch(bits_per_pixel)
            {
            case 16: /* Read 16-bit bitfields data */
                success = read_bmp_16bpp(*image, *width, *height, file,
                    redmask, greenmask, bluemask, alphamask);
                break;
            case 32: /* Read 32-bit bitfields data */
                success = read_bmp_32bpp(*image, *width, *height, file,
                    redmask, greenmask, bluemask, alphamask);
                break;
            }
            break;
    }
    
    if(!success)
        ErrorMessage("Error reading BMP data.\n");
    
fail:   /* There was a problem, clean up and exit */
    if(palette)
        Free(palette);
    
    if(!success && *image)
        Free(*image);
    
    return success;
}


static int write_bmp(const uint32_t *image, int width, int height, FILE *file)
{
    const uint8_t *image_ptr = (uint8_t *)image;
    uint32_t *palette = NULL;
    uint32_t pixel;
    long int imageSize;
    int use_palette, num_colors, use_color, use_alpha;
    int x, y, i, row_padding, success = 0;

    
    if(!image)
        return 0;
    
    palette = get_image_palette(&num_colors, &use_color, &use_alpha,
        image, width, height);
    
    /* Decide whether to use 8-bit palette or 24-bit RGB format */
    if(palette && 2*num_colors < width*height)
        use_palette = 1;
    else
        use_palette = num_colors = 0;
    
    /* Tell file to use buffering */
    setvbuf(file, 0, _IOFBF, FILE_BUFFER_CAPACITY);
    
    if(use_palette)
    {
        row_padding = (-width)&3;
        imageSize = (width + row_padding)*((long int)height);
    }
    else
    {
        row_padding = (-3*width)&3;
        imageSize = (3*width + row_padding)*((long int)height);
    }
    
    /*** Write the header ***/
    
    /* Write the BMP header */
    putc(0x42, file);                       /* magic numbers             */
    putc(0x4D, file);
    
    /* filesize */
    write_u32_le(54 + 4*num_colors + imageSize, file);
    
    write_u32_le(0, file);                  /* Reserved fields           */
    write_u32_le(54 + 4*num_colors, file);   /* image data offset */
    
    /* Write the infoheader */
    write_u32_le(40, file);                 /* infoheader size           */
    write_u32_le(width, file);              /* image width               */
    write_u32_le(height, file);             /* image height              */
    write_u16_le(1, file);                   /* Number of colorplanes     */
    write_u16_le((use_palette) ? 8:24, file); /* bits per pixel */
    write_u32_le(0, file);                  /* compression method (none) */
    write_u32_le(imageSize, file);          /* image size                */
    write_u32_le(2835, file);               /* HResolution (2835=72dpi)  */
    write_u32_le(2835, file);               /* VResolution               */
    
    /* Number of colors */
    write_u32_le((!use_palette || num_colors == 256) ? 0:num_colors, file);
    
    write_u32_le(0, file);                  /* Important colors          */
    
    if(ferror(file))
    {
        ErrorMessage("Error during write to file.\n");
        goto fail;
    }
    
    if(use_palette)
    {   /* Write the palette */
        for(i = 0; i < num_colors; i++)
        {
            pixel = palette[i];
            putc(((uint8_t *)&pixel)[2], file);     /* blue   */
            putc(((uint8_t *)&pixel)[1], file);     /* green  */
            putc(((uint8_t *)&pixel)[0], file);     /* red    */
            putc(0, file);                          /* Unused */
        }
    }
    
    /* Write the image data */
    width <<= 2;
    image_ptr += ((long int)width)*((long int)height - 1);
            
    for(y = height; y; y--, image_ptr -= width)
    {
        if(use_palette)
        {   /* 8-bit palette image data */
            for(x = 0; x < width; x += 4)
            {
                pixel = *((uint32_t *)(image_ptr + x));
                
                for(i = 0; i < num_colors; i++)
                    if(pixel == palette[i])
                        break;
                
                putc(i, file);
            }
        }
        else
        {   /* 24-bit RGB image data */
            for(x = 0; x < width; x += 4)
            {
                putc(image_ptr[x+2], file);  /* Write blue component  */
                putc(image_ptr[x+1], file);  /* Write green component */
                putc(image_ptr[x+0], file);  /* Write red component   */
            }
        }
        
        for(x = row_padding; x; x--)         /* Write row padding */
            putc(0, file);
    }
    
    if(ferror(file))
    {
        ErrorMessage("Error during write to file.\n");
        goto fail;
    }
    
    success = 1;
fail:
    if(palette)
        Free(palette);
    return success;
}


#ifdef USE_LIBJPEG

typedef struct{
    struct jpeg_error_mgr pub;
    jmp_buf jmpbuf;
} hooked_jerr;


METHODDEF(void) jerr_exit(j_common_ptr cinfo)
{
    hooked_jerr *jerr = (hooked_jerr *) cinfo->err;
    (*cinfo->err->output_message)(cinfo);
    longjmp(jerr->jmpbuf, 1);
}


static int read_jpeg(uint32_t **image, int *width, int *height, FILE *file)
{
    struct jpeg_decompress_struct cinfo;
    hooked_jerr jerr;
    JSAMPARRAY buffer;
    uint8_t *image_ptr;
    unsigned i, row_size;
    
    *image = 0;
    *width = *height = 0;
    cinfo.err = jpeg_std_error(&jerr.pub);
    jerr.pub.error_exit = jerr_exit;
    
    if(setjmp(jerr.jmpbuf))
        goto fail;  /* If this code is reached, libjpeg has signaled an error. */
    
    jpeg_create_decompress(&cinfo);
    jpeg_stdio_src(&cinfo, file);
    jpeg_read_header(&cinfo, 1);
    cinfo.out_color_space = JCS_RGB;   /* Ask for RGB image data */
    jpeg_start_decompress(&cinfo);
    *width = (int)cinfo.output_width;
    *height = (int)cinfo.output_height;
    
    if(*width > MAX_IMAGE_SIZE || *height > MAX_IMAGE_SIZE)
    {
        ErrorMessage("image dimensions exceed MAX_IMAGE_SIZE.\n");
        jpeg_abort_decompress(&cinfo);
        goto fail;
    }
    
    /* Allocate image memory */
    if(!(*image = (uint32_t *)Malloc(sizeof(uint32_t)
        *((size_t)*width)*((size_t)*height))))
    {
        jpeg_abort_decompress(&cinfo);
        goto fail;
    }
    
    /* Allocate a one-row-high array that will go away when done */
    row_size = cinfo.output_width * cinfo.output_components;
    buffer = (*cinfo.mem->alloc_sarray) ((j_common_ptr) &cinfo,
        JPOOL_IMAGE, row_size, 1);
    image_ptr = (uint8_t *)*image;
    
    while(cinfo.output_scanline < cinfo.output_height)
        for(jpeg_read_scanlines(&cinfo, buffer, 1), i = 0;
            i < row_size; i += 3)
        {
            *(image_ptr++) = buffer[0][i];   /* red   */
            *(image_ptr++) = buffer[0][i+1]; /* green */
            *(image_ptr++) = buffer[0][i+2]; /* blue  */
            *(image_ptr++) = 0xFF;
        }

    jpeg_finish_decompress(&cinfo);
    jpeg_destroy_decompress(&cinfo);
    return 1;
    
fail:
    if(*image)
        Free(*image);
    
    *width = *height = 0;
    jpeg_destroy_decompress(&cinfo);
    return 0;
}


static int write_jpeg(const uint32_t *image, int width, int height,
    FILE *file, int quality)
{
    struct jpeg_compress_struct cinfo;
    hooked_jerr jerr;
    uint8_t *buffer = 0, *image_ptr;
    unsigned i, row_size;
    
    if(!image)
        return 0;
    
    cinfo.err = jpeg_std_error(&jerr.pub);
    jerr.pub.error_exit = jerr_exit;
    
    if(setjmp(jerr.jmpbuf))
        goto fail;  /* libjpeg has signaled an error */
    
    jpeg_create_compress(&cinfo);
    jpeg_stdio_dest(&cinfo, file);
    cinfo.image_width = width;
    cinfo.image_height = height;
    cinfo.input_components = 3;
    cinfo.in_color_space = JCS_RGB;
    jpeg_set_defaults(&cinfo);
    jpeg_set_quality(&cinfo, (quality < 100) ? quality : 100, 1);
    jpeg_start_compress(&cinfo, 1);

    row_size = 3*width;
    image_ptr = (uint8_t *)image;

    if(!(buffer = (uint8_t *)Malloc(row_size)))
        goto fail;
    
    while(cinfo.next_scanline < cinfo.image_height)
    {
        for(i = 0; i < row_size; i += 3)
        {
            buffer[i] = image_ptr[0];   /* red   */
            buffer[i+1] = image_ptr[1]; /* green */
            buffer[i+2] = image_ptr[2]; /* blue  */
            image_ptr += 4;
        }

        jpeg_write_scanlines(&cinfo, &buffer, 1);
    }

    if(buffer)
        Free(buffer);

    jpeg_finish_compress(&cinfo);
    jpeg_destroy_compress(&cinfo);
    return 1;
fail:
    if(buffer)
        Free(buffer);

    jpeg_destroy_compress(&cinfo);
    return 0;
}
#endif /* USE_LIBJPEG */


#ifdef USE_LIBPNG

static int read_png(uint32_t **image, int *width, int *height, FILE *file)
{
    png_bytep *row_pointers;
    png_byte header[8];
    png_structp png;
    png_infop info;
    png_uint_32 png_width, png_height;
    int bit_depth, color_type, interlace_type;
    unsigned row;
    
    *image = 0;
    *width = *height = 0;
    
    /* Check that file is a PNG file */
    if(fread(header, 1, 8, file) != 8 || png_sig_cmp(header, 0, 8))
        return 0;
    
    /* Read the info header */
    if(!(png = png_create_read_struct(
        PNG_LIBPNG_VER_STRING, NULL, NULL, NULL))
        || !(info = png_create_info_struct(png)))
    {
        if(png)
            png_destroy_read_struct(&png, (png_infopp)NULL, (png_infopp)NULL);
        
        return 0;
    }
        
    if(setjmp(png_jmpbuf(png)))
        goto fail; /* libpng has signaled an error */
    
    png_init_io(png, file);
    png_set_sig_bytes(png, 8);
    png_set_user_limits(png, MAX_IMAGE_SIZE, MAX_IMAGE_SIZE);
    png_read_info(png, info);
    png_get_IHDR(png, info, &png_width, &png_height, &bit_depth, &color_type,
        &interlace_type, (int*)NULL, (int*)NULL);
    *width = (int)png_width;
    *height = (int)png_height;
    
    /* Tell libpng to convert everything to 32-bit RGBA */
    if(color_type == PNG_COLOR_TYPE_PALETTE)
        png_set_palette_to_rgb(png);
    if(color_type == PNG_COLOR_TYPE_GRAY && bit_depth < 8)
        png_set_expand_gray_1_2_4_to_8(png);
    if(color_type == PNG_COLOR_TYPE_GRAY
        || color_type == PNG_COLOR_TYPE_GRAY_ALPHA)
        png_set_gray_to_rgb(png);
    if(png_get_valid(png, info, PNG_INFO_tRNS))
        png_set_tRNS_to_alpha(png);
    
    png_set_strip_16(png);
    png_set_filler(png, 0xFF, PNG_FILLER_AFTER);

    png_set_interlace_handling(png);
    png_read_update_info(png, info);
    
    /* Allocate image memory and row pointers */
    if(!(*image = (uint32_t *)Malloc(sizeof(uint32_t)
        *((size_t)*width)*((size_t)*height)))
        || !(row_pointers = (png_bytep *)Malloc(sizeof(png_bytep)
        *png_height)))
        goto fail;

    for(row = 0; row < png_height; row++)
        row_pointers[row] = (png_bytep)(*image + png_width*row);
    
    /* Read the image data */
    png_read_image(png, row_pointers);
    Free(row_pointers);
    png_destroy_read_struct(&png, &info, (png_infopp)NULL);
    return 1;
    
fail:
    if(*image)
        Free(*image);
    
    *width = *height = 0;
    png_destroy_read_struct(&png, &info, (png_infopp)NULL);
    return 0;
}


static int write_png(const uint32_t *image, int width, int height, FILE *file)
{
    const uint32_t *image_ptr;
    uint32_t *palette = NULL;
    uint8_t *row_buffer;
    png_structp png;
    png_infop info;
    png_color png_palette[256];
    png_byte png_trans[256];
    uint32_t pixel;
    int png_color_type, num_colors, use_color, use_alpha;
    int x, y, i, success = 0;

    
    if(!image)
        return 0;
    
    if(!(row_buffer = (uint8_t *)Malloc(4*width)))
        return 0;
    
    if(!(png = png_create_write_struct(PNG_LIBPNG_VER_STRING,
        NULL, NULL, NULL))
        || !(info = png_create_info_struct(png)))
    {
        if(png)
            png_destroy_write_struct(&png, (png_infopp)NULL);
    
        Free(row_buffer);
        return 0;
    }
        
    if(setjmp(png_jmpbuf(png)))
    {   /* If this code is reached, libpng has signaled an error. */
        goto fail;
    }

    /* Configure PNG output */
    png_init_io(png, file);
    png_set_compression_level(png, Z_BEST_COMPRESSION);
    
    palette = get_image_palette(&num_colors, &use_color, &use_alpha,
        image, width, height);
        
    /* The image is written according to the analysis of get_image_palette */
    if(palette && use_color)
        png_color_type = PNG_COLOR_TYPE_PALETTE;
    else if(use_alpha)
        png_color_type = PNG_COLOR_TYPE_RGB_ALPHA;
    else if(use_color)
        png_color_type = PNG_COLOR_TYPE_RGB;
    else
        png_color_type = PNG_COLOR_TYPE_GRAY;
    
    png_set_IHDR(png, info, width, height, 8, png_color_type,
        PNG_INTERLACE_NONE, PNG_COMPRESSION_TYPE_BASE, PNG_FILTER_TYPE_BASE);
        
    if(png_color_type == PNG_COLOR_TYPE_PALETTE)
    {
        for(i = 0; i < num_colors; i++)
        {
            pixel = palette[i];
            png_palette[i].red = ((uint8_t *)&pixel)[0];
            png_palette[i].green = ((uint8_t *)&pixel)[1];
            png_palette[i].blue = ((uint8_t *)&pixel)[2];
            png_trans[i] = ((uint8_t *)&pixel)[3];
        }
        
        png_set_PLTE(png, info, png_palette, num_colors);
        
        if(use_alpha)
            png_set_tRNS(png, info, png_trans, num_colors, NULL);
    }
    
    png_write_info(png, info);
    
    for(y = 0, image_ptr = image; y < height; y++, image_ptr += width)
    {
        switch(png_color_type)
        {
        case PNG_COLOR_TYPE_RGB_ALPHA:
            png_write_row(png, (png_bytep)image_ptr);
            break;
        case PNG_COLOR_TYPE_RGB:
            for(x = 0; x < width; x++)
            {
                pixel = image_ptr[x];
                row_buffer[3*x + 0] = ((uint8_t *)&pixel)[0];
                row_buffer[3*x + 1] = ((uint8_t *)&pixel)[1];
                row_buffer[3*x + 2] = ((uint8_t *)&pixel)[2];
            }
            
            png_write_row(png, (png_bytep)row_buffer);
            break;
        case PNG_COLOR_TYPE_GRAY:
            for(x = 0; x < width; x++)
            {
                pixel = image_ptr[x];
                row_buffer[x] = ((uint8_t *)&pixel)[0];
            }
            
            png_write_row(png, (png_bytep)row_buffer);
            break;
        case PNG_COLOR_TYPE_PALETTE:
            for(x = 0; x < width; x++)
            {
                pixel = image_ptr[x];
                
                for(i = 0; i < num_colors; i++)
                    if(pixel == palette[i])
                        break;
                                    
                row_buffer[x] = i;
            }
            
            png_write_row(png, (png_bytep)row_buffer);
            break;
        }
    }

    png_write_end(png, info);
    success = 1;
fail:
    if(palette)
        Free(palette);
    png_destroy_write_struct(&png, &info);
    Free(row_buffer);
    return success;
}
#endif /* USE_LIBPNG */


#ifdef USE_LIBTIFF

static int read_tiff(uint32_t **image, int *width, int *height,
    const char *filename, unsigned directory)
{
    TIFF *tiff;
    uint32 image_width, image_height;

    *image = 0;
    *width = *height = 0;
    
    if(!(tiff = TIFFOpen(filename, "r")))
    {
        ErrorMessage("TIFFOpen failed to open file.\n");
        return 0;
    }
    
    TIFFSetDirectory(tiff, directory);
    TIFFGetField(tiff, TIFFTAG_IMAGEWIDTH, &image_width);
    TIFFGetField(tiff, TIFFTAG_IMAGELENGTH, &image_height);
    *width = (int)image_width;
    *height = (int)image_height;
    
    if(*width > MAX_IMAGE_SIZE || *height > MAX_IMAGE_SIZE)
    {
        ErrorMessage("Image dimensions exceed MAX_IMAGE_SIZE.\n");
        goto fail;
    }
    
    if(!(*image = (uint32_t *)Malloc(
        sizeof(uint32_t)*image_width*image_height)))
        goto fail;
    
    if(!TIFFReadRGBAImageOriented(tiff, image_width, image_height,
        (uint32 *)*image, ORIENTATION_TOPLEFT, 1))
        goto fail;
    
    TIFFClose(tiff);
    return 1;
    
fail:
    if(*image)
        Free(*image);
    
    *width = *height = 0;
    TIFFClose(tiff);
    return 0;
}


static int write_tiff(const uint32_t *image, int width, int height,
    const char *filename)
{
    TIFF *tiff;
    uint16 alpha = EXTRASAMPLE_ASSOCALPHA;

    if(!image)
        return 0;
    
    if(!(tiff = TIFFOpen(filename, "w")))
    {
        ErrorMessage("TIFFOpen failed to open file.\n");
        return 0;
    }
    
    if(TIFFSetField(tiff, TIFFTAG_IMAGEWIDTH, width) != 1
        || TIFFSetField(tiff, TIFFTAG_IMAGELENGTH, height) != 1
        || TIFFSetField(tiff, TIFFTAG_SAMPLESPERPIXEL, 4) != 1
        || TIFFSetField(tiff, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_RGB) != 1
        || TIFFSetField(tiff, TIFFTAG_EXTRASAMPLES, 1, &alpha) != 1
        || TIFFSetField(tiff, TIFFTAG_BITSPERSAMPLE, 8) != 1
        || TIFFSetField(tiff, TIFFTAG_ORIENTATION, ORIENTATION_TOPLEFT) != 1
        || TIFFSetField(tiff, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG) != 1
        /* compression can be COMPRESSION_NONE, COMPRESSION_DEFLATE,
        COMPRESSION_LZW, or COMPRESSION_JPEG */
        || TIFFSetField(tiff, TIFFTAG_COMPRESSION, COMPRESSION_LZW) != 1)
    {
        ErrorMessage("TIFFSetField failed.\n");
        TIFFClose(tiff);
        return 0;
    }
    
    if(TIFFWriteEncodedStrip(tiff, 0, (tdata_t)image,
        4*((size_t)width)*((size_t)height)) < 0)
    {
        ErrorMessage("Error writing data to file.\n");
        TIFFClose(tiff);
        return 0;
    }

    TIFFClose(tiff);
    return 1;
}
#endif /* USE_LIBTIFF */


static void *convert_to_format(uint32_t *src, int width, int height,
    unsigned format)
{
    const int num_pixels = width*height;
    const int num_channels = (format & IMAGEIO_GRAYSCALE) ?
        1 : ((format & IMAGEIO_STRIP_ALPHA) ? 3 : 4);
    const int channel_stride = (format & IMAGEIO_PLANAR) ? num_pixels : 1;
    const int channel_stride2 = 2*channel_stride;
    const int channel_stride3 = 3*channel_stride;
    double *dest_f64;
    float *dest_f32;
    uint8_t *dest_u8;
    uint32_t pixel;
    int order[4] = {0, 1, 2, 3};
    int i, x, y, pixel_stride, row_stride;
    
    
    pixel_stride = (format & IMAGEIO_PLANAR) ? 1 : num_channels;
    
    if(format & IMAGEIO_COLUMNMAJOR)
    {
        row_stride = pixel_stride;
        pixel_stride *= height;
    }
    else
        row_stride = width*pixel_stride;
    
    if(format & IMAGEIO_BGRFLIP)
    {
        order[0] = 2;
        order[2] = 0;
    }
    
    if((format & IMAGEIO_AFLIP) && !(format & IMAGEIO_STRIP_ALPHA))
    {
        order[3] = order[2];
        order[2] = order[1];
        order[1] = order[0];
        order[0] = 3;
    }
    
    switch(format & (IMAGEIO_U8 | IMAGEIO_SINGLE | IMAGEIO_DOUBLE))
    {
    case IMAGEIO_U8:  /* destination type is uint8_t */
        if(!(dest_u8  = (uint8_t *)Malloc(
            sizeof(uint8_t)*num_channels*num_pixels)))
            return NULL;
        
        switch(num_channels)
        {
        case 1: /* Convert RGBA U8 to grayscale U8 */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_u8[i] = (uint8_t)(0.299f*((uint8_t *)&pixel)[0]
                        + 0.587f*((uint8_t *)&pixel)[1]
                        + 0.114f*((uint8_t *)&pixel)[2] + 0.5f);
                }
            break;
        case 3: /* Convert RGBA U8 to RGB (or BGR) U8 */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_u8[i] =
                        ((uint8_t *)&pixel)[order[0]];
                    dest_u8[i + channel_stride] =
                        ((uint8_t *)&pixel)[order[1]];
                    dest_u8[i + channel_stride2] =
                        ((uint8_t *)&pixel)[order[2]];
                }
            break;
        case 4: /* Convert RGBA U8 to RGBA (or BGRA, ARGB, or ABGR) U8 */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_u8[i] =
                        ((uint8_t *)&pixel)[order[0]];
                    dest_u8[i + channel_stride] =
                        ((uint8_t *)&pixel)[order[1]];
                    dest_u8[i + channel_stride2] =
                        ((uint8_t *)&pixel)[order[2]];
                    dest_u8[i + channel_stride3] =
                        ((uint8_t *)&pixel)[order[3]];
                }
            break;
        }
        return dest_u8;
    case IMAGEIO_SINGLE:  /* destination type is float */
        if(!(dest_f32 = (float *)Malloc(
            sizeof(float)*num_channels*num_pixels)))
            return NULL;
        
        switch(num_channels)
        {
        case 1: /* Convert RGBA U8 to grayscale float */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_f32[i] =
                        1.172549019607843070675535e-3f*((uint8_t *)&pixel)[0]
                        + 2.301960784313725357840079e-3f*((uint8_t *)&pixel)[1]
                        + 4.470588235294117808150007e-4f*((uint8_t *)&pixel)[2];
                }
            break;
        case 3: /* Convert RGBA U8 to RGB (or BGR) float */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y; x < width;
                    x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_f32[i] =
                        ((uint8_t *)&pixel)[order[0]]/255.0f;
                    dest_f32[i + channel_stride] =
                        ((uint8_t *)&pixel)[order[1]]/255.0f;
                    dest_f32[i + channel_stride2] =
                        ((uint8_t *)&pixel)[order[2]]/255.0f;
                }
            break;
        case 4: /* Convert RGBA U8 to RGBA (or BGRA, ARGB, or ABGR) float */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_f32[i] =
                        ((uint8_t *)&pixel)[order[0]]/255.0f;
                    dest_f32[i + channel_stride] =
                        ((uint8_t *)&pixel)[order[1]]/255.0f;
                    dest_f32[i + channel_stride2] =
                        ((uint8_t *)&pixel)[order[2]]/255.0f;
                    dest_f32[i + channel_stride3] =
                        ((uint8_t *)&pixel)[order[3]]/255.0f;
                }
            break;
        }
        return dest_f32;
    case IMAGEIO_DOUBLE:  /* destination type is double */
        if(!(dest_f64 = (double *)Malloc(sizeof(double)*num_channels*num_pixels)))
            return NULL;
        
        switch(num_channels)
        {
        case 1: /* Convert RGBA U8 to grayscale double */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y; x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_f64[i] = 1.172549019607843070675535e-3*((uint8_t *)&pixel)[0]
                        + 2.301960784313725357840079e-3*((uint8_t *)&pixel)[1]
                        + 4.470588235294117808150007e-4*((uint8_t *)&pixel)[2];
                }
            break;
        case 3: /* Convert RGBA U8 to RGB (or BGR) double */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y; x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_f64[i] =
                        ((uint8_t *)&pixel)[order[0]]/255.0;
                    dest_f64[i + channel_stride] =
                        ((uint8_t *)&pixel)[order[1]]/255.0;
                    dest_f64[i + channel_stride2] =
                        ((uint8_t *)&pixel)[order[2]]/255.0;
                }
            break;
        case 4: /* Convert RGBA U8 to RGBA (or BGRA, ARGB, or ABGR) double */
            for(y = 0; y < height; y++, src += width)
                for(x = 0, i = row_stride*y; x < width; x++, i += pixel_stride)
                {
                    pixel = src[x];
                    dest_f64[i] =
                        ((uint8_t *)&pixel)[order[0]]/255.0;
                    dest_f64[i + channel_stride] =
                        ((uint8_t *)&pixel)[order[1]]/255.0;
                    dest_f64[i + channel_stride2] =
                        ((uint8_t *)&pixel)[order[2]]/255.0;
                    dest_f64[i + channel_stride3] =
                        ((uint8_t *)&pixel)[order[3]]/255.0;
                }
            break;
        }
        return dest_f64;
    default:
        return NULL;
    }
}


static uint32_t *convert_from_format(void *src, int width, int height,
    unsigned format)
{
    const int num_pixels = width*height;
    const int num_channels = (format & IMAGEIO_GRAYSCALE) ?
        1 : ((format & IMAGEIO_STRIP_ALPHA) ? 3 : 4);
    const int channel_stride = (format & IMAGEIO_PLANAR) ? num_pixels : 1;
    const int channel_stride2 = 2*channel_stride;
    const int channel_stride3 = 3*channel_stride;
    double *src_f64 = (double *)src;
    float *src_f32 = (float *)src;
    uint8_t *src_u8 = (uint8_t *)src;
    uint8_t *dest, *dest_ptr;
    int order[4] = {0, 1, 2, 3};
    int i, x, y, pixel_stride, row_stride;
    
    if(!(dest = (uint8_t *)Malloc(sizeof(uint32_t)*num_pixels)))
        return NULL;
    
    dest_ptr = dest;
    pixel_stride = (format & IMAGEIO_PLANAR) ? 1 : num_channels;
    
    if(format & IMAGEIO_COLUMNMAJOR)
    {
        row_stride = pixel_stride;
        pixel_stride *= height;
    }
    else
        row_stride = width*pixel_stride;
    
    if(format & IMAGEIO_BGRFLIP)
    {
        order[0] = 2;
        order[2] = 0;
    }
    
    if((format & IMAGEIO_AFLIP) && !(format & IMAGEIO_STRIP_ALPHA))
    {
        order[3] = order[2];
        order[2] = order[1];
        order[1] = order[0];
        order[0] = 3;
    }
    
    switch(format & (IMAGEIO_U8 | IMAGEIO_SINGLE | IMAGEIO_DOUBLE))
    {
    case IMAGEIO_U8:  /* Source type is uint8_t */
        switch(num_channels)
        {
        case 1: /* Convert grayscale U8 to RGBA U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x] =
                    dest_ptr[4*x + 1] =
                    dest_ptr[4*x + 2] = src_u8[i];
                    dest_ptr[4*x + 3] = 255;
                }
            break;
        case 3: /* Convert RGB (or BGR) U8 to RGBA U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x + order[0]] = src_u8[i];
                    dest_ptr[4*x + order[1]] = src_u8[i + channel_stride];
                    dest_ptr[4*x + order[2]] = src_u8[i + channel_stride2];
                    dest_ptr[4*x + 3] = 255;
                }
            break;
        case 4: /* Convert RGBA U8 to RGBA (or BGRA, ARGB, or ABGR) U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x + order[0]] = src_u8[i];
                    dest_ptr[4*x + order[1]] = src_u8[i + channel_stride];
                    dest_ptr[4*x + order[2]] = src_u8[i + channel_stride2];
                    dest_ptr[4*x + order[3]] = src_u8[i + channel_stride3];
                }
            break;
        }
        break;
    case IMAGEIO_SINGLE:  /* Source type is float */
        switch(num_channels)
        {
        case 1: /* Convert grayscale float to RGBA U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x] =
                    dest_ptr[4*x + 1] =
                    dest_ptr[4*x + 2] = ROUNDCLAMPF(src_f32[i]);
                    dest_ptr[4*x + 3] = 255;
                }
            break;
        case 3: /* Convert RGBA U8 to RGB (or BGR) float */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x + order[0]] =
                        ROUNDCLAMPF(src_f32[i]);
                    dest_ptr[4*x + order[1]] =
                        ROUNDCLAMPF(src_f32[i + channel_stride]);
                    dest_ptr[4*x + order[2]] =
                        ROUNDCLAMPF(src_f32[i + channel_stride2]);
                    dest_ptr[4*x + 3] = 255;
                }
            break;
        case 4: /* Convert RGBA U8 to RGBA (or BGRA, ARGB, or ABGR) float */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x + order[0]] =
                        ROUNDCLAMPF(src_f32[i]);
                    dest_ptr[4*x + order[1]] =
                        ROUNDCLAMPF(src_f32[i + channel_stride]);
                    dest_ptr[4*x + order[2]] =
                        ROUNDCLAMPF(src_f32[i + channel_stride2]);
                    dest_ptr[4*x + order[3]] =
                        ROUNDCLAMPF(src_f32[i + channel_stride3]);
                }
            break;
        }
        break;
    case IMAGEIO_DOUBLE:  /* Source type is double */
        switch(num_channels)
        {
        case 1: /* Convert grayscale double to RGBA U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x] =
                    dest_ptr[4*x + 1] =
                    dest_ptr[4*x + 2] = ROUNDCLAMP(src_f64[i]);
                    dest_ptr[4*x + 3] = 255;
                }
            break;
        case 3: /* Convert RGB (or BGR) double to RGBA U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y;
                    x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x + order[0]] =
                        ROUNDCLAMP(src_f64[i]);
                    dest_ptr[4*x + order[1]] =
                        ROUNDCLAMP(src_f64[i + channel_stride]);
                    dest_ptr[4*x + order[2]] =
                        ROUNDCLAMP(src_f64[i + channel_stride2]);
                    dest_ptr[4*x + 3] = 255;;
                }
            break;
        case 4: /* Convert RGBA (or BGRA, ARGB, or ABGR) double to RGBA U8 */
            for(y = 0; y < height; y++, dest_ptr += 4*width)
                for(x = 0, i = row_stride*y; x < width; x++, i += pixel_stride)
                {
                    dest_ptr[4*x + order[0]] =
                        ROUNDCLAMP(src_f64[i]);
                    dest_ptr[4*x + order[1]] =
                        ROUNDCLAMP(src_f64[i + channel_stride]);
                    dest_ptr[4*x + order[2]] =
                        ROUNDCLAMP(src_f64[i + channel_stride2]);
                    dest_ptr[4*x + order[3]] =
                        ROUNDCLAMP(src_f64[i + channel_stride3]);
                }
            break;
        }
        break;
    default:
        return NULL;
    }
    
    return (uint32_t *)dest;
}


int identify_image_type(char *type, const char *filename)
{
    FILE *file;
    uint32_t magic;
    
    type[0] = '\0';
    
    if(!(file = fopen(filename, "rb")))
        return 0;
    
    /* Determine the file format by reading the first 4 bytes */
    magic = ((uint32_t)getc(file));
    magic |= ((uint32_t)getc(file)) << 8;
    magic |= ((uint32_t)getc(file)) << 16;
    magic |= ((uint32_t)getc(file)) << 24;
    
    /* Test for errors */
    if(ferror(file))
    {
        fclose(file);
        return 0;
    }
    
    fclose(file);
    
    if((magic & 0x0000FFFFL) == 0x00004D42L)                /* BMP */
        strcpy(type, "BMP");
    else if((magic & 0x00FFFFFFL) == 0x00FFD8FFL)           /* JPEG/JFIF */
        strcpy(type, "JPEG");
    else if(magic == 0x474E5089L)                           /* PNG */
        strcpy(type, "PNG");
    else if(magic == 0x002A4949L || magic == 0x2A004D4DL)   /* TIFF */
        strcpy(type, "TIFF");
    else if(magic == 0x38464947L)                           /* GIF */
        strcpy(type, "GIF");
    else if(magic == 0x474E4D8AL)                           /* MNG */
        strcpy(type, "MNG");
    else if((magic & 0xF0FF00FFL) == 0x0001000AL            /* PCX */
        && ((magic >> 8) & 0xFF) < 6)
        strcpy(type, "PCX");
    else
        return 0;
    
    return 1;
}


void *read_image(int *width, int *height,
    const char *filename, unsigned format)
{
    void *image = NULL;
    uint32_t *image_u8 = NULL;
    FILE *file;
    char type[8];
    
    
    identify_image_type(type, filename);
    
    if(!(file = fopen(filename, "rb")))
    {
        ErrorMessage("Unable to open file \"%s\".\n", filename);
        return 0;
    }
    
    if(!strcmp(type, "BMP"))
    {
        if(!read_bmp(&image_u8, width, height, file))
            ErrorMessage("Failed to read \"%s\".\n", filename);
    }
    else if(!strcmp(type, "JPEG"))
    {
#ifdef USE_LIBJPEG
        if(!(read_jpeg(&image_u8, width, height, file)))
            ErrorMessage("Failed to read \"%s\".\n", filename);
#else
        ErrorMessage("file \"%s\" is a JPEG image.\n"
                     "Compile with USE_LIBJPEG to enable JPEG reading.\n",
                     filename);
#endif
    }
    else if(!strcmp(type, "PNG"))
    {
#ifdef USE_LIBPNG
        if(!(read_png(&image_u8, width, height, file)))
            ErrorMessage("Failed to read \"%s\".\n", filename);
#else
        ErrorMessage("file \"%s\" is a PNG image.\n"
                     "Compile with USE_LIBPNG to enable PNG reading.\n",
                     filename);
#endif
    }
    else if(!strcmp(type, "TIFF"))
    {
#ifdef USE_LIBTIFF
        fclose(file);
        
        if(!(read_tiff(&image_u8, width, height, filename, 0)))
            ErrorMessage("Failed to read \"%s\".\n", filename);
        
        file = NULL;
#else
        ErrorMessage("file \"%s\" is a TIFF image.\n"
                     "Compile with USE_LIBTIFF to enable TIFF reading.\n",
                     filename);
#endif
    }
    else
    {
        /* file format is unsupported. */
        if(type[0])
            ErrorMessage("file \"%s\" is a %s image.", filename, type);
        else
            ErrorMessage("file \"%s\" is an unrecognized format.", filename);
        fprintf(stderr, "\nSorry, only "
            READIMAGE_FORMATS_SUPPORTED " reading is supported.\n");
    }
    
    if(file)
        fclose(file);
    
    if(image_u8 && format)
    {
        image = convert_to_format(image_u8, *width, *height, format);
        Free(image_u8);
    }
    else
        image = image_u8;
    
    return image;
}


int write_image(void *image, int width, int height,
    const char *filename, unsigned format, int quality)
{
    FILE *file;
    uint32_t *image_u8;
    enum {BMP_FORMAT, JPEG_FORMAT, PNG_FORMAT, TIFF_FORMAT} fileformat;
    int success = 0;
    
    if(!image || width <= 0 || height <= 0)
    {
        ErrorMessage("Null image.\n");
        ErrorMessage("Failed to write \"%s\".\n", filename);
        return 0;
    }
    
    if(string_ends_with(filename, ".bmp"))
        fileformat = BMP_FORMAT;
    else if(string_ends_with(filename, ".jpg")
        || string_ends_with(filename, ".jpeg"))
    {
        fileformat = JPEG_FORMAT;
#ifndef USE_LIBJPEG
        ErrorMessage("Failed to write \"%s\".\n", filename);
        ErrorMessage("Compile with USE_LIBJPEG to enable JPEG writing.\n");
        return 0;
#endif
    }
    else if(string_ends_with(filename, ".png"))
    {
        fileformat = PNG_FORMAT;
#ifndef USE_LIBPNG
        ErrorMessage("Failed to write \"%s\".\n", filename);
        ErrorMessage("Compile with USE_LIBPNG to enable PNG writing.\n");
        return 0;
#endif
    }
    else if(string_ends_with(filename, ".tif")
        || string_ends_with(filename, ".tiff"))
    {
        fileformat = TIFF_FORMAT;
#ifndef USE_LIBTIFF
        ErrorMessage("Failed to write \"%s\".\n", filename);
        ErrorMessage("Compile with USE_LIBTIFF to enable TIFF writing.\n");
        return 0;
#endif
    }
    else
    {
        ErrorMessage("Failed to write \"%s\".\n", filename);
        
        if(string_ends_with(filename, ".gif"))
            ErrorMessage("GIF is not supported.  ");
        else if(string_ends_with(filename, ".mng"))
            ErrorMessage("MNG is not supported.  ");
        else if(string_ends_with(filename, ".pcx"))
            ErrorMessage("PCX is not supported.  ");
        else
            ErrorMessage("Unable to determine format from extension.\n");
        
        ErrorMessage("Sorry, only "
            WRITEIMAGE_FORMATS_SUPPORTED " writing is supported.\n");
        return 0;
    }
    
    if(!(file = fopen(filename, "wb")))
    {
        ErrorMessage("Unable to write to file \"%s\".\n", filename);
        return 0;
    }
    
    if(!(image_u8 = convert_from_format(image, width, height, format)))
        return 0;
    
    switch(fileformat)
    {
    case BMP_FORMAT:
        success = write_bmp(image_u8, width, height, file);
        break;
    case JPEG_FORMAT:
#ifdef USE_LIBJPEG
        success = write_jpeg(image_u8, width, height, file, quality);
#else
        /* Dummy operation to avoid unused variable warning if compiled
        without libjpeg.  Note that execution returns above if
        format == JPEG_FORMAT and USE_LIBJPEG is undefined. */
        success = quality;
#endif
        break;
    case PNG_FORMAT:
#ifdef USE_LIBPNG
        success = write_png(image_u8, width, height, file);
#endif
        break;
    case TIFF_FORMAT:
#ifdef USE_LIBTIFF
        fclose(file);
        success = write_tiff(image_u8, width, height, filename);
        file = 0;
#endif
        break;
    }
    
    if(!success)
        ErrorMessage("Failed to write \"%s\".\n", filename);
    
    Free(image_u8);
    
    if(file)
        fclose(file);
    
    return success;
}