---

mosaic_byte_export.c

Go to the documentation of this file.

#include        <stdio.h>
#include        "sadie.h"
#include        "mosaic.h"

/*-Copyright Information------------------------------------------------------*/
/* Copyright (c) 1999 by the University of Arizona Digital Image Analysis Lab */
/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function writes a TREES mosaic  image from main memory to disk      */
/*   in BIP format.                                                           */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void MOSAIC_BYTE_EXPORT (in,name)
IMAGE_BYTE *in;     /*  I   Pointer to the image.                             */
char *name;         /*  I   String, containing the name of the disk file.     */
/*----------------------------------------------------------------------------*/
{   
    static char array[2] = { 0x00, 0x01 };
    char msg[SLEN];
    short nbnd, nbit;
    unsigned short nlin, npix;
    unsigned short j,k;
    long nbytes;
    PIXEL_BYTE *column=NULL;
    PIXEL gmin, gmax;
    FILE *fp=NULL;

    /* if (TIMES) TIMING(T_MOSAIC_EXPORT); */

    /* check input */

    if (!CHECKIMG_BYTE(in)) {
        MESSAGE('E'," Can't identify image.");
        goto the_end;
    }

    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"MOSAIC_BYTE_EXPORT");
        MESSAGE('I',"");
        sprintf(msg," Input image:               %s",in->text);
        MESSAGE('I',msg);
        sprintf(msg," Output disk file:          %s",name);
        MESSAGE('I',msg);
        sprintf(msg," Output format:             8-bit BIP in MOSAIC_IMAGE format");
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
    }


    if (in->nbnd > 1) {
        MESSAGE('I'," Only first band of multi-band image will be written to disk.");
    }


    in->nbit = 8;

    /* swap bytes, if necessary */
    nbnd = 1;
    nlin = in->nlin;
    npix = in->npix;
    nbit = in->nbit;
    if (*(short *)array != 1) {

        //printf("MOSAIC_EXPORT: Swapping bytes...\n"); 

    
        nbnd = SWAP((unsigned char *)&nbnd);
        nlin = SWAP((unsigned char *)&nlin);
        npix = SWAP((unsigned char *)&npix);
        nbit = SWAP((unsigned char *)&nbit);
    } 
    //else {
    //    //printf("MOSAIC_EXPORT: No need to swap bytes...\n");
    //}

    printf("Exporting mosaic of size %dx%d\n",in->npix,in->nlin);

    /* open image file */
    if (!(fp=fopen(name,"w"))) {
        sprintf(msg," Can't open file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }



    /*
     * write image header to file 
     */

    /* number of bands */
    if (fwrite(&nbnd,sizeof(nbnd),1,fp) != 1) {
        sprintf(msg," Can't write number of bands to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }

    /* number of lines/band */
    if (fwrite(&nlin,sizeof(nlin),1,fp) != 1) {
        sprintf(msg," Can't write number of lines/band to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }

    /* number of PIXEL_BYTEs/line */
    if (fwrite(&npix,sizeof(npix),1,fp) != 1) {
        sprintf(msg," Can't write number of PIXEL_BYTEs/line to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }

    /* number of bits/PIXEL_BYTE */
    if (fwrite(&nbit,sizeof(nbit),1,fp) != 1) {
        sprintf(msg," Can't write number of bits/PIXEL_BYTE to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }

    gmin = (PIXEL)in->gmin;
    if (fwrite(&gmin,sizeof(gmin),1,fp) != 1) {
        sprintf(msg," Can't write gmin to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }
    gmax = (PIXEL)in->gmax;
    if (fwrite(&gmax,sizeof(gmax),1,fp) != 1) {
        sprintf(msg," Can't write gmax to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }

    /* rest of header */
    nbytes = (long)(sizeof(in->pstr)+sizeof(in->text));
    if (fwrite(&(in->pstr),nbytes,1,fp) != 1) {
        sprintf(msg," Can't write image header to file %s.",name);
        MESSAGE('E',msg);
        goto the_end;
    }



    /* write the image data to file */
    nbytes = (long)in->nlin;
    if ((column=(PIXEL_BYTE *)malloc((long)in->nlin * sizeof(PIXEL_BYTE))) == NULL) {
        MESSAGE('E'," Memory request failed.");
        goto the_end;
    }
    
    for (k = 0; k < in->npix; k++) {
        for (j=0; j < in->nlin; j++) {
            column[j] = in->data[0][j][k];
        }
        
        if (fwrite(column,nbytes,1,fp) != 1) {
            sprintf(msg," Can't write image data to file %s.",name);
            MESSAGE('E',msg);
        }
    }


    /* close image file */
    the_end:

    if (column) free(column);
    
    if (fp) {
        if (fclose(fp)) {
            sprintf(msg," Can't close file %s.",name);
            MESSAGE('W',msg);
        }
    }

    /* if (TIMES) TIMING(T_EXIT); */

}

---